Shallow Ice-Sheet Composite Structure Revealed by Seismic Imaging Near the West Antarctic Ice Sheet (WAIS) Divide Camp

• Published in: Journal of Geophysical Research: Earth Surface
• Main Authors: Zhang, Zhendong, Nakata, Nori, Karplus, Marianne, Kaip, Galen, Yi, Jia
• Other Authors: Massachusetts Institute of Technology, Cambridge, MA, USA, Lawrence Berkeley National Laboratory, Berkeley, CA, USA, University of Texas at El Paso, El Paso, TX, USA, China Earthquake Disaster Prevention Center, China Earthquake Administration, Beijing, China
• Format: Article in Journal/Newspaper
• Language: unknown
• Published: American Geophysical Union (AGU) 2022
• Subjects: Antarctic; The Antarctic; West Antarctic Ice Sheet; Antarc*; Ice Sheet; Ice Shelves
• Online Access: http://hdl.handle.net/10754/689942; https://doi.org/10.1029/2022JF006777

Small-scale polar ice-sheet composite structures can influence coarse-grained climate models. For example, surface melting can hydro-fracture the ice and may eventually lubricate the base. The drainage process will change the petrophysical properties of the subsurface and is predictable from high-resolution seismic imaging. Ray-based imaging methods, though limited by the approximations made, have mapped seismic velocities in 1D, but have yet to offer definitive information about the lateral heterogeneity. Here we use the wave-equation-based seismic imaging method to estimate shear-wave velocities and to image an englacial reflector in 2D. We use 7-day ambient-noise data recorded by a linear array near the West Antarctic Ice Sheet (WAIS) Divide drilling site to perform the study. We obtain 2D vertical and horizontal shear-wave velocity models and observe a lateral variation of the radial anisotropy along the acquisition line. The inverted velocity model is further used to estimate the thicknesses of firn-air and firn layers. We also observe evident SH reflections and calculate the reflector image using robust zero-offset imaging and a more precise reverse time migration imaging method. The imaged reflector is at about 1,700 m depth and is dipping to the southwest. We anticipate that a porosity change at that depth may cause this englacial seismic discontinuity near WAIS Divide. Our results demonstrate the feasibility of estimating the petrophysical properties of polar ice using seismic methods. We anticipate our work to be further used for understanding ice dynamics and thermodynamics, such as the stability of ice shelves, the retreat of the Antarctic ice sheet, and restoring the paleo-sedimentary environment. The authors thank Tariq Alkhalifah for helpful discussions. The authors thank the HPC team at KAUST for providing guidance on using IBEX and Shaheen clusters. The computing for this project was partly performed at the OU Supercomputing Center for Education and Research (OSCER) at the University of ...