Search NASA⌕ Search

SEARCH · Search NASA

Results for “model analyses”

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 55 records · Page 3

Elements of complexity in subsurface modeling, exemplified with three case studies

There are complexity elements to consider when applying subsurface flow and transport models to support environmental analyses. Modelers balance the benefits and costs of modeling along the spectrum of complexity, taking into account the attributes of more simple models (e.g., lower cost, faster execution, easier to explain, less mechanistic) and the attributes of more complex models (higher cost, slower execution, harder to explain, more mechanistic and technically defensible). In this paper, modeling complexity is examined with respect to considering this balance. The discussion of modeling complexity is organized into three primary elements: 1) modeling approach, 2) description of process, and 3) description of heterogeneity. Three examples are used to examine these complexity elements. Two of the examples use simulations generated from a complex model to develop simpler models for efficient use in model applications. The first example is designed to support performance evaluation of soil vapor extraction remediation in terms of groundwater protection. The second example investigates the importance of simulating different categories of geochemical reactions for carbon sequestration and selecting appropriate simplifications for use in evaluating sequestration scenarios. In the third example, the modeling history for a uranium-contaminated site demonstrates that conservative parameter estimates were inadequate surrogates for complex, critical processes and there is discussion on the selection of more appropriate model complexity for this application. All three examples highlight how complexity considerations are essential to create scientifically defensible models that achieve a balance between model simplification and complexity.

conceptual models↗

Nonnegative tensor factorization for contaminant source identification

Unsupervised Machine Learning (ML) is becoming increasingly popular for solving various types of data analytics problems including feature extraction, blind source separation, exploratory analyses, model diagnostics, etc. In this work, we have developed a new unsupervised ML method based on Nonnegative Tensor Factorization (NTF) for identification of the original groundwater types (including contaminant sources) present in geochemical mixtures observed in an aquifer. Frequently, groundwater types with different geochemical signatures are related to different background and/or contamination sources. The characterization of groundwater mixing processes is a challenging but very important task critical for any environmental management project aiming to characterize the fate and transport of contaminants in the subsurface and perform contaminant remediation. This task typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. Additionally, the application of inverse methods may introduce biases in the analyses through the assumptions made in the model development process. Here, we substitute the model inversion with unsupervised ML analysis. The ML analysis does not make any assumptions about underlying physical and geochemical processes occurring in the aquifer. Our ML methodology, called NTFk, is capable of identifying (1) the unknown number of groundwater types (contaminant sources) present in the aquifer, (2) the original geochemical concentrations (signatures) of these groundwater types and (3) spatial and temporal dynamics in the mixing of these groundwater types. These results are obtained only from the measured geochemical data without any additional site information. In general, the NTFk methodology allows for interpretation of large high-dimensional datasets representing diverse spatial and temporal components such as state variables and velocities. NTFk has been tested on synthetic and real-world site three-dimensional datasets. Finally, the NTFk algorithm is designed to work with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).

58 GEOSCIENCES↗

LDEF geometry/mass model for radiation analyses

A three-dimensional geometry/mass model of LDEF is under development for ionizing radiation analyses. This model, together with ray tracing algorithms, is being programmed for use both as a stand alone code in determining three-dimensional shielding distributions at dosimetry locations and as a geometry module that can be interfaced with radiation transport codes.

Colborn, B. L.↗

Utah FORGE Phase Native State FALCON Model Files

The submission includes FALCON input file and mesh for the an initial pressure-temperature simulation, and a second set for pressure-temperature-displacement simulation. All simulations are steady state. Data and input for the FORGE Phase 2 native state model were compiled from historic and FORGE specific characterization efforts. The simulations look to find the in-situ pressure, temperature, and stress distributions in the vicinity of the Utah FORGE site. The purpose of the native state model was to establish a baseline for all other modeling analyses—so that the whole modeling team can start from a consistent set of conditions which makes for comparison of results more straightforward and meaningful. Requires the FALCON simulator which can be found in the "FALCON Simulator" link in the submission. Instructions for compiling and running FALCON can be found on the FORGE branch of the FALCON git repository in the README.

15 GEOTHERMAL ENERGY↗

Air Quality Modeling Using the NASA GEOS-5 Multispecies Data Assimilation System

The NASA Goddard Earth Observing System (GEOS) data assimilation system (DAS) has been expanded to include chemically reactive tropospheric trace gases including ozone (O3), nitrogen dioxide (NO2), and carbon monoxide (CO). This system combines model analyses from the GEOS-5 model with detailed atmospheric chemistry and observations from MLS (O3), OMI (O3 and NO2), and MOPITT (CO). We show results from a variety of assimilation test experiments, highlighting the improvements in the representation of model species concentrations by up to 50% compared to an assimilation-free control experiment. Taking into account the rapid chemical cycling of NO2 when applying the assimilation increments greatly improves assimilation skills for NO2 and provides large benefits for model concentrations near the surface. Analysis of the geospatial distribution of the assimilation increments suggest that the free-running model overestimates biomass burning emissions but underestimates lightning NOx emissions by 5-20%. We discuss the capability of the chemical data assimilation system to improve atmospheric composition forecasts through improved initial value and boundary condition inputs, particularly during air pollution events. We find that the current assimilation system meaningfully improves short-term forecasts (1-3 day). For longer-term forecasts more emphasis on updating the emissions instead of initial concentration fields is needed.

Keller, Christoph A.↗

Space Shuttle filament wound case compressive strength study. II - Analysis

In order to explain the results of the testing conducted in Part I and relate the subscale test specimen performance to case performance during prelaunch loading, a comprehensive analysis program was conducted. Finite element models of both full-scale case and subscale test coupons were made. Validated axisymmetric models were used for the case analyses. Modeling of coupons was accomplished using a generalized plane strain finite element model in most cases, although some three-dimensional finite element modeling was used to investigate the importance of three-dimensional effects on the coupon performance. The large number of plies and complexity of the test coupons used limited the finite element grids to a through-the-thickness discretization of one element per ply. Delaminations observed in the helical ply dropoff region were modeled using frictionless sliders. Comparisons of predicted strains and displacements were made with the test data from both cases and coupons to validate the finite element models. Based on observations of the test results and finite element predictions, a failure criterion for both coupons and cases was postulated and then validated by comparing performance predictions to test results. Finally, the validated failure criterion was used to predict the performance of the first two flight cases under prelaunch loads.

Messick, M. J.↗

Diagnosing AIRS Sampling with CloudSat Cloud Classes

AIRS yield and sampling vary with cloud state. Careful utilization of collocated multiple satellite sensors is necessary. Profile differences between AIRS and ECMWF model analyses indicate that AIRS has high sampling and excellent accuracy for certain meteorological conditions. Cloud-dependent sampling biases may have large impact on AIRS L2 and L3 data in climate research. MBL clouds / lower tropospheric stability relationship is one example. AIRS and CloudSat reveal a reasonable climatology in the MBL cloud regime despite limited sampling in stratocumulus. Thermodynamic parameters such as EIS derived from AIRS data map these cloud conditions successfully. We are working on characterizing AIRS scenes with mixed cloud types.

model analyses↗

4-D Cloud Water Content Fields Derived from Operational Satellite Data

In order to improve operational safety and efficiency, the transportation industry, including aviation, has an urgent need for accurate diagnoses and predictions of clouds and associated weather conditions. Adverse weather accounts for 70% of all air traffic delays within the U.S. National Airspace System. The Federal Aviation Administration has determined that as much as two thirds of weather-related delays are potentially avoidable with better weather information and roughly 20% of all aviation accidents are weather related. Thus, it is recognized that an important factor in meeting the goals of the Next Generation Transportation System (NexGen) vision is the improved integration of weather information. The concept of a 4-D weather cube is being developed to address that need by integrating observed and forecasted weather information into a shared 4-D database, providing an integrated and nationally consistent weather picture for a variety of users and to support operational decision support systems. Weather analyses and forecasts derived using Numerical Weather Prediction (NWP) models are a critical tool that forecasters rely on for guidance and also an important element in current and future decision support systems. For example, the Rapid Update Cycle (RUC) and the recently implemented Rapid Refresh (RR) Weather Research and Forecast (WRF) models provide high frequency forecasts and are key elements of the FAA Aviation Weather Research Program. Because clouds play a crucial role in the dynamics and thermodynamics of the atmosphere, they must be adequately accounted for in NWP models. The RUC, for example, cycles at full resolution five cloud microphysical species (cloud water, cloud ice, rain, snow, and graupel) and has the capability of updating these fields from observations. In order to improve the models initial state and subsequent forecasts, cloud top altitude (or temperature, T(sub c)) derived from operational satellite data, surface observations of cloud base altitude, radar reflectivity, and lightning data are used to help build and remove clouds in the models assimilation system. Despite this advance and the many recent advances made in our understanding of cloud physical processes and radiative effects, many problems remain in adequately representing clouds in models. While the assimilation of cloud top information derived from operational satellite data has merit, other information is available that has not yet been exploited. For example, the vertically integrated cloud water content (CWC) or cloud water path (CWP) and cloud geometric thickness (delta Z) are standard products being derived routinely from operational satellite data. These and other cloud products have been validated under a variety of conditions. Since the uncertainties have generally been found to be less than those found in model analyses and forecasts, the satellite products should be suitable for data assimilation, provided an appropriate strategy can be developed that links the satellite-derived cloud parameters with cloud parameters specified in the model. In this paper, we briefly outline such a strategy and describe a methodology to retrieve cloud water content profiles from operational satellite data. Initial results and future plans are presented. It is expected that the direct assimilation of this new product will provide the most accurate depiction of the vertical distribution of cloud water ever produced at the high spatial and temporal resolution needed for short term weather analyses and forecasts.

Smith, William L., Jr.↗

LDEF geometry/mass model for radiation analyses

A 3-D geometry/mass model of Long Duration Exposure Facility (LDEF) was generated to aid in the interpretation of ionizing radiation experiments relative to the influence of varying shielding distributions around the dosimetry, and to allow more definitive calculations and comparisons with the measured data. This model takes into account the major individual structural members of the LDEF spacecraft, the mass in each experiment tray, and, for selected trays containing ionizing radiation dosimetry, major parts within the tray are modeled. The geometry/mass model, together with a ray tracing algorithms, was programmed for use both as a stand alone code in determining 3-D shielding distributions at dosimetry locations, and as a geometry module that can be interfaced with radiation transport codes.

Colborn, B. L.↗

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗

Air–Sea Turbulent Heat Fluxes in Climate Models and Observational Analyses: What Drives Their Variability?

A traditional view is that the ocean outside of the tropics responds passively to atmosphere forcing, which implies that air–sea heat fluxes are mainly driven by atmosphere variability. This paper tests this viewpoint using state-of-the-art air–sea turbulent heat flux observational analyses and a climate model run at different resolutions. It is found herein that in midlatitude ocean frontal zones the variability of air–sea heat fluxes is not predominantly driven by the atmosphere variations but instead is forced by sea surface temperature (SST) variations arising from intrinsic oceanic variability. Meanwhile in most of the tropics and subtropics wind is the dominant driver of heat flux variability, and atmosphere humidity is mainly important in higher latitudes. The predominance of ocean forcing of heat fluxes found in frontal regions occurs on scales of around 700 km or less. Spatially smoothing the data to larger scales results in the traditional atmosphere-driving case, while filtering to retain only small scales of 5° or less leads to ocean forcing of heat fluxes over most of the globe. All observational analyses examined (1° OAFlux; 0.25° J-OFURO3; 0.25° SeaFlux) show this general behavior. A standard resolution (1°) climate model fails to reproduce the midlatitude, small-scale ocean forcing of heat flux: refining the ocean grid to resolve eddies (0.1°) gives a more realistic representation of ocean forcing but the variability of both SST and of heat flux is too high compared to observational analyses.

54 ENVIRONMENTAL SCIENCES↗

Pajarito models

Pajarito models are developed to perform various types of model analyses related to the LANL-site groundwater contamination problems. Pajarito models include several models representing the LANL Chromium (Cr) and RDX sites. Pajarito models also include a model that is at the scale of the Pajarito Plateau. The models are designed to perform model prediction of groundwater flow and contaminant transport. The models are bundled in repositories. The model repositories include scripts, input control files, and drivers. A detailed description of the various files included in the repositories is provided in electronic form within the repositories.

Vesselinov, Velimir↗

Surface temperature comparison of the Arctic winter MOSAiC observations, ERA5 reanalysis, and MODIS satellite retrieval

Atmospheric model systems, such as those used for weather forecast and reanalysis production, often have significant and systematic errors in their representation of the Arctic surface energy budget and its components. The newly available observation data of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (2019/2020) enable a range of model analyses and validation in order to advance our understanding of potential model deficiencies. In the present study, we analyze deficiencies in the surface radiative energy budget over Arctic sea ice in the ERA5 global atmospheric reanalysis by comparing against the winter MOSAiC campaign data, as well as, a pan-Arctic level-2 MODIS ice surface temperature remote sensing product. We find that ERA5 can simulate the timing of radiatively clear periods, though it is not able to distinguish the two observed radiative Arctic winter states, radiatively clear and opaquely cloudy, in the distribution of the net surface radiative budget. The ERA5 surface temperature over Arctic sea ice has a conditional error with a positive bias in radiatively clear conditions and a negative bias in opaquely cloudy conditions. The mean surface temperature error is 4°C for radiatively clear situations at MOSAiC and up to 15°C in some parts of the Arctic. The spatial variability of the surface temperature, given by 4 observation sites at MOSAiC, is not captured by ERA5 due to its spatial resolution but represented in the level-2 satellite product. The sensitivity analysis of possible error sources, using satellite products of snow depth and sea ice thickness, shows that the positive surface temperature errors during radiatively clear events are, to a large extent, caused by insufficient sea ice thickness and snow depth representation in the reanalysis system. A positive bias characterizes regions with ice thickness greater than 1.5 m, while the negative bias for thinner ice is partly compensated by the effect of snow.

54 ENVIRONMENTAL SCIENCES↗

Modeling Solar Lyman Alpha Irradiance

Solar Lyman alpha irradiance is estimated from various solar indices using linear regression analyses. Models developed with multiple linear regression analysis, including daily values and 81-day running means of solar indices, predict reasonably well both the short- and long-term variations observed in Lyman alpha. It is shown that the full disk equivalent width of the He line at 1083 nm offers the best proxy for Lyman alpha, and that the total irradiance corrected for sunspot effect also has a high correlation with Lyman alpha.

Pap, J.↗

Updated global and regional trends of stratospheric ozone profiles

We present updated evaluation of stratospheric ozone profile trends in the 60° S–60°N latitude range using long-term ground-based and satellite climate data records, as well as simulations by chemistry-climate models. The trends are evaluated using the LOTUS (Long-term Ozone Trends and Uncertainties in the Stratosphere) regression model. Analyses of satellite data confirm the statistically significant positive ozone trends in the period 2000–2024 in the upper stratosphere of ~1–3% per decade, with larger trends at mid-latitudes compared to the tropics. The trends are slightly positive or close to zero in the middle stratosphere, and mostly negative, -1 to -2% per decade, in the lower stratosphere, but they are not statistically significant. The morphology and magnitude of ozone trends are similar to previous analyses (2000–2020 trends). Ozone trends in 2000–2024 predicted by chemistry-climate model simulations are in good agreement with combined satellite trends. In the upper stratosphere, models predict a slightly stronger ozone recovery than observations. In the lower stratosphere, both models and satellite observations report negative trends in the tropics, while modelled ozone trends are slightly positive at mid-latitudes. Ozone profile trends over several stations estimated from ground-based records capture the same overall vertical pattern of ozone trends as merged gridded satellite datasets. Analyses of regional ozone profile trends in 2003–2024 using merged satellite datasets confirmed the previous observations of a longitudinal structure in ozone trends in the NH mid-latitude stratosphere, with positive trends over Scandinavia and negative trends over Siberia. However, the magnitude of this dipole-like structure is reduced compared to previous analyses.

trends↗