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.

Ocean heat uptake and its consequences for the magnitude of sea level rise and climate change

Under increasing greenhouse gas concentrations, ocean heat uptake moderates the rate of climate change, and thermal expansion makes a substantial contribution to sea level rise. In this paper we quantify the differences in projections among atmosphere-ocean general circulation models of the Coupled Model Intercomparison Project in terms of transient climate response, ocean heat uptake efficiency and expansion efficiency of heat. The CMIP3 and CMIP5 ensembles have statistically indistinguishable distributions in these parameters. The ocean heat uptake efficiency varies by a factor of two across the models, explaining about 50% of the spread in ocean heat uptake in CMIP5 models with CO2 increasing at 1%/year. It correlates with the ocean global-mean vertical profiles both of temperature and of temperature change, and comparison with observations suggests the models may overestimate ocean heat uptake and underestimate surface warming, because their stratification is too weak. The models agree on the location of maxima of shallow ocean heat uptake (above 700 m) in the Southern Ocean and the North Atlantic, and on deep ocean heat uptake (below 2000 m) in areas of the Southern Ocean, in some places amounting to 40% of the top-to-bottom integral in the CMIP3 SRES A1B scenario. The Southern Ocean dominates global ocean heat uptake; consequently the eddy-induced thickness diffusivity parameter, which is particularly influential in the Southern Ocean, correlates with the ocean heat uptake efficiency. The thermal expansion produced by ocean heat uptake is 0.12 m YJ−1, with an uncertainty of about 10% (1 YJ = 1024 J).

The effect of climate change on the variability of the Northern Hemisphere stratospheric polar vortex

With extreme variability of the Arctic polar vortex being a key link for stratosphere–troposphere influences, its evolution into the twenty-first century is important for projections of changing surface climate in response to greenhouse gases. Variability of the stratospheric vortex is examined using a state-of-the-art climate model and a suite of specifically developed vortex diagnostics. The model has a fully coupled ocean and a fully resolved stratosphere. Analysis of the standard stratospheric zonal mean wind diagnostic shows no significant increase over the twenty-first century in the number of major sudden stratospheric warmings (SSWs) from its historical value of 0.7 events per decade, although the monthly distribution of SSWs does vary, with events becoming more evenly dispersed throughout the winter. However, further analyses using geometric-based vortex diagnostics show that the vortex mean state becomes weaker, and the vortex centroid is climatologically more equatorward by up to 2.5°, especially during early winter. The results using these diagnostics not only characterize the vortex structure and evolution but also emphasize the need for vortex-centric diagnostics over zonally averaged measures. Finally, vortex variability is subdivided into wave-1 (displaced) and -2 (split) components, and it is implied that vortex displacement events increase in frequency under climate change, whereas little change is observed in splitting events.

A verification framework for interannual-to-decadal predictions experiments

Decadal predictions have a high profile in the climate science community and beyond, yet very little is known about their skill. Nor is there any agreed protocol for estimating their skill. This paper proposes a sound and coordinated framework for verification of decadal hindcast experiments. The framework is illustrated for decadal hindcasts tailored to meet the requirements and specifications of CMIP5 (Coupled Model Intercomparison Project phase 5). The chosen metrics address key questions about the information content in initialized decadal hindcasts. These questions are: (1) Do the initial conditions in the hindcasts lead to more accurate predictions of the climate, compared to un-initialized climate change projections? and (2) Is the prediction model’s ensemble spread an appropriate representation of forecast uncertainty on average? The first question is addressed through deterministic metrics that compare the initialized and uninitialized hindcasts. The second question is addressed through a probabilistic metric applied to the initialized hindcasts and comparing different ways to ascribe forecast uncertainty. Verification is advocated at smoothed regional scales that can illuminate broad areas of predictability, as well as at the grid scale, since many users of the decadal prediction experiments who feed the climate data into applications or decision models will use the data at grid scale, or downscale it to even higher resolution. An overall statement on skill of CMIP5 decadal hindcasts is not the aim of this paper. The results presented are only illustrative of the framework, which would enable such studies. However, broad conclusions that are beginning to emerge from the CMIP5 results include (1) Most predictability at the interannual-to-decadal scale, relative to climatological averages, comes from external forcing, particularly for temperature; (2) though moderate, additional skill is added by the initial conditions over what is imparted by external forcing alone; however, the impact of initialization may result in overall worse predictions in some regions than provided by uninitialized climate change projections; (3) limited hindcast records and the dearth of climate-quality observational data impede our ability to quantify expected skill as well as model biases; and (4) as is common to seasonal-to-interannual model predictions, the spread of the ensemble members is not necessarily a good representation of forecast uncertainty. The authors recommend that this framework be adopted to serve as a starting point to compare prediction quality across prediction systems. The framework can provide a baseline against which future improvements can be quantified. The framework also provides guidance on the use of these model predictions, which differ in fundamental ways from the climate change projections that much of the community has become familiar with, including adjustment of mean and conditional biases, and consideration of how to best approach forecast uncertainty.

Aerosol contribution to the rapid warming of near-term climate under RCP 2.6

The importance of aerosol emissions for near term climate projections is investigated by analysing simulations with the HadGEM2-ES model under two different emissions scenarios: RCP2.6 and RCP4.5. It is shown that the near term warming projected under RCP2.6 is greater than under RCP4.5, even though the greenhouse gas forcing is lower. Rapid and substantial reductions in sulphate aerosol emissions due to a reduction of coal burning in RCP2.6 lead to a reduction in the negative shortwave forcing due to aerosol direct and indirect effects. Indirect effects play an important role over the northern hemisphere oceans, especially the subtropical northeastern Pacific where an anomaly of 5-10\,Wm$^{-2}$ develops. The pattern of surface temperature change is consistent with the expected response to this surface radiation anomaly, whilst also exhibiting features that reflect redistribution of energy, and feedbacks, within the climate system. These results demonstrate the importance of aerosol emissions as a key source of uncertainty in near term projections of global and regional climate.

Revisiting Asian monsoon formation and change associated with Tibetan Plateau forcing: II. Change

Data analysis based on station observations reveals that many meteorological variables averaged over the Tibetan Plateau (TP) are closely correlated, and their trends during the past decades are well correlated with the rainfall trend of the Asian summer monsoon. However, such correlation does not necessarily imply causality. Further diagnosis confirms the existence of a weakening trend in TP thermal forcing, characterized by weakened surface sensible heat flux in spring and summer during the past decades. This weakening trend is associated with decreasing summer precipitation over northern South Asia and North China and increasing precipitation over northwestern China, South China, and Korea. An atmospheric general circulation model, the HadAM3, is employed to elucidate the causality between the weakening TP forcing and the change in the Asian summer monsoon rainfall. Results demonstrate that a weakening in surface sensible heating over the TP results in reduced summer precipitation in the plateau region and a reduction in the associated latent heat release in summer. These changes in turn result in the weakening of the near-surface cyclonic circulation surrounding the plateau and the subtropical anticyclone over the subtropical western North Pacific, similar to the results obtained from the idealized TP experiment in Part I of this study. The southerly that normally dominates East Asia, ranging from the South China Sea to North China, weakens, resulting in a weaker equilibrated Sverdrup balance between positive vorticity generation and latent heat release. Consequently, the convergence of water vapor transport is confined to South China, forming a unique anomaly pattern in monsoon rainfall, the so-called “south wet and north dry.” Because the weakening trend in TP thermal forcing is associated with global warming, the present results provide an effective means for assessing projections of regional climate over Asia in the context of global warming.

Warming caused by cumulative carbon emissions towards the trillionth tonne

Global efforts to mitigate climate change are guided by projections of future temperatures1. But the eventual equilibrium global mean temperature associated with a given stabilization level of atmospheric greenhouse gas concentrations remains uncertain1, 2, 3, complicating the setting of stabilization targets to avoid potentially dangerous levels of global warming4, 5, 6, 7, 8. Similar problems apply to the carbon cycle: observations currently provide only a weak constraint on the response to future emissions9, 10, 11. Here we use ensemble simulations of simple climate-carbon-cycle models constrained by observations and projections from more comprehensive models to simulate the temperature response to a broad range of carbon dioxide emission pathways. We find that the peak warming caused by a given cumulative carbon dioxide emission is better constrained than the warming response to a stabilization scenario. Furthermore, the relationship between cumulative emissions and peak warming is remarkably insensitive to the emission pathway (timing of emissions or peak emission rate). Hence policy targets based on limiting cumulative emissions of carbon dioxide are likely to be more robust to scientific uncertainty than emission-rate or concentration targets. Total anthropogenic emissions of one trillion tonnes of carbon (3.67 trillion tonnes of CO2), about half of which has already been emitted since industrialization began, results in a most likely peak carbon-dioxide-induced warming of 2 °C above pre-industrial temperatures, with a 5–95% confidence interval of 1.3–3.9 °C.

A probabilistic analysis of the future potential of evaporative cooling systems in a temperate climate

The probabilistic projections of climate change for the United Kingdom (UK Climate Impacts Programme) show a trend towards hotter and drier summers. This suggests an expected increase in cooling demand for buildings – a conflicting requirement to reducing building energy needs and related CO2 emissions. Though passive design is used to reduce thermal loads of a building, a supplementary cooling system is often necessary. For such mixed-mode strategies, indirect evaporative cooling is investigated as a low energy option in the context of a warmer and drier UK climate. Analysis of the climate projections shows an increase in wet-bulb depression; providing a good indication of the cooling potential of an evaporative cooler. Modelling a mixed-mode building at two different locations, showed such a building was capable of maintaining adequate thermal comfort in future probable climates. Comparing the control climate to the scenario climate, an increase in the median of evaporative cooling load is evident. The shift is greater for London than for Glasgow with a respective 71.6% and 3.3% increase in the median annual cooling load. The study shows evaporative cooling should continue to function as an effective low-energy cooling technique in future, warming climates.

The neurochemically diverse intermedius nucleus of the medulla as a source of excitatory and inhibitory synaptic input to the nucleus tractus solitarii

Sensory afferent signals from neck muscles have been postulated to influence central cardiorespiratory control as components of postural reflexes, but neuronal pathways for this action have not been identified. The intermedius nucleus of the medulla (InM) is a target of neck muscle spindle afferents and is ideally located to influence such reflexes but is poorly investigated. To aid identification of the nucleus, we initially produced three-dimensional reconstructions of the InM in both mouse and rat. Neurochemical analysis including transgenic reporter mice expressing green fluorescent protein in GABA-synthesizing neurons, immunohistochemistry, and in situ hybridization revealed that the InM is neurochemically diverse, containing GABAegric and glutamatergic neurons with some degree of colocalization with parvalbumin, neuronal nitric oxide synthase, and calretinin. Projections from the InM to the nucleus tractus solitarius (NTS) were studied electrophysiologically in rat brainstem slices. Electrical stimulation of the NTS resulted in antidromically activated action potentials within InM neurons. In addition, electrical stimulation of the InM resulted in EPSPs that were mediated by excitatory amino acids and IPSPs mediated solely by GABA(A) receptors or by GABA(A) and glycine receptors. Chemical stimulation of the InM resulted in (1) a depolarization of NTS neurons that were blocked by NBQX (2,3-dioxo-6-nitro-1,2,3,4-tetrahydrobenzo[f]quinoxaline-7-sulfonoamide) or kynurenic acid and (2) a hyperpolarization of NTS neurons that were blocked by bicuculline. Thus, the InM contains neurochemically diverse neurons and sends both excitatory and inhibitory projections to the NTS. These data provide a novel pathway that may underlie possible reflex changes in autonomic variables after neck muscle spindle afferent activation.

Predicting the risks from climate change to forage and crop production for animal feed

Climate change is expected to bring warmer temperatures, changes to rainfall patterns, and increased frequency of extreme weather. Projections of climate impacts on feed crops show that there will likely be opportunities for increased productivity as well as considerable threats to crop productivity in different parts of the world over the next 20 to 50 years. On balance, we anticipate substantial risks to the volume, volatility, and quality of animal feed supply chains from climate change. Adaptation strategies and investment informed by high quality research at the interface of crop and animal science will be needed, both to respond to climate change and to meet the increasing demand for animal products expected over the coming decades.

Increasing influence of heat stress on French maize yields from the 1960s to the 2030s

Improved crop yield forecasts could enable more effective adaptation to climate variability and change. Here, we explore how to combine historical observations of crop yields and weather with climate model simulations to produce crop yield projections for decision relevant timescales. Firstly, the effects on historical crop yields of improved technology, precipitation and daily maximum temperatures are modelled empirically, accounting for a nonlinear technology trend and interactions between temperature and precipitation, and applied specifically for a case study of maize in France. The relative importance of precipitation variability for maize yields in France has decreased significantly since the 1960s, likely due to increased irrigation. In addition, heat stress is found to be as important for yield as precipitation since around 2000. A significant reduction in maize yield is found for each day with a maximum temperature above 32 °C, in broad agreement with previous estimates. The recent increase in such hot days has likely contributed to the observed yield stagnation. Furthermore, a general method for producing near-term crop yield projections, based on climate model simulations, is developed and utilized. We use projections of future daily maximum temperatures to assess the likely change in yields due to variations in climate. Importantly, we calibrate the climate model projections using observed data to ensure both reliable temperature mean and daily variability characteristics, and demonstrate that these methods work using retrospective predictions. We conclude that, to offset the projected increased daily maximum temperatures over France, improved technology will need to increase base level yields by 12% to be confident about maintaining current levels of yield for the period 2016–2035; the current rate of yield technology increase is not sufficient to meet this target.

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.

A method for incorporating climate variability in climate change impact assessments: sensitivity of river flows in the Eden catchment to precipitation scenarios

Interest in the impacts of climate change is ever increasing. This is particularly true of the water sector where understanding potential changes in the occurrence of both floods and droughts is important for strategic planning. Climate variability has been shown to have a significant impact on UK climate and accounting for this in future climate cahgne projections is essential to fully anticipate potential future impacts. In this paper a new resampling methodology is developed which includes the variability of both baseline and future precipitation. The resampling methodology is applied to 13 CMIP3 climate models for the 2080s, resulting in an ensemble of monthly precipitation change factors. The change factors are applied to the Eden catchment in eastern Scotland with analysis undertaken for the sensitivity of future river flows to the changes in precipitation. Climate variability is shown to influence the magnitude and direction of change of both precipitation and in turn river flow, which are not apparent without the use of the resampling methodology. The transformation of precipitation changes to river flow changes display a degree of non-linearity due to the catchment's role in buffering the response. The resampling methodology developed in this paper provides a new technique for creating climate change scenarios which incorporate the important issue of climate variability.

Transport impacts on atmosphere and climate: shipping

Emissions of exhaust gases and particles from oceangoing ships are a significant and growing contributor to the total emissions from the transportation sector. We present an assessment of the contribution of gaseous and particulate emissions from oceangoing shipping to anthropogenic emissions and air quality. We also assess the degradation in human health and climate change created by these emissions. Regulating ship emissions requires comprehensive knowledge of current fuel consumption and emissions, understanding of their impact on atmospheric composition and climate, and projections of potential future evolutions and mitigation options. Nearly 70% of ship emissions occur within 400 km of coastlines, causing air quality problems through the formation of ground-level ozone, sulphur emissions and particulate matter in coastal areas and harbours with heavy traffic. Furthermore, ozone and aerosol precursor emissions as well as their derivative species from ships may be transported in the atmosphere over several hundreds of kilometres, and thus contribute to air quality problems further inland, even though they are emitted at sea. In addition, ship emissions impact climate. Recent studies indicate that the cooling due to altered clouds far outweighs the warming effects from greenhouse gases such as carbon dioxide (CO2) or ozone from shipping, overall causing a negative present-day radiative forcing (RF). Current efforts to reduce sulphur and other pollutants from shipping may modify this. However, given the short residence time of sulphate compared to CO2, the climate response from sulphate is of the order decades while that of CO2 is centuries. The climatic trade-off between positive and negative radiative forcing is still a topic of scientific research, but from what is currently known, a simple cancellation of global mean forcing components is potentially inappropriate and a more comprehensive assessment metric is required. The CO2 equivalent emissions using the global temperature change potential (GTP) metric indicate that after 50 years the net global mean effect of current emissions is close to zero through cancellation of warming by CO2 and cooling by sulphate and nitrogen oxides.

Describing seasonal variability in the distribution of daily effective temperatures for 1985-2009 compared to 1904-1984 for De Bilt, Holland

This paper proposes a method for describing the distribution of observed temperatures on any day of the year such that the distribution and summary statistics of interest derived from the distribution vary smoothly through the year. The method removes the noise inherent in calculating summary statistics directly from the data thus easing comparisons of distributions and summary statistics between different periods. The method is demonstrated using daily effective temperatures (DET) derived from observations of temperature and wind speed at De Bilt, Holland. Distributions and summary statistics are obtained from 1985 to 2009 and compared to the period 1904–1984. A two-stage process first obtains parameters of a theoretical probability distribution, in this case the generalized extreme value (GEV) distribution, which describes the distribution of DET on any day of the year. Second, linear models describe seasonal variation in the parameters. Model predictions provide parameters of the GEV distribution, and therefore summary statistics, that vary smoothly through the year. There is evidence of an increasing mean temperature, a decrease in the variability in temperatures mainly in the winter and more positive skew, more warm days, in the summer. In the winter, the 2% point, the value below which 2% of observations are expected to fall, has risen by 1.2 °C, in the summer the 98% point has risen by 0.8 °C. Medians have risen by 1.1 and 0.9 °C in winter and summer, respectively. The method can be used to describe distributions of future climate projections and other climate variables. Further extensions to the methodology are suggested.