Search NASA⌕ Search

Engineering topics

Zhang, Haijiang

Publications and source records attributed to Zhang, Haijiang.

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗