Beyond a deterministic representation of the temperature dependence of soil respiration
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Direct normal irradiance (DNI) is often interpreted differently in ground measurements and forecasts using numerical weather prediction (NWP) models, leading to substantial bias in DNI forecasting especially under cloudy-sky conditions. To mitigate the bias, we use the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) to conduct physics-based simulations of solar radiation in the circumsolar region. The DNI is quantified using the sum of the radiation along the sun direction and the scattered radiation within the circumsolar region. FARMS-DNI is coupled with the Weather Research and Forecasting model with solar extensions (WRF-Solar) to predict day-ahead DNI in 9-km pixels over an extended domain covering the contiguous United States. The DNI forecasts for the entire year of 2018 are validated using surface-based observations at the Atmospheric Radiation Measurement (ARM) and Surface Radiation Budget Network (SURFRAD) sites and a satellite-based solar radiation product developed by the procedure of the National Solar Radiation Data Base (NSRDB). The results show that the WRF-Solar coupled with the conventional Beer-Bouguer-Lambert law underestimes the overall DNI forecasts by 30%-60%, which is a joint effect of the uncertainties in computing the circumsolar radiation and forecasting cloud properties. This bias is significantly reduced by using FARMS-DNI.
Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.
Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.
Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan−Wigner transformation (JW). The Hartree−Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > $\frac{1}{2}$.
Not provided.
In this study, we present a graph neural network-based learning approach using an autoencoder setup to derive low-dimensional variables from features observed in experimental crystal structures. These variables are then biased in enhanced sampling to observe state-to-state transitions and reliable thermodynamic weights. Our approach uses simple convolution and pooling methods. To verify the effectiveness of our protocol, we examined the nucleation of various allotropes and polymorphs of iron and glycine from their molten states. Our graph latent variables when biased in well-tempered metadynamics consistently show transitions between states and achieve accurate free energy calculations in agreement with experiments, both of which are indicators of dependable sampling. This underscores the strength and promise of our graph neural net variables for improved sampling. Furthermore, the protocol shown here should be applicable for other systems and with other sampling methods.
A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.
Abstract Days of extreme precipitation over California are evaluated in Coupled Model Intercomparison Project Phase 6 (CMIP6) models and the ERA‐Interim reanalysis. In the current climate, the model spread in composited precipitation on extreme precipitation days is closely related to the magnitude of composited integrated vapor transport (IVT) across models, a proxy for the intensity of atmospheric rivers. Most models underestimate the magnitude of IVT associated with extreme precipitation, according to ERA‐Interim. This is due mostly to the contribution of moisture, which almost all models overestimate, while the contribution of lower‐tropospheric wind speed is generally closer to the reanalyses. Moreover, most models underestimate the variance in the latitude of maxima of numerous variables among days of extreme California precipitation. That is, in the general circulation models there is a lack of diversity in the latitude of the disturbances bringing winter precipitation to California. In the future climate, most models project a decrease in the frequency of southward‐displaced disturbances among California extreme precipitation days. Hence, the greatest increases in extreme precipitation are over northern California. However, the historical underestimate of the latitudinal variance of disturbances calls into question the reliability of these projections. This bias should be especially considered for dynamical downscaling efforts over the region.
Here, we introduce a physics-based aerosol wet removal scheme with unified treatments of aerosol transport and removal by convective clouds into the Community Atmosphere Model version 6. Since several important physical processes are still neglected or poorly represented in this new physics-based scheme, we develop secondary improvements to the parameterizations of aerosol activation, resuspension, and cloud-borne aerosol detrainment in this new scheme. Changes in the aerosol wet removal scheme cause tropospheric aerosol concentrations to decrease to different extents: compared to the control run, the physics-based scheme significantly decreases aerosol burdens by up to 60% over the southern Pacific Ocean, whereas the secondary improvements mitigate the decreasing tendency. The burden changes also depend on aerosol chemical components: the sulfate mass decrease is compensated by secondary production, black carbon (BC) is effectively removed via increasing the hygroscopicity of particulate organic matter from 0 to 0.2, and dust shows the most spatially heterogeneous changes. Simulated aerosol profiles are evaluated against aircraft-based observations over the Pacific and Atlantic Oceans. The secondary-improved scheme reduces the overestimations of upper tropospheric BC and sea salt concentrations by a factor of 10 and 1,000, respectively, and reproduces the dependence of BC mass decrease rates on cloud types. Consideration of convective cloud-borne aerosol detrainment plays the most important role in enhancing the aerosol wet removal and decreasing the positive biases of tropospheric BC and sea salt concentrations. We also summarize unresolved issues related to convective cloud genesis and microphysics, cloud-borne aerosol evolution, and BC and dust emissions.
The warm Gulf Stream sea surface temperatures (SSTs) strongly impact the evolution of winter clouds behind atmospheric cold fronts. Such cloud evolution remains challenging to model. Here, the Gulf Stream is too wide within the ERA5 and MERRA2 reanalyses, affecting the turbulent surface fluxes. Known problems within the ERA5 boundary layer (too-dry and too-cool with too strong westerlies), ascertained primarily from ACTIVATE 2020 campaign aircraft dropsondes and secondarily from older buoy measurements, reinforce surface flux biases. In contrast, MERRA2 winter surface winds and air-sea temperature/humidity differences are slightly too weak, producing surface fluxes that are too low. Reanalyses boundary layer heights in the strongly-forced winter cold-air-outbreak regime are realistic, whereas late-summer quiescent stable boundary layers are too shallow. Nevertheless, the reanalysis biases are small, and reanalyses adequately support their use for initializing higher-resolution cloud process modeling studies of cold-air outbreaks.
Abstract The ability of an observationally‐constrained cloud‐system resolving model (Weather Research and Forecasting; WRF, 4‐km grid spacing) and a global climate model (Energy Exascale Earth System Model; E3SM, 1‐degree grid spacing) to represent the precipitation diurnal cycle over the Amazon basin during the 2014 wet season is assessed. The WRF model coupled with a 3‐D variational data assimilation scheme reproduces the spatial variability of the precipitation diurnal cycle over the basin and the lifecycle of westward propagating MCSs initiated by the coastal sea‐breeze front. In contrast, a single morning peak in rainfall is produced by E3SM for simulations despite the nudging of large‐scale winds toward global reanalysis, indicating precipitation in E3SM is largely controlled by local convection associated with diurnal heating. The role of propagating MCS on the environment are discussed by using a multivariate perturbation analysis. We also find that the advection of moisture perturbations from ocean to inland regions have a higher correlation with the occurrence of MCSs in the Amazon than the intensity of colder air intrusion associated with sea breezes along the coast. Moreover, the presence of large cold pools over the central Amazon basin are responsible for the maintenance of propagating deep convection.
Abstract Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land‐surface model (LSM) simulations during drought. However, systematic evaluations of both optimality‐based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf‐level observations of gas‐exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr −1 . We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf‐gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses.
Understanding the water balance of a catchment in relation to its regional climate forcings and catchment characteristics is critical for predicting current and future water resources amid changing climate and land cover. This study intends to improve Zhang's monthly water balance model (a physics-based conceptual hydrologic model) that reflects the physical partitioning process of the hydrological cycle at the basin level based on regional climate and catchment characteristics. The existing model does not include snow process and has confronted evident limitations in snow-affected areas, which is a critical aspect since snowmelt water has been a significant source of water resources for many regions, especially in the temperate and frigid zones. We introduce a snow module based on surface energy balance and Budyko-limits on melting and combine it with the existing water balance equations. Moreover, monthly parameterization is applied to the model to better explain the time-varying hydrological characteristics of a catchment. The proposed model involves five different monthly parameters, which determine the physical partitioning process of the hydrological cycle, and they are regionally calibrated and validated under Budyko-type constraints. The model is applied to 1,210 basins across the continental United States (CONUS), and the simulated streamflow is compared to the observed data. The proposed model significantly outperformed the original model, improving the median NSE by 31% (from 0.51 to 0.67) and increasing the number of catchments with an acceptable NSE by 58%. The spatial variability of the basin characteristics across the CONUS is also investigated based on the calibrated parameters.
The significance of phloem hydrodynamics to plant mortality and survival, which impacts ecosystem-scale carbon and water cycling, is not in dispute. The phloem provides the conduits for products of photosynthesis to be transported to different parts of the plant for consumption or storage. The M$\ddot{u}$nch pressure flow hypothesis (PFH) is the leading framework to mathematically represent this transport. It assumes that osmosis provides the necessary pressure differences to drive water and sucrose within the phloem. Mathematical models utilizing the PFH approximate the phloem by a relatively rigid slender semi-permeable tube. However, the phloem consists of living cells that contract and expand in response to pressure fluctuations. The effect of membrane elasticity on osmotically driven sucrose front speed has rarely been considered and frames the scope here. Laboratory experiments were conducted to elucidate the elastic-to-plastic pressure-deformation relation in membranes and their effect on sucrose front speeds. It is demonstrated that membrane elasticity acts to retard the sucrose front speed. The retardation emerges because of two effects: (a) part of the osmotic pressure is diverted to perform mechanical work to expand the membrane instead of pressurizing water, and (b) expansion of the membrane reduces the sucrose concentration driving osmotic potential due to volume increases and concomitant dilution effects. Finally, these results offer a novel perspective about the much discussed presence of sieve plates throughout the phloem acting as structural expansion dampers.