Atmospheric rivers (ARs) are filamentary structures within the atmosphere that account for a substantial portion of poleward moisture transport and play an important role in Earth’s hydroclimate. However, there is no one quantitative definition for what constitutes an atmospheric river, leading to uncertainty in quantifying how these systems respond to global change. This study seeks to better understand how different AR detection tools (ARDTs) respond to changes in climate states utilizing single-forcing climate model experiments under the aegis of the Atmospheric River Tracking Method Intercomparison Project (ARTMIP). We compare a simulation with an early Holocene orbital configuration and another with CO2 levels of the Last Glacial Maximum to a pre-industrial control simulation to test how the ARDTs respond to changes in seasonality and mean climate state, respectively. We find good agreement among the algorithms in the AR response to the changing orbital configuration, with a poleward shift in AR frequency that tracks seasonal poleward shifts in atmospheric water vapor and zonal winds. In the low CO2 simulation, the algorithms generally agree on the sign of AR changes but there is substantial spread in their magnitude, indicating that mean-state changes lead to larger uncertainty. This disagreement likely arises primarily from differences between algorithms in their thresholds for water vapor and its transport used for identifying ARs. These findings warrant caution in ARDT selection for paleoclimate and climate change studies in which there is a change to the mean climate state, as ARDT selection contributes substantial uncertainty in such cases.

Global kilometer-scale models are the future of Earth system models as they can explicitly simulate organized convective storms and their associated extreme weather. Here, we comprehensively examined tropical mesoscale convective system (MCS) characteristics in the DYAMOND (DYnamics of the atmospheric general circulation modeled on non-hydrostatic domains) models for both summer and winter phases by applying eight different feature trackers to the simulations and satellite observations. Although different trackers produce substantial differences (a factor of 2-3) in observed MCS frequency and their contribution to total precipitation, model-observation differences in MCS statistics are more consistent among the trackers. DYAMOND models are generally skillful in simulating tropical mean MCS frequency, with multi-model mean biases of 2.9% over land and -0.5% over ocean. However, most models underestimate the MCS precipitation amount (23%) and their contribution to total precipitation (17%) relative to observations. These biases show large inter-model variability, but are generally smaller over land (13%) than over ocean (21%) on average. MCS diurnal cycle and cloud shield characteristics are better simulated than precipitation. Most models overestimate MCS precipitation intensity and underestimate stratiform rain contribution (up to a factor of 2), particularly over land. Models also predict a wide range of precipitable water in the tropics compared to reanalysis and satellite observations, and many models simulate a greater sensitivity of MCS precipitation intensity to precipitable water. The MCS metrics developed in this work provide process-oriented diagnostics for future model development efforts.

Mesoscale convective systems (MCSs) play an important role in modulating the global hydrological cycle, general circulation, and radiative energy budget. In this study, we evaluate MCS simulations in the second version of U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv2). E3SMv2 atmosphere model (EAMv2) is run at the uniform 0.25° horizontal resolution. We track MCSs consistently in the model and observations using the PyFLEXTRKR algorithm, which defines MCS based on both cloud-top brightness temperature (Tb) and surface precipitation. Results from using Tb only to define MCS, commonly used in previous studies, are also discussed. Furthermore, sensitivity experiments are performed to examine the impact of new cloud and convection parameterizations developed for EAMv3 on simulated MCSs. Our results show that EAMv2 simulated MCS precipitation is largely underestimated in the tropics and contiguous United States. This is mainly attributed to the underestimated precipitation intensity in EAMv2. In contrast, the simulated MCS frequency becomes more comparable to observations if MCSs are defined only based on cloud-top Tb. The Tb-based MCS tracking method, however, includes many cloud systems with very weak precipitation which conflicts with the MCS definition. This result illustrates the importance of accounting for precipitation in evaluating simulated MCSs. We also find that the new physics parameterizations help increase the relative contribution of convective precipitation to total precipitation in the tropics, but the simulated MCS properties are overall not significantly improved. This suggests that simulating MCSs will remain a challenge for the next version of E3SM.

Radiation-topography interaction plays an important role in the surface energy balance over the Tibetan Plateau (TP). However, the impacts of such interaction over the TP on climate locally and in the Asian regions remain unclear. This study uses the Energy Exascale Earth System Model (E3SM) to evaluate the regional and teleconnected impacts of radiation-topography interaction over the TP. Land-atmosphere coupled experiments show that topography regulates the surface energy balance, snow processes, and surface climate over the TP across seasons. Accounting for radiation-topography interaction overall improves E3SM’s performance in simulating surface climate. The winter cold bias in air temperature decreases from -4.48 K to -3.70 K, and the wet bias in summer precipitation is mitigated in southern TP. The TP’s radiation-topography interaction further reduces the South and East Asian summer precipitation biases. Our results demonstrate the topographic roles in regional climate over the TP and highlight its teleconnected climate impacts.

Flow routing is a fundamental process of Earth System Models’ (ESMs) river component. Traditional flow routing models rely on Cartesian rectangular meshes, which exhibit limitations, particularly when coupled with unstructured mesh-based ocean components. They also lack the support for regionally refined models (RRMs). While previous studies have highlighted the potential benefits of unstructured meshes for flow routing, their widespread application and comprehensive evaluation within ESMs remain limited. This study extends the river component of the Energy Exascale Earth System Model (E3SM) to unstructured Voronoi meshes. We evaluated the model’s performance in simulating river discharge and water depth across three watersheds spanning the Arctic, temperate, and tropical regions. The results show that while providing several benefits, unstructured mesh-based flow routing can achieve comparable performance to structured mesh-based routing, and their difference is often less than 10%. Although the unstructured mesh-based method could address several existing limitations, this research also shows that additional improvements in the numerical method are needed to fully exploit the advantages of unstructured mesh for hydrologic and ESMs.

The linear relationship between gross primary productivity (GPP) and evapotranspiration (ET), evidenced by site-scale observations, is well recognized as an indicator of the close interactions between carbon and hydrologic processes in terrestrial ecosystems. However, it is not clear whether this relationship holds at the catchment scale, and if so, what are the controlling factors of its slope and intercept. This study proposes and examines a generalized GPP-ET relationship at 380 near-natural catchments across various climatic and landscape conditions in the contiguous U.S., based on monthly remote sensing-based GPP data, vegetation phenology, and several hydrometeorological variables. We demonstrate the validity of this GPP-ET relationship at the catchment scale, with Pearson’s r ≥ 0.6 for 97% of the 380 catchments. Furthermore, we propose a regionalization strategy for estimating the slope and intercept of the generalized GPP-ET relationship at the catchment scale by linking the parameter values a priori with hydrometeorological data. We validate the monthly GPP predicted from the relationship and regionalized parameters against remote-sensing based GPP product, yielding Kling-Gupta Efficient (KGE) values ≥ 0.5 for 92% of the catchments. Finally, we verify the relationship and its parameter regionalization at 35 AmeriFlux sites with KGE ≥ 0.5 for 25 sites, demonstrating that the new relationship is transferable across the site, catchment, and regional scales. The relationship will be valuable for diagnosing coupled water–carbon simulations in land surface and Earth system models and constraining remote-sensing based estimation of monthly ET.

This work describes the implementation and evaluation of the Slab Ocean Model com16 ponent of the Energy Exascale Earth System Model version 2 (E3SMv2-SOM) and its application to understanding the climate sensitivity to ocean heat transports (OHTs) and CO2 forcing. E3SMv2-SOM reproduces the baseline climate and Equilibrium Climate Sensitivity (ECS) of the fully coupled E3SMv2 experiments reasonably well, with a pattern correlation close to 1 and a global mean bias of less than 1% of the fully coupled surface temperature and precipitation. Sea ice extent and volume are also well reproduced in the SOM. Consistent with general model behavior, the ECS estimated from the SOM (4.5K) exceeds the effective climate sensitivity obtained from extrapolation to equilibrium in the fully coupled model (4.0K). The E3SMv2 baseline climate also shows a large sensitivity to OHT strengths, with a global surface temperature difference of about 4.0◦ C between high-/low-OHT experiments with prescribed forcings derived from fully coupled experiments with realistic/weak ocean circulation strengths. Similar to their forc ng pattern, the surface temperature response occurs mainly over the subpolar regions in both hemispheres. However, the Southern Ocean shows more surface temperature sensitivity to high/low-OHT forcing due to a positive/negative shortwave cloud radiative effect caused by decreases/increases in mid-latitude marine low-level clouds. This large temperature sensitivity also causes an overcompensation between the prescribed OHTs and atmosphere heat transports. The SOM’s ECS estimate is also sensitive to the prescribed OHT and the associated baseline climate it is initialized from; the high-OHT ECS is 0.5K lower than the low-OHT ECS.

In the face of a changing climate, the understanding, predictions and projections of natural and human systems are increasingly crucial to prepare and cope with extremes and cascading hazards, determine unexpected feedbacks and potential tipping points, inform long-term adaptation strategies, and guide mitigation approaches. Increasingly complex socio-economic systems require enhanced predictive information to support advanced practices. Such new predictive challenges drive the need to fully capitalize on ambitious scientific and technological opportunities. These include the unrealized potential for very high-resolution modeling of global-to-local Earth system processes across timescales, a reduction of model biases, enhanced integration of human systems and the Earth Systems, better quantification of predictability and uncertainties; expedited science-to-service pathways and co-production of actionable information with stakeholders. Enabling technological opportunities include exascale computing, advanced data storage, novel observations and powerful data analytics, including artificial intelligence and machine learning. Looking to generate community discussions on how to accelerate progress on U.S. climate predictions and projections, representatives of Federally-funded U.S. modeling groups outline here perspectives on a six-pillar national approach grounded in climate science that builds on the strengths of the U.S. modeling community and agency goals. This calls for an unprecedented level of coordination to capitalize on transformative opportunities, augmenting and complementing current modeling center capabilities and plans to support agency missions. Tangible outcomes include projections with horizontal spatial resolutions finer than 10 km, representing extremes and associated risks in greater detail, reduced model errors, better predictability estimates, and more customized projections to support the next generation of climate services.

Mesoscale convective systems (MCSs) produce over 50% of tropical precipitation and account for the majority of extreme rainfall and flooding events. MCSs are considered the building blocks of larger-scale convectively coupled equatorial waves (CCEWs). While CCEWs can provide favorable environments for convection, how CCEWs can systematically impact organized convection and thereby MCS characteristics is less clear. We examine this question by analyzing a global MCS tracking dataset. During the active phase of CCEWs, MCS frequency increases and MCSs rain harder, produce more lifetime total rain, and grow larger in size. The probability of extreme MCSs also elevates. These changes are most pronounced when MCSs are associated with Kelvin waves and tropical depression-type waves while least so with the Madden-Julian Oscillation. These results can be benchmarks to improve model representation of MCS interactions with large-scale circulations and can be leveraged operationally for extended forecasts of high-impact MCSs.

Forecasting rapid intensification (RI) of tropical cyclones (TC) is a mission known for large errors. One under-researched factor that affects TC intensification is sea surface salinity (SSS), which is an important proxy for density stratification in certain ocean regions and can affect the surface enthalpy flux under a strengthening hurricane. To investigate the effects of SSS on TC RI, we use a previously built statistical model consisting of a variety of machine learning (ML) methods. A calibrator was trained on top of the ML models to correct probability forecasts. The ML model performance is improved with the addition of SSS in the Eastern North Pacific and the Caribbean subregion of the North Atlantic. Limited improvement is found in the Western North Pacific. In the Indian Ocean, SSS is also notably correlated with RI occurrence, but the TC samples are not sufficient to train ML models.

Large-scale dynamical and thermodynamical processes are common environmental drivers of extreme weather events. However, such large-scale environmental conditions often display systematic biases in climate simulations, posing challenges to evaluating extreme weather events and associated risks in current and future climate. In this paper, a machine learning (ML) approach was employed to bias correct the large-scale wind, temperature, and humidity simulated by the E3SM atmosphere model at $\sim 1^\circ$ resolution. The usefulness of the proposed ML approach for extreme weather analysis was demonstrated with a focus on three extreme weather events, including tropical cyclones (TCs), extratropical cyclones (ETCs), and atmospheric rivers (ARs). We show that the ML model can effectively reduce climate bias in large-scale wind, temperature, and humidity while preserving their responses to imposed climate change perturbations. The bias correction is found to directly improve the water vapor transport associated with ARs, and the representations of thermodynamical flows associated with ETCs. When the bias-corrected large-scale winds are used to drive a synthetic TC track forecast model over the Atlantic basin, the resulting TC track density agrees better with that of the TC track model driven by observed winds. In addition, the ML model insignificantly interferes with the mean climate change signals of large-scale storm environments as well as the occurrence and intensity of three extreme events. This study suggests that the proposed ML approach can be used to improve the downscaling of extreme weather events by providing more realistic large-scale storm environments simulated by low-resolution climate models.

Coastal wetlands play an important role in the global water and biogeochemical cycles. Climate change is making them more difficult to adapt to the fluctuation of sea levels and other environment changes. Given the importance of eco-geomorphological processes for coastal wetland resilience, many eco-geomorphology models differing in complexity and numerical schemes have been developed in recent decades. But their divergent estimates on the response of coastal wetlands to climate change indicate that substantial structural uncertainties exist in these models. To investigate the structural uncertainty of coastal wetland eco-geomorphology models, we developed a multi-algorithm model framework of eco-geomorphological processes, such as mineral accretion and organic matter accretion, within a single hydrodynamics model. The framework is designed to explore possible ways to represent coastal wetland eco-geomorphology in Earth system models and reduce the related uncertainties in global applications. We tested this model framework at three representative coastal wetland sites: two saltmarsh wetland (Venice Lagoon and Plum Island Estuary) and a mangrove wetland (Hunter Estuary). Through the model-data comparison, we showed the importance to use a multi-algorithm ensemble approach for more robust predictions of the evolution of coastal wetlands. We also find that more observations of mineral and organic matter accretion at different elevations of coastal wetlands and evaluation of the coastal wetland models at different sites of diverse environments can help reduce the model uncertainty.