Journal of the Persian Gulf
(Marine Science)/Vol. 6/No. 22/December 2015/13/35-48
Numerical Modeling of Tsunami Waves Associated
With Worst Earthquake Scenarios of the Makran
Subduction Zone in the Jask Port, Iran
Akbarpour Jannat, Mahmood Reza*1; Rastgoftar, Ehsan1
Iranian National Institute for Oceanography and Atmospheric Science (INIOAS), Tehran, IR Iran
Received: June 2015
Accepted: November 2015
© 2015 Journal of the Persian Gulf. All rights reserved.
Abstract
The recent studies show that the past researches may have significantly underestimated
earthquake and tsunami hazard in the Makran Subduction Zone (MSZ) and this region is
potentially capable of producing major earthquakes. In this study, the worst case possible
earthquake scenarios of the MSZ are simulated using fully nonlinear boussinesq model to
investigate tsunami hazards on the Jask Port, one of the most important southern coastal areas
of Iran. The 9.1Mw and 8.7Mw earthquakes which respectively, represent the rupture of about
full and half length of the MSZ plate boundary, have been considered in the modeling. The
global and local numerical simulation were conducted based on 1-minute and 3-arc sec
resolution bathymetry data, respectively, in order to capture tsunami generation, propagation
and inundation at the Jask Port. Results show the water subsidence is firstly obsereved along
the Jask Port coastlines just after the earthquake occurrence, which can alert residents as a
natural warning sign. Model results also reveal that due to the 9.1Mw earthquake and the
8.7Mw western Makran earthquake, the most coastal areas around the Jask Port are inundated
by tsunami waves. However, the tsunami run-up is not observed when the 8.7Mw earthquake
occurs at the eastern segment of the MSZ because the tsunami waves almost propagate in the
north-south direction, perpendicular to the Makran fault and it causes a little effect of the
tsunami attack on the western coast of Makran. Therefore, the most crucial factor that
determines the tsunami risk at the Jask Port is the location of earthquake focal points in the
Makran region.
Keywords: Makran Subduction Zone, Tsunami, Jask Port, Numerical model, Inundation
1. Introduction
indescribable damage and casualties, paying more
serious attention to this phenomenon and to assessing
the vulnerability of different coasts around the world,
especially for marginal Indian Ocean countries like
Iran, is required. Definitely, Makran tsunami in 1945
and powerful tsunami of Indonesia in 2004 are not
the first or the last such occurrences in the region.
Indeed, the behaviour of three Eurasian, Indian and
Following the catastrophic Indian Ocean tsunami
in December of 2004, which resulted in fulmination
to more than 225,000 humans live and homelessness
of more than one million people (Geist, Titov et al.
2006) and 2011 Japan tsunami which caused
*
Email: akbarpour@inio.ac.ir
35
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
Arabian tectonic plates indicates that the recurrence
of tsunami in this zone is quite conceivable
(Ambraseys and Melville 1982). Tsunami hazard
assessment in any particular region and its irreparable
losses requires the compilation and analysis of past
tsunami records. Naturally, a better understanding of
the tsunami vulnerability of any region will not only
be crucial for all aspects of regional development,
including land-use and urban activities, but also for
estimating associated social and economical losses in
the aftermath of such events (Heidarzadeh et al. 2008).
Therefore, cataloguing the history of tsunami
occurrences around the world is a well-heeded
international practice.
Although tsunamis can be generated by variety of
different sources like earthquakes, submarine
landslides, volcanic eruptions and meteorite impact, but
subduction zones earthquakes are known as the most
common source of tsunamis. However, it should be
noted that usually earthquakes that take place along the
interface between the tectonic plates can lead to
tsunami formation. For example, the 2012 huge
earthquake occurred close to epicentre of the 2004
Indian Ocean earthquake was out of subduction zone
and did not cause considerable disturbance at ocean
surface; while the 2004 Indian Ocean earthquake
occurred inside the subduction zone, between Sunda
and Indian plates, generated huge tsunami that
devastated the surrounding region and caused about
230,000 losses of lives.
The Makran Subduction Zone (MSZ) is a
tsunamigenic narrow region extended from the Strait
of Hormuz at the end of Zindan-Minab fault in Iran
to the Ornach-Nal fault near Karachi in Pakistan. It is
stretching about 900 km (Heidarzadeh, et al. 2008)
due to reconstructions occurring among three
tectonic Eurasian, Indian and Arabian plates (Fig.1).
The rate of convergence along the Makran boundary
increases slightly from the west towards the east
(DeMets, Gordon et al. 1990), and the mean
subduction rate at the MSZ is about 19.5 mm/y
(Vernant, Nilforoushan et al. 2004). According to
other subduction rates related to different zones, such
as the Tonga Subduction Zone between Fiji and
Pacific plates equalling 160 mm/y (Bevis, Taylor et al.
2002), Japan Subduction Zone among Pacific and
Eurasian and Philippine plates at 80 mm/y (Kawasaki,
Asai et al. 2001), Sumatra Subduction Zone among
Australian, Indian and Eurasian plates amounting to
65 mm/y (Gahalaut and Catherine 2006), the MSZ is
characterised as a relatively slow moving subduction
zone. From the seismicity point of view, in the MSZ
large and shallow earthquakes are common which can
be regarded as a partial evidence for the tsunami
generation potential associated with order of
magnitude 8 earthquakes with an approximate average
return period of 100 to 250 years (Ambraseys and
Melville 1982; Heidarzadeh et al. 2008).
Fig. 1: Tectonic map of the MSZ at the northwestern Indian Ocean (Heidarzadeh et al., 2008)
36
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
Few references based on historical reports
document that quite a number of tsunamis have
occurred associated with different sources until 1945
A.D. in the MSZ. Heck (1947) provides a list of
global tsunamis generated from 479 B.C. to 1946
A.D., some of which were related to tsunami
occurrences in the Arabian Sea and adjacent regions.
Berninghausen (1966) presents another list of tsunami
events including three from the West coast of India
and the Arabian Sea. As an important tsunami
resource, Page (1979) describes existing evidence for
the recurrence of large-magnitude earthquakes along
the MSZ and predicts that the 1945-type earthquakes
can occur with a return period of about 150 to 250
years within the Eastern Makran. Ambraseys and
Melville (1982) submitted the results of studies related
to the historical earthquakes inside Iran. Another list
compiled by Murty and Rafiq (1991) used a variety of
sources in the Indian Ocean region from 326 B.C. to
1974 A.D., as well as Byrne et al. (1992), which
treated great thrust earthquakes along the plate
boundary of the MSZ. The threshold tsunami
generation potential of the region was studied by
Pararas-Carayannis (2006) along the MSZ.
The results of recent studies carried out based on
thermal modeling of the MSZ by Smith et al (2013),
shows that past assumptions may have significantly
underestimated the earthquake and tsunami hazard in
the MSZ; and it is potentially capable of producing
major earthquakes, up to Mw 8.7-9.2. Therefore, the
probability of recurrence of earthquake and tsunami
in the Makran zone is relatively high and the tsunami
can be considered as an obvious hazard for
neighbouring countries, like Iran and Oman.
Among the existing approaches to study tsunamis,
the use of Paleotsunami data, historical tsunami
records, online DART systems, local tide gages and
numerical modeling, using historical data and
numerical modeling due to being economically
advantageous are predominant to investigate tsunami
behaviour after occurrence. Numerical modeling of
tsunami in the MSZ provides insight into past events
and simulates possible future scenarios. Recently,
some numerical works have established in the MSZ to
study the generation and propagation of the tsunami
waves (Heidarzadeh et al., 2009; Heidarzadeh and
Kijko, 2011). In particular, the determination of the
potential run-up values and areas prone to inundation
via numerical modeling from a local or distant tsunami
is recognized useful, since data from past events is
usually insufficient (Rafi and Mahmood 2010).
The Jask Port is one of the most important coastal
areas of Iran, regarding its resident population and
strategic situation. It is located at the west of Makran
zone and can be affected by possible tsunamis. In
recent years, within the national development plans
of Iran, major ongoing port construction activities
have been undertaken for expansion of the Jask Port
as a portal on the transit corridor. Hence, the tsunami
hazard assessment is important in this region for
future development. In order to assess tsunami
hazard, the impact and run-up of the MSZ seismic
tsunamis on the Jask Port is simulated in this study.
However, the certainty of tsunami occurrence is
undeniable in the MSZ, and therefore, any
substantial investment in the development of the port
facilities is contingent upon a thorough investigation
of potential tsunami impacts.
The current study present the results of the
performed numerical modeling of tsunami propagation
as well as coastal zone inundation around the Jask due
to proximity to the MSZ based on evaluating three
plausible tsunami scenarios of varying consequences.
In the defined scenarios for tsunami modeling,
earthquakes having moment magnitude of 9.1Mw
and 8.7Mw are considered; which are equivalent to
the rupture of about full and half length of the plate
boundary in the MSZ. For this purpose, the open
source numerical model, GEOWAVE is applied as a
comprehensive numerical model. GEOWAVE is a
combination of TOPICS and FUNWAVE models, in
which the tsunami generation is simulated by
TOPICS and then tsunami propagation and run-up
are simulated by FUNWAVE. Unlike previous
simulations of Makran tsunamis (Rajendran et al,
2008; Neetu et al., 2011), containing just one largescale global model, in this study one high resolution
local model is defined around the Jask Port in order
to capture tsunami run up accurately.
37
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
2. Numerical Model
-4500
3500
-3000
2500
-2000
1500
-1000
-4.5-4000
-4 --3.5
-3 --2.5
-2--1.5
-1 -500
-0.5 00 5001000
0.5 1 1500
1.5 2000
2 2500
2.53000
3
28 N
-1000 -800
-600 -400
-200
0
200
400
(b)
(a)
25.7
Jask
26 N
25.6
25.5
24 N
25.4
25.3
22 N
56 E
58 E
60 E
62 E
64 E
66 E
68 E
57.6
57.7
57.8
Fig. 2: Bathymetry of studying domains for Global (a) and Local (b) models
38
-500
-300
-400
0
-100
Bathymetry (m)
-200
200
100
400
300
600
500
800
Bathymetry (km)
700
1000
Numerical modeling is an impressive approach to
studying the events of tsunami, particularly for the
future hazards. Numerical modeling of seismically
induced tsunamis may be organized into the following
three phases: generation, propagation and run-up
(inundation). The generation phase involves modeling
an underwater earthquake source and translates its
effects into an initial free surface water-level
displacement and an associated velocity field. The
second phase consists of propagating the generated
tsunami across Deep Ocean/Sea. Finally, inundation
associated with the run-up phase simulates the shallow
ocean zone behaviour of the tsunami in order to predict
coastal flooding and inundation. The interaction
between the waves and the shoreline is considered for
calculating wave propagation over dry land.
The achievement of tsunami studies depends on
the combination of precise tsunami sources, an
advanced tsunami propagation and finally inundation
model. As yet, many numerical models have been
modified or lately developed to impressively and
carefully study the generation, propagation, run-up
and inundation of tsunami. In this study, one of the
tsunami models, GEOWAVE, is adopted to study the
generation, propagation and inundation of potential
tsunami generated from MSZ. GEOWAVE is formed
by merging the Tsunami Open and Progressive Initial
Conditions System (TOPICS) with the Fully Nonlinear
Boussinesq water wave model FUNWAVE. An
900
accurate run-up and inundation by using a slot
method and wave breaking with a dissipative
breaking model were implemented in the model.
This model has capability of simulation of the whole
life range of tsunami from generation, propagation
till inundation. This model has been validated with
analytical solutions and experimental studies (Wei
and Kirby 1995; Wei et al. 1995; Chen et al. 2000;
Kennedy et al. 2000), and also benchmark problems
related to several historical events including the 1946
Unimak, Alaska; the 1946 Skagway, Alaska and the
1998 Papua New Guinea (Watts et al., 1999;Watts et
al. 2003, 2005; Waythomas and Watts 2003; Day et
al. 2005; Ioualalen et al. 2006; Grilli et al. 2007) and
indicated satisfactory achievements in comparison
with either field or satellite observations and tidal
gauge measurements. Figure 2 shows the global and
local domain of tsunami modeling. The data of the
Makran 1945 tsunami event authorize us to validate
the numerical model of the tsunami propagation.
Because of the lack of persistent sea-level
measurements of the tsunami, the sea-level
observations of the Makran 1945 tsunami from the
tide gauge at Karachi in Pakistan is persumed for
numerical model validation. The maximum wave
height of tsunami recorded at Karachi was 44 cm in
which compatible with the existent numerical model
that calculated 50 cm wave height at the same station.
Compatible with the reports by Neetu et al. (2011) and
Rajendran et al. (2008), the present model result
indicates the wave height of 250 cm in Pasni.
57.9
58
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
2.1. Model Grids
breaking, bottom friction and shoreline moving
boundary effects have been applied in this study. The
According to the long wavelength of tsunami waves,
especially in deep water which is usually in order of
tens of kilometers, generation and propagation of
waves can be simulated based on coarse hydrographic
data, typically in the order of kilometers; but accurate
calculation of tsunami run-up requires much finer
topography and hydrographic data since the behaviour
of tsunami waves strongly depends on the bathymetric
and topographic features near the coastal area.
governing volume conservation and momentum
equations are as follows,
(1)
distinguished
grids
of
bathymetric
(2)
in which (
) is water surface elevation, h(x,
y) is the local mean water depth;
is horizontal
gradient and M is depth-integrated horizontal volume
flux which is given by Eq. (3).
The present numerical modeling is performed in
two
)
(
In Eq.(3),
and
is horizontal velocity at an elevation
= −0.531h. Eq.(4) and Eq.(5)
taken to be
topographic data; The global type grid (refers to
represent dispersive effects,
global model) which covers the entire region of
related to the wave breaking and bottom friction are
Oman Sea where the fault plane is located, with 1-
explained by
minute grid resolution, and the local type grid (refers
functions of velocity given as Eq. (6) and Eq. (7).
to local model) covers the entire Jask port and its
(Grilli et al., 2007; Kennedy et al., 2000)
surrounding areas with 90 m grid resolution. The
Μ=
high resolution bathymetric and topographical data
was achieved from Community Model Interface for
2
α
(
V2 =
run-up and inundation in the Jask area for possible
tsunami events. The initial wave of tsunami is
1 2
(h
6
2
h At
2
Tsunami (ComMIT) database which is approved by
NOAA, and used to obtain the applicable details of
2
α
uα +
V1 =
and
+
1
+
2
.
generated in the global model and then it propagates
2)
h +
h
)2
+
α
2
2
1
) + (
2
( +
hA
1
+ (h
2
)
hB
At + Bt
2 )(u
α.
2
α
(3)
(4)
)
(5)
(6)
)
.(
from MSZ to the coastlines. The local model used
. The forces
respectively. A and B are
α h Bt
)(uα .
α
and
(7)
For moving boundary condition in the shoreline, a
when the waves reach to areas that have high
porous slot method was implemented by inserting
significance for calculation of tsunami run-up. When
factors
the waves enter to the vicinity of local models, the
and
as follows;
1 .
=
𝛿 + (1 𝛿)𝑒 𝜆(
global model is stopped and the local model is
begun. Free surface elevation and water velocities
which are achieved at the final time step of global
=
model, are constrained to the local model as the
initial conditions.
(
) + 𝛿(
𝛿( +
Here,
0)
+
+
(1
0)
+
𝛿)
𝜆
(1
0
𝑒 𝜆(
)/ 0
𝛿)
𝜆
.
0
)
1
0
1
𝑒
≥
≤
𝑒
𝜆(1+
𝜆(1+
(8)
0)
0)
.
.
≥
≤
(9)
is the slot layer depth which must be
deeper than the depth of maximum wave rundown
2.2. Governing Equations
during a simulation. The selection of
is explained
by Kennedy et al. (2000). According to the number
According to the Wei et al. (1995), the fully
of tsunami run-up events investigated by Watts et al.
nonlinear Boussinesq equations with considering
(2003), Day et al. (2005) and Grilli et al. (2007),
39
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
parameter values
=0.08 and
Wells and Coppersmith (1994) based on source
=25 are used for run-
up calculation.
parameters of 244 historical earthquakes. After initial
The major advantage of a Boussinesq wave
estimation, the source parameters are modified
propagation model over nonlinear shallow water
slightly to satisfy the relation between moment
wave model is considering the horizontal velocities
magnitude and seismic moment (Eqs. 10 and 11).
and also vertical acceleration variations along depth
2
Log10 ( M 0 ) 6
3
M 0 LW
Mw
(Watts et al., 2003). Moreover, the equations include
(10)
(11)
Where Mw is moment magnitude (dimensionless
eddy viscosity terms to model the dissipation caused
by wave breaking.
number), M0 is seismic moment in N m, is rigidity
of the earth (assumed 3×1010 Nm-2), L and W are
3. Tsunami Simulation for Different Scenario
Magnitudes
fault length and width respectively and is the
displacement on the fault surface.
3.1. Earthquake Scenario
3.2 Source Models
Earthquake scenarios (magnitude and source zone)
must be defined for any tsunami simulations. For the
The determined source parameters are inserted to
9.1Mw earthquake, the MSZ should experience the
TOPICS, tsunami generator model of GEOWAVE, and
full rupture to obtain relevant seismic moment; while
it calculates initial disturbance of ocean surface using
in order to produce the moment magnitude of 8.7Mw
Okada's elastic deformation source model (Okada,
about half of the plate boundary should be faulted. In
1985). Calculated free surface elevations actually are
the latter case, the western or eastern part of the MSZ
seabed vertical displacements transferred to ocean
can be faulted, according to the earthquake epicenter.
surface identically, because of instantaneous generation
Since estimation of future earthquake place is
of seismic tsunamis. As expected, maximum values of
impossible, the 8.7Mw earthquake is considered in
the initial wave crest and trough are observed during
intended scenarios both at the western and eastern
scenario 1 in which the 9.1Mw earthquake leads to
segments of the MSZ. So, for modeling of Makran
about +8m uplift and -6m subsidence of the ocean
tsunamis three scenarios are considered totally, as
surface (Fig.3a). Due to variation of the fault strike
presented in Table 1. In each scenario the source
between western and eastern part of the MSZ (Byrne et
parameters, defining dimensions of faulted area, are
al, 1992), the full rupture of the MSZ in scenario 1 is
estimated using empirical relationships presented by
divided to two separated segments.
Table 1: Characteristics and source parameters of scenarios for MSZ earthquakes
Scenario Earthquake moment
number
Faulted area
Source parameters
magnitude
Fault length (km)
Fault width
Fault displacement
(km)
(m)
1
9.1
All part of the MSZ
900
70
25
2
8.7
Western part of the MSZ
500
50
14
3
8.7
Eastern part of the MSZ
500
50
14
40
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
Free Surface Elevation (m)
-6 -4 -2 0 2 4 6 8 10 12 14
28 N
28 N
t=0
t=6min
26 N
26 N
24 N
24 N
22 N
22 N
(b)
(a)
58 E
60 E
62 E
64 E
66 E
68 E
28 N
58 E
60 E
62 E
64 E
66 E
64 E
66 E
68 E
28 N
t=12min
t=18min
26 N
26 N
24 N
24 N
22 N
22 N
(d)
(c)
58 E
60 E
62 E
64 E
66 E
68 E
58 E
60 E
62 E
68 E
Fig. 3: Snapshots of global tsunami simulations at times t = 0, 6, 12 and 18 min after the earthquake due to the Mw 9.1 magnitude
3.3. Global Model
but they amplify again approaching the coast due to
water depth reduction and shoaling effect. The waves
Free surface elevations obtained from TOPICS
propagate perpendicular to the faulted area mainly in
feature initial wave of tsunami and provide initial
condition
of
the
tsunami
propagation
the north-south direction, unlike landslide tsunamis
model,
propagate radial and identically in all directions
FUNWAVE; which is based on fully nonlinear
(Soltanpour and Rastgoftar, 2011). Figure 3 also shows
Boussinesq equations developed by Wei et al. (1995).
The
major
advantage
of
a
Boussinesq
that along the coastlines locating north of the MSZ
wave
falling water is first observed and then main waves of
propagation model over nonlinear shallow water wave
tsunami arrive; while southern coastlines of the MSZ,
models is that horizontal velocity variations along depth
like Oman, first receive tsunami positive waves.
and also vertical acceleration are considered (Watts et
al., 2003). Moreover, FUNWAVE equations include
3.4. Local Model
eddy viscosity terms to model the dissipation caused by
wave breaking. The governing equations are solved
Free surface elevation and water horizontal
first at global model grid points having 1-arc min
velocities of global model are transferred to the local
resolution. Thus, the propagation of tsunami waves
model as the initial conditions (Fig. 4). It should be
from source to near the coasts is simulated, as
noted that the initial conditions of global model also
illustrated in Fig. 3. It is observed that the tsunami
include water velocities, but they are set to zero;
waves dissipate firstly, at the beginning of propagation;
because no initial velocities are assume for the
41
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
earthquake tsunamis.
data bank. The time step of the local model is
After providing the required initial conditions, the
determined 0.23 s to satisfy the stability conditions; it
Boussinesq equations of FUNWAVE are solved at
is clear that the global model has greater time step
this time in the local model domain. The total
(2.43 s) according to larger mesh size. Finally, the free
number of grid points in the local model area is
surface elevation and horizontal velocities are
554,238 (751×738) according to 3-arc sec (about 90
calculated for all time steps and the run-up of tsunami
m) bathymetric data extracted from ComMIT model
on coasts of the Jask Port is computed (Fig. 5).
Free Surface Elevation (m)
-2
-1
0
1
2
3
4
V Velocity (m/s)
U Velocity (m/s)
5
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
-0.8 -0.4
25.75
25.75
25.65
25.65
25.55
25.55
25.45
25.45
25.35
25.35
0
0.4
0.8
1.2
1.6
2
2.4
25.7
25.6
25.5
25.4
25.3
25.25
25.25
57.6
57.7
57.8
57.9
58.0
(c)
(b)
(a)
57.6
57.7
57.8
57.9
58.0
57.6
57.7
57.8
57.9
58.0
Fig. 4: Local model initial conditions for scenario-1; (a) free surface elevation, (b) longitudinal velocity and (c) latitudinal velocity
Free Surface Elevation (m)
-10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16 18 20 22 24 26
t=20min
t=18min
t=22min
25.7
25.6
(a)
(b)
(c)
t=26min
t=24min
t=28min
25.7
25.6
(e)
(d)
57.7
57.8
57.9
57.7
57.8
57.9
(f)
57.7
57.8
57.9
Fig. 5: Tsunami induced inundation at Jask coasts for scenario-1 at times t= 18, 20, 22, 24, 26 and 28 min after the earthquake
42
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
4. Results and Discussion
coast of the port around the harbor, while the gauges
'S' and 'A' are at the eastern coast near the soccer
stadium and the airport, respectively (Fig. 6). Figure 7
shows the wave heights time series extracted during
the scenarios 1 and 2; while the time series of the
scenario 3 are presented in Fig. 8.
Three virtual gauges were deployed along the Jask
Port coastline to extract the time series of tsunami.
The location and water depth of the gauges are
presented in Table 2. The gauge 'H' located at western
Table 2: Locations of virtual tidal gauges
Gauge No.
Longitude (°E)
Latitude (°N)
Depth (m)
A
57.8000
25.6516
0.5
H
57.7672
25.6460
0.6
S
57.7720
25.6439
0.7
Fig. 6: Location of virtual tide gauges 'A', 'H' and 'S' for wave height time series analysis (filled circles)
Fig. 7: Modeled time series of water surface elevations at various locations along the Jask coast for the Mw 9.1 (solid line) and the Mw
8.7 western Makran (dashed line) tsunamis
43
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
Fig. 8: Modeled time series of water surface elevations at various locations along the Jask coast for the Mw 8.7 eastern Makran tsunami
Figures 7 and 8 demonstrate, as expected, the
maximum wave heights are observed in scenario 1 that
represents the largest earthquake. During scenario 3
tsunami waves appear at all the gauges much smaller
than two other scenarios. For example, the gauge 'S' has
recorded about 14 and 10-meter height waves in
scenario 1 and 2, respectively; but in scenario 3 the
maximum wave height observed in this gauge is less
than 0.4 meter. The reason is that when earthquake
occurs at the eastern part of the MSZ, the tsunami
waves reaches slightly along the western coasts of the
Makran zone, like Jask Port; because the tsunamis
caused by earthquakes spread perpendicular to the
source fault in the north-south direction in the MSZ.
Figure 7 reveals that in scenarios 1 and 2 the time
between earthquake occurrence and arrival of tsunami
waves is about 20 minutes. This shortness of time
makes the tsunami warning process for the Jask Port
very difficult. On the other hand, the negative wave of
tsunami is first observed along the Jask Port coastline
and the sea bottom becomes visible at shallow areas
when the negative wave height is more than the water
depth. This situation can be seen in Fig. 7 where there
are no data displayed, mostly at the beginning of the
time series. Residents of the port can use this natural
warning sign and be aware of the event before arrival
of the destructive waves. The figure also shows the
gauges 'S' and 'A' have recorded the waves having
approximately the same height, whereas the gauge 'H'
has experienced waves with smaller heights. This can
be attributed to the location of gauge 'H' that is at the
western coast of the port and so receives the diffracted
tsunami waves.
Maximum tsunami heights recorded during
scenarios 1 and 3 are presented in Fig. 9. This figure
also displays the tsunami inundation extent for all the
scenarios (The tsunami inundation in scenario 2 has
been specified by bold line). It is observed that the
waves due to the 8.7Mw eastern Makran earthquake
(scenario 3) cannot run-up on the Jask Port at all.
However, the waves caused by the 9.1Mw earthquake
(scenario 1) and the 8.7Mw western Makran
earthquake (scenario 2), not only inundate the port
completely, but also capture on farther inland areas.
Therefore, it can be concluded that the most crucial
factor determines the tsunami risk in the Jask Port is
the location of the Makran earthquakes. In the other
words, the huge earthquakes occur at east of the
MSZ can be considered as small hazard for the Jask
Port, while the occurrence of earthquakes having less
magnitudes at west of the MSZ leads to the high
waves devastate the port seriously. While it was
previously thought that the tsunami effect depends
essentially on the earthquake magnitude. According
to Fig. 9, the inundation extent of tsunami is a
function of the maximum tsunami height. In scenario
3, where no inundation takes place, the maximum
water level around the coastline of the Jask Port is
less than 1m; whereas it is between 10 m and 15 m in
scenario 1.
Although the most coastal areas around the Jask
Port are inundated by the tsunami waves in scenarios
1 and 2, but there are regions where are safe against
the tsunami waves. These regions are shown by circle
in Fig. 9. As adjacent areas of the Jask Port are
considered for development in future plans of the
Makran zone, the estimated safe regions seem
advisable regarding to tsunami threat.
44
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
Free Surface Elevation (m)
Latitude (deg.)
0
5
10
15
20
Free Surface Elevation (m)
0
25
25.75
25.75
25.70
25.70
25.65
25.65
(a)
25.60
57.75
57.8
57.85
57.9
0.2
0.4
0.6
0.8
1
(b)
25.60
57.95
1.2
57.75
57.8
57.85
57.9
57.95
Longitude
Longitude
Fig. 9: Maximum computed free surface elevation in the studying domain produced by (a) Mw 9.1 and (b) Mw 8.7 eastern Makran
earthquakes. Bold line represents the inundation line of Mw 8.7 western Makran earthquakes. The circles mark the regions are safe
against the tsunami waves
5. Conclusions
inundated by the tsunami waves. However, the
tsunami run-up is not observed at all when the
8.7Mw earthquake occurs at the eastern segment of
the MSZ. When 8.7Mw eastern Makran earthquake
occurs, the tsunami waves reach slightly along
western coasts of Makran zone; because seismic
tsunamis propagate perpendicular to the source fault,
in north-south direction in the MSZ. Therefore,
unlike previously thought, the tsunami risk at the
Jask Port depends more on the location of focal
points MSZ earthquakes than their magnitudes.
Largest possible earthquakes of the MSZ were
simulated to investigate tsunami hazards on the Jask
Port. The 9.1Mw and 8.7Mw earthquakes, featuring
respectively the rupture of about full and half length
of the MSZ's plate boundary, were assumed in defined
scenarios. The smaller earthquake was considered
both at the western and eastern segments of the MSZ.
GEOWAVE numerical model was employed to
simulate tsunami generation, propagation and
inundation. Unlike previous studies, containing just
one large-scale global model, in this study a smallscale local model was also defined based on 3-arc
sec resolution bathymetry data, in order to capture
tsunami run-up on the Jask Port.
The results show that during critical tsunamis
there is a short time between earthquake occurrence
and arrival of tsunami waves (about 20 minutes).
But, along the coastlines locating north of the MSZ,
like the Jask Port, falling water of tsunami was first
observed. Residents of the Jask Port can use this
natural warning sign and vacate the port before
arrival of the destructive waves.
Model results also reveal that due to the 9.1Mw
earthquake and the 8.7Mw western Makran earthquake
the most coastal areas around the Jask Port are
References
Ambraseys N.N. and Melville C.P., 1982. A History
of Persian Earthquakes, Cambridge University
Press, Cambridge, 219P
Berninghausen W.H., 1966. Tsunamis and Seismic
seiches reported from regions adjacent to the Indian
Ocean, Bulletin of Seismological Society of America
56(1):69–74
Bevis M., Taylor F. W., Schutz B. E., Recy J., Isacks B.
L., Helu S., Singh R., Kendrick E., Stowell J., Taylor
B., and Calmantli S., 2002. Geodetic observations of
very rapid convergence and back-arc extension at the
tonga arc, Nature, 374, 249 – 251
Byrne D., Sykes L.R. and Davis D.M., 1992. Great
45
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
thrust earthquakes and aseismic slip along the
plate boundary of the Makran subduction zone,
Journal of Geophysical Research. 97:449–478
Chen Q., Kirby J.T., Dalrymple R.A., Kennedy A.B.
and Chawla A., 2000 Boussinesq modeling of wave
transformation, breaking, and run up. II: 2D, Journal
of Waterway Port Coastal and Ocean Engineering
126(1): 48–56
Day S.J., Watts P., Grilli S.T. and Kirby J.T., 2005
Mechanical models of the 1975 Kalapana, Hawaii
earthquake and tsunami. Marine Geology 215(1–2):
59–92
DeMets C., Gordon R. G., Argus D. F. and Stein S.,
1990. Current plate motions, Geophysical Journal
International, 101, 425-478
Geist E.L., Titov V.V. and Synolakis C.E., 2006
Tsunami: Wave of change, Scientific American
56-63
Gahalaut V. K. and Catherine J. K., 2006, Rupture
characteristics of 28 march 2005 Sumatra earthquake
from GPS measurements and its implication for
tsunami generation, Earth Planet Science Letter, 249:
39-46
Grilli S.T., Ioualalen M., Asavanant J., Shi F., Kirby
J.T. and Watts P., 2007. Source Constraints and
Model Simulation of the December 26, 2004,
Indian Ocean Tsunami, J. Waterway Port Coastal
Ocean Eng 133(6):414-428
Heck N.H., 1947. List of seismic sea wave, Bull
Seismol. Soc. Am 37(4):269–286
Heidarzadeh M., Dolatshahi P., Hajizadeh Zaker N. and
Yalciner A.C., 2008 Historical tsunamis in the
Makran subduction zone off Southern Coasts of Iran
and Pakistan and results of preliminary numerical
modeling, Ocean Engineering 35:774–786
Heidarzadeh M., Dolatshahi P., Hajizadeh Zaker N.,
2009. Modeling the near-field effects of the worstcase tsunami in the Makran subduction zone,
Ocean Engineering. 36(5): 368–376
Heidarzadeh M. and Kijko A., 2011. A probabilistic
tsunami hazard assessment for the Makran
subduction zone at the northwestern Indian Ocean,
Natural Hazards 56(3):577-593
Ioualalen M., Pelletier J., Watts P. and Regnier M.,
2006 Numerical modeling of the 26th November
1999 Vanuatu tsunami, Journal Geophysical
Research 111 C06030
Kawasaki I., Asai Y., and Tamura Y, 2001. Space–
time distribution of interplate moment release
including slow earthquakes and the seismo-geodetic
coupling in the sanriku-oki region along the japan
trench, Tectonophysics, 330, 267-283
Kennedy A.B., Chen Q., Kirby J.T. and Dalrymple
R.A., 2000 Boussinesq modeling of wave
transformation, breaking, and run up. I: 1D, Journal
of Waterway Port Coastal Ocean Eng 126(1):39–47
Murty T. and Rafiq M., 1991. A tentative list of
tsunamis in the marginal seas of the north indian
ocean, Natural Hazards. 4:81-83
Neetu S., Suresh I., Shankar R., et al., 2011 Trapped
waves of the 27 November 1945 Makran tsunami:
observations and numerical modeling, Natural
Hazards DOI 10.1007/s11069-011-9854-0
Okada Y., 1985. Surface deformation due to shear
and tensile faults in a half-space. Bull Seismol Soc
Am 75(4):1135–1154
Page W.D., Alt J.N., Cluff L.S. and Plafker G., 1979.
Evidence for the recurrence of large magnitude
earthquakes along the Makran Coast of Iran and
Pakistan, Tectonophysics 52:533–547
Pararas-Carayannis G., 2006. The potential for
tsunami generation along the Makran subduction
zone in the northern Arabian Sea. Case study: the
earthquake and tsunami of November 28, 1945,
Sci Tsunami Hazards 24(5):358–384
Rafi Z., Mahmood N., 2010. Numerical modeling of
tsunami inundation for potential earthquake at
Makran subduction zone - a case study for Gwadar
coastal area, Pakistan Meteorological Department,
Islamabad PMD-45
Rajendran C.P., Ramanamurthy M.V., Reddy N.T.
and Rajendran K., 2008. Hazard implications of
the late arrival of the 1945 Makran tsunami,
Current Science, 95(12):1739-1743
46
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
IUGG Birmingham UK JSS42 Abstract
Watts P., Grilli S.T., Kirby J.T., Fryer G.J. and
Tappin D.R., 2003 Landslide tsunami case studies
using Boussinesq model and a fully nonlinear
tsunami generation model, Natural Hazards and
Earth Science 3:391–402
Watts P., Ioulalen M., Grilli S.T., Shi F., Kirby J.T.,
2005 Numerical simulation of the December 26,
2004 Indian Ocean tsunami using a higher-order
Boussinesq model, Proc 5th Int Conf on Ocean
Wave Measurement and Analysis, WAVES 2005
Madrid Spain Paper No. 221
Waythomas C.F. and Watts P., 2003, Numerical
simulation of tsunami generation by pyroclastic
flow at Aniakchak Volcano, Alaska, Geophys Res
Lett 30(14):1751–1755
Wei G. and Kirby J.T., 1995 Time-dependent
numerical code for extended Boussinesq equations,
Journl of Waterway Port Coastal and Ocean
Engineering 21 (5): 251–261
Wei G., Kirby J.T., Grilli S.T. and Subramanya R.,
1995. A fully nonlinear Boussinesq model for free
surface waves. Part 1: Highly nonlinear unsteady
waves, Journal of Fluid Mechanics 294: 71–92
Soltanpour, M. and Rastgoftar, E., 2011. Study of
tsunami attacks on neighboring countries of
Caspian Sea caused by a probable submarine
landslide", Journal of Coastal Research, SI 64: 1195
– 1199. Szczecin, Poland
Smith G.L., McNeill L.C., Wang K., He J. and
Henstock T.J., 2013. Thermal structure and
megathrust seismogenic potential of the Makran
subduction zone, Geophys Res Lett. doi: 10. 1002/
grl. 50374
Vernant P.h., Nilforoushan F. and Hatzfeld D. et al.,
2004. Present-day crustal deformation and plate
kinematics in the middle east constrained by GPS
measurements in Iran and Northern Oman,
Geophysical Journal International 157:381–398
Wells D.L. and Coppersmith K.J., 1994 New
empirical relationships among magnitude, rupture
length, rupture width, rupture area, and surface
displacement, Bulletin of the Seismological
Society of America 84 (4), 974–1002
Watts P., Borrero J.C., Tappin D.R., Bardet J.P. and
Grilli S.T., et al., 1999 Novel simulation technique
employed on the 1998 Papua New Guinea
tsunami, Proceedings of 22nd General Assembly
Akbarpour Jannat and Rastgoftar / Numerical Modeling of Tsunami Waves Associated…
Journal of the Persian Gulf (Marine Science)/Vol .6/No .22/December 2015/13/35-48
Journal of the Persian Gulf
(Marine Science)/Vol. 6/No. 22/December 2015/13/35-48
47
48