Global aviation operations contribute to anthropogenic climate change via a complex set of processes that lead to a net surface warming. Of importance are aviation emissions of carbon dioxide (CO2), nitrogen oxides (NOx), water vapor, soot and sulfate aerosols, and increased cloudiness due to contrail formation. Aviation grew strongly over the past decades (1960-2018) in terms of activity, with revenue passenger kilometers increasing from 109 to 8269 billion km yr-1, and in terms of climate change impacts, with CO2 emissions increasing by a factor of 6.8 to 1034 Tg CO2 yr-1. Over the period 2013-2018, the growth rates in both terms show a marked increase. Here, we present a new comprehensive and quantitative approach for evaluating aviation climate forcing terms. Both radiative forcing (RF) and effective radiative forcing (ERF) terms and their sums are calculated for the years 2000-2018. Contrail cirrus, consisting of linear contrails and the cirrus cloudiness arising from them, yields the largest positive net (warming) ERF term followed by CO2 and NOx emissions. The formation and emission of sulfate aerosol yields a negative (cooling) term. The mean contrail cirrus ERF/RF ratio of 0.42 indicates that contrail cirrus is less effective in surface warming than other terms. For 2018 the net aviation ERF is +100.9 milliwatts (mW) m-2 (5-95% likelihood range of (55, 145)) with major contributions from contrail cirrus (57.4 mW m-2), CO2 (34.3 mW m-2), and NOx (17.5 mW m-2). Non-CO2 terms sum to yield a net positive (warming) ERF that accounts for more than half (66%) of the aviation net ERF in 2018. Using normalization to aviation fuel use, the contribution of global aviation in 2011 was calculated to be 3.5 (4.0, 3.4) % of the net anthropogenic ERF of 2290 (1130, 3330) mW m-2. Uncertainty distributions (5%, 95%) show that non-CO2 forcing terms contribute about 8 times more than CO2 to the uncertainty in the aviation net ERF in 2018. The best estimates of the ERFs from aviation aerosol-cloud interactions for soot and sulfate remain undetermined. CO2-warming-equivalent emissions based on global warming potentials (GWP* method) indicate that aviation emissions are currently warming the climate at approximately three times the rate of that associated with aviation CO2 emissions alone. CO2 and NOx aviation emissions and cloud effects remain a continued focus of anthropogenic climate change research and policy discussions.

In the context of a changing climate, the acceleration of the Brewer-Dobson circulation [Butchart 2014] leads to a decrease of the ozone total column in the tropics. This decrease affects directly surface ultraviolet radiation, which are already very high in this area. Following the work of (Bais et al., 2011), (Butchart, 2014) and (Hegglin & Shepherd, 2009) on the future evolution of surface irradiance derived from Chemistry Climate Models (CCM) projections, we projected here surface irradiance from 2010 to 2100 with focus on the tropics. We used the latest chemistry climate projection exercise ; Chemistry Climate Model Initiative (CCMI) coupled with a radiative transfer model (TUV (Madronich, 1993)) to calculate the evolution of surface Ultraviolet radiation throughout the 21st century. Ultraviolet Index (UVi) has been specifically considered (McKenzie, Matthews, & Johnston, 1991). At first, simulation from RefC2 Chemistry Climate Model Initiative have been coupled with a radiative transfer model, in order to obtained modeled UV index (UVi-M). UVi-M is then compared against available satellite ultraviolet radiation observations (OMI OMUVbd product) between 2005 and 2016. Statistical difference and variance have been analysed versus different parameters: geographical location, model or ensemble of model outputs used in the radiative transfer calculation.

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137Cs and 133Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τe, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τe lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

A new Geoengineering Model Intercomparison Project (GeoMIP) experiment "G4 specified stratospheric aerosols" (short name: G4SSA) is proposed to investigate the impact of stratospheric aerosol geoengineering on atmosphere, chemistry, dynamics, climate, and the environment. In contrast to the earlier G4 GeoMIP experiment, which requires an emission of sulfur dioxide (SO2) into the model, a prescribed aerosol forcing file is provided to the community, to be consistently applied to future model experiments between 2020 and 2100. This stratospheric aerosol distribution, with a total burden of about 2 Tg S has been derived using the ECHAM5-HAM microphysical model, based on a continuous annual tropical emission of 8 Tg SO2 yr−1. A ramp-up of geoengineering in 2020 and a ramp-down in 2070 over a period of 2 years are included in the distribution, while a background aerosol burden should be used for the last 3 decades of the experiment. The performance of this experiment using climate and chemistry models in a multi-model comparison framework will allow us to better understand the impact of geoengineering and its abrupt termination after 50 years in a changing environment. The zonal and monthly mean stratospheric aerosol input data set is available at https://www2.acd.ucar.edu/gcm/geomip-g4-specified-stratospheric-aerosol-data-set.

We have diagnosed the lifetimes of long‐lived source gases emitted at the surface and removed in the stratosphere using six three‐dimensional chemistry‐climate models and a two‐dimensional model. The models all used the same standard photochemical data. We investigate the effect of different definitions of lifetimes, including running the models with both mixing ratio (MBC) and flux (FBC) boundary conditions. Within the same model, the lifetimes diagnosed by different methods agree very well. Using FBCs versus MBCs leads to a different tracer burden as the implied lifetime contained in the MBC value does not necessarily match a model's own calculated lifetime. In general, there are much larger differences in the lifetimes calculated by different models, the main causes of which are variations in the modeled rates of ascent and horizontal mixing in the tropical midlower stratosphere. The model runs have been used to compute instantaneous and steady state lifetimes. For chlorofluorocarbons (CFCs) their atmospheric distribution was far from steady state in their growth phase through to the 1980s, and the diagnosed instantaneous lifetime is accordingly much longer. Following the cessation of emissions, the resulting decay of CFCs is much closer to steady state. For 2100 conditions the model circulation speeds generally increase, but a thicker ozone layer due to recovery and climate change reduces photolysis rates. These effects compensate so the net impact on modeled lifetimes is small. For future assessments of stratospheric ozone, use of FBCs would allow a consistent balance between rate of CFC removal and model circulation rate.

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng g−1 for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng g−1 for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng g−1. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model–measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60–90° N) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07–0.25) W m−2 and 0.18 (0.06–0.28) W m−2 in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W m−2 for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

In this study we examine the performance of 31 global model radiative transfer schemes in cloud-free conditions with prescribed gaseous absorbers and no aerosols (Rayleigh atmosphere), with prescribed scattering-only aerosols, and with more absorbing aerosols. Results are compared to benchmark results from high-resolution, multi-angular line-by-line radiation models. For purely scattering aerosols, model bias relative to the line-by-line models in the top-of-the atmosphere aerosol radiative forcing ranges from roughly −10 to 20%, with over- and underestimates of radiative cooling at lower and higher solar zenith angle, respectively. Inter-model diversity (relative standard deviation) increases from ~10 to 15% as solar zenith angle decreases. Inter-model diversity in atmospheric and surface forcing decreases with increased aerosol absorption, indicating that the treatment of multiple-scattering is more variable than aerosol absorption in the models considered. Aerosol radiative forcing results from multi-stream models are generally in better agreement with the line-by-line results than the simpler two-stream schemes. Considering radiative fluxes, model performance is generally the same or slightly better than results from previous radiation scheme intercomparisons. However, the inter-model diversity in aerosol radiative forcing remains large, primarily as a result of the treatment of multiple-scattering. Results indicate that global models that estimate aerosol radiative forcing with two-stream radiation schemes may be subject to persistent biases introduced by these schemes, particularly for regional aerosol forcing.

We use observations of N2O and mean age to identify realistic transport in models in order to explain their ozone predictions. The results are applied to 15 chemistry climate models (CCMs) participating in the 2010 World Meteorological Organization ozone assessment. Comparison of the observed and simulated N2O, mean age and their compact correlation identifies models with fast or slow circulations and reveals details of model ascent and tropical isolation. This process-oriented diagnostic is more useful than mean age alone because it identifies models with compensating transport deficiencies that produce fortuitous agreement with mean age. The diagnosed model transport behavior is related to a model's ability to produce realistic lower stratosphere (LS) O-3 profiles. Models with the greatest tropical transport problems compare poorly with O-3 observations. Models with the most realistic LS transport agree more closely with LS observations and each other. We incorporate the results of the chemistry evaluations in the Stratospheric Processes and their Role in Climate (SPARC) CCMVal Report to explain the range of CCM predictions for the return-to-1980 dates for global (60 degrees S-60 degrees N) and Antarctic column ozone. Antarctic O-3 return dates are generally correlated with vortex Cl-y levels, and vortex Cly is generally correlated with the model's circulation, although model Cl chemistry and conservation problems also have a significant effect on return date. In both regions, models with good LS transport and chemistry produce a smaller range of predictions for the return-to-1980 ozone values. This study suggests that the current range of predicted return dates is unnecessarily broad due to identifiable model deficiencies.

Monthly averaged surface erythemal solar irradiance (UV-Ery) for local noon from 1960 to 2100 has been derived using radiative transfer calculations and projections of ozone, temperature and cloud change from 14 chemistry climate models (CCM), as part of the CCMVal-2 activity of SPARC. Our calculations show the influence of ozone depletion and recovery on erythemal irradiance. In addition, we investigate UV-Ery changes caused by climate change due to increasing greenhouse gas concentrations. The latter include effects of both stratospheric ozone and cloud changes. The derived estimates provide a global picture of the likely changes in erythemal irradiance during the 21st century. Uncertainties arise from the assumed scenarios, different parameterizations – particularly of cloud effects on UV-Ery – and the spread in the CCM projections. The calculations suggest that relative to 1980, annually mean UV-Ery in the 2090s will be on average ~12 % lower at high latitudes in both hemispheres, ~3 % lower at mid latitudes, and marginally higher (~1 %) in the tropics. The largest reduction (~16 %) is projected for Antarctica in October. Cloud effects are responsible for 2–3 % of the reduction in UV-Ery at high latitudes, but they slightly moderate it at mid-latitudes (~1 %). The year of return of erythemal irradiance to values of certain milestones (1965 and 1980) depends largely on the return of column ozone to the corresponding levels and is associated with large uncertainties mainly due to the spread of the model projections. The inclusion of cloud effects in the calculations has only a small effect of the return years. At mid and high latitudes, changes in clouds and stratospheric ozone transport by global circulation changes due to greenhouse gases will sustain the erythemal irradiance at levels below those in 1965, despite the removal of ozone depleting substances. At northern high latitudes (60°–90°), the projected decreases in cloud transmittance towards the end of the 21st century will reduce the yearly average surface erythemal irradiance by ~5 % with respect to the 1960s.

The response of stratospheric climate and circulation to increasing amounts of greenhouse gases (GHGs) and ozone recovery in the twenty-first century is analyzed in simulations of 11 chemistry-climate models using near-identical forcings and experimental setup. In addition to an overall global cooling of the stratosphere in the simulations (0.59 +/- 6 0.07 K decade(-1) at 10 hPa), ozone recovery causes a warming of the Southern Hemisphere polar lower stratosphere in summer with enhanced cooling above. The rate of warming correlates with the rate of ozone recovery projected by the models and, on average, changes from 0.8 to 0.48 K decade(-1) at 100 hPa as the rate of recovery declines from the first to the second half of the century. In the winter northern polar lower stratosphere the increased radiative cooling from the growing abundance of GHGs is, in most models, balanced by adiabatic warming from stronger polar downwelling. In the Antarctic lower stratosphere the models simulate an increase in low temperature extremes required for polar stratospheric cloud (PSC) formation, but the positive trend is decreasing over the twenty-first century in all models. In the Arctic, none of the models simulates a statistically significant increase in Arctic PSCs throughout the twenty-first century. The subtropical jets accelerate in response to climate change and the ozone recovery produces a westward acceleration of the lower-stratospheric wind over the Antarctic during summer, though this response is sensitive to the rate of recovery projected by the models. There is a strengthening of the Brewer-Dobson circulation throughout the depth of the stratosphere, which reduces the mean age of air nearly everywhere at a rate of about 0.05 yr decade(-1) in those models with this diagnostic. On average, the annual mean tropical upwelling in the lower stratosphere (similar to 70 hPa) increases by almost 2% decade(-1), with 59% of this trend forced by the parameterized orographic gravity wave drag in the models. This is a consequence of the eastward acceleration of the subtropical jets, which increases the upward flux of (parameterized) momentum reaching the lower stratosphere in these latitudes.

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an "ensemble of opportunity" of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21st century, up-to and after the time when ozone concentrations return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

Simulations of 15 coupled chemistry climate models, for the period 1960–2100, are presented. The models include a detailed stratosphere, as well as including a realistic representation of the tropospheric climate. The simulations assume a consistent set of changing greenhouse gas concentrations, as well as temporally varying chlorofluorocarbon concentrations in accordance with observations for the past and expectations for the future. The ozone results are analyzed using a nonparametric additive statistical model. Comparisons are made with observations for the recent past, and the recovery of ozone, indicated by a return to 1960 and 1980 values, is investigated as a function of latitude. Although chlorine amounts are simulated to return to 1980 values by about 2050, with only weak latitudinal variations, column ozone amounts recover at different rates due to the influence of greenhouse gas changes. In the tropics, simulated peak ozone amounts occur by about 2050 and thereafter total ozone column declines. Consequently, simulated ozone does not recover to values which existed prior to the early 1980s. The results also show a distinct hemispheric asymmetry, with recovery to 1980 values in the Northern Hemisphere extratropics ahead of the chlorine return by about 20 years. In the Southern Hemisphere midlatitudes, ozone is simulated to return to 1980 levels only 10 years ahead of chlorine. In the Antarctic, annually averaged ozone recovers at about the same rate as chlorine in high latitudes and hence does not return to 1960s values until the last decade of the simulations.

The impact of stratospheric ozone on the tropospheric general circulation of the Southern Hemisphere (SH) is examined with a set of chemistry-climate models participating in the Stratospheric Processes and their Role in Climate (SPARC)/Chemistry-Climate Model Validation project phase 2 (CCMVal-2). Model integrations of both the past and future climates reveal the crucial role of stratospheric ozone in driving SH circulation change: stronger ozone depletion in late spring generally leads to greater poleward displacement and intensification of the tropospheric midlatitude jet, and greater expansion of the SH Hadley cell in the summer. These circulation changes are systematic as poleward displacement of the jet is typically accompanied by intensification of the jet and expansion of the Hadley cell. Overall results are compared with coupled models participating in the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), and possible mechanisms are discussed. While the tropospheric circulation response appears quasi-linearly related to stratospheric ozone changes, the quantitative response to a given forcing varies considerably from one model to another. This scatter partly results from differences in model climatology. It is shown that poleward intensification of the westerly jet is generally stronger in models whose climatological jet is biased toward lower latitudes. This result is discussed in the context of quasi-geostrophic zonal mean dynamics.