Search NASASearch

SEARCH · Search NASA

Results for “model parameterization”

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

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals

Evaluating mesoscale model predictions of diurnal speedup events in the Altamont Pass Wind Resource Area of California

Mesoscale model predictions of wind, turbulence, and wind energy capacity factors are evaluated in the Altamont Pass Wind Resource Area of California (APWRA), where the diurnal regional sea breeze and associated terrain-driven speedup flows drive wind energy production during the summer months. Results from the Weather Research and Forecasting model version 4.4 using a novel three-dimensional planetary boundary layer (3D PBL) scheme, which treats both vertical and horizontal turbulent mixing, are compared to those using a well-established one-dimensional (1D) scheme that treats only vertical turbulent mixing. Each configuration is evaluated over a nearly 3-month-long period during the Hill Flow Study, and due to the recurring nature of the observed speedup flows, diurnal composite averaging is used to capture robust trends in model performance. Both model configurations showed similar overall skill. The general timing and direction of the speedup flows is captured, but their magnitude is overestimated within a typical wind turbine rotor layer. Both also fail to capture a persistent observed near-surface jet-like flow, likely due to the limited grid resolution that is typical of mesoscale models. However, the 3D PBL configuration shows several minor improvements over the 1D PBL configuration, including improved wind speed and turbulence kinetic energy profiles during the accelerating phase of the speedup events, as well as reduced positive wind speed bias at surface stations across the APWRA region. Using a mesoscale wind farm parameterization, modeled capacity factors are also compared to monthly data reported to the US Energy Information Administration (EIA) during the study period. Although the monthly trend in the data is captured, both model configurations overestimate capacity factors by roughly 7 %–11 %. Through model evaluation, this study provides confidence in the 3D PBL scheme for wind energy applications in complex terrain and provides guidance for future testing.

17 WIND ENERGY

Diskoseismology - Signatures of black hole accretion disks

General relativity requires the existence of a spectrum of oscillations which are trapped near the inner edge of accretion disks around black holes. We have developed a general formalism for analyzing the normal modes of such acoustic perturbations of arbitrary thin disk models, approximating the dominant relativistic effects via a modified Newtonian potential (these modes do not exist in Newtonian gravity). The eigenfunctions and eigenfrequencies of a variety of disk models are found to fall in to two main classes, which are analogous to the p-modes and g-modes in the sun. In this work, we compute the eigenfunctions and eigenfrequencies of isothermal disks. The (relatively small) rates of growth or damping of these oscillations due to gravitational radiation and parameterized models of viscosity are also computed.

Nowak, Michael

Global High Resolution Crustal Magnetic Field at the Surface of the Moon from Low-Altitude Lunar Prospector Magnetic Gradient Data

We derived new vector gradients based models of crustal magnetic field at the lunar surface with data from the Lunar Prospector (LP) satellite using two model parameterization approaches: a global set of 35820 1° spaced (~30 km) equal area monopoles at 20 km below the surface (O’Brien and Parker, 1994; Olsen et al., 2017) and combined results of subsets of 100000 0.66° spaced monopoles at the same depth. We use the scheme of iteratively reweighted least-squares inversion to compute the initial model. Then the amplitudes of these monopoles are determined by minimizing the misfit to the components together with the global average of |Br| at the ellipsoid surface (i.e. applying a L1 model regularization of Br). In previous approaches using vector fields for modeling, we found that external field contamination leads to spurious anomalies in the downward continued field models even with stringent data selection criteria and ad-hoc noise removal techniques (e.g., satellite’s position in the Moon’s wake w.r.t. the solar wind and in the Earth’s magnetotail, internal/external dipoles fields removal, low-order polynomial removal, joint equivalent source cross-validation technique and visually removing remaining anomalous segments). On the other hand, with the use of gradients-only data (along-track first finite differences), we were able to completely bypass the ad-hoc techniques. Similar processing of Kaguya magnetic data, which have only higher altitude coverage over the most of the Moon except in the region of the South Pole–Aitken (SPA) basin, completely misses some of the anomalies seen in the Lunar Prospector data. The combined Lunar Prospector and Kaguya gradient-based models also severely degrade the derivation of the anomaly fields in many regions. Euler analysis of isolated anomaly features from the Reiner Gamma swirl suggests top depths of about 0.3 to 1.5 km and center depths of 10-14 km; in the region between Stein and Vallier craters north of the SPA basin our analysis suggests top depths of around 1.5 km and center depths of 13-15 km. With the spectral depth determination techniques, the SPA basin region yields depths to the base of magnetization ranging between 15 and 40 km. Three-dimensional modeling and the bulk magnetization determinations of the sources constrained by the Euler and spectral methods is underway.

Gradient

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES

Examining the Impacts of High-Resolution Land Surface Initialization on Model Predictions of Convection in the Southeastern U.S.

One of the most challenging weather forecast problems in the southeastern U.S. is daily summertime pulse convection. During the summer, atmospheric flow and forcing are generally weak in this region; thus, convection typically initiates in response to local forcing along sea/lake breezes, and other discontinuities often related to horizontal gradients in surface heating rates. Numerical simulations of pulse convection usually have low skill, even in local predictions at high resolution, due to the inherent chaotic nature of these precipitation systems. Forecast errors can arise from assumptions within physics parameterizations, model resolution limitations, as well as uncertainties in both the initial state of the atmosphere and land surface variables such as soil moisture and temperature. For this study, it is hypothesized that high-resolution, consistent representations of surface properties such as soil moisture and temperature, ground fluxes, and vegetation are necessary to better simulate the interactions between the land surface and atmosphere, and ultimately improve predictions of local circulations and summertime pulse convection. The NASA Short-term Prediction Research and Transition (SPORT) Center has been conducting studies to examine the impacts of high-resolution land surface initialization data generated by offline simulations of the NASA Land Informatiot~ System (LIS) on subsequent numerical forecasts using the Weather Research and Forecasting (WRF) model (Case et al. 2008, to appear in the Journal of Hydrometeorology). Case et al. presents improvements to simulated sea breezes and surface verification statistics over Florida by initializing WRF with land surface variables from an offline LIS spin-up run, conducted on the exact WRF domain and resolution. The current project extends the previous work over Florida, focusing on selected case studies of typical pulse convection over the southeastern U.S., with an emphasis on improving local short-term WRF simulations over the Mobile, AL and Miami, FL NWS county warning areas. Future efforts may involve examining the impacts of assimilating remotely-sensed soil moisture data, and/or introducing weekly greenness vegetation fraction composites (as opposed to monthly climatologies) into ol'fline NASA LIS runs. Based on positive impacts, the offline LIS runs could be transitioned into an operational mode, providing land surface initialization data to NWS forecast offices in real time.

Case, Jonathan L.

NASA Downscaling Project: Final Report

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA - 2 reanalyses were used to drive the NU - WRF regional climate model and a GEOS - 5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000 - 2010. The results of these experiments were compared to observational datasets to evaluate the output.

dynamical downsizing

Superclustering in the explosion scenario

A simple toy model is used to study the spatial distribution of rich clusters in a generic type of explosion scenario. The model, parameterized by the distribution of shell radii and the filling factor, places spherical shells at random and identifies each 'knot' as a cluster. The resulting cluster correlation function is close to a power law extending to the diameter of the largest spheres. Richer clusters form at the intersections of bigger shells and so have stronger correlations. Typical shell radii and filling factors are required to produce the observed number density of clusters. Models with a power-law radius distribution also reproduce the richness distribution of clusters in the Abell catalog. Supercluster multiplicity functions, void probabilities, number counts, topology statistics, and velocity correlations confirm the presence of strong superclustering and quantify the non-Gaussian nature of the model.

Weinberg, David H.

Potential Link Between Ice Nucleation and Climate Model Spread in Arctic Amplification

Arctic amplification (AA) is simulated by all global climate models, however the spread in the degree of projected warming is large and the underlying mechanisms driving it are poorly understood. The impact of the temperature dependence of immersion freezing on cloud feedbacks and AA is studied using NASA's GEOS-5 model. Parameterizations that exhibit low ice-nucleating particle (INP) concentrations in the high Arctic during summer are found to weaken the cloud-phase feedback. This allows sunlight to readily melt sea-ice in the summer, which decreases the stability of the lower troposphere, causing a decrease in wintertime cloud fraction over open ocean. Arctic amplification was found to span from ∼1.4 to >2.6, which spans 30% of the spread in AA in the coupled model intercomparison project models, depending on the temperature dependence of immersion freezing. These results suggest that summertime INP concentrations may provide an observational constraint on AA.

Ice Nucleation

Implications of using MERRA-2 cloud water in aerosol wet scavenging for Pb-210 surface concentrations and vertical profiles in GEOS-Chem

Precipitation scavenging is the dominant loss process for a whole suite of aerosols but model parameterizations of this process are highly uncertain, substantially contributing to large uncertainties in the simulated loadings and radiative forcing of aerosols. Lead-210 (210Pb, radioactive half-life of 22.3 years) is produced by radioactive decay of soil-emitted gaseous 222Rn. It attaches to ambient submicron aerosols and is subject to precipitation scavenging processes. Liu et al. [2001] estimated the global mean lifetime of tropospheric 210Pb aerosols to be ~9 days using the GEOS-Chem model. Luo et al. [2019, 2020] provided an option to use MERRA-2 cloud water in the model aerosol scavenging scheme. In this study, we examine the impact of this revision on 210Pb simulations and compare results with climatological observations of 210Pb surface and UT/LS concentrations, as well as aircraft profile measurements from eleven NASA field missions.

Bo Zhang

NASA/GSFC Research Activities for the Global Ocean Carbon Cycle: A Prospectus for the 21st Century

There are increasing concerns that anthropogenic inputs of carbon dioxide into the Earth system have the potential for climate change. In response to these concerns, the GSFC Laboratory for Hydrospheric Processes has formed the Ocean Carbon Science Team (OCST) to contribute to greater understanding of the global ocean carbon cycle. The overall goals of the OCST are to: 1) detect changes in biological components of the ocean carbon cycle through remote sensing of biooptical properties, 2) refine understanding of ocean carbon uptake and sequestration through application of basic research results, new satellite algorithms, and improved model parameterizations, 3) develop and implement new sensors providing critical missing environmental information related to the oceanic carbon cycle and the flux of CO2 across the air-sea interface. The specific objectives of the OCST are to: 1) establish a 20-year time series of ocean color, 2) develop new remote sensing technologies, 3) validate ocean remote sensing observations, 4) conduct ocean carbon cycle scientific investigations directly related to remote sensing data, emphasizing physiological, empirical and coupled physical/biological models, satellite algorithm development and improvement, and analysis of satellite data sets. These research and mission objectives are intended to improve our understanding of global ocean carbon cycling and contribute to national goals by maximizing the use of remote sensing data.

Gregg, W. W.

FIRE Arctic Clouds Experiment

An overview is given of the First ISCCP Regional Experiment (FIRE) Arctic Clouds Experiment that was conducted in the Arctic during April through July, 1998. The principal goal of the field experiment was to gather the data needed to examine the impact of arctic clouds on the radiation exchange between the surface, atmosphere, and space, and to study how the surface influences the evolution of boundary layer clouds. The observations will be used to evaluate and improve climate model parameterizations of cloud and radiation processes, satellite remote sensing of cloud and surface characteristics, and understanding of cloud-radiation feedbacks in the Arctic. The experiment utilized four research aircraft that flew over surface-based observational sites in the Arctic Ocean and Barrow, Alaska. In this paper we describe the programmatic and science objectives of the project, the experimental design (including research platforms and instrumentation), conditions that were encountered during the field experiment, and some highlights of preliminary observations, modelling, and satellite remote sensing studies.

Curry, J. A.

Regional-Scale Modeling at NASA Goddard Space Flight Center

Over the past decade, the Goddard Mesoscale Modeling and Dynamics Group has used a popular regional scale model, MM5, to study precipitation processes. Our group is making contributions to the MM5 by incorporating the following physical and numerical packages: improved Goddard cloud processes, a land processes model (Parameterization for Land-Atmosphere-Cloud Exchange - PLACE), efficient but sophisticated radiative processes, conservation of hydrometeor mass (water budget), four-dimensional data assimilation for rainfall, and better computational methods for trace gas transport. At NASA Goddard, the MM5 has been used to study: (1) the impact of initial conditions, assimilation of satellite-derived rainfall, and cumulus parameterizations on rapidly intensifying oceanic cyclones, hurricanes and typhoons, (2) the dynamic and thermodynamic processes associated with the development of narrow cold frontal rainbands, (3) regional climate and water cycles, (4) the impact of vertical transport by clouds and lightning on trace gas distributiodproduction associated with South and North American mesoscale convective systems, (5) the development of a westerly wind burst (WWB) that occurred during the TOGA COARE and the diurnal variation of precipitation in the tropics, (6) a Florida sea breeze convective event and a Mid-US flood event using a sophisticated land surface model, (7) the influence of soil heterogeneity on land surface energy balance in the southwest GCIP region, (8) explicit simulations (with 1.33 to 4 km horizontal resolution) of hurricanes Bob (1991) and Bonnie (1998), (9) a heavy precipitation event over Taiwan, and (10) to make real time forecasts for a major NASA field program. In this paper, the modifications and simulated cases will be described and discussed.

Tao, W.-K.

Risk analysis: divergent models and convergent interpretations

Material presented at a NASA-sponsored workshop on risk models for exposure conditions relevant to prolonged space flight are described in this paper. Analyses used mortality data from experiments conducted at Argonne National Laboratory on the long-term effects of external whole-body irradiation on B6CF1 mice by 60Co gamma rays and fission neutrons delivered as a single exposure or protracted over either 24 or 60 once-weekly exposures. The maximum dose considered was restricted to 1 Gy for neutrons and 10 Gy for gamma rays. Proportional hazard models were used to investigate the shape of the dose response at these lower doses for deaths caused by solid-tissue tumors and tumors of either connective or epithelial tissue origin. For protracted exposures, a significant mortality effect was detected at a neutron dose of 14 cGy and a gamma-ray dose of 3 Gy. For single exposures, radiation-induced mortality for neutrons also occurred within the range of 10-20 cGy, but dropped to 86 cGy for gamma rays. Plots of risk relative to control estimated for each observed dose gave a visual impression of nonlinearity for both neutrons and gamma rays. At least for solid-tissue tumors, male and female mortality was nearly identical for gamma-ray exposures, but mortality risks for females were higher than for males for neutron exposures. As expected, protracting the gamma-ray dose reduced mortality risks. Although curvature consistent with that observed visually could be detected by a model parameterized to detect curvature, a relative risk term containing only a simple term for total dose was usually sufficient to describe the dose response. Although detectable mortality for the three pathology end points considered typically occurred at the same level of dose, the highest risks were almost always associated with deaths caused by tumors of epithelial tissue origin.

Non-NASA Center

Modeling and Evaluation of Canted Coil Springs as High Temperature Seal Preloading Devices

Future reusable launch vehicles will require advanced structural seals. This includes propulsion seals along edges and hinge lines in hypersonic engines, and control surface seals for movable flaps and elevons on proposed reentry vehicles. Seals must remain in sealing engagement with opposing surfaces, for multiple missions, even though the seal gap may be opening and closing due to thermal and structural loads. To meet this requirement either the seals themselves must be resilient or there must be a resilient structural element behind the seals. Case Western Reserve University is working with NASA s Glenn Research Center to develop more resilient high temperature seal components and preloading devices. Results are presented for a finite element analysis of a canted coil spring that is being considered as a high temperature seal preloading device. This type of spring is a leading candidate due to its ability to provide nearly constant force over a large deflection. The finite element analyses were verified by comparing them to experimental results of canted coil springs of three different stiffnesses, measured at Glenn Research Center. Once validated the parameterized model was combined with a scripting algorithm to assess the effects of key spring design variables (wire diameter, coils per inch, cant amplitude, eccentricity, and spring width) on spring stiffness and maximum Von Mises stress to aid in subsequent design.

Oswald, Jay J.

Evaluation of the Sensitivity of the Amazonian Diurnal Cycle to Convective Intensity in Reanalyses

Model parameterizations of tropical deep convection are unable to reproduce the observed diurnal and spatial variability of convection in the Amazon, which contributes to climatological biases in the water cycle and energy budget. Convective intensity regimes are defined using percentiles of daily minimum 3-hourly averaged outgoing longwave radiation (OLR) from Clouds and the Earth's Radiant Energy System (CERES). This study compares the observed spatial variability of convective diurnal cycle statistics for each regime to MERRA-2 and ERA-Interim (ERA) reanalysis data sets. Composite diurnal cycle statistics are computed for daytime hours (06:00-21:00 local time) in the wet season (December-January-February). MERRA-2 matches observations more closely than ERA for domain averaged composite diurnal statistics-specifically precipitation. However, ERA reproduces mesoscale features of OLR and precipitation phase associated with topography and the propagation of the coastal squall line. Both reanalysis models are shown to underestimate extreme convection.

Itterly, Kyle F.

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

Soil Moisture Active Passive (SMAP) mission brightness temperature (T(b) ) observations are assimilated into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update simulations of surface and root-zone soil moisture. Different time series components of the T(b) observations are assimilated including anomalies, inter-annual variations, and high frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the uncertainties of modeled and observed T(b) is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. In a benchmark experiment, T(b) observations are assimilated using a spatially constant 4 Kelvin (K) observation uncertainty, as in the operational SMAP Level-4 algorithm. All DA experiments exhibit notable skill improvements in most regions. Improvements are largest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for surface soil moisture simulations (0.06). No notable difference in skill - neither improvement nor deterioration - is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. This may be explained by the presence of large observation operator errors, which are analytically shown to have the potential to render post-update uncertainty insensitive to inaccuracies in estimates of the Kalman gain. These results have important implications for the design of soil moisture DA systems, in particular for parameterizing model and observation uncertainties.

Hydrology