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 January 2021 | Supplemental Material + Published
Journal Article Open

Direct P-wave anisotropy measurements at Homestake Mine: implications for wave propagation in continental crust

Abstract

We measured anisotropic seismic properties of schists of the Homestake Formation located at a depth of 1478 m in the Sanford Underground Research Facility (SURF) in the Black Hills of South Dakota, USA. We deployed a 24-element linear array of three-component geophones in an area in the Homestake Mine called 19-ledge. An airless jackhammer source was used to shoot two profiles: (1) a walkaway survey to appraise any distance dependence and (2) a fan shot profile to measure variations with azimuth. Slowness estimates from the fan shot profile show a statistically significant deviation with azimuth with the expected 180° variation with azimuth. We measured P-wave particle motion deviations from data rotated to ray coordinates using three methods: (1) a conventional principal component method, (2) a novel grid search method that maximized longitudinal motion over a range of search angles and (3) the multiwavelet method. The multiwavelet results were computed in two frequency bands of 200–600 and 100–300 Hz. Results were binned by azimuth and averaged with a robust estimation method with error bars estimated by a bootstrap method. The particle motion results show large, statistically significant variations with azimuth with a 180° cyclicity. We modelled the azimuthal variations in compressional wave speed and angular deviation from purely longitudinal particle motion of P-waves using an elastic tensor method to appraise the relative importance of crystalline fabric relative to fracturing parallel to foliation. The model used bulk averages of crystal fabric measured for an analogous schist sample from southeast Vermont rotated to the Homestake Formation foliation directions supplied by SURF from old mine records. We found with average crustal crack densities crack induced anisotropy had only a small effect on the observables. We found strong agreement in the traveltime data. The observed amplitudes of deviations of P particle motion showed significantly larger variation than the model predictions and a 20° phase shift in azimuth. We attribute the inadequacies of the model fit to the particle motion data to inadequacies in the analogue rock and/or near receiver distortions from smaller scale heterogeneity. We discuss the surprising variability of signals recorded in this experimental data. We show clear examples of unexplained resonances and unexpected variations on a scale much smaller than a wavelength that has broad implications for wave propagation in real rocks.

Additional Information

© The Author(s) 2020. 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 2020 August 14. Received 2020 July 22; in original form 2019 November 17. This paper would not have been possible without the amazing support the staff at SURF provided to collect this experimental data. We especially thank Tom Regan and Neil Engle for keeping us safe underground and being great helpers. Jaret Heise was invaluable in coordinating the experimental work with SURF personnel and getting us access to geologic data from Homestake. Sincere thanks to the DUGL team who fielded the broad-band experiment and to which this experiment was an added component. We also thank Terry Stigall of Indiana University who was a key player in keeping instrumentation working and coming up with her usual creative solutions to problems. Sincere thanks to Jeffrey Park and an anonymous reviewer for constructive comments that helped improve this paper significantly. This paper emerged from the B.S. thesis of JA. The text was extended and revised by all authors led by GLP. The multiwavelet code is the legacy of a body of work by GLP spanning two decades. SJB did the model calculations. The seismic instruments were provided by the Incorporated Research Institutions for Seismology (IRIS) through the PASSCAL Instrument Center at New Mexico Tech. Data used are available through the IRIS Data Management Center. The facilities of the IRIS Consortium are supported by the National Science Foundation under Cooperative Agreement EAR-1261681 and the DOE National Nuclear Security Administration. This research was supported by the National Science Foundation under EAR-1526894.

Attached Files

Published - ggaa416.pdf

Supplemental Material - ggaa416_supplemental_files.zip

Files

ggaa416.pdf
Files (14.3 MB)
Name Size Download all
md5:0692f42d4b85185892ac36d93355f94a
2.1 MB Preview Download
md5:b9f2abe883686bb7f64cf5e8a83bcf10
12.1 MB Preview Download

Additional details

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