Numerical analysis of impacts of 2011 Japan Tohoku tsunami on China Coast

Ren, Zhiyuan; Wang, Ben-long; Fan, Tingting; Liu, Hua

doi:10.1016/s1001-6058(11)60399-6

Cited by 31 publications

(14 citation statements)

References 13 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Again, we use the linear high‐order Boussinesq model as reference, but this time we also include computations based on the nondispersive linear shallow water model [ Ren et al ., ] in order to illustrate the influence of dispersion. Both models are set up to cover the first quadrant of the problem using 2000 by 2000 grid points with grid size dx = dy = 1.0 km and 3000 time steps with step size dt = 5.0 s. A perspective snapshot of the computed surface elevation is shown in Figure .…”

Section: Applications On a Constant Depthmentioning

confidence: 99%

“…Weak dispersion can typically be handled by any standard Boussinesq model and these are now being used for tsunami modeling on a regular basis [see e.g., Løvholt et al ., ; Glimsdal et al ., ; Grilli et al ., ]. It is, however, still common to use models based on the nondispersive nonlinear shallow water (NSW) equations [see e.g., Tang et al ., ; Ren et al ., ]. In principle, the NSW equations cannot correctly predict the temporal evolution of the transient waves, but it has to be admitted that these models have been relatively successful in simulating geophysical tsunamis.…”

Section: Introductionmentioning

confidence: 99%

“…This includes an initial 1‐D rectangular disturbance, an initial 2‐D square disturbance, and an initial 2‐D Gaussian disturbance. The convolution solutions are compared with numerical results from models solving linear high‐order Boussinesq equations [ Madsen et al ., ] as well as linear shallow water equations [ Ren et al ., ].…”

Section: Introductionmentioning

confidence: 99%

“…Section 7 covers the extension to uneven bottom and a special convolution procedure suited for this purpose. As a test example, we focus on the 2011 Tohoku tsunami, and compare the convolution results with DART measurements and with numerical results obtained from solving the nonlinear shallow water equations [ Ren et al ., ].…”

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

Uniform asymptotic approximations for transient waves due to an initial disturbance

et al. 2016

View full text Add to dashboard Cite

In this work, we first present a semianalytical method for the evolution of linear fully dispersive transient waves generated by an initial surface displacement and propagating over a constant depth. The procedure starts from Fourier and Hankel transforms and involves a combination of the method of stationary phase, the method of uniform asymptotic approximations and various Airy integral formulations. Second, we develop efficient convolution techniques expressed as single and double summations over the source area. These formulations are flexible, extremely fast, and highly accurate even for the dispersive tail of the transient waves. To verify the accuracy of the embedded dispersion properties, we consider test cases with sharp‐edged disturbances in 1‐D and 2‐D. Furthermore, we consider the case of a relatively blunt Gaussian disturbance in 2‐D. In all cases, the agreement between the convolution results and simulations with a high‐order Boussinesq model is outstanding. Finally, we make an attempt to extend the convolution methods to geophysical tsunami problems taking into account, e.g., uneven bottom effects. Unfortunately, refraction/diffraction effects cannot easily be incorporated, so instead we focus on the incorporation of linear shoaling and its effect on travel time and temporal evolution of the surface elevation. The procedure is tested on data from the 2011 Japan tsunami. Convolution results are likewise compared to model simulations based on the nonlinear shallow water equations and both are compared with field observations from 10 deep water DART buoys. The near‐field results are generally satisfactory, while the far‐field results leave much to be desired.

show abstract

Section: Applications On a Constant Depthmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Uniform asymptotic approximations for transient waves due to an initial disturbance

et al. 2016

View full text Add to dashboard Cite

show abstract

“…In this paper we present a GPU‐accelerated version of the open source GeoClaw code distributed as part of Clawpack (Clawpack Development Team, ), which is widely used for modeling tsunamis and has undergone a benchmarking process in 2011 (Horrillo et al, ) as part of the National Tsunami Hazard Mitigation Program, allowing its use on tsunami hazard modeling projects supported by this program (e.g., González et al, ; Titov et al, ). It has been used for many other projects involving tsunami hazard assessment and/or the study of historic or paleotsunamis, (e.g., Amir et al, ; Arcos & LeVeque, ; Cienfuegos et al, ; Galanti et al, ; Hayes & Furlong, ; Johnstone & Lence, ; MacInnes et al, ; Ren et al, ). GeoClaw has also been used in modeling storm surge (e.g., Mandli & Dawson, ) and dam failures (e.g., George, ; Turzewski et al, ), and so the GPU‐accelerated version should prove useful beyond tsunami modeling.…”

Section: Introductionmentioning

confidence: 99%

Accelerating an Adaptive Mesh Refinement Code for Depth‐Averaged Flows Using GPUs

Qin

LeVeque

Motley

2019

J Adv Model Earth Syst

View full text Add to dashboard Cite

Solving the shallow water equations efficiently is critical to the study of natural hazards induced by tsunami and storm surge, since it provides more response time in an early warning system and allows more runs to be done for probabilistic assessment where thousands of runs may be required. Using adaptive mesh refinement speeds up the process by greatly reducing computational demands while accelerating the code using the graphics processing unit (GPU) does so through using faster hardware. Combining both, we present an efficient CUDA implementation of GeoClaw, an open source Godunov‐type high‐resolution finite volume numerical scheme on adaptive grids for shallow water system with varying topography. The use of adaptive mesh refinement and spherical coordinates allows modeling transoceanic tsunami simulation. Numerical experiments on the 2011 Japan tsunami and a local tsunami triggered by a hypothetical Mw 7.3 earthquake on the Seattle Fault illustrate the correctness and efficiency of the code, which implements a simplified dimensionally split version of the algorithms. Both numerical simulations are conducted on subregions on a sphere with adaptive grids that adequately resolve the propagating waves. The implementation is shown to be accurate and faster than the original when using Central Processing Units (CPUs) alone. The GPU implementation, when running on a single GPU, is observed to be 3.6 to 6.4 times faster than the original model running in parallel on a 16‐core CPU. Three metrics are proposed to evaluate relative performance of the model, which shows efficient usage of hardware resources.

show abstract

A Self‐Consistent Fault Slip Model for the 2011 Tohoku Earthquake and Tsunami

2018

View full text Add to dashboard Cite

The unprecedented geophysical and hydrographic data sets from the 2011 Tohoku earthquake and tsunami have facilitated numerous modeling and inversion analyses for a wide range of dislocation models. Significant uncertainties remain in the slip distribution as well as the possible contribution of tsunami excitation from submarine slumping or anelastic wedge deformation. We seek a self‐consistent model for the primary teleseismic and tsunami observations through an iterative approach that begins with downsampling of a finite fault model inverted from global seismic records. Direct adjustment of the fault displacement guided by high‐resolution forward modeling of near‐field tsunami waveform and runup measurements improves the features that are not satisfactorily accounted for by the seismic wave inversion. The results show acute sensitivity of the runup to impulsive tsunami waves generated by near‐trench slip. The adjusted finite fault model is able to reproduce the DART records across the Pacific Ocean in forward modeling of the far‐field tsunami as well as the global seismic records through a finer‐scale subfault moment‐ and rake‐constrained inversion, thereby validating its ability to account for the tsunami and teleseismic observations without requiring an exotic source. The upsampled final model gives reasonably good fits to onshore and offshore geodetic observations albeit early after‐slip effects and wedge faulting that cannot be reliably accounted for. The large predicted slip of over 20 m at shallow depth extending northward to ~39.7°N indicates extensive rerupture and reduced seismic hazard of the 1896 tsunami earthquake zone, as inferred to varying extents by several recent joint and tsunami‐only inversions.

show abstract

Numerical analysis of impacts of 2011 Japan Tohoku tsunami on China Coast

Cited by 31 publications

References 13 publications

Uniform asymptotic approximations for transient waves due to an initial disturbance

Uniform asymptotic approximations for transient waves due to an initial disturbance

Accelerating an Adaptive Mesh Refinement Code for Depth‐Averaged Flows Using GPUs

A Self‐Consistent Fault Slip Model for the 2011 Tohoku Earthquake and Tsunami

Contact Info

Product

Resources

About