Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Stable isotope record, Mg/Ca ratios and age model of St. Stephens Quarry, Alabama

In the largest global cooling event of the Cenozoic Era, between 33.8 and 33.5 Myr ago, warm, high-CO2 conditions gave way to the variable 'icehouse' climates that prevail today. Despite intense study, the history of cooling versus ice-sheet growth and sea-level fall reconstructed from oxygen isotope values in marine sediments at the transition has not been resolved. Here, we analyse oxygen isotopes and Mg/Ca ratios of benthic foraminifera, and integrate the results with the stratigraphic record of sea-level change across the Eocene-Oligocene transition from a continental-shelf site at Saint Stephens Quarry, Alabama. Comparisons with deep-sea (Sites 522 (South Atlantic) and 1218 (Pacific)) d18O and Mg/Ca records enable us to reconstruct temperature, ice-volume and sea-level changes across the climate transition. Our records show that the transition occurred in at least three distinct steps, with an increasing influence of ice volume on the oxygen isotope record as the transition progressed. By the early Oligocene, ice sheets were ~25% larger than present. This growth was associated with a relative sea-level decrease of approximately 105 m, which equates to a 67 m eustatic fall.

Lithology and organic geochemistry at DSDP Hole 76-534

The organic facies of Early and middle Cretaceous sediments drilled at DSDP Site 534 is dominated by terrestrially derived plant remains and charcoal. Marine organic matter is mixed with the terrestrial components, but through much of this period was diluted by the terrestrial material. The supply of terrestrial organic matter was high here because of the nearness of the shore and high runoff promoted by a humid temperate coastal climate. Reducing conditions favored preservation of both marine and terrestrial organic matter, the terrestrial materials having reached the site mostly in turbidity currents or in the slow-moving, near-bottom nepheloid layer. An increase in the abundance of terrestrial organic matter occurred when the sea level dropped in the Valanginian and again in the Aptian-Albian, because rivers dumped more terrigenous elastics into the Basin and marine productivity was lower at these times than when sea level was high. A model is proposed to explain the predominance of reducing conditions in the Valanginian-Aptian, of oxidizing conditions in the late Aptian, and of reducing conditions in the Albian-Cenomanian. The model involves influx of oxygen-poor subsurface waters from the Pacific at times of high or rising sea level (Valanginian-Aptian, and Albian- Cenomanian) and restriction of that influx at times of low sea level (late Aptian). In the absence of a supply of oxygenpoor deep water, the bottom waters of the North Atlantic became oxidizing in the late Aptian, probably in response to development of a Mediterranean type of circulation. The influx of nutrients from the Pacific led to an increase in productivity through time, accounting for an increase in the proportion of marine organic matter from the Valanginian into the Aptian and from the Albian to the Cenomanian. Conditions were dominantly oxidizing through the Middle Jurassic into the Berriasian, with temporary exceptions when bottom waters became reducing, as in the Callovian. Mostly terrestrial and some marine organic matter accumulated during the Callovian reducing episode. When Jurassic bottom waters were oxidizing, only terrestrial organic matter was buried in the sediments, in very small amounts.

Freshwater runoff for Taylor Slough to Florida Bay for 1965 to 2000 for the FATHOM model

The Florida Bay ecosystem supports a number of economically important ecosystem services, including several recreational fisheries, which may be affected by changing salinity and temperature due to climate change. In this paper, we use a combination of physical models and habitat suitability index models to quantify the effects of potential climate change scenarios on a variety of juvenile fish and lobster species in Florida Bay. The climate scenarios include alterations in sea level, evaporation and precipitation rates, coastal runoff, and water temperature. We find that the changes in habitat suitability vary in both magnitude and direction across the scenarios and species, but are on average small. Only one of the seven species we investigate (Lagodon rhomboides, i.e., pinfish) sees a sizable decrease in optimal habitat under any of the scenarios. This suggests that the estuarine fauna of Florida Bay may not be as vulnerable to climate change as other components of the ecosystem, such as those in the marine/terrestrial ecotone. However, these models are relatively simplistic, looking only at single species effects of physical drivers without considering the many interspecific interactions that may play a key role in the adjustment of the ecosystem as a whole. More complex models that capture the mechanistic links between physics and biology, as well as the complex dynamics of the estuarine food web, may be necessary to further understand the potential effects of climate change on the Florida Bay ecosystem.

Grain size composition and age of ODP Hole 167-1017E sediments

Siliciclastic sedimentation at Ocean Drilling Program Site 1017 on the southern slope of the Santa Lucia Bank, central California margin, responded closely to oceanographic and climatic change over the past ~130 ka. Variation in mean grain-size and sediment sorting within the ~25-m-thick succession from Hole 1017E show Milankovitch-band to submillenial-scale variation. Mean grain size of the "sortable silt" fraction (10-63 µm) ranges from 17.6 to 33.9 µm (average 24.8 µm) and is inversely correlated with the degree of sorting. Much of the sediment has a bimodal or trimodal grain-size distribution that is composed of distinct fine silt, coarse silt to fine sand, and clay-size components. The position of the mode and the sorting of each component changes through the succession, but the primary variation is in the presence or abundance of the coarse silt fraction that controls the overall mean grain size and sorting of the sample. The occurrence of the best-sorted, finest grained sediment at high stands of sea level (Holocene, marine isotope Substages 5c and 5e) reflect the linkage between global climate and the sedimentary record at Site 1017 and suggest that the efficiency of off-shelf transport is a key control of sedimentation on the Santa Lucia Slope. It is not clear what proportion of the variation in grain size and sorting may also be caused by variations in bottom current strength and in situ hydrodynamic sorting.

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.

Glacial-interglacial record of ODP Site 181-1119 in the Canterbury Bright, Southwest Pacific Ocean

Ocean Drilling Program (ODP) Site 1119 is located at water depth 395 m near the subtropical front (STF; here represented by the Southland Front), just downslope from the shelf edge of eastern South Island, New Zealand. The upper 86.19 metres composite depth (mcd) of Site 1119 sediment was deposited at an average sedimentation rate of 34 cm/kyr during Marine Isotope Stages (MIS) 1-8 (0-252 ka), and is underlain across a ~25 kyr intra-MIS 8 unconformity by MIS 8.5-11 (277-367 ka) and older sediment deposited at ~14 cm/kyr. A time scale is assigned to Site 1119 using radiocarbon dates for the period back to ~39 ka, and, prior to then, by matching its climatic record with that of the Vostok ice core, which it closely resembles. Four palaeoceanographic proxy measures for surface water masses vary together with the sandy-muddy, glacial-interglacial (G/I) cyclicity at the site. Interglacial intervals are characterised by heavy delta13C, high colour reflectance (a proxy for carbonate content), low Q-ray (a proxy for clay content) and light delta18O; conversely, glacial intervals exhibit light delta13C, low reflectance, high Q-ray and heavy delta18O signatures. Early interglacial intervals are represented by silty clays with 10-105-cm-thick beds of sharp-based (Chondrites-burrowed), shelly, graded, fine sand. The sands are rich in foraminifera, and were deposited distant from the shoreline under the influence of longitudinal flow in relatively deep water. Glacial intervals comprise mostly micaceous silty clay, though with some thin (2-10 cm thick) sands present also at peak cold periods, and contain the cold-water scallop Zygochlamys delicatula. Interglacial sandy intervals are characterised by relatively low sedimentation rates of 5-32 cm/kyr; cold climate intervals MIS 10, 6 and 2 have successively higher sedimentation rates of 45, 69 and 140 cm/kyr. Counter-intuitively,and forced by the bathymetric control of a laterally-moving shoreline during G/I and I/G transitions, the 1119 core records a southeasterly (seaward) movement of the STF during early glacial periods, accompanied by the incursion of subtropical water (STW) above the site, and northwesterly (landward) movement during late glacial and interglacial times, resulting in a dominant influence then of subantarctic surface water (SAW). The history of passage of these different water masses at the site is clearly delineated by their characteristic delta13C values. The intervals of thin, graded sands-muds which occur within MIS 2-3, 6, 7.4 and 10 indicate the onset at times of peak cold of intermittent bottom currents caused by strengthened and expanded frontal flows along the STF, which at such times lay near Site 1119 in close proximity to seaward-encroaching subantarctic waters within the Bounty gyre. In common with other nearby Southern Hemisphere records, the cold period which represents the last glacial maximum lasted between ~23-18 ka at Site 1119, during which time the STF and Subantarctic Front (SAF) probably merged into a single intense frontal zone around the head of the adjacent Bounty Trough.

Core alignment and composite depth scale of DSDP Site 86-577

Deep Sea Drilling Project Site 577 on Shatsky Rise (North Pacific Ocean) recovered a series of cores at three holes that contain calcareous nannofossil ooze of latest Cretaceous (late Maastrichtian) through early Eocene age. Several important records have been generated using samples from these cores, but the stratigraphy has remained outdated and confusing. Here we revise the stratigraphy at Site 577. This includes refining several age datums, realigning cores in the depth domain, and placing all stratigraphic markers on a current time scale. The work provides a template for appropriately bringing latest Cretaceous and Paleogene data sets at old drill sites into current paleoceanographic literature for this time interval. While the Paleocene Eocene Thermal Maximum (PETM) lies within core gaps at Holes 577* and 577A, the sedimentary record at the site holds other important events and remains crucially relevant to understanding changes in oceanographic conditions from the latest Cretaceous through early Paleogene.

Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, links model results in MATLAB format

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

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.

Average fluorescence and dissolved iron and Fe-binding ligand characteristics during POLARSTERN expedition ANT-XXIV/3 in 2008

Organic complexation of dissolved iron (dFe) was investigated in the Atlantic sector of the Southern Ocean in order to understand the distribution of Fe over the whole water column. The total concentration of dissolved organic ligands ([Lt]) measured by voltammetry ranged between 0.54 and 1.84 nEq of M Fe whereas the conditional binding strength (K') ranged between 10**21.4 and 10**22.8. For the first time, trends in Fe-organic complexation were observed in an ocean basin by examining the ratio ([Lt]/[dFe]), defined as the organic ligand concentration divided by the dissolved Fe concentration. The [Lt]/[dFe] ratio indicates the saturation state of the natural ligands with Fe; a ratio near 1 means saturation of the ligands leading to precipitation of Fe. Reversely, high ratios mean Fe depletion and show a high potential for Fe solubilisation. In surface waters where phytoplankton is present low dissolved Fe and high variable ligand concentrations were found. Here the [Lt]/[dFe] ratio was on average 4.4. It was especially high (5.6-26.7) in the HNLC (High Nutrient, Low Chlorophyll) regions, where Fe was depleted. The [Lt]/[dFe] ratio decreased with depth due to increasing dissolved Fe concentrations and became constant below 450 m, indicating a steady state between ligand and Fe. Relatively low [Lt]/[dFe] ratios (between 1.1 and 2.7) existed in deep water north of the Southern Boundary, facilitating Fe precipitation. The [Lt]/[dFe] ratio increased southwards from the Southern Boundary on the Zero Meridian and from east to west in the Weddell Gyre due to changes both in ligand characteristics and in dissolved iron concentration. High [Lt]/[dFe] ratio expresses Fe depletion versus ligand production in the surface. The decrease with depth reflects the increase of [dFe] which favours scavenging and (co-) precipitation, whereas a horizontal increase in the deep waters results from an increasing distance from Fe sources. This increase in the [Lt]/[dFe] ratio at depth shows the very resistant nature of the dissolved organic ligands.

Tide amplitude at Ilulissat, and surface height and flow velocities of Jakobshavn Isbrae in August 2004

During the summer of 2004, the front area of the Jakobshavn Isbræ was monitored using a geodetic-photogrammetric survey with temporarily coincident precise observations of local ocean tides in the Disko Bay close to Ilulissat. The geodetic and photogrammetric observations were conducted at the southern margin of the glacier front. The largest observed horizontal flow velocities are in the central part of the front with values up to 45 m/d. This is a factor of 2 greater than the average velocities at the front area observed in the last century. Our new observations confirm previous estimates of an acceleration of glacier flow during the last decade. The photogrammetric survey provided flow trajectories for 4000 surface points with a time resolution of 30 min. These flow trajectories were used to compare the vertical motion of the glacier with the observed tides. The existence of a free-floating glacier tongue in 2004 was confirmed by these data. However, it occupied only a small belt, of at most a few 100 m width, in the central part of the glacier front. Horizontal motion did not appear to depend on the tidal phase, unlike some of the fast-moving ice streams of West Antarctica.

Model of interactive visualization of diffusion processes from different strains, and source routes and spatial diffusion of capuchin monkeys over the past 6 million years

Source routes and Spatial Diffusion of capuchin monkeys over the past 6 million years, rebuilt in the SPREAD 1.0.6 from the MCC tree. The map shows the 10 different regions to which distinctive samples were associated. The different transmission routes have been calculated from the average rate over time. Only rates with Bayes factor> 3 were considered as significantly different from zero. Significant diffusion pathways are highlighted with color varying from dark brown to red, being the dark brown less significant rates and deep red the most significant rates.

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.