A 5 million year reconstruction of sea level, temperature, and sea water d18O

Marine sediment records from the Oligocene and Miocene reveal clear 400,000-year (400-kyr) climate cycles related to variations in orbital eccentricity. These cycles are also observed in the Plio-Pleistocene records of the global carbon cycle. However they are absent in the Late Pleistocene ice-age record over the past 1.5 million years. Here, we present a simulation of global ice volume over the past 5 million years with a coupled system of four 3-D ice-sheet models. Our simulation shows that the 400-kyr long eccentricity cycles of Antarctica vary coherently with d13C records during the Pleistocene suggesting that they drive the long-term carbon cycle changes throughout the past 35 million years. The 400-kyr response of Antarctica is eventually suppressed by the dominant 100-kyr glacial cycles of the large ice sheets in the Northern Hemisphere (NH).

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.

Sedimentology of surface samples from the Weddell Sea, Antarctica

Mineralogical and granulometric properties of glacial-marine surface sediments of the Weddell Sea and adjoining areas were studied in order to decipher spatial variations of provenance and transport paths of terrigenous detritus from Antarctic sources. The silt fraction shows marked spatial differences in quartz contents. In the sand fractions heavy-mineral assemblages display low mineralogical maturity and are dominated by garnet, green hornblende, and various types of clinopyroxene. Cluster analysis yields distinct heavy-mineral assemblages, which can be attributed to specific source rocks of the Antarctic hinterland. The configuration of modern mineralogical provinces in the near-shore regions reflects the geological variety of the adjacent hinterland. In the distal parts of the study area, sand-sized heavy minerals are good tracers of ice-rafting. Granulometric characteristics and the distribution of heavy-mineral provinces reflect maxima of relative and absolute accumulation of ice-rafted detritus in accordance with major iceberg drift tracks in the course of the Weddell Gyre. Fine-grained and coarse-grained sediment fractions may have different origins. In the central Weddell Sea, coarse ice-rafted detritus basically derives from East Antarctic sources, while the fine-fraction is discharged from weak permanent bottom currents and/or episodic turbidity currents and shows affinities to southern Weddell Sea sources. Winnowing of quartz-rich sediments through intense bottom water formation in the southern Weddell Sea provides muddy suspensions enriched in quartz. The influence of quartz-rich suspensions moving within the Weddell Gyre contour current can be traced as far as the continental slope in the northwestern Weddell Sea. In general, the focusing of mud by currents significantly exceeds the relative and absolute contribution of ice-rafted detritus beyond the shelves of the study area.

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.

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.

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.

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.