Assessment of Intraspecific Variability in the Forest Dormouse (Dryomys nitedula) and Woolly Dormouse (Dryomys laniger) from Türkiye and Adjacent Regions Based on Mitochondrial DNA

This study aimed to reveal intraspecific variations in two Dryomys species distributed in Türkiye, based on mitochondrial DNA cytochrome b gene sequences, and to discuss the factors driving these variations in the context of phylogeography and genetic species concepts. As a result of Maximum Likelihood, Bayesian Inference, and Network analyses, which included haplogroups or lineages from Italy, Russia, the Caucasus, and Iran identified in previous studies, along with Turkish haplotypes, three major clades (MC1, MC2, and MC3) were identified within Dryomys nitedula. These clades began to diverge evolutionarily in the middle of the Late Miocene (8.82 million years ago) and exhibit significant genetic differences from one another. The Turkish haplotypes were divided into five distinct lineages (N1-N5), each within five subclades (SC1-SC5), which were nested within these MCs. These lineages, their geographical distributions, and the subspecies defined in previous studies that correspond to these lineages are as follows: N1 from the Thrace region (Dryomys nitedula wingei), N2 from the Black Sea region (potentially a new subspecies), N3 from western and central Anatolia (Dryomys nitedula phrygius), N4 from northeastern Anatolia (Dryomys nitedula tichomirowi), and N5 from eastern Anatolia (Dryomys nitedula pictus). The N2 lineage, distributed in areas close to the coastal side of the Eastern Black Sea region and with a range close to both N3 (D. n. phrygius) and N4 (D. n. tichomirowi), exhibited high genetic differentiation from these two lineages and was a candidate to be treated as a new subspecies of Dryomys nitedula in Türkiye. The N5 lineage, which includes haplotypes from the distribution areas of the populations initially classified as Dryomys pictus and later as Dryomys nitedula pictus in previous studies, was found to be more closely related to Dryomys nitedula kurdistanicus from the Zagros Mountains than to D. n. pictus from the central regions of Iran. Combining the results of this study with previous research, it is clear that the D. nitedula lineages in Türkiye, along with haplogroups or subspecies in neighboring regions diverged between the middle Late Miocene and Middle Pleistocene. This divergence is believed to have been driven by climatic cycles and geomorphological processes that shaped the topography of their distribution range. The high genetic diversity observed in the lineages of Anatolia suggests that the region may have served as a glacial refuge for D. nitedula. Similarly to the processes and factors shaping the evolution of D. nitedula, Dryomys laniger was found to have diverged into two lineages, western (L1) and eastern (L2 or Dryomys anatolicus), within its distribution range during the Late Pliocene (2.94 Mya). To make a more accurate taxonomic assessment of D. laniger, a larger number of samples is needed, and the distribution limits should be more clearly defined.

Reduced North Pacific Deep Water formation across the Northern Hemisphere Glaciation

Today, deep waters do not form in the northern high latitudes of the Pacific Ocean, but this may not have been the case during the Pliocene. Evidence suggests there was a Pacific Meridional Overturning Circulation during the warm Late Pliocene, similar to the modern Atlantic Ocean with a weak halocline in the subpolar North Pacific resulting in North Pacific Deep Water (NPDW) formation. However, much of this evidence comes from proxies that can be biased by changes in local productivity. We present a coupled Mg/Ca-δ18O record from the North Pacific which shows two distinct water masses in the Pliocene North Pacific Ocean, with NPDW colder and fresher than the underlying deeper water. Here, we show a decline in NPDW formation during glacials from 2.7 million years ago, which we suggest demonstrates the strong sensitivity of ocean gateways to sea level and ice volume change in shaping deep water circulation, and thus the climate system.

Climate and Human Evolution: Insights from Marine Records

The relationship between climate and human evolution is complex, and the causal mechanisms remain unknown. Here, we review and synthesize what is currently known about climate forcings on African landscapes, focusing mainly on the last 4 million years. We use information derived from marine sediment archives and data-numerical climate model comparisons and integration. There exists a heterogeneity in pan-African hydroclimate changes, forced by a combination of orbitally paced, low-latitude fluctuations in insolation; polar ice volume changes; tropical sea surface temperature gradients linked to the Walker circulation; and possibly greenhouse gases. Pan-African vegetation changes do not follow the same pattern, which is suggestive of additional influences, such as CO2 and temperature. We caution against reliance on temporal correlations between global or regional climate, environmental changes, and human evolution and briefly proffer some ideas on how pan-African climate trends could help create novel conceptual frameworks to determine the causal mechanisms of associations between climate/habitat change and hominin evolution.

A new discovery of cryptic species in barred sand burrower Limnichthys fasciatus (Pisces: Creediidae) from the western Pacific, with evolutionary perspectives of anti-equatorial species

Although the barred sand burrower, Limnichthys fasciatus, is widely distributed throughout the western Pacific, including Japan, Taiwan, and Australia, its morphology and genetics are poorly known. We discovered four cryptic species of Limnichthys from the western Pacific based on mtDNA cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (rRNA) sequences. Genetic distances showed remarkably large differences (12.7%-24.3% in COI and 7.9%-19.6% in 16S rRNA) between true L. fasciatus (type locality: southeastern Australia) and the others. A relaxed clock model with optimized selected substitution models showed that their deep divergence began in the middle Miocene epoch and subsequently diverged into the current cryptic species in the Plio-Pleistocene. A eurythermal common ancestor may have evolved independently in each region due to geographical events and paleoclimatic fluctuations, which made it possible for L. fasciatus complex to be an anti-equatorial species. Despite their deep genetic divergence, they showed marked phenotypic similarity, suggesting that they have experienced similar selective pressures related to their specific behavior.

Revisiting the mid-Pleistocene transition ocean circulation crisis

The mid-Pleistocene transition (MPT) [~1.25 to 0.85 million years ago (Ma)] marks a shift in the character of glacial-interglacial climate (1, 2). One prevailing hypothesis for the origin of the MPT is that glacial deep ocean circulation fundamentally changed, marked by a circulation "crisis" at ~0.90 Ma (marine isotope stages 24 to 22) (3). Using high-resolution paired neodymium, carbon, and oxygen isotope data from the South Atlantic Ocean (Cape Basin) across the MPT, we find no evidence of a substantial change in deep ocean circulation. Before and during the early MPT (~1.30 to 1.12 Ma), the glacial deep ocean variability closely resembled that of the most recent glacial cycle. The carbon storage facilitated by developing deep ocean stratification across the MPT required only modest circulation adjustments.

Recent and Recurrent Autopolyploidization Fueled Diversification of Snow Carp on the Tibetan Plateau

Whole-genome duplication (WGD), or polyploidization, is a major contributor to biodiversity. However, the establishment and survival of WGDs are often considered to be stochastic, since elucidating the processes of WGD establishment remains challenging. In the current study, we explored the processes leading to polyploidy establishment in snow carp (Cyprinidae: Schizothoracinae), a predominant component of the ichthyofauna of the Tibetan Plateau and its surrounding areas. Using large-scale genomic data from isoform sequencing, we analyzed ohnolog genealogies and divergence in hundreds to thousands of gene families across major snow carp lineages. Our findings demonstrated that independent autopolyploidization subsequent to speciation was prevalent, while autopolyploidization followed by speciation also occurred in the diversification of snow carp. This was further supported by matrilineal divergence and drainage evolution evidence. Contrary to the long-standing hypothesis that ancient polyploidization preceded the diversification of snow carp, we determined that polyploidy in extant snow carp was established by recurrent autopolyploidization events during the Pleistocene. These findings indicate that the diversification of extant snow carp resembles a coordinated duet: first, the uplift of the Tibetan Plateau orchestrated the biogeography and diversification of their diploid progenitors; then, the extensive Pliocene-Pleistocene climate changes acted as relay runners, further fueling diversification through recurrent autopolyploidization. Overall, this study not only reveals a hitherto unrecognized recent WGD lineage in vertebrates but also advances current understanding of WGD processes, emphasizing that WGD establishment is a nonstochastic event, emerging from numerous adaptations to environmental challenges and recurring throughout evolutionary history rather than merely in plants.

Geologically constrained 2-million-year-long simulations of Antarctic Ice Sheet retreat and expansion through the Pliocene

Pliocene global temperatures periodically exceeded modern levels, offering insights into ice sheet sensitivity to warm climates. Ice-proximal geologic records from this period provide crucial but limited glimpses of Antarctic Ice Sheet behavior. We use an ice sheet model driven by climate model snapshots to simulate transient glacial cyclicity from 4.5 to 2.6 Ma, providing spatial and temporal context for geologic records. By evaluating model simulations against a comprehensive synthesis of geologic data, we translate the intermittent geologic record into a continuous reconstruction of Antarctic sea level contributions, revealing a dynamic ice sheet that contributed up to 25 m of glacial-interglacial sea level change. Model grounding line behavior across all major Antarctic catchments exhibits an extended period of receded ice during the mid-Pliocene, coincident with proximal geologic data around Antarctica but earlier than peak warmth in the Northern Hemisphere. Marine ice sheet collapse is triggered with 1.5 °C model subsurface ocean warming.

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Ice core evidence for atmospheric oxygen decline since the Mid-Pleistocene transition

The history of atmospheric oxygen (PO₂) and the processes that act to regulate it remain enigmatic because of difficulties in quantitative reconstructions using indirect proxies. Here, we extend the ice-core record of PO₂ using 1.5-million-year-old (Ma) discontinuous ice samples drilled from Allan Hills Blue Ice Area, East Antarctica. No statistically significant difference exists in PO₂ between samples at 1.5 Ma and 810 thousand years (ka), suggesting that the Late-Pleistocene imbalance in O₂ sources and sinks began around the time of the transition from 40- to 100-ka glacial cycles in the Mid-Pleistocene between ~1.2 Ma and 700 ka. The absence of a coeval secular increase in atmospheric CO₂ over the past ~1 Ma requires negative feedback mechanisms such as Pco₂-dependent silicate weathering. Fast processes must also act to suppress the immediate Pco₂ increase because of the imbalance in O₂ sinks over sources beginning in the Mid-Pleistocene.

Evidence for a Northern Hemispheric trigger of the 100,000-y glacial cyclicity

The causes of the Mid-Pleistocene Transition, the shift from ∼41-ky to 100-ky interglacial-glacial cycles and more intense ice ages, remain intensely debated, as this fundamental change occurred between ∼1,250 and 650 ka without substantial changes in astronomical climate forcings. Recent studies disagree about the relative importance of events and processes in the Northern and Southern Hemispheres, as well as whether the shift occurred gradually over several interglacial-glacial cycles or abruptly at ∼900 ka. We address these issues using a north-to-south reconstruction of the Atlantic arm of the global meridional overturning ocean circulation, a primary means for distributing heat around the globe, using neodymium (Nd) isotopes. Results reveal a period of intense erosion affecting the cratonic shields surrounding the North Atlantic between Marine Isotope Stages (MIS) 27 and 25 (∼980 and 950 ka), reflected by unusually low Nd isotope ratios in deep North Atlantic seawater. This episode preceded a major ocean circulation weakening between MIS 25 and 21 (950 and 860 ka) that coincided with the first ∼100-ky-long interglacial-glacial onset of Northern Hemisphere glaciation at around 2.4 to 2.8 Ma. The data point to a Northern Hemisphere-sourced initiation for the transition, possibly induced through regolith loss and increased exposure of the crystalline bedrock, which would lead to increased friction, enabling larger ice sheets that are characteristic of the 100-ky interglacial-glacial cycles.

Historical biogeography of the genus Rhadinaea (Squamata: Dipsadinae)

Multiple geological and climatic events have created geographical or ecological barriers associated with speciation events, playing a role in biological diversification in North and Central America. Here, we evaluate the influence of the Neogene and Quaternary geological events, as well as the climatic changes in the diversification of the colubrid snake genus Rhadinaea using molecular dating and ancestral area reconstruction. A multilocus sequence dataset was generated for 37 individuals of Rhadinaea from most of the biogeographical provinces where the genus is distributed, representing 19 of the 21 currently recognized species, and two undescribed species. Our analyses show that the majority of the Rhadinaea species nest in two main clades, herein identified as "Eastern" and "Southern". These clades probably diverged from each other in the early Miocene, and their divergence was followed by 11 divergences during the middle to late Miocene, three divergences during the Pliocene, and six divergences in the Pleistocene. The ancestral distribution of Rhadinaea was reconstructed across the Sierra Madre del Sur. Our phylogenetic analyses do not support the monophyly of Rhadinaea. The Miocene and Pliocene geomorphology, perhaps in conjunction with climate change, appears to have triggered the diversification of the genus, while the climatic changes during the Miocene probably induced the diversification of Rhadinaea in the Sierra Madre del Sur. Our analysis suggests that the uplifting of the Trans-Mexican Volcanic Belt and Chiapan-Guatemalan highlands in this same period resulted in northward and southward colonization events. This was followed by more recent, independent colonization events in the Pliocene and Pleistocene involving the Balsas Basin, Chihuahuan Desert, Pacific Coast, Sierra Madre Occidental, Sierra Madre Oriental, Sierra Madre del Sur, Trans-Mexican Volcanic Belt, and Veracruz provinces, probably driven by the climatic fluctuations of the time.

Causes and timing of recurring subarctic Pacific hypoxia

Several North Pacific studies of the last deglaciation show hypoxia throughout the ocean margins and attribute this phenomenon to the effects of abrupt warming and meltwater inputs. Yet, because of the lack of long records spanning multiple glacial cycles and deglaciation events, it is unclear whether deoxygenation was a regular occurrence of warming events and whether deglaciation and/or other conditions promoted hypoxia throughout time. Here, subarctic Pacific laminated sediments from the past 1.2 million years demonstrate that hypoxic events recurred throughout the Pleistocene as episodes of highly productive phytoplankton growth and were generally associated with interglacial climates, high sea levels, and enhanced nitrate utilization-but not with deglaciations. We suggest that hypoxia was typically stimulated by high productivity from iron fertilization facilitated by redox-remobilized iron from flooded continental shelves.

Paleoceanography of the northwestern Pacific across the Early-Middle Pleistocene boundary (Marine Isotope Stages 20-18)

The fluctuating position of the boundary between the Kuroshio (warm) and Oyashio (cold) currents in the mid-latitude western North Pacific affects both heat transport and air-ocean interactions and has significant consequences for the East Asian climate. We reconstruct the paleoceanography of Marine Isotope Stages (MIS) 20-18, MIS 19 being one of the closest astronomical analogues to the present interglacial, through multiple proxies including microfossil assemblage data, planktonic foraminiferal isotopes (δ18O and δ13C), and foraminiferal Mg/Ca-based temperature records, from the Chiba composite section (CbCS) exposed on the Boso Peninsula, east-central Japan. Principal component analysis (PCA) is used to capture dominant patterns of the temporal variation in these marine records, and shows that the relative abundances of calcareous nannofossil and radiolarian taxa are consistent with the water mass types inferred from geochemical proxies. The leading mode (36.3% of total variance) mirrors variation in the terrestrial East Asian winter monsoon (EAWM), reflecting seasonal trends dominated by the winter monsoon system. In the CbCS, this mode is interpreted as reflecting the interplay between the warm Kuroshio and cold Oyashio waters, which is likely related to the latitudinal shift of the subtropical-subarctic gyre boundary in the North Pacific. The second mode (15.4% of total variance) is closely related to subsurface conditions. The leading mode indicates that MIS 19b and 19a are represented by millennial-scale stadial/interstadial oscillations. Northerly positions for the gyre boundary during late MIS 19c, the interstadials of MIS 19a, and early MIS 18 are inferred from the leading mode, which is consistent with a weak EAWM and consequent mild winter climate in East Asia. Nonetheless, the northerly positions for the gyre boundary during late MIS 19c and early MIS 19a were not associated with subsurface warming presumably due to the suppressed gyre circulation itself caused by the weak Aleutian Low. Intermittent southerly positions for the gyre boundary are inferred for the stadials of MIS 19b and 19a. Regional sea surface temperature (SST) comparisons in the western North Pacific reveal that the moderate SSTs during MIS 19a through early MIS 18 were restricted to the mid- to high latitudes, influenced by the weak EAWM. Comparison between MIS 20-18 and MIS 2-1 suggests that glacial MIS 20 and 18 had significantly milder winters than MIS 2, likely related to the relatively weak EAWM.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40645-020-00395-3.

Sea-level stands from the Western Mediterranean over the past 6.5 million years

Sea-level reconstructions are important for understanding past ice sheet variability and its response to past and future warming. Here we present Neogene and Quaternary sea-level snapshots using phreatic overgrowths on speleothems (POS) from caves on Mallorca, Spain. POS are excellent sea level index points because of their clear relationship to sea level and precise U–Pb chronology. We find that local sea-level before and at the onset of the Messinian Salinity Crisis was at 33.3 ± 0.25 m (6.54 ± 0.37 Ma) and 31.8 ± 0.25 m (5.86 ± 0.60 Ma) above present levels, respectively. We further present global mean sea level (GMSL) estimates, i.e. local sea level corrected for glacial isostatic adjustment and long-term uplift, for three other POS. The results show that GMSL during the Pliocene–Pleistocene Transition was 6.4 m (− 2.0–8.8 m) at 2.63 ± 0.11 Ma and during the beginning and the end of the Mid-Pleistocene Transition was − 1.1 m (− 5.6–2.4 m) and 5 m (1.5–8.1 m), respectively. These estimates provide important constraints for the past evolution of sea level and show that local sea level prior to the MSC was similar to the highest stand during the Pliocene, with markedly lower position afterwards.

Combining Modern and Paleoceanographic Perspectives on Ocean Heat Uptake

Monitoring Earth's energy imbalance requires monitoring changes in the heat content of the ocean. Recent observational estimates indicate that ocean heat uptake is accelerating in the twenty-first century. Examination of estimates of ocean heat uptake over the industrial era, the Common Era of the last 2,000 years, and the period since the Last Glacial Maximum, 20,000 years ago, permits a wide perspective on modern-day warming rates. In addition, this longer-term focus illustrates how the dynamics of the deep ocean and the cryosphere were active in the past and are still active today. The large climatic shifts that started with the melting of the great ice sheets have involved significant ocean heat uptake that was sustained over centuries and millennia, and modern-ocean heat content changes are small by comparison.

A new sea-level record for the Neogene/Quaternary boundary reveals transition to a more stable East Antarctic Ice Sheet

Sea-level rise resulting from the instability of polar continental ice sheets represents a major socioeconomic hazard arising from anthropogenic warming, but the response of the largest component of Earth's cryosphere, the East Antarctic Ice Sheet (EAIS), to global warming is poorly understood. Here we present a detailed record of North Atlantic deep-ocean temperature, global sea-level, and ice-volume change for ∼2.75 to 2.4 Ma ago, when atmospheric partial pressure of carbon dioxide (pCO₂) ranged from present-day (>400 parts per million volume, ppmv) to preindustrial (<280 ppmv) values. Our data reveal clear glacial-interglacial cycles in global ice volume and sea level largely driven by the growth and decay of ice sheets in the Northern Hemisphere. Yet, sea-level values during Marine Isotope Stage (MIS) 101 (∼2.55 Ma) also signal substantial melting of the EAIS, and peak sea levels during MIS G7 (∼2.75 Ma) and, perhaps, MIS G1 (∼2.63 Ma) are also suggestive of EAIS instability. During the succeeding glacial-interglacial cycles (MIS 100 to 95), sea levels were distinctly lower than before, strongly suggesting a link between greater stability of the EAIS and increased land-ice volumes in the Northern Hemisphere. We propose that lower sea levels driven by ice-sheet growth in the Northern Hemisphere decreased EAIS susceptibility to ocean melting. Our findings have implications for future EAIS vulnerability to a rapidly warming world.

Interglacials of the Quaternary defined by northern hemispheric land ice distribution outside of Greenland

Glacial/interglacial dynamics during the Quaternary were suggested to be mainly driven by obliquity (41-kyr periodicity), including irregularities during the last 1 Myr that resulted in on average 100-kyr cycles. Here, we investigate this so-called Mid-Pleistocene Transition via model-based deconvolution of benthic δ¹⁸O, redefining interglacials by lack of substantial northern hemispheric land ice outside of Greenland. We find that in 67%, 88% and 52% of the obliquity cycles during the early, middle and late Quaternary, respectively, a glacial termination is realized leading to irregular appearances of new interglacials during various parts of the last 2.6 Myr. This finding suggests that the proposed idea of terminations leading to new interglacials in the Quaternary as obliquity driven with growing influence of land ice volume on the timing of deglaciations during the last 1 Myr might be too simple. Alternatively, the land ice-based definition of interglacials needs revision if applied to the entire Quaternary.

This study presents a new definition of interglacials during the Quaternary. The authors find the appearance of interglacials is in general following the 41-kyr cycle of obliquity with various exceptions, suggesting a more complex physical mechanism triggering glacial terminations.

Predominant east to west colonizations across major oceanic barriers: Insights into the phylogeographic history of the hydroid superfamily Plumularioidea, suggested by a mitochondrial DNA barcoding marker

We provide preliminary insights into the global phylogeographic and evolutionary patterns across species of the hydrozoan superfamily Plumularioidea (Cnidaria: Hydrozoa). We analyzed 1,114 16S sequences of 198 putative species of Plumularioidea collected worldwide. We investigated genetic connections and divergence in relation to present-day and ancient biogeographic barriers, climate changes and oceanic circulation. Geographical distributions of most species are generally more constrained than previously assumed. Some species able to raft are dispersed widely. Human-mediated dispersal explains some wide geographical ranges. Trans-Atlantic genetic connections are presently unlikely for most of the tropical-temperate species, but were probably more frequent until the Miocene-Pliocene transition, before restriction of the Tethys Sea and the Central American Seaway. Trans-Atlantic colonizations were predominantly directed westwards through (sub)tropical waters. The Azores were colonized multiple times and through different routes, mainly from the east Atlantic, at least since the Pliocene. Extant geminate clades separated by the Isthmus of Panama have predominantly Atlantic origin. Various ancient colonizations mainly directed from the Indian Ocean to the Atlantic occurred through the Tethys Sea and around South Africa in periods of lower intensity of the Benguela upwelling. Thermal tolerance, population sizes, dispersal strategies, oceanic currents, substrate preference, and land barriers are important factors for dispersal and speciation of marine hydroids.

Extensive marine-terminating ice sheets in Europe from 2.5 million years ago

Geometries of Early Pleistocene [2.58 to 0.78 million years (Ma) ago] ice sheets in northwest Europe are poorly constrained but are required to improve our understanding of past ocean-atmosphere-cryosphere coupling. Ice sheets are believed to have changed in their response to orbital forcing, becoming, from about 1.2 Ma ago, volumetrically larger and longer-lived. We present a multiproxy data set for the North Sea, extending to over a kilometer below the present-day seafloor, which demonstrates spatially extensive glaciation of the basin from the earliest Pleistocene. Ice sheets repeatedly entered the North Sea, south of 60°N, in water depths of up to ~250 m from 2.53 Ma ago and subsequently grounded in the center of the basin, in deeper water, from 1.87 Ma ago. Despite lower global ice volumes, these ice sheets were near comparable in spatial extent to those of the Middle and Late Pleistocene but possibly thinner and moving over slippery (low basal resistance) beds.

Glacial Steady State Topography Controlled by the Coupled Influence of Tectonics and Climate

Glaciers and rivers are the main agents of mountain erosion. While in the fluvial realm empirical relationships and their mathematical description, such as the stream power law, improved the understanding of fundamental controls on landscape evolution, simple constraints on glacial topography and governing scaling relations are widely lacking. We present a steady state solution for longitudinal profiles along eroding glaciers in a coupled system that includes tectonics and climate. We combined the shallow ice approximation and a glacial erosion rule to calculate ice surface and bed topography from prescribed glacier mass balance gradient and rock uplift rate. Our approach is inspired by the classic application of the stream power law for describing a fluvial steady state but with the striking difference that, in the glacial realm, glacier mass balance is added as an altitude-dependent variable. From our analyses we find that ice surface slope and glacial relief scale with uplift rate with scaling exponents indicating that glacial relief is less sensitive to uplift rate than relief in most fluvial landscapes. Basic scaling relations controlled by either basal sliding or internal deformation follow a power law with the exponent depending on the exponents for the glacial erosion rule and Glen's flow law. In a mixed scenario of sliding and deformation, complicated scaling relations with variable exponents emerge. Furthermore, a cutoff in glacier mass balance or cold ice in high elevations can lead to substantially larger scaling exponents which may provide an explanation for high relief in high latitudes.

100- kyr cyclicity in volcanic ash emplacement: evidence from a 1.1 Myr tephra record from the NW Pacific

It is a longstanding observation that the frequency of volcanism periodically changes at times of global climate change. The existence of causal links between volcanism and Earth’s climate remains highly controversial, partly because most related studies only cover one glacial cycle. Longer records are available from marine sediment profiles in which the distribution of tephras records frequency changes of explosive arc volcanism with high resolution and time precision. Here we show that tephras of IODP Hole U1437B (northwest Pacific) record a cyclicity of explosive volcanism within the last 1.1 Myr. A spectral analysis of the dataset yields a statistically significant spectral peak at the ~100 kyr period, which dominates the global climate cycles since the Middle Pleistocene. A time-domain analysis of the entire eruption and δ¹⁸O record of benthic foraminifera as climate/sea level proxy shows that volcanism peaks after the glacial maximum and ∼13 ± 2 kyr before the δ¹⁸O minimum right at the glacial/interglacial transition. The correlation is especially good for the last 0.7 Myr. For the period 0.7–1.1 Ma, during the Middle Pleistocene Transition (MPT), the correlation is weaker, since the 100 kyr periodicity in the δ¹⁸O record diminishes, while the tephra record maintains its strong 100 kyr periodicity.

Sea ice dynamics across the Mid-Pleistocene transition in the Bering Sea

Sea ice and associated feedback mechanisms play an important role for both long- and short-term climate change. Our ability to predict future sea ice extent, however, hinges on a greater understanding of past sea ice dynamics. Here we investigate sea ice changes in the eastern Bering Sea prior to, across, and after the Mid-Pleistocene transition (MPT). The sea ice record, based on the Arctic sea ice biomarker IP₂₅ and related open water proxies from the International Ocean Discovery Program Site U1343, shows a substantial increase in sea ice extent across the MPT. The occurrence of late-glacial/deglacial sea ice maxima are consistent with sea ice/land ice hysteresis and land-glacier retreat via the temperature-precipitation feedback. We also identify interactions of sea ice with phytoplankton growth and ocean circulation patterns, which have important implications for glacial North Pacific Intermediate Water formation and potentially North Pacific abyssal carbon storage.

Active Pacific meridional overturning circulation (PMOC) during the warm Pliocene

An essential element of modern ocean circulation and climate is the Atlantic meridional overturning circulation (AMOC), which includes deep-water formation in the subarctic North Atlantic. However, a comparable overturning circulation is absent in the Pacific, the world's largest ocean, where relatively fresh surface waters inhibit North Pacific deep convection. We present complementary measurement and modeling evidence that the warm, ~400-ppmv (parts per million by volume) CO₂ world of the Pliocene supported subarctic North Pacific deep-water formation and a Pacific meridional overturning circulation (PMOC) cell. In Pliocene subarctic North Pacific sediments, we report orbitally paced maxima in calcium carbonate accumulation rate, with accompanying pigment and total organic carbon measurements supporting deep-ocean ventilation-driven preservation as their cause. Together with high accumulation rates of biogenic opal, these findings require vigorous bidirectional communication between surface waters and interior waters down to ~3 km in the western subarctic North Pacific, implying deep convection. Redox-sensitive trace metal data provide further evidence of higher Pliocene deep-ocean ventilation before the 2.73-Ma (million years) transition. This observational analysis is supported by climate modeling results, demonstrating that atmospheric moisture transport changes, in response to the reduced meridional sea surface temperature gradients of the Pliocene, were capable of eroding the halocline, leading to deep-water formation in the western subarctic Pacific and a strong PMOC. This second Northern Hemisphere overturning cell has important implications for heat transport, the ocean/atmosphere cycle of carbon, and potentially the equilibrium response of the Pacific to global warming.

Impact of glacial/interglacial sea level change on the ocean nitrogen cycle

The continental shelves are the most biologically dynamic regions of the ocean, and they are extensive worldwide, especially in the western North Pacific. Their area has varied dramatically over the glacial/interglacial cycles of the last million years, but the effects of this variation on ocean biological and chemical processes remain poorly understood. Conversion of nitrate to N₂ by denitrification in sediments accounts for half or more of the removal of biologically available nitrogen ("fixed N") from the ocean. The emergence of continental shelves during ice ages and their flooding during interglacials have been hypothesized to drive changes in sedimentary denitrification. Denitrification leads to the occurrence of phosphorus-bearing, N-depleted surface waters, which encourages N₂ fixation, the dominant N input to the ocean. An 860,000-y record of foraminifera shell-bound N isotopes from the South China Sea indicates that N₂ fixation covaried with sea level. The N₂ fixation changes are best explained as a response to changes in regional excess phosphorus supply due to sea level-driven variations in shallow sediment denitrification associated with the cyclic drowning and emergence of the continental shelves. This hypothesis is consistent with a glacial ocean that hosted globally lower rates of fixed N input and loss and a longer residence time for oceanic fixed N-a "sluggish" ocean N budget during ice ages. In addition, this work provides a clear sign of sea level-driven glacial/interglacial oscillations in biogeochemical fluxes at and near the ocean margins, with implications for coastal organisms and ecosystems.

A persistent and dynamic East Greenland Ice Sheet over the past 7.5 million years

Climate models show that ice-sheet melt will dominate sea-level rise over the coming centuries, but our understanding of ice-sheet variations before the last interglacial 125,000 years ago remains fragmentary. This is because terrestrial deposits of ancient glacial and interglacial periods are overrun and eroded by more recent glacial advances, and are therefore usually rare, isolated and poorly dated. In contrast, material shed almost continuously from continents is preserved as marine sediment that can be analysed to infer the time-varying state of major ice sheets. Here we show that the East Greenland Ice Sheet existed over the past 7.5 million years, as indicated by beryllium and aluminium isotopes (¹⁰Be and ²⁶Al) in quartz sand removed by deep, ongoing glacial erosion on land and deposited offshore in the marine sedimentary record. During the early Pleistocene epoch, ice cover in East Greenland was dynamic; in contrast, East Greenland was mostly ice-covered during the mid-to-late Pleistocene. The isotope record we present is consistent with distinct signatures of changes in ice sheet behaviour coincident with major climate transitions. Although our data are continuous, they are from low-deposition-rate sites and sourced only from East Greenland. Consequently, the signal of extensive deglaciation during short, intense interglacials could be missed or blurred, and we cannot distinguish between a remnant ice sheet in the East Greenland highlands and a diminished continent-wide ice sheet. A clearer constraint on the behaviour of the ice sheet during past and, ultimately, future interglacial warmth could be produced by ¹⁰Be and ²⁶Al records from a coring site with a higher deposition rate. Nonetheless, our analysis challenges the possibility of complete and extended deglaciation over the past several million years.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.

Determining the distribution loss of brown eared-pheasant (Crossoptilon mantchuricum) using historical data and potential distribution estimates

We analyzed the synchronous relationship between forest cover and species distribution to explain the contraction in the distribution range of the brown eared-pheasant (Crossoptilon mantchuricum) in China. Historical resources can provide effective records for reconstructing long-term distribution dynamics. The brown eared-pheasant's historical distribution from 25 to 1947 CE, which included the three provinces of Shaanxi, Shanxi, and Hebei based on this species' habitat selection criteria, the history of the forests, ancient climate change records, and fossil data. The current species distribution covers Shaanxi, Shanxi, and Hebei provinces, as well as Beijing city, while Shanxi remains the center of the distribution area. MaxEnt model indicated that the suitable conditions of the brown eared-pheasant had retreated to the western regions of Shanxi and that the historical distribution area had reduced synchronously with the disappearance of local forest cover in Shanxi. We built a correlative relationship between the presence/absence of brown eared-pheasants and forest coverage and found that forest coverage in the north, northeast, central, and southeast areas of the Shanxi province were all less than 10% in 1911. Wild brown eared-pheasants are stable in the Luliang Mountains, where forest coverage reached 13.2% in 2000. Consequently, we concluded that the distribution of this species is primarily determined by vegetation conditions and that forest cover was the most significant determining factor.

Comparative phylogeography and evolutionary history of schizothoracine fishes in the Changtang Plateau and their implications for the lake level and Pleistocene climate fluctuations

The water level oscillation of endorheic lakes and extent change of glaciers associated with the Asian monsoon are known as prominent representatives of climatic and environmental events in the Tibetan Plateau during the Quaternary. However, details process in spatial and temporal changes are still debated. We use the schizothoracines as a palaeoclimatic proxy to test two hypotheses concerning the evolution of Quaternary glaciations and lakes of the Changtang Plateau: (1) the Tibetan glaciations generally tended to decrease since the middle Pleistocene; (2) the lakes expansion was driven by summer monsoon rainfall. Based on a wide range‐wide sampling throughout in the Changtang Plateau and its adjacent drainages, we constructed phylogeny and demographic histories of schizothoracines in the Changtang Plateau. Our results showed that the populations of the exorheic rivers and lakes in southern Tibet possessed higher genetic variability, earlier coalescent and expansion times than those of the endorheic lakes in the Changtang Plateau. Population expansions are highly consistent with phases of strong summer monsoon and high lake level during interglacial stages. The maximum growth rate intervals showed three pulses from 64.7 to 54.8, 39.6 to 31.0, and 14.9 to 2.4 kya respectively. The significant positive correlations were found between regional precipitation and genetic diversity, as well as coalescence time of populations in the endorheic lakes. We suggested that the demographic history of the schizothoracines reflects the spatial and temporal changes in climate and lake level, in particular, in regional precipitation gradients associated with changes of the South Asian monsoon, and supports the climatic hypothesis of a general diminishing tend in Tibetan glaciations in the Tibetan Plateau since the middle Pleistocene.

SEA-LEVEL RISE. Sea-level rise due to polar ice-sheet mass loss during past warm periods

Interdisciplinary studies of geologic archives have ushered in a new era of deciphering magnitudes, rates, and sources of sea-level rise from polar ice-sheet loss during past warm periods. Accounting for glacial isostatic processes helps to reconcile spatial variability in peak sea level during marine isotope stages 5e and 11, when the global mean reached 6 to 9 meters and 6 to 13 meters higher than present, respectively. Dynamic topography introduces large uncertainties on longer time scales, precluding robust sea-level estimates for intervals such as the Pliocene. Present climate is warming to a level associated with significant polar ice-sheet loss in the past. Here, we outline advances and challenges involved in constraining ice-sheet sensitivity to climate change with use of paleo-sea level records.

Exposure age and ice-sheet model constraints on Pliocene East Antarctic ice sheet dynamics

The Late Pliocene epoch is a potential analogue for future climate in a warming world. Here we reconstruct Plio-Pleistocene East Antarctic Ice Sheet (EAIS) variability using cosmogenic nuclide exposure ages and model simulations to better understand ice sheet behaviour under such warm conditions. New and previously published exposure ages indicate interior-thickening during the Pliocene. An ice sheet model with mid-Pliocene boundary conditions also results in interior thickening and suggests that both the Wilkes Subglacial and Aurora Basins largely melted, offsetting increased ice volume. Considering contributions from West Antarctica and Greenland, this is consistent with the most recent IPCC AR5 estimate, which indicates that the Pliocene sea level likely did not exceed +20 m on Milankovitch timescales. The inception of colder climate since ∼3 Myr has increased the sea ice cover and inhibited active moisture transport to Antarctica, resulting in reduced ice sheet thickness, at least in coastal areas.

Temperature impacts on deep-sea biodiversity

Temperature is considered to be a fundamental factor controlling biodiversity in marine ecosystems, but precisely what role temperature plays in modulating diversity is still not clear. The deep ocean, lacking light and in situ photosynthetic primary production, is an ideal model system to test the effects of temperature changes on biodiversity. Here we synthesize current knowledge on temperature-diversity relationships in the deep sea. Our results from both present and past deep-sea assemblages suggest that, when a wide range of deep-sea bottom-water temperatures is considered, a unimodal relationship exists between temperature and diversity (that may be right skewed). It is possible that temperature is important only when at relatively high and low levels but does not play a major role in the intermediate temperature range. Possible mechanisms explaining the temperature-biodiversity relationship include the physiological-tolerance hypothesis, the metabolic hypothesis, island biogeography theory, or some combination of these. The possible unimodal relationship discussed here may allow us to identify tipping points at which on-going global change and deep-water warming may increase or decrease deep-sea biodiversity. Predicted changes in deep-sea temperatures due to human-induced climate change may have more adverse consequences than expected considering the sensitivity of deep-sea ecosystems to temperature changes.

Sea-level and deep-sea-temperature variability over the past 5.3 million years

Ice volume (and hence sea level) and deep-sea temperature are key measures of global climate change. Sea level has been documented using several independent methods over the past 0.5 million years (Myr). Older periods, however, lack such independent validation; all existing records are related to deep-sea oxygen isotope (δ(18)O) data that are influenced by processes unrelated to sea level. For deep-sea temperature, only one continuous high-resolution (Mg/Ca-based) record exists, with related sea-level estimates, spanning the past 1.5 Myr. Here we present a novel sea-level reconstruction, with associated estimates of deep-sea temperature, which independently validates the previous 0-1.5 Myr reconstruction and extends it back to 5.3 Myr ago. We find that deep-sea temperature and sea level generally decreased through time, but distinctly out of synchrony, which is remarkable given the importance of ice-albedo feedbacks on the radiative forcing of climate. In particular, we observe a large temporal offset during the onset of Plio-Pleistocene ice ages, between a marked cooling step at 2.73 Myr ago and the first major glaciation at 2.15 Myr ago. Last, we tentatively infer that ice sheets may have grown largest during glacials with more modest reductions in deep-sea temperature.

Long-term cycles in the carbon reservoir of the Quaternary ocean: a perspective from the South China Sea

In recent years, long-term, high-resolution records from the deep sea and ice-cores have offered new research opportunities for Quaternary science. Paleoclimate studies are no longer restricted to individual glacial cycles, but extend to long-term (≥105 yr) processes across those cycles. Ocean Drilling Program Leg 184 of the South China Sea in 1999 uncovered well-preserved sediment sections, in which three long-term cycles in Pleistocene carbon isotope (δ13C) sequence have been found and demonstrated to be common in the global ocean. Subsequent discoveries confirm the existence of long-term processes of 105 yr in both the hydrologic (ice-sheet changes) and carbon (biogeochemical changes) cycles, posing the question whether the two processes are related. The present review shows that the long-eccentricity cycles prevail throughout the δ13C and other biogeochemical records in geologic history, and 400-kyr cycles in the oceanic δ13C sequence before the Quaternary can be hypothetically explained by changes in ratio between particulate and dissolved organic carbon (POC/DOC) in the ocean, depending on the monsoon-controlled nutrient supply. This is a ‘DOC hypothesis’. However, ocean restructuring at 1.6 Ma marked by the isolation of a sluggish abyss under the Southern Ocean has obscured the long-eccentricity 400-kyr signal in oceanic δ13C. The last million-year period has experienced two major changes in the climate regime, namely the mid-Pleistocene transition (MPT) centered at 0.9 Ma and the mid-Brunhes event (MBE) around 0.4 Ma. The MPT and MBE were preluded by δ13C maxima-III (δ13Cmax-III) ∼ 1.0 Ma and δ13Cmax-II ∼ 0.5 Ma, respectively. Together with similar hydroclimatic phenomena over corresponding glacial cycles, the two groups of hydrologic and biogeochemical events appear to have been driven largely by oceanographic changes in the Southern Ocean. Therefore, we interpret that the long-term biogeochemical processes originating from the Southern Ocean must have played a crucial role in Quaternary ice-sheet waxing and waning.

On the identification of a Pliocene time slice for data-model comparison

The characteristics of the mid-Pliocene warm period (mPWP: 3.264-3.025 Ma BP) have been examined using geological proxies and climate models. While there is agreement between models and data, details of regional climate differ. Uncertainties in prescribed forcings and in proxy data limit the utility of the interval to understand the dynamics of a warmer than present climate or evaluate models. This uncertainty comes, in part, from the reconstruction of a time slab rather than a time slice, where forcings required by climate models can be more adequately constrained. Here, we describe the rationale and approach for identifying a time slice(s) for Pliocene environmental reconstruction. A time slice centred on 3.205 Ma BP (3.204-3.207 Ma BP) has been identified as a priority for investigation. It is a warm interval characterized by a negative benthic oxygen isotope excursion (0.21-0.23‰) centred on marine isotope stage KM5c (KM5.3). It occurred during a period of orbital forcing that was very similar to present day. Climate model simulations indicate that proxy temperature estimates are unlikely to be significantly affected by orbital forcing for at least a precession cycle centred on the time slice, with the North Atlantic potentially being an important exception.

The multimillennial sea-level commitment of global warming

Global mean sea level has been steadily rising over the last century, is projected to increase by the end of this century, and will continue to rise beyond the year 2100 unless the current global mean temperature trend is reversed. Inertia in the climate and global carbon system, however, causes the global mean temperature to decline slowly even after greenhouse gas emissions have ceased, raising the question of how much sea-level commitment is expected for different levels of global mean temperature increase above preindustrial levels. Although sea-level rise over the last century has been dominated by ocean warming and loss of glaciers, the sensitivity suggested from records of past sea levels indicates important contributions should also be expected from the Greenland and Antarctic Ice Sheets. Uncertainties in the paleo-reconstructions, however, necessitate additional strategies to better constrain the sea-level commitment. Here we combine paleo-evidence with simulations from physical models to estimate the future sea-level commitment on a multimillennial time scale and compute associated regional sea-level patterns. Oceanic thermal expansion and the Antarctic Ice Sheet contribute quasi-linearly, with 0.4 m °C(-1) and 1.2 m °C(-1) of warming, respectively. The saturation of the contribution from glaciers is overcompensated by the nonlinear response of the Greenland Ice Sheet. As a consequence we are committed to a sea-level rise of approximately 2.3 m °C(-1) within the next 2,000 y. Considering the lifetime of anthropogenic greenhouse gases, this imposes the need for fundamental adaptation strategies on multicentennial time scales.

Patterns and mechanisms of early Pliocene warmth

About five to four million years ago, in the early Pliocene epoch, Earth had a warm, temperate climate. The gradual cooling that followed led to the establishment of modern temperature patterns, possibly in response to a decrease in atmospheric CO2 concentration, of the order of 100 parts per million, towards preindustrial values. Here we synthesize the available geochemical proxy records of sea surface temperature and show that, compared with that of today, the early Pliocene climate had substantially lower meridional and zonal temperature gradients but similar maximum ocean temperatures. Using an Earth system model, we show that none of the mechanisms currently proposed to explain Pliocene warmth can simultaneously reproduce all three crucial features. We suggest that a combination of several dynamical feedbacks underestimated in the models at present, such as those related to ocean mixing and cloud albedo, may have been responsible for these climate conditions.

Islands beneath islands: phylogeography of a groundwater amphipod crustacean in the Balearic archipelago

Background

Metacrangonyctidae (Amphipoda, Crustacea) is an enigmatic continental subterranean water family of marine origin (thalassoid). One of the species in the genus, Metacrangonyx longipes, is endemic to the Balearic islands of Mallorca and Menorca (W Mediterranean). It has been suggested that the origin and distribution of thalassoid crustaceans could be explained by one of two alternative hypotheses: (1) active colonization of inland freshwater aquifers by a marine ancestor, followed by an adaptative shift; or (2) passive colonization by stranding of ancestral marine populations in coastal aquifers during marine regressions. A comparison of phylogenies, phylogeographic patterns and age estimations of clades should discriminate in favour of one of these two proposals.

Results

Phylogenetic relationships within M. longipes based on three mitochondrial DNA (mtDNA) and one nuclear marker revealed five genetically divergent and geographically structured clades. Analyses of cytochrome oxidase subunit 1 (cox1) mtDNA data showed the occurrence of a high geographic population subdivision in both islands, with current gene flow occurring exclusively between sites located in close proximity. Molecular-clock estimations dated the origin of M. longipes previous to about 6 Ma, whereas major cladogenetic events within the species took place between 4.2 and 2.0 Ma.

Conclusions

M. longipes displayed a surprisingly old and highly fragmented population structure, with major episodes of cladogenesis within the species roughly correlating with some of the major marine transgression-regression episodes that affected the region during the last 6 Ma. Eustatic changes (vicariant events) -not active range expansion of marine littoral ancestors colonizing desalinated habitats-explain the phylogeographic pattern observed in M. longipes.

Phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs have been determined by geological processes and climate change in the Late Cenozoic

AIM: Our aims were to assess the phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs and to estimate divergence times using different geological scenarios. We related divergence times to past geological events and discuss the relevance of our data for the systematics of eastern Mediterranean water frogs. LOCATION: The eastern Mediterranean region. METHODS: Genetic diversity and divergence were calculated using sequences of two protein-coding mitochondrial (mt) genes: ND2 (1038 bp, 119 sequences) and ND3 (340 bp, 612 sequences). Divergence times were estimated in a Bayesian framework under four geological scenarios representing alternative possible geological histories for the eastern Mediterranean. We then compared the different scenarios using Bayes factors and additional geological data. RESULTS: Extensive genetic diversity in mtDNA divides eastern Mediterranean water frogs into six main haplogroups (MHG). Three MHGs were identified on the Anatolian mainland; the most widespread MHG with the highest diversity is distributed from western Anatolia to the northern shore of the Caspian Sea, including the type locality of Pelophylax ridibundus. The other two Anatolian MHGs are restricted to south-eastern Turkey, occupying localities west and east of the Amanos mountain range. One of the remaining three MHGs is restricted to Cyprus; a second to the Levant; the third was found in the distribution area of European lake frogs (P. ridibundus group), including the Balkans. MAIN CONCLUSIONS: Based on geological evidence and estimates of genetic divergence we hypothesize that the water frogs of Cyprus have been isolated from the Anatolian mainland populations since the end of the Messinian salinity crisis (MSC), i.e. since c. 5.5-5.3 Ma, while our divergence time estimates indicate that the isolation of Crete from the mainland populations (Peloponnese, Anatolia) most likely pre-dates the MSC. The observed rates of divergence imply a time window of c. 1.6-1.1 million years for diversification of the largest Anatolian MHG; divergence between the two other Anatolian MHGs may have begun about 3.0 Ma, apparently as a result of uplift of the Amanos Mountains. Our mtDNA data suggest that the Anatolian water frogs and frogs from Cyprus represent several undescribed species.

The Great American Biotic Interchange: Dispersals, Tectonics, Climate, Sea Level and Holding Pens

The biotic and geologic dynamics of the Great American Biotic Interchange are reviewed and revised. Information on the Marine Isotope Stage chronology, sea level changes as well as Pliocene and Pleistocene vegetation changes in Central and northern South America add to a discussion of the role of climate in facilitating trans-isthmian exchanges. Trans-isthmian land mammal exchanges during the Pleistocene glacial intervals appear to have been promoted by the development of diverse non-tropical ecologies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10914-010-9144-8) contains supplementary material, which is available to authorized users.

Tropical ocean temperatures over the past 3.5 million years

Determining the timing and amplitude of tropical sea surface temperature (SST) change is an important part of solving the puzzle of the Plio-Pleistocene ice ages. Alkenone-based tropical SST records from the major ocean basins show coherent glacial-interglacial temperature changes of 1 degrees to 3 degrees C that align with (but slightly lead) global changes in ice volume and deep ocean temperature over the past 3.5 million years. Tropical temperatures became tightly coupled with benthic delta18O and orbital forcing after 2.7 million years. We interpret the similarity of tropical SST changes, in dynamically dissimilar regions, to reflect "top-down" forcing through the atmosphere. The inception of a strong carbon dioxide-greenhouse gas feedback and amplification of orbital forcing at approximately 2.7 million years ago connected the fate of Northern Hemisphere ice sheets with global ocean temperatures since that time.