Matteo Spagnolo

Morphometric Analysis on the Size, Shape and Areal Distribution of Glacial Cirques in the Maritime Alps (Western French‐Italian Alps)

Abstract The morphometry of 432 glacial cirques in the Maritime Alps (Western French‐Italian Alps), studied over several years of fieldwork, was analysed with the use of a geographical information system. Some of the parameters automatically evaluated from digital elevation models required an objective and relatively new definition. In particular, cirque length was measured along a line that, from the threshold midpoint, splits the cirque into two equivalent surfaces; cirque width was automatically drawn as the longest line inscribed in the cirque and perpendicular to the length line. Significant correlations were found among the different factors and parameters analysed. In particular, cirque shape analysis showed that cirques develop allometrically in the three dimensions, i.e. more in length and width than in altitudinal range. Nevertheless cirques of the Maritime Alps have a regular, almost circular shape (mean L/W value = 1.07). The correlations among length, width and area are all very high (r2= 0.8–0.9). In terms of size, cirques show a wide range in area from 0.06 to 5.2 km2 with a mean value of 0.4 km2. The largest cirques are found on SSW‐facing slopes and at high elevations. Small cirques can be found at all altitudes but all those at high elevation are part of compound cirques at the main head valleys. Most cirques (37%) are characterized by a northern aspect; NE and SW are also frequent directions.

A GIS tool for automatic calculation of glacier equilibrium-line altitudes

A toolbox for the automated calculation of glacier equilibrium-line altitudes (ELAs) using the Accumulation Area Ratio, Area-Altitude Balance Ratio, Area-Altitude and Kurowski methods is presented. These are the most commonly-used methods of ELA calculation in palaeo-glacier reconstructions. The toolbox has been coded in Python and runs in ArcGIS requiring only the reconstructed surface of the palaeo-glacier (a DEM) as input. Through fast and automatic calculation this toolbox simplifies the process of ELA determination and can successfully work both for a single glacier and for large datasets of multiple glaciers. We describe Equilibrium Line Altitude calculation methods for palaeoglaciers.We examine suitability of each method for different glacier types.We present a toolbox for automatic ELA calculation on AAR and AABR methods.Toolbox is coded in Python and runs in ArcGIS.Toolbox use and operation are described.

Ice stream motion facilitated by a shallow-deforming and accreting bed

Ice streams drain large portions of ice sheets and play a fundamental role in governing their response to atmospheric and oceanic forcing, with implications for sea-level change. The mechanisms that generate ice stream flow remain elusive. Basal sliding and/or bed deformation have been hypothesized, but ice stream beds are largely inaccessible. Here we present a comprehensive, multi-scale study of the internal structure of mega-scale glacial lineations (MSGLs) formed at the bed of a palaeo ice stream. Analyses were undertaken at macro- and microscales, using multiple techniques including X-ray tomography, thin sections and ground penetrating radar (GPR) acquisitions. Results reveal homogeneity in stratigraphy, kinematics, granulometry and petrography. The consistency of the physical and geological properties demonstrates a continuously accreting, shallow-deforming, bed and invariant basal conditions. This implies that ice stream basal motion on soft sediment beds during MSGL formation is accommodated by plastic deformation, facilitated by continuous sediment supply and an inefficient drainage system.

A geomorphological overview of glacial landforms on the Icelandic continental shelf

Abstract Please click here to download the map associated with this article. The availability of a bathymetric database that covers about 80% of the Icelandic shelf has made it possible to produce a geomorphological map of the glacial landforms. The digital elevation model of the bathymetry was analyzed as a series of shaded relief images. Trough edges, bulging trough mouths, moraines, eskers, melt water channels, streamlined bedrock and streamlined drift, mostly hitherto unmapped, distributed all around the island have been identified. Moraines are found on the shelf, within troughs and inside fjords. Streamlined landforms are always confined to the bottom of troughs. Troughs appear to have been cut by ice streams draining an ice sheet that likely covered the entire shelf. At the shelf break, most troughs terminate with contours that bulge in a convex-outwards fashion. This suggests that an ice stream eroded, transported and finally deposited large amounts of sediment at the trough mouth. Overall, the glacial morphology of the shelf highlights a radial pattern that indicates a main ice divide near the centre of Iceland.

Manual mapping of drumlins in synthetic landscapes to assess operator effectiveness

Mapped topographic features are important for understanding processes that sculpt the Earths surface. This paper presents maps that are the primary product of an exercise that brought together 27 researchers with an interest in landform mapping wherein the efficacy and causes of variation in mapping were tested using novel synthetic DEMs containing drumlins. The variation between interpreters (e.g. mapping philosophy, experience) and across the study region (e.g. woodland prevalence) opens these factors up to assessment. A priori known answers in the synthetics increase the number and strength of conclusions that may be drawn with respect to a traditional comparative study. Initial results suggest that overall detection rates are relatively low (34–40%), but reliability of mapping is higher (72–86%). The maps form a reference dataset.

Glacial history of the Maritime Alps from the Last Glacial Maximum to the Little Ice Age

Abstract A complete sequence of glacial deposits and moraines within the same valley system in the Maritime Alps, spanning from the Last Glacial Maximum (LGM) to the Little Ice Age is presented. The sequence is geomorphologically and morphostratigraphically coherent and most stadials have been chronologically constrained by their cosmogenic exposure ages, lichenometry and by correlation with radiocarbon-dated moraines in neighbouring valleys. The shape, extent and thickness of the palaeoglaciers at each stadial have also been reconstructed and their equilibrium line altitude calculated. The LGM moraine of the Gesso Basin bears a similar equilibrium line altitude and age to that of other LGM moraines across the Alps. The recognized Late-glacial stadials show strong similarities with the corresponding stadials of the central–eastern Alpine valleys, such as Gschnitz, Bühl, Daun and Egesen. The recalculation of the exposure ages of moraine boulders with a new production rate better defines the LGM (24.0 ka) and the Egesen Stadial (13.0 ka), while the Bühl Stadial (18.5 ka) is dated for the first time in the Alps. Three early Holocene glacial advances are defined and correlated to the Kartell, Kromer and Göschenen I stadials, widely recognized in other Alpine sectors. Lichenometric dates indicate a three-fold oscillation during the Little Ice Age (thirteenth, seventeenth and nineteenth centuries).

GlaRe, a GIS tool to reconstruct the 3D surface of palaeoglaciers

Glacier reconstructions are widely used in palaeoclimatic studies and this paper presents a new semi-automated method for generating glacier reconstructions: GlaRe, is a toolbox coded in Python and operating in ArcGIS. This toolbox provides tools to generate the ice thickness from the bed topography along a palaeoglacier flowline applying the standard flow law for ice, and generates the 3D surface of the palaeoglacier using multiple interpolation methods. The toolbox performance has been evaluated using two extant glaciers, an icefield and a cirque/valley glacier from which the subglacial topography is known, using the basic reconstruction routine in GlaRe. Results in terms of ice surface, ice extent and equilibrium line altitude show excellent agreement that confirms the robustness of this procedure in the reconstruction of palaeoglaciers from glacial landforms such as frontal moraines. GlaRe is a tool for palaeoglacier 3D surface reconstruction from bed topography.GlaRe is coded in Python and runs in ArcGIS as a toolbox.The theory of perfect plasticity equilibrium glacier profile is considered.Interpolation methods for glacier 3D surface creation are described and discussed.GlaRe is tested with two extant glaciers, showing a small, acceptable error.

Using UAV acquired photography and structure from motion techniques for studying glacier landforms: application to the glacial flutes at Isfallsglaciären

© 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd. Glacier and ice sheet retreat exposes freshly deglaciated terrain which often contains small-scale fragile geomorphological features which could provide insight into subglacial or submarginal processes. Subaerial exposure results in potentially rapid landscape modification or even disappearance of the minor-relief landforms as wind, weather, water and vegetation impact on the newly exposed surface. Ongoing retreat of many ice masses means there is a growing opportunity to obtain high resolution geospatial data from glacier forelands to aid in the understanding of recent subglacial and submarginal processes. Here we used an unmanned aerial vehicle to capture close-range aerial photography of the foreland of Isfallsglaciaren, a small polythermal glacier situated in Swedish Lapland. An orthophoto and a digital elevation model with ~2 cm horizontal resolution were created from this photography using structure from motion software. These geospatial data was used to create a geomorphological map of the foreland, documenting moraines, fans, channels and flutes. The unprecedented resolution of the data enabled us to derive morphological metrics (length, width and relief) of the smallest flutes, which is not possible with other data products normally used for glacial landform metrics mapping. The map and flute metrics compare well with previous studies, highlighting the potential of this technique for rapidly documenting glacier foreland geomorphology at an unprecedented scale and resolution. The vast majority of flutes were found to have an associated stoss-side boulder, with the remainder having a likely explanation for boulder absence (burial or erosion). Furthermore, the size of this boulder was found to strongly correlate with the width and relief of the lee-side flute. This is consistent with the lee-side cavity infill model of flute formation. Whether this model is applicable to all flutes, or multiple mechanisms are required, awaits further study.

Exploring Explanations of Subglacial Bedform Sizes Using Statistical Models

Sediments beneath modern ice sheets exert a key control on their flow, but are largely inaccessible except through geophysics or boreholes. In contrast, palaeo-ice sheet beds are accessible, and typically characterised by numerous bedforms. However, the interaction between bedforms and ice flow is poorly constrained and it is not clear how bedform sizes might reflect ice flow conditions. To better understand this link we present a first exploration of a variety of statistical models to explain the size distribution of some common subglacial bedforms (i.e., drumlins, ribbed moraine, MSGL). By considering a range of models, constructed to reflect key aspects of the physical processes, it is possible to infer that the size distributions are most effectively explained when the dynamics of ice-water-sediment interaction associated with bedform growth is fundamentally random. A ‘stochastic instability’ (SI) model, which integrates random bedform growth and shrinking through time with exponential growth, is preferred and is consistent with other observations of palaeo-bedforms and geophysical surveys of active ice sheets. Furthermore, we give a proof-of-concept demonstration that our statistical approach can bridge the gap between geomorphological observations and physical models, directly linking measurable size-frequency parameters to properties of ice sheet flow (e.g., ice velocity). Moreover, statistically developing existing models as proposed allows quantitative predictions to be made about sizes, making the models testable; a first illustration of this is given for a hypothesised repeat geophysical survey of bedforms under active ice. Thus, we further demonstrate the potential of size-frequency distributions of subglacial bedforms to assist the elucidation of subglacial processes and better constrain ice sheet models.

A New Approach for the Study of the Coast Indentation Index

Abstract The indentation index, which is the ratio between the real length of a coast and its Euclidean length, is a parameter applied to characterize rock coasts and to study their evolution. Rather than subjectively selecting two or more sectors of a rock coast, the method proposed in this paper considers analyzing the indentation index on the same coastline previously split into several adjacent tracts with equal Euclidean length. By digitizing the coastline in a GIS environment, it becomes possible to test several Euclidean length values on the same coastline, obtaining a different spatial variability of the indentation index with each trial. The best length values that maximize the spatial variability of the indentation index are those that determine an indentation index pattern characterized by high variance and low spatial autocorrelation. The spatial distribution of the indentation index can eventually be analyzed considering known littoral forces acting on the studied coast. When more than one Euclidean length value is found to maximize the variability of indentation index within the same coast, it is likely that there are one or more littoral forces acting or interacting differently at different scales.