海底地震有限断层破裂模型对近场海啸数值预报的影响

快速准确的海啸源模型是近场海啸精确预警的关键.尽管目前还没有办法直接对其进行正演定量计算,但是可以通过多源地震、海啸观测数据进行反演或联合反演推算.不同的海啸源可能导致不同的预警结论,了解不同类型海啸源适用性、评估海啸源特征差异对近场海啸的影响,无论对于海啸预警还是海啸模拟研究尤为重要.本文评估分析了6种不同同震断层模型对2011年3月11日日本东北地震海啸近场数值预报的影响,重点对比分析了有限断层模型与均一滑动场模型对近场海啸产生、传播、淹没特征的影响及各自的误差.研究表明:近场海啸波能量分布主要取决于海啸源分布特征,特别是走向角的差异对海啸能量分布影响较大;有限断层模型对海啸灾害最为严重的39°N以南沿岸地区的最大海啸爬坡高度明显优于均一滑动场模型结果;综合对比DART浮标、GPS浮标及近岸潮位站共32个站次的海啸波幅序列结果发现有限断层模型整体平均绝对/相对误差比均一滑动场模型平均误差要低,其中Fujii海啸源的平均绝对/相对误差最小,分别是0.56m和26.71%.UCSB海啸源的平均绝对/相对误差次之.3个均一滑动场模型中USGSCMT海啸源模拟精度最高.相对于深海、浅海观测站,有限断层模型比均一滑动场模型对近岸观测站计算精度更高.海啸源误差具有显著的方向性,可能与反演所采用的波形数据的代表性有关;谱分析结果表明Fujii海啸源对在12至60min主频波谱的模拟要优于UCSB海啸源.海啸源中很难真实反映海底地震破裂过程,然而通过联合反演海啸波形数据推算海啸源的方法可以快速确定海啸源,并且最大限度的降低地震破裂过程与海啸产生的不确定性带来的误差.     

Modeling of the 2011 Japan Tsunami: Lessons for Near-Field Forecast

During the devastating 11 March 2011 Japanese tsunami, data from two tsunami detectors were used to determine the tsunami source within 1.5 h of earthquake origin time. For the first time, multiple near-field tsunami measurements of the 2011 Japanese tsunami were used to demonstrate the accuracy of the National Oceanic and Atmospheric Administration (NOAA) real-time flooding forecast system in the far field. To test the accuracy of the same forecast system in the near field, a total of 11 numerical models with grids telescoped to 2 arcsec (~60 m) were developed to hindcast the propagation and coastal inundation of the 2011 Japanese tsunami along the entire east coastline of Japan. Using the NOAA tsunami source computed in near real-time, the model results of tsunami propagation are validated with tsunami time series measured at different water depths offshore and near shore along Japan’s coastline. The computed tsunami runup height and spatial distribution are highly consistent with post-tsunami survey data collected along the Japanese coastline. The computed inundation penetration also agrees well with survey data, giving a modeling accuracy of 85.5 % for the inundation areas along 800 km of coastline between Ibaraki Prefecture (north of Kashima) and Aomori Prefecture (south of Rokkasho). The inundation model results highlighted the variability of tsunami impact in response to different offshore bathymetry and flooded terrain. Comparison of tsunami sources inferred from different indirect methods shows the crucial importance of deep-ocean tsunami measurements for real-time tsunami forecasts. The agreement between model results and observations along Japan’s coastline demonstrate the ability and potential of NOAA’s methodology for real-time near-field tsunami flooding forecasts. An accurate tsunami flooding forecast within 30 min may now be possible using the NOAA forecast methodology with carefully placed tsunameters and large-scale high-resolution inundation models with powerful computing capabilities.     

Combined Effects of Tectonic and Landslide-Generated Tsunami Runup at Seward, Alaska During the M W 9.2 1964 Earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the $M_{\rm W}$ 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5?min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211?million m³ (Haeussler et?al. in Submarine mass movements and their consequences, pp 269?C278, 2007). The first tectonic tsunami wave arrived in Resurrection Bay about 30?min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et?al. in Pure Appl Geophys 166:131?C152, 2009). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559?C572, 1994) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.     

Field Survey and Geological Effects of the 15 November 2006 Kuril Tsunami in the Middle Kuril Islands

The near-field expression of the tsunami produced by the 15 November 2006 Kuril earthquake (M_w 8.1–8.4) in the middle Kuril Islands, Russia, including runup of up to 20 m, remained unknown until we conducted a post-tsunami survey in the summer of 2007. Because the earthquake occurred between summer field expeditions in 2006 and 2007, we have observations, topographic profiles, and photographs from three months before and nine months after the tsunami. We thoroughly surveyed portions of the islands of Simushir and Matua, and also did surveys on parts of Ketoi, Yankicha, Ryponkicha, and Rasshua. Tsunami runup in the near-field of the middle Kuril Islands, over a distance of about 200 km, averaged 10 m over 130 locations surveyed and was typically between 5 and 15 m. Local topography strongly affected inundation and somewhat affected runup. Higher runup generally occurred along steep, protruding headlands, whereas longer inundation distances occurred on lower, flatter coastal plains. Sediment transport was ubiquitous where sediment was available—deposit grain size was typically sand, but ranged from mud to large boulders. Wherever there were sandy beaches, a more or less continuous sand sheet was present on the coastal plain. Erosion was extensive, often more extensive than deposition in both space and volume, especially in areas with runup of more than 10 m. The tsunami eroded the beach landward, stripped vegetation, created scours and trim lines, cut through ridges, and plucked rocks out of the coastal plain.     

Evaluation of the 15 July 2009 Fiordland,New Zealand Tsunami in the Source Region

We undertake detailed near-field numerical modelling of the tsunami generated by the 15 July 2009 earthquake (Mw 7.8) in Fiordland, New Zealand. High resolution bathymetry and topography data at Breaksea and Dusky Sounds, and Chalky and Preservation Inlets are derived mostly from digitised New Zealand nautical charts, Shuttle Radar Topographic Mission (SRTM) 3 arc-second data, and General Bathymetric Chart of the Ocean (GEBCO) 30 s data. A combination of continuous and campaign Global Positioning System (GPS), satellite radar (ALOS/PALSAR InSAR images) and seismology data are used to constrain the seafloor deformation for the initial tsunami condition. This source model, derived independently of DART observations, provides an excellent fit to observed tsunami elevations recorded by DART buoy 55015. The model results in the near field show maximum tsunami elevations in the range 0.5–2.0 m inside the sounds and inlets with maximum flow speeds of 3.0 m/s. Along the open coast, maximum tsunami elevations reach 2.0 m. The high flow speeds through the inlets may change the inlet stratifications and water mass inside the sounds. Media reports and field reconnaissance data show some tsunami evidence at Cormorant Cove, Duck and Goose Coves, and Passage Point.     

Tsunami Modeling to Validate Slip Models of the 2007 M w 8.0 Pisco Earthquake, Central Peru

Following the 2007, August 15th, M _w 8.0, Pisco earthquake in central Peru, Sladen et al. (J Geophys Res 115: B02405, 2010) have derived several slip models of this event. They inverted teleseismic data together with geodetic (InSAR) measurements to look for the co-seismic slip distribution on the fault plane, considering those data sets separately or jointly. But how close to the real slip distribution are those inverted slip models? To answer this crucial question, the authors generated some tsunami records based on their slip models and compared them to DART buoys, tsunami records, and available runup data. Such an approach requires a robust and accurate tsunami model (non-linear, dispersive, accurate bathymetry and topography, etc.) otherwise the differences between the data and the model may be attributed to the slip models themselves, though they arise from an incomplete tsunami simulation. The accuracy of a numerical tsunami simulation strongly depends, among others, on two important constraints: (i) A fine computational grid (and thus the bathymetry and topography data sets used) which is not always available, unfortunately, and (ii) a realistic tsunami propagation model including dispersion. Here, we extend Sladen’s work using newly available data, namely a tide gauge record at Callao (Lima harbor) and the Chilean DART buoy record, while considering a complete set of runup data along with a more realistic tsunami numerical that accounts for dispersion, and also considering a fine-resolution computational grid, which is essential. Through these accurate numerical simulations we infer that the InSAR-based model is in better agreement with the tsunami data, studying the case of the Pisco earthquake indicating that geodetic data seems essential to recover the final co-seismic slip distribution on the rupture plane. Slip models based on teleseismic data are unable to describe the observed tsunami, suggesting that a significant amount of co-seismic slip may have been aseismic. Finally, we compute the runup distribution along the central part of the Peruvian coast to better understand the wave amplification/attenuation processes of the tsunami generated by the Pisco earthquake.     

Reconstruction of Tsunami Inland Propagation on December 26, 2004 in Banda Aceh, Indonesia, through Field Investigations

This paper presents the results from an extensive field data collection effort following the December 26, 2004 earthquake and tsunami in Banda Aceh, Sumatra. The data were collected under the auspices of TSUNARISQUE, a joint French-Indonesian program dedicated to tsunami research and hazard mitigation, which has been active since before the 2004 event. In total, data from three months of field investigations are presented, which detail important aspects of the tsunami inundation dynamics in Banda Aceh. These include measurements of runup, tsunami wave heights, flow depths, flow directions, event chronology and building damage patterns. The result is a series of detailed inundation maps of the northern and western coasts of Sumatra including Banda Aceh and Lhok Nga. Among the more important findings, we obtained consistent accounts that approximately ten separate waves affected the region after the earthquake; this indicates a high-frequency component of the tsunami wave energy in the extreme near-field. The largest tsunami wave heights were on the order of 35 m with a maximum runup height of 51 m. This value is the highest runup value measured in human history for a seismically generated tsunami. In addition, our field investigations show a significant discontinuity in the tsunami wave heights and flow depths along a line approximately 3 km inland, which the authors interpret to be the location of the collapse of the main tsunami bore caused by sudden energy dissipation. The propagating bore looked like a breaking wave from the landward side although it has distinct characteristics. Patterns of building damage are related to the location of the propagating bore with overall less damage to buildings beyond the line where the bore collapsed. This data set was built to be of use to the tsunami community for the purposes of calibrating and improving existing tsunami inundation models, especially in the analysis of extreme near-field events.     

Modeling the seismic source and tsunami generation of the December 12, 1992 Flores Island,Indonesia, earthquake

On December 12, 1992 a large earthquake (M _s 7.5) occurred just north of Flores Island, Indonesia which, along with the tsunami it generated, killed more than 2,000 people. In this study, teleseismicP andSH waves, as well asPP waves from distances up to 123°, are inverted for the orientations and time histories of multiple point sources. By repeating the inversion for reasonable values of depth, time separation and spatial separation, a 2-fault model is developed. Next, the vertical deformation of the seafloor is estimated from this fault model. Using a detailed bathymetric model, linear and nonlinear tsunami propagation models are tested. The data consist of a single tide gauge record at Palopo (650 km to the north), as well as tsunami runup height measurements from Flores Island and nearby islands. Assuming a tsunami runup amplification factor of two, the two-fault model explains the tide gauge record and the tsunami runup heights on most of Flores Island. It cannot, however, explain the large tsunami runup heights observed near Leworahang (on Hading Bay) and Riangkroko (on the northeast peninsula). Massive coastal slumping was observed at both of these locations. A final model, which in addition to the two faults, includes point sources of large vertical displacement at these two locations explains the observations quite well.     

Tsunami Source of the 2010 Mentawai, Indonesia Earthquake Inferred from Tsunami Field Survey and Waveform Modeling

The 2010 Mentawai earthquake (magnitude 7.7) generated a destructive tsunami that caused more than 500 casualties in the Mentawai Islands, west of Sumatra, Indonesia. Seismological analyses indicate that this earthquake was an unusual “tsunami earthquake,” which produces much larger tsunamis than expected from the seismic magnitude. We carried out a field survey to measure tsunami heights and inundation distances, an inversion of tsunami waveforms to estimate the slip distribution on the fault, and inundation modeling to compare the measured and simulated tsunami heights. The measured tsunami heights at eight locations on the west coasts of North and South Pagai Island ranged from 2.5 to 9.3 m, but were mostly in the 4–7 m range. At three villages, the tsunami inundation extended more than 300 m. Interviews of local residents indicated that the earthquake ground shaking was less intense than during previous large earthquakes and did not cause any damage. Inversion of tsunami waveforms recorded at nine coastal tide gauges, a nearby GPS buoy, and a DART station indicated a large slip (maximum 6.1 m) on a shallower part of the fault near the trench axis, a distribution similar to other tsunami earthquakes. The total seismic moment estimated from tsunami waveform inversion was 1.0 × 10²¹ Nm, which corresponded to M_w 7.9. Computed coastal tsunami heights from this tsunami source model using linear equations are similar to the measured tsunami heights. The inundation heights computed by using detailed bathymetry and topography data and nonlinear equations including inundation were smaller than the measured ones. This may have been partly due to the limited resolution and accuracy of publically available bathymetry and topography data. One-dimensional run-up computations using our surveyed topography profiles showed that the computed heights were roughly similar to the measured ones.     

Numerical Modelling of the 26th December 2004 Indian Ocean Tsunami for the Southeastern Coast of India

A numerical simulation of the 26th December, 2004 Indian Ocean tsunami of the Tamil Nadu coastal zone is presented. The simulation approach is based on a fully nonlinear Boussinesq tsunami propagation model and included an accurate computational domain and a robust coseismic source. The simulation is first confronted to available tide gauge and runup observations. The agreement between observations and the predicted wave heights allowed a reasonable validation of the simulation. As a result, a full picture of the tsunami impact is provided over the entire coastal zone Tamil Nadu. The processes responsible for coastal vulnerability are discussed.     

A macroseismic model of a tsunami source and estimation of its efficiency by numerical modeling

A typical model of the source of a tsunami (“macroseismic source”) is suggested for use in approximate estimation of maximum tsunami height using straightforward numerical modeling. In this paper the model is tested using three actual events: the 1952 North Kuril Is., 1971 Moneron, and 1994 Shikotan earthquakes, which excited considerable tsunamis at Russia’s Far East coasts. Comparison of the maximum tsunami runup values as obtained in numerical experiments with observations of actual tsunamis showed that the numerical model proposed here is suitable for crude estimation of tsunami runup and tsunami waiting times for coastal population centers in the near zone of a tsunami source.     

Evaluation of the use of paleotsunami deposits to reconstruct inundation distance and runup heights associated with prehistoric inundation events,Crescent City,southern Cascadia margin

Historic‐ and prehistoric‐tsunami sand deposits are used to independently establish runup records for tsunami hazard mitigation and modeled runup verification in Crescent City, California, located in the southern Cascadia Subduction Zone. Inundation from historic (1964) farfield tsunami (~5–6 m runup height) left sand sheet deposits (100–200 m width) in wetlands located behind a low beach ridge [3–4 m elevation of the National Geodetic Vertical Datum of 1988 (NAVD88)]. The most landward flooding lines (4·5–5 m elevation) in high‐gradient alluvial wetlands exceed the 1964 sand sheet records of inundation by 1–2 m in elevation. The most landward flooding in low‐gradient alluvial wetlands exceed the corresponding sand sheet record of inundation distance by 1000 m. Nevertheless, the sand sheet record is an important proxy for high‐velocity inundation. Sand sheet deposition from the 1964 historic tsunami closely corresponds to the landward extent of large debris transport and structural damage in the Crescent City waterfront. The sand sheet deposits provide a proxy for maximum hazard or ‘kill zone’ in the study area. Six paleotsunami sand sheets (0·3–3 ka) are recorded in the back‐ridge marshes in Crescent City, yielding a ~450 year mean recurrence interval for nearfield Cascadia tsunami. Two paleotsunami sand deposit records, likely correlated to Cascadia ruptures between 1·0 and 1·5 ka, are traced to 1·2 km distance and 9–10 m elevation, as adjusted for paleo‐sea level. The paleotsunami sand deposits demonstrate at least twice the runup height, and four times the inundation distance of the farfield 1964 tsunami sand sheet in the same marsh system. The preserved paleotsunami deposits in Crescent City are compared to the most landward flooding, as modeled by other investigators from a predicted Cascadia (~ Mw 9) rupture. The short geologic record (~1·5 ka) yields slightly lower runup records than those predicted for the modeled Mw 9 rupture scenario in the same marsh, but it generally verifies predicted maximum tsunami runup for use in the planning of emergency response and rapid evacuation. Copyright © 2011 John Wiley & Sons, Ltd.     

2004年苏门答腊地震的几个断层滑动模型的全球同震位移对比

2004年苏门答腊大地震后,不同作者根据地震波和/或GPS观测,提出了不同的断层错动模型.在利用同震位移观测资料反演断层滑动模型时,由于使用半无限空间均匀介质模型或半无限空间分层介质模型,一般只能利用近场位移GPS观测约束,无法利用远场资料,这些模型有时差别颇大,如何区别这些模型的优劣是一个仍尚未解决的问题.本文采用等效体力有限元方法,在考虑地球球形和分层的条件下,对四个不同作者提供的2004年苏门答腊地震的断层滑动模型计算全球同震位移.由于采用了球形模型,所以不仅可以把四个模型的近场位移计算结果与GPS数据进行对比,而且可以把远场位移计算结果与GPS数据进行对比.我们发现,垂直位移对断层滑动模型的依赖性小于水平位移.四个模型计算的近场位移与GPS位移符合程度均较好,但是四个模型计算的远场位移与GPS位移符合情况有很大不同,其中Chlieh等（2007）模型在近场与远场符合程度均很好,是四个模型中最好的.另外还探讨了断层反演数据资料、断层几何模型以及地球模型对计算结果的影响.对于特大地震,全球同震位移观测与计算值吻合程度的好坏是衡量断层滑动模型的合理性的一个重要依据.     

Validation and Verification of a Numerical Model for Tsunami Propagation and Runup

A robust numerical model to simulate propagation and runup of tsunami waves in the framework of non-linear shallow water theory is developed. The numerical code adopts a staggered leapfrog finite-difference scheme to solve the shallow water equations formulated for depth-averaged water fluxes in spherical coordinates. A temporal position of the shoreline is calculated using a free-surface moving boundary algorithm. For large scale problems, the developed algorithm is efficiently parallelized employing a domain decomposition technique. The developed numerical model is benchmarked in an exhaustive series of tests suggested by NOAA. We conducted analytical and laboratory benchmarking for the cases of solitary wave runup on simple beaches, runup of a solitary wave on a conically-shaped island, and the runup in the Monai Valley, Okushiri Island, Japan, during the 1993 Hokkaido-Nansei-Oki tsunami. In all conducted tests the calculated numerical solution is within an accuracy recommended by NOAA standards. We summarize results of numerical benchmarking of the model, its strengths and limits with regards to reproduction of fundamental features of coastal inundation, and also illustrate some possible improvements.     

Satellite Data for a Rapid Assessment of Tsunami Inundation Areas after the 2011 Tohoku Tsunami

The M _w = 9.0 earthquake that occurred off the coast of Japan’s Tohoku region produced a great tsunami causing catastrophic damage and loss of life. Within hours of the tsunami event, satellite data were readily available and massive media coverage immediately circulated thousands of photographs and videos of the tsunami. Satellite data allow a rapid assessment of inundated areas where access can be difficult either as a result of damaged infrastructure (e.g., roads, bridges, ports, airports) or because of safety issues (e.g., the hazard at Nuclear Power Plant at Fukushima). In this study, we assessed in a day tsunami inundation distances and runup heights using satellite data (very high-resolution satellite images from the GeoEye1 satellite and from the DigitalGlobe worldview, SRTM and ASTER GDEM) of the Tohoku region, Northeast Japan. Field survey data by Japanese and other international scientists validated our results. This study focused on three different locations. Site selection was based on coastal morphologies and the distance to the tsunami source (epicenter). Study sites are Rikuzentakata, Oyagawahama, and Yagawahama in the Oshika Peninsula, and the Sendai coastal plain (Sendai City to Yamamoto City). Maximum inundation distance (6 km along the river) and maximum runup (39 m) at Rikuzentakata estimated from satellite data agree closely with the 39.7 m inundation reported in the field. Here the ria coastal morphology and horn shaped bay enhanced the tsunami runup and effects. The Sendai coastal plain shows large inundation distances (6 km) and lower runup heights. Natori City and Wakabayashi Ward, on the Sendai plain, have similar runup values (12 and 16 m, respectively) obtained from SRTM data; these are comparable to those obtained from field surveys (12 and 9.5 m). However, at Yagawahama and Oyagawahama, Miyagi Prefecture, both SRTM and ASTER data provided maximum runup heights (41 to 45 m and 33 to 34 m, respectively), which are higher than those measured in the field (about 27 m). This difference in DEM and field data is associated with ASTER and SRTM DEM’s pixel size and vertical accuracy, the latter being dependent on ground coverage, slope, aspect and elevation. Countries with less access to technology and infrastructure can benefit from the use of satellite imagery and freely available DEMs for an initial, pre-field surveys, rapid estimate of inundated areas, distances and runup, and for assisting in hazard management and mitigation after a natural disaster.     

Coastal Sedimentation Associated with the Tohoku Tsunami of 11 March 2011 in South Kuril Islands, NW Pacific Ocean

Sediment deposited by the Tohoku tsunami of March 11, 2011 in the Southern Kurils (Kunashir, Shikotan, Zeleniy, Yuri, Tanfiliev islands) was radically different from sedimentation during local strong storms and from tsunamis with larger runup at the same location. Sediments from the 2011 Tohoku tsunami were surveyed in the field, immediately and 6 months after the event, and analyzed in the laboratory for sediment granulometry, benthos Foraminifa assemblages, and diatom algae. Run-up elevation and inundation distance were calculated from the wrackline (accumulations of driftwood, woody debris, grass, and seaweed) marking the distal edge of tsunami inundation. Run-up of the tsunami was 5 m at maximum, and 3–4 m on average. Maximum distance of inundation was recorded in river mouths (up to 630 m), but was generally in the range of 50–80 m. Although similar to the local strong storms in runup height, the tsunami generally did not erode the coast, nor leave a deposit. However, deposits uncharacteristic of tsunami, described as brown aleuropelitic (silty and clayey) mud rich in organic matter, were found in closed bays facing the South Kuril Strait. These closed bays were covered with sea ice at the time of tsunami. As the tsunami waves broke the ice, the ice floes enhanced the bottom erosion on shoals and destruction of low-lying coastal peatland even at modest ranges of runup. In the muddy tsunami deposits, silt comprised up to 64 % and clay up to 41.5 %. The Foraminifera assemblages displayed features characteristic of benthic microfauna in the near-shore zone. Deep-sea diatoms recovered from tsunami deposits in two closely situated bays, namely Krabovaya and Otradnaya bays, had different requirements for environmental temperature, suggesting these different diatoms were brought to the bays by the tsunami wave entraining various water masses when skirting the island from the north and from the south.     

Tsunami Probability in the Caribbean Region

We calculated tsunami runup probability (in excess of 0.5 m) at coastal sites throughout the Caribbean region. We applied a Poissonian probability model because of the variety of uncorrelated tsunami sources in the region. Coastlines were discretized into 20 km by 20 km cells, and the mean tsunami runup rate was determined for each cell. The remarkable ~500-year empirical record compiled by O’Loughlin and Lander (2003) was used to calculate an empirical tsunami probability map, the first of three constructed for this study. However, it is unclear whether the 500-year record is complete, so we conducted a seismic moment-balance exercise using a finite-element model of the Caribbean-North American plate boundaries and the earthquake catalog, and found that moment could be balanced if the seismic coupling coefficient is c = 0.32. Modeled moment release was therefore used to generate synthetic earthquake sequences to calculate 50 tsunami runup scenarios for 500-year periods. We made a second probability map from numerically-calculated runup rates in each cell. Differences between the first two probability maps based on empirical and numerical-modeled rates suggest that each captured different aspects of tsunami generation; the empirical model may be deficient in primary plate-boundary events, whereas numerical model rates lack backarc fault and landslide sources. We thus prepared a third probability map using Bayesian likelihood functions derived from the empirical and numerical rate models and their attendant uncertainty to weight a range of rates at each 20 km by 20 km coastal cell. Our best-estimate map gives a range of 30-year runup probability from 0–30% regionally.     

2010年智利地震海啸数值模拟及其对我国沿海的影响分析

2010年2月27日06时34分(北京时间14时34分),在智利中南部近岸(36.1°S,72.6°W)发生Mw8.8级地震,并引发了泛太平洋范围的海啸,太平洋沿岸多个国家的验潮站和海啸监测浮标均监测到了强震引发的海啸;海啸波传播25 h后到达我国沿海.本文利用海啸数值模型对此次地震海啸进行了数值模拟.重点模拟了我国沿...     

Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

Field Measurements and Numerical Simulations of the 2004 Tsunami Impact on the East Coast of Sri Lanka

This paper employs a numerical model of tsunami propagation together with documented observations and field measurements of the evidence left behind by the tsunami in December 2004, to identify and interpret the factors that have contributed to the significant spatial variability of the level of tsunami impact along the coastal belt of the eastern province of Sri Lanka. The model results considered in the present analysis include the distribution of the amplitude of the tsunami and the pattern of wave propagation over the continental shelf off the east coast, while the field data examined comprise the maximum water levels measured at or near the shoreline, the horizontal inundation distances and the number of housing and other buildings damaged. The computed maximum amplitude of the tsunami at water points nearest the shoreline along the east coast shows considerable variation ranging from 2.2 m to 11.4 m with a mean value of 5.7 m; moreover, the computed amplitudes agree well with the available field measurements. We also show that the shelf bathymetry off the east coast, particularly the submarine canyons at several locations, significantly influences the near-shore transformation of tsunami waves, and consequently, the spatial variation of the maximum water levels along the coastline. The measured values of inundation also show significant variation along the east coast and range from 70 m to 4560 m with a median value of 700 m. Our analyses of field data also show the dominant influence of the coastal topography and geomorphology on the extent of tsunami inundation. Furthermore, the measured inundation distances indicate no apparent correlation with the computed tsunami heights at the respective locations. We also show that both the computed tsunami heights and the measured inundation distances for the east coast closely follow the log-normal statistical distribution.