We’re smarter together. Learn from this collection of community knowledge and add your expertise.

North Canterbury, New Zealand Earthquake

by Respected Advisor on ‎11-30-2016 12:59 AM - edited on ‎12-01-2016 09:25 AM by Community Manager (1,182 Views)

On 14 November 2016 at just after midnight local time a major earthquake struck in the North Canterbury region of New Zealand. This isn’t surprising as the country sits on the boundary of the Australian and Pacific tectonic plates, part of the Earth’s crust, and it has a long history of earthquake activity.

 

Some of the impacts of this 7.8 magnitude quake make amazing reading. Scientists have estimated that the two biggest land masses of New Zealand, the North and South Islands moved closer together by 2 metres!  Around the other side of the world the Large Hadron Collider (LHC) picked up the earthquake waves as they distorted the 27km in circumference particle accelerator ring by 0.1mm:

 

LHC.jpg

 

The LHC is undoubtedly the world’s most expensive (and unintended) earthquake detector. However they need to adjust the alignment of the sub-atomic particle beams even for movements like this.

 

So why am I talking about earthquakes? Well, inspired in part by @DarthPathos I wanted to experiment with using SAS to analyse earthquake data. Fortunately there is a great website run by Geonet, a public organisation in New Zealand that freely provides this data: http://www.geonet.org.nz/

 

I extracted earthquake data into a CSV file starting from 5 days prior to the big earthquake, to 5 days after. I then tried the SAS SGPLOT procedure to get this:

 

North Canterbury New Zealand Earthquake: Magnitude 7.8 Depth 15km

 

Quake1.jpg

 

As you can see the magnitude 7.8 earthquake stands out like a sore thumb. Leading up to the big event there are very few quakes happening and those that do happen are quite small. You wouldn’t normally feel these unless they were very close by.

However contrast that with the quake storm that follows. There are literally hundreds of them, with dozens exceeding magnitude 4 or 5. These you would definitely feel even if you were some distance away.

 

I then decided to enhance the SGPLOT to tell more of a story.

 

North Canterbury New Zealand Earthquake: Magnitude 7.8 Depth 15km

 

Quake2.jpg

 

I used the DROPLINE statement to identify the big shake and I’ve overlaid the depths of the quakes using the COLORRESPONSE option. It is interesting to note that most of the aftershocks are quite shallow – less than 30km. Generally speaking the shallower the quake, the more violent the shaking so the quake depth is an important measure. What is also interesting about the aftershocks is the trend in the magnitudes. Initially the magnitudes are high but they decay over the days and weeks following in an exponential trend.

 

Here is the complete program for this graph:

 

 

proc sgplot data= WORK.EARTHQUAKES ;
  title "North Canterbury New Zealand Earthquake: Magnitude 7.8 Depth 15km";
  where depth <=50;
  scatter x = origin_dt y = magnitude / group = Quakegroup colorresponse = depth colormodel = (red green blue black);
  dropline x = "13NOV2016:11:02:56"dt y = 7.82 / dropto = both lineattrs=(color=red pattern=dot) ; 
run;

 

Fortunately the 7.8 earthquake happened in sparsely-populated area and at a time when most people were at home so there were only 2 deaths and a small number of injuries. Here is a map from SAS VA that shows the area. This was very quick and easy to do. I loaded the earthquake dataset into the VA LASR library. I then used the earthquake ID to define a custom geographic location using latitude and longitude. I then created a Geo Bubble Map report and defined the size of the dot to show the earthquake magnitude and the colour to show the depth. This screenshot really doesn't do justice to the amazing level of detail. I've also only shown the magnitude 5 or higher quakes otherwise the map would just be obliterated with dots.

 

Quake3.jpg

 

On a personal note, I live in the city of Wellington which is about 200km from the epicentre of the 7.8 earthquake. There was significant damage in the city even at that distance including a lot of buildings that are now having to be demolished or to have major repairs. The building I work in needs some minor repair work so I could be working offsite for a month or two! Meanwhile the aftershocks are continuing and we are still feeling them on most days.

 

So there you have it. A real life experience that motivated me to use SAS to explain what living through earthquakes is like. 

Comments
by Super User
on ‎11-30-2016 05:00 AM

Thanks for sharing your experiences and analysis @SASKiwi. Sounds like there was quite a bit of damage. Hopefully the Kiwi spirit remains strong. I wasn't aware the aftershocks are still happening. Keep safe!

by Respected Advisor
on ‎11-30-2016 02:21 PM

Thanks @MichelleHomes! So far so good. This is proving to be a good test for the resilience of our SAS systems and so far they are holding up well.

by Regular Contributor
on ‎11-30-2016 08:28 PM

I'm honoured that I inspired such an awesome use of Free Data!  @BeverlyBrown has been my source of inspiration, motivation and challenging me - glad it's infectious!

(Thanks @MichelleHomes for pointing me to this!)

Keep it up and looking forward to what else you come up with other free datasets 8-)

Chris

by Super User
on ‎11-30-2016 08:43 PM

My pleasure @DarthPathos!

 

As @SASKiwi had referenced you I thought I'd let you know. It seems member tagging in Article documents isn't possible - maybe one of the administrators can look into this.

by Respected Advisor
on ‎12-01-2016 03:12 AM

Hey thanks @DarthPathos! I've been following your articles with interest. It was excellent to discover such a great source of free data in my own backyard that is so topical with recent local events.

by Regular Contributor
on ‎12-01-2016 07:40 AM

@SASKiwi Wow, thank you!  I always love finding new datasets.  This week's article is another "fun" one, but then I am going to get back to "serious" stuff next week - found some interesting datasets from Africa I'd like to use.

 

Looking forward to whatever else you find!

Chris

by Community Manager
on ‎12-01-2016 07:52 AM

@MichelleHomes, et al, -- at-mentions are not supported in the library articles.  I dunno why, really, but it's by design.  I think we're more than making up for it by using them in the comments though!

by Super User
on ‎12-01-2016 09:12 AM
by Frequent Contributor
on ‎12-01-2016 07:50 PM

Thanks, @SASKiwi, not only a great teaching example, but also a comfort for all the other New Zealand based SAS users in the quake-affected zone who browse the community postings to see for themselves that the magnitude of the aftershocks generally decrease, albeit exponentially. As seen by other New Zealanders, Wellingtonians seem to take quakes in their stride, and almost downplay the many smaller quakes that Wellington experiences, but this set was pretty severe.  Are you one of those, with nerves of steel?

 

 

by Respected Advisor
on ‎12-02-2016 04:16 PM

Nerves of steel? Perhaps not so much @Damien_Mather! Thanks for the kind words. It keeps me motivated to continue the contributions to the community. BTW I know you from the early days of SAS conferences in NZ. I always found your presentations highly entertaining as well as informing. Great to see you are still using SAS!

by Regular Contributor
on ‎12-12-2016 02:14 AM

I'm with @SASKiwi - there's no nerves of steel in Wellington. There was another 4.8 magnitude shock this evening, about 50km away. Fortunately where we are, we often don't feel them, or very much (although the 7.8 monster scared the **bleep** out of us). But many others in Wellington were nervous again. And this is nothing compared with the poor sods in Christchurch who went through this for years, and are still rudely reminded of what they went through.

 

You never get blasé about earthquakes.

Contributors
Your turn
Sign In!

Want to write an article? Sign in with your profile.