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').

Point clouds of measurements in the Dechen Cave near Iserlohn, Germany

Modeling natural phenomena from 3D information enhances our understanding of the environment. Dense 3D point clouds are increasingly used as highly detailed input datasets. In addition to the capturing techniques of point clouds with LiDAR, low-cost sensors have been released in the last few years providing access to new research fields and facilitating 3D data acquisition for a broader range of applications. This letter presents an analysis of different speleothem features using 3D point clouds acquired with the gaming device Microsoft® Kinect. We compare the Kinect sensor with terrestrial LiDAR reference measurements using the KinFu pipeline for capturing complete 3D objects (< 4m**3). The results demonstrate the suitability of the Kinect to capture flowstone walls and to derive morphometric parameters of cave features. Although the chosen capturing strategy (KinFu) reveals a high correlation (R2=0.92) of stalagmite morphometry along the vertical object axis, a systematic overestimation (22% for radii and 44% for volume) is found. The comparison of flowstone wall datasets predominantly shows low differences (mean of 1 mm with 7 mm standard deviation) of the order of the Kinect depth precision. For both objects the major differences occur at strongly varying and curved surface structures (e.g. with fine concave parts).

Global 2009 1km MODIS (MOD35/MOD09) cloud frequency and MOD35 processing path

Identifying cloud interference in satellite-derived data is a critical step toward developing useful remotely sensed products. Most MODIS land products use a combination of the MODIS (MOD35) cloud mask and the 'internal' cloud mask of the surface reflectance product (MOD09) to mask clouds, but there has been little discussion of how these masks differ globally. We calculated global mean cloud frequency for both products, for 2009, and found that inflated proportions of observations were flagged as cloudy in the Collection 5 MOD35 product. These erroneously categorized areas were spatially and environmentally non-random and usually occurred over high-albedo land-cover types (such as grassland and savanna) in several regions around the world. Additionally, we found that spatial variability in the processing path applied in the Collection 5 MOD35 algorithm affects the likelihood of a cloudy observation by up to 20% in some areas. These factors result in abrupt transitions in recorded cloud frequency across landcover and processing-path boundaries impeding their use for fine-scale spatially contiguous modeling applications. We show that together, these artifacts have resulted in significantly decreased and spatially biased data availability for Collection 5 MOD35-derived composite MODIS land products such as land surface temperature (MOD11) and net primary productivity (MOD17). Finally, we compare our results to mean cloud frequency in the new Collection 6 MOD35 product, and find that landcover artifacts have been reduced but not eliminated. Collection 6 thus increases data availability for some regions and land cover types in MOD35-derived products but practitioners need to consider how the remaining artifacts might affect their analysis.

Ground-based electromagnetic (EM) and drill-hole ice and snow thickness and melt pond depth measurements during Polarstern cruise ARK-IX/4, ARK-XI/1, and ARK-XII

Drillhole-determined sea-ice thickness was compared with values derived remotely using a portable small-offset loop-loop steady state electromagnetic (EM) induction device during expeditions to Fram Strait and the Siberian Arctic, under typical winter and summer conditions. Simple empirical transformation equations are derived to convert measured apparent conductivity into ice thickness. Despite the extreme seasonal differences in sea-ice properties as revealed by ice core analysis, the transformation equations vary little for winter and summer. Thus, the EM induction technique operated on the ice surface in the horizontal dipole mode yields accurate results within 5 to 10% of the drillhole determined thickness over level ice in both seasons. The robustness of the induction method with respect to seasonal extremes is attributed to the low salinity of brine or meltwater filling the extensive pore space in summer. Thus, the average bulk ice conductivity for summer multiyear sea ice derived according to Archie's law amounts to 23 mS/m compared to 3 mS/m for winter conditions. These mean conductivities cause only minor differences in the EM response, as is shown by means of 1-D modeling. However, under summer conditions the range of ice conductivities is wider. Along with the widespread occurrence of surface melt ponds and freshwater lenses underneath the ice, this causes greater scatter in the apparent conductivity/ice thickness relation. This can result in higher deviations between EM-derived and drillhole determined thicknesses in summer than in winter.

Pore water and solid-phase measurements on sediment cores from the Black Sea

The surface sediments in the Black Sea are underlain by extensive deposits of iron (Fe) oxide-rich lake sediments that were deposited prior to the inflow of marine Mediterranean Sea waters ca. 9000 years ago. The subsequent downward diffusion of marine sulfate into the methane-bearing lake sediments has led to a multitude of diagenetic reactions in the sulfate-methane transition zone (SMTZ), including anaerobic oxidation of methane (AOM) with sulfate. While the sedimentary cycles of sulfur (S), methane and Fe in the SMTZ have been extensively studied, relatively little is known about the diagenetic alterations of the sediment record occurring below the SMTZ. Here we combine detailed geochemical analyses of the sediment and pore water with multicomponent diagenetic modeling to study the diagenetic alterations below the SMTZ at two sites in the western Black Sea. We focus on the dynamics of Fe, S and phosphorus (P) and demonstrate that diagenesis has strongly overprinted the sedimentary burial records of these elements. Our results show that sulfate-mediated AOM substantially enhances the downward diffusive flux of sulfide into the deep limnic deposits. During this downward sulfidization, Fe oxides, Fe carbonates and Fe phosphates (e.g. vivianite) are converted to sulfide phases, leading to an enrichment in solid phase S and the release of phosphate to the pore water. Below the sulfidization front, high concentrations of dissolved ferrous Fe (Fe2+) lead to sequestration of downward diffusing phosphate as authigenic vivianite, resulting in a transient accumulation of total P directly below the sulfidization front. Our model results further demonstrate that downward migrating sulfide becomes partly re-oxidized to sulfate due to reactions with oxidized Fe minerals, fueling a cryptic S cycle and thus stimulating slow rates of sulfate-driven AOM (~ 1-100 pmol/cm**3/d) in the sulfate-depleted limnic deposits. However, this process is unlikely to explain the observed release of dissolved Fe2+ below the SMTZ. Instead, we suggest that besides organoclastic Fe oxide reduction, AOM coupled to the reduction of Fe oxides may also provide a possible mechanism for the high concentrations of Fe2+ in the pore water at depth. Our results reveal that methane plays a key role in the diagenetic alterations of Fe, S and P records in Black Sea sediments. The downward sulfidization into the limnic deposits is enhanced through sulfate-driven AOM with sulfate and AOM with Fe oxides may provide a deep source of dissolved Fe2+ that drives the sequestration of P in vivianite below the sulfidization front.

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.

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.