Detection and attribution of climate change: a regional perspective

The Intergovernmental Panel on Climate Change fourth assessment report, published in 2007 came to a more confident assessment of the causes of global temperature change than previous reports and concluded that ‘it is likely that there has been significant anthropogenic warming over the past 50 years averaged over each continent except Antarctica.’ Since then, warming over Antarctica has also been attributed to human influence, and further evidence has accumulated attributing a much wider range of climate changes to human activities. Such changes are broadly consistent with theoretical understanding, and climate model simulations, of how the planet is expected to respond. This paper reviews this evidence from a regional perspective to reflect a growing interest in understanding the regional effects of climate change, which can differ markedly across the globe. We set out the methodological basis for detection and attribution and discuss the spatial scales on which it is possible to make robust attribution statements. We review the evidence showing significant human-induced changes in regional temperatures, and for the effects of external forcings on changes in the hydrological cycle, the cryosphere, circulation changes, oceanic changes, and changes in extremes. We then discuss future challenges for the science of attribution. To better assess the pace of change, and to understand more about the regional changes to which societies need to adapt, we will need to refine our understanding of the effects of external forcing and internal variability

Decadal to multidecadal variability and the climate change background

Three prominent quasi-global patterns of variability and change are observed using the Met Office's sea surface temperature (SST) analysis and almost independent night marine air temperature analysis. The first is a global warming signal that is very highly correlated with global mean SST. The second is a decadal to multidecadal fluctuation with some geographical similarity to the El Niño–Southern Oscillation (ENSO). It is associated with the Pacific Decadal Oscillation (PDO), and its Pacific-wide manifestation has been termed the Interdecadal Pacific Oscillation (IPO). We present model investigations of the relationship between the IPO and ENSO. The third mode is an interhemispheric variation on multidecadal timescales which, in view of climate model experiments, is likely to be at least partly due to natural variations in the thermohaline circulation. Observed climatic impacts of this mode also appear in model simulations. Smaller-scale, regional atmospheric phenomena also affect climate on decadal to interdecadal timescales. We concentrate on one such mode, the winter North Atlantic Oscillation (NAO). This shows strong decadal to interdecadal variability and a correspondingly strong influence on surface climate variability which is largely additional to the effects of recent regional anthropogenic climate change. The winter NAO is likely influenced by both SST forcing and stratospheric variability. A full understanding of decadal changes in the NAO and European winter climate may require a detailed representation of the stratosphere that is hitherto missing in the major climate models used to study climate change.

The determination of total nitrogen and total phosphorus concentrations in freshwaters from land use, stock headage and population data: testing of a model for use in conservation and water quality management

1. Nutrient concentrations (particularly N and P) determine the extent to which water bodies are or may become eutrophic. Direct determination of nutrient content on a wide scale is labour intensive but the main sources of N and P are well known. This paper describes and tests an export coefficient model for prediction of total N and total P from: (i) land use, stock headage and human population; (ii) the export rates of N and P from these sources; and (iii) the river discharge. Such a model might be used to forecast the effects of changes in land use in the future and to hindcast past water quality to establish comparative or baseline states for the monitoring of change. 2. The model has been calibrated against observed data for 1988 and validated against sets of observed data for a sequence of earlier years in ten British catchments varying from uplands through rolling, fertile lowlands to the flat topography of East Anglia. 3. The model predicted total N and total P concentrations with high precision (95% of the variance in observed data explained). It has been used in two forms: the first on a specific catchment basis; the second for a larger natural region which contains the catchment with the assumption that all catchments within that region will be similar. Both models gave similar results with little loss of precision in the latter case. This implies that it will be possible to describe the overall pattern of nutrient export in the UK with only a fraction of the effort needed to carry out the calculations for each individual water body. 4. Comparison between land use, stock headage, population numbers and nutrient export for the ten catchments in the pre-war year of 1931, and for 1970 and 1988 show that there has been a substantial loss of rough grazing to fertilized temporary and permanent grasslands, an increase in the hectarage devoted to arable, consistent increases in the stocking of cattle and sheep and a marked movement of humans to these rural catchments. 5. All of these trends have increased the flows of nutrients with more than a doubling of both total N and total P loads during the period. On average in these rural catchments, stock wastes have been the greatest contributors to both N and P exports, with cultivation the next most important source of N and people of P. Ratios of N to P were high in 1931 and remain little changed so that, in these catchments, phosphorus continues to be the nutrient most likely to control algal crops in standing waters supplied by the rivers studied.

Modelling the impact of land use change on water quality in agricultural catchments

Export coefficient modelling was used to model the impact of agriculture on nitrogen and phosphorus loading on the surface waters of two contrasting agricultural catchments. The model was originally developed for the Windrush catchment where the highly reactive Jurassic limestone aquifer underlying the catchment is well connected to the surface drainage network, allowing the system to be modelled using uniform export coefficients for each nutrient source in the catchment, regardless of proximity to the surface drainage network. In the Slapton catchment, the hydrological path-ways are dominated by surface and lateral shallow subsurface flow, requiring modification of the export coefficient model to incorporate a distance-decay component in the export coefficients. The modified model was calibrated against observed total nitrogen and total phosphorus loads delivered to Slapton Ley from inflowing streams in its catchment. Sensitivity analysis was conducted to isolate the key controls on nutrient export in the modified model. The model was validated against long-term records of water quality, and was found to be accurate in its predictions and sensitive to both temporal and spatial changes in agricultural practice in the catchment. The model was then used to forecast the potential reduction in nutrient loading on Slapton Ley associated with a range of catchment management strategies. The best practicable environmental option (BPEO) was found to be spatial redistribution of high nutrient export risk sources to areas of the catchment with the greatest intrinsic nutrient retention capacity.

Evaluation and management of the impact of land use change on the nitrogen and phosphorus load delivered to surface waters: the export coefficient modelling approach

A manageable, relatively inexpensive model was constructed to predict the loss of nitrogen and phosphorus from a complex catchment to its drainage system. The model used an export coefficient approach, calculating the total nitrogen (N) and total phosphorus (P) load delivered annually to a water body as the sum of the individual loads exported from each nutrient source in its catchment. The export coefficient modelling approach permits scaling up from plot-scale experiments to the catchment scale, allowing application of findings from field experimental studies at a suitable scale for catchment management. The catchment of the River Windrush, a tributary of the River Thames, UK, was selected as the initial study site. The Windrush model predicted nitrogen and phosphorus loading within 2% of observed total nitrogen load and 0.5% of observed total phosphorus load in 1989. The export coefficient modelling approach was then validated by application in a second research basin, the catchment of Slapton Ley, south Devon, which has markedly different catchment hydrology and land use. The Slapton model was calibrated within 2% of observed total nitrogen load and 2.5% of observed total phosphorus load in 1986. Both models proved sensitive to the impact of temporal changes in land use and management on water quality in both catchments, and were therefore used to evaluate the potential impact of proposed pollution control strategies on the nutrient loading delivered to the River Windrush and Slapton Ley

Steady state and dynamic modelling of nitrogen in the River Kennet: impacts of land use change since the 1930s

Steady state and dynamic models have been developed and applied to the River Kennet system. Annual nitrogen exports from the land surface to the river have been estimated based on land use from the 1930s and the 1990s. Long term modelled trends indicate that there has been a large increase in nitrogen transport into the river system driven by increased fertiliser application associated with increased cereal production, increased population and increased livestock levels. The dynamic model INCA Integrated Nitrogen in Catchments. has been applied to simulate the day-to-day transport of N from the terrestrial ecosystem to the riverine environment. This process-based model generates spatial and temporal data and reproduces the observed instream concentrations. Applying the model to current land use and 1930s land use indicates that there has been a major shift in the short term dynamics since the 1930s, with increased river and groundwater concentrations caused by both non-point source pollution from agriculture and point source discharges. �

Wildfire responses to abrupt climate change in North America

It is widely accepted, based on data from the last few decades and on model simulations, that anthropogenic climate change will cause increased fire activity. However, less attention has been paid to the relationship between abrupt climate changes and heightened fire activity in the paleorecord. We use 35 charcoal and pollen records to assess how fire regimes in North America changed during the last glacial–interglacial transition (15 to 10 ka), a time of large and rapid climate changes. We also test the hypothesis that a comet impact initiated continental-scale wildfires at 12.9 ka; the data do not support this idea, nor are continent-wide fires indicated at any time during deglaciation. There are, however, clear links between large climate changes and fire activity. Biomass burning gradually increased from the glacial period to the beginning of the Younger Dryas. Although there are changes in biomass burning during the Younger Dryas, there is no systematic trend. There is a further increase in biomass burning after the Younger Dryas. Intervals of rapid climate change at 13.9, 13.2, and 11.7 ka are marked by large increases in fire activity. The timing of changes in fire is not coincident with changes in human population density or the timing of the extinction of the megafauna. Although these factors could have contributed to fire-regime changes at individual sites or at specific times, the charcoal data indicate an important role for climate, and particularly rapid climate change, in determining broad-scale levels of fire activity.

Pollen-based continental climate reconstructions at 6 and 21 ka: a global synthesis

Subfossil pollen and plant macrofossil data derived from 14C-dated sediment profiles can provide quantitative information on glacial and interglacial climates. The data allow climate variables related to growing season warmth, winter cold, and plant-available moisture to be reconstructed. Continental-scale reconstructions have been made for the mid-Holocene (MH, around 6 ka) and Last Glacial Maximum (LGM, around 21 ka), allowing comparison with palaeoclimate simulations currently being carried out as part of the fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change. The synthesis of the available MH and LGM climate reconstructions and their uncertainties, obtained using modern analogue, regression and model-inversion techniques, is presented for four temperature variables and two moisture variables. Reconstructions of the same variables based on surface-pollen assemblages are shown to be accurate and unbiased. Reconstructed LGM and MH climate anomaly patterns are coherent, consistent between variables, and robust with respect to the choice of technique. They support a conceptual model of the controls of Late Quaternary climate change whereby the first-order effects of orbital variations and greenhouse forcing on the seasonal cycle of temperature are predictably modified by responses of the atmospheric circulation and surface energy balance.

Sensitivity of an Earth system climate model to idealized radiative forcing

We diagnose forcing and climate feedbacks in benchmark sensitivity experiments with the new Met Office Hadley Centre Earth system climate model HadGEM2-ES. To identify the impact of newly-included biogeophysical and chemical processes, results are compared to a parallel set of experiments performed with these processes switched off, and different couplings with the biogeochemistry. In abrupt carbon dioxide quadrupling experiments we find that the inclusion of these processes does not alter the global climate sensitivity of the model. However, when the change in carbon dioxide is uncoupled from the vegetation, or when the model is forced with a non-carbon dioxide forcing – an increase in solar constant – new feedbacks emerge that make the climate system less sensitive to external perturbations. We identify a strong negative dust-vegetation feedback on climate change that is small in standard carbon dioxide sensitivity experiments due to the physiological/fertilization effects of carbon dioxide on plants in this model.

Global and regional temperature-change potentials for near-term climate forcers

We examine the climate effects of the emissions of near-term climate forcers (NTCFs) from 4 continental regions (East Asia, Europe, North America and South Asia) using radiative forcing from the task force on hemispheric transport of air pollution source-receptor global chemical transport model simulations. These simulations model the transport of 3 aerosol species (sulphate, particulate organic matter and black carbon) and 4 ozone precursors (methane, nitric oxides (NOx), volatile organic compounds and carbon monoxide). From the equilibrium radiative forcing results we calculate global climate metrics, global warming potentials (GWPs) and global temperature change potentials (GTPs) and show how these depend on emission region, and can vary as functions of time. For the aerosol species, the GWP(100) values are −37±12, −46±20, and 350±200 for SO2, POM and BC respectively for the direct effects only. The corresponding GTP(100) values are −5.2±2.4, −6.5±3.5, and 50±33. This analysis is further extended by examining the temperature-change impacts in 4 latitude bands. This shows that the latitudinal pattern of the temperature response to emissions of the NTCFs does not directly follow the pattern of the diagnosed radiative forcing. For instance temperatures in the Arctic latitudes are particularly sensitive to NTCF emissions in the northern mid-latitudes. At the 100-yr time horizon the ARTPs show NOx emissions can have a warming effect in the northern mid and high latitudes, but cooling in the tropics and Southern Hemisphere. The northern mid-latitude temperature response to northern mid-latitude emissions of most NTCFs is approximately twice as large as would be implied by the global average.

The benefits of a compute cluster approach to high spatial resolution biodiversity richness modelling: projecting the impact of climate change on Mediterranean flora

High spatial resolution environmental data gives us a better understanding of the environmental factors affecting plant distributions at fine spatial scales. However, large environmental datasets dramatically increase compute times and output species model size stimulating the need for an alternative computing solution. Cluster computing offers such a solution, by allowing both multiple plant species Environmental Niche Models (ENMs) and individual tiles of high spatial resolution models to be computed concurrently on the same compute cluster. We apply our methodology to a case study of 4,209 species of Mediterranean flora (around 17% of species believed present in the biome). We demonstrate a 16 times speed-up of ENM computation time when 16 CPUs were used on the compute cluster. Our custom Java ‘Merge’ and ‘Downsize’ programs reduce ENM output files sizes by 94%. The median 0.98 test AUC score of species ENMs is aided by various species occurrence data filtering techniques. Finally, by calculating the percentage change of individual grid cell values, we map the projected percentages of plant species vulnerable to climate change in the Mediterranean region between 1950–2000 and 2020.

Real-time multi-model decadal climate predictions

We present the first climate prediction of the coming decade made with multiple models, initialized with prior observations. This prediction accrues from an international activity to exchange decadal predictions in near real-time, in order to assess differences and similarities, provide a consensus view to prevent over-confidence in forecasts from any single model, and establish current collective capability. We stress that the forecast is experimental, since the skill of the multi-model system is as yet unknown. Nevertheless, the forecast systems used here are based on models that have undergone rigorous evaluation and individually have been evaluated for forecast skill. Moreover, it is important to publish forecasts to enable open evaluation, and to provide a focus on climate change in the coming decade. Initialized forecasts of the year 2011 agree well with observations, with a pattern correlation of 0.62 compared to 0.31 for uninitialized projections. In particular, the forecast correctly predicted La Niña in the Pacific, and warm conditions in the north Atlantic and USA. A similar pattern is predicted for 2012 but with a weaker La Niña. Indices of Atlantic multi-decadal variability and Pacific decadal variability show no signal beyond climatology after 2015, while temperature in the Niño3 region is predicted to warm slightly by about 0.5 °C over the coming decade. However, uncertainties are large for individual years and initialization has little impact beyond the first 4 years in most regions. Relative to uninitialized forecasts, initialized forecasts are significantly warmer in the north Atlantic sub-polar gyre and cooler in the north Pacific throughout the decade. They are also significantly cooler in the global average and over most land and ocean regions out to several years ahead. However, in the absence of volcanic eruptions, global temperature is predicted to continue to rise, with each year from 2013 onwards having a 50 % chance of exceeding the current observed record. Verification of these forecasts will provide an important opportunity to test the performance of models and our understanding and knowledge of the drivers of climate change.

How will organic carbon stocks in mineral soils evolve under future climate? Global projections using RothC for a range of climate change scenarios

We use a soil carbon (C) model (RothC), driven by a range of climate models for a range of climate scenarios to examine the impacts of future climate on global soil organic carbon (SOC) stocks. The results suggest an overall global increase in SOC stocks by 2100 under all scenarios, but with a different extent of increase among the climate model and emissions scenarios. The impacts of projected land use changes are also simulated, but have relatively minor impacts at the global scale. Whether soils gain or lose SOC depends upon the balance between C inputs and decomposition. Changes in net primary production (NPP) change C inputs to the soil, whilst decomposition usually increases under warmer temperatures, but can also be slowed by decreased soil moisture. Underlying the global trend of increasing SOC under future climate is a complex pattern of regional SOC change. SOC losses are projected to occur in northern latitudes where higher SOC decomposition rates due to higher temperatures are not balanced by increased NPP, whereas in tropical regions, NPP increases override losses due to higher SOC decomposition. The spatial heterogeneity in the response of SOC to changing climate shows how delicately balanced the competing gain and loss processes are, with subtle changes in temperature, moisture, soil type and land use, interacting to determine whether SOC increases or decreases in the future. Our results suggest that we should stop looking for a single answer regarding whether SOC stocks will increase or decrease under future climate, since there is no single answer. Instead, we should focus on improving our prediction of the factors that determine the size and direction of change, and the land management practices that can be implemented to protect and enhance SOC stocks.

Precipitation, radiative forcing and global temperature change

Radiative forcing is a useful tool for predicting equilibrium global temperature change. However, it is not so useful for predicting global precipitation changes, as changes in precipitation strongly depend on the climate change mechanism and how it perturbs the atmospheric and surface energy budgets. Here a suite of climate model experiments and radiative transfer calculations are used to quantify and assess this dependency across a range of climate change mechanisms. It is shown that the precipitation response can be split into two parts: a fast atmospheric response that strongly correlates with the atmospheric component of radiative forcing, and a slower response to global surface temperature change that is independent of the climate change mechanism, ∼2-3% per unit of global surface temperature change. We highlight the precipitation response to black carbon aerosol forcing as falling within this range despite having an equilibrium response that is of opposite sign to the radiative forcing and global temperature change.

Updated estimate of aerosol direct radiative forcing from satellite observations and comparison against the Hadley Centre climate model

The fourth assessment report of the Intergovernmental Panel on Climate Change (IPCC) includes a comparison of observation-based and modeling-based estimates of the aerosol direct radiative forcing. In this comparison, satellite-based studies suggest a more negative aerosol direct radiative forcing than modeling studies. A previous satellite-based study, part of the IPCC comparison, uses aerosol optical depths and accumulation-mode fractions retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) at collection 4. The latest version of MODIS products, named collection 5, improves aerosol retrievals. Using these products, the direct forcing in the shortwave spectrum defined with respect to present-day natural aerosols is now estimated at −1.30 and −0.65 Wm−2 on a global clear-sky and all-sky average, respectively, for 2002. These values are still significantly more negative than the numbers reported by modeling studies. By accounting for differences between present-day natural and preindustrial aerosol concentrations, sampling biases, and investigating the impact of differences in the zonal distribution of anthropogenic aerosols, good agreement is reached between the direct forcing derived from MODIS and the Hadley Centre climate model HadGEM2-A over clear-sky oceans. Results also suggest that satellite estimates of anthropogenic aerosol optical depth over land should be coupled with a robust validation strategy in order to refine the observation-based estimate of aerosol direct radiative forcing. In addition, the complex problem of deriving the aerosol direct radiative forcing when aerosols are located above cloud still needs to be addressed.