SEARCH · Engineering Papers

Results for “mean”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Accuracy of HVAC Load Predictions: Validation of EnergyPlus and DOE-2 using FLEXLAB Measurements

The aim of the project reported here was to better understand the level of accuracy of three building energy simulation (BES) engines (‘engines’) — EnergyPlus™, DOE-2.1e, and DOE-2.2 — by identifying and investigating significant deviations between the performance predicted by these engines and actual performance as measured in the FLEXLAB® test facility at Lawrence Berkeley National Laboratory (LBNL). The specific test conditions included some of those prescribed in ANSI/ASHRAE Standard 140 - Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs. Detailed measurements of FLEXLAB performance, including indoor temperatures and heat fluxes and air-flow and water flow rates and temperatures in the Heating, Ventilating and Air Conditioning (HVAC) system, together with hourly weather data, were recorded and used in analyzing the simulation results from EnergyPlus v8.8, DOE-2.2 v3.65 and DOE-2.1e v127. These engines are commonly used in the United States for building energy code compliance, federal, state, and utility incentives programs, as well as energy efficient design of new buildings and energy retrofit of existing buildings. Seven conventional overhead mixing ventilation scenarios were tested and each engine was found to have a similar level of agreement with the measurements of space-level heating and sensible cooling loads. These results provide useful information regarding the accuracy of these engines in predicting the cooling and heating load elements of whole building energy performance. This information is intended for practitioners who are concerned about transitioning between simulation tools with different engines and for managers of utility programs leveraging these tools for evaluating and/or projecting measure savings to be incentivized under their programs. The results of the comparisons of simulated and measured performance indicate that the predictions from all three engines are not significantly different. The 24-hour average value of the absolute mean bias indicates the likely magnitude of the error in any particular case. The average mean bias is reduced by cancelation of overprediction in one case by underprediction in another. The daytime absolute mean biases, which may be more important for both energy performance and occupant comfort, are ~6%, presumably because of the greater complexity involved in simulating in the presence of solar radiation. EnergyPlus typically overpredicts the cooling load and/or underpredicts the heating load by ~1.5% and the DOE-2 engines typically underpredict the cooling load by approximately the same amount. The Root Mean Square Error is relatively more sensitive to shorter term variations in the difference between predicted and measured loads; the three engines have similar values, ~10%, suggesting that the uncertainties in their predictions of peak loads may also be similar in magnitude. The implication of these results is that users, both designers and program analysts, can use EnergyPlus, DOE-2.1e, or DOE-2.2 to model conventional commercial buildings equipped with overhead mixing ventilation with a similar level of confidence. Further work is required to better understand the variability in the level of agreement between the engine predictions and FLEXLAB measurements, where a particular engine will agree well with FLEXLAB in some cases and not so well in others and another engine will agree or disagree in different cases. As the sources of this variability are identified and eliminated or reduced significantly, it is recommended that the experimental capabilities and methods developed in the study reported here should be applied to validating heating and cooling load calculations for spaces with different types of furniture and miscellaneous loads. These methods should then be applied to low energy space conditioning systems in EnergyPlus including, in particular, radiant slab and radiant ceiling panel cooling and heating systems and ‘mixed mode’ systems that combine mechanical cooling and natural ventilation systems, focusing on controls, including control of thermal mass. The work reported here addresses the conventional method of heating and cooling occupied spaces; other methods, such as the use of radiant heating and cooling systems have the potential to provide equivalent occupant comfort, or better, with lower energy consumption. These systems are addressed more explicitly in EnergyPlus but there is a need for empirical validation to give users the same level of confidence in modeling these systems that they have, or should have, in modeling conventional systems, based on the results presented here.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

A computational study of turbulent separated flow over a wall-mounted cube at two different Reynolds numbers and incoming velocity profiles

Simulations of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at Reynolds numbers (Re) where the bulk flow characteristics are known to become relatively Re-insensitive using both direct numerical simulation (DNS) and wall-resolved large eddy simulation (WRLES). The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions are selected to allow for easy comparisons for model development purposes and avoid confounding effects from freestream turbulence. In addition to mean velocity profiles, select turbulence statistics are presented in detail. We find that both mean velocity profiles and mean wake turbulence are Re-insensitive over the range examined here.

42 ENGINEERING↗

DOE OSTI2021

Online Monitoring System for Concrete Structures Affected by Alkali-Silica Reaction

This report presents a comprehensive study of developing ultrasonic wave and acoustic emission (AE) techniques for long-term monitoring of alkali-silica reaction (ASR) development in concrete. Small, medium, and full-scale concrete specimens were cast, conditioned, and monitored for various periods, from at least one year to 2.4 years. The reactive concrete specimens contain coarse or fine reactive aggregates to study the effects of different types of reactive aggregates. Confinements were also designed to simulate the 2-dimensional confinement effects of reinforcement in the shield building of nuclear power plants. The concrete specimens were stored in a environmental chamber with high humidity and high temperature to accelerate the ASR development. The ultrasonic monitoring data shows high sensitivity to ASR development and could detect cracking initiation well before visible surface cracks occurred. However, the linear ultrasonic analysis based on wave velocity is strongly affected by temperature variation. Therefore, a nonlinear ultrasonic method was proposed to measure thermally induced nonlinear acoustic responses of concrete (thermal modulation of ultrasonic wave). The measured nonlinear acoustic parameter shows a high correlation with ASR expansion across specimens with different reactive aggregates and confinement conditions. The same conclusion was obtained from nonlinear resonance tests on small concrete prisms. Compared to the linear acoustic methods, the nonlinear tests show high sensitivities to ASR damage from internal microcracking initiation at the early stage to visible cracks at the late stage of ASR. The attributes of the thermal modulation of nonlinear ultrasonic method include high sensitivity, immunity to temperature effects, and strong correlation with ASR expansion, which demonstrate great potentials of the nonlinear ultrasonic method for diagnosis of ASR damage and prediction of concrete deterioration process. Acoustic emission is a passive sensing technique for damage assessment, and access to only one surface is needed even for thick and heavily reinforced elements such as the walls utilized for nuclear shield building. Additionally, relatively few sensors are required to monitor the progression of the damage process. Results indicate that the AE data can be related to the damage rating index, which is a petrography-based means of assessing damage due to ASR in reinforced concrete. Furthermore, AE is capable of detecting ASR damage long before surface cracking is visually noticeable. Boundary conditions play in important role in the progression of ASR damage and differences in boundary conditions are reflected in the AE data. Entropy based data assessment methods provide a means to assess the damage state, and convolutional neural network based data assessment provided a means to assess the damage state in real-time. Results indicate that artificial neural network models may be used as a means to predict volumetric expansion based on AE data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Release of Evaluated 235 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 235 U from 200 keV to 15 MeV that is a potential release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation had to be re-done from "scratch", as the input to the $\overline{v}_p$ evaluation of the previous library, ENDF/B-VIII.0, was lost. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. Another major difference to ENDF/B-VIII.0 is that this evaluation includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to perform two types of validation for the new 235 U $\overline{v}_p$: On the one hand, one can employ evaluated CGMF parameters obtained from fitting to experimental 235 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data indicating the evaluated 235 U(n,f) $\overline{v}_p$ are physical. On the other hand, one can validate 235 U $\overline{v}_p$ with respect to integral responses such as fast ICSBEP critical assemblies or LLNL pulsed spheres. LLNL pulsed-sphere neutron-leakage spectra are minimally impacted by the new 235 U $\overline{v}_p$ as these experimental data are shape data and the $\overline{v}_p$ would mostly lead to a change in normalization of the data as the spheres are relatively thin (0.7 and 1.5 mean-free path) and, thus, mostly depend on 235 U $\overline{v}_p$ from 12-15 MeV. The change in the predicted effective neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, however, is large compared to values using ENDF/B-VIII.0 and experimental k eff : The average bias is 108 pcm across all studied k eff values versus 12 pcm for ENDF/B-VIII.0. A reasonable performance in simulating keff (mean bias of 14 pcm) can be retained by tweaking 235 U $\overline{v}_p$ from 3-5 MeV, and combining it with a recent 235 U PFNS evaluation that is also a ENDF/B-VIII.1 release candidate.

235U↗

DOE OSTI2023

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a high-order spectral element code Nek5000/RS, which is supported under the DOE's Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms. In our earlier ANL report, we demonstrated our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗

DOE OSTI2025

Front-of-Meter Model Results

These files contains aggregations of key variables from the NREL Distributed Wind Futures Study using full parcel level data. These variables describe total technical and economic potential for distributed wind turbine deployment. Aggregations are available at the (1) county, (2) zipcode (zip code tabulation area or zcta), and (3) US Census block group level. Each scenario is coded with the scenario name (e.g., baseline) and year (e.g., 2022). Those files postfixed with 'econpot' contain results for only those parcels that are economically viable while the files postfixed with 'techpot' include results for all parcels that are technically feasible. Hence these correspond to technoeconomic and technical potential respectively. The data are available as CSV or Geopackage. Columns in the files are as follows: * geoid: geographic identifier (FIPS code or similar) * min_techpot_sum_kw: technical potential for all parcels in kW using turbines downsized to demand when appropriate * max_techpot_sum_kw: technical potential for all parcels in kW without downsizing turbines * aep_sum_kwh: annual energy production estimate in kWh * cf_mean_ratio: mean capacity factor * lcoe_mean_cents_per_kwh: mean levelized cost of energy for parcels in geography in cents per kWh * lcoe_std_cents_per_kwh: standard deviation of the above * parcel_area_sum_acres: total area of viable parcels in acres * n_turbines: number of cited turbines (one per viable parcel currently) Note: These are preliminary results from the full-parcel 2024 update of the Distributed Wind Energy Futures study. Please take care when making use of the data, and feel free to contact the team with any questions. Full documentation in support of these data is in progress and will follow.

17 WIND ENERGY↗

DOE OSTI2025

Behind-the-Meter Model Results

17 WIND ENERGY↗

DOE OSTI2023

Tropical Tree Crop Simulation with a Process-Based, Daily Timestep Simulation Model (ALMANAC): Description of Model Adaptation and Examples with Coffee and Cocoa Simulations

Coffee (Coffea species) and Cocoa (Theobroma cacao) are important cash crops grown in the tropics but traded globally. This study was conducted to apply the ALMANAC model to these crops for the first time, and to test its ability to simulate them under agroforestry management schemes and varying precipitation amounts. To create this simulation, coffee was grown on a site in Kaua’i, Hawai’i, USA, and cocoa was grown on a site in Sefwi Bekwai, Ghana. A stand-in for a tropical overstory tree was created for agroforestry simulations using altered parameters for carob, a common taller tropical tree for these regions. For both crops, ALMANAC was able to realistically simulate yields when compared to the collected total yield data. On Kaua’i, the mean simulated yield was 2% different from the mean measured yield, and in all three years, the simulated values were within 10% of the measured values. For cocoa, the mean simulated yield was 3% different from the mean measured yield and the simulated yield was within 10% of measured yields for all four available years. When precipitation patterns were altered, in Ghana, the wetter site showed lower percent changes in yield than the drier site in Hawai’i. When agroforestry-style management was simulated, a low Leaf Area Index (LAI) of the overstory showed positive or no effect on yields, but when LAI climbed too high, the simulation was able to show the detrimental effect this competition had on crop yields. These simulation results are supported by other literature documenting the effects of agroforestry on tropical crops. This research has applied ALMANAC to new crops and demonstrated its simulation of different management and environmental conditions. The results show promise for ALMANAC’s applicability to these scenarios as well as its potential to be further tested and utilized in new circumstances.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Coreset Clustering on Small Quantum Computers

Many quantum algorithms for machine learning require access to classical data in superposition. However, for many natural data sets and algorithms, the overhead required to load the data set in superposition can erase any potential quantum speedup over classical algorithms. Recent work by Harrow introduces a new paradigm in hybrid quantum-classical computing to address this issue, relying on coresets to minimize the data loading overhead of quantum algorithms. We investigated using this paradigm to perform k-means clustering on near-term quantum computers, by casting it as a QAOA optimization instance over a small coreset. We used numerical simulations to compare the performance of this approach to classical k-means clustering. We were able to find data sets with which coresets work well relative to random sampling and where QAOA could potentially outperform standard k-means on a coreset. However, finding data sets where both coresets and QAOA work well—which is necessary for a quantum advantage over k-means on the entire data set—appears to be challenging.

42 ENGINEERING↗

DOE OSTI2022

Antiviral Activity and Crystal Structures of HIV-1 gp120 Antagonists

As part of our effort to discover drugs that target HIV-1 entry, we report the antiviral activity and crystal structures of two novel inhibitors in a complex with a gp120 core. NBD-14204 showed similar antiviral activity against all the clinical isolates tested. The IC50 values were in the range of 0.24–0.9 µM with an overall mean of 0.47 ± 0.03 µM, showing slightly better activity against the clinical isolates than against the lab-adapted HIV-1 HXB2 (IC 50 = 0.96 ± 0.1 µM). Moreover, the antiviral activity of NBD-14208 was less consistent, showing a wider range of IC 50 values (0.66–5.7 µM) with an overall mean of 3 ± 0.25 µM and better activity against subtypes B and D (Mean IC 50 2.2–2.5 µM) than the A, C and Rec viruses (Mean IC 50 2.9–3.9 µM). SI of NBD-14204 was about 10-fold higher than NBD-14208, making it a better lead compound for further optimization. In addition, we tested these compounds against S375Y and S375H mutants of gp120, which occurred in some clades and observed these to be sensitive to NBD-14204 and NBD-14208. These inhibitors also showed modest activity against HIV-1 reverse transcriptase. Furthermore, we determined the crystal structures of both inhibitors in complexes with gp120 cores. As expected, both NBD-14204 and NBD-14208 bind primarily within the Phe43 cavity. It is noteworthy that the electron density of the thiazole ring in both structures was poorly defined due to the flexibility of this scaffold, suggesting that these compounds maintain substantial entropy, even when bound to the Phe43 cavity.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

The Evolution and Mass Dependence of Galaxy Cluster Pressure Profiles at 0.05 ≤ z ≤ 0.60 and 4 × 10 14 M ⊙ ≤ M 500 ≤ 30 × 10 14 M ⊙

Abstract We have combined X-ray observations from Chandra with Sunyaev–Zel’dovich effect data from Planck and Bolocam to measure intracluster medium pressure profiles from 0.03 R 500 ≤ R ≤ 5 R 500 for a sample of 21 low- z galaxy clusters with a median redshift of 〈 z 〉 = 0.08 and a median mass of 〈 M 500 〉 = 6.1 × 10 14 M ⊙ and a sample of 19 mid- z galaxy clusters with 〈 z 〉 = 0.50 and 〈 M 500 〉 = 10.6 × 10 14 M ⊙ . The mean scaled pressure in the low- z sample is lower at small radii and higher at large radii, a trend that is accurately reproduced in similarly selected samples from The Three Hundred simulations. This difference appears to be primarily due to dynamical state at small radii, evolution at intermediate radii, and a combination of evolution and mass dependence at large radii. Furthermore, the overall flattening of the mean scaled pressure profile in the low- z sample compared to the mid- z sample is consistent with expectations due to differences in the mass accretion rate and the fractional impact of feedback mechanisms. In agreement with previous studies, the fractional scatter about the mean scaled pressure profile reaches a minimum of ≃20% near 0.5 R 500 . This scatter is consistent between the low- z and mid- z samples at all radii, suggesting it is not strongly impacted by sample selection, and this general behavior is reproduced in The Three Hundred simulations. Finally, analytic functions that approximately describe the mass and redshift trends in mean pressure profile shape are provided.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The Chemodynamics of the Stellar Populations in M31 from APOGEE Integrated-light Spectroscopy

We present an analysis of nearly 1000 near-infrared, integrated-light spectra from APOGEE in the inner ~7 kpc of M31. We utilize full-spectrum fitting with A-LIST simple stellar population spectral templates that represent a population of stars with the same age, [M/H], and [α/M]. With this, we determine the mean kinematics, metallicities, α abundances, and ages of the stellar populations of M31's bar, bulge, and inner disk (~4–7 kpc). We find a nonaxisymmetric velocity field in M31 resulting from the presence of a bar. The bulge of M31 is less metal-rich (mean [M/H] = $0.149$$^{+0.067}_{-0.081}$ dex) than the disk, features minima in metallicity on either side of the bar ([M/H] ~ -0.2), and is enhanced in α abundance (mean [α/M] = $0.281$$^{+0.035}_{-0.038}$). The disk of M31 within ~7 kpc is enhanced in both metallicity ([M/H] = $-0.023$$^{+0.050}_{-0.052}$) and α abundance ([α/M] = $0.274$$^{+0.020}_{-0.025}$). Both of these structural components are uniformly old at ≃12 Gyr. We find the mean metallicity increases with distance from the center of M31, with the steepest gradient along the disk major axis (0.043 ± 0.021 dex kpc -1 ). This gradient is the result of changing light contributions from the bulge and disk. The chemodynamics of stellar populations encodes information about a galaxy's chemical enrichment, star formation history, and merger history, allowing us to discuss new constraints on M31's formation. Our results provide a stepping stone between our understanding of the Milky Way and other external galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Data from: Comparing traditional and Bayesian approaches to ecological meta-analysis

Despite the wide application of meta-analysis in ecology, some of the traditional methods used for meta-analysis may not perform well given the type of data characteristic of ecological meta-analyses. We reviewed published meta-analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies (ni) and the number of studies or records (k) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta-analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that ni and k were highly variable, distributions were right-skewed, and were generally small (median ni =5, median k=44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval based on the t-distribution that uses an adjusted standard error (the Hartung-Knapp-Sidik-Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z-distribution confidence intervals had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta-analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random-effects meta-analysis, the performance of the best frequentist and Bayesian methods were similar for the same combinations of factors (k and mean replication), though the Bayesian approaches had higher than nominal (>95%) coverage for the mean effect when k was very low (k<15). Our literature review suggests that many meta-analyses that used z-distribution or bootstrapping confidence intervals may have over-estimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta-analyses.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Identifying the sources of uncertainty in climate model simulations of solar radiation modification with the G6sulfur and G6solar Geoengineering Model Intercomparison Project (GeoMIP) simulations

We present here results from the Geoengineering Model Intercomparison Project (GeoMIP) simulations for the experiments G6sulfur and G6solar for six Earth system models participating in the Climate Model Intercomparison Project (CMIP) Phase 6. The aim of the experiments is to reduce the warming that results from a high-tier emission scenario (Shared Socioeconomic Pathways SSP5-8.5) to that resulting from a medium-tier emission scenario (SSP2-4.5). These simulations aim to analyze the response of climate models to a reduction in incoming surface radiation as a means to reduce global surface temperatures, and they do so either by simulating a stratospheric sulfate aerosol layer or, in a more idealized way, through a uniform reduction in the solar constant in the model. We find that over the final two decades of this century there are considerable inter-model spreads in the needed injection amounts of sulfate (29 ± 9 Tg-SO 2 /yr between 2081 and 2100), in the latitudinal distribution of the aerosol cloud and in the stratospheric temperature changes resulting from the added aerosol layer. Even in the simpler G6solar experiment, there is a spread in the needed solar dimming to achieve the same global temperature target (1.91 ± 0.44 %). The analyzed models already show significant differences in the response to the increasing CO 2 concentrations for global mean temperatures and global mean precipitation (2.05 K ± 0.42 K and 2.28 ± 0.80 %, respectively, for SSP5-8.5 minus SSP2-4.5 averaged over 2081–2100). With aerosol injection, the differences in how the aerosols spread further change some of the underlying uncertainties, such as the global mean precipitation response (–3.79 ± 0.76 % for G6sulfur compared to –2.07 ± 0.40 % for G6solar against SSP2-4.5 between 2081 and 2100). These differences in the behavior of the aerosols also result in a larger uncertainty in the regional surface temperature response among models in the case of the G6sulfur simulations, suggesting the need to devise various, more specific experiments to single out and resolve particular sources of uncertainty. The spread in the modeled response suggests that a degree of caution is necessary when using these results for assessing specific impacts of geoengineering in various aspects of the Earth system. However, all models agree that compared to a scenario with unmitigated warming, stratospheric aerosol geoengineering has the potential to both globally and locally reduce the increase in surface temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Multi-century dynamics of the climate and carbon cycle under both high and net negative emissions scenarios

Abstract. Future climate projections from Earth system models (ESMs) typically focus on the timescale of this century. We use a set of five ESMs and one Earth system model of intermediate complexity (EMIC) to explore the dynamics of the Earth's climate and carbon cycles under contrasting emissions trajectories beyond this century to the year 2300. The trajectories include a very-high-emissions, unmitigated fossil-fuel-driven scenario, as well as a mitigation scenario that diverges from the first scenario after 2040 and features an “overshoot”, followed by a decrease in atmospheric CO2 concentrations by means of large net negative CO2 emissions. In both scenarios and for all models considered here, the terrestrial system switches from being a net sink to either a neutral state or a net source of carbon, though for different reasons and centered in different geographic regions, depending on both the model and the scenario. The ocean carbon system remains a sink, albeit weakened by carbon cycle feedbacks, in all models under the high-emissions scenario and switches from sink to source in the overshoot scenario. The global mean temperature anomaly is generally proportional to cumulative carbon emissions, with a deviation from proportionality in the overshoot scenario that is governed by the zero emissions commitment. Additionally, 23rd century warming continues after the cessation of carbon emissions in several models in the high-emissions scenario and in one model in the overshoot scenario. While ocean carbon cycle responses qualitatively agree in both globally integrated and zonal mean dynamics in both scenarios, the land models qualitatively disagree in zonal mean dynamics, in the relative roles of vegetation and soil in driving C fluxes, in the response of the sink to CO2, and in the timing of the sink–source transition, particularly in the high-emissions scenario. The lack of agreement among land models on the mechanisms and geographic patterns of carbon cycle feedbacks, alongside the potential for lagged physical climate dynamics to cause warming long after CO2 concentrations have stabilized, points to the possibility of surprises in the climate system beyond the 21st century time horizon, even under relatively mitigated global warming scenarios, which should be taken into consideration when setting global climate policy.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production

Terrestrial photosynthesis is the basis for vegetation growth and drives the land carbon cycle. Accurately simulating gross primary production (GPP, ecosystem-level apparent photosynthesis) is key for satellite monitoring and Earth system model predictions under climate change. While robust models exist for describing leaf-level photosynthesis, predictions diverge due to uncertain photosynthetic traits and parameters which vary on multiple spatial and temporal scales. Here, we describe and evaluate a GPP (photosynthesis per unit ground area) model, the P-model, that combines the Farquhar–von Caemmerer–Berry model for C 3 photosynthesis with an optimality principle for the carbon assimilation–transpiration trade-off, and predicts a multi-day average light use efficiency (LUE) for any climate and C 3 vegetation type. The model builds on the theory developed in Prentice et al. (2014) and Wang et al. (2017a) and is extended to include low temperature effects on the intrinsic quantum yield and an empirical soil moisture stress factor. The model is forced with site-level data of the fraction of absorbed photosynthetically active radiation (fAPAR) and meteorological data and is evaluated against GPP estimates from a globally distributed network of ecosystem flux measurements. Although the P-model requires relatively few inputs, the R 2 for predicted versus observed GPP based on the full model setup is 0.75 (8 d mean, 126 sites) – similar to comparable satellite-data-driven GPP models but without predefined vegetation-type-specific parameters. The R 2 is reduced to 0.70 when not accounting for the reduction in quantum yield at low temperatures and effects of low soil moisture on LUE. The R 2 for the P-model-predicted LUE is 0.32 (means by site) and 0.48 (means by vegetation type). Applying this model for global-scale simulations yields a total global GPP of 106–122 Pg C yr –1 (mean of 2001–2011), depending on the fAPAR forcing data. The P-model provides a simple but powerful method for predicting – rather than prescribing – light use efficiency and simulating terrestrial photosynthesis across a wide range of conditions. The model is available as an R package (rpmodel).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The GGCMI Phase 2 emulators: global gridded crop model responses to changes in CO 2 , temperature, water, and nitrogen (version 1.0)

Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: that growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project

Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.

Geology↗