Land surface feedbacks and palaeomonsoons in northern Africa

We ran a sequence of climate model experiments for 6000 years ago, with land-surface conditions based on a realistic map of palaeovegetation, lakes and wetlands, to quantify the effects of land-surface feedbacks in the Saharan region. Vegetation-induced albedo and moisture flux changes produced year-round warming, forced the monsoon to 17°–25°N two months earlier, and shifted the precipitation belt ≈300 km northwards compared to the effects of orbital forcing alone. The addition of lakes and wetlands produced localised changes in evaporation and precipitation, but caused no further extension of the monsoon belt. Diagnostic analyses with biome and continental hydrology models showed that the combined land-surface feedbacks, although substantial, could neither maintain grassland as far north as observed (≈26°N) nor maintain Lake “MegaChad” (330,000 km²).

The ocean's role in polar climate change: asymmetric Arctic and Antarctic responses to greenhouse gas and ozone forcing

In recent decades, the Arctic has been warming and sea ice disappearing. By contrast, the Southern Ocean around Antarctica has been (mainly) cooling and sea-ice extent growing. We argue here that interhemispheric asymmetries in the mean ocean circulation, with sinking in the northern North Atlantic and upwelling around Antarctica, strongly influence the sea-surface temperature (SST) response to anthropogenic greenhouse gas (GHG) forcing, accelerating warming in the Arctic while delaying it in the Antarctic. Furthermore, while the amplitude of GHG forcing has been similar at the poles, significant ozone depletion only occurs over Antarctica. We suggest that the initial response of SST around Antarctica to ozone depletion is one of cooling and only later adds to the GHG-induced warming trend as upwelling of sub-surface warm water associated with stronger surface westerlies impacts surface properties. We organize our discussion around ‘climate response functions’ (CRFs), i.e. the response of the climate to ‘step’ changes in anthropogenic forcing in which GHG and/or ozone-hole forcing is abruptly turned on and the transient response of the climate revealed and studied. Convolutions of known or postulated GHG and ozone-hole forcing functions with their respective CRFs then yield the transient forced SST response (implied by linear response theory), providing a context for discussion of the differing warming/cooling trends in the Arctic and Antarctic. We speculate that the period through which we are now passing may be one in which the delayed warming of SST associated with GHG forcing around Antarctica is largely cancelled by the cooling effects associated with the ozone hole. By mid-century, however, ozone-hole effects may instead be adding to GHG warming around Antarctica but with diminished amplitude as the ozone hole heals. The Arctic, meanwhile, responding to GHG forcing but in a manner amplified by ocean heat transport, may continue to warm at an accelerating rate.

The importance of wind and buoyancy forcing for the boundary density variations and the geostrophic component of the AMOC at 26°N

It is widely thought that changes in both the surface buoyancy fluxes and wind stress drive variability in the Atlantic meridional overturning circulation (AMOC), but that they drive variability on different time scales. For example, wind forcing dominates short-term variability through its effects on Ekman currents and coastal upwelling, whereas buoyancy forcing is important for longer time scales (multiannual and decadal). However, the role of the wind forcing on multiannual to decadal time scales is less clear. Here the authors present an analysis of simulations with the Nucleus for European Modelling of the Ocean (NEMO) ocean model with the aim of explaining the important drivers of the zonal density gradient at 26°N, which is directly related to the AMOC. In the experiments, only one of either the wind stress or the buoyancy forcing is allowed to vary in time, whereas the other remains at its seasonally varying climatology. On subannual time scales, variations in the density gradient, and in the AMOC minus Ekman, are driven largely by local wind-forced coastal upwelling at both the western and eastern boundaries. On decadal time scales, buoyancy forcing related to the North Atlantic Oscillation dominates variability in the AMOC. Interestingly, however, it is found that wind forcing also plays a role at longer time scales, primarily impacting the interannual variability through the excitation of Rossby waves in the central Atlantic, which propagate westward to interact with the western boundary, but also by modulating the decadal time-scale response to buoyancy forcing.

Evaluating adjusted forcing and model spread for historical and future scenarios in the CMIP5 generation of climate models

We utilize energy budget diagnostics from the Coupled Model Intercomparison Project phase 5 (CMIP5) to evaluate the models' climate forcing since preindustrial times employing an established regression technique. The climate forcing evaluated this way, termed the adjusted forcing (AF), includes a rapid adjustment term associated with cloud changes and other tropospheric and land-surface changes. We estimate a 2010 total anthropogenic and natural AF from CMIP5 models of 1.9 ± 0.9 W m−2 (5–95% range). The projected AF of the Representative Concentration Pathway simulations are lower than their expected radiative forcing (RF) in 2095 but agree well with efficacy weighted forcings from integrated assessment models. The smaller AF, compared to RF, is likely due to cloud adjustment. Multimodel time series of temperature change and AF from 1850 to 2100 have large intermodel spreads throughout the period. The intermodel spread of temperature change is principally driven by forcing differences in the present day and climate feedback differences in 2095, although forcing differences are still important for model spread at 2095. We find no significant relationship between the equilibrium climate sensitivity (ECS) of a model and its 2003 AF, in contrast to that found in older models where higher ECS models generally had less forcing. Given the large present-day model spread, there is no indication of any tendency by modelling groups to adjust their aerosol forcing in order to produce observed trends. Instead, some CMIP5 models have a relatively large positive forcing and overestimate the observed temperature change.

Regional and monthly and clear-sky aerosol direct radiative effect (and forcing) derived from the GlobAEROSOL-AATSR satellite aerosol product

Using the GlobAEROSOL-AATSR dataset, estimates of the instantaneous, clear-sky, direct aerosol radiative effect and radiative forcing have been produced for the year 2006. Aerosol Robotic Network sun-photometer measurements have been used to characterise the random and systematic error in the GlobAEROSOL product for 22 regions covering the globe. Representative aerosol properties for each region were derived from the results of a wide range of literature sources and, along with the de-biased GlobAEROSOL AODs, were used to drive an offline version of the Met Office unified model radiation scheme. In addition to the mean AOD, best-estimate run of the radiation scheme, a range of additional calculations were done to propagate uncertainty estimates in the AOD, optical properties, surface albedo and errors due to the temporal and spatial averaging of the AOD fields. This analysis produced monthly, regional estimates of the clear-sky aerosol radiative effect and its uncertainty, which were combined to produce annual, global mean values of (−6.7±3.9)Wm−2 at the top of atmosphere (TOA) and (−12±6)Wm−2 at the surface. These results were then used to give estimates of regional, clear-sky aerosol direct radiative forcing, using modelled pre-industrial AOD fields for the year 1750 calculated for the AEROCOM PRE experiment. However, as it was not possible to quantify the uncertainty in the pre-industrial aerosol loading, these figures can only be taken as indicative and their uncertainties as lower bounds on the likely errors. Although the uncertainty on aerosol radiative effect presented here is considerably larger than most previous estimates, the explicit inclusion of the major sources of error in the calculations suggest that they are closer to the true constraint on this figure from similar methodologies, and point to the need for more, improved estimates of both global aerosol loading and aerosol optical properties.

ERA-Interim/Land: a global land surface reanalysis data set

ERA-Interim/Land is a global land surface reanalysis data set covering the period 1979–2010. It describes the evolution of soil moisture, soil temperature and snowpack. ERA-Interim/Land is the result of a single 32-year simulation with the latest ECMWF (European Centre for Medium-Range Weather Forecasts) land surface model driven by meteorological forcing from the ERA-Interim atmospheric reanalysis and precipitation adjustments based on monthly GPCP v2.1 (Global Precipitation Climatology Project). The horizontal resolution is about 80 km and the time frequency is 3-hourly. ERA-Interim/Land includes a number of parameterization improvements in the land surface scheme with respect to the original ERA-Interim data set, which makes it more suitable for climate studies involving land water resources. The quality of ERA-Interim/Land is assessed by comparing with ground-based and remote sensing observations. In particular, estimates of soil moisture, snow depth, surface albedo, turbulent latent and sensible fluxes, and river discharges are verified against a large number of site measurements. ERA-Interim/Land provides a global integrated and coherent estimate of soil moisture and snow water equivalent, which can also be used for the initialization of numerical weather prediction and climate models.

Global analysis of radiative forcing from fire-induced shortwave albedo change

Land surface albedo, a key parameter to derive Earth's surface energy balance, is used in the parameterization of numerical weather prediction, climate monitoring and climate change impact assessments. Changes in albedo due to fire have not been fully investigated on a continental and global scale. The main goal of this study, therefore, is to quantify the changes in instantaneous shortwave albedo produced by biomass burning activities and their associated radiative forcing. The study relies on the MODerate-resolution Imaging Spectroradiometer (MODIS) MCD64A1 burned-area product to create an annual composite of areas affected by fire and the MCD43C2 bidirectional reflectance distribution function (BRDF) albedo snow-free product to compute a bihemispherical reflectance time series. The approximate day of burning is used to calculate the instantaneous change in shortwave albedo. Using the corresponding National Centers for Environmental Prediction (NCEP) monthly mean downward solar radiation flux at the surface, the global radiative forcing associated with fire was computed. The analysis reveals a mean decrease in shortwave albedo of −0.014 (1σ = 0.017), causing a mean positive radiative forcing of 3.99 Wm−2 (1σ = 4.89) over the 2002–20012 time period in areas affected by fire. The greatest drop in mean shortwave albedo change occurs in 2002, which corresponds to the highest total area burned (378 Mha) observed in the same year and produces the highest mean radiative forcing (4.5 Wm−2). Africa is the main contributor in terms of burned area, but forests globally give the highest radiative forcing per unit area and thus give detectable changes in shortwave albedo. The global mean radiative forcing for the whole period studied (~0.0275 Wm−2) shows that the contribution of fires to the Earth system is not insignificant.

Sensitivity of terrestrial precipitation trends to the structural evolution of sea surface temperatures

Pronounced intermodel differences in the projected response of land surface precipitation (LSP) to future anthropogenic forcing remain in the Coupled Model Intercomparison Project Phase 5 model integrations. A large fraction of the intermodel spread in projected LSP trends is demonstrated here to be associated with systematic differences in simulated sea surface temperature (SST) trends, especially the representation of changes in (i) the interhemispheric SST gradient and (ii) the tropical Pacific SSTs. By contrast, intermodel differences in global mean SST, representative of differing global climate sensitivities, exert limited systematic influence on LSP patterns. These results highlight the importance to regional terrestrial precipitation changes of properly simulating the spatial distribution of large-scale, remote changes as reflected in the SST response to increasing greenhouse gases. Moreover, they provide guidance regarding which region-specific precipitation projections may be potentially better constrained for use in climate change impact assessments.

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.

Direct and ozone-mediated forcing of the Southern Annular Mode by greenhouse gases

We assess the roles of long-lived greenhouse gases and ozone depletion in driving meridional surface pressure gradients in the southern extratropics; these gradients are a defining feature of the Southern Annular Mode. Stratospheric ozone depletion is thought to have caused a strengthening of this mode during summer, with increasing long-lived greenhouse gases playing a secondary role. Using a coupled atmosphere-ocean chemistry-climate model, we show that there is cancelation between the direct, radiative effect of increasing greenhouse gases by the also substantial indirect—chemical and dynamical—feedbacks that greenhouse gases have via their impact on ozone. This sensitivity of the mode to greenhouse gas-induced ozone changes suggests that a consistent implementation of ozone changes due to long-lived greenhouse gases in climate models benefits the simulation of this important aspect of Southern Hemisphere climate.

Southern Hemisphere jet latitude biases in CMIP5 models linked to shortwave cloud forcing

Substantial biases in shortwave cloud forcing (SWCF) of up to ±30 W m−2are found in the midlatitudes of the Southern Hemisphere in the historical simulations of 34 CMIP5 coupled general circulation models. The SWCF biases are shown to induce surface temperature anomalies localized in the midlatitudes, and are significantly correlated with the mean latitude of the eddy-driven jet, with a negative SWCF bias corresponding to an equatorward jet latitude bias. Aquaplanet model experiments are performed to demonstrate that the jet latitude biases are primarily induced by the midlatitude SWCF anomalies, such that the jet moves toward (away from) regions of enhanced (reduced) temperature gradients. The results underline the necessity of accurately representing cloud radiative forcings in state-of-the-art coupled models.

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.

Global climatology of surface water temperatures of large lakes by remote sensing

Lake surface water temperatures (LSWTs) of 246 globally distributed large lakes were derived from Along-Track Scanning Radiometers (ATSR) for the period 1991–2011. The climatological cycles of mean LSWT derived from these data quantify on a global scale the responses of large lakes' surface temperatures to the annual cycle of forcing by solar radiation and the ambient meteorological conditions. LSWT cycles reflect the twice annual peak in net solar radiation for lakes between 1°S to 12°N. For lakes without a lake-mean seasonal ice cover, LSWT extremes exceed air temperatures by 0.5–1.7 °C for maximum and 0.7–1.9 °C for minimum temperature. The summer maximum LSWTs of lakes from 25°S to 35°N show a linear decrease with increasing altitude; −3.76 ± 0.17 °C km−1 (inline image = 0.95), marginally lower than the corresponding air temperature decrease with altitude −4.15 ± 0.24 °C km−1 (inline image = 0.95). Lake altitude of tropical lakes account for 0.78–0.83 (inline image) of the variation in the March to June LSWT–air temperature differences, with differences decreasing by 1.9 °C as the altitude increases from 500 to 1800 m above sea level (a.s.l.) We define an ‘open water phase’ as the length of time the lake-mean LSWT remains above 4 °C. There is a strong global correlation between the start and end of the lake-mean open water phase and the spring and fall 0 °C air temperature transition days, (inline image = 0.74 and 0.80, respectively), allowing for a good estimation of timing and length of the open water phase of lakes without LSWT observations. Lake depth, lake altitude and distance from coast further explain some of the inter-lake variation in the start and end of the open water phase.

Contrasting fast precipitation responses to tropospheric and stratospheric ozone forcing

The precipitation response to radiative forcing (RF) can be decomposed into a fast precipitation response (FPR), which depends on the atmospheric component of RF, and a slow response, which depends on surface temperature change. We present the first detailed climate model study of the FPR due to tropospheric and stratospheric ozone changes. The FPR depends strongly on the altitude of ozone change. Increases below about 3 km cause a positive FPR; increases above cause a negative FPR. The FPR due to stratospheric ozone change is, per unit RF, about 3 times larger than that due to tropospheric ozone. As historical ozone trends in the troposphere and stratosphere are opposite in sign, so too are the FPRs. Simple climate model calculations of the time-dependent total (fast and slow) precipitation change, indicate that ozone's contribution to precipitation change in 2011, compared to 1765, could exceed 50% of that due to CO2 change.

The impact of variable sea ice roughness on changes in Arctic Ocean surface stress: a model study

The Arctic sea ice cover is thinning and retreating, causing changes in surface roughness that in turn modify the momentum flux from the atmosphere through the ice into the ocean. New model simulations comprising variable sea ice drag coefficients for both the air and water interface demonstrate that the heterogeneity in sea ice surface roughness significantly impacts the spatial distribution and trends of ocean surface stress during the last decades. Simulations with constant sea ice drag coefficients as used in most climate models show an increase in annual mean ocean surface stress (0.003 N/m2 per decade, 4.6%) due to the reduction of ice thickness leading to a weakening of the ice and accelerated ice drift. In contrast, with variable drag coefficients our simulations show annual mean ocean surface stress is declining at a rate of -0.002 N/m2 per decade (3.1%) over the period 1980-2013 because of a significant reduction in surface roughness associated with an increasingly thinner and younger sea ice cover. The effectiveness of sea ice in transferring momentum does not only depend on its resistive strength against the wind forcing but is also set by its top and bottom surface roughness varying with ice types and ice conditions. This reveals the need to account for sea ice surface roughness variations in climate simulations in order to correctly represent the implications of sea ice loss under global warming.