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.

Bulk rock magnetic measurements and thermomagnetic analyses of sediment core GeoB4901-8 from the southeastern flank of the Niger deep-sea fan

Diagenesis has extensively affected the magnetic mineral inventory of organic-rich late Quaternary sediments in the Niger deep-sea fan. Changes in concentration, grain size, and coercivity document modifications of the primary magnetic mineral assemblages at two horizons. The first front, the modern iron redox boundary, is characterized by a drastic decline in magnetic mineral content, coarsening of the grain size spectrum, and reduction in coercivity. Beneath a second front, the transition from the suboxic to the sulfidic anoxic domain, a further but less pronounced decrease in concentration and bulk grain size occurs. Finer grains and higher coercive magnetic constituents substantially increase in the anoxic environment. Low- and high-temperature experiments were performed on bulk sediments and on extracts which have also been examined by X-ray diffraction. Thermomagnetic analyses proved ferrimagnetic titanomagnetites of terrigenous provenance as the principal primary magnetic mineral components. Their broad range of titanium contents reflects the volcanogenic traits of the Niger River drainage areas. Diagenetic alteration is not only a grain size selective process but also critically depends on titanomagnetite composition. Low-titanium compounds are less resistant to diagenetic dissolution. Intermediate titanium content titanomagnetite thus persists as the predominant magnetic mineral fraction in the sulfidic anoxic sediments. At the Fe redox boundary, precipitation of authigenic, possibly bacterial, magnetite is documented. The presence of hydrogen sulfide in the pore water suggests a formation of secondary magnetic iron sulfides in the anoxic domain. Grain size-specific data argue for a gradual development of a superparamagnetic and single-domain iron sulfide phase in this milieu, most likely greigite.

Carbonate and terrigenous content and chemical composition of sediments in the Central Equatorial Pacific

This study addresses changes in the absolute magnitude and spatial geometry of particle flux and export production in a meridional transect across the central equatorial Pacific Ocean's upwelling system during oxygen isotope Stage 11 and Stage 12 and compares these time periods to the current Holocene interglacial system. Temporal and spatial variability in several chemical proxies of export production, and in particular the distributions of Ba, scavenged Al, and P, are studied in a suite of sediment cores gathered along a cross-equator transect at 5°S, 2°S, 0°, 2°N, and 4°N. Because this latitudinal range preserves strong gradients in biogenic particle flux in the modern equatorial Pacific Ocean, we are able to assess variations in the relative magnitude of export production as well as the meridional width of the equatorial system through the late Quaternary glacial/interglacial cycles. During interglacial oxygen isotope Stage 11 the chemical proxies each indicate lower particle flux and export production than during Stage 12. These records are consistent throughout the transect during this time period, but geographic narrowing (during the interglacial) and widening (during the glacial) of the meridional gradient also occurs. Although carbonate concentration varies dramatically through glacial/interglacial cycles at all latitudes studied, the productivity proxies record only minimal glacial/interglacial change at 5°S and 4°N, indicating that the carbonate minima at these latitudes is controlled dominantly by dissolution rather than production. The chemical data indicate that although the spatial geometry of the system during Stages 11 and 12 indicates maximum productivity at the equator during both glacial and interglacial conditions, the absolute magnitude of export production integrated from 5°S to 4°N during Stage 11 was 25-50% less than during Stage 12, and also was 25-50% less than it is now.

Isopollen maps for 18,000 years B.P. of the offshore of Northwest Africa: Evidence for paleowind circulation

Distribution patterns of the most important pollen types from southern European and northwest African source areas for the 18,000 years B.P. time slice are reconstructed from pollen records of 14 well-dated deep-sea cores located between 37° and 9°N and compared with the modern pollen distribution in this area. It is concluded that the belt with maximum African Easterly Jet transport did not shift latitudinally during the last glacial-interglacial transition but remained at about 20°N. Furthermore, it is substantiated that the trade winds did not shift latitudinally during the last glacial-interglacial transition. This evidence is not compatible with an atmospheric circulation model that assumes a zone of surface westerlies in the northern part of northwest Africa. Trade winds during glacial episodes did, however, intensify, especially from about 36° to 24° N.

Eolian records from sediment core RC27-61 in the Arabian Sea

The modern Indian Ocean summer monsoon is driven by differential heating between the Asian continent and the Indian Ocean to the south. This differential heating produces a strong pressure gradient which drives southwest monsoon winds during June, July, and August. Satellite and meteorological observations, aerosol measurements, sediment trap studies, and mineralogical studies indicate an atmospheric mode of transport for modern lithogenic sediments in the northwest Arabian Sea. Analyses of lithogenic grain size and mass accumulation rate (MAR) records from the Owen Ridge indicate that eolian transport has been the primary mode of transport for the past 370 kyr. Visual inspection shows that the MAR record is positively correlated with global ice volume as indicated by the marine delta18O record. In contrast, the grain-size record varies at a much higher frequency, showing little correlation to either the MAR or the delta18O records. Spectral analyses confirm these relationships, indicating that the lithogenic grain-size and MAR records are coherent only over the precession band whereby the grain size leads the MAR by 124° (~8 kyr). We conclude that an eolian transport mechanism is the only mechanism that allows for this phase difference and at the same time is supported by comparison of the grain size and MAR with independent eolian records. We use lithogenic grain size as a paleoclimatic indicator of summer monsoon wind strength and lithogenic MAR as a paleoclimatic indicator of source-area aridity. These interpretations are supported by comparison of the lithogenic records to independent indicators of wind strength (Globigerina bulloides upwelling record) and aridity (a loess record from central China). Such comparisons indicate high coherence and zero phase relationships. Our work supports the findings of previous studies which have documented the link between monsoon strength and the Earth's axial precession cycles. Both the lithogenic MAR and the grain-size records have high coherency with precessional insolation. Maximum lithogenic MAR (source-area aridity) is in phase with delta18O (global ice volume) and leads maximum precessional insolation by 88° (~6 kyr). We attribute this lead to the influence of glacial conditions on the aridity, and therefore the deflation potential, of the source areas. Maximum lithogenic grain size (summer monsoon wind strength) lags maximum precession by 148° (~9 kyr). We attribute this lag both to the influence of global and/or local ice volume and to the availability of latent heat from the southern hemisphere Indian Ocean, the two of which combine to determine the strength of the Indian Ocean monsoon.

Nitrogen isotope variations along the northeast Pacific margin

Glacial-interglacial changes in sedimentary d15N over the last 120 kyr display a remarkably similar pattern in timing and amplitude in core records extending from the denitrification zone in the eastern tropical North Pacific (ETNP), where subsurface denitrification is active, to the Oregon margin, where no denitrification occurs today. Low d15N values (4-6 per mil) generally characterize glacial stages 2 and 4, and higher d15N values (7-10 per mil) are representative of the Holocene, millennial-scale periods within stage 3, and stage 5. The inferred synchroneity of d15N variations along the entire margin implies that the nitrate isotopic signal produced in the oxygen-poor subsurface waters in the ETNP is rapidly advected northward and recorded at sites far beyond the boundaries of the modern denitrification zone. Similar to d15N, primary production indicators (percent Corg, Ba/Al, and percent opal) show glacial-interglacial as well as millennial-scale variations along the NE Pacific margin, with higher primary production during warm periods. However, the relative phasing between d15N and paleoproduction tracers within individual records changes latitudinally. Whereas d15N and primary production vary approximately synchronously in the midlatitudes, production lags d15N in the ETNP by several kiloyears. This lag calls for a new understanding of the processes driving denitrification in the ETNP. We suggest that oxygen input by the Equatorial Undercurrent as well as local organic matter flux controls denitrification rates in the ETNP. Moreover, the differences in relative timing point to a time-transgressive development of upwelling-favorable winds along the NE Pacific margin after the last glaciation, with those in the north developing several kiloyears earlier.

Geochemistry of basement carbonates from ODP Site 801 in the western Pacific Ocean

Many studies argue, based partly on Pb isotopic evidence, that recycled, subducted slabs reside in the mantle source of ocean island basalts (OIB) (Hofmann and White, 1982, doi:10.1016/0012-821X(82)90161-3; Weaver, 1991 doi:10.1016/0012-821X(91)90217-6; Lassiter, and Hauri, 1998, doi:10.1016/S0012-821X(98)00240-4). Such models, however, have remained largely untested against actual subduction zone inputs, due to the scarcity of comprehensive measurements of both radioactive parents (Th and U) and radiogenic daughter (Pb) in altered oceanic crust (AOC). Here, we discuss new, comprehensive measurements of U, Th, and Pb concentrations in the oldest AOC, ODP Site 801, and consider the effect of subducting this crust on the long-term Pb isotope evolution of the mantle. The upper 500 m of AOC at Site 801 shows >4-fold enrichment in U over pristine glass during seafloor alteration, but no net change to Pb or Th. Without subduction zone processing, ancient AOC would evolve to low 208Pb/206Pb compositions unobserved in the modern mantle (Hart and Staudigel, 1989 [Isotopic characterization and identification of recycled components, in: Crust/Mantle Recycling at Convergence Zones, Eds. S.R. Hart, L. Gqlen, NATO ASI Series. Series C: Mathematical and Physical Sciences 258, pp. 15-28, D. Reidel Publishing Company, Dordrecht-Boston, 1989]). Subduction, however, drives U-Th-Pb fractionation as AOC dehydrates in the earth's interior. Pacific arcs define mixing trends requiring 8-fold enrichment in Pb over U in AOC-derived fluid. A mass balance across the Mariana subduction zone shows that 44-75% of Pb but <10% of U is lost from AOC to the arc, and a further 10-23% of Pb and 19-40% of U is lost to the back-arc. Pb is lost shallow and U deep from subducted AOC, which may be a consequence of the stability of phases binding these elements during seafloor alteration: U in carbonate and Pb in sulfides. The upper end of these recycling estimates, which reflect maximum arc and back-arc growth rates, remove enough Pb and U from the slab to enable it to evolve rapidly (<<0.5 Ga) to sources suitable to explain the 208Pb/206Pb isotopic array of OIB, although these conditions fail to simultaneously satisfy the 207Pb/206Pb system. Lower growth rates would require additional U loss (29%) at depths beyond the zones of arc and back-arc magmagenesis, which would decrease upper mantle kappa (232Th/238U) over time, consistent with one solution to the "kappa conundrum" (Elliott et al., 1999, doi:10.1016/S0012-821X(99)00077-1). The net effects of alteration (doubling of l [238U/204Pb]) and subduction (doubling of omega [232Th/204Pb]) are sufficient to create the Pb isotopic signatures of oceanic basalts.

Age determination of sediment cores from the Skagerrak

Paleoceanographic and paleoenvironmental interpretations based on foraminifera, sedimentary data, radiocarbon dates, and stable isotope measurements were derived from two sections in the Skagerrak: a 115-m-thick Holocene marine section drilled onshore at Skagen near the northernmost tip of Jutland, Denmark, and a 9-m piston core from the Skagerrak, north of Skagen. The foraminiferal data show that arctic-subarctic environments in the deep Skagerrak-Kattegat area were succeeded by boreal conditions at 9.6 ka. This was a result of northward migration of the Atlantic polar front and inflow of warm Atlantic water into the area through the Norwegian Channel. A gradual warming of the water masses after 9.6 ka is indicated by the data. Rare foraminifera and high sedimentation rates are found between approximately 8.6 ka and 7.6 ka at both core locations. The modern foraminiferal assemblages of the area were fully established at 7.6 ka indicating that the modern circulation pattern in the Skagerrak-Kattegat after the opening of the English Channel and the Danish Straits was not established before this date. At 5.5 ka a sudden change to coarser sediments (higher-energy environments) and the appearance of the foraminifer Eoeponidella laesoeensis is recorded in the Skagen core. This indicates a rapid change in the hydrography reflecting altered meteorological and hydrographic conditions in the Skagerrak-Kattegat, including a strengthening of the Jutland Current and increased inflow of North Sea water into the Kattegat. The event is interpreted as a response to cooling at the end of the Holocene climatic optimum in late Atlantic time and possibly reflects a rapid cooling event of North Atlantic surface water masses.

Grain size analysis and sedimentation rates of ODP Hole 175-1085A off the Namibian coast

Sediments from the ODP Site 1085A were studied to investigate the impacts of global cooling in the Middle and Late Miocene on the climate in Southwestern Africa. The size composition of the sediment was analysed emphasising the silt fraction. A comparison with the modern grain size distribution and suitable transport processes made it possible to assign specific transport processes to the grain size composition. Three processes are considered for transport of terrigeneous silt: while there was no evidence found for (1) transport by ocean currents, the analyses showed signals of (2) wind transport indicating dry conditions associated with a cool climate and (3) fluvial transport that points to humid and warm conditions. Three climatic phases were defined. The first phase from 13.8 to 11.8 Myr reveals a stable humid climate in Southwest Africa independent of the Antarctic glaciations. During the second phase from 11.8 to 10.4 Myr the regional climate cooled considerably but was not drier. Additionally, the climate during this phase reacted to the Antarctic glaciations. This cooling-trend continued during phase 3 from 10.4 to 9.0 Myr with a significant increase in dust input, pointing to overall drier conditions. However, fluvial transport still remained as the main source.

Late Cretaceous planktonic foraminiferal assemblages from the Exmouth Plateau, ODP Sites 122-762 and 122-763, eastern Indian Ocean

The evolution of planktonic foraminifera during the Late Cretaceous is marked in the Santonian by the disappearance of complex morphotypes (the marginotruncanids), and the contemporary increasing importance and diversification of another group of complex taxa, the globotruncanids. Upper Turonian to lower Campanian planktonic foraminiferal assemblages from Holes 762C and 763B (Ocean Drilling Program, Leg 122, Exmouth Plateau, 47°S palaeolatitude) were studied in detail to evaluate the compositional variations at the genus and species level based on the assumption that, in the Cretaceous oceans as in the modern, any faunal change was associated with changes in the characteristics and the degree of stability of the oceanic surface waters. Three major groups were recognised based on gross morphology, and following the assumption that Cretaceous planktonic foraminifera, although extinct, had life-history strategies comparable to those of modern planktonics: 1 - r-selected opportunists; 2 - k-selected specialists; 3 - r/k intermediate morphotypes which include all genera that display a range of trophic strategies in-between opportunist and specialist taxa. Although planktonic foraminiferal assemblages are characterised by a progressive appearance of complex taxa, this trend is discontinuous. Variation in number of species and specimens within genera has allowed recognition of five discrete intervals each of them reflecting different oceanic conditions based on fluctuations in diversity and abundance of the major morphotypes. Planktonic forms show cyclical fluctuations in diversity and abundance of cold (r-strategists) and warm taxa (k-strategists), perhaps representing alternating phases of unstable conditions (suggesting a weakly stratified upper water column in a mesotrophic environment), and well-stratified surface and near-surface waters (indicating a more oligotrophic environment). Interval 1, middle Turonian to early Coniacian in age, is dominated by the r/k intermediate morphotypes which alternate with r-strategists. These cyclical alternations are used to identify three additional subintervals. Interval 2, aged middle to late Coniacian, is characterised by the increasing number of species and relative abundance of k-strategists. After this maximum diversification the k-strategists show a progressive decrease reaching a minimum value in Interval 3 (early to late Santonian), which corresponds to the extinction of the genus Marginotruncana. In the Interval 4, latest Santonian in age, the k-strategists, represented mainly by the genera Globotruncana, increase again in diversity and abundance. The last Interval 5 (early Campanian) is dominated by juvenile globotruncanids and r-strategists which fluctuate in opposite phase. The positive peak (Interval 2) related to the maximum diversification of warm taxa (k-strategists) in the Coniacian seems to correspond to a warmer episode. It is followed by a marked decrease in the relative abundance of warm taxa (k-strategists crisis) with a minimum in the late Santonian (Interval 3), reflecting a decrease in temperature. Detailed analysis of faunal variations allows the Santonian faunal turnover to be ascribed to a cooling event strong enough to cause the extinction of the marginotruncanids.

Glacial-interglacial sediment record from the North Atlantic

Studies from the subtropical western and eastern Atlantic Ocean, using the 231Pa/230Th ratio as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of meridional overturning circulation (MOC) over the last deglaciation. In this study, we present a compilation of existing and new sedimentary 231Pa/230Th records from North Atlantic cores between 1710 and 4550 m water depth. Comparing sedimentary 231Pa/230Th from different depths provides new insights into the evolution of the geometry and rate of deep water formation in the North Atlantic during the last 20,000 years. The 231Pa/230Th ratio measured in upper Holocene sediments indicates slow water renewal above ?2500 m and rapid flushing below, consistent with our understanding of modern circulation. In contrast, during the Last Glacial Maximum (LGM), Glacial North Atlantic Intermediate Water (GNAIW) drove a rapid overturning circulation to a depth of at least ?3000 m depth. Below ~4000 m, water renewal was much slower than today. At the onset of Heinrich event 1, transport by the overturning circulation declined at all depths. GNAIW shoaled above 3000 m and significantly weakened but did not totally shut down. During the Bølling-Allerød (BA) that followed, water renewal rates further decreased above 2000 m but increased below. Our results suggest for the first time that ocean circulation during that period was quite distinct from the modern circulation mode, with a comparatively higher renewal rate above 3000 m and a lower renewal rate below in a pattern similar to the LGM but less accentuated. MOC during the Younger Dryas appears very similar to BA down to 2000 m and slightly slower below.

Alkenone accumulation and SST reconstruction for the Benguela upwelling zone

A feature of Pliocene climate is the occurrence of "permanent El Niño-like" or "El Padre" conditions in the Pacific Ocean. From the analysis of sediment cores in the modern northern Benguela upwelling, we show that the mean oceanographic state off Southwest Africa during the warm Pliocene epoch was also analogous to that of a persistent Benguela "El Niño". At present these events occur when massive southward flows of warm and nutrient-poor waters extend along the coasts of Angola and Namibia, with dramatic effects on regional marine ecosystems and rainfall. We propose that the persistent warmth across the Pliocene in the Benguela upwelling ended synchronously with the narrowing of the Indonesian seaway, and the early intensification of the Northern Hemisphere Glaciations around 3.0-3.5 Ma. The emergence of obliquity-related cycles in the Benguela sea surface temperatures (SST) after 3 Ma highlights the development of strengthened links to high latitude orbital forcing. The subsequent evolution of the Benguela upwelling system was characterized by the progressive intensification of the meridional SST gradients, and the emergence of the 100 ky cycle, until the modern mean conditions were set at the end of the Mid Pleistocene transition, around 0.6 Ma. These findings support the notion that the interplay of changes in the depth of the global thermocline, atmospheric circulation and tectonics preconditioned the climate system for the end of the warm Pliocene epoch and the subsequent intensification of the ice ages.

Age models and summer sea surface temperature and winter sea ice concentration for the EPILOG-LGM time slice in the Pacific Southern Ocean

Sea surface temperatures and sea-ice extent are the most critical variables to evaluate the Southern Ocean paleoceanographic evolution in relation to the development of the global carbon cycle, atmospheric CO2 variability and ocean-atmosphere circulation. In contrast to the Atlantic and the Indian sectors, the Pacific sector of the Southern Ocean has been insufficiently investigated so far. To cover this gap of information we present diatom-based estimates of summer sea surface temperature (SSST) and winter sea-ice concentration (WSI) from 17 sites in the polar South Pacific to study the Last Glacial Maximum (LGM) at the EPILOG time slice (19,000-23,000 cal. years BP). Applied statistical methods are the Imbrie and Kipp Method (IKM) and the Modern Analog Technique (MAT) to estimate temperature and sea-ice concentration, respectively. Our data display a distinct LGM east-west differentiation in SSST and WSI with steeper latitudinal temperature gradients and a winter sea-ice edge located consistently north of the Pacific-Antarctic Ridge in the Ross sea sector. In the eastern sector of our study area, which is governed by the Amundsen Abyssal Plain, the estimates yield weaker latitudinal SSST gradients together with a variable extended winter sea-ice field. In this sector, sea-ice extent may have reached sporadically the area of the present Subantarctic Front at its maximum LGM expansion. This pattern points to topographic forcing as major controller of the frontal system location and sea-ice extent in the western Pacific sector whereas atmospheric conditions like the Southern Annular Mode and the ENSO affected the oceanographic conditions in the eastern Pacific sector. Although it is difficult to depict the location and the physical nature of frontal systems separating the glacial Southern Ocean water masses into different zones, we found a distinct temperature gradient in latitudes straddled by the modern Southern Subtropical Front. Considering that the glacial temperatures north of this zone are similar to the modern, we suggest that this represents the Glacial Southern Subtropical Front (GSSTF), which delimits the zone of strongest glacial SSST cooling (>4K) to its North. The southern boundary of the zone of maximum cooling is close to the glacial 4°C isotherm. This isotherm, which is in the range of SSST at the modern Antarctic Polar Front (APF), represents a circum-Antarctic feature and marks the northern edge of the glacial Antarctic Circumpolar Current (ACC). We also assume that a glacial front was established at the northern average winter sea ice edge, comparable with the modern Southern Antarctic Circumpolar Current Front (SACCF). During the glacial, this front would be located in the area of the modern APF. The northward deflection of colder than modern surface waters along the South American continent leads to a significant cooling of the glacial Humboldt Current surface waters (4-8K), which affects the temperature regimes as far north as into tropical latitudes. The glacial reduction of ACC temperatures may also result in the significant cooling in the Atlantic and Indian Southern Ocean, thus may enhance thermal differentiation of the Southern Ocean and Antarctic continental cooling. Comparison with temperature and sea ice simulations for the last glacial based on numerical simulations show that the majority of modern models overestimate summer and winter sea ice cover and that there exists few models that reproduce our temperature data rather well.