Jimenez, C. et al.

Paper:

Seismic Source of 1746 Callao Earthquake fromTsunami Numerical Modeling
Cesar Jimenez∗1,∗2 , Nabilt Moggiano∗2 , Erick Mas∗3 , Bruno Adriano∗3 ,Shunichi Koshimura∗3 , Yushiro Fujii∗4 , and Hideaki Yanagisawa∗5∗1 Fenlab, Universidad Nacional Mayor de San Marcos (UNMSM)

Av Venezuela s/n, Lima, Peru
E-mail: cjimenezt@unmsm.edu.pe∗2 Direccion de Hidrografia y Navegacion (DHN)
Calle Roca N 116, Chucuito-Callao, Peru
∗3 Laboratory of Remote Sensing and Geoinformatics for Disaster Management,International Research Institute of Disaster Science, Tohoku University
Aoba 6-6-03, Sendai 980-8579, Japan∗4 International Institute of Seismology and Earthquake Engineering, Building Research InstituteTatehara, Tsukuba, Ibaraki 305-0802, Japan∗5 Department of Regional Management, Faculty of Liberal Arts, Tohoku Gakuin University
2-1-1 Tenjinzawa, Izumi-ku, Sendai, Miyagi 981-3193, Japan[Received November 2, 2012; accepted February 8, 2013]

In this paper a model of slip distribution is proposed
for the 1746 Callao earthquake and tsunami based on
macroseismic observations written in historical documents. This is done using computational tools such as
tsunami numerical simulation through a forward process by trial and error. The idea is to match historicalobservations with numerical simulation results to obtain a plausible seismic source model. Results showa high asperity from Cañete to Huacho, which would
explain the great destruction in this area. The rupture
directivity of the seismic source, from north to south,
would explain the value of the arrival time of the firsttsunami wave at Callao. A kinematic seismic source
model was used as a first approximation of the event.
The estimated magnitude was Mw9.0.
Keywords: seismic source, tsunami, simulation

1. Introduction
Many researchers have tried to infer seismic sources ofhistorical events such as the 1755 Lisbon earthquake [1,2], the 1693 Catania earthquake [3], and the 1700 Cascadia earthquake [4], etc., from macroseismic observations
recorded in historical documents and from observationsof marine geology and paleo-tsunami studies [5, 6], using
tsunami numerical modeling techniques. Other authors
such as Dorbath [7], using only macroseismic information, have qualitatively estimated the magnitude and rupture area size.
At the present time, we have suitable numerical tools
and geophysical instrumentation such as broadband seismic stations, digital tide gauge stations, geodetic GPS andradar interferometry INSAR for studying seismic sources266

and crustal deformation. This is not the case for historical events such as the 1746 earthquake, however. In
this sense, we can only infer or estimate a seismic source
model from macroseismic and tsunami descriptions ofhistorical documents found in the literature.

2. Historical Aspects
Silgado [8], among other authors, performed a compilation of the historical aspects of this event. Below is a summary of the historical account provided by
Soloviev [9]:The tragedy began on October 28, 1746, at 22:30 pm.
To the north, the quake was felt as far as Guayaquil,
1100 km from Lima, and a Jesuit mission locatednear the confluence of the Marañón and Huallaga
rivers, 750 km from Lima. In Huancavelica, southsouth-east of Lima, there was severe thunder accompanying movement. Buildings south of Lima collapsed all the way to Cañete, and north to Huacho
about 120 km north of Lima, where a new bridgeover the river collapsed at Huaura. Roads leading
to the interior were blocked by landslides. The valleys of Supe, Barranca and Pativilca were seriously
damaged by the earthquake and tsunami. Movement
was also felt in Cusco, Tacna and border towns. Aftershocks were felt in Lima overnight and within
24 hours, no fewer than 200 events were counted.
The port of Callao suffered the most damage fromthis catastrophe. The city was enclosed by a wallwhose base was sometimes beaten by the tide andwaves. The population reached 5,000 people. Halfan hour after the quake, the sea rose 10 m and, in
Journal of Disaster Research Vol.8 No.2, 2013

Seismic Source of 1746 Callao Earthquake from Tsunami Numerical Modelingits advance, broke the walls and inundated the city.
When the sea receded, most of the houses and buildings were taken from their foundations and washedaway. Most of the walls of the city, including thegates, were washed away. Almost all of the inhabitants of the city perished in this disaster. The water rose about 4 km inland, reaching even those who
tried to escape to Lima. Only 200 people were saved,
and these by clinging to wooden objects, but werethrown between the coast and San Lorenzo island, adistance of up to 8 km. Of the 23 ships anchoredin the harbor, 19 were sunk and 4 were carried inland. At around 04:00 of the next day, Callao wasinundated again by another tsunami wave. The maximum inundation height was estimated at 24 m.
In the port of Santa (see Fig. 1), waves hit the ship
Concepcion with such force that it sprang a leak andsank. The crew of the Soledad, which was near
Nazca, noted that the sea had retreated, and took precautions, so the ship was saved. The tsunami destroyed the ports of Caballas (Ica), Pisco, Chancay
and Guañape (La Libertad). Near Huacho, the roadwas completely flooded and vehicles that stayed on
the road, together mules, were swept away when thesea receded. In Huaura (Salinas), the sea flooded areas around 4 km inland, drowning mules and drivers.
Near Callao, after the earthquake, part of the coastline suffered subsidence, to where a new bay wasformed. In the space of six hours after the inundation
at Callao, the tsunami reached Concepción, Chile;and at Acapulco, Mexico, a ship was washed ashore.
From historical records of authors cited in the literature [10], we found the following important facts, which
will be useful in defining seismic source through a forward process based on trial and error. Note that in doingso, parameters of the simulated tsunami will need to fit
these historical observations.
a) The arrival time in Callao of the first tsunami wave
was at around 30 minutes.
b) The maximum height of the first tsunami wave wasabout 10 m at Callao.
c) The maximum inundation height (runup) was 24 m
(at the cliffs of the Costa Verde zone), see paragraph
5.3 and Fig. 6.
d) After the earthquake, a permanent subsidence at
Callao was reported.
e) Total destruction stretched from Cañete in the south
to Barranca and Pativilca in the north. This implies
large asperity or slip between these locations.
f) Maximum horizontal inundation was 4 km probablyat the old Rimac River bed.
g) Tsunami waves arrived at Callao from the northwest.
See Fig. 2, from a picture of the period [11].
Journal of Disaster Research Vol.8 No.2, 2013

Fig. 1. Location of ports and resorts affected by the tsunami.

Fig. 2. 1746 Tsunami, from a picture of the period [11].

h) According to Silgado [8], the intensity was X-XI
MM in Lima and Callao. It is possible that the occurrence of some cracks in the main square in Lima
would imply very great acceleration.
Due to the lack of more historical data to thoroughly
confirm the testimonies in the record, we will, at thisstage, use available information as the main constraint tobe satisfied by the seismic source model proposed here.

3. Data Acquisition
3.1. Topographic and Bathymetric Data
To prepare digital elevation model data, satellite andsurvey data topographic and bathymetric was acquired.
These data were digitally processed to obtain a digital elevation model with a grid resolution of 30 m (inundationgrid) and 900 m (generation and propagation grid), as wellas intermediate grids. Table 1 specifies the resolution of267

Jimenez, C. et al.

Table 1. Characteristics of digital elevation models.
Model
SRTM 90
DHN survey
Gebco 30

Type
Topography
BathymetryTopo-bathymetry

Resolution
90 m
100 m
900 m

each topographic and bathymetric model.
The DHN (Dirección de Hidrografı́a y Navegación)bathymetry and SRTM topography have been processed(through two-dimensional interpolation) to obtain a fineresolution digital elevation model of 30 m in raster or matrix format to be used in the inundation phase for the numerical simulation. To develop codes and routines, weused Fortran and Matlab programming language, libraries
and Mapping Toolbox utilities as well as Surfer 7.0 software. For application of the linear model of tsunami propagation, we used Gebco30 bathymetry.

4. Methodology
4.1. Forward Method by Trial and Error
If we know the parameters that characterize a system
and its initial and boundary conditions, then we can predict the behavior of the system at any instant in time. This
problem of predicting future behavior, is called the “forward problem.” The “inverse problem” is to use the results of certain observations: tidal, field data, and paleotsunami deposits, etc. to infer values of parameters characterizing a system [12].
In this paper we use information from historical documents to estimate parameters of a seismic source (slip
distribution). The solution of the forward problem is notunique, but chooses that set of values providing the best
correlation between observed and simulated data. The
trial and error method involves varying parameters in thiscase, slip and depth and performing the numerical modeland comparing the output with historical data, then repeating the process as often as necessary to achieve a roughcorrelation.
4.2. Constraints on the Size of the Source and Magnitude of the Earthquake
According to the seismotectonic pattern, the size ofseismic source in the central region of Peru is restrictedby the Nazca ridge (15◦ S) to the south and the Mendaña’s
fracture (10◦ S) to the north. This region covers a distancearound 600 km. According to the empirical relation of
Papazachos [13]:log(L) = 0.55M − 2.19

. . . . . . . . . (1)where L is in km. Magnitude would have a maximum
limit of Mw9.0. Furthermore, the fact that it caused cracks
in the main square in Lima indicates that there was strong
acceleration due to an event of great magnitude, i.e., ofat least Mw9.0. Dorbath [7] estimated a magnitude ofMw 8.6 and 9.0 Mt, Beck [14] estimated Mw8.8, and268

Table 2. Focal mechanism parameters.
Focal mechanism
Strike angle
Dip angle
Slip angle

Angle
329◦
20◦
90◦

Table 3. Geometry parameters of the fault.
Parameter
Longitude
Width
Mean slip
Depth

550 km140 km11.5 m
8.0 kmOcola [15] calculated Mw 9.2 based on macroseismic intensity.
This paper proposes a magnitude of Mw9.0 as an initialassumption. Later, it will be shown that this magnitude
fits results from tsunami numerical simulation.

4.3. Focal Mechanism and Fault Parameters
Focal mechanism parameters (see Table 2) are estimated according to following assumptions: the azimuthal
angle should be parallel to the marine trench in central
Peru (about 329◦ ). The dip angle for the central regionof Peru is 20◦ (average dip for CMT solutions) and thedisplacement angle (rake or slip angle) is estimated in 90◦
for maximum vertical deformation. These are the fixedparameters in the proposed method.
According to empirical equations relating fault dimensions to magnitude, fault dimensions would be for a magnitude Mw9.0: L = 575 km and W = 144 km. However,in this paper we set these dimensions to: L = 550 km and
W = 140 km. According to the relation between scalar
seismic moment and moment magnitude:log(M0 ) = 1.5Mw + 9.1. . . . . . . . . (2)= 3.98 × 1022 Nm, using the following def-

We have: M0
inition of scalar seismic moment:
M0 = μ LW D

. . . . . . . . . . . . . (3)where: μ = 4.5 × 1010 N/m2 , is the average rigidity of theelastic medium. After calculating the dislocation value ormean slip D, we have: D = 11.5 m. The geometry of the
rupture area must be heterogeneous composed of several
fragments to represent the maximum asperity in the areafrom Cañete to Barranca. The depth of the top of the fault
is taken initially as: H = 5 km, and after the trial and error
method used, is set finally to H = 8 km (see Table 3).

5. Tsunami Numerical Modeling5.1. Seismic Source Model
According to Perfettini [16], the Nazca ridge appears to
be a barrier aseismic to earthquake rupture. Because no
Journal of Disaster Research Vol.8 No.2, 2013

Seismic Source of 1746 Callao Earthquake from Tsunami Numerical ModelingFig. 3. Fragments of the seismic source (slip distribution).
Table 4. Slip distribution for each of 5 fragments. Coordinates correspond to the lower left corner and depth to theupper side along the vertical axis faults.
N◦
1
2
3
4
5

Lat (◦ )-14.2000-13.3520-12.5040-11.6561-10.8082

Lon (◦ )-77.4200-77.9295-78.4390-78.9485-79.4580

H j (km)
8.00
8.00
8.00
8.00
8.00

Slip (m)
09.00
17.50
17.00
07.00
07.00

great historical earthquakes have been reported in northern Peru since at least the 17th century, it is reasonable
to assume that great earthquakes at Lima correspond to arupture between latitudes 15◦ S to 9◦ S. We have chosen afault length L = 550 km, a distance encompassing to theselatitudes.
We have assumed that Peruvian seismicity distribution
is located on a fault plane from the trench. We have chosen a dip angle for subduction plane δ = 20◦ corresponding to mean dip of CMT solutions in the selected area.
Width W of the fault in the direction of the dip angle ischosen as W = 140 km and corresponds to depth limit
W sin(δ ) + H of 58 km, a depth consistent with the extension of seismicity with depth. The lower limit of the horizontal projection of the fault corresponds approximately
to the position of the coast, often seen to be the limit ofgreat subduction earthquakes.
To simulate historical tsunami information we havedivided the total rupture area into five fragments (seeFig. 3). Slip distribution is obtained by the forward
method by trial and error until the model agrees with historical observations. The present model may be appropriate for explaining available tsunami information, but notnecessarily the observed intensities due to shaking of the
Journal of Disaster Research Vol.8 No.2, 2013

Fig. 4. Initial deformation of the seabed [19].

ground. This slip distribution is given in Table 4, wherein
the mean slip value corresponds to 11.5 m (see Figs. 3 and4).
To calculate initial coseismic deformation of our source
model (see Fig. 4) we used the model of Okada [17],which will be used to set the initial condition of thetsunami propagation phase.

5.2. Tsunami Propagation Model
We have used the numerical model TUNAMI [18],which numerically solves using a finite difference scheme(the leapfrog method) equations of conservation of linear
momentum and the continuity equation of fluid dynamics
to simulate the propagation (linear model in spherical coordinates) and inundation (nonlinear model in Cartesian
coordinates for inundation grid) of the tsunami, with theinitial condition of coseismic deformation given by theseismic source model.
To fit the requirement of “tsunami waves arriving at
Callao after half an hour from the northwest,” we musttake into account a source rupture directivity from northto south for a total rupture time of about 3 minutes(vr = 3 km/s). This implies that the rupture could startfrom fault 01 (north) consecutively until fault 05 (south),(Fig. 5). This is in agreement with the North-West to
South-East unilateral rupture propagation that has characterized most great earthquakes in Peru, such as the 1996
Nazca, 2001 Camaná and 2007 Pisco earthquakes. This
probably means that the rupture began north of the rupture area and then propagated south-east.

269

Jimenez, C. et al.

Fig. 5. Kinematic model of the seismic source for the 1746 Callao earthquake and tsunami.

Fig. 6. Virtual tidal gauges. Notice lower saturation due tothe proximity of virtual gauges to the shoreline. The heightof the first wave in Callao was 10 m while the maximum
runup was 22 m at Costa Verde.

of La Punta, the first wave was less than 10 m, while in
the area of the Costa Verde, the height of the first wave
was greater than 22 m. These differences are due to thebathymetry, topography and morphology of each particular area.
We also observe a peculiarity of the coseismic subsidence phenomenon that started after the earthquake, i.e.,a decrease in the sea level as evidenced by virtual gauges
before the arrival of the first tsunami wave. This is due tothe geometry of the subduction zone and the focal mechanism of the earthquake.
Figure 7 shows an inundation map for Callao. Inundation would arrive at 4 km from La Punta to the landmarkof Plaza Vea. In the airport zone, inundation would arrive
at 1.8 km inland, the airport is located 3.6 km inland, so
the airport would not be affected. It is also noted that at
Costa Verde, waves may be more than 22 m high, whichis consistent with historical information.

6. Conclusions
5.3. Tsunami Inundation Model
For simulation of this phase, we have used a non-linear
model to take into account friction between fluid andground through the Manning roughness coefficient. As
a result of non-linear simulation of tsunami dynamics, wehave obtained parameters such as arrival time and wave
height.
Figure 6 shows that in a virtual gauge, adjacent to the
Lima port area (Callao DHN), the first wave had a heightof 10 m with an arrival time of 23 minutes. In the area
270

The seismic source parameters modeled in this study
are derived from historical documents and informationfrom seismic catalogs based on macroseismic parameters.
The magnitude of 1746 Callao earthquake is estimatedat Mw9.0, a value obtained quantitatively by the forward
process of trial and error through numerical modeling oftsunamis.
Heterogeneous initial deformation or slip distribution
indicates greater asperity from Cañete to the south to Barranca to the north. This agrees with the total destruction
Journal of Disaster Research Vol.8 No.2, 2013

Seismic Source of 1746 Callao Earthquake from Tsunami Numerical ModelingFig. 7. Tsunami inundation map for Callao and Miraflores Bay. Note the height (red) reached by the tsunami
waves in Costa Verde. Blue triangles correspond to virtual gauges: M1 = Callao-DHN, M2 = La Punta, M3 =Costa Verde.

reported in this region.
For great earthquakes or for very slow rupture velocity,
it is necessary to apply a kinematic source model. Thisaffects the values of tsunami arrival time, depending on
the directivity of the rupture process.
For a near field tsunami at Lima and Callao, the arrival
time of the first wave would be 23 minutes (at Callao)
according to the nonlinear inundation model.
According to the numerical model, the height of thefirst wave at virtual gauge M1 (Callao DHN) is 10 m. The
maximum vertical inundation height or run-up is 22 m at
virtual gauge M3 (Costa Verde), (Fig. 7).
Results obtained by this paper correspond to a numerical model and, to some extent, are validated by data from
historical documents.
Acknowledgements
The authors acknowledge support from: the project of the Scienceand Technology Research Partnership for Sustainable Development supported by JST-JICA, Dirección de Hidrografı́a y Navegación of the Peruvian Navy and the Universidad Nacional Mayorde San Marcos - Fenlab.

References:[1] M. Baptista, P. Miranda, J. Miranda, and V. Mendes, “Constrainson the source of the 1755 Lisbon Tsunami inferred from numerical modelling of historical data on the source of the 1755 Lisbon Tsunami,” Journal of Geodynamics, Vol.25, No.2, pp. 159-174,
1998.
[2] A. Santos, S. Koshimura, and F. Imamura, “The 1755 Lisbon
Tsunami: Tsunami Source determination and its validation,” Journal of Disaster Research, Vol.4, No.1, pp. 41-52. 2009.
[3] M. Gutscher, J. Roger, M. Baptista, J. Miranda, and S. Tinti,“Source of the 1693 Catania earthquake and tsunami (southern
Italy): new evidence from tsunami modeling of a locked subduc-

Journal of Disaster Research Vol.8 No.2, 2013tion fault plane,” Geophysical Research Letters, Vol.33, L08309,
2006.
[4] K. Satake, K. Wang, and B. Atwater, “Fault slip and seismic
moment of the 1700 Cascadia earthquake inferred from Japanese
tsunami descriptions,” Journal of Geophysical Research, Vol.108,
No.B11, 2003.
[5] D. Sugawara, F. Imamura, H. Matsumoto, K. Goto, and K. Minoura,“Reconstruction of the AD 869 Jogan earthquake induced tsunami
by using the geological data,” Journal of Japanese Society of Natural Disaster Science, Vol.29, pp. 501-516, 2011.
[6] K. Goto, C. Chagué-Goff, S. Fujino, J. Goff, and B. Jaffe, Y.
Nishimura, B. Richmond, et al., “New insights of tsunami hazard
from the 2011 Tohoku-oki event,” Marine Geology, Vol.290, pp. 4650, 2011.
[7] L. Dorbath, A. Cisternas, and C. Dorbath, “Assessment of the size oflarge and great historical earthquakes in Peru,” Bulletin of Seismological Society of America, Vol.80, No.3, pp. 551-576, June 1990.
[8] E. Silgado, “Historia de los Sismos más notables ocurridos en el
Perú (1513-1974),” Boletı́n, No.3, Instituto de Geologı́a y Minerı́a,Lima, Perú.
[9] S. Soloviev and Ch. Go, “Catalogue of Tsunamis on the Eastern
Shore of the Pacific Ocean (1513-1968),” Nauka Publishing House,p. 204, Moscow, USSR, 1975.
[10] L. Seiner, “Historia de los Sismos en el Perú: Catálogo Siglos
XVIII-XIX,” Fondo Editorial Universidad de Lima, Perú 2011.
[11] E. Ortega, C. Jiménez, and N. Moggiano, Revista de Tsunamis, Dirección de Hidrografı́a y Navegación, Callao, Perú, 2012.
[12] A. Tarantola (Ed.), “Inverse problem theory and methods for model
parameter estimation,” Society for Industrial and Applied Mathematics, Philadelphia, 2002.
[13] B. Papazachos, “Global relations between seismic fault parametersand moment magnitude of Earthquakes,” Bulletin of the GeologicalSociety of Greece, Vol.XXXVI, pp. 1482-1489, 2004.
[14] S. Beck and S. Nishenko, “Variations in the mode of great earthquake rupture along the central Peru subduction zone,” GeophysicalResearch Letters, Vol.57, pp. 1969-1972, 1990.
[15] L. Ocola, “Catálogo Sı́smico de la Republica del Perú,” SISAN, IV,
Lima, Perú, p. 466, 1984.
[16] H. Perfettini, J. P. Avouac, H. Tavera, A. Kositsky, J. M. Nocquet,
F. Bondoux, M. Chlieh, A. Sladen, L. Audin, D. L. Farber, and P.
Soler, “Seismic and aseismic slip on the Central Peru megathrust,”
Nature, Vol.465, pp. 78-81, doi:10.1038/nature09062, 2010.
[17] Y. Okada, “Surface deformation due to shear and tensile faults in a
half space,” Bulletin of Seismological Society of America, Vol.75,
No.4, pp. 1135-1154, 1985.

271

Jimenez, C. et al.

[18] F. Imamura, “Review of Tsunami Simulation with a Finite Difference Method, Long Waves Runup Models,” World Scientific Publishing Co. Pte. Ltd., Singapore, 1996.
[19] A. Bruno, E. Mas, S. Koshimura, Y. Fujii, S. Yauri, C. Jimenez, andH. Yanagisawa, “Tsunami Inundation Mapping in Lima, for Two
Tsunami Source Scenarios,” Journal of Disaster Research, Vol.8,
No.2, pp. 274-284, 2013 (this number).

Name:
Erick Mas Samanez

Affiliation:
Assistant Professor, Laboratory of Remote Sensing and Geoinformatics for Disaster Management (ReGiD), International Research Instituteof Disaster Science (IRIDeS), Tohoku University

Name:
Cesar Jimenez Tintaya

Affiliation:

Address:

Fenlab, Universidad Nacional Mayor de San
Marcos (UNMSM)
Dirección de Hidrografı́a y Naveagacı́on (DHN)

Aoba 6-6-3, Sendai 980-8579, Japan

Address:Av. Venezuela s/n, Lima, PerúCalle Roca 116, Chucuito-Callao, Perú

Brief Career:
2000-2007 Assistant Research, Instituto Geofisico del Perú IGP
2008-2012 Research Scientist, Dirección de Hidrografı́a y Naveagacı́on
2009-2012 Assistant Professor, UNMSM, Perú
2012- Graduate Studies in Geophysics, UNMSM, Perú

Selected Publications:• M. Ioualalen, H. Perfettini, S. Yauri, C. Jimenez and H. Tavera,“Tsunami Modeling to Validate Slip Models of the 2007 Mw8.0 Pisco
Earthquake, Central Peru,” Pure and Applied Geophysics. 2012.
• E. Mas, S. Koshimura, A. Suppasri, M. Matsuoka, M. Matsuyama, T.
Yoshii, C. Jimenez, et al. (2012). “Developing Tsunami fragility curves
using remote sensing and survey data of the 2010 Chilean Tsunami in
Dichato,” Natural Hazards and Earth System Science, 12, 2689-2697.
doi:10.5194/nhess-12-2689-2012.
• T. Yoshii, M. Imamura, M. Matsuyama, S. Koshimura, M. Matsuoka, E.
Mas and C. Jimenez. “Salinity in Soils and Tsunami Deposits in Areas
Affected by the 2010 Chile and 2011 Japan Tsunamis,” Pure and AppliedGeophysics, 2012.

Brief Career:
1999-2004 B.S. Civil Engineering, National University of Engineering,Peru
2006-2009 M.Sc. Disaster Risk Management, National University ofEngineering, Peru
2009-2012 PhD Civil Engineering, Tsunami Engineering, Tohoku
University, Japan
2012- Assistant Professor, ReGiD, IRIDeS, Tohoku University, Japan

Selected Publications:• E. Mas, A. Suppasri, F. Imamura, and S. Koshimura, “Agent Basedsimulation of the 2011 Great East Japan Earthquake Tsunami evacuation.
An integrated model of tsunami inundation and evacuation,” Journal of
Natural Disaster Science, Vol.34, Iss.1, pp. 41-57, 2012.
• E. Mas, S. Koshimura, A. Suppasri, M. Matsuoka, M. Matsuyama, T.
Yoshii, C. Jimenez, F. Yamazaki, and F. Imamura, “Developing Tsunami
fragility curves using remote sensing and survey data of the 2010 Chilean
Tsunami in Dichato,” Natural Hazards and Earth System Science, Vol.12,
pp. 2689-2697, 2012, doi:10.5194/nhess-12-2689-2012.

Academic Societies & Scientific Organizations:• Japan Geoscience Union (JpGU)• European Geosciences Union (EGU)• American Geophysical Union (AGU)

Name:

Academic Societies & Scientific Organizations:

Bruno Adriano Ortega• Peruvian Society of Physics (SOPERFI)

Affiliation:
Research Student, Laboratory of Remote Sensing and Geoinformatics for Disaster Management (ReGiD), International Research Instituteof Disaster Science (IRIDeS), Tohoku University

Name:
Nabilt Moggiano

Affiliation:
Geophysical Division, Oceanography Department, Dirección de Hidrografı́a y Navegación(DHN), Peruvian NavyLaboratory FENLAB, Universidad Nacional
Mayor de San Marcos (UNMSM)
Participant of Tsunami Group and Intergovemmental Oceanographic Commision Tsunami
Programme (IOC ICG/PTWS)

Address:
Calle Roca N 116, Chucuito-Callao, Peru

Brief Career:
2012-2014 M.Sc. Student, Faculty of Physical Sciences, National
University of San Marcos, Peru
2010 B.Sc. Physics, Faculty of Physical Sciences, National University of
San Marcos, Peru

Selected Publications:• C. Jiménez and N. Moggiano, “Seismic and source of Pisco 2007
Earthquake,” Science of Tsunami Hazards (to be published).
• N. Moggiano and C. Jimenez, “The 15 August 2007 Peru tsunami:
Numerical modeling, observation and validation,” XIV Giambiagi WinterSchool “Applied and Environmental Geophysics,” held in Universidad de
Buenos Aires (UBA) 2012 (poster), Buenos Aires, Argentina, 2012, Website http://giambiagi.df.uba.ar/272

Address:
Aoba 6-6-3, Sendai 980-8579, Japan

Brief Career:
2009-2010 Master of Disaster Management, National Graduate Institutefor Policy Studies, Japan
2010-2012 Adjunct Professor, Faculty of Civil Engineering, National
University of Engineering, Peru
2012-present Research Student, ReGiD, IRIDeS, Tohoku University, Japan

Selected Publications:• B. Adriano, E. Mas, S. Koshimura, and Y. Fujii, “Remote Sensing-based
Assessment of Tsunami Vulnerability in the coastal area of Lima, Peru,” inThe 10th International Workshop on Remote Sensing for Disaster
Management, Japan, 2012.
• E. Mas, B. Adriano, S. Koshimura, F. Imamura, J. Kuroiwa, F.
Yamazaki, C. Zavala, and M. Estrada, “Evaluation of Tsunami Evacuation
Building Demand through the Multi-Agent System Simulation of
Residents’ Behavior,” in Proceedings of International Sessions in Coastal
Engineering, JSCE, Vol.3, 2012.
• B. Adriano, S. Koshimura, and Y. Fujii, “Tsunami Source and Inundation
Modeling of the June 2001 Peru Earthquake,” in Joint Conference
Proceedings 9CUEE/4ACEE, Tokyo, pp. 2061-2065, 2012.

Academic Societies & Scientific Organizations:• Peruvian Engineering College• Japan Geoscience Union (JpGU)

Journal of Disaster Research Vol.8 No.2, 2013

Seismic Source of 1746 Callao Earthquake from Tsunami Numerical Modeling

Name:

Name:

Shunichi Koshimura

Hideaki Yanagisawa

Affiliation:

Affiliation:

Professor, Laboratory of Remote Sensing
and Geoinformatics for Disaster Management(ReGiD), International Research Institute ofDisaster Science (IRIDeS), Tohoku University

Department of Regional Management, Facultyof Liberal Arts, Tohoku Gakuin University

Address:

Address:

Aoba 6-6-3, Sendai 980-8579, Japan

2-1-1 Tenjinzawa, Izumi-ku, Sendai, Miyagi 981-3193, Japan

Brief Career:

Brief Career:

2000-2002 JSPS Research Fellow, National Oceanic and Atmospheric
Administration
2002-2005 Research Scientist, Disaster Reduction and Human Renovation
Institute
2005-2012 Associate Professor, Tohoku University
2012- Professor, ReGiD, IRIDeS, Tohoku University, Japan

2008-2009 Post-Doctoral Research Fellow, Graduate School ofEngineering, Tohoku University
2009- Company Member, Tokyo Electric Power Services Company
Limited2012- Lecturer, Department of Regional Management, Faculty of Liberal
Arts, Tohoku Gakuin University

Selected Publications:

Selected Publications:• S. Koshimura, T. Oie, H. Yanagisawa, and F. Imamura, “Developing
fragility functions for tsunami damage estimation using numerical modeland post-tsunami data from Banda Aceh, Indonesia,” Coastal Engineering
Journal, No.3, pp. 243-273, 2009.
• S. Koshimura, Y. Hayashi, K. Munemoto, and F. Imamura, “Effect of theEmperor seamounts on trans-oceanic propagation of the 2006 Kuril Island
earthquake tsunami,” Geophysical Research letters, Vol.35, L02611,
doi:10.1029/2007GL032129, 24, 2008.
• S. Koshimura, T. Katada, H. O.Mofjeld, and Y. Kawata, “A method for
estimating casualties due to the tsunami inundation flow,” Natural Hazard
Vol.39, pp. 265-274, 2006.

Academic Societies & Scientific Organizations:• Japan Society of Civil Engineers (JSCE)• Institute of Social Safety Science• Japan Association for Earthquake Engineering (JAEE)• Japan Society for Computational Engineering and Science (JSCES)• American Geophysical Union (AGU)• H. Yanagisawa, S. Koshimura, K. Goto, T. Miyagi, F. Imamura, A.
Ruangrassamee, and C. Tanavud, “Damage of mangrove forest by the2004 Indian Ocean tsunami at Pakarang Cape and Namkem, Thailand,”
Polish Journal of Environmental Studies, Vol.18, No.1, pp. 35-42, 2009.
• H. Yanagisawa, S. Koshimura, K. Goto, T. Miyagi, F. Imamura, A.
Ruangrassamee, and C. Tanavud, “The reduction effects of mangrove
forest on a tsunami based on field surveys at Pakarang Cape, Thailand andnumerical analysis, Estuarine,” Coastal and Shelf Science, Vol.81, pp.
27-37, 2009.

Academic Societies & Scientific Organizations:• Japan Society of Civil Engineers (JSCE)• Japan Society for Mangroves

Name:
Yushiro Fujii

Affiliation:
Senior Research Scientist, International Institute of Seismology and Earthquake Engineering,
Building Research Institute

Address:
1 Tachihara, Tsukuba, Ibaraki 305-0802, Japan

Brief Career:
2004 Ph.D. Earth and Planetary Sciences, Faculty of Sciences, Kyushu
University, Japan
2005 Researcher, Research Institute for Information Technology, Kyushu
University, Japan
2006 Researcher, AFRC, Advanced Industrial Science and Technology
2006- Research Scientist, IISEE, Building Research Institute, Japan

Selected Publications:• Y. Fujii and K. Satake, “Slip Distribution and Seismic Moment of the2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms
and Coastal Geodetic Data,” Pure and Applied Geophysics, 2012.
• Y. Fujii, K. Satake, S. Sakai, M. Shinohara, and T. Kanazawa, “Tsunami
source of the 2011 off the Pacific coast of Tohoku Earthquake,” Earth,Planets and Space, Vol.63, pp. 815-20, 2011.
• Y. Fujii and K. Satake, “Tsunami Source of the 2004 Sumatra-Andaman
Earthquake Inferred from Tide Gauge and Satellite Data,” Bulletin of the
Seismological Society of America, Vol.97, pp. S192-S207, 2007.

Academic Societies & Scientific Organizations:• Seismological Society of Japan (SSJ)• Japan Geoscience Union (JpGU)• American Geophysical Union (AGU)

Journal of Disaster Research Vol.8 No.2, 2013273