SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This is a CSV spreadsheet containing UTM and Latitude and Longitude coordinates and elevations for Wells 78-32, 58-32, and 16A(78)-32 and BOR1, BOR2, BOR3, FOR2, FOR5, FORK, FORU, and FORW seismic stations. These are from a GPS survey conducted by the Utah Geological Survey in June, 2021.
This catalog describes the seismicity associated with the 2019 stimulation at Utah FORGE. Containing both matched-filter detections (Dzubay et al., 2022) and Schlumberger-recorded events (detected with a 12-level geophone string), the final combined catalog contains a total of 534 microseismic events spanning -2.0 Mw to -0.1 Mw. Users may differentiate between SLB and MF events using the fact that SLB event magnitudes are recorded to a higher level of precision (MF mags determined using relative amplitude ratios). Users should be wary of locations and depths (measured from sea level) for MF events, as all detections were assigned the same locations as their template events.
Abstract We investigate the merits of the more recently developed finite‐frequency approach to tomography against the more traditional and approximate ray theoretical approach for state of the art seismic models developed for western North America. To this end, we employ the spectral element method to assess the agreement between observations on real data and measurements made on synthetic seismograms predicted by the models under consideration. We check for phase delay agreement as well as waveform cross‐correlation values. Based on statistical analyses on S wave phase delay measurements, finite frequency shows an improvement over ray theory. Random sampling using cross‐correlation values identifies regions where synthetic seismograms computed with ray theory and finite‐frequency models differ the most. Our study suggests that finite‐frequency approaches to seismic imaging exhibit measurable improvement for pronounced low‐velocity anomalies such as mantle plumes.
Using an updated data set of ballistic PKIKP travel time data at antipodal distances, we test different models of anisotropy in the Earth's innermost inner core (IMIC) and obtain significantly better fits for a fast axis aligned with Earth's rotation axis, rather than a quasi-equatorial direction, as proposed recently. Reviewing recent results on the single crystal structure and elasticity of iron at core conditions, we find that an hcp structure with the fast c axis parallel to Earth's rotation is more likely but a body-centered cubic structure with the [111] axis aligned in that direction results in very similar predictions for seismic anisotropy. These models are therefore not distinguishable based on current seismological data. In addition, to match the seismological observations, the inferred strength of anisotropy in the IMIC (6–7%) implies almost perfect alignment of iron crystals, an intriguing, albeit unlikely situation, especially in the presence of heterogeneity, which calls for further studies. Fast axis of anisotropy in the central part of the inner core aligned with Earth's axis of rotation Lastly, the structure of iron in the inner core is most likely hcp, not bcc Not currently possible to distinguish between hcp and bcc structures from seismic observations
Subsea ice-bearing permafrost (IBPF) and associated gas hydrate in the Arctic have been subject to a warming climate and saline intrusion since the last transgression at the end of the Pleistocene. The consequent degradation of IBPF is potentially associated with significant degassing of dissociating gas hydrate deposits. Previous studies interpreted the distribution of subsea permafrost on the U.S. Beaufort continental shelf based on geographically sparse data sets and modeling of expected thermal history. The most cited work projects subsea permafrost to the shelf edge (~100 m isobath). This study uses a compilation of stacking velocity analyses from ~100,000 line-km of industry-collected multichannel seismic reflection data acquired over 57,000 km 2 of the U.S. Beaufort shelf to delineate continuous subsea IBPF. Gridded average velocities of the uppermost 750 ms two-way travel time range from 1475 to 3110 m s -1 . The monotonic, cross-shore pattern in velocity distribution suggests that the seaward extent of continuous IBPF is within 37 km of the modern shoreline at water depths < 25 m. These interpretations corroborate recent Beaufort seismic refraction studies and provide the best, margin-scale evidence that continuous subsea IBPF does not currently extend to the northern limits of the continental shelf.
Here, we conducted an experiment in Pahrump, Nevada, in June 2017, where artificial seismic signals were created using a seismic hammer, and the possibility of detecting them from their acoustic signature was examined. In this work, we analyze the pressure signals recorded by highly sensitive barometers deployed on the ground and on tethers suspended from balloons. Our signal processing results show that wind noise experienced by a barometer on a free-flying balloon is lower compared to one on a moored balloon. This has never been experimentally demonstrated in the lower troposphere. While seismoacoustic signals were not recorded on the hot air balloon platform owing to operational challenges, we demonstrate the detection of seismoacoustic signals on our moored balloon platform. Our results have important implications for performing seismology in harsh surface environments such as Venus through atmospheric remote sensing.
Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.
For characterizing geothermal systems, it is important to have clear images of steeply-dipping fault zones because they may confine the boundaries of geothermal reservoirs and influence hydrothermal flow. Elastic reverse-time migration (ERTM) is the most promising tool for subsurface imaging with multicomponent seismic data. However, conventional ERTM usually generates significant artifacts caused by the cross correlation of undesired wavefields and the polarity reversal of shear waves. In addition, it is difficult for conventional ERTM to directly image steeply-dipping fault zones. We develop a new ERTM imaging method in this paper to reduce these artifacts and directly image steeply-dipping fault zones. In our new ERTM method, forward-propagated source wavefields and backward-propagated receiver wavefields are decomposed into compressional (P) and shear (S) components. Furthermore, each component of these wavefields is separated into left- and right-going, or downgoing and upgoing waves. The cross correlation imaging condition is applied to the separated wavefields along opposite propagation directions. For converted waves (P-to-S or S-to-P), the polarity correction is applied to the separated wavefields based on the analysis of Poynting vectors. Numerical imaging examples of synthetic seismic data demonstrate that our new ERTM method produces high-resolution images of steeply-dipping fault zones.
Relying primarily on 3D-seismic data, an assessment of the major structural elements of the Miocene section, northern Gulf of Mexico (Texas and Louisiana), was initiated here to identify prospective subsurface areas for permanent geologic storage of anthropogenic CO 2 . Identifying and mapping fault planes and key stratigraphic surfaces in the seismic data helped identify several areas that may be suitable for Carbon Capture and Storage (CCS) projects of scale. Prospective sites that included proximity to anthropogenic sources of CO 2 , RMS amplitude intervals that indicate confining zones that greatly retard vertical migration of buoyant CO 2 , and structural closures were ranked on a regional scale. Local site (Site 1) assessment used hydrocarbon-based and volumetric-based methods to estimate the amount of CO 2 that can be safely injected into a prospective site’s permanent storage reservoirs.
Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.
Wind Energy is the one of the most promising renewable energy. Suction bucket foundation is considered to be a viable type of wind turbine foundation. Soil liquefaction caused by earthquakes at offshore seismic active area may lead to a significant degradation of soil strength and stiffness. In this study, nine centrifuge tests were carried out to investigate the seismic response of suction bucket foundation under earthquake loading. Both dry and saturate soil conditions were considered in tests. The geometric design of five suction bucket models considered the bucket diameter, penetration depth, and modified buckets with inside compartments. It was found that soil underlying and near the bucket foundation shown a better ability to resist liquefaction in saturated tests comparing to free field while no significant differences were observed in dry tests. The five bucket models performed quite differently, which demonstrated the aspect ratio effects and inside-bucket compartment effects. The results provide insight into optimized design of suction bucket foundation for wind turbine.
Ambient seismic recordings taken at broad locations across Ross Ice Shelf and a dense array near West Antarctic Ice Sheet (WAIS) Divide, Antarctica, show pervasive temporally variable resonance peaks associated with trapped seismic waves in near-surface firn layers. These resonance peaks feature splitting on the horizontal components, here interpreted as frequency-dependent anisotropy in the firn and underlying ice due to several overlapping mechanisms driven by ice flow. Frequency peak splitting magnitudes and fast/slow axes were systematically estimated at single stations using a novel algorithm and compared with good agreement with active source anisotropy measurements at WAIS Divide determined via active sources recorded on a 1 km circular array. The approach was further applied to the broad Ross Ice Shelf (RIS) array, where anisotropy axes were directly compared with visible surface features and ice shelf flow lines. The near-surface firn, depicted by anisotropy above 30 Hz, was shown to exhibit a novel plastic stretching mechanism of anisotropy, whereby the fast direction in snow aligns with accelerating ice shelf flow.
Omphacite is a major mineral phase of eclogite, which provides the main driving force for the slab subduction into the Earth's interior. We have measured the single-crystal elastic moduli of omphacite at high pressures for the first time up to 18 GPa at ambient temperature using Brillouin spectroscopy. A least squares fit of the velocity-pressure data to the third-order finite strain equation of state yields K S0 ' = 4.5 (3), G 0 ' = 1.6 (1) with ρ 0 = 3.34 (1) g/cm 3 , K S0 = 123 (3) GPa, and G 0 = 74 (2) GPa. In addition, the synchrotron single-crystal X-ray diffraction data have been collected up to 18 GPa and 700 K. The fitting to Holland-Powell thermal-pressure equation of state yields K T0 ' = 4.6 (5) and α 0 = 2.7 (8) × 10 -5 K -1 . Based on the obtained thermoelastic parameters of omphacite, the anisotropic seismic velocities of eclogite are modeled and compared with pyrolite between 200 and 500 km. The largest contrast between the eclogite and pyrolite in terms of seismic properties is observed between ~310 and 410 km.
Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.
There are an infinite number of locations where seismic waves interact with Earth's surface, yet seismologists can only record a finite set of ground motions. This requires assumptions about how known measurements or predictions can be extrapolated to nearby locations. We test to what extent extrapolation is viable using an array of 1,108 geophones installed 10–20 m apart within a 0.36 km 2 study area for 1 month atop the San Jacinto Fault in Southern California. Measuring ground motions from 38 small local earthquakes, we observe peak ground velocity (PGV) variations up to 22% across the array. We find the basin structure consistently traps energy, amplifying the seismic velocities, and that the damaged region within the fault reduces the PGVs. Comparing surface and borehole measurements (148 m deep), we find PGVs recorded at the surface are 3–10 times greater than the measurements at depth. We attribute these differences to amplification from soft sediments within the basin. We show that by leveraging information from >1,000 measurements per earthquake the ground motion velocities can be tightly constrained. These results suggest the spatial extrapolation of PGVs in regions of known faults and basins, even across a small area, should be done with caution.
Abstract Peterie et al. question one observation in our paper: associating pressure increases to injection volumes at distances of up to 25 km from an injection well. In this reply, we show that the comment misunderstands our analysis and the evidence that led to this conclusion. We also show that gauge‐depth‐corrected pressures, used by the authors to produce statewide pressure maps, are discrepant with the static fluid level data, provided in our original compilation and analysis. The discrepancies are a result of the pressure correction method employed, which naïvely substitutes formation pressure for bottomhole pressure to calculate wellbore fluid density. Their linearly interpolated pressure maps, based on sparse data, contain interpolation and extrapolation artifacts that contradict injection trends in the state, the Theis solution, and the superposition principle. We reiterate that pressure and static fluid level increases in Class I wells existed prior to 2013, most notably in central Kansas, where recent earthquakes are cited in the comment as evidence of a pressure plume emanating from the Kansas‐Oklahoma border, 90 km away. We show that the space‐time pattern of seismicity in this area is inconsistent with a northward propagating pressure plume and, instead, seismicity appears to be centered on and near a cluster of high‐rate injection wells, two of which are among the highest rate wells in the state. These observations, along with recent M4 + earthquakes during continued decreases in wastewater injection in southern Kansas and northern Oklahoma, question the usefulness of the comment for understanding and managing societally significant earthquakes.
Abstract On September 12, 2018 and April 14, 2019 shallow (<2.5 km from the surface, <1 km relative to sea level) earthquakes ( M w 4.0 and 4.1, respectively) struck a relatively aseismic area in south‐central Utah, within the Black Rock Desert volcanic field. Seismic waveforms are dominated by long‐period energy, no clear direct S‐wave arrivals, and signal durations of ∼100 s. Using co‐seismic interferograms, we detect clear and distinguishable surface deformation confirming the shallowness of the two mainshocks, while full moment tensor solutions rule out an isotropic source, such as an explosion or collapse. For the April 2019 sequence, we build an earthquake catalog using a dense temporary geophone array. The enhanced earthquake catalog shows a rapidly decaying aftershock sequence, with 85% of aftershocks occurring during the first three hours. Combining these observations, we argue that these earthquake sequences are related to the Black Rock Desert volcanic field.
Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.