Search NASA⌕ Search

SEARCH · Search NASA

Results for “geophysical methods”

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 145 records · Page 8

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗

Novel CCS Monitoring of Approved Class VI Storage in North Dakota

Presentation at Dallas Geophysical Society Meeting, Dallas, Texas, April 24, 2025. This high-level overview of focuses on novel and sustainable monitoring methods at various stages of planning or demonstration to accelerate the deployment at future carbon capture, utilization, and storage (CCUS) sites.

02 PETROLEUM↗

Measurement of the Vertical Gradient of the Semidiurnal Tidal Wind Phase in Winter at the 95 Km Level

When supplemented by absolute reflection height measurements, low frequency wind measurements in the 90-100 km height range become truly competitive in comparison with the more widely used radar meteor wind observations. For example, height profiles of the wind parameters in the so-called meteor zone can be obtained due to the considerable interdiurnal variability of the average nighttime reflection heights controlled by geomagnetic activity. The phase of the semidiurnal tidal wind is particularly height-dependent. The measured vertical gradient of 1/4 h/km in winter corresponds to a vertical wavelength of about 50 km. Wind measurements in the upper atmosphere, at heights between 90 and 100 km, were carried out at the Collm Geophysical Observatory of Karl Marx University Leipzig for a number of years. These measurements use the closely-spaced receiver method and three measuring paths, on 179, 227, and 272 kHz. They take place every day between sunset and sunrise, i.e., nightly. A night in this sense may last as long as 18 hours in winter. Both the measurements and their evaluation are completely automatic, and the prevailing winds and tides are separated.

Schminder, R.↗

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

Sensitivity of Spacebased Microwave Radiometer Observations to Ocean Surface Evaporation

Ocean surface evaporation and the latent heat it carries are the major components of the hydrologic and thermal forcing on the global oceans. However, there is practically no direct in situ measurements. Evaporation estimated from bulk parameterization methods depends on the quality and distribution of volunteer-ship reports which are far less than satisfactory. The only way to monitor evaporation with sufficient temporal and spatial resolutions to study global environment changes is by spaceborne sensors. The estimation of seasonal-to-interannual variation of ocean evaporation, using spacebased measurements of wind speed, sea surface temperature (SST), and integrated water vapor, through bulk parameterization method,s was achieved with reasonable success over most of the global ocean, in the past decade. Because all the three geophysical parameters can be retrieved from the radiance at the frequencies measured by the Scanning Multichannel Microwave Radiometer (SMMR) on Nimbus-7, the feasibility of retrieving evaporation directly from the measured radiance was suggested and demonstrated using coincident brightness temperatures observed by SMMR and latent heat flux computed from ship data, in the monthly time scale. However, the operational microwave radiometers that followed SMMR, the Special Sensor Microwave/Imager (SSM/I), lack the low frequency channels which are sensitive to SST. This low frequency channels are again included in the microwave imager (TMI) of the recently launched Tropical Rain Measuring Mission (TRMM). The radiance at the frequencies observed by both TMI and SSM/I were simulated through an atmospheric radiative transfer model using ocean surface parameters and atmospheric temperature and humidity profiles produced by the reanalysis of the European Center for Medium Range Weather Forecast (ECMWF). From the same ECMWF data set, coincident evaporation is computed using a surface layer turbulent transfer model. The sensitivity of the radiance to evaporation over various seasons and geographic locations are examined. The microwave frequencies with radiance that are significant correlated with evaporation are identify and capability of estimating evaporation directly from TMI will be discussed.

Liu, Timothy W.↗

New Earth Science Data and Access Methods

NASA's Earth Science Enterprise, working with its domestic and international partners, provides scientific data and analysis to improve life here on Earth. NASA provides science data products that cover a wide range of physical, geophysical, biochemical and other parameters, as well as services for interdisciplinary Earth science studies. Management and distribution of these products is administered through the Earth Observing System Data and Information System (EOSDIS) Distributed Active Archive Centers (DAACs), which all hold data within a different Earth science discipline. This paper will highlight selected EOS datasets and will focus on how these observations contribute to the improvement of essential services such as weather forecasting, climate prediction, air quality, and agricultural efficiency. Emphasis will be placed on new data products derived from instruments on board Terra, Aqua and ICESat as well as new regional data products and field campaigns. A variety of data tools and services are available to the user community. This paper will introduce primary and specialized DAAC-specific methods for finding, ordering and using these data products. Special sections will focus on orienting users unfamiliar with DAAC resources, HDF-EOS formatted data and the use of desktop research and application tools.

Moses, John F.↗

Global analyses of water vapor, cloud and precipitation derived from a diagnostic assimilation of SSM/I geophysical retrievals

An analytical approach is described for diagnostically assimilating moisture data from Special Sensor Microwave Imager (SSM/I) into a global analysis of water vapor, cloud content, and precipitation. In this method, 3D fields of wind and temperature values taken from ECMWF gridded analysis are used to drive moisture conservation equations with parameterized microphysical treatment of vapor, liquid, and ice; the evolving field of water vapor is periodically updated or constrained by SSM/I retrievals of precipitable water. Initial results indicate that this diagnostic model can produce realistic large-scale fields of cloud and precipitation. The resulting water vapor analyses agree well with SSM/I and have an additional advantage of being synoptic.

Robertson, Franklin R.↗

Investigation of finite element: ABC methods for electromagnetic field simulation

The mechanics of wave propagation in the presence of obstacles is of great interest in many branches of engineering and applied mathematics like electromagnetics, fluid dynamics, geophysics, seismology, etc. Such problems can be broadly classified into two categories: the bounded domain or the closed problem and the unbounded domain or the open problem. Analytical techniques have been derived for the simpler problems; however, the need to model complicated geometrical features, complex material coatings and fillings, and to adapt the model to changing design parameters have inevitably tilted the balance in favor of numerical techniques. The modeling of closed problems presents difficulties primarily in proper meshing of the interior region. However, problems in unbounded domains pose a unique challenge to computation, since the exterior region is inappropriate for direct implementation of numerical techniques. A large number of solutions have been proposed but only a few have stood the test of time and experiment. The goal of this thesis is to develop an efficient and reliable partial differential equation technique to model large three dimensional scattering problems in electromagnetics.

Chatterjee, A.↗

Carbon Storage Core Characterization Efforts at NETL

The multi-scale Computed Tomography (CT) and core flow facility in the Geocharacterization Laboratory at NETL, Morgantown yields porosity, permeability, and fracture properties of rock core samples obtained from the subsurface while maintaining the integrity of the sample. Additionally, geophysical bulk rock properties are analyzed with the laboratory’s GeoTEK multi-sensor core logger in a comparable fashion to downhole methods. NETL researchers collaborate with stakeholders within the carbon storage, oil and gas, and critical minerals sectors. Since 2017, over 1.88 miles of core have been analyzed within the laboratory and all data is publicly available through the Technical Report Series (TRS) on the Energy Data eXchange (EDX). Additionally, the website, RokBase, was curated to extrapolate and visualize the high-resolution data from field operations. The characterization of the Lively Grove #1 (LG#1) well provides a case study into the full capabilities of the Geocharacterization Laboratory. During the comprehensive study of LG#1 ~1-2 mm in diameter, vertical to bedding, cylindrical structures were identified throughout the St. Peter Formation. These structures are pervasive throughout the St. Peter Formation at depth and are characterized as the trace fossil, Skolithos.

Isom, Shelby L↗

The strong equivalence principle and its violation

Theoretical and observational aspects of an SEP violation are discussed. A two-times theory is presented as a possible framework to handle an SEP violation, and tests performed to check the compatibility of such violation with a host of data ranging from nucleosynthesis to geophysics are summarized. Also discussed are the dynamical equations needed to analyze radar ranging data to reveal an SEP violation and in particular the method employed by Shapiro and Reasenberg (1976).

Canuto, V. M.↗

Methods of training the graduate level and professional geologist in remote sensing technology

Requirements for a basic course in remote sensing to accommodate the needs of the graduate level and professional geologist are described. The course should stress the general topics of basic remote sensing theory, the theory and data types relating to different remote sensing systems, an introduction to the basic concepts of computer image processing and analysis, the characteristics of different data types, the development of methods for geological interpretations, the integration of all scales and data types of remote sensing in a given study, the integration of other data bases (geophysical and geochemical) into a remote sensing study, and geological remote sensing applications. The laboratories should stress hands on experience to reinforce the concepts and procedures presented in the lecture. The geologist should then be encouraged to pursue a second course in computer image processing and analysis of remotely sensed data.

Kolm, K. E.↗

The scatterometer - Data and applications

The techniques used to extract oceanographic and meteorological parameters from satellite radar scatterometer data are described and demonstrated, primarily on the basis of Seasat data. The principles of backscatter measurement and the general relationships between wind and wave fields are summarized, and the construction of geophysical data records from scatterometer swaths is explained. Particular attention is given to surface-layer dynamics and stresses; planetary-boundary-layer dynamics; averaging methods; the estimation of heat, moisture, and CO2 flux coefficients; the determination of pressure fields; storm monitoring; synoptic-scale analysis; climate analysis; and ocean dynamics. The capabilities and potential applications of the planned ERS-1, NSCAT, and SCANSCAT satellite scatterometers are briefly indicated.

Brown, R. A.↗

Passive Microwave Remote Sensing of Soil Moisture

Microwave remote sensing provides a unique capability for direct observation of soil moisture. Remote measurements from space afford the possibility of obtaining frequent, global sampling of soil moisture over a large fraction of the Earth's land surface. Microwave measurements have the benefit of being largely unaffected by cloud cover and variable surface solar illumination, but accurate soil moisture estimates are limited to regions that have either bare soil or low to moderate amounts of vegetation cover. A particular advantage of passive microwave sensors is that in the absence of significant vegetation cover soil moisture is the dominant effect on the received signal. The spatial resolutions of passive Microwave soil moisture sensors currently considered for space operation are in the range 10-20 km. The most useful frequency range for soil moisture sensing is 1-5 GHz. System design considerations include optimum choice of frequencies, polarizations, and scanning configurations, based on trade-offs between requirements for high vegetation penetration capability, freedom from electromagnetic interference, manageable antenna size and complexity, and the requirement that a sufficient number of information channels be available to correct for perturbing geophysical effects. This paper outlines the basic principles of the passive microwave technique for soil moisture sensing, and reviews briefly the status of current retrieval methods. Particularly promising are methods for optimally assimilating passive microwave data into hydrologic models. Further studies are needed to investigate the effects on microwave observations of within-footprint spatial heterogeneity of vegetation cover and subsurface soil characteristics, and to assess the limitations imposed by heterogeneity on the retrievability of large-scale soil moisture information from remote observations.

Njoku, Eni G.↗

PRELIMINARY MODEL FOR SOIL MOISTURE RETRIEVAL USING P-BAND RADIOMETER OBSERVATIONS

Soil Moisture is an important geophysical variable that needs reliable quantification for applications in hydrology, meteorology and agriculture. L-band radiometry has proved to be one of the best methods in soil moisture estimation using microwave signals. However, they provide measurements that correspond to a shallow depth of 5 cm and are also affected by the presence of overlaying vegetation and roughness. In contrast, P-band radiometry is expected to provide moisture information on a deeper layer of soil. Moreover, these lower frequency measurements are expected to be less affected by soil roughness and vegetation contributions. Consequently, this pilot study uses the Polarimetric P-band Multi beam Radiometer(PPMR) at 740 MHz to evaluate the response of theP-band radiometer over a realistic range of surface condition sat the field scale. A preliminary framework of P-band MicrowaveEmission of theBiosphere (P-MEB) has been developed as a forward model that simulates brightness temperature from soil moisture and other ancillary data collected from the field. This paper presents the model for the bare soil condition observed during June 2018 to August 2018. The results show that H-polarized PPMR data has better correlation to the soil moisture over a depth of 10 cm than the V-polarized PPMR data. A model is under improvement by incorporating a more suitable effective temperature formulation.

Edward J Kim↗

Supporting Responsible Machine Learning in Heliophysics

Over the last decade, Heliophysics researchers have increasingly adopted a variety of machine learning methods such as artificial neural networks, decision trees, and clustering algorithms into their workflow. Adoption of these advanced data science methods had quickly outpaced institutional response, but many professional organizations such as the European Commission, the National Aeronautics and Space Administration (NASA), and the American Geophysical Union have now issued (or will soon issue) standards for artificial intelligence and machine learning that will impact scientific research. These standards add further (necessary) burdens on the individual researcher who must now prepare the public release of data and code in addition to traditional paper writing. Support for these is not reflected in the current state of institutional support, community practices, or governance systems. We examine here some of these principles and how our institutions and community can promote their successful adoption within the Heliophysics discipline.

Machine learning↗

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗

Application of Heat Capacity Mapping Mission data to regional geologic analysis for mineral and energy resource evaluation

Heat Capacity Mapping Mission thermal-inertia images of a diversity of terrains and geologic settings were examined in conjunction with topographic, geologic, geophysical, and LANDSAT data. The images were found to have attributes similar to bedrock maps. In the Cascades region, two new features were identified and a method was developed to characterize regional terranes using linear feature data. Two northeast-trending Lineaments were discovered in the Overthrust Belt of Montana and Idaho. The longer of the two extends from the Idaho-Oregon border, through the Idaho batholith and across the Lewis thrust. It coincides, along segments, with mapped faults and an aeromagnetic pattern change. A major lineament crossing the Colorado Plateau and the Southern Rocky Mountians was detected on several thermal-inertial images and evidence was found for the existence of a geologic discontinuity. Vegetation-covered areas in Richfield and the Silver City quadrangle (Arizona and New Mexico) displayed thermal-inertia differences within heavily vegetation areas although no apreciable correlation was found between vegetation cover and thermal inertia. Resistant ridges and knolls have high thermal inertias and thermal-inertia contrasts occurred at lithologic and fault contacts. In the heavy vegetated Pinaleno Mountains, Arizona, a Lithologic unit obscured on LANDSAT MSS data due to the vegetation cover, exhibited a thermal-inertia contrast with its surroundings.

Watson, K.↗