Color data, carbon and oxygen isotopes, and offsets from BBCP Polecat Bench cores 2A and 2B

The Earth's climate abruptly warmed by 5-8 °C during the Palaeocene-Eocene thermal maximum (PETM), about 55.5 million years ago**1,2. This warming was associated with a massive addition of carbon to the ocean-atmosphere system, but estimates of the Earth systemresponse to this perturbation are complicated by widely varying estimates of the duration of carbon release, which range from less than a year to tens of thousands of years. In addition the source of the carbon, and whether it was released as a single injection or in several pulses, remains the subject of debate**2-4. Here we present a new high-resolution carbon isotope record from terrestrial deposits in the Bighorn Basin (Wyoming, USA) spanning the PETM, and interpret the record using a carbon-cycle boxmodel of the ocean-atmosphere-biosphere system.Our record shows that the beginning of the PETMis characterized by not one but two distinct carbon release events, separated by a recovery to background values. To reproduce this pattern, our model requires two discrete pulses of carbon released directly to the atmosphere, at average rates exceeding 0.9 Pg C yr**-1, with the first pulse lasting fewer than 2,000 years.

Clay mineraology and grain size composition of surface sediments from South Pacific Ocean

The reconstruction of low-latitude ocean-atmosphere interactions is one of the major issues of (paleo-)environmental studies. The trade winds, extending over 20° to 30° of latitude in both hemispheres, between the subtropical highs and the intertropical convergence zone, are major components of the atmospheric circulation and little is known about their long-term variability on geological time-scales, in particular in the Pacific sector. We present the modern spatial pattern of eolian-derived marine sediments in the eastern equatorial and subtropical Pacific (10°N to 25°S) as a reference data set for the interpretation of SE Pacific paleo-dust records. The terrigenous silt and clay fractions of 75 surface sediment samples have been investigated for their grain-size distribution and clay-mineral compositions, respectively, to identify their provenances and transport agents. Dust delivered to the southeast Pacific from the semi- to hyper-arid areas of Peru and Chile is rather fine-grained (4-8 µm) due to low-level transport within the southeast trade winds. Nevertheless, wind is the dominant transport agent and eolian material is the dominant terrigenous component west of the Peru-Chile Trench south of ~ 5°S. Grain-size distributions alone are insufficient to identify the eolian signal in marine sediments due to authigenic particle formation on the sub-oceanic ridges and abundant volcanic glass around the Galapagos Islands. Together with the clay-mineral compositions of the clay fraction, we have identified the dust lobe extending from the coasts of Peru and Chile onto Galapagos Rise as well as across the equator into the doldrums. Illite is a very useful parameter to identify source areas of dust in this smectite-dominated study area.

Radiocarbon ages and plateau definitions of sediments from the abyssal ocean

Ice core records demonstrate a glacial-interglacial atmospheric CO2 increase by ~100 ppm, while 14C calibration efforts document a strong decrease in atmospheric 14C concentration during this period. A calculated transfer of ~530 Gt of 14C depleted carbon is required to produce the deglacial coeval rise of carbon in the atmosphere and terrestrial biosphere. This amount is usually ascribed to oceanic carbon release, although the actual mechanisms remained elusive, since an adequately old and carbon-enriched deep-ocean reservoir seemed unlikely. Here we present a new, though still fragmentary, ocean-wide d14C dataset showing that during the Last Glacial Maximum (LGM) and Heinrich Stadial 1 (HS-1) the maximum 14C age difference between ocean deep waters and the atmosphere exceeded the modern values by up to 1500 14C yr, in the extreme reaching 5100 14C yr. Below 2000 m depth the 14C ventilation age of modern ocean waters is directly linked to the concentration of dissolved inorganic carbon (DIC). We propose as working hypothesis that the modern regression of DIC vs d14C also applies for LGM times, which implies that a mean LGM aging by ~600 14C yr corresponded to a global rise of ~85-115 µmol DIC/kg in the deep ocean. Thus, the prolonged residence time of ocean deep waters may indeed have made it possible to absorb an additional ~730-980 Gt DIC, one third of which possibly originated from intermediate waters. We also infer that LGM deep-water O2 dropped to suboxic values of <10µmol/kg in the Atlantic sector of the Southern Ocean, possibly also in the subpolar North Pacific. The outlined deglacial transfer of the extra aged, deep-ocean carbon to the atmosphere via the dynamic ocean-atmosphere carbon exchange would be sufficient to account for two trends observed, (1) for the increase in atmospheric CO2 and (2) for the 190-permil drop in atmospheric d14C during the so-called HS-1 'Mystery Interval', when atmospheric 14C production rates were largely constant.

Coral oxygen isotope and Sr/Ca data from the Northern Gulf of Aqaba (Red Sea)

The last interglacial period (about 125,000 years ago) is thought to have been at least as warm as the present climate (Kukla et al., 2002, doi:10.1006/qres.2001.2316). Owing to changes in the Earth's orbit around the Sun, it is thought that insolation in the Northern Hemisphere varied more strongly than today on seasonal timescales (Berger, 1987, doi:10.1175/1520-0469(1978)035<2362:LTVODI>2.0.CO;2), which would have led to corresponding changes in the seasonal temperature cycle (Montoya et al., 2000, doi:10.1175/1520-0442(2000)013<1057:CSFKBW>2.0.CO;2). Here we present seasonally resolved proxy records using corals from the northernmost Red Sea, which record climate during the last interglacial period, the late Holocene epoch and the present. We find an increased seasonality in the temperature recorded in the last interglacial coral. Today, climate in the northern Red Sea is sensitive to the North Atlantic Oscillation (Felis et al., 2000 doi:10.1029/1999PA000477; Rimbu et al., 2001, doi:10.1029/2001GL013083), a climate oscillation that strongly influences winter temperatures and precipitation in the North Atlantic region. From our coral records and simulations with a coupled atmosphere-ocean circulation model, we conclude that a tendency towards the high-index state of the North Atlantic Oscillation during the last interglacial period, which is consistent with European proxy records (Zagwijn, 1996, doi:10.1016/0277-3791(96)00011-X; Aalbersberg and Litt, 1998, doi:10.1002/(SICI)1099-1417(1998090)13:5<367::AID-JQS400>3.0.CO;2-I; Klotz et al., 2003, doi:10.1016/S0921-8181(02)00222-9), contributed to the larger amplitude of the seasonal cycle in the Middle East.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Age model and grain size distribution of the fine fraction of ODP Site 114-704

The physical properties of sediments beneath an upwelling area in the southern part of the Atlantic Ocean (ODP Hole 704A) were investigated. Highly significant correlations characterize the relationship of carbonate content to bulk density (R = 0.85), carbonate content to porosity (R = 0.84), and carbonate content to impedance (R = 0.84). No relationship exists between carbonate content and compressional-wave velocity (R = 0.24), indicating that amplitude variations in impedance are primarily controlled by variations in bulk density, which, in turn, are controlled by climatically driven biogenic opal and carbonate deposition. In general, maxima in impedance correspond to maxima in carbonate content (minima in opal content). The impedance record exhibits its most drastic change at about 2.4 Ma, marking dramatic increases in the average content of biogenic opal and the beginning of large-amplitude fluctuations. Between 0.7 and 0.4 Ma carbonate content, bulk density, and grain density decrease while opal content drastically increases. Similar changes have been observed in sediments beneath an upwelling cell off northwest Africa, indicating an oceanwide enhancement in upwelling or in the calcite corrosiveness of bottom water that appears to be isochronous.

Model experiment on optimal size of a fish stock with environmental uncertainties

We analyze the effect of environmental uncertainties on optimal fishery management in a bio-economic fishery model. Unlike most of the literature on resource economics, but in line with ecological models, we allow the different biological processes of survival and recruitment to be affected differently by environmental uncertainties. We show that the overall effect of uncertainty on the optimal size of a fish stock is ambiguous, depending on the prudence of the value function. For the case of a risk-neutral fishery manager, the overall effect depends on the relative magnitude of two opposing effects, the 'convex-cost effect' and the 'gambling effect'. We apply the analysis to the Baltic cod and the North Sea herring fisheries, concluding that for risk neutral agents the net effect of environmental uncertainties on the optimal size of these fish stocks is negative, albeit small in absolute value. Under risk aversion, the effect on optimal stock size is positive for sufficiently high coefficients of constant relative risk aversion.

Benthic foraminiferal assemblage changes during interstadial 4 to 8 in ODP Site 146-893A, Santa Barbara Basin

Benthic foraminiferal assemblages from Santa Barbara Basin exhibit major faunal and ecological switches associated with late Quaternary millennial- to decadal-scale global climate oscillations. Repeated turnovers of entire faunas occurred rapidly (<40-400 yr) without extinction or speciation in conjunction with Dansgaard-Oeschger shifts in thermohaline circulation, ventilation, and climate, confirming evolutionary model predictions of Roy et al. Consistent faunal successions of dysoxic taxa during successive interstadials reflect the extreme sensitivity and adaptation of the benthic ecosystem to the rapid environmental changes that marked the late Quaternary and possibly other transitional intervals in the history of the Earth's ocean-atmosphere-cryosphere system. These data support the hypothesis that broad segments of the biosphere are well adapted to rapid climate change.

Age model and stable isotope record of DSDP holes 48-401 and 86-577

Early Paleogene warm climates may have been linked to different modes and sources of deepwater formation. Warm polar temperatures of the Paleocene and Eocene may have resulted from either increased atmospheric trace gases or increased heat transport through deep and intermediate waters. The possibility of increasing ocean heat transport through the production of warm saline deep waters (WSDW) in the Tethyan region has generated considerable interest. In addition, General Circulation Model results indicate that deepwater source regions may be highly sensitive to changing basin configurations. To decipher deepwater changes, we examined detailed benthic foraminiferal faunal and isotopic records of the late Paleocene through the early Eocene (~60 to 50 Ma) from two critical regions: the North Atlantic (Bay of Biscay Site 401) and the Pacific (Shatsky Rise Site 577). These records are compared with published data from the Southern Ocean (Maud Rise Site 690, Islas Orcadas Rise Site 702). During the late Paleocene, similar benthic foraminiferal delta18O values were recorded at all four sites. This indicates uniform deepwater temperatures, consistent with a single source of deep water. The highest delta13C values were recorded in the Southern Ocean and were 0.5 per mil more positive than those of the Pacific. We infer that the Southern Ocean was proximal to a source of nutrient-depleted deep water during the late Paleocene. Upper Paleocene Reflector Ab was cut on the western Bermuda Rise by cyclonically circulating bottom water, also suggesting a vigorous source of bottom water in the Southern Ocean. A dramatic negative excursion in both carbon and oxygen isotopes occurred in the latest Paleocene in the Southern Ocean. This is a short-term (<100 kyr), globally synchronous event which also is apparent in both the Atlantic and Pacific records as a carbon isotopic excursion of approximately 1 per mil. Faunal analyses from the North Atlantic and Pacific sites indicate that the largest benthic foraminiferal faunal turnover of the Cenozoic was synchronous with the isotopic excursion, lending support to the hypothesis that the extinctions were caused by a change in deepwater circulation. We speculate that the Southern Ocean deepwater source was reduced or eliminated at the time of the excursion. During the early Eocene, Southern Ocean delta13C values remained enriched relative to the North Atlantic and Pacific. However, the Southern Ocean was also enriched in delta18O relative to these basins. We interpret that these patterns indicate that although the Southern Ocean was proximal to a source of cool, nutrient-depleted water, the intermediate to upper deep water sites of the North Atlantic and Pacific were ventilated by a different source that probably originated in low latitudes, i.e., WSDW.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Barium intensities and tie points between holes of Leg 208

The Paleocene - Eocene thermal maximum (PETM) is one of the best known examples of a transient climate perturbation, associated with a brief, but intense, interval of global warming and a massive perturbation of the global carbon cycle from injection of isotopically light carbon into the ocean-atmosphere system. One key to quantifying the mass of carbon released, identifying the source(s), and understanding the ultimate fate of this carbon is to develop high-resolution age models. Two independent strategies have been employed, cycle stratigraphy and analysis of extraterrestrial Helium (HeET), both of which were first tested on Ocean Drilling Program (ODP) Site 690. Both methods are in agreement for the onset of the PETM and initial recovery, or the clay layer ("main body"), but seem to differ in the final recovery phase of the event above the clay layer, where the carbonate contents rise and carbon isotope values return toward background values. Here we present a state-of-the-art age model for the PETM derived from a new orbital chronology developed with cycle stratigraphic records from sites drilled during ODP Leg 208 (Walvis Ridge, Southeastern Atlantic) integrated with published records from Site 690 (Weddell Sea, Southern Ocean, ODP Leg 113). During Leg 208, five Paleocene - Eocene (P-E) boundary sections (Sites 1262 to 1267) were recovered in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge yielding the first stratigraphically complete P-E deep-sea sequence with moderate to relatively high sedimentation rates (1 to 3 cm/kyr). A detailed chronology was developed with non-destructive X-ray fluorescence (XRF) core scanning records on the scale of precession cycles, with a total duration of the PETM now estimated to be ~ 170 kyr. The revised cycle stratigraphic record confirms original estimates for the duration of the onset and initial recovery, but suggests a new duration for the final recovery that is intermediate to the previous estimates by cycle stratigraphy and HeET.

Coral Sr/Ca records from Bonaire Holocene corals

Proxy reconstructions of tropical Atlantic sea surface temperature (SST) that extend beyond the period of instrumental observations have primarily focused on centennial to millennial variability rather than on seasonal to multidecadal variability. Here we present monthly-resolved records of Sr/Ca (a proxy of SST) from fossil annually-banded Diploria strigosa corals from Bonaire (southern Caribbean Sea). The individual corals provide time-windows of up to 68 years length, and the total number of 295 years of record allows for assessing the natural range of seasonal to multidecadal SST variability in the western tropical Atlantic during snapshots of the mid- to late Holocene. Comparable to modern climate, the coral Sr/Ca records reveal that mid- to late Holocene SST was characterised by clear seasonal cycles, persistent quasi-biennial and prominent interannual as well as inter- to multidecadal-scale variability. However, the magnitude of SST variations on these timescales has varied over the last 6.2 ka. The coral records show increased seasonality during the mid-Holocene consistent with climate model simulations indicating that southern Caribbean SST seasonality is induced by insolation changes on orbital timescales, whereas internal dynamics of the climate system play an important role on shorter timescales. Interannual SST variability is linked to ocean-atmosphere interactions of Atlantic and Pacific origin. Pronounced interannual variability in the western tropical Atlantic is indicated by a 2.35 ka coral, possibly related to a strengthening of the variability of the El Niño/Southern Oscillation throughout the Holocene. Prominent inter- to multidecadal SST variability is evident in the coral records and slightly more pronounced in the mid-Holocene. We finally argue that our coral data provide a target for studying Holocene climate variability on seasonal and interannual to multidecadal timescales, when using further numerical models and high-resolution proxy data.

Negative effects of ocean acidification on calcification vary within the coccolithophore genus Calcidiscus

A large percentage of CO2 emitted into the atmosphere is absorbed by the oceans, causing chemical changes in surface waters known as ocean acidification (OA). Despite the high interest and increased pace of OA research to understand the effects of OA on marine organisms, many ecologically important organisms remain unstudied. Calcidiscus is a heavily calcified coccolithophore genus that is widespread and genetically and morphologically diverse. It contributes substantially to global calcium carbonate production, organic carbon production, oceanic carbon burial, and ocean-atmosphere CO2 exchange. Despite the importance of this genus, relatively little work has examined its responses to OA. We examined changes in growth, morphology, and carbon allocation in multiple strains of Calcidiscus leptoporus in response to ocean acidification. We also, for the first time, examined the OA response of Calcidiscus quadriperforatus, a larger and more heavily calcified Calcidiscus congener. All Calcidiscus coccolithophores responded negatively to OA with impaired coccolith morphology and a decreased ratio of particulate inorganic to organic carbon (PIC:POC). However, strains responded variably; C. quadriperforatus showed the most sensitivity, while the most lightly calcified strain of C. leptoporus showed little response to OA. Our findings suggest that calcium carbonate production relative to organic carbon production by Calcidiscus coccolithophores may decrease in future oceans and that Calcidiscus distributions may shift if more resilient strains and species become dominant in assemblages. This study demonstrates that variable responses to OA may be strain or species specific in a way that is closely linked to physiological traits, such as cellular calcite quota.

Age determination of sediment cores of the California Current System

Alkenone sea surface temperature (SST) records were generated from the Ocean Drilling Program's (ODP) Sites 1014 and 1016 to examine the response of the California Current System to global climate change during the last 136 ka. The temperature differences between these sites (Delta SST(NEP)=SST(ODP1014)-SST(ODP1016)) reflected the intensity of the California Current and varied between 0.4 and 6.1 °C. A high Delta SST(NEP) (weaker California Current) was found for late marine isotope stage (MIS) 2 and early MIS 5e, while a low Delta SST(NEP) (stronger California Current) was detected for mid-MIS 5e and MIS 1. Spectral analysis indicated that this variation pattern dominated 23- (precession) and 30-ka periods. Comparison of the Delta SST(NEP) and SST based on data from core MD01-2421 at the Japan margin revealed anti-phase variation; the high Delta SST(NEP) (weakening of the California Current) corresponded to the low SST at the Japan margin (the southward displacement of the NW Pacific subarctic boundary), and vice versa. This variation was synchronous with a model prediction of the tropical El Niño-Southern Oscillation behavior. These findings suggest that the intensity of the North Pacific High varied in response to precessional forcing, and also that the response has been linked with the changes of tropical ocean-atmosphere interactions.