Search for: All records

In this study, I demonstrate that distributed acoustic sensing (DAS) raw strain rate data can directly be used to estimate spectral source parameters through an Empirical Green's Function (EGF) deconvolution analysis. Previously, DAS had been widely used in passive seismology to image the subsurface and analyze ground motion variations by converting strain or strain rate to particle velocity or acceleration prior to analysis. In this study, spectral analysis is applied to the PoroTomo joint DAS and seismic Nodal array in the Brady Hot Springs geothermal field to obtain source parameters for two M4 earthquakes via EGF analysis, where nearly collocated smaller events are used as an EGF to remove path and site effects. The EGF workflow is applied to raw DAS strain rate data without conversion to particle velocities and raw Nodal seismic data. The DAS and Nodal results are very consistent with similar features of spectral ratios, corner frequencies and moment ratios for the same event pairs. The uncertainty due to stacked spectral measurement is much lower on the DAS array, suggesting better stability of spectral shape measurement, possibly due to the much denser spatial sampling. The uncertainty due to model fitting is similar between DAS and Nodal arrays with slightly lower uncertainty on the DAS array. These observations demonstrate potential for directly using the strain rate measurements from DAS arrays for earthquake source characterizations.

 

Abstract High-resolution passive seismic imaging of shallow subsurface structures is often challenged by the scarcity of coherent body-wave energy in ambient noise recorded at surface stations. We show that the autocorrelation (AC) of teleseismic P-wave coda extracted from just one month of continuous recording at 5 Hz geophones can overcome this limitation. We apply this method to investigate the longitudinal subsurface bedrock structure of Unaweep Canyon—a paleovalley in western Colorado (United States) with complex evolution. Both fluvial and glacial processes have been proposed to explain the canyon’s genesis and morphology. The teleseismic P-wave coda AC retrieves zero-offset reflections from the shallow (200–500 m depth) basement interface at 120 stations along a 5 km long profile. In addition, we invert interferometrically retrieved surface-wave dispersion for the shear-wave structure of the sedimentary fill. Combined interpretation of these results and other geophysical and well data suggests an overdeepened basement geometry most consistent with glacial processes. 

It is well known that large earthquakes often exhibit significant rupture complexity such as well separated subevents. With improved recording and data processing techniques, small earthquakes have been found to exhibit rupture complexity as well. Studying these small earthquakes offers the opportunity to better understand the possible causes of rupture complexities. Specifically, if they are random or are related to fault properties. We examine microearthquakes (M < 3) in the Parkfield, California, area that are recorded by a high-resolution borehole network. We quantify earthquake complexity by the deviation of source time functions and source spectra from simple circular (omega-square) source models. We establish thresholds to declare complexity, and find that it can be detected in earthquakes larger than magnitude 2, with the best resolution above M2.5. Comparison between the two approaches reveals good agreement (>90 per cent), implying both methods are characterizing the same source complexity. For the two methods, 60–80 per cent (M 2.6–3) of the resolved events are complex depending on the method. The complex events we observe tend to cluster in areas of previously identified structural complexity; a larger fraction of the earthquakes exhibit complexity in the days following the Mw 6 2004 Parkfield earthquake. Ignoring the complexity of these small events can introduce artefacts or add uncertainty to stress drop measurements. Focusing only on simple events however could lead to systematic bias, scaling artefacts and the lack of measurements of stress in structurally complex regions.

 

Earthquake stress drop is an important source parameter that directly links to strong ground motion and fundamental questions in earthquake physics. Stress drop estimations may contain significant uncertainties due to such factors as variations in material properties and data limitations, which limit the applications of stress drop interpretations. Using a high‐resolution borehole network, we estimate stress drop for 4551 (M0‐4) earthquakes on the San Andreas Fault at Parkfield, California, between 2001 and 2016 using spectral decomposition and an improved stacking method. To evaluate the influence of spatiotemporal variations of material properties on stress drop estimations, we apply different strategies to account for spatial variations of velocity and attenuation changes, and divide earthquakes into three separate time periods to correct temporal variations of attenuation. These results show that appropriate corrections can significantly reduce the scatter in stress drop estimates, and decrease apparent depth and magnitude dependence. We find that insufficient bandwidth can cause systematic underestimation of stress drop estimates and increased scatter. The stress drop measurements from the high‐frequency borehole recordings exhibit complex stable spatial patterns with no clear correlation with the nature of fault slip, or the slip distribution of the 2004 M6 earthquake. Temporal variations are significantly smaller, less well resolved and varying spatially. They do not affect the long‐term stress drop spatial variations, suggesting local material properties may control the spatial heterogeneity of stress drop.

 

A labial salivary gland biopsy (LSGB) plays an essential role in diagnosing Sjögren's syndrome (SS), but its clinical application is limited due to its invasiveness. Here, we present a handheld single snapshot multiple-frequency demodulation-spatial frequency domain imaging (SSMD-SFDI) device for a rapid optical biopsy of labial salivary glands noninvasively. The structural and physiological parameters of lower lip mucosa were obtained from the light reflectance of the layered oral mucosa. The recovered parameters were found to correlate strongly with the progression of SS. In our pilot study on 15 healthy subjects and 183 SS patients, a support vector machine (SVM) classifier using the measured parameters distinguished healthy subjects, LSGB I, II, III, and IV patients in sequence with AUCs of 0.979, 0.898, 0.906, and 0.978, respectively. Critical structural and physiological alterations in the mucosa due to SS were further identified and used to assess its risk using an explainable neural network. The handheld spatial frequency domain imager may serve as a valuable label-free and noninvasive tool for early diagnosing and surveying SS.

 

To better quantify how injection, prior seismicity, and fault properties control rupture growth and propagation of induced earthquakes, we perform a finite‐fault slip inversion on aM_w4.0 earthquake that occurred in April 2015, the largest earthquake in an induced sequence near Guthrie, Oklahoma. The slip inversion reveals a rupture with slip patches that are anti‐correlated to the locations of prior seismicity. The prior seismicity driven by low pore pressure changes and static stress changes occurred on weaker portions of the fault, while theM_w4.0 earthquake likely ruptured relatively stronger portions of the fault. To resolve if pore pressure changes or the initial underlying stress distribution and fault strength controlled the final slip distribution of the GuthrieM_w4.0 earthquake, we compare strike‐slip events of similar magnitude from tectonically active regions and previously inactive regions. Earthquakes on reactivated faults exhibit different slip distributions than active regions, they have more prominent and well separated slip patches, a behavior often associated with faults of lower fault maturity. Pore pressure shows little effect on the distributions. These observations suggest that the initial underlying stress distribution and fault strength of reactivated faults in low deformation regions is the primary controlling factor of the slip distribution with pore pressure perturbations and earthquake interactions being secondary. Therefore, GuthrieM_w4.0 earthquakes slip distribution was enhanced by pore‐pressure perturbations and earthquake interactions by creating an optimal stress state for its failure, but the slip distribution itself is controlled by its fault's initial stress and strength state.

 

We calculate rupture directivity and velocity for earthquakes in three well‐recorded repeating sequences (2001–2016) on the San Andreas Fault at Parkfield usingPwaves from borehole recordings and the empirical Green's function method. The individual events in each sequence all show the same directivity; the largest magnitude sequence (M ~ 2.7, 8 events) ruptures unilaterally NW (at ~0.8Vs), the second sequence (M ~ 2.3, 9 events) ruptures unilaterally SE, and the smallest magnitude sequence (M ~ 2, 11 events) is less well resolved. The highly repetitive rupture suggests that geometry or material properties might control nucleation of small locked patches. The source spectra of theM ~ 2.7 sequence exhibit no detectable temporal variation. The smallerMsequences both exhibit a decrease in high‐frequency energy following theM6 earthquake that recovers with time. This could indicate a decrease in stress drop, an increase in attenuation, or a combination of the two, followed by gradual healing.

 

We compare source parameter estimates for earthquakes in the 2011 Prague Mw 5.7, Oklahoma, sequence to investigate random uncertainty and systematic bias, and resolve reliable relative variations in stress drop. Source parameters provide insight into the earthquake rupture processes but large variations between studies occur. The Prague earthquake sequence is a prime example of this, with different studies reaching contrasting interpretations of the effects of injection on source parameters. We examine the Prague earthquake sequence using a single coherent catalog for all the events detected by the Oklahoma Geological Survey (OGS) and McMahon et al. (2017). We use three principal approaches to estimate stress drop in order to understand the biases of each: a spectral decomposition method based on stacking, individual event spectral modeling, and a spectral ratio method based on highly correlated events. We also compare our results with previous studies for the Prague sequences aftershocks, as well as past results for the Mw 4.8 foreshock and Mw 4.8 aftershock and Mw 5.7 mainshock. The absolute values of stress drop vary significantly between methods, but the relative patterns remain consistent, except when low quality or low bandwidth data are included. The consistent relative patterns reveal that the stress drops of aftershocks are dependent on the fault orientation and the proximity of the events to the mainshocks slip. These results indicate that fault structure as well as past events play an important role in stress drop patterns.