Abstract:The Main objective of this study is to obtain a model that can explain the surface wave data and Bouguer gravity anomaly simultaneously and also to alleviate the non-uniqueness of surface wave inversion. We applied a novel non-linear joint inversion technique using both surface wave and gravity anomalies to obtain high resolution 3D shear-wave velocity crustal model of Central North China Craton including Shanxi graben. The crustal model helps us to understand the crust and mantle dynamics and the evolution history of the Shanxi rift. We bandpass filtered WGM2012 global Bouguer gravity anomaly data to remove the data with wavelength less than 50 km and greater than 200 km. We then projected the gravity data and surface wave data from ambient noise tomography to a Cartesian coordinate system with grid spacing 50 km x 50 km. The joint inversion was performed in this coordinate system. The inversion volume is much larger than the target region. This parameterization strategy was chosen to minimize the likelihood edge effects from gravity modelling. We used Crust 2.0 global model as our initial model in the inversion. An approximate non-linear P velocity and density relationship is derived from combing Birch's law and Nafe-Drake's curve. For each iteration of inversion, the sensitivity kernel of surface wave and gravity data were recalculated and updated until the inversion converged. We did several experiments to obtain an optimized weighting parameter to balance the influences of both datasets. A surface-wave-only inversion was also performed to compare the results of joint inversion. Both shear wave velocity models from joint inversion and surface-wave-only inversion can fit surface wave data well though only the joint inversion provides an acceptable fit to the gravity data. Gravity misfits decrease from >30 mGal to similar to 3.4 mGal after joint inversion. The velocity anomalies from both inversions all show a remarkable correlation with surface geology and upper crustal structures beneath Shanxi rift. The Luliang mountain region and Datong volcanic region all show relatively high shear wave velocity in the upper crust. Meanwhile, the basins in the southern rift show slow velocities. The southern part of the Shanxi rift is slower than the northern part in the upper crust. The Cenozoic extension started from the northern rift and then passed on to the southern rift. However, in the mid-and-lower crust beneath the Datong volcanic region, the surface-wave-only inversion resolves much slower shear wave velocity compared to the joint inversion model. The most pronounced difference occurs in the mid-and-lower crust beneath the Luliang mountain region. The joint inversion reveals a slow velocity anomaly with amplitude about 2%, while the surface-wave-only inversion model shows a fast velocity. The Datong volcanic region, north part of the Shanxi rift, might have experienced mantle upwelling since the Cenozoic time. The hot and buoyant mantle heated the upper most mantle and the crust, which caused the large-scale mid and lower crust slow velocity and negative gravity anomaly in the Datong volcanic region. The slow velocity anomaly beneath the Luliang mountain region is also related to the mantle upwelling beneath the Datong volcanic region. The Cenozoic mantle upwelling might control the formation and evolution of the Shanxi rift.

Rayleigh wave phase velocity and azimuthal anisotropy of central North China Craton derived from ambient noise tomography

High-resolution Lithospheric Structure Beneath Mainland China from Ambient Noise and Earthquake Surface-Wave Tomography

Lithospheric structure of the Ordos Block and its boundary areas inferred from Rayleigh wave dispersion

Crust and upper mantle structure of the North China Craton and the NE Tibetan Plateau and its tectonic implications

Three‐dimensional Lithospheric S Wave Velocity Model of the NE Tibetan Plateau and Western North China Craton

Rayleigh wave azimuthal anisotropy in the Southeastern Tibetan Plateau from Eikonal tomography

Rayleigh wave phase velocity tomography in Shaanxi and its adjacent regions

High resolution 3-D crustal structure beneath NE China from joint inversion of ambient noise and receiver functions using NECESSArray data

High-resolution crustal velocity structure in the Shanxi rift zone and its tectonic implications

Rayleigh phase velocity and azimuthal anisotropy from high dense seismic arrays reveal strong control of block boundary faults on crustal material migration in SW China

Crustal Azimuthal Anisotropy Beneath the Central North China Craton Revealed by Receiver Functions

3-D S-wave Velocity Structure of Crust and Upper Mantle Beneath North China

Sedimentary and Crustal Velocity Structure of Trans-North China Orogen from Joint Inversion of Rayleigh Wave Phase Velocity and Ellipticity and Some Implication for Syn-rift Volcanism

Ambient noise tomography in north China craton

Upper-mantle seismic anisotropy structure and dynamics beneath NE China inferred from SKS splitting analysis

Three-dimensional crustal model of Shanxi graben from 3D joint inversion of ambient noise surface wave and Bouguer gravity anomalies

Mantle Seismic Anisotropy Beneath NE China and Implications for the Lithospheric Delamination Hypothesis Beneath the Southern Great Xing'an Range

Three Dimensional Shear Wave Velocity Structure of the Crust and Upper Mantle Beneath China from Ambient Noise Surface Wave Tomography

Crustal Shear Wave Velocity and Radial Anisotropy in the Xiaojiang Fault Zone System (SE Tibet) Revealed by Ambient Noise Interferometry

Lithosphere Structure and Thickness Beneath the North China Craton from Joint Inversion of Ambient Noise and Surface Wave Tomography

Seismic Azimuthal Anisotropy of Northeastern Tibetan Plateau From Ambient Noise Double Beamforming Tomography: Implications for Crustal Deformation