(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Downscaled land surface temperatures obtained by a particle smoothing approach, with links to NetCDF files

This dataset contains continuous time series of land surface temperature (LST) at spatial resolution of 300m around the 12 experimental sites of the PAGE21 project (grant agreement number 282700, funded by the EC seventh Framework Program theme FP7-ENV-2011). This dataset was produced from hourly LST time series at 25km scale, retrieved from SSM/I data (André et al., 2015, doi:10.1016/j.rse.2015.01.028) and downscaled to 300m using a dynamic model and a particle smoothing approach. This methodology is based on two main assumptions. First, LST spatial variability is mostly explained by land cover and soil hydric state. Second, LST is unique for a land cover class within the low resolution pixel. Given these hypotheses, this variable can be estimated using a land cover map and a physically based land surface model constrained with observations using a data assimilation process. This methodology described in Mechri et al. (2014, doi:10.1002/2013JD020354) was applied to the ORCHIDEE land surface model (Krinner et al., 2005, doi:10.1029/2003GB002199) to estimate prior values of each land cover class provided by the ESA CCI-Land Cover product (Bontemps et al., 2013) at 300m resolution . The assimilation process (particle smoother) consists in simulating ensemble of LST time series for each land cover class and for a large number of parameter sets. For each parameter set, the resulting temperatures are aggregated considering the grid fraction of each land cover and compared to the coarse observations. Miniminizing the distance between the aggregated model solutions and the observations allow us to select the simulated LST and the corresponding parameter sets which fit the observations most closely. The retained parameter sets are then duplicated and randomly perturbed before simulating the next time window. At the end, the most likely LST of each land cover class are estimated and used to reconstruct LST maps at 300m resolution using ESA CCI-Land Cover. The resulting temperature maps on which ice pixels were masked, are provided at daily time step during the nine-year analysis period (2000-2009).

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

(Table 1) Sea surface temperature reconstruction for eastern equatorial Pacific surface sediment samples

Significant uncertainties persist in the reconstruction of past sea surface temperatures in the eastern equatorial Pacific, especially regarding the amplitude of the glacial cooling and the details of the post-glacial warming. Here we present the first regional calibration of alkenone unsaturation in surface sediments versus mean annual sea surface temperatures (maSST). Based on 81 new and 48 previously published data points, it is shown that open ocean samples conform to established global regressions of Uk'37 versus maSST and that there is no systematic bias from seasonality in the production or export of alkenones, or from surface ocean nutrient concentrations or salinity. The flattening of the regression at the highest maSSTs is found to be statistically insignificant. For the near-coastal Peru upwelling zone between 11-15°S and 76-79°W, however, we corroborate earlier observations that Uk'37 SST estimates significantly over-estimate maSSTs at many sites. We posit that this is caused either by uncertainties in the determination of maSSTs in this highly dynamic environment, or by biasing of the alkenone paleothermometer toward El Niño events as postulated by Rein et al. (2005).

Mean dynamic topography of the North Indian Ocean

A satellite-only Mean Dynamic Topography (MDT) of the North Indian Ocean is estimated from the DIRR5 geoid and CNES_CLS11 mean sea surface (Schaffer et al. 2012). DIRR5 geoid is estimated from the latest release (Release 5) of GOCE gravity data according to previous studies (e.g., Johannessen et al. 2003; Raj, 2014). Note that this MDT estimated is referenced to a time period of 7 years (1993-1999). A correction data obtained from AVISO is later used to convert the MDT to a time reference of 20 years (1993-2012). More details are given in Raj (2016).

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.

Long chain diols and diol indices of surface sediments of the Iberian Atlantic margin

Long chain diols are lipids that have gained interest over the last years due to their high potential to serve as biomarkers and diol indices have been proposed to reconstruct upwelling conditions and sea surface temperature (SST). However, little is known about the sources of the diols and the mechanisms impacting their distribution. Here we studied the factors controlling diol distributions in the Iberian Atlantic margin, which is characterized by a dynamic continental shelf under the influence of upwelling of nutrient-rich cold deep waters, and fluvial input. We analyzed suspended particulate matter (SPM) of the Tagus river, marine SPM and marine surface sediments along five transects off the Iberian margin, as well as riverbank sediments and soil from the catchment area of the Tagus river. Relatively high fractional abundances of the C32 1,15-diol (normalized with respect to the 1,13- and 1,15-diols) were observed in surface sediments in front of major river mouths and this abundance correlates strongly with the BIT index, a tracer for continental input of organic carbon. Together with an even higher fractional abundance of the C32 1,15-diol in the Tagus river SPM, and the absence of long chain diols in the watershed riverbank sediments and soils, we suggest that this long chain diol is produced in-situ in the river. Further support for this hypothesis comes from the small but distinct stable carbon isotopic difference of 1.3? with the marine C28 1,13-diol. The 1,14-diols are relatively abundant in surface sediments directly along the northern part of the coast, close to the upwelling zone, suggesting that Diol Indices based on 1,14-diols would work well as upwelling tracers in this region. Strikingly, we observed a significant difference in stable carbon isotopic composition between the monounsaturated C30:1 1,14- and the saturated C28 1,14-diol (3.8±0.7 per mil), suggesting different sources, in accordance with their different distributions. In addition, the Long chain Diol Index (LDI), a proxy for sea surface temperature, was applied for the surface sediments. The results correlate well with satellite SSTs offshore but reveal a significant discrepancy with satellite-derived SSTs in front of the Tagus and Sado rivers. This suggests that river outflow might compromise the applicability of this proxy.

Age and alkenone-derived Holocene sea-surface temperature records of sediment core M40-4-SL78

This Special Issue of The Holocene contains 16 research papers based on a symposium at the 11th International Meeting of the European Union of Geosciences held in Strasbourg in April 2001. The aim of the symposium was a state-of-the-art assessment of empirical studies of postglacial marine and terrestrial climatic archives and their integration with numerical climate models. This editorial places the individual papers in the broader context of natural climate variability and anthropogenic impacts on the global climate system, regional differences in climate between maritime and continental areas, and the need for an improved theoretical basis for understanding the underlying causes of environmental change. The focus of the Special Issue is the dynamic and relatively well-understood climate of the North Atlantic and the European realm, where, in relation to the steepest offshore temperature gradient on Earth, observational data are abundant and many recent advances have been made in climate reconstruction from proxy archives. The editorial also contains a summary and overview of the papers included in the four main sections of the Special Issue, which emphasize: (1) numerical modelling experiments; (2) models of glacier buildup and equilibrium-line altitude; (3) marine and terrestrial proxy records of climatic change; and (4) multiproxy palaeoenvironmental reconstruction of a Portuguese lagoonal system.

The distribution of meiofauna in surface sediments of the Fladen Ground, North Sea

A general study of structure, biomass, and dynamic estimates on meiofauna was carried out during PREFLEX (1975) and FLEX (1976), in 117- 141 m water depth. On the basis of these data an attempt was made to estimate meiofauna production, and this is discussed in relation to the energy input from the spring phytoplankton bloom. Sampling was performed at five stations, but only the stations 1, 4, and 5 were covered by a complete series from August 1975 to July 1976. At each station, from four replicate box core samples, two were withdrawn to study the abundance, distribution, and biomass of meiofauna, the content of chloroplastic pigment equivalents (CPE), and chemical and grain size analyses. At all stations grain size fell in the range of fine sand having median diameters (MD) of < 125 µm. From station 1 to 5 an increase in MD was observed. Highest values of CPE (7.81 µg m l**-1) and organic matter (4.7 %) were obtained in June and July (1976)/ August (1975), respectively. Meiofauna abundance was fairly uniform at all stations examined. Station 1 displayed maximal numbers during the whole investigation period. The abundance per 100 cm**2 varied between 15,550 and 34,900 organisms. All meiofauna studied both in total and as separate taxa showed annual cycles of abundance. Low abundance values were recorded during early summer, and maximum values during winter. High numbers of Foraminifera were obtained for August 1975 (9,460 per 100 cm**2) and July 1976 (9,710 per 100 cm**2). From December to June the values decreased from 3,280 to 1,030 per 100 cm**2. At station 1 maximum values of meiofauna biomass were recorded ranging from 1.5 to 2.7 g DWT m**-2. The mean meiofauna dry weight amounted to 2.1 g DWT m**-2. Based on minimum production, the P/B ratio for the area of station 1 might have a mean of 1.4. Taking into consideration generation times we believe that a turnover ratio of 2 is a conservative value for the Fladen Ground meiofauna. The annual production would amount to 4.2 g DWT m**-2 yr**-1. This is 27.5 % of the energy supply during the spring phytoplankton bloom, which is channelled into the meiofauna. The hypothesis is put forward that the energetic strategy of deep offshore meiofauna differs distinctively from that of shallow inshore meiofauna. While the shallow inshore meiofauna show a relatively fast response to organic matter input, the deep offshore meiofauna reacts much more slowly, the food energy consumption seems to be spread out over a longer period.