Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model intercomparison”

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 163 records · Page 9

The Terrestrial Biosphere Model Farm

Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901–2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.

54 ENVIRONMENTAL SCIENCES↗

CONUS-wide Projected Flood Frequency and Uncertainty Estimates, Version 1.0

This dataset presents a large-ensemble of CONUS-wide projected flood frequency and uncertainty estimates across ~2.7 million NHDPlusV2 river reaches over the CONUS. The framework producing this dataset leverages a multi-model, uncertainty-aware modeling framework that allows evaluating shifts in flood frequences at the stream reach level across the CONUS. CONUS-wide ensemble streamflow projections generated from hydrologic simulations driven by downscaled and bias-corrected Coupled Model Intercomparison Project Phase 6 (CMIP6) outputs are used to derive these flood frequency and uncertainty estimates over the period 1980 - 2099. A spatially consistent regional L-moment algorithm is applied across clusters defined by the US Hydrologic Unit Code Subregions (HUC4s and HUC8s) and NHDPlusV2 stream orders to estimate flood frequencies. The dataset also includes at-site based flood estimates that allow for the comparison between local and regional approach-based estimates, assess projected changes, and characterize their uncertainties. For more reliable estimation of rare flood frequencies such as 500 and 1000-year return periods, super-ensemble based estimates are also included in the dataset. This dataset is derived to support the "Impact-Informed Dam Safety Risk Assessment for Securing Hydropower Assests" project for the US Department of Energy (DOE) Hydropower and Hydrokinetic Office (H2O). For further details, refer to Kao et al. (2022), Ghimire et al. (2023), Ghimire et al. (2025), and Hosking and Wallis (1997).

Ghimire, Ganesh [ORNL] (ORCID:0000000242843941)↗

The Surface Albedo of Sea Ice in CMIP6 and the Implications for the Surface Albedo Feedback

The Arctic has experienced rapid sea ice loss and a substantial decline in surface albedo, significantly impacting its radiation budget. Climate models from the Coupled Model Intercomparison Project (CMIP6) reproduce these changes. However, inconsistencies remain among models regarding the magnitude, spatial distribution, and seasonal patterns of Arctic surface albedo evolution. This study investigates these discrepancies by comparing model outputs with observation from the Clouds and the Earth's Radiant Energy System (CERES). We develop a decomposition method to assess the contributions of sea ice albedo, sea ice concentration, and sea ice extent to Arctic surface albedo. Over land, differences in snow cover account for the substantial inter-model spread in surface albedo, while over the ocean, sea ice albedo, concentration, and extent all contribute. Comparisons between CMIP6 and the Atmospheric Model Intercomparison Project (AMIP) simulations, which use identically prescribed sea ice concentrations, reveals considerable inter-model spread in Arctic Ocean surface albedo due to differences in sea ice albedo. Applying the decomposition method to projections shows that models predicting larger decreases in sea ice concentration and extent, especially in the Central Arctic, exhibit lower surface albedo and stronger sea ice albedo feedback. Beyond 2050, Arctic Ocean surface albedo decline is mainly influenced by sea ice extent indicating that the retreat of the ice edge is the most important process to constrain the surface albedo feedback. This study provides insights into factors contributing to the spread and changes in Arctic surface albedo and the associated sea ice albedo feedback.

Patrick C Taylor↗

Effective radiative forcing and adjustments in CMIP6 models

Abstract. The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

54 ENVIRONMENTAL SCIENCES↗

Effective radiative forcing and adjustments in CMIP6 models

The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W/sq. m, comprised of 1.81 (±0.09) W/sq. m from CO2, 1.08 (± 0.21) W/sq. m from other well-mixed greenhouse gases, −1.01 (± 0.23) W/sq. m from aerosols and −0.09 (±0.13) W/sq. m from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W/sq. m. The majority of the remaining 0.21 W/sq. m is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W/sq. m, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

effective radiative forcing↗

The GGCMI Phase 2 Experiment: Global Gridded Crop Model Simulations Under Uniform Changes in CO2, Temperature, Water, and Nitrogen Levels (Protocol Version 1.0)

Concerns about food security under climate change motivate efforts to better understand future changes in crop yields. Process-based crop models, which represent plant physiological and soil processes, are necessary tools for this purpose since they allow representing future climate and management conditions not sampled in the historical record and new locations to which cultivation may shift. However, process-based crop models differ in many critical details, and their responses to different interacting factors remain only poorly understood. The Global Gridded Crop Model Intercomparison (GGCMI) Phase 2 experiment, an activity of the Agricultural Model Intercomparison and Improvement Project (AgMIP), is designed to provide a systematic parameter sweep focused on climate change factors and their interaction with overall soil fertility, to allow both evaluating model behavior and emulating model responses in impact assessment tools. In this paper we describe the GGCMI Phase 2 experimental protocol and its simulation data archive. A total of 12 crop models simulate five crops with systematic uniform perturbations of historical climate, varying CO2, temperature, water supply, and applied nitrogen (“CTWN”) for rainfed and irrigated agriculture, and a second set of simulations represents a type of adaptation by allowing the adjustment of growing season length. We present some crop yield results to illustrate general characteristics of the simulations and potential uses of the GGCMI Phase 2 archive. For example, in cases without adaptation, modeled yields show robust decreases to warmer temperatures in almost all regions, with a nonlinear dependence that means yields in warmer baseline locations have greater temperature sensitivity. Inter-model uncertainty is qualitatively similar across all the four input dimensions but is largest in high-latitude regions where crops may be grown in the future.

GGCMI↗

A multi-model CMIP6-PMIP4 study of Arctic sea ice at 127 ka: sea ice data compilation and model differences

The Last Interglacial period (LIG) is a period with increased summer insolation at high northern latitudes, which results in strong changes in the terrestrial and marine cryosphere. Understanding the mechanisms for this response via climate modelling and comparing the models' representation of climate reconstructions is one of the objectives set up by the Paleoclimate Modelling Intercomparison Project for its contribution to the sixth phase of the Coupled Model Intercomparison Project. Here we analyse the results from 16 climate models in terms of Arctic sea ice. The multi-model mean reduction in minimum sea ice area from the pre industrial period (PI) to the LIG reaches 50 % (multi-model mean LIG area is 3.20×10^6 sq.km, compared to 6.46×10^6 sq.km for the PI). On the other hand, there is little change for the maximum sea ice area (which is 15–16×10^6 sq.km for both the PI and the LIG. To evaluate the model results we synthesise LIG sea ice data from marine cores collected in the Arctic Ocean, Nordic Seas and northern North Atlantic. The reconstructions for the northern North Atlantic show year-round ice-free conditions, and most models yield results in agreement with these reconstructions. Model–data disagreement appear for the sites in the Nordic Seas close to Greenland and at the edge of the Arctic Ocean. The northernmost site with good chronology, for which a sea ice concentration larger than 75 % is reconstructed even in summer, discriminates those models which simulate too little sea ice. However, the remaining models appear to simulate too much sea ice over the two sites south of the northernmost one, for which the reconstructed sea ice cover is seasonal. Hence models either underestimate or overestimate sea ice cover for the LIG, and their bias does not appear to be related to their bias for the pre-industrial period. Drivers for the inter-model differences are different phasing of the up and down short-wave anomalies over the Arctic Ocean, which are associated with differences in model albedo; possible cloud property differences, in terms of optical depth; and LIG ocean circulation changes which occur for some, but not all, LIG simulations. Finally, we note that inter-comparisons between the LIG simulations and simulations for future climate with moderate (1 %/yr) CO2 increase show a relationship between LIG sea ice and sea ice simulated under CO2 increase around the years of doubling CO2. The LIG may therefore yield insight into likely 21st century Arctic sea ice changes using these LIG simulations.

Arctic sea ice↗

Global and northern-high-latitude net ecosystem production in the 21st century from CMIP6 experiments

Climate warming is accelerating the changes in the global terrestrial ecosystems and particularly those in the northern high latitudes (NHLs; poleward of 50° N) and rendering the land–atmosphere carbon exchange highly uncertain. The Coupled Model Intercomparison Project Phase 6 (CMIP6) employs the most updated climate models to estimate terrestrial ecosystem carbon dynamics driven by a new set of socioeconomic and climate change pathways. By analyzing the future (2015–2100) carbon fluxes estimated by 10 CMIP6 models, we quantitatively evaluated the projected magnitudes, trends, and uncertainties in the global and NHL carbon fluxes under four scenarios plus the role of NHLs in the global terrestrial ecosystem carbon dynamics. Overall, the models suggest that the global and NHL terrestrial ecosystems will be consistent carbon sinks in the future, and the magnitude of the carbon sinks is projected to be larger under scenarios with higher radiative forcing. By the end of this century, the models on average estimate the NHL net ecosystem productivity (NEP) as 0.54 ± 0.77, 1.01 ± 0.98, 0.97 ± 1.62, and 1.05 ± 1.83 Pg C yr –1 under SSP126, SSP245, SSP370, and SSP585, respectively. The uncertainties are not substantially reduced compared with earlier results, e.g., the Coupled Climate–Carbon Cycle Model Intercomparison Project (C4MIP). Although NHLs contribute a small fraction of the global carbon sink (~13 %), the relative uncertainties in NHL NEP are much larger than the global level. Our results provide insights into future carbon flux evolutions under future scenarios and highlight the urgent need to constrain the large uncertainties associated with model projections for making better climate mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

54 ENVIRONMENTAL SCIENCES↗

ISMIP6 Antarctica: A Multi-Model Ensemble of the Antarctic Ice Sheet Evolution Over the 21st Century

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between -7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica ass change varies between -6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

Hélène Seroussi↗

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗

helios: An R package to process heating and cooling degrees for GCAM

helios is an open-source R package that estimates population-weighted heating and cooling degree-hours (HDH and CDH) and degree-days (HDD and CDD) at various temporal (e.g., energy dispatch segments, monthly, yearly) and spatial scales (e.g., U.S. states, global political regions, countries). The degree hour and degree day outputs from helios are used to inform electricity demand load in the Global Change Analysis Model (GCAM) as well as in GCAM-USA (which is the version of GCAM with U.S. state-level details). helios uses a workflow with four steps: processing raw data; calculating heating and cooling degrees; visualizing performance diagnostics; and outputing results in various formats. There are two sources of widely-used climate data compatible with helios: (1) hourly climate data with 12-km resolution that are dynamically downscaled with the Weather Research and Forecasting (WRF) model and projected using a thermal global warming (TGW) approach; and (2) daily climate data with 0.5-degree resolution from the Coupled Model Intercomparison Project (CMIP) that is bias-adjusted and statistical downscaled by the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP). In summary, helios is a model that standardizes methodology of heating and cooling degrees-hours and degree-days using publicly available data and advance the understanding of the impact of spatial and temporal temperature variability on building energy services.

97 MATHEMATICS AND COMPUTING↗

Carbon-Temperature-Water Change Analysis for Peanut Production Under Climate Change: A Prototype for the AgMIP Coordinated Climate-Crop Modeling Project (C3MP)

Climate change is projected to push the limits of cropping systems and has the potential to disrupt the agricultural sector from local to global scales. This article introduces the Coordinated Climate-Crop Modeling Project (C3MP), an initiative of the Agricultural Model Intercomparison and Improvement Project (AgMIP) to engage a global network of crop modelers to explore the impacts of climate change via an investigation of crop responses to changes in carbon dioxide concentration ([CO2]), temperature, and water. As a demonstration of the C3MP protocols and enabled analyses, we apply the Decision Support System for Agrotechnology Transfer (DSSAT) CROPGRO-Peanut crop model for Henry County, Alabama, to evaluate responses to the range of plausible [CO2], temperature changes, and precipitation changes projected by climate models out to the end of the 21st century. These sensitivity tests are used to derive crop model emulators that estimate changes in mean yield and the coefficient of variation for seasonal yields across a broad range of climate conditions, reproducing mean yields from sensitivity test simulations with deviations of ca. 2% for rain-fed conditions. We apply these statistical emulators to investigate how peanuts respond to projections from various global climate models, time periods, and emissions scenarios, finding a robust projection of modest (<10%) median yield losses in the middle of the 21st century accelerating to more severe (>20%) losses and larger uncertainty at the end of the century under the more severe representative concentration pathway (RCP8.5). This projection is not substantially altered by the selection of the AgMERRA global gridded climate dataset rather than the local historical observations, differences between the Third and Fifth Coupled Model Intercomparison Project (CMIP3 and CMIP5), or the use of the delta method of climate impacts analysis rather than the C3MP impacts response surface and emulator approach.

climate change↗

Regional-scale, sector-specific evaluation of global CO2 inversion models using aircraft data from the ACT-America project

We use vertical profiles of airborne measurements of CO2 from frontal cases during the summer 2016 Atmospheric Carbon and Transport – America (ACT-America) campaign to evaluate the skill of a set of ten global CO2 inversion models participating in the Orbiting Carbon Observatory – 2 (OCO-2) Model Intercomparison Project (MIP). Model errors and biases (model minus observation) were categorized by region (Mid-Atlantic, Midwest, and South), frontal sector (warm or cold), and transport model (predominantly Tracer Model 5 (TM5) and Goddard Earth Observing System – Chemistry (GEOS-Chem)). Overall, the inversions reproduce the general structures of the observed vertical profiles and the enhanced / depleted low-level CO2 in warm / cold sectors, but tend to underestimate the magnitude of the sector difference in each region. In the Midwest and South warm sectors, inversion biases were about 1 ppm above 1500 m AGL, though model spread (quantified by interquartile range) was often even smaller; below 1500 m AGL, model biases were negative and about -2 ppm near the surface. For the Midwest and South cold sectors, models had +2-3 ppm biases below 1500 m AGL, but with comparable model spread. Uniquely, in the Mid-Atlantic there was a consistent difference between TM5 and GEOS-Chem mole fractions for both sectors up to 3000 m AGL (TM5 lower by 2 ppm), on the order of the observation-relative biases. In the MidAtlantic TM5 inversions had negative biases in warm sectors, while GEOS-Chem inversions had positive biases in cold sectors. Possible reasons for the regional variability are discussed.

Gaudet, Brian J.↗

Snow Reconciles Observed and Simulated Phase Partitioning and Increases Cloud Feedback

The surprising increase of Earth's climate sensitivity in the most recent coupled model intercomparison project (CMIP) models has been largely attributed to extratropical cloud feedback, which is thought to be driven by greater supercooled water in present-day cloud phase partitioning (CPP). Here, we report that accounting for precipitation in the Goddard Institute for Space Studies ModelE3 radiation scheme, neglected in more than 60% of CMIP6 and 90% of CMIP5 models, systematically changes its apparent CPP and substantially increases its cloud feedback, consistent with results using CMIP models. Including precipitation in the comparison with cloud–aerosol lidar and infrared pathfinder satellite observations (CALIPSO) measurements and in model radiation schemes is essential to faithfully constrain cloud amount and phase partitioning, and simulate cloud feedbacks. Our findings suggest that making radiation schemes precipitation-aware (missing in most CMIP6 models) should strengthen their positive cloud feedback and further increase their already high mean climate sensitivity.

54 ENVIRONMENTAL SCIENCES↗

Evaluating Tropical Precipitation Relations in CMIP6 Models with ARM Data

Abstract A set of diagnostics based on simple, statistical relationships between precipitation and the thermodynamic environment in observations is implemented to assess phase 6 of the Coupled Model Intercomparison Project (CMIP6) model behavior with respect to precipitation. Observational data from the Atmospheric Radiation Measurement (ARM) permanent field observational sites are augmented with satellite observations of precipitation and temperature as an observational baseline. A robust relationship across observational datasets between column water vapor (CWV) and precipitation, in which conditionally averaged precipitation exhibits a sharp pickup at some critical CWV value, provides a useful convective onset diagnostic for climate model comparison. While a few models reproduce an appropriate precipitation pickup, most models begin their pickup at too low CWV and the increase in precipitation with increasing CWV is too weak. Convective transition statistics compiled in column relative humidity (CRH) partially compensate for model temperature biases—although imperfectly since the temperature dependence is more complex than that of column saturation. Significant errors remain in individual models and weak pickups are generally not improved. The conditional-average precipitation as a function of CRH can be decomposed into the product of the probability of raining and mean precipitation during raining times (conditional intensity). The pickup behavior is primarily dependent on the probability of raining near the transition and on the conditional intensity at higher CRH. Most models roughly capture the CRH dependence of these two factors. However, compensating biases often occur: model conditional intensity that is too low at a given CRH is compensated in part by excessive probability of precipitation.

54 ENVIRONMENTAL SCIENCES↗

Interdisciplinary Workshop on the Impacts of Land Use and Land Management on Earth System Evolution, Biogeochemical Cycles, Extremes, and Inter-Sectoral Dynamics (Final Report)

The Aspen Global Change Institute (AGCI) convened a week-long interdisciplinary science workshop from September 15-20, 2019 entitled, “Impacts of Land Use and Land Management on Earth System Evolution, Biogeochemical Cycles, Extremes, and Inter-Sectoral Dynamics.” The purpose of the workshop was to advance research on the impacts of land use and land management on climate, carbon emissions, and extremes, as well as the impact of land-use change on important societal issues such as food, water, and energy security. The workshop was co-chaired by Kate Calvin (PNNL), George Hurtt (University of Maryland), and Dave Lawrence (National Center for Atmospheric Research). Co-funders of the workshop included NASA, NOAA, and travel support for European colleagues through CRESCENDO. This workshop brought together 31 professionals with expertise across biogeochemistry, climate impacts, mitigation, land use datasets, and societal impacts to coordinate on upcoming research and publications on the role of biogeochemical and biophysical properties of the Earth surface and their implications for climate, as conducted through the Land Use model Intercomparison Project (LUMIP), which is part of Phase 6 of the Coupled Model Intercomparison Project (CMIP6). A primary outcome of the meeting was the definition and planning of a series of LUMIP papers that will rely on the LUMIP and CMIP6 simulations. By the end of the meeting more than twenty papers and paper lead authors had been identified and scoped out. Around ten of these papers were submitted by the December 31, 2019 IPCC deadline.

54 ENVIRONMENTAL SCIENCES↗

Relative Impacts of Mitigation, Temperature, and Precipitation on 21st-Century Megadrought Risk in the American Southwest

Megadroughts are comparable in severity to the worst droughts of the 20th century but are of much longer duration. A megadrought in the American Southwest would impose unprecedented stress on the limited water resources of the area, making it critical to evaluate future risks not only under different climate change mitigation scenarios but also for different aspects of regional hydroclimate. We find that changes in the mean hydroclimate state, rather than its variability, determine megadrought risk in the American Southwest. Estimates of megadrought probabilities based on precipitation alone tend to underestimate risk. Furthermore, business-as-usual emissions of greenhouse gases will drive regional warming and drying, regardless of large precipitation uncertainties. We find that regional temperature increases alone push megadrought risk above 70, 90, or 99% by the end of the century, even if precipitation increases moderately, does not change, or decreases, respectively. Although each possibility is supported by some climate model simulations, the latter is the most common outcome for the American Southwest in Coupled Model Intercomparison 5 generation models. An aggressive reduction in global greenhouse gas emissions cuts megadrought risks nearly in half.

risk↗