Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published October 2021 | Supplemental Material + Published
Journal Article Open

Seismic wavefield reconstruction using a pre-conditioned wavelet–curvelet compressive sensing approach

Abstract

The proliferation of large seismic arrays have opened many new avenues of geophysical research; however, most techniques still fundamentally treat regional and global scale seismic networks as a collection of individual time-series rather than as a single unified data product. Wavefield reconstruction allows us to turn a collection of individual records into a single structured form that treats the seismic wavefield as a coherent 3-D or 4-D entity. We propose a split processing scheme based on a wavelet transform in time and pre-conditioned curvelet-based compressive sensing in space to create a sparse representation of the continuous seismic wavefield with smooth second-order derivatives. Using this representation, we illustrate several applications, including surface wave gradiometry, Helmholtz–Hodge decomposition of the wavefield into irrotational and solenoidal components, and compression and denoising of seismic records.

Additional Information

© The Author(s) 2021. Published by Oxford University Press on behalf of The Royal Astronomical Society. This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model). Accepted 2021 June 3. Received 2021 May 23; in original form 2021 February 21. Published: 07 June 2021. The authors would like to thank the editor, Andrew Valentine, and an anonymous reviewer, for their insightful reviews that clarified several points. They would also like to thank Charles Langston for a particularly enthusiastic review. JBM acknowledges the financial support of the Origin Energy Foundation and the General Sir John Monash Foundation during his PhD studies. ZZ thanks the support from NSF CAREER award EAR 1848166 and the NSF/IUCRC GMG Center. Seismic data were processed using Obspy (Beyreuther et al. 2010), and figures were created using Matplotlib (Hunter 2007) using Cartopy for mapping (Met Office 2010). Data Availability: Data for the Mw 7.0 Loyalty Islands Earthquake were obtained through the obspy FDSN service using the Southern California Earthquake Data Center as a provider (SCEDC 2013). CSN data for the Mw 7.1 Ridgecrest Earthquake are available through the CSN website at http://csn.caltech.edu/data/. Code and data to recreate the Helmholtz–Hodge decomposition are available at https://github.com/jbmuir/HelmholtzHodgeCSN.

Attached Files

Published - ggab222.pdf

Supplemental Material - ggab222_supplemental_files.zip

Files

ggab222.pdf
Files (24.6 MB)
Name Size Download all
md5:97d83f93833b3c7b9f1aee80a1bd45e3
4.1 MB Preview Download
md5:118c586b67317d07a2d44d4a0a27384c
20.5 MB Preview Download

Additional details

Created:
August 20, 2023
Modified:
October 20, 2023