DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression

Abstract

Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress themmore » before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.« less

Authors:
 [1];  [2];  [3]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Univ. of Wyoming, Laramie, WY (United States)
  3. National Cheng Kung Univ., Tainan City (Taiwan)
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
OSTI Identifier:
1258513
Alternate Identifier(s):
OSTI ID: 1358859
Report Number(s):
LLNL-JRNL-678678
Journal ID: ISSN 0098-3004
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Computers and Geosciences
Additional Journal Information:
Journal Volume: 93; Journal ID: ISSN 0098-3004
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING

Citation Formats

Lindstrom, Peter, Chen, Po, and Lee, En-Jui. Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression. United States: N. p., 2016. Web. doi:10.1016/j.cageo.2016.04.009.
Lindstrom, Peter, Chen, Po, & Lee, En-Jui. Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression. United States. https://doi.org/10.1016/j.cageo.2016.04.009
Lindstrom, Peter, Chen, Po, and Lee, En-Jui. Thu . "Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression". United States. https://doi.org/10.1016/j.cageo.2016.04.009. https://www.osti.gov/servlets/purl/1258513.
@article{osti_1258513,
title = {Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression},
author = {Lindstrom, Peter and Chen, Po and Lee, En-Jui},
abstractNote = {Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.},
doi = {10.1016/j.cageo.2016.04.009},
journal = {Computers and Geosciences},
number = ,
volume = 93,
place = {United States},
year = {Thu May 05 00:00:00 EDT 2016},
month = {Thu May 05 00:00:00 EDT 2016}
}

Journal Article:

Citation Metrics:
Cited by: 25 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Seismic Tomography
journal, October 1984


Adaptive Local Cosine transform for Seismic Image Compression
conference, December 2006


Bandwidth Reduction through Multithreaded Compression of Seismic Images
conference, May 2011

  • Aqrawi, Ahmed A.; Elster, Anne C.
  • Distributed Processing, Workshops and Phd Forum (IPDPSW), 2011 IEEE International Symposium on Parallel and Distributed Processing Workshops and Phd Forum
  • DOI: 10.1109/IPDPS.2011.330

Low bit-rate efficient compression for seismic data
journal, December 2001

  • Averbuch, A. Z.; Meyer, F.; Stromberg, J. -O.
  • IEEE Transactions on Image Processing, Vol. 10, Issue 12
  • DOI: 10.1109/83.974565

Finite-Moment Tensor of the 3 September 2002 Yorba Linda Earthquake
journal, June 2005

  • Chen, P.
  • Bulletin of the Seismological Society of America, Vol. 95, Issue 3
  • DOI: 10.1785/0120040094

Full 3D Tomography for the Crustal Structure of the Los Angeles Region
journal, August 2007

  • Chen, P.; Zhao, L.; Jordan, T. H.
  • Bulletin of the Seismological Society of America, Vol. 97, Issue 4
  • DOI: 10.1785/0120060222

Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods
journal, July 2007


The TeraShake Computational Platform for Large-Scale Earthquake Simulations
book, January 2009


Seismic Data Compression Using 2D Lifting-Wavelet Algorithms
journal, January 2015


Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint methods
journal, December 2009

  • Fichtner, Andreas; Kennett, Brian L. N.; Igel, Heiner
  • Geophysical Journal International, Vol. 179, Issue 3, p. 1703-1725
  • DOI: 10.1111/j.1365-246X.2009.04368.x

Group testing for image compression
journal, August 2002


Seismic data compression
conference, January 1990

  • Jonsson, S. B.; Spanias, A. S.
  • Ninth Annual International Phoenix Conference on Computers and Communications. 1990 Conference Proceedings
  • DOI: 10.1109/PCCC.1990.101632

Full-3-D tomography for crustal structure in Southern California based on the scattering-integral and the adjoint-wavefield methods: F3DT FOR SOUTHERN CALIFORNIA CRUST
journal, August 2014

  • Lee, En-Jui; Chen, Po; Jordan, Thomas H.
  • Journal of Geophysical Research: Solid Earth, Vol. 119, Issue 8
  • DOI: 10.1002/2014JB011346

Testing Waveform Predictions of 3D Velocity Models against Two Recent Los Angeles Earthquakes
journal, October 2014

  • Lee, E. -J.; Chen, P.; Jordan, T. H.
  • Seismological Research Letters, Vol. 85, Issue 6
  • DOI: 10.1785/0220140093

An optimized parallel LSQR algorithm for seismic tomography
journal, December 2013


Fast and Efficient Compression of Floating-Point Data
journal, September 2006

  • Lindstrom, Peter; Isenburg, Martin
  • IEEE Transactions on Visualization and Computer Graphics, Vol. 12, Issue 5
  • DOI: 10.1109/TVCG.2006.143

Fixed-Rate Compressed Floating-Point Arrays
journal, December 2014

  • Lindstrom, Peter
  • IEEE Transactions on Visualization and Computer Graphics, Vol. 20, Issue 12
  • DOI: 10.1109/TVCG.2014.2346458

i PACKMAN: high-quality, low-complexity texture compression for mobile phones
conference, January 2005

  • Ström, Jacob; Akenine-Möller, Tomas
  • Proceedings of the ACM SIGGRAPH/EUROGRAPHICS conference on Graphics hardware - HWWS '05
  • DOI: 10.1145/1071866.1071877

Seismic tomography of the southern California crust based on spectral-element and adjoint methods
journal, January 2010


Comparison of wavelet image coding schemes for seismic data compression
conference, October 1997

  • Vassiliou, Anthony A.; Wickerhouser, Mladen V.
  • Optical Science, Engineering and Instrumentation '97, SPIE Proceedings
  • DOI: 10.1117/12.279684

Seismic data compression based on integer wavelet transform
journal, November 2004


Seismic data Compression Methods
journal, August 1974


Adaptive Seismic Data Compression Using Wavelet Packets
conference, July 2006

  • Wu, W.; Yang, Z.; Qin, Q.
  • 2006 IEEE International Symposium on Geoscience and Remote Sensing
  • DOI: 10.1109/IGARSS.2006.202

Fast Lossless Compression of Seismic Floating-Point Data
conference, May 2009

  • Xie, Xing; Qin, Qianqing
  • 2009 International Forum on Information Technology and Applications (IFITA)
  • DOI: 10.1109/IFITA.2009.556

Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models
journal, December 2005

  • Zhao, L.
  • Bulletin of the Seismological Society of America, Vol. 95, Issue 6
  • DOI: 10.1785/0120050081

Strain Green's Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies
journal, October 2006

  • Zhao, L.; Chen, P.; Jordan, T. H.
  • Bulletin of the Seismological Society of America, Vol. 96, Issue 5
  • DOI: 10.1785/0120050253

A fast compression algorithm for seismic data from non-cable seismographs
conference, October 2012

  • Zheng, Fan; Liu, Shufen
  • 2012 World Congress on Information and Communication Technologies (WICT)
  • DOI: 10.1109/WICT.2012.6409260

Works referencing / citing this record:

Data Reduction Techniques for Simulation, Visualization and Data Analysis: Survey on Scientific Data Reduction Techniques
journal, March 2018

  • Li, S.; Marsaglia, N.; Garth, C.
  • Computer Graphics Forum, Vol. 37, Issue 6
  • DOI: 10.1111/cgf.13336

Compression Challenges in Large Scale Partial Differential Equation Solvers
journal, September 2019

  • Götschel, Sebastian; Weiser, Martin
  • Algorithms, Vol. 12, Issue 9
  • DOI: 10.3390/a12090197

An Efficient Encoding Algorithm Using Local Path on Huffman Encoding Algorithm for Compression
journal, February 2019

  • Erdal, Erdal; Ergüzen, Atilla
  • Applied Sciences, Vol. 9, Issue 4
  • DOI: 10.3390/app9040782