SEARCH · Search NASA

Results for “SYSTEMS MODELING”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Evaluation of the interactive stratospheric ozone (O3v2) module in the E3SM version 1 Earth system model

Stratospheric ozone affects climate directly as the predominant heat source in the stratosphere and indirectly through chemical reactions controlling other greenhouse gases. The U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) implemented a new ozone chemistry module that improves the simulation of the sharp tropopause gradients, replacing a version based partly on long-term average climatologies that poorly represented heating rates in the lowermost stratosphere. The new O3v2 module extends seamlessly into the troposphere and preserves the naturally sharp cross-tropopause gradient, with 20 %–40 % less ozone in this region. Additionally, O3v2 enables the diagnosis of stratosphere–troposphere exchange flux of ozone, a key budget term lacking in E3SMv1. Here, we evaluate key features in ozone abundance and other closely related quantities in atmosphere-only E3SMv1 simulations driven by observed sea surface temperatures (SSTs, years 1990–2014), comparing them with satellite observations of ozone and also with the University of California, Irvine chemistry transport model (UCI CTM) using the same stratospheric chemistry scheme but driven by European Centre forecast fields for the same period. In terms of stratospheric column ozone, O3v2 shows reduced mean bias and improved northern midlatitude variability, but it is not quite as good as the UCI CTM. As expected, SST-forced E3SMv1 simulations cannot synchronize with observed quasi-biennial oscillations (QBOs), but they do show the typical QBO pattern seen in column ozone. This new O3v2 E3SMv1 model mostly retains the same climate state and climate sensitivity as the previous version, and we recommend its use for other climate models that still use ozone climatologies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Predicting Near-Term Changes in the Earth System: A Large Ensemble of Initialized Decadal Prediction Simulations Using the Community Earth System Model

The objective of near-term climate prediction is to improve our fore-knowledge, from years to a decade or more in advance, of impactful climate changes that can in general be attributed to a combination of internal and externally forced variability. Predictions initialized using observations of past climate states are tested by comparing their ability to reproduce past climate evolution with that of uninitialized simulations in which the same radiative forcings are applied. A new set of decadal prediction (DP) simulations has recently been completed using the Community Earth System Model (CESM) and is now available to the community. This new large-ensemble (LE) set (CESM-DPLE) is composed of historical simulations that are integrated forward for 10 years following initialization on 1 November of each year between 1954 and 2015. CESM-DPLE represents the “initialized” counterpart to the widely studied CESM Large Ensemble (CESM-LE); both simulation sets have 40-member ensembles, and they use identical model code and radiative forcings. Comparing CESM-DPLE to CESM-LE highlights the impacts of initialization on prediction skill and indicates that robust assessment and interpretation of DP skill may require much larger ensembles than current protocols recommend. CESM-DPLE exhibits significant and potentially useful prediction skill for a wide range of fields, regions, and time scales, and it shows widespread improvement over simpler benchmark forecasts as well as over a previous initialized system that was submitted to phase 5 of the Coupled Model Intercomparison Project (CMIP5). Finally, the new DP system offers new capabilities that will be of interest to a broad community pursuing Earth system prediction research.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Variable C : N : P stoichiometry of dissolved organic matter cycling in the Community Earth System Model

Dissolved organic matter (DOM) plays an important role in the ocean's biological carbon pump by providing an advective/mixing pathway for ~ 20% of export production. DOM is known to have a stoichiometry depleted in nitrogen (N) and phosphorus (P) compared to the particulate organic matter pool, a fact that is often omitted from biogeochemical ocean general circulation models. However the variable C : N : P stoichiometry of DOM becomes important when quantifying carbon export from the upper ocean and linking the nutrient cycles of N and P with that of carbon. Here we utilize recent advances in DOM observational data coverage and offline tracer-modeling techniques to objectively constrain the variable production and remineralization rates of the DOM C : N : P pools in a simple biogeochemical-ocean model of DOM cycling. The optimized DOM cycling parameters are then incorporated within the Biogeochemical Elemental Cycling (BEC) component of the Community Earth System Model (CESM) and validated against the compilation of marine DOM observations. The optimized BEC simulation including variable DOM C : N : P cycling was found to better reproduce the observed DOM spatial gradients than simulations that used the canonical Redfield ratio. Global annual average export of dissolved organic C, N, and P below 100 m was found to be 2.28 Pg C yr -1 (143 Tmol C yr -1 , 16.4 Tmol N yr -1 , and 1 Tmol P yr -1 , respectively, with an average export C : N : P stoichiometry of 225 : 19 : 1 for the semilabile (degradable) DOM pool. Dissolved organic carbon (DOC) export contributed ~ 25% of the combined organic C export to depths greater than 100 m.

Letscher, R. T.↗

DOE OSTI2016

Zymomonas mobilis as a model system for production of biofuels and biochemicals

Summary Zymomonas mobilis is a natural ethanologen with many desirable industrial biocatalyst characteristics. In this review, we will discuss work to develop Z. mobilis as a model system for biofuel production from the perspectives of substrate utilization, development for industrial robustness, potential product spectrum, strain evaluation and fermentation strategies. This review also encompasses perspectives related to classical genetic tools and emerging technologies in this context.

09 BIOMASS FUELS↗

DOE OSTI2019

Regional grid refinement in an Earth system model: impacts on the simulated Greenland surface mass balance

In this study, the resolution dependence of the simulated Greenland ice sheet surface mass balance (GrIS SMB) in the variable-resolution Community Earth System Model (VR-CESM) is investigated. Coupled atmosphere–land simulations are performed on two regionally refined grids over Greenland at 0.5° (~55 km) and 0.25° (~28 km), maintaining a quasi-uniform resolution of 1° (~111 km) over the rest of the globe. On the refined grids, the SMB in the accumulation zone is significantly improved compared to airborne radar and in situ observations, with a general wetting (more snowfall) at the margins and a drying (less snowfall) in the interior GrIS. Total GrIS precipitation decreases with resolution, which is in line with best-available regional climate model results. In the ablation zone, CESM starts developing a positive SMB bias with increased resolution in some basins, notably in the east and the north. The mismatch in ablation is linked to changes in cloud cover in VR-CESM, and a reduced effectiveness of the elevation classes subgrid parametrization in CESM. Overall, our pilot study introduces VR-CESM as a new tool in the cryospheric sciences, which could be used to dynamically downscale SMB in scenario simulations and to force dynamical ice sheet models through the CESM coupling framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Convex Relaxation of Grid-Connected Energy Storage System Models With Complementarity Constraints in DC OPF

Including complementarity constraints in energy storage system (ESS) models in optimization problems ensure an optimal solution will not produce a physically unrealizable control strategy where there is simultaneous charging and discharging. However, the current approaches to impose complementarity constraints require the use of non-convex optimization methods. Here, we propose a convex relaxation for a common ESS model that has terms for both charging and discharging based on a penalty reformulation for use in a model predictive control (MPC) based optimal power flow (DC OPF) problem. In this approach, the complementarity constraints are omitted and a penalty term is added to the optimization objective function. For the DC OPF problem, we provide analysis for the conditions under which the convex relaxation of the complementarity constraint ensures that a solution with simultaneous ESS charging and discharging operation is suboptimal. Simulation results demonstrating ESS behavior with and without the penalty reformulation are provided for an MPC-based DC OPF problem on multiple IEEE test systems.

25 ENERGY STORAGE↗

NASA NTRSDate not supplied

Combining Earth System Modeling and Machine Learning to Investigate Volcanic Sulfate Deposition in Polar Ice Cores

Volcanic eruptions emit large amounts of sulfur dioxide (SO2), water, and other chemicals into the atmosphere, both in the troposphere and the stratosphere. Most of the SO2 is converted to sulfate aerosol, which is eventually deposited following long-range transport. The deposits from large eruptions are potentially detectable in ice cores, but there are many cases in which sulfate layers have not been linked to their source volcanoes. As volcanoes can act as significant shocks to the global climate system, we are interested in locating these eruptions in order to increase understanding of the volcanic record. To narrow down the search, we performed 140 simulations of volcanic eruptions using the GISS ModelE Earth system model. We varied the latitude, longitude, Julian day, plume top, plume bottom, and injected SO2 and H2O amounts using a Latin hypercube sampling approach, and analyzed correlations between these parameters and sulfate depositions at ice core sites in Antarctica and Greenland. Using machine learning and parameter estimation, we generated probability distributions and maximum likelihood estimates for the parameters given sulfate deposition data, which can predict latitude with some skill. We find that the volcano latitude and SO2 content are best correlated with sulfate depositions at each pole, while longitude, Julian day, and H2O have small or insignificant effects. Plume altitude and thickness are important because they determine how much of the SO2 is injected into the stratosphere, which has implications for sulfur transport and lifetimes.

Earth system models↗

DOE OSTI2022

A Macroalgal Cultivation Modeling System (MACMODS): Evaluating the Role of Physical-Biological Coupling on Nutrients and Farm Yield

Offshore aquaculture has the potential to expand the macroalgal industry. However, moving into deeper waters requires suspended structures that will present novel farm-environment interactions. Here, we present a computational modeling framework, the Macroalgal Cultivation Modeling System (MACMODS), to explore within-farm modifications to light, seawater flow, and nutrient fields across time and space scales relevant to macroalgae. A regional ocean model informs the site-specific setting, the Santa Barbara Channel in the Southern California Bight. A fine-scale hydrodynamic model predicts modified flows and turbulent mixing within the farm. A spatially resolved macroalgal growth model, parameterized for giant kelp, Macrocystis pyrifera , predicts kelp biomass. Key findings from model integration are that regional ocean conditions set overall farm performance, while fine-scale within-farm circulation and nutrient delivery are important to resolve variation in within-farm macroalgal performance. Therefore, we conclude that models resolving within-farm dynamics can provide benefit to farmers with insight on how farm design and regional ocean conditions interact to influence overall yield. Here, the presence of repeating longlines aligned with the mean current generate flow diversions around the farm as well as attached Langmuir circulations and increased turbulence intensity. These flow-induced phenomena lead to less biomass in the interior portion of the farm relative to the edges. We also find that there is an effluent “footprint” that extends as much as 20 km beyond the farm. In this regard, MACMODS can be used to not only evaluate farm design and cultivation practices that maximize yield but also explore interactions between the farm and ecosystem in order to minimize impacts.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1971

Binder model system to be used for determination of prepolymer functionality

Development of a method for determining the functionality distribution of prepolymers used for rocket binders is discussed. Research has been concerned with accurately determining the gel point of a model polyester system containing a single trifunctional crosslinker, and the application of these methods to more complicated model systems containing a second trifunctional crosslinker, monofunctional ingredients, or a higher functionality crosslinker. Correlations of observed with theoretical gel points for these systems would allow the methods to be applied directly to prepolymers.

Martinelli, F. J.↗

DOE OSTI2023

Climate, Variability, and Climate Sensitivity of “Middle Atmosphere” Chemistry Configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6))

Abstract Simulating whole atmosphere dynamics, chemistry, and physics is computationally expensive. It can require high vertical resolution throughout the middle and upper atmosphere, as well as a comprehensive chemistry and aerosol scheme coupled to radiation physics. An unintentional outcome of the development of one of the most sophisticated and hence computationally expensive model configurations is that it often excludes a broad community of users with limited computational resources. Here, we analyze two configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6)) with simplified “middle atmosphere” chemistry at nominal 1 and 2° horizontal resolutions. Using observations, a reanalysis, and direct model comparisons, we find that these configurations generally reproduce the climate, variability, and climate sensitivity of the 1° nominal horizontal resolution configuration with comprehensive chemistry. While the background stratospheric aerosol optical depth is elevated in the middle atmosphere configurations as compared to the comprehensive chemistry configuration, it is comparable among all configurations during volcanic eruptions. For any purposes other than those needing an accurate representation of tropospheric organic chemistry and secondary organic aerosols, these simplified chemistry configurations deliver reliable simulations of the whole atmosphere that require 35% and 86% fewer computational resources at nominal 1 and 2° horizontal resolution, respectively.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Reducing Long‐Standing Surface Ozone Overestimation in Earth System Modeling by High‐Resolution Simulation and Dry Deposition Improvement

The overestimation of surface ozone concentration in low‐resolution global atmospheric chemistry and climate models has been a long‐standing issue. We first update the ozone dry deposition scheme in both high‐ (0.25°) and low‐resolution (1°) Community Earth System Model (CESM) version 1.3 runs, by adding the effects of leaf area index and correcting the sunlit and shaded fractions of stomatal resistances. With this update, 5‐year‐long summer simulations (2015–2019) using the low‐resolution CESM still exhibit substantial ozone overestimation (by 6.0–16.2 ppbv) over the U.S., Europe, eastern China, and ozone pollution hotspots. The ozone dry deposition scheme is further improved by adjusting the leaf cuticle conductance, reducing the mean ozone bias by 19%, and increasing the model resolution further reduces the ozone overestimation by 43%. We elucidate the mechanism by which model grid spacing influences simulated ozone, revealing distinctive pathways in urban versus rural areas. In rural areas, grid spacing mainly affects daytime ozone levels, where additional NO x emissions from nearby urban areas result in an ozone boost and overestimation in low‐resolution simulations. In contrast, over urban areas, daytime ozone overestimation follows a similar mechanism due to the influence of volatile organic compounds from surrounding rural areas. However, nighttime ozone overestimation is closely linked to weakened NO titration owing to the redistribution of urban NO x to rural areas. Additionally, stratosphere‐troposphere exchange may also contribute to reducing ozone bias in high‐resolution simulations, warranting further investigation. This optimized high‐resolution CESM may enhance understanding of ozone formation mechanisms, sources, and changes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Forced Patterns of Sea Level Rise in the Community Earth System Model Large Ensemble From 1920 to 2100

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid- to late-20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near-surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. Finally, the relative roles of surface flux and ocean convergence anomalies are examined.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Thermodynamic versus Dynamic Controls on Extreme Precipitation in a Warming Climate from the Community Earth System Model Large Ensemble

The moisture budget is evaluated as a function of the probability distribution of precipitation for the end of the twentieth century and projected end of the twenty-first century in the Community Earth System Model Large Ensemble. For a given precipitation percentile, a conditional moisture budget equation relates precipitation minus evaporation ( P − E) to vertical moisture transport, horizontal moisture advection, and moisture storage. At high percentiles, moisture advection and moisture storage cancel and evaporation is negligible, so that precipitation is approximately equal to vertical moisture transport, and likewise for projected changes. Therefore, projected changes to extreme precipitation are approximately equal to the sum of thermodynamic and dynamic tendencies, representing changes to the vertical profiles of moisture content and mass convergence, respectively. The thermodynamic tendency is uniform across percentiles and regions as an intensification of the hydrological cycle, but the dynamic tendency is more complex. For extreme events, per degree of warming, in the mid-to-high latitudes the dynamic tendency is small, so that precipitation approximately scales by the Clausius–Clapeyron 7% K −1 increase. In the subtropics, a drying tendency originating from dynamics offsets the thermodynamic wetting tendency, with the net effect on precipitation varying among regions. The effect of this dynamic drying decreases with increasing percentile. In the deep tropics, a positive dynamic tendency occurs with magnitude similar to or greater than the positive thermodynamic tendency, resulting in generally a 10%–15% K −1 precipitation increase, and with a >25% K −1 increase over the tropical east Pacific. This reinforcing dynamical tendency increases rapidly for high percentiles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Response of Simple, Model Systems to Extreme Conditions

The focus of the research was on the application of high-pressure/high-temperature techniques, together with intense energetic ion beams, to the study of the behavior of simple oxide systems (e.g., SiO 2 , GeO 2 , CeO 2 , TiO 2 , HfO 2 , SnO 2 , ZnO and ZrO 2 ) under extreme conditions. These simple stoichiometries provide unique model systems for the analysis of structural responses to pressure up to and above 1 Mbar, temperatures of up to several thousands of kelvin, and the extreme energy density generated by energetic heavy ions (tens of keV/atom). The investigations included systematic studies of radiation- and pressure-induced amorphization of high P-T polymorphs. By studying the response of simple stoichiometries that have multiple structural “outcomes”, we have established the basic knowledge required for the prediction of the response of more complex structures to extreme conditions. We especially focused on the amorphous state and characterized the different non-crystalline structure-types that result from the interplay of radiation and pressure. For such experiments, we made use of recent technological developments, such as the perforated diamond-anvil cell and in situ investigation using synchrotron x-ray sources. We have been particularly interested in using extreme pressures to alter the electronic structure of a solid prior to irradiation. We expected that the effects of modified band structure would be evident in the track structure and morphology, information which is much needed to describe theoretically the fundamental physics of track-formation. Finally, we investigated the behavior of different simple-oxide, composite nanomaterials (e.g., uncoated nanoparticles vs. core/shell systems) under coupled, extreme conditions. This provided insight into surface and boundary effects on phase stability under extreme conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Photovoltaic System Modeling: Uncertainty and Sensitivity Analyses

We report an uncertainty and sensitivity analysis for modeling AC energy from photovoltaic systems. Output from a PV system is predicted by a sequence of models. We quantify uncertainty in the output of each model using empirical distributions of each model’s residuals. We propagate uncertainty through the sequence of models by sampling these distributions to obtain an empirical distribution of a PV system’s output. We consider models that: (1) translate measured global horizontal, direct and global diffuse irradiance to plane-of-array irradiance; (2) estimate effective irradiance; (3) predict cell temperature; (4) estimate DC voltage, current and power; (5) reduce DC power for losses due to inefficient maximum power point tracking or mismatch among modules; and (6) convert DC to AC power. Our analysis considers a notional PV system comprising an array of FirstSolar FS-387 modules and a 250 kW AC inverter; we use measured irradiance and weather at Albuquerque, NM. We found the uncertainty in PV system output to be relatively small, on the order of 1% for daily energy. We found that uncertainty in the models for POA irradiance and effective irradiance to be the dominant contributors to uncertainty in predicted daily energy. Our analysis indicates that efforts to reduce the uncertainty in PV system output predictions may yield the greatest improvements by focusing on the POA and effective irradiance models.

14 SOLAR ENERGY↗

DOE OSTI2018

Workshop on the Initialization of High-Resolution Earth System Models

From weather to climate timescales, weeks to decades, societal needs require detailed predictions of the atmosphere, ocean, and land-surface to facilitate informed planning. Central to this scientific endeavor is the ability to accurately initialize the earth system models (ESMs), since frequently the skill in forecasts, predictions, and projections depends on the quality of the initialization. However, initialization of high-resolution ESMs remains an outstanding challenge, particularly for the climate modeling community where initialization has received relatively less attention and where longer-term information is needed. This need, after a review of previous workshop outcomes, led the U.S. Department of Energy (DOE) and National Oceanic and Atmospheric Administration (NOAA) to host a workshop that focused exclusively on the opportunities that data assimilation (DA) and, specifically, coupled data assimilation (CDA) might provide to the ESM enterprise with particular attention toward initialization. The workshop was held April 9-10, 2018, at the Rockville Hilton in Rockville, Maryland, with 41 participants. It was composed of seven sessions, each with a distinct topical theme that related to the broad challenge of initializing high-resolution ESMs. The first session was used to frame the workshop challenge. Two of the sessions (2 and 7) focused the current state of the art of CDA within laboratory and operational centers at both the short-term, subseasonal and long-term, climate time scales. Three of the sessions (3, 4, and 5) addressed new algorithmic approaches to CDA at high resolution. Session 6 focused on the computational science solutions for the management of large input/output streams required to carry out the assimilation procedure.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

STITCHES: a Python package to amalgamate existing Earth system model output into new scenario realizations

Understanding the interaction between humans and the Earth system is a computationally daunting task, with many possible approaches depending on resources available and questions of interest. For example, state-of-the-art impact models require decade-long time series of relatively high frequency, spatially resolved and often multiple variables representing climatic impact-drivers (Ruane et al., 2022). Most commonly these are derived from the outputs of detailed, computationally expensive Earth System Models (ESMs) run according to a standard, limited set of future scenarios, the latest being the SSP-RCPs run under CMIP6/ScenarioMIP (Eyring et al., 2016; O’Neill et al., 2016). At the time of writing, O’Neill et al. (2016) has been cited more than 1750 times and Eyring et al. (2016) more than 5000 times, highlighting the broad, general applications of this data. Often, however, impact modeling seeks to explore new scenarios that were not part of the ScenarioMIP protocol, and/or needs a larger set of initial condition ensemble members than are typically available to quantify the effects of ESM internal variability. In addition, the recognition that the human and Earth systems are fundamentally intertwined, and may feature potentially significant feedback loops, is making integrated, simultaneous modeling of the coupled human-Earth system increasingly necessary, if computationally challenging with most existing tools (Thornton et al., 2017). For impact modelers, climate model emulators can be the answer to meet both the needs of: 1) creating realizations for novel scenarios and 2) achieving a simplified, computationally tractable representation of ESM behavior in a coupled human-Earth system modeling framework. We proposed a new, comprehensive approach to such emulation of gridded, multivariate ESM outputs for novel scenarios without the computational cost of a full ESM, STITCHES (Tebaldi et al., 2022). The approach outlined in Tebaldi et al. (2022) should be extensible to future CMIP eras, although the STITCHES software at present is strictly focused on CMIP6/ScenarioMIP data hosted on Pangeo (https://gallery.pangeo.io/repos/pangeo-gallery/cmip6/). The corresponding STITCHES Python package uses existing archives of ESMs’ scenario experiments from CMIP6/ScenarioMIP to construct gridded, multivariate realizations of new scenarios provided by reduced complexity climate models (Hartin et al., 2015; Meinshausen et al., 2011; Smith et al., 2018), or to enrich existing initial condition ensembles. Its output provides the same characteristics as the emulated ESM output: multivariate (spanning potentially all variables that the ESM has saved), spatially resolved (down to the native grid of the ESM), and preserving the same high frequency as the original data. A new realization of multiple variables can be generated on the order of minutes with STITCHES, rather than the hours or sometimes days that ESMs require.

97 MATHEMATICS AND COMPUTING↗