Multiple Glacial-interglacial Cycles Research Articles

Effective diffusivity of sulfuric acid in Antarctic ice cores

Abstract. Deposition of sulfuric acid in ice cores is important both for understanding past volcanic activity and for synchronizing ice core timescales. Sulfuric acid has a low eutectic point, so it can potentially exist in liquid at grain boundaries and veins, accelerating chemical diffusion. A high effective diffusivity would allow post-depositional diffusion to obscure the climate history and the peak matching among older portions of ice cores. Here, we use records of sulfate from the European Project of Ice Coring in Antarctica (EPICA) Dome C (EDC) ice core to estimate the effective diffusivity of sulfuric acid in ice. We focus on EDC because multiple glacial–interglacial cycles are preserved, allowing analysis for long timescales and deposition in similar climates. We calculate the mean concentration gradient and the width of prominent volcanic events, and analyze the evolution of each with depth and age. We find the effective diffusivities for interglacial and glacial maximums to be 5±3×10-9 m2 a−1, an order of magnitude lower than a previous estimate derived from the Holocene portion of EDC (Barnes et al., 2003). The effective diffusivity may be even smaller if the bias from artificial smoothing from the sampling is accounted for. Effective diffusivity is not obviously affected by the ice temperature until about −10 ∘C, 3000 m depth, which is also where anomalous sulfate peaks begin to be observed (Traversi et al., 2009). Low effective diffusivity suggests that sulfuric acid is not readily diffusing in liquid-like veins in the upper portions of the Antarctic Ice Sheet and that records may be preserved in deep, old ice if the ice temperature remains well below the pressure melting point.

Elevated Southern Hemisphere moisture availability during glacial periods.

Late Pleistocene ice-age climates are routinely characterized as having imposed moisture stress on low- to mid-latitude ecosystems1-5. This idea is largely based on fossil pollen evidence for widespread, low-biomass glacial vegetation, interpreted as indicating climatic dryness6. However, woody plant growth is inhibited under low atmospheric CO2(refs. 7,8), so understanding glacial environments requires the development of new palaeoclimate indicators that are independent of vegetation9. Here we show that, contrary to expectations, during the past 350 kyr, peaks in southern Australian climatic moisture availability were largely confined to glacial periods, including the Last Glacial Maximum, whereas warm interglacials were relatively dry. By measuring the timing of speleothem growth in the Southern Hemisphere subtropics, which today has a predominantly negative annual moisture balance, we developed a record of climatic moisture availability that is independent of vegetation and extends through multiple glacial-interglacial cycles. Our results demonstrate that a cool-moist response is consistent across the austral subtropics and, in part, may result from reduced evaporation under cool glacial temperatures. Insofar as cold glacial environments in the Southern Hemisphere subtropics have been portrayed as uniformly arid3,10,11, our findings suggest that their characterization as evolutionary or physiological obstacles to movement and expansion of animal, plant and, potentially, human populations10 should be reconsidered.

The Coupled Ice Sheet–Earth System Model Bern3D v3.0

Abstract Understanding climate variability from millennial to glacial–interglacial time scales remains challenging due to the complex and nonlinear feedbacks between ice, ocean, sediments, biosphere, and atmosphere. Complex climate models generally struggle to dynamically and comprehensively simulate such long time periods as a result of the large computational costs. Here, we therefore coupled a dynamical ice sheet model to the Bern3D Earth system model of intermediate complexity, which allows for simulating multiple glacial–interglacial cycles. The performance of the model is first validated against modern observations and its response to abrupt perturbations, such as atmospheric CO2 changes and North Atlantic freshwater hosing, is investigated. To further test the fully coupled model, the climate evolution over the entire last glacial cycle is explored in a transient simulation forced by variations in the orbital configuration and greenhouse gases and aerosols. The model simulates global mean surface temperature in fair agreement with reconstructions, exhibiting a gradual cooling trend since the last interglacial that is interrupted by two more rapid cooling events during the early Marine Isotope Stage (MIS) 4 and Last Glacial Maximum (LGM). Simulated Northern Hemispheric ice sheets show pronounced variability on orbital time scales, and ice volume more than doubles from MIS3 to the LGM in good agreement with recent sea level reconstructions. At the LGM, the Atlantic overturning has a strength of about 14 Sv (1 Sv ≡ 106 m3 s−1), which is a reduction by about one-quarter compared to the preindustrial. We thus demonstrate that the new coupled model is able to simulate large-scale aspects of glacial–interglacial cycles.

Glacial–Interglacial Cycles and Early Human Evolution in China

China is a crucial region for investigating the relationship between climate change and hominin evolution across diverse terrestrial ecosystems. With the continuous development of palaeoclimatology, chronology, and archaeology, the environmental and hominin record of the Early and Middle Pleistocene in China is steadily accumulating, shedding light on the effects of climate change on the distribution of early human settlements and population dynamics. However, the migration and dispersal of these early humans within long-term climate fluctuations and their underlying mechanisms remain to be clarified. Based on the spatial-temporal distribution of 95 Early to Middle Pleistocene archaeological sites in China, we found that intensified hominin activities gradually shifted southward under the influence of multiple glacial–interglacial cycles. The frequent bidirectional movements of these early humans between north and south were assumed as follows. During glacial periods, hominins living in North China migrated to southern areas, while inter-glacial periods witnessed the northward expansion of hominins inhabiting South China. Among all the potential driving mechanisms, we suggest that the available resources in terrestrial ecosystems may be the most fundamental factor. Combined with paleoenvironmental and archaeological records, we provide an Asian perspective for a better understanding of how the glacial–interglacial cycles shaped early human evolution.

Valley morphology and Quaternary seismic stratigraphy of the Manicouagan impact crater lake (Eastern Canada)

ABSTRACTLakes in formerly glaciated regions can provide valuable paleoclimate archives. Lake Manicouagan (Manikuakan, according to Innu toponymy), formed in the basin of the ~214‐Ma Manicouagan impact crater of eastern Québec, is a key area for reconstructing long‐term environmental change, as it was directly affected by the Pleistocene glaciations and the waxing and waning of the Laurentide Ice Sheet. Here, we present high‐resolution seismic data revealing an overdeepened bedrock valley filled with a sedimentary sequence. We assess its potential to serve as a paleoclimate archive. The varying shape of the overdeepened valley indicates complex erosional processes. A lower narrow V‐shaped gorge is indicative of either pressurized subglacial meltwater or pre‐Quaternary fluvial erosion, or a combination of both. Three scenarios are discussed regarding deposition of the sedimentary sequence: (i) deposition only during and after retreat of the last glacial episode; (ii) deposition during multiple glacial–interglacial cycles; and (iii) deposition mainly during a subglacial lake stage. We suggest subglacial followed by proglacial sedimentation as the most probable scenario for deposition of the sedimentary succession. We recommend considering the sediments of Lake Manicouagan as a paleoclimate archive reaching back at least to 7.5 ka, but the lake also probably contains sediments deposited before the last deglacial period.

Environmental and climate evolution in the Southwest USA since the last interglacial deduced from the pollen record from Stoneman lake, Arizona

Long and continuous lake sedimentary records offer enormous potential for interpreting the paleoenvironmental histories of the past and for understanding how terrestrial environments might adapt in the context of current global warming. However, sedimentary records that contain multiple glacial-interglacial cycles are scarce in continental basins. An ∼80 m sediment core was recently obtained from Stoneman Lake (STL), Arizona, containing a unique record of the last ∼1.3 Ma. Here we show a detailed pollen study of the topmost ∼10 m of the record, covering the last climatic cycle since the Last Interglacial period (MIS5-MIS1; last ∼130,000 years = 130 kyr), with the goal of broadening our knowledge of the paleoenvironmental history of the arid North American Southwest in the past. The STL pollen record shows that the MIS5e interglacial was the warmest period of the last 130 kyr. This is deduced by the abundance of pollen types from plants that today exist at lower elevations that occurred around the STL at that time. These include Pinus edulis and other associated low elevation thermophilous plants such as Juniperus, Ambrosia, Amaranthaceae, Asteraceae and Artemisia. Climate cooled rapidly and dramatically at the MIS5-4 boundary, which triggered a displacement of forest species towards lower elevation, causing P. ponderosa to occupy the study area. MIS3 was characterized by relatively warmer climate conditions with 3 prominent climatic oscillations (MIS3a, b and c). The coldest conditions were reached during MIS2 (LGM), when a ∼1000 m displacement towards lower elevations of the subalpine forest species relative to present is observed. This is deduced by the highest abundance of Picea (∼20–25%) and Abies in the STL record, indicating their occurrence in the study area. Warming during the last deglaciation is evidenced by a shift of vegetation towards higher altitudes and the development of a montane forest composed mainly of Pinus ponderosa and Quercus replacing the LGM subalpine species. This montane forest remained abundant throughout the Holocene. This study shows that orbital-scale climate changes (mainly precession and eccentricity changes) forced vegetation and lake-level oscillations, documenting that insolation had a main role in controlling environmental change in this area. Climate projections of enhanced warming predict that P. edulis and Juniperus forest species will occupy the study area in the near future.

Ancient marine sediment DNA reveals diatom transition in Antarctica

Antarctica is one of the most vulnerable regions to climate change on Earth and studying the past and present responses of this polar marine ecosystem to environmental change is a matter of urgency. Sedimentary ancient DNA (sedaDNA) analysis can provide such insights into past ecosystem-wide changes. Here we present authenticated (through extensive contamination control and sedaDNA damage analysis) metagenomic marine eukaryote sedaDNA from the Scotia Sea region acquired during IODP Expedition 382. We also provide a marine eukaryote sedaDNA record of ~1 Mio. years and diatom and chlorophyte sedaDNA dating back to ~540 ka (using taxonomic marker genes SSU, LSU, psbO). We find evidence of warm phases being associated with high relative diatom abundance, and a marked transition from diatoms comprising <10% of all eukaryotes prior to ~14.5 ka, to ~50% after this time, i.e., following Meltwater Pulse 1A, alongside a composition change from sea-ice to open-ocean species. Our study demonstrates that sedaDNA tools can be expanded to hundreds of thousands of years, opening the pathway to the study of ecosystem-wide marine shifts and paleo-productivity phases throughout multiple glacial-interglacial cycles.

Middle-Late Pleistocene Eastern Mediterranean nutricline depth and coccolith preservation linked to Monsoon activity and Atlantic Meridional Overturning Circulation

The eastern Mediterranean Sea lies under the influence of high- and low-latitude climatic systems. The northern part of the basin is affected by Atlantic depressions and continental and polar air masses that promote intermediate and deep-water formation. The southern part is influenced by subtropical conditions and monsoon activity. Monsoon intensification results in enhanced freshwater discharge from the Nile River and other (now dry) systems along the North African margin. This freshwater influx into the Mediterranean Sea reduces surface water buoyancy loss. Disentangling the influences of these diverse climatic forcings is hindered by inherent proxy data limitations and by interactions between the climatic forcings. Here we use a wealth of published and new paleoclimate records across Termination II to understand the impacts of the higher latitude and subtropical/monsoon climate influences on coccolithophore ecology and holococcolith preservation in Aegean Sea sediment core LC21. We then use these findings to interpret coccolith assemblage variations at Ocean Drilling Program Site 967 (located nearby LC21, at the Eratosthenes Seamount) during multiple glacial-interglacial cycles across the Middle Pleistocene (marine isotopic stages 14–9). The LC21 analysis suggests that holococcolith preservation was enhanced during Heinrich Stadial 11 (∼133 ka) and cold spell C26 (∼119 ka). These two events have been previously linked to cold conditions in the North Atlantic and Atlantic Meridional Overturning Circulation weakening. We propose that associated atmospheric perturbations over the Mediterranean Sea promoted deep-water formation, and thus holococcolith preservation. Similarly, in the Middle Pleistocene (MIS 14-9) of Site 967, we observe temporal coincidence between ten episodes of enhanced holococcolith preservation and episodes of Atlantic Meridional Overturning Circulation slowdown. In Site 967, we also identified repeated fluctuations in placoliths and in Florisphaera profunda, which indicate nutricline depth variations. The development of a deep chlorophyll maximum is associated with the North Africa and wet phases, as recently observed using elemental proxy records at Site 967, during the deposition of sapropel layers. A further deep chlorophyll maximum development is identified during MISs 12 and 10, as a result of pycnocline and nutricline shoaling within the lower part of the photic zone due to glacial sea-level lowering and water mass transport reduction at both the Gibraltar and Sicily Straits. Finally, enhanced holococcolith preservation during cold/dry events is clearly correlated to weakened monsoon activity in both Africa and Asia.

Evaluating refugia in recent human evolution in Africa.

Homo sapiens have adapted to an incredible diversity of habitats around the globe. This capacity to adapt to different landscapes is clearly expressed within Africa, with Late Pleistocene Homo sapiens populations occupying savannahs, woodlands, coastlines and mountainous terrain. As the only area of the world where Homo sapiens have clearly persisted through multiple glacial-interglacial cycles, Africa is the only continent where classic refugia models can be formulated and tested to examine and describe changing patterns of past distributions and human phylogeographies. The potential role of refugia has frequently been acknowledged in the Late Pleistocene palaeoanthropological literature, yet explicit identification of potential refugia has been limited by the patchy nature of palaeoenvironmental and archaeological records, and the low temporal resolution of climate or ecological models. Here, we apply potential climatic thresholds on human habitation, rooted in ethnographic studies, in combination with high-resolution model datasets for precipitation and biome distributions to identify persistent refugia spanning the Late Pleistocene (130–10 ka). We present two alternate models suggesting that between 27% and 66% of Africa may have provided refugia to Late Pleistocene human populations, and examine variability in precipitation, biome and ecotone distributions within these refugial zones.This article is part of the theme issue ‘Tropical forests in the deep human past’.

A Long-Term Record of Quaternary Facies Patterns and Palaeonvironmental Trends from the Po Plain (NE Italy) as Revealed by Bio-Sedimentary Data

Understanding Quaternary dynamics of delta-coastal plains across multiple glacial-interglacial cycles in the Milankovitch band (~100 kyrs) is crucial to achieve a robust evaluation of possible environmental response to future climate-change scenarios. In this work, we document the long-term bio-sedimentary record of core 204 S16 (~205 m long), which covers a wide portion of the post-MPR (Mid-Pleistocene Revolution) interval, taking advantage of the highly subsiding context of the SE Po Plain (NE Italy). Detailed facies characterization through an integrated sedimentological and meiofauna (benthic foraminifers and ostracods) approach allowed for the identification of a repetitive pattern of alluvial deposits alternating with four fossiliferous, paralic to shallow-marine units (Units 1–4). The transgressive surfaces identified at the base of these units mark major flooding events, forced by Holocene (Unit 4), Late Pleistocene (Unit 3) and Middle Pleistocene (Units 1, 2) interglacials. Distinct stratigraphic patterns typify the Middle Pleistocene interval, which includes coastal-marine (tidal inlet and bay) deposits. In contrast, lagoonal sediments record the maximum marine influence in the Late Pleistocene-Holocene succession. As a whole, the meiofauna tracks a regressive trend, with the deepest conditions recorded by the oldest Unit 1 (MIS 9/11 age?).

Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches

Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Magnesium in subaqueous speleothems as a potential palaeotemperature proxy

Few palaeoclimate archives beyond the polar regions preserve continuous and datable palaeotemperature proxy time series over multiple glacial-interglacial cycles. This hampers efforts to develop a more coherent picture of global patterns of past temperatures. Here we show that Mg concentrations in a subaqueous speleothem from an Italian cave track regional sea-surface temperatures over the last 350,000 years. The Mg shows higher values during warm climate intervals and converse patterns during cold climate stages. In contrast to previous studies, this implicates temperature, not rainfall, as the principal driver of Mg variability. The depositional setting of the speleothem gives rise to Mg partition coefficients that are more temperature dependent than other calcites, enabling the effect of temperature change on Mg partitioning to greatly exceed the effects of changes in source-water Mg/Ca. Subaqueous speleothems from similar deep-cave environments should be capable of providing palaeotemperature information over multiple glacial-interglacial cycles.

Late Pleistocene Landscape and Settlement Dynamics of Portuguese Estremadura

ABSTRACTHere we report the results of an integrated geoarchaeological survey to study Palaeolithic human settlement dynamics in the coastal region of Portuguese Estremadura. The region has been an important focus of human occupation across multiple glacial-interglacial cycles, including periods of well-documented abrupt climate instability during MIS 3 and 2. The pedestrian survey covered a roughly 10 km wide strip of land between São Pedro de Muel and Peniche. The survey intensively targeted three landscape settings with Pleistocene-age surfaces and sediments: the coastal bluffs with exposed aeolian, fluvial, and colluvial sands; the Caldas da Rainha diapiric valley and associated fluvial/estuarine fills; and, Cretaceous chert-rich limestone uplands that bound the inland margin of the study area. We discovered dozens of new Palaeolithic sites, analyzed numerous Pleistocene sedimentary sections, and applied widespread OSL-dating to establish age control that allowed us to build a regional geomorphic history to contextualize Late Pleistocene human settlement across the region.

Biodiversity gradients emerge.

A simulation of expansion, fragmentation and extirpation of species ranges over multiple glacial–interglacial cycles matches empirical biodiversity gradients and shows that high levels of biodiversity in the tropics can emerge from temporally variable but spatially patchy precipitation regimes, driven by allopatric speciation.

Quantitative reconstruction of early Holocene and last glacial climate on the Balkan Peninsula using coupled hydrological and isotope mass balance modelling

We investigate the modern hydrology of Lake Ohrid (Macedonia/Albania) using a combined hydrological and isotope-based modelling approach and present a new evaluation of contemporary water balance and palaeoclimate estimates. The combined model is able to estimate hydrological components that cannot be directly measured, and indicates that sublacustrine spring inflow is in the order of 50% higher than previous estimates and groundwater outflow comprises approximately a third of overall water outflow. In combination with sediment core oxygen isotope data, we used the combined model to quantitatively reconstruct past climate, in particular precipitation, during the early Holocene and last glacial period. Calculated precipitation in the early Holocene was higher than the value for present day and was approximately 44% lower than present during the last glacial, assuming the majority of precipitation fell as snow. The estimated amount of precipitation in the last glacial would have been high enough to provide refugial conditions at Lake Ohrid and to support the continuous existence of arboreal vegetation in the catchment. The improved understanding of the modern isotope hydrology of Lake Ohrid is fundamental for explaining the systematics of past isotope variation and providing context for extended sediment records from the lake, which will provide longer-term palaeoclimate reconstructions covering multiple glacial-interglacial cycles.

Late Quaternary multiple incised valley systems: An unusually well‐preserved stratigraphic record of two interglacial valley‐fill successions from the Arno Plain (northern Tuscany, Italy)

AbstractUn‐fragmented stratigraphic records of late Quaternary multiple incised valley systems are rarely preserved in the subsurface of alluvial‐delta plains due to older valley reoccupation. The identification of a well‐preserved incised valley fill succession beneath the southern interfluve of the Last Glacial Maximum Arno palaeovalley (northern Italy) represents an exceptional opportunity to examine in detail evolutionary trends of a Mediterranean system over multiple glacial–interglacial cycles. Through sedimentological and quantitative meiofauna (benthic foraminifera and ostracods) analyses of two reference cores (80 m and 100 m long) and stratigraphic correlations, a mid‐Pleistocene palaeovalley, 5 km wide and 50 m deep, was reconstructed. Whereas valley filling is chronologically constrained to the penultimate interglacial (Marine Isotope Stage 7) by four electron spin resonance ages on bivalve shells (Cerastoderma glaucum), its incision is tentatively correlated with the Marine Isotope Stage 8 sea‐level fall. Above basal fluvial‐channel gravels, the incised valley fill is formed by a mud‐prone succession, up to 44 m thick, formed by a lower floodplain unit and an upper unit with brackish meiofauna that reflects the development of a wave‐dominated estuary. Subtle meiofauna changes towards less confined conditions record two marine flooding episodes, chronologically linked to the internal Marine Isotope Stage 7 climate‐eustatic variability. After the maximum transgressive phase, recorded by coastal sands, the interfluves were flooded around 200 ka (latest Marine Isotope Stage 7). The subsequent shift in river incision patterns, possibly driven by neotectonic activity, prevented valley reoccupation guiding the northward formation of the Last Glacial Maximum palaeovalley. The applied multivariate approach allowed the sedimentological characterization of the Marine Isotope Stage 7 and Marine Isotope Stage 1 palaeovalley fills, including shape, size and facies architecture, which revealed a consistent river‐coastal system response over two non‐consecutive glacial–interglacial cycles (Marine Isotope Stages 8 to 7 and Marine Isotope Stages 2 to 1). The recurring stacking pattern of facies documents a predominant control exerted on stratigraphy by Milankovitch and sub‐Milankovitch glacio‐eustatic oscillations across the late Quaternary period.

Phylogeography of postglacial range expansion in Juglans mandshurica (Juglandaceae) reveals no evidence of bottleneck, loss of genetic diversity, or isolation by distance in the leading-edge populations

The past studies of postglacial recolonization patterns in high latitude regions have revealed a significant role of dispersal capacity in shaping the genetic diversity and population structure of temperate trees. However, most of these studies have focused on species with long-distance dispersal followed by exponential population growth and were therefore unable to reveal the patterns in the case of a gradual expansion. Here we studied the impacts of postglacial range expansions on the distribution of genetic diversity in the Manchurian walnut (Juglans mandshurica), a common tree of East Asian cool-temperate deciduous forests that apparently lacks long-distance seed dispersal ability. The genetic diversity and structure of 19 natural walnut populations in Northeast China and the Korean Peninsula were examined using 17 nuclear simple sequence repeat (SSR) loci. Potential habitats under current and past climatic conditions were predicted using the ecological niche modelling (ENM) method. Bayesian clustering analysis revealed three groups, which were inferred to have diverged through multiple glacial-interglacial cycles in multiple refugia during the Quaternary Period. ENM estimated a southward range shift at the LGM, but high suitability scores still occurred in the western parts of the Changbai Mountains (Northeast China), the Korean peninsula and the exposed seafloor of the Yellow Sea. In contrast to most other cool-temperate trees co-occurring in the same region, the Manchurian walnut did not show any evidence of a population bottleneck, loss of genetic diversity or isolation by distance during the postglacial expansion. Our study clearly indicates that current northern populations originated from one glacial lineage and recolonization via a gradually advancing front due to the lack of a long-distance seed dispersal mechanism led to no latitudinal decrease in genetic diversity.

Multiple glacial refugia for cool-temperate deciduous trees in northern East Asia: the Mongolian oak as a case study.

In East Asia, temperate forests are predicted to have retracted southward to c. 30° N during the last glacial maximum (LGM) based on fossil pollen data, whereas phylogeographic studies have often suggested glacial in situ survival of cool-temperate deciduous trees in their modern northern ranges. Here we report a study of the genetic diversity and structure of 29 natural Mongolian oak (Quercus mongolica) populations using 19 nuclear simple sequence repeat (nSSR) loci and four chloroplast DNA fragments. Bayesian clustering analysis with nSSRs revealed five groups, which were inferred by approximate Bayesian computation (ABC) to have diverged in multiple refugia through multiple glacial-interglacial cycles. Analysis of chloroplast DNA variation revealed four lineages that were largely but incompletely geographically disjunct. Ecological niche modelling (ENMs) indicated a southward range shift of the oak's distribution at the LGM, although high suitability scores were also evident in the Changbai Mts. (Northeast China), the Korean Peninsula, areas surrounding the Bohai Sea, and along the coast of the Russian Far East. In addition, endemic chloroplast DNA haplotypes and nuclear lineages occurred in high-latitude northern areas where the ENM predicted no suitable habitat. The combined evidence from nuclear and chloroplast DNA, and the results of the ENM clearly demonstrate that multiple northern refugia, including cryptic ones, were maintained across the current distributional range of the Mongolian oak during the LGM or earlier glacial periods. Though spatially limited, postglacial expansions from these refugia have led to a pattern of decreased genetic diversity with increasing latitude.

Population genetic structure of Oryza rufipogon and Oryza nivara: implications for the origin of O.nivara.

Ecological speciation plays a primary role in driving species divergence and adaptation. Oryza rufipogon and Oryza nivara are two incipient species at the early stage of speciation with distinct differences in morphology, life history traits and habitat preference, and therefore provide a unique model for the study of ecological speciation. However, the population genetic structure of the ancestral O.rufipogon has been controversial despite substantial study, and the origin of the derivative O.nivara remains unclear. Here, based on sequences of 10 nuclear and two chloroplast loci from 26 wild populations across the entire geographic ranges of the two species, we conducted comprehensive analyses using population genetics, phylogeography and species distribution modelling (SDM) approaches. In addition to supporting the two previously reported major subdivisions, we detected four genetically distinct groups within O.rufipogon and found no correlation between the genetic groups and either species identity or geographical regions. The SDM clearly showed substantial change in the distribution range of O.rufipogon in history, demonstrating that the repeated extinction and colonization of local populations due to multiple glacial-interglacial cycles during the Quaternary was most likely the main factor shaping the confounding population genetic structure of O.rufipogon. Moreover, we found significant differences between the two species in climate preferences, suggestive of an important role for climatic factors in the adaptation, persistence and expansion of O.nivara. Finally, based on the genetic pattern and dynamics of the O. nivara populations, we hypothesize that O.nivara might have independently originated multiple times from different O.rufipogon populations.

Coupled micromorphological and stable isotope analysis of Quaternary calcrete development

Pedogenic calcretes are widespread in arid and semi-arid regions. Using calcrete profiles from four river terraces of the Rio Alias in southeast Spain, this study explores the potential of using detailed micromorphological and stable isotopic analysis to more fully understand the impacts of Quaternary environmental change on calcrete development. The four profiles increase in carbonate complexity with progressive age, reflecting calcretisation over multiple glacial–interglacial cycles since MIS 9 (c. 300ka). Calcrete profiles contain a mixture of Alpha (non-biogenic) and Beta (biogenic) microfabrics. Alpha fabrics have higher δ13C and δ18O values. The profiles contain a range of crystal textures, but there is little difference between the δ13C and δ18O values of spar, microspar, and micrite cements. Strong positive covariance between δ13C and δ18O suggests that both isotopes are responding to the same environmental parameter, which is inferred to be relative aridity. The study reveals that the detailed co-analysis of calcrete micromorphology and stable isotope signatures can allow patterns of calcrete formation to be placed into a wider palaeoclimatic context. This demonstrates the potential of this technique to more reliably constrain the palaeoenvironmental significance of secondary carbonates in dryland settings where other proxy records may be poorly preserved.