Early onset and tropical forcing of 100,000-year Pleistocene glacial cycles

Between 1.5 and 0.6 Myr ago, the period of the Earth's glacial cycles changed from 41 kyr, the period of the Earth's obliquity cycles, to 100 kyr, the period of the Earth's orbital eccentricity, which has a much smaller effect on global insolation. The timing of this transition and its causes pose one of the most perplexing problems in palaeoclimate research. Here we use complex demodulation to examine the phase evolution of precession and semiprecession cycles--the latter of which are phase-coupled to both precession and eccentricity--in the tropical and extratropical Atlantic Ocean. We find that about 1.5 Myr ago, tropical semiprecession cycles (with periods of about 11.5 kyr) started to propagate to higher latitudes, coincident with a growing amplitude envelope of the 100-kyr cycles. Evidence from numerical models suggests that cycles of about 10 kyr in length may be required to explain the high amplitude of the 100-kyr cycles. Combining our results with consideration of a modern analogue, we conclude that increased heat flow across the equator or from the tropics to higher latitudes around 1.5 Myr ago strengthened the semiprecession cycle in the Northern Hemisphere, and triggered the transition to sustained 100-kyr glacial cycles.     

Astronomical pacing of late Palaeocene to early Eocene global warming events

At the boundary between the Palaeocene and Eocene epochs, about 55 million years ago, the Earth experienced a strong global warming event, the Palaeocene-Eocene thermal maximum. The leading hypothesis to explain the extreme greenhouse conditions prevalent during this period is the dissociation of 1,400 to 2,800 gigatonnes of methane from ocean clathrates, resulting in a large negative carbon isotope excursion and severe carbonate dissolution in marine sediments. Possible triggering mechanisms for this event include crossing a threshold temperature as the Earth warmed gradually, comet impact, explosive volcanism or ocean current reorganization and erosion at continental slopes, whereas orbital forcing has been excluded. Here we report a distinct carbonate-poor red clay layer in deep-sea cores from Walvis ridge, which we term the Elmo horizon. Using orbital tuning, we estimate deposition of the Elmo horizon at about 2 million years after the Palaeocene-Eocene thermal maximum. The Elmo horizon has similar geochemical and biotic characteristics as the Palaeocene-Eocene thermal maximum, but of smaller magnitude. It is coincident with carbon isotope depletion events in other ocean basins, suggesting that it represents a second global thermal maximum. We show that both events correspond to maxima in the approximately 405-kyr and approximately 100-kyr eccentricity cycles that post-date prolonged minima in the 2.25-Myr eccentricity cycle, implying that they are indeed astronomically paced.     

Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary

Between 34 and 15 million years (Myr) ago, when planetary temperatures were 3-4 degrees C warmer than at present and atmospheric CO2 concentrations were twice as high as today, the Antarctic ice sheets may have been unstable. Oxygen isotope records from deep-sea sediment cores suggest that during this time fluctuations in global temperatures and high-latitude continental ice volumes were influenced by orbital cycles. But it has hitherto not been possible to calibrate the inferred changes in ice volume with direct evidence for oscillations of the Antarctic ice sheets. Here we present sediment data from shallow marine cores in the western Ross Sea that exhibit well dated cyclic variations, and which link the extent of the East Antarctic ice sheet directly to orbital cycles during the Oligocene/Miocene transition (24.1-23.7 Myr ago). Three rapidly deposited glacimarine sequences are constrained to a period of less than 450 kyr by our age model, suggesting that orbital influences at the frequencies of obliquity (40 kyr) and eccentricity (125 kyr) controlled the oscillations of the ice margin at that time. An erosional hiatus covering 250 kyr provides direct evidence for a major episode of global cooling and ice-sheet expansion about 23.7 Myr ago, which had previously been inferred from oxygen isotope data (Mi1 event).     

High-latitude influence on the eastern equatorial Pacific climate in the early Pleistocene epoch

Many records of tropical sea surface temperature and marine productivity exhibit cycles of 23 kyr (orbital precession) and 100 kyr during the past 0.5 Myr (refs 1-5), whereas high-latitude sea surface temperature records display much more pronounced obliquity cycles at a period of about 41 kyr (ref. 6). Little is known, however, about tropical climate variability before the mid-Pleistocene transition about 900 kyr ago, which marks the change from a climate dominated by 41-kyr cycles (when ice-age cycles and high-latitude sea surface temperature variations were dictated by changes in the Earth's obliquity) to the more recent 100-kyr cycles of ice ages. Here we analyse alkenones from marine sediments in the eastern equatorial Pacific Ocean to reconstruct sea surface temperatures and marine productivity over the past 1.8 Myr. We find that both records are dominated by the 41-kyr obliquity cycles between 1.8 and 1.2 Myr ago, with a relatively small contribution from orbital precession, and that early Pleistocene sea surface temperatures varied in the opposite sense to local annual insolation in the eastern equatorial Pacific Ocean. We conclude that during the early Pleistocene epoch, climate variability at our study site must have been determined by high-latitude processes that were driven by orbital obliquity forcing.     

Obliquity pacing of the late Pleistocene glacial terminations

The 100,000-year timescale in the glacial/interglacial cycles of the late Pleistocene epoch (the past approximately 700,000 years) is commonly attributed to control by variations in the Earth's orbit. This hypothesis has inspired models that depend on the Earth's obliquity (approximately 40,000 yr; approximately 40 kyr), orbital eccentricity (approximately 100 kyr) and precessional (approximately 20 kyr) fluctuations, with the emphasis usually on eccentricity and precessional forcing. According to a contrasting hypothesis, the glacial cycles arise primarily because of random internal climate variability. Taking these two perspectives together, there are currently more than thirty different models of the seven late-Pleistocene glacial cycles. Here we present a statistical test of the orbital forcing hypothesis, focusing on the rapid deglaciation events known as terminations. According to our analysis, the null hypothesis that glacial terminations are independent of obliquity can be rejected at the 5% significance level, whereas the corresponding null hypotheses for eccentricity and precession cannot be rejected. The simplest inference consistent with the test results is that the ice sheets terminated every second or third obliquity cycle at times of high obliquity, similar to the original proposal by Milankovitch. We also present simple stochastic and deterministic models that describe the timing of the late-Pleistocene glacial terminations purely in terms of obliquity forcing.     

Impacts of orbital forcing and atmospheric carbon dioxide on Miocene ice-sheet expansion

The processes causing the middle Miocene global cooling, which marked the Earth's final transition into an 'icehouse' climate about 13.9 million years ago (Myr ago), remain enigmatic. Tectonically driven circulation changes and variations in atmospheric carbon dioxide levels have been suggested as driving mechanisms, but the lack of adequately preserved sedimentary successions has made rigorous testing of these hypotheses difficult. Here we present high-resolution climate proxy records, covering the period from 14.7 to 12.7 million years ago, from two complete sediment cores from the northwest and southeast subtropical Pacific Ocean. Using new chronologies through the correlation to the latest orbital model, we find relatively constant, low summer insolation over Antarctica coincident with declining atmospheric carbon dioxide levels at the time of Antarctic ice-sheet expansion and global cooling, suggesting a causal link. We surmise that the thermal isolation of Antarctica played a role in providing sustained long-term climatic boundary conditions propitious for ice-sheet formation. Our data document that Antarctic glaciation was rapid, taking place within two obliquity cycles, and coincided with a striking transition from obliquity to eccentricity as the drivers of climatic change.     

A complete insect from the Late Devonian period

After terrestrialization, the diversification of arthropods and vertebrates is thought to have occurred in two distinct phases, the first between the Silurian and the Frasnian stages (Late Devonian period) (425-385?million years (Myr) ago), and the second characterized by the emergence of numerous new major taxa, during the Late Carboniferous period (after 345?Myr ago). These two diversification periods bracket the depauperate vertebrate Romer's gap (360-345?Myr ago) and arthropod gap (385-325?Myr ago), which could be due to preservational artefact. Although a recent molecular dating has given an age of 390?Myr for the Holometabola, the record of hexapods during the Early-Middle Devonian (411.5-391?Myr ago, Pragian to Givetian stages) is exceptionally sparse and based on fragmentary remains, which hinders the timing of this diversification. Indeed, although Devonian Archaeognatha are problematic, the Pragian of Scotland has given some Collembola and the incomplete insect Rhyniognatha, with its diagnostic dicondylic, metapterygotan mandibles. The oldest, definitively winged insects are from the Serpukhovian stage (latest Early Carboniferous period). Here we report the first complete Late Devonian insect, which was probably a terrestrial species. Its 'orthopteroid' mandibles are of an omnivorous type, clearly not modified for a solely carnivorous diet. This discovery narrows the 45-Myr gap in the fossil record of Hexapoda, and demonstrates further a first Devonian phase of diversification for the Hexapoda, as in vertebrates, and suggests that the Pterygota diversified before and during Romer's gap.     

Mantle thermal pulses below the Mid-Atlantic Ridge and temporal variations in the formation of oceanic lithosphere

A 20-Myr record of creation of oceanic lithosphere is exposed along a segment of the central Mid-Atlantic Ridge on an uplifted sliver of lithosphere. The degree of melting of the mantle that is upwelling below the ridge, estimated from the chemistry of the exposed mantle rocks, as well as crustal thickness inferred from gravity measurements, show oscillations of approximately 3-4 Myr superimposed on a longer-term steady increase with time. The time lag between oscillations of mantle melting and crustal thickness indicates that the mantle is upwelling at an average rate of approximately 25 mm x yr(-1), but this appears to vary through time. Slow-spreading lithosphere seems to form through dynamic pulses of mantle upwelling and melting, leading not only to along-axis segmentation but also to across-axis structural variability. Also, the central Mid-Atlantic Ridge appears to have become steadily hotter over the past 20 Myr, possibly owing to north-south mantle flow.     

Atlantic overturning responses to Late Pleistocene climate forcings

The factors driving glacial changes in ocean overturning circulation are not well understood. On the basis of a comparison of 20 climate variables over the past four glacial cycles, the SPECMAP project proposed that summer insolation at high northern latitudes (that is, Milankovitch forcing) drives the same sequence of ocean circulation and other climate responses over 100-kyr eccentricity cycles, 41-kyr obliquity cycles and 23-kyr precession cycles. SPECMAP analysed the circulation response at only a few sites in the Atlantic Ocean, however, and the phase of circulation response has been shown to vary by site and orbital band. Here we test the SPECMAP hypothesis by measuring the phase of orbital responses in benthic delta(13)C (a proxy indicator of ocean nutrient content) at 24 sites throughout the Atlantic over the past 425 kyr. On the basis of delta(13)C responses at 3,000-4,010 m water depth, we find that maxima in Milankovitch forcing are associated with greater mid-depth overturning in the obliquity band but less overturning in the precession band. This suggests that Atlantic overturning is strongly sensitive to factors beyond ice volume and summer insolation at high northern latitudes. A better understanding of these processes could lead to improvements in model estimates of overturning rates, which range from a 40 per cent increase to a 40 per cent decrease at the Last Glacial Maximum and a 10-50 per cent decrease over the next 140 yr in response to projected increases in atmospheric CO(2) (ref. 4).     

A 300-million-year record of atmospheric carbon dioxide from fossil plant cuticles

To understand better the link between atmospheric CO2 concentrations and climate over geological time, records of past CO2 are reconstructed from geochemical proxies. Although these records have provided us with a broad picture of CO2 variation throughout the Phanerozoic eon (the past 544 Myr), inconsistencies and gaps remain that still need to be resolved. Here I present a continuous 300-Myr record of stomatal abundance from fossil leaves of four genera of plants that are closely related to the present-day Ginkgo tree. Using the known relationship between leaf stomatal abundance and growing season CO2 concentrations, I reconstruct past atmospheric CO2 concentrations. For the past 300 Myr, only two intervals of low CO2 (<1,000 p.p.m.v.) are inferred, both of which coincide with known ice ages in Neogene (1-8 Myr) and early Permian (275-290 Myr) times. But for most of the Mesozoic era (65-250 Myr), CO2 levels were high (1,000-2,000 p.p.m.v.), with transient excursions to even higher CO2 (>2,000 p.p.m.v.) concentrations. These results are consistent with some reconstructions of past CO2 (refs 1, 2) and palaeotemperature records, but suggest that CO2 reconstructions based on carbon isotope proxies may be compromised by episodic outbursts of isotopically light methane. These results support the role of water vapour, methane and CO2 in greenhouse climate warming over the past 300 Myr.     

Eocene global warming events driven by ventilation of oceanic dissolved organic carbon

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (～65-34 million years (Myr) ago). The most extreme hyperthermal was the ～170 thousand year (kyr) interval of 5-7 °C global warming during the Palaeocene-Eocene Thermal Maximum (PETM, 56?Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs, and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon. Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (～40?kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was re-sequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM. Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources, but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.     

Using the fossil record to estimate the age of the last common ancestor of extant primates

Divergence times estimated from molecular data often considerably predate the earliest known fossil representatives of the groups studied. For the order Primates, molecular data calibrated with various external fossil dates uniformly suggest a mid-Cretaceous divergence from other placental mammals, some 90 million years (Myr) ago, whereas the oldest known fossil primates are from the basal Eocene epoch (54-55 Myr ago). The common ancestor of primates should be earlier than the oldest known fossils, but adequate quantification is needed to interpret possible discrepancies between molecular and palaeontological estimates. Here we present a new statistical method, based on an estimate of species preservation derived from a model of the diversification pattern, that suggests a Cretaceous last common ancestor of primates, approximately 81.5 Myr ago, close to the initial divergence time inferred from molecular data. It also suggests that no more than 7% of all primate species that have ever existed are known from fossils. The approach unites all the available palaeontological methods of timing evolutionary events: the fossil record, extant species and clade diversification models.     

Warm tropical ocean surface and global anoxia during the mid-Cretaceous period

The middle of the Cretaceous period (about 120 to 80 Myr ago) was a time of unusually warm polar temperatures, repeated reef-drowning in the tropics and a series of oceanic anoxic events (OAEs) that promoted both the widespread deposition of organic-carbon-rich marine sediments and high biological turnover. The cause of the warm temperatures is unproven but widely attributed to high levels of atmospheric greenhouse gases such as carbon dioxide. In contrast, there is no consensus on the climatic causes and effects of the OAEs, with both high biological productivity and ocean 'stagnation' being invoked as the cause of ocean anoxia. Here we show, using stable isotope records from multiple species of well-preserved foraminifera, that the thermal structure of surface waters in the western tropical Atlantic Ocean underwent pronounced variability about 100 Myr ago, with maximum sea surface temperatures 3-5 degrees C warmer than today. This variability culminated in a collapse of upper-ocean stratification during OAE-1d (the 'Breistroffer' event), a globally significant period of organic-carbon burial that we show to have fundamental, stratigraphically valuable, geochemical similarities to the main OAEs of the Mesozoic era. Our records are consistent with greenhouse forcing being responsible for the warm temperatures, but are inconsistent both with explanations for OAEs based on ocean stagnation, and with the traditional view (reviewed in ref. 12) that past warm periods were more stable than today's climate.     

Cooler winters as a possible cause of mass extinctions at the Eocene/Oligocene boundary

The Eocene/Oligocene boundary, at about 33.7 Myr ago, marks one of the largest extinctions of marine invertebrates in the Cenozoic period. For example, turnover of mollusc species in the US Gulf coastal plain was over 90% at this time. A temperature change across this boundary--from warm Eocene climates to cooler conditions in the Oligocene--has been suggested as a cause of this extinction event, but climate reconstructions have not provided support for this hypothesis. Here we report stable oxygen isotope measurements of aragonite in fish otoliths--ear stones--collected across the Eocene/Oligocene boundary. Palaeo-temperatures reconstructed from mean otolith oxygen isotope values show little change through this interval, in agreement with previous studies. From incremental microsampling of otoliths, however, we can resolve the seasonal variation in temperature, recorded as the otoliths continue to accrete new material over the life of the fish. These seasonal data suggest that winters became about 4 degrees C colder across the Eocene/Oligocene boundary. We suggest that temperature variability, rather than change in mean annual temperature, helped to cause faunal turnover during this transition.     

Combined obliquity and precession pacing of late Pleistocene deglaciations

Milankovitch proposed that Earth resides in an interglacial state when its spin axis both tilts to a high obliquity and precesses to align the Northern Hemisphere summer with Earth's nearest approach to the Sun. This general concept has been elaborated into hypotheses that precession, obliquity or combinations of both could pace deglaciations during the late Pleistocene. Earlier tests have shown that obliquity paces the late Pleistocene glacial cycles but have been inconclusive with regard to precession, whose shorter period of about 20,000 years makes phasing more sensitive to timing errors. No quantitative test has provided firm evidence for a dual effect. Here I show that both obliquity and precession pace late Pleistocene glacial cycles. Deficiencies in time control that have long stymied efforts to establish orbital effects on deglaciation are overcome using a new statistical test that focuses on maxima in orbital forcing. The results are fully consistent with Milankovitch's proposal but also admit the possibility that long Southern Hemisphere summers contribute to deglaciation.     

Ecosystem remodelling among vertebrates at the Permian-Triassic boundary in Russia

The mass extinction at the Permian-Triassic boundary, 251 million years (Myr) ago, is accepted as the most profound loss of life on record. Global data compilations indicate a loss of 50% of families or more, both in the sea and on land, and these figures scale to a loss of 80-96% of species, based on rarefaction analyses. This level of loss is confirmed by local and regional-scale studies of marine sections, but the terrestrial record has been harder to analyse in such close detail. Here we document the nature of the event in Russia in a comprehensive survey of 675 specimens of amphibians and reptiles from 289 localities spanning 13 successive geological time zones in the South Urals basin. These changes in diversity and turnover cannot be explained simply by sampling effects. There was a profound loss of genera and families, and simplification of ecosystems, with the loss of small fish-eaters and insect-eaters, medium and large herbivores and large carnivores. Faunal dynamics also changed, from high rates of turnover through the Late Permian period to greater stability at low diversity through the Early Triassic period. Even after 15 Myr of ecosystem rebuilding, some guilds were apparently still absent-small fish-eaters, small insect-eaters, large herbivores and top carnivores.     

珠江口盆地中中新世韩江组天文调谐地质年代表

利用多种频谱分析方法,对珠江口盆地Langhian阶韩江组三角洲相硅质碎屑岩及浅海相碳酸盐岩韵律性交替地层进行旋回地层学分析,然后选择最新天文周期解决方案中的北纬65°夏半年日照量曲线及斜率曲线作为目标曲线,对自然伽马测井曲线进行天文调谐,得到高分辨率天文调谐地质年代表。结果表明:珠江口盆地中中新世韩江组地层旋回受天文轨道周期变化的驱动;利用天文调谐地质年代表进行的定年结果与前人利用古地磁对南中国海相邻盆地的定年一致;通过天文调谐地质年代表求得的研究层段高分辨率沉积速率变化与天文轨道周期偏心率曲线可以对应,推测沉积速率变化的主要控制因素是405 ka长周期偏心率。     

基于米兰柯维奇理论的地层划分新方法

提出了以米兰柯维奇理论为基础的地层划分新方法 ,该方法利用频谱中峰值频率的位置进行旋回识别。定义了岁差和斜率的旋回发育状况函数 ,利用该函数来反映地层中米兰柯维奇旋回的发育程度 ,并据此分析了地层的连续性、判断出沉积间断的位置及进行了地层划分。研究结果表明 ,利用地层中存在的米兰柯维奇旋回可以进行精细地层划分 ,划分结果可进行远距离对比。该方法在大港油田、胜利油田、塔河油田的实际应用中取得了较好的效果 ,具有准确性高、效率高、应用范围广的特点。     

Iron meteorite evidence for early formation and catastrophic disruption of protoplanets

In our Solar System, the planets formed by collisional growth from smaller bodies. Planetesimals collided to form Moon-to-Mars-sized protoplanets in the inner Solar System in 0.1-1 Myr, and these collided more energetically to form planets. Insights into the timing and nature of collisions during planetary accretion can be gained from meteorite studies. In particular, iron meteorites offer the best constraints on early stages of planetary accretion because most are remnants of the oldest bodies, which accreted and melted in <1.5 Myr, forming silicate mantles and iron-nickel metallic cores. Cooling rates for various groups of iron meteorites suggest that if the irons cooled isothermally in the cores of differentiated bodies, as conventionally assumed, these bodies were 5-200 km in diameter. This picture is incompatible, however, with the diverse cooling rates observed within certain groups, most notably the IVA group, but the large uncertainties associated with the measurements do not preclude it. Here we report cooling rates for group IVA iron meteorites that range from 100 to 6,000 K Myr(-1), increasing with decreasing bulk Ni. Improvements in the cooling rate model, smaller error bars, and new data from an independent cooling rate indicator show that the conventional interpretation is no longer viable. Our results require that the IVA meteorites cooled in a 300-km-diameter metallic body that lacked an insulating mantle. This body probably formed approximately 4,500 Myr ago in a 'hit-and-run' collision between Moon-to-Mars-sized protoplanets. This demonstrates that protoplanets of approximately 10(3) km size accreted within the first 1.5 Myr, as proposed by theory, and that fragments of these bodies survived as asteroids.     

Oxidation of the Ediacaran ocean

Oxygenation of the Earth's surface is increasingly thought to have occurred in two steps. The first step, which occurred approximately 2,300 million years (Myr) ago, involved a significant increase in atmospheric oxygen concentrations and oxygenation of the surface ocean. A further increase in atmospheric oxygen appears to have taken place during the late Neoproterozoic period ( approximately 800-542 Myr ago). This increase may have stimulated the evolution of macroscopic multicellular animals and the subsequent radiation of calcified invertebrates, and may have led to oxygenation of the deep ocean. However, the nature and timing of Neoproterozoic oxidation remain uncertain. Here we present high-resolution carbon isotope and sulphur isotope records from the Huqf Supergroup, Sultanate of Oman, that cover most of the Ediacaran period (approximately 635 to approximately 548 Myr ago). These records indicate that the ocean became increasingly oxygenated after the end of the Marinoan glaciation, and they allow us to identify three distinct stages of oxidation. When considered in the context of other records from this period, our data indicate that certain groups of eukaryotic organisms appeared and diversified during the second and third stages of oxygenation. The second stage corresponds with the Shuram excursion in the carbon isotope record and seems to have involved the oxidation of a large reservoir of organic carbon suspended in the deep ocean, indicating that this event may have had a key role in the evolution of eukaryotic organisms. Our data thus provide new insights into the oxygenation of the Ediacaran ocean and the stepwise restructuring of the carbon and sulphur cycles that occurred during this significant period of Earth's history.