This paper provides an overview of the United States (U.S.) Department of Energy’s (DOE’s) Energy Exascale Earth System Model version 2.1 with an Arctic regionally refined mesh (RRM), hereafter referred to as E3SMv2.1-Arctic, for the atmosphere (25 km), land (25 km), and ocean/ice (10 km) components. We evaluate the atmospheric component and its interactions with land, ocean, and cryosphere by comparing the RRM (E3SM2.1-Arctic) historical simulations (1950-2014) with the uniform low-resolution (LR) counterpart, reanalysis products, and observational datasets. The RRM generally reduces biases in the LR model, improving simulations of Arctic large-scale mean fields, such as precipitation, atmospheric circulation, clouds, atmospheric river frequency, and sea ice dynamics. However, the RRM introduces a seasonally dependent surface air temperature bias, reducing the LR cold bias in summer but enhancing the LR warm bias in winter. The RRM also underestimates winter sea ice area and volume, consistent with its strong winter warm bias. Radiative feedback analysis shows similar climate feedback strengths in both RRM and LR, with the RRM exhibiting a more positive surface albedo feedback and contributing to a stronger surface warming than LR. These findings underscore the importance of high-resolution modeling for advancing our understanding of Arctic climate changes and their broader global impacts, although some persistent biases appear to be independent of model resolution at 10-100 km scales.

Hydroclimate and terrestrial hydrology greatly influence the local community, ecosystem, and economy in Alaska and Yukon River Basin. A high-resolution re-simulation of the historical climate in Alaska can provide an important benchmark for climate change studies. In this study, we utilized the Regional Arctic Systems Model (RASM) and conducted coupled land-atmosphere modeling for Alaska and Yukon River Basin at 4-km grid spacing. In RASM, the land model was replaced with the Community Terrestrial Systems Model (CTSM) given its comprehensive process representations for cold regions. The microphysics schemes in the Weather Research and Forecast (WRF) atmospheric model were manually tuned for optimal model performance. This study aims to maintain good model performance for both hydroclimate and terrestrial hydrology, especially streamflow, which was rarely a priority in coupled models. Therefore, we implemented a strategy of iterative testing and re-optimization of CTSM. A multi-decadal climate dataset (1990-2021) was generated using RASM with optimized land parameters and manually tuned WRF microphysics. When evaluated against multiple observational datasets, this dataset well captures the climate statistics and spatial distributions for five key weather variables and hydrologic fluxes, including precipitation, air temperature, snow fraction, evaporation-to-precipitation ratios, and streamflow. The simulated precipitation shows wet bias during the spring season and simulated air temperatures exhibit dampened seasonality with warm biases in winter and cold biases in summer. We used transfer entropy to investigate the discrepancy in connectivity of hydrologic fluxes between the offline CTSM and coupled models, which contributed to their discrepancy in streamflow simulations.

A set of decadal simulations has been completed and evaluated for gains using the Regional Arctic System Model (RASM) to dynamically downscale data from a global Earth system model (ESM) and two atmospheric reanalyses. RASM is a fully coupled atmosphere - land - ocean - sea ice regional Earth system model. Nudging to the forcing data is applied to approximately the top half of the atmosphere. RASM simulations were also completed with a modification to the atmospheric physics for evaluating changes to the modeling system. The results show that for the top half of the atmosphere, the RASM simulations follow closely to that of the forcing data, regardless of the forcing data. The results for the lower half of the atmosphere, as well as the surface, show a clustering of atmospheric state and surface fluxes based on the modeling system. At all levels of the atmosphere the imprint of the weather from the forcing data is present as indicated in the pattern of the monthly and annual means. Biases, in comparison to reanalyses, are evident in the ESM forced simulations for the top half of the atmosphere but are not present in the lower atmosphere. This suggests that bias correction is not needed for fully-coupled dynamical downscaling simulations. While the RASM simulations tended to go to the same mean state for the lower atmosphere, there is a different evolution of the weather across the ensemble of simulations. These differences in the weather result in variances in the sea ice and oceanic states.