Mechanistic forecasts of species responses to climate change: The promise of biophysical ecology

A core challenge in global change biology is to predict how species will respond to future environmental change and to manage these responses. To make such predictions and management actions robust to novel futures, we need to accurately characterize how organisms experience their environments and the biological mechanisms by which they respond. All organisms are thermodynamically connected to their environments through the exchange of heat and water at fine spatial and temporal scales and this exchange can be captured with biophysical models. Although mechanistic models based on biophysical ecology have a long history of development and application, their use in global change biology remains limited despite their enormous promise and increasingly accessible software. We contend that greater understanding and training in the theory and methods of biophysical ecology is vital to expand their application. Our review shows how biophysical models can be implemented to understand and predict climate change impacts on species' behavior, phenology, survival, distribution, and abundance. It also illustrates the types of outputs that can be generated, and the data inputs required for different implementations. Examples range from simple calculations of body temperature at a particular site and time, to more complex analyses of species' distribution limits based on projected energy and water balances, accounting for behavior and phenology. We outline challenges that currently limit the widespread application of biophysical models relating to data availability, training, and the lack of common software ecosystems. We also discuss progress and future developments that could allow these models to be applied to many species across large spatial extents and timeframes. Finally, we highlight how biophysical models are uniquely suited to solve global change biology problems that involve predicting and interpreting responses to environmental variability and extremes, multiple or shifting constraints, and novel abiotic or biotic environments.     

Incorporating spatial autocorrelation into species distribution models alters forecasts of climate‐mediated range shifts

Species distribution models (SDMs) are widely used to forecast changes in the spatial distributions of species and communities in response to climate change. However, spatial autocorrelation (SA) is rarely accounted for in these models, despite its ubiquity in broad‐scale ecological data. While spatial autocorrelation in model residuals is known to result in biased parameter estimates and the inflation of type I errors, the influence of unmodeled SA on species' range forecasts is poorly understood. Here we quantify how accounting for SA in SDMs influences the magnitude of range shift forecasts produced by SDMs for multiple climate change scenarios. SDMs were fitted to simulated data with a known autocorrelation structure, and to field observations of three mangrove communities from northern Australia displaying strong spatial autocorrelation. Three modeling approaches were implemented: environment‐only models (most frequently applied in species' range forecasts), and two approaches that incorporate SA; autologistic models and residuals autocovariate (RAC) models. Differences in forecasts among modeling approaches and climate scenarios were quantified. While all model predictions at the current time closely matched that of the actual current distribution of the mangrove communities, under the climate change scenarios environment‐only models forecast substantially greater range shifts than models incorporating SA. Furthermore, the magnitude of these differences intensified with increasing increments of climate change across the scenarios. When models do not account for SA, forecasts of species' range shifts indicate more extreme impacts of climate change, compared to models that explicitly account for SA. Therefore, where biological or population processes induce substantial autocorrelation in the distribution of organisms, and this is not modeled, model predictions will be inaccurate. These results have global importance for conservation efforts as inaccurate forecasts lead to ineffective prioritization of conservation activities and potentially to avoidable species extinctions.     

Predicting future coexistence in a North American ant community

Global climate change will remodel ecological communities worldwide. However, as a consequence of biotic interactions, communities may respond to climate change in idiosyncratic ways. This makes predictive models that incorporate biotic interactions necessary. We show how such models can be constructed based on empirical studies in combination with predictions or assumptions regarding the abiotic consequences of climate change. Specifically, we consider a well‐studied ant community in North America. First, we use historical data to parameterize a basic model for species coexistence. Using this model, we determine the importance of various factors, including thermal niches, food discovery rates, and food removal rates, to historical species coexistence. We then extend the model to predict how the community will restructure in response to several climate‐related changes, such as increased temperature, shifts in species phenology, and altered resource availability. Interestingly, our mechanistic model suggests that increased temperature and shifts in species phenology can have contrasting effects. Nevertheless, for almost all scenarios considered, we find that the most subordinate ant species suffers most as a result of climate change. More generally, our analysis shows that community composition can respond to climate warming in nonintuitive ways. For example, in the context of a community, it is not necessarily the most heat‐sensitive species that are most at risk. Our results demonstrate how models that account for niche partitioning and interspecific trade‐offs among species can be used to predict the likely idiosyncratic responses of local communities to climate change.     

Sensitivity analysis of the tree distribution model Phenofit to climatic input characteristics: implications for climate impact assessment

Species distributions are already affected by climate change. Forecasting their long‐term evolution requires models with thoroughly assessed validation. Our aim here is to demonstrate that the sensitivity of such models to climate input characteristics may complicate their validation and introduce uncertainties in their predictions. In this study, we conducted a sensitivity analysis of a process‐based tree distribution model Phenofit to climate input characteristics. This analysis was conducted for two North American trees which differ greatly in their distribution and eight different types of climate input for the historic period which differ in their spatial (local or gridded data) and temporal (daily vs. monthly) resolution as well as their type (locally recorded, extrapolated or simulated by General Circulation Models). We show that the climate data resolution (spatial and temporal) and their type, highly affect the model predictions. The sensitivity analysis also revealed, the importance, for global climate change impact assessment, of (i) the daily variability of temperatures in modeling the biological processes shaping species distribution, (ii) climate data at high latitudes and elevations and (iii) climate data with high spatial resolution.     

Range position and climate sensitivity: The structure of among‐population demographic responses to climatic variation

Species’ distributions will respond to climate change based on the relationship between local demographic processes and climate and how this relationship varies based on range position. A rarely tested demographic prediction is that populations at the extremes of a species’ climate envelope (e.g., populations in areas with the highest mean annual temperature) will be most sensitive to local shifts in climate (i.e., warming). We tested this prediction using a dynamic species distribution model linking demographic rates to variation in temperature and precipitation for wood frogs (Lithobates sylvaticus) in North America. Using long‐term monitoring data from 746 populations in 27 study areas, we determined how climatic variation affected population growth rates and how these relationships varied with respect to long‐term climate. Some models supported the predicted pattern, with negative effects of extreme summer temperatures in hotter areas and positive effects on recruitment for summer water availability in drier areas. We also found evidence of interacting temperature and precipitation influencing population size, such as extreme heat having less of a negative effect in wetter areas. Other results were contrary to predictions, such as positive effects of summer water availability in wetter parts of the range and positive responses to winter warming especially in milder areas. In general, we found wood frogs were more sensitive to changes in temperature or temperature interacting with precipitation than to changes in precipitation alone. Our results suggest that sensitivity to changes in climate cannot be predicted simply by knowing locations within the species’ climate envelope. Many climate processes did not affect population growth rates in the predicted direction based on range position. Processes such as species‐interactions, local adaptation, and interactions with the physical landscape likely affect the responses we observed. Our work highlights the need to measure demographic responses to changing climate.     

Riparian plant species differ in sensitivity to both the mean and variance in groundwater stores

确定气候的长期缓慢变化与放大的气候变异性之间相互作用的生态后果是植物生态学的一个重要前沿研究。我们将最近的气候敏感性函数方法与修订的水文“水桶模型”相结合，以改进对植物物种如何响应地下水资源平均值和方差变化的预测。利用河岸植被覆盖和地下水位长期数据集的时空变化，构建了旱地河岸廊道常见植物种类的首个地下水敏感性函数。我们的研究结果展示了这种方法的重要性，它可以识别在日益变化的气候层，随着地下水储量的下降，哪些植物物种能够繁茂的生长(或衰落)。研究结果表明，河岸带植物物种对地下水位平均值和方差的敏感性各不相同。杨树(Populus deltoides ssp. wislizenii)随着地下水年际方差的增大，植被覆盖度预计将下降，而狼柳树 (Salix exigua)和其他原生湿地物种预计将受益于更大的年际方差。非本地物种对地下水变化不敏感，但俄罗斯橄榄(Elaeagnus angustifolia)在更深的平均地下水位条件下敏感性将降低。暖空气温度调节了杨树对地下水的敏感性，在最高温度更高的年份/地点，杨树对地下水的变化比在低温地点/时期更敏感。地下水年内变异系数越大，杨树覆盖度下降幅度越大，但与年际变异系数相关性不显著，这可能因为相对于杨树的寿命，这个时间(16 a)太短。与此相反，非本地柳树(Tamarix chinensis)覆盖度随地下水年内和年际变异系数的增加而增加。总之，我们的研究结果预测，地下水变异性和平均值的变化将通过单个植物物种对地下水储量中的平均值和方差的不同敏感性影响河岸带植物群落。     

Evidence that implicit assumptions of ‘no evolution’ of disease vectors in changing environments can be violated on a rapid timescale

Projected impacts of climate change on vector-borne disease dynamics must consider many variables relevant to hosts, vectors and pathogens, including how altered environmental characteristics might affect the spatial distributions of vector species. However, many predictive models for vector distributions consider their habitat requirements to be fixed over relevant time-scales, when they may actually be capable of rapid evolutionary change and even adaptation. We examine the genetic signature of a spatial expansion by an invasive vector into locations with novel temperature conditions compared to its native range as a proxy for how existing vector populations may respond to temporally changing habitat. Specifically, we compare invasions into different climate ranges and characterize the importance of selection from the invaded habitat. We demonstrate that vector species can exhibit evolutionary responses (altered allelic frequencies) to a temperature gradient in as little as 7–10 years even in the presence of high gene flow, and further, that this response varies depending on the strength of selection. We interpret these findings in the context of climate change predictions for vector populations and emphasize the importance of incorporating vector evolution into models of future vector-borne disease dynamics.     

Inclusion of local environmental conditions alters high-latitude vegetation change predictions based on bioclimatic models

Current predictions of how species will respond to climate change are typically based on coarse-grained climate surfaces utilizing bioclimate envelope modelling. However, the suitability of environmental conditions for a given species might result from a variety of factors including some unrelated to climate. To address this issue, we investigated whether the inclusion of topographical and soil information in bioclimatic envelope models would significantly alter predictions of climate change—induced fine-scale tree and shrub species range size changes at the tree-limit in subarctic Europe. Using generalized additive models and data on current climate and species distributions and three different climate scenarios for the period 2040–2069, we developed predictions of the currently suitable area and potential range size changes of seven tree and shrub species in an area of 1,100 km² at a resolution of 1-ha. The inclusion of topography and soil information increased the predictive accuracy of climate-only models for all studied species. The predicted changes in species distribution volumes were contradictory, and the predicted occurrences varied greatly depending on the model used. Our results therefore support the arguments that vegetation responses to climate change can be influenced by local environmental conditions and that attention should be paid to the combined effects of these factors. We conclude that disregarding local topography and soil conditions in bioclimatic models may result in biased projections of range expansions and the associated colonization, extinction and turnover assessments.     

Improving assessment and modelling of climate change impacts on global terrestrial biodiversity

Understanding how species and ecosystems respond to climate change has become a major focus of ecology and conservation biology. Modelling approaches provide important tools for making future projections, but current models of the climate-biosphere interface remain overly simplistic, undermining the credibility of projections. We identify five ways in which substantial advances could be made in the next few years: (i) improving the accessibility and efficiency of biodiversity monitoring data, (ii) quantifying the main determinants of the sensitivity of species to climate change, (iii) incorporating community dynamics into projections of biodiversity responses, (iv) accounting for the influence of evolutionary processes on the response of species to climate change, and (v) improving the biophysical rule sets that define functional groupings of species in global models.     

Uncertainty in assessing the impacts of global change with coupled dynamic species distribution and population models

Concern over rapid global changes and the potential for interactions among multiple threats are prompting scientists to combine multiple modelling approaches to understand impacts on biodiversity. A relatively recent development is the combination of species distribution models, land‐use change predictions, and dynamic population models to predict the relative and combined impacts of climate change, land‐use change, and altered disturbance regimes on species' extinction risk. Each modelling component introduces its own source of uncertainty through different parameters and assumptions, which, when combined, can result in compounded uncertainty that can have major implications for management. Although some uncertainty analyses have been conducted separately on various model components – such as climate predictions, species distribution models, land‐use change predictions, and population models – a unified sensitivity analysis comparing various sources of uncertainty in combined modelling approaches is needed to identify the most influential and problematic assumptions. We estimated the sensitivities of long‐run population predictions to different ecological assumptions and parameter settings for a rare and endangered annual plant species (Acanthomintha ilicifolia, or San Diego thornmint). Uncertainty about habitat suitability predictions, due to the choice of species distribution model, contributed most to variation in predictions about long‐run populations.     

Climate change at the landscape scale: predicting fine-grained spatial heterogeneity in warming and potential refugia for vegetation

Current predictions of how species will respond to climate change are based on coarse‐grained climate surfaces or idealized scenarios of uniform warming. These predictions may erroneously estimate the risk of extinction because they neglect to consider spatially heterogenous warming at the landscape scale or identify refugia where species can persist despite unfavourable regional climate. To address this issue, we investigated the heterogeneity in warming that has occurred in a 10 km × 10 km area from 1972 to 2007. We developed estimates by combining long‐term daily observations from a limited number of weather stations with a more spatially comprehensive dataset (40 sites) obtained during 2005–2006. We found that the spatial distribution of warming was greater inland, at lower elevations, away from streams, and at sites exposed to the northwest (NW). These differences corresponded with changes in weather patterns, such as an increasing frequency of hot, dry NW winds. As plant species were biased in the topographic and geographic locations they occupied, these differences meant that some species experienced more warming than others, and are at greater risk from climate change. This species bias could not be detected at coarser scales. The uneven seasonal nature of warming (e.g. more warming in winter, minimums increased more than maximums) means that climate change predictions will vary according to which predictors are selected in species distribution models. Models based on a limited set of predictors will produce erroneous predictions when the correct limiting factor is not selected, and this is difficult to avoid when temperature predictors are correlated because they are produced using elevation‐sensitive interpolations. The results reinforce the importance of downscaling coarse‐grained (∼50 km) temperature surfaces, and suggest that the accuracy of this process could be improved by considering regional weather patterns (wind speed, direction, humidity) and topographic exposure to key wind directions.     

Multi‐model comparison highlights consistency in predicted effect of warming on a semi‐arid shrub

A number of modeling approaches have been developed to predict the impacts of climate change on species distributions, performance, and abundance. The stronger the agreement from models that represent different processes and are based on distinct and independent sources of information, the greater the confidence we can have in their predictions. Evaluating the level of confidence is particularly important when predictions are used to guide conservation or restoration decisions. We used a multi‐model approach to predict climate change impacts on big sagebrush (Artemisia tridentata), the dominant plant species on roughly 43 million hectares in the western United States and a key resource for many endemic wildlife species. To evaluate the climate sensitivity of A. tridentata, we developed four predictive models, two based on empirically derived spatial and temporal relationships, and two that applied mechanistic approaches to simulate sagebrush recruitment and growth. This approach enabled us to produce an aggregate index of climate change vulnerability and uncertainty based on the level of agreement between models. Despite large differences in model structure, predictions of sagebrush response to climate change were largely consistent. Performance, as measured by change in cover, growth, or recruitment, was predicted to decrease at the warmest sites, but increase throughout the cooler portions of sagebrush's range. A sensitivity analysis indicated that sagebrush performance responds more strongly to changes in temperature than precipitation. Most of the uncertainty in model predictions reflected variation among the ecological models, raising questions about the reliability of forecasts based on a single modeling approach. Our results highlight the value of a multi‐model approach in forecasting climate change impacts and uncertainties and should help land managers to maximize the value of conservation investments.     

Incorporating local adaptation into forecasts of species’ distribution and abundance under climate change

Populations of many species are genetically adapted to local historical climate conditions. Yet most forecasts of species’ distributions under climate change have ignored local adaptation (LA), which may paint a false picture of how species will respond across their geographic ranges. We review recent studies that have incorporated intraspecific variation, a potential proxy for LA, into distribution forecasts, assess their strengths and weaknesses, and make recommendations for how to improve forecasts in the face of LA. The three methods used so far (species distribution models, response functions, and mechanistic models) reflect a trade‐off between data availability and the ability to rigorously demonstrate LA to climate. We identify key considerations for incorporating LA into distribution forecasts that are currently missing from many published studies, including testing the spatial scale and pattern of LA, the confounding effects of LA to nonclimatic or biotic drivers, and the need to incorporate empirically based dispersal or gene flow processes. We suggest approaches to better evaluate these aspects of LA and their effects on species‐level forecasts. In particular, we highlight demographic and dynamic evolutionary models as promising approaches to better integrate LA into forecasts, and emphasize the importance of independent model validation. Finally, we urge closer examination of how LA will alter the responses of central vs. marginal populations to allow stronger generalizations about changes in distribution and abundance in the face of LA.     

Implications of Climate Change for Bird Conservation in the Southwestern U.S. under Three Alternative Futures

Future expected changes in climate and human activity threaten many riparian habitats, particularly in the southwestern U.S. Using Maximum Entropy (MaxEnt3.3.3) modeling, we characterized habitat relationships and generated spatial predictions of habitat suitability for the Lucy’s warbler (Oreothlypis luciae), the Southwestern willow flycatcher (Empidonax traillii extimus) and the Western yellow-billed cuckoo (Coccyzus americanus). Our goal was to provide site- and species-specific information that can be used by managers to identify areas for habitat conservation and/or restoration along the Rio Grande in New Mexico. We created models of suitable habitat for each species based on collection and survey samples and climate, biophysical, and vegetation data. We projected habitat suitability under future climates by applying these models to conditions generated from three climate models for 2030, 2060 and 2090. By comparing current and future distributions, we identified how habitats are likely to change as a result of changing climate and the consequences of those changes for these bird species. We also examined whether land ownership of high value sites shifts under changing climate conditions. Habitat suitability models performed well. Biophysical characteristics were more important that climate conditions for predicting habitat suitability with distance to water being the single most important predictor. Climate, though less important, was still influential and led to declines of suitable habitat of more than 60% by 2090. For all species, suitable habitat tended to shrink over time within the study area leaving a few core areas of high importance. Overall, climate changes will increase habitat fragmentation and reduce breeding habitat patch size. The best strategy for conserving bird species within the Rio Grande will include measures to maintain and restore critical habitat refugia. This study provides an example of a presence-only habitat model that can be used to inform the management of species at intermediate scales.     

Untangling climate and water chemistry to predict changes in freshwater macrophyte distributions

Forecasting changes in the distributions of macrophytes is essential to understanding how aquatic ecosystems will respond to climate and environmental changes. Previous work in aquatic ecosystems has used climate data at large scales and chemistry data at small scales; the consequence of using these different data types has not been evaluated. This study combines a survey of macrophyte diversity and water chemistry measurements at a large regional scale to demonstrate the feasibility and necessity of including ecological measurements, in addition to climate data, in species distribution models of aquatic macrophytes. A survey of 740 water bodies stratified across 327,000 square kilometers was conducted to document Characeae (green macroalgae) species occurrence and water chemistry data. Chemistry variables and climate data were used separately and in concert to develop species distribution models for ten species across the study area. The impacts of future environmental changes on species distributions were modeled using a range of global climate models (GCMs), representative concentration pathways (RCPs), and pollution scenarios. Models developed with chemistry variables generally gave the most accurate predictions of species distributions when compared with those using climate variables. Calcium and conductivity had the highest total relative contribution to models across all species. Habitat changes were most pronounced in scenarios with increased road salt and deicer influences, with two species predicted to increase in range by >50% and four species predicted to decrease in range by >50%. Species of Characeae have distinct habitat ranges that closely follow spatial patterns of water chemistry. Species distribution models built with climate data alone were insufficient to predict changes in distributions in the study area. The development and implementation of standardized, large‐scale water chemistry databases will aid predictions of habitat changes for aquatic ecosystems.     

Biotic interactions influence the projected distribution of a specialist mammal under climate change

Aim

To measure the effects of including biotic interactions on climate‐based species distribution models (SDMs) used to predict distribution shifts under climate change. We evaluated the performance of distribution models for an endangered marsupial, the northern bettong (Bettongia tropica), comparing models that used only climate variables with models that also took into account biotic interactions.

Location

North‐east Queensland, Australia.

Methods

We developed separate climate‐based distribution models for the northern bettong, its two main resources and a competitor species. We then constructed models for the northern bettong by including climate suitability estimates for the resources and competitor as additional predictor variables to make climate + resource and climate + resource + competition models. We projected these models onto seven future climate scenarios and compared predictions of northern bettong distribution made by these differently structured models, using a ‘global’ metric, the I similarity statistic, to measure overlap in distribution and a ‘local’ metric to identify where predictions differed significantly.

Results

Inclusion of food resource biotic interactions improved model performance. Over moderate climate changes, up to 3.0 °C of warming, the climate‐only model for the northern bettong gave similar predictions of distribution to the more complex models including interactions, with differences only at the margins of predicted distributions. For climate changes beyond 3.0 °C, model predictions diverged significantly. The interactive model predicted less contraction of distribution than the simpler climate‐only model.

Main conclusions

Distribution models that account for interactions with other species, in particular direct resources, improve model predictions in the present‐day climate. For larger climate changes, shifts in distribution of interacting species cause predictions of interactive models to diverge from climate‐only models. Incorporating interactions with other species in SDMs may be needed for long‐term prediction of changes in distribution of species under climate change, particularly for specialized species strongly dependent on a small number of biotic interactions.     

Spatial scale,topography and thermoregulatory behaviour interact when modelling species’ thermal niches

The spatial scale at which climate and species’ occupancy data are gathered, and the resolution at which ecological models are run, can strongly influence predictions of species performance and distributions. Running model simulations at coarse rather than fine spatial resolutions, for example, can determine if a model accurately predicts the distribution of a species. The impacts of spatial scale on a model's accuracy are particularly pronounced across mountainous terrain. Understanding how these discrepancies arise requires a modelling approach in which the underlying processes that determine a species’ distribution are explicitly described. Here we use a process‐based model to explore how spatial resolution, topography and behaviour alter predictions of a species thermal niche, which in turn constrains its survival and geographic distribution. The model incorporates biophysical equations to predict the operative temperature (T_e), thermal‐dependent performance and survival of a typical insect, with a complex life‐cycle, in its microclimate. We run this model with geographic data from a mountainous terrain in South Africa using climate data at three spatial resolutions. We also explore how behavioural thermoregulation affects predictions of a species performance and survival by allowing the animal to select the optimum thermal location within each coarse‐grid cell. At the regional level, coarse‐resolution models predicted lower T_e at low elevations and higher T_e at high elevations than models run at fine‐resolutions. These differences were more prominent on steep, north‐facing slopes. The discrepancies in T_e in turn affected estimates of the species thermal niche. The modelling framework revealed how spatial resolution and topography influence predictions of species distribution models, including the potential impacts of climate change. These systematic biases must be accounted for when interpreting the outputs of future modelling studies, particularly when species distributions are predicted to shift from uniform to topographically heterogeneous landscapes.     

Simple models for stomatal conductance derived from a process model: cross-validation against sap flux data

Representation of stomatal physiology in models of plant-atmosphere gas exchange is minimal, and direct application of process-based models is limited by difficulty of parameter estimation. We derived simple models of stomatal conductance from a recent process-based model, and cross-validated them against measurements of sap flux (176-365 d in length) in 36 individual trees of two age classes for two Eucalyptus species across seven sites in the mountains of southeastern Australia. The derived models - which are driven by irradiance and evaporative demand and have two to four parameters that represent sums and products of biophysical parameters in the process model - reproduced a median 83-89% of observed variance in half-hourly and diurnally averaged sap flux, and performed similarly whether fitted using a random sample of all data or using 1 month of data from spring or autumn. Our simple models are an advance in predicting plant water use because their parameters are transparently related to reduced processes and properties, enabling easy accommodation of improved knowledge about how those parameters respond to environmental change and differ among species.     

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species‐based Dynamic Bioclimate Envelope Model (DBEM) with a size‐based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness‐of‐fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter‐specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.     

Introduction: food crops in a changing climate

Changes in both the mean and the variability of climate, whether naturally forced, or due to human activities, pose a threat to crop production globally. This paper summarizes discussions of this issue at a meeting of the Royal Society in April 2005. Recent advances in understanding the sensitivity of crops to weather, climate and the levels of particular gases in the atmosphere indicate that the impact of these factors on crop yields and quality may be more severe than previously thought. There is increasing information on the importance to crop yields of extremes of temperature and rainfall at key stages of crop development. Agriculture will itself impact on the climate system and a greater understanding of these feedbacks is needed. Complex models are required to perform simulations of climate variability and change, together with predictions of how crops will respond to different climate variables. Variability of climate, such as that associated with El Ni?o events, has large impacts on crop production. If skilful predictions of the probability of such events occurring can be made a season or more in advance, then agricultural and other societal responses can be made. The development of strategies to adapt to variations in the current climate may also build resilience to changes in future climate. Africa will be the part of the world that is most vulnerable to climate variability and change, but knowledge of how to use climate information and the regional impacts of climate variability and change in Africa is rudimentary. In order to develop appropriate adaptation strategies globally, predictions about changes in the quantity and quality of food crops need to be considered in the context of the entire food chain from production to distribution, access and utilization. Recommendations for future research priorities are given.