The Community Velocity Model V.1.0 of Southwest China, Constructed from Joint Body- and Surface-Wave Travel-Time Tomography

Author(s):  
Ying Liu ◽  
Huajian Yao ◽  
Haijiang Zhang ◽  
Hongjian Fang

Abstract Southwest China, located at the southeastern margin of the Tibetan plateau, plays an important role for the plateau growth and its material extrusion. It has complicated tectonic environment and strong seismic activities including the 2008 Wenchuan great earthquake. Numerous geophysical studies have been conducted in southwest China. However, a community velocity model (CVM) in this region is still not available, which makes it difficult to have a consistent catalog of earthquake locations and focal mechanisms and a consistent velocity model for simulating strong ground motions and evaluating earthquake hazards. In this study, we aim at building a high-resolution CVM (both VP and VS) of the crust and uppermost mantle in southwest China along with earthquake locations by joint inversion of body- and surface-wave travel-time data. In total, we have assembled 386,958 P- and 372,662 S-wave first arrival times and nearly 8100 Rayleigh-wave dispersion curves in the period band of 5–50 s. A multigrid strategy is adopted in the joint inversion. A coarser horizontal grid interval of 0.5° is first used and then a finer grid interval of 0.25° is used with initial models interpolated from the coarser-grid inverted velocity models. The spatial resolution of both VP and VS models can reach up to 0.5° horizontally and 10 km vertically according to the checkerboard tests. The comparisons of our inverted VP and VS models with those from other studies show general consistency in large-scale features. The inverted models are further validated by P-wave arrival times from active sources and Rayleigh-wave data. In general, our velocity models show two low-velocity zones in the middle-lower crust and a prominent high-velocity region in between them. Our new models have been served as the first version of the CVM in southwest China (SWChinaCVM-1.0) for future studies.

Author(s):  
Shoucheng Han ◽  
Haijiang Zhang ◽  
Hailiang Xin ◽  
Weisen Shen ◽  
Huajian Yao

Abstract Xin et al. (2019) presented 3D seismic velocity models (VP and VS) of crust and uppermost mantle of continental China using seismic body-wave travel-time tomography, which are referred to as Unified Seismic Tomography Models for Continental China Lithosphere 1.0 (USTClitho1.0). Compared with previous models of continental China, the VP and VS models of USTClitho1.0 have the highest spatial resolution of 0.5°–1.0° in the horizontal direction and are useful for better understanding the complex tectonics of continental China. Although USTClitho1.0 is implicitly constrained by surface-wave data by using the VS model from surface-wave tomography and the converted VP model as initial models for body-wave travel-time tomography, the predicted surface-wave dispersion curves from USTClitho1.0 do not fit the observed data well. Here, we present updated 3D VP and VS models of the continental China lithosphere (USTClitho2.0) by joint inversion of body-wave arrival times and surface-wave dispersion data. Compared with the previous joint inversion scheme of Zhang et al. (2014), similar to Fang et al. (2016), it is further improved by including the sensitivity of surface-wave dispersion data to VP in the new joint inversion system. As a result, the shallow VP structure is also better imaged. In addition, the new joint inversion scheme considers the large topography variations between the eastern and western parts of China. Thus, USTClitho2.0 better resolves the upper-crustal structure of the Tibetan plateau. Compared with USTClitho1.0, USTClitho2.0 fits both body-wave arrival times and surface-wave dispersion data. Thus, the new velocity models are more accurate and can serve as a better reference model for regional-scale tomography and geodynamic studies in continental China.


Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1062-1065 ◽  
Author(s):  
Thomas Gruber ◽  
Stewart A. Greenhalgh

Rectangular grid velocity models and their derivatives are widely used in geophysical inversion techniques. Specifically, seismic tomographic reconstruction techniques, whether they be based on raypath methods (Bregman et al., 1989; Moser, 1991; Schneider et al., 1992; Cao and Greenhalgh, 1993; Zhou, 1993) or full wave equation methods (Vidale, 1990; Qin and Schuster, 1993; Cao and Greenhalgh, 1994) for calculating synthetic arrival times, involve propagation through a grid model. Likewise, migration of seismic reflection data, using asymptotic ray theory or finite difference/pseudospectral methods (Stolt and Benson, 1986; Zhe and Greenhalgh, 1997) involve assigning traveltimes to upward and downward propagating waves at every grid point in the model. The traveltimes in both cases depend on the grid specification. However, the precision level of such numerical models and their dependence on the model parameters is often unknown. In this paper, we describe a two‐dimensional velocity model and derive an error bound for first‐break times calculated with such a model. The analysis provides clear guidelines for grid specifications.


Geophysics ◽  
1992 ◽  
Vol 57 (1) ◽  
pp. 9-14 ◽  
Author(s):  
Gérard C. Herman

A nonlinear inversion method is presented, especially suited for the determination of global velocity models. In a certain sense, it can be considered as a generalization of methods based on traveltimes of reflections, with the requirement of accurately having to determine traveltimes replaced by the (less stringent and less subjective) requirement of having to define time windows around main reflections (or composite reflections) of interest. It is based on an error norm, related to the phase of the wavefield, which is directly computed from wavefield measurements. Therefore, the cumbersome step of interpreting arrivals and measuring arrival times is avoided. The method is applied to the reconstruction of a depth‐dependent global velocity model from a set of plane‐wave responses and is compared to other methods. Despite the fact that the new error norm only makes use of data having a temporal bandwidth of a few Hz, its behavior is very similar to the behavior of the error norm used in traveltime inversion.


1990 ◽  
Vol 80 (2) ◽  
pp. 395-410 ◽  
Author(s):  
Glenn D. Nelson ◽  
John E. Vidale

Abstract We present a new method for locating earthquakes in a region with arbitrarily complex three-dimensional velocity structure, called QUAKE3D. Our method searches a gridded volume and finds the global minimum travel-time residual location within the volume. Any minimization criterion may be employed. The L1 criterion, which minimizes the sum of the absolute values of travel-time residuals, is especially useful when the station coverage is sparse and is more robust than the L2 criterion (which minimizes the RMS sum) employed by most earthquake location programs. On a UNIX workstation with 8 Mbytes memory, travel-time grids of size 150 by 150 by 50 are reasonably employed, with the actual geographic coverage dependent on the grid spacing. Location precision is finer than the grid spacing. Earthquake recordings at six stations in Bear Valley are located as an example, using various layered and laterally varying velocity models. Locations with QUAKE3D are nearly identical to HYPOINVERSE locations when the same flat-layered velocity model is used. For the examples presented, the computation time per event is approximately 4 times slower than HYPOINVERSE, but the computation time for QUAKE3D is dependent only on the grid size and number of stations, and independent of the velocity model complexity. Using QUAKE3D with a laterally varying velocity model results in locations that are physically more plausible and statistically more precise. Compared to flat-layered solutions, the earthquakes are more closely aligned with the surface fault trace, are more uniform in depth distribution, and the event and station travel-time residuals are much smaller. Hypocentral error bars computed by QUAKE3D are more realistic in that the trade-off of depth versus origin time is implicit in our error estimation, but ignored by HYPOINVERSE.


Author(s):  
Jiayan Tan ◽  
Charles A. Langston ◽  
Sidao Ni

ABSTRACT Ambient noise cross-correlations, used to obtain fundamental-mode Rayleigh-wave group velocity estimates, and teleseismic P-wave receiver functions are jointly modeled to obtain a 3D shear-wave velocity model for the crust and upper mantle of Oklahoma. Broadband data from 82 stations of EarthScope Transportable Array, the U.S. National Seismic Network, and the Oklahoma Geological Survey are used. The period range for surface-wave ambient noise Green’s functions is from 4.5 to 30.5 s constraining shear-wave velocity to a depth of 50 km. We also compute high-frequency receiver functions at these stations from 214 teleseismic earthquakes to constrain individual 1D velocity models inferred from the surface-wave tomography. Receiver functions reveal Ps conversions from the Moho, intracrustal interfaces, and shallow sedimentary basins. Shallow low-velocity zones in the model correlate with the large sedimentary basins of Oklahoma. The velocity model significantly improves the agreement of synthetic and observed seismograms from the 6 November 2011 Mw 5.7 Prague, Oklahoma earthquake suggesting that it can be used to improve earthquake location and moment tensor inversion of local and regional earthquakes.


2017 ◽  
Vol 54 (2) ◽  
pp. 163-172 ◽  
Author(s):  
Shutian Ma ◽  
Pascal Audet

Models of the seismic velocity structure of the crust in the seismically active northern Canadian Cordillera remain poorly constrained, despite their importance in the accurate location and characterization of regional earthquakes. On 29 August 2014, a moderate earthquake with magnitude 5.0, which generated high-quality Rayleigh wave data, occurred in the Northwest Territories, Canada, ∼100 km to the east of the Cordilleran Deformation Front. We carefully selected 23 seismic stations that recorded the Rayleigh waves and divided them into 13 groups according to the azimuth angle between the earthquake and the stations; these groups mostly sample the Cordillera. In each group, we measured Rayleigh wave group velocity dispersion, which we inverted for one-dimensional shear-wave velocity models of the crust. We thus obtained 13 models that consistently show low seismic velocities with respect to reference models, with a slow upper and lower crust surrounding a relatively fast mid crustal layer. The average of the 13 models is consistent with receiver function data in the central portion of the Cordillera. Finally, we compared earthquake locations determined by the Geological Survey of Canada using a simple homogenous crust over a mantle half space with those estimated using the new crustal velocity model, and show that estimates can differ by as much as 10 km.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. U13-U33 ◽  
Author(s):  
Alexey Stovas

A proposed power-gradient velocity model incorporates several well-known velocity models as special cases. The model covers a wide range of possible velocity distributions and has four parameters, giving more flexibility in velocity-model manipulation. For this nonlinear model, the kinematic characteristics — offset-traveltime parameteric equations, traveltime parameters, relative geometric spreading, and propagator phase — are computed. The characteristics are investigated with respect to a parameter responsible for nonlinearity of velocity distribution. The inversion of travel-time parameters was studied in three- and four-parameter frameworks.


2009 ◽  
Vol 47 (6) ◽  
Author(s):  
L. De Luca ◽  
R. De Franco ◽  
G. Biella ◽  
A. Corsi ◽  
R. Tondi

We performed an analysis of refraction data recorded in Italy since 1968 in the frame of the numerous deep seismic sounding and wide-angle reflection/refraction projects. The aims of this study are to construct a parametric database including the recording geometric information relative to each profile, the phase pickings and the results of some kinematic analyses performed on the data, and to define a reference 1D velocity model for the Italian territory from all the available refraction data. As concerns the first goal, for each seismic section we picked the P-wave first-arrival-times, evaluated the uncertainties of the arrival-times pickings and determined from each travel time-offset curve the 1D velocity model. The study was performed on 419 seismic sections. Picking was carried out manually by an algorithm which includes the computation of three picking functions and the picking- error estimation. For each of the travel time-offset curves a 1D velocity model has been calculated. Actually, the 1D velocity-depth functions were estimated in three different ways which assume: a constant velocitygradient model, a varying velocity-gradient model and a layered model. As regards the second objective of this work, a mean 1D velocity model for the Italian crust was defined and compared with those used for earthquake hypocentre locations and seismic tomographic studies by different institutions operating in the Italian area, to assess the significance of the model obtained. This model can be used in future works as input for a next joint tomographic inversion of active and passive seismic data.


Sign in / Sign up

Export Citation Format

Share Document