Search NASA⌕ Search

SEARCH · Search NASA

Results for “Global Reference Atmospheric Models”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

End-of-Century Changes in Orographic Precipitation with the Intermediate Complexity Atmospheric Research Model over the Western United States

Abstract Downscaled precipitation projections were created using the Intermediate Complexity Atmospheric Research (ICAR) model over the western United States to increase the physical realism in orographic precipitation changes. End-of-century simulations from eight models in phase 5 of the Coupled Model Intercomparison Project (CMIP5) were downscaled with ICAR and compared to the widely utilized statistically downscaled dataset, localized constructed analogs (LOCAs), to understand where and why projections of cool-season (September–May) precipitation differed. ICAR and LOCA precipitation projections were similar, but their sign differed in hydrologically relevant regions likely due to ICAR’s simulation of microphysics and mesoscale dynamics with high-resolution topography (6 km). In the Pacific Northwest, cool-season precipitation projections from ICAR showed an increase on the windward side of the Cascades and no significant change within the lee. This difference between the windward and leeward side was attributed to reduced zonal wind speeds, allowing more time for microphysical processes within ICAR. This contrast is enhanced by rain’s faster fall speed compared to snow, limiting transport into the lee. Meanwhile, LOCA projected an increase in precipitation across the Cascades. In the Upper Colorado River basin, LOCA projected an increase in precipitation in high elevation regions (>3000 m), but ICAR projected no significant change or a decrease in precipitation. High elevation differences were most evident in the spring and fall and were also attributed to a snow-to-rain transition and dynamical processes that impacted orographic enhancement within ICAR. Idealized, controlled studies are needed to better isolate individual processes, but these results underscore the importance of including microphysics and mesoscale dynamics within regional-scale precipitation projections. Significance Statement A set of global climate model simulations was downscaled using an atmospheric model that contains key physical equations, referred to as Intermediate Complexity Atmospheric Research (ICAR). ICAR was used to examine projected changes in end-of-century cool-season precipitation over mountains in the western United States. Precipitation projections from ICAR were similar to projections that used statistical relationships to downscale climate projections. However, projections differed between ICAR and statistically downscaled datasets in whether they increased, decreased, or stayed the same in specific, hydrologically relevant regions such as the eastern Cascades and high elevation areas of the Upper Colorado River basin. These differences were attributed to the simulation of physical processes in ICAR. The results highlight the importance of kilometer-scale atmospheric processes in regional climate projections.

Currier, William Ryan [NOAA/Physical Sciences Labo↗
DOE OSTI2024

Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations

Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Investigating the Atmospheric Sources and Sinks of Perfluorooctanoic Acid Using a Global Chemistry Transport Model

Perfluorooctanoic acid, PFOA, is one of the many concerning pollutants in our atmosphere; it is highly resistant to environmental degradation processes, which enables it to accumulate biologically. With direct routes of this chemical to the environment decreasing, as a consequence of the industrial phase out of PFOA, it has become more important to accurately model the effects of indirect production routes, such as environmental degradation of precursors; e.g., fluorotelomer alcohols (FTOHs). The study reported here investigates the chemistry, physical loss and transport of PFOA and its precursors, FTOHs, throughout the troposphere using a 3D global chemical transport model, STOCHEM-CRI. Furthermore, this investigation includes an important loss process of PFOA in the atmosphere via the addition of the stabilised Criegee intermediates, hereby referred to as the “Criegee Field.” Whilst reaction with Criegee intermediates is a significant atmospheric loss process of PFOA, it does not result in its permanent removal from the atmosphere. Additionally, the atmospheric fate of the resultant hydroperoxide product from the reaction of PFOA and Criegee intermediates resulted in a ≈0.04 Gg year-1 increase in the production flux of PFOA. Furthermore, the physical loss of the hydroperoxide product from the atmosphere (i.e., deposition), whilst decreasing the atmospheric concentration, is also likely to result in the reformation of PFOA in environmental aqueous phases, such as clouds, precipitation, oceans and lakes. As such, removal facilitated by the “Criegee Field” is likely to simply result in the acceleration of PFOA transfer to the surface (with an expected decrease in PFOA atmospheric lifetime of ≈10 h, on average from ca. 80 h without Criegee loss to 70 h with Criegee loss).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Technical descriptions of the experimental dynamical downscaling simulations over North America by the CAM–MPAS variable-resolution model

Abstract. Comprehensive assessment of climate datasets is important for communicating model projections and associated uncertainties to stakeholders. Uncertainties can arise not only from assumptions and biases within the model but also from external factors such as computational constraint and data processing. To understand sources of uncertainties in global variable-resolution (VR) dynamical downscaling, we produced a regional climate dataset using the Model for Prediction Across Scales (MPAS; dynamical core version 4.0) coupled to the Community Atmosphere Model (CAM; version 5.4), which we refer to as CAM–MPAS hereafter. This document provides technical details of the model configuration, simulations, computational requirements, post-processing, and data archive of the experimental CAM–MPAS downscaling data. The CAM–MPAS model is configured with VR meshes featuring higher resolutions over North America as well as quasi-uniform-resolution meshes across the globe. The dataset includes multiple uniform- (240 and 120 km) and variable-resolution (50–200, 25–100, and 12–46 km) simulations for both the present-day (1990–2010) and future (2080–2100) periods, closely following the protocol of the North American Coordinated Regional Climate Downscaling Experiment. A deviation from the protocol is the pseudo-warming experiment for the future period, using the ocean boundary conditions produced by adding the sea surface temperature and sea-ice changes from the low-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-LR) in the Coupled Model Intercomparison Project Phase 5 to the present-day ocean state from a reanalysis product. Some unique aspects of global VR models are evaluated to provide background knowledge to data users and to explore good practices for modelers who use VR models for regional downscaling. In the coarse-resolution domain, strong resolution sensitivity of the hydrological cycles exists over the tropics but does not appear to affect the midlatitude circulations in the Northern Hemisphere, including the downscaling target of North America. The pseudo-warming experiment leads to similar responses of large-scale circulations to the imposed radiative and boundary forcings in the CAM–MPAS and MPI-ESM-LR models, but their climatological states in the historical period differ over various regions, including North America. Such differences are carried to the future period, suggesting the importance of the base state climatology. Within the refined domain, precipitation statistics improve with higher resolutions, and such statistical inference is verified to be negligibly influenced by horizontal remapping during post-processing. Limited (≈50 % slower) throughput of the current code is found on a recent many-core/wide-vector high-performance computing system, which limits the lengths of the 12–46 km simulations and indirectly affects sampling uncertainty. Our experience shows that global and technical aspects of the VR downscaling framework require further investigations to reduce uncertainties for regional climate projection.

58 GEOSCIENCES↗
DOE OSTI2026

Convective Biases in the US DOE Global Storm‐Resolving Model: Insights From Regionally Refined Simulations During the CACTI Campaign

Accurately simulating convective processes in complex terrain remains a critical challenge for global storm-resolving models (GSRMs). This study systematically evaluates moist convective biases in the Regionally Refined Mesh configuration of the U.S. Department of Energy Simple Cloud-Resolving E3SM Atmosphere Model (RRM-SCREAM) using comprehensive observations and large-eddy simulations from the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) campaign in the mountainous area of central Argentina. Comparisons of simulations with high-resolution observations and reanalysis data indicate that RRM-SCREAM effectively captures large-scale meteorological patterns, including regional atmospheric gradients and diurnal variability. However, RRM-SCREAM disproportionately produces smaller precipitation clusters referred to as “popcorn convection,” and exaggerated rainfall intensities compared to observations and reference models. Detailed examination of a representative orographic shallow-to-deep convective transition case shows that RRM-SCREAM delays initial shallow convection growth due to lower-tropospheric dryness and sustained convective inhibition, but once triggered, deep convection becomes overly vigorous with excessively strong vertical velocities and elevated cloud ice content, linked to a thermodynamic structure characterized by suppressed low-level moistening and excessive upper-level moisture retention. Our results highlight specific deficiencies in the model representation of convective vertical velocity, cloud microphysical processes, and convective precipitation organization within RRM-SCREAM. Addressing these biases is essential for improving the predictions of convective clouds and precipitation in the global high-resolution atmospheric models.

Su, Tianning [Lawrence Livermore National Laborato↗
DOE OSTI2024

Expanding the representation of aerosol, cloud, and precipitation processes with graph network-based simulators

We explored a novel framework for simulating the small-scale processes that drive the evolution of aerosol, cloud, and precipitation particles, which are a critical gap in the predictive understanding of weather and climate. Particle-based methods have emerged as an effective tool for modeling aerosol-cloud-precipitation interactions, but existing particle-based models are computationally too expensive to simulate the large domains relevant for the atmosphere or to represent the full suite of relevant processes. The lack of a comprehensive and efficient reference model is a critical bottleneck in our understanding of cloud and precipitation processes and our ability to parameterize these processes for regional- and global-scale simulations. To address this need, we explored an approach to accelerate and expand particle-based models using a new machine learning approach, graph network-based simulators (GNS). Rather than modeling the evolution of the system by numerically integrating continuity equations, the GNS represents dynamics through learned message passing. Our aim was to develop fast and accurate surrogate models for particle-based simulations. We explored applying GNS to simulate cloud droplet transport, growth, and evaporation under turbulent conditions, but we found the GNS over-smoothed the simulations. We then applied the GNS to simulate aerosol dynamics through gas condensation and found the GNS was able to reproduce the benchmark, physics-based simulation with high accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Responses in the Subpolar North Atlantic in Two Climate Model Sensitivity Experiments with Increased Stratospheric Aerosols

Abstract The subpolar North Atlantic (SPNA) shows contrasting responses in two sensitivity experiments with increased stratospheric aerosols, offering insight into the physical processes that may impact the Atlantic meridional overturning circulation (AMOC) in a warmer climate. In one, the upper ocean becomes warm and salty, but in the other it becomes cold and fresh. The changes are accompanied by diverging AMOC responses. The first experiment strengthens the AMOC, opposing the weakening trend in the reference simulation. The second experiment shows a much smaller impact. Both simulations use the Community Earth System Model with the Whole Atmosphere Community Climate Model component (CESM-WACCM) but differ in model versions and stratospheric aerosol specifications. Despite both experiments using similar approaches to increase stratospheric aerosols to counteract the rising global temperature, the contrasting SPNA and AMOC responses indicate a considerable dependency on model physics, climate states, and model responses to forcings. This study focuses on examining the physical processes involved with the impact of stratospheric aerosols on the SPNA salinity changes and their potential connections with the AMOC and the Arctic. We find that in both cases, increased stratospheric aerosols act to enhance the SPNA upper-ocean salinity by reducing freshwater export from the Arctic, which is closely tied to the Arctic sea ice changes. The impact on AMOC is primarily through the thermal component of the surface buoyancy fluxes, with negligible contributions from the freshwater component. These experiments shed light on the physical processes that dictate the important connections between the SPNA, the Arctic, the AMOC, and their subsequent feedbacks on the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Root Foraging Alters Global Patterns of Ecosystem Legacy From Climate Perturbations

The response of terrestrial ecosystems to climate perturbations typically persist longer than the timescale of the forcing, a phenomenon broadly referred to as legacy. Understanding the strength of legacy is critical for predicting ecosystem sensitivity to climate extremes and the extent that persistent changes in surface-atmosphere exchange might feedback onto climate. The cause of ecosystem legacy has been associated with myriad factors such as changes in aboveground biomass, however, few studies have tested how changes in the depth distribution of fine roots in response to perturbation might alter an ecosystem's recovery time. We explore this question using an Earth System Model that includes a dynamic root module where vegetation can forage for water and nutrients by altering their root profiles. The global simulations presented here show that in response to water stress events most ecosystems deepen their root profiles. In semi-arid ecosystems, this response expedites recovery (i.e., less legacy) relative to simulations without dynamics roots because access to deeper water pools after the initial event remains favorable. In wetter ecosystems, the development of deeper root profiles slows down the recovery timescale (i.e., more legacy) because the deeper root profile reduces access to nutrients and is thus unfavorable. Here, the simulations show that while the inclusion of dynamic roots might only minimally affect global patterns of Gross Primary Productivity and Evapotranspiration, the shift in root profile alters the timescale of recovery. Studies interested in the sustained response of land surfaces fluxes to climate disturbances should consider models that include dynamic root capability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The GFDL‐CM4X Climate Model Hierarchy, Part I: Model Description and Thermal Properties

We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of 0.25° (referred to as CM4X‐p25) and the other that uses a grid 0.125° (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of -0.02 W m -2 within $\mathcal{O}$ (150) years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly 400 ZJ less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughly 1100 ZJ more heat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer the mesoscale dominance hypothesis for why CM4X‐p125 shows such favorable thermal equilibration properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Prescribing the aerosol effective radiative forcing in the Simple Cloud-Resolving E3SM Atmosphere Model v1

Aerosol effective radiative forcing critically influences climate projections but remains poorly constrained. Using the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) v1 configuration, we assess the sensitivity of the radiative forcing due to anthropogenic aerosol changes using a simplified prescribed aerosol scheme (SPA) derived from E3SM v3. Nudged simulations at 3 and 12 km horizontal grid spacings reveal a more negative aerosol forcing than the reference 100 km E3SM v3 whence the SPA properties are derived. The resulting globally averaged aerosol forcing signal is largely due to aerosol–cloud interactions and exhibits little overall resolution sensitivity, while hints of resolution sensitivity appear regionally between the 3 and 12 km runs. While the default SPA scheme overestimates cloud droplet dependence on aerosols, parameterization adjustments in the activation process reconcile forcing estimates with the reference model. Our results demonstrate the ability to use a prescribed aerosol scheme to hold aerosol forcing to a desired strength across resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Causes of and Solutions to Wind Speed Bias in NREL's 2020 Offshore Wind Resource Assessment for the California Pacific Outer Continental Shelf

This report provides the results of a detailed analysis into the causes of high wind speed bias in the 20-year wind resource data set for offshore California the National Renewable Energy Laboratory (NREL) released in 2020, herein called CA20. The data set was developed using the state-of-the-art Weather Research and Forecasting (WRF) model. Notably, no floating lidars were available at the time in offshore California to validate offshore hub-height wind speeds. In late 2020, the Pacific Northwest National Laboratory (PNNL) deployed two floating lidars in the California outer continental shelf (OCS), near the Bureau of Ocean Energy Management (BOEM) call areas of Humboldt and Morro Bay. Using these observations through 2021, NREL found considerable bias in modeled hub-height winds at both locations: up to +2 m/s at Humboldt over a 6-month period, and up to +1 m/s at Morro Bay over a one-year period. Upon the discovery of this bias, the Department of Energy (DOE) and BOEM funded NREL and PNNL to investigate the causes of, impacts of, and solutions to the bias in the CA20 data set. This report summarizes the findings of this research. We first investigated whether different WRF model setups could lead to reduced bias. We found that the choice of planetary boundary layer (PBL) scheme - which controls the vertical turbulent mixing of momentum, heat, and moisture in the lowermost part of the atmosphere - greatly affected hub-height wind speeds in the region. Specifically, switching from the Mellor-Yamada-Nakanishi-Niino (MYNN) scheme used in CA20 (and widely used across a range of operational and research weather models) to the less common Yonsei University (YSU) scheme nearly eliminated the bias at both the Humboldt and Morro Bay lidar locations. The large discrepancy between the MYNN- and YSU-modeled hub-height winds pointed towards the role of atmospheric stability. In general, PBL schemes agree well in conditions of high turbulence and mixing, normally referred to as "unstable" conditions. By contrast, PBL schemes start to diverge in "stable" conditions, where turbulence is low and thermal stratification (i.e., higher temperature air sitting on top of colder air) greatly suppresses vertical mixing. Under such conditions, winds aloft can decouple from surface effects and greatly accelerate, causing high wind speeds at hub-height and frequent low-level jets (LLJs). We determined that these stable conditions are in fact dominant in offshore California. The region is characterized by moderate-to-extreme stable stratification with a LLJ on average around 200 meters above sea-level. To our knowledge, no wind energy area globally has as strongly stable stratification as offshore California. Under these extreme conditions, we determined that the MYNN scheme models higher stability than YSU, resulting in less vertical turbulent mixing than YSU, allowing for the acceleration of hub-height winds, more intense LLJs, and higher-amplitude inertial oscillations. Using surface observations, we found that MYNN overestimates near-surface stability, whereas YSU tends to model stability better. We then considered several short-term case studies to assess additional meteorological drivers of the bias at Humboldt. We found that during synoptic scale northerly flows driven by the North Pacific High and inland thermal low, a coastal warm bias in the MYNN case studies contributes to the modeled wind speed bias by altering the boundary layer thermodynamics via a thermal wind mechanism. Given the strong performance of the YSU-based runs in offshore California, NREL has produced and published an updated version of the CA20 data set with YSU as the PBL scheme. This updated data set is now part of NREL's 2023 National Offshore Wind (NOW-23) data set, which covers all the U.S. offshore waters. The development and final validation of the NOW-23 data set in offshore California is documented in this report.

17 WIND ENERGY↗
DOE OSTI2023

A Strengthened Teleconnection of the Quasi‐Biennial Oscillation and Tropical Easterly Jet in the Past Decades in E3SMv1

The Energy Exascale Earth System Model version 1 (E3SMv1) is a new global model that can generate a weakened tropical easterly jet (TEJ) during the easterly quasi-biennial oscillation (EQBO) and a strengthened TEJ during the westerly QBO (WQBO), which is referred to as the QBO-TEJ teleconnection. Although the longitudinal structure of the QBO-TEJ teleconnection is different between reanalysis and simulations, we hypothesize that simulations can still provide insights into the observed long-term changes in the QBO-TEJ teleconnection. The Atmospheric Model Intercomparison Project historical simulations (AMIP-historical) generate a strengthened trend in the correlation between the QBO and TEJ indices in the period of 1950–2014 as in reanalysis. The comparison of the AMIP-historical and the Coupled Model Intercomparison Project coupled historical simulations (coupled-historical) shows that the warming of the Indian Ocean has played a vital role in the strengthened QBO-TEJ teleconnection over the past six decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Accuracy of Vaisala RS41 and RS92 Upper Tropospheric Humidity Compared to Satellite Hyperspectral Infrared Measurements

Radiosondes are important for calibrating satellite sensors and assessing sounding retrievals. Vaisala RS41 radiosondes have mostly replaced RS92 in the Global Climate Observing System (GCOS) Reference Upper Air Network (GRUAN) and the conventional network. This study assesses RS41 and RS92 upper tropospheric humidity (UTH) accuracy by comparing with Infrared Atmospheric Sounding Interferometer (IASI) upper tropospheric water vapor absorption spectrum measurements. Using single RS41 and RS92 soundings at three GRUAN and DOE Atmospheric Radiation Measurement (ARM) sites and dual RS92/RS41 launches at three additional GRUAN sites, collocated with cloud-free IASI radiances (OBS), we compute Line-by-Line Radiative Transfer Model radiances for radiosonde profiles (CAL). We analyze OBS-CAL differences from 2015 to 2020, for daytime, nighttime, and dusk/dawn separately if data is available, for standard (STD) RS92 and RS41 processing, and RS92 GRUAN Data Processing (GDP; RS41 GDP is in development). We find that daytime RS41 (even without GDP) has ~1% smaller UTH errors than GDP RS92. RS41 may still have a dry bias of 1–1.5% for both daytime and nighttime, and a similar error for nighttime RS92 GDP, while standard RS92 may have a dry bias of 3–4%. These sonde humidity biases are probably upper limits since “cloud-free” scenes could still be cloud contaminated. Radiances computed from European Centre for Medium-Range Weather Forecasts (ECMWF) analyses match better than radiosondes with IASI measurements, perhaps because ECMWF assimilates IASI measurements. Relative differences between RS41 STD and RS92 GDP, or between radiosondes and ECMWF humidity profiles obtained from the radiance analysis, are consistent with their differences obtained directly from the RH measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantifying Carbon Cycle Extremes and Attributing Their Causes Under Climate and Land Use and Land Cover Change From 1850 to 2300

The increasing atmospheric carbon dioxide (CO 2 ) mole fraction affects global climate through radiative (trapping longwave radiation) and physiological effects (reduction of plant transpiration). We use the simulations of the Community Earth System Model (CESM1-BGC) forced with Representative Concentration Pathway 8.5 to investigate climate-vegetation feedbacks from 1850 to the year 2300. Human-induced land use and land cover change (LULCC), through biogeochemical and biogeophysical processes, alter the climate and modify photosynthetic activity. The changing characteristics of extreme anomalies in photosynthesis, referred to as carbon cycle extremes, increase the uncertainty of terrestrial ecosystems to act as a net carbon sink. However, the role of LULCC in altering carbon cycle extremes under business-as-usual (continuously rising) CO 2 emissions is unknown. Here we show that LULCC magnifies the intensity, frequency, and extent of carbon cycle extremes, resulting in a net reduction in expected photosynthetic activity in the future. We found that large temporally contiguous negative carbon cycle extremes are due to a persistent decrease in soil moisture, which is triggered by declines in precipitation. With LULCC and global warming, vegetation exhibits increased vulnerability to hot and dry environmental conditions, increasing the frequency of fire events and resulting in considerable losses in photosynthetic activity. While most regions show strengthening of negative carbon cycle extremes, a few locations show a weakening effect driven by declining vegetation cover or benign climate conditions for photosynthesis. Increasing hot, dry, and fire-driven carbon cycle extremes are essential for improving carbon cycle modeling and estimation of ecosystem responses to LULCC and rising CO 2 mole fractions. Moreover, large aberrations in vegetation productivity represent potential and growing threats to human lives, wildlife, and food security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding Recent Global Hydroclimate Change Using Multivariate Detection and Attribution Techniques and GCM Experiments

The observed slowdown in the rate of global average surface temperature (T) warming during the post-1998 decade (referred to at the time as the “global warming hiatus”) posed a challenge for climate researchers because it was not captured in global climate models projection. Observational and modeling evidence hinted that that leveling of the T trend was linked with the simultaneous appearance of cold surface temperature in the eastern equatorial Pacific that formed a protracted La Niña-like state, a phenomenon known as the Inter-decadal Pacific Oscillation (IPO). It was however not clear whether the concurrent equatorial Pacific cooling was a regional response to the rise in atmospheric greenhouse gas levels or anthropogenic aerosols that is not captured by the climate models or whether it was an unrelated, serendipitous manifestation of a natural phenomenon, internal to the climate system. To resolve this issue, we proposed in this project to develop and apply fingerprinting methods to differentiate between the observed climate response to various external forcing agents and the manifestations of internal climate variability. The project also sought to find out if climate models are able to reproduce the change signals detected in key observational phenomena and if they can be used to investigate the underlying processes. As such we examined changes in global hydroclimate variability and in an outstanding example of regional hydroclimate change, namely the multi-decadal 20 th century change in Sahel rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Historical and Future Windstorms in the Northeastern United States

Large-scale windstorms represent an important atmospheric hazard in the Northeastern US (NE) and are associated with substantial socioeconomic losses. Regional simulations performed with the Weather Research and Forecasting (WRF) model using lateral boundary conditions from three Earth System Models (ESMs: Geophysical Fluid Dynamics Laboratory (GFDL), Hadley Centre Global Environment Model (HadGEM) and Max Planck Institute (MPI)) are used to quantify possible future changes in windstorm characteristics and/or changes in the parent cyclone types responsible for windstorms. WRF nested within MPI ESM best represents important aspects of historical windstorms and the cyclone types responsible for generating windstorms compared with a reference simulation performed with the ERA-Interim reanalysis for the historical climate. The spatial scale and frequency of the largest windstorms in each simulation defined using the greatest extent of exceedance of local 99.9th percentile wind speeds (U > U999) plus 50-year return period wind speeds (U50,RP) do not exhibit secular trends. Projections of extreme wind speeds and windstorm intensity/frequency/geolocation and dominant parent cyclone type associated with windstorms vary markedly across the simulations. Only the MPI nested simulations indicate statistically significant differences in windstorm spatial scale, frequency and intensity over the NE in the future and historical periods. This model chain, which also exhibits the highest fidelity in the historical climate, yields evidence of future increases in 99.9th percentile 10 m height wind speeds, the frequency of simultaneous U > U999 over a substantial fraction (5–25%) of the NE and the frequency of maximum wind speeds above 22.5 ms−1. These geophysical changes, coupled with a projected doubling of population, leads to a projected tripling of a socioeconomic loss index, and hence risk to human systems, from future windstorms.

Pryor, Sara C. (ORCID:0000000348473440)↗
DOE OSTI2023

A Consistent Representation of Cloud Overlap and Cloud Subgrid Vertical Heterogeneity

Many global climate models underestimate the cloud cover and overestimate the cloud albedo, especially for low-level clouds. We determine how a correct representation of the vertical structure of clouds can fix part of this bias. We use the 1D McICA framework and focus on low-level clouds. Using Large Eddy Simulations results as reference, we propose a method based on exponential-random overlap that represents the cloud overlap between layers and the subgrid cloud properties over several vertical scales, with a single value of the overlap parameter. Starting from a coarse vertical grid, representative of atmospheric models, this algorithm is used to generate the vertical profile of the cloud fraction with a finer vertical resolution, or to generate it on the coarse grid but with subgrid heterogeneity and cloud overlap that ensures a correct cloud cover. Doing so we find decorrelation lengths are dependent on the vertical resolution, except if the vertical subgrid heterogeneity and interlayer overlap are taken into account coherently. We confirm that the frequently used maximum-random overlap leads to a significant error by underestimating the low-level cloud cover with a relative error of about 50%, that can lead to an error of SW cloud albedo as big as 70%. Not taking into account the subgrid vertical heterogeneity of clouds can cause a relative error of 20% in brightness, assuming the cloud cover is correct.

54 ENVIRONMENTAL SCIENCES↗