Coulomb stress
interaction (click
here for a
recent review
paper)
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.
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.
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)
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
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.