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 a presentation for the 2024 Signal Processing Applied to Nonproliferation (SPAN) 2024 workshop. It discusses seismic research conducted at the Redmond Mine as part of the NA-22 STILGAR project.
The subject of neotectonics, covering the structures and structural activity of the last 5 million years (i.e., post-Miocene) is a well-recognized field, including "active tectonics," focussed on the last 500,000 years in a 1986 National Research Council report of that title. However, there is a cartographic gap between tectonic maps, generally showing all features regardless of age, and maps of current seismic or volcanic activity. We have compiled a map intended to bridge this gap, using modern data bases and computer-aided cartographic techniques. The maps presented here are conceptually descended from an earlier map showing tectonic and volcanic activity of the last one million years. Drawn by hand with the National Geographic Society's 1975 "The Physical World" map as a base, the 1981 map in various revisions has been widely reproduced in textbooks and various technical publications. However, two decades of progress call for a completely new map that can take advantage of new knowledge and cartographic techniques. The digital tectonic activity map (DTM), presented in shaded relief (Fig. 1) and schematic (Fig. 2) versions, is the result. The DTM is intended to show tectonism and volcanism of the last one million years, a period long enough to be representative of global activity, but short enough that features such as fault scarps and volcanos are still geomorphically recognizable. Data Sources and Cartographic Methods The DTM is based on a wide range of sources, summarized in Table 1. The most important is the digital elevation model, used to construct a shaded relief map. The bathymetry is largely from satellite altimetry, specifically the marine gravity compilations by Smith and Sandwell (1996). The shaded relief map was designed to match the new National Geographic Society world physical map (1992), although drawn independently, from the digital elevation model. The Robinson Projection is used instead of the earlier Van der Grinten one. Although neither conformal nor equal-area, the Robinson Projection provides a reasonable compromise and retains useful detail at high latitudes.
Connecting the geophysical and geochemical signals recorded at and above the surface of volcanoes with source mechanisms is fundamental to understand transitions from quiescence to eruption, and to integrate the behavior of volcanoes with their regional seismotectonic and hydrological context. Some of these signals include shallow volcanic tremor, a long-lasting (minutes-to-years) ground vibration detected in volcanic areas during unrest; volcanic outgassing, which shows intriguing periodicities over multiple scales, as revealed by space-borne instruments and by ground-based, high-frequency (~1 Hz), sampling methods; inflation and deflation of volcanic edifices, commonly associated with subsurface volume changes; and large-scale (from a few to tens of km2), low-temperature, thermal anomalies, a newly-discovered signal that has been observed to emerge on volcanic flanks from months-to-years prior to gas- and magma-driven eruptions. Traditionally, these geophysical and geochemical signals have been modelled independently, which limits our interpretation of subsurface processes and thus our assessment of unrest. In this work, we show through lumped-parameter models, numerical simulations, and preliminary laboratory experiments that many of the signals recorded around volcanoes can be integrated into a common “stop-and-go” mechanism. In particular, we found that they can arise spontaneously when taking into account the feedbacks between the temporary accumulation of gas and/or magma in the crust (“stop”), and the diffusive transfer of gas and heat towards the surface (“go”). For example, we find that the “stop-and-go” mechanism can explain: (i) the emergence of monochromatic, broadband, and harmonic tremor; (ii) the periodic components commonly observed in volcanic outgassing time series; (iii) the link between large-scale thermal anomalies and shallow hydrothermal systems; (iv) the reported time lags between large-scale thermal anomalies and deformation; and (v) the breathing of potentially destructive volcanic calderas. The “stop-and-go” mechanism provides a consistent and realistic framework to link multiple geophysical and geochemical signals with the processes leading to volcanic unrest and eruption.
Like all operations in which fluids are injected into the subsurface, geologic carbon storage (GCS) presents inherit risks that must be assessed and mitigated to successfully deploy such technology in the field. For example, alterations to pore pressure and the stress state due to fluid injection may trigger certain geohazards, such as potential fault activation and induced seismicity or unintentional hydraulic fracturing. Here we present capabilities for assessing these potential risks using the State of Stress Analysis Tool (SOSAT): a Python library and web-based application capable of probabilistically estimating the subsurface state of stress informed by various field observations. Using posterior distributions of principal stress components, SOSAT can assess the risk of fault activation within the injection formation for either a critically oriented fault or a fault with a user-prescribed orientation. Additionally, SOSAT can estimate the probability of unintentional hydraulic fracturing of the intact reservoir rock. Here, we demonstrate these SOSAT capabilities using a hypothetical GCS site. This type of risk assessment can assist with the effective deployment of GCS technology in the field by informing safe design and management practices.
The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.
Satellite magnetometer observations of the Greenland-Iceland region compare quite well with lower altitude data. The satellite magnetic data suggest magnetically enhanced crust was emplaced by the Iceland Plume. Crustal thicknesses, which may be more than 30 km for the Greenland-Scotland Ridge, were obtained from inversion of the compensating terrain gravity effects that were estimated by spectral correlation analysis of the free-air gravity anomalies and terrain gravity effects. Regional magnetic anomaly maxima overlie possible thickened crust from eastern Iceland to the Greenland Coast. The Iceland-Faroe Ridge may involve thinner crust than the Greenland-Iceland portion of the Greenland-Scotland Ridge. The gravity derived crustal model exceeds a 0.7 correlation with available seismic estimates. In thermally active areas our gravity Moho estimates are systematically deeper than the seismic estimates suggesting local density reductions of the underlying lower crust/upper mantle. In south central Greenland, on the other hand, the gravity Moho estimates are shallower than seismic estimates to suggest a local enhancement of the lower crust/upper mantle density. The dichotomous crust of the Greenland-Iceland and Iceland-Faroe Ridges suggests unequal crustal development by the Iceland Plume and the Mid-Atlantic Ridge, where more crustal material may have been contributed to the North Atlantic Plate than the Eurasian Plate. A new thermal modeling scheme based on Poisson's relation between point pole gravity and thermal potentials allows estimation of magnetic crustal thicknesses. Subsequent magnetic anomaly inversion for susceptibility contrasts infers crustal development of the Greenland-Scotland Ridge by temporally variable pulses in plume strength.
Subsurface environments host diverse microorganisms in fluid-filled fractures; however, little is known about how geological and hydrological processes shape the subterranean biosphere. Here, we sampled three flowing boreholes weekly for 10 mo in a 1478-m-deep fractured rock aquifer to study the role of fracture activity (defined as seismically or aseismically induced fracture aperture change) and advection on fluid-associated microbial community composition. We found that despite a largely stable deep-subsurface fluid microbiome, drastic community-level shifts occurred after events signifying physical changes in the permeable fracture network. The community-level shifts include the emergence of microbial families from undetected to over 50% relative abundance, as well as the replacement of the community in one borehole by the earlier community from a different borehole. Null-model analysis indicates that the observed spatial and temporal community turnover was primarily driven by stochastic processes (as opposed to deterministic processes). We, therefore, conclude that the observed community-level shifts resulted from the physical transport of distinct microbial communities from other fracture(s) that outpaced environmental selection. Given that geological activity is a major cause of fracture activity and that geological activity is ubiquitous across space and time on Earth, our findings suggest that advection induced by geological activity is a general mechanism shaping the microbial biogeography and diversity in deep-subsurface habitats across the globe.
On the six Apollo landed missions, the Astronauts deployed the Apollo Lunar Surface Experiments Package (ALSEP) science stations which measured active and passive seismic events, magnetic fields, charged particles, solar wind, heat flow, the diffuse atmosphere, meteorites and their ejecta, lunar dust, etc. Today's scientists are able to extract new information and make new discoveries from the old ALSEP data utilizing recent advances in computer capabilities and new analysis techniques. However, current-day investigators are encountering problems trying to use the ALSEP data. In 2007 archivists from NASA Goddard Space Flight Center (GSFC) National Space Science Data Center (NSSDC) estimated only about 50 percent of the processed ALSEP lunar surface data-of-interest to current lunar science investigators were in the NSSDC archives. The current-day lunar science investigators found most of the ALSEP data, then in the NSSDC archives. were extremely difficult to use. The data were in forms often not well described in the published reports and rerecording anomalies existed in the data which could only be resolved by tape experts. To resolve this problem, the DPS Lunar Data Node was established in 2008 at NSSDC and is in the process of successfully making the existing archived ALSEP data available to current-day investigators in easily useable forms. In July of 2010 the NASA Lunar Science Institute (NLSI) at Ames Research Center established the Recovery of Missing ALSEP Data Focus Group in recognition of the importance of the current activities to find the raw and processed ALSEP data missing from the NSSDC archives.
What: Thermal moonquakes are small amplitude events that are produced by diurnal temperature changes. Why: Finding the locations of thermal moonquakes will lead to information about lunar surface processes. Where: Apollo 17 Lunar Seismic Profiling Experiment (LSPE) is able to locate thermal moonquake unlike Apollo 12-16. The primary science goal was an active source experiment to study the detailed structure of the lunar crust using 8 explosive packages (EPs). The secondary science goal was to passively listen for lunar seismic activity.
Manipulation of subsurface stresses through fluid injection/extraction may induce seismicity through modulation of overpressure along faults. Globally, most industrial operations do not induce significant earthquakes. However, the possibility of inducing moderate to large earthquakes is recognized as a significant risk for basin-scale carbon dioxide (CO 2 ) storage operations (Niass et al., 2017). One possible method to reduce the occurrence of induced seismicity is to maintain only small (or zero) fluid overpressure changes in the subsurface through extraction of formation brine. For active pressure management (APM), operators may employ different procedures including extraction of various brine volumes and different configurations of injection/extraction well locations. However, the impact each these operational procedures on induced seismicity must be understood in order to develop management strategies. Here, we describe a numerical modeling workflow to evaluate the efficacy of APM techniques, that when applied to a specific site may be used to inform operational management strategies. We couple the multiphase reservoir model, NUFT, with the earthquake simulator, RSQSim and simulate seismicity induced by CO 2 injection into a hypothetical brine saturated reservoir. Reservoir overpressure is modified by co-extraction of brine at various fractions of the injected CO 2 volume from five potential wells. Our goal is to understand how the well-pair configuration and the extracted brine volume impact the total number and the maximum event magnitude. Additionally, the performance of each scenario is also evaluated in terms of the volume of CO 2 that can be stored before an earthquake above a tolerance threshold is induced.
Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.
Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.
Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.
Abstract The large‐scale implementation of geological carbon sequestration has raised concerns about potential fault activation and induced seismicity, which could compromise storage integrity and pose seismic risks. We theoretically compared two failure criteria, Mohr‐Coulomb (MC) and Modified Cam‐Clay criteria (MCC), to assess fault rupture during CO 2 storage. Both criteria characterize fault behavior in a specific stress regime but differ in reducing the complexity of fault rupture to a few key mechanisms. Using a coupled hydromechanical model, we demonstrate that the choice of a failure criterion and the physics of fault weakening associated with these criteria strongly condition fault response to a given fluid injection. MC mainly relates rupture to friction, while MCC relates rupture to fault poro‐plasticity. Our findings highlight that the selection of a failure criterion, being inherently subjective, can significantly alter the predicted fault behavior during CO 2 storage, thereby impacting the reliability of geomechanical risk assessments.
A compilation of maps of global geophysical and geological data plotted on a common scale and projection is presented. The maps include satellite gravity, magnetic, seismic, volcanic, tectonic activity, and mantle velocity anomaly data. The Bibliographic references for all maps are included.
A fault investigation was performed at the Mars and Uranus antenna sites at the Goldstone Deep Space Communications Complex in the Mojave desert. The Mars/Uranus Station consists of two large-diameter reflector antennas used for communication and control of deep-space probes and other missions. The investigation included interpretation of Landsat thematic mapper scenes, side-looking airborne radar transparencies, and both color-infrared and black-and-white aerial photography. Four photolineaments suggestive of previously undocumented faults were identified. Three generally discrete morphostratigraphic alluvial-fan deposits were also recognized and dated using geomorphic and soil stratigraphic techniques. Fourteen trenches were excavated across the four lineaments; the trenches show that three of the photolineaments coincide with faults. The last displacement of two of the faults occurred between about 12,000 and 35,000 years ago. The third fault was judged to be older than 12,000 years before present (ybp), although uncertainty remains. None of the surface traces of the three faults crosses under existing antennas or structures; however, their potential activity necessitates appropriate seismic retrofit designs and loss-prevention measures to mitigate potential earthquake damage to facilities and structures.
We analyze over a century of geodetic data to study crustal deformation and plate motion around the Japanese Islands, using the block-fault model for crustal deformation developed by Matsu'ura et al. (1986). We model the area including the Japanese Islands with 19 crustal blocks and 104 faults based on the distribution of active faults and seismicity. Geodetic data are used to obtain block motions and average slip rates of faults. This geodetic model predicts that the Pacific plate moves N deg 69 +/- 2 deg W at about 80 +/- 3 mm/yr relative to the Eurasian plate which is much lower than that predicted in geologic models. Substantial aseismic slip occurs on the subduction boundaries. The block containing the Izu Peninsula may be separated from the rigid part of the Philippine Sea plate. The faults on the coast of Japan Sea and the western part of the Median Tectonic Line have slip rates exceeding 4 mm/yr, while the Fossa Magna does not play an important role in the tectonics of the central Japan. The geodetic model requires the division of northeastern Japan, contrary to the hypothesis that northeastern Japan is a part of the North American plate. Owing to rapid convergence, the seismic risk in the Nankai trough may be larger than that of the Tokai gap.
Radiometric (K-Ar and Ar-40/Ar-39) age determinations of volcanic and volcaniclastic rocks, combined with structural, gravity, and seismic reflection data, are used to constrain the age of sedimentary strata contained within the seismically and volcanically active northern Malawi (Nyasa) rift and to characterize changes in basin and flank morphologies with time. Faulting and volcanism within the Tukuyu-Karonga basin began at approximately 8.6 Ma, when sediments were deposited in abroad, initially asymmetric lake basin bounded on its northeastern side by a border fault system with minor topographic relief. Extensions, primarily by a slip along the border fault, and subsequent regional isostatic compensation led to the development of a 5-km-deep basin bounded by broad uplifted flanks. Along the low-relief basin margin opposite border fault, younger stratigraphic sequences commonly onlap older wedge-shaped sequences, although their internal geometry is often progradational. Intrabasinal faulting, flankuplift, and basaltic and felsic volcanism from centers at the northern end of the basin became more important at about 2.5 Ma when cross-rift transfer faults developed to link the Tukuyu-Karonga basin to the Rukwa basin. Local uplift and volcanic construction at the northern end of the basin led to a southeastward shift in the basin's depocenter. Sequence boundaries are commonly erosional along this low-relief (hanging wall) margin and conformable in the deep lake basin. The geometry of stratigraphic sequences and the distribution of the erosion indicate that horizontal and vertical crustal movements both across and along the length of the rift basin led to changes in levels of the lake, irrespective of paleoclimatic fluctuations.