Geology Reports⌕ Search
Home/Archive/USGS

USGS · 70000320

Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method

J. Xia·R. D. Miller·Y. Xu
View original source ↗Export metadata ↓

Abstract

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (>2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data. ?? Birkhaueser 2008.

Explore related subjects

Pure and Applied GeophysicsArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J. Xia, R. D. Miller, Y. Xu. 2008-07-31. Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method. https://doi.org/10.1007/s00024-008-0364-2

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

KEEP EXPLORING

Related USGS reports

USGS2021

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↗
USGS2021

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↗