New insights into Late Pleistocene glacial and postglacial history of northernmost Ungava (Canada) from Pingualuit Crater Lake sediments

The Pingualuit Crater was formed by a meteoritic impact ca. 1.4 million years ago in northernmost Ungava (Canada). Due to its geographical position near the center of successive North American ice sheets and its favorable morphometry, the Pingualuit Crater Lake (water depth = 246 m) promises to yield a unique continuous sedimentary sequence covering several glacial/interglacial cycles in the terrestrial Canadian Arctic. In this paper, we suggest the existence of a subglacial lake at least during the Last Glacial Maximum (LGM) by hydraulic potential modeling using LGM ice-surface elevation and bed topography derived from a digital elevation model. These results support the hypothesis that the bottom sediments of the Crater Lake escaped glacial erosion and may contain a long-term continental sedimentary sequence. We also present the stratigraphy of a 9 m-long core retrieved from the deep basin of the lake as well as a multiproxy reconstruction of its deglacial and postglacial history. The base of the core is formed by very dense diamicton reflecting basal melt-out environments marking the end of subglacial conditions at the coring site. The overlying finely laminated silt are related to the onset of proglacial conditions characterized by extremely low lacustrine productivity. Infra Red Stimulated Luminescence and AMS ¹⁴C dating, as well as biostratigraphic data indicate sediment mixing between recent (e.g. Holocene) and much older (pre- to mid-Wisconsinan) material reworked by glacier activity. This process prevents the precise dating of these sediments that we interpret as being deposited just before the final deglaciation of the lake. Two finer grained and organic-rich intervals reflect the inception of lacustrine productivity resulting from the cessation of glacial meltwater inputs and ice-free periods. The lower organic interval corresponds to the early postglacial period (6850–5750 cal BP) and marks the transition between proglacial and postglacial conditions during the Holocene Thermal Maximum, while the uppermost organic-rich core section represents late Holocene sediments (～4200–600 cal BP). The organic intervals are separated by a basin-scale erosive slide occurring around 4200 cal BP and likely related to 1) a seismic event due to the glacio-isostatic rebound following the last deglaciation or 2) slope instabilities associated with rapid discharge events of the lake.     

Subaqueous morphology of Lake Lucerne (Central Switzerland): implications for mass movements and glacial history

Bathymetric data available for Swiss lakes have typically only low to moderate resolution and variable quality, making them insufficient for detailed underwater geomorphological studies. This article presents results of a new bathymetric survey in perialpine Lake Lucerne using modern hydrographic equipment. A digital terrain model (DTM) of the lake floor (raster dataset with 1 m cell size) covering the Chrüztrichter and Vitznau basins documents signatures of major Holocene mass movements and relics from the glacial history of the lake. Combining the bathymetry data with reflection seismic profiles and an existing event chronology allows investigating the morphology in its geological context. Subaqueous sediment slide scars with sharp headwalls cover large areas on moderately inclined slopes. The particularly large Weggis slide complex, correlated with an historical earthquake (ad 1601), features a ~9 km long and 4–7 m high headwall and covers an area of several square kilometers. Large debris cones of prehistoric rockfalls and the deposits of recent rockfall events imaged on the almost flat basin plain document mass-movement activity on steep slopes above the lake. Six transverse moraines, visible as subaqueous ridges, as lake-floor lineaments, or only imaged on reflection seismic profiles, indicate a complex glacial-inherited morphology. As many of the documented features result from potentially catastrophic events, high-resolution bathymetry can significantly improve natural hazard assessment for lakeshore communities by extending classical hazard maps to the subaqueous domain.     

Micromechanical aspects of the shear strength of wet granular soils

This paper presents a micromechanical model for the analysis of wet granular soils at low saturation (below 30%). The discrete element method is employed to model the solid particles. The capillary water is assumed to be in a pendular state and thus exists in the form of liquid bridges at the particle‐to‐particle contacts. The resulting inter‐particle adhesion is accounted for using the toroidal approximation of the bridge. Hydraulic hysteresis is accounted for based on the possible mechanism of the formation and breakage of the liquid bridges during wetting and drying phases. Shear test computational simulations were conducted at different water contents under relatively low net normal stresses. The results of these simulations suggest that capillary‐induced attractive forces and hydraulic hysteresis play an important role in affecting the shear strength of the soil. These attractive forces produce a tensile stress that contributes to the apparent cohesion of the soil and increases its stiffness. During a drying phase, capillary‐induced tensile stresses, and hence shear strength, tend to be larger than those during a wetting phase. The proposed model appears to capture the macroscopic response of wet granular materials and revealed a number of salient micromechanical mechanisms and response patterns consistent with theoretical considerations. Copyright © 2008 John Wiley & Sons, Ltd.     

<Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">V</Emphasis>–<Emphasis Type="Italic">T</Emphasis> equation of state of CaAl<Subscript>4</Subscript>Si<Subscript>2</Subscript>O<Subscript>11</Subscript> CAS phase

The thermoelastic parameters of the CAS phase (CaAl₄Si₂O₁₁) were examined by in situ high-pressure (up to 23.7 GPa) and high-temperature (up to 2,100 K) synchrotron X-ray diffraction, using a Kawai-type multi-anvil press. P–V data at room temperature fitted to a third-order Birch–Murnaghan equation of state (BM EOS) yielded: V _0,300 = 324.2 ± 0.2 Å³ and K _0,300 = 164 ± 6 GPa for K′ _0,300 = 6.2 ± 0.8. With K′ _0,300 fixed to 4.0, we obtained: V _0,300 = 324.0 ± 0.1 Å³ and K _0,300 = 180 ± 1 GPa. Fitting our P–V–T data with a modified high-temperature BM EOS, we obtained: V _0,300 = 324.2 ± 0.1 Å³, K _0,300 = 171 ± 5 GPa, K′ _0,300 = 5.1 ± 0.6 (?K _0,T /?T) _P  = ?0.023 ± 0.006 GPa K^?1, and α_0,T  = 3.09 ± 0.25 × 10^?5 K^?1. Using the equation of state parameters of the CAS phase determined in the present study, we calculated a density profile of a hypothetical continental crust that would contain ~10 vol% of CaAl₄Si₂O₁₁. Because of the higher density compared with the coexisting minerals, the CAS phase is expected to be a plunging agent for continental crust subducted in the transition zone. On the other hand, because of the lower density compared with lower mantle minerals, the CAS phase is expected to remain buoyant in the lowermost part of the transition zone.     

Exploring controls of the early and stepped deglaciation on the western margin of the British Irish Ice Sheet

New optically stimulated luminescence dating and Bayesian models integrating all legacy and BRITICE-CHRONO geochronology facilitated exploration of the controls on the deglaciation of two former sectors of the British–Irish Ice Sheet, the Donegal Bay (DBIS) and Malin Sea ice-streams (MSIS). Shelf-edge glaciation occurred ~27 ka, before the global Last Glacial Maximum, and shelf-wide retreat began 26–26.5 ka at a rate of ~18.7–20.7 m a^–1. MSIS grounding zone wedges and DBIS recessional moraines show episodic retreat punctuated by prolonged still-stands. By ~23–22 ka the outer shelf (~25 000 km²) was free of grounded ice. After this time, MSIS retreat was faster (~20 m a^–1 vs. ~2–6 m a^–1 of DBIS). Separation of Irish and Scottish ice sources occurred ~20–19.5 ka, leaving an autonomous Donegal ice dome. Inner Malin shelf deglaciation followed the submarine troughs reaching the Hebridean coast ~19 ka. DBIS retreat formed the extensive complex of moraines in outer Donegal Bay at 20.5–19 ka. DBIS retreated on land by ~17–16 ka. Isolated ice caps in Scotland and Ireland persisted until ~14.5 ka. Early retreat of this marine-terminating margin is best explained by local ice loading increasing water depths and promoting calving ice losses rather than by changes in global temperatures. Topographical controls governed the differences between the ice-stream retreat from mid-shelf to the coast.     

Increasing concentrations of arsenic and vanadium in (southern) Swedish streams

The aim of this study was to investigate temporal trends and controlling factors of As and V in running waters throughout Sweden. For this purpose, data on stream water chemistry from 62 streams of varying catchment size and characteristics, included in the Swedish environmental monitoring programmes were evaluated. The geochemical software Visual MINTEQ was used to model the speciation and trend analyses were performed on total concentrations of As and V as well as modelled fractions (dissolved species as well as arsenate and vanadate adsorbed to ferrihydrite). The trend analyses showed increasing total concentrations of As and V in southern Sweden. Concentrations of As and V correlated significantly to Fe concentrations in 59 and 60 of the 62 streams respectively, indicating that Fe is an important determining factor for As and V concentrations in Swedish streams. This was confirmed by the geochemical modelling that indicated that the adsorbed fraction is the dominant form of As and V and that the concentrations of As and V in Swedish streams are thus highly determined by concentrations of colloidal or particulate Fe. It is therefore suggested that the increasing trends of As and V are to a large extent due to increasing concentrations of colloidal Fe, which is stabilised by increasing concentrations of DOC. Further the geochemical modelling indicates that the dissolved fraction of As and V generally is small, with the exception of a few streams with high pH and/or phosphate concentrations.     

Recuperated optically stimulated luminescence dating of middle-size quartz grains from the Palaeolithic site of Bonneval (Eure-et-Loir,France)

The site of Bonneval (France) exhibits a 6 m thick stratigraphical sequence in which Middle and Lower Palaeolithic artefacts were discovered. The stratigraphy of the site consists of an alternation of loess-like silt layers deposited along glacial stages, and paleosoil horizons corresponding to the vegetal stabilization of the area during interglacials. These layers hence provide past environmental informations and can be used as stratigraphical markers and compared with well dated sequences of Northern France. Eight sediment samples from this Middle and Upper Pleistocene sequence were dated using the recuperated optically stimulated luminescence (Re-OSL) signal of quartz grains. The middle-size fraction (40–63 μm) shows a good reproducibility of the OSL signals which allows the application of a multiple-aliquot additive-dose approach for D_e determination. For all samples, the growth curves are linear up to an accumulated dose of 3000 Gy. The Re-OSL ages range from ～100 to ～450 ka and are in good agreement with the independent chronological information; they establish when handaxe industries have been made during the Middle and Upper Pleistocene in the Paris Basin.     

Application-driven ground motion prediction equation for seismic hazard assessments in non-cratonic moderate-seismicity areas

We present a ground motion prediction equation (GMPE) for probabilistic seismic hazard assessments (PSHA) in low-to-moderate seismicity areas, such as Germany. Starting from the NGA-West2 flat-file (Ancheta et al. in Earthquake Spectra 30:989–1005, 2014), we develop a model tailored to the hazard application in terms of data selection and implemented functional form. In light of such hazard application, the GMPE is derived for hypocentral distance (along with the Joyner-Boore one), selecting recordings at sites with vs30 ≥ 360 m/s, distances within 300 km, and magnitudes in the range 3 to 8 (being 7.4 the maximum magnitude for the PSHA in the target area). Moreover, the complexity of the considered functional form is reflecting the availability of information in the target area. The median predictions are compared with those from the NGA-West2 models and with one recent European model, using the Sammon’s map constructed for different scenarios. Despite the simplification in the functional form, the assessed epistemic uncertainty in the GMPE median is of the order of those affecting the NGA-West2 models for the magnitude range of interest of the hazard application. On the other hand, the simplification of the functional form led to an increment of the apparent aleatory variability. In conclusion, the GMPE developed in this study is tailored to the needs for applications in low-to-moderate seismic areas and for short return periods (e.g., 475 years); its application in studies where the hazard is involving magnitudes above 7.4 and for long return periods is not advised.     

Pattern,style and timing of British–Irish Ice Sheet retreat: Shetland and northern North Sea sector

The offshore sector around Shetland remains one of the least well-studied parts of the former British–Irish Ice Sheet with several long-standing scientific issues unresolved. These key issues include (i) the dominance of a locally sourced ‘Shetland ice cap’ vs an invasive Fennoscandian Ice Sheet; (ii) the flow configuration and style of glaciation at the Last Glacial Maximum (i.e. terrestrial vs marine glaciation); (iii) the nature of confluence between the British–Irish and Fennoscandian Ice Sheets; (iv) the cause, style and rate of ice sheet separation; and (v) the wider implications of ice sheet uncoupling on the tempo of subsequent deglaciation. As part of the Britice-Chrono project, we present new geological (seabed cores), geomorphological, marine geophysical and geochronological data from the northernmost sector of the last British–Irish Ice Sheet (north of 59.5°N) to address these questions. The study area covers ca. 95 000 km², an area approximately the size of Ireland, and includes the islands of Shetland and the surrounding continental shelf, some of the continental slope, and the western margin of the Norwegian Channel. We collect and analyse data from onshore in Shetland and along key transects offshore, to establish the most coherent picture, so far, of former ice-sheet deglaciation in this important sector. Alongside new seabed mapping and Quaternary sediment analysis, we use a multi-proxy suite of new isotopic age assessments, including 32 cosmogenic-nuclide exposure ages from glacially transported boulders and 35 radiocarbon dates from deglacial marine sediments, to develop a synoptic sector-wide reconstruction combining strong onshore and offshore geological evidence with Bayesian chronosequence modelling. The results show widespread and significant spatial fluctuations in size, shape and flow configuration of an ice sheet/ice cap centred on, or to the east of, the Orkney–Shetland Platform, between ~30 and ~15 ka BP. At its maximum extent ca. 26–25 ka BP , this ice sheet was coalescent with the Fennoscandian Ice Sheet to the east. Between ~25 and 23 ka BP the ice sheet in this sector underwent a significant size reduction from ca. 85 000 to <50 000 km², accompanied by several ice-margin oscillations. Soon after, connection was lost with the Fennoscandian Ice Sheet and a marine corridor opened to the east of Shetland. This triggered initial (and unstable) re-growth of a glaciologically independent Shetland Ice Cap ca. 21–20 ka BP with a strong east–west asymmetry with respect to topography. Ice mass growth was followed by rapid collapse, from an area of ca. 45 000 km² to ca. 15 000 km² between 19 and 18 ka BP , stabilizing at ca. 2000 km² by ~17 ka BP. Final deglaciation of Shetland occurred ca. 17–15 ka BP , and may have involved one or more subsidiary ice centres on now-submerged parts of the continental shelf. We suggest that the unusually dynamic behaviour of the northernmost sector of the British–Irish Ice Sheet between 21 and 18 ka BP – characterized by numerous extensive ice sheet/ice mass readvances, rapid loss and flow redistributions – was driven by significant changes in ice mass geometry, ice divide location and calving flux as the glaciologically independent ice cap adjusted to new boundary conditions. We propose that this dynamism was forced to a large degree by internal (glaciological) factors specific to the strongly marine-influenced Shetland Ice Cap.