Search USGSSearch

USGS · 70037726

A trade-off solution between model resolution and covariance in surface-wave inversion

Abstract

Regularization is necessary for inversion of ill-posed geophysical problems. Appraisal of inverse models is essential for meaningful interpretation of these models. Because uncertainties are associated with regularization parameters, extra conditions are usually required to determine proper parameters for assessing inverse models. Commonly used techniques for assessment of a geophysical inverse model derived (generally iteratively) from a linear system are based on calculating the model resolution and the model covariance matrices. Because the model resolution and the model covariance matrices of the regularized solutions are controlled by the regularization parameter, direct assessment of inverse models using only the covariance matrix may provide incorrect results. To assess an inverted model, we use the concept of a trade-off between model resolution and covariance to find a proper regularization parameter with singular values calculated in the last iteration. We plot the singular values from large to small to form a singular value plot. A proper regularization parameter is normally the first singular value that approaches zero in the plot. With this regularization parameter, we obtain a trade-off solution between model resolution and model covariance in the vicinity of a regularized solution. The unit covariance matrix can then be used to calculate error bars of the inverse model at a resolution level determined by the regularization parameter. We demonstrate this approach with both synthetic and real surface-wave data. ?? 2010 Birkh??user / Springer Basel AG.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J. Xia, Y. Xu, R. D. Miller, C. Zeng. 2010-03-16. A trade-off solution between model resolution and covariance in surface-wave inversion. https://doi.org/10.1007/s00024-010-0107-z

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Local variations in broadband sensor installations: Orientations, sensitivities, and noise levels

As seismologists continue to place more stringent demands on data quality, accurately described metadata are becoming increasingly important. In order to better constrain the orientation and sensitivities of seismometers deployed in U.S. Geological Survey networks, the Albuquerque Seismological Laboratory (ASL) has recently begun identifying true north with a fiber optic gyroscope (FOG) and has developed methodologies to constrain mid-band, vertical component sensitivity levels to less than 1% in a controlled environment. However, questions remain regarding the accuracy of this new alignment technique as well as if instrument sensitivities and background noise levels are stable when the seismometers are installed in different environmental settings. In this study, we examine the stability and repeatability of these parameters by reinstalling two high-quality broadband seismometers (Streckeisen STS-2.5 and Nanometrics T-360 Global Seismographic Network (GSN) version) at different locations around the ASL and comparing them to each other and a reference STS-6 seismometer that stayed stationary for the duration of the experiment. We find that even in different environmental conditions, the sensitivities of the two broadband seismometers stayed stable to within 0.1% and that orientations attained using the FOG are generally accurate to within a degree. However, one install was off by 5° due to a mistake made by the installation team. These results indicate that while technology and methodologies are now in place to calibrate and orient a seismometer to within 1°, human error both during the installation and while producing the metadata is often a limiting factor. Finally, we find that background noise levels at short periods (0.1–1 s) become noisier when the sensors are emplaced in unconsolidated materials, whereas the noise levels at long periods (30–100 s) are not sensitive to local geological structure on the vertical components.

Pure and Applied Geophysics

Inter-source interferometry of seismic body waves: Required conditions and examples

Seismic interferometry is widely applied to retrieve wavefields propagating between receivers. Another version of seismic interferometry, called inter-source interferometry, uses the principles of seismic reciprocity and expands interferometric applications to retrieve waves that propagate between two seismic sources. Previous studies of inter-source interferometry usually involve surface-wave and coda-wave estimations. We use inter-source interferometry to estimate the P-waves propagating between two sources rather than the estimation of surface waves and coda waves. We show that the recovered arrival times are dependent on the accuracy of the earthquake catalog of the two sources. Using inter-source interferometry, one can recover the waveform of the direct body waves and potentially reconstruct the waveform of coda waves, depending on the source-receiver geometry. The retrieval of these waveforms is accurate only when the wavefield is sampled with approximately 4 receivers per wavelength in the stationary phase zone. We show that using only receivers inside the stationary phase region for inter-source interferometry introduces the phase error of approximately 0.3 radians. In our study, we show an example of the P-wavefield reconstruction between two earthquakes using the seismic records from an array along San Andreas Fault. The retrieved P waves give a qualitative estimation of the thickness of the low-velocity zone of San Andreas Fault of approximately 4 km.

Pure and Applied Geophysics