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.

Age model and stable oxygen isotope mean values of ice core SZ99 from Akademii Nauk ice cap on Severnaya Zemlya

Between 1999 and 2001, a 724 m long ice core was drilled on Akademii Nauk, the largest glacier on Severnaya Zemlya, Russian Arctic. The drilling site is located near the summit. The core is characterized by high melt-layer content. The melt layers are caused by melting and even by rain during the summer. We present high-resolution data of density, electrical conductivity (dielectrical profiling), stable water isotopes and melt-layer content for the upper 136 m (120 m w.e.) of the ice core. The dating by isotopic cycles and electrical conductivity peak identification suggests that this core section covers approximately the past 275 years. Singularities of volcanogenic and anthropogenic origin provide well-defined additional time markers. Long-term temperatures inferred from 12 year running mean averages of d18O reach their lowest level in the entire record around 1790. Thereafter the d18O values indicate a continuously increasing mean temperature on the Akademii Nauk ice cap until 1935, interrupted only by minor cooling episodes. The 20th century is found to be the warmest period in this record.

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.

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.

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.

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.

(Table 1) Age model of sediment core MD03-2705

A reconstruction of northwest African summer monsoon strength during the cold marine isotopic stage (MIS) 6 indicates a link to the seasonal migration of the Intertropical Convergence Zone (ITCZ). High-resolution studies of eolian dust supply and sea surface temperature recorded in marine core MD03-2705, on the Mauritanian margin, provide a better understanding about the penultimate glacial history of northwestern African aridity/humidity and upwelling coastal activity. Today, site MD03-2705 experiences increased upwelling and dust flux during the winter months, when the ITCZ is in a southerly position. Analyses of foraminifera isotopic composition suggest that during MIS 6.5 (180-168 ka) the average position of the ITCZ migrated north, marked by an increase in the strength of the summer monsoon, which decreased eolian dust transport and the coastal upwelling activity. The northward migration is in phase with a specific orbital combination of a low precessional index with a high obliquity signal. High-resolution analysis of stable isotopes (d18O and d13C) and microscale resolution geochemical (Ti/Al and quartz grain counts) determinations reveal that the transition between monsoonal humid (MIS 6.5) and dry (MIS 6.4) conditions has occurred in less than 1.3 ka. Such rapid changes suggest a nonlinear link between the African monsoonal rainfall system and environmental changes over the continent. This study provides new insights about the influence of vegetation and oceanic temperature feedbacks on the onset of African summer monsoon and demonstrates that, during the penultimate glacial period, changes in tropical dynamics had regional and global impacts.

Age model and stable isotopes of sediment cores from north-central Chile and south-central Chile

Terrigenous sediment supply, marine transport, and depositional processes along tectonically active margins are key to decoding turbidite successions as potential archives of climatic and seismic forcings. Sequence stratigraphic models predict coarse-grained sediment delivery to deep-marine sites mainly during sea-level fall and lowstand. Marine siliciclastic deposition during transgressions and highstands has been attributed to sustained connectivity between terrigenous sources and marine sinks facilitated by narrow shelves. To decipher the controls on Holocene highstand turbidite deposition, we analyzed 12 sediment cores from spatially discrete, coeval turbidite systems along the Chile margin (29° - 40°S) with changing climatic and geomorphic characteristics but uniform changes in sea level. Sediment cores from intraslope basins in north-central Chile (29° - 33°S) offshore a narrow to absent shelf record a shut-off of turbidite deposition during the Holocene due to postglacial aridification. In contrast, core sites in south-central Chile (36° - 40°S) offshore a wide shelf record frequent turbidite deposition during highstand conditions. Two core sites are linked to the Biobío river-canyon system and receive sediment directly from the river mouth. However, intraslope basins are not connected via canyons to fluvial systems but yield even higher turbidite frequencies. High sediment supply combined with a wide shelf and an undercurrent moving sediment toward the shelf edge appear to control Holocene turbidite sedimentation and distribution. Shelf undercurrents may play an important role in lateral sediment transport and supply to the deep sea and need to be accounted for in sediment-mass balances.

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.

Sea ice thickness distribution (ITD) Model for the Arctic, links to NetCDF files

Sea ice models contain many different parameterizations of which one of the most commonly used is a subgrid-scale ice thickness distribution (ITD). The effect of this model component and the associated ice strength formulation on the reproduction of observed Arctic sea ice is assessed. To this end the model's performance in reproducing satellite observations of sea ice concentration, thickness and drift is evaluated. For an unbiased comparison, different model configurations with and without an ITD are tuned with an automated parameter optimization. The original combination of ITD and ice strength parameterization does not lead to better results than a simple single category model. Yet changing to a simpler ice strength formulation, which depends linearly on the mean ice thickness across all thickness categories, allows to clearly improve the model-data misfit when using an ITD. In the original formulation, the ice strength depends strongly on the number of thickness categories, so that introducing more categories can lead to thicker albeit weaker ice on average.