Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim  Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location  Europe. Methods  We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000 yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results  Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions  The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Palaeobiology of the medieval population of Albano (Rome, Italy): a combined morphological and biomolecular approach

Recent years have seen increased interest in skeletal populations from the Imperial Roman Age in Italy, but much less is known about diet and standards of living in the subsequent medieval period. To fill this gap, we conducted a morphological analysis of human remains from Albano, an Italian town near Rome, as well as a stable isotope analysis of bone collagen to reconstruct diet. The sample was recovered from a Medieval cemetery (1040–1220 cal. yr. BP) located in the gardens of the historical Palazzo Doria Pamphili in Albano. A minimum number of 40 individuals (31 adults and 9 sub-adults) were examined using standard methods. Though the general health status of the population was good, signs of cribra orbitalia and diffuse enthesopathies were noted during the morphological examination. Stable carbon and nitrogen isotope analyses of the bone collagen from 24 adult humans and three faunal bones indicate that the diet of the population may be described as predominantly terrestrial and C3-plant based although the data for some of the individuals suggest a modest consumption of C4-(millet) based or aquatic proteins. No evidence of significant dietary differences between the sexes was found. The comparison of the isotope data from Albano with those from populations recovered in the same region is consistent with a shift from a terrestrial, possibly plant foods-dominated subsistence in the Early Middle Ages to a diet with a higher contribution from animal proteins, both terrestrial and aquatic, in the Later Middle Ages.

Performance of a BP c-Si building Integrated PV system after 15years of installation

In this study, the performance, yield and characteristics of a 15 year old photovoltaic system installation has been investigated. The technology, BP Saturn modules which were steel-blue polycrystalline silicon cells are no longer in production. A bespoke monitoring system was designed and purpose built to monitor the characteristics of 6 strings, of 18 modules connected in series. The total output of the system is configured to 6.5kWp (series to parallel configuration). The PV system is demonstrating system outputs to be inferior by 0.7% per year. However,efficiency values in comparison to standard test conditions have remained practically the same. This output though very relevant to the possible performance and stability of aging cells, requires additional parametric studies to develop a more robust argument. The result presented in this paper is part of an on going investigation into PV system aging effects.

Evaluating the Bank of England density forecasts of inflation

We consider evaluating the UK Monetary Policy Committee's inflation density forecasts using probability integral transform goodness-of-fit tests. These tests evaluate the whole forecast density. We also consider whether the probabilities assigned to inflation being in certain ranges are well calibrated, where the ranges are chosen to be those of particular relevance to the MPC, given its remit of maintaining inflation rates in a band around per annum. Finally, we discuss the decision-based approach to forecast evaluation in relation to the MPC forecasts

Evaluating the Survey of Professional Forecasters probability distributions of expected inflation based on derived event probability forecasts

Techniques are proposed for evaluating forecast probabilities of events. The tools are especially useful when, as in the case of the Survey of Professional Forecasters (SPF) expected probability distributions of inflation, recourse cannot be made to the method of construction in the evaluation of the forecasts. The tests of efficiency and conditional efficiency are applied to the forecast probabilities of events of interest derived from the SPF distributions, and supplement a whole-density evaluation of the SPF distributions based on the probability integral transform approach.

Combining probability forecasts

We consider different methods for combining probability forecasts. In empirical exercises, the data generating process of the forecasts and the event being forecast is not known, and therefore the optimal form of combination will also be unknown. We consider the properties of various combination schemes for a number of plausible data generating processes, and indicate which types of combinations are likely to be useful. We also show that whether forecast encompassing is found to hold between two rival sets of forecasts or not may depend on the type of combination adopted. The relative performances of the different combination methods are illustrated, with an application to predicting recession probabilities using leading indicators.

Probability distributions or point predictions? Survey forecasts of US output growth and inflation

We consider whether survey respondents’ probability distributions, reported as histograms, provide reliable and coherent point predictions, when viewed through the lens of a Bayesian learning model. We argue that a role remains for eliciting directly-reported point predictions in surveys of professional forecasters.

Aerosol impacts on drizzle properties in warm clouds from ARM Mobile Facility maritime and continental deployments

We have extensively evaluated the response of cloud-base drizzle rate (Rcb; mm day–1) in warm clouds to liquid water path (LWP; g m–2) and to cloud condensation nuclei (CCN) number concentration (NCCN; cm–3), an aerosol proxy. This evaluation is based on a 19-month long dataset of Doppler radar, lidar, microwave radiometers and aerosol observing systems from the Atmospheric Radiation Measurement (ARM) Mobile Facility deployments at the Azores and in Germany. Assuming 0.55% supersaturation to calculate NCCN, we found a power law , indicating that Rcb decreases by a factor of 2–3 as NCCN increases from 200 to 1000 cm–3 for fixed LWP. Additionally, the precipitation susceptibility to NCCN ranges between 0.5 and 0.9, in agreement with values from simulations and aircraft measurements. Surprisingly, the susceptibility of the probability of precipitation from our analysis is much higher than that from CloudSat estimates, but agrees well with simulations from a multi-scale high-resolution aerosol-climate model. Although scale issues are not completely resolved in the intercomparisons, our results are encouraging, suggesting that it is possible for multi-scale models to accurately simulate the response of LWP to aerosol perturbations.

Dynamic Devensian ice flow in NE England: a sedimentological reconstruction

This study reconstructs the depositional environments that accompanied both ice advance and ice retreat of the last British–Irish Ice Sheet in NE England during the Last Glacial Maximum, and proposes three regional ice-flow phases. The Late Devensian (29–22 cal. ka BP) Tyne Gap Ice Stream initially deposited the Blackhall Till Formation during shelf-edge glaciation (Phase I). This subglacial traction till comprises several related facies, including stratified and laminated diamictons, tectonites, and sand and gravel beds deposited both in subglacial canals and in proglacial streams. Eventually, stagnation of the Tyne Gap Ice Stream led to ice-marginal sedimentation in County Durham (Phase II). During the Dimlington Stadial (21 cal. ka BP), the North Sea Lobe advanced towards the coastline of N Norfolk. This resulted initially in sandur deposition (widespread, tabular sand and gravel; the Peterlee Sand and Gravel Formation; Phase II) and ultimately in deposition of the Horden Till Formation (Phase III), a massive subglacial till. As the North Sea Lobe overrode previous formations, it thrusted and stacked sediments in County Durham, and dammed proglacial lakes between the east-coast ice, the Pennine uplands and the remaining Pennine ice. The North Sea Lobe retreated after Heinrich Event 1 (16 ka). This study highlights the complexity of ice flow during the Late Devensian glaciation of NE England, with changing environmental and oceanic conditions forcing a mobile and sensitive ice sheet.

The water balance of northern Africa during the mid-Holocene: an evaluation of the 6ka BP PMIP simulations

Runoff fields over northern Africa (10–25°N, 20°W–30°E) derived from 17 atmospheric general circulation models driven by identical 6 ka BP orbital forcing, sea surface temperatures, and CO2 concentration have been analyzed using a hydrological routing scheme (HYDRA) to simulate changes in lake area. The AGCM-simulated runoff produced six-fold differences in simulated lake area between models, although even the largest simulated changes considerably underestimate the observed changes in lake area during the mid-Holocene. The inter-model differences in simulated lake area are largely due to differences in simulated runoff (the squared correlation coefficient, R2, is 0.84). Most of these differences can be attributed to differences in the simulated precipitation (R2=0.83). The higher correlation between runoff and simulated lake area (R2=0.92) implies that simulated differences in evaporation have a contributory effect. When runoff is calculated using an offline land-surface scheme (BIOME3), the correlation between runoff and simulated lake area is (R2=0.94). Finally, the spatial distribution of simulated precipitation can exert an important control on the overall response.

Pollen-based reconstructions of Japanese biomes at 0, 6000 and 18,000 14C yr bp

A biomization method, which objectively assigns individual pollen assemblages to biomes ( Prentice et al., 1996 ), was tested using modern pollen data from Japan and applied to fossil pollen data to reconstruct palaeovegetation patterns 6000 and 18,000 14C yr bp Biomization started with the assignment of 135 pollen taxa to plant functional types (PFTs), and nine possible biomes were defined by specific combinations of PFTs. Biomes were correctly assigned to 54% of the 94 modern sites. Incorrect assignments occur near the altitudinal limits of individual biomes, where pollen transport from lower altitudes blurs the local pollen signals or continuous changes in species composition characterizes the range limits of biomes. As a result, the reconstructed changes in the altitudinal limits of biomes at 6000 and 18,000 14C yr bp are likely to be conservative estimates of the actual changes. The biome distribution at 6000 14C yr bp was rather similar to today, suggesting that changes in the bioclimate of Japan have been small since the mid-Holocene. At 18,000 14C yr bp the Japanese lowlands were covered by taiga and cool mixed forests. The southward expansion of these forests and the absence of broadleaved evergreen/warm mixed forests reflect a pronounced year-round cooling.

How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

Intercomparison of simulated global vegetation distribution in response to 6 kyr BP orbital forcing

The response of ten atmospheric general circulation models to orbital forcing at 6 kyr BP has been investigated using the BIOME model, which predicts equilibrium vegetation distribution, as a diagnostic. Several common features emerge: (a) reduced tropical rain forest as a consequence of increased aridity in the equatorial zone, (b) expansion of moisture-demanding vegetation in the Old World subtropics as a consequence of the expansion of the Afro–Asian monsoon, (c) an increase in warm grass/shrub in the Northern Hemisphere continental interiors in response to warming and enhanced aridity, and (d) a northward shift in the tundra–forest boundary in response to a warmer growing season at high northern latitudes. These broadscale features are consistent from model to model, but there are differences in their expression at a regional scale. Vegetation changes associated with monsoon enhancement and high-latitude summer warming are consistent with palaeoenvironmental observations, but the simulated shifts in vegetation belts are too small in both cases. Vegetation changes due to warmer and more arid conditions in the midcontinents of the Northern Hemisphere are consistent with palaeoenvironmental data from North America, but data from Eurasia suggests conditions were wetter at 6 kyr BP than today. The models show quantitatively similar vegetation changes in the intertropical zone, and in the northern and southern extratropics. The small differences among models in the magnitude of the global vegetation response are not related to differences in global or zonal climate averages, but reflect differences in simulated regional features. Regional-scale analyses will therefore be necessary to identify the underlying causes of such differences among models.

The climate and biomes of Europe at 6000 yr BP: comparison of model simulations and pollen-based reconstructions

14C-dated pollen and lake-level data from Europe are used to assess the spatial patterns of climate change between 6000 yr BP and present, as simulated by the NCAR CCM1 (National Center for Atmospheric Research, Community Climate Model, version 1) in response to the change in the Earth’s orbital parameters during this perod. First, reconstructed 6000 yr BP values of bioclimate variables obtained from pollen and lake-level data with the constrained-analogue technique are compared with simulated values. Then a 6000 yr BP biome map obtained from pollen data with an objective biome reconstruction (biomization) technique is compared with BIOME model results derived from the same simulation. Data and simulations agree in some features: warmer-than-present growing seasons in N and C Europe allowed forests to extend further north and to higher elevations than today, and warmer winters in C and E Europe prevented boreal conifers from spreading west. More generally, however, the agreement is poor. Predominantly deciduous forest types in Fennoscandia imply warmer winters than the model allows. The model fails to simulate winters cold enough, or summers wet enough, to allow temperate deciduous forests their former extended distribution in S Europe, and it incorrectly simulates a much expanded area of steppe vegetation in SE Europe. Similar errors have also been noted in numerous 6000 yr BP simulations with prescribed modern sea surface temperatures. These errors are evidently not resolved by the inclusion of interactive sea-surface conditions in the CCM1. Accurate representation of mid-Holocene climates in Europe may require the inclusion of dynamical ocean–atmosphere and/or vegetation–atmosphere interactions that most palaeoclimate model simulations have so far disregarded.

The neolithisation of Liguria (NW Italy): an environmental archaeological and palaeoenvironmental perspective

The archaeological evidence compiled for Liguria has enabled the formulation of a comprehensive model of Neolithic social, technological and economic development (∼7800–5700 cal yrs BP). The model indicates that during the Early and Middle Neolithic (∼7800–6300 cal yrs BP; ‘Impressed Ware’ and ‘Square Mouthed’ pottery cultures) human activity mainly focussed on low (coastal) and mid-altitude areas. By the Late Neolithic (∼6300–5700 cal yrs BP; ‘Chassey’ culture) farming practices were taking place over a wider range of altitudes and involved transhumant pastoralism. Complementary environmental archaeological and palaeoecological records from caves, open-air sites, lakes and mires indicate that human activities had a more significant impact on the environment than previously thought. This included clearance, especially Abies, Ulmus, Fraxinus and Tilia, and woodland utilisation and management (e.g. leaf foddering), as well as cereal cultivation and animal husbandry. The influence of Middle Holocene climatic changes, especially from ∼7800 cal yrs BP, on the direction of vegetation changes and socio-economic developments during the Neolithic remain uncertain.