Risk assessment of UK skylark populations using life-history and individual-based landscape models

Abstract: Following a workshop exercise, two models, an individual-based landscape model (IBLM) and a non-spatial life-history model were used to assess the impact of a fictitious insecticide on populations of skylarks in the UK. The chosen population endpoints were abundance, population growth rate, and the chances of population persistence. Both models used the same life-history descriptors and toxicity profiles as the basis for their parameter inputs. The models differed in that exposure was a pre-determined parameter in the life-history model, but an emergent property of the IBLM, and the IBLM required a landscape structure as an input. The model outputs were qualitatively similar between the two models. Under conditions dominated by winter wheat, both models predicted a population decline that was worsened by the use of the insecticide. Under broader habitat conditions, population declines were only predicted for the scenarios where the insecticide was added. Inputs to the models are very different, with the IBLM requiring a large volume of data in order to achieve the flexibility of being able to integrate a range of environmental and behavioural factors. The life-history model has very few explicit data inputs, but some of these relied on extensive prior modelling needing additional data as described in Roelofs et al.(2005, this volume). Both models have strengths and weaknesses; hence the ideal approach is that of combining the use of both simple and comprehensive modeling tools.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.

Physically consistent responses of the global atmospheric hydrological cycle in models and observations

Robust and physically understandable responses of the global atmospheric water cycle to a warming climate are presented. By considering interannual responses to changes in surface temperature (T), observations and AMIP5 simulations agree on an increase in column integrated water vapor at the rate 7 %/K (in line with the Clausius­Clapeyron equation) and of precipitation at the rate 2-­3 %/K (in line with energetic constraints). Using simple and complex climate models, we demonstrate that radiative forcing by greenhouse gases is currently suppressing global precipitation (P) at ~ -0.15 %/decade. Along with natural variability, this can explain why observed trends in global P over the period 1988-2008 are close to zero. Regional responses in the global water cycle are strongly constrained by changes in moisture fluxes. Model simulations show an increased moisture flux into the tropical wet region at 900 hPa and an enhanced outflow (of smaller magnitude) at around 600 hPa with warming. Moisture transport explains an increase in P in the wet tropical regions and small or negative changes in the dry regions of the subtropics in CMIP5 simulations of a warming climate. For AMIP5 simulations and satellite observations, the heaviest 5-day rainfall totals increase in intensity at ~15 %/K over the ocean with reductions at all percentiles over land. The climate change response in CMIP5 simulations shows consistent increases in P over ocean and land for the highest intensities, close to the Clausius-Clapeyron scaling of 7 %/K, while P declines for the lowest percentiles, indicating that interannual variability over land may not be a good proxy for climate change. The local changes in precipitation and its extremes are highly dependent upon small shifts in the large-scale atmospheric circulation and regional feedbacks.

The use of land-sea warming contrast under climate change to improve impact metrics

A favoured method of assimilating information from state-of-the-art climate models into integrated assessment models of climate impacts is to use the transient climate response (TCR) of the climate models as an input, sometimes accompanied by a pattern matching approach to provide spatial information. More recent approaches to the problem use TCR with another independent piece of climate model output: the land-sea surface warming ratio (φ). In this paper we show why the use of φ in addition to TCR has such utility. Multiple linear regressions of surface temperature change onto TCR and φ in 22 climate models from the CMIP3 multi-model database show that the inclusion of φ explains a much greater fraction of the inter-model variance than using TCR alone. The improvement is particularly pronounced in North America and Eurasia in the boreal summer season, and in the Amazon all year round. The use of φ as the second metric is beneficial for three reasons: firstly it is uncorrelated with TCR in state-of-the-art climate models and can therefore be considered as an independent metric; secondly, because of its projected time-invariance, the magnitude of φ is better constrained than TCR in the immediate future; thirdly, the use of two variables is much simpler than approaches such as pattern scaling from climate models. Finally we show how using the latest estimates of φ from climate models with a mean value of 1.6—as opposed to previously reported values of 1.4—can significantly increase the mean time-integrated discounted damage projections in a state-of-the-art integrated assessment model by about 15 %. When compared to damages calculated without the inclusion of the land-sea warming ratio, this figure rises to 65 %, equivalent to almost 200 trillion dollars over 200 years.

Realism of rainfall in a very high-resolution regional climate model

The realistic representation of rainfall on the local scale in climate models remains a key challenge. Realism encompasses the full spatial and temporal structure of rainfall, and is a key indicator of model skill in representing the underlying processes. In particular, if rainfall is more realistic in a climate model, there is greater confidence in its projections of future change. In this study, the realism of rainfall in a very high-resolution (1.5 km) regional climate model (RCM) is compared to a coarser-resolution 12-km RCM. This is the first time a convection-permitting model has been run for an extended period (1989–2008) over a region of the United Kingdom, allowing the characteristics of rainfall to be evaluated in a climatological sense. In particular, the duration and spatial extent of hourly rainfall across the southern United Kingdom is examined, with a key focus on heavy rainfall. Rainfall in the 1.5-km RCM is found to be much more realistic than in the 12-km RCM. In the 12-km RCM, heavy rain events are not heavy enough, and tend to be too persistent and widespread. While the 1.5-km model does have a tendency for heavy rain to be too intense, it still gives a much better representation of its duration and spatial extent. Long-standing problems in climate models, such as the tendency for too much persistent light rain and errors in the diurnal cycle, are also considerably reduced in the 1.5-km RCM. Biases in the 12-km RCM appear to be linked to deficiencies in the representation of convection.

Climate change projections and stratosphere–troposphere interaction

Climate change is expected to increase winter rainfall and flooding in many extratropical regions as evaporation and precipitation rates increase, storms become more intense and storm tracks move polewards. Here, we show how changes in stratospheric circulation could play a significant role in future climate change in the extratropics through an additional shift in the tropospheric circulation. This shift in the circulation alters climate change in regional winter rainfall by an amount large enough to significantly alter regional climate change projections. The changes are consistent with changes in stratospheric winds inducing a change in the baroclinic eddy growth rate across the depth of the troposphere. A change in mean wind structure and an equatorward shift of the tropospheric storm tracks relative to models with poor stratospheric resolution allows coupling with surface climate. Using the Atlantic storm track as an example, we show how this can double the predicted increase in extreme winter rainfall over Western and Central Europe compared to other current climate projections

A robust mechanism for strengthening of the Brewer–Dobson circulation in response to climate change: critical-layer control of subtropical wave breaking

Climate models consistently predict a strengthened Brewer–Dobson circulation in response to greenhouse gas (GHG)-induced climate change. Although the predicted circulation changes are clearly the result of changes in stratospheric wave drag, the mechanism behind the wave-drag changes remains unclear. Here, simulations from a chemistry–climate model are analyzed to show that the changes in resolved wave drag are largely explainable in terms of a simple and robust dynamical mechanism, namely changes in the location of critical layers within the subtropical lower stratosphere, which are known from observations to control the spatial distribution of Rossby wave breaking. In particular, the strengthening of the upper flanks of the subtropical jets that is robustly expected from GHG-induced tropospheric warming pushes the critical layers (and the associated regions of wave drag) upward, allowing more wave activity to penetrate into the subtropical lower stratosphere. Because the subtropics represent the critical region for wave driving of the Brewer–Dobson circulation, the circulation is thereby strengthened. Transient planetary-scale waves and synoptic-scale waves generated by baroclinic instability are both found to play a crucial role in this process. Changes in stationary planetary wave drag are not so important because they largely occur away from subtropical latitudes.

Projections of UV radiation changes in the 21st century: impact of ozone recovery and cloud effects

Monthly averaged surface erythemal solar irradiance (UV-Ery) for local noon from 1960 to 2100 has been derived using radiative transfer calculations and projections of ozone, temperature and cloud change from 14 chemistry climate models (CCM), as part of the CCMVal-2 activity of SPARC. Our calculations show the influence of ozone depletion and recovery on erythemal irradiance. In addition, we investigate UV-Ery changes caused by climate change due to increasing greenhouse gas concentrations. The latter include effects of both stratospheric ozone and cloud changes. The derived estimates provide a global picture of the likely changes in erythemal irradiance during the 21st century. Uncertainties arise from the assumed scenarios, different parameterizations – particularly of cloud effects on UV-Ery – and the spread in the CCM projections. The calculations suggest that relative to 1980, annually mean UV-Ery in the 2090s will be on average 12% lower at high latitudes in both hemispheres, 3% lower at mid latitudes, and marginally higher (1 %) in the tropics. The largest reduction (16 %) is projected for Antarctica in October. Cloud effects are responsible for 2–3% of the reduction in UV-Ery at high latitudes, but they slightly moderate it at mid-latitudes (1 %). The year of return of erythemal irradiance to values of certain milestones (1965 and 1980) depends largely on the return of column ozone to the corresponding levels and is associated with large uncertainties mainly due to the spread of the model projections. The inclusion of cloud effects in the calculations has only a small effect of the return years. At mid and high latitudes, changes in clouds and stratospheric ozone transport by global circulation changes due to greenhouse gases will sustain the erythemal irradiance at levels below those in 1965, despite the removal of ozone depleting substances.

Chemistry–climate model simulations of twenty-first century stratospheric climate and circulation changes

The response of stratospheric climate and circulation to increasing amounts of greenhouse gases (GHGs) and ozone recovery in the twenty-first century is analyzed in simulations of 11 chemistry–climate models using near-identical forcings and experimental setup. In addition to an overall global cooling of the stratosphere in the simulations (0.59 6 0.07 K decade21 at 10 hPa), ozone recovery causes a warming of the Southern Hemisphere polar lower stratosphere in summer with enhanced cooling above. The rate of warming correlates with the rate of ozone recovery projected by the models and, on average, changes from 0.8 to 0.48 Kdecade21 at 100 hPa as the rate of recovery declines from the first to the second half of the century. In the winter northern polar lower stratosphere the increased radiative cooling from the growing abundance of GHGs is, in most models, balanced by adiabatic warming from stronger polar downwelling. In the Antarctic lower stratosphere the models simulate an increase in low temperature extremes required for polar stratospheric cloud (PSC) formation, but the positive trend is decreasing over the twenty-first century in all models. In the Arctic, none of the models simulates a statistically significant increase in Arctic PSCs throughout the twenty-first century. The subtropical jets accelerate in response to climate change and the ozone recovery produces awestward acceleration of the lower-stratosphericwind over theAntarctic during summer, though this response is sensitive to the rate of recovery projected by the models. There is a strengthening of the Brewer–Dobson circulation throughout the depth of the stratosphere, which reduces the mean age of air nearly everywhere at a rate of about 0.05 yr decade21 in those models with this diagnostic. On average, the annual mean tropical upwelling in the lower stratosphere (;70 hPa) increases by almost 2% decade21, with 59% of this trend forced by the parameterized orographic gravity wave drag in the models. This is a consequence of the eastward acceleration of the subtropical jets, which increases the upward flux of (parameterized) momentum reaching the lower stratosphere in these latitudes.

Chemistry‐climate model simulations of spring Antarctic ozone

Coupled chemistry‐climate model simulations covering the recent past and continuing throughout the 21st century have been completed with a range of different models. Common forcings are used for the halogen amounts and greenhouse gas concentrations, as expected under the Montreal Protocol (with amendments) and Intergovernmental Panel on Climate Change A1b Scenario. The simulations of the Antarctic ozone hole are compared using commonly used diagnostics: the minimum ozone, the maximum area of ozone below 220 DU, and the ozone mass deficit below 220 DU. Despite the fact that the processes responsible for ozone depletion are reasonably well understood, a wide range of results is obtained. Comparisons with observations indicate that one of the reasons for the model underprediction in ozone hole area is the tendency for models to underpredict, by up to 35%, the area of low temperatures responsible for polar stratospheric cloud formation. Models also typically have species gradients that are too weak at the edge of the polar vortex, suggesting that there is too much mixing of air across the vortex edge. Other models show a high bias in total column ozone which restricts the size of the ozone hole (defined by a 220 DU threshold). The results of those models which agree best with observations are examined in more detail. For several models the ozone hole does not disappear this century but a small ozone hole of up to three million square kilometers continues to occur in most springs even after 2070.

What is the climate system able to do ‘on its own’?

The climate of the Earth, like planetary climates in general, is broadly controlled by solar irradiation, planetary albedo and emissivity as well as its rotation rate and distribution of land (with its orography) and oceans. However, the majority of climate fluctuations that affect mankind are internal modes of the general circulation of the atmosphere and the oceans. Some of these modes, such as El Nino-Southern Oscillation (ENSO), are quasi-regular and have some longer-term predictive skill; others like the Arctic and Antarctic Oscillation are chaotic and generally unpredictable beyond a few weeks. Studies using general circulation models indicate that internal processes dominate the regional climate and that some like ENSO events have even distinct global signatures. This is one of the reasons why it is so difficult to separate internal climate processes from external ones caused, for example, by changes in greenhouse gases and solar irradiation. However, the accumulation of the warmest seasons during the latest two decades is lending strong support to the forcing of the greenhouse gases. As models are getting more comprehensive, they show a gradually broader range of internal processes including those on longer time scales, challenging the interpretation of the causes of past and present climate events further.

Towards understanding the climate of Venus

ESA’s Venus Express Mission has monitored Venus since April 2006, and scientists worldwide have used mathematical models to investigate its atmosphere and model its circulation. This book summarizes recent work to explore and understand the climate of the planet through a research program under the auspices of the International Space Science Institute (ISSI) in Bern, Switzerland. Some of the unique elements that are discussed are the anomalies with Venus’ surface temperature (the huge greenhouse effect causes the surface to rise to 460°C, without which would plummet as low as -40°C), its unusual lack of solar radiation (despite being closer to the Sun, Venus receives less solar radiation than Earth due to its dense cloud cover reflecting 76% back) and the juxtaposition of its atmosphere and planetary rotation (wind speeds can climb up to 200 m/s, much faster than Venus’ sidereal day of 243 Earth-days).

Origin of regional climate differences: role of boundary conditions and model formulation in two GCMs

Model differences in projections of extratropical regional climate change due to increasing greenhouse gases are investigated using two atmospheric general circulation models (AGCMs): ECHAM4 (Max Planck Institute, version 4) and CCM3 (National Center for Atmospheric Research Community Climate Model version 3). Sea-surface temperature (SST) fields calculated from observations and coupled versions of the two models are used to force each AGCM in experiments based on time-slice methodology. Results from the forced AGCMs are then compared to coupled model results from the Coupled Model Intercomparison Project 2 (CMIP2) database. The time-slice methodology is verified by showing that the response of each model to doubled CO2 and SST forcing from the CMIP2 experiments is consistent with the results of the coupled GCMs. The differences in the responses of the models are attributed to (1) the different tropical SST warmings in the coupled simulations and (2) the different atmospheric model responses to the same tropical SST warmings. Both are found to have important contributions to differences in implied Northern Hemisphere (NH) winter extratropical regional 500 mb height and tropical precipitation climate changes. Forced teleconnection patterns from tropical SST differences are primarily responsible for sensitivity differences in the extratropical North Pacific, but have relatively little impact on the North Atlantic. There are also significant differences in the extratropical response of the models to the same tropical SST anomalies due to differences in numerical and physical parameterizations. Differences due to parameterizations dominate in the North Atlantic. Differences in the control climates of the two coupled models from the current climate, in particular for the coupled model containing CCM3, are also demonstrated to be important in leading to differences in extratropical regional sensitivity.

Midlatitude Forcing Mechanisms for Glacier Mass Balance Investigated Using General Circulation Models

A process-oriented modeling approach is applied in order to simulate glacier mass balance for individual glaciers using statistically downscaled general circulation models (GCMs). Glacier-specific seasonal sensitivity characteristics based on a mass balance model of intermediate complexity are used to simulate mass balances of Nigardsbreen (Norway) and Rhonegletscher (Switzerland). Simulations using reanalyses (ECMWF) for the period 1979–93 are in good agreement with in situ mass balance measurements for Nigardsbreen. The method is applied to multicentury integrations of coupled (ECHAM4/OPYC) and mixed-layer (ECHAM4/MLO) GCMs excluding external forcing. A high correlation between decadal variations in the North Atlantic oscillation (NAO) and mass balance of the glaciers is found. The dominant factor for this relationship is the strong impact of winter precipitation associated with the NAO. A high NAO phase means enhanced (reduced) winter precipitation for Nigardsbreen (Rhonegletscher), typically leading to a higher (lower) than normal annual mass balance. This mechanism, entirely due to internal variations in the climate system, can explain observed strong positive mass balances for Nigardsbreen and other maritime Norwegian glaciers within the period 1980–95. It can also partly be responsible for recent strong negative mass balances of Alpine glaciers.

Dynamics, stratospheric ozone, and climate change

Dynamics affects the distribution and abundance of stratospheric ozone directly through transport of ozone itself and indirectly through its effect on ozone chemistry via temperature and transport of other chemical species. Dynamical processes must be considered in order to understand past ozone changes, especially in the northern hemisphere where there appears to be significant low-frequency variability which can look “trend-like” on decadal time scales. A major challenge is to quantify the predictable, or deterministic, component of past ozone changes. Over the coming century, changes in climate will affect the expected recovery of ozone. For policy reasons it is important to be able to distinguish and separately attribute the effects of ozone-depleting substances and greenhouse gases on both ozone and climate. While the radiative-chemical effects can be relatively easily identified, this is not so evident for dynamics — yet dynamical changes (e.g., changes in the Brewer-Dobson circulation) could have a first-order effect on ozone over particular regions. Understanding the predictability and robustness of such dynamical changes represents another major challenge. Chemistry-climate models have recently emerged as useful tools for addressing these questions, as they provide a self-consistent representation of dynamical aspects of climate and their coupling to ozone chemistry. We can expect such models to play an increasingly central role in the study of ozone and climate in the future, analogous to the central role of global climate models in the study of tropospheric climate change.