Uncertainty in the relationship between climate forcing and hydrological response in UK catchments

This paper assesses the relationship between amount of climate forcing – as indexed by global mean temperature change – and hydrological response in a sample of UK catchments. It constructs climate scenarios representing different changes in global mean temperature from an ensemble of 21 climate models assessed in the IPCC AR4. The results show a considerable range in impact between the 21 climate models, with – for example - change in summer runoff at a 2oC increase in global mean temperature varying between -40% and +20%. There is evidence of clustering in the results, particularly in projected changes in summer runoff and indicators of low flows, implying that the ensemble mean is not an appropriate generalised indicator of impact, and that the standard deviation of responses does not adequately characterise uncertainty. The uncertainty in hydrological impact is therefore best characterised by considering the shape of the distribution of responses across multiple climate scenarios. For some climate model patterns, and some catchments, there is also evidence that linear climate change forcings produce non-linear hydrological impacts. For most variables and catchments, the effects of climate change are apparent above the effects of natural multi-decadal variability with an increase in global mean temperature above 1oC, but there are differences between catchments. Based on the scenarios represented in the ensemble, the effect of climate change in northern upland catchments will be seen soonest in indicators of high flows, but in southern catchments effects will be apparent soonest in measures of summer and low flows. The uncertainty in response between different climate model patterns is considerably greater than the range due to uncertainty in hydrological model parameterisation.

A comparative analysis of projected impacts of climate change on river runoff from global and catchment-scale hydrological models

We present a comparative analysis of projected impacts of climate change on river runoff from two types of distributed hydrological model, a global hydrological model (GHM) and catchment-scale hydrological models (CHM). Analyses are conducted for six catchments that are global in coverage and feature strong contrasts in spatial scale as well as climatic and development conditions. These include the Liard (Canada), Mekong (SE Asia), Okavango (SW Africa), Rio Grande (Brazil), Xiangu (China) and Harper's Brook (UK). A single GHM (Mac-PDM.09) is applied to all catchments whilst different CHMs are applied for each catchment. The CHMs typically simulate water resources impacts based on a more explicit representation of catchment water resources than that available from the GHM, and the CHMs include river routing. Simulations of average annual runoff, mean monthly runoff and high (Q5) and low (Q95) monthly runoff under baseline (1961-1990) and climate change scenarios are presented. We compare the simulated runoff response of each hydrological model to (1) prescribed increases in global mean temperature from the HadCM3 climate model and (2)a prescribed increase in global-mean temperature of 2oC for seven GCMs to explore response to climate model and structural uncertainty. We find that differences in projected changes of mean annual runoff between the two types of hydrological model can be substantial for a given GCM, and they are generally larger for indicators of high and low flow. However, they are relatively small in comparison to the range of projections across the seven GCMs. Hence, for the six catchments and seven GCMs we considered, climate model structural uncertainty is greater than the uncertainty associated with the type of hydrological model applied. Moreover, shifts in the seasonal cycle of runoff with climate change are presented similarly by both hydrological models, although for some catchments the monthly timing of high and low flows differs.This implies that for studies that seek to quantify and assess the role of climate model uncertainty on catchment-scale runoff, it may be equally as feasible to apply a GHM as it is to apply a CHM, especially when climate modelling uncertainty across the range of available GCMs is as large as it currently is. Whilst the GHM is able to represent the broad climate change signal that is represented by the CHMs, we find, however, that for some catchments there are differences between GHMs and CHMs in mean annual runoff due to differences in potential evaporation estimation methods, in the representation of the seasonality of runoff, and in the magnitude of changes in extreme monthly runoff, all of which have implications for future water management issues.

Bistability of the Atlantic overturning circulation in a global climate model and links to ocean freshwater transport

The possibility of a rapid collapse in the strength of the Atlantic meridional overturning circulation (AMOC), with associated impacts on climate, has long been recognized. The suggested basis for this risk is the existence of two stable regimes of the AMOC (‘on’ and ‘off’), and such bistable behaviour has been identified in a range of simplified climate models. However, up to now, no state-of-the-art atmosphere-ocean coupled global climate model (AOGCM) has exhibited such behaviour, leading to the interpretation that the AMOC is more stable than simpler models indicate. Here we demonstrate AMOC bistability in the response to freshwater perturbations in the FAMOUS AOGCM - the most complex AOGCM to exhibit such behaviour to date. The results also support recent suggestions that the direction of the net freshwater transport at the southern boundary of the Atlantic by the AMOC may be a useful physical indicator of the existence of bistability. We also present new estimates for this net freshwater transport by the AMOC from a range of ocean reanalyses which suggest that the Atlantic AMOC is currently in a bistable regime, although with large uncertainties. More accurate observational constraints, and an improved physical understanding of this quantity, could help narrow uncertainty in the future evolution of the AMOC and to assess the risk of a rapid AMOC collapse.

Regime dependent changes in global precipitation

Assessment of changes in precipitation (P) as a function of percentiles of surface temperature (T) and 500 hPa vertical velocity (ω) are presented, considering present-day simulations and observational estimates from the Global Precipitation Climatology Project (GPCP) combined with the European Centre for Medium-range Weather Forecasts Interim reanalysis (ERA Interim). There is a tendency for models to overestimate P in the warm, subsiding regimes compared to GPCP, in some cases by more than 100%, while many models underestimate P in the moderate temperature regimes. Considering climate change projections between 1980–1999 and 2080–2099, responses in P are characterised by dP/dT ≥ 4%/K over the coldest 10–20% of land points and over warm, ascending ocean points while P declines over the warmest, descending regimes (dP/dT ∼ − 4%/K for model ensemble means). The reduced Walker circulation limits this contrasting dP/dT response in the tropical wet and dry regimes only marginally. Around 70% of the global surface area exhibits a consistent sign for dP/dT in at least 6 out of a 7-member model ensemble when considering P composites in terms of dynamic regime.

Spin-up and adjustment of the Antarctic Circumpolar Current and global pycnocline

A theory is presented for the adjustment of the Antarctic Circumpolar Current (ACC) and global pycnocline to a sudden and sustained change in wind forcing. The adjustment timescale is controlled by the mesoscale eddy diffusivity across the ACC, the mean width of the ACC, the surface area of the ocean basins to the north, and deep water formation in the North Atlantic. In particular, northern sinking may have the potential to shorten the timescale and reduce its sensitivity to Southern Ocean eddies, but the relative importance of northern sinking and Southern Ocean eddies cannot be determined precisely, largely due to limitations in the parameterization of northern sinking. Although it is clear that the main processes that control the adjustment timescale are those which counteract the deepening of the global pycnocline, the theory also suggests that the timescale can be subtly modified by wind forcing over the ACC and global diapycnal mixing. Results from calculations with a reduced-gravity model compare well with the theory. The multidecadal-centennial adjustment timescale implies that long observational time series will be required to detect dynamic change in the ACC due to anthropogenic forcing. The potential role of Southern Ocean mesoscale eddy activity in determining both the equilibrium state of the ACC and the timescale over which it adjusts suggests that the response to anthropogenic forcing may be different in coupled ocean-atmosphere climate models that parameterize and resolve mesoscale eddies.

Kinetic energy analysis of the response of the Atlantic meridional overturning circulation to CO2-forced climate change

Atmosphere–ocean general circulation models (AOGCMs) predict a weakening of the Atlantic meridional overturning circulation (AMOC) in response to anthropogenic forcing of climate, but there is a large model uncertainty in the magnitude of the predicted change. The weakening of the AMOC is generally understood to be the result of increased buoyancy input to the north Atlantic in a warmer climate, leading to reduced convection and deep water formation. Consistent with this idea, model analyses have shown empirical relationships between the AMOC and the meridional density gradient, but this link is not direct because the large-scale ocean circulation is essentially geostrophic, making currents and pressure gradients orthogonal. Analysis of the budget of kinetic energy (KE) instead of momentum has the advantage of excluding the dominant geostrophic balance. Diagnosis of the KE balance of the HadCM3 AOGCM and its low-resolution version FAMOUS shows that KE is supplied to the ocean by the wind and dissipated by viscous forces in the global mean of the steady-state control climate, and the circulation does work against the pressure-gradient force, mainly in the Southern Ocean. In the Atlantic Ocean, however, the pressure-gradient force does work on the circulation, especially in the high-latitude regions of deep water formation. During CO2-forced climate change, we demonstrate a very good temporal correlation between the AMOC strength and the rate of KE generation by the pressure-gradient force in 50–70°N of the Atlantic Ocean in each of nine contemporary AOGCMs, supporting a buoyancy-driven interpretation of AMOC changes. To account for this, we describe a conceptual model, which offers an explanation of why AOGCMs with stronger overturning in the control climate tend to have a larger weakening under CO2 increase

Solar influence on global and regional climates

The literature relevant to how solar variability influences climate is vast—but much has been based on inadequate statistics and non-robust procedures. The common pitfalls are outlined in this review. The best estimates of the solar influence on the global mean air surface temperature show relatively small effects, compared with the response to anthropogenic changes (and broadly in line with their respective radiative forcings). However, the situation is more interesting when one looks at regional and season variations around the global means. In particular, recent research indicates that winters in Eurasia may have some dependence on the Sun, with more cold winters occurring when the solar activity is low. Advances in modelling ‘‘top-down’’ mechanisms, whereby stratospheric changes influence the underlying troposphere, offer promising explanations of the observed phenomena. In contrast, the suggested modulation of low-altitude clouds by galactic cosmic rays provides an increasingly inadequate explanation of observations.

What influence will future solar activity changes over the 21st century have on projected global near-surface temperature changes?

During the 20th century, solar activity increased in magnitude to a so-called grand maximum. It is probable that this high level of solar activity is at or near its end. It is of great interest whether any future reduction in solar activity could have a significant impact on climate that could partially offset the projected anthropogenic warming. Observations and reconstructions of solar activity over the last 9000 years are used as a constraint on possible future variations to produce probability distributions of total solar irradiance over the next 100 years. Using this information, with a simple climate model, we present results of the potential implications for future projections of climate on decadal to multidecadal timescales. Using one of the most recent reconstructions of historic total solar irradiance, the likely reduction in the warming by 2100 is found to be between 0.06 and 0.1 K, a very small fraction of the projected anthropogenic warming. However, if past total solar irradiance variations are larger and climate models substantially underestimate the response to solar variations, then there is a potential for a reduction in solar activity to mitigate a small proportion of the future warming, a scenario we cannot totally rule out. While the Sun is not expected to provide substantial delays in the time to reach critical temperature thresholds, any small delays it might provide are likely to be greater for lower anthropogenic emissions scenarios than for higher-emissions scenarios.

Tropical cyclone climatology in a 10-km global atmospheric GCM: toward weather-resolving climate modeling

Northern Hemisphere tropical cyclone (TC) activity is investigated in multiyear global climate simulations with theECMWFIntegrated Forecast System (IFS) at 10-km resolution forced by the observed records of sea surface temperature and sea ice. The results are compared to analogous simulationswith the 16-, 39-, and 125-km versions of the model as well as observations. In the North Atlantic, mean TC frequency in the 10-km model is comparable to the observed frequency, whereas it is too low in the other versions. While spatial distributions of the genesis and track densities improve systematically with increasing resolution, the 10-km model displays qualitatively more realistic simulation of the track density in the western subtropical North Atlantic. In the North Pacific, the TC count tends to be too high in thewest and too low in the east for all resolutions. These model errors appear to be associated with the errors in the large-scale environmental conditions that are fairly similar in this region for all model versions. The largest benefits of the 10-km simulation are the dramatically more accurate representation of the TC intensity distribution and the structure of the most intense storms. The model can generate a supertyphoon with a maximum surface wind speed of 68.4 m s21. The life cycle of an intense TC comprises intensity fluctuations that occur in apparent connection with the variations of the eyewall/rainband structure. These findings suggest that a hydrostatic model with cumulus parameterization and of high enough resolution could be efficiently used to simulate the TC intensity response (and the associated structural changes) to future climate change.

How much information is lost by using global-mean climate metrics? an example using the transport sector

Metrics are often used to compare the climate impacts of emissions from various sources, sectors or nations. These are usually based on global-mean input, and so there is the potential that important information on smaller scales is lost. Assuming a non-linear dependence of the climate impact on local surface temperature change, we explore the loss of information about regional variability that results from using global-mean input in the specific case of heterogeneous changes in ozone, methane and aerosol concentrations resulting from emissions from road traffic, aviation and shipping. Results from equilibrium simulations with two general circulation models are used. An alternative metric for capturing the regional climate impacts is investigated. We find that the application of a metric that is first calculated locally and then averaged globally captures a more complete and informative signal of climate impact than one that uses global-mean input. The loss of information when heterogeneity is ignored is largest in the case of aviation. Further investigation of the spatial distribution of temperature change indicates that although the pattern of temperature response does not closely match the pattern of the forcing, the forcing pattern still influences the response pattern on a hemispheric scale. When the short-lived transport forcing is superimposed on present-day anthropogenic CO2 forcing, the heterogeneity in the temperature response to CO2 dominates. This suggests that the importance of including regional climate impacts in global metrics depends on whether small sectors are considered in isolation or as part of the overall climate change.

Variation in the global-scale impacts of climate change on crop productivity due to climate model uncertainty and adaptation

Crop production is inherently sensitive to fluctuations in weather and climate and is expected to be impacted by climate change. To understand how this impact may vary across the globe many studies have been conducted to determine the change in yield of several crops to expected changes in climate. Changes in climate are typically derived from a single to no more than a few General Circulation Models (GCMs). This study examines the uncertainty introduced to a crop impact assessment when 14 GCMs are used to determine future climate. The General Large Area Model for annual crops (GLAM) was applied over a global domain to simulate the productivity of soybean and spring wheat under baseline climate conditions and under climate conditions consistent with the 2050s under the A1B SRES emissions scenario as simulated by 14 GCMs. Baseline yield simulations were evaluated against global country-level yield statistics to determine the model's ability to capture observed variability in production. The impact of climate change varied between crops, regions, and by GCM. The spread in yield projections due to GCM varied between no change and a reduction of 50%. Without adaptation yield response was linearly related to the magnitude of local temperature change. Therefore, impacts were greatest for countries at northernmost latitudes where warming is predicted to be greatest. However, these countries also exhibited the greatest potential for adaptation to offset yield losses by shifting the crop growing season to a cooler part of the year and/or switching crop variety to take advantage of an extended growing season. The relative magnitude of impacts as simulated by each GCM was not consistent across countries and between crops. It is important, therefore, for crop impact assessments to fully account for GCM uncertainty in estimating future climates and to be explicit about assumptions regarding adaptation.

Tree species diversity interacts with elevated CO2 to induce a greater root system response

As a consequence of land use change and the burning of fossil fuels, atmospheric concentrations of CO2 are increasing and altering the dynamics of the carbon cycle in forest ecosystems. In a number of studies using single tree species, fine root biomass has been shown to be strongly increased by elevated CO2. However, natural forests are often intimate mixtures of a number of co-occurring species. To investigate the interaction between tree mixture and elevated CO2, Alnus glutinosa, Betula pendula and Fagus sylvatica were planted in areas of single species and a three species polyculture in a free-air CO2 enrichment study (BangorFACE). The trees were exposed to ambient or elevated CO2 (580 µmol mol-1) for four years. Fine and coarse root biomass, together with fine root turnover and fine root morphological characteristics were measured. Fine root biomass, and morphology responded differentially to elevated CO2 at different soil depths in the three species when grown in monocultures. In polyculture, a greater response to elevated CO2 was observed in coarse roots to a depth of 20 cm, and fine root area index to a depth of 30 cm. Total fine root biomass was positively affected by elevated CO2 at the end of the experiment, but not by species diversity. Our data suggest that existing biogeochemical cycling models parameterised with data from species grown in monoculture may be underestimating the belowground response to global change.

Support for global climate reorganization during the "Medieval Climate Anomaly"

Widely distributed proxy records indicate that the Medieval Climate Anomaly (MCA; *900–1350 AD) was characterized by coherent shifts in large-scale Northern Hemisphere atmospheric circulation patterns. Although cooler sea surface temperatures in the central and eastern equatorial Pacific can explain some aspects of medieval circulation changes, they are not sufficient to account for other notable features, including widespread aridity through the Eurasian sub-tropics, stronger winter westerlies across the North Atlantic and Western Europe, and shifts in monsoon rainfall patterns across Africa and South Asia. We present results from a full-physics coupled climate model showing that a slight warming of the tropical Indian and western Pacific Oceans relative to the other tropical ocean basins can induce a broad range of the medieval circulation and climate changes indicated by proxy data, including many of those not explained by a cooler tropical Pacific alone. Important aspects of the results resemble those from previous simulations examining the climatic response to the rapid Indian Ocean warming during the late twentieth century, and to results from climate warming simulations—especially in indicating an expansion of the Northern Hemisphere Hadley circulation. Notably, the pattern of tropical Indo-Pacific sea surface temperature (SST) change responsible for producing the proxy-model similarity in our results agrees well with MCA-LIA SST differences obtained in a recent proxy-based climate field reconstruction. Though much remains unclear, our results indicate that the MCA was characterized by an enhanced zonal Indo-Pacific SST gradient with resulting changes in Northern Hemisphere tropical and extra-tropical circulation patterns and hydroclimate regimes, linkages that may explain the coherent regional climate shifts indicated by proxy records from across the planet. The findings provide new perspectives on the nature and possible causes of the MCA—a remarkable, yet incompletely understood episode of Late Holocene climatic change.