Category Archives: Earthquakes

How and why does COMET respond to earthquakes?

A personal view from COMET* director, Prof Tim Wright 

14 February 2023 

We have all been shocked and horrified by the pictures coming from Türkiye and Syria over the past week following the two large earthquakes that happened on 6 February. COMET aims to generate and freely distribute robust scientific satellite deformation data sets following all moderate-large earthquakes on land. These can assist the humanitarian and scientific response.  

Here I want to briefly explain COMET’s approach to responding to an event like this from my perspective not just as COMET director, but also as a scientist who has been involved in responding to many earthquakes using satellite deformation data and who has been working on active faulting in Türkiye for my entire career. 

Deformation data from satellites are now a fundamental data set for understanding what happened in an earthquake, complementing the information that is available from seismology and from field studies. Deformation data shows, very precisely, which faults moved in an earthquake and by how much. Earthquakes are not point sources – they involve slip on faults. As earthquakes get larger, they involve larger/longer faults with increasing amounts of slip. Understanding how much and where slip occurred on a fault helps us understand the amount of shaking that is likely to have been experienced in different locations – distance to the fault is a primary control on the intensity of shaking. 

In COMET, we have been developing a service to automatically produce ground movement data sets from radar satellites following all earthquakes bigger than Magnitude 5.5 that are likely to impact the land surface, and are now a trusted source for these data sets. This is possible because the Sentinel-1 satellites, part of the European Commission’s Copernicus programme, have a systematic acquisition strategy over tectonic and volcanic areas, something that COMET worked with the European Space Agency to help define. We produce results from satellite data as fast as possible following the satellite acquisition, make them available to the public via a dedicated portal, and typically let people know about them using the @NERC_COMET and @COMET_database accounts on twitter.  

In the case of the Türkiye earthquakes, the data sets that we produced in COMET, which was also analysed by other international groups, show that the fault that ruptured in the initial Magnitude 7.8 earthquake was a 300 km long section of the East Anatolian Fault, and the fault that ruptured in the second large earthquake, a magnitude 7.5, was over 100 km long and occurred along a different branch of the same fault system. In UK terms, the length of the first fault is about the same as the distance from Bristol to Hull; and is close to the distance from London to Paris. This great extent is one of the reasons why the devastation has been so widespread and horrendous. Both earthquakes occurred on mapped faults and in areas where national maps of seismic hazard required the strictest building codes. 

We believe it is important to analyse satellite data as quickly as possible and to be open and transparent with the results, sharing them with the wider community and allowing people to use them in whatever way they wish. For most earthquakes, we now have data within a few days of the event, and we post results within a few hours of the satellite data acquisition, whether there is any media interest in an event or not. We would like to be able to do this within a few hours of every event, and with new satellite systems from the Europe, the US/India, Japan, Canada, China, and others coming on stream, there is good prospect of being able to provide results within a day of most events in the second half of this decade. 

However, even if the results are not yet always available within the 72-hour window required for initial search and rescue efforts, the analysis of satellite deformation data can still be useful. Various groups around the world are using these types of data (provided by COMET or processed by other groups) to build models of the event – these can in turn be used to help understand how an earthquake influences the activity on nearby faults. More directly, the data, and information derived from it such as the location of individual fault ruptures, can help responders understand the potential impact on key infrastructure such as bridges and roads, which are vital for the relief effort but might have been impacted by the surface rupture, or other slope failures resulting from the earthquake. The data also help guide field geologists to sites in the field where they can study the earthquake fault rupture up close – such features can degrade very quickly in poor weather conditions; documenting them is important for understanding what palaeoseismologists see when they are investigating records of ancient earthquake ruptures and for understanding the detail of what happened in the causative earthquake. In turn, detailed surface rupture mapping can help to add constraints to models of the earthquake based on geodetic data. In the longer-term, forensically understanding every earthquake helps us prepare for future events in different regions.  

As well as providing data sets, COMET scientists, including myself, often discuss the meaning of these data sets with colleagues in open fora like twitter, and respond to queries from colleagues and the media. Leading on from comments by American seismologist and science communicator Dr Lucy Jones, my view also is that this information sharing is vital so that scientists involved directly in the response, or those involved in communication to the public, have access to the latest data and can understand what happened during an event. Often this understanding evolves quite quickly during the hours and days that follow the earthquake. Having reliable information communicated directly in platforms like twitter, and indirectly via the media, is vital for combatting misinformation and conspiracy theories that can unfortunately proliferate in the absence of reliable scientific comment. Twitter has become a key platform for many of those discussions. 

Behind the scenes, we also are passing on information to local and international partners directly, so that they have the data they need to respond, and we are responding to queries to help people use the data. We are also providing data and information to help assist with the UK government’s emergency response via our partners in the British Geological Survey. 

We are of course acutely conscious that our scientific responses to events like the recent earthquakes are happening in the context of humanitarian tragedies. If you are able, I encourage you to give generously to organisations such as the Disasters Emergency Committee in the UK, who are helping survivors in urgent need of aid and assistance. 

 *COMET is the UK Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics, a partnership between scientists in UK Universities and the British Geological Survey tasked with delivering National Capability science on behalf of the Natural Environment Research Council. Our focus is on using Earth Observation data alongside other data sets to help understand earthquake and volcanic processes and hazards. 

Thanks to Wendy Bohon, Ruth Lawford-Rolfe, and Laura Gregory for providing input on early drafts of this text. 

Türkiye-Syria Earthquakes, February 2023 

On 6 February, a 7.8-magnitude earthquake struck the East Anatolian Fault affecting large areas of Southern Türkiye and Northern Syria. This was followed by a 7.5-magnitude event approximately 9 hours later, around 60 miles to the north. To date more than 37,000 people are confirmed to have died, large numbers of people are affected across the region and the damage to buildings and infrastructure is significant.   

Images from ESA’s Sentinel-1A satellite captured on 9/10 February clearly showed the physical effects of the earthquake on the ground, including deformation of up to 6 metres along a 300km section of the fault, and the second event causing a second ~125km rupture. Many population centres sit close to these zones, explaining the significant human impact of the event.  

By combining Sentinel-1A imagery from before and after the earthquake, COMET scientists have been able to measure surface deformation that is clearly visible in InSAR and pixel offset tracking data sets shown below: 

In addition to the results from the satellite radar data, we have also used the pre- and post-event optical images from Sentinel-2 to estimate ground movement in the earthquakes also using pixel tracking: 

The processing outputs from Sentinel-1A data are available for download at our LiCSAR system event page. The results from Sentinel-2 are available here. 

The images above contain modified Copernicus Sentinel-1 and Sentinel-2 data analysed by the Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET). Data processing uses JASMIN, the UK’s collaborative data analysis environment (http://jasmin.ac.uk). 

For more information on why and how COMET responds to earthquake events read this article by our Director – Professor Tim Wright. 

COMET – 14 February 2023 

COMET Central Asia Fault Database

COMET researcher Tamarah King, based at the University of Oxford, has recently written a blog providing a research update on the COMET Central Asia Fault Database; progress report.

The COMET Central Asia Fault Database integrates decades of fault mapping and field-studies by researchers from the UK NERC Centre for Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET), and global collaborators.

Central Asia is home to one of the world’s great mountain ranges–the Tien Shan–which is formed by vigorous crustal convergence across a multitude of tectonic faults. Here we describe the motivation to assemble the database and the choices that we have made in its structure, which are based on utility, necessity, and limitations in available information. We are working towards a full public release of the dataset, so keep an eye out!

Key points:

  • COMET researchers have assembled a comprehensive database of active faults and associated attributes within Central Asia.
  • The database is comprised of structures identified by COMET researchers from both remote and field mapping (rather than a digitization of all published maps).
  • Faults are represented at three scales to suit various applications, e.g., geotechnical site exposure, geomorphic and neotectonic science, structural continuity for regional deformation models, and an inventory of seismogenic sources.
Image from field-work along the Dzhungarian Fault, Tien Shan, Kazakhstan (Credit: Austin Elliot)

Background Motivation

COMET researchers have been investigating active tectonic structures across the Central Asia region since the early 2000s through programs such as Earthquakes without Frontiers and Looking Inside the Continents from Space, along with local partner institutions in the region such as the Kyrgyz Institute of Seismology. Alongside remote mapping, field campaigns with collaborators have produced a large amount of tectonic and earthquake related information. The COMET Central Asia Fault Database assembles these data to produce a cohesive fault database of use to a wide range of geoscientists, as has been done recently across other regions of the planet.

Figure 2. Trenching with collaborators along the Dzhungarian Fault, Tien Shan, Kazakhstan (Credit: Austin Elliot)

Individual contributors had mapped faults at variable times (~2000 to present), variable resolutions (field-mapping to coarse satellite imagery), and for variable purposes (field site to continental-tectonic scale studies). Rather than reduce this variability to a single representation of the fault network, we produce a database that contains three resolution levels, to increase suitability for various applications.

The blog continues at Blog – Earthquakes in Central Asia.

In the meantime, if you’d like to be involved or would like more information, please get in touch with Tamarah King via tamarah.king@earth.ox.ac.uk

 

Harmony: Mission Candidate for the Earth Explorer 10

COMET scientists Professor Juliet Biggs (University of Bristol) and Professor Andy Hooper (University of Leeds) both serve on the Harmony Mission Advisory Group and are delighted to have been chosen to develop the concept further.

On February 18-19, ESA’s Programme Board for Earth Observation (PB-EO) decided on the continuation of the three Earth Explorer (EE) mission candidates towards the next phase in the path to their implementation. The three missions, namely, Daedalus, Hydroterra and Harmony, were selected in 2018 for a Phase-0 feasibility study out of 21 submitted proposals. The PB-EO has made now the unprecedented decision of selecting only one mission for Phase A, namely Harmony, instead of more than one as done in previous EE calls.

The Harmony mission is dedicated to the observation and quantification of small-scale motion and deformation fields at the air-sea interface (winds, waves, surface currents), of solid Earth (tectonic strain and height changes at volcanoes), and in the cryosphere (glacier flows and height changes). In order to achieve the different mission goals, the Harmony mission shall deploy two companion satellites following one of ESA’s Copernicus Sentinel-1 satellites. The companions will be flying in two different formations (see Figure 1): the stereo formation, with one Harmony satellite placed in front and one behind Sentinel-1, in both cases at a distance of about 350 km from it; and the cross-track formation, with both Harmony units flying close to each other (~200-500 m) also at 350 km from Sentinel-1. Each Harmony satellite carries as main payload a receive-only synthetic aperture radar (SAR), which shall acquire the reflected signals transmitted by Sentinel-1 towards the Earth. A multi-view thermal infra-red payload is also included to measure cloud height and cloud motion vectors. The angular diversity provided by the Harmonies in combination with Sentinel-1 will allow the retrieval of deformation measurements of the sea and earth surface with unprecedented accuracy (see Figure 2), while the cross-track configuration will allow the accurate measurement of elevation changes for land-ice and volcanic applications.

 

 

 

 

 

Figure 1: Representation of the (left) stereo and (right) cross-track flying formations for Harmony. The Sentinel-1 satellite is depicted in black color. Sentinel-1 transmits a signal and acquires the backscattered echoes (represented with magenta arrows), while the Harmony satellites receive part of the energy that bounces towards them (represented with the green arrows). Copyright: Harmony Mission Advisory Group.

Dr. Paco López-Dekker from the Delft University of Technology and principal investigator of the Harmony mission, comments “It is very exciting that our multi static-SAR concept, which combines many ideas that were matured during my years at HR, has made it to this final stage. During Phase-0 we have drafted a beautiful and elegant mission concept promising an unprecedented view at Earth System processes. Now we have the responsibility to look at it from all sides and be sure that it will work. Challenging and fun.”

Professor Juliet Biggs from the University of Bristol and member of Harmony’s Mission Advisory Group at ESA adds “The Harmony mission is remarkable in that it promises new scientific discoveries across an astonishing breadth of topics: from the gradual motion of tectonic plates to small-scale processes on the ocean surface. I’m delighted that we have been selected to develop the concept further and that Harmony is one step closer to becoming a reality”

Dr. Pau Prats, from the German Aerospace Centre, DLR and member of Harmony’s Mission Advisory Group at ESA, is convinced of the benefits a mission like Harmony will bring to the community: “The unique configuration of the Harmony satellites in combination with Sentinel-1 will allow us to literally add a new dimension to SAR observations, a fact that will foster SAR technology and its applications during the next two decades.”

Figure 2: Coloured areas show regions straining at greater than 10 nanostrain per year (the threshold above which 95% of earthquake fatalities occur). Blue regions are those that have a small component of north-south strain and can be imaged by Sentinel-1 alone. Red regions indicate the extra area that will be constrained by Harmony. From Harmony Report for Assessment. 2020.  

So, what’s next? Even though Harmony is currently the only EE-10 mission candidate it does not mean it will be implemented. The industry and science teams have one and a half years of hard work ahead to demonstrate the mission has reached the technological and scientific level of maturity required to enter into the next phase, that will ultimately result in the launch of the Harmony satellites by the end of this decade.

Announcement can be found on ESA website: https://www.esa.int/Applications/Observing_the_Earth/ESA_moves_forward_with_Harmony

Title figure for the Harmony mission. Copyright: ESA.

 

Earthquake research could improve seismic forecasts

The timing and size of three deadly earthquakes that struck Italy in 2016 may have been pre-determined, according to new research published in Earth and Planetary Science Letters that could improve future earthquake forecasts.

Credit: L.Gregory

A joint British-Italian team of geologists and seismologists have shown that the clustering of the three quakes might have been caused by the arrangement of a cross-cutting network of underground faults.

The findings show that although all three earthquakes occurred on the same major fault, several smaller faults prevented a single massive earthquake from occurring instead and also acted as pathways for naturally occurring fluids that triggered later earthquakes.

The cluster of three earthquakes, termed a “seismic sequence” by seismologists, each had magnitudes greater than six and killed more than 300 people in Italy’s Apennine mountains between 24 August and 30 October 2016.

Earthquake sequences

The research, led by COMET scientist Richard Walters from Durham University, comes on the second anniversary of the start of the earthquake sequence.

The researchers say the findings could have wider implications for the study of seismic hazards, enabling scientists to better understand potential earthquake sequences following a quake.

Dr Walters said: “These results address a long-standing mystery in earthquake science – why a major fault system sometimes fails in a single large earthquake that ruptures its entire length, versus failing in multiple smaller earthquakes drawn-out over months or years.

“Our results imply that even though we couldn’t have predicted when the earthquake sequence would start, once it got going, both the size and timing of the major earthquakes may have been pre-determined by the arrangement of faults at depth.

“This is all information we could hypothetically know before the event, and therefore, this could be a hugely important avenue for improving future earthquake forecasts.”

Location of survey site at rupture across a road near Castelluccio. The rupture occurred during the third earthquake in the seismic sequence and gives researchers a record of the deformation. Credit: L.Gregory

Thousands of aftershocks

Dr Walters and the team used satellite data to estimate which part of the fault failed in each earthquake, and compared this pattern with the location and timing of thousands of tiny aftershocks throughout the seismic sequence.

They found that intersections of small faults with the main fault system separated each of the three largest earthquakes, strongly suggesting these intersections stop the growth of each earthquake and prevent the faults failing in a single large event.

In addition, the scientists also found that after the first earthquake, thousands of aftershocks crept northwards along these same fault intersections at a rate of around 100 metres per day, in a manner consistent with naturally occurring water and gas being pumped along the faults by the first earthquake on 24 August, 2016.

The second earthquake, on the 26 October, occurred exactly when these fluids reached its location, therefore controlling the relative timing of failure.

Dr Walters added: “It was a big surprise that these relatively small faults were having such a huge influence over the whole sequence.

“They stop the first earthquake in its tracks, and then they channel the fluids that start the sequence up again months later. No-one’s ever seen this before.”

Seismic hazard

Co-author and COMET associate Dr Laura Gregory from the University of Leeds said it was important to understand whether or not a fault fails in a seismic sequence, and that the team’s results were only made possible by combining a varied array of different datasets.

Dr Gregory said: “A seismic sequence has vastly different implications for seismic hazard compared to a single large earthquake. If the faults in Italy in 2016 had failed together in one big event, the impact on the local population would have been much worse.

“This is the first time we’ve ever had this quality of modern data over one of these earthquake sequences, and bringing together a range of specialists was key for unpicking how the earthquakes related to one another.

“I was scrambling over the mountainside immediately after each earthquake with British and Italian colleagues, measuring the metre-high cliffs that had suddenly formed. Meanwhile, other members of our team were analysing data from seismometers stationed around the world, or were mapping the tiny bending of the ground around the faults using satellites orbiting the planet at 500 miles altitude.”

The research was partly supported by both COMET and a NERC Urgency Grant.

 

 

Turkish fault reveals seismic steadiness

Satellite data has shed new light on seismic hazard in one of the world’s most deadly earthquake zones.

Published today in Nature Communications, the COMET study describes how tectonic strain builds up along Turkey’s North Anatolian Fault at a remarkably steady rate.

This means that present-day measurements can not only reflect past and future strain accumulation, but also provide vital information on events still to come.

The strain, which builds up as Turkey is squeezed between three major tectonic plates, has caused almost the entire length of the fault to rupture since 1939 in a series of major earthquakes gradually migrating east-west towards Istanbul.

Strain rates along Turkey’s North Anatolian Fault (past ruptures shown in purple/yellow) alongside westward progression of earthquakes since 1939

Led by COMET PhD student Ekbal Hussain[1], the team used satellite images from the European Space Agency’s Envisat mission to identify tiny ground movements at earthquake locations along the fault.

Dr Hussain explained: “Because we know so much about the fault’s recent history, we could look at the strain build up at specific places knowing how much time had passed since the last earthquake.”

The 600-plus satellite images, taken between 2002 and 2010, provided insights into the equivalent of 250 years of the fault’s earthquake repeat cycle.

Remarkably, apart from the ten years immediately after an earthquake, strain rates levelled out at about 0.5 microstrain per year, equivalent to 50mm over a 100km region, regardless of where or when the last earthquake took place.

Dr Hussain added: “This means that the strain rates we measure over the short term can also reflect what’s happening in the longer term, telling us how much energy is being stored on the fault and could eventually be released in an earthquake.”

Istanbul, straddling the Bosphorus Strait at the centre of the picture, and the surrounding area in northwestern Turkey captured by Envisat.  Turkey’s location makes it vulnerable to earthquakes, with the North Anatolian fault lying just 15 km south of Istanbul. Credit: ESA

Until the satellite era, it was difficult to get a clear picture of how strain built up on the fault.  Now, satellites like Envisat, alongside the newer Sentinel-1 mission, can detect ground movements of less than a millimetre, indicating how and where strain is accumulating.

The findings suggest that some existing hazard assessment models, which presume that strain rates vary over time, need to be rethought.  This is especially true for regions where there are long gaps between earthquakes, such as the Himalayas.

Co-author and COMET Director Tim Wright said: “Discovering this consistent strain accumulation will help us to reassess how we model seismic hazards, as well as improving understanding of the earthquake cycle worldwide.”

The full paper is: Hussain et al. (2018) Constant strain accumulation rate between major earthquakes on the North Anatolian Fault, Nature Communications, doi:10.1038/s41467-018-03739-2

[1] Now Remote Sensing Scientist at BGS Keyworth.  Dr Hussain is available for comment (ekhuss@bgs.ac.uk).

Sentinel-1 satellites reveal pre-event movements and source areas of the Maoxian landslides, China

At about 5:38am local time on 24 June 2017 (21:38 on 23 June 2017 UTC), a massive landslide struck Xinmo Village, Maoxian County, Sichuan Province in China. Sichuan province is prone to earthquakes, including the 2008 Mw 7.9 Wenchuan earthquake that killed over 70,000 people, as well as the 1933 Mw 7.3 Diexi earthquake with a death toll of up to 9300.

Authorities have confirmed that the Maoxian landslide was triggered by heavy rain. The Maoxian landslide swept 64 homes in Xinmo village, blocking a 2km section of river and burying 1,600 meters of road. The collapsed rubble was estimated to be about eight million cubic meters (Figure 1a). Three days later (on 27 June 2017), a second landslide hit Xinmo Village; almost in the same time, another landslide occurred in Shidaguan Town, 20km away from Xinmo Village (Figure 1b).

Figure 1. The Maoxian landslides: (a) Xinmo Village (24 June 2017 local time); (b) Shidaguan  Town (27 June 2017 local time)

As part of the ESA’s Copernicus Program, the Sentinel-1 mission comprises a constellation of two polar-orbiting satellites, operating day and night performing C-band synthetic aperture radar imaging in all-weather conditions. Sentinel-1 images acquired before and after events such as landslides, earthquakes or volcanoes offer information on the extent and surface displacements of affected areas, which can be used for damage and future hazard assessment.

A joint team from Newcastle University (UK), Chengdu University of Technology, Tongji University, China Academy of Space Technology and Wuhan University (China) have been racing against time to respond these two events by combining ESA’s Sentinel-1, Chinese Gaofen-2/3 with field observations.

“Sentinel-1 acquired a post-event image thirteen and half hours after the Xinmo event, and provided us the first interferogram for the Xinmo landslide.” said Professor Zhenhong Li, Professor of Imaging Geodesy at Newcastle University, “This first Sentinel-1 interferogram, together with its corresponding coherence and amplitude maps, not only helped us identify the source area of this massive landslide, but also assisted with mapping the landslide boundary (Figure 2).”

Figure 2. The Xinmo landslide: (a) the first Sentinel-1 interferogram (20170612-20170624);  (b) UAV imagery. Note: (1) Red lines represent the landslide boundary derived from InSAR observations; (2) Yellow lines indicate relatively stable areas; and (3) Green and while circles imply the source area of this massive landslide.

“More importantly, through the analysis of the archived Sentinel-1 data, we found that pre-event movements exhibited in the source area during the period from 14 May to 19 June 2017 for the Xinmo event.”

GIF

Pre-event signals are even clearer for the Shidaguan landslide, suggesting it had been sliding for a while.

GIF

“It is well known that landslides are hard to predict” said Professor Qiang Xu, Geologist with Chengdu University of Technology, “This study convincingly demonstrates that InSAR can be used to detect and map active landslides, which is a great achievement.”

After presenting the joint research findings in the Dragon-4 symposium at Copenhagen on 27 June 2017, Professor Deren Li, Academician of Chinese Academy of Sciences (CAS) and Chinese Academy of Engineering (CAE) at Wuhan University concluded, “This joint effort suggests that landslide Early Warning System (EWS) might be encouraging. There are a range of factors we should consider for EWS, such as real time, automatic, spatial and temporal resolutions, and data uncertainty.”

 

Massive quake shifts perception of earthquake hazard models

Dramatic data from last year’s major earthquake in Kaikoura, New Zealand, will change the way scientists think about earthquake hazards in tectonic plate boundary zones.

The South Island earthquake was one of the most comprehensively recorded earthquakes in history.  Satellite images of the earthquake and its aftermath have enabled scientists to analyse the quake in an unprecedented level of detail.

COMET Director Tim Wright, study co-author, said: “We’ve never seen anything like the Kaikoura quake before, it was one of the most complex ever recorded.

“An earthquake commonly ruptures across a single fault line or faults that are closely grouped; Kaikoura ruptured at least 12 major faults. This challenges many assumptions about how individual faults control earthquake ruptures.”

Immediately following the start of Kaikoura’s earthquake in November, the Sentinel-1 and ALOS-2 satellites were tasked with gathering images and data.  COMET researchers at the University of Leeds used the data to provide a rapid analysis of what was happening during the quake.

They found that seismic readings of the earthquake were not giving accurate assessments of where the ruptures were occurring.

Professor Wright said: “Seismic readings are currently the fastest method of gathering earthquake data as they can quickly read shockwaves sent through the earth – but they paint a crude picture. The complexity of the Kaikoura earthquake caused seismologists to completely misinterpret the earthquake based on these seismic waves.

“Satellites such as the European Space Agency’s Sentinel-1 measure the way the ground deforms at very high resolution. The detail in the images showed us that ruptures took place across many separate faults.

“One of the aims of COMET is to expand satellite capability to provide rapid response earthquake data to ensure residents and rescue workers have access to fast and accurate information.”

The satellites provided pre- and post-earthquake images to measure the extent of land movement. Kaikaoura’s earthquake caused sections of earth to move up to 25 metres and created surface ruptures measuring 12 metres. This caused large scale landslides and triggered a tsunami.

Map of the South Island of New Zealand showing the ground deformation and location of the surface rupture from the November 2016 Kaikoura earthquake. The rainbow contours of ground motion are in 20 cm intervals whilst the red/blue colours show up to +/-6 m of motion near the fault. The sharp change in colour nearer the coast is where the fault has ruptured through to the surface.

In a study published today in Science, the team of researchers reports the full range of data analysed from Kaikaoura’s earthquake, including satellite imagery, field observation, GPS data and coastal uplift data.

The research will prompt reassessment of how many different faults can be involved in a single earthquake and could potentially feed into revaluations of seismic hazard models.

Lead researcher Dr Ian Hamling, a natural hazards geodesist from New Zealand research institute GNS Science, said: “There was growing evidence internationally that conventional seismic hazard models are too simple and restrictive.

“The message from Kaikoura is that earthquake science should be more open to a wider range of possibilities when rupture models are being developed. It underlines the importance of re-evaluating how rupture scenarios are defined for seismic hazard models.”

New Zealand’s complex network of faults are similar to those found in western United States, Japan and central Asia.

Leeds co-author and COMET Associate Dr John Elliott said: “While earthquakes like Kaikoura’s do not commonly occur, the data we’ve gathered from this event will expand our understanding of similar boundary zones around the world.

“Not only could the data help inform us for the future but it may change how we’ve interpreted ancient earthquakes. If an earthquake like Kaikoura’s took place thousands of years ago, current methods of paleo-seismology would possibly see it as a series of earthquakes over a long period of time, rather than as one large single quake.”

The full paper, ‘Complex multi-fault rupture during the 2016 M7.8 Kaikoura earthquake, New Zealand’ was published on 23 March in Science First Release.

Earthquake resilience in the East African Rift: the challenge of making baseline measurements

By COMET Scientist Juliet Biggs

Malawi is a small country at the tip of the East African Rift. Small, that is, until you spend several weeks driving around in loops over challenging roads.

Ake Fagereng (Cardiff), Hassan Mdala (Geological Survey) and I covered nearly 5000 km in less than three weeks. Why would we do this to ourselves? To establish a baseline geodetic network at the southern tip of the East African Rift.

We want to know how strain is distributed, and which are the active structures. This latest field season is just another piece in the tectonic jigsaw that we hope will ultimately lead to improving earthquake resilience in developing countries.

Our partner for the last few years has been the Geological Survey Department of Malawi – in particular, Hassan, who operates a small mobile office in the back seat of the car. From here, he fixes all our problems, often before we even realize they exist; making use of a boundless network of friends and family, all of whom seem to be just waiting to help with our latest peculiar request.

Each GPS station is composed of an antenna, a receiver, a car battery, a regulator, a solar panel and a metal pin that we drill into the rock. Most of the kit only stays for 4 days, but we’ll need to re-measure the position of the pin for many years to come. In our search for suitable sites, we met many enthusiastic head teachers, police chiefs and district commissioners and one rather terrifying headmistress, all of whom helped with security arrangements in schools or telecommunications masts.

At Mua Primary school, with the Head Teacher, Deputy Head and our GPS antenna
At Mua Primary school, with the Head Teacher, Deputy Head and our GPS antenna

Everywhere we went, there were spectacular landscapes and interesting geology, most of which we were forced to drive straight past. But so much of my research is done from satellites orbiting at 600 km that I can get a bit over-excited about getting my nose to the rock again, and we took advantage of what opportunity we had.

Malawi is largely composed of various types of gneiss, a high-grade metamorphic rock, with lots of migmatites, blobs that were once partially molten rock. Perhaps the most spectacular landscapes are the ‘inselbergs’, large granite intrusions, which now stick several thousand metres above the surrounding landscape. On our only day off, squeezed somehow into a hectic schedule, we hiked up the largest of these, Mulanje, for some truly spectacular views. We also found some fascinating faults, but you’ll have to wait for the paper to hear about those!

Hiking on Mulanje Massif. According to the Official Malawi Tourism site, the third largest mountain in Africa – a rather dubious statistic, but still steep enough to trouble my knees
Hiking on Mulanje Massif. According to the Official Malawi Tourism site, the third largest mountain in Africa – a rather dubious statistic, but still steep enough to trouble my knees

Many people ask about the personal challenges to fieldwork in Africa, and there are indeed plenty of those: insect bites, the side-effects of malaria medication, power cuts, cold showers, endless meetings getting permission, long days driving over roads covered with pot-holes, goats, dogs, chickens and bicycles, and on this occasion a car break-in (I’m wondering if there’s such a thing as a frequent flyer card for local police stations…?).

But none of that matters much; I love every minute (well, almost, I could skip the police stations). A great trip requires great company: Hassan, for whom every minor disaster is ‘a learning opportunity’, and Ake whose response to the next challenge is usually ‘well, that’ll be fun’.

Ake and Hassan contemplating a job well done. ‘Til next time.
Ake and Hassan contemplating a job well done. ‘Til next time

All I can say, is I learned a lot and had fun doing it, so I guess they are both right.  I simply feel privileged to spend time in a country with such interesting scientific questions, beautiful landscapes and warm, friendly people. I suspect we’ll find a way to come back before too long!

Our fieldwork was funded by Global Challenges Research Fund through an EPSRC Institutional Sponsorship Award to the University of Bristol; a BGS-COMET fieldwork award and a Researcher Mobility Grant from the World Universities Network.

Detecting geohazards with GPS

COMET’s Marek Ziebart and Chris Atkins, both at UCL, have been using sidereal filtering to detect geohazards in near-real time.  This article describes how they have applied this to earthquakes in the Aleutian Islands, Alaska.

Geodetic quality GPS (Global Positioning System) data – fundamentally high precision range measurements made between orbiting spacecraft and tracking stations, often rigidly attached to bedrock – is one of the key observables available to us to analyse plate tectonics and the earthquake cycle. Such tracking stations now often make measurements at 1 Hz or higher. This enables us to calculate high rate position time series for the tracking stations, revealing various kinds of motion related to seismic waves and other forms of deformation.

 Figure 1: Deep braced GNSS tracking antenna at the Aleutian Islands, Cape Sarichef, Alaska (AV24-WestdahlNWAK2008). Photo: UNAVCO.org

Figure 1: Deep braced GNSS tracking antenna at the Aleutian Islands, Cape Sarichef, Alaska (AV24-WestdahlNWAK2008). Credit: UNAVCO.org

One such time series is shown below, calculated by a standard processing method called PPP (precise point positioning).  The plot shows the change in position in the north-south direction. When do you think the earthquake seismic waves begin to arrive?

Marek2

Figure 2: 25 minutes of the 1Hz GPS time series for AV24 – a deep braced GNSS antenna on the Aleutian Islands, Cape Sarichef, Unimak Island, Alaska at the time of the 2011 Tohoku earthquake (WestdahlNWAK2008)

GPS Seismology

How can we use this type of data? The time series in figure 2 shows an estimate of the GPS antenna movement in the North-South direction during the 2011 magnitude 9.0 Tohoku earthquake over a period of 25 minutes. This is a measurement of the ground shaking – the data is used to understand the nature of the earthquake and to compare it to other earthquakes in different places around the world. Earthquakes often cause a change in the shape of the ground that is retained after the earthquake is over, and which can be determined from a time series. This kind of information can be used to estimate the size of the earthquake, how much energy was released and tells us about the earthquake’s epicentre location.

Another kind of instrument, a seismometer, can be used to determine similar information. The drawback of seismometers is that they can ‘clip’ when the seismic waves become too large to handle. A GPS antenna, on the other hand, has no such limitation – no matter how big the amplitude of the seismic waves, they can still be measured using GPS technology. Because of this, GPS antennas such as the ones shown in figures 1 and 6 below are sometimes called broadband GPS seismometers. In practice we combine information from both GPS antennas and seismometers – they are complimentary. There are now thousands of such GPS antennas installed around the world, connected to the internet and taking measurements 24 hours a day.

Marek3

Figure 3: an example of a GPS tracking station network – this one is in the Gulf of Alaska and forms part of the US Plate Boundary Observatory (PBO)

So what’s the problem? The problem is Multipath – this is akin to radio interference where the signals transmitted by the orbiting satellites have more than one way of getting to the antenna. The strongest, direct signal – which is the one we want to measure – is distorted by signals bouncing off adjacent reflecting surfaces in the antenna’s vicinity. This causes spurious patterns of motion in the receiver time series. These can vary rapidly over a few seconds in a noise-like way (with amplitudes at the few mm level), or they can appear as a random walk – a slowly varying position error (with changes at the level of a few centimetres). Both these effects limit our ability to use GPS data for seismometry. To get some feel for the effect consider Figure 4.

On the day of the earthquake (shown as the black time series) the onset of the arrival of the largest seismic waves is shown clearly at around 6:03:00. A not dissimilar signature is also visible around 5:50:00. However, that same signature (an apparent ground movement of some 25mm over a few minutes, shaking and then reducing in magnitude) is visible at almost the same time on the previous day (in the grey time series). It doesn’t mean the earthquake was repeating – these are multipath distortions. They appear indistinguishable from seismic waves. Our problem, then, is how to remove the spurious, apparent motion from the time series, whilst retaining the actual signature of the earthquake. So what’s the solution?

Marek4

Figure 4: the GPS time series (North-South component) for station AV24 on the day of the earthquake and on the day before. The time series on the vertical axis are offset by 60mm for clarity

Sidereal Filtering

Sidereal filtering is a technique used to reduce errors caused by multipath in the positioning of GPS receivers. It relies upon the receiver remaining static from one day to the next relative to its surrounding environment and takes advantage of the ground track repeat time of the GPS satellites, which is just less than a sidereal day (usually around 23 hours, 55 minutes, 55 seconds). The repeating multipath error can thus be identified and largely removed from the following day by subtraction.

A conventional position-domain sidereal filter (PDSF) identifies and removes this repeating pattern from a position time series. However, the ground track repeat time of individual GPS satellites can differ from each other by a few seconds, whereas a PDSF has to assume that all satellites have the same repeat time. This can cause problems for the PDSF, especially when the oscillations caused by multipath interference are particularly high. We have developed an observation-domain sidereal filter (ODSF) that identifies and removes multipath errors from the GPS phase measurements themselves.

Unlike the PDSF, it can account for the fact that the ground track repeat times of the GPS satellites differ by a few seconds from one another. This means that it is more effective at removing the effects of high-frequency multipath error and is less sensitive to satellite outages. For each phase measurement, the ODSF algorithm searches for an appropriate correction based on the azimuth and elevation of the relevant satellite. That correction is derived from the measurement residuals on the previous day that most closely correspond to that azimuth and elevation. The precision of these measures of azimuth and elevation needs to be high – a hundredth of a degree or better. This is readily achievable in our algorithm.

Marek5

Figure 5: The same GPS time series for station AV24, but with the time series resulting from the two types of sidereal filter also shown. The time series on the vertical axis are offset by 50mm for clarity

Figure 5 shows the same 25-minute time series of 1 Hz displacements in northing as shown in Figure 4. Also plotted in red and green are the time series that result after applying the two types of sidereal filter: the PDSF and ODSF respectively. Notice that both the PDSF and the ODSF are largely successful at removing these oscillating errors. They reveal a surface wave arrival time of about 06:00 with the largest Love waves arriving at around 06:03. These events would not have been so easy to distinguish in the original unfiltered time series.

Marek6

Figure 6: Deep braced GNSS tracking antenna on the Aleutian Islands, Cape Sarichef, Unimak Island, Alaska (AV27-WestdahlSWAK2008). Credit: UNAVCO.org

Marek7

Figure 7: GPS time series for station AV27 including those resulting from the two types of sidereal filter. The time series on the vertical axis are offset by 50mm for clarity

Figure 7 similarly shows 1 Hz displacements at station AV27 which is only 11 km from AV24. One would therefore expect a very similar displacement signal. However, it is clear that the standard PPP time series, plotted in black and grey, are severely affected by strong short-period (~11 s) multipath error. In this case, the PDSF was unable to remove such a high-frequency error and instead increased the amplitude of these oscillations. The ODSF on the other hand was far more effective because it could take into account the differing repeat times of the satellites.

Conclusions

The precision of GPS technology has improved dramatically since 1995 when the system reached FOC (full operational capability). It is an effective tool for measuring geohazards and other geophysical phenomena in near real-time. However, data processing and interpretation must be handled carefully. Multipath effects on position time series can introduce spurious signatures that resemble seismic events and hamper our ability to exploit the technology. Sidereal filtering offers a very effective way to clean up GPS positioning time series to reveal more clearly geophysical events, effectively in real-time.