Explore the vast range of titles available.

The main advancements of the Beijing Climate Center (BCC) climate system model from phase 5 of the Coupled Model Intercomparison Project (CMIP5) to phase 6 (CMIP6) are presented, in terms of physical parameterizations and model performance. BCC-CSM1.1 and BCC-CSM1.1m are the two models involved in CMIP5, whereas BCC-CSM2-MR, BCC-CSM2-HR, and BCC-ESM1.0 are the three models configured for CMIP6. Historical simulations from 1851 to 2014 from BCC-CSM2-MR (CMIP6) and from 1851 to 2005 from BCC-CSM1.1m (CMIP5) are used for models assessment. The evaluation matrices include the following: (a) the energy budget at top-of-atmosphere; (b) surface air temperature, precipitation, and atmospheric circulation for the global and East Asia regions; (c) the sea surface temperature (SST) in the tropical Pacific; (d) sea-ice extent and thickness and Atlantic Meridional Overturning Circulation (AMOC); and (e) climate variations at different timescales, such as the global warming trend in the 20th century, the stratospheric quasi-biennial oscillation (QBO), the Madden–Julian Oscillation (MJO), and the diurnal cycle of precipitation. Compared with BCC-CSM1.1m, BCC-CSM2-MR shows significant improvements in many aspects including the tropospheric air temperature and circulation at global and regional scales in East Asia and climate variability at different timescales, such as the QBO, the MJO, the diurnal cycle of precipitation, interannual variations of SST in the equatorial Pacific, and the long-term trend of surface air temperature.

In the late twentieth century, global mean surface air temperature especially on land is continuously warming. Our analyses show that the global mean of dust increased since 1980, using the Modern-Era Retrospective Analysis version 2 for Research and Applications (MERRA-2) reanalysis data. This variation of global dust is mainly contributed by the dust increase outside of dust core areas (i.e. high dust mass concentration region). The causes to result in global dust variations are explored. In dust core areas, surface wind is the primary driving factor for surface dust, both of which show no remarkable trends of increase or decrease since 1980. In areas outside of the core areas, especially in arid and semi-arid areas in North and Middle Asia, surface air temperature warming is the primary impact factor causing the dust increase. An increase in surface air temperature is accompanied by enhancement of atmospheric instability which can trigger more upward motion and bring more dust. All 9 Earth System Models (ESMs) for the Aerosol Chemistry Model Intercomparison Project (AerChemMIP) reproduce the reasonable spatial distribution and seasonal cycle of dust in the present day. But only a few models such as BCC-ESM1 and GFDL-ESM4 simulate the increasing trend of dust similar to MERRA-2. While the primary impact of wind in dust core areas, and surface temperature outside of the core areas, especially in middle to high latitudes in Eurasian continent, are presented in most ESMs.

Climate models participated in the High Resolution Model Intercomparison Project (HighResMIP) of Coupled Model Intercomparison Project 6 (CMIP6) are evaluated to reveal the impact of enhanced resolution in simulating the climatological distribution of precipitation over China. The multi-model mean (MME) of five models with 30–50 km horizontal resolution in the atmosphere (MME-50) had a better performance in reproducing the observed spatial patterns of precipitation over Northwest China and Southwest China than the MME of their lower-resolution (70–140 km) models (MME-100). Such an improvement is mainly attributed to the topographical rainfall reproduced by the higher-resolution models, which have the superiority of reproducing the local vertical circulation over the complex terrain. The MME-50 also improves the skill score of the 850-hPa southwesterly of the Indian-Burma trough relative to the MME-100, which may contribute to better simulation skill of precipitation over Southwest China. The MME-50 (0.92) has a close skill score to the MME-100 (0.91) in the simulation of East Asian summer monsoon, which explains why the MME-50 cannot improve the simulation skill of the precipitation over Southeast China and Northern China. The skill score of the precipitation over the Tibetan Plateau (TP) simulated by the MME-50 is even lower than the MME-100, revealing that simulating the climate over the high plateau is still a challenge for the high-resolution models.

The Beijing Climate Center Earth System Model version 1 (BCC-ESM1) is the first version of a fully coupled Earth system model with interactive atmospheric chemistry and aerosols developed by the Beijing Climate Center, China Meteorological Administration. Major aerosol species (including sulfate, organic carbon, black carbon, dust, and sea salt) and greenhouse gases are interactively simulated with a whole panoply of processes controlling emission, transport, gas-phase chemical reactions, secondary aerosol formation, gravitational settling, dry deposition, and wet scavenging by clouds and precipitation. Effects of aerosols on radiation, cloud, and precipitation are fully treated. The performance of BCC-ESM1 in simulating aerosols and their optical properties is comprehensively evaluated as required by the Aerosol Chemistry Model Intercomparison Project (AerChemMIP), covering the preindustrial mean state and time evolution from 1850 to 2014. The simulated aerosols from BCC-ESM1 are quite coherent with Coupled Model Intercomparison Project Phase 5 (CMIP5)-recommended data, in situ measurements from surface networks (such as IMPROVE in the US and EMEP in Europe), and aircraft observations. A comparison of modeled aerosol optical depth (AOD) at 550 nm with satellite observations retrieved from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging SpectroRadiometer (MISR) and surface AOD observations from the AErosol RObotic NETwork (AERONET) shows reasonable agreement between simulated and observed AOD. However, BCC-ESM1 shows weaker upward transport of aerosols from the surface to the middle and upper troposphere, likely reflecting the deficiency of representing deep convective transport of chemical species in BCC-ESM1. With an overall good agreement between BCC-ESM1 simulated and observed aerosol properties, it demonstrates a success of the implementation of interactive aerosol and atmospheric chemistry in BCC-ESM1.

By conducting several sets of hindcast experiments using the Beijing Climate Center Climate System Model, which participates in the Sub-seasonal to Seasonal (S2S) Prediction Project, we systematically evaluate the model’s capability in forecasting MJO and its main deficiencies. In the original S2S hindcast set, MJO forecast skill is about 16 days. Such a skill shows significant seasonal-to-interannual variations. It is found that the model-dependent MJO forecast skill is more correlated with the Indian Ocean Dipole (IOD) than with the El Niño–Southern Oscillation. The highest skill is achieved in autumn when the IOD attains its maturity. Extended skill is found when the IOD is in its positive phase. MJO forecast skill’s close association with the IOD is partially due to the quickly strengthening relationship between MJO amplitude and IOD intensity as lead time increases to about 15 days, beyond which a rapid weakening of the relationship is shown. This relationship transition may cause the forecast skill to decrease quickly with lead time, and is related to the unrealistic amplitude and phase evolutions of predicted MJO over or near the equatorial Indian Ocean during anomalous IOD phases, suggesting a possible influence of exaggerated IOD variability in the model. The results imply that the upper limit of intraseasonal predictability is modulated by large-scale external forcing background state in the tropical Indian Ocean. Two additional sets of hindcast experiments with improved atmosphere and ocean initial conditions (referred to as S2S_IEXP1 and S2S_IEXP2, respectively) are carried out, and the results show that the overall MJO forecast skill is increased to 21–22 days. It is found that the optimization of initial sea surface temperature condition largely accounts for the increase of the overall MJO forecast skill, even though the improved initial atmosphere conditions also play a role. For the DYNAMO/CINDY field campaign period, the forecast skill increases to 27 days in S2S_IEXP2. Nevertheless, even with improved initialization, it is still difficult for the model to predict MJO propagation across the western hemisphere–western Indian Ocean area and across the eastern Indian Ocean–Maritime Continent area. Especially, MJO prediction is apparently limited by various interrelated deficiencies (e.g., overestimated IOD, shorter-than-observed MJO life cycle, Maritime Continent prediction barrier), due possibly to the model bias in the background moisture field over the eastern Indian Ocean and Maritime Continent. Thus, more efforts are needed to correct the deficiency in model physics in this region, in order to overcome the well-known Maritime Continent predictability barrier.

Gravity wave parameterizations contribute to uncertainties in middle atmosphere modeling. To investigate the potential for using machine learning to represent atmospheric gravity waves and the impact of implementing such schemes in a general circulation model (GCM), we train a random forest (RF) emulator on outputs from an existing complex parameterization scheme for orographic gravity wave drag (GWD). The performance of the RF emulator is then evaluated with a focus on stratospheric climatology and variability in climate simulations from the middle atmosphere resolving Beijing Climate Center Atmospheric Circulation Model. In offline tests, the predicted orographic GWD by the RF agrees remarkably well with the target GWD throughout the troposphere and the middle atmosphere. The RF emulator can reproduce the observed climatology of zonal‐mean wind and air temperature in the GCM simulation, as well as its target scheme. Compared to the target orographic GWD parameterization scheme, the RF emulator can reproduce the breakdown of the polar vortex in the Southern Hemisphere stratosphere. This study demonstrates the feasibility of using machine learning to emulate parameterized orographic GWD for modeling the stratosphere with a GCM. Plain Language Summary Machine learning has been utilized to learn atmospheric physical parameterizations such as moist convection, cloud microphysics, radiation, nonorographic gravity waves, etc. However, it remains unknown whether it is feasible to apply machine learning algorithms to parameterize orographic gravity waves and how such schemes perform when interactively coupled in atmospheric general circulation models. Here we employ the widely used random forest algorithm to learn from an existing complex parameterization scheme for orographic gravity wave drag. Then, we implement the random forest scheme in a high‐top atmospheric general circulation model to test its performance. We find that the random forest scheme acts as well as the target parameterization scheme. This study demonstrates the potential of using machine learning to emulate parameterized orographic gravity waves for modeling the middle atmosphere. In future work, it is recommended that machine learning may be used to develop new parameterization schemes directly from observations or high‐resolution model outputs to overcome critical deficiencies in current orographic gravity wave parameterizations, helping to reduce the common biases in simulating the stratosphere with general circulation models. Key Points Parameterized gravity wave drag from unresolved orography can be emulated with a random forest The random forest emulator simulates vertical distributions of the zonal wind and air temperature quite well in a general circulation model (GCM) The GCM configured with the random forest scheme reproduces a reasonable breakdown of the austral polar vortex in the stratosphere

This study aimed to examine the spatio-temporal variations in the atmospheric boundary layer height (ABLH) over the Tarim Basin (TB). Monthly ABLH data from the ERA-Interim dataset from January 1979 to December 2018 were used. Periodicity analysis and the Mann–Kendall Abrupt Changes test were employed to identify the change cycle and abrupt change year of the boundary layer height. The Empirical Orthogonal Function (EOF) method was utilized to determine the spatial distribution of the boundary layer height, and the RF method was used to establish the relationship between the ABLH and influencing factors. The results demonstrated that the highest values of ABLH (over 1900 m) were observed in the middle parts of the study area in June, and the ABLH exhibited a significant increase over the TB throughout the study period. Abrupt changes in the ABLH were also identified in 2004, as well as in 2-, 5-, 9-, and 15-year changing cycles. The first EOF ABLH mode indicated that the middle and northeast regions are relatively high ABLH areas within the study area. Additionally, the monthly variations in ABLH show a moderately positive correlation with air temperature, while exhibiting a negative correlation with air pressure and relative humidity.

This paper provides a comprehensive assessment of Asian summer monsoon prediction skill as a function of lead time and its relationship to sea surface temperature prediction using the seasonal hindcasts of the Beijing Climate Center Climate System Model, BCC CSM1.1(m). For the South and Southeast Asian summer monsoon, reasonable skill is found in the model’s forecasting of certain aspects of monsoon climatology and spatiotemporal variability. Nevertheless, deficiencies such as significant forecast errors over the tropical western North Pacific and the eastern equatorial Indian Ocean are also found. In particular, overestimation of the connections of some dynamical monsoon indices with large-scale circulation and precipitation patterns exists in most ensemble mean forecasts, even for short lead-time forecasts. Variations of SST, measured by the first mode over the tropical Pacific and Indian oceans, as well as the spatiotemporal features over the Ni˜no3.4 region, are overall well predicted. However, this does not necessarily translate into successful forecasts of the Asian summer monsoon by the model. Diagnostics of the relationships between monsoon and SST show that difficulties in predicting the South Asian monsoon can be mainly attributed to the limited regional response of monsoon in observations but the extensive and exaggerated response in predictions due partially to the application of ensemble average forecasting methods. In contrast, in spite of a similar deficiency, the Southeast Asian monsoon can still be forecasted reasonably, probably because of its closer relationship with large-scale circulation patterns and El Ni˜no–Southern Oscillation.

This paper provides a systematic evaluation of the ability of 12 Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) to simulate the spatial inhomogeneity of the atmospheric carbon dioxide (CO 2 ) concentration. The multi-model ensemble mean (MME) can reasonably simulate the increasing trend of CO 2 concentration from 1850 to 2014, compared with the observation data from the Scripps CO 2 Program and CMIP6 prescribed data, and improves upon the CMIP5 MME CO 2 concentration (which is overestimated after 1950). The growth rate of CO 2 concentration in the northern hemisphere (NH) is higher than that in the southern hemisphere (SH), with the highest growth rate in the mid-latitudes of the NH. The MME can also reasonably simulate the seasonal amplitude of CO 2 concentration, which is larger in the NH than in the SH and grows in amplitude after the 1950s (especially in the NH). Although the results of the MME are reasonable, there is a large spread among ESMs, and the difference between the ESMs increases with time. The MME results show that regions with relatively large CO 2 concentrations (such as northern Russia, eastern China, Southeast Asia, the eastern United States, northern South America, and southern Africa) have greater seasonal variability and also exhibit a larger inter-model spread. Compared with CMIP5, the CMIP6 MME simulates an average spatial distribution of CO 2 concentration that is much closer to the site observations, but the CMIP6-inter-model spread is larger. The inter-model differences of the annual means and seasonal cycles of atmospheric CO 2 concentration are both attributed to the differences in natural sources and sinks of CO 2 between the simulations.

This study isolates the role of the 2019 Southern Hemisphere (SH) sudden stratospheric warming (SSW) for surface climate in eight Stratospheric Nudging And Predictable Surface Impacts project (SNAPSI) models. The novel nudging experiments allow for a clearer disentanglement of the SSW’s impact on surface climate than the standard free-running forecasts, in which the entire atmospheric system evolves freely and deviates from observed reality. SNAPSI models capture the downward propagation of the negative Southern Annular Mode (SAM) from the stratosphere when the zonally symmetric stratospheric circulation is nudged toward the observations. In addition, zonally asymmetric stratospheric variations amplify warm anomalies and dry conditions over eastern Australia, enhancing the Australian local high-pressure system. The composite analysis reveals that the zonally asymmetric stratospheric variations driven by the westward-shifted vortex influence the eastern Australian precipitation forecast by strengthening the downward coupling of the SAM and enhancing the meridional ridge-trough structure over Australia. Westward-shifted vortex events are identified in over 60% more ensemble members in the nudged run than in the free run, associated with precipitation closer to observed values in the nudged run. The wildfire weather potential risk in Australia is assessed using the Fraction Attributable Risk (FAR) of the Hot Dry Windy (HDW) wildfire weather index. The positive FAR of the HDW wildfire weather index indicates that the nudged stratosphere contributes up to ~30% to the increased wildfire weather risk along eastern and southern Australia, highlighting the contribution of the 2019 SSW event to Australian extreme weather from mid-October to mid-November.