Search NASA⌕ Search

SEARCH · Search NASA

Results for “CODAS”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗
DOE OSTI2024

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗
DOE OSTI2025

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Probabilistic Seismic Hazard Analysis for Iraq Based on the Updated Earthquake Catalog (1900-2021) and Ground Motion Characteristics

Onur et al. (2017) compiled the first comprehensive earthquake catalog for Iraq, covering 1900 to 2009 within 26°–40°N latitude and 36°–51°E longitude. This catalog was utilized in a probabilistic seismic hazard assessment (PSHA) by Abdulnaby et al. (2020) to aid in updating Iraq’s building code seismic provisions. Recently, we have updated the earthquake catalog for Iraq by adding earthquakes recorded from 2010 to 2021 and directly calculating moment magnitude (Mw) for about 2,800 earthquakes using the coda envelope methodology and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq.

58 GEOSCIENCES↗
NASA NTRS1977

Elastic wave propagation in a highly scattering medium - A diffusion approach

The propagation of elastic waves in the moon, where the first seismograms were characterized by the presence of a long coda attributed to strongly scattered waves, is modeled with the aid of the time-dependent equation of radiative transfer. The average energy density as a function of time and space is described by the diffusion equation with linear dissipation on the assumption that all the energy present has been scattered many times and the time and distance scales of the problem are long compared to the scales of the scattering process. Ultrasonic experiments in the laboratory confirm the applicability of the formalism.

Dainty, A. M.↗
NASA NTRS1981

Structure of the lunar crust at highland site Apollo Station 16

The seismic crustal structure of the moon is known in the region of Oceanus Procellarum from the analysis of artificial impact data. To extend this knowledge we have used data from natural lunar seismic events to search for secondary seismic wave arrivals in the form of peg-leg multiples caused by reflections at crustal interfaces and converted arrivals caused by refractions at crustal interfaces. A polarization filter has been applied to the data to enhance the rectilinear particle motion expected for the onset of these body wave arrivals in the scattered coda of lunar seismograms. The results of this work tentatively indicate that the highland site at station 16 has a 75 km thick crust with an intermediate 20 km interface, compared to the 60 km crust and 20 km interface in Oceanus Procellarum, a mare region. Since the 20 km upper crust appears to exist at both highland and mare sites, it probably does not represent a mare basalt layer but rather a more general feature of the crust. Crustal thickness may partially control elevation by isostasy.

Goins, N. R.↗
NASA NTRS1993

Using Transparent Informed Prefetching (TIP) to reduce file read latency

As processor performance gains continue to outstrip Input/Output gains, I/O performance is becoming critical to overall system performance. File read latency is the most significant bottleneck for high performance I/O. Other aspects of I/O performance benefit from recent advances in disk bandwidth and throughput resulting from disk arrays, and in write performance derived from buffered write behind and the Log-structured File System. The access gap problem limiting improvements in read latency is exacerbated by distributed file systems operating over networks with diverse bandwidth. Focus is on extending the power of caching and prefetching to reduce file read latencies by exploiting hints from high-levels of a system. Such Transparent Informed Prefetching, TIP, and its benefits are described. It is argued that hints that disclose high level knowledge are a means for transferring optimization information across, without violating, module boundaries. How TIP can be used to convert the high throughput of new technologies such as disk arrays and log-structured file systems into low latency for applications is discussed. Our preliminary experiments show reductions in wall - clock execution time of 13 percent and 20 percent for a multiple module compilation tool (make) accessing data on a local disk and remote Coda file server, respectively, and a reduction of 30 percent for a text search (grep) remotely accessing many small files.

Patterson, R. H.↗
NASA NTRS2001

Recent Developments in DAO's Finite-Volume Data Assimilation System

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. The data assimilation method implemented in CODAS incorporates a simplified version of the model bias estimation and correction algorithm, as described by Dee and da Silva (1998). In this talk we will briefly describe the general system formulation, and focus on the impact of 3 data types recently introduced, namely: 1) cloud tracks winds from the Multi-angle Imaging Spectrometer by the US Air Force, and 3) temperature and moisture information derived from GPS refractivity occultation measurements. The impact of these data types on observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills will be discussed. In addition we will assess the impact of cloud assimilation on top of the atmosphere radiation fields estimated from CERES measurements.

daSilva, Arlindo↗
NASA NTRS2012

Recovery of Deep Moonquake Focal Mechanisms

Deep moonquakes are clustered not only in space but also in time: their recurrence times correspond to the durations of the anomalistic and draconic months, with some clusters preferring one of the two periods, while others are active with both periods. A key constraint for the understanding of the connection between the orbital motion of the Moon and its seismic activity is the focal mechanism: the orientation of the fault surface on which failure occurs during the quake. Due to the small aperture of the Apollo seismic network and the strong scattering of seismic waves within the lunar crust, the evaluation of P wave first motions to constrain the strike and dip of the fault planes is not feasible. Instead we evaluate the amplitude ratios of P and S waves. Seismograms are rotated into the P-SV-SH coordinate frame and amplitudes are determined as averages over short time windows after the arrival to reduce the impact of the scattering coda, which is independent of the source orientation. We allow for reversals of the fault motion, as observed for some clusters in previous studies, by taking into account the absolute amplitude only, without sign. An empirical site correction factor is applied to correct for amplitude distortions in the crust. We construct ensembles of fault plane solutions using an exhaustive grid search by accepting all orientations that reproduce the measured amplitude ratios within the observed standard deviations. Since all events of a given cluster are supposed to share the same fault plane, the combination of the individual inversion results further constrains the orientation. We evaluate 106 events from 25 different moonquake clusters. The most active cluster A001 contributes 37 events, while others contribute 1 to 9 events per cluster. Comparison of fault orientations with the variation of the tidal stress results in preferred orientations.

Weber, Renee C.↗
NASA NTRSDate not supplied

Near real-time air quality forecasts using the NASA GEOS model

NASA's Global Modeling and Assimilation Office (GMAO) produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time5-day forecasts of atmospheric composition at unprecedented horizontal resolution of 0.25 degrees (~25 km). This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module (version 12) to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). Satellite observations are assimilated into the system for improved representation of weather and smoke. The assimilation system is being expanded to include chemically reactive trace gases. We discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions, notably incorporating new observations (TROPOMI, geostationary satellites) and machine learning techniques. We show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.

Air Quality Forecast↗
NASA NTRSDate not supplied

NASA's High-Resolution GEOS Forecasting and Reanalysis Products: A Unified Tool from Local to Global Scales

NASA's GMAO produces high-resolution global forecasts for weather, aerosols, and air quality. The NASA Global Earth Observing System (GEOS) model has been expanded to provide global near-real-time 5-day forecasts of atmospheric chemical composition at unprecedented horizontal resolution of 0.25 degrees (~25 km), freely available to the public. This composition forecast system (GEOS-CF) combines the operational GEOS weather forecasting model with the state-of-the-science GEOS-Chem chemistry module to provide detailed analysis of a wide range of air pollutants such as ozone, carbon monoxide, nitrogen oxides, and fine particulate matter (PM2.5). GEOS-CF also assimilated satellite observations into the system for improved representation of weather and smoke. The assimilation system is currently being expanded to include chemically reactive trace gases. While the main focus of this new product is on tropospheric air quality information, the GEOS-Chem chemistry model used in this system includes the unified tropospheric stratospheric chemistry mechanism for improved forecasts of total column ozone during anomalous dynamical and chemical events. I will discuss current capabilities of the GEOS Constituent Data Assimilation System (CoDAS) to improve atmospheric composition modeling and possible future directions for GEOS-CF and reanalysis products. In addition, I will show how machine learning techniques can be used to correct for sub-grid-scale variability, which further improves model estimates at a given observation site.

Co-DAS↗
NASA NTRS2021

Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, H II Bubble Size, and Self-shielded Systems

Using the CoDa II simulation, we study the Lyα transmissivity of the intergalactic medium (IGM) during reionization. At z>6, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Lyα, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger H II regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The H II regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within H II regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68% range of the transmissivity is similar to or greater than the median for galaxies with M(sub UV) ≽ -21 at z ≽ 7, implying that more than a hundred galaxies would be needed to measure the transmission to 10% accuracy.

Reionization↗
NASA NTRSDate not supplied

Generating Global CH4 NASA GEOS Product by Assimilating TROPOMI

Radiative properties of methane (CH4) significantly contribute to climate change and the recent acceleration of global CH4 growth requires thorough investigation into its causes. Examination of temporal and spatial CH4 variability is crucial for better understanding the shifts that are currently taking place. Integrating different CH4 observations to create a reliable, easily to use global atmospheric CH4 product could contribute to assessment of emissions processes, but remains challenging due to the lack of long satellite data records. Here we present the NASA Goddard Earth Observing System (GEOS) based global CH4 product constrained by atmospheric transport from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and assimilated CH4 from the TROPOspheric Monitoring Instrument (TROPOMI). This product uses the Constituent Data Assimilation System (CoDAS) component of the GEOS capability, which has been expanded and generalized to assimilate any satellite, ground-based, or in situ observation of atmospheric composition that can be expressed as a vertical sounding or point sample. The product is similar to the one generated by assimilating carbon dioxide (CO2) data into GEOS from Orbiting Carbon Observatory 2 (OCO-2): OCO-2 GEOS Level 3 daily and monthly, 0.5x0.625 assimilated CO2 V10r. The preliminary product is investigated using a variety of quality check approaches with the help of in situ observations to examine bias correction, error inflation, and performance of the assimilation system. Potential applications of this approach include support for interpretation of high-resolution point source detection approaches, climate and greenhouse gas reanalyses, and boundary conditions for regional modeling approaches.

Nikolay Balashov↗
NASA NTRSDate not supplied

Low Latency Flux and Concentration Datasets in Support of Greenhouse Gas Monitoring Based on NASA's GEOS Modeling and Data Assimilation System

We present efforts to develop space-based greenhouse gas monitoring systems that can provide low latency information and traceability to independent observations. Through support from its Carbon Monitoring System program, NASA has developed the capability to assimilate XCO2 retrievals from the Orbiting Carbon Observatory, 2 (OCO-2) into the Goddard Earth Observing System (GEOS) Constituent Data Assimilation System (CoDAS) to create gap-filled, three-dimensional (3D) estimates of CO2 mixing ratio. When OCO-2 data are not available, concentration fields are further informed by a bottom-up flux package based on remotely sensed fire radiative power, nighttime lights, and vegetation reflectance combined with estimates of atmospheric growth rate based on surface in situ data. The 3D nature of this dataset supports evaluation with independent aircraft data, helping to ensure transparency of remotely sensed data products. These quasi-operational data are currently produced 2-3 months behind real time and are distributed via NASA and international dashboard services to a variety of end users. In this presentation, we provide an overview of the system as well as remaining data gaps and modeling challenges. We also highlight the application of this dataset for detecting emissions anomalies associated with COVID-19 and comparing against independent emissions estimates. Finally, we highlight a new NASA initiative called the Earth Information System (EIS), which aims to support open science and applications by leveraging emerging cloud computing capabilities to increase access to NASA’s greenhouse gas datasets, opportunities for co-development, and transparency in methods for analysis and flux attribution.

Lesley Ott↗
NASA NTRSDate not supplied

GEOS Constituent Data Assimilation beyond Aura MLS: Assimilating NASA SAGE III/ISS profiles of stratospheric water vapor

Water vapor in the lower stratosphere is an important trace gas for atmospheric chemistry and radiative budget with direct impact on climate. During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the stratospheric ozone loss and the subsequent recovery as well as trends in water vapor linked to surface temperature trends. The SAGE III instrument aboard the International Space Station (ISS) continues the SAGE mission record, with high vertical resolution profiles of water vapor (and other constituents) since mid-2017. The NASA GEOS Earth system model has the capability to assimilate multi-constituents from ground and space-based instruments using the GEOS Constituent Data Assimilation System (CoDAS). Reanalysis products of stratospheric water vapor without data constraints are historically poor, and two recent reanalyses which assimilated stratospheric constituents observed by the Microwave Limb Sounder (MLS) improve the representation of stratospheric composition, including water vapor, when compared against independent observations. The MLS instrument is on NASA's Aura satellite which is expected to be decommissioned in the coming years. Here we demonstrate that while the number of solar occultation observations a day from SAGE III/ISS is about 1 % of the total number of profiles observed globally by MLS, the chemical timescales of water vapor in the lower stratosphere are long enough that the SAGE III/ISS data can provide a useful constraint on the assimilated product. Analyzed water vapor fields from assimilating only SAGE III/ISS water vapor profiles agree well with independent stratospheric observations, including when compared against the frost-point hygrometers at the three NOAA stations. Over the five-year period from 2017 through 2022, the analyzed water vapor using SAGE III/ISS captures many of the features seen with the assimilation of MLS observations. Note, this ability is reduced over regions where the SAGE III/ISS instrument provides infrequent or no observations (associated with the orbit of the ISS), as highlighted in the case study period following the eruption of the tropical underwater volcano Hunga Tonga in early 2022. Nevertheless, there is a clear benefit to the assimilation of stratospheric water vapor from SAGE III/ISS observations, allowing us to continue to monitor stratospheric composition for climate assessments following the impending loss of Aura.

K. Emma Knowland↗
NASA NTRSDate not supplied

GEOS Constituent Data Assimilation Beyond Aura MLS: Assimilating NASA SAGE III/ISS Profiles of Stratospheric Water Vapor

Water vapor in the lower stratosphere is an important trace gas for atmospheric chemistry and radiative budget with direct impact on climate. During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the stratospheric ozone loss and the subsequent recovery as well as trends in water vapor linked to surface temperature trends. The SAGE III instrument aboard the International Space Station (ISS) continues the SAGE mission record, with high vertical resolution profiles of water vapor (and other constituents) since mid-2017. The NASA GEOS Earth system model has the capability to assimilate multi-constituents from ground and space-based instruments using the GEOS Constituent Data Assimilation System (CoDAS). Reanalysis products of stratospheric water vapor without data constraints are historically poor, and two recent reanalyses which assimilated stratospheric constituents observed by the Microwave Limb Sounder (MLS) improve the representation of stratospheric composition, including water vapor, when compared against independent observations. The MLS instrument is on NASA's Aura satellite which is expected to be decommissioned in the coming years. Here we demonstrate that while the number of solar occultation observations a day from SAGE III/ISS is about 1 % of the total number of profiles observed globally by MLS, the chemical timescales of water vapor in the lower stratosphere are long enough that the SAGE III/ISS data can provide a useful constraint on the assimilated product. Analyzed water vapor fields from assimilating only SAGE III/ISS water vapor profiles agree well with independent stratospheric observations, including when compared against the frost-point hygrometers at the three NOAA stations. Over the five-year period from 2017 through 2022, the analyzed water vapor using SAGE III/ISS captures many of the features seen with the assimilation of MLS observations. Note, this ability is reduced over regions where the SAGE III/ISS instrument provides infrequent or no observations (associated with the orbit of the ISS), as highlighted in the case study period following the eruption of the tropical underwater volcano Hunga Tonga in early 2022. Nevertheless, there is a clear benefit to the assimilation of stratospheric water vapor from SAGE III/ISS observations, allowing us to continue to monitor stratospheric composition for climate assessments following the impending loss of Aura.

K. Emma Knowland↗