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,099 records · Page 61

Hemispherically symmetric strategies for stratospheric aerosol injection

Stratospheric aerosol injection (SAI) comes with a wide range of possible design choices, such as the location and timing of the injection. Different stratospheric aerosol injection strategies can yield different climate responses; therefore, understanding the range of possible climate outcomes is crucial to making informed future decisions on SAI, along with the consideration of other factors. Yet, to date, there has been no systematic exploration of a broad range of SAI strategies. This limits the ability to determine which effects are robust across different strategies and which depend on specific injection choices. This study systematically explores how the choice of SAI strategy affects climate responses in one climate model. Here, we introduce four hemispherically symmetric injection strategies, all of which are designed to maintain the same global mean surface temperature: an annual injection at the Equator (EQ), an annual injection of equal amounts of SO 2 at 15° N and 15° S (15N+15S), an annual injection of equal amounts of SO 2 at 30° N and 30° S (30N+30S), and a polar injection strategy that injects equal amounts of SO 2 at 60° N and 60° S only during spring in each hemisphere (60N+60S). We compare these four hemispherically symmetric SAI strategies with a more complex injection strategy that injects different quantities of SO 2 at 30° N, 15° N, 15° S, and 30° S in order to maintain not only the global mean surface temperature but also its large-scale horizontal gradients. All five strategies are simulated using version 2 of the Community Earth System Model with the middle atmosphere version of the Whole Atmosphere Community Climate model, version 6, as the atmospheric component, CESM2(WACCM6-MA), with the global warming scenario, Shared Socioeconomic Pathway (SSP)2-4.5. We find that the choice of SAI strategy affects the spatial distribution of aerosol optical depths, injection efficiency, and various surface climate responses. In addition, injecting in the subtropics produces more global cooling per unit injection, with the EQ and the 60N+60S cases requiring, respectively, 59 % and 50 % more injection than the 30N+30S case to meet the same global mean temperature target. Injecting at higher latitudes results in larger Equator-to-pole temperature gradients. While all five strategies restore Arctic September sea ice, the high-latitude injection strategy is more effective due to the SAI-induced cooling occurring preferentially at higher latitudes. These results suggest trade-offs wherein different strategies appear better or worse, depending on which metrics are deemed important.

54 ENVIRONMENTAL SCIENCES↗

Calcium is associated with specific soil organic carbon decomposition products at Blodgett Forest Research Center, Georgetown, California as analysed with scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy

This data is from the paper calcium is associated with specific soil organic carbon decomposition products, published in SOIL. DOI: https://doi.org/10.5194/soil-11-381-2025, 2025.This file contains CSVs with spectral data and bulk soil data and there is no specific program required to open this data. The data includes Scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy. data from the measurement of samples from the Whole-soil Warming project, run by the Belowground Biogeochemistry team at Blodgett Forest Research Center, Georgetown, California run by the University of California, Berkeley. It also includes bulk soil chemical properties. The University of California's Blodgett Forest Research Station (Forest) is situated in the Sierra Nevada foothills (1370 m a.s.l.) near Georgetown, California. The samples were collected from here: 38.912013, -120.661469, https://maps.app.goo.gl/291bCJ1zVqUhgktz6. The Forest soils were characterised as Alfisols, which are equivalent to Dystric Cambisols (IUSS Working Group WRB, 2015), and formed in granitic parent materials, in a temperate climate, under thinned, mixed-coniferous forest (Fig. S3; Gaudinski et al., 2009). With these analyses we aimed to answer the question, is calcium associated with a specific type of organic matter enriched in aromatic and phenolic carbon at the microscale in samples from Blodgett Forest Research Center? and how does this specific type of carbon respond to experiments targetted at removing and adding calcium to the soils, specifically cation exchange and incubation after calcium addition? Abstract from the paper can be found below: Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil-types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca-SOC is removed through cation-exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, highlighting a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca, and in turn, this association inhibits dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots on the microscale shifting local decomposition processes and thereby explaining the colocation of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth System Models could improve our understanding, predictions, and management of carbon dynamics in soils, and account for their response to Ca-rich amendments.

54 ENVIRONMENTAL SCIENCES↗

Pervasive shifts in forest dynamics in a changing world

Forest dynamics arise from the interplay of environmental drivers and disturbances with the demographic processes of recruitment, growth, and mortality, subsequently driving biomass and species composition. However, forest disturbances and subsequent recovery are shifting with global changes in climate and land use, altering these dynamics. Changes in environmental drivers, land use, and disturbance regimes are forcing forests toward younger, shorter stands. Rising carbon dioxide, acclimation, adaptation, and migration can influence these impacts. Recent developments in Earth system models support increasingly realistic simulations of vegetation dynamics. In parallel, emerging remote sensing datasets promise qualitatively new and more abundant data on the underlying processes and consequences for vegetation structure. In combination, these advances hold promise for improving the scientific understanding of changes in vegetation demographics and disturbances.

54 ENVIRONMENTAL SCIENCES↗

Phosphorus limitation of early growth differs between nitrogen‐fixing and nonfixing dry tropical forest tree species

Summary Tropical forests are often characterized by low soil phosphorus (P) availability, suggesting that P limits plant performance. However, how seedlings from different functional types respond to soil P availability is poorly known but important for understanding and modeling forest dynamics under changing environmental conditions. We grew four nitrogen (N)‐fixing Fabaceae and seven diverse non‐N‐fixing tropical dry forest tree species in a shade house under three P fertilization treatments and evaluated carbon (C) allocation responses, P demand, P‐use, investment in P acquisition traits, and correlations among P acquisition traits. Nitrogen fixers grew larger with increasing P addition in contrast to non‐N fixers, which showed fewer responses in C allocation and P use. Foliar P increased with P addition for both functional types, while P acquisition strategies did not vary among treatments but differed between functional types, with N fixers showing higher root phosphatase activity (RPA) than nonfixers. Growth responses suggest that N fixers are limited by P, but nonfixers may be limited by other resources. However, regardless of limitation, P acquisition traits such as mycorrhizal colonization and RPA were nonplastic across a steep P gradient. Differential limitation among plant functional types has implications for forest succession and earth system models.

59 BASIC BIOLOGICAL SCIENCES↗

Can Biogeochemical Tracer Observations Constrain Southern Ocean Diapycnal Mixing Rates?

Abstract Direct observations of background diapycnal mixing rates in the Southern Ocean (SO) are limited spatially and temporally, making the choice of an appropriate value to parameterize this mixing in Earth system models a challenge. However, the deployment of Argo floats throughout the SO has provided an extensive range of observations of both physical and biogeochemical parameters. We use an ocean state estimate run with various background diapycnal mixing coefficients to assess if biogeochemical tracer observations can be used to better constrain SO diapycnal mixing rates. We find that vertical tracer distributions in the SO are highly sensitive to the rate of background diapycnal mixing and can provide an upper limit on background mixing rates. This demonstrates the importance of biogeochemical tracer observations throughout the full depth of the water column to validate ocean models.

Ellison, Elizabeth [The Australian National Univer↗

Advancing Understanding of Deep Convective Anvil Clouds

In this project we use ground-based DOE ARM observations, satellite observations, and Weather Research and Forecasting (WRF) model simulations to arrive at a more complete understanding (or, new theory) of how atmospheric heating profiles in thunderstorm cores drives how quickly air rises in the storms. The rate at which air motion increases or decreases with altitude directly relates to how much cloud mass is ejected laterally into the adjacent atmosphere, which then influences the development and maintenance of “stratiform raining cloud anvils”. Our work has led to new equations connecting the heating in these thunderstorms to these stratiform clouds, with the connection modulated by the stability structure of the atmosphere (i.e., if more stable, then a more rapid change in thunderstorm heating with altitude results in more stratiform cloud, independent of whether over land or ocean, or how organized the thunderstorms are). Improving knowledge of stratiform clouds is important, since they contribute to substantial accumulated rainfall (important to agriculture and local hydrology) as well as radiation, given their large areal extents and lifetime. The radiation impacts likely impact simulated cloud feedback in all Earth System models (ESMs).

54 ENVIRONMENTAL SCIENCES↗

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗

Integrating Maximum Entropy Production Theory and Machine Learning to Improve Global Evapotranspiration Modeling

Accurate estimation of terrestrial evapotranspiration (ET) is vital for understanding global water and energy cycles. However, current global ET estimations are not well constrained. This study introduces an integrated framework combining the Maximum Entropy Production (MEP) theory with Random Forest (RF) model to improve global ET estimation. Specifically, in contrast to direct ET estimation by the RF model, the integrated framework (MEP‐RF) trains to predict error of MEP‐simulated ET. MEP‐RF outperforms RF in spatiotemporal extrapolation. Attribution analysis with in situ observations reveals that the inputs of MEP are the most critical variables for the ET process, including net radiation, vegetated area, soil moisture, and surface temperature. We further drive MEP‐RF with global reanalysis and satellite data sets of these four inputs, yielding a global mean terrestrial ET of 548 mm/year, with 77% attributed to transpiration. The global ET increased at a rate of 0.85 mm/year per year during 2003–2021, primarily due to vegetation greening rather than rising temperature, while decreasing soil moisture led to decreasing regional ET. The integrated framework provides a novel approach for the estimation of global ET without the need for hard‐to‐obtain and thus uncertain inputs, such as wind speed, surface roughness, aerodynamic and canopy stomatal resistance. Therefore, MEP‐RF offers an independent method on existing global ET products. It represents a promising physically based approach that can be incorporated into Earth System Models to enhance water and energy cycle simulations.

54 ENVIRONMENTAL SCIENCES↗

Upscaling Soil Organic Carbon Measurements at the Continental Scale Using Multivariate Clustering Analysis and Machine Learning

Abstract Estimates of soil organic carbon (SOC) stocks are essential for many environmental applications. However, significant inconsistencies exist in SOC stock estimates for the U.S. across current SOC maps. We propose a framework that combines unsupervised multivariate geographic clustering (MGC) and supervised Random Forests regression, improving SOC maps by capturing heterogeneous relationships with SOC drivers. We first used MGC to divide the U.S. into 20 SOC regions based on the similarity of covariates (soil biogeochemical, bioclimatic, biological, and physiographic variables). Subsequently, separate Random Forests models were trained for each SOC region, utilizing environmental covariates and SOC observations. Our estimated SOC stocks for the U.S. (52.6 ± 3.2 Pg for 0–30 cm and 108.3 ± 8.2 Pg for 0–100 cm depth) were within the range estimated by existing products like Harmonized World Soil Database, HWSD (46.7 Pg for 0–30 cm and 90.7 Pg for 0–100 cm depth) and SoilGrids 2.0 (45.7 Pg for 0–30 cm and 133.0 Pg for 0–100 cm depth). However, independent validation with soil profile data from the National Ecological Observatory Network showed that our approach ( R 2 = 0.51) outperformed the estimates obtained from Harmonized World Soil Database ( R 2 = 0.23) and SoilGrids 2.0 ( R 2 = 0.39) for the topsoil (0–30 cm). Uncertainty analysis (e.g., low representativeness and high coefficients of variation) identified regions requiring more measurements, such as Alaska and the deserts of the U.S. Southwest. Our approach effectively captures the heterogeneous relationships between widely available predictors and the current SOC baseline across regions, offering reliable SOC estimates at 1 km resolution for benchmarking Earth system models.

58 GEOSCIENCES↗

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

Evapotranspiration depletes groundwater under warming over the contiguous United States

A warmer climate increases evaporative demand. However, response to warming depends on water availability. Existing earth system models represent soil moisture but simplify groundwater connections, a primary control on soil moisture. Here we apply an integrated surface-groundwater hydrologic model to evaluate the sensitivity of shallow groundwater to warming across the majority of the US. We show that as warming shifts the balance between water supply and demand, shallow groundwater storage can buffer plant water stress; but only where shallow groundwater connections are present, and not indefinitely. As warming persists, storage can be depleted and connections lost. Similarly, in the arid western US warming does not result in significant groundwater changes because this area is already largely water limited. The direct response of shallow groundwater storage to warming demonstrates the strong and early effect that low to moderate warming may have on groundwater storage and evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗

Labrador Sea freshening linked to Beaufort Gyre freshwater release

Abstract The Beaufort Gyre (BG), the largest Arctic Ocean freshwater reservoir, has drastically increased its liquid freshwater content by 40% in the past two decades. If released within a short period, the excess freshwater could potentially impact the large-scale ocean circulation by freshening the upper subpolar North Atlantic. Here, we track BG-sourced freshwater using passive tracers in a global ocean sea-ice model and show that this freshwater exited the Arctic mostly through the Canadian Arctic Archipelago, rather than Fram Strait, during an historical release event in 1983–1995. The Labrador Sea is the most affected region in the subpolar North Atlantic, with a freshening of 0.2 psu on the western shelves and 0.4 psu in the Labrador Current. Given that the present BG freshwater content anomaly is twice the historical analog studied here, the impact of a future rapid release on Labrador Sea salinity could be significant, easily exceeding similar fluxes from Greenland meltwater.

54 ENVIRONMENTAL SCIENCES↗

High-resolution leaf area index maps generated from unoccupied aerial system, Teller Mile 27, Seward Peninsula, Alaska

Leaf area index (LAI), a measure of the amount of one-side leaf area per ground unit, is an important indicator of plant carbon, energy, and water cycle. In the heterogeneous Arctic landscapes, it has been challenging to accurately measure LAI across species and space needed for Earth system model validation. Here, we use multispectral unoccupied aerial systems (UASs) to scale up and map leaf area index (LAI) , in a low-Arctic tundra landscape on the Seward Peninsula, Alaska. We linked previous published LAI measurements with high-resolution, UAS-collected multispectral data collected over the region of Next Generation Ecosystem Experiments in the Arctic (NGEE Arctic)’s Teller Mile Maker 27 site in 2022 to develop random forest (RF) machine learning models to predict and map LAI. 100 RF models were developed to account for uncertainties in ground LAI plot measurements and process scaling. This dataset includes a raster (*.tif) map of the mean LAI value of the 100 RF models, a raster (*.tif) map of the standard deviation of the RF-modeled LAI data, and a user guide (*.pdf).

54 ENVIRONMENTAL SCIENCES↗

Irrigation-induced land water depletion aggravated by climate change

Agricultural irrigation has experienced rapid expansion, and its growing freshwater consumption is potentially exacerbating water scarcity issues. Previous studies predominantly relied on observations or land-only simulations, often neglecting land–atmosphere interactions or failing to capture long-term evolution. We therefore analyse the effects of historical irrigation expansion on water fluxes and resources using seven Earth system models. Here we show that irrigation expansion in many regions substantially decreases the net water influx from the atmosphere to land, further aggravating the existing drying trends caused by climate change. For example, irrigation expansion changed the trend of this net influx from −0.664 ( ± 0.283) to −1.461 ( ± 0.261) mm yr −2 in South Asia after 1960. Consequently, the local terrestrial water storage depletion rate is substantially enlarged by irrigation expansion (for example, from −2.559 ( ± 0.094) to −16.008 ( ± 0.557) mm yr −1 ). Our results attribute the land water loss to irrigation expansion and climate change, calling for immediate solutions to tackle the negative trends.

irrigation↗

Leaf isoprene and monoterpene emission distribution across hyperdominant tree genera in the Amazon basin

Tropical forests are acknowledged to be the largest global source of isoprene (C 5 H 8 ) and monoterpenes (C 10 H 16 ) emissions, with current synthesis studies suggesting few tropical species emit isoprenoids (20-38%) and do so with highly variable emission capacities, including within the same genera. This apparent lack of a clear phylogenetic thread has created difficulties both in linking isoprenoid function with evolution and for the development of accurate biosphere-atmosphere models. In this work, we present a systematic emission study of "hyperdominant" tree species in the Amazon Basin. Across 162 individuals, distributed among 25 botanical families and 113 species, isoprenoid emissions were widespread among both early and late successional species (isoprene: 61.9% of the species; monoterpenes: 15.0%; both isoprene and monoterpenes: 9.7%). The hyperdominant species (69) across the top five most abundant genera, which make up about 50% of all individuals in the Basin, had a similar abundance of isoprenoid emitters (isoprene: 63.8%; monoterpenes: 17.4%; both 11.6%). Among the abundant genera, only Pouteria had a low frequency of isoprene emitting species (15.8% of 19 species). In contrast, Protium, Licania, Inga, and Eschweilera were rich in isoprene emitting species (83.3% of 12 species, 61.1% of 18 species, 100% of 8 species, and 100% of 12 species, respectively). Light response curves of individuals in each of the five genera showed light-dependent, photosynthesis-linked emission rates of isoprene and monoterpenes. Importantly, in every genus, we observed species with light-dependent isoprene emissions together with monoterpenes including β-ocimene. These observations support the emerging view of the evolution of isoprene synthases from β-ocimene synthases. Our results have important implications for understanding isoprenoid function-evolution relationships and the development of more accurate Earth System Models.

59 BASIC BIOLOGICAL SCIENCES↗

Wildfire Smoke Demonstrates Significant and Predictable Black Carbon Light Absorption Enhancements

Abstract Black carbon (BC) is estimated to have the second largest anthropogenic radiative forcing in earth‐systems models (ESMs), but there is significant uncertainty in its impact due to complex mixing with organics. Laboratory‐generated particles show that co‐mixed non‐absorbing material enhances absorption by BC by a factor of 2–3.5 as predicted by optical models. However, weak or no enhancements are often reported for field studies. The cause of lower‐than‐expected absorption is not well understood and implies a lower radiative impact of BC compared to how many ESMs currently treat aerosols. By analyzing BC aerosol particle‐by‐particle we reconcile observed and expected absorption for ambient smoke plumes varying in geographic origin, fuel types, burn conditions, atmospheric age and transport. Although particle‐by‐particle tracking is computationally prohibitive for sophisticated ESMs we show that realistic BC absorption is reliably estimated by bulk properties of the plume providing a suitable parameterization to constrain black carbon radiative forcing.

54 ENVIRONMENTAL SCIENCES↗

Concurrent Changes in Extreme Hydroclimate Events in the Colorado River Basin

Extreme events resulting in catastrophic damage have more than doubled in the last five years, costing hundreds of lives and thousands of homes, and heavily undermining regional economic stability. At present, most of these hydroclimatic extreme events are documented by the media as individual events; however, in scientific terms, many are better understood as concurrent events—concurrent extremes of both temperature and precipitation (e.g., drought, floods). This paper considers concurrent changes in hydroclimate extremes, including heatwaves, drought, flooding, and low flows, in six historical-to-future (1970–1999, 2070–2099) Earth System Model (ESM) climate scenarios for the Colorado River basin. Results indicate that temperature-driven Impacts (heatwaves, drought) have the strongest responses while precipitation-driven Impacts have weaker responses. All Impacts exhibit an increase in magnitude from synoptic to annual time scales, with heatwaves increasing in strength about three times at the annual time scale versus the synoptic, while low flows only increase slightly. Critical watersheds in the Colorado were identified, highlighting the Blue River basin, Uncompahgre, East Taylor, Salt/Verde watersheds, locations of important water infrastructures, water resources, and hydrological research. Our results indicate that concurrent extreme hydroclimate events are projected to increase in the future and intensify within critical regions of the Colorado River basin. Considering extreme hydroclimate events concurrently is an important step towards linking economic and social effects of these events and their associated instabilities on a regional scale.

54 ENVIRONMENTAL SCIENCES↗