Earthquake physics

Coulomb stress interaction (click here for a recent review paper)

Slip_840-1003.jpg (731518 bytes)

SeabedMov_840_1003.jpg (507548 bytes)

Tsunami hazard in the Indian Ocean

In light of the research described in the section below, we have formed a collaboration with Kerry Sieh of Caltech and Massimo Cocco, Alessio Piatanesi, Carlo Giunchi, and Spina Concetti of INGV - Rome to model the next likely tsunami in the Indian Ocean. This NERC funded project involves the following steps:

1) We simulate a range of possible slip distributions for different potential earthquakes. As the length of any future large earthquake is unconstrained, we consider 4 possible lengths based on the geology - 210, 330, 630, and 840 km. For each rupture length, we generate a suite of synthetic slip distributions with realistic moments and heterogeneity and without unphysical effects such as large slips at the edges. 2) These slip distributions are sent to Kerry who ranks them based on his geologic knowledge of past ruptures. This filter leads to a range of potential earthquakes that are somewhat characteristic but we also retain a small number of less likely slip distributions. 3) Andrea Antonioli has developed a 3D finite element model of the Sunda Trench. He uses the chosen slip distributions as inputs to his model and calculates the sea-floor displacements resulting from each. 4) These sea-floor displacements are sent to Alessio in Rome who computes tsunami wave heights and wave energies along the coast of Sumatra and in the Indian ocean. We have digitized the bathymetry just offshore Sumatra as well as its coastline in order to improve the
resolution in this vitally important area. 5) Finally, we have a range of possible tsunamis that may affect the region. We are exhibiting this work at the Royal Society Summer Science Exhibition in London and Glasgow. The exhibit is supported by EPSRC and the Royal Society.
eqs1&2_stress_sumatra.jpg (33968 bytes)

coseismic_sumatra2.jpg (614679 bytes)

Seismic hazard in Sumatra

Following the great 26 Dec.2004 Sumatra earthquake, we calculated the stress changes from it on the Sunda Trench immediately to the south and on the Sumatra Fault that runs down the middle of the island. We found that stresses had increased significantly on both structures and in a paper published in Nature on 17 March 2005 we warned that the seismic hazard had consequently increased in the area. On 28 March 2005 an M=8.7 earthquake occurred in the area of concern along the Sunda Trench although fortunately it did not trigger a large tsunami.

We have since recalculated the stresses from the two earthquakes and in a subsequent paper highlight the continued hazard along the Sumatra Fault as well as the increased hazard of a further earthquake and tsunami in the Mentawai Islands region, a portion of which has not experienced a large earthquake since 1797. We are now seeking funding from NERC to investigate stress interaction and the timing of future large earthquakes in a global catalog.

sc4_sections2.jpg (889341 bytes)

Aftershock hazard

The PRESAP project (EVG1-1999-00001) was led by the Geophysics Group and funded under Framework 5. It was designed to test the feasibility of making near real time Coulomb stress calculations in order to estimate the likely spatial distribution of aftershocks following moderate to large earthquakes in Europe. Although previous retrospective research had shown a clear causal link between stress changes and aftershocks, no one had investigated whether it would be possible to assess aftershock hazard in the hours to days following an earthquake. The aims of PRESAP were to determine whether this was both scientifically and technologically feasible.

The results and methodologies developed during the project were tested in a pseudo real time scenario in which data (disguised in location and time) were fed to participants in the time sequence in which it would have been available in reality. Only 4 hours and 25 minutes after the main shock, a Coulomb stress map was “released”; this release highlighted the danger of a particular fault that ruptured in a large magnitude event approximately 6.5 hours later. A subsequent release, 14 minutes before the large aftershock, increased the emphasis on this potential region of high hazard.

The results from PRESAP clearly show that it is both scientifically and technologically feasible to make real time estimations of the likely spatial distribution of aftershocks. The challenge now is to move the science from “feasible” to “practical” in order to reduce aftershock hazard in Europe.

Click for some of our papers from this project ( pap1, pap2, pap3)

home.h2.jpg (17612 bytes) Stress interaction in Turkey

In 2002, we published a paper on Coulomb triggering and earthquake hazard along the East Anatolian Fault in Turkey. We identified two regions as being highly loaded and suggested that detailed paleoseismological work be carried out to investigate previous earthquakes along specific segments. In 2003, an M=6.4 earthquake occurred near the city of Bingol. This event was not, however, on the EAF but rather on a small conjugate fault. As our subsequent paper shows, the relationship between this earthquake and activity on the EAF is complex and has implications for commonly used methods of secular loading. This work was funded by NERC and a follow-up project currently underway is funded by an EC Intra-European Fellowship. In a related project, Emre Evren is presently doing Ph.D. research on earthquake hazard in the Istanbul area.

Earthquake modeling

home.h1.jpg (56555 bytes) We have investigated the effects of fault zone heterogeneity on the nucleation and growth of large earthquakes ( link) and the magnitude-frequency distribution in 2D cellular automata. In the latter paper we developed a new set of stress transfer rules that allow for realistic stress concentrations in a nearest-neighbor model.

We recently completed an EPSRC funded project to investigate the conditions under which the magnitude-frequency distribution regionally and on individual faults obeys a power-law. The research was carried out in a complex 3-D network model that combines cellular automata with realistic nearest-neighbor stress transfer on individual faults and 3-D elastic stress transfer between faults. As part of the project, the model was parallelized and we acquired a 24 processor high performance computing system. The model is presently being modified as part of Alison Hetherington's Ph.D. project.