Search NASA⌕ Search

SEARCH · Search NASA

Results for “fault seismicity”

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.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Analysis of the 2021 Milford, Utah earthquake swarm: Enhanced earthquake catalog and migration patterns

An earthquake swarm occurred in Spring 2021 in south-central Utah near the town of Milford. The University of Utah Seismograph Stations located 125 earthquakes between March 19 and May 10 with magnitudes ranging from 0.5 to 3.2. We implement a matched-filter technique in order to identify additional earthquakes that went undetected during the routine network location. The 125 network-located earthquakes are used as templates and are cross-correlated with continuous data for the dates Feb 17—June 10. This time period corresponds to approximately 1 month before the earthquake swarm began through 1 month after it ended. For the matched-filter analysis, we rely heavily on station FOR1, which is located within 5 km of most template events. Four other stations within 20–30 km of most template events provide a supplement to the closest station. The matched-filter implementation results in the detection of over 600 earthquakes in addition to the original 125 catalog events. This is one of the largest swarms ever recorded in Utah, and no previous large swarms have been recorded in this location. We use HypoDD to obtain relative double-difference locations of the catalog events. Both routine locations and HypoDD relocations of the catalog events suggest a fault dipping west beneath the mountain range, opposite of typical Basin and Range normal faults that dip beneath alluvial valleys. Moment tensors for the largest five events show normal faulting consistent with the west-dipping fault seen in the seismicity. Hydrothermal features in the area, including a geothermal power plant, suggest that fluids may be a contributing factor to the earthquake swarm triggering. We examine the role of fluids by exploring bounds on diffusion parameters and investigating spatial migration characteristics of the swarm seismicity. We conclude that this swarm is the result of heterogenous stress conditions in a prefractured region.

58 GEOSCIENCES↗
NASA NTRS1976

Strainmeters and tiltmeters in geophysics

Several types of sensitive strainmeters and tiltmeters have been developed, and it is now becoming clear which geophysical applications are most suitable for these instruments. In general, strainmeters and tiltmeters are used for observing ground deformation at periods of minutes to days. Small-scale lateral inhomogeneities at the instrument sites distort signals by a few percent, although the effects of large structures can be calculated. In earth tide work these lateral inhomogeneities and unknown ocean loading signals prevent accurate values of the regional tide from being obtained. This limits tidal investigations to looking for temporal variations, possibly associated with pre-earthquake dilatancy, and spatial variations caused by gross elasticity contrasts in the local geological structure. Strainmeters and tiltmeters are well suited for observing long-period seismic waves, seismic slip events on faults and volcano tumescence, where small site-induced distortions in the measured signals are seldom important.

Goulty, N. R.↗
NASA NTRS1983

Two-dimensional Finite Element Modeling for Modeling Tectonic Stress and Strain

Techniques of finite element analysis in two dimensional plane strain were applied to problems of geophysics and tectonics. More specifically, the flexibility of the finite element method was employed to address problems involving geological complexity and fault interactions. The modeling of effective anisotropy in material elastic properties proved useful in describing the deformation of faulted crustal blocks. The applications of this modeling work to problems of actual tectonics in southern California was explored. Preliminary models show encouraging agreement with measured tectonic strain in this region, and modeling work was done to gain an understanding of the stress state in a locked fault region with future seismic potential.

Lyzenga, G. A.↗
NASA NTRS1993

GPS measurements of crustal deformation associated with the 22 April 1991, Valle de la Estrella , Costa Rica earthquake

Results are presented for analysis of Global Positioning System (GPS) measurements made at sites in Costa Rica in February and July 1991. Significant horizontal and vertical displacements relative to February positions were observed. Differences were found in GPS derived vertical and horizontal displacements as compared to other types of geodetic measurements of uplift in the coastal regions. A slip dislocation model which fits the GPS measured displacement was computed. Differences between these data sets and their associated models were investigated in light of unmodeled slip heterogeneity on the fault and post-seismic displacements.

Lundgren, Paul R.↗
DOE OSTI2021

Machine-Learning-Based High-Resolution Earthquake Catalog Reveals How Complex Fault Structures Were Activated during the 2016–2017 Central Italy Sequence

The 2016–2017 central Italy seismic sequence occurred on an 80 km long normal-fault system. The sequence initiated with the Mw 6.0 Amatrice event on 24 August 2016, followed by the Mw 5.9 Visso event on 26 October and the Mw 6.5 Norcia event on 30 October. We analyze continuous data from a dense network of 139 seismic stations to build a high-precision catalog of ~900,000 earthquakes spanning a 1 yr period, based on arrival times derived using a deep-neural-network-based picker. Our catalog contains an order of magnitude more events than the catalog routinely produced by the local earthquake monitoring agency. Aftershock activity reveals the geometry of complex fault structures activated during the earthquake sequence and provides additional insights into the potential factors controlling the development of the largest events. Activated fault structures in the northern and southern regions appear complementary to faults activated during the 1997 Colfiorito and 2009 L’Aquila sequences, suggesting that earthquake triggering primarily occurs on critically stressed faults. Delineated major fault zones are relatively thick compared to estimated earthquake location uncertainties, and a large number of kilometer-long faults and diffuse seismicity were activated during the sequence. These properties might be related to fault age, roughness, and the complexity of inherited structures. The rich details resolvable in this catalog will facilitate continued investigation of this energetic and well-recorded earthquake sequence.

58 GEOSCIENCES↗
DOE OSTI2026

One Earth Energy Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, 2D seismic data was acquired at both sites. This report summarizes the findings from the 2D and 3D seismic interpretation at the One Earth Energy site near Gibson City, Illinois. The seismic data confirms the stratigraphic continuity of the Mt. Simon Arkose Zone storage interval and the Eau Claire confining unit across the project area. The seismic data also indicates that there are faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir within the modeled CO 2 plume (for more detailed information, see Faults and Fractures section of One Earth Energy Class VI Permit applications). However, the seismic data also shows that there are no faults within the modeled CO 2 plume that transect the confining unit Eau Claire Formation. The faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir all tip out in the Lower Mt. Simon Formation and do not reach the overlying Eau Claire confining unit. A small 3D survey acquired around the One Earth Energy #1 characterization well confirms these findings.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
NASA NTRS1983

Focal depths and fault plane solutions of earthquakes under the Tibetan plateau

Seismic data are used to place constraints on the focal depths and fault plane solutions of 16 crustal earthquakes beneath the highest parts of the Tibetan plateau. Synthetic seismograms were computed for large events from 1962 to 1976, and compared with the observed waveforms. All of the events studied are found to have occurred at depths less than 15 km and probably between 5 and 10 km, suggesting that the crust below 10-15 km is essentially aseismic, as in most of the Western United States. Fault-plane solutions show combinations of normal and strike slip faulting with T axes consistently oriented approximately east-west, supporting the inference that the upper crust of Tibet is actively extending in an east-west direction. In addition, solutions for two intermediate-depth events at depths of 90 and 85 km show primarily normal faulting with east-west T axes. Results show that brittle deformation occurs at shallow depths in the crust and in the uppermost mantle, but the lower crust seems to be aseismic.

Molnar, P.↗
NASA NTRS2012

Plate Margin Deformation and Active Tectonics Along the Northern Edge of the Yakutat Terrane in the Saint Elias Orogen, Alaska and Yukon, Canada

The northwest directed motion of the Pacific plate is accompanied by migration and collision of the Yakutat terrane into the cusp of southern Alaska. The nature and magnitude of accretion and translation on upper crustal faults and folds is poorly constrained, however, due to pervasive glaciation. In this study we used high-resolution topography, geodetic imaging, seismic, and geologic data to advance understanding of the transition from strike-slip motion on the Fairweather fault to plate margin deformation on the Bagley fault, which cuts through the upper plate of the collisional suture above the subduction megathrust. The Fairweather fault terminates by oblique-extensional splay faulting within a structural syntaxis, allowing rapid tectonic upwelling of rocks driven by thrust faulting and crustal contraction. Plate motion is partly transferred from the Fairweather to the Bagley fault, which extends 125 km farther west as a dextral shear zone that is partly reactivated by reverse faulting. The Bagley fault dips steeply through the upper plate to intersect the subduction megathrust at depth, forming a narrow fault-bounded crustal sliver in the obliquely convergent plate margin. Since . 20 Ma the Bagley fault has accommodated more than 50 km of dextral shearing and several kilometers of reverse motion along its southern flank during terrane accretion. The fault is considered capable of generating earthquakes because it is linked to faults that generated large historic earthquakes, suitably oriented for reactivation in the contemporary stress field, and locally marked by seismicity. The fault may generate earthquakes of Mw <= 7.5.

Bruhn, Ronald L.↗
NASA NTRS1986

Correlation of data on strain accumulation adjacent to the San Andreas Fault with available models

Theoretical and numerical studies of deformation on strike slip faults were performed and the results applied to geodetic observations performed in the vicinity of the San Andreas Fault in California. The initial efforts were devoted to an extensive series of finite element calculations of the deformation associated with cyclic displacements on a strike-slip fault. Measurements of strain accumulation adjacent to the San Andreas Fault indicate that the zone of strain accumulation extends only a few tens of kilometers away from the fault. There is a concern about the tendency to make geodetic observations along the line to the source. This technique has serious problems for strike slip faults since the vector velocity is also along the fault. Use of a series of stations lying perpendicular to the fault whose positions are measured relative to a reference station are suggested to correct the problem. The complexity of faulting adjacent to the San Andreas Fault indicated that the homogeneous elastic and viscoelastic approach to deformation had serious limitations. These limitation led to the proposal of an approach that assumes a fault is composed of a distribution of asperities and barriers on all scales. Thus, an earthquake on a fault is treated as a failure of a fractal tree. Work continued on the development of a fractal based model for deformation in the western United States. In order to better understand the distribution of seismicity on the San Andreas Fault system a fractal analog was developed. The fractal concept also provides a means of testing whether clustering in time or space is a scale-invariant process.

Turcotte, Donald L.↗
DOE OSTI2020

The Value of Information From Horizontal Distributed Acoustic Sensing Compared to Multicomponent Geophones Via Machine Learning

Abstract Faults play an important role in recharging many geothermal reservoirs, and seismic information can image the locations of these faults. The value of information (VOI) metric is used to objectively quantify and compare the value of two types of seismic receiver data via a machine learning approach. The demonstrated VOI methodology is novel by including spatial models from seismic data and obtaining the information statistics from machine learning. Our two-dimensional numerical experiments compare images created from sparsely spaced (80 m), two-component geophone sampling to high spatial resolution (1 m), single-component DAS. We used a three-fold cross validation of a U-Net convolutional neural networks to achieve average classification statistics. The results suggest that when horizontal sources are utilized, geophones and DAS identify reflectors and non-reflectors at roughly the same rate. The average F1 score for horizontal DAS is 0.939 and 0.931 for geophones. For images created from a vertical source, DAS performed marginally better (F1 = 0.919) than geophones (F1 = 0.877). Our transferrable methodology can provide guidance on which acquisition scenarios can improve images of important structures in the subsurface and present an efficient method for obtaining reliability statistics from high-dimensional, spatial data.

Energy & Fuels↗
NASA NTRS1985

Low-angle normal faulting and isostatic response in the Gulf of Suez: Evidence from seismic interpretation and geometric reconstruction

Tectonic extension within continental crust creates a variety of major features best classed as extensional orogens. These features have come under increasing attention in recent years, with the welding of field observation and theoretical concepts. Most recent advances have come from the Basin and Range Province of the southwestern United States and from the North Sea. Application of these geometric and isostatic concepts, in combination with seismic interpretation, to the southern Gulf of Suez, an active extensional orogen, allows generation of detailed structural maps and geometrically balanced sections which suggest a regional structural model. Geometric models which should prove to be a valuable adjunct to numerical and thermal models for the rifting process are discussed.

Perry, S. K.↗
DOE OSTI2020

Stable Rate of Slip Along the Karakax Section of the Altyn Tagh Fault from Observation of Interglacial and Postglacial Offset Morphology and Surface Dating

Digital elevation maps obtained using TanDEM-X and Pleiades data combined with newly obtained surface age estimates using cosmogenic radionuclide (CRN) and optically simulated luminescence (OSL) methods are used to quantify the slip rate along the western section of the Altyn Tagh Fault in southern Xinjiang. The reconstruction of the conical shape of massive alluvial fans inferred to be from the Eemian (115±7 ka) from CRN dating shows consistent left-lateral offsets of 300±20 m, yielding a slip rate of 2.6±0.3 mm/year. Successive episodes of incision have left cut terraces inset in wide canyons, 10–25 m below the fans' surface. The incision was followed by the deposition of a broad terrace of early Holocene age, which is reincised by modern stream channels. Near the village of Shanxili, a 200 m wide valley is partially dammed by a shutter ridge displaced by the fault. A fill terrace deposited upstream from the ridge has an OSL age of 8.8±0.6 ka. The 23±2 m offset of the riser incising the terrace indicates a minimum postdepositional movement on the fault, yielding a Holocene rate of 2.6±0.5 mm/year, consistent with the 115-ka average slip rate. Furthermore, scarp degradation analysis using mass diffusion reveals a nonlinear relationship between fault displacement and degradation coefficient along the progressively exposed fault scarp, a pattern suggesting either seismic clustering or variable diffusion rate since the Eemian. Together with the Gozha Co-Longmu Co fault to the south, the Karakax section of the Altyn Tagh Fault contributes to the eastward movement of the western corner of Tibet.

58 GEOSCIENCES↗
NASA NTRS1977

Faulting parameters derived from computer simulation of earthquakes

Seismic source parameters, average displacement, rupture length, and strain energy release are investigated by computer simulation using a coupled massive block model of the sliding along an active fault. Average displacements and energy release vary considerably with the degree of heterogeneity in the friction and elastic parameters used in the model. Strain energy release is determined primarily by the product of dynamic friction, rupture length, and average displacement. Interrelationships among the faulting parameters are consistent with theoretical arguments and experimental data. The variation in the frequency of occurrence of simulation events with strain energy release is different from the variation in the frequency of naturally occurring events with seismic energy.

Cohen, S. C.↗
NASA NTRS2008

Rapid Ice Mass Loss: Does It Have an Influence on Earthquake Occurrence in Southern Alaska?

The glaciers of southern Alaska are extensive, and many of them have undergone gigatons of ice wastage on time scales on the order of the seismic cycle. Since the ice loss occurs directly above a shallow main thrust zone associated with subduction of the Pacific-Yakutat plate beneath continental Alaska, the region between the Malaspina and Bering Glaciers is an excellent test site for evaluating the importance of recent ice wastage on earthquake faulting potential. We demonstrate the influence of cumulative glacial mass loss following the 1899 Yakataga earthquake (M=8.1) by using a two dimensional finite element model with a simple representation of ice fluctuations to calculate the incremental stresses and change in the fault stability margin (FSM) along the main thrust zone (MTZ) and on the surface. Along the MTZ, our results indicate a decrease in FSM between 1899 and the 1979 St. Elias earthquake (M=7.4) of 0.2 - 1.2 MPa over an 80 km region between the coast and the 1979 aftershock zone; at the surface, the estimated FSM was larger but more localized to the lower reaches of glacial ablation zones. The ice-induced stresses were large enough, in theory, to promote the occurrence of shallow thrust earthquakes. To empirically test the influence of short-term ice fluctuations on fault stability, we compared the seismic rate from a reference background time period (1988-1992) against other time periods (1993-2006) with variable ice or tectonic change characteristics. We found that the frequency of small tectonic events in the Icy Bay region increased in 2002-2006 relative to the background seismic rate. We hypothesize that this was due to a significant increase in the rate of ice wastage in 2002-2006 instead of the M=7.9, 2002 Denali earthquake, located more than 100km away.

Sauber, Jeanne M.↗
NASA NTRS1984

Geodynamics of crustal deformation and seismotectonic block movements in central Europe

Geological observations reveal the style of neotectonic near-surface stresses and deformations in central Europe. Seismic activity, focal depths and fault plane solutions of earthquakes indicate kinematic reactions within the crust. A crustal deformation model which may account for the Rhine graben systems and the associated seismotectonic block movements in Europe is presented. A computer aided tomography to gravity anomalies is used in determining the crustal stresses in central Europe. Tomographical interpretations of gravity data with respect to seismic stresses are discussed. Kinematics and dynamics are integrated to show that the measured regional stresses in central Europe are derivable from the convection generated traction on the boundary of the elastic spherical shell of the crust as inferred from satellite derived gravity data.

Liu, H. S.↗
DOE OSTI2022

Limited Dynamic Triggering in the Utah Region, USA

The state of Utah, USA, experiences around 3800 catalogued earthquakes per year, highlighting that the region is seismically active and susceptible to earthquakes. Following the 2002 Denali Fault (M7.9) earthquake in Alaska, the region showed an elevated seismicity rate for 3 weeks following the passage of high amplitude surface waves, suggesting that the region may be particularly susceptible to dynamic triggering. With over 23,396 faults and each fault presenting a potential fault for triggering, we systematically search for dynamic triggering throughout the state of Utah caused by large, global earthquakes with M ≥ 7. Specifically, we analyse earthquake catalogues and all available waveform data to determine statistical increases of seismicity rate following the passage of seismic arrivals. While we find instances of dynamic triggering, our results show that these events occur sparsely in the region. In total, less than 20 per cent of the 273 main shocks that occur from 2000 to the end of 2017 show a statistical indication of dynamic triggering throughout the Utah region, highlighting that dynamic triggering is limited for stresses created by transient signals from global M ≥ 7 earthquakes, with the exception being the Denali Fault (M7.9), Alaska earthquake (i.e. an instance of significant triggering).

58 GEOSCIENCES↗
DOE OSTI2022

Permeability-controlled migration of induced seismicity to deeper depths near Venus in North Texas

Abstract Migration of seismic events to deeper depths along basement faults over time has been observed in the wastewater injection sites, which can be correlated spatially and temporally to the propagation or retardation of pressure fronts and corresponding poroelastic response to given operation history. The seismicity rate model has been suggested as a physical indicator for the potential of earthquake nucleation along faults by quantifying poroelastic response to multiple well operations. Our field-scale model indicates that migrating patterns of 2015–2018 seismicity observed near Venus, TX are likely attributed to spatio-temporal evolution of Coulomb stressing rate constrained by the fault permeability. Even after reducing injection volumes since 2015, pore pressure continues to diffuse and steady transfer of elastic energy to the deep fault zone increases stressing rate consistently that can induce more frequent earthquakes at large distance scales. Sensitivity tests with variation in fault permeability show that (1) slow diffusion along a low-permeability fault limits earthquake nucleation near the injection interval or (2) rapid relaxation of pressure buildup within a high-permeability fault, caused by reducing injection volumes, may mitigate the seismic potential promptly.

58 GEOSCIENCES↗