Search

The control of the production of ozone-depleting substances through the Montreal Protocol means that the stratospheric ozone layer is recovering.sup.1 and that consequent increases in harmful surface ultraviolet radiation are being avoided.sup.2,3. The Montreal Protocol has co-benefits for climate change mitigation, because ozone-depleting substances are potent greenhouse gases.sup.4-7. The avoided ultraviolet radiation and climate change also have co-benefits for plants and their capacity to store carbon through photosynthesis.sup.8, but this has not previously been investigated. Here, using a modelling framework that couples ozone depletion, climate change, damage to plants by ultraviolet radiation and the carbon cycle, we explore the benefits of avoided increases in ultraviolet radiation and changes in climate on the terrestrial biosphere and its capacity as a carbon sink. Considering a range of strengths for the effect of ultraviolet radiation on plant growth.sup.8-12, we estimate that there could have been 325-690 billion tonnes less carbon held in plants and soils by the end of this century (2080-2099) without the Montreal Protocol (as compared to climate projections with controls on ozone-depleting substances). This change could have resulted in an additional 115-235 parts per million of atmospheric carbon dioxide, which might have led to additional warming of global-mean surface temperature by 0.50-1.0 degrees. Our findings suggest that the Montreal Protocol may also be helping to mitigate climate change through avoided decreases in the land carbon sink. Modelling suggests that the Montreal Protocol may be mitigating climate change by protecting the land carbon sink, as well as by protecting the ozone layer and reducing greenhouse gas emissions., Author(s): Paul J. Young [sup.1] [sup.2] [sup.3] , Anna B. Harper [sup.4] [sup.5] , Chris Huntingford [sup.6] , Nigel D. Paul [sup.1] [sup.7] , Olaf Morgenstern [sup.8] , Paul A. [...]

An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH3CCl3) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom-up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long-lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH-SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long-term trend and emissions derived from the measured hemispheric gradient, the combination of HFC-32 (CH2F2), HFC-134a (CH2FCF3, HFC-152a (CH3CHF2), and HCFC-22 (CHClF2), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF.

The most recent generation of climate models that has informed the Sixth Assessment Report (AR6) of the Intergovernmental Panel on Climate Change (IPCC) is characterized by the presence of several models with larger equilibrium climate sensitivities (ECSs) and transient climate responses (TCRs) than exhibited by the previous generation. Partly as a result, AR6 did not use any direct quantifications of ECSs and TCRs based on the 4 × CO2 and 1pctCO2 simulations and relied on other evidence when assessing the Earth's actual ECS and TCR. Here I use historical observed global-mean temperature and simulations produced under the Detection and Attribution Model Intercomparison Project to constrain the ECS, TCR, and historical aerosol-related cooling. I introduce additivity criteria that disqualify 8 of the participating 16 models from consideration in multi-model averaging calculations. Based on the remaining eight models, I obtain an average adjusted ECS of 3.5 ± 0.4 K and a TCR of 1.8 ± 0.3 K (both at 68 % confidence). Both are consistent with the AR6 estimates but with substantially reduced uncertainties. Furthermore, importantly I find that the optimal cooling due to short-lived climate forcers consistent with the observed temperature record should, on average, be about 47 % ± 39 % of what these models simulate in their aerosol-only simulations, yielding a multi-model mean, global-mean, and annual-mean cooling due to near-term climate forcers for 2000–2014, relative to 1850–1899, of 0.24 ± 0.11 K (at 68 % confidence). This is consistent with but at the lower end of the very likely uncertainty range of the IPCC's AR6. There is a correlation between the models' ECSs and their aerosol-related cooling, whereby large-ECS models tend to be associated also with strong aerosol-related cooling. The results imply that a reduction in the aerosol-related cooling, along with a more moderate adjustment of the greenhouse-gas-related warming for most models, would bring the historical global-mean temperature simulated by these models into better agreement with observations. [ABSTRACT FROM AUTHOR]

This study presents an updated and homogenised ozone time series covering 34 years (1987–2020) of ozonesonde measurements at Lauder, New Zealand, and attributes vertically resolved ozone trends using a multiple linear regression (MLR) analysis and a chemistry–climate model (CCM). Homogenisation of the time series leads to a marked difference in ozone values before 1997, in which the ozone trends are predominantly negative from the surface to ∼ 30 km, ranging from ∼ - 2 % per decade to - 13 % per decade, maximising at around 12–13 km, in contrast to the uncorrected time series which shows no clear trends for this period. For the post-2000 period, ozone at Lauder shows negative trends in the stratosphere, maximising just below 20 km (∼ - 5 % per decade) despite the fact that stratospheric chlorine and bromine from ozone-depleting substances (ODSs) have both been declining since 1997. However, the ozone trends change from negative for 1987–1999 to positive in the post-2000 period in the free troposphere. The post-2000 ozone trends calculated from the ozonesonde measurements compare well with those derived from the co-located low-vertical-resolution Fourier-transform infrared spectroscopy (FTIR) ozone time series. The MLR analysis identifies that the increasing tropopause height, associated with CO2 -driven dynamical changes, is the leading factor driving the continuous negative trend in lower-stratospheric ozone at Lauder over the whole observational period, whilst the ozone-depleting substances (ODSs) only contribute to the negative ozone trend in the lower stratosphere over the pre-1999 period. Meanwhile, stratospheric temperature changes contribute significantly to the negative ozone trend above 20 km over the post-2000 period. Furthermore, the chemistry–climate model (CCM) simulations that separate the effects of individual forcings show that the predominantly negative modelled trend in ozone for the 1987–1999 period is driven not only by ODSs but also by increases in greenhouse gases (GHGs), with large but opposing impacts from methane (positive) and CO2 (negative), respectively. Over the 2000–2020 period, the model results show that the CO2 increase is the dominant driver for the negative trend in the lower stratosphere, in agreement with the MLR analysis. Although the model underestimates the observed negative ozone trend in the lower stratosphere for both periods, it clearly shows that CO2 -driven dynamical changes have played an increasingly important role in driving the lower-stratospheric ozone trends in the vicinity of Lauder. [ABSTRACT FROM AUTHOR]

Large uncertainty exists in the sign of long‐term changes in regional scale mean precipitation across the current generation of global climate models. To explore the physical drivers of this uncertainty for New Zealand, here we adopt a storyline approach applying cluster analysis to spatial patterns of future projected seasonal mean precipitation change across CMIP6 models (n = 43). For the winter precipitation change signal, the models split roughly into two main groups: both groups have a very robust wet signal across the west coast of the South Island but differ notably in terms of the sign of precipitation change across the north of the North Island. These far north winter precipitation differences appear related to how far the Hadley cell edge and regional eddy‐driven jet shift across the models relative to their historical positions. In contrast, for summer, most models have a markedly weaker and spatially non‐uniform response, where internal variability often plays a large role. However, a small group of models predict a robust wet signal across most of the country in summer. This "wet model" group is characterized by a regional La Niña‐like increase in high pressure shifted further to the south‐east of New Zealand, associated with more frequent north‐easterly flow over the country and accompanied by significant warming of local sea surface temperatures. This regional circulation response appears related to changes in stationary Rossby wave paths as opposed to changes in La Niña occurrence frequency itself. Plain Language Summary: At the regional scale, future changes to mean precipitation under climate change could carry large societal consequences. Unfortunately, large uncertainties still exist on regional scales which may hinder climate change adaptation efforts. Here we explore and characterize these uncertainties across the latest generation of global climate models for the New Zealand region. Across the models, winter precipitation changes are shown to be much more consistent compared to summer precipitation changes. In winter, changes in the jet stream and Hadley cell edge positions in the models are important for determining the regional spatial patterns of precipitation change. In summer, internal variability uncertainty plays a larger role, models that predict robust wet changes across the country are associated with more north‐easterly flow conditions in the future period. Changes to Rossby wave pathways appear important for setting up this regional circulation response in summer. Key Points: Storylines are used to characterize and explain the main precipitation change patterns across modelsSpatial patterns of precipitation change are more robust in winter, inter‐model differences relate to Hadley cell and jet changesSpatial patterns of precipitation change are less robust in summer, internal variability and Rossby wave pathway changes are important [ABSTRACT FROM AUTHOR]

In this work we analyse the impact of including a regional, high-resolution ocean model on simulated atmospheric climate in a coupled earth system model. The resolution of the regional, nested ocean model is approximately 0.2° compared to the ~1° resolution of the global ocean model within which it is embedded and this work complements previously published work on ocean circulation and marine heatwaves using this setup, referred to as the New Zealand Earth System Model, NZESM. After a brief discussion of the wider model setup, the persistent Southern Ocean warm bias in climate models and the validation data sets used, we show the effects of the altered ocean physics on air temperature, precipitation and evaporation, latent and sensible surface heat balances, westerly winds, the storm track and the effect on total cloud amount. Overall we find that the NZESM provides a better representation of regional atmospheric climate compared to its parent model - UKESM1 - although this improvement is not universal. For example, although the NZESM shows better agreement in surface air temperature within the nested ocean region, there is also some deterioration in the agreement at higher southern latitudes where the seasonal sea ice edge coincides with a transition from negative to positive correlation between air temperature and cloud amount. The lack of additional model tuning in the NZESM after the nested ocean model's inclusion largely accounts for the presence of these improvement-deterioration pairs with respect to observations. The reader is encouraged to read the paper of Behrens et al. (Behrens et al, 2020) before this one since it provides much additional information which will aid understanding. This study aims to provide a high-level reference ontology for how changing one aspect of the ocean physics in a coupled model can impact simulated atmospheric climate. [ABSTRACT FROM AUTHOR]

The most recent generation of climate models that has informed the 6th Assessment Report (AR6) of IPCC is characterized by the presence of several models with anomalously large equilibrium climate sensitivities (ECSs) relative to the previous generation. Partly as a result, AR6 did not use any direct quantifications of ECSs using 4×CO2 simulations and relied on other evidence when assessing the Earth's actual ECS. Here I use the historical observed global-mean surface air temperature and simulations produced under the Detection and AttributionModel Intercomparison Project to constrain the ECS and historical aerosol-related cooling. Based on 15 largely independent models I obtain an average adjusted ECS of 3.4±0.8 K (at 68% confidence), which is very consistent with the AR6 estimate. Furthermore, importantly I find that the optimal cooling due to anthropogenic aerosols consistent with the observed temperature record should on average be about 34±31% of what these models simulate, yielding an aerosol-related global-mean cooling for 2000-2014, relative to 1850-1899, of -0.19±0.14 K (at 68% confidence), when these models simulate on average -0.63±0.28 K. For 12 models the reduction in aerosol-related cooling equals or exceeds 50%. There is a correlation between the models' ECS and their aerosol-related cooling, whereby large-ECS models tend to be associated also with large aerosol-related cooling. The results imply that a large reduction of the aerosol-related cooling, along with a more moderate adjustment of the greenhouse-gas related warming, for most models would bring the historical global mean temperature simulated by these models into better agreement with observations. [ABSTRACT FROM AUTHOR]

This study presents an updated and homogenised ozone time series covering 34 years (1987-2020) of ozonesonde measurements at Lauder, New Zealand, and derived vertically resolved ozone trends. Over the period of 1987-1999, the ozone trends in the homogenised ozone data are predominantly negative from the surface to ~30 km, ranging from ~ -2 to -12% decade-1, maximising at around 12-13 km. These negative trends are statistically significant at 95% confidence level below 5 km and above 17 km. For the post-2000 period, ozone at Lauder shows negative trends in the stratosphere (but the trends are only statistically significant above 17 km), maximising just below 20 km (~ -5% decade-1), despite stratospheric chlorine and bromine from ozone-depleting substances (ODSs) both declining in this period. In the troposphere, the ozone trends change from negative for 1987-1999 to positive in the post-2000 period. The post-2000 ozone trends from the ozonesonde measurements compare well with those from a low-vertical resolution Fourier-transform infrared spectroscopy (FTIR) ozone time series. A multiple-linear regression analysis indicates that anthropogenic forcing plays a significant role in driving the significant negative trend in the stratospheric ozone at Lauder, in which the effect of greenhouse gas (GHG)-driven dynamical and chemical changes is reflected in the significant positive trends in tropopause height and tropospheric temperature, and significant negative trends of stratospheric temperature observed at Lauder. The interannual variation in lower stratospheric ozone is largely explained by the variation in tropopause height at Lauder, which is highly anti-correlated with stratospheric temperature and correlated with tropospheric temperature. Furthermore, the impact of ODSs and GHGs on ozone over Lauder is assessed in a chemistry-climate model using a series of single forcing simulations. The model simulations show that the predominantly negative modelled trend in ozone for the 1987-1999 period is driven not only by ODSs, but also by increases in GHGs with large but opposing impacts from methane (positive) and CO2 (negative), respectively. Over the 2000-2020 period, although the model underestimates the observed negative ozone trend in the lower stratosphere but clearly shows that CO2-driven dynamical changes have had an increasingly important role in driving ozone trends in this region. [ABSTRACT FROM AUTHOR]

This study presents an updated and homogenised ozone time series covering 34 years (1987–2020) of ozonesonde measurements at Lauder, New Zealand, and derived vertically resolved ozone trends. Over the period of 1987–1999, the ozone trends in the homogenised ozone data are predominantly negative from the surface to ∼30 km, ranging from −2 to −10 % decade−1, maximising at around 12–13 km. These negative trends are statistically significant at 95 % confidence level below 5 km and above 17 km. For the post-2000 period, ozone at Lauder shows negative trends in the stratosphere (but the trends are only statistically significant above 17 km), maximising just below 20 km (∼ −5 % decade−1), despite stratospheric chlorine and bromine from ozone-depleting substances (ODSs) both declining in this period. In the troposphere, the ozone trends change from negative for 1987–1999 to positive in the post-2000 period. The post-2000 ozone trends from the ozonesonde measurements compare well with those from a low-vertical resolution Fourier-transform infrared spectroscopy (FTIR) ozone time series. A multiple-linear regression analysis indicates that anthropogenic forcing plays a significant role in driving the significant negative trend in the stratospheric ozone at Lauder, in which the effect of greenhouse gas (GHG)-driven dynamical and chemical changes is reflected in the significant positive trends in tropopause height and tropospheric temperature, and significant negative trends of stratospheric temperature observed at Lauder. The interannual variation in lower stratospheric ozone is largely explained by the variation in tropopause height at Lauder, which is highly anti-correlated with stratospheric temperature and correlated with tropospheric temperature. Furthermore, the impact of ODSs and GHGs on ozone over Lauder is assessed in a chemistry-climate model using a series of single forcing simulations. The model simulations show that the predominantly negative modelled trend in ozone for the 1987–1999 period is driven not only by ODSs, but also by increases in GHGs with large but opposing impacts from methane (positive) and CO2 (negative), respectively. Over the 2000–2020 period, although the model underestimates the observed negative ozone trend in the lower stratosphere but clearly shows that CO2-driven dynamical changes have had an increasingly important role in driving ozone trends in this region. [ABSTRACT FROM AUTHOR]

In the Arctic stratosphere, the combination of chemical ozone depletion by halogenated ozone-depleting substances (hODSs) and dynamic fluctuations can lead to severe ozone minima. These Arctic ozone minima are of great societal concern due to their health and climate impacts. Owing to the success of the Montreal Protocol, hODSs in the stratosphere are gradually declining, resulting in a recovery of the ozone layer. On the other hand, continued greenhouse gas (GHG) emissions cool the stratosphere, possibly enhancing the formation of polar stratospheric clouds (PSCs) and, thus, enabling more efficient chemical ozone destruction. Other processes, such as the acceleration of the Brewer–Dobson circulation, also affect stratospheric temperatures, further complicating the picture. Therefore, it is currently unclear whether major Arctic ozone minima will still occur at the end of the 21st century despite decreasing hODSs. We have examined this question for different emission pathways using simulations conducted within the Chemistry-Climate Model Initiative (CCMI-1 and CCMI-2022) and found large differences in the models' ability to simulate the magnitude of ozone minima in the present-day climate. Models with a generally too-cold polar stratosphere (cold bias) produce pronounced ozone minima under present-day climate conditions because they simulate more PSCs and, thus, high concentrations of active chlorine species (ClOx). These models predict the largest decrease in ozone minima in the future. Conversely, models with a warm polar stratosphere (warm bias) have the smallest sensitivity of ozone minima to future changes in hODS and GHG concentrations. As a result, the scatter among models in terms of the magnitude of Arctic spring ozone minima will decrease in the future. Overall, these results suggest that Arctic ozone minima will become weaker over the next decades, largely due to the decline in hODS abundances. We note that none of the models analysed here project a notable increase of ozone minima in the future. Stratospheric cooling caused by increasing GHG concentrations is expected to play a secondary role as its effect in the Arctic stratosphere is weakened by opposing radiative and dynamical mechanisms. [ABSTRACT FROM AUTHOR]

Tropical cyclones – TCs – affecting the South Pacific region are studied using coupled atmosphere-ocean earth system models and (offline) storm tracking software which tracks the position of simulated pressure lows through time. The models used are the United Kingdom Earth System Model, version 1 – UKESM1 – and the related New Zealand Earth System Model, the NZESM. The model pair considered here differ only in their treatment of the ocean and the NZESM has a nominal resolution of 0.2° in the region surrounding New Zealand and 1° elsewhere; UKESM1 has a a uniform 1° resolution everywhere. After validating the storm tracking algorithm against the track of cyclone Giselle from 1968 and cyclone Gabrielle from 2023 we use the Saffir-Simpson scale to split the tracked systems into categories based on their severity. For systems formed in the vicinity of New Zealand (and globally) the overall number is overestimated but stronger (category 2 and 3) storms are underestimated. We also see a general decrease in the total number of storms as radiative forcing, F , increases although there is some evidence of a small increase at extreme levels of warming. In the metrics studied here we find no difference between the ensembles of UKESM1 and NZESM simulations and going forward use the UKESM1, which has larger available ensembles. The power dissipation index, PDI, gives a first order measure of TC strength and we find that the average PDI per storm increases with F by up to 26 % under a 'fossil-fuelled development' scenario. Although the physical mechanisms behind the increase in average PDI with F are relatively simple to understand, those governing the frequency of occurrence are not. In the results shown here, vertical wind shear increases with F which tends to reduce TC numbers but the effect of the tropospheric relative humidity is much less clear. The increase in the area of the tropics bounded by the 26.5° isotherm should, on its own, increase the number of TCs, in opposition to the general behaviour observed, except perhaps at extreme levels of future warming. [ABSTRACT FROM AUTHOR]

Using nine chemistry‐climate and eight associated no‐chemistry models, we investigate the persistence and timing of cold episodes occurring in the Arctic and Antarctic stratosphere during the period 1980–2014. We find systematic differences in behavior between members of these model pairs. In a first group of chemistry models whose dynamical configurations mirror their no‐chemistry counterparts, we find an increased persistence of such cold polar vortices, such that these cold episodes often start earlier and last longer, relative to the times of occurrence of the lowest temperatures. Also the date of occurrence of the lowest temperatures, both in the Arctic and the Antarctic, is often delayed by 1–3 weeks in chemistry models, versus their no‐chemistry counterparts. This behavior exacerbates a widespread problem occurring in most or all models, a delayed occurrence, in the median, of the most anomalously cold day during such cold winters. In a second group of model pairs there are differences beyond just ozone chemistry. In particular, here the chemistry models feature more levels in the stratosphere, a raised model top, and differences in non‐orographic gravity wave drag versus their no‐chemistry counterparts. Such additional dynamical differences can completely mask the above influence of ozone chemistry. The results point toward a need to retune chemistry‐climate models versus their no‐chemistry counterparts. Plain Language Summary: Ozone is a chemical constituent of the atmosphere acting as an absorber of both solar ultraviolet light and infrared radiation emitted by the Earth. It therefore needs to be considered in climate models. Explicit ozone chemistry is a computationally challenging addition to a climate model; hence in most cases ozone is simply prescribed. Especially during relatively cold stratospheric winter/spring seasons, Antarctic and Arctic ozone depletion can be considerable. Such anomalous ozone loss is not reflected in the imposed ozone field, and hence differences in behavior are expected for such situations between chemistry‐ and no‐chemistry models. Indeed for such cold winters/springs, we find an enhanced persistence of such cold spells in a set of chemistry‐climate models, versus their no‐chemistry counterparts; such enhanced persistence generally makes the chemistry model less realistic than its no‐chemistry counterpart. However, if there are substantial further differences between the members of these model pairs, such as regarding their grid configuration or physical processes beyond chemistry, these can obscure the effect of ozone chemistry. We thus claim that adding stratospheric ozone chemistry to a climate model necessitates retuning to counteract a deterioration of the simulated stratospheric climate that can otherwise occur. Key Points: Coupling in ozone chemistry causes an increase in persistence of low temperature anomalies over both polesIn the Antarctic, coupling in chemistry amplifies pre‐existing stratospheric cold biasesThese effects can be masked by other dynamical differences present in some models [ABSTRACT FROM AUTHOR]

We quantify the impacts of halogenated ozone‐depleting substances (ODSs), greenhouse gases (GHGs), and short‐lived ozone precursors on ozone changes between 1850 and 2014 using single‐forcing perturbation simulations from several Earth system models with interactive chemistry participating in the Coupled Model Intercomparison Project Aerosol and Chemistry Model Intercomparison Project. We present the responses of ozone to individual forcings and an attribution of changes in ozone columns and vertically resolved stratospheric and tropospheric ozone to these forcings. We find that whilst substantial ODS‐induced ozone loss dominates the stratospheric ozone changes since the 1970s, in agreement with previous studies, increases in tropospheric ozone due to increases in short‐lived ozone precursors and methane since the 1950s make increasingly important contributions to total column ozone (TCO) changes. Increases in methane also lead to substantial extra‐tropical stratospheric ozone increases. Impacts of nitrous oxide and carbon dioxide on stratospheric ozone are significant but their impacts on TCO are small overall due to several opposing factors and are also associated with large dynamical variability. The multi‐model mean (MMM) results show a clear change in the stratospheric ozone trends after 2000 due to now declining ODSs, but the trends are generally not significantly positive, except in the extra‐tropical upper stratosphere, due to relatively small changes in forcing over this period combined with large model uncertainty. Although the MMM ozone compares well with the observations, the inter‐model differences are large primarily due to the large differences in the models' representation of ODS‐induced ozone depletion. Plain Language Summary: Overhead ozone absorbs harmful solar ultraviolet light, protecting life on Earth. Due to human activities since the nineteenth century, emissions of greenhouse gases (GHGs) and ozone‐depleting substances (ODSs) containing chlorine and bromine have profoundly affected stratospheric ozone. Near the Earth' surface, ozone has increased substantially leading to worsening air quality. In this study, we use Earth system models to interactively assess the roles of ODSs, ozone‐forming pollutants, and GHGs including methane, carbon dioxide (CO2), and nitrous oxide (N2O) on ozone changes from the surface to the upper stratosphere. Whilst substantial reductions in stratospheric ozone due to ODSs occurred since the 1970s, the lower‐atmospheric ozone increases due to anthropogenic pollution have counteracted this decrease. Increases in GHGs lead to various positive and negative effects on stratospheric ozone in different regions, and their impacts vary with ODS levels in the atmosphere. Amongst the GHGs assessed here, the increase in methane leads to overwhelming positive trends in both stratospheric and tropospheric ozone through mainly chemical effects. The impact of changes in N2O and CO2 on total column ozone is more uncertain due to large inter‐model differences, although their overall impact is small during the simulation period. Key Points: New multi‐model results show significant positive effects of ozone precursors on near‐global ozone offsetting the negative effects of ozone‐depleting substances (ODSs)ODS and greenhouse gases dominate stratospheric ozone changes but with large inter‐model differences due to uncertainties in responses to ODS changesIncreases in carbon dioxide and nitrous oxide significantly impact stratospheric ozone, but their net effects on total columns are small due to cancellations [ABSTRACT FROM AUTHOR]

We evaluate total‐column ozone trends using a piecewise linear regression approach and maximizing usage of five gridded total‐column ozone data sets. The new approach yields more consistent estimates of observed ozone loss during 1979–2000, when halocarbon concentrations were increasing, and consequently, using CMIP6 simulations, an increased effective radiative forcing estimate of ozone‐depleting substances with a substantially reduced uncertainty range versus an earlier evaluation. At more than 84% confidence, it is now larger than zero and compares more favorably with three previous evaluations. We furthermore find significant positive post‐1997 global‐ and Southern‐Hemisphere‐mean trends, respectively, in these five data sets. For the extrapolar region (60°S–60°N) and for the Northern Hemisphere, the assessment whether there is a positive trend is inconclusive and depends on which observational data set is included in the calculation. Plain Language Summary: Changes in overhead ozone amounts reflect the impact of the Montreal Protocol, designed to protect the ozone layer, and several other influences. Here, we assess five different ozone data sets using satellite and ground‐based observations as well as fields generated by present‐generation climate models. For the period 1979–2000, during which stratospheric ozone depletion got established, we find good agreement for the whole globe and for selected subregions across the observational data sets. For 1997–2020, in the global and Southern Hemisphere means, we find significant, positive ozone trends. For a region excluding both poles and for the Northern Hemisphere, the uncertainty ranges still include zero. Using observational and modeled ozone trends for 1979–2000, we recalculate the impact of ozone‐depleting substances, accounting for ozone depletion, on the Earth' radiation balance. We find a slightly larger net impact than a previous evaluation, which within the uncertainty bounds is much more likely to be positive and is more consistent with three other literature references. Key Points: We evaluate total‐column ozone trends for 1979–2000 and 1997–2020For 1997–2020, we find significant global‐ and Southern‐Hemisphere‐mean positive trendsThe effective radiative forcing of ozone‐depleting substances is now more consistent with three previous evaluations [ABSTRACT FROM AUTHOR]

Insufficient reflection of short-wave radiation especially over the Southern Ocean region is still a leading issue in many present-day global climate models. One of the potential reasons for this observed bias is an inadequate representation of clouds. In a previous study, we modified the cloud micro-physics scheme in the Unified Model and showed that choosing a more realistic value for the capacitance or shape parameter of atmospheric ice-crystals, in better agreement with theory and observations, benefits the simulation of short-wave radiation over the Southern Ocean by brightening the clouds. However, attempts to modify the cloud phase by directly adjusting the micro-physics process rates like capacitance tend to affect both the hemispheres symmetrically whereas we seek to brighten only the high-latitude Southern Hemisphere clouds. In this study we implement a simple prognostic parametrisation whereby the heterogeneous ice nucleation temperature is made to vary three-dimensionally as a function of the mineral dust distribution in the model. As a result, those regions with less dust number density would have lower nucleation temperature compared to the default global value of -10° C. By using mineral dust as an indicator for ice nucleating particles in the model, this parametrisation thus captures the impact of ice nucleating particles on the cloud distribution due to its general paucity over the Southern Ocean region. This approach thus improves the physics of the model with minimal complexity. [ABSTRACT FROM AUTHOR]

Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼ 300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer–Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions. [ABSTRACT FROM AUTHOR]

The connection between the dominant mode of interannual variability in the tropical troposphere, the El Niño–Southern Oscillation (ENSO), and the entry of stratospheric water vapor is analyzed in a set of model simulations archived for the Chemistry-Climate Model Initiative (CCMI) project and for Phase 6 of the Coupled Model Intercomparison Project. While the models agree on the temperature response to ENSO in the tropical troposphere and lower stratosphere, and all models and observations also agree on the zonal structure of the temperature response in the tropical tropopause layer, the only aspect of the entry water vapor response with consensus in both models and observations is that La Niña leads to moistening in winter relative to neutral ENSO. For El Niño and for other seasons, there are significant differences among the models. For example, some models find that the enhanced water vapor for La Niña in the winter of the event reverses in spring and summer, some models find that this moistening persists, and some show a nonlinear response, with both El Niño and La Niña leading to enhanced water vapor in both winter, spring, and summer. A moistening in the spring following El Niño events, the signal focused on in much previous work, is simulated by only half of the models. Focusing on Central Pacific ENSO vs. East Pacific ENSO, or temperatures in the mid-troposphere compared with temperatures near the surface, does not narrow the inter-model discrepancies. Despite this diversity in response, the temperature response near the cold point can explain the response of water vapor when each model is considered separately. While the observational record is too short to fully constrain the response to ENSO, it is clear that most models suffer from biases in the magnitude of the interannual variability of entry water vapor. This bias could be due to biased cold-point temperatures in some models, but others appear to be missing forcing processes that contribute to observed variability near the cold point. [ABSTRACT FROM AUTHOR]

Quantifying forcings from anthropogenic perturbations to the Earth system (ES) is important for understanding changes in climate since the pre-industrial (PI) period. Here, we quantify and analyse a wide range of present-day (PD) anthropogenic effective radiative forcings (ERFs) with the UK's Earth System Model (ESM), UKESM1, following the protocols defined by the Radiative Forcing Model Intercomparison Project (RFMIP) and the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP). In particular, quantifying ERFs that include rapid adjustments within a full ESM enables the role of various chemistry–aerosol–cloud interactions to be investigated. Global mean ERFs for the PD (year 2014) relative to the PI (year 1850) period for carbon dioxide (CO 2), nitrous oxide (N 2 O), ozone-depleting substances (ODSs), and methane (CH 4) are 1.89 ± 0.04, 0.25 ± 0.04, - 0.18 ± 0.04, and 0.97 ± 0.04 W m -2 , respectively. The total greenhouse gas (GHG) ERF is 2.92 ± 0.04 W m -2. UKESM1 has an aerosol ERF of - 1.09 ± 0.04 W m -2. A relatively strong negative forcing from aerosol–cloud interactions (ACI) and a small negative instantaneous forcing from aerosol–radiation interactions (ARI) from sulfate and organic carbon (OC) are partially offset by a substantial forcing from black carbon (BC) absorption. Internal mixing and chemical interactions imply that neither the forcing from ARI nor ACI is linear, making the aerosol ERF less than the sum of the individual speciated aerosol ERFs. Ozone (O 3) precursor gases consisting of volatile organic compounds (VOCs), carbon monoxide (CO), and nitrogen oxides (NO x), but excluding CH 4 , exert a positive radiative forcing due to increases in O 3. However, they also lead to oxidant changes, which in turn cause an indirect aerosol ERF. The net effect is that the ERF from PD–PI changes in NO x emissions is negligible at 0.03 ± 0.04 W m -2 , while the ERF from changes in VOC and CO emissions is 0.33 ± 0.04 W m -2. Together, aerosol and O 3 precursors (called near-term climate forcers (NTCFs) in the context of AerChemMIP) exert an ERF of - 1.03 ± 0.04 W m -2 , mainly due to changes in the cloud radiative effect (CRE). There is also a negative ERF from land use change (- 0.17 ± 0.04 W m -2). When adjusted from year 1850 to 1700, it is more negative than the range of previous estimates, and is most likely due to too strong an albedo response. In combination, the net anthropogenic ERF (1.76 ± 0.04 W m -2) is consistent with other estimates. By including interactions between GHGs, stratospheric and tropospheric O 3 , aerosols, and clouds, this work demonstrates the importance of ES interactions when quantifying ERFs. It also suggests that rapid adjustments need to include chemical as well as physical adjustments to fully account for complex ES interactions. [ABSTRACT FROM AUTHOR]

Automatic lidars and ceilometers (ALCs) provide valuable information on cloud and aerosols but have not been systematically used in the evaluation of general circulation models (GCMs) and numerical weather prediction (NWP) models. Obstacles associated with the diversity of instruments, a lack of standardisation of data products and open processing tools mean that the value of large ALC networks worldwide is not being realised. We discuss a tool, called the Automatic Lidar and Ceilometer Framework (ALCF), that overcomes these problems and also includes a ground-based lidar simulator, which calculates the radiative transfer of laser radiation and allows one-to-one comparison with models. Our ground-based lidar simulator is based on the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP), which has been extensively used for spaceborne lidar intercomparisons. The ALCF implements all steps needed to transform and calibrate raw ALC data and create simulated attenuated volume backscattering coefficient profiles for one-to-one comparison and complete statistical analysis of clouds. The framework supports multiple common commercial ALCs (Vaisala CL31, CL51, Lufft CHM 15k and Droplet Measurement Technologies MiniMPL), reanalyses (JRA-55, ERA5 and MERRA-2) and models (the Unified Model and AMPS – the Antarctic Mesoscale Prediction System). To demonstrate its capabilities, we present case studies evaluating cloud in the supported reanalyses and models using CL31, CL51, CHM 15k and MiniMPL observations at three sites in New Zealand. We show that the reanalyses and models generally underestimate cloud fraction. If sufficiently high-temporal-resolution model output is available (better than 6-hourly), a direct comparison of individual clouds is also possible. We demonstrate that the ALCF can be used as a generic evaluation tool to examine cloud occurrence and cloud properties in reanalyses, NWP models, and GCMs, potentially utilising the large amounts of ALC data already available. This tool is likely to be particularly useful for the analysis and improvement of low-level cloud simulations which are not well monitored from space. This has previously been identified as a critical deficiency in contemporary models, limiting the accuracy of weather forecasts and future climate projections. While the current focus of the framework is on clouds, support for aerosol in the lidar simulator is planned in the future. [ABSTRACT FROM AUTHOR]

We assess the effective radiative forcing due to ozone‐depleting substances using models participating in the Aerosols and Chemistry and Radiative Forcing Model Intercomparison Projects (AerChemMIP, RFMIP). A large intermodel spread in this globally averaged quantity necessitates an "emergent constraint" approach whereby we link the radiative forcing to ozone declines measured and simulated during 1979–2000, excluding two volcanically perturbed periods. During this period, ozone‐depleting substances were increasing, and several merged satellite‐based climatologies document the ensuing decline of total‐column ozone. Using these analyses, we find an effective radiative forcing of −0.05 to 0.13 W m−2. Our best estimate (0.04 W m−2) is on the edge of the "likely" range given by the Fifth Assessment Report of IPCC of 0.03 to 0.33 W m−2 but is in better agreement with two other literature results. Plain Language Summary: Chloroflourocarbons and other compounds involved in ozone depletion are also powerful greenhouse gases, but their contribution to global warming is reduced due to the cooling effect of the ozone loss which they induce. Models informing an upcoming climate report disagree on the ozone loss and thus on the climate influence of these gases. Here we use observed ozone loss to reduce the resultant uncertainty in their overall climate influence and infer a smaller warming influence of these substances than was considered likely in a 2013 climate report. The result implies a smaller benefit to climate due to their phase‐out, mandated under the Montreal Protocol, than would have been the case under previous understanding. Key Points: Effective radiative forcing of ozone‐depleting substances, as discerned from CMIP6 simulations, spans a large rangeUsing an emergent constraint approach, our new estimate is consistent with observational climatologies of total‐column ozoneThis range implies a smaller forcing than the estimate provided by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [ABSTRACT FROM AUTHOR]

The connection between the dominant mode of interannual variability in the tropical troposphere, El Niño Southern Oscillation (ENSO), and entry of stratospheric water vapor, is analyzed in a set of the model simulations archived for the Chemistry-Climate Model Initiative (CCMI) project and for phase 6 of the Coupled Model Intercomparison Project. While the models agree on the temperature response to ENSO in the tropical troposphere and lower stratosphere, and all models also agree on the zonal structure of the response in the tropical tropopause layer, the only aspect of the entry water vapor with consensus is that La Niña leads to moistening in winter relative to neutral ENSO. For El Niño and for other seasons there are significant differences among the models. For example, some models find that the enhanced water vapor for La Niña in the winter of the event reverses in spring and summer, other models find that this moistening persists, while some show a nonlinear response with both El Niño and La Niña leading to enhanced water vapor in both winter, spring, and summer. Focusing on Central Pacific ENSO versus East Pacific ENSO, or temperatures in the mid-troposphere as compared to temperatures near the surface, does not narrow the inter-model discrepancies. Despite this diversity in response, the temperature response near the cold point can explain the response of water vapor when each model is considered separately. While the observational record is too short to fully constrain the response to ENSO, it is clear that most models suffer from biases in the magnitude of interannual variability of entry water vapor. This bias could be due to missing forcing processes that contribute to observed variability in cold point temperatures. [ABSTRACT FROM AUTHOR]

The present generation of global climate models is characterised by insufficient reflection of short-wave radiation over the Southern Ocean due to a misrepresentation of clouds. This is a significant concern as it leads to excessive heating of the ocean surface, sea surface temperature biases and subsequent problems with atmospheric dynamics. In this study, we modify cloud microphysics in a recent version of the Met Office's Unified Model and show that choosing a more realistic value for the shape parameter of atmospheric ice crystals, in better agreement with theory and observations, benefits the simulation of short-wave radiation. In the model, for calculating the growth rate of ice crystals through deposition, the default assumption is that all ice particles are spherical in shape. We modify this assumption to effectively allow for oblique shapes or aggregates of ice crystals. Along with modified ice nucleation temperatures, we achieve a reduction in the annual-mean short-wave cloud radiative effect over the Southern Ocean by up to ∼4 W m -2 and seasonally much larger reductions compared to the control model. By slowing the growth of the ice phase, the model simulates substantially more supercooled liquid cloud. [ABSTRACT FROM AUTHOR]

This paper describes the development of New Zealand's Earth System Model (NZESM) and evaluates its performance against its parent model (United Kingdom Earth System Model, UKESM) and observations. The main difference between the two earth system models is an embedded high‐resolution (1/5°) nested region over the oceans around New Zealand in the NZESM. Due to this finer ocean model mesh, currents such as the East Australian Current, East Australian Current Extension, Tasman Front, and Tasman Leakage, and their volume and heat transports are better simulated in the NZESM. The improved oceanic transports have led to a reduction in upper ocean temperature and salinity biases over the nested region. In addition, net transports through the Tasman Sea of volume, heat and salt in the NZESM agree better with previously reported estimates. A consequence of the increased cross‐Tasman Sea transports in the NZESM is increased temperatures and salinity west of Australia and in the Southern Ocean reducing the meridional sea surface temperature gradient between the subtropics and sub‐Antarctic. This also leads to a weakening of the westerly winds between 60°S and 45°S over large parts of the Southern Ocean, which reduces the northward Ekman transport, reduces the formation of Antarctic Intermediate Water, and allows for a southward expansion of the Super‐Gyre in all ocean basins. Connecting an improved oceanic circulation around New Zealand to a basin‐wide Super‐Gyre response is an important step forward in our current understanding of how local scales can influence global scales in a fully coupled earth system model. Plain Language Summary: We describe the model development of the New Zealand Earth System Model and assess its performance against the model on which it is based (United Kingdom Earth System Model) and observations. The New Zealand Earth System Model is a fully coupled earth system model, which aims to model all relevant bio‐physical processes in and between the atmosphere, land, ocean, and sea‐ice. The main difference between both models is that the oceans around New Zealand in the New Zealand Earth System Model are more precisely modeled, due to a refined ocean model mesh in this region. That results in a more accurate oceanic circulation around New Zealand in the New Zealand Earth System Model compared to the United Kingdom Earth System Model and reduced model biases of temperature and salinity. These oceanic changes have implications beyond the oceans around New Zealand, causing a warming in the Southern Ocean and a related weakening of the westerly winds over the Southern Ocean. This weakening of the winds allows subtropical waters to reach further south into the Southern Ocean. It is notable that regional changes in the ocean circulation can have implications on the global scale. Key Points: NZESM is a nested fully coupled ESM based on UKESM with a high‐resolution ocean grid of 1/5° around New ZealandThe oceanic circulation is improved in NZESM over the nested domain and model biases of temperature and salinity are reducedThe Super‐Gyre intensifies and expands southward due to wind changes triggered by changes in the large‐scale heat transport [ABSTRACT FROM AUTHOR]

Ozone profile measurements collected at L'Aquila (Italy, 42.4° N) during seventeen years of radio-sounding (2000-2016) are presented here, with an analysis of derived trends. Model results from the SPARC-CCMI exercise are used in parallel to highlight the physical and chemical mechanisms regulating mid-latitude ozone trends. The statistically significant trends highlighted in time series at L'Aquila are those in the mid-upper stratosphere (+5.9 ± 4.2), mid troposphere (+5.9 ± 2.4) and upper troposphere (+2.5 ± 0.9), all in percent/decade. The upper stratospheric positive trend was already well documented in recent WMO assessments and attributed to the starting decline of stratospheric Cly and Bry and to the stratospheric cooling induced by increasing well mixed greenhouse gases, thus slowing down gas-phase reactions that destroy ozone in the upper stratosphere. The ozone increase in the mid-upper troposphere is largely regulated by the increasing strength of the Brewer-Dobson circulation, which moves more ozone from the tropics to the extratropics and enhances the tropospheric influx from the lowermost stratosphere. This climate feedback mechanism on tropospheric ozone is only partially compensated by the increasing chemical ozone loss associated to higher H2O values in response to the tropospheric warming. We also note that ozone trends obtained in the lower stratosphere are negative (-2.2 percent/decade), but do not result to be statistically significant in our analyses. [ABSTRACT FROM AUTHOR]

Southern Ocean (SO) shortwave (SW) radiation biases are a common problem in contemporary general circulation models (GCMs), with most models exhibiting a tendency to absorb too much incoming SW radiation. These biases have been attributed to deficiencies in the representation of clouds during the austral summer months, either due to cloud cover or cloud albedo being too low. The problem has been the focus of many studies, most of which utilised satellite datasets for model evaluation. We use multi-year ship-based observations and the CERES spaceborne radiation budget measurements to contrast cloud representation and SW radiation in the atmospheric component Global Atmosphere (GA) version 7.1 of the HadGEM3 GCM and the MERRA-2 reanalysis. We find that the prevailing bias is negative in GA7.1 and positive in MERRA-2. GA7.1 performs better than MERRA-2 in terms of absolute SW bias. Significant errors of up to 21 W m -2 (GA7.1) and 39 W m -2 (MERRA-2) are present in both models in the austral summer. Using ship-based ceilometer observations, we find low cloud below 2 km to be predominant in the Ross Sea and the Indian Ocean sectors of the SO. Utilising a novel surface lidar simulator developed for this study, derived from an existing Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) – active remote sensing simulator (ACTSIM) spaceborne lidar simulator, we find that GA7.1 and MERRA-2 both underestimate low cloud and fog occurrence relative to the ship observations on average by 4 %–9 % (GA7.1) and 18 % (MERRA-2). Based on radiosonde observations, we also find the low cloud to be strongly linked to boundary layer atmospheric stability and the sea surface temperature. GA7.1 and MERRA-2 do not represent the observed relationship between boundary layer stability and clouds well. We find that MERRA-2 has a much greater proportion of cloud liquid water in the SO in austral summer than GA7.1, a likely key contributor to the difference in the SW radiation bias. Our results suggest that subgrid-scale processes (cloud and boundary layer parameterisations) are responsible for the bias and that in GA7.1 a major part of the SW radiation bias can be explained by cloud cover underestimation, relative to underestimation of cloud albedo. [ABSTRACT FROM AUTHOR]

One of the key questions in the air quality and climate sciences is how tropospheric ozone concentrations will change in the future. This will depend on two factors: changes in stratosphere-to-troposphere transport (STT) and changes in tropospheric chemistry. Here we aim to identify robust changes in STT using simulations from the Chemistry Climate Model Initiative (CCMI) under a common climate change scenario (RCP6.0). We use two idealized stratospheric tracers to isolate changes in transport: stratospheric ozone (O3S), which is exactly like ozone but has no chemical sources in the troposphere, and st80, a passive tracer with fixed volume mixing ratio in the stratosphere. We find a robust increase in the tropospheric columns of these two tracers across the models. In particular, stratospheric ozone in the troposphere is projected to increase 10 %–16 % by the end of the 21st century in the RCP6.0 scenario. Future STT is enhanced in the subtropics due to the strengthening of the shallow branch of the Brewer–Dobson circulation (BDC) in the lower stratosphere and of the upper part of the Hadley cell in the upper troposphere. The acceleration of the deep branch of the BDC in the Northern Hemisphere (NH) and changes in eddy transport contribute to increased STT at high latitudes. These STT trends are caused by greenhouse gas (GHG) increases, while phasing out of ozone-depleting substances (ODS) does not lead to robust transport changes. Nevertheless, the decline of ODS increases the reservoir of ozone in the lower stratosphere, which results in enhanced STT of O3S at middle and high latitudes. A higher emission scenario (RCP8.5) produces stronger STT trends, with increases in tropospheric column O3S more than 3 times larger than those in the RCP6.0 scenario by the end of the 21st century. [ABSTRACT FROM AUTHOR]

Automatic lidars and ceilometers provide valuable information on cloud and aerosols, but have not been used systematically in the evaluation of GCMs and NWP models. Obstacles associated with the diversity of instruments, a lack of standardisation of data products and open processing tools mean that the value of the large ALC networks worldwide is not being realised. We discuss a tool, called the Automatic Lidar and Ceilometer Framework (ALCF), that overcomes these problems and also includes a ground-based lidar simulator, which calculates the radiative transfer of laser radiation, and allows one-to-one comparison with models. Our ground-based lidar simulator is based on the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) which has been used extensively for spaceborne lidar intercomparisons. The ALCF implements all steps needed to transform and calibrate raw ALC data and create simulated backscatter profiles for one-to-one comparison and complete statistical analysis of cloud. The framework supports multiple common commercial ALCs (Vaisala CL31, CL51, Lufft CHM 15k and Sigma Space MiniMPL), reanalyses (JRA-55, ERA5 and MERRA-2) and models (AMPS and the Unified Model). To demonstrate its capabilities, we present case studies evaluating cloud in the supported reanalyses and models using CL31, CL51, CHM 15k and MiniMPL observations at three sites in New Zealand. We show that the reanalyses and models generally underestimate cloud fraction and overestimate cloud albedo, the common too few too bright problem. If sufficiently high temporal resolution model output is available (better than 6 hourly), a direct comparison of individual clouds is also possible. We demonstrate that the ALCF can be used as a generic evaluation tool to examine cloud occurrence and cloud properties in reanalyses, NWP models and GCMs, potentially utilising the large amounts of ALC data already available. This tool is likely to be particularly useful for the analysis and improvement of low-level cloud simulations which are not well monitored from space. This has previously been identified as a critical deficiency in contemporary models, limiting the accuracy of weather forecasts and future climate projections. [ABSTRACT FROM AUTHOR]

Here we present a description of the UKCA StratTrop chemical mechanism, which is used in the UKESM1 Earth system model for CMIP6. The StratTrop chemical mechanism is a merger of previously well-evaluated tropospheric and stratospheric mechanisms, and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained for tropospheric ozone and its budget terms from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address. [ABSTRACT FROM AUTHOR]

Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here we evaluate long-term changes in these species from the pre-industrial (1850) to the end of the 21st century in CMIP6 models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations, although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global total column ozone (TCO) has increased from ∼300 DU in 1850 to ∼305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0 and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer-Dobson Circulation under other SSPs. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼0.5 ppmv from the pre-industrial to the present day and are projected to increase further by the end of the 21st century. The largest increases (∼2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Both TCO and tropical lower stratospheric water vapour show large variability following explosive volcanic eruptions. [ABSTRACT FROM AUTHOR]

Quantifying forcings from anthropogenic perturbations to the Earth System (ES) is important for understanding changes in climate since the pre-industrial period. In this paper, we quantify and analyse a wide range of present-day (PD) anthropogenic climate forcings with the UK's Earth System Model (ESM), UKESM1, following the protocols defined by the Radiative Forcing Model Intercomparison Project (RFMIP) and the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP). In particular, by quantifying effective radiative forcings (ERFs) that include rapid adjustments within a full ESM, it enables the role of various climate-chemistry-aerosol-cloud feedbacks to be quantified. Global mean ERFs are 1.83, 0.13, −0.33, and 0.93 W m−2 at the PD (Year 2014) relative to the pre-industrial (PI; Year 1850) for carbon dioxide, nitrous oxide, ozone-depleting substances, and methane, respectively. The PD total greenhouse gas ERF is 2.89 W m−2, larger than the sum of the individual GHG ERFs. UKESM1 has an aerosol forcing of −1.13 W m−2. A relatively strong negative forcing from aerosol-cloud interactions and a small negative instantaneous forcing from aerosol-radiation interactions are partially offset by a substantial forcing from black carbon absorption. Internal mixing and chemical interactions mean that neither the forcing from aerosol-radiation interactions nor aerosol-cloud interactions are linear, making the total aerosol ERF less than the sum of the individual speciated aerosol ERFs. Tropospheric ozone precursors, in addition to exerting a positive forcing due to ozone, lead to oxidant changes which in turn cause an indirect aerosol ERF, altering the sign of the net ERF from nitrogen oxide emissions. Together, aerosol and tropospheric ozone precursors (near-term climate forcers, NTCFs) exert a global mean ERF of −1.12 W m−2, mainly due to changes in the cloud radiative effect. There is also a negative PD ERF from land use (−0.32 W m−2). It is outside the range of previous estimates, and is most likely due to too strong an albedo response. In combination, the net anthropogenic ERF is potentially biased low (1.61 W m−2) relative to other estimates, due to the inclusion of non-linear feedbacks and ES interactions. By including feedbacks between greenhouse gases, stratospheric and tropospheric ozone, aerosols, and clouds, some of which act non-linearly, this work demonstrates the importance of ES interactions when quantifying climate forcing. It also suggests that rapid adjustments need to include chemical as well as physical adjustments to fully account for complex ES interactions. [ABSTRACT FROM AUTHOR]

The hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting methane lifetime (τCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of τCH4 differences among 10 models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D), the mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapour, carbon monoxide (CO), the ratio of NO:NOx , and formaldehyde (HCHO) explain moderate differences in τCH4 , while isoprene, methane, the photolysis frequency of NO2 by visible light (JNO2), overhead ozone column, and temperature account for little to no model variation in τCH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in τCH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3 , JO1D , NOx , and H2O , with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in τCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapour on OH and τCH4 , imparting an increasing and decreasing trend of about 0.5 % decade -1 , respectively. The responses due to NOx , ozone column, and temperature are also in reasonably good agreement between the two studies. [ABSTRACT FROM AUTHOR]

With low concentrations of tropospheric aerosol, the Southern Ocean offers a "natural laboratory" for studies of aerosol–cloud interactions. Aerosols over the Southern Ocean are produced from biogenic activity in the ocean, which generates sulfate aerosol via dimethylsulfide (DMS) oxidation, and from strong winds and waves that lead to bubble bursting and sea spray emission. Here, we evaluate the representation of Southern Ocean aerosols in the Hadley Centre Global Environmental Model version 3, Global Atmosphere 7.1 (HadGEM3-GA7.1) chemistry–climate model. Compared with aerosol optical depth (AOD) observations from two satellite instruments (the Moderate Resolution Imaging Spectroradiometer, MODIS-Aqua c6.1, and the Multi-angle Imaging Spectroradiometer, MISR), the model simulates too-high AOD during winter and too-low AOD during summer. By switching off DMS emission in the model, we show that sea spray aerosol is the dominant contributor to AOD during winter. In turn, the simulated sea spray aerosol flux depends on near-surface wind speed. By examining MODIS AOD as a function of wind speed from the ERA-Interim reanalysis and comparing it with the model, we show that the sea spray aerosol source function in HadGEM3-GA7.1 overestimates the wind speed dependency. We test a recently developed sea spray aerosol source function derived from measurements made on a Southern Ocean research voyage in 2018. In this source function, the wind speed dependency of the sea spray aerosol flux is less than in the formulation currently implemented in HadGEM3-GA7.1. The new source function leads to good agreement between simulated and observed wintertime AODs over the Southern Ocean; however, it reveals partially compensating errors in DMS-derived AOD. While previous work has tested assumptions regarding the seawater climatology or sea–air flux of DMS, we test the sensitivity of simulated AOD, cloud condensation nuclei and cloud droplet number concentration to three atmospheric sulfate chemistry schemes. The first scheme adds DMS oxidation by halogens and the other two test a recently developed sulfate chemistry scheme for the marine troposphere; one tests gas-phase chemistry only, while the second adds extra aqueous-phase sulfate reactions. We show how simulated sulfur dioxide and sulfuric acid profiles over the Southern Ocean change as a result and how the number concentration and particle size of the soluble Aitken, accumulation and coarse aerosol modes are affected. The new DMS chemistry scheme leads to a 20 % increase in the number concentration of cloud condensation nuclei and cloud droplets, which improves agreement with observations. Our results highlight the importance of atmospheric chemistry for simulating aerosols and clouds accurately over the Southern Ocean. [ABSTRACT FROM AUTHOR]

The present generation of global climate models is characterized by insufficient reflection of short-wave radiation over the Southern Ocean due to a misrepresentation of clouds. This is a significant concern as it leads to excessive heating of the ocean surface, sea surface temperature biases, and subsequent problems with atmospheric dynamics. In this study we modify cloud micro-physics in a recent version of the Met Office's Unified Model and show that choosing a more realistic value for the shape parameter of atmospheric ice-crystals, in better agreement with theory and observations, benefits the simulation of short-wave radiation. In the model, for calculating the growth rate of ice crystals through deposition, the default assumption is that all ice particles are spherical in shape. We modify this assumption to effectively allow for oblique shapes or aggregates of ice crystals. Along with modified ice nucleation temperatures, we achieve a reduction in the annual-mean short-wave cloud radiative effect over the Southern Ocean by up to 4 W/m2, and seasonally much larger reductions. By slowing the growth of the ice phase, the model simulates substantially more supercooled liquid cloud. We hypothesize that such abundant supercooled liquid cloud is the result of a paucity of ice nucleating particles in this part of the atmosphere. [ABSTRACT FROM AUTHOR]

Here we present a description of the UKCA StratTrop chemical mechanism which is used in the UKESM1 Earth System Model for CMIP6. The StratTrop chemical mechanism is a merger of previously well evaluated tropospheric and stratospheric mechanisms and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address. [ABSTRACT FROM AUTHOR]

Hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting CH4 lifetime (τCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of τCH4 differences among ten models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO¹D) due mostly to clouds, mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapor, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in τCH4, while isoprene, CH4, the photolysis frequency of NO2 by visible light (JNO2), overhead O3 column, and temperature account for little-to-no model variation in τCH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in τCH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO¹D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in τCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapor on OH and τCH4, imparting an increasing and decreasing trend of about 0.5 % decade−1, respectively. The responses due to NOx, O3 column, and temperature are also in reasonably good agreement between the two studies, though a discrepancy in the CH4 response highlights a need for further examination of the CH4 feedback on the abundance of OH. [ABSTRACT FROM AUTHOR]

We have derived values of the ultraviolet index (UVI) at solar noon using the Tropospheric Ultraviolet Model (TUV) driven by ozone, temperature and aerosol fields from climate simulations of the first phase of the Chemistry-Climate Model Initiative (CCMI-1). Since clouds remain one of the largest uncertainties in climate projections, we simulated only the clear-sky UVI. We compared the modelled UVI climatologies against present-day climatological values of UVI derived from both satellite data (the OMI-Aura OMUVBd product) and ground-based measurements (from the NDACC network). Depending on the region, relative differences between the UVI obtained from CCMI/TUV calculations and the ground-based measurements ranged between -5.9 % and 10.6 %. We then calculated the UVI evolution throughout the 21st century for the four Representative Concentration Pathways (RCPs 2.6, 4.5, 6.0 and 8.5). Compared to 1960s values, we found an average increase in the UVI in 2100 (of 2 %–4 %) in the tropical belt (30 ∘ N–30 ∘ S). For the mid-latitudes, we observed a 1.8 % to 3.4 % increase in the Southern Hemisphere for RCPs 2.6, 4.5 and 6.0 and found a 2.3 % decrease in RCP 8.5. Higher increases in UVI are projected in the Northern Hemisphere except for RCP 8.5. At high latitudes, ozone recovery is well identified and induces a complete return of mean UVI levels to 1960 values for RCP 8.5 in the Southern Hemisphere. In the Northern Hemisphere, UVI levels in 2100 are higher by 0.5 % to 5.5 % for RCPs 2.6, 4.5 and 6.0 and they are lower by 7.9 % for RCP 8.5. We analysed the impacts of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) on UVI from 1960 by comparing CCMI sensitivity simulations (1960–2100) with fixed GHGs or ODSs at their respective 1960 levels. As expected with ODS fixed at their 1960 levels, there is no large decrease in ozone levels and consequently no sudden increase in UVI levels. With fixed GHG, we observed a delayed return of ozone to 1960 values, with a corresponding pattern of change observed on UVI, and looking at the UVI difference between 2090s values and 1960s values, we found an 8 % increase in the tropical belt during the summer of each hemisphere. Finally we show that, while in the Southern Hemisphere the UVI is mainly driven by total ozone column, in the Northern Hemisphere both total ozone column and aerosol optical depth drive UVI levels, with aerosol optical depth having twice as much influence on the UVI as total ozone column does. [ABSTRACT FROM AUTHOR]

The Northern Hemisphere and tropical circulation response to interannual variability in Arctic stratospheric ozone is analyzed in a set of the latest model simulations archived for the Chemistry-Climate Model Initiative (CCMI) project. All models simulate a connection between ozone variability and temperature/geopotential height in the lower stratosphere similar to that observed. A connection between Arctic ozone variability and polar cap surface air pressure is also found, but additional statistical analysis suggests that it is mediated by the dynamical variability that typically drives the anomalous ozone concentrations. While the CCMI models also show a connection between Arctic stratospheric ozone and the El Niño–Southern Oscillation (ENSO), with Arctic stratospheric ozone variability leading to ENSO variability 1 to 2 years later, this relationship in the models is much weaker than observed and is likely related to ENSO autocorrelation rather than any forced response to ozone. Overall, Arctic stratospheric ozone is related to lower stratospheric variability. Arctic stratospheric ozone may also influence the surface in both polar and tropical latitudes, though ozone is likely not the proximate cause of these impacts and these impacts can be masked by internal variability if data are only available for ∼40 years. [ABSTRACT FROM AUTHOR]

Climate model simulations show an acceleration of the Brewer–Dobson circulation (BDC) in response to climate change. While the general mechanisms for the BDC strengthening are widely understood, there are still open questions concerning the influence of the details of the wave driving. Mean age of stratospheric air (AoA) is a useful transport diagnostic for assessing changes in the BDC. Analyzing AoA from a subset of Chemistry–Climate Model Initiative part 1 climate projection simulations, we find a remarkable agreement between most of the models in simulating the largest negative AoA trends in the extratropical lower to middle stratosphere of both hemispheres (approximately between 20 and 25 geopotential kilometers (gpkm) and 20–50 ∘ N and S). We show that the occurrence of AoA trend minima in those regions is directly related to the climatological AoA distribution, which is sensitive to an upward shift of the circulation in response to climate change. Also other factors like a reduction of aging by mixing (AbM) and residual circulation transit times (RCTTs) contribute to the AoA distribution changes by widening the AoA isolines. Furthermore, we analyze the time evolution of AbM and RCTT trends in the extratropics and examine the connection to possible drivers focusing on local residual circulation strength, net tropical upwelling and wave driving. However, after the correction for a vertical shift of pressure levels, we find only seasonally significant trends of residual circulation strength and zonal mean wave forcing (resolved and unresolved) without a clear relation between the trends of the analyzed quantities. This indicates that additional causative factors may influence the AoA, RCTT and AbM trends. In this study, we postulate that the shrinkage of the stratosphere has the potential to influence the RCTT and AbM trends and thereby cause additional AoA changes over time. [ABSTRACT FROM AUTHOR]

One of the key questions in the air quality and climate sciences is how will tropospheric ozone concentrations change in the future. This will depend on two factors: changes in stratosphere-to-troposphere transport (STT) and changes in tropospheric chemistry. Here we aim to identify robust changes in STT using simulations from the Chemistry Climate Model Initiative (CCMI) under a common climate change scenario (RCP6.0). We use two idealized stratospheric tracers to isolate changes in transport: stratospheric ozone (O3S), which is exactly like ozone but has no chemical sources in the troposphere, and st80, a passive tracer with fixed volume mixing ratio in the stratosphere. We find a robust increase in the tropospheric columns of these two tracers across the models. In particular, stratospheric ozone in the troposphere is projected to increase 10–16 % by the end of the 21st century in the RCP6.0 scenario. Future STT is enhanced in the subtropics due to the strengthening of the shallow branch of the Brewer-Dobson circulation (BDC) in the lower stratosphere and of the upper part of the Hadley cell in the upper troposphere. The acceleration of the deep branch of the BDC and changes in eddy transport contribute to increase STT at high latitudes. The idealized tracer st80 shows that these STT changes are dominated by greenhouse gas (GHG) increases, while phasing out of ozone depleting substances (ODS) does not lead to robust STT changes. Nevertheless, the increase of O3S concentrations in the troposphere is attributed to GHG only in the subtropics. At middle and high latitudes it is due to stratospheric ozone recovery linked to ODS decline. A higher emission scenario (RCP8.5) produces qualitatively similar but stronger STT trends, with changes in tropospheric column O3S more than three times larger than those in the RCP6.0 scenario by the end of the 21st century. [ABSTRACT FROM AUTHOR]

Transport from the Northern Hemisphere (NH) midlatitudes to the Arctic plays a crucial role in determining the abundance of trace gases and aerosols that are important to Arctic climate via impacts on radiation and chemistry. Here we examine this transport using an idealized tracer with a fixed lifetime and predominantly midlatitude land-based sources in models participating in the Chemistry Climate Model Initiative (CCMI). We show that there is a 25 %–45 % difference in the Arctic concentrations of this tracer among the models. This spread is correlated with the spread in the location of the Pacific jet, as well as the spread in the location of the Hadley Cell (HC) edge, which varies consistently with jet latitude. Our results suggest that it is likely that the HC-related zonal-mean meridional transport rather than the jet-related eddy mixing is the major contributor to the inter-model spread in the transport of land-based tracers into the Arctic. Specifically, in models with a more northern jet, the HC generally extends further north and the tracer source region is mostly covered by surface southward flow associated with the lower branch of the HC, resulting in less efficient transport poleward to the Arctic. During boreal summer, there are poleward biases in jet location in free-running models, and these models likely underestimate the rate of transport into the Arctic. Models using specified dynamics do not have biases in the jet location, but do have biases in the surface meridional flow, which may result in differences in transport into the Arctic. In addition to the land-based tracer, the midlatitude-to-Arctic transport is further examined by another idealized tracer with zonally uniform sources. With equal sources from both land and ocean, the inter-model spread of this zonally uniform tracer is more related to variations in parameterized convection over oceans rather than variations in HC extent, particularly during boreal winter. This suggests that transport of land-based and oceanic tracers or aerosols towards the Arctic differs in pathways and therefore their corresponding inter-model variabilities result from different physical processes. [ABSTRACT FROM AUTHOR]

Southern Ocean (SO) shortwave (SW) radiation biases are a common problem in contemporary general circulation models (GCMs), with most models exhibiting a tendency to absorb too much incoming SW radiation. These biases have been attributed to deficiencies in the representation of clouds during the austral summer months, either due to cloud cover or cloud optical thickness being too low. The problem has been the focus of many studies, most of which utilised satellite datasets for model evaluation. We use multi-year ship based observations and the CERES spaceborne radiation budget measurements to contrast cloud representation and SW radiation in the atmospheric component Global Atmosphere (GA) version 7.0 and 7.1 of the HadGEM3 GCM and the MERRA-2 reanalysis. We find that MERRA-2 is biased in the opposite direction to GA (reflects too much SW radiation). In addition, MERRA-2 performs better in terms of absolute SW bias than nudged runs of GA7.0 and GA7.1 in the 60-70°S latitude band. GA7.1 reduces the SO SW radiation biases relative to GA7.0, but significant errors remain at up to 20Wm-2 between 60 and 70°S in the austral summer months. Using ship-based ceilometer observations, we find low cloud below 2km to be predominant in the Ross Sea and the Indian Ocean sector of the SO. Utilising a novel surface lidar simulator developed for this study, derived from an existing COSP-ACTSIM spaceborne lidar simulator, we find that GA7.0 and MERRA-2 both underestimate low cloud occurrence relative to the ship observations by 18-25% on average, though the cloud cover in MERRA-2 is closer to observations by about 7%. Based on radiosonde observations, we also find the low cloud to be strongly linked to boundary-layer atmospheric stability and the sea surface temperature. GA7.0 and MERRA-2 agree well with observations in terms of boundary-layer stability, suggesting that subgrid-scale parametrisations do not generate enough cloud in response to the thermodynamic profile of the atmosphere and the surface temperature. Our analysis shows that MERRA-2 has a much greater proportion of cloud liquid water in the SO in January than GA7.0, a likely key contributor to the difference in SW radiation. We show that boundary-layer stability and relative humidity fields are very similar in GA7.0 and MERRA-2, and unlikely to be the cause of the different cloud representation, suggesting that subgrid-scale parametrisations are responsible for the difference between the models. [ABSTRACT FROM AUTHOR]

Model evaluation is difficult and generally relies on analysis that can mask compensating errors. This paper defines new metrics, using clusters generated from a machine learning algorithm, to estimate mean and compensating errors in different model runs. As a test case, we investigate the Southern Ocean shortwave radiative bias using clusters derived by applying self‐organizing maps to satellite data. In particular, the effects of changing cloud phase parameterizations in the MetOffice Unified Model are examined. Differences in cluster properties show that the regional radiative biases are substantially different than the global bias, with two distinct regions identified within the Southern Ocean, each with a different signed bias. Changing cloud phase parameterizations can reduce errors at higher latitudes but increase errors at lower latitudes of the Southern Ocean. Ranking the parameterizations often shows a contrast in mean and compensating errors, notably in all cases large compensating errors remain. Key Points: A novel method for climate model evaluation is used to identify both mean errors and potential compensating errorsAs an example, we apply this methodology to investigate the quality of different cloud parameterizations over the Southern OceanChanges to the cloud phase parameterizations can reduce shortwave radiative bias regionally, but large compensating errors remain [ABSTRACT FROM AUTHOR]