Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 541 records · Page 30

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗

Decadal Variability of Southeast US Rainfall in an Eddying Global Coupled Model

Abstract Ocean variability is a dominant source of remote rainfall predictability, but in many cases the physical mechanisms driving this predictability are not fully understood. This study examines how ocean mesoscales (i.e., the Gulf Stream SST front) affect decadal Southeast US (SEUS) rainfall, arguing that the local imprint of large‐scale teleconnections is sensitive to resolved mesoscale features. Based on global coupled model experiments with eddying and eddy‐parameterizing ocean, we find that a resolved Gulf Stream improves localized rainfall and remote circulation response in the SEUS. The eddying model generally improves the air‐sea interactions in the Gulf Stream and the North Atlantic Subtropical High that modulate SEUS rainfall over decadal timescales. The eddy‐parameterizing simulation fails to capture the sharp SST gradient associated with the Gulf Stream and overestimates the role of tropical Pacific SST anomalies in the SEUS rainfall.

54 ENVIRONMENTAL SCIENCES↗

Baroclinic Control of Southern Ocean Eddy Upwelling Near Topography

In the Southern Ocean, mesoscale eddies contribute to the upwelling of deep waters along sloping isopycnals, helping to close the upper branch of the meridional overturning circulation. Eddy energy (EE) is not uniformly distributed along the Antarctic Circumpolar Current (ACC). Instead, “hotspots” of EE that are associated with enhanced eddy-induced upwelling exist downstream of topographic features. This study shows that, in idealized eddy-resolved simulations, a topographic feature in the ACC path can enhance and localize eddy-induced upwelling. However, the upwelling systematically occurs in regions where eddies grow through baroclinic instability, rather than in regions where EE is large. Across a range of parameters, along-stream eddy growth rate is a more reliable indicator of eddy upwelling than traditional parameterizations such as eddy kinetic energy, eddy potential energy, or isopycnal slope. Ocean eddy parameterizations should consider metrics specific to the growth of baroclinic instability to accurately model eddy upwelling near topography.

54 ENVIRONMENTAL SCIENCES↗

Mesoscale Convective Systems Simulated by a High–Resolution Global Nonhydrostatic Model Over the United States and China

Mesoscale convective systems (MCSs) contribute a large fraction of warm-season precipitation and generate hazardous weather with substantial socio-economic impacts. Uncertainties in convection parameterizations in climate models limit our understanding of MCS characteristics and reliability of future projection. We examine the MCS simulation from the global 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) without cumulus parameterization against satellite observation from Global Precipitation Measurement (GPM) during 2001-2008. We focus on MCSs over the central U.S. and eastern China where MCSs prevail from March to August. A process-oriented tracking method incorporating both cloud and precipitation criteria is used to identify and track MCSs. About 140/100 MCSs initiate in the central U.S./eastern China per warm season and most of them initiate east of high mountains and in coastal regions. The frequency distribution of MCS lifetime is well captured in NICAM. But the simulated MCSs have stronger precipitation, smaller precipitation area, and larger cold cloud system than observed in both regions, which may be caused by weak entrainment as it is not well resolved at 14 km resolution. Here, the simulated MCS number is also underestimated in summer. By examining the climatological and MCS large-scale environments, the significant underestimation of MCS number in summer over the central U.S. may be attributed to the large climatological dry bias in the atmosphere. For China, mean moisture in summer is well simulated but deficiency in capturing the dynamic condition related to the coastal topography for triggering convection may have contributed to underestimation of MCS even in a sufficiently moist environment.

54 ENVIRONMENTAL SCIENCES↗

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗

Comparing Model Representations of Physiological Limits on Transpiration at a Semi-arid Ponderosa Pine Site

Mechanistic representations of biogeochemical processes in ecosystem models are rapidly advancing, requiring advancements in model evaluation approaches. Here we quantify multiple aspects of model functional performance to evaluate improved process representations in ecosystem models. We compare semi-empirical stomatal models with hydraulic constraints against more mechanistic representations of stomatal and hydraulic functioning at a semi-arid pine site using a suite of metrics and analytical tools. We find that models generally perform similarly under unstressed conditions, but performance diverges under atmospheric and soil drought. The more empirical models better capture synergistic information flows between soil water potential and vapor pressure deficit to transpiration, while the more mechanistic models are overly deterministic. Although models can be parameterized to yield similar functional performance, alternate parameterizations could not overcome structural model constraints that underestimate the unique information contained in soil water potential about transpiration. Additionally, both multilayer canopy and big-leaf models were unable to capture the magnitude of canopy temperature divergence from air temperature, and we demonstrate that errors in leaf temperature can propagate to considerable error in simulated transpiration. This study demonstrates the value of merging underutilized observational data streams with emerging analytical tools to characterize ecosystem function and discriminate among model process representations.

54 ENVIRONMENTAL SCIENCES↗

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗

Simulating Southern Ocean Aerosol and Ice Nucleating Particles in the Community Earth System Model Version 2

Southern Ocean (SO) low-level mixed phase clouds have been a long-standing challenge for Earth system models to accurately represent. While improvements to the Community Earth System Model version 2 (CESM2) resulted in increased supercooled liquid in SO clouds and improved model radiative biases, simulated SO clouds in CESM2 now contain too little ice. Previous observational studies have indicated that marine particles are major contributor to SO low-level cloud heterogeneous ice nucleation, a process that initiates a number of cloud processes that govern cloud radiative properties. In this study, we utilize detailed aerosol and ice nucleating particle (INP) measurements from two recent measurement campaigns to assess simulated aerosol abundance, number size distributions, and composition and INP parameterizations for use in CESM2. Our results indicate that CESM2 has a positive bias in simulated surface-level total aerosol surface area at latitudes north of 58°S. Measured INP populations were dominated by marine INPs and we present evidence of refractory INPs present over the SO assumed here to be mineral dust INPs. Results highlight a critical need to assess simulated mineral dust number and size distributions in CESM2 in order to adequately represent SO INP populations and their response to long-term changes in atmospheric transport patterns and land use change. Furthermore, we also discuss important cautions and limitations in applying a commonly used mineral dust INP parameterization to remote regions like the pristine SO.

58 GEOSCIENCES↗

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗

Bulk Transfer Coefficients Estimated From Eddy-Covariance Measurements Over Lakes and Reservoirs

The drag coefficient, Stanton number and Dalton number are of particular importance for estimating the surface turbulent fluxes of momentum, heat and water vapor using bulk parameterization. Although these bulk transfer coefficients have been extensively studied over the past several decades in marine and large-lake environments, there are no studies analyzing their variability for smaller lakes. Here, we evaluated these coefficients through directly measured surface fluxes using the eddy-covariance technique over more than 30 lakes and reservoirs of different sizes and depths. Our analysis showed that the transfer coefficients (adjusted to neutral atmospheric stability) were generally within the range reported in previous studies for large lakes and oceans. All transfer coefficients exhibit a substantial increase at low wind speeds (<3 m s -1 ), which was found to be associated with the presence of gusts and capillary waves (except Dalton number). Stanton number was found to be on average a factor of 1.3 higher than Dalton number, likely affecting the Bowen ratio method. At high wind speeds, the transfer coefficients remained relatively constant at values of 1.6·10 -3 , 1.4·10 -3 , 1.0·10 -3 , respectively. We found that the variability of the transfer coefficients among the lakes could be associated with lake surface area. In flux parameterizations at lake surfaces, it is recommended to consider variations in the drag coefficient and Stanton number due to wind gustiness and capillary wave roughness while Dalton number could be considered as constant at all wind speeds.

54 ENVIRONMENTAL SCIENCES↗

Contrasting the Biophysical and Radiative Effects of Rising CO 2 Concentrations on Ozone Dry Deposition Fluxes

Abstract The dry deposition of ozone from the atmosphere to ecosystems is an important coupling mechanism between atmospheric chemistry and terrestrial biogeochemical processes. In most Earth system models, dry deposition is simulated using a resistor‐in‐series approach that aims to parameterize the governing biological, chemical, and physical processes through a series of functional approximations. Here, we evaluate the influence of carbon cycle‐climate responses on this parameterization using the results of the Energy Exascale Earth System Model v1.1 Biogeochemistry simulation campaign. This simulation campaign was designed in part to explore the biophysical and radiative effects of rising historical CO 2 concentrations on the Earth system. We find that while the global annual ozone dry deposition is relatively insensitive to these effects, regionally the influence can be up to 10%. The strongest regional sensitivities in ozone dry deposition are predominantly in higher latitudes over land in the northern hemisphere and are dominated by the radiative effect of CO 2 , with little net influence of biophysical responses. Of all the impacts of the radiative effect of CO 2 , we point to the potential importance of accurately representing ozone deposition to snow in Earth System Models and provide recommendations for future simulation campaigns.

54 ENVIRONMENTAL SCIENCES↗

On the Role of Macrophysics and Microphysics in Km-Scale Simulations of Mixed-Phase Clouds During Cold Air Outbreaks

Regional atmospheric models struggle to maintain supercooled liquid in mixed-phase clouds during polar cold-air outbreaks (CAOs). Previous studies focused on the parameterization of aerosol, microphysics and turbulence to understand the origin of this widespread model bias. This study investigates the role of sub-grid cloud fraction parameterizations (CFP) in the simulation of mixed-phase cloud. Km-scale simulations are performed for a large number of CAO cases over Norway, for which continuous ground observations were collected at one site over 6 months. We use a novel analysis that attributes the cloud-radiative errors to deficiencies in specific cloud regimes. We show that the CFP matters for cloud-radiative effects in CAOs, but that it is probably not the primary cause of the lack of liquid water in simulated mixed-phase clouds. Of all the CFP sensitivities explored in this study, the prognostic representation of both liquid and ice shows most promise in increasing the liquid water path. A newly proposed hybrid CFP with prognostic frozen and diagnostic liquid cloud fraction reproduces some of the benefits of the prognostic scheme at reduced cost and complexity. All simulations appear to produce too large precipitation particles. As a result, reducing the snow deposition rate increases the liquid water path and improves the precipitation size distributions, while simulations are less sensitive to reduced riming rates.

54 ENVIRONMENTAL SCIENCES↗

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Differentiable, Learnable, Regionalized Process-Based Models With Multiphysical Outputs can Approach State-Of-The-Art Hydrologic Prediction Accuracy

Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Winter Heat Conduction Through the Sea Ice System During MOSAiC

Abstract Models struggle to accurately simulate observed sea ice thickness changes, which could be partially due to inadequate representation of thermodynamic processes. We analyzed co‐located winter observations of the Arctic sea ice from the Multidisciplinary Drifting Observatory for the Study of the Arctic Climate for evaluating and improving thermodynamic processes in sea ice models, aiming to enable more accurate predictions of the warming climate system. We model the sea ice and snow heat conduction for observed transects forced by realistic boundary conditions to understand the impact of the non‐resolved meter‐scale snow and sea ice thickness heterogeneity on horizontal heat conduction. Neglecting horizontal processes causes underestimating the conductive heat flux of 10% or more. Furthermore, comparing model results to independent temperature observations reveals a ∼5 K surface temperature overestimation over ice thinner than 1 m, attributed to shortcomings in parameterizing surface turbulent and radiative fluxes rather than the conduction. Assessing the model deficiencies and parameterizing these unresolved processes is required for improved sea ice representation.

Geology↗

Basal Melting and Oceanic Observations Beneath Central Fimbulisen, East Antarctica

Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year‐round melt record (2017–2021) from a phase‐sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long‐term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr −1 , marked by substantial sub‐weekly variability ranging from 0.4 to 3.5 m yr −1 . 36‐h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September–March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability ( r = 0.84), highlighting the dominant role of shear‐driven turbulence. This relationship enables parameterization of melt rates for the decade‐long ocean record (2010–2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite‐derived melt rates but with a tenfold lower seasonal amplitude and a 3‐month delay in seasonality. These detailed concurrent ice–ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice‐shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice‐shelf stability and sea‐level change.

54 ENVIRONMENTAL SCIENCES↗