Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth system 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.

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Estimating the Likelihood of GHG Concentration Scenarios From Probabilistic Integrated Assessment Model Simulations

The climate scenarios that form the basis for current climate risk assessments have no assigned probabilities, and this impedes the analysis of future climate risks. This paper proposes an approach to estimate the probability of carbon dioxide (CO 2 ) concentration scenarios used in key climate change modeling experiments. It computes the CO 2 emissions compatible with the concentrations prescribed by Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 experiments. The distribution of these compatible cumulative emissions is interpreted as the likelihood of future emissions given a concentration pathway. Using Bayesian analysis, the probability of each pathway can be estimated from a probabilistic sample of future emissions. The approach is demonstrated with five probabilistic CO 2 emission simulation ensembles from four Integrated Assessment Models (IAM), leading to independent estimates of the likelihood of the CO 2 concentration of Representative Concentration Pathways (RCP) and Shared Socioeconomic Pathways (SSP). Results suggest that SSP5-8.5 is unlikely for the second half of the 21st century, but offer no clear consensus on which of the remaining scenarios is most likely. Estimates of likelihoods of CO 2 concentrations associated with RCP and SSP scenarios are affected by sampling errors, differences in emission sources simulated by the IAMs, and a lack of a common experimental framework for IAM simulations. These shortcomings, along with a small IAM ensemble size, limit the applicability of the results presented here. Novel joint IAM and the Earth System Model experiments are needed to deliver actionable probabilistic climate risk assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global gridded crop harvested area, production, yield, and monthly physical area data circa 2015

Abstract Here we provide an update to global gridded annual and monthly crop datasets. This new dataset uses the crop categories established by the Global Agro-Ecological Zones (GAEZ) Version 3 model, which is based on the Food and Agricultural Organization of the United Nations (FAO) crop production data. We used publicly available data from the FAOSTAT database as well as GAEZ Version 4 global gridded dataset to generate circa 2015 annual crop harvested area, production, and yields by crop production system (irrigated and rainfed) for 26 crops and crop categories globally at 5-minute resolution. We additionally used available data on crop rotations, cropping intensity, and planting and harvest dates to generate monthly gridded cropland data for physical areas for the 26 crops by production system. These data are in standard georeferenced gridded format, and can be used by any global hydrology, land surface, or other earth system model that requires gridded annual or monthly crop data inputs.

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 OSTI2023

East African Monsoon as a Drawbridge for the Circumnavigation of Madden–Julian Oscillation Events

The seasonality of circumnavigation of Madden-Julian Oscillation (MJO) events is examined using 42 years of observations and global model simulations. MJO events are most likely to propagate across the Atlantic Ocean, Africa and back to equatorial Indian Ocean relatively unweakened during March-April-May and to a lesser extent in October-November-December. These seasons correspond to the "long-rain" and "short-rain" periods of East African monsoon when precipitation over Africa is near the equator and extend farthest east toward Indian Ocean. Moisture budget analysis of a circumnavigating MJO event under real and hypothetical seasonally varying background conditions indicates that the seasonal zonal movement of zonal component of background circulation regulates the moistening of MJO crossing from Eastern Africa to Indian Ocean. Ensembles of Energy Exa-scale Earth System Model (E3SM) simulations of the MJO event under idealized seasonally varying background state confirm the observed drawbridge-like role of East African monsoon in the propagation of MJO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Subseasonal Tropical Cyclone Prediction and Modulations by MJO and ENSO in CESM2

Subseasonal tropical cyclone (TC) reforecasts from the Community Earth System Model version 2 (CAM6) subseasonal prediction system are examined in this study. Here, we evaluate the modeled TC climatology and the probabilistic forecast skill of basin-wide TC genesis at weekly temporal resolution. Prediction skill is calculated using the Brier skill score relative to a constant annual mean climatology and to a monthly varying seasonal climatology during TC season. The model captures the observed basin-wide climatological TC seasonality and spatial distributions at weeks 1–6, but TC genesis is largely underestimated from Week 2 onward. For some basins and lead times, the predicted TC genesis is primarily controlled by the number of TC “seeds” and the mean-state climate condition. The model has good prediction skill relative to the constant climatology across all the basins and lead times, but is only skillful in the eastern Pacific, North Indian Ocean, and Southern Hemisphere at Week 1 when compared to the seasonal climatology, indicating limited skill in predicting deviations from the seasonal cycle. We find strong modulations of the predicted TC genesis at up to 3 weeks of forecast lead time by the Madden-Julian Oscillation. The interannual variability of predicted TC genesis and accumulated cyclone energy are skillfully predicted in the North Atlantic and the Northwestern Pacific, with a strong modulation by the El Nino-Southern Oscillation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geospatial Diffusion for Land Cover Imperviousness Change Forecasting

Land-use and land-cover (LULC) has a significant effect on several Earth system processes. For example, impervious surfaces reduce infiltration and speed water flow, impacting regional hydrology and flood risk. While Earth System models have improved forecasting hydrologic and atmospheric processes at higher resolutions, the ability to forecast LULC change has lagged behind. In this paper, we propose a new paradigm exploiting Generative AI (GenAI) for land cover change forecasting by framing it as a data synthesis problem conditioned on historical and auxiliary data-sources. To demonstrate the feasibility of our methodology, we perform experiments where a diffusion model is trained for decadal forecasting of imperviousness change across the entire United States. We find that our model yields MAE lower than a no-change baseline for resolutions ≥ 0.7 X 0.7km2 on average, demonstrating its ability to capture and project accurate spatiotemporal patterns. Finally, we discuss future research to incorporate Earth's physical properties and enabling scenario simulations via driver variables.

Varshney, Debvrat [ORNL] (ORCID:0000000188981736)↗
DOE OSTI2020

Global biogeography of fungal and bacterial biomass carbon in topsoil

Bacteria and fungi, representing two major soil microorganism groups, play an important role in global nutrient biogeochemistry. Biogeographic patterns of bacterial and fungal biomass are of fundamental importance for mechanistically understanding nutrient cycling. We synthesized 1323 data points of phospholipid fatty acid-derived fungal biomass C (FBC), bacterial biomass C (BBC), and fungi:bacteria (F:B) ratio in topsoil, spanning 11 major biomes. The FBC, BBC, and F:B ratio display clear biogeographic patterns along latitude and environmental gradients including mean annual temperature, mean annual precipitation, net primary productivity, root C density, soil temperature, soil moisture, and edaphic factors. At the biome level, tundra has the highest FBC and BBC densities at 3684 (95% confidence interval: 1678–8084) mg kg -1 and 428 (237–774) mg kg -1 , respectively; desert has the lowest FBC and BBC densities at 16.92 (14.4–19.89) mg kg -1 and 6.83 (6.1–7.65) mg kg -1 , respectively. The F:B ratio varies dramatically, ranging from 1.8 (1.6–2.1) in savanna to 8.6 (6.7–11.0) in tundra. An empirical model was developed for the F:B ratio and it is combined with a global dataset of soil microbial biomass C to produce global maps for FBC and BBC in 0–30 cm topsoil. Across the globe, the highest FBC is found in boreal forest and tundra while the highest BBC is in boreal forest and tropical/subtropical forest, the lowest FBC and BBC are in shrub and desert. Global stocks of living microbial biomass C were estimated to be 12.6 (6.6–16.4) Pg C for FBC and 4.3 (0.5–10.3) Pg C for BBC in topsoil. These findings advance our understanding of the global distribution of fungal and bacterial biomass, which facilitates the incorporation of fungi and bacteria into Earth system models. The global maps of bacterial and fungal biomass serve as a benchmark for validating microbial models in simulating the global C cycle under a changing climate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A unified ensemble soil moisture dataset across the continental United States

Abstract A unified ensemble soil moisture (SM) package has been developed over the Continental United States (CONUS). The data package includes 19 products from land surface models, remote sensing, reanalysis, and machine learning models. All datasets are unified to a 0.25-degree and monthly spatiotemporal resolution, providing a comprehensive view of surface SM dynamics. The statistical analysis of the datasets leverages the Koppen-Geiger Climate Classification to explore surface SM’s spatiotemporal variabilities. The extracted SM characteristics highlight distinct patterns, with the western CONUS showing larger coefficient of variation values and the eastern CONUS exhibiting higher SM values. Remote sensing datasets tend to be drier, while reanalysis products present wetter conditions. In-situ SM observations serve as the basis for wavelet power spectrum analyses to explain discrepancies in temporal scales across datasets facilitating daily SM records. This study provides a comprehensive soil moisture data package and an analysis framework that can be used for Earth system model evaluations and uncertainty quantification, quantifying drought impacts and land–atmosphere interactions and making recommendations for drought response planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large emissions of CO 2 and CH 4 due to active-layer warming in Arctic tundra

Climate warming may accelerate decomposition of Arctic soil carbon, but few controlled experiments have manipulated the entire active layer. To determine surface-atmosphere fluxes of carbon dioxide and methane under anticipated end-of-century warming, here we used heating rods to warm (by 3.8 °C) to the depth of permafrost in polygonal tundra in $Utqia\dot{g}vik$ (formerly Barrow), Alaska and measured fluxes over two growing seasons. We show that ecosystem respiration is ~30% higher in warmed plots than in control plots (0.99 μmol m -2 s -1 versus 0.67 μmol m -2 s -1 , p < 0.0001, n = 79). Additionally, the observed temperature sensitivity (Q 10 of 2.8) is higher than that imposed for soil in Earth system models or reported by arctic experiments warming only the surface. A shoulder-season warming experiment revealed that rapid snow melt, which is becoming a more common event, can result in large methane emissions that may have otherwise been oxidized to carbon dioxide. Thus, warming promotes greenhouse gas emissions from the whole, deepening active layer and may contribute to climate change amplification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cooling outweighs warming across phenological transitions in the Northern Hemisphere

Vegetation phenology, i.e., seasonal biological events such as leaf-out and leaf-fall, regulates local climate through biophysical processes like evapotranspiration (ET) and albedo. However, the net surface temperature impact of these processes—whether ET cooling or albedo-induced warming predominates—and how the dominance changes across phenological transitions and regions remains poorly understood. Here, we investigated the effects of vegetation foliage on daytime land surface temperature (LST) following six phenological transitions, spanning from the start of season to end of season, in deciduous and mixed forests across the mid- to high-latitude Northern Hemisphere during 2013–2021 using multiple satellite products and ground observations. We quantified vegetation effect as the difference between observed LST and LST estimates from the Annual Temperature Cycle (ATC) model, representing a no-foliage scenario. We found that vegetation-induced cooling consistently outweighs warming following all phenological transitions except for the end of the season. Cooling intensity increased with vegetation greenness, ranging from 1.0 ± 0.5 °C (mean ± 0.15 SD) in 59% of forests after the start of the season (SOS) to 6.1 ± 0.8 °C in 89% of forests following the onset of maturity, before declining toward the end of the season. Over half of the regions experiencing cooling showed intensification of surface cooling with climate warming, suggesting an amplified vegetation-mediated cooling under future climate change. The findings provide a more precise understanding of the role of vegetation in modulating climate at the intraseasonal scale, highlighting the importance of integrating phenological impacts into climate adaptation strategies and Earth system modeling.

albedo↗
DOE OSTI2022

What explains the year-to-year variation in growing season timing of boreal black spruce forests?

Amplified climate warming in high latitudes is expected to affect growing season timing of the vast boreal biome. It is unclear whether the presence of permafrost (perennially frozen ground) might have an influence on changes in growing season timing. Here, this study examined how different environmental variables explained, either directly or indirectly, the variation in growing season timing of boreal forest stands with and without permafrost. We expected that environmental variables explaining the variation in growing season timing differed or had different explanatory power depending on permafrost presence or absence. The growing season was delineated from daily gross primary productivity (GPP) time series derived from 40 site-year data of net ecosystem carbon dioxide exchange measured with eddy covariance techniques over five black spruce (Picea mariana [Mill.])-dominated boreal forest stands in North America. In permafrost-free forest stands, a combination of start in canopy ‘green-up’ in spring and the timing of air and soil temperature increasing above freezing explained the start-of-season (SOS GPP ). Results from commonality analysis and structural equation modeling suggest that canopy ‘green-up’ and air temperature directly affected SOS GPP in permafrost-free forest stands. In addition, soil temperature acted as mediator for an indirect effect of air temperature on SOS GPP . In contrast, none of the environmental variables, or their combination, explained the variation in SOS GPP in forest stands with permafrost. The explanatory power of environmental variables was more consistent regarding the end-of-season (EOS GPP ). In both, forest stands with and without permafrost, EOS GPP was directly explained by mean soil water content in the fall and the first day of continuous snowpack formation. A better understanding how environmental variables control SOS GPP and EOS GPP in forest stands with and without permafrost will help to refine parameterizations of the boreal biome in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dampening of the Precipitation Response to Aerosol Pollution From Turbulence in Cumulus Clouds

Using aircraft observations and a cloud model with particle‐based microphysics, it is shown that the enhancement of drop collision‐coalescence from turbulence in clouds not only leads to earlier onset of rain in warm cumulus clouds as past studies have suggested, but also significantly dampens the precipitation susceptibility to aerosol loading. Enhanced drop coalescence from turbulence substantially increases the production of drizzle embryos just above cloud base, which in turn act as seeds that accelerate rain drop growth at mid and upper cloud levels even in highly polluted conditions. In contrast, pollution aerosols strongly inhibit rainfall when the commonly assumed gravitational‐only collision kernel is used and turbulent coalescence is neglected. Overall, turbulence‐enhanced drop coalescence strongly influences the response of warm cumulus clouds and precipitation to aerosol loading, suggesting that the effects of turbulent coalescence should be included in Earth system model representations of aerosol‐cloud‐precipitation interactions and aerosol indirect radiative forcing.

Chandrakar, Kamal Kant [NSF National Center for At↗
DOE OSTI2024

Surrogate construction via weight parameterization of residual neural networks

Surrogate model development is a critical step for uncertainty quantification or other sample-intensive tasks for complex computational models. Here, in this work, we develop a multi-output surrogate form using a class of neural networks (NNs) that employ shortcut connections, namely Residual NNs (ResNets). ResNets are known to regularize the surrogate learning problem and improve the efficiency and accuracy of the resulting surrogate. Inspired by the continuous, Neural ODE analogy, we augment ResNets with weight parameterization strategy with respect to ResNet depth. Weight-parameterized ResNets regularize the NN surrogate learning problem and allow better generalization with a drastically reduced number of learnable parameters. We demonstrate that weight-parameterized ResNets are more accurate and efficient than conventional feed-forward multi-layer perceptron networks. We also compare various options for parameterization of the weights as functions of ResNet depth. We demonstrate the results on both synthetic examples and a large scale earth system model of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A mountain-induced moist baroclinic wave test case for the dynamical cores of atmospheric general circulation models

Idealized test cases for the dynamical cores of atmospheric general circulation models are informative tools to assess the accuracy of the numerical designs and investigate the general characteristics of atmospheric motions. A new test case is introduced that is built upon a baroclinically unstable base state with an added orographic barrier. The topography is analytically prescribed and acts as a trigger of both baroclinic Rossby waves and inertia–gravity waves on a rotating, regular-sized planet. Both dry and idealized moist configurations are suggested. The latter utilizes the Kessler warm-rain precipitation scheme. The test case enhances the complexity of the existing test suite hierarchy and focuses on the impacts of two midlatitudinal mountain ridges on the circulation. Selected simulation examples from four dynamical cores are shown. These are the Spectral Element and Finite Volume dynamical cores, which are part of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM), versions 2.1.3 and 2.2, and the Cubed-Sphere Finite Volume dynamical cores, which is new to CESM version 2.2. In addition, the Model for Prediction Across Scales (MPAS) is tested. The overall flow patterns agree well in the four dynamical cores, but the details can vary greatly. The examples highlight the broad palette of use cases for the test case and reveal physics–dynamics coupling issues.

58 GEOSCIENCES↗
DOE OSTI2020

Expanding the design space of stratospheric aerosol geoengineering to include precipitation-based objectives and explore trade-offs

Previous climate modeling studies demonstrate the ability of feedback-regulated, stratospheric aerosol geoengineering with injection at multiple independent latitudes to meet multiple simultaneous temperature-based objectives in the presence of anthropogenic climate change. However, the impacts of climate change are not limited to rising temperatures but also include changes in precipitation, loss of sea ice, and many more; knowing how a given geoengineering strategy will affect each of these climate metrics is vital to understanding the limits and trade-offs of geoengineering. Here, we first introduce a new method of visualizing the design space in which desired climate outcomes are represented by 2-D surfaces on a 3-D graph. Surface orientations represent how different injection choices influence that objective, and intersecting surfaces represent objectives which can be met simultaneously. Using this representation as a guide, we present simulations of two new strategies for feedback-regulated aerosol injection, using the Community Earth System Model with the Whole Atmosphere Community Climate Model – CESM1(WACCM). The first simultaneously manages global mean temperature, tropical precipitation centroid, and Arctic sea ice extent, while the second manages global mean precipitation, tropical precipitation centroid, and Arctic sea ice extent. Both simulations control the tropical precipitation centroid to within 5 % of the goal, and the latter controls global mean precipitation to within 1% of the goal. Additionally, the first simulation overcompensates sea ice, while the second undercompensates sea ice; all of these results are consistent with the expectations of our design space model. In addition to showing that precipitation-based climate metrics can be managed using feedback alongside other goals, our simulations validate the utility of our design space visualization in predicting our climate model behavior under a given geoengineering strategy, and together they help illustrate the fundamental limits and trade-offs of stratospheric aerosol geoengineering.

58 GEOSCIENCES↗