Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Gamma-Ray Observations of the Supernova Remnant RX J0852.0-4622 with the Fermi Large Area Telescope

We report on gamma-ray observations of the supernova remnant (SNR) RX J0852.04622 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. In the Fermi-LAT data, we find a spatially extended source at the location of the SNR. The extension is consistent with the SNR size seen in other wavelengths such as X-rays and TeV gamma rays, leading to the identification of the gamma-ray source with the SNR. The spectrum is well described as a power law with a photon index of = 1.85 0.06 (stat)+0.18 0.19 (sys), which smoothly connects to the H.E.S.S. spectrum in the TeV energy band. We discuss the gamma-ray emission mechanism based on multiwavelength data. The broadband data can be fit well by a model in which the gamma rays are of hadronic origin. We also consider a scenario with inverse Compton scattering of electrons as the emission mechanism of the gamma rays. Although the leptonic model predicts a harder spectrum in the Fermi-LAT energy range, the model can fit the data considering the statistical and systematic errors.

Tanaka, T.↗
NASA NTRS2012

A Review of Low Frequency Electromagnetic Wave Phenomena Related to Tropospheric-Ionospheric Coupling Mechanisms

Investigation of coupling mechanisms between the troposphere and the ionosphere requires a multidisciplinary approach involving several branches of atmospheric sciences, from meteorology, atmospheric chemistry, and fulminology to aeronomy, plasma physics, and space weather. In this work, we review low frequency electromagnetic wave propagation in the Earth-ionosphere cavity from a troposphere-ionosphere coupling perspective. We discuss electromagnetic wave generation, propagation, and resonance phenomena, considering atmospheric, ionospheric and magnetospheric sources, from lightning and transient luminous events at low altitude to Alfven waves and particle precipitation related to solar and magnetospheric processes. We review in situ ionospheric processes as well as surface and space weather phenomena that drive troposphere-ionosphere dynamics. Effects of aerosols, water vapor distribution, thermodynamic parameters, and cloud charge separation and electrification processes on atmospheric electricity and electromagnetic waves are reviewed. We also briefly revisit ionospheric irregularities such as spread-F and explosive spread-F, sporadic-E, traveling ionospheric disturbances, Trimpi effect, and hiss and plasma turbulence. Regarding the role of the lower boundary of the cavity, we review transient surface phenomena, including seismic activity, earthquakes, volcanic processes and dust electrification. The role of surface and atmospheric gravity waves in ionospheric dynamics is also briefly addressed. We summarize analytical and numerical tools and techniques to model low frequency electromagnetic wave propagation and solving inverse problems and summarize in a final section a few challenging subjects that are important for a better understanding of tropospheric-ionospheric coupling mechanisms.

Simoes, Fernando↗
NASA NTRS2014

A Numerical Testbed for Remote Sensing of Aerosols, and its Demonstration for Evaluating Retrieval Synergy from a Geostationary Satellite Constellation of GEO-CAPE and GOES-R

We present a numerical testbed for remote sensing of aerosols, together with a demonstration for evaluating retrieval synergy from a geostationary satellite constellation. The testbed combines inverse (optimal-estimation) software with a forward model containing linearized code for computing particle scattering (for both spherical and non-spherical particles), a kernel-based (land and ocean) surface bi-directional reflectance facility, and a linearized radiative transfer model for polarized radiance. Calculation of gas absorption spectra uses the HITRAN (HIgh-resolution TRANsmission molecular absorption) database of spectroscopic line parameters and other trace species cross-sections. The outputs of the testbed include not only the Stokes 4-vector elements and their sensitivities (Jacobians) with respect to the aerosol single scattering and physical parameters (such as size and shape parameters, refractive index, and plume height), but also DFS (Degree of Freedom for Signal) values for retrieval of these parameters. This testbed can be used as a tool to provide an objective assessment of aerosol information content that can be retrieved for any constellation of (planned or real) satellite sensors and for any combination of algorithm design factors (in terms of wavelengths, viewing angles, radiance and/or polarization to be measured or used). We summarize the components of the testbed, including the derivation and validation of analytical formulae for Jacobian calculations. Benchmark calculations from the forward model are documented. In the context of NASA's Decadal Survey Mission GEOCAPE (GEOstationary Coastal and Air Pollution Events), we demonstrate the use of the testbed to conduct a feasibility study of using polarization measurements in and around the O2 A band for the retrieval of aerosol height information from space, as well as an to assess potential improvement in the retrieval of aerosol fine and coarse mode aerosol optical depth (AOD) through the synergic use of two future geostationary satellites, GOES-R (Geostationary Operational Environmental Satellite R-series) and TEMPO (Tropospheric Emissions: Monitoring of Pollution). Strong synergy between GEOS-R and TEMPO are found especially in their characterization of surface bi-directional reflectance, and thereby, can potentially improve the AOD retrieval to the accuracy required by GEO-CAPE.

Aerosols↗
NASA NTRS2017

GES DISC Greenhouse Gas Data Sets and Associated Services

NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC) archives and distributes rich collections of data on atmospheric greenhouse gases from multiple missions. Hosted data include those from the Atmospheric Infrared Sounder (AIRS) mission (which has observed CO2, CH4, ozone, and water vapor since 2002); legacy water vapor and ozone retrievals from TIROS Operational Vertical Sounder (TOVS); and Upper Atmosphere Research Satellite (UARS) going back to the early 1980s. GES DISC also archives and supports data from seven projects of the Making Earth System Data Records for Use in Research Environments (MEaSUREs) program that have ozone and water vapor records. Greenhouse gases data from the A-Train satellite constellation is also available: (1) Aura-Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) ozone, nitrous oxide, and water vapor since 2004; (2) Greenhouse Gases Observing Satellite (GOSAT) CO2 observations since 2009 from the Atmospheric CO2 Observations from Space (ACOS) task; and (3) Orbiting Carbon Observatory-2 (OCO-2) CO2 data since 2014. The most recent related data set that the GES DISC archives is methane flux for North America, as part of NASAs Carbon Monitoring System (CMS) project. This dataset contains estimates of methane emission in North America based on an inversion of the GEOS-Chem chemical transport model constrained by GOSAT observations (Turner et al., 2015). Along with data stewardship, an important focus area of the GES DISC is to enhance the usability of its data and broaden its user base. Users have unrestricted access to a new user-friendly search interface, which includes many services such as variable subsetting, format conversion, quality screening, and quick browse. The majority of the GES DISC data sets are also accessible through Open-source Project for a Network Data Access Protocol (OPeNDAP) and Web Coverage Service (WCS). The latter two services provide more options for specialized subsetting, format conversion, and image viewing. Additional data exploration, data preview, and preliminary analysis capabilities are available via NASA Giovanni, which obviates the need forusers to download the data (Acker and Leptoukh, 2007). Giovanni provides a bridge between the data and science and has been very successful in extending GES DISC data to educational users and to users with limited resources.

data ordering↗
NASA NTRS2022

Evaluation of remote sensing-based evapotranspiration products at low-latitude eddy covariance sites

Remote sensing-based evapotranspiration (ET) products have been evaluated primarily using data from northern middle latitudes; therefore, little is known about their performance at low latitudes. To address this bias, an evaluation dataset was compiled using eddy covariance data from 40 sites between latitudes 30° S and 30° N. The flux data were obtained from the emerging network in Mexico (MexFlux) and from openly available databases of FLUXNET, AsiaFlux, and OzFlux. This unique reference dataset was then used to evaluate remote sensing-based ET products in environments that have been underrepresented in earlier studies. The evaluated products were: MODIS ET (MOD16, both the discontinued collection 5 (C5) and the latest collection (C6)), Global Land Evaporation Amsterdam Model (GLEAM) ET, and Atmosphere-Land Exchange Inverse (ALEXI) ET. Products were compared with unadjusted fluxes (ETorig) and with fluxes corrected for the lack of energy balance closure (ETebc). Three common statistical metrics were used: coefficient of determination (R2), root mean square error (RMSE), and percent bias (PBIAS). The effect of a vegetation mismatch between pixel and site on product evaluation results was investigated by examining the relationship between the statistical metrics and product-specific vegetation match indexes. Evaluation results of this study and those published in the literature were used to examine the performance of the products across latitudes. Differences between the MOD16 collection 5 and 6 datasets were generally smaller than differences with the other products. Performance and ranking of the evaluated products depended on whether ETorig or ETebc was used. When using ETorig, GLEAM generally had the highest R2, smallest PBIAS, and best RMSE values across the studied land cover types and climate zones. Neither MOD16 nor ALEXI performed consistently better than the other. When using ETebc, none of the products stood out in terms of both low bias and strong correlations. The use of ETebc instead of ETorig affected the biases more than the correlations. The product evaluation results showed no significant relationship with the degree of match between the vegetation at the pixel and site scale. The latitudinal comparison showed tendencies of lower R2 (all products) but better PBIAS and normalized RMSE values (MOD16 and GLEAM) for forests at low latitudes than for forests at northern middle latitudes. For non-forest vegetation, the products showed no clear latitudinal differences in performance.

Diego Salazar-Martínez↗
NASA NTRSDate not supplied

The PACE-MAPP Polarized Ocean Bio-optical Model

The MAPP polarimeter retrieval algorithm, originally developed for the airborne NASA GISS Research Scanning Polarimeter, is being adapted to the PACE observing system to allow for a simultaneous retrieval of aerosol, thin cirrus, and ocean products using both PACE polarimeters and the OCI shortwave infrared channels. Here, we describe the polarized ocean bio-optical model we have developed for PACE-MAPP, which is designed to use a combination of homogenous and coated spheres to more accurately model the absorption and backscattering properties of algal and non-algal particles, and which also includes absorption by colored dissolved organic matter. We first use Lorenz-Mie computations to create lookup tables to accurately and efficiently characterize the optical properties of homogeneous and coated spheres embedded in water for PACE. We then use these lookup tables to develop a bio-optical model that consists of an external mixture of homogeneous (non-algal) and coated (algal) spheres. Lastly, we use the bio-optical model in forward radiative transfer computations to perform inversions of MODIS remote sensing reflectance for a variety of ocean regions. We explore the sensitivity of these inversions to ocean constituent inherent optical properties.

James Allen↗
NASA NTRS2023

Development of High-Resolution Soil Hydraulic Parameters with Use of Earth Observations for Enhancing Root Zone Soil Moisture Product

Regional quantification of energy and water balance fluxes depends inevitably on the estimation of surface and rootzone soil moisture. The simulation of soil moisture depends on the soil retention characteristics, which are difficult to estimate at a regional scale. Thus, the present study proposes a new method to estimate high-resolution Soil Hydraulic Parameters (SHPs) which in turn help to provide high-resolution (spatial and temporal) rootzone soil moisture (RZSM) products. The study is divided into three phases—(I) involves the estimation of finer surface soil moisture (1 km) from the coarse resolution satellite soil moisture. The algorithm utilizes MODIS 1 km Land Surface Temperature (LST) and 1 km Normalized difference vegetation Index (NDVI) for downscaling 25 km C-band derived soil moisture from AMSR-2 to 1 km surface soil moisture product. At one of the test sites, soil moisture is continuously monitored at 5, 20, and 50 cm depth, while at 44 test sites data were collected randomly for validation. The temporal and spatial correlation for the downscaled product was 70% and 83%, respectively. (II) In the second phase, downscaled soil moisture product is utilized to inversely estimate the SHPs for the van Genuchten model (1980) at 1 km resolution. The numerical experiments were conducted to understand the impact of homogeneous SHPs as compared to the three-layered parameterization of the soil profile. It was seen that the SHPs estimated using the downscaled soil moisture (I-d experiment) performed with similar efficiency as compared to SHPs estimated from the in-situ soil moisture data (I-b experiment) in simulating the soil moisture. The normalized root mean square error (nRMSE) for the two treatments was 0.37 and 0.34, respectively. It was also noted that nRMSE for the treatment with the utilization of default SHPs (I-a) and AMSR-2 soil moisture (I-c) were found to be 0.50 and 0.43, respectively. (III) Finally, the derived SHPs were used to simulate both surface soil moisture and RZSM. The final product, RZSM which is the daily 1 km product also showed a nearly 80% correlation at the test site. The estimated SHPs are seen to improve the mean NSE from 0.10 (I-a experiment) to 0.50 (I-d experiment) for the surface soil moisture simulation. The mean nRMSE for the same was found to improve from 0.50 to 0.31.

soil moisture downscaling↗
NASA NTRSDate not supplied

Geophysical Retrievals and Cloud Analyses from Merged Airborne Radiometer Datasets Covering 10–684 GHz

Airborne microwave radiometers provide insight about numerous aspects of Earth’s atmosphere and yield critical validation datasets for spaceborne radiometers. Three radiometers that are important to NASA’s airborne remote-sensing arsenal include: the Advanced Microwave Precipitation Radiometer (AMPR), covering 10–85 GHz; the Conical Scanning Millimeter-wave Imaging Radiometer (CoSMIR), covering 50–183 GHz; and the Compact Scanning Submillimeter-wave Imaging Radiometer (CoSSIR), covering 170–684 GHz. The NASA field campaigns of interest to this study include: the Integrated Precipitation and Hydrology Experiment (IPHEx) in 2014, the Olympic Mountains Experiment and Radar Definition Experiment (OLYMPEX/RADEX) in 2015–2016, the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) in 2020–2023, and the Airborne Lightning Observatory for FEGS and TGFs (ALOFT) in 2023. To provide a more comprehensive perspective on clouds and precipitation observed during these airborne field campaigns, AMPR data were merged spatiotemporally with CoSMIR data for IPHEx, OLYMPEX/RADEX, and IMPACTS (2020 and 2022), while considering differences in instrument characteristics and operations, providing brightness temperature (Tb) values from 10–183 GHz in a common background grid throughout each flight. AMPR and CoSSIR data were similarly merged for IMPACTS (2023) and ALOFT, providing a common background grid with Tb values covering 10–684 GHz throughout each flight. These merged Tb data were employed in geophysical retrievals using the Community Radiative Transfer Model (CRTM), an Eddington radiative transfer model, and a one-dimensional variational (1DVAR) inversion method. Retrievals of cloud liquid water path were of primary interest. This presentation will include an overview of the methods for the radiometer data mergers, the radiative transfer methods, the geophysical retrievals, and detailed results from examining trends in Tb and cloud liquid water path in clouds, precipitation, and cloud-to-precipitation transition zones.

Corey G Amiot↗
NASA NTRS1996

Temperature, Pressure, and Infrared Image Survey of an Axisymmetric Heated Exhaust Plume

The focus of this research is to numerically predict an infrared image of a jet engine exhaust plume, given field variables such as temperature, pressure, and exhaust plume constituents as a function of spatial position within the plume, and to compare this predicted image directly with measured data. This work is motivated by the need to validate computational fluid dynamic (CFD) codes through infrared imaging. The technique of reducing the three-dimensional field variable domain to a two-dimensional infrared image invokes the use of an inverse Monte Carlo ray trace algorithm and an infrared band model for exhaust gases. This report describes an experiment in which the above-mentioned field variables were carefully measured. Results from this experiment, namely tables of measured temperature and pressure data, as well as measured infrared images, are given. The inverse Monte Carlo ray trace technique is described. Finally, experimentally obtained infrared images are directly compared to infrared images predicted from the measured field variables.

Nelson, Edward L.↗
DOE OSTI2025

Advancements and opportunities to improve bottom–up estimates of global wetland methane emissions

Wetlands are the single largest natural source of atmospheric methane (CH 4 ), contributing approximately 30% of total surface CH 4 emissions, and they have been identified as the largest source of uncertainty in the global CH 4 budget based on the most recent Global Carbon Project CH 4 report. High uncertainties in the bottom–up estimates of wetland CH 4 emissions pose significant challenges for accurately understanding their spatiotemporal variations, and for the scientific community to monitor wetland CH 4 emissions from space. In fact, there are large disagreements between bottom–up estimates versus top–down estimates inferred from inversion of atmospheric CH 4 concentrations. To address these critical gaps, we review recent development, validation, and applications of bottom–up estimates of global wetland CH 4 emissions, as well as how they are used in top–down inversions. These bottom–up estimates, using (1) empirical biogeochemical modeling (e.g. WetCHARTs: 125–208 TgCH 4 yr -1 ); (2) process-based biogeochemical modeling (e.g. WETCHIMP: 190 ± 39 TgCH 4 yr -1 ); and (3) data-driven machine learning approach (e.g. UpCH4: 146 ± 43 TgCH 4 yr -1 ). Bottom–up estimates are subject to significant uncertainties (~80 Tg CH 4 yr -1 ), and the ranges of different estimates do not overlap, further amplifying the overall uncertainty when combining multiple data products. These substantial uncertainties highlight gaps in our understanding of wetland CH 4 biogeochemistry and wetland inundation dynamics. Major tropical and arctic wetland complexes are regional hotspots of CH 4 emissions. However, the scarcity of satellite data over the tropics and northern high latitudes offer limited information for top–down inversions to improve bottom–up estimates. Recent advances in surface measurements of CH 4 fluxes (e.g. FLUXNET-CH 4 ) across a wide range of ecosystems including bogs, fens, marshes, and forest swamps provide an unprecedented opportunity to improve existing bottom–up estimates of wetland CH 4 estimates. We suggest that continuous long-term surface measurements at representative wetlands, high fidelity wetland mapping, combined with an appropriate modeling framework, will be needed to significantly improve global estimates of wetland CH 4 emissions. There is also a pressing unmet need for fine-resolution and high-precision satellite CH 4 observations directed at wetlands.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) and 10 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become trapped and radiated from an axion sphere. The trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a diffusion approximation for axion-energy transport. The axion opacity due to inverse nucleon-nucleon, axion bremsstrahlung is computed; and then the numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The estimate of the axion mass is confirmed above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration - the most sensitive barometer of axion cooling - it is concluded that for an axion mass greater than about 3 eV axion emission would not have had a significant effect on the neutrino bursts detected by KII and IMB. It is strongly suggested that an axion with mass in the interval 10(exp -3) to 3 eV is excluded by the observation of neutrinos from SN 1987A.

Burrows, Adam↗
NASA NTRS1998

Inversion of Crater Morphometric Data to Gain Insight on the Cratering Process

In recent years, morphometric data for Venus and several outer planet satellites have been collected, so we now have observational data of complex Craters formed in a large range of target properties. We present general inversion techniques that can utilize the morphometric data to quantitatively test various models of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic pressure on a planet, or the factor c/pg.

Herrick, Robert R.↗
NASA NTRS1971

Vibrational population distributions in nonequilibrium nozzle expansion flows

Experimental measurements and theoretical calculations of the vibrational population distribution in nonequilibrium nozzle expansion flows of gas mixtures are reported. These studies were directed toward determining whether vibrational energy exchange pumping could lead to laser action on the vibrational bands of a diatomic molecule. Three different types of experiments were conducted. These showed (1) that vibrational energy was preferentially transferred from N2 to CO in supersonic nozzle flows containing these gases; (2) that under some conditions this vibrational energy exchange pumping mechanism created population inversions in the vibrational levels of CO; and (3) that at large expansion ratios the magnitude of these population inversions was sufficient to sustain lasing in the nozzle. A theoretical model was developed to calculate vibrational state population distributions in gas dynamic expansions of a mixture of diatomic gases. Although only isothermal calculations have been completed, these data indicate that population inversions are predicted for conditions similar to those obtained in the nozzle expansion flows.

Watt, W. S.↗
NASA NTRS1979

Direct and inverse problems in radiation of sound from discrete random sources on two coaxial rings

An analytical model consisting of two ring sources of sound is developed to study the direct radiation in terms of correlation, coherence, and phase and also to aid in solving the inverse-radiation problem of determining the noise source in terms of farfield measurements. The rings consist of discrete sources which are either monopoles or quadrupoles with Gaussian autocorrelations. Only adjacent sources, both within and between the rings, are correlated. Results show that from the farfield information one can determine when the sources are compact or noncompact with respect to the acoustic wavelength and distinguish between the types of sources. In addition, from the inverse-radiation approach one can recover the center of mass, the location and separation distance of the ring, and the respective diameters.

Maestrello, L.↗
DOE OSTI2021

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗