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.

Dead Sea lake level for the last 260 ka

The rain regime of the Levant during the late Quaternary was controlled primarily by Mediterranean cyclonic systems associated with North Atlantic climate shifts. Lake levels in the Dead Sea basin have been robust recorders of the regional hydrology and generally indicate highstand (wet) conditions throughout glacial intervals and lowstands (dry) during interglacials. However, sporadic deposition of travertines and speleothems occurred in the Negev Desert and Arava Valley during past interglacials, suggesting intrusions of humidity from southern sources probably in association with enhanced activity of mid-latitude Red Sea synoptic troughs and/or low-latitude tropical plumes. The southerly incursions of wetness were superimposed on the long-term interglacial Levantine arid conditions, as reflected by the current prevailing hyperaridity, and could have had an important impact on human migra- tion through the Red Sea-Dead Sea corridor.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

Discharge data derived from five water level gauges and discharge measurements in the Aguima and Niaou catchment, Benin, West Africa

This paper analyzes the hydrological processes and the impact of soil properties and land use on these processes in tropical headwater catchment in the sub-humid part of Benin (West-Africa), the Aguima catchment. The presented study is integrated in the GLOWA IMPETUS project, which investigates the effects of global change on the water cycle and water availability on a regional scale in Morocco and Benin. The lack of field investigations concerning soil and surface hydrology in the Benin research area necessitates detailed field measurements including measurements of discharge, soil water dynamics, soil physical properties etc. on the local scale in order to understand the dominant runoff generation processes and its influencing factors. This is a pre-requisite to be able to forecast the effects which global change has on hydrological processes and water availability in the region. The paper gives an overview over the hydrologic measuring concept of the IMPETUS-Benin project focusing on measurements concerning the soil saturated conductivity ksat and discharge behaviour of two different sub-catchment of the Aguima catchment. The results of ksat measurements revealed that interflow is the dominant runoff process on the hillslopes of the investigated catchment. Concerning the impact of land use on the hydrological processes infiltration experiments showed that infiltration rates were reduced on cultivated land compared to natural land cover. This results in significant differences in runoff behaviour and runoff ratios while comparing natural and agricultural used catchments.

Lake level reconstructions at Laguna Potrok Aike (southern Patagonia, Argentina) for the last 49 ka

Four seismic surveys and a stratigraphic record from southernmost Patagonia (Argentina) based on 51 AMS-14C dates obtained in the framework of ICDP expedition 5022 "Potrok Aike Maar Lake Sediment Archive Drilling Project" (PASADO) provide a database to compare the 106 m composite profile from the lake centre with piston cores from the littoral and outcrops in the catchment area. Based on event correlation using distinct volcanic ash layers with unique geochemical composition and optically stimulated luminescence (OSL) dates on feldspars, sediment records are firmly linked. This approach allows to match the sediment record with water levels during the past ca. 49 ka providing evidence for lake level variations. Reconstructed lake levels were 20 m higher than today during the last Glacial until the early Holocene. With the migration of the Southern Hemispheric Westerlies over this site the lake level dropped ca. 55 m for a period of two millennia. Thereupon the water balance was more positive again causing a stepwise rise of the lake level until the maximum was reached during the Little Ice Age with a subsequent lowering since the 20th century. We suggest that the mid- to late-Holocene lake level variation is caused by intensity changes of the Southern Hemispheric Westerlies.

Properties and faunal abundances of level sea ice and sea-ice ridges in the Chukchi/Beaufort Seas and Canada Basin

The abundances and distribution of metazoan within-ice meiofauna (13 stations) and under-ice fauna (12 stations) were investigated in level sea ice and sea-ice ridges in the Chukchi/Beaufort Seas and Canada Basin in June/July 2005 using a combination of ice coring and SCUBA diving. Ice meiofauna abundance was estimated based on live counts in the bottom 30 cm of level sea ice based on triplicate ice core sampling at each location, and in individual ice chunks from ridges at four locations. Under-ice amphipods were counted in situ in replicate (N=24-65 per station) 0.25 m**2 quadrats using SCUBA to a maximum water depth of 12 m. In level sea ice, the most abundant ice meiofauna groups were Turbellaria (46%), Nematoda (35%), and Harpacticoida (19%), with overall low abundances per station that ranged from 0.0 to 10.9 ind/l (median 0.8 ind/l). In level ice, low ice algal pigment concentrations (<0.1-15.8 µg Chl a /l), low brine salinities (1.8-21.7) and flushing from the melting sea ice likely explain the low ice meiofauna concentrations. Higher abundances of Turbellaria, Nematoda and Harpacticoida also were observed in pressure ridges (0-200 ind/l, median 40 ind/l), although values were highly variable and only medians of Turbellaria were significantly higher in ridge ice than in level ice. Median abundances of under-ice amphipods at all ice types (level ice, various ice ridge structures) ranged from 8 to 114 ind/m**2 per station and mainly consisted of Apherusa glacialis (87%), Onisimus spp. (7%) and Gammarus wilkitzkii (6%). Highest amphipod abundances were observed in pressure ridges at depths >3 m where abundances were up to 42-fold higher compared with level ice. We propose that the summer ice melt impacted meiofauna and under-ice amphipod abundance and distribution through (a) flushing, and (b) enhanced salinity stress at thinner level sea ice (less than 3 m thickness). We further suggest that pressure ridges, which extend into deeper, high-salinity water, become accumulation regions for ice meiofauna and under-ice amphipods in summer. Pressure ridges thus might be crucial for faunal survival during periods of enhanced summer ice melt. Previous estimates of Arctic sea ice meiofauna and under-ice amphipods on regional and pan-Arctic scales likely underestimate abundances at least in summer because they typically do not include pressure ridges.

Total V9 rDNA information organized at the OTU level

The present data set provides a tab separated text file compressed in a zip archive. The file includes metadata for each TaraOceans V9 rDNA OTU including the following fields: md5sum = identifier of the representative (most abundant) sequence of the swarm; cid = identifier of the OTU; totab = total abundance of barcodes in this OTU; TARA_xxx = number of occurrences of barcodes in this OTU in each of the 334 samples;rtotab = total abundance of the representative barcode; pid = percentage identity of the representative barcode to the closest reference sequence from V9_PR2; lineage = taxonomic path assigned to the representative barcode ; refs = best hit reference sequence(s) with respect to the representative barcode ; taxogroup = high-taxonomic level assignation of the representative barcode. The file also includes three categories of functional annotations: (1) Chloroplast: yes, presence of permanent chloroplast; no, absence of permanent chloroplast ; NA, undetermined. (2) Symbiont (small partner): parasite, the species is a parasite; commensal, the species is a commensal; mutualist, the species is a mutualist symbiont, most often a microalgal taxon involved in photosymbiosis; no the species is not involved in a symbiosis as small partner; NA, undetermined. (3) Symbiont (host): photo, the host species relies on a mutualistic microalgal photosymbiont to survive (obligatory photosymbiosis); photo_falc, same as photo, but facultative relationship; photo_klep, the host species maintains chloroplasts from microalgal prey(s) to survive; photo_klep_falc, same as photo_klep, but facultative; Nfix, the host species must interact with a mutualistic symbiont providing N2 fixation to survive; Nfix_falc, same as Nfix, but facultative; no, the species is not involved in any mutualistic symbioses; NA, undetermined.

Benthic foraminifera and sea level analyses from IODP Hole 310-M0005D

The course of sea-level fluctuations during Termination II (TII; the penultimate deglaciation), which is critical for understanding ice-sheet dynamics and suborbital climate variability, has yet to be established. This is partly because most shallow-water sequences encompassing TII were eroded during sea-level lowstands of the last glacial period or were deposited below the present sea level. Here we report a new sequence recording sea-level changes during TII in the Pleistocene sequence at Hole M0005D (water depth: 59.63 m below sea level [mbsl]) off Tahiti, French Polynesia, which was drilled during Integrated Ocean Drilling Program Expedition 310. Lithofacies variations and stratigraphic changes in the taxonomic composition, preservation states, and intraspecific test morphology of large benthic foraminifers indicate a deepening-upward sequence in the interval from Core 310-M0005D-26R (core depth: 134 mbsl) through -16R (core depth: 106 mbsl). Reconstruction of relative sea levels, based on paleodepth estimations using large benthic foraminifers, indicated a rise in sea level of about 90 m during this interval, suggesting its correlation with one of the terminations. Assuming that this rise in sea level corresponds to that during TII, after correcting for subsidence since the time of deposition, a highstand sea-level position would be 2 ± 15 m above present sea level (masl), which is generally consistent with highstand sea-level positions in MIS 5e (4 ± 2 masl). If this rise in sea level corresponds to that during older terminations, the subsidence-corrected highstand sea-level positions (30 ± 15 masl for Termination III and 54 ± 15 masl for Termination IV) are not consistent with reported ranges of interglacial sea-level highstands (-18 to 15 masl). Therefore, the studied interval likely records the rise in sea level and associated environmental changes during TII. In particular, the intervening cored materials between the two episodes of sea-level rise found in the studied interval might record the sea-level reversal event during TII. This conclusion is consistent with U/Th ages of around 133 ka, which were obtained from slightly diagenetically altered (i.e., < 1% calcite) in situ corals in the studied interval (Core 310-M0005D-20R [core depth: 118 mbsl]). This study also suggests that our inverse approach to correlate a stratigraphic interval with an approximate time frame could be useful as an independent check on the accuracy of uranium-series dating, which has been applied extensively to fossil corals in late Quaternary sea-level studies.

Calculated sea level records and deep ocean temperatures from DSDP and ODP sites

During the mid-Pleistocene transition the dominant 41 ka periodicity of glacial cycles transitioned to a quasi-100 ka periodicity for reasons not yet known. This study investigates the potential role of deep ocean hydrography by examining oxygen isotope ratios in benthic foraminifera. Oxygen isotope records from the Atlantic, Pacific and Indian Ocean basins are separated into their ice volume and local temperature/hydrography components using a piece-wise linear transfer function and a temperature calibration. Although our method has certain limitations, the deep ocean hydrography reconstructions show that glacial deep ocean temperatures approached freezing point as the mid-Pleistocene transition progressed. Further analysis suggests that water mass reorganisation could have been responsible for these temperature changes, leading to such stable conditions in the deep ocean that some obliquity cycles were skipped until precessional forcing triggered deglaciation, creating the apparent quasi-100 ka pattern. This study supports previous work that suggests multiples of obliquity cycles dominate the quasi-100 ka glacial cycles with precession components driving deglaciations.

Response of the calcifying coccolithophore Emiliania huxleyi to low pH/high pCO2: from physiology to molecular level, 2011

The emergence of ocean acidification as a significant threat to calcifying organisms in marine ecosystems creates a pressing need to understand the physiological and molecular mechanisms by which calcification is affected by environmental parameters. We report here, for the first time, changes in gene expression induced by variations in pH/pCO2 in the widespread and abundant coccolithophore Emiliania huxleyi. Batch cultures were subjected to increased partial pressure of CO2 (pCO2; i.e. decreased pH), and the changes in expression of four functional gene classes directly or indirectly related to calcification were investigated. Increased pCO2 did not affect the calcification rate and only carbonic anhydrase transcripts exhibited a significant down-regulation. Our observation that elevated pCO2 induces only limited changes in the transcription of several transporters of calcium and bicarbonate gives new significant elements to understand cellular mechanisms underlying the early response of E. huxleyi to CO2-driven ocean acidification.

Temperature, water level and bathymetry of thermokarst lakes in the continuous permafrost zone of northern Siberia - Lena River Delta, Siberia

Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009-2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths > 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W/m**2 in summer and with heat released back into the water column at a rate of less than 1 W/m**2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Hydrochemistry, water level, and discharge of water from Samoylov Island, Lena Delta, northeastern Siberia, in 2008

Although ponds make up roughly half of the total area of surface water in permafrost landscapes, their relevance to carbon dioxide emissions on a landscape scale has, to date, remained largely unknown. We have therefore investigated the inflows and outflows of dissolved organic and inorganic carbon from lakes, ponds, and outlets on Samoylov Island, in the Lena Delta of northeastern Siberia in September 2008, together with their carbon dioxide emissions. Outgassing of carbon dioxide (CO2) from these ponds and lakes, which cover 25% of Samoylov Island, was found to account for between 74 and 81% of the calculated net landscape-scale CO2 emissions of 0.2-1.1 g C/m**2/d during September 2008, of which 28-43% was from ponds and 27-46% from lakes. The lateral export of dissolved carbon was negligible compared to the gaseous emissions due to the small volumes of runoff. The concentrations of dissolved inorganic carbon in the ponds were found to triple during freezeback, highlighting their importance for temporary carbon storage between the time of carbon production and its emission as CO2. If ponds are ignored the total summer emissions of CO2-C from water bodies of the islands within the entire Lena Delta (0.7-1.3 Tg) are underestimated by between 35 and 62%.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.