Table 1: Technetium-99 in seawater in the West Spitsbergen Current and adjacent areas, sorted from east to west

99Tc levels were measured in seawater samples collected between 2000 and 2002 in the West Spitsbergen Current (WSC) and along the western coast of Svalbard or Spitzbergen and compared with available oceanographic 3-D modelling results for the late 1990s. Additional data from related regions are also presented in order to support the data interpretation. The seawater in the Arctic fjord Kongsfjorden on the western coast of Svalbard is influenced by the WSC, as shown by the 99Tc levels in surface water. By means of the WSC, 99Tc reaches the Eastern Fram Strait, where one branch of the WSC turns west into the East Greenland Current (EGC), and another branch continues northwards into the Arctic Ocean. Surface seawater collected in the central part of the WSC during a cruise on board the R/V "Polarstern" in the summer of 2000, showed higher levels of 99Tc than samples measured in Kongsfjorden in the spring of 2000. However, all levels measured in surface water are of the same order of magnitude. Data from sampling of deeper water in the WSC and EGC provide information pertaining to the lateral distribution of 99Tc. In all vertical profiling surveys (conducted in spring and summer), the highest levels of 99Tc were found in surface water. Comparison with oceanographic 3-D modelling indicates both significant seasonal variations in the lateral stratification of the WSC and variations with depth over shorter vertical distances. This information can be applied in sampling strategies, environmental monitoring, long-range transport of pollutants and physical oceanography.

3D models created from ROV videos and their corresponding subtransect dimensions

Underwater video transects have become a common tool for quantitative analysis of the seafloor. However a major difficulty remains in the accurate determination of the area surveyed as underwater navigation can be unreliable and image scaling does not always compensate for distortions due to perspective and topography. Depending on the camera set-up and available instruments, different methods of surface measurement are applied, which make it difficult to compare data obtained by different vehicles. 3-D modelling of the seafloor based on 2-D video data and a reference scale can be used to compute subtransect dimensions. Focussing on the length of the subtransect, the data obtained from 3-D models created with the software PhotoModeler Scanner are compared with those determined from underwater acoustic positioning (ultra short baseline, USBL) and bottom tracking (Doppler velocity log, DVL). 3-D model building and scaling was successfully conducted on all three tested set-ups and the distortion of the reference scales due to substrate roughness was identified as the main source of imprecision. Acoustic positioning was generally inaccurate and bottom tracking unreliable on rough terrain. Subtransect lengths assessed with PhotoModeler were on average 20% longer than those derived from acoustic positioning due to the higher spatial resolution and the inclusion of slope. On a high relief wall bottom tracking and 3-D modelling yielded similar results. At present, 3-D modelling is the most powerful, albeit the most time-consuming, method for accurate determination of video subtransect dimensions.

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.

Snow and ice surfaces on King George Island (Antarctic Peninsula) with high-resolution TerraSAR-X time series, with links to shapfiles and GeoTIFFs

Changes of glaciers and snow cover in polar regions affect a wide range of physical and ecosystem processes on land and in the adjacent marine environment. In this study, we investigate the potential of 11-day repeat high-resolution satellite image time series from the TerraSAR-X mission to derive glaciological and hydrological parameters on King George Island, Antarctica during the period Oct/25/2010 to Apr/19/2011. The spatial pattern and temporal evolution of snow cover extent on ice-free areas can be monitored using multi-temporal coherence images. SAR coherence is used to map glacier extent of land terminating glaciers with an average accuracy of 25 m. Multi-temporal SAR color composites identify the position of the late summer snow line at about 220 m above sea level. Glacier surface velocities are obtained from intensity feature-tracking. Surface velocities near the calving front of Fourcade Glacier were up to 1.8 ± 0.01 m/d. Using an intercept theorem based on fundamental geometric principles together with differential GPS field measurements, the ice discharge of Fourcade Glacier was estimated to 20700 ± 5500 m**3/d (corresponding to ~19 ± 5 kt/d). The rapidly changing surface conditions on King George Island and the lack of high-resolution digital elevation models for the region remain restrictions for the applicability of SAR data and the precision of derived products.

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.

Micro-granulometry of sedimentary dykes, Kathmandu basin, Nepal

Soft-sediment deformation structures have been analyzed at six sites of the Kathmandu valley. Microgranulometric study (this Supplement and Fig. 3B of Mugnier et al., Tectonophysics, 2011) reveals that silty levels (60 to 80% silt) favor the development of soft-sediment deformation structures, while sandy levels (60 to 80% sand) are passively deformed. Nonetheless well sorted sand levels (more than 80% sand) generate over-fluid pressure during compaction if located beneath a silty cap, leading to fluidization and dike development. 3-D geometry of seismites indicates a very strong horizontal shearing during their development. Using a physical approach based on soil liquefaction during horizontal acceleration, we show that the fluidization zone progressively grows down-section during the shaking, but does not exactly begin at the surface. The comparison of bed-thickness and strength/depth evolution indicates three cases: i) no soft-sediment deformation occurs for thin (few centimeters) silty beds; ii) the thickness of soft-sediment deformation above sandy beds is controlled by the lithological contrast; iii) the thickness of soft-sediment deformation depends on the shaking intensity for very thick silty beds. These 3 cases are evidenced in the Kathmandu basin. We use the 30 cm-thick soft-sediment deformation level formed during the 1833 earthquake as a reference: the 1833 earthquake rupture zone extended very close to Kathmandu, inducing there MMI IX-X damages. A 90 cm-thick sediment deformation has therefore to be induced by an event greater than MMI X. From a compilation of paleo and historic seismology studies, it is found that the great (M ~ 8.1) historical earthquakes are not characteristic of the greatest earthquakes of Himalaya; hence earthquakes greater than M ~ 8.6 occurred. Kathmandu is located above one of the asperities that laterally limits the extent of mega-earthquake ruptures and two successive catastrophic events already affected Kathmandu, in 1255 located to the west of this asperity and in ~ 1100 to the east.

Dominance of facultative and obligate oligotrophic bacteria in surface waters above Gunnerus and Astrid Ridge, Antarctic Ocean

Facultative and obligate oligotrophs have been enumerated in March/April 1990 by the MPN-method with 14C-protein hydrolysate as tracer substrate. Obligate (10-3360 cells/ml) and facultative (110-9000 cells/ml) oligotrophs revealed to be the dominant population above Gunnerus Ridge (65°30'-68°S; 31-35°E) at a depth of 25 m compared with eutrophic bacteria (5 to 260 CFU/ml). Above Astrid Ridge (65-68°S; 8-18°E), obligate (0-1100 cells/ml) and facultative oligotrophs (300-9000 cells/ml) were also abundant but not always dominant. Bacterial biomass above Gunnerus Ridge was only between 7.3 and 43.6% of particulate biomass, but biomass of bacteria above Astrid Ridge amounted from 56.9 to >100% of particulate biomass; an exception was station no. PS16/552 with only 22.2% of bacterial biomass. Ratio of bacterial biomass to particulate biomass was negatively correlated with maximal primary production, complementing the view that phytoplankton was the dominant population above Gunnerus Ridge, whereas bacteria predominated above Astrid Ridge. Eutrophic bacteria were also more abundant above Astrid Ridge, with 3 to 6380 CFU/ml. Total bacteria by acridine orange direct counts amounted from 1 x 10**4 to 34.2 x 10**4 cells/ml. Bacterial biomass above Gunnerus Ridge was 1.8 to 10.7, and above Astrid Ridge 5.7 to 13.6 mg C/m*3. Maximal primary production above Gunnerus Ridge was 4.5 to 11.0, and above Astrid Ridge 2.3 to 3.5 mg C/m**3/d.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Numerical modeling of marine radiocarbon reservoir ages during the last deglaciation

A critical problem in radiocarbon dating is the spatial and temporal variability of marine reservoir ages (MRAs). We assessed the MRA evolution during the last deglaciation by numerical modeling, applying a self-consistent iteration scheme in which an existing radiocarbon chronology (derived by Hughen et al., Quat. Sci. Rev., 25, pp. 3216-3227, 2006) was readjusted by transient, 3-D simulations of marine and atmospheric Delta14C. To estimate the uncertainties regarding the ocean ventilation during the last deglaciation, we considered various ocean overturning scenarios which are based on different climatic background states (PD: modern climate, GS: LGM climate conditions). Minimum and maximum MRAs are included in file 'MRAminmax_21-14kaBP.nc'. Three further files include MRAs according to equilibrium simulations of the preindustrial ocean (file 'C14age_preindustrial.nc'; this is an update of our results published in 2005) and of the glacial ocean (files 'C14age_spinupLGM_GS.nc' and 'C14age_spinupLGM_PD.nc').

Production, oxygen uptake, and sinking velocity of copepod fecal pellets originated from the central North Sea

Production, oxygen uptake, and sinking velocity of copepod fecal pellets egested by Temora longicornis were measured using a nanoflagellate (Rhodomonas sp.), a diatom (Thalassiosira weissflogii), or a coccolithophorid (Emiliania huxleyi) as food sources. Fecal pellet production varied between 0.8 pellets ind**-1 h**-1 and 3.8 pellets ind**-1 h**-1 and was significantly higher with T. weissflogii than with the other food sources. Average pellet size varied between 2.2 x 10**5 µm**3 and 10.0 x 10**5 µm**3. Using an oxygen microsensor, small-scale oxygen fluxes and microbial respiration rates were measured directly with a spatial resolution of 2 µm at the interface of copepod fecal pellets and the surrounding water. Averaged volume-specific respiration rates were 4.12 fmol O2 µm**-3 d**-1, 2.86 fmol O2 µm**-3 d**-1, and 0.73 fmol O2 µm**-3 d**-1 in pellets produced on Rhodomonas sp., T. weissflogii, and E. huxleyi, respectively. The average carbon-specific respiration rate was 0.15 d**-1 independent on diet (range: 0.08-0.21 d**-1). Because of ballasting of opal and calcite, sinking velocities were significantly higher for pellets produced on T. weissflogii (322 +- 169 m d**-1) and E. huxleyi (200 +- 93 m d**-1) than on Rhodomonas sp. (35 +- 29 m d**-1). Preservation of carbon was estimated to be approximately 10-fold higher in fecal pellets produced when T. longicornis was fed E. huxleyi or T. weissflogii rather than Rhodomonas sp. Our study directly demonstrates that ballast increases the sinking rate of freshly produced copepod fecal pellets but does not protect them from decomposition.

Consumption rates of the periwinkle Littorina obtusata on the brown seaweed Ascophyllum nodosum

Anti-herbivory defenses support persistence of seaweeds. Little is known, however, about temporal dynamics in the induction of grazer-deterrent seaweed traits. In two induction experiments, consumption rates of the periwinkle Littorina obtusata (L.) on the brown seaweed Ascophyllum nodosum (L.) Le Jolis were measured in 3-d intervals. Changes in palatability of directly grazed A. nodosum were tested every 3 d with feeding assays using fresh and reconstituted seaweed pieces. Likewise, assays with fresh A. nodosum assessed changes in seaweed palatability in response to water-borne cues from nearby grazed conspecifics. Consumption rates of L. obtusata varied significantly during the 27-d induction phase of each experiment. Direct grazing by L. obtusata lowered palatability of fresh and reconstituted A. nodosum pieces to conspecific grazers after 15 d as well as after 6 and 12 d, respectively. After 12, 18, and 24 d, fresh A. nodosum located downstream of L. obtusata-grazed conspecifics was significantly less palatable than A. nodosum located downstream of ungrazed conspecifics. Changes in L. obtusata consumption rates and A. nodosum palatability during both induction experiments suggest temporal variation of grazer-deterrent responses, which may complicate experimental detection of inducible anti-herbivory defenses.