Journal
BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA
Volume 107, Issue 5, Pages 2183-2193Publisher
SEISMOLOGICAL SOC AMER
DOI: 10.1785/0120170077
Keywords
-
Categories
Funding
- Southern California Earthquake Center (SCEC) [7911]
- National Science Foundation (NSF) [ACI-1450451]
- NSF Cooperative Agreement [EAR-1033462]
- U.S. Geological Survey (USGS) [G12AC20038]
- Direct For Computer & Info Scie & Enginr
- Office of Advanced Cyberinfrastructure (OAC) [1450451] Funding Source: National Science Foundation
Ask authors/readers for more resources
In a realistic geological structure with a large contrast in seismic wavespeed between shallow and deep regions, simulation of seismic wave propagation using a spatially uniform grid can be computationally very demanding, due to over-discretization of the high-speed material. Thus, numerical methods that allow for coarser discretization of the faster regions have the potential to be much more efficient. Discontinuous mesh (DM) methods, operating by exchanging wavefield information between media partitions discretized with two different grid spacings, provide a convenient way to improve such efficiency issues. Unfortunately, discontinuous staggered-grid finite-difference (FD) methods typically suffer from inherent stability problems, in particular in strongly heterogeneous media, arising from numerical noise generated at the overlap of the two regions with different grid spacing. We have developed a 3D fourth-order velocity-stress staggered-grid FD DM anelastic wave propagation method (AWP-DM) for seismic wavefield estimation using a discontinuous mesh interface (WEDMI) between fine and coarse meshes. Benchmarks in models with realistic 3D velocity variations and finite-fault sources across the grid interface show stable results for a number of timesteps, exceeding the need dictated by current high-frequency ground-motion simulations. In the case of a factor-of-three ratio between the coarse and fine grid sizes, this method is capable of producing a level of accuracy comparable to that from the uniform fine-grid scheme, using at least 7-8 grid points per minimum S wavelength inside the mesh overlap zone.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available