Deforesting Europe: pollen-inferred Holocene land cover change using the 'pseudobiomization' approach

Maps of continental-scale land cover are utilized by a range of diverse users but whilst a range of products exist that describe present and recent land cover in Europe, there are currently no datasets that describe past variations over long time-scales. User groups with an interest in past land cover include the climate modelling community, socio-ecological historians and earth system scientists. Europe is one of the continents with the longest histories of land conversion from forest to farmland, thus understanding land cover change in this area is globally significant. This study applies the pseudobiomization method (PBM) to 982 pollen records from across Europe, taken from the European Pollen Database (EPD) to produce a first synthesis of pan-European land cover change for the period 9000 BP to present, in contiguous 200 year time intervals. The PBM transforms pollen proportions from each site to one of eight land cover classes (LCCs) that are directly comparable to the CORINE land cover classification. The proportion of LCCs represented in each time window provides a spatially aggregated record of land cover change for temperate and northern Europe, and for a series of case study regions (western France, the western Alps, and the Czech Republic and Slovakia). At the European scale, the impact of Neolithic food producing economies appear to be detectable from 6000 BP through reduction in broad-leaf forests resulting from human land use activities such as forest clearance. Total forest cover at a pan-European scale moved outside the range of previous background variability from 4000 BP onwards. From 2200 BP land cover change intensified, and the broad pattern of land cover for preindustrial Europe was established by 1000 BP. Recognizing the timing of anthropogenic land cover change in Europe will further the understanding of land cover-climate interactions, and the origins of the modern cultural landscape.

Abrupt vegetation change after the Late Quaternary megafaunal extinction in southeastern Australia

A substantial extinction of megafauna occurred in Australia between 50 and 45 kyr ago, a period that coincides with human colonization of Australia. Large shifts in vegetation also occurred around this time, but it is unclear whether the vegetation changes were driven by the human use of fire-and thus contributed to the extinction event-or were a consequence of the loss of megafaunal grazers. Here we reconstruct past vegetation changes in southeastern Australia using the stable carbon isotopic composition of higher plant wax n-alkanes and levels of biomass burning from the accumulation rates of the biomarker levoglucosan from a well-dated sediment core offshore from the Murray-Darling Basin. We find that from 58 to 44 kyr ago, the abundance of plants with the C-4 carbon fixation pathway was generally high-between 60 and 70%. By 43 kyr ago, the abundance of C-4 plants dropped to 30% and biomass burning increased. This transient shift lasted for about 3,000 years and came after the period of human arrival and directly followed megafauna extinction at 48.9-43.6 kyr ago. We conclude that the vegetation shift was not the cause of the megafaunal extinction in this region. Instead, our data are consistent with the hypothesis that vegetation change was the consequence of the extinction of large browsers and led to the build-up of fire-prone vegetation in the Australian landscape.

A Late Pleistocene sea level stack

Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0-430 ka and five records from 0-798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the d18O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic d18O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the d18O of seawater.

Late Weichselian fluvial evolution of the southern Kara Sea

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1-12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network. These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel-levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land. The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses. Asymmetrical submarine channel-levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel-levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.

Physical/chemical properties and stable oxygen and carbon isotope records of Holocene to late Pleistocene sediments in the Gela Basin, central Mediterranean Sea

The Holocene Twin Slides form the most recent of recurrent mass wasting events along the NE portion of Gela Basin within the Sicily Channel, central Mediterranean Sea. Here, we present new evidence on the morphological evolution and stratigraphic context of this coeval slide complex based on deepdrilled sediment sequences providing a >100 ka paleo-oceanographic record. Both Northern (NTS) and Southern Twin Slide (STS) involve two failure stages, a debris avalanche and a translational slide, but are strongly affected by distinct preconditioning factors linked to the older and buried Father Slide. Core-acoustic correlations suggest that sliding occurred along sub-horizontal weak layers reflecting abrupt physical changes in lithology or mechanical properties. Our results show further that headwall failure predominantly took place along sub-vertical normal faults, partly through reactivation of buried Father Slide headscarps.

Holocene flood frequency across the Central Alps

The frequency of large-scale heavy precipitation events in the European Alps is expected to undergo substantial changes with current climate change. Hence, knowledge about the past natural variability of floods caused by heavy precipitation constitutes important input for climate projections. We present a comprehensive Holocene (10,000 years) reconstruction of the flood frequency in the Central European Alps combining 15 lacustrine sediment records. These records provide an extensive catalog of flood deposits, which were generated by flood-induced underflows delivering terrestrial material to the lake floors. The multi-archive approach allows suppressing local weather patterns, such as thunderstorms, from the obtained climate signal. We reconstructed mainly late spring to fall events since ice cover and precipitation in form of snow in winter at high-altitude study sites do inhibit the generation of flood layers. We found that flood frequency was higher during cool periods, coinciding with lows in solar activity. In addition, flood occurrence shows periodicities that are also observed in reconstructions of solar activity from 14C and 10Be records (2500-3000, 900-1200, as well as of about 710, 500, 350, 208 (Suess cycle), 150, 104 and 87 (Gleissberg cycle) years). As atmospheric mechanism, we propose an expansion/shrinking of the Hadley cell with increasing/decreasing air temperature, causing dry/wet conditions in Central Europe during phases of high/low solar activity. Furthermore, differences between the flood patterns from the Northern Alps and the Southern Alps indicate changes in North Atlantic circulation. Enhanced flood occurrence in the South compared to the North suggests a pronounced southward position of the Westerlies and/or blocking over the northern North Atlantic, hence resembling a negative NAO state (most distinct from 4.2 to 2.4 kyr BP and during the Little Ice Age). South-Alpine flood activity therefore provides a qualitative record of variations in a paleo-NAO pattern during the Holocene. Additionally, increased South Alpine flood activity contrasts to low precipitation in tropical Central America (Cariaco Basin) on the Holocene and centennial time scale. This observation is consistent with a Holocene southward migration of the Atlantic circulation system, and hence of the ITCZ, driven by decreasing summer insolation in the Northern hemisphere, as well as with shorter-term fluctuations probably driven by solar activity.

Accumulation of particulate organic carbon at the Eurasian continental margin during late Quaternary times

Data on the amount and composition of organic carbon were determined in sediment cores from the Kara and Laptev Sea continental margin, representing oxygen isotope stages 1-6. The characterization of organic matter is based on hydrogen index (HI) values, n-alkanes and maceral composition, indicating the predominance of terrigenous organic matter through space and time. The variations in the amount and composition of organic carbon are mainly influenced by changes in fluvial sediment supply, Atlantic water inflow, and continental ice sheets. During oxygen isotope stage (OIS) 6, high organic carbon contents in sediments from the Laptev Sea and western East Siberian Sea continental margin were probably caused by the increased glacial erosion and further transport in the eastward-flowing boundary current along the continental margin. During OIS 5 and early OIS 3, some increased amounts of marine organic matter were preserved in sediments east of the Lomonosov Ridge, suggesting an influence of nutrient-rich Pacific waters. During OIS 2, terrigenous organic carbon supply was increased along the Barents and western Kara Sea continental margin caused by extended continental ice sheets in the Barents Sea (Svalbard to Franz Josef Land) area and increased glacial erosion. Along the Laptev Sea continental margin, on the other hand, the supply of terrigenous (organic) matter was significantly reduced due to the lack of major ice sheets and reduced river discharge. Towards the Holocene, the amount of total organic carbon (TOC) increased along the Kara and Laptev Sea continental margin, reaching average values of up to 0.5 g C/cm**2/ky. Between about 8 and 10 ka (9 and 11 Cal ka), i.e., during times when the inner shallow Kara and Laptev seas became largely flooded for the first time after the Last Glacial Maximum, maximum supply of terrigenous organic carbon occurred, which is related to an increase in coastal erosion and Siberian river discharge. During the last 8000 years, the increased amount of marine organic carbon preserved in the sediments from the Kara and Laptev Sea continental margin is interpreted as a result of the intensification of Atlantic water inflow along the Eurasian continental margin.

Late pleistocene and holocene sedimentation in the central and eastern Persian Gulf

The sandfraction of the sediment was analysed in five cores, taken from 65 m water depth in the central and eastern part of the Persian Gulf. The holocene marls are underlayn by aragonite muds, which are probably 10-11,000 years old. 1. The cores could be subdivided into coarse grained and fine grained layers. Sorting is demonstrated by the following criteria: With increasing median values of the sandfraction - the fine grained fraction decreases within each core; - the median of each biogenic component, benthonic as well as planktonic, increases; - the median of the relict sediment, which in core 1179 was carried upward into the marl by bioturbation, increases; - the percentages of pelecypods, gastropods, decapods and serpulid worms in the sandfraction increase, the percentages of foraminifera and ostracods decrease; - the ratios of pelecypods to foraminifera and of decapods to ostracods increase; - the ratios of benthonic molluscs to planktonic molluscs (pteropods) and of benthonic foraminifera to planktonic foraminifera increase (except in core 1056 and 1179); - the ratio of planktonic molluscs (pteropods) to planktonic foraminifera increases; - the globigerinas without orbulinas increase, the orbulinas decrease in core 1056. Different settling velocities of these biogenic particles help in better understanding the results : the settling velocities, hence the equivalent hydrodynamic diameters, of orbulinas are smaller than those of other globigerinas, those of planktonic foraminifera are smaller than those of planktonic molluscs, those of planktonic molluscs are smaller than those of benthonic molluscs, those of pelecypods are smaller than those of gastropods. Bioturbation could not entirely distroy this "grain-size-stratification". Sorting has been stronger in the coarse layers than in the finer ones. As a cause variations in the supply of terrigenous material at constant strength of tidal currents is suggested. When much terrigenous material is supplied (large contents of fine grained fraction) the sedimentation rates are high: the respective sediment surface is soon covered and removed from the influence of tidal currents. When, however, the supply of terrigenous material is small, more sandy material is taken away in all locations within the influence of terrigenous supply. Thus the biogenic particles in the sediment do not only reflect the organic production, but also the influence of currents. 2. There is no parameter present in all cores that is independently variable from grain size and can be used for stratigraphic correlation. The two cores from the Strait of Hormus were correlated by their sequences of coarse and fine grained layers. 3. The sedimentation rates of terrigenous material, of total planktonic and benthonic organisms and of molluscs, foraminifera, echinoids and ophiuroids are shown in table 1 (total sediment 6.3-75.5 cm/1000 yr, biogenic carbonate 1.9-3.6 cm/1000 yr). The sedimentation rates of benthonic organisms are nearly the same in the cores of the Strait of Hormus, whereas near the Central Swell they are smaller. In the upper parts of the two cores of the Strait of Hormus sedimentation rates are higher than in the deeper parts, where higher median values point to stronger reworking. 4. The sequence of coarse and fine grained intervals in the two cores of the Hormus Strait, attributed to variations in climate, as well as the increase of terrigenous supply from the deeper to the upper parts of the cores, agrees with the descriptions in the literature of the post Pleistocene climate as becoming more humid. The rise of sea level is sedimentologically not measurable in the marly sediments - except perhaps for the higher content of echinoids in the lower part of core 1056. These may be attributed to the influence of a migrating wave-base. 5. The late Pleistocene aragonite mud is very fine grained (> 50%< 2 p) and poor in fossils (0.5-1.8%) biogenic particles of total sediment. The sand fraction consists almost entirely of white clumps, c. 0.1 mm in diameter (1177), composed of aragonite needles and of detrital minerals with the same size (1201). The argonite mud was probably not formed in situ, because the water depth at time of formation was at most 35 m at least 12 m. The sorting of the sediment (predominance of the fine grained sand), the absence of larger biogenic components and of pellets, c. 0.2-0.5 mm in diameter, which are typical for Recent and Pleistocene locations of aragonite formation, as well as the sedimentological conditions near the sampling points, indicate rather a transport of aragonite mud from an area of formation in very shallow waters. Sorting as well as lenticular fabric in core 1201 point to sedimentation within the influence of currents. During alternating sedimentation - and reworking processes the aragonitic matrix was separated from the silt - and sand-sized minerals. The lenses grade into touches because of bioturbation. 6. In core 1056 D2 from Hormus Bay the percentages of organic carbon, total nitrogen and total carbonate were determined. With increasing amounts of smaller grain sizes the content of organic matter increases, whereas the amount of carbonate decreases. The amounts of organic carbon and of nitrogen decrease with increasing depth, probably due to early-diagenetic decomposition processes. Most of the total nitrogen is of organic origin, only about 10% may well be inorganically fixed as ammonium-nitrogen. In the upper part of the core the C/N-ratio increases with increasing depth. This may be connected with a stronger decomposition of nitrogen-containing organic compounds. The general decrease of the C/N-ratios in the lower part of the core may be explained by the relative increase of inorganically fixed ammonium-nitrogen with decreasing content of organic matter.

Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.

Late Glacial and Holocene data from sediment cores MD99-2236 and MD99-2256

The Laurentide Ice Sheet (LIS) was a large, dynamic ice sheet in the early Holocene. The glacial events through Hudson Strait leading to its eventual demise are recorded in the well-dated Labrador shelf core, MD99-2236 from the Cartwright Saddle. We develop a detailed history of the timing of ice-sheet discharge events from the Hudson Strait outlet of the LIS during the Holocene using high-resolution detrital carbonate, ice rafted detritus (IRD), d18O, and sediment color data. Eight detrital carbonate peaks (DCPs) associated with IRD peaks and light oxygen isotope events punctuate the MD99-2236 record between 11.5 and 8.0 ka. We use the stratigraphy of the DCPs developed from MD99-2236 to select the appropriate DeltaR to calibrate the ages of recorded glacial events in Hudson Bay and Hudson Strait such that they match the DCPs in MD99-2236. We associate the eight DCPs with H0, Gold Cove advance, Noble Inlet advance, initial retreat of the Hudson Strait ice stream (HSIS) from Hudson Strait, opening of the Tyrrell Sea, and drainage of glacial lakes Agassiz and Ojibway. The opening of Foxe Channel and retreat of glacial ice from Foxe Basin are represented by a shoulder in the carbonate data. DeltaR of 350 years applied to the radiocarbon ages constraining glacial events H0 through the opening of the Tyrell Sea provided the best match with the MD99-2236 DCPs; DeltaR values and ages from the literature are used for the younger events. A very close age match was achieved between the 8.2 ka cold event in the Greenland ice cores, DCP7 (8.15 ka BP), and the drainage of glacial lakes Agassiz and Ojibway. Our stratigraphic comparison between the DCPs in MD99-2236 and the calibrated ages of Hudson Strait/Bay deglacial events shows that the retreat of the HSIS, the opening of the Tyrell Sea, and the catastrophic drainage of glacial lakes Agassiz and Ojibway at 8.2 ka are separate events that have been combined in previous estimates of the timing of the 8.2 ka event from marine records. SW Iceland shelf core MD99-2256 documents freshwater entrainment into the subpolar gyre from the Hudson Strait outlet via the Labrador, North Atlantic, and Irminger currents. The timing of freshwater release from the LIS Hudson Strait outlet in MD99-2236 matches evidence for freshwater forcing and LIS icebergs carrying foreign minerals to the SW Iceland shelf between 11.5 and 8.2 ka. The congruency of these records supports the conclusion of the entrainment of freshwater from the retreat of the LIS through Hudson Strait into the subpolar gyre and provides specific time periods when pulses of LIS freshwater were present to influence climate.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Clay-mineral composition and biogenic opal content of Late Neogene to Quaternary sediments of the Antarctic Peninsula region

Clay-mineral composition and biogenic opal content in upper Miocene to Quaternary drift sediments recovered at two Ocean Drilling Program (ODP) sites from the continental rise in the Bellingshausen Sea had been analyzed in order to reconstruct the climatic and glacial history of the Antarctic Peninsula. The clay mineral composition at both sites is dominated by smectite, illite, and chlorite, and alternates between a smectite-enriched and a chlorite-enriched assemblage throughout the last 9.3 my. The spatial distribution of clay minerals in Holocene sediments west of the Antarctic Peninsula facilitates the identification of particular source areas, and thus the reconstruction of transport pathways. The similarity to clay mineral variations reported from upper Quaternary sequences suggests that the short-term clay-mineralogical fluctuations in the ODP cores reflect glacial-interglacial cyclicity. Thus, repeated ice advances and retreats in response to a varying size of the Antarctic Peninsula ice cap are likely to have occurred throughout the late Neogene and Quaternary. The clay minerals in the drift sediments exhibit only slight long-term variations, which are caused by local changes in glacial erosion and in supply of source rocks, rather than by major climatic changes. The opal records at the ODP sites are dominated by long-term variations since the late Miocene. We infer that the opal content in the drift sediments, although it is influenced by dissolution in the water column and the sediment column and by the burial with lithogenic detritus, provides a signal of paleoproductivity. Because the annual sea-ice coverage is regarded as the main factor controlling biological productivity, the opal signal helps to reconstruct paleoceanographic changes in the Bellingshausen Sea. Slightly enhanced opal deposition during the late Miocene indicates slightly warmer climatic conditions in the Antarctic Peninsula area than at present. During the early Pliocene, enhanced opal deposition in the Pacific sector of the Southern Ocean and coinciding high opal concentrations in sedimentary sequences from the Atlantic and Indian sectors document a strong reduction of sea-ice cover and relatively warm climatic conditions. Thereby, the early onset of the Pliocene warmth in the Bellingshausen Sea points to a positive feedback of regional Antarctic climate on the global thermohaline circulation. A decrease of opal deposition between 3.1 and 2.6 Ma likely reflects sea-ice expansion in response to reduced supply of northern-sourced deep-waters to the Southern Ocean, caused by the onset of Northern Hemisphere glaciation. Throughout the Quaternary, a relatively constant level of opal deposition on the Antarctic continental margin indicates relatively stable climatic conditions.