Why the Early Paleozoic was intrinsically prone to marine extinction

The geological record of marine animal biodiversity reflects the interplay between changing rates of speciation versus extinction. Compared to mass extinctions, background extinctions have received little attention. To disentangle the different contributions of global climate state, continental configuration, and atmospheric oxygen concentration (pO2) to variations in background extinction rates, we drive an animal physiological model with the environmental outputs from an Earth system model across intervals spanning the past 541 million years. We find that climate and continental configuration combined to make extinction susceptibility an order of magnitude higher during the Early Paleozoic than during the rest of the Phanerozoic, consistent with extinction rates derived from paleontological databases. The high extinction susceptibility arises in the model from the limited geographical range of marine organisms. It stands even when assuming present-day pO2, suggesting that increasing oxygenation through the Paleozoic is not necessary to explain why extinction rates apparently declined with time.

Why supercontinents became shorter lived as the Earth evolved

Periodic assembly and break-up of supercontinents since at least two billion years ago (Ga), dubbed the supercontinent cycle, provides the first-order tectonic control on the evolution of the Earth System including episodic orogenic events, mineralization, the formation and closure of oceans and superoceans, and even the evolution of life. However, the lifespan of the supercontinents appears to decrease with time, from ∼300 million years (Myr) for Nuna/Columbia, to 200-250 Myr for Rodinia and ∼150 Myr for the youngest supercontinent Pangaea. To understand what caused such a secular decrease in supercontinental lifespan, we conduct 3-D geodynamic modeling using realistic tectonic settings. The results show that the yield stress of newly formed orogens during the assembly of a supercontinent provides the dominant control on the lifespan of the supercontinent, implying that the yield stress of young orogens becomes lower with time. We hypothesize that the decreasing mantle temperature due to Earth's secular cooling might have caused new orogens to become weaker.

Combining palaeontological and neontological data shows a delayed diversification burst of carcharhiniform sharks likely mediated by environmental change

Estimating deep-time species-level diversification processes remains challenging. Both the fossil record and molecular phylogenies allow the estimation of speciation and extinction rates, but each type of data may still provide an incomplete picture of diversification dynamics. Here, we combine species-level palaeontological (fossil occurrences) and neontological (molecular phylogenies) data to estimate deep-time diversity dynamics through process-based birth-death models for Carcharhiniformes, the most speciose shark order today. Despite their abundant fossil record dating back to the Middle Jurassic, only a small fraction of extant carcharhiniform species is recorded as fossils, which impedes relying only on the fossil record to study their recent diversification. Combining fossil and phylogenetic data, we recover a complex evolutionary history for carcharhiniforms, exemplified by several variations in diversification rates with an early low diversity period followed by a Cenozoic radiation. We further reveal a burst of diversification in the last 30 million years, which is partially recorded with fossil data only. We also find that reef expansion and temperature change can explain variations in speciation and extinction through time. These results pinpoint the primordial importance of these environmental variables in the evolution of marine clades. Our study also highlights the benefit of combining the fossil record with phylogenetic data to address macroevolutionary questions.

Continental configuration controls ocean oxygenation during the Phanerozoic

The early evolutionary and much of the extinction history of marine animals is thought to be driven by changes in dissolved oxygen concentrations ([O₂]) in the ocean^1-3. In turn, [O₂] is widely assumed to be dominated by the geological history of atmospheric oxygen (pO₂)^4,5. Here, by contrast, we show by means of a series of Earth system model experiments how continental rearrangement during the Phanerozoic Eon drives profound variations in ocean oxygenation and induces a fundamental decoupling in time between upper-ocean and benthic [O₂]. We further identify the presence of state transitions in the global ocean circulation, which lead to extensive deep-ocean anoxia developing in the early Phanerozoic even under modern pO₂. Our finding that ocean oxygenation oscillates over stable thousand-year (kyr) periods also provides a causal mechanism that might explain elevated rates of metazoan radiation and extinction during the early Palaeozoic Era⁶. The absence, in our modelling, of any simple correlation between global climate and ocean ventilation, and the occurrence of profound variations in ocean oxygenation independent of atmospheric pO₂, presents a challenge to the interpretation of marine redox proxies, but also points to a hitherto unrecognized role for continental configuration in the evolution of the biosphere.

Dataset of Phanerozoic continental climate and Köppen–Geiger climate classes

This article describes a suite of global climate model output files that provide continental climatic conditions (monthly temperatures, precipitation, evaporation, precipitation minus evaporation balance, runoff) together with the calculated Köppen–Geiger climate classes and topography, for 28 evenly spaced time slices through the Phanerozoic (Cambrian to Quaternary, 540 Ma to 0 Ma). Climatic variables were simulated with the Fast Ocean Atmosphere Model (FOAM), using a recent set of open-access continental reconstructions with paleotopography and recent atmospheric CO₂ and solar luminosity estimates. FOAM is a general circulation model frequently used in paleoclimate studies, especially in the Palaeozoic. Köppen–Geiger climate classes were calculated based on simulated temperature and precipitation fields using Wong Hearing et al.’s [1] implementation of Peel et al.’s [2] updated classification. This dataset provides a unique window onto changing continental climate throughout the Phanerozoic that accounts for the simultaneous evolution of paleogeography (continental configuration and topography), atmospheric composition and greenhouse gas forcing, and solar luminosity.

Post-extinction recovery of the Phanerozoic oceans and biodiversity hotspots

The fossil record of marine invertebrates has long fuelled the debate as to whether or not there are limits to global diversity in the sea1-5. Ecological theory states that, as diversity grows and ecological niches are filled, the strengthening of biological interactions imposes limits on diversity6,7. However, the extent to which biological interactions have constrained the growth of diversity over evolutionary time remains an open question1-5,8-11. Here we present a regional diversification model that reproduces the main Phanerozoic eon trends in the global diversity of marine invertebrates after imposing mass extinctions. We find that the dynamics of global diversity are best described by a diversification model that operates widely within the exponential growth regime of a logistic function. A spatially resolved analysis of the ratio of diversity to carrying capacity reveals that less than 2% of the global flooded continental area throughout the Phanerozoic exhibits diversity levels approaching ecological saturation. We attribute the overall increase in global diversity during the Late Mesozoic and Cenozoic eras to the development of diversity hotspots under prolonged conditions of Earth system stability and maximum continental fragmentation. We call this the 'diversity hotspots hypothesis', which we propose as a non-mutually exclusive alternative to the hypothesis that the Mesozoic marine revolution led this macroevolutionary trend12,13.

Cenozoic evolution of deep ocean temperature from clumped isotope thermometry

Characterizing past climate states is crucial for understanding the future consequences of ongoing greenhouse gas emissions. Here, we revisit the benchmark time series for deep ocean temperature across the past 65 million years using clumped isotope thermometry. Our temperature estimates from the deep Atlantic Ocean are overall much warmer compared with oxygen isotope-based reconstructions, highlighting the likely influence of changes in deep ocean pH and/or seawater oxygen isotope composition on classical oxygen isotope records of the Cenozoic. In addition, our data reveal previously unrecognized large swings in deep ocean temperature during early Eocene acute greenhouse warmth. Our results call for a reassessment of the Cenozoic history of ocean temperatures to achieve a more accurate understanding of the nature of climatic responses to tectonic events and variable greenhouse forcing.

Evidence for active upper mantle flow in the Atlantic and Indo-Australian realms since the Upper Jurassic from hiatus maps and spreading rate changes

Histories of large-scale horizontal and vertical lithosphere motion hold important information on mantle convection. Here, we compare continent-scale hiatus maps as a proxy for mantle flow induced dynamic topography and plate motion variations in the Atlantic and Indo-Australian realms since the Upper Jurassic, finding they frequently correlate, except when plate boundary forces may play a significant role. This correlation agrees with descriptions of asthenosphere flow beneath tectonic plates in terms of Poiseuille/Couette flow, as it explicitly relates plate motion changes, induced by evolving basal shear forces, to non-isostatic vertical motion of the lithosphere. Our analysis reveals a timescale, on the order of a geological series, between the occurrence of continent-scale hiatus and plate motion changes. This is consistent with the presence of a weak upper mantle. It also shows a spatial scale for interregional hiatus, on the order of 2000-3000 km in diameter, which can be linked by fluid dynamic analysis to active upper mantle flow regions. Our results suggest future studies should pursue large-scale horizontal and vertical lithosphere motion in combination, to track the expressions of past mantle flow. Such studies would provide powerful constraints for adjoint-based geodynamic inverse models of past mantle convection.

Eutrophication and Deoxygenation Forcing of Marginal Marine Organic Carbon Burial During the PETM

The Paleocene-Eocene Thermal Maximum (PETM) is recognized globally by a negative excursion in stable carbon isotope ratios (δ¹³C) in sedimentary records, termed the carbon isotope excursion (CIE). Based on the CIE, the cause, duration, and mechanisms of recovery of the event have been assessed. Here, we focus on the role of increased organic carbon burial on continental margins as a key driver of CO₂ drawdown and global exogenic δ¹³C during the recovery phase. Using new and previously published sediment proxy data, we show evidence for widespread enhanced primary production, low oxygen waters, and high organic carbon (C_org) burial in marginal and restricted environments throughout the δ¹³C excursion. With a new biogeochemical box model for deep and marginal environments, we show that increased phosphorus availability and water column stratification on continental margins can explain the increased C_org burial during the PETM. Deoxygenation and recycling of phosphorus relative to C_org were relatively mild, compared to modern day anoxic marine systems. Our model reproduces the conditions reconstructed by field data, resulting in a burial of 6,000 Pg across the PETM, in excess of late Paleocene burial, and ∼3,300 Pg C for the critical first 40 kyr of the recovery, primarily located on continental margins. This value is consistent with prior data and model estimates (∼2,000-3,000 Pg C). To reproduce global exogenic δ¹³C patterns, this C_org burial implies an injection of 5,000-10,000 Pg C during the first ∼100-150 kyr of the PETM, depending on the source's δ¹³C (-11‰ to -55‰).

Running across the Silurian/Devonian Boundary along Northern Gondwana: A Conodont Perspective

The Global Stratotype Section and Point (GSSP) of the Silurian/Devonian boundary, Lower Devonian Series and Lochkovian Stage was formally placed in 1977 at Klonk, in the Czech Republic, at the first appearance of the graptolite Uncinatograptus uniformis uniformis (Přibyl). However, since then, correlation of this limit has been often hampered in carbonate facies where graptolites are uncommon or totally absent. A large calcareous deposition occurred at the Silurian/Devonian boundary along the northern and peri-Gondwana margin, thus representing an ideal location to select and test a possible additional biostratigraphic marker of the limit among conodonts. The first appearance of Caudicriodus hesperius almost simultaneously at the base of the Devonian in Bohemia, the Carnic Alps, Sardinia, Morocco and elsewhere indicates that this taxon is the conodont that best approximates the beginning of the Period. The first or last appearance of other species (e.g., Ozarkodina confluens, Zieglerodina klonkensis, Z. remscheidensis and Caudicriodus woschmidti) may help to recognise the boundary as well.

The Deep Rocky Biosphere: New Geomicrobiological Insights and Prospects

Rocks that react with liquid water are widespread but spatiotemporally limited throughout the solar system, except for Earth. Rock-forming minerals with high iron content and accessory minerals with high amounts of radioactive elements are essential to support rock-hosted microbial life by supplying organics, molecular hydrogen, and/or oxidants. Recent technological advances have broadened our understanding of the rocky biosphere, where microbial inhabitation appears to be difficult without nutrient and energy inputs from minerals. In particular, microbial proliferation in igneous rock basements has been revealed using innovative geomicrobiological techniques. These recent findings have dramatically changed our perspective on the nature and the extent of microbial life in the rocky biosphere, microbial interactions with minerals, and the influence of external factors on habitability. This study aimed to gather information from scientific and/or technological innovations, such as omics-based and single-cell level characterizations, targeting deep rocky habitats of organisms with minimal dependence on photosynthesis. By synthesizing pieces of rock-hosted life, we can explore the evo-phylogeny and ecophysiology of microbial life on Earth and the life’s potential on other planetary bodies.

Long-Term Earth-Moon Evolution With High-Level Orbit and Ocean Tide Models

Tides and Earth-Moon system evolution are coupled over geological time. Tidal energy dissipation on Earth slows E a r t h ' s rotation rate, increases obliquity, lunar orbit semi-major axis and eccentricity, and decreases lunar inclination. Tidal and core-mantle boundary dissipation within the Moon decrease inclination, eccentricity and semi-major axis. Here we integrate the Earth-Moon system backwards for 4.5 Ga with orbital dynamics and explicit ocean tide models that are "high-level" (i.e., not idealized). To account for uncertain plate tectonic histories, we employ Monte Carlo simulations, with tidal energy dissipation rates (normalized relative to astronomical forcing parameters) randomly selected from ocean tide simulations with modern ocean basin geometry and with 55, 116, and 252 Ma reconstructed basin paleogeometries. The normalized dissipation rates depend upon basin geometry and E a r t h ' s rotation rate. Faster Earth rotation generally yields lower normalized dissipation rates. The Monte Carlo results provide a spread of possible early values for the Earth-Moon system parameters. Of consequence for ocean circulation and climate, absolute (un-normalized) ocean tidal energy dissipation rates on the early Earth may have exceeded t o d a y ' s rate due to a closer Moon. Prior to ∼ 3 Ga , evolution of inclination and eccentricity is dominated by tidal and core-mantle boundary dissipation within the Moon, which yield high lunar orbit inclinations in the early Earth-Moon system. A drawback for our results is that the semi-major axis does not collapse to near-zero values at 4.5 Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation.

Influence of historical changes in tropical reef habitat on the diversification of coral reef fishes

Past environmental changes are expected to have profoundly impacted diversity dynamics through time. While some previous studies showed an association between past climate changes or tectonic events and important shifts in lineage diversification, it is only recently that past environmental changes have been explicitly integrated in diversification models to test their influence on diversification rates. Here, we used a global reconstruction of tropical reef habitat dynamics during the Cenozoic and phylogenetic diversification models to test the influence of (i) major geological events, (ii) reef habitat fragmentation and (iii) reef area on the diversification of 9 major clades of tropical reef fish (Acanthuridae, Balistoidea, Carangoidea, Chaetodontidae, Haemulinae, Holocentridae, Labridae, Pomacentridae and Sparidae). The diversification models revealed a weak association between paleo-habitat changes and diversification dynamics. Specifically, the fragmentation of tropical reef habitats over the Cenozoic was found to be a driver of tropical reef fish diversification for 2 clades. However, overall, our approach did not allow the identification of striking associations between diversification dynamics and paleo-habitat fragmentation in contrast with theoretical model's predictions.

A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins

Living turtles are characterized by extraordinarily low species diversity given their age. The clade's extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene-Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene-Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossil-bearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities.

Investigating Mesozoic Climate Trends and Sensitivities With a Large Ensemble of Climate Model Simulations

The Mesozoic era (∼252 to 66 million years ago) was a key interval in Earth's evolution toward its modern state, witnessing the breakup of the supercontinent Pangaea and significant biotic innovations like the early evolution of mammals. Plate tectonic dynamics drove a fundamental climatic transition from the early Mesozoic supercontinent toward the Late Cretaceous fragmented continental configuration. Here, key aspects of Mesozoic long-term environmental changes are assessed in a climate model ensemble framework. We analyze so far the most extended ensemble of equilibrium climate states simulated for evolving Mesozoic boundary conditions covering the period from 255 to 60 Ma in 5 Myr timesteps. Global mean temperatures are generally found to be elevated above the present and exhibit a baseline warming trend driven by rising sea levels and increasing solar luminosity. Warm (Triassic and mid-Cretaceous) and cool (Jurassic and end-Cretaceous) anomalies result from pCO2 changes indicated by different reconstructions. Seasonal and zonal temperature contrasts as well as continental aridity show an overall decrease from the Late Triassic-Early Jurassic to the Late Cretaceous. Meridional temperature gradients are reduced at higher global temperatures and less land area in the high latitudes. With systematic sensitivity experiments, the influence of paleogeography, sea level, vegetation patterns, pCO2, solar luminosity, and orbital configuration on these trends is investigated. For example, long-term seasonality trends are driven by paleogeography, but orbital cycles could have had similar-scale effects on shorter timescales. Global mean temperatures, continental humidity, and meridional temperature gradients are, however, also strongly affected by pCO2.

Reconciling atmospheric CO2, weathering, and calcite compensation depth across the Cenozoic

The Cenozoic era (66 to 0 million years) is marked by long-term aberrations in carbon cycling and large climatic shifts, some of which challenge the current understanding of carbon cycle dynamics. Here, we investigate possible mechanisms responsible for the observed long-term trends by using a novel approach that features a full-fledged ocean carbonate chemistry model. Using a compilation of pCO2, pH, and calcite compensation depth (CCD) observational evidence and a suite of simulations, we reconcile long-term Cenozoic climate and CCD trends. We show that the CCD response was decoupled from changes in silicate and carbonate weathering rates, challenging the continental uplift hypothesis. The two dominant mechanisms for decoupling are shelf-basin carbonate burial fractionation combined with proliferation of pelagic calcifiers. The temperature effect on remineralization rates of marine organic matter also plays a critical role in controlling the carbon cycle dynamics, especially during the warmer periods of the Cenozoic.

Continent-scale Hiatus Maps for the Atlantic Realm and Australia since the Upper Jurassic and links to mantle flow induced dynamic topography

Interregional geological maps hold important information for geodynamic models. Here, we use such maps to visualize major conformable and unconformable contacts at interregional scales and at the level of geologic series from the Upper Jurassic onward across North and South America, Europe, Africa and Australia. We extract hiatus information from these paleogeological maps, which we plot in a paleogeographical reference frame to link the maps to the plate and plume modes of mantle convection. We assume that interregional patterns of hiatus surfaces are proxy records of continent-scale mantle-induced vertical motion of the lithosphere. We find significant differences in the distribution of hiatus across and between continents at the timescale of geologic series, that is ten to a few tens of millions of years (Myrs). This is smaller than the mantle transit time, which, as the timescale of convection, is about 100-200 Myrs. Our results imply that different timescales for convection and topography in convective support must be an integral component of time-dependent geodynamic Earth models, consistent with the presence of a weaker upper mantle relative to the lower mantle. Additional geological constraints together with interregional geological maps at the resolution of stages (1-2 Myrs), are needed to assist in future geodynamic interpretations of interregional geologic hiatus.

A possible link between coral reef success, crustose coralline algae and the evolution of herbivory

Crustose coralline red algae (CCA) play a key role in the consolidation of many modern tropical coral reefs. It is unclear, however, if their function as reef consolidators was equally pronounced in the geological past. Using a comprehensive database on ancient reefs, we show a strong correlation between the presence of CCA and the formation of true coral reefs throughout the last 150 Ma. We investigated if repeated breakdowns in the potential capacity of CCA to spur reef development were associated with sea level, ocean temperature, CO₂ concentration, CCA species diversity, and/or the evolution of major herbivore groups. Model results show that the correlation between the occurrence of CCA and the development of true coral reefs increased with CCA diversity and cooler ocean temperatures while the diversification of herbivores had a transient negative effect. The evolution of novel herbivore groups compromised the interaction between CCA and true reef growth at least three times in the investigated time interval. These crises have been overcome by morphological adaptations of CCA.

Marine Os isotopic evidence for multiple volcanic episodes during Cretaceous Oceanic Anoxic Event 1b

The Aptian-Albian boundary is marked by one of the major oceanic perturbations during the Cretaceous, called Oceanic Anoxic Event (OAE) 1b. Extensive volcanic episodes at the Southern Kerguelen Plateau has been suggested as the trigger of OAE1b, but compelling evidence remains lacking. Here, we reconstructed the temporal variations of marine Os isotopic ratios across the Aptian-Albian boundary in the Tethyan and Pacific pelagic sedimentary records to elucidate the causal links between OAE1b, the biotic turnover, and volcanic episodes. Our new Os isotopic records show two negative spikes that correlate with a period of planktonic foraminiferal turnover across the Aptian-Albian boundary during OAE1b and suggest multiple submarine volcanic events. By comparing our Os isotopic profile with carbon isotopic compositions of carbonate, CaCO₃ content, and the relative abundances of agglutinated foraminifera, we conclude that ocean acidification caused by the massive release of CO₂ through extensive volcanic episodes could have promoted the major planktonic foraminiferal turnover during OAE1b.

Temperate rainforests near the South Pole during peak Cretaceous warmth

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years^1-5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume⁶. In the near absence of proximal geological records from south of the Antarctic Circle, it is disputed whether polar ice could exist under such environmental conditions. Here we use a sedimentary sequence recovered from the West Antarctic shelf-the southernmost Cretaceous record reported so far-and show that a temperate lowland rainforest environment existed at a palaeolatitude of about 82° S during the Turonian-Santonian age (92 to 83 million years ago). This record contains an intact 3-metre-long network of in situ fossil roots embedded in a mudstone matrix containing diverse pollen and spores. A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric carbon dioxide concentrations of 1,120-1,680 parts per million by volume and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo under high levels of atmospheric carbon dioxide.

A record of seafloor methane seepage across the last 150 million years

Seafloor methane seepage is a significant source of carbon in the marine environment. The processes and temporal patterns of seafloor methane seepage over multi-million-year time scales are still poorly understood. The microbial oxidation of methane can store carbon in sediments through precipitation of carbonate minerals, thus providing a record of past methane emission. In this study, we compiled data on methane-derived carbonates to build a proxy time series of methane emission over the last 150 My and statistically compared it with the main hypothesised geological controllers of methane emission. We quantitatively demonstrate that variations in sea level and organic carbon burial are the dominant controls on methane leakage since the Early Cretaceous. Sea level controls methane seepage variations by imposing smooth trends on timescales in the order of tens of My. Organic carbon burial is affected by the same cyclicities, and instantaneously controls methane release because of the geologically rapid generation of biogenic methane. Both the identified fundamental (26–27 My) and higher (12 My) cyclicities relate to global phenomena. Temporal correlation analysis supports the evidence that modern expansion of hypoxic areas and its effect on organic carbon burial may lead to higher seawater methane concentrations over the coming centuries.

Ecological constraints coupled with deep-time habitat dynamics predict the latitudinal diversity gradient in reef fishes

We develop a spatially explicit model of diversification based on palaeohabitat to explore the predictions of four major hypotheses potentially explaining the latitudinal diversity gradient (LDG), namely, the 'time-area', 'tropical niche conservatism', 'ecological limits' and 'evolutionary speed' hypotheses. We compare simulation outputs to observed diversity gradients in the global reef fish fauna. Our simulations show that these hypotheses are non-mutually exclusive and that their relative influence depends on the time scale considered. Simulations suggest that reef habitat dynamics produced the LDG during deep geological time, while ecological constraints shaped the modern LDG, with a strong influence of the reduction in the latitudinal extent of tropical reefs during the Neogene. Overall, this study illustrates how mechanistic models in ecology and evolution can provide a temporal and spatial understanding of the role of speciation, extinction and dispersal in generating biodiversity patterns.

The herophile species group of Calisto (Lepidoptera : Nymphalidae : Satyrinae), new taxa and historical biogeography

The genus Calisto is endemic tothe West Indiesand the only representative there of the Satyrinae. Here wereconstruct the evolutionary relationshipsof the herophile group and describe five new species from Cuba: Calisto gundlachi sp. nov., Calisto siguanensis sp. nov., Calisto disjunctus sp. nov., Calisto sharkeyae sp. nov. and Calisto lastrai sp. nov.We employ one mitochondrial and four nuclear markers to assess the phylogenetic position, Maximum Likelihood and Bayesian Inference approaches, of the new taxa. Our phylogenetic trees yielded two strongly supported main clades with four of the new species included within them and C. sharkeyae as sister group to the rest of the major main clade. We conduct time-divergence estimations and ancestral area reconstructions using BEAST and BioGeoBEARS. The group originated 12.15 million years ago during the middle Miocene in north-eastern Cuba, Nipe-Sagua-Baracoa Massif. After 6 million years of in situ evolution most lineages started to colonise other Cuban territories and the Bahamas. This scenario is consistent with key geological events, including the closure of the western Havana–Matanzas channel 8–6 million years ago, the uplift of the Sierra Maestra 6–5 million years ago, and the land connections among Cuban regions during the Miocene–Pleistocene sea level drops. Dispersal and vicariance processes may have occurred, with populations surviving floodings on the major and minor mountain ranges, which remained as ‘islands’. http://zoobank.org/urn:lsid:zoobank.org:act:03690F79-F938-42A0-B234-4A228D5C1913

The Latitudinal Diversity Gradient: Novel Understanding through Mechanistic Eco-evolutionary Models

The latitudinal diversity gradient (LDG) is one of the most widely studied patterns in ecology, yet no consensus has been reached about its underlying causes. We argue that the reasons for this are the verbal nature of existing hypotheses, the failure to mechanistically link interacting ecological and evolutionary processes to the LDG, and the fact that empirical patterns are often consistent with multiple explanations. To address this issue, we synthesize current LDG hypotheses, uncovering their eco-evolutionary mechanisms, hidden assumptions, and commonalities. Furthermore, we propose mechanistic eco-evolutionary modeling and an inferential approach that makes use of geographic, phylogenetic, and trait-based patterns to assess the relative importance of different processes for generating the LDG.

Major intensification of Atlantic overturning circulation at the onset of Paleogene greenhouse warmth

During the Late Cretaceous and early Cenozoic the Earth experienced prolonged climatic cooling most likely caused by decreasing volcanic activity and atmospheric CO₂ levels. However, the causes and mechanisms of subsequent major global warming culminating in the late Paleocene to Eocene greenhouse climate remain enigmatic. We present deep and intermediate water Nd-isotope records from the North and South Atlantic to decipher the control of the opening Atlantic Ocean on ocean circulation and its linkages to the evolution of global climate. The marked convergence of Nd-isotope signatures 59 million years ago indicates a major intensification of deep-water exchange between the North and South Atlantic, which coincided with the turning point of deep-water temperatures towards early Paleogene warming. We propose that this intensification of Atlantic overturning circulation in concert with increased atmospheric CO₂ from continental rifting marked a climatic tipping point contributing to a more efficient distribution of heat over the planet.

Evidence for Basement Reactivation during the Opening of the Labrador Sea from the Makkovik Province, Labrador, Canada: Insights from Field Data and Numerical Models

The onshore exposures adjacent to modern, offshore passive continental margins may preserve evidence of deformation from the pre-, syn-, and post-rift phases of continental breakup that allow us to investigate the processes associated with and controlling rifting and breakup. Here, we characterize onshore brittle deformation and pre-rift basement metamorphic mineral fabric from onshore Labrador in Eastern Canada in the Palaeoproterozoic Aillik Domain of the Makkovik Province. Stress inversion (1) was applied to these data and then compared to (2) numerical models of hybrid slip and dilation tendency, (3) independent calculations of the regional geopotential stress field, and (4) analyses of palaeo-stress in proximal regions from previous work. The stress inversion shows well-constrained extensional deformation perpendicular to the passive margin, likely related to pre-breakup rifting in the proto-Labrador Sea. Hybrid slip and dilatation analysis indicates that inherited basement structures were likely oriented in a favorable orientation to be reactivated during rifting. Reconstructed geopotential stresses illuminate changes of the ambient stress field over time and confirm the present paleo-stress estimates. The new results and numerical models provide a consistent picture of the late Mesozoic-Cenozoic lithospheric stress field evolution in the Labrador Sea region. The proto-Labrador Sea region was characterized by a persistent E–W (coast-perpendicular) extensional stress regime, which we interpret as the pre-breakup continental rifting that finally led to continental breakup. Later, the ridge push of the Labrador Sea spreading ridge maintained this general direction of extension. We see indications for anti-clockwise rotation of the direction of extension along some of the passive margins. However, extreme persistent N–S-oriented extension as indicated by studies further north in West Greenland cannot be confirmed.

Volcanic ash as a driver of enhanced organic carbon burial in the Cretaceous

On greater than million year timescales, carbon in the ocean-atmosphere-biosphere system is controlled by geologic inputs of CO₂ through volcanic and metamorphic degassing. High atmospheric CO₂ and warm climates in the Cretaceous have been attributed to enhanced volcanic emissions of CO₂ through more rapid spreading at mid-ocean ridges and, in particular, to a global flare-up in continental arc volcanism. Here, we show that global flare-ups in continental arc magmatism also enhance the global flux of nutrients into the ocean through production of windblown ash. We show that up to 75% of Si, Fe and P is leached from windblown ash during and shortly after deposition, with soluble Si, Fe and P inputs from ash alone in the Cretaceous being higher than the combined input of dust and rivers today. Ash-derived nutrient inputs may have increased the efficiency of biological productivity and organic carbon preservation in the Cretaceous, possibly explaining why the carbon isotopic signature of Cretaceous seawater was high. Variations in volcanic activity, particularly continental arcs, have the potential of profoundly altering carbon cycling at the Earth’s surface by increasing inputs of CO₂ and ash-borne nutrients, which together enhance biological productivity and burial of organic carbon, generating an abundance of hydrocarbon source rocks.

Zircon evidence for incorporation of terrigenous sediments into the magma source of continental basalts

Crustal components may be incorporated into continental basalts by either shallow contamination or deep mixing. While the former proceeds at crustal depths with common preservation of refractory minerals, the latter occurs at mantle depths with rare survival of relict minerals. Discrimination between the two mechanisms has great bearing to subcontinental mantle geochemistry. Here we report the occurrence of relict zircons in Cenozoic continental basalts from eastern China. A combined study of zircon U-Pb ages and geochemistry indicates that detrital zircons were carried by terrigenous sediments into a subcontinental subduction zone, where the zircon were transferred by fluids into the magma sources of continental basalts. The basalts were sampled from three petrotectonic units with distinct differences in their magmatic and metamorphic ages, making the crustal contamination discernible. The terrigenous sediments were carried by the subducting oceanic crust into the asthenospheric mantle, producing both soluble and insoluble materials at the slab-mantle interface. These materials were served as metasomatic agents to react with the overlying mantle wedge peridotite, generating a kind of ultramafic metasomatites that contain the relict zircons. Therefore, the occurrence of relict zircons in continental basalts indicates that this refractory mineral can survive extreme temperature-pressure conditions in the asthenospheric mantle.

Trophic and tectonic limits to the global increase of marine invertebrate diversity

The marine invertebrate fossil record provides the most comprehensive history of how the diversity of animal life has evolved through time. One of the main features of this record is a modest rise in diversity over nearly a half-billion years. The long-standing view is that ecological interactions such as resource competition and predation set upper limits to global diversity, which, in the absence of external perturbations, is maintained indefinitely at equilibrium. However, the effect of mechanisms associated with the history of the seafloor, and their influence on the creation and destruction of marine benthic habitats, has not been explored. Here we use statistical methods for causal inference to investigate the drivers of marine invertebrate diversity dynamics through the Phanerozoic. We find that diversity dynamics responded to secular variations in marine food supply, substantiating the idea that global species richness is regulated by resource availability. Once diversity was corrected for changes in food resource availability, its dynamics were causally linked to the age of the subducting oceanic crust. We suggest that the time elapsed between the formation (at mid-ocean ridges) and destruction (at subduction zones) of ocean basins influences the diversity dynamics of marine invertebrates and may have contributed to constrain their diversification.

Macrostratigraphy and Its Promise for Paleobiology

Macrostratigraphy is the study and statistical analysis of sediment packages that formed continuously at a specified scale of temporal resolution and that are bound by gaps recognizable at that same scale. The temporal ranges of gap-bound packages, compiled separately for different geographic locations, permit area-weighted, survivorship-based measures of rock quantity and spatio-temporal environmental continuity to be measured. Analytical basin fill models suggest that the parameters controlling sedimentation and sequence stratigraphic architecture, such as base level and sediment supply, can be detected quantitatively by macrostratigraphy.

Macrostratigraphic analysis of the marine sedimentary rock record in the United States at a temporal resolution of ~106 years reproduces most of the well-known Sloss sequences, but it also identifies two prominent megasequences, the Paleozoic and Modern megasequences, which are separated by a Permian-Triassic discontinuity and Phanerozoic minimum in rock quantity. Many short- and long-term features of the macroevolutionary history of marine animals are reproduced by macrostratigraphy, including 1) many patterns in genus richness, 2) patterns in rates of genus extinction and, to a lesser degree, rates of origination, and 3) patterns of extinction selectivity and the shifting relative richness of Sepkoski's Paleozoic and Modern evolutionary faunas. The extent to which macrostratigraphy reproduces the macroevolutionary history of marine animals transcends what is expected by geologically-controlled sampling biases. Instead, the processes which control the spatio-temporal dynamics of shelf sedimentation, including expansions and contractions of shallow epicontinental seas, have probably exerted a consistent influence on the macroevolutionary history of marine animals. Exploring the common cause hypothesis by putting fossils back into rocks and rocks into a new quantitative framework for physical environmental change holds considerable promise for paleobiology.

Constraining climate sensitivity and continental versus seafloor weathering using an inverse geological carbon cycle model

The relative influences of tectonics, continental weathering and seafloor weathering in controlling the geological carbon cycle are unknown. Here we develop a new carbon cycle model that explicitly captures the kinetics of seafloor weathering to investigate carbon fluxes and the evolution of atmospheric CO₂ and ocean pH since 100 Myr ago. We compare model outputs to proxy data, and rigorously constrain model parameters using Bayesian inverse methods. Assuming our forward model is an accurate representation of the carbon cycle, to fit proxies the temperature dependence of continental weathering must be weaker than commonly assumed. We find that 15–31 °C (1σ) surface warming is required to double the continental weathering flux, versus 3–10 °C in previous work. In addition, continental weatherability has increased 1.7–3.3 times since 100 Myr ago, demanding explanation by uplift and sea-level changes. The average Earth system climate sensitivity is  K (1σ) per CO₂ doubling, which is notably higher than fast-feedback estimates. These conclusions are robust to assumptions about outgassing, modern fluxes and seafloor weathering kinetics.

The influence of tectonics, continental weathering, and seafloor weathering in the geological carbon cycle remain unclear. Here, the authors develop a new carbon cycle model and, through comparison with proxy data, critically evaluate the influence of these components on carbon fluxes since 100 Ma.

Australian shelf sediments reveal shifts in Miocene Southern Hemisphere westerlies

Global climate underwent a major reorganization when the Antarctic ice sheet expanded ~14 million years ago (Ma) (1). This event affected global atmospheric circulation, including the strength and position of the westerlies and the Intertropical Convergence Zone (ITCZ), and, therefore, precipitation patterns (2-5). We present new shallow-marine sediment records from the continental shelf of Australia (International Ocean Discovery Program Sites U1459 and U1464) providing the first empirical evidence linking high-latitude cooling around Antarctica to climate change in the (sub)tropics during the Miocene. We show that Western Australia was arid during most of the Middle Miocene. Southwest Australia became wetter during the Late Miocene, creating a climate gradient with the arid interior, whereas northwest Australia remained arid throughout. Precipitation and river runoff in southwest Australia gradually increased from 12 to 8 Ma, which we relate to a northward migration or intensification of the westerlies possibly due to increased sea ice in the Southern Ocean (5). Abrupt aridification indicates that the westerlies shifted back to a position south of Australia after 8 Ma. Our midlatitude Southern Hemisphere data are consistent with the inference that expansion of sea ice around Antarctica resulted in a northward movement of the westerlies. In turn, this may have pushed tropical atmospheric circulation and the ITCZ northward, shifting the main precipitation belt over large parts of Southeast Asia (4).

Kinematics and dynamics of the East Pacific Rise linked to a stable, deep-mantle upwelling

Earth's tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth's dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. The mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

Tracking the spatiotemporal variations of statistically independent components involving enrichment of rare-earth elements in deep-sea sediments

Deep-sea sediments have attracted much attention as a promising resource for rare-earth elements and yttrium (REY). In this study, we show statistically independent components characterising REY-enrichment in the abyssal ocean that are decoded by Independent Component Analysis of a multi-elemental dataset of 3,968 bulk sediment samples from 101 sites in the Pacific and Indian oceans. This study for the first time reconstructs the spatiotemporal variations of the geochemical signatures, including hydrothermal, hydrogenous, and biogenic calcium phosphate components that were closely involved in the formation of REY-rich mud over the past 65 million years. An underlying key factor of significant REY-enrichment is a sufficiently low sedimentation rate that enables the mud to accumulate REY from seawater. In the early Cenozoic, a remarkably small supply of aeolian dust, compared with any other time and region, facilitated the deposition of very high-grade REY-rich mud in the South Pacific. This indicates an important link between the genesis of the seafloor mineral resources and Earth’s dynamic phenomena such as climate change and plate tectonics.

On the enigmatic birth of the Pacific Plate within the Panthalassa Ocean

The oceanic Pacific Plate started forming in Early Jurassic time within the vast Panthalassa Ocean that surrounded the supercontinent Pangea, and contains the oldest lithosphere that can directly constrain the geodynamic history of the circum-Pangean Earth. We show that the geometry of the oldest marine magnetic anomalies of the Pacific Plate attests to a unique plate kinematic event that sparked the plate's birth at virtually a point location, surrounded by the Izanagi, Farallon, and Phoenix Plates. We reconstruct the unstable triple junction that caused the plate reorganization, which led to the birth of the Pacific Plate, and present a model of the plate tectonic configuration that preconditioned this event. We show that a stable but migrating triple junction involving the gradual cessation of intraoceanic Panthalassa subduction culminated in the formation of an unstable transform-transform-transform triple junction. The consequent plate boundary reorganization resulted in the formation of a stable triangular three-ridge system from which the nascent Pacific Plate expanded. We link the birth of the Pacific Plate to the regional termination of intra-Panthalassa subduction. Remnants thereof have been identified in the deep lower mantle of which the locations may provide paleolongitudinal control on the absolute location of the early Pacific Plate. Our results constitute an essential step in unraveling the plate tectonic evolution of "Thalassa Incognita" that comprises the comprehensive Panthalassa Ocean surrounding Pangea.

Plate tectonics drive tropical reef biodiversity dynamics

The Cretaceous breakup of Gondwana strongly modified the global distribution of shallow tropical seas reshaping the geographic configuration of marine basins. However, the links between tropical reef availability, plate tectonic processes and marine biodiversity distribution patterns are still unknown. Here, we show that a spatial diversification model constrained by absolute plate motions for the past 140 million years predicts the emergence and movement of diversity hotspots on tropical reefs. The spatial dynamics of tropical reefs explains marine fauna diversification in the Tethyan Ocean during the Cretaceous and early Cenozoic, and identifies an eastward movement of ancestral marine lineages towards the Indo-Australian Archipelago in the Miocene. A mechanistic model based only on habitat-driven diversification and dispersal yields realistic predictions of current biodiversity patterns for both corals and fishes. As in terrestrial systems, we demonstrate that plate tectonics played a major role in driving tropical marine shallow reef biodiversity dynamics.

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.

Biogeography of Drosophila (Diptera: Drosophilidae) in East and Southeast Asia

The causes of high biological diversity in biodiversity hotspots have long been a major subject of study in conservation biology. To investigate this matter, we conducted a phylogeographic study of five Drosophila (Diptera: Drosophilidae) species from East and Southeast Asia: Drosophila albomicans Duda, D. formosana Duda, D. immigrans Sturtevant, D. melanogaster Meigen, and D. simulans Sturtevant. We collected 185 samples from 28 localities in eight countries. From each collected individual, we sequenced the autosomal extra sex comb gene (esc) and seven mitochondrial genes, including nicotinamide adenine dinucleotide hydrate-reductase dehydrogenase subunit 4 (ND4), ND4L, tRNA-His, tRNA-Pro, tRNA-Thr, partial ND5, and partial ND6. Phylogenetic analyses using maximum- likelihood and Bayesian methods revealed interesting population structure and identified the existence of two distinct D. formosana lineages (Southeast Asian and Taiwanese populations). Genetic differentiation among groups of D. immigrans suggests the possibility of endemic speciation in Taiwan. In contrast, D. melanogaster remained one extensively large population throughout East and Southeast Asia, including nearby islets. A molecular clock was used to estimate divergence times, which were compared with past geographical events to infer evolutionary scenarios. Our findings suggest that interglacial periods may have caused population isolation, thus enhancing population differentiation more strongly for some of the Drosophila species. The population structure of each Drosophila species in East and Southeast Asia has been influenced by past geographic events.

Seawater Incursion Events in a Cretaceous Paleo-lake Revealed by Specific Marine Biological Markers

Many large paleo-lakes in North China were formed after the Triassic Era. Seawater incursion events (SWIEs) in these lakes have been extensively discussed in the literature, yet lack reliable methodology and solid evidence, which are essential for reconstructing and confirming SWIEs. The present study employs specific marine biological markers (24-n-propyl and 24-isopropyl cholestanes) to trace SWIEs in a dated core taken from the Songliao Basin (SLB). Two SWIEs were identified. The first SWIE from 91.37 to 89.00 Ma, was continuous and variable but not strong, while the second SWIE from 84.72 to 83.72 Ma was episodic and strong. SWIEs caused high total organic carbon (TOC) and negative δ¹³C_org values in the sediments, which were interpreted as an indication of high productivity in the lake, due to the enhancement of nutrient supplies as well as high levels of aqueous CO₂, due to the mixing of alkaline seawater and acidic lake water. The SWIEs in SLB were controlled by regional tectonic activity and eustatic variation. Movement direction changes of the Izanagi/Kula Plate in 90 Ma and 84 Ma created faults and triggered SWIEs. A high sea level, from 90 to 84 Ma, also facilitated the occurrence of SWIEs in SLB.

Biological evidence supports an early and complex emergence of the Isthmus of Panama

The formation of the Isthmus of Panama, which linked North and South America, is key to understanding the biodiversity, oceanography, atmosphere, and climate in the region. Despite its importance across multiple disciplines, the timing of formation and emergence of the Isthmus and the biological patterns it created have been controversial. Here, we analyze molecular and fossil data, including terrestrial and marine organisms, to show that biotic migrations across the Isthmus of Panama began several million years earlier than commonly assumed. An earlier evolution of the Isthmus has broad implications for the mechanisms driving global climate (e.g., Pleistocene glaciations, thermohaline circulation) as well as the rich biodiversity of the Americas.

The linking of North and South America by the Isthmus of Panama had major impacts on global climate, oceanic and atmospheric currents, and biodiversity, yet the timing of this critical event remains contentious. The Isthmus is traditionally understood to have fully closed by ca. 3.5 million years ago (Ma), and this date has been used as a benchmark for oceanographic, climatic, and evolutionary research, but recent evidence suggests a more complex geological formation. Here, we analyze both molecular and fossil data to evaluate the tempo of biotic exchange across the Americas in light of geological evidence. We demonstrate significant waves of dispersal of terrestrial organisms at approximately ca. 20 and 6 Ma and corresponding events separating marine organisms in the Atlantic and Pacific oceans at ca. 23 and 7 Ma. The direction of dispersal and their rates were symmetrical until the last ca. 6 Ma, when northern migration of South American lineages increased significantly. Variability among taxa in their timing of dispersal or vicariance across the Isthmus is not explained by the ecological factors tested in these analyses, including biome type, dispersal ability, and elevation preference. Migration was therefore not generally regulated by intrinsic traits but more likely reflects the presence of emergent terrain several millions of years earlier than commonly assumed. These results indicate that the dramatic biotic turnover associated with the Great American Biotic Interchange was a long and complex process that began as early as the Oligocene–Miocene transition.

Sea level change. Inherited landscapes and sea level change

Enabled by recently gained understanding of deep-seated and surficial Earth processes, a convergence of views between geophysics and sedimentary geology has been quietly taking place over the past several decades. Surface topography resulting from lithospheric memory, retained at various temporal and spatial scales, has become the connective link between these two methodologically diverse geoscience disciplines. Ideas leading to the hypothesis of plate tectonics originated largely with an oceanic focus, where dynamic and mostly horizontal movements of the crust could be envisioned. But when these notions were applied to the landscapes of the supposedly rigid plate interiors, there was less success in explaining the observed anomalies in terrestrial topography. Solid-Earth geophysics has now reached a developmental stage where vertical movements can be measured and modeled at meaningful scales and the deep-seated structures can be imaged with increasing resolution. Concurrently, there have been advances in quantifying mechanical properties of the lithosphere (the solid outer skin of Earth, usually defined to include both the crust and the solid but elastic upper mantle above the asthenosphere). The lithosphere acts as the intermediary that transfers the effects of mantle dynamics to the surface. These developments have allowed us to better understand the previously puzzling topographic features of plate interiors and continental margins. On the sedimentary geology side, new quantitative modeling techniques and holistic approaches to integrating source-to-sink sedimentary systems have led to clearer understanding of basin evolution and sediment budgets that allow the reconstruction of missing sedimentary records and past geological landscapes.

Plate tectonic controls on atmospheric CO2 levels since the Triassic

Climate trends on timescales of 10s to 100s of millions of years are controlled by changes in solar luminosity, continent distribution, and atmosphere composition. Plate tectonics affect geography, but also atmosphere composition through volcanic degassing of CO2 at subduction zones and midocean ridges. So far, such degassing estimates were based on reconstructions of ocean floor production for the last 150 My and indirectly, through sea level inversion before 150 My. Here we quantitatively estimate CO2 degassing by reconstructing lithosphere subduction evolution, using recent advances in combining global plate reconstructions and present-day structure of the mantle. First, we estimate that since the Triassic (250-200 My) until the present, the total paleosubduction-zone length reached up to ∼200% of the present-day value. Comparing our subduction-zone lengths with previously reconstructed ocean-crust production rates over the past 140 My suggests average global subduction rates have been constant, ∼6 cm/y: Higher ocean-crust production is associated with longer total subduction length. We compute a strontium isotope record based on subduction-zone length, which agrees well with geological records supporting the validity of our approach: The total subduction-zone length is proportional to the summed arc and ridge volcanic CO2 production and thereby to global volcanic degassing at plate boundaries. We therefore use our degassing curve as input for the GEOCARBSULF model to estimate atmospheric CO2 levels since the Triassic. Our calculated CO2 levels for the mid Mesozoic differ from previous modeling results and are more consistent with available proxy data.

Phylogeography and biogeography of the lower Central American Neotropics: diversification between two continents and between two seas

Lower Central America (LCA) provides a geologically complex and dynamic, richly biodiverse model for studying the recent assembly and diversification of a Neotropical biota. Here, we review the growing literature of LCA phylogeography studies and their contribution to understanding the origins, assembly, and diversification of the LCA biota against the backdrop of regional geologic and climatic history, and previous biogeographical inquiry. Studies to date reveal that phylogeographical signal within taxa of differing distributions reflects a diversity of patterns and processes rivalling the complexities of LCA landscapes themselves. Even so, phylogeography is providing novel insights into regional diversification (e.g. cryptic lineage divergences), and general evolutionary patterns are emerging. Congruent multi-taxon phylogeographic breaks are found across the Nicaraguan depression, Chorotega volcanic front, western and central Panama, and the Darién isthmus, indicating that a potentially shared history of responses to regional-scale (e.g. geological) processes has shaped the genetic diversity of LCA communities. By contrast, other species show unique demographic histories in response to overriding historical events, including no phylogeographic structure at all. These low-structure or incongruent patterns provide some evidence for a role of local, ecological factors (e.g. long-distance dispersal and gene flow in plants and bats) in shaping LCA communities. Temporally, comparative phylogeographical structuring reflects Pliocene-Pleistocene dispersal and vicariance events consistent with the timeline of emergence of the LCA isthmus and its major physiographic features, e.g. cordilleras. We emphasise the need to improve biogeographic inferences in LCA through in-depth comparative phylogeography projects capitalising on the latest statistical phylogeographical methods. While meeting the challenges of reconstructing the biogeographical history of this complex region, phylogeographers should also take up the critical service to society of applying their work to the conservation of its fascinating biodiversity.

Earth's punctuated tectonic evolution: cause and effect

Peaks in the Precambrian preserved crustal record are associated with major volcanic, tectonic and climatic events. These include addition of juvenile continental crust, voluminous high-temperature volcanism, massive mantle depletion, widespread orogeny and mineralization, large apparent polar wander velocity spikes, and subsequent palaeomagnetic intensity increases. These events impinge on the glaciation record, atmospheric and ocean chemistry, and on the rise of oxygen. Here we summarize and assess a number of geodynamic models that have been proposed to explain the observed episodicity in the Precambrian record. We find that episodic behaviour from nonlinear slab-driven models best explains the observed record. Examples of such slab-driven systems include mantle avalanches or episodic subduction. In these cases, rapid descent of slabs into the mantle drives fast plate motions and convergence at the surface. This is accompanied by large-scale upwellings of deep hot mantle, which contribute to voluminous volcanism. Further modelling will determine the relative importance of each mechanism, and reinforce the fundamental contribution of the mantle to the evolution of Earth's surface systems.

The Early Eocene equable climate problem: can perturbations of climate model parameters identify possible solutions?

Geological data for the Early Eocene (56-47.8 Ma) indicate extensive global warming, with very warm temperatures at both poles. However, despite numerous attempts to simulate this warmth, there are remarkable data-model differences in the prediction of these polar surface temperatures, resulting in the so-called 'equable climate problem'. In this paper, for the first time an ensemble with a perturbed climate-sensitive model parameters approach has been applied to modelling the Early Eocene climate. We performed more than 100 simulations with perturbed physics parameters, and identified two simulations that have an optimal fit with the proxy data. We have simulated the warmth of the Early Eocene at 560 ppmv CO2, which is a much lower CO2 level than many other models. We investigate the changes in atmospheric circulation, cloud properties and ocean circulation that are common to these simulations and how they differ from the remaining simulations in order to understand what mechanisms contribute to the polar warming. The parameter set from one of the optimal Early Eocene simulations also produces a favourable fit for the last glacial maximum boundary climate and outperforms the control parameter set for the present day. Although this does not 'prove' that this model is correct, it is very encouraging that there is a parameter set that creates a climate model able to simulate well very different palaeoclimates and the present-day climate. Interestingly, to achieve the great warmth of the Early Eocene this version of the model does not have a strong future climate change Charney climate sensitivity. It produces a Charney climate sensitivity of 2.7(°)C, whereas the mean value of the 18 models in the IPCC Fourth Assessment Report (AR4) is 3.26(°)C±0.69(°)C. Thus, this value is within the range and below the mean of the models included in the AR4.

Evidence for a rapid release of carbon at the Paleocene-Eocene thermal maximum

The Paleocene/Eocene thermal maximum (PETM) and associated carbon isotope excursion (CIE) are often touted as the best geologic analog for the current anthropogenic rise in pCO2. However, a causal mechanism for the PETM CIE remains unidentified because of large uncertainties in the duration of the CIE's onset. Here, we report on a sequence of rhythmic sedimentary couplets comprising the Paleocene/Eocene Marlboro Clay (Salisbury Embayment). These couplets have corresponding δ(18)O cycles that imply a climatic origin. Seasonal insolation is the only regular climate cycle that can plausibly account for δ(18)O amplitudes and layer counts. High-resolution stable isotope records show 3.5‰ δ(13)C decrease over 13 couplets defining the CIE onset, which requires a large, instantaneous release of (13)C-depleted carbon. During the CIE, a clear δ(13)C gradient developed on the shelf with the largest excursions in shallowest waters, indicating atmospheric δ(13)C decreased by ~20‰. Our observations and revised release rate are consistent with an atmospheric perturbation of 3,000-gigatons of carbon (GtC).