skip to main content


Title: Borehole fibre-optic seismology inside the Northeast Greenland Ice Stream
SUMMARY

Ice streams are major contributors to ice sheet mass loss and sea level rise. Effects of their dynamic behaviour are imprinted into seismic properties, such as wave speeds and anisotropy. Here, we present results from a distributed acoustic sensing (DAS) experiment in a deep ice-core borehole in the onset region of the Northeast Greenland Ice Stream, with focus on phenomenological and methodological aspects. A series of active seismic surface sources produced clear recordings of the P and S wavefield, including internal reflections, along a 1500 m long fibre-optic cable that was placed into the borehole. The combination of nonlinear traveltime tomography with a firn model constrained by multimode surface wave data, allows us to invert for P and S wave speeds with depth-dependent uncertainties on the order of only 10 m s−1, and vertical resolution of 20–70 m. The wave speed model in conjunction with the regularly spaced DAS data enable a straightforward separation of internal upward reflections followed by a reverse-time migration that provides a detailed reflectivity image of the ice. While the differences between P and S wave speeds hint at anisotropy related to crystal orientation fabric, the reflectivity image seems to carry a pronounced climatic imprint caused by rapid variations in grain size. Further improvements in resolution do not seem to be limited by the DAS channel spacing. Instead, the maximum frequency of body waves below ∼200 Hz, low signal-to-noise ratio caused by poor coupling, and systematic errors produced by the ray approximation, appear to be the leading-order issues. Among these, only the latter has a simple existing solution in the form of full-waveform inversion. Improving signal bandwidth and quality, however, will likely require a significantly larger effort in terms of both sensing equipment and logistics.

 
more » « less
NSF-PAR ID:
10469702
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
235
Issue:
3
ISSN:
0956-540X
Format(s):
Medium: X Size: p. 2430-2441
Size(s):
["p. 2430-2441"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Full waveform inversion (FWI) and distributed acoustic sensing (DAS) are powerful tools with potential to improve how seismic site characterization is performed. FWI is able to provide true 2D or 3D images of the subsurface by inverting stress wave recordings collected over a wide variety of scales. DAS can be used to efficiently collect high-resolution stress wave recordings from long and complex fiber optic arrays and is well-suited for large-scale site characterization projects. Due to the relative novelty of combining FWI and DAS, there is presently little published literature regarding the application of FWI to DAS data for near-surface (depths < 30 m) site characterization. We perform 2D FWI on DAS data collected at a well-characterized site using four different, site-specific 1D and 2D starting models. We discuss the unique benefits and challenges associated with inverting DAS data compared to traditional geophone data. We examine the impacts of using the various starting models on the final 2D subsurface images. We demonstrate that while the inversions performed using all four starting models are able to fit the major features of the DAS waveforms with similar misfit values, the final subsurface images can be quite different from one another at depths greater than about 10 m. As such, the best representation(s) of the subsurface are evaluated based on: (1) their agreement with borehole lithology logs that were not used in the development of the starting models, and (2) consistency at shallow depths between the final inverted images derived from multiple starting models. Our results demonstrate that FWI applied to DAS data has significant potential as a tool for near-surface site characterization while also emphasizing the significant impact that starting model selection can have on FWI results. 
    more » « less
  2. Abstract

    We conducted a vertical seismic profile (VSP) in the borehole of International Ocean Discovery Program/International Continental Scientific Drilling Program Expedition 364 Site M0077 to better understand the nature of the seismic reflectivity and the in situ seismic properties associated with the Chicxulub impact structure peak ring. Extraction of the up‐going wavefield from the VSP shows that a strong seismic reflection event imaged in seismic reflection data results from discontinuities in the elastic impedanceZ(the product of density and wave speed) at the top and bottom of a zone of hydrothermally altered melt‐bearing polymict breccia (suevite) that are characterized by anomalously lowZ. Below this strong carbonate/suevite reflection event, the upgoing seismic wavefield is chaotic, indicating high levels of scattering from the suevites and underlying melt rocks and shocked granitoids of the peak ring, in contrast to the clear coherent reflections throughout the overlying Cenozoic sediments. We extract shear wave speeds, which, together with those provided from the complementary sonic log and densities from core scanning, allowed determination ofVP/VSand Poisson's ratiov. These values are anomalously high relative to comparable terrestrial lithologies. We also calculate a variety of damage parameters for the disrupted peak ring granitoids. These values may assist in linking seismic observations to shock levels that are necessary to calibrate current impact models and may also be useful in assessing levels of fracturing within major fault zones.

     
    more » « less
  3. ABSTRACT

    Distributed acoustic sensing (DAS) technology is an emerging field of seismic sensing that enables recording ambient noise seismic data along the entire length of a fiber-optic cable at meter-scale resolution. Such a dense spatial resolution of recordings over long distances has not been possible using traditional methods because of limited hardware resources and logistical concerns in an urban environment. The low spatial resolution of traditional passive seismic acquisition techniques has limited the accuracy of the previously generated velocity profiles in many important urban regions, including the Reno-area basin, to the top 100 m of the underlying subsurface. Applying the method of seismic interferometry to ambient noise strain rate data obtained from a dark-fiber cable allows for generating noise cross correlations, which can be used to infer shallow and deep subsurface properties and basin geometry. We gathered DAS ambient noise seismic data for this study using a 12 km portion of a dark-fiber line in Reno, Nevada. We used gathered data to generate and invert dispersion curves to estimate the near-surface shear-wave velocity structure. Comparing the generated velocity profiles with previous regional studies shows good agreement in determining the average depth to bedrock and velocity variations in the analyzed domain. A synthetic experiment is also performed to verify the proposed framework further and better understand the effect of the infrastructural cover along the cable. The results obtained from this research provide insight into the application of DAS using dark-fiber lines in subsurface characterization in urban environments. It also discusses the potential effects of the conduit that covers such permanent fiber installations on the produced inversion results.

     
    more » « less
  4. SUMMARY

    Sedimentary basins can strongly amplify seismic waves from earthquakes. To better predict future ground motions, detailed knowledge of the sediment thickness and internal structure of basins is required. We image the sediment-to-bedrock interface of the Kanto Basin in Japan using the P-wave reflectivity response from earthquake and ambient seismic noise autocorrelation functions (ACFs) at 286 shallow borehole stations. Earthquake ACFs are computed using P-wave records from 50 Mw 6+ teleseismic events. Noise ACFs are obtained using 1 month of continuous data. Both methods are used to retrieve P-wave traveltimes between the surface and the bedrock interface and map the basin basement geometry. Our prediction of the basement depth agrees generally well with that from a reference velocity model, except for smoother variations in the central part of the basin. Using full-wavefield simulations, we show that the nature of the autocorrelated wavefield has a significant impact on the shape of the ACF waveforms and that earthquake ACFs yield more accurate results in the Kanto Basin.

     
    more » « less
  5. SUMMARY

    Distributed acoustic sensing (DAS) is a relatively new technology for recording the propagation of seismic waves, with promising applications in both engineering and geophysics. DAS's ability to simultaneously collect high spatial resolution waveforms over long arrays suggests that it is well-suited for near-surface imaging applications such as 2-D multichannel analysis of surface waves (MASWs), which require, at a minimum, long, linear arrays of single-component receivers. The 2-D MASW method uses a large number of sensor subarrays deployed along a linear alignment to produce 1-D shear-wave velocity (VS) profiles beneath each subarray. The 1-D VS profiles are then combined to form a pseudo-2-D VS image beneath the entire linear alignment that can be used for the purpose of identifying and characterizing lateral variations in subsurface layering. Traditionally, 2-D MASW is conducted using arrays consisting of either 24 or 48 geophones. While additional receivers could easily be incorporated into the testing configuration, it is rare for researchers and practitioners to have access to greater numbers of seismographs and geophones. When a limited number of geophones are available for deployment, there is a need to pre-determine the geophone spacing and subarray length prior to field data acquisition. Studies examining how the choice of subarray geometry impacts the resulting pseudo-2-D VS cross-sections have been largely limited to synthetic data. In response, this study utilizes DAS data to examine the effects of using various subarray lengths by comparing pseudo-2-D VS cross-sections derived from active-source waveforms collected at a well-characterized field site. DAS is particularly useful for 2-D MASW applications because the subarray geometry does not need to be determined prior to field data acquisition. We organize the DAS waveforms into multiple sets of overlapping MASW subarrays of differing lengths, ranging from 11 to 47 m, along the same alignment, allowing for direct comparison of the derived pseudo-2-D VS results at the site. We show that the length of the individual MASW subarrays has a significant effect on the resulting VS cross-sections, including the resolved location of a strong impedance contrasts at our study site, and evaluate the results relative to ground truth from invasive testing. Our results suggest that the choice of subarray length is important and should be carefully chosen to meet project-specific goals. Furthermore, analysts may consider using multiple subarray geometries during the data processing stage, as is made possible by DAS, to properly evaluate the uncertainty of 2-D MASW results. This study demonstrates the potential of using DAS to collect data for 2-D MASW in a manner that is efficient and flexible, and can be easily scaled up for use with very long arrays.

     
    more » « less