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.

Poor air quality is currently responsible for large impacts on human health across the world. In addition, the air pollutants ozone (O₃) and particulate matter less than 2.5 µm in diameter (PM_2.5) are also radiatively active in the atmosphere and can influence Earth's climate. It is important to understand the effect of air quality and climate mitigation measures over the historical period and in different future scenarios to ascertain any impacts from air pollutants on both climate and human health. The Coupled Model Intercomparison Project Phase 6 (CMIP6) presents an opportunity to analyse the change in air pollutants simulated by the current generation of climate and Earth system models that include a representation of chemistry and aerosols (particulate matter). The shared socio-economic pathways (SSPs) used within CMIP6 encompass a wide range of trajectories in precursor emissions and climate change, allowing for an improved analysis of future changes to air pollutants. Firstly, we conduct an evaluation of the available CMIP6 models against surface observations of O₃ and PM_2.5. CMIP6 models consistently overestimate observed surface O₃ concentrations across most regions and in most seasons by up to 16 ppb, with a large diversity in simulated values over Northern Hemisphere continental regions. Conversely, observed surface PM_2.5 concentrations are consistently underestimated in CMIP6 models by up to 10 µg m⁻³, particularly for the Northern Hemisphere winter months, with the largest model diversity near natural emission source regions. The biases in CMIP6 models when compared to observations of O₃ and PM_2.5 are similar to those found in previous studies. Over the historical period (1850–2014) large increases in both surface O₃ and PM_2.5 are simulated by the CMIP6 models across all regions, particularly over the mid to late 20th century, when anthropogenic emissions increase markedly. Large regional historical changes are simulated for both pollutants across East and South Asia with an annual mean increase of up to 40 ppb for O₃ and 12 µg m⁻³ for PM_2.5. In future scenarios containing strong air quality and climate mitigation measures (ssp126), annual mean concentrations of air pollutants are substantially reduced across all regions by up to 15 ppb for O₃ and 12 µg m⁻³ for PM_2.5. However, for scenarios that encompass weak action on mitigating climate and reducing air pollutant emissions (ssp370), annual mean increases in both surface O₃ (up 10 ppb) and PM_2.5 (up to 8 µg m⁻³) are simulated across most regions, although, for regions like North America and Europe small reductions in PM_2.5 are simulated due to the regional reduction in precursor emissions in this scenario. A comparison of simulated regional changes in both surface O₃ and PM_2.5 from individual CMIP6 models highlights important regional differences due to the simulated interaction of aerosols, chemistry, climate and natural emission sources within models. The projection of regional air pollutant concentrations from the latest climate and Earth system models used within CMIP6 shows that the particular future trajectory of climate and air quality mitigation measures could have important consequences for regional air quality, human health and near-term climate. Differences between individual models emphasise the importance of understanding how future Earth system feedbacks influence natural emission sources, e.g. response of biogenic emissions under climate change.

This paper quantifies the pre-industrial (1850) to present-day (2014) effective radiative forcing (ERF) of anthropogenic emissions of NO_X, volatile organic compounds (VOCs; including CO), SO₂, NH₃, black carbon, organic carbon, and concentrations of methane, N₂O and ozone-depleting halocarbons, using CMIP6 models. Concentration and emission changes of reactive species can cause multiple changes in the composition of radiatively active species: tropospheric ozone, stratospheric ozone, stratospheric water vapour, secondary inorganic and organic aerosol, and methane. Where possible we break down the ERFs from each emitted species into the contributions from the composition changes. The ERFs are calculated for each of the models that participated in the AerChemMIP experiments as part of the CMIP6 project, where the relevant model output was available. The 1850 to 2014 multi-model mean ERFs (± standard deviations) are −1.03 ± 0.37 W m⁻² for SO₂ emissions, −0.25 ± 0.09 W m⁻² for organic carbon (OC), 0.15 ± 0.17 W m⁻² for black carbon (BC) and −0.07 ± 0.01 W m⁻² for NH₃. For the combined aerosols (in the piClim-aer experiment) it is −1.01 ± 0.25 W m⁻². The multi-model means for the reactive well-mixed greenhouse gases (including any effects on ozone and aerosol chemistry) are 0.67 ± 0.17 W m⁻² for methane (CH₄), 0.26 ± 0.07 W m⁻² for nitrous oxide (N₂O) and 0.12 ± 0.2 W m⁻² for ozone-depleting halocarbons (HC). Emissions of the ozone precursors nitrogen oxides (NO_x), volatile organic compounds and both together (O₃) lead to ERFs of 0.14 ± 0.13, 0.09 ± 0.14 and 0.20 ± 0.07 W m⁻² respectively. The differences in ERFs calculated for the different models reflect differences in the complexity of their aerosol and chemistry schemes, especially in the case of methane where tropospheric chemistry captures increased forcing from ozone production.

Quantification and attribution of long-term tropospheric ozone trends are critical for understanding the impact of human activity and climate change on atmospheric chemistry but are also challenged by the limited coverage of long-term ozone observations in the free troposphere where ozone has higher production efficiency and radiative potential compared to that at the surface. In this study, we examine observed tropospheric ozone trends, their attributions, and radiative impacts from 1995–2017 using aircraft observations from the In-service Aircraft for a Global Observing System database (IAGOS), ozonesondes, and a multi-decadal GEOS-Chem chemical model simulation. IAGOS observations above 11 regions in the Northern Hemisphere and 19 of 27 global ozonesonde sites have measured increases in tropospheric ozone (950–250 hPa) by 2.7 ± 1.7 and 1.9 ± 1.7 ppbv per decade on average, respectively, with particularly large increases in the lower troposphere (950–800 hPa) above East Asia, the Persian Gulf, India, northern South America, the Gulf of Guinea, and Malaysia/Indonesia by 2.8 to 10.6 ppbv per decade. The GEOS-Chem simulation driven by reanalysis meteorological fields and the most up-to-date year-specific anthropogenic emission inventory reproduces the overall pattern of observed tropospheric ozone trends, including the large ozone increases over the tropics of 2.1–2.9 ppbv per decade and above East Asia of 0.5–1.8 ppbv per decade and the weak tropospheric ozone trends above North America, Europe, and high latitudes in both hemispheres, but trends are underestimated compared to observations. GEOS-Chem estimates an increasing trend of 0.4 Tg yr⁻¹ of the tropospheric ozone burden in 1995–2017. We suggest that uncertainties in the anthropogenic emission inventory in the early years of the simulation (e.g., 1995–1999) over developing regions may contribute to GEOS-Chem's underestimation of tropospheric ozone trends. GEOS-Chem sensitivity simulations show that changes in global anthropogenic emission patterns, including the equatorward redistribution of surface emissions and the rapid increases in aircraft emissions, are the dominant factors contributing to tropospheric ozone trends by 0.5 Tg yr⁻¹. In particular, we highlight the disproportionately large, but previously underappreciated, contribution of aircraft emissions to tropospheric ozone trends by 0.3 Tg yr⁻¹, mainly due to aircraft emitting NO_x in the mid-troposphere and upper troposphere where ozone production efficiency is high. Decreases in lower-stratospheric ozone and the stratosphere–troposphere flux in 1995–2017 contribute to an ozone decrease at mid-latitudes and high latitudes. We estimate the change in tropospheric ozone radiative impacts from 1995–1999 to 2013–2017 is +18.5 mW m⁻², with 43.5 mW m⁻² contributed by anthropogenic emission changes (20.5 mW m⁻² alone by aircraft emissions), highlighting that the equatorward redistribution of emissions to areas with strong convection and the increase in aircraft emissions are effective for increasing tropospheric ozone's greenhouse effect.

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.

BCC-CSM2-HR is a high-resolution version of the Beijing Climate Center (BCC) Climate System Model (T266 in the atmosphere and $\mathrm{1}/\mathrm{4}$ ^∘ latitude × $\mathrm{1}/\mathrm{4}$ ^∘ longitude in the ocean). Its development is on the basis of the medium-resolution version BCC-CSM2-MR (T106 in the atmosphere and 1^∘ latitude × 1^∘ longitude in the ocean) which is the baseline for BCC participation in the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study documents the high-resolution model, highlights major improvements in the representation of atmospheric dynamical core and physical processes. BCC-CSM2-HR is evaluated for historical climate simulations from 1950 to 2014, performed under CMIP6-prescribed historical forcing, in comparison with its previous medium-resolution version BCC-CSM2-MR. Observed global warming trends of surface air temperature from 1950 to 2014 are well captured by both BCC-CSM2-MR and BCC-CSM2-HR. Present-day basic atmospheric mean states during the period from 1995 to 2014 are then evaluated at global scale, followed by an assessment on climate variabilities in the tropics including the tropical cyclones (TCs), the El Niño–Southern Oscillation (ENSO), the Madden–Julian Oscillation (MJO), and the quasi-biennial oscillation (QBO) in the stratosphere. It is shown that BCC-CSM2-HR represents the global energy balance well and can realistically reproduce the main patterns of atmospheric temperature and wind, precipitation, land surface air temperature, and sea surface temperature (SST). It also improves the spatial patterns of sea ice and associated seasonal variations in both hemispheres. The bias of the double intertropical convergence zone (ITCZ), obvious in BCC-CSM2-MR, almost disappears in BCC-CSM2-HR. TC activity in the tropics is increased with resolution enhanced. The cycle of ENSO, the eastward propagative feature and convection intensity of MJO, and the downward propagation of QBO in BCC-CSM2-HR are all in a better agreement with observations than their counterparts in BCC-CSM2-MR. Some imperfections are, however, noted in BCC-CSM2-HR, such as the excessive cloudiness in the eastern basin of the tropical Pacific with cold SST biases and the insufficient number of tropical cyclones in the North Atlantic.