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.

Molecular differentiation by PCR-RFLP technique shows that Antarctic fishes of the genus Notothenia, claimed to be two species, are phenotypic varieties with the same genetic pattern (Table 1)

The taxonomy of Antarctic fishes has been predominantly based on morphological characteristics rather than on genetic criteria. A typical example is the Notothenia group, which includes N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii Richardson, 1844. The Polymerase Chain Reaction and Restriction Fragment Length Polymorphism (PCR-RFLP) technique was used to determine whether N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951 are different or whether they are the same species with morphological, physiological and behavioural variability. N. rossii was used as control. Mitochondrial DNA (mtDNA) was isolated from muscle specimens of N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii, which were collected in Admiralty Bay, King George Island. The DNA was used to amplify a fragment (690 base pairs) of the mitochondrial gene coding region of NADH dehydrogenase subunit 2. Further, the amplicon was digested with the following restriction enzymes: DdeI, HindIII and RsaI. The results showed a variation of the digestion pattern of the fragment amplified between N. rossii, and N. coriiceps Richardson, 1844 or N. neglecta Nybelin, 1951. However, no differences were found between N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951, on the grounds of the same genetic pattern shown by the two fish.

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Age model for ODP Leg 175 sites

Upwelling along the western coast of Africa south of the equator may be partitioned into three major areas, each having its own dynamics and history: (1) the eastern equatorial region, comprising the Congo Fan and the area of Mid-Angola; (2) the Namibia upwelling system, extending from the Walvis Ridge to Lüderitz; and (3) the Cape Province region, where upwelling is subdued. The highest nutrient contents in thermocline waters are in the northern region, the lowest in the southern one. Wind effects are at a maximum near the southern end of the Namibia upwelling system, and maximum productivity occurs near Walvis Bay, where the product between upwelling rate and nutrient content of upwelled waters is at a maximum. In the Congo/Angola region, opal tends to follow organic carbon quite closely in the Quaternary record. However, organic carbon has a strong precessional component, while opal does not. Despite relatively low opal content, sediments off Angola show the same patterns as those off the Congo; thus, they are part of the same regime. The spectrum shows nonlinear interference patterns between high- and low-latitude forcing, presumably tied to thermocline fertility and wind. On Walvis Ridge, as in the Congo-Angola region, the organic matter record behaves normally; that is, supply is high during glacial periods. In contrast, interglacial periods are favorable for opal deposition. The pattern suggests reduction in silicate content of the thermocline during glacial periods. The reversed phase (opal abundant during interglacials) persists during the entire Pleistocene and can be demonstrated deep into the Pliocene, not just on Walvis Ridge but all the way to the Oranje River and off the Cape Province. From comparison with other regions, it appears that silicate is diminished in the global thermocline, on average, whenever winds become strong enough to substantially shorten the residence time of silicate in upper waters (Walvis Hypothesis, solving the Walvis Paradox of reversed phase in opal deposition). The central discovery during Leg 175 was the documentation of a late Pliocene opal maximum for the entire Namibia upwelling system (early Matuyama Diatom Maximum [MDM]). The maximum is centered on the period between the end of the Gauss Chron and the beginning of the Olduvai Chron. A rather sharp increase in both organic matter deposition and opal deposition occurs near 3 Ma in the middle of the Gauss Chron, in association with a series of major cooling steps. As concerns organic matter, high production persists at least to 1 Ma, when there are large changes in variability, heralding subsequent pulsed production periods. From 3 to 2 Ma, organic matter and opal deposition run more or less parallel, but after 2 Ma opal goes out of phase with organic matter. Apparently, this is the point when silicate becomes limiting to opal production. Thus, the MDM conundrum is solved by linking planetary cooling to increased mixing and upwelling (ramping up to the MDM) and a general removal of silicate from the upper ocean through excess precipitation over global supply (ramping down from the MDM). The hypothesis concerning the origin of the Namibia opal acme or MDM is fundamentally the same as the Walvis Hypothesis, stating that glacial conditions result in removal of silicate from the thermocline (and quite likely from the ocean as a whole, given enough time). The Namibia opal acme, and other opal maxima in the latest Neogene in other regions of the ocean, marks the interval when a cooling ocean selectively removes the abundant silicate inherited from a warm ocean. When the excess silicate is removed, the process ceases. According to the data gathered during Leg 175, major upwelling started in the late part of the late Miocene. Presumably, this process contributed to the drawing down of carbon dioxide from the atmosphere, helping to prepare the way for Northern Hemisphere glaciation.

Results of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958-2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Revised age model of DSDP Site 95-612

Alkenone-based Cenozoic records of the partial pressure of atmospheric carbon dioxide (pCO2) are founded on the carbon isotope fractionation that occurred during marine photosynthesis (epsilon [p37:2]). However, the magnitude of epsilon [p37:2] is also influenced by phytoplankton cell size - a consideration lacking in previous alkenone-based CO2 estimates. In this study, we reconstruct cell size trends in ancient alkenone-producing coccolithophores (the reticulofenestrids) to test the influence that cell size variability played in determining epsilon [p37:2] trends and pCO2 estimates during the middle Eocene to early Miocene. At the investigated deep-sea sites, the reticulofenestrids experienced high diversity and largest mean cell sizes during the late Eocene, followed by a long-term decrease in maximum cell size since the earliest Oligocene. Decreasing haptophyte cell sizes do not account for the long-term increase in the stable carbon isotopic composition of alkenones and associated decrease in epsilon [p37:2] values during the Paleogene, supporting the conclusion that the secular pattern of epsilon [p37:2] values is primarily controlled by decreasing CO2 concentration since the earliest Oligocene. Further, given the physiology of modern alkenone producers, and considering the timings of coccolithophorid cell size change, extinctions, and changes in reconstructed pCO2 and temperature, we speculate that the selection of smaller reticulofenestrid cells during the Oligocene primarily reflects an adaptive response to increased [CO2(aq)] limitation.

Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.

Age model, radiocarbon dates and dinoflagellate cyst measurements from sediment core DA06-139G

A marine sediment core from Vaigat in Disko Bugt, West Greenland, has been analysed in terms of lithology, dinoflagellate cysts and foraminifera in order to evaluate the influence of oceanographic variability on West Greenland glacier stability. The data show that during the past 5200 years the Atlantic foraminiferal abundance in the subsurface waters of the West Greenland Current (WGC) episodically increased, indicating periods of increases in the inflow of subsurface warm Atlantic water at 2000 - 1500 cal. yr BP and 1300 cal. yr BP as well as periods of less pronounced increased bottom-water temperatures around 4700 - 4000 cal. yr BP, 3100 - 2800, 2600, 1000 - 800, 500 - 400, and at 200 cal. yr. The sedimentological and dinoflagellate cyst data indicate that these episodes with enhanced advection of Irminger Sea-derived waters are accompanied by increased iceberg rafting, which we link to increased iceberg calving in relation to destabilization of the Jakobshavn Isbrae. The long-term trend in the data documents the end of a late-Holocene Thermal Maximum between 5200 and 4300 cal. yr BP and a final onset of the Neoglaciation at 3500 cal. yr BP. Increased responses of the iceberg rafting after 3500 cal. yr BP, reflects a westward/seaward advance of the glacier margin in relation to onset of Neoglaciation and a development of the glacier into a floating tongue after 2000 cal. yr BP. A comparison of our record with a record from the eastern North Atlantic indicates that a NAO-like anomaly pattern between subsurface waters in West Greenland and atmospheric temperature in the Eastern North Atlantic may have been operating during most of the late Holocene. However, during the past 1000 years the NAO signal may have weakened as some other mode of climate variability overprints the anti-phase climate signal in this region.

Age determinaton and age model of sediments from the Bay of Biscay

Continuous high-resolution mass accumulation rates (MAR) and X-ray fluorescence (XRF) measurements from marine sediment records in the Bay of Biscay (NE Atlantic) have allowed the determination of the timing and the amplitude of the 'Fleuve Manche' (Channel River) discharges during glacial stages MIS 10, MIS 8, MIS 6 and MIS 4-2. These results have yielded detailed insight into the Middle and Late Pleistocene glaciations in Europe and the drainage network of the western and central European rivers over the last 350 kyr. This study provides clear evidence that the 'Fleuve Manche' connected the southern North Sea basin with the Bay of Biscay during each glacial period and reveals that 'Fleuve Manche' activity during the glaciations MIS 10 and MIS 8 was significantly less than during MIS 6 and MIS 2. We correlate the significant 'Fleuve Manche' activity, detected during MIS 6 and MIS 2, with the extensive Saalian (Drenthe Substage) and the Weichselian glaciations, respectively, confirming that the major Elsterian glaciation precedes the glacial MIS 10. In detail, massive 'Fleuve Manche' discharges occurred at ca 155 ka (mid-MIS 6) and during Termination I, while no significant discharges are found during Termination II. It is assumed that a substantial retreat of the European ice sheet at ca 155 kyr, followed by the formation of ice-free conditions between the British Isles and Scandinavia until Termination II, allowed meltwater to flow northwards through the North Sea basin during the second part of the MIS 6. We assume that this glacial pattern corresponds to the Warthe Substage glacial maximum, therefore indicating that the data presented here equates to the Drenthe and the Warthe glacial advances at ca 175-160 ka and ca 150-140 ka, respectively. Finally, the correlation of our records with ODP site 980 reveals that massive 'Fleuve Manche' discharges, related to partial or complete melting of the European ice masses, were synchronous with strong decreases in both the rate of deep-water formation and the strength of the Atlantic thermohaline circulation. 'Fleuve Manche' discharges over the last 350 kyr probably participated, with other meltwater sources, in the collapse of the thermohaline circulation by freshening the northern Atlantic surface water.

(Table 2) Age model and sedimentation rate of sediment core GeoB5559-2

The Canary Basin lies in a region of strong interaction between the atmospheric and ocean circulation systems: Trade winds drive seasonal coastal upwelling and dust storm outbreaks from the neighbouring Sahara desert are the major source of terrigenous sediment. To investigate the forcing mechanisms for dust input and wind strength in the North Canary Basin, the temporal pattern of variability of sedimentological and geochemical proxy records has been analysed in two sediment cores between latitudes 30°30'N and 31°40'N. Spectral analysis of the dust proxy records indicates that insolation changes related to eccentricity and precession are the main periods of temporal variation in the record. Si/Al and grain-size of the terrigenous fraction show an increase in glacial-interglacial transitions while Al concentration and Fe/Al ratio are both in phase with minima in the precessional index. Hence, the results obtained show that the wind strength was intensified at Terminations. At times of maxima of Northern Hemisphere seasonal insolation, when the African monsoon was enhanced, the North Canary Basin also received higher dust input. This result suggests that the moisture brought by the monsoon may have increased the availability of dust in the source region.