We use P wave receiver functions from the western U.S. and adjacent regions to construct a receiver function wavefield interpolation scheme that helps to equalize the lateral sampling of the receiver functions and the surface wave dispersion and to greatly simplify the receiver functions. Spatial interpolation and smoothing suppress poorly sampled and difficult to interpret back azimuthal variations and allow the extraction of the first‐order features in the receiver function wavefield, including observations from several ray parameter ranges. We combine the interpolated receiver functions with Rayleigh wave dispersion estimates and surface gravity observations to estimate the 3‐D shear wave speed beneath the region. Speed variations in the 3‐D model correlate strongly with expected geologic variations and illuminate broad‐scale features of the western U.S. crust and upper mantle. The model is smooth, self‐consistent, and demonstrates the compatibility of the interpolated receiver functions and dispersion observations.

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. Here we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

The important task of tracking seismic activity requires both sensitive detection and accurate earthquake location. Approximate earthquake locations can be estimated promptly and automatically; however, accurate locations depend on precise seismic phase picking, which is a laborious and time‐consuming task. We adapted a deep neural network (DNN) phase picker trained on local seismic data to mesoscale hydraulic fracturing experiments. We designed a novel workflow, transfer learning‐aided double‐difference tomography, to overcome the 3 orders of magnitude difference in both spatial and temporal scales between our data and data used to train the original DNN. Only 3,500 seismograms (0.45% of the original DNN data) were needed to retrain the original DNN model successfully. The phase picks obtained with transfer‐learned model are at least as accurate as the analyst's and lead to improved event locations. Moreover, the effort required for picking once the DNN is trained is a small fraction of the analyst's.

The final version of the above article was posted prematurely on 16 July 2021, owing to a technical error. The final, corrected version of record will be made fully available at a later date.

• Mesoscale hydraulic fracturing in crystalline rock observed with multi-geophysical sensor array at close proximity • Created fracture network consists of multi-strand hydraulic fractures and reactivated preexisting structures • Hydraulic fracture growth is strongly influenced by rock fabric, pre-existing fractures,