Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

The evolutionary history of size variation of planktic foraminiferal assemblages in the Cenozoic

The iterative evolutionary radiation of planktic foraminifers is a well-documented macroevolutionary process. Here we document the accompanying size changes in entire planktic foraminiferal assemblages for the past 70 My and their relationship to paleoenvironmental changes. After the size decrease at the Cretaceous/Paleogene (K/P) boundary, high latitude assemblages remained consistently small. Size evolution in low latitudes can be divided into three major phases: the first is characterized by dwarfs (65-42 Ma), the second shows moderate size fluctuations (42-14 Ma), and in the third phase, planktic foraminifers have grown to the unprecedented sizes observed today. Our analyses of size variability with paleoproxy records indicate that periods of size increase coincided with phases of global cooling (Eocene and Neogene). These periods were characterized by enhanced latitudinal and vertical temperature gradients in the oceans and high diversity (polytaxy). In the Paleocene and during the Oligocene, the observed (minor) size changes of the largely low-diversity (oligotaxic) assemblages seem to correlate with productivity changes. However, polytaxy per se was not responsible for larger test sizes.

Radiolarian assemblages of the late Neogene North Pacific

Radiolarian census and abundance data were collected from three deep-sea cores drilled by the Ocean Drilling Program Sites 884, 887 and 1151 to investigate patterns of ecologic changes in space and time during the last 16 million years for the mid-latitude to subarctic North Pacific. High concentrations of radiolarians occurred between 9.0 and 2.7 Ma. Radiolarian species richness was highest in the early middle Miocene at each site and gradually decreased up to about 7 Ma, coinciding with a well-established global cooling trend. A degree of overlap index calculated for radiolarian assemblages revealed 11 faunal change events, of which 8 corresponded to global cooling events and expansions of polar ice sheets. Three of the faunal change events were observed within the peak of radiolarian accumulation rate and were ascribed to changes in primary productivity in the North Pacific rather than global climatic changes. Our assemblage analyses revealed that north-south differentiation in radiolarian assemblages in the northwestern Pacific has existed since 16 Ma and became more distinct via major steps at 6.8 Ma and 2.7 Ma, coinciding with major glaciation events, and that east-west faunal contrasts in the subarctic region became obvious beginning at 11.7 Ma and changed to a different mode around 6.8 Ma. The observed east-west faunal differences possibly reflect east to west climate differences that were characterized by cooler temperatures in the east than the west during the late Miocene (11.7-6.8 Ma) and then by the opposite temperature trend (6.8 Ma-Recent). A severe glaciation at 2.7 Ma played a large role, particularly in temporal changes in radiolarian accumulation rate and assemblage composition.

Occurrence and isotopic composition (d18O, d13C) of benthic foraminifera in sediments from DSDP Leg 72

Changes in the vertical water mass structure of the Vema Channel during the Pliocene have been inferred from benthic foraminiferal assemblages and stable isotopic analyses from three sites of DSDP Leg 72 (South Atlantic). Faunal and isotopic results from Sites 516A and 518 suggest that a major change occurred in deep-water circulation patterns in the late Pliocene near 3.2 Ma. Benthic oxygen isotopic records from Sites 516A and 518 show a characteristic increase in d18O values near 3.2 Ma. This has been documented in numerous Pliocene isotopic records. The magnitude of the oxygen isotopic enrichment near 3.2 Ma appears to increase with water depth from an average enrichment of 0.34 per mil in Site 516A (1313 m) to an average enrichment of 0.58 per mil in Site 518 (3944 m). We suggest that this enrichment resulted partly from a change in deep-water circulation patterns which included a decrease in bottom-water temperatures. Planktonic d18O values near 3.2 Ma show no evidence of an enrichment which would be indicative of an increase in global ice volume. On the contrary, d18O values in Sites 517 and 518 become more depleted near 3.2 Ma, indicating a surface-water warming perhaps due to a change in the strength and/or position of the Brazil Current. An increase in the relative abundance of the benthic foraminifer Nuttalides umbonifera, which is associated with Antarctic Bottom Water (AABW) in the modern ocean, coincides with the benthic 18O enrichment in Site 518. At 3.2 Ma, oxygen and carbon isotopic gradients between Sites 518 (3944 m) and 516A (1313 m) show a marked increase such that Site 518 becomes enriched in 18O and depleted in 13C relative to Site 516A. This enrichment in d18O is interpreted as partly representing a temperature decrease at Site 518; the depletion in d13C indicates a corrosive water mass which is high in metabolic CO2. We suggest that benthic foraminiferal and stable isotopic changes in Site 518 resulted from a pulse-like increase in the formation of AABW near 3.2 Ma. The cause of this circulation event may have been linked to global cooling and/or the final closure of the Central American Seaway.

Stable isotope record, opal accumulation rate and geochemistry of middle Miocene sediments of IODP Site 321-U1338

During the middle Miocene, Earth's climate transitioned from a relatively warm phase (Miocene climatic optimum) into a colder mode with re-establishment of permanent ice sheets on Antarctica, thus marking a fundamental step in Cenozoic cooling. Carbon sequestration and atmospheric CO2 drawdown through increased terrestrial and/or marine productivity have been proposed as the main drivers of this fundamental transition. We integrate high-resolution (1-3 k.y.) benthic stable isotope data with XRF-scanner derived biogenic silica and carbonate accumulation estimates in an exceptionally well-preserved sedimentary archive, recovered at Integrated Ocean Drilling Program Site U1338, to reconstruct eastern equatorial Pacific productivity variations and to investigate temporal linkages between high- and low-latitude climate change over the interval 16-13 Ma. Our records show that the climatic optimum (16.8-14.7 Ma) was characterized by high amplitude climate variations, marked by intense perturbations of the carbon cycle. Episodes of peak warmth at (southern hemisphere) insolation maxima coincided with transient shoaling of the carbonate compensation depth and enhanced carbonate dissolution in the deep ocean. A switch to obliquity-paced climate variability after 14.7 Ma concurred with a general improvement in carbonate preservation and the onset of stepwise global cooling, culminating with extensive ice growth over Antarctica at ~13.8 Ma. We find that two massive increases in opal accumulation at ~14.0 and ~13.8 Ma occurred just before and during the final and most prominent cooling step, supporting the hypothesis that enhanced siliceous productivity in the eastern equatorial Pacific contributed to CO2 drawdown.

Astronomically modulated Neogene sediment records of ODP Site 184-1148

General global cooling over the Neogene has been modulated by changes in Earth's orbital parameters. Investigations of deep-sea sediment sequences show that various orbital cycles can dominate climate records for different latitudes or for different time intervals. However, a comprehensive understanding of astronomical imprints over the entire Neogene has been elusive because of the general absence of long, continuous records extending beyond the Pliocene. We present benthic foraminiferal d18O and d13C records over the past 23 Ma at Ocean Drilling Program Site 1148 in the northern South China Sea and construct an astronomically tuned timescale (TJ08) for these records based on natural gamma radiation and color reflectance data at this site. Our results show that a 41 ka cycle has dominated sediment records at this location over the Neogene, displaying a linear response to orbital forcing. A 100 ka cycle has also been significant. However, it is correlated nonlinearly with Earth's orbital variations at the 100 ka band. The sediment records also display a prominent 405 ka cycle. Although this cycle was coherent with orbital forcing during the Oligocene and the early Miocene, it was not coherent with Earth's orbital variations at the 405 ka band over the whole Neogene. Amplification of Northern Hemisphere and Southern Hemisphere glaciation since the middle Miocene may be responsible for this change in sedimentary response. Our benthic foraminifera d18O and d13C records further exhibit amplitude variations with longer periods of 600, 1000, 1200, and 2400 ka. Apparently, these cycles are nonlinear responses to insolation forcing.

(Appendix) Analysis of benthic foraminifera from ODP Leg 181 and DSDP Leg 90

Canonical correspondence analysis indicates that the distribution of Neogene benthic foraminiferal faunas (>63 µm) in seven DSDP and ODP sites (500-4500 m water depth) east of New Zealand (38-51°S, 170°E-170°W) is most strongly influenced by depth (water mass stratification), and secondly by age (palaeoceanographic changes influencing faunal composition and biotic evolution). Stratigraphic faunal changes are interpretted in terms of the pulsed sequential development of southern, and later northern, polar glaciation and consequent cooling of bottom waters, increased vertical and lateral stratification of ocean water masses, and increased overall and seasonal surface water productivity. Oligocene initiation of the Antarctic Circumpolar Current and Deep Western Boundary Current (DWBC), flowing northwards past New Zealand, resulted in extensive hiatuses throughout the Southwest Pacific, some extending through into the Miocene. Planktic foraminiferal fragmentation index values indicate that carbonate dissolution was significant at abyssal depths throughout most of the Neogene, peaking at upper abyssal depths in the late Miocene (11-7 Ma), with the lysocline progressively deepened thereafter. Miocene abyssal faunas are dominated by Globocassidulina subglobosa and Oridorsalis umbonatus, with increasing Epistominella exigua after 16 Ma at upper abyssal depths. Peak abundances of Epistominella umbonifera indicate increased input of cold Southern Component Water to the DWBC at 7-6 Ma. Faunal association changes imply establishment of the modern Oxygen Minimum Zone (upper Circumpolar Deep Water) in the latest Miocene. Significant latitudinal differences between the benthic foraminiferal faunas at lower bathyal depths indicate the existence of an oceanic front along the Chatham Rise (location of present Subtropical Front), since the early late Miocene at least, with more pulsed productivity (higher E. exigua) along the south side. Modern Antarctic Intermediate Water faunal associations were established north of the Chatham Rise at 10-9 Ma, and south of it at 3-1.5 Ma. Middle-upper bathyal faunas on the Campbell Plateau are dominated by reticulate bolivinids during the early and middle Miocene, indicative of sustained productivity above relatively sluggish, suboxic bottom waters. Faunal changes and hiatuses indicate increased current vigour over the Campbell Plateau from the latest Miocene on. Surface water productivity (food supply) appears to have increased in three steps (at times of enhanced global cooling) marked by substantially increased relative abundance of: (1) Abditodentrix pseudothalmanni, Alabaminella weddellensis, Cassidulina norvangi (16-15 Ma, increased pulsed productivity); (2) Bulimina marginata f. aculeata, Nonionella auris, Trifarina angulosa, Uvigerina peregrina (3-1.5 Ma, increased overall productivity); and (3) Cassidulina carinata (1-0.5 Ma, increased overall productivity). Three intervals of deep-sea benthic foraminiferal taxonomic turnover are recognised (16-15, 11.5-10, 2-0.5 Ma) corresponding to intervals of enhanced global cooling and possible productivity changes. The late Pliocene-middle Pleistocene extinction, associated with increasing Northern Hemisphere glaciation, culminating in the middle Pleistocene climatic transition, was more significant in the study area than the earlier Neogene turnovers.

Geochemistry, lithology, and bulk density of Cricket Flat paleosol samples

The high-resolution marine isotope climate record indicates pronounced global cooling during the Langhian (16-13.8 Ma), beginning with the warm middle Miocene climatic optimum and ending with significant Antarctic ice sheet expansion and the transition to "icehouse" conditions. Terrestrial paleoclimate data from this interval is sparse and sometimes conflicting. In particular, there are gaps in the terrestrial record in the Pacific Northwest during the late Langhian and early Serravallian between about 14.5 and 12.5 Ma. New terrestrial paleoclimate data from this time and region could reconcile these conflicting records. Paleosols are particularly useful for reconstructing paleoenvironment because the rate and style of pedogenesis is primarily a function of surface environmental conditions; however, complete and well-preserved paleosols are uncommon. Most soils form in erosive environments that are not preserved, or in environments such as floodplains that accumulate in small increments; the resulting cumulic soils are usually thin, weakly developed, and subject to diagenetic overprinting from subsequent soils. The paleosol at Cricket Flat in northeastern Oregon is an unusually complete and well-preserved paleosol from a middle Miocene volcanic sequence in the Powder River Volcanic Field. An olivine basalt flow buried the paleosol at approximately 13.8 ± 0.6 Ma, based on three 40Ar/39Ar dates on the basalt. We described the Cricket Flat paleosol and used its physical and chemical profile and micromorphology to assess pedogenesis. The Cricket Flat paleosol is an Ultisol-like paleosol, chemically consistent with a high degree of weathering. Temperature and rainfall proxies suggest that Cricket Flat received 1120 ± 180 mm precipitation y-1 and experienced a mean annual temperature of 14.5 ± 2.1 °C during the formation of the paleosol, significantly warmer and wetter than today. This suggests slower cooling after the middle Miocene climatic optimum than is seen in the existing paleosol record.

Middle Miocene mineralogical events of Marion Plateau, northeast of Australia

Three ODP sites located on the Marion Plateau, Northeast Australian margin, were investigated for clay mineral and bulk mineralogy changes through the early to middle Miocene. Kaolinite to smectite (K/S) ratios, as well as mass accumulation rates of clays, point to a marked decrease in accumulation of smectite associated with an increase in accumulation of kaolinite starting at ~15.6 Ma, followed by a second increase in accumulation of kaolinite at ~13.2 Ma. Both of these increases are correlative to an increase in the calcite to detritus ratio. Comparison of our record with published precipitation proxies from continental Queensland indicates that increases in kaolinite did not correspond to more intense tropical-humid conditions, but instead to periods of greater aridity. Three mechanisms are explored to explain the temporal trends in clay on the Marion Plateau: sea-level changes, changes in oceanic currents, and denudation of the Australian continent followed by reworking and eolian transport of clays. Though low mass accumulation rates of kaolinite are compatible with a possible contribution of eolian material after 14 Ma, when Australia became more arid, the lateral distribution of kaolinite along slope indicates mainly fluvial input for all clays and thus rules out this mechanism as well as oceanic current transport as the main controls behind clay accumulation on the plateau. We propose a model explaining the good correlation between long-term sea-level fall, decrease in smectite accumulation, increase in kaolinite accumulation and increase in carbonate input to the distal slope locations. We hypothesize that during low sea level and thus periods of drier continental climate in Queensland, early Miocene kaolinite-rich lacustrine deposits were being reworked, and that the progradation of the heterozoan carbonate platforms towards the basin center favored input of carbonate to the distal slope sites. The major find of our study is that increase kaolinite fluxes on the Queensland margin during the early and middle Miocene did not reflect the establishment of a tropical climate, and this stresses that care must be taken when reconstructing Australian climate based on deep-sea clay records alone.

Oxygen and Carbon stable isotope ratios for Site U1334

The Oligocene-Miocene transition (OMT) (~23 Ma) is interpreted as a transient global cooling event, associated with a large-scale Antarctic ice sheet expansion. Here we present a 2.23 Myr long high-resolution (~3 kyr) benthic foraminiferal oxygen and carbon isotope (d18O and d13C) record from Integrated Ocean Drilling Program Site U1334 (eastern equatorial Pacific Ocean), covering the interval from 21.91 to 24.14 Ma. To date, five other high-resolution benthic foraminiferal stable isotope stratigraphies across this time interval have been published, showing a ~1 per mil increase in benthic foraminiferal d18O across the OMT. However, these records are still few and spatially limited and no clear understanding exists of the global versus local imprints. We show that trends and the amplitudes of change are similar at Site U1334 as in other high-resolution stable isotope records, suggesting that these represent global deep water signals. We create a benthic foraminiferal stable isotope stack across the OMT by combining Site U1334 with records from ODP Sites 926, 929, 1090, 1264, and 1218 to best approximate the global signal. We find that isotopic gradients between sites indicate interbasinal and intrabasinal variabilities in deep water masses and, in particular, note an offset between the equatorial Atlantic and the equatorial Pacific, suggesting that a distinct temperature gradient was present during the OMT between these deep water masses at low latitudes. A convergence in the d18O values between infaunal and epifaunal species occurs between 22.8 and 23.2 Ma, associated with the maximum d18O excursion at the OMT, suggesting climatic changes associated with the OMT had an effect on interspecies offsets of benthic foraminifera. Our data indicate a maximum glacioeustatic sea level change of ~50 m across the OMT.

Palynology and sea surface temperature reconstruction for the mid-Cretaceous Oceanic Anoxic Event 2 at ODP Hole 210-1276A

The mid-Cretaceous is thought to be a greenhouse world with significantly higher atmospheric pCO2 and sea-surface temperatures as well as a much flatter latitudinal thermal gradient compared to the present. This time interval was punctuated by the Cenomanian/Turonian Oceanic Anoxic Event (OAE-2, ~ 93.5 Myr ago), an episode of global, massive organic carbon burial that likely resulted in a large and abrupt pCO2 decline. However, the climatic consequences of this pCO2 drop are yet poorly constrained. We determined the first, high-resolution sea-surface temperature (SST) record across OAE-2 from a deep-marine sedimentary sequence at Ocean Drilling Program (ODP) Site 1276 in the mid-latitudinal Newfoundland Basin, NW Atlantic. By employing the organic palaeothermometer TEX86, we found that SSTs across the OAE-2 interval were extremely high, but were punctuated by a remarkably large cooling (5-11 °C), which is synchronous with the 2.5-5.5 °C cooling in SST records from equatorial Atlantic sites, and the "Plenus Cold Event". Because this global cooling event is concurrent with increased organic carbon burial, it likely acted in response to the associated pCO2 drop. Our findings imply a substantial increase in the latitudinal SST gradient in the proto-North Atlantic during this period of global cooling and reduced atmospheric pCO2, suggesting a strong coupling between pCO2 and latitudinal thermal gradients under greenhouse climate conditions.

Nd isotope data for ODP Leg 208 holes

The flow of deep-water masses is a key component of heat transport in the modern climate system, yet the role of deep-ocean heat transport during periods of extreme warmth is poorly understood. The present mode of meridional overturning circulation is characterized by deep-water formation in both the North Atlantic and the Southern Ocean. However, a different mode of meridional overturning circulation operated during the extreme greenhouse warmth of the early Cenozoic, during which time the Southern Ocean was the dominant region of deep-water formation. The combination of general global cooling and tectonic evolution of the Atlantic basins over the past ~55 m.y. ultimately led to the development of a mode of overturning circulation characterized by both Southern Ocean and North Atlantic deep-water sources. The change in deep-water circulation mode may, in turn, have affected global climate; however, unraveling the causes and consequences of this transition requires a better understanding of the timing of the transition. New Nd isotope data from the southeastern Atlantic Ocean indicate that the initial transition to a bipolar mode of deep-water circulation occurred in the early Oligocene, ca. 33 Ma. The likely cause of significant deep-water production in the North Atlantic was tectonic deepening of the sill separating the Greenland-Norwegian Sea from the North Atlantic.

Grain size analysis and sedimentation rates of ODP Hole 175-1085A off the Namibian coast

Sediments from the ODP Site 1085A were studied to investigate the impacts of global cooling in the Middle and Late Miocene on the climate in Southwestern Africa. The size composition of the sediment was analysed emphasising the silt fraction. A comparison with the modern grain size distribution and suitable transport processes made it possible to assign specific transport processes to the grain size composition. Three processes are considered for transport of terrigeneous silt: while there was no evidence found for (1) transport by ocean currents, the analyses showed signals of (2) wind transport indicating dry conditions associated with a cool climate and (3) fluvial transport that points to humid and warm conditions. Three climatic phases were defined. The first phase from 13.8 to 11.8 Myr reveals a stable humid climate in Southwest Africa independent of the Antarctic glaciations. During the second phase from 11.8 to 10.4 Myr the regional climate cooled considerably but was not drier. Additionally, the climate during this phase reacted to the Antarctic glaciations. This cooling-trend continued during phase 3 from 10.4 to 9.0 Myr with a significant increase in dust input, pointing to overall drier conditions. However, fluvial transport still remained as the main source.

Brachiopods in Early Carboniferous

Changes in Mississippian global paleogeography derived from the reconfiguration of the continents, a reversal in ocean currents and global cooling. Although the tectonic and climatic changes are well-documented, their effects on the distribution of brachiopod fauna are poorly documented. Here we present systematic quantitative analyses on global paleobiogeography based on a global brachiopod database from the Mississippian (i.e., Tournaisian, Visean, and Serpukhovian). The dataset consists of 2123 species of 344 brachiopod genera from 1156 localities. Our results reveal that global provincialism was not evident during the Tournaisian and Visean Stages. Two realms, i.e., the Gondwanan and Paleoequatorial Realms, are recognized during the Tournaisian. The Paleoequatorial Realm dominates during the Visean Stage, whereas the Gondwanan Realm is not documented due to the absence of data points. In contrast to the early and middle Mississippian stages, faunal provincialism is greatly enhanced in the Serpukhovian Stage with Paleotethyan and North American realms easily distinguished. This indicates that the Rheic Ocean was closed before the Serpukhovian due to the collision between Gondwana and Laurussia, that disrupted faunal interchange between the Paleotethys and North America. In addition, the paleolatitude-related thermal gradient was enhanced and the Boreal Realm was distinguished from the Paleotethyan Realm during the onset of the Late Palaeozoic Ice Age (LPIA) in the Serpukhovian. The paleolatitude diversity gradient pattern further shows a distinct shift of diversity center from the southern tropic zone in the Tournaisian and Visean to the northern tropic zone in the Serpukhovian.