Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2019

Geophysical Monitoring of Moisture-Induced Landslides: A Review

Geophysical monitoring of landslides can provide insights into spatial and temporal variations of subsurface properties associated with slope failure. Recent improvements in equipment, data analysis, and field operations have led to a significant increase in the use of such techniques in monitoring. Geophysical methods complement intrusive approaches, which sample only a very small proportion of the subsurface, and walk-over or remotely sensed data, which principally provide information only at the ground surface. In particular, recent studies show that advances in geophysical instrumentation, data processing, modeling, and interpretation in the context of landslide monitoring are significantly improving the characterization of hillslope hydrology and soil and rock hydrology and strength and their dynamics over time. This review appraises the state of the art of geophysical monitoring, as applied to moisture-induced landslides. Here we focus on technical and practical uses of time-lapse methods in geophysics applied to monitoring moisture-induced landslide. The case studies identified in this review show that several geophysical techniques are currently used in the monitoring of subsurface landslide processes. These geophysical contributions to monitoring and predicting the evolution of landslide processes are currently underrealized. Hence, the further integration of multiple-parametric and geotechnically coupled geophysical monitoring systems has considerable potential. The complementary nature of certain methods to map the distribution of subsurface moisture and elastic moduli will greatly increase the predictive and monitoring capacity of early warning systems in moisture-induced landslide settings.

58 GEOSCIENCES↗
DOE OSTI2020

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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 OSTI2013

Source Characterization and Temporal Variation of Methane Seepage from Thermokarst Lakes on the Alaska North Slope in Response to Arctic Climate Change

The goals of this research were to characterize the source, magnitude and temporal variability of methane seepage from thermokarst lakes (TKL) within the Alaska North Slope gas hydrate province, assess the vulnerability of these areas to ongoing and future arctic climate change and determine if gas hydrate dissociation resulting from permafrost melting is contributing to the current lake emissions. Analyses were focused on four main lake locations referred to in this report: Lake Qalluuraq (referred to as Lake Q) and Lake Teshekpuk (both on Alaska's North Slope) and Lake Killarney and Goldstream Bill Lake (both in Alaska's interior). From analyses of gases coming from lakes in Alaska, we showed that ecological seeps are common in Alaska and they account for a larger source of atmospheric methane today than geologic subcap seeps. Emissions from the geologic source could increase with potential implications for climate warming feedbacks. Our analyses of TKL sites showing gas ebullition were complemented with geophysical surveys, providing important insight about the distribution of shallow gas in the sediments and the lake bottom manifestation of seepage (e.g., pockmarks). In Lake Q, Chirp data were limited in their capacity to image deeper sediments and did not capture the thaw bulb. The failure to capture the thaw bulb at Lake Q may in part be related to the fact that the present day lake is a remnant of an older, larger, and now-partially drained lake. These suggestions are consistent with our analyses of a dated core of sediment from the lake that shows that a wetland has been present at the site of Lake Q since approximately 12,000 thousand years ago. Chemical analyses of the core indicate that the availability of methane at the site has changed during the past and is correlated with past environmental changes (i.e. temperature and hydrology) in the Arctic. Discovery of methane seeps in Lake Teshekpuk in the northernmost part of the lake during 2009 reconnaissance surveys provided a strong impetus to visit this area in 2010. The seismic methods applied in Lake Teshekpuk were able to image pockmarks, widespread shallow gas in the sediments, and the relationship among different sediment packages on the lake's bottom, but even boomer seismics did not detect permafrost beneath the northern part of the lake. By characterizing the biogeochemistry of shallow TKL with methane seeps we showed that the radical seasonal shifts in ice cover and temperature. These seasonal environmental differences result in distinct consumption and production processes of biologically-relevant compounds. The combined effects of temperature, ice-volume and other lithological factors linked to seepage from the lake are manifest in the distribution of sedimentary methane in Lake Q during ice-covered and ice-free conditions. The biogeochemistry results illustrated very active methanotrophy in TKLs. Substantial effort was subsequently made to characterize the nature of methanotrophic communities in TKLs. We applied stable isotope probing approaches to genetically characterize the methanotrophs most active in utilizing methane in TKLs. Our study is the first to identify methane oxidizing organisms active in arctic TKLs, and revealing that type I methanotrophs and type II methanotrophs are abundant and active in assimilating methane in TKLs. These organisms play an important role in limiting the flux of methane from these sites. Our investigations indicate that as temperatures increase in the Arctic, oxidation rates and active methanotrophic populations will also shift. Whether these changes can offset predicted increases in methanogenesis is an important question underlying models of future methane flux and resultant climate change. Overall our findings indicate that TKLs and their ability to act as both source and sink of methane are exceedingly sensitive to environmental change.

03 NATURAL GAS↗
DOE OSTI2018

Interseismic Ground Deformation and Fault Slip Rates in the Greater San Francisco Bay Area From Two Decades of Space Geodetic Data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. We combine interferometric synthetic aperture radar data from the ERS-1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike-slip rates on major locked faults. We use the estimated deep rates to simulate the long-wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long-wavelength signal from the InSAR horizontal mean velocity field, we estimate fault-parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Finally, nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment-compaction processes, respectively.

58 GEOSCIENCES↗
DOE OSTI2017

Multiobjective Optimization of Regional and Teleseismic Data to Constrain the Source of the 12 September 2016 Mw 5.4 Earthquake in South Korea

In this paper, we propose a new approach for source inversion with automatic criteria for the station selection process and using a multiobjective optimization scheme combining teleseismic body-wave and regional surface-wave observations. To illustrate the approach, we determine the source parameters of the 12 September 2016 M w 5.4 South Korea earthquake, as well as its foreshock (same day) and its main aftershock on 19 September 2016. Available teleseismic (short period and broadband) and regional (broadband) seismograms from the Incorporated Research Institutions for Seismology International–Federation of Digital Seismograph Networks (IRIS-FDSN) webservice are selected (see Data and Resources). First, all teleseismic data with signal-to-noise ratio (SNR) above 2 in the 0.3–1.5 Hz frequency band are selected and a cepstrum function is blindly applied to help detect reflected wave arrivals inside a P coda. This method allows the identification of records with weak P-coda and clear depth phase (pP, sP) arrivals and can be used to select an optimal set of teleseismic seismograms to be used for the inversion. Second, regional seismograms are selected according to the fit between their observed and theoretical dispersion curves extracted from local 1D velocity model. Depending on the epicentral distance, two frequency bands are proposed for regional data: 15–50 s for epicentral distances smaller than 2° and 30–80 s otherwise. Finally, we explore the space parameters from exhaustive grid searches for the two objective functions Tf and Rf for teleseismic and regional data, respectively. We show that the seismic moment parameters can be exclusively derived from Rf and can be fixed for Tf. However, the depth parameter is not sensitive to Rf and can be simultaneously inverted from both. The selection of the optimal solution satisfying both cost functions uses the optimum Pareto front that allows a decrease in epistemic uncertainty. Lastly, we discuss the advantage of this approach for improving the characterization of moderate earthquake sources in areas with limited regional instrumentation.

58 GEOSCIENCES↗
NASA NTRS1998

GPS Monitoring of Surface Change During and Following the Fortuitous Occurrence of the M(sub w) = 7.3 Landers Earthquake in our Network

Accomplishments: (1) Continues GPS monitoring of surface change during and following the fortuitous occurrence of the M(sub w) = 7.3 Landers earthquake in our network, in order to characterize earthquake dynamics and accelerated activity of related faults as far as 100's of kilometers along strike. (2) Integrates the geodetic constraints into consistent kinematic descriptions of the deformation field that can in turn be used to characterize the processes that drive geodynamics, including seismic cycle dynamics. In 1991, we installed and occupied a high precision GPS geodetic network to measure transform-related deformation that is partitioned from the Pacific - North America plate boundary northeastward through the Mojave Desert, via the Eastern California shear zone to the Walker Lane. The onset of the M(sub w) = 7.3 June 28, 1992, Landers, California, earthquake sequence within this network poses unique opportunities for continued monitoring of regional surface deformation related to the culmination of a major seismic cycle, characterization of the dynamic behavior of continental lithosphere during the seismic sequence, and post-seismic transient deformation. During the last year, we have reprocessed all three previous epochs for which JPL fiducial free point positioning products available and are queued for the remaining needed products, completed two field campaigns monitoring approx. 20 sites (October 1995 and September 1996), begun modeling by development of a finite element mesh based on network station locations, and developed manuscripts dealing with both the Landers-related transient deformation at the latitude of Lone Pine and the velocity field of the whole experiment. We are currently deploying a 1997 observation campaign (June 1997). We use GPS geodetic studies to characterize deformation in the Mojave Desert region and related structural domains to the north, and geophysical modeling of lithospheric behavior. The modeling is constrained by our existing and continued GPS measurements, which will provide much needed data on far-field strain accumulation across the region and on the deformational response of continental lithosphere during and following a large earthquake, forming the basis for kinematic and dynamic modeling of secular and seismic-cycle deformation. GPS geodesy affords both regional coverage and high precision that uniquely bear on these problems.

Miller, M. Meghan↗
NASA NTRS2010

Observing Traveling Ionospheric Disturbances Caused by Tsunamis Using GPS TEC Measurements

Ground-based Global Positioning System (GPS) measurements of ionospheric Total Electron Content (TEC) show variations consistent with atmospheric internal gravity waves caused by ocean tsunamis following two recent seismic events: the American Samoa earthquake of September 29, 2009, and the Chile earthquake of February 27, 2010. Fluctuations in TEC correlated in time, space, and wave properties with these tsunamis were observed in TEC estimates processed using JPL's Global Ionospheric Mapping Software. These TEC estimates were band-pass filtered to remove ionospheric TEC variations with wavelengths and periods outside the typical range of internal gravity waves caused by tsunamis. Observable variations in TEC appear correlated with the tsunamis in certain locations, but not in others. Where variations are observed, the typical amplitude tends to be on the order of 1% of the background TEC value. Variations with amplitudes ~ 0.1 - 0.2 TECU are observable with periods and timing affiliated with the tsunami. These observations are compared to estimates of expected tsunami-driven TEC variations produced by Embry Riddle Aeronautical University's Spectral Full Wave Model, an atmosphere-ionosphere coupling model, and found to be in good agreement in some locations, though there are cases when the model predicts an observable tsunami-driven signature and none is observed. These TEC variations are not always seen when a tsunami is present, but in these two events the regions where a strong ocean tsunami was observed did coincide with clear TEC observations, while a lack of clear TEC observations coincided with smaller tsunami amplitudes. There exists the potential to apply these detection techniques to real-time GPS TEC data, providing estimates of tsunami speed and amplitude that may be useful for early warning systems.

Galvan, David A.↗
DOE OSTI2023

Economic and operational investigation of CO 2 sequestration through enhanced oil recovery in unconventional reservoirs in Colorado, USA

The ongoing CCUS commercial projects are highly relied on the support of government incentives due to massive capital investment. Here, this study analyzes the economics of carbon capture utilization and sequestration (CCUS) projects, shows a state-wide CCUS deployment exercise, followed by simulation results of enhanced oil recovery (EOR) based CO 2 storage in unconventional reservoirs. The comprehensive economic analysis of capture, transportation, sequestration costs, enhanced 45Q tax credits, and EOR revenue implies the practicality of CO 2 -EOR to offset the high CCUS costs. With the economics understanding, we study the top CO 2 sources, existing CO 2 pipelines, and sequestration sinks in the state of Colorado, USA. This paper next presents results from EOR simulation in one section of the unconventional Denver-Julesburg (DJ) Basin Niobrara and Codell reservoirs. The simulation model is based on a geological static model, incorporated with hydraulic fracture stimulation, history matched to production, and calibrated to the microseismic and time-lapse surface seismic data. The CO 2 -EOR simulation results show that oil production can be increased and more CO 2 stored with: a longer primary production period; the presence of a shut-in period; higher injection rates; and multi-well injectors. The modeling results show that about 7–10 Mscf of CO 2 will be stored when recovering 1 stb of EOR oil. By adding the enhanced oil revenue and the carbon credits together, it is estimated that the most economic case can generate $\$13$ MM when oil price is assumed to be $\$80$/stb, and the EOR oil revenue is 3.4 times greater than that generated from 45Q incentives. It corresponds to the scenario that a five-year primary production is followed by CO 2 injection into four wells with the sequence of injection (4 MMscf/day for 6 months), shut-in (6 months) and production (12 months). The best practices in this study will provide valuable insights for similar CCUS projects in other unconventional fields. Furthermore, this study defines a term named “Carbon Neutrality Index (CNI)” by comparing the amount of CO 2 stored with that burned by EOR oil. The CNI value of 0 indicates the enhanced oil is carbon neutral; a negative CNI value implies there is a net reduction in carbon emission. The 4-year huff-n-puff (HnP) simulation leads to a positive CNI value, indicating that the EOR oil generated in this process is not carbon neutral yet.

03 NATURAL GAS↗
NASA NTRS2020

Neighbors on the Red Planet: Mars Science Laboratory Relay Coordination Post InSight Arrival

The Mars Science Laboratory (MSL) strategic communications planning toolset selects orbiter relay opportunities that the rover will use for transmitting data back to Earth. The toolset was extensively reworked in preparation for the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport's (InSight) arrival at Mars in November 2018 and the regular use of Mars Atmosphere and Volatile Evolution (MAVEN) and Trace Gas Orbiter (TGO) orbiters as relay assets in addition to Mars Reconnaissance Orbiter (MRO) and 2001 Mars Odyssey (ODY). Overflight selection criteria was automated in the new toolset to take into account overflight deconfliction and down-selection as well as tactical timeline planning impacts and total data return. This was done while remaining flexible and configurable for changing mission priorities. As the Curiosity rover ages, the MSL planning team must overcome issues such as reduced budgets, memory bank anomalies, and reduced power availability. These are some examples of factors that affect the strategic communications planning toolset. In addition to adapting to evolving internal mission needs, the toolset must also be flexible to changes in the relay planning interface with other landers. The concept of shared relay or “split passes” was introduced when InSight began operating on Mars just 600 kilometers away and at the same longitude as MSL. This proximity meant that orbiters could now communicate with two landers during the same relay session. Over time, this new operational use-case became common practice as InSight and MSL settled into relay planning negotiations. Today, many TGO relay sessions are shared between MSL and InSight. As more orbiters pursue integrating this capability into nominal relay operations, and more landers arrive on Mars, the strategic communications process and toolsets increases in complexity. Development of a single tool that schedules relay sessions for all orbiters and landers simultaneously may become necessary as the number of Martian spacecraft increases. Examples of how toolset selection criteria and capabilities have helped or hindered MSL planning will be presented. Ongoing improvements to MSL toolsets and processes, as well as to shared relay tools such as Mars Relay Operations Service (MaROS) and General Telecom Predictor (GTP) are also discussed in this paper. Examples of how toolset selection criteria and capabilities have helped or hindered MSL planning will be presented. Ongoing improvements to MSL toolsets and processes, as well as to shared relay tools such as Mars Relay Operations Service (MaROS) and General Telecom Predictor (GTP) are also discussed in this paper.

Vandermey, Nancy↗
DOE OSTI2021

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗
DOE OSTI2019

High frequency (4th order) sequence stratigraphy of Early Miocene deltaic shorelines, offshore Texas and Louisiana

Well logs and 3D seismic reflection data are integrated to image Lower Miocene depocenters of the upper Texas and westernmost Louisiana coastal and offshore areas. Although previously interpreted at a large scale (3rd order cycles), the detailed stratigraphy and depositional history of the early Miocene succession has not been fully developed. The upper lower Miocene interval from Robulus L (MFS10) to Amphistegina B (MFS9) has been further subdivided into five 4th order cycles (each 200 to 300 m thick) to provide finer scale stratigraphic detail and interpretation of depositional environments. Results of the finer scale correlations are presented in a series of sandstone percent maps. Overall, the maps display a strike-elongate (sub-parallel to the present day coastline) orientation for the sandstone bodies which thin to the southeast (basinward). Net sandstone thickness is relatively low in older cycles and increases in younger cycles suggesting an overall progradation of the deltaic system into the mapped area through time. Continued marine regression and southward movement of the shoreline are also confirmed by seismic amplitude 2 maps. After a major transgression associated with Robulus L deltaic progradation occurred under rising sea level conditions. Deltaic sediments prograded southward and sandstone brought to the upper Texas coast was distributed laterally by longshore currents to form strandplain coast lines. Maximum regression occurred during the interval MFS9_2 to MFS9_3 when a deltaic depocenter formed offshore Texas in east High Island area. However, sandstone thickness progressively decreased laterally to the east (offshore Louisiana) where more shale was present and marine processes reworked deltaic derived sandstones into shore parallel bars. There was a clear retreat of the shoreline during the youngest cycle marking the beginning of transgression associated with Amphistegina B . Detailed stratigraphic interpretation (at a 4th order scale) shows more variability in the dominant shoreline processes during early Miocene than previously thought. Recognition of sedimentary systems variability at a 4th order scale is critical for improved hydrocarbon exploration and in understanding potential future CO 2 storage in the area.

58 GEOSCIENCES↗
DOE OSTI2015

An explosion aftershock model with application to on-site inspection

An estimate of aftershock activity due to a theoretical underground nuclear explosion is produced using an aftershock rate model. The model is developed with data from the Nevada National Security Site, formerly known as the Nevada Test Site, and the Semipalatinsk Test Site, which we take to represent soft-rock and hard-rock testing environments, respectively. Estimates of expected magnitude and number of aftershocks are calculated using the models for different testing and inspection scenarios. These estimates can help inform the Seismic Aftershock Monitoring System (SAMS) deployment in a potential Comprehensive Test Ban Treaty On-Site Inspection (OSI), by giving the OSI team a probabilistic assessment of potential aftershocks in the Inspection Area (IA). The aftershock assessment, combined with an estimate of the background seismicity in the IA and an empirically derived map of threshold magnitude for the SAMS network, could aid the OSI team in reporting. Here, we apply the hard-rock model to a M5 event and combine it with the very sensitive detection threshold for OSI sensors to show that tens of events per day are expected up to a month after an explosion measured several kilometers away.

58 GEOSCIENCES↗
DOE OSTI2025

Southwest Regional Partnership on Carbon Sequestration: Phase III (Final Scientific/Technical Report)

The Southwest Regional Partnership on Carbon Sequestration (SWP) is one of 7 regional partnerships formed in 2003 under the U.S. Department of Energy’s (DOE) Regional Carbon Sequestration Partnerships (RCSPs) initiative. The overall purpose of the initiative was to help determine and implement the technology, infrastructure, and regulations most appropriate to promote carbon storage in different regions of the country. Covering Arizona, Colorado, New Mexico, Oklahoma, Utah, and parts of Texas, Wyoming, and Kansas, the SWP evaluated regional carbon storage and utilization potential and focused on technologies and sites that could complement the region’s strong position in energy production. The project progressed through three phases: • Phase I (2003–2005): Characterized regional geologic formations and CO 2 sources, assessed sequestration potential, and identified pilot test sites. • Phase II (2005–2013): Conducted small-scale field tests to validate sequestration methods, including geologic and terrestrial projects. • Phase III (2008–2022): Demonstrated large-scale CO 2 injection at a commercial oil field to test monitoring, verification, and long-term storage strategies. This report covers Phase III. The final project site, the Farnsworth Unit (FWU) in Texas, provided real-world testing of reservoir characterization, monitoring, and risk evaluation tools and processes that could be used in any commercial scale carbon capture, utilization, and storage (CCUS) project. Extensive data collection and analysis helped refine best practices for reservoir characterization, injection monitoring, and storage verification. The SWP contributed to national databases, DOE best practice manuals, and regional geological assessments to support future sequestration efforts. Key lessons learned include the importance of robust data management, strategic site selection, regulatory navigation, and effective industry collaboration. The project’s findings will inform ongoing and future carbon storage initiatives. Task 1 (Regional Characterization) • The SWP continued to participate in national outreach efforts and NATCARB. • The SWP evaluated multiple potential sites before selecting the FWU as the primary field test location. Task 2 (Public Outreach and Education) • The SWP contributed to national databases, DOE best practice manuals, and regional geological assessments to support future sequestration efforts. Task 3 (Permitting and Regulatory Compliance) • The SWP ensured compliance with federal and state regulations, including National Environmental Policy Act (NEPA) requirements. • The SWP obtained all necessary permits for drilling, injection, and monitoring activities. Task 4 (Site Characterization and Planning) • The SWP developed work plans for four key activities: characterization, simulation, monitoring and verification, and risk evaluation. • The SWP collected and synthesized legacy data from multiple sources to build initial static geological models and dynamic reservoir models demonstrating project feasibility. • The SWP conducted an initial risk evaluation and developed mitigation plans. Task 5 (Field Operations and Data Collection) • The SWP drilled, logged, and cored three characterization wells to gather critical subsurface data. • The SWP conducted multiple geophysical surveys, including 3D seismic, crosswell seismic, and vertical seismic profiling, to improve reservoir characterization. Task 6 (Monitoring and Verification) • The SWP performed extensive geological characterization using data from characterization wells and seismic surveys. • The SWP established a surface monitoring network to track CO 2 flux in soil gas, groundwater chemistry, and near-surface atmospheric CO 2 levels. • The SWP built and refined reservoir models to study the effects of relative permeability on simulation behavior and improve calibration with experimental data. Task 7 (Risk Assessment and Model Refinement) • The SWP conducted multiple studies to evaluate reservoir integrity, predict CO 2 plume behavior and improve predictive modeling capabilities. • The SWP refined geological models and used them to enhance the accuracy of simulation models. • The SWP continued quantitative risk assessment of top-ranked risks and strengthened the link between qualitative and quantitative risk methodologies.

02 PETROLEUM↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2021

Oregon Rivers Characterization Tool

DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.

58 GEOSCIENCES↗
DOE OSTI2016

Development of a CO 2 Chemical Sensor for Downhole CO 2 Monitoring in Carbon Sequestration

Geologic storage of carbon dioxide (CO 2 ) has been proposed as a viable means for reducing anthropogenic CO 2 emissions. The means for geological sequestration of CO 2 is injection of supercritical CO 2 underground, which requires the CO 2 to remain either supercritical, or in solution in the water/brine present in the underground formation. However, there are aspects of geologic sequestration that need further study, particularly in regards to safety. To date, none of the geologic sequestration locations have been tested for storage integrity under the changing stress conditions that apply to the sequestration of very large amounts of CO 2 . Establishing environmental safety and addressing public concerns require widespread monitoring of the process in the deep subsurface. In addition, studies of subsurface carbon sequestration such as flow simulations, models of underground reactions and transports require a comprehensive monitoring process to accurately characterize and understand the storage process. Real-time information about underground CO 2 movement and concentration change is highly helpful for: (1) better understanding the uncertainties present in CO 2 geologic storage; (2) improvement of simulation models; and (3) evaluation of the feasibility of geologic CO 2 storage. Current methods to monitor underground CO 2 storage include seismic, geoelectric, isotope and tracer methods, and fluid sampling analysis. However, these methods commonly resulted low resolution, high cost, and the inability to monitor continuously over the long time scales of the CO 2 storage process. A preferred way of monitoring in-situ underground CO 2 migration is to continuous measure CO 2 concentration change in brine during the carbon storage process. An approach to obtain the real time information on CO 2 concentration change in formation solution is highly demanded in carbon storage to understand the CO 2 migration subsurface and to answer the public safety problem. The objective of the study is to develop a downhole CO 2 sensor that can in-situ, continuously monitor CO 2 concentration change in deep saline. The sensor is a Severinghaus-type CO 2 sensor with small size, which renders it can be embedded in monitoring well casing or integrated with pressure/temperature transducers, enabling the development of “smart” wells. The studies included: (1) prepare and characterize metal-oxide electrodes. Test the electrodes response to pH change. Investigate different ions and brine concentration effects on the electrode’s performance. Study the stability of the electrode in brine solution; (2) fabricate a downhole CO 2 sensor with the metal-oxide electrodes prepared in the laboratory. Test the performance of the CO 2 sensor in brine solutions. Study high pressure effects on the performance of the sensor; (3) design and conduct CO 2 /brine coreflooding experiments with the CO2 sensor. Monitor CO 2 movement along the core and test the performance of the sensor in coreflooding tests. Develop a data acquisition system that can digitize the sensor’s output voltage. Our completed research has resulted in deep understanding of downhole CO 2 sensor development and CO 2 monitoring in CO 2 storage process. The developed downhole CO 2 sensor included a metal-oxide electrode, a gas-permeable membrane, a porous steel cup, and a bicarbonate-based internal electrolyte solution. Iridium oxide-based electrode was prepared and used for preparation the CO 2 sensor. The prepared iridium oxide-based electrode displayed a linearly response to pH change. Different factors such as different ions and ions concentration, temperature, and pressure effects on the electrode performance on pH response were investigated. The results indicated that the electrode exhibited a good performance even in high salt concentration of produced water. To improve the electrode performance under high pressure, IrO 2 nanoparticles with the particle size in the range of 1-2 nm were prepared and electrodeposited on stainless steel substrate by cyclic voltammetry. It was observed that the thin film of iridium oxide was formed on the substrate surface and such iridium oxide-based electrode displayed excellent performance under high pressure for longer term. A downhole CO 2 sensor with the iridium oxide-based electrode was prepared. The working principle of the CO 2 sensor is based on the measurement of the pH change of the internal electrolyte solution caused by the hydrolysis of CO 2 and then determination of the CO 2 concentration in water. The prepared downhole CO 2 sensor had the size of diameter of 0.7 in. and length of 1.5 in. The sensor was tested under the pressures of 500 psi, 2,000 psi, and 3,000 psi. A linear correlation was observed between the sensor potential change and dissolved CO 2 concentration in water. The response time of the CO 2 sensor was in the range of 60-100 minutes. Further tests indicated that the CO 2 sensor exhibited good reproducibility under high pressure. A CO 2 /brine coreflooding system was constructed to simulate the real-world CO 2 storage process. The prepared downhole CO 2 sensor was loaded in the system to monitor CO 2 movement during CO 2 /brine coreflooding test. The results indicated that the sensor could detect CO 2 movement in the tests. Further studies showed that the sensor could be recovered by brine flooding after CO 2 /brine flushed the core. The results of the coreflooding tests demonstrated that the sensor had potential application for CO 2 monitoring in carbon sequestration. A data acquisition system for the downhoe CO 2 sensor was developed and coded. The system converted the sensor output signal into digital data and transported the data from downhole to wellhead surface. The data acquisition system was tested and evaluated in the laboratory with the prepared sensor for data collection.

36 MATERIALS SCIENCE↗