1 Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand Polina Berezina 1 Institute of Geology, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine Supervisor: Prof. Kenji Satake Earthquake Research Institute, The University of Tokyo, Tokyo, Japan (Dated: August 20, 2017) Abstract Computer simulation of tsunami generated by the 2009 Dusky Sound earthquake was done to investigate tsunami characteristics and energy. A single-fault model was created and further improved to generate a seafloor displacement with Okada’s solution and a matching sea surface displacement. The simulation was run up to seven hours with a bathymetry grid resolution of 60 arc-seconds and time interval of 5 seconds. It was found that the generated tsunami had a maximum height of around 0.6 m, which was characteristic of the near-fault area. In addition, computed synthetic waveforms closely matched observed DART and tide gauge records at the time of the tsunami arrival. As time passed, the computation results deviated from the gauge data due to an increased error, as well as noisiness of observed waveforms at bays. 1. INTRODUCTION The Dusky Sound earthquake occurred on July 15 2009 at 09:22 am GMT. The epicenter was located in the Fiordland region, South Island of New Zealand. Having a magnitude of 7.8 (Mw), it was one of the biggest earthquakes in the history of New Zealand, on par with the 2016 Kaikoura earthquake [1]. Multiple aftershocks were recorded in the following 20 days with a magnitude smaller than 5.7, as illustrated in figure 2. The geological setting of the quake is characterized by oblique subduction of the Pacific oceanic plate under Indo-Australian continental plate. The earthquake struck in the northern part of the Puysegur trench, which represents one of the most seismically active areas in New Zealand. The convergence rates in this region is about 35- 45 mm/year [2]. Figure 1 shows earthquakes with a magnitude larger than 4.5 from July 2009 to June 2017, according to existing CMT solutions. The “beach ball” representations of focal mechanisms, as shown in figure 2, describe a strike-slip faulting with a general strike orientation along the subduction zone. On the right of the image, a total frequency of these earthquakes is presented in a logarithmic scale. Evidently, earthquakes with a magnitude smaller than 6.0 are predominant in the region. The earthquake was felt in the South Island and southern part of the North Island, causing a minor damage. It triggered several landslides in the Fiordland National Park near the Dusky Sound. As a precaution, the Pacific Tsunami Warning Center Potential issued tsunami warnings that were soon diminished [1]. There were only some tsunami deposits preserved, vegetation disturbances, and little coastal deformation due to a gradual rather than sudden motion during the earthquake [3]. Figure 1. Epicenter (a star) and aftershocks, according to USGS data. The diameter and color intensity correspond to the magnitude
6
Embed
Tsunami Simulation of 2009 Dusky Sound Earthquake in … · Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand ... The grid resolution of 60 arc ... obtained through
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
1
Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand
Polina Berezina
1 Institute of Geology, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
Supervisor: Prof. Kenji Satake
Earthquake Research Institute, The University of Tokyo, Tokyo, Japan
(Dated: August 20, 2017)
Abstract
Computer simulation of tsunami generated by the 2009 Dusky Sound earthquake was done to
investigate tsunami characteristics and energy. A single-fault model was created and further improved to
generate a seafloor displacement with Okada’s solution and a matching sea surface displacement. The
simulation was run up to seven hours with a bathymetry grid resolution of 60 arc-seconds and time interval
of 5 seconds. It was found that the generated tsunami had a maximum height of around 0.6 m, which was
characteristic of the near-fault area. In addition, computed synthetic waveforms closely matched observed
DART and tide gauge records at the time of the tsunami arrival. As time passed, the computation results
deviated from the gauge data due to an increased error, as well as noisiness of observed waveforms at bays.
1. INTRODUCTION
The Dusky Sound earthquake occurred on July
15 2009 at 09:22 am GMT. The epicenter was
located in the Fiordland region, South Island of
New Zealand. Having a magnitude of 7.8 (Mw), it
was one of the biggest earthquakes in the history
of New Zealand, on par with the 2016 Kaikoura
earthquake [1]. Multiple aftershocks were
recorded in the following 20 days with a
magnitude smaller than 5.7, as illustrated in
figure 2.
The geological setting of the quake is
characterized by oblique subduction of the Pacific
oceanic plate under Indo-Australian continental
plate. The earthquake struck in the northern part
of the Puysegur trench, which represents one of
the most seismically active areas in New Zealand.
The convergence rates in this region is about 35-
45 mm/year [2]. Figure 1 shows earthquakes
with a magnitude larger than 4.5 from July 2009
to June 2017, according to existing CMT solutions.
The “beach ball” representations of focal
mechanisms, as shown in figure 2, describe a
strike-slip faulting with a general strike
orientation along the subduction zone. On the
right of the image, a total frequency of these
earthquakes is presented in a logarithmic scale.
Evidently, earthquakes with a magnitude smaller
than 6.0 are predominant in the region.
The earthquake was felt in the South Island
and southern part of the North Island, causing a
minor damage. It triggered several landslides in
the Fiordland National Park near the Dusky
Sound. As a precaution, the Pacific Tsunami
Warning Center Potential issued tsunami
warnings that were soon diminished [1]. There
were only some tsunami deposits preserved,
vegetation disturbances, and little coastal
deformation due to a gradual rather than sudden
motion during the earthquake [3].
Figure 1. Epicenter (a star) and aftershocks,
according to USGS data. The diameter and color
intensity correspond to the magnitude
2
Figure 2. Earthquakes with magnitude > 4.5 since 2009 until 2017 and their CMT solutions
2. METHODS
In case of a shallow-water wave or a long wave,
water surface displacement can be approximated
to the seafloor displacement [4]. That is why
seafloor displacement is calculated first to further