A team effort to develop a Community Integrated Earth System Model (CIESM) was initiated in China in 2012. The model was based on NCAR Community Earth System Model (Version 1.2.1) with several novel developments and modifications aimed to overcome some persistent systematic biases, such as the double Intertropical Convergence Zone problem and underestimated marine boundary layer clouds. Aerosols' direct and indirect effects are prescribed using the MACv2‐SP approach and data sets. The spin‐up of a 500‐year preindustrial simulation and three historical simulations are described and evaluated. Prominent improvements include alleviated double Intertropical Convergence Zone problem, increased marine boundary layer clouds, and better El Niño Southern Oscillation amplitude and periods. One deficiency of the model is the significantly underestimated Arctic and Antarctic sea ice in warm seasons. The historical warming is about 0.55 °C greater than observations toward 2014. CIESM has an equilibrium climate sensitivity of 5.67 K, mainly resulted from increased positive shortwave cloud feedback. Our efforts on porting and redesigning CIESM for the heterogeneous Sunway TaihuLight supercomputer are also introduced, including some ongoing developments toward a future version of the model.

In this study, we compare eight recently developed snow depth products over Arctic sea ice, which use satellite observations, modeling, or a combination of satellite and modeling approaches. These products are further compared against various ground-truth observations, including those from ice mass balance observations and airborne measurements. Large mean snow depth discrepancies are observed over the Atlantic and Canadian Arctic sectors. The differences between climatology and the snow products early in winter could be in part a result of the delaying in Arctic ice formation that reduces early snow accumulation, leading to shallower snowpacks at the start of the freeze-up season. These differences persist through spring despite overall more winter snow accumulation in the reanalysis-based products than in the climatologies. Among the products evaluated, the University of Washington (UW) snow depth product produces the deepest spring (March–April) snowpacks, while the snow product from the Danish Meteorological Institute (DMI) provides the shallowest spring snow depths. Most snow products show significant correlation with snow depths retrieved from Operational IceBridge (OIB) while correlations are quite low against buoy measurements, with no correlation and very low variability from University of Bremen and DMI products. Inconsistencies in reconstructed snow depth among the products, as well as differences between these products and in situ and airborne observations, can be partially attributed to differences in effective footprint and spatial–temporal coverage, as well as insufficient observations for validation/bias adjustments. Our results highlight the need for more targeted Arctic surveys over different spatial and temporal scales to allow for a more systematic comparison and fusion of airborne, in situ and remote sensing observations.

With semiconductor technology gradually approaching its physical and thermal limits, recent supercomputers have adopted major architectural changes to continue increasing the performance through more power-efficient heterogeneous many-core systems. Examples include Sunway TaihuLight that has four management processing elements (MPEs) and 256 computing processing elements (CPEs) inside one processor and Summit that has two central processing units (CPUs) and six graphics processing units (GPUs) inside one node. Meanwhile, current high-resolution Earth system models that desperately require more computing power generally consist of millions of lines of legacy code developed for traditional homogeneous multicore processors and cannot automatically benefit from the advancement of supercomputer hardware. As a result, refactoring and optimizing the legacy models for new architectures become key challenges along the road of taking advantage of greener and faster supercomputers, providing better support for the global climate research community and contributing to the long-lasting societal task of addressing long-term climate change. This article reports the efforts of a large group in the International Laboratory for High-Resolution Earth System Prediction (iHESP) that was established by the cooperation of Qingdao Pilot National Laboratory for Marine Science and Technology (QNLM), Texas A&M University (TAMU), and the National Center for Atmospheric Research (NCAR), with the goal of enabling highly efficient simulations of the high-resolution (25 km atmosphere and 10 km ocean) Community Earth System Model (CESM-HR) on Sunway TaihuLight. The refactoring and optimizing efforts have improved the simulation speed of CESM-HR from 1 SYPD (simulation years per day) to 3.4 SYPD (with output disabled) and supported several hundred years of pre-industrial control simulations. With further strategies on deeper refactoring and optimizing for remaining computing hotspots, as well as redesigning architecture-oriented algorithms, we expect an equivalent or even better efficiency to be gained on the new platform than traditional homogeneous CPU platforms. The refactoring and optimizing processes detailed in this paper on the Sunway system should have implications for similar efforts on other heterogeneous many-core systems such as GPU-based high-performance computing (HPC) systems.

To increase our capability to predict Arctic sea ice and climate, we have developed a coupled atmosphere‐sea ice‐ocean model configured for the pan‐Arctic with sufficient flexibility. The Los Alamos Sea Ice Model is coupled with the Weather Research and Forecasting Model and the Regional Ocean Modeling System in the Coupled Ocean‐Atmosphere‐Wave‐Sediment Transport modeling system. It is well known that dynamic models used to predict Arctic sea ice at short‐term periods strongly depend on model initial conditions. Parallel Data Assimilation Framework is implemented into the new modeling system to assimilate sea ice observations and generate skillful model initialization, which aid in the prediction procedures. The Special Sensor Microwave Imager/Sounder sea ice concentration, the CyroSat‐2, and Soil Moisture and Ocean Salinity sea ice thickness are assimilated with the localized error subspace transform ensemble Kalman filter. We conduct Arctic sea ice prediction for the melting seasons of 2017 and 2018. Predictions with improved initial sea ice conditions show reasonable sea ice evolution and small biases in the minimum sea ice extent, although the ice refreezing is delayed. Our prediction experiments suggest that the use of appropriate uncertainty for the observed sea ice thickness can lead to improved spatial distribution of the initial ice thickness and thus the predicted sea ice distribution. Our new modeling system initialized by the output of the National Centers for Environmental Prediction Climate Forecast System seasonal forecasts with data assimilation can significantly increase the sea ice prediction skills in sea ice extent for the entire Arctic as well as in the Northern Sea Route compared with the predictions by the National Centers for Environmental Prediction Climate Forecast System.

With the aid of the newly developed ‘Sunway’ heterogeneous-architecture supercomputer, which has world-leading HPC (high-performance computer) capability, a series of high-resolution coupled Earth system models (SW-HRESMs) with up to 5 km of atmosphere and 3 km of ocean have been developed. These models can meet the needs of multiscale interaction studies with different computational costs. Here we describe the progress of SW-HRESMs development, with an overview of the major advancements made by the international Earth science community in HR-ESMs. We also show the preliminary results of SW-HRESMs with regard to capturing major weather-climate extremes in the atmosphere and ocean, stressing the importance of permitted clouds and ocean submesoscale eddies in modeling tropical cyclones and eddy-mean flow interactions, and paving the way for further model development to resolve finer scales with even higher resolution and more realistic physics. Finally, in addition to increasing model resolution, the development procedure for a non-hydrostatic cloud and ocean submesoscale resolved ESM is discussed, laying out the major scientific directions of such a huge modeling advancement.

The Soil Moisture Active Passive (SMAP) mission was designed to provide a global mapping of soil moisture (SM) measured by L-band passive and active microwave sensors. In this study, we evaluate the newly released SMAP enhanced SM products over the Tibetan Plateau by performing comparisons among SMAP standard products, in-situ observations and Community Land Model (CLM) simulations driven by high-resolution meteorological forcing. At local scales, the enhanced SMAP products, the standard products and CLM simulations all generally compare well with the in-situ observations. The SMAP products show stronger correlations (0.64–0.88) but slightly larger unbiased root mean square errors (ubRMSE, ~0.06) relative to the CLM simulations (0.58–0.79 and 0.037–0.047, for correlation and ubRMSE, respectively). At the regional scale, both SMAP products show similar spatial distributions of SM on the TP (Tibetan Plateau), although, as expected, the enhanced product provides more fine details. The SMAP enhanced product is in good agreement with model simulations with respect to temporal and spatial variations in SM over most of the TP. Regions with low correlation between SMAP enhanced products and model simulations are mainly located in the northwestern TP and regions of complex topography, where meteorological stations are sparse and non-existent or elevation is highly variable. In such remote regions, CLM simulations may be problematic due to inaccurate land cover maps and/or uncertainties in meteorological forcing. The independent, high-resolution observations provided by SMAP could help to constrain the model simulation and, ultimately, improve the skill of models in these problematic regions.