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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Country-level land carbon sink and its causing components by the middle of the twenty-first century

Abstract Background Countries have long been making efforts by reducing greenhouse-gas emissions to mitigate climate change. In the agreements of the United Nations Framework Convention on Climate Change, involved countries have committed to reduction targets. However, carbon (C) sink and its involving processes by natural ecosystems remain difficult to quantify. Methods Using a transient traceability framework, we estimated country-level land C sink and its causing components by 2050 simulated by 12 Earth System Models involved in the Coupled Model Intercomparison Project Phase 5 (CMIP5) under RCP8.5. Results The top 20 countries with highest C sink have the potential to sequester 62 Pg C in total, among which, Russia, Canada, USA, China, and Brazil sequester the most. This C sink consists of four components: production-driven change, turnover-driven change, change in instantaneous C storage potential, and interaction between production-driven change and turnover-driven change. The four components account for 49.5%, 28.1%, 14.5%, and 7.9% of the land C sink, respectively. Conclusion The model-based estimates highlight that land C sink potentially offsets a substantial proportion of greenhouse-gas emissions, especially for countries where net primary production (NPP) likely increases substantially and inherent residence time elongates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

Abstract. This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of CMIP6 Streamflow in the Arctic

Earth system models (ESMs) serve as the primary basis for projecting future streamflow changes, but they are biased in terms of the their ability to reproduce historical observations of streamflow. Rigorous downscaling and bias-correction procedures, which are time-consuming and introduce uncertainties, are commonly used to address these biases in streamflow or runoff. However, the applicable limits of streamflow projections remain largely untested through direct comparisons with observations. In this work, we compare historical streamflow time series observations with a suite of Coupled Model Intercomparison Project phase 6 (CMIP6) ESMs to assess the space and time scales at which each ESM represents historic streamflow observations. We focus our analysis on streamflow across the Arctic given the higher degree of warming experienced over this region, underscoring its importance for understanding the future changes. A series of metrics—volume, seasonality, extreme capturing streamflow variability, drought, and flooding is tested to assess ESM simulation skill for total events and overall distributions. Our results indicate that although improvements are necessary in ESM runoff and streamflow projections, the use of temporal averaging and targeting the best-fit model for comparison to observations when a suite of ESMs is considered can provide robust streamflow projections. Moreover, this work provides a basis for contextualizing the limits of ESM applicability for studies aiming to assess streamflow in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Blocking Simulations in GFDL GCMs for CMIP5 and CMIP6

Abstract The frequency of atmospheric blocking has been largely underestimated by general circulation models (GCMs) participating in the Coupled Model Intercomparison Project (CMIP). Errors in the onset, persistence, barotropicity, geographical preference, seasonality, intensity, and moving speed of global blocking were diagnosed in 10 Geophysical Fluid Dynamics Laboratory (GFDL) GCMs for recent CMIP5 and CMIP6 using a detection approach that combines zonal eddies and the reversal of zonal winds. The blocking frequency, similar at 500 and 250 hPa, is underestimated by 50% in the Atlantic–Europe region during December–February but is overestimated by 60% in the Pacific–North America region during that season and by 70% in the southwest Pacific during July–August. These blocking biases at 500 hPa were investigated in the five CMIP6 models that showed improvements over the CMIP5 versions. The Atlantic–Europe underestimate corresponds to lower instantaneous blocking rates, lower persistent blocking rates, and higher persistent stationary ridge rates; the number of blocks with a duration of 4–5 days is only 40%–65% of that in observations. In contrast, the overestimate consists of excessive blocks with a duration longer than 12 days in the Pacific–North America and up to twice as many 4–6-day events in the southwest Pacific. Simulated December–February blocks up to 12 days in the Pacific–North America region tend to be stronger and to move more slowly than those in observations. Diagnostic sensitivity tests indicated that the zonal mean and zonal eddy components of the mean state play a key role, as replacing each with that of observations substantially reduced many of the outstanding biases in these GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Improving Projections of Antarctic Ice Sheet Contribution to Sea‐Level Change Through 2300 by Capturing Gravitational, Rotational, and Deformational Effects

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

58 GEOSCIENCES↗
DOE OSTI2025

Integrated Modeling Driven Evaluation of Opportunities for Climate‐Resilient Perennial Biomass Crop Plantings in Flood‐Prone Agricultural Landscapes

Adapting to future climate change in flood-prone landscapes will require climate-resilient agricultural systems. Planting perennial crops, like switchgrass and willow, along river corridors can mitigate future flooding while supporting bioenergy markets. We developed an integrated assessment linking climate, hydrologic, and inundation model results to assess future flood risk to river-adjacent agricultural lands in the Mid-Atlantic Region (MAR) and explore this opportunity. We produced ensemble streamflow projections for every MAR stream using a hydrologic model driven by a suite of downscaled and bias-corrected Coupled Model Intercomparison Project Phase 6 climate projections. We then conducted high-resolution inundation mapping based on projected flood frequencies for baseline and future periods. Results show that in the near-term future, at least two-thirds of the streams will experience 100-year floods more severe than the baseline 200-year floods. Riparian zones are projected to face a median rise of inundation by 9.5%–24.1%. Results show that there is an opportunity to mitigate flooding in over half of MAR's counties with the quantities of switchgrass and willow plantings anticipated for mature bioenergy markets, even under the most extreme (200-year) flood events. Our integrated modeling framework can guide similar regions to evaluate opportunities for flood-resilient agricultural systems under climate change.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

The Southern hemisphere midlatitude circulation response to rapid adjustments and sea surface temperature driven feedbacks

Rapid adjustments – the response of meteorology to external forcing while sea surface temperatures (SST) and sea ice are held fixed – can affect the midlatitude circulation and contribute to long-term forced circulation responses in climate simulations. This study examines rapid adjustments in the Southern hemisphere (SH) circulation using nine models from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), which perform fixed SST and coupled ocean experiments for five perturbations: a doubling of carbon dioxide (2xCO2), a tripling of methane (3xCH4), a five-fold increase in sulphate aerosol (5xSO4), a ten-fold increase in black carbon aerosol (10xBC) and a 2% increase in solar constant (2%Sol). In the coupled experiments, the SH eddy-driven jet shifts poleward and strengthens for forcings that produce global warming (and vice versa for 5xSO4), with the strongest response found in austral summer. In austral winter, the responses project more strongly onto a change in jet strength. For 10xBC, which induces strong shortwave absorption, the multi-model mean (MMM) rapid adjustment in DJF jet latitude is ~75% of the change in the coupled simulations. For the other forcings, which induce larger SST changes, the effect of SST-mediated feedbacks on the SH circulation is larger than the rapid adjustment. Nevertheless, for these perturbations the magnitude of the MMM jet shift due to the rapid adjustment is still around 25-33% of that in the coupled experiments. The results demonstrate the need to understand the mechanisms for rapid adjustments in the midlatitude circulation, in addition to the effect of changing SSTs.

Wood, Thomas↗
DOE OSTI2023

Evaluating Enhanced Reservoir Evaporation Losses From CMIP6‐Based Future Projections in the Contiguous United States

Abstract Enhanced reservoir evaporation has become an emerging concern regarding water loss, especially when compounded with the ever‐increasing water demand. In this study, we evaluated the evaporation rates and losses for 678 major reservoirs (representing nearly 90% of total storage capacity) in the Contiguous United States over historical baseline (1980–2019), near‐term (2020–2039), and mid‐term (2040–2059) future periods. The evaporation rate was estimated using the Lake Evaporation Model (LEM), an advanced lake evaporation model that addresses both heat storage and fetch effects, driven by multi‐ensemble downscaled Coupled Model Intercomparison Projects 6 (CMIP6) projections under the SSP585 emission scenario. The results project that the evaporation loss may increase by 2.5 × 10 7 m 3 /yr through the research period (1980–2059). Among all regions, the Rio Grande is projected to have the largest increasing rate in the near‐term and mid‐term future, with values of 7.11% of 10.25%, respectively. At the seasonal scale, the most significant increase in the evaporation rate is projected during the fall. The evaporation is projected to increase faster than the streamflow over many of the regions in the southwestern US during the summer/fall, suggesting that the shortage of water will be further exacerbated. The climate models contribute the most to the variance, as compared to the other components related to the projection of evaporation losses (e.g., hydrological model, downscaling method, and historical meteorological data set). These findings demonstrate the need to consider accelerated water loss through open water evaporation in long‐term water resources planning across various spatiotemporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

On the Correspondence Between Atmosphere-Only and Coupled Simulations for Radiative Feedbacks and Forcing From CO 2

Atmosphere-only experiments are widely used to investigate climate feedbacks simulated in more computationally expensive fully coupled global climate model simulations. We confirm that this remains a valid approach by comparing the radiative feedbacks and forcing between coupled and atmosphere-only simulations for the latest models taking part in the 6th phase of the Coupled Model Intercomparison Project (CMIP6). For global-mean cloud feedbacks, we find a better than previously known correspondence between these experiments, which applies even to the response of individual cloud properties (amount, altitude, and optical depth) and holds even when considering atmosphere-only simulations of only 1 yr duration. For regional cloud feedbacks, the correspondence between the two experiments is generally present at every geographic location except for the tropical Pacific but takes longer experiments to reveal. For the lapse rate and surface albedo feedbacks, the correspondence between the two experiments is weaker due to the non-uniform warming pattern and loss of sea ice in the coupled experiment. For the across-model relationship between 4xCO 2 effective radiative forcing and feedback, we find a different behavior across experiments in CMIP6 than in CMIP5, casting doubt on the physical significance of previous results that highlighted an anti-correlation between the two quantities. Overall, these results confirm the utility of atmosphere-only experiments particularly to study cloud feedbacks, which are the dominant source of inter-model spread in climate sensitivity.

58 GEOSCIENCES↗
DOE OSTI2025

High-Resolution Fire Weather Index Data for the Conterminous US (1980–2099), Version 1

This dataset presents a suite of high-resolution fire weather index datasets calculated from observation (gridMet, Livneh, Daymet V4), reanalysis (AgERA5), downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Aside from the daily FWI datasets, we also include a set of FWI extreme indicators at annual, seasonal, and monthly scales, including 1) fwixx: maximum FWI; 2) fwisa: mean FWI. For annual fwisa, it refers to the season with the maximum seasonal average; 3) fwils: Length of fire season over a specified period, where fire season is defined as the days exceeding the median value of the normalized FWI during the reference period (1980-1984); 4) fwixd: Number of extreme fire weather days over a specified period, where extreme day is defined as the day with FWI > the 95th percentile of the FWI during the reference period (1980-1984). All FWI datasets cover 1980-2020 baseline and the model simulated products including the downscaled products additionally include 2021-2099 near-future periods under the high-end (SSP585) emission scenario.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Model spread in tropical low cloud feedback tied to overturning circulation response to warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of Global Climate Models for Use in Energy Analysis

The interplay between energy, climate, and weather is becoming more complex due to increasing contributions of renewable energy generation, energy storage, electrified end uses, and the increasing frequency of extreme weather events. Energy system analyses commonly rely on meteorological inputs to estimate renewable energy generation and energy demand; however, these inputs rarely represent the estimated impacts of future climate change. Climate models and publicly available climate change datasets can be used for this purpose, but the selection of inputs from the myriad of available models and datasets is a nuanced and subjective process. In this work, we assess datasets from various global climate models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). We present evaluations of their skills with respect to the historical climate and comparisons of their future projections of climate change for two climate change scenarios. We present the results for different climatic and energy system regions and include interactive figures in the accompanying software repository. Previous work has presented similar GCM evaluations, but none have presented variables and metrics specifically intended for comprehensive energy systems analysis including impacts on energy demand, thermal cooling, hydropower, water availability, solar energy generation, and wind energy generation. We focus on GCM output meteorological variables that directly affect these energy system components including the representation of extreme values that can drive grid resilience events. The objective of this work is not to recommend the best climate model and dataset for a given analysis, but instead to provide a reference to facilitate the selection of climate models and scenarios in subsequent work.

14 SOLAR ENERGY↗
DOE OSTI2022

Warming Pattern over the Northern Hemisphere Midlatitudes in Boreal Summer 1979–2020

Abstract Significant surface air temperature warming during summer 1979–2020 is not uniformly distributed in the northern midlatitudes over land but rather is confined to several longitudinal sectors including Europe, central Siberia and Mongolia, and both coasts of North America. These hot spots are accompanied by a chain of high pressure ridges from an anomalous, circumglobal Rossby wave train in the upper troposphere. From reanalysis data and several baseline experiments from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we find that the circulation trend pattern is associated with fluctuations of the Atlantic multidecadal variability (AMV) and the interdecadal Pacific oscillation. The phase shift of AMV in the 1990s is particularly noteworthy for accelerating warming averaged over the northern midlatitude land. The amplitude of the observed trend in both surface air temperature and the upper-level geopotential height generally falls beyond the range of multidecadal trends simulated by the CMIP6 preindustrial control runs, supporting the likelihood that anthropogenic forcing played a critical role in the observed trend. On the other hand, the fidelity of the simulated low-frequency modes of variability and their teleconnections, especially on multidecadal time scales, is difficult to assess because of the relatively short observational records. Our mechanistic modeling results indicate that synoptic eddy–mean flow interaction is a key to the formation of the anomalous wave train but how the multidecadal modes can modulate the synoptic eddies through atmosphere–ocean and atmosphere–land interactions remains poorly understood. This gap in our knowledge makes it challenging to quantify the roles of the low-frequency modes and external forcings in causing the observed multidecadal trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The GFDL Variable-Resolution Global Chemistry-Climate Model for Research at the Nexus of US Climate and Air Quality Extremes

We present a variable-resolution global chemistry-climate model (AM4VR) developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) for research at the nexus of US climate and air quality extremes. AM4VR has a horizontal resolution of 13 km over the US, allowing it to resolve urban-to-rural chemical regimes, mesoscale convective systems, and land-surface heterogeneity. With the resolution gradually reducing to 100 km over the Indian Ocean, we achieve multi-decadal simulations driven by observed sea surface temperatures at 50% of the computational cost for a 25-km uniform-resolution grid. In contrast with GFDL's AM4.1 contributing to the sixth Coupled Model Intercomparison Project at 100 km resolution, AM4VR features much improved US climate mean patterns and variability. In particular, AM4VR shows improved representation of: precipitation seasonal-to-diurnal cycles and extremes, notably reducing the central US dry-and-warm bias; western US snowpack and summer drought, with implications for wildfires; and the North American monsoon, affecting dust storms. AM4VR exhibits excellent representation of winter precipitation, summer drought, and air pollution meteorology in California with complex terrain, enabling skillful prediction of both extreme summer ozone pollution and winter haze events in the Central Valley. AM4VR also provides vast improvements in the process-level representations of biogenic volatile organic compound emissions, interactive dust emissions from land, and removal of air pollutants by terrestrial ecosystems. We highlight the value of increased model resolution in representing climate–air quality interactions through land-biosphere feedbacks. AM4VR offers a novel opportunity to study global dimensions to US air quality, especially the role of Earth system feedbacks in a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modular hybrid modeling to increase efficiency, explore structural uncertainty, and allow multidimensional complexity scaling in land surface models

Land surface models (LSMs) are indispensable tools for predicting hydrologic extremes, as well as a particularly uncertain component of Earth system models that has stubbornly resisted the convergence in projections over several successive generations of model intercomparisons. This uncertainty in LSMs is poorly quantified and poorly attributed to specific processes, which has hampered efforts to focus research in reducing uncertainty. This has resulted from sparse sampling of the possible uncertainty space—which is high-dimensional and has contributions from parametric, structural, initial, and boundary condition uncertainties—as an artifact of CMIP-type ensembles of opportunity and limitations inherent in observational benchmarks. A new approach is needed to understand and reduce this uncertainty, based around individual LSMs that can represent the breadth of assumptions represented in current CMIP-type efforts, while at the same time exploring that uncertainty in a systematic way, confronting multiple types of observations, and where justified, replacing process representations with ML-driven emulators. We propose an approach of modular hybrid modeling to address these challenges.

58 GEOSCIENCES↗
DOE OSTI2024

Estimating Uncertainty in Simulated ENSO Statistics

Abstract Large ensembles of model simulations are frequently used to reduce the impact of internal variability when evaluating climate models and assessing climate change induced trends. However, the optimal number of ensemble members required to distinguish model biases and climate change signals from internal variability varies across models and metrics. Here we analyze the mean, variance and skewness of precipitation and sea surface temperature in the eastern equatorial Pacific region often used to describe the El Niño–Southern Oscillation (ENSO), obtained from large ensembles of Coupled model intercomparison project phase 6 climate simulations. Leveraging established statistical theory, we develop and assess equations to estimate, a priori, the ensemble size or simulation length required to limit sampling‐based uncertainties in ENSO statistics to within a desired tolerance. Our results confirm that the uncertainty of these statistics decreases with the square root of the time series length and/or ensemble size. Moreover, we demonstrate that uncertainties of these statistics are generally comparable when computed using either pre‐industrial control or historical runs. This suggests that pre‐industrial runs can sometimes be used to estimate the expected uncertainty of statistics computed from an existing historical member or ensemble, and the number of simulation years (run duration and/or ensemble size) required to adequately characterize the statistic. This advance allows us to use existing simulations (e.g., control runs that are performed during model development) to design ensembles that can sufficiently limit diagnostic uncertainties arising from simulated internal variability. These results may well be applicable to variables and regions beyond ENSO.

54 ENVIRONMENTAL SCIENCES↗