New age estimates for the Palaeolithic assemblages and Pleistocene succession of Casablanca, Morocco

Marine and aeolian Quaternary sediments from Casablanca, Morocco were dated using the optically stimulated luminescence (OSL) signal of quartz grains. These sediments form part of an extensive succession spanning the Pleistocene, and contain a rich faunal and archaeological record, including an Acheulian lithic assemblage from before the Brunhes–Matayama boundary, and a Homo erectus jaw from younger cave deposits. Sediment samples from the sites of Reddad Ben Ali, Oulad J’mel, Sidi Abderhamane and Thomas Quarries have been dated, in order to assess the upper limits of OSL. The revision of previously measured mammalian tooth enamel electron spin resonance (ESR) dates from the Grotte des Rhinocéros, Oulad Hamida Quarry 1, incorporating updated environmental dose rate measurements and attenuation calculations, also provide chronological constraint for the archaeological material preserved at Thomas Quarries. Several OSL age estimates extend back to around 500,000 years, with a single sample providing an OSL age close to 1 Ma in magnetically reversed sediments. These luminescence dates are some of the oldest determined, and their reliability is assessed using both internal criteria based on stratigraphic consistency, and external lithostratigraphic, morphostratigraphic and independent chronological constraints. For most samples, good internal agreement is observed using single aliquot regenerative-dose OSL measurements, while multiple aliquot additive-dose measurements generally have poorer resolution and consistency. Novel slow-component and component-resolved OSL approaches applied to four samples provide significantly enhanced dating precision, and an examination of the degree of signal zeroing at deposition. A comparison of the OSL age estimates with the updated ESR dates and one U-series date demonstrate that this method has great potential for providing reliable age estimates for sediments of this antiquity. We consider the cause of some slight age inversion observed at Thomas Quarries, and provide recommendations for further luminescence dating within this succession.

Assessing the completeness of optical resetting of quartz OSL in the natural environment

Resetting of previously accumulated optically stimulated luminescence (OSL) signals during transport of sediment is a fundamental requirement for reliable optical dating. The completeness of optical resetting of 46 modern-age quartz samples from a variety of depositional environments was examined. All equivalent dose (De) estimates were View the MathML source, with the majority of aeolian samples View the MathML source, and fluvial samples View the MathML source. The OSL signal of quartz originates from several trap types with different rates of charge loss during illumination. As such, incomplete bleaching may be identifiable as an increase in De from easy-to-bleach through to hard-to-bleach components. For all modern fluvial samples with non-zero De values, SAR De(t) analysis and component-resolved linearly modulated OSL (LM OSL) De estimates showed this to be the case, implying incomplete resetting of previously accumulated charge. LM OSL measurements were also made to investigate the extent of bleaching of the slow components in the natural environment. In aeolian sediments examined, the natural LM OSL was effectively zero (i.e. all components were fully reset). The slow components of modern fluvial samples displayed measurable residual signals up to 15 Gy.

Identification of partial resetting using De as a function of illumination time

Modern age samples from various depositional environments were examined for signal resetting. For 19 modern aeolian/beach samples all De values obtained were View the MathML source, with ∼70% having View the MathML source. For 21 fluvial/colluvial samples, all De values were View the MathML source with ∼80% being View the MathML source. De as a function of illumination (OSL measurement) time (De(t)) plots were examined for all samples. Based on previous laboratory experiments, increases in De(t) were expected for partially reset samples, and constant De(t) for fully reset samples. All aeolian samples, both modern age and additional ‘young’ samples (<1000 years), showed constant (flat) De(t) while all modern, non-zero De, fluvial/colluvial samples showed increasing De(t). ‘Replacement plots’, where a regenerated signal is substituted for the natural, yielded constant (flat) De(t). These findings support strongly the use of De(t) as a method of identifying incomplete resetting in fluvial samples. Potential complicating factors, such as illumination (bleaching) spectrum, thermal instability and component composition are discussed and a series of internal checks on the applicability of the De(t) for each individual aliquot/grain level are outlined.

Potential of the slow component of quartz OSL for age determination of sedimentary samples

The slow component of quartz OSL exhibits a high thermal stability, and, in many of the samples studied, a high dose saturation level (several hundreds or, even thousands, of Grays). These properties suggest that the slow component has potential as a long-range dating tool. Initial attempts have been made to obtain equivalent doses for a number of sedimentary samples. Single- and multiple-aliquot techniques were modified for use with the slow component. The results indicate that there is a good potential for sediment dating, particularly for samples of significant age. Experiments concerning the optical resetting of the slow component suggest that, given its slow optical depletion rate, dating may be restricted to aeolian sediments.

Pollen, plant macrofossil and charcoal records for palaeovegetation reconstruction in the Mediterranean-Black Sea Corridor since the Last Glacial Maximum

New reconstructions of changing vegetation patterns in the Mediterranean-Black Sea Corridor since the Last Glacial Maximum are being produced by an improved biomisation scheme that uses both pollen and plant macrofossil data, in conjunction. Changes in fire regimes over the same interval will also be reconstructed using both microscopic and macroscopic charcoal remains. These reconstructions will allow a diagnosis of the interactions between climate, fire and vegetation on millennial timescales, and will also help to clarify the role of coastline and other geomorphic changes, salinity and impacts of human activities in this region. These new data sets are being produced as a result of collaboration between the Palynology Working Group (WG-2) within the IGCP-521 project and the international Palaeovegetation Mapping Project (BIOME 6000). The main objective of this paper is to present the goals of this cooperation, methodology, including limitations and planned improvements, and to show the initial results of some applications.

Glacial-interglacial changes in dust deposition on the Chinese Loess Plateau

The Chinese Loess Plateau (CLP) contains an extensive record of aeolian deposition through multiple glacial–interglacial cycles. Independent chronologies based on pedostratigraphy, magnetic susceptibility, radiocarbon and luminescence dating were developed for 79 sites and used to estimate aeolian mass accumulation rates (MARs) for marine isotope stages 1–5. The regional median value of MAR for Stage 2 is 310 g/m2/yr compared to an estimate of 65 g/m2/yr for Stage 5. Estimated MARs from individual sites for Stage 2 are approximately 4.3 times greater than MARs for Stage 5 and 2.1 times greater than for Stage 1. MAR values at individual sites are consistently highest in the northwest and lowest in the southwest of the CLP during all marine isotope stages. MARs estimated on sections through loess terraces are consistently higher than MAR estimates at other sites, indicating that local recycling of loess material from exposed river valley deposits has been significant throughout the last 130 kyr. Although the spatial and temporal patterns in MAR are robust, there are uncertainties about the magnitude of these changes due to (a) lack of bulk density measurements and uncertainties in the chronologies for individual sites, (b) site and chronological biases in the sampling used to derive regional estimates, and (c) the unquantified nature of human impact on accumulation rates during the late Holocene. Nevertheless, the records from the CLP pose a number of challenges which could be addressed by numerical models of the palaeo-dust cycle.

Mid-Holocene climates of the Americas: a dynamical response to changed seasonality

Simulations of the climatic response to mid-Holocene (6 ka BP) orbital forcing with two coupled ocean–atmosphere models (FOAM and CSM) show enhancement of monsoonal precipitation in parts of the American Southwest, Central America and northernmost South America during Northern Hemisphere summer. The enhanced onshore flow that brings precipitation into Central America is caused by a northward displacement of the inter-tropical convergence zone, driven by cooling of the equatorial and warming of the northern subtropical and mid-latitude ocean. Ocean feedbacks also enhance precipitation over the American Southwest, although the increase in monsoon precipitation there is largely driven by increases in land-surface temperature. The northward shift in the equatorial precipitation band that causes enhanced precipitation in Central America and the American Southwest has a negative feedback effect on monsoonal precipitation in northern South America. The simulations demonstrate that mid-Holocene aridity in the mid-continent of North America is dynamically linked to the orbitally induced enhancement of the summer monsoon in the American Southwest, with a spatial structure (wet in the Southwest and dry in the mid-continent) similar to that found in strong monsoon years today. Changes in winter precipitation along the west coast of North America, in Central America and along the Gulf Coast, caused by southward-displacement of the westerly storm tracks, indicate that changes in the Northern Hemisphere winter monsoon also play a role in regional climate changes during the mid-Holocene. Although the simulations with FOAM and CSM differ in detail, the general mechanisms and patterns are common to both. The model results thus provide a coherent dynamical explanation for regional patterns of increased or decreased aridity shown by vegetation, lake status and aeolian data from the Americas

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.

Palaeovegetation of China: a pollen data-based synthesis for the mid-Holocene and last glacial maximum.

Pollen data from China for 6000 and 18,000 14C yr bp were compiled and used to reconstruct palaeovegetation patterns, using complete taxon lists where possible and a biomization procedure that entailed the assignment of 645 pollen taxa to plant functional types. A set of 658 modern pollen samples spanning all biomes and regions provided a comprehensive test for this procedure and showed convincing agreement between reconstructed biomes and present natural vegetation types, both geographically and in terms of the elevation gradients in mountain regions of north-eastern and south-western China. The 6000 14C yr bp map confirms earlier studies in showing that the forest biomes in eastern China were systematically shifted northwards and extended westwards during the mid-Holocene. Tropical rain forest occurred on mainland China at sites characterized today by either tropical seasonal or broadleaved evergreen/warm mixed forest. Broadleaved evergreen/warm mixed forest occurred further north than today, and at higher elevation sites within the modern latitudinal range of this biome. The northern limit of temperate deciduous forest was shifted c. 800 km north relative to today. The 18,000 14C yr bp map shows that steppe and even desert vegetation extended to the modern coast of eastern China at the last glacial maximum, replacing today’s temperate deciduous forest. Tropical forests were excluded from China and broadleaved evergreen/warm mixed forest had retreated to tropical latitudes, while taiga extended southwards to c. 43°N.

Pollen-based biomes for Beringia 18,000, 6000 and 0 14C yr B.P.

The objective biomization method developed by Prentice et al. (1996) for Europe was extended using modern pollen samples from Beringia and then applied to fossil pollen data to reconstruct palaeovegetation patterns at 6000 and 18,000 14C yr bp. The predicted modern distribution of tundra, taiga and cool conifer forests in Alaska and north-western Canada generally corresponds well to actual vegetation patterns, although sites in regions characterized today by a mosaic of forest and tundra vegetation tend to be preferentially assigned to tundra. Siberian larch forests are delimited less well, probably due to the extreme under-representation of Larix in pollen spectra. The biome distribution across Beringia at 6000 14C yr bp was broadly similar to today, with little change in the northern forest limit, except for a possible northward advance in the Mackenzie delta region. The western forest limit in Alaska was probably east of its modern position. At 18,000 14C yr bp the whole of Beringia was covered by tundra. However, the importance of the various plant functional types varied from site to site, supporting the idea that the vegetation cover was a mosaic of different tundra types.

Mid-Holocene and glacial-maximum vegetation geography of the northern continents and Africa.

BIOME 6000 is an international project to map vegetation globally at mid-Holocene (6000 14C yr bp) and last glacial maximum (LGM, 18,000 14C yr bp), with a view to evaluating coupled climate-biosphere model results. Primary palaeoecological data are assigned to biomes using an explicit algorithm based on plant functional types. This paper introduces the second Special Feature on BIOME 6000. Site-based global biome maps are shown with data from North America, Eurasia (except South and Southeast Asia) and Africa at both time periods. A map based on surface samples shows the method’s skill in reconstructing present-day biomes. Cold and dry conditions at LGM favoured extensive tundra and steppe. These biomes intergraded in northern Eurasia. Northern hemisphere forest biomes were displaced southward. Boreal evergreen forests (taiga) and temperate deciduous forests were fragmented, while European and East Asian steppes were greatly extended. Tropical moist forests (i.e. tropical rain forest and tropical seasonal forest) in Africa were reduced. In south-western North America, desert and steppe were replaced by open conifer woodland, opposite to the general arid trend but consistent with modelled southward displacement of the jet stream. The Arctic forest limit was shifted slighly north at 6000 14C yr bp in some sectors, but not in all. Northern temperate forest zones were generally shifted greater distances north. Warmer winters as well as summers in several regions are required to explain these shifts. Temperate deciduous forests in Europe were greatly extended, into the Mediterranean region as well as to the north. Steppe encroached on forest biomes in interior North America, but not in central Asia. Enhanced monsoons extended forest biomes in China inland and Sahelian vegetation into the Sahara while the African tropical rain forest was also reduced, consistent with a modelled northward shift of the ITCZ and a more seasonal climate in the equatorial zone. Palaeobiome maps show the outcome of separate, independent migrations of plant taxa in response to climate change. The average composition of biomes at LGM was often markedly different from today. Refugia for the temperate deciduous and tropical rain forest biomes may have existed offshore at LGM, but their characteristic taxa also persisted as components of other biomes. Examples include temperate deciduous trees that survived in cool mixed forest in eastern Europe, and tropical evergreen trees that survived in tropical seasonal forest in Africa. The sequence of biome shifts during a glacial-interglacial cycle may help account for some disjunct distributions of plant taxa. For example, the now-arid Saharan mountains may have linked Mediterranean and African tropical montane floras during enhanced monsoon regimes. Major changes in physical land-surface conditions, shown by the palaeobiome data, have implications for the global climate. The data can be used directly to evaluate the output of coupled atmosphere-biosphere models. The data could also be objectively generalized to yield realistic gridded land-surface maps, for use in sensitivity experiments with atmospheric models. Recent analyses of vegetation-climate feedbacks have focused on the hypothesized positive feedback effects of climate-induced vegetation changes in the Sahara/Sahel region and the Arctic during the mid-Holocene. However, a far wider spectrum of interactions potentially exists and could be investigated, using these data, both for 6000 14C yr bp and for the LGM.

Monsoon changes for 6000 years ago: results of 18 simulations from the Paleoclimate Modeling Intercomparison Project (PMIP).

Amplification of the northern hemisphere seasonal cycle of insolation during the mid-Holocene causes a northward shift of the main regions of monsoon precipitation over Africa and India in all 18 simulations conducted for the Paleoclimate Modeling Intercomparison Project (PMIP). Differences among simulations are related to differences in model formulation. Despite qualitative agreement with paleoecological estimates of biome shifts, the magnitude of the monsoon increases over northern Africa are underestimated by all the models.

Land surface feedbacks and palaeomonsoons in northern Africa

We ran a sequence of climate model experiments for 6000 years ago, with land-surface conditions based on a realistic map of palaeovegetation, lakes and wetlands, to quantify the effects of land-surface feedbacks in the Saharan region. Vegetation-induced albedo and moisture flux changes produced year-round warming, forced the monsoon to 17°–25°N two months earlier, and shifted the precipitation belt ≈300 km northwards compared to the effects of orbital forcing alone. The addition of lakes and wetlands produced localised changes in evaporation and precipitation, but caused no further extension of the monsoon belt. Diagnostic analyses with biome and continental hydrology models showed that the combined land-surface feedbacks, although substantial, could neither maintain grassland as far north as observed (≈26°N) nor maintain Lake “MegaChad” (330,000 km²).

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.