Pollen-based quantitative reconstructions of Holocene regional vegetation cover (plant-functional types and land-cover types) in Europe suitable for climate modelling

We present quantitative reconstructions of regional vegetation cover in north-western Europe, western Europe north of the Alps, and eastern Europe for five time windows in the Holocene around 6k, 3k, 0.5k, 0.2k, and 0.05k calendar years before present (bp)] at a 1 degrees x1 degrees spatial scale with the objective of producing vegetation descriptions suitable for climate modelling. The REVEALS model was applied on 636 pollen records from lakes and bogs to reconstruct the past cover of 25 plant taxa grouped into 10 plant-functional types and three land-cover types evergreen trees, summer-green (deciduous) trees, and open land]. The model corrects for some of the biases in pollen percentages by using pollen productivity estimates and fall speeds of pollen, and by applying simple but robust models of pollen dispersal and deposition. The emerging patterns of tree migration and deforestation between 6k bp and modern time in the REVEALS estimates agree with our general understanding of the vegetation history of Europe based on pollen percentages. However, the degree of anthropogenic deforestation (i.e. cover of cultivated and grazing land) at 3k, 0.5k, and 0.2k bp is significantly higher than deduced from pollen percentages. This is also the case at 6k in some parts of Europe, in particular Britain and Ireland. Furthermore, the relationship between summer-green and evergreen trees, and between individual tree taxa, differs significantly when expressed as pollen percentages or as REVEALS estimates of tree cover. For instance, when Pinus is dominant over Picea as pollen percentages, Picea is dominant over Pinus as REVEALS estimates. These differences play a major role in the reconstruction of European landscapes and for the study of land cover-climate interactions, biodiversity and human resources.

Temperate Mountain Forest Biodiversity under Climate Change: Compensating Negative Effects by Increasing Structural Complexity

Species adapted to cold-climatic mountain environments are expected to face a high risk of range contractions, if not local extinctions under climate change. Yet, the populations of many endothermic species may not be primarily affected by physiological constraints, but indirectly by climate-induced changes of habitat characteristics. In mountain forests, where vertebrate species largely depend on vegetation composition and structure, deteriorating habitat suitability may thus be mitigated or even compensated by habitat management aiming at compositional and structural enhancement. We tested this possibility using four cold-adapted bird species with complementary habitat requirements as model organisms. Based on species data and environmental information collected in 300 1-km2 grid cells distributed across four mountain ranges in central Europe, we investigated (1) how species’ occurrence is explained by climate, landscape, and vegetation, (2) to what extent climate change and climate-induced vegetation changes will affect habitat suitability, and (3) whether these changes could be compensated by adaptive habitat management. Species presence was modelled as a function of climate, landscape and vegetation variables under current climate; moreover, vegetation-climate relationships were assessed. The models were extrapolated to the climatic conditions of 2050, assuming the moderate IPCC-scenario A1B, and changes in species’ occurrence probability were quantified. Finally, we assessed the maximum increase in occurrence probability that could be achieved by modifying one or multiple vegetation variables under altered climate conditions. Climate variables contributed significantly to explaining species occurrence, and expected climatic changes, as well as climate-induced vegetation trends, decreased the occurrence probability of all four species, particularly at the low-altitudinal margins of their distribution. These effects could be partly compensated by modifying single vegetation factors, but full compensation would only be achieved if several factors were changed in concert. The results illustrate the possibilities and limitations of adaptive species conservation management under climate change.

The role of climate-driven chemical weathering on soil production

Climate plays an important role in controlling rates of weathering and weathered regolith production. Regolith production functions, however, seldom take climate parameters into account. Based on a climate-dependent weathered regolith production model, at low denudation rates, relative regolith thicknesses are less sensitive to changes in precipitation rates, while at high denudation rates, small changes in climatic parameters can result in complete stripping of hillslopes. This pattern is compounded by the long residence times and system response times associated with low denudation rates, and vice versa. As others have shown, the transition between regolith-mantled and bedrock slopes is dependent on the ratio of denudation to production. Here, we further suggest that this is itself a function of precipitation rate and temperature. We suggest that climatic parameters can be easily incorporated into existing soil production models and that such additions improve the predictive power of soil production models. (C) 2013 Elsevier B.V. All rights reserved.

Simple measures of climate, soil properties and plant traits predict national-scale grassland soil carbon stocks

Soil carbon (C) storage is a key ecosystem service. Soil C stocks play a vital role in soil fertility and climate regulation, but the factors that control these stocks at regional and national scales are unknown, particularly when their composition and stability are considered. As a result, their mapping relies on either unreliable proxy measures or laborious direct measurements. Using data from an extensive national survey of English grasslands, we show that surface soil (0–7 cm) C stocks in size fractions of varying stability can be predicted at both regional and national scales from plant traits and simple measures of soil and climatic conditions. Soil C stocks in the largest pool, of intermediate particle size (50–250 μm), were best explained by mean annual temperature (MAT), soil pH and soil moisture content. The second largest C pool, highly stable physically and biochemically protected particles (0·45–50 μm), was explained by soil pH and the community abundance-weighted mean (CWM) leaf nitrogen (N) content, with the highest soil C stocks under N-rich vegetation. The C stock in the small active fraction (250–4000 μm) was explained by a wide range of variables: MAT, mean annual precipitation, mean growing season length, soil pH and CWM specific leaf area; stocks were higher under vegetation with thick and/or dense leaves. Testing the models describing these fractions against data from an independent English region indicated moderately strong correlation between predicted and actual values and no systematic bias, with the exception of the active fraction, for which predictions were inaccurate. Synthesis and applications. Validation indicates that readily available climate, soils and plant survey data can be effective in making local- to landscape-scale (1–100 000 km2) soil C stock predictions. Such predictions are a crucial component of effective management strategies to protect C stocks and enhance soil C sequestration.

Comparison of co-located independent ground-based middle atmospheric wind and temperature measurements with numerical weather prediction models

High-resolution, ground-based and independent observations including co-located wind radiometer, lidar stations, and infrasound instruments are used to evaluate the accuracy of general circulation models and data-constrained assimilation systems in the middle atmosphere at northern hemisphere midlatitudes. Systematic comparisons between observations, the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analyses including the recent Integrated Forecast System cycles 38r1 and 38r2, the NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalyses, and the free-running climate Max Planck Institute–Earth System Model–Low Resolution (MPI-ESM-LR) are carried out in both temporal and spectral dom ains. We find that ECMWF and MERRA are broadly consistent with lidar and wind radiometer measurements up to ~40 km. For both temperature and horizontal wind components, deviations increase with altitude as the assimilated observations become sparser. Between 40 and 60 km altitude, the standard deviation of the mean difference exceeds 5 K for the temperature and 20 m/s for the zonal wind. The largest deviations are observed in winter when the variability from large-scale planetary waves dominates. Between lidar data and MPI-ESM-LR, there is an overall agreement in spectral amplitude down to 15–20 days. At shorter time scales, the variability is lacking in the model by ~10 dB. Infrasound observations indicate a general good agreement with ECWMF wind and temperature products. As such, this study demonstrates the potential of the infrastructure of the Atmospheric Dynamics Research Infrastructure in Europe project that integrates various measurements and provides a quantitative understanding of stratosphere-troposphere dynamical coupling for numerical weather prediction applications.

Independent at heart: persistent association of altitude with ischaemic heart disease mortality after consideration of climate, topography and built environment.

BACKGROUND Living at higher altitude was dose-dependently associated with lower risk of ischaemic heart disease (IHD). Higher altitudes have different climatic, topographic and built environment properties than lowland regions. It is unclear whether these environmental factors mediate/confound the association between altitude and IHD. We examined how much of the altitude-IHD association is explained by variations in exposure at place of residence to sunshine, temperature, precipitation, aspect, slope and distance to main road. METHODS We included 4.2 million individuals aged 40-84 at baseline living in Switzerland at altitudes 195-2971 m above sea level (ie, full range of residence), providing 77 127 IHD deaths. Mortality data 2000-2008, sociodemographic/economic information and coordinates of residence were obtained from the Swiss National Cohort, a longitudinal, census-based record linkage study. Environment information was modelled to residence level using Weibull regression models. RESULTS In the model not adjusted for other environmental factors, IHD mortality linearly decreased with increasing altitude resulting in a lower risk (HR, 95% CI 0.67, 0.60 to 0.74) for those living >1500 m (vs<600 m). This association remained after adjustment for all other environmental factors 0.74 (0.66 to 0.82). CONCLUSIONS The benefit of living at higher altitude was only partially confounded by variations in climate, topography and built environment. Rather, physical environment factors appear to have an independent effect and may impact on cardiovascular health in a cumulative way. Inclusion of additional modifiable factors as well as individual information on traditional IHD risk factors in our combined environmental model could help to identify strategies for the reduction of inequalities in IHD mortality.

Variability of daily winter wind speed distribution over Northern Europe during the past millennium in regional and global climate simulations

We analyse the variability of the probability distribution of daily wind speed in wintertime over Northern and Central Europe in a series of global and regional climate simulations covering the last centuries, and in reanalysis products covering approximately the last 60 years. The focus of the study lies on identifying the link of the variations in the wind speed distribution to the regional near-surface temperature, to the meridional temperature gradient and to the North Atlantic Oscillation. Our main result is that the link between the daily wind distribution and the regional climate drivers is strongly model dependent. The global models tend to behave similarly, although they show some discrepancies. The two regional models also tend to behave similarly to each other, but surprisingly the results derived from each regional model strongly deviates from the results derived from its driving global model. In addition, considering multi-centennial timescales, we find in two global simulations a long-term tendency for the probability distribution of daily wind speed to widen through the last centuries. The cause for this widening is likely the effect of the deforestation prescribed in these simulations. We conclude that no clear systematic relationship between the mean temperature, the temperature gradient and/or the North Atlantic Oscillation, with the daily wind speed statistics can be inferred from these simulations. The understand- ing of past and future changes in the distribution of wind speeds, and thus of wind speed extremes, will require a detailed analysis of the representation of the interaction between large-scale and small-scale dynamics.

North Pacific Carbon Cycle Response to Climate Variability on Seasonal to Decadal Timescales

Climate variability drives significant changes in the physical state of the North Pacific, and there may be important impacts of this variability on the upper ocean carbon balance across the basin. We address this issue by considering the response of seven biogeochemical ocean models to climate variability in the North Pacific. The models' upper ocean pCO(2) and air-sea CO(2) flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO(2) and its temperature- and non-temperature-driven components at three contrasting oceanographic sites capture the basic features found in observations (Takahashi et al., 2002, 2006; Keeling et al., 2004; Brix et al., 2004). However, particularly in the Western Subarctic Gyre, the models have difficulty representing the temporal structure of the total pCO(2) seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO(2) flux interannual variability (1 sigma) in the North Pacific is smaller ( ranges across models from 0.03 to 0.11 PgC/yr) than in the Tropical Pacific ( ranges across models from 0.08 to 0.19 PgC/yr), and the time series of the first or second EOF of the air-sea CO(2) flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO(2) flux anomalies are correlated with the PDO, their magnitudes are small ( up to +/- 0.025 PgC/yr ( 1 sigma)). Flux anomalies are damped because anomalies in the key drivers of pCO(2) ( temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO(2) anomalies.

Age models, color measurements, and calculated sea surface temperatures of sediment cores from the Bering Sea and the Sea of Okhotsk

Past changes in North Pacific sea surface temperatures and sea-ice conditions are proposed to play a crucial role in deglacial climate development and ocean circulation but are less well known than from the North Atlantic. Here, we present new alkenone-based sea surface temperature records from the subarctic northwest Pacific and its marginal seas (Bering Sea and Sea of Okhotsk) for the time interval of the last 15 kyr, indicating millennial-scale sea surface temperature fluctuations similar to short-term deglacial climate oscillations known from Greenland ice-core records. Past changes in sea-ice distribution are derived from relative percentage of specific diatom groups and qualitative assessment of the IP25 biomarker related to sea-ice diatoms. The deglacial variability in sea-ice extent matches the sea surface temperature fluctuations. These fluctuations suggest a linkage to deglacial variations in Atlantic meridional overturning circulation and a close atmospheric coupling between the North Pacific and North Atlantic. During the Holocene the subarctic North Pacific is marked by complex sea surface temperature trends, which do not support the hypothesis of a Holocene seesaw in temperature development between the North Atlantic and the North Pacific.

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Age models and stabile isotope ratios from SO35 cores

Stable isotope measurements on the planktonic foraminifer Globigerinoides ruber (white) have been carried out on a number of selected deep-seas sediment cores from the South Lau and Norlh Fiji Basins. The d18O-curves show good correlation with the inter-ocean oraphic correlation composite d18O-record of the standard reference section (Prell et al. 1986), which, in combination with the chronostratigraphic classifications of Herterich & Sarnthein (1984, modified) and Imbrie et al. 1984), allows a detailed dating of the sedimentary sequences. The deepest layers in core no. 119 (southern Lau Basin) could be assigned to Isotope Stage 24. Measurements made on bulk carbonate in two cores show a much higher glacial-interglacial amplitude, allowing the general identification of the conventional oxygen isotope stages. The d13C-values of the benthic foraminifer Cibicidoides wuellerstorfi show progressively lighter values northwards reflecting an increasing contribution of the isotopically lighter CO2 from the remineralisation of organic matter during the general northward movement of the deep water masses. Cyclicities in the sedimentation rates were observed in core nos. 117 and 119 (both southern Lau Basin) where the interglacials exhibit higher levels than the glacials. Calculated new or export paleoproductivity show that the glacials had higher productivity in the euphotic zone. From the oxygen isotope stratigraphy, the five ash layers in core nos. 117 and 119 could be dated as about 530 ka B.P. in Stage 14, 695 ka B.P. in Stage 18, 775 ka B.P. in Stage 21, 790 ka B.P. and 825 ka B.P. in Stage 22. Carbonate dissolution occurred during stages 5, 8 and 10 to 12.

Revised age models and foraminiferal records obtained for the Last Interglacial period in six marine sediment cores

The dataset contains the revised age models and foraminiferal records obtained for the Last Interglacial period in six marine sediment cores: - the Southern Ocean core MD02-2488 (age model, sea surface temperatures, benthic d18O and d13C for the period 136-108 ka), - the North Atlantic core MD95-2042 (age model, planktic d18O, benthic d18O and d13C for the period 135-110 ka), - the North Atlantic core ODP 980 (age model, planktic d18O, sea surface temperatures, seawater d18O, benthic d18O and d13C, ice-rafted detritus for the period 135-110 ka), - the North Atlantic core CH69-K09 (age model, planktic d18O, sea surface temperatures, seawater d18O, benthic d18O and d13C, ice-rafted detritus for the period 135-110 ka), - the Norwegian Sea core MD95-2010 (age model, percentage of Neogloboquadrina pachyderma sinistral, sea surface temperatures, benthic d18O, ice-rafted detritus for the period 134-110 ka), - the Labrador Sea core EW9302-JPC2 (age model, percentage of Neogloboquadrina pachyderma sinistral, sea surface temperatures, benthic d18O for the period 134-110 ka).

Revised composite depth scale, Fe content and orbital tuning of Middle Eocene Climate Optimum samples from ODP Site 207-1260

A high-resolution stratigraphy is essential toward deciphering climate variability in detail and understanding causality arguments of events in earth history. Because the highly dynamic middle to late Eocene provides a suitable testing ground for carbon cycle models for a waning warm world, an accurate time scale is needed to decode climate-driving mechanisms. Here we present new results from ODP Site 1260 (Leg 207) which covers a unique expanded middle Eocene section (magnetochrons C18r to C20r, late Lutetian to early Bartonian) of the tropical western Atlantic including the chron C19r transient hyperthermal event and the Middle Eocene Climate Optimum (MECO). To establish a detailed cyclostratigraphy we acquired a distinctive iron intensity records by XRF scanning Site 1260 cores. We revise the shipboard composite section, establish a cyclostratigraphy and use the exceptional eccentricity modulated precession cycles for orbital tuning. The new astrochronology revises the age of magnetic polarity chrons C19n to C20n, validates the position of very long eccentricity minima at 40.2 and 43.0 Ma in the orbital solutions, and extends the Astronomically Tuned Geological Time Scale back to 44 Ma. For the first time the new data provide clear evidence for an orbital pacing of the chron C19r event and a likely involvement of the very long eccentricity cycle contributing to the evolution of the MECO.

Age models of ODP Site 184-1143

Here we present an improved astronomical timescale since 5 Ma as recorded in the ODP Site 1143 in the southern South China Sea, using a recently published Asian summer monsoon record (hematite to goethite content ratio, Hm/Gt) and a parallel benthic d18O record. Correlation of the benthic d18O record to the stack of 57 globally distributed benthic d18O records (LR04 stack) and the Hm/Gt curve to the 65°N summer insolation curve is a particularly useful approach to obtain refined timescales. Hence, it constitutes the basis for our effort. Our proposed modifications result in a more accurate and robust chronology than the existing astronomical timescale for the ODP Site 1143. This updated timescale further enables a detailed study of the orbital variability of low-latitude Asian summer monsoon throughout the Plio-Pleistocene. Comparison of the Hm/Gt record with the d18O record from the same core reveals that the oscillations of low-latitude Asian summer monsoon over orbital scales differed considerably from the glacial-interglacial climate cycles. The popular view that summer monsoon intensifies during interglacial stages and weakens during glacial stages appears to be too simplistic for low-latitude Asia. In low-latitude Asia, some strong summer monsoon intervals appear to have also occurred during glacial stages in addition to their increased occurrence during interglacial stages. Vice versa, some notably weak summer monsoon intervals have also occurred during interglacial stages next to their anticipated occurrence during glacial stages. The well-known mid-Pleistocene transition (MPT) is only identified in the benthic d18O record but not in the Hm/Gt record from the same core. This suggests that the MPT may be a feature of high- and middle-latitude climates, possibly determined by high-latitude ice sheet dynamics. For low-latitude monsoonal climate, its orbital-scale variations respond more directly to insolation and are little influenced by high-latitude processes, thus the MPT is likely not recorded. In addition, the Hm/Gt record suggests that low-latitude Asian summer monsoon intensity has a long-term decreasing trend since 2.8 Ma with increased oscillation amplitude. This long-term variability is presumably linked to the Northern Hemisphere glaciation since then.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.