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.

Wetland maps including open water extent dynamics based on ENVISAT ASAR WS for Siberia, 2007 and 2008, links to GeoTIFFs

Wetlands store large amounts of carbon, and depending on their status and type, they release specific amounts of methane gas to the atmosphere. The connection between wetland type and methane emission has been investigated in various studies and utilized in climate change monitoring and modelling. For improved estimation of methane emissions, land surface models require information such as the wetland fraction and its dynamics over large areas. Existing datasets of wetland dynamics present the total amount of wetland (fraction) for each model grid cell, but do not discriminate the different wetland types like permanent lakes, periodically inundated areas or peatlands. Wetland types differently influence methane fluxes and thus their contribution to the total wetland fraction should be quantified. Especially wetlands of permafrost regions are expected to have a strong impact on future climate due to soil thawing. In this study ENIVSAT ASAR Wide Swath data was tested for operational monitoring of the distribution of areas with a long-term SW near 1 (hSW) in northern Russia (SW = degree of saturation with water, 1 = saturated), which is a specific characteristic of peatlands. For the whole northern Russia, areas with hSW were delineated and discriminated from dynamic and open water bodies for the years 2007 and 2008. The area identified with this method amounts to approximately 300,000 km**2 in northern Siberia in 2007. It overlaps with zones of high carbon storage. Comparison with a range of related datasets (static and dynamic) showed that hSW represents not only peatlands but also temporary wetlands associated with post-forest fire conditions in permafrost regions. Annual long-term monitoring of change in boreal and tundra environments is possible with the presented approach. Sentinel-1, the successor of ENVISAT ASAR, will provide data that may allow continuous monitoring of these wetland dynamics in the future complementing global observations of wetland fraction.

Gotland Deep (Baltic Sea) independent age model based on Palaeomagnetic Secular Variation (PSV) and lead (Pb) deposition profiles

Dating of sediment cores from the Baltic Sea has proven to be difficult due to uncertainties surrounding the 14C reservoir age and a scarcity of macrofossils suitable for dating. Here we present the results of multiple dating methods carried out on cores in the Gotland Deep area of the Baltic Sea. Particular emphasis is placed on the Littorina stage (8 ka ago to the present) of the Baltic Sea and possible changes in the 14C reservoir age of our dated samples. Three geochronological methods are used. Firstly, palaeomagnetic secular variations (PSV) are reconstructed, whereby ages are transferred to PSV features through comparison with varved lake sediment based PSV records. Secondly, lead (Pb) content and stable isotope analysis are used to identify past peaks in anthropogenic atmospheric Pb pollution. Lastly, 14C determinations were carried out on benthic foraminifera (Elphidium spec.) samples from the brackish Littorina stage of the Baltic Sea. Determinations carried out on smaller samples (as low as 4 µg C) employed an experimental, state-of-the-art method involving the direct measurement of CO2 from samples by a gas ion source without the need for a graphitisation step - the first time this method has been performed on foraminifera in an applied study. The PSV chronology, based on the uppermost Littorina stage sediments, produced ten age constraints between 6.29 and 1.29 cal ka BP, and the Pb depositional analysis produced two age constraints associated with the Medieval pollution peak. Analysis of PSV data shows that adequate directional data can be derived from both the present Littorina saline phase muds and Baltic Ice Lake stage varved glacial sediments. Ferrimagnetic iron sulphides, most likely authigenic greigite (Fe3S4), present in the intermediate Ancylus Lake freshwater stage sediments acquire a gyroremanent magnetisation during static alternating field (AF) demagnetisation, preventing the identification of a primary natural remanent magnetisation for these sediments. An inferred marine reservoir age offset (deltaR) is calculated by comparing the foraminifera 14C determinations to a PSV & Pb age model. This deltaR is found to trend towards younger values upwards in the core, possibly due to a gradual change in hydrographic conditions brought about by a reduction in marine water exchange from the open sea due to continued isostatic rebound.

Age model of ODP Hole 107-654A

It has long been speculated that glacio-eustatic sea level drop may have caused or contributed to the isolation and consequent desiccation of the Mediterranean in the late Miocene (the 'Messinian salinity crisis'). Ocean Drilling Program site 654 on the Sardinia margin of the Tyrrhenian Sea is the first deep-sea drill site to penetrate through upper Messinian evaporites into lower Messinian/upper Tortonian open marine sediments, and thus offers a unique opportunity to date the onset of the salinity crisis. A reexamination of the magnetostratigraphic, biostratigraphic, and stable isotope-stratigraphic constraints on the preevaporite sediments of site 654 has yielded two possible ages for the contact between salinity crisis sediments and the underlying normal marine sediments. One magnetostratigraphic interpretation plus the biostratigraphically determined position of the Tortonian/Messinian boundary imply a date of about 6.2 Ma for the youngest presalinity crisis sediments. An alternative magnetostratigraphic interpretation plus the carbon isotope stratigraphy imply a date of about 5.2 Ma. The younger of these dates coincides with a delta18O spike in open ocean sediments [Keigwin, 1987 doi:10.1029/PA002i006p00639], which is attributed to increased ice volume.

Model parameter

Interannual environmental variability in Peru is dominated by the El Niño Southern Oscillation (ENSO). The most dramatic changes are associated with the warm El Niño (EN) phase (opposite the cold La Niña phase), which disrupts the normal coastal upwelling and affects the dynamics of many coastal marine and terrestrial resources. This study presents a trophic model for Sechura Bay, located at the northern extension of the Peruvian upwelling system, where ENSO-induced environmental variability is most extreme. Using an initial steady-state model for the year 1996, we explore the dynamics of the ecosystem through the year 2003 (including the strong EN of 1997/98 and the weaker EN of 2002/03). Based on support from literature, we force biomass of several non-trophically-mediated 'drivers' (e.g. Scallops, Benthic detritivores, Octopus, and Littoral fish) to observe whether the fit between historical and simulated changes (by the trophic model) is improved. The results indicate that the Sechura Bay Ecosystem is a relatively inefficient system from a community energetics point of view, likely due to the periodic perturbations of ENSO. A combination of high system productivity and low trophic level target species of invertebrates (i.e. scallops) and fish (i.e. anchoveta) results in high catches and an efficient fishery. The importance of environmental drivers is suggested, given the relatively small improvements in the fit of the simulation with the addition of trophic drivers on remaining functional groups' dynamics. An additional multivariate regression model is presented for the scallop Argopecten purpuratus, which demonstrates a significant correlation between both spawning stock size and riverine discharge-mediated mortality on catch levels. These results are discussed in the context of the appropriateness of trophodynamic modeling in relatively open systems, and how management strategies may be focused given the highly environmentally influenced marine resources of the region.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

Simulated Methyl iodide emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by OMIP data

Production pathways of the prominent volatile organic halogen compound methyl iodide (CH3I) are not fully understood. Based on observations, production of CH3I via photochemical degradation of organic material or via phytoplankton production has been proposed. Additional insights could not be gained from correlations between observed biological and environmental variables or from biogeochemical modeling to identify unambiguously the source of methyl iodide. In this study, we aim to address this question of source mechanisms with a three-dimensional global ocean general circulation model including biogeochemistry (MPIOM-HAMOCC (MPIOM - Max Planck Institute Ocean Model HAMOCC - HAMburg Ocean Carbon Cycle model)) by carrying out a series of sensitivity experiments. The simulated fields are compared with a newly available global data set. Simulated distribution patterns and emissions of CH3I differ largely for the two different production pathways. The evaluation of our model results with observations shows that, on the global scale, observed surface concentrations of CH3I can be best explained by the photochemical production pathway. Our results further emphasize that correlations between CH3I and abiotic or biotic factors do not necessarily provide meaningful insights concerning the source of origin. Overall, we find a net global annual CH3I air-sea flux that ranges between 70 and 260 Gg/yr. On the global scale, the ocean acts as a net source of methyl iodide for the atmosphere, though in some regions in boreal winter, fluxes are of the opposite direction (from the atmosphere to the ocean).

Digital elevation model (DEM), surveyed stream network (SSN), mask and stream sampling points (SSP)

These data are provided to allow users for reproducibility of an open source tool entitled 'automated Accumulation Threshold computation and RIparian Corridor delineation (ATRIC)'

Stable isotope record and age model for sapropel stratigraphy of ODP Site 160-967

Stable oxygen isotope data from four holes drilled at the Ocean Drilling Program Site 967, which is located on the lower northern slope of the Eratosthenes Seamount, provide a continuous record of Eastern Mediterranean surface-water conditions during the last 3.2 Ma. A high-resolution stratigraphy for the Pliocene-Pleistocene sequence was established by using a combination of astronomical calibration of sedimentary cycles, nannofossil stratigraphy, and stable oxygen isotope fluctuations. Sapropels and color cycles are present throughout the last 3.2 Ma at Site 967, and their ages, as determined by calibration against the precessional component of the astronomical record, are consistent with those estimated for the sapropels of the classical land-based marine sequences of the Punta Piccola, San Nicola, Singa, and Vrica sections (southern Italy). The Site 967 oxygen isotope record shows large amplitude fluctuations mainly caused by variations in surface water salinity throughout the entire period. Spectral analysis shows that fluctuations in the d18O record were predominantly influenced by orbital obliquity and precessional forcing from 3.2 to 1 Ma, and all main orbital frequencies characterize the d18O record for the last million years. The start of sapropel formation at 3.2 Ma indicates a possible link between sapropel formation and the build up of northern hemisphere ice sheets. The dominance of the obliquity cycle in the interval from 3.2-1 Ma further points to the sensitivity of Eastern Mediterranean climate to the fluctuations in the volume of Arctic ice sheets. An intensification of negative isotope anomalies at Site 967, relative to the open ocean, supports a link between high run-off (during warm periods) and sapropel formation. freshwater input would have inhibited deep-water formation, which led to stagnation of deeper waters. Comparison with the land sections also confirms that differential preservation and diagenesis play a key role in sapropel occurrence.

Age model of ODP Site 113-697 (Table 2)

We examined diatom preservation patterns in Pliocene age sediments of Jane Basin (ODP Site 697) and compared them with diatom distribution in more northerly sites at various sectors of the Southern Ocean. Our data from Site 697, as well as other sites from around the Southern Ocean, support the view that there was significant ice growth on Antarctica during the late Pliocene. DSDP Site 514 in the Atlantic sector shows increased relative abundance of Eucampia antarctica, an ice-related form, in the upper part of the Gauss Chron with a larger increase just above it. With one exception, all sites included in the present study show increased relative abundance of E. antarctica in the upper part of the Gauss. Our view that there was ice growth on Antarctica during the late Gauss Chron is supported by the results from ODP Site 697. While diatoms are present and percent opal is high in the early and middle Gauss Chron (suggesting more open-ocean conditions), late Gauss sediments contain low percentages of opal and few or no diatoms. This is also true for the early Matuyama Chron. If we accept spring and summer sea-ice cover as the major suppressant of diatom productivity in the Southern Ocean, then we conclude that sea-ice covered the region around Site 697 through much of the year during this interval. Further, the absence of diatoms and the low percentages of opal in middle and late Matuyama chron sediments suggests increased sea-ice cover over the Jane Basin during this time. Although warmer openocean intervals are inferred for intervals near the Olduvai and Jaramillo Subchrons, most of the Matuyama Chron was marked by extensive sea-ice cover with low seasonal contrast. Our results for the early part of the Brunhes Chron are similar, at least for the Jane Basin. During this time, sea-ice cover over the basin apparently extended well into the growing season. In contrast, the later Brunhes Chron is marked by alternating open water (during the growing season) and extensive, almost year-round, sea-ice.

Surface wind stress modeled for the Panama Seaway with the Community Climate System Model 3

The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) - induced by the closure of the gateway - was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling.

Age model and pollen analysis from site GIK16867 in the northern Angola Basin

Sporomorphs and dinoflagellate cysts from site GIK16867 in the northern Angola Basin record the vegetation history of the West African forest during the last 700 ka in relation to changes in salinity and productivity of the eastern Gulf of Guinea. During most cool and cold periods, the Afromontane forest, rather than the open grass-rich dry forest, expanded to lower altitudes partly replacing the lowland rain forest of the borderlands east of the Gulf of Guinea. Except in Stage 3, when oceanic productivity was high during a period of decreased atmospheric circulation, high oceanic productivity is correlated to strong winds. The response of marine productivity in the course of a climatic cycle, however, is earlier than that of wind vigour and makes wind-stress-induced oceanic upwelling in the area less likely. Monsoon variation is well illustrated by the pollen record of increased lowland rain forest that is paired to the dinoflagellate cyst record of decreased salinity forced by increased precipitation and run-off.

Model simulation output data (CCSM3) related to uplift of Africa and Himalaya/Tibet Plateau and the influence on tropical African vegetation

Hominid evolution in the late Miocene has long been hypothesized to be linked to the retreat of the tropical rainforest in Africa. One cause for the climatic and vegetation change often considered was uplift of Africa, but also uplift of the Himalaya and the Tibetan Plateau was suggested to have impacted rainfall distribution over Africa. Recent proxy data suggest that in East Africa open grassland habitats were available to the common ancestors of hominins and apes long before their divergence and do not find evidence for a closed rainforest in the late Miocene. We used the coupled global general circulation model CCSM3 including an interactively coupled dynamic vegetation module to investigate the impact of topography on African hydro-climate and vegetation. We performed sensitivity experiments altering elevations of the Himalaya and the Tibetan Plateau as well as of East and Southern Africa. The simulations confirm the dominant impact of African topography for climate and vegetation development of the African tropics. Only a weak influence of prescribed Asian uplift on African climate could be detected. The model simulations show that rainforest coverage of Central Africa is strongly determined by the presence of elevated African topography. In East Africa, despite wetter conditions with lowered African topography, the conditions were not favorable enough to maintain a closed rainforest. A discussion of the results with respect to other model studies indicates a minor importance of vegetation-atmosphere or ocean-atmosphere feedbacks and a large dependence of the simulated vegetation response on the land surface/vegetation model.