The impact of salinity perturbations on the future uptake of heat by the Atlantic Ocean

Anthropogenic ocean heat uptake is a key factor in determining climate change and sea-level rise. There is considerable uncertainty in projections of freshwater forcing of the ocean, with the potential to influence ocean heat uptake. We investigatethis by adding either -0.1 Sv or +0.1 Sv freshwater to the Atlantic in global climate model simulations, simultaneously imposing an atmospheric CO2 increase. The resulting changes in the Atlantic meridional overturning circulation are roughly equal and opposite (±2Sv). The impact of the perturbation on ocean heat content is more complex, although it is relatively small (~5%) compared to the total anthropogenic heat uptake. Several competing processes either accelerate or retard warming at different depths. Whilst positive freshwater perturbations cause an overall heating of the Atlantic, negative perturbations produce insignificant net changes in heat content. The processes active in our model appear robust, although their net result is likely model- and experiment-dependent.

Will Arctic sea ice thickness initialization improve seasonal forecast skill?

Arctic sea ice thickness is thought to be an important predictor of Arctic sea ice extent. However, coupled seasonal forecast systems do not generally use sea ice thickness observations in their initialization and are therefore missing a potentially important source of additional skill. To investigate how large this source is, a set of ensemble potential predictability experiments with a global climate model, initialized with and without knowledge of the sea ice thickness initial state, have been run. These experiments show that accurate knowledge of the sea ice thickness field is crucially important for sea ice concentration and extent forecasts up to 8 months ahead, especially in summer. Perturbing sea ice thickness also has a significant impact on the forecast error in Arctic 2 m temperature a few months ahead. These results suggest that advancing capabilities to observe and assimilate sea ice thickness into coupled forecast systems could significantly increase skill.

Counteracting the climate effects of volcanic eruptions using short-lived greenhouse gases

A large volcanic eruption might constitute a climate emergency, significantly altering global temperature and precipitation for several years. Major future eruptions will occur, but their size or timing cannot be predicted. We show, for the first time, that it may be possible to counteract these climate effects through deliberate emissions of short-lived greenhouse gases, dampening the abrupt impact of an eruption. We estimate an emission pathway countering a hypothetical eruption 3 times the size of Mount Pinatubo in 1991. We use a global climate model to evaluate global and regional responses to the eruption, with and without counteremissions. We then raise practical, financial, and ethical questions related to such a strategy. Unlike the more commonly discussed geoengineering to mitigate warming from long-lived greenhouse gases, designed emissions to counter temporary cooling would not have the disadvantage of needing to be sustained over long periods. Nevertheless, implementation would still face significant challenges.

The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Aerosol direct radiative effect of smoke over clouds over the southeast Atlantic Ocean from 2006 to 2009

The aerosol direct radiative effect (DRE) of African smoke was analyzed in cloud scenes over the southeast Atlantic Ocean, using Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) satellite observations and Hadley Centre Global Environmental Model version 2 (HadGEM2) climate model simulations. The observed mean DRE was about 30–35 W m−2 in August and September 2006–2009. In some years, short episodes of high-aerosol DRE can be observed, due to high-aerosol loadings, while in other years the loadings are lower but more prolonged. Climate models that use evenly distributed monthly averaged emission fields will not reproduce these high-aerosol loadings. Furthermore, the simulated monthly mean aerosol DRE in HadGEM2 is only about 6 W m−2 in August. The difference with SCIAMACHY mean observations can be partly explained by an underestimation of the aerosol absorption Ångström exponent in the ultraviolet. However, the subsequent increase of aerosol DRE simulation by about 20% is not enough to explain the observed discrepancy between simulations and observations.

Stratospheric dynamics and midlatitude jets under geoengineering with space mirrors and sulfate and titania aerosols

The impact on the dynamics of the stratosphere of three approaches to geoengineering by solar radiation management is investigated using idealized simulations of a global climate model. The approaches are geoengineering with sulfate aerosols, titania aerosols, and reduction in total solar irradiance (representing mirrors placed in space). If it were possible to use stratospheric aerosols to counterbalance the surface warming produced by a quadrupling of atmospheric carbon dioxide concentrations, tropical lower stratospheric radiative heating would drive a thermal wind response which would intensify the stratospheric polar vortices. In the Northern Hemisphere this intensification results in strong dynamical cooling of the polar stratosphere. Northern Hemisphere stratospheric sudden warming events become rare (one and two in 65 years for sulfate and titania, respectively). The intensification of the polar vortices results in a poleward shift of the tropospheric midlatitude jets in winter. The aerosol radiative heating enhances the tropical upwelling in the lower stratosphere, influencing the strength of the Brewer-Dobson circulation. In contrast, solar dimming does not produce heating of the tropical lower stratosphere, and so there is little intensification of the polar vortex and no enhanced tropical upwelling. The dynamical response to titania aerosol is qualitatively similar to the response to sulfate.

Potential role of vegetation feedback in the climate sensitivity of high-latitude regions: a case study at 6000 years B.P.

Previous climate model simulations have shown that the configuration of the Earth's orbit during the early to mid-Holocene (approximately 10–5 kyr) can account for the generally warmer-than-present conditions experienced by the high latitudes of the northern hemisphere. New simulations for 6 kyr with two atmospheric/mixed-layer ocean models (Community Climate Model, version 1, CCMl, and Global ENvironmental and Ecological Simulation of Interactive Systems, version 2, GENESIS 2) are presented here and compared with results from two previous simulations with GENESIS 1 that were obtained with and without the albedo feedback due to climate-induced poleward expansion of the boreal forest. The climate model results are summarized in the form of potential vegetation maps obtained with the global BIOME model, which facilitates visual comparisons both among models and with pollen and plant macrofossil data recording shifts of the forest-tundra boundary. A preliminary synthesis shows that the forest limit was shifted 100–200 km north in most sectors. Both CCMl and GENESIS 2 produced a shift of this magnitude. GENESIS 1 however produced too small a shift, except when the boreal forest albedo feedback was included. The feedback in this case was estimated to have amplified forest expansion by approximately 50%. The forest limit changes also show meridional patterns (greatest expansion in central Siberia and little or none in Alaska and Labrador) which have yet to be reproduced by models. Further progress in understanding of the processes involved in the response of climate and vegetation to orbital forcing will require both the deployment of coupled atmosphere-biosphere-ocean models and the development of more comprehensive observational data sets

Dominant role of greenhouse-gas forcing in the recovery of Sahel rainfall

Sahelian summer rainfall, controlled by the West African monsoon, exhibited large-amplitude multidecadal variability during the twentieth century. Particularly important was the severe drought of the 1970s and 1980s, which had widespread impacts1–6. Research into the causes of this drought has identified anthropogenic aerosol forcing3,4,7 and changes in sea surface temperatures (SSTs; refs 1,2,6,8–11) as the most important drivers. Since the 1980s, there has been some recovery of Sahel rainfall amounts2–6,11–14, although not to the pre-drought levels of the 1940s and 1950s. Here we report on experiments with the atmospheric component of a state-of-the-art global climate model to identify the causes of this recovery. Our results suggest that the direct influence of higher levels of greenhouse gases in the atmosphere was the main cause, with an additional role for changes in anthropogenic aerosol precursor emissions. We find that recent changes in SSTs, although substantial, did not have a significant impact on the recovery. The simulated response to anthropogenic greenhouse-gas and aerosol forcing is consistent with a multivariate fingerprint of the observed recovery, raising confidence in our findings. Although robust predictions are not yet possible, our results suggest that the recent recovery in Sahel rainfall amounts is most likely to be sustained or amplified in the near term.

The solsticial pause on Mars: 2 modelling and investigation of causes

The martian solsticial pause, presented in a companion paper (Lewis et al., this issue), was investigated further through a series of model runs using the UK version of the LMD/UK Mars Global Climate Model. It was found that the pause could not be adequately reproduced if radiatively active water ice clouds were omitted from the model. When clouds were used, along with a realistic time-dependent dust opacity distribution, a substantial minimum in near-surface transient eddy activity formed around solstice in both hemispheres. The net effect of the clouds in the model is, by altering the thermal structure of the atmosphere, to decrease the vertical shear of the westerly jet near the surface around solstice, and thus reduce baroclinic growth rates. A similar effect was seen under conditions of large dust loading, implying that northern midlatitude eddy activity will tend to become suppressed after a period of intense flushing storm formation around the northern cap edge. Suppression of baroclinic eddy generation by the barotropic component of the flow and via diabatic eddy dissipation were also investigated as possible mechanisms leading to the formation of the solsticial pause but were found not to make major contributions. Zonal variations in topography were found to be important, as their presence results in weakened transient eddies around winter solstice in both hemispheres, through modification of the near-surface flow. The zonal topographic asymmetry appears to be the primary reason for the weakness of eddy activity in the southern hemisphere relative to the northern hemisphere, and the ultimate cause of the solsticial pause in both hemispheres. The meridional topographic gradient was found to exert a much weaker influence on near-surface transient eddies.

Mechanisms of the negative shortwave cloud feedback in mid to high latitudes

Increases in cloud optical depth and liquid water path (LWP) are robust features of global warming model simulations in high latitudes, yielding a negative shortwave cloud feedback, but the mechanisms are still uncertain. We assess the importance of microphysical processes for the negative optical depth feedback by perturbing temperature in the microphysics schemes of two aquaplanet models, both of which have separate prognostic equations for liquid water and ice. We find that most of the LWP increase with warming is caused by a suppression of ice microphysical processes in mixed-phase clouds, resulting in reduced conversion efficiencies of liquid water to ice and precipitation. Perturbing the temperature-dependent phase partitioning of convective condensate also yields a small LWP increase. Together, the perturbations in large-scale microphysics and convective condensate partitioning explain more than two-thirds of the LWP response relative to a reference case with increased SSTs, and capture all of the vertical structure of the liquid water response. In support of these findings, we show the existence of a very robust positive relationship between monthly-mean LWP and temperature in CMIP5 models and observations in mixed-phase cloud regions only. In models, the historical LWP sensitivity to temperature is a good predictor of the forced global warming response poleward of about 45°, although models appear to overestimate the LWP response to warming compared to observations. We conclude that in climate models, the suppression of ice-phase microphysical processes that deplete cloud liquid water is a key driver of the LWP increase with warming and of the associated negative shortwave cloud feedback.

Mixed-phase cloud physics and Southern Ocean cloud feedback in climate models

Increasing optical depth poleward of 45° is a robust response to warming in global climate models. Much of this cloud optical depth increase has been hypothesized to be due to transitions from ice-dominated to liquid-dominated mixed-phase cloud. In this study, the importance of liquid-ice partitioning for the optical depth feedback is quantified for 19 Coupled Model Intercomparison Project Phase 5 models. All models show a monotonic partitioning of ice and liquid as a function of temperature, but the temperature at which ice and liquid are equally mixed (the glaciation temperature) varies by as much as 40 K across models. Models that have a higher glaciation temperature are found to have a smaller climatological liquid water path (LWP) and condensed water path and experience a larger increase in LWP as the climate warms. The ice-liquid partitioning curve of each model may be used to calculate the response of LWP to warming. It is found that the repartitioning between ice and liquid in a warming climate contributes at least 20% to 80% of the increase in LWP as the climate warms, depending on model. Intermodel differences in the climatological partitioning between ice and liquid are estimated to contribute at least 20% to the intermodel spread in the high-latitude LWP response in the mixed-phase region poleward of 45°S. It is hypothesized that a more thorough evaluation and constraint of global climate model mixed-phase cloud parameterizations and validation of the total condensate and ice-liquid apportionment against observations will yield a substantial reduction in model uncertainty in the high-latitude cloud response to warming.

Modelled and observed changes in aerosols and surface solar radiation over Europe between 1960 and 2009

Substantial changes in anthropogenic aerosols and precursor gas emissions have occurred over recent decades due to the implementation of air pollution control legislation and economic growth. The response of atmospheric aerosols to these changes and the impact on climate are poorly constrained, particularly in studies using detailed aerosol chemistry–climate models. Here we compare the HadGEM3-UKCA (Hadley Centre Global Environment Model-United Kingdom Chemistry and Aerosols) coupled chemistry–climate model for the period 1960–2009 against extensive ground-based observations of sulfate aerosol mass (1978–2009), total suspended particle matter (SPM, 1978–1998), PM10 (1997–2009), aerosol optical depth (AOD, 2000–2009), aerosol size distributions (2008–2009) and surface solar radiation (SSR, 1960–2009) over Europe. The model underestimates observed sulfate aerosol mass (normalised mean bias factor (NMBF) = −0.4), SPM (NMBF = −0.9), PM10 (NMBF = −0.2), aerosol number concentrations (N30 NMBF = −0.85; N50 NMBF = −0.65; and N100 NMBF = −0.96) and AOD (NMBF = −0.01) but slightly overpredicts SSR (NMBF = 0.02). Trends in aerosol over the observational period are well simulated by the model, with observed (simulated) changes in sulfate of −68 % (−78 %), SPM of −42 % (−20 %), PM10 of −9 % (−8 %) and AOD of −11 % (−14 %). Discrepancies in the magnitude of simulated aerosol mass do not affect the ability of the model to reproduce the observed SSR trends. The positive change in observed European SSR (5 %) during 1990–2009 ("brightening") is better reproduced by the model when aerosol radiative effects (ARE) are included (3 %), compared to simulations where ARE are excluded (0.2 %). The simulated top-of-the-atmosphere aerosol radiative forcing over Europe under all-sky conditions increased by > 3.0 W m−2 during the period 1970–2009 in response to changes in anthropogenic emissions and aerosol concentrations.

The solsticial pause on Mars: 1. A planetary wave reanalysis

Large-scale planetary waves are diagnosed from an analysis of profiles retrieved from the Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft during its scientific mapping phase. The analysis is conducted by assimilating thermal profiles and total dust opacity retrievals into a Mars global circulation model. Transient waves are largest throughout the northern hemisphere autumn, winter and spring period and almost absent during the summer. The southern hemisphere exhibits generally weaker transient wave behaviour. A striking feature of the low-altitude transient waves in the analysis is that they show a broad subsidiary minimum in amplitude centred on the winter solstice, a period when the thermal contrast between the summer hemisphere and the winter pole is strongest and baroclinic wave activity might be expected to be strong. This behaviour, here called the ‘solsticial pause,’ is present in every year of the analysis. This strong pause is under-represented in many independent model experiments, which tend to produce relatively uniform baroclinic wave activity throughout the winter. This paper documents and diagnoses the transient wave solsticial pause found in the analysis; a companion paper investigates the origin of the phenomenon in a series of model experiments.

Irreducible uncertainty in near-term climate projections

Model simulations of the next few decades are widely used in assessments of climate change impacts and as guidance for adaptation. Their non-linear nature reveals a level of irreducible uncertainty which it is important to understand and quantify, especially for projections of near-term regional climate. Here we use large idealised initial condition ensembles of the FAMOUS global climate model with a 1 %/year compound increase in CO2 levels to quantify the range of future temperatures in model-based projections. These simulations explore the role of both atmospheric and oceanic initial conditions and are the largest such ensembles to date. Short-term simulated trends in global temperature are diverse, and cooling periods are more likely to be followed by larger warming rates. The spatial pattern of near-term temperature change varies considerably, but the proportion of the surface showing a warming is more consistent. In addition, ensemble spread in inter-annual temperature declines as the climate warms, especially in the North Atlantic. Over Europe, atmospheric initial condition uncertainty can, for certain ocean initial conditions, lead to 20 year trends in winter and summer in which every location can exhibit either strong cooling or rapid warming. However, the details of the distribution are highly sensitive to the ocean initial condition chosen and particularly the state of the Atlantic meridional overturning circulation. On longer timescales, the warming signal becomes more clear and consistent amongst different initial condition ensembles. An ensemble using a range of different oceanic initial conditions produces a larger spread in temperature trends than ensembles using a single ocean initial condition for all lead times. This highlights the potential benefits from initialising climate predictions from ocean states informed by observations. These results suggest that climate projections need to be performed with many more ensemble members than at present, using a range of ocean initial conditions, if the uncertainty in near-term regional climate is to be adequately quantified.

Selective environmental stress caused by magmatic sulfur emissions from continental flood basalt eruptions

Several biotic crises during the past 300 million years have been linked to episodes of continental flood basalt volcanism, and in particular to the release of massive quantities of magmatic sulphur gas species. Flood basalt provinces were typically formed by numerous individual eruptions, each lasting years to decades. However, the environmental impact of these eruptions may have been limited by the occurrence of quiescent periods that lasted hundreds to thousands of years. Here we use a global aerosol model to quantify the sulphur-induced environmental effects of individual, decade-long flood basalt eruptions representative of the Columbia River Basalt Group, 16.5–14.5 million years ago, and the Deccan Traps, 65 million years ago. For a decade-long eruption of Deccan scale, we calculate a decadal-mean reduction in global surface temperature of 4.5 K, which would recover within 50 years after an eruption ceased unless climate feedbacks were very different in deep-time climates. Acid mists and fogs could have caused immediate damage to vegetation in some regions, but acid-sensitive land and marine ecosystems were well-buffered against volcanic sulphur deposition effects even during century-long eruptions. We conclude that magmatic sulphur from flood basalt eruptions would have caused a biotic crisis only if eruption frequencies and lava discharge rates had been high and sustained for several centuries at a time.