Uncertainty in simulating wheat yields under climate change

Projections of climate change impacts on crop yields are inherently uncertain1. Uncertainty is often quantified when projecting future greenhouse gas emissions and their influence on climate2. However, multi-model uncertainty analysis of crop responses to climate change is rare because systematic and objective comparisons among process-based crop simulation models1, 3 are difficult4. Here we present the largest standardized model intercomparison for climate change impacts so far. We found that individual crop models are able to simulate measured wheat grain yields accurately under a range of environments, particularly if the input information is sufficient. However, simulated climate change impacts vary across models owing to differences in model structures and parameter values. A greater proportion of the uncertainty in climate change impact projections was due to variations among crop models than to variations among downscaled general circulation models. Uncertainties in simulated impacts increased with CO2 concentrations and associated warming. These impact uncertainties can be reduced by improving temperature and CO2 relationships in models and better quantified through use of multi-model ensembles. Less uncertainty in describing how climate change may affect agricultural productivity will aid adaptation strategy development andpolicymaking.

Tackling regional climate change by leaf albedo bio-geoengineering

The likelihood that continuing greenhouse-gas emissions will lead to an unmanageable degree of climate change [1] has stimulated the search for planetary-scale technological solutions for reducing global warming [2] (“geoengineering”), typically characterized by the necessity for costly new infrastructures and industries [3]. We suggest that the existing global infrastructure associated with arable agriculture can help, given that crop plants exert an important influence over the climatic energy budget 4 and 5 because of differences in their albedo (solar reflectivity) compared to soils and to natural vegetation [6]. Specifically, we propose a “bio-geoengineering” approach to mitigate surface warming, in which crop varieties having specific leaf glossiness and/or canopy morphological traits are specifically chosen to maximize solar reflectivity. We quantify this by modifying the canopy albedo of vegetation in prescribed cropland areas in a global-climate model, and thereby estimate the near-term potential for bio-geoengineering to be a summertime cooling of more than 1°C throughout much of central North America and midlatitude Eurasia, equivalent to seasonally offsetting approximately one-fifth of regional warming due to doubling of atmospheric CO2[7]. Ultimately, genetic modification of plant leaf waxes or canopy structure could achieve greater temperature reductions, although better characterization of existing intraspecies variability is needed first.

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.

Centennial-scale climate change in Ireland during the Holocene

We examine mid- to late Holocene centennial-scale climate variability in Ireland using proxy data from peatlands, lakes and a speleothem. A high degree of between-record variability is apparent in the proxy data and significant chronological uncertainties are present. However, tephra layers provide a robust tool for correlation and improve the chronological precision of the records. Although we can find no statistically significant coherence in the dataset as a whole, a selection of high-quality peatland water table reconstructions co-vary more than would be expected by chance alone. A locally weighted regression model with bootstrapping can be used to construct a ‘best-estimate’ palaeoclimatic reconstruction from these datasets. Visual comparison and cross-wavelet analysis of peatland water table compilations from Ireland and Northern Britain show that there are some periods of coherence between these records. Some terrestrial palaeoclimatic changes in Ireland appear to coincide with changes in the North Atlantic thermohaline circulation and solar activity. However, these relationships are inconsistent and may be obscured by chronological uncertainties. We conclude by suggesting an agenda for future Holocene climate research in Ireland.

Global vegetation and terrestrial carbon cycle changes after the last ice age

•In current models, the ecophysiological effects of CO2 create both woody thickening and terrestrial carbon uptake, as observed now, and forest cover and terrestrial carbon storage increases that took place after the last glacial maximum (LGM). Here, we aimed to assess the realism of modelled vegetation and carbon storage changes between LGM and the pre-industrial Holocene (PIH). •We applied Land Processes and eXchanges (LPX), a dynamic global vegetation model (DGVM), with lowered CO2 and LGM climate anomalies from the Palaeoclimate Modelling Intercomparison Project (PMIP II), and compared the model results with palaeodata. •Modelled global gross primary production was reduced by 27–36% and carbon storage by 550–694 Pg C compared with PIH. Comparable reductions have been estimated from stable isotopes. The modelled areal reduction of forests is broadly consistent with pollen records. Despite reduced productivity and biomass, tropical forests accounted for a greater proportion of modelled land carbon storage at LGM (28–32%) than at PIH (25%). •The agreement between palaeodata and model results for LGM is consistent with the hypothesis that the ecophysiological effects of CO2 influence tree–grass competition and vegetation productivity, and suggests that these effects are also at work today.

Evaluation of a photosynthesis-based biogenic isoprene emission scheme in JULES and simulation of isoprene emissions under present-day climate conditions

We have incorporated a semi-mechanistic isoprene emission module into the JULES land-surface scheme, as a first step towards a modelling tool that can be applied for studies of vegetation – atmospheric chemistry interactions, including chemistry-climate feedbacks. Here, we evaluate the coupled model against local above-canopy isoprene emission flux measurements from six flux tower sites as well as satellite-derived estimates of isoprene emission over tropical South America and east and south Asia. The model simulates diurnal variability well: correlation coefficients are significant (at the 95 % level) for all flux tower sites. The model reproduces day-to-day variability with significant correlations (at the 95 % confidence level) at four of the six flux tower sites. At the UMBS site, a complete set of seasonal observations is available for two years (2000 and 2002). The model reproduces the seasonal pattern of emission during 2002, but does less well in the year 2000. The model overestimates observed emissions at all sites, which is partially because it does not include isoprene loss through the canopy. Comparison with the satellite-derived isoprene-emission estimates suggests that the model simulates the main spatial patterns, seasonal and inter-annual variability over tropical regions. The model yields a global annual isoprene emission of 535 ± 9 TgC yr−1 during the 1990s, 78 % of which from forested areas.

Nature of vegetation and building morphology characteristics across a city: influence on shadow patterns and mean radiant temperatures in London

Vegetation and building morphology characteristics are investigated at 19 sites on a north-south LiDAR transect across the megacity of London. Local maxima of mean building height and building plan area density at the city centre are evident. Surprisingly, the mean vegetation height (zv3) is also found to be highest in the city centre. From the LiDAR data various morphological parameters are derived as well as shadow patterns. Continuous images of the effects of buildings and of buildings plus vegetationon sky view factor (Ψ) are derived. A general reduction of Ψ is found, indicating the importance of including vegetation when deriving Ψ in urban areas. The contribution of vegetation to the shadowing at ground level is higher during summer than in autumn. Using these 3D data the influence on urban climate and mean radiant temperature (T mrt ) is calculated with SOLWEIG. The results from these simulations highlight that vegetation can be most effective at reducing heat stress within dense urban environments in summer. The daytime average T mrt is found to be lowest in the densest urban environments due to shadowing; foremost from buildings but also from trees. It is clearly shown that this method could be used to quantify the influence of vegetation on T mrt within the urban environment. The results presented in this paper highlight a number of possible climate sensitive planning practices for urban areas at the local scale (i.e. 102- 5 × 103 m).

Satellite-based estimate of the direct and indirect aerosol climate forcing

The main uncertainty in anthropogenic forcing of the Earth’s climate stems from pollution aerosols, particularly their ‘‘indirect effect’’ whereby aerosols modify cloud properties. We develop a new methodology to derive a measurement-based estimate using almost exclusively information from an Earth radiation budget instrument (CERES) and a radiometer (MODIS). We derive a statistical relationship between planetary albedo and cloud properties, and, further, between the cloud properties and column aerosol concentration. Combining these relationships with a data set of satellite-derived anthropogenic aerosol fraction, we estimate an anthropogenic radiative forcing of �-0.9 ± 0.4 Wm�-2 for the aerosol direct effect and of �-0.2 ± 0.1 Wm�-2 for the cloud albedo effect. Because of uncertainties in both satellite data and the method, the uncertainty of this result is likely larger than the values given here which correspond only to the quantifiable error estimates. The results nevertheless indicate that current global climate models may overestimate the cloud albedo effect.

The sensitivity of the tropical circulation and Maritime Continent precipitation to climate model resolution

The dependence of the annual mean tropical precipitation on horizontal resolution is investigated in the atmospheric version of the Hadley Centre General Environment Model (HadGEM1). Reducing the grid spacing from about 350 km to 110 km improves the precipitation distribution in most of the tropics. In particular, characteristic dry biases over South and Southeast Asia including the Maritime Continent as well as wet biases over the western tropical oceans are reduced. The annual-mean precipitation bias is reduced by about one third over the Maritime Continent and the neighbouring ocean basins associated with it via the Walker circulation. Sensitivity experiments show that much of the improvement with resolution in the Maritime Continent region is due to the specification of better resolved surface boundary conditions (land fraction, soil and vegetation parameters) at the higher resolution. It is shown that in particular the formulation of the coastal tiling scheme may cause resolution sensitivity of the mean simulated climate. The improvement in the tropical mean precipitation in this region is not primarily associated with the better representation of orography at the higher resolution, nor with changes in the eddy transport of moisture. Sizeable sensitivity to changes in the surface fields may be one of the reasons for the large variation of the mean tropical precipitation distribution seen across climate models.

A comprehensive benchmarking system for evaluating global vegetation models

We present a benchmark system for global vegetation models. This system provides a quantitative evaluation of multiple simulated vegetation properties, including primary production; seasonal net ecosystem production; vegetation cover; composition and height; fire regime; and runoff. The benchmarks are derived from remotely sensed gridded datasets and site-based observations. The datasets allow comparisons of annual average conditions and seasonal and inter-annual variability, and they allow the impact of spatial and temporal biases in means and variability to be assessed separately. Specifically designed metrics quantify model performance for each process, and are compared to scores based on the temporal or spatial mean value of the observations and a "random" model produced by bootstrap resampling of the observations. The benchmark system is applied to three models: a simple light-use efficiency and water-balance model (the Simple Diagnostic Biosphere Model: SDBM), the Lund-Potsdam-Jena (LPJ) and Land Processes and eXchanges (LPX) dynamic global vegetation models (DGVMs). In general, the SDBM performs better than either of the DGVMs. It reproduces independent measurements of net primary production (NPP) but underestimates the amplitude of the observed CO2 seasonal cycle. The two DGVMs show little difference for most benchmarks (including the inter-annual variability in the growth rate and seasonal cycle of atmospheric CO2), but LPX represents burnt fraction demonstrably more accurately. Benchmarking also identified several weaknesses common to both DGVMs. The benchmarking system provides a quantitative approach for evaluating how adequately processes are represented in a model, identifying errors and biases, tracking improvements in performance through model development, and discriminating among models. Adoption of such a system would do much to improve confidence in terrestrial model predictions of climate change impacts and feedbacks.

Isoprene emissions and climate

Biogenic volatile organic compounds (BVOCs) play an important role in atmospheric chemistry and the carbon cycle. Isoprene is quantitatively the most important of the non-methane BVOCs (NMBVOCs), with an annual emission of about 400–600 TgC; about 90% of this is emitted by terrestrial plants. Incorporating a mechanistic treatment of isoprene emissions within land-surface schemes has recently become a focus for the modelling community, the aim being to quantify the potential magnitude of associated climate feedbacks. However, these efforts are hampered by major uncertainties about why plants emit isoprene and the relative importance of different environmental controls on isoprene emission. The availability and reliability of observations of isoprene fluxes from different types of vegetation is limited, and this also imposes constraints on model development. Nevertheless, progress is being made towards the development of mechanistic models of isoprene emission which, in conjunction with atmospheric chemistry models, will ultimately allow improved quantification of the feedbacks between the terrestrial biosphere and climate under past and future climate states.

Synergistic feedbacks between ocean and vegetation on mid- and high-latitude climates during the mid-Holocene

Simulations with the IPSL atmosphere–ocean model asynchronously coupled with the BIOME1 vegetation model show the impact of ocean and vegetation feedbacks, and their synergy, on mid- and high-latitude (>40°N) climate in response to orbitally-induced changes in mid-Holocene insolation. The atmospheric response to orbital forcing produces a +1.2 °C warming over the continents in summer and a cooling during the rest of the year. Ocean feedback reinforces the cooling in spring but counteracts the autumn and winter cooling. Vegetation feedback produces warming in all seasons, with largest changes (+1 °C) in spring. Synergy between ocean and vegetation feedbacks leads to further warming, which can be as large as the independent impact of these feedbacks. The combination of these effects causes the high northern latitudes to be warmer throughout the year in the ocean–atmosphere-vegetation simulation. Simulated vegetation changes resulting from this year-round warming are consistent with observed mid-Holocene vegetation patterns. Feedbacks also impact on precipitation. The atmospheric response to orbital-forcing reduces precipitation throughout the year; the most marked changes occur in the mid-latitudes in summer. Ocean feedback reduces aridity during autumn, winter and spring, but does not affect summer precipitation. Vegetation feedback increases spring precipitation but amplifies summer drying. Synergy between the feedbacks increases precipitation in autumn, winter and spring, and reduces precipitation in summer. The combined changes amplify the seasonal contrast in precipitation in the ocean–atmosphere-vegetation simulation. Enhanced summer drought produces an unrealistically large expansion of temperate grasslands, particularly in mid-latitude Eurasia.

The impact of sea ice dynamics on the Arctic climate system

Five paired global climate model experiments, one with an ice pack that only responds thermodynamically (TI) and one including sea-ice dynamics (DI), were used to investigate the sensitivity of Arctic climates to sea-ice motion. The sequence of experiments includes situations in which the Arctic was both considerably colder (Glacial Inception, ca 115,000 years ago) and considerably warmer (3 × CO2) than today. Sea-ice motion produces cooler anomalies year-round than simulations without ice dynamics, resulting in reduced Arctic warming in warm scenarios and increased Arctic cooling in cold scenarios. These changes reflect changes in atmospheric circulation patterns: the DI simulations favor outflow of Arctic air and sea ice into the North Atlantic by promoting cyclonic circulation centered over northern Eurasia, whereas the TI simulations favor southerly inflow of much warmer air from the North Atlantic by promoting cyclonic circulation centered over Greenland. The differences between the paired simulations are sufficiently large to produce different vegetation cover over >19% of the land area north of 55°N, resulting in changes in land-surface characteristics large enough to have an additional impact on climate. Comparison of the DI and TI experiments for the mid-Holocene (6000 years ago) with paleovegetation reconstructions suggests the incorporation of sea-ice dynamics yields a more realistic simulation of high-latitude climates. The spatial pattern of sea-ice anomalies in the warmer-than-modern DI experiments strongly resembles the observed Arctic Ocean sea-ice dipole structure in recent decades, consistent with the idea that greenhouse warming is already impacting the high-northern latitudes.

The role of dust in climate changes today, at the last glacial maximum and in the future.

Natural mineral aerosol (dust) is an active component of the climate system and plays multiple roles in mediating physical and biogeochemical exchanges between the atmosphere, land surface and ocean. Changes in the amount of dust in the atmosphere are caused both by changes in climate (precipitation, wind strength, regional moisture balance) and changes in the extent of dust sources caused by either anthropogenic or climatically induced changes in vegetation cover. Models of the global dust cycle take into account the physical controls on dust deflation from prescribed source areas (based largely on soil wetness and vegetation cover thresholds), dust transport within the atmospheric column, and dust deposition through sedimentation and scavenging by precipitation. These models successfully reproduce the first-order spatial and temporal patterns in atmospheric dust loading under modern conditions. Atmospheric dust loading was as much as an order-of-magnitude larger than today during the last glacial maximum (LGM). While the observed increase in emissions from northern Africa can be explained solely in terms of climate changes (colder, drier and windier glacial climates), increased emissions from other regions appear to have been largely a response to climatically induced changes in vegetation cover and hence in the extent of dust source areas. Model experiments suggest that the increased dust loading in tropical regions had an effect on radiative forcing comparable to that of low glacial CO2 levels. Changes in land-use are already increasing the dust loading of the atmosphere. However, simulations show that anthropogenically forced climate changes substantially reduce the extent and productivity of natural dust sources. Positive feedbacks initiated by a reduction of dust emissions from natural source areas on both radiative forcing and atmospheric CO2 could substantially mitigate the impacts of land-use changes, and need to be considered in climate change assessments.