Search

The glacial landscape of the Alps has fascinated generations of explorers, artists, mountaineers, and scientists with its diversity, including erosional features of all scales from high-mountain cirques to steep glacial valleys and large overdeepened basins. Using previous glacier modelling results and empirical inferences of bedrock erosion under modern glaciers, we compute a distribution of potential glacier erosion in the Alps over the last glacial cycle from 120 000 years ago to the present. Despite large uncertainties pertaining to the climate history of the Alps and unconstrained glacier erosion processes, the resulting modelled patterns of glacier erosion include persistent features. The cumulative imprint of the last glacial cycle shows a very strong localization of erosion potential with local maxima at the mouths of major Alpine valleys and some other upstream sections where glaciers are modelled to have flowed with the highest velocity. The potential erosion rates vary significantly through the glacial cycle but show paradoxically little relation to the total glacier volume. Phases of glacier advance and maximum extension see a localization of rapid potential erosion rates at low elevation, while glacier erosion at higher elevation is modelled to date from phases of less extensive glaciation. The modelled erosion rates peak during deglaciation phases, when frontal retreat results in steeper glacier surface slopes, implying that climatic conditions that result in rapid glacier erosion might be quite transient and specific. Our results depict the Alpine glacier erosion landscape as a time-transgressive patchwork, with different parts of the range corresponding to different glaciation stages and time periods. [ABSTRACT FROM AUTHOR]

Modelling paleo-glacier networks in mountain ranges on the millennial timescales requires ice flow approximations. Hybrid models calculating ice flow by combining vertical shearing (shallow ice approximation) and longitudinal stretching (shallow shelf approximation) have been applied to model paleo-glacier networks on steep terrain, yet their validity has not yet been assessed quantitatively. Moreover, hybrid models consistently yield higher ice thicknesses than Last Glacial Maximum geomorphological reconstructions in the European Alps. Here, we compare results based on the hybrid Parallel Ice Sheet Model (PISM) and the Stokes model Elmer/Ice on the Rhine Glacier, a catchment of the former European Alpine Icefield. For PISM, we also test two magnitudes of flux limitation in a scheme that reduces shearing velocities. We find that the flux limitation typically used in PISM yields significantly reduced shearing speeds and increases ice thicknesses by up to 500 m, partly explaining previous overestimations. However, reducing the ice flux limitation allows the hybrid model to minimize this mismatch and captures sliding speeds, ice thicknesses, ice extent and basal temperatures in close agreement with those obtained with the Stokes model. [ABSTRACT FROM AUTHOR]

The European Alps, the cradle of pioneering glacial studies, are one of the regions where geological markers of past glaciations are most abundant and well-studied. Such conditions make the region ideal for testing numerical glacier models based on simplified ice flow physics against field-based reconstructions and vice versa. Here, we use the Parallel Ice Sheet Model (PISM) to model the entire last glacial cycle (120-0 ka) in the Alps, using horizontal resolutions of 2 and 1 km. Climate forcing is derived using two sources: present-day climate data from WorldClim and the ERA-Interim reanalysis; time-dependent temperature offsets from multiple palaeo-climate proxies. Among the latter, only the European Project for Ice Coring in Antarctica (EPICA) ice core record yields glaciation during marine oxygen isotope stages 4 (69-62 ka) and 2 (34-18 ka). This is spatially and temporally consistent with the geological reconstructions, while the other records used result in excessive early glacial cycle ice cover and a late Last Glacial Maximum. Despite the low variability of this Antarcticbased climate forcing, our simulation depicts a highly dynamic ice sheet, showing that Alpine glaciers may have advanced many times over the foreland during the last glacial cycle. Ice flow patterns during peak glaciation are largely governed by subglacial topography but include occasional transfluences through the mountain passes. Modelled maximum ice surface is on average 861m higher than observed trimline elevations in the upper Rhône Valley, yet our simulation predicts little erosion at high elevation due to cold-based ice. Finally, despite the uniform climate forcing, differences in glacier catchment hypsometry produce a timetransgressive Last Glacial Maximum advance, with some glaciers reaching their modelled maximum extent as early as 27 ka and others as late as 21 ka. [ABSTRACT FROM AUTHOR]

After more than a century of geological research, the Cordilleran ice sheet of North America remains among the least understood in terms of its former extent, volume, and dynamics. Because of the mountainous topography on which the ice sheet formed, geological studies have often had only local or regional relevance and shown such a complexity that ice-sheet-wide spatial reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (62.2-56.9 ka) and 2 (23.2-16.9 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.7 ka). [ABSTRACT FROM AUTHOR]

The article discusses the seasonal and spatial effects of temperature variability in a mass-balance model of positive degree-day (PDD) glacier surface. It mentions that the model is a parameterization of surface melt which is presumed to be proportional to the number of PDD defines and the integral of Celsius temperature over the interval of time. An overview of modelling glaciers is also presented.

Discrimination of volcanic and glacial seismic signals remains one of the key challenges for studies on volcano-glacier interactions. Here we show that even in the absence of any active volcano, interpretation of so-called Long-Period (LP) seismic events of glacial origin is truly challenging.Discoveries of LP seismic events in volcanic, tectonic and glacial environments have led to a decades-long debate on their source mechanisms, which could be either fluid-driven rupture or slow-slip. LP events in volcanic environments are of pivotal importance for eruption forecasting because they can provide information on fluid migration beneath the volcanic edifice. LP events in environments of fast glacier flow may also correspond to the resonance of fluid-filled cracks or represent sliding. However, other sources have to be ruled out before making such an assertion. Here we present observations from three comprehensive experiments at a tidewater Bowdoin Glacier, northwest Greenland, where thousands of apparently hybrid LP events were recorded. We conducted a detailed analysis of their properties by applying advanced seismological approaches and using a novel-to-glaciology GPS-stacking technique. The events are found to be highly repetitive, year-round events that are localized near the calving front, with tide-modulated amplitudes that increased with waiting time and an occurrence rate that exhibited a temporal dependence on ice speed. The emergent onset of the events is polarized, with a mixed polarity that is primarily in the horizontal plane. Taken together, the records suggest that the source of the LP events, which is sensitive to ice speed and meltwater supply, may correspond to localized slip with fluid resonance. However, further research is necessary to fully resolve the source mechanism of these LP events since we cannot dismiss the possibility of pure fluid-filled cracking due to a complex link between basal-water pressure and ice speed. [ABSTRACT FROM AUTHOR]