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.

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.

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.

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.

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.

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.

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.

A document by Chang Liao. Click on the document to view its contents.

In recent decades, Arctic sea ice coverage experienced a drastic decline in winter, when sea ice is expected to recover following the melting season. Using observations and climate model simulations, we found a robust frequency increase in atmospheric rivers (ARs, intense corridors of moisture transport) over Barents-Kara Seas and the neighboring central Arctic (ABK) in early winter. The extensive moisture carried by more frequent ARs has intensified surface downward longwave radiation and liquid rainfall, caused stronger melting of thin, fragile ice cover, and slowed the seasonal recovery of sea ice, contributing to the sea ice cover decline in ABK. A series of model ensemble experiments suggests that, in addition to a uniform AR increase in response to anthropogenic forcing, the contribution of tropical Pacific variability is indispensable in the observed Arctic AR changes. These findings have significant implications for understanding the rapidly changing Arctic hydroclimate and the cryosphere.

The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low-resolution (roughly 110 km) and high-resolution (roughly 25 km) configurations — as established by the High Resolution Model Intercomparison Project protocol — to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes to the water cycle come from reductions in precipitation and evapotranspiration, the latter of which results from decreases in evaporative fraction. While the HR simulation is warmer than the LR, moisture convergence decreases despite the increased atmospheric water vapor, suggesting circulation biases are an important factor. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.

Flow direction modeling consists of (1) an accurate representation of the river network and (2) digital elevation model (DEM) processing to preserve characteristics with hydrological significance. In part 1 of our study, we presented a mesh-independent approach to representing river networks on different types of meshes. This follow-up part 2 study presents a novel DEM processing approach for flow direction modeling. This approach consists of (1) a topological relationship-based hybrid breaching-filling method to conduct stream burning for the river network and (2) a modified depression removal method for rivers and hillslopes. Our methods minimize modifications to surface elevations and provide a robust two-step procedure to remove local depressions in DEM. They are mesh-independent and can be applied to both structured and unstructured meshes. We applied our new methods to the Susquehanna River Basin with different model configurations. The results show that topological relationship-based stream burning and depression-filling methods can reproduce the correct river networks, providing high-quality flow direction and other characteristics for hydrologic and Earth system models.