Search NASA⌕ Search

SEARCH · Search NASA

Results for “evaporative fraction regimes”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A simple framework to characterize land aridity based on surface energy partitioning regimes

Land aridity is often characterized by the aridity index (AI), which does not account for land surface water-energy interactions that are crucially important in determining regional climate. Such interactions can be captured by the evaporative fraction (EF, ratio of evapotranspiration to available energy) regimes. As EF is subject to energy and water limitations in humid and dry areas, respectively, EF regimes may be used to characterize land aridity to account for the influence of complex land characteristics and their impact on water availability. Here, we propose a simple framework to characterize land aridity by statistically ranking the coupling strength between EF and surface energy and water terms. The framework is demonstrated using gridded data and compared with AI over the U.S. and China. Results show that regionalization of aridity zones based on EF regimes and a two-tiered classification scheme may provide information such as surface energy and water variability complementary to the background aridity depicted by AI, with implications for extreme events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Direct numerical simulation of turbulent dispersion of evaporative aerosol clouds produced by an intense expiratory event

Airborne particles are a major route for transmission of COVID-19 and many other infectious diseases. When a person talks, sings, coughs, or sneezes, nasal and throat secretions are spewed into the air. After a short initial fragmentation stage, the expelled material is mostly composed of spherical particles of different sizes. While the dynamics of the largest droplets are dominated by gravitational effects, the smaller aerosol particles, mostly transported by means of hydrodynamic drag, form clouds that can remain afloat for long times. In subsaturated air environments, the dependence of pathogen-laden particle dispersion on their size is complicated due to evaporation of the aqueous fraction. Particle dynamics can significantly change when ambient conditions favor rapid evaporation rates that result in a transition from buoyancyto-drag dominated dispersion regimes. To investigate the effect of particle size and evaporation on pathogen-laden cloud evolution, a direct numerical simulation of a mild cough was coupled with an evaporative Lagrangian particle advection model. The results suggest that while the dispersion of cough particles in the tails of the size distribution are unlikely to be disrupted by evaporative effects, preferential aerosol diameters (30-40 μm) may exhibit significant increases in the residence time and horizontal range under typical ambient conditions. Using estimations of the viral concentration in the spewed fluid and the number of ejected particles in a typical respiratory event, we obtained a map of viral load per volume of air at the end of the cough and the number of virus copies per inhalation in the emitter vicinity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Land–Atmosphere Coupling at the U.S. Southern Great Plains: A Comparison on Local Convective Regimes between ARM Observations, Reanalysis, and Climate Model Simulations

Using the 9-yr warm-season observations at the Atmospheric Radiation Measurement Southern Great Plains site, in this work we assess the land–atmosphere (LA) coupling in the North American Regional Reanalysis (NARR) and two climate models: hindcasts with the Community Atmosphere Model version 5.1 by Cloud-Associated Parameterizations Testbed (CAM5-CAPT) and nudged runs with the Energy Exascale Earth System Model Atmosphere Model version 1 Regionally Refined Model (EAMv1-RRM). We focus on three local convective regimes and diagnose model behaviors using the local coupling metrics. NARR agrees well with observations except a slightly warmer and drier surface with higher downwelling shortwave radiation and lower evaporative fraction. On clear-sky days, it shows warmer and drier early-morning conditions in both models with significant underestimates in surface evaporation by EAMv1-RRM. On the majority of the ARM-observed shallow cumulus days, there is no or little low-level clouds in either model. When captured in models, the simulated shallow cumulus shows much less cloud fraction and lower cloud bases than observed. On the days with late-afternoon deep convection, models tend to present a stable early-morning lower atmosphere more frequently than the observations, suggesting that the deep convection is triggered more often by elevated instabilities. Generally, CAM5-CAPT can reproduce the local LA coupling processes to some extent due to the constrained early-morning conditions and large-scale winds. EAMv1-RRM exhibits large precipitation deficits and warm and dry biases toward mid-to-late summers, which may be an amplification through a positive LA feedback among initial atmosphere and land states, convection triggering and large-scale circulations.

58 GEOSCIENCES↗
DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Effects of Wind Shear and Aerosol Conditions on the Organization of Precipitating Marine Stratocumulus Clouds

Here this study examines how wind shear affects precipitating marine stratocumulus clouds under different cloud droplet number concentrations (N d ). We performed a series of large eddy simulations (LES) of nocturnal marine stratocumulus clouds using Cloud Model 1 (CM1). The simulations show that N d is the dominant factor for cloud cellular organization in this cloud regime rather than wind shear. Low N d characterizes the open cellular structure with a high in-cloud liquid water path (LWP). When wind shear is increased, the cloud fraction tends to decrease along with LWP, suggesting the cloud top region is significantly influenced by the entrainment and mixing of dry air from the free troposphere. We also examine cold pools in open and closed cellular clouds. Open-cell clouds produce larger and deeper cold pools compared to closed-cell clouds. Interestingly, cold pools can exist without surface precipitation and are produced by evaporation of light precipitation (drizzle) below the cloud base with weak downdrafts. The evaporation of raindrops and drizzle play an important role in initiating new convection, particularly where colliding outflows occur downstream of the cloud. This secondary convection contributes to the development and maintenance of the cloud cellular organization and formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimates of global recycling coefficients for LTX- β discharges

We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Microphysics Dominates Sub‐Cloud Rain Evaporation in Trade Cumuli Over Barbados

This study investigates the microphysical and thermodynamic influences on North Atlantic trade-cumulus rain evaporation in the sub-cloud layer over Barbados from January to November 2020. Using radar observations and a 1-D model, the results reveal that microphysical properties, namely geometric mean diameter (D g ) and raindrop concentration, primarily control evaporation processes under conditions of limited thermodynamic variability. Thermodynamic factors like cloud base height and relative humidity exert a smaller influence due to their narrow observed ranges. Small D g leads to rapid increases in rain evaporation fraction (REF), as small drops evaporate completely, while larger D g yields slower, more gradual REF changes. Frequent weak-rain cases reduce the mean evaporation flux, whereas infrequent intense-rain events substantially enhance it. These findings highlight the pivotal role of microphysical variability in shaping sub-cloud evaporation in this trade cumulus regime, while thermodynamic effects may be more important in other environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mixed-state entanglement and information recovery in thermalized states and evaporating black holes

We study the universal behavior of quantum information-theoretic quantities in thermalized isolated quantum many-body systems and evaporating black holes. In particular, we study a genuine mixed-state entanglement measure called the logarithmic negativity, other correlation measures including the Renyi negativities and the mutual information, and a signature of multipartite entanglement called the reflected entropy. We also probe the feasibility of recovering quantum information from subsystems of a thermalized quantum many-body system or from the radiation of an evaporating black hole, using quantities such as relative entropy and Petz map fidelity. A recently developed technique called the equilibrium approximation allows us to probe these quantities at finite temperature. We find striking qualitative differences from the infinite temperature case, which has been the topic of previous studies using Haar-random states. In particular, we find regimes where the logarithmic negativity is extensive but the mutual information is sub-extensive, indicating a large amount of undistillable, bound entanglement in thermalized states. For evaporating black holes at finite temperature, both the logarithmic negativity and the Petz map fidelity reveal an important new time scale t b , which is earlier than the Page time t p by a finite fraction of the total evaporation time. We find that t b , as opposed to t p , is the time scale at which quantum entanglement between different parts of the radiation becomes extensive, and the fidelity of information recovery for a large diary thrown into the black hole starts to grow.

2D Gravity↗
DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗
DOE OSTI2022

Disentangling the Relative Drivers of Seasonal Evapotranspiration Across a Continental-Scale Aridity Gradient

Evapotranspiration (ET) is a significant ecosystem flux, governing the partitioning of energy at the land surface. Understanding the seasonal pattern and magnitude of ET is critical for anticipating a range of ecosystem impacts, including drought, heat-wave events, and plant mortality. Here, in this study, we identified the relative controls of seasonal variability in ET, and how these controls vary among ecosystems. We used overlapping AmeriFlux and PhenoCam time series at a daily timestep from 20 sites to explore these linkages (# site-years >100), and our study area covered a broad climatological aridity gradient in the U.S. and Canada. We focused on disentangling the most important controls of bulk surface conductance (G s ) and evaporative fraction (EF = LE/[H + LE]), where LE and H represent latent and sensible heat fluxes, respectively. Specifically, we investigated how vegetation phenology varied in importance relative to meteorological variables (vapor pressure deficit and antecedent precipitation) as a driver of G s and EF using path analysis, a framework for quantifying and comparing the causal linkages among multiple response and explanatory variables. Our results revealed that the drivers of G s and EF seasonality varied significantly between energy- and water-limited ecosystems. Specifically, precipitation had a much higher effect in water-limited ecosystems, while seasonal patterns in canopy greenness emerged as a stronger control in energy-limited ecosystems. Given that phenology is expected to shift under future climate, our findings provide key information for understanding and predicting how phenology may impact 21st-century hydroclimate regimes and the surface-energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗