Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (98)

Search Parameters:
Keywords = snow avalanche

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
33 pages, 15029 KiB  
Article
Coupling Different Machine Learning and Meta-Heuristic Optimization Techniques to Generate the Snow Avalanche Susceptibility Map in the French Alps
by Enes Can Kayhan and Ömer Ekmekcioğlu
Water 2024, 16(22), 3247; https://doi.org/10.3390/w16223247 - 12 Nov 2024
Viewed by 393
Abstract
The focus of this study is to introduce a hybrid predictive framework encompassing different meta-heuristic optimization and machine learning techniques to identify the regions susceptible to snow avalanches. To accomplish this aim, the present research sought to acquire the best-performed model among nine [...] Read more.
The focus of this study is to introduce a hybrid predictive framework encompassing different meta-heuristic optimization and machine learning techniques to identify the regions susceptible to snow avalanches. To accomplish this aim, the present research sought to acquire the best-performed model among nine different hybrid scenarios encompassing three different meta-heuristics, namely particle swarm optimization (PSO), gravitational search algorithm (GSA), and Cuckoo Search (CS), and three different ML approaches, i.e., support vector classification (SVC), stochastic gradient boosting (SGB), and k-nearest neighbors (KNN), pertaining to different predictive families. According to diligent analysis performed with regard to the blinded testing set, the PSO-SGB illustrated the most satisfactory predictive performance with an accuracy of 0.815, while the precision and recall were found to be 0.824 and 0.821, respectively. The F1-score of the predictions was found to be 0.821, and the area under the receiver operating curve (AUC) was obtained to be 0.9. Despite attaining similar predictive success via the CS-SGB model, the time-efficiency analysis underscored the PSO-SGB, as the corresponding process consumed considerably less computational time compared to its counterpart. The SHapley Additive exPlanations (SHAP) implementation further informed that slope, elevation, and wind speed are the most contributing attributes to detecting snow avalanche susceptibility in the French Alps. Full article
Show Figures

Figure 1

Figure 1
<p>Research flowchart.</p>
Full article ">Figure 2
<p>Study Domain.</p>
Full article ">Figure 3
<p>Generated layers for utilized factors. (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) LULC, (<b>g</b>) TPI, (<b>h</b>) TWI, (<b>i</b>) TRI, (<b>j</b>) lithology, (<b>k</b>) rainfall, (<b>l</b>) wind speed, (<b>m</b>) minimum temperature, (<b>n</b>) maximum temperature, (<b>o</b>) solar radiation, (<b>p</b>) snow depth, (<b>q</b>) distance to faults.</p>
Full article ">Figure 3 Cont.
<p>Generated layers for utilized factors. (<b>a</b>) elevation, (<b>b</b>) slope, (<b>c</b>) aspect, (<b>d</b>) profile curvature, (<b>e</b>) plan curvature, (<b>f</b>) LULC, (<b>g</b>) TPI, (<b>h</b>) TWI, (<b>i</b>) TRI, (<b>j</b>) lithology, (<b>k</b>) rainfall, (<b>l</b>) wind speed, (<b>m</b>) minimum temperature, (<b>n</b>) maximum temperature, (<b>o</b>) solar radiation, (<b>p</b>) snow depth, (<b>q</b>) distance to faults.</p>
Full article ">Figure 4
<p>Convergence graph of PSO with respect to the validation set (<b>a</b>) SVC, (<b>b</b>) SGB, and (<b>c</b>) KNN.</p>
Full article ">Figure 5
<p>Convergence graph of GSA with respect to the validation set (<b>a</b>) SVC, (<b>b</b>) SGB, and (<b>c</b>) KNN.</p>
Full article ">Figure 6
<p>Convergence graph of CS with respect to the validation set (<b>a</b>) SVC, (<b>b</b>) SGB, and (<b>c</b>) KNN.</p>
Full article ">Figure 7
<p>Confusion matrices for ML models with regard to the testing set (<b>a</b>) PSO-SVC, (<b>b</b>) PSO-SGB, (<b>c</b>) PSO-KNN, (<b>d</b>) GSA-SVC, (<b>e</b>) GSA-SGB, (<b>f</b>) GSA-KNN, (<b>g</b>) CS-SVC, (<b>h</b>) CS-SGB, and (<b>i</b>) CS-KNN.</p>
Full article ">Figure 8
<p>ROC plots of the ML outcomes based on the testing set (<b>a</b>) PSO, (<b>b</b>) GSA, and (<b>c</b>) CS.</p>
Full article ">Figure 9
<p>Avalanche susceptibility map for testing set based on the best-performed model.</p>
Full article ">Figure 10
<p>SHAP summary plot.</p>
Full article ">
17 pages, 2807 KiB  
Article
Anomalous Diffusion by Ocean Waves and Eddies
by Joey J. Voermans, Alexander V. Babanin, Alexei T. Skvortsov, Cagil Kirezci, Muhannad W. Gamaleldin, Henrique Rapizo, Luciano P. Pezzi, Marcelo F. Santini and Petra Heil
J. Mar. Sci. Eng. 2024, 12(11), 2036; https://doi.org/10.3390/jmse12112036 - 11 Nov 2024
Viewed by 493
Abstract
Understanding the dispersion of floating objects and ocean properties at the ocean surface is crucial for various applications, including oil spill management, debris tracking and search and rescue operations. While mesoscale turbulence has been recognized as a primary driver of dispersion, the role [...] Read more.
Understanding the dispersion of floating objects and ocean properties at the ocean surface is crucial for various applications, including oil spill management, debris tracking and search and rescue operations. While mesoscale turbulence has been recognized as a primary driver of dispersion, the role of submesoscale processes is poorly understood. This study investigates the largely unexplored mechanism of dispersion by refracted wave fields. In situ observations demonstrate significantly faster and distinct dispersion patterns for objects influenced by wind, waves and currents compared to those solely driven by ocean currents. Numerical simulations of wave fields refracted by ocean eddies corroborate these findings, revealing diffusivities that exceed those of turbulent diffusion at scales up to 10 km during energetic sea states. Our results highlight the importance of ocean waves in dispersing surface material, suggesting that refracted wave fields may play a significant role in submesoscale spreading. As atmospheric forcing at the ocean surface will only strengthen due to anthropogenic contributions, additional research into wave refraction is necessary. This requires concurrent high-resolution measurements of wind, waves and currents to inform the revisions of large-scale coupled models to better include the submesoscale physics. Full article
Show Figures

Figure 1

Figure 1
<p>Mean square separation distance <math display="inline"><semantics> <msup> <mi>s</mi> <mn>2</mn> </msup> </semantics></math> of drogued (blue) and undrogued (red) drifter pairs derived from the Global Drifter Program dataset. Shaded area gives the 95% confidence interval of the mean, determined using bootstrap algorithm (1000 samples with replacements). Only pairs with minimum separation distance of up to 1 km were considered. Geographical bias was limited by restricting the number of observations per area per unit time. Fits to the undrogued drifter data are given by <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>50</mn> <msup> <mi>t</mi> <mrow> <mn>9</mn> <mo>/</mo> <mn>5</mn> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>300</mn> <msup> <mi>t</mi> <mrow> <mn>8</mn> <mo>/</mo> <mn>5</mn> </mrow> </msup> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>5</mn> <msup> <mi>t</mi> <mn>3</mn> </msup> </mrow> </semantics></math> for the drogued drifter data, with <span class="html-italic">t</span> in days.</p>
Full article ">Figure 2
<p>Root mean square separation distance of three wave buoy clusters deployed in the Southern Ocean (red). Simplified model was fitted to the observations (green, Equation (<a href="#FD10-jmse-12-02036" class="html-disp-formula">10</a>) with <math display="inline"><semantics> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.03</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>u</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> as given in <a href="#app1-jmse-12-02036" class="html-app">Figure S2</a>). Black and gray lines represent the simulated tracer dispersion by an eddy refracted wave field, with eddy radius of 5 and 2.5 km, respectively. Four different values of the eddy velocity scale <span class="html-italic">U</span> are used, namely 0.05, 0.1, 0.2 and 0.4 m/s.</p>
Full article ">Figure 3
<p>(<b>a</b>) Pattern of refracted wave rays by an idealized ocean eddy. (<b>b</b>) Tracer particle trajectories following the Stokes drift velocity field of the refracted wave field.</p>
Full article ">Figure 4
<p>Horizontal diffusivity of drogued (blue) and undrogued (red) drifter pairs against root mean square separation distance as derived from the GDP dataset. Simulated wave-induced diffusivity of a refracted wave field is given in gray for different properties of the wave field and turbulent eddy, where the solid line represents a fit to the simulations and dashed lines are an extrapolation thereof. Best fits to GDP data when <math display="inline"><semantics> <mrow> <mn>2</mn> <mo>≤</mo> <mi>s</mi> <mo>≤</mo> <mn>20</mn> </mrow> </semantics></math> km are <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>≈</mo> <mn>1.3</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> <msup> <mi>s</mi> <mrow> <mn>1.4</mn> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>K</mi> <mo>≈</mo> <mn>7.5</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>2</mn> </mrow> </msup> <msup> <mi>s</mi> <mrow> <mn>1.0</mn> </mrow> </msup> </mrow> </semantics></math> for the drogued and undrogued drifters, respectively.</p>
Full article ">Figure 5
<p>Mean square separation distance <math display="inline"><semantics> <msup> <mi>s</mi> <mn>2</mn> </msup> </semantics></math> of drogued (blue) and undrogued (red) drifter pairs derived from the Global Drifter Program dataset. Shaded area encloses the 10th and 90th percentiles of the undrogued (red) and drogued (blue) datasets, where the top of the shaded area represents the 10th percentile and the bottom represents the 90th percentile. Only pairs with minimum separation distance of up to 1 km were considered. Geographical bias was limited by restricting the number of observations per area per unit time. Fits to the undrogued drifter data are given by <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>50</mn> <msup> <mi>t</mi> <mrow> <mn>9</mn> <mo>/</mo> <mn>5</mn> </mrow> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>300</mn> <msup> <mi>t</mi> <mrow> <mn>8</mn> <mo>/</mo> <mn>5</mn> </mrow> </msup> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>≈</mo> <mn>5</mn> <msup> <mi>t</mi> <mn>3</mn> </msup> </mrow> </semantics></math> for the drogued drifter data, with <span class="html-italic">t</span> in days.</p>
Full article ">
24 pages, 2680 KiB  
Review
Remote Sensing Techniques for Assessing Snow Avalanche Formation Factors and Building Hazard Monitoring Systems
by Natalya Denissova, Serik Nurakynov, Olga Petrova, Daniker Chepashev, Gulzhan Daumova and Alena Yelisseyeva
Atmosphere 2024, 15(11), 1343; https://doi.org/10.3390/atmos15111343 - 9 Nov 2024
Viewed by 496
Abstract
Snow avalanches, one of the most severe natural hazards in mountainous regions, pose significant risks to human lives, infrastructure, and ecosystems. As climate change accelerates shifts in snowfall and temperature patterns, it is increasingly important to improve our ability to monitor and predict [...] Read more.
Snow avalanches, one of the most severe natural hazards in mountainous regions, pose significant risks to human lives, infrastructure, and ecosystems. As climate change accelerates shifts in snowfall and temperature patterns, it is increasingly important to improve our ability to monitor and predict avalanches. This review explores the use of remote sensing technologies in understanding key geomorphological, geobotanical, and meteorological factors that contribute to avalanche formation. The primary objective is to assess how remote sensing can enhance avalanche risk assessment and monitoring systems. A systematic literature review was conducted, focusing on studies published between 2010 and 2025. The analysis involved screening relevant studies on remote sensing, avalanche dynamics, and data processing techniques. Key data sources included satellite platforms such as Sentinel-1, Sentinel-2, TerraSAR-X, and Landsat-8, combined with machine learning, data fusion, and change detection algorithms to process and interpret the data. The review found that remote sensing significantly improves avalanche monitoring by providing continuous, large-scale coverage of snowpack stability and terrain features. Optical and radar imagery enable the detection of crucial parameters like snow cover, slope, and vegetation that influence avalanche risks. However, challenges such as limitations in spatial and temporal resolution and real-time monitoring were identified. Emerging technologies, including microsatellites and hyperspectral imaging, offer potential solutions to these issues. The practical implications of these findings underscore the importance of integrating remote sensing data with ground-based observations for more robust avalanche forecasting. Enhanced real-time monitoring and data fusion techniques will improve disaster management, allowing for quicker response times and more effective policymaking to mitigate risks in avalanche-prone regions. Full article
Show Figures

Figure 1

Figure 1
<p>Flow chart of the literature search strategy.</p>
Full article ">Figure 2
<p>Geographic distribution of study areas where relevant literature was found.</p>
Full article ">Figure 3
<p>Number of publications per year.</p>
Full article ">Figure 4
<p>Word cloud illustrating the frequency of terms in titles of reviewed articles.</p>
Full article ">Figure 5
<p>Clustered co-occurrence map of most relevant terms from titles of the compiled articles.</p>
Full article ">
19 pages, 9849 KiB  
Article
Unraveling Information from Seismic Signals Generated by Gravitational Mass Movements
by Emma Suriñach and Elsa Leticia Flores-Márquez
Geosciences 2024, 14(11), 294; https://doi.org/10.3390/geosciences14110294 - 1 Nov 2024
Viewed by 489
Abstract
A practical analysis of the spectrograms of the seismic data generated by gravitational mass movements (GMMs), such as snow avalanches, landslides, lahars, and debris flows recorded on one sensor, is presented. The seismic signal produced by these movements is analyzed in terms of [...] Read more.
A practical analysis of the spectrograms of the seismic data generated by gravitational mass movements (GMMs), such as snow avalanches, landslides, lahars, and debris flows recorded on one sensor, is presented. The seismic signal produced by these movements is analyzed in terms of the shape of the initial section of the spectrogram, which corresponds to the start of the movement of the gravitational mass. The shape of the envelope of the spectrogram is a consequence of the progressive reception of high-frequency energy in the signal as the gravitational mass (GM) approaches the sensor because of the attenuation properties of the seismic waves in the ground. An exponential law was used to fit this envelope of the onset signal. The proposed methodology allows us to obtain the propagation characteristics of different types of GMM. The analysis of the adjusted parameters for different types of GMM allows us to assert that differences of one order of magnitude exist in the values of these parameters depending on the type of event. In addition, differences in the values of the exponent were obtained between the events of each type of the analyzed GMM. We present a template of different curves for each type of GMM with the corresponding parameter values that can help professionals characterize a GMM with only one seismic record (one seismic sensor) whenever the mass movement approaches the recording sensor or passes over it. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Time series (seismogram) recorded from a mass movement. (<b>b</b>) Spectrogram of the time series. The amplitude values ((m s<sup>−1</sup>)<sup>2</sup> s) are on a log<sub>10</sub> scale. SON and SOB sections are separated by a solid black line. Time (s) in horizontal axes. (<b>c</b>) Amplitude frequency transect for t = 10 s indicated as (<b>c</b>) in (<b>b</b>). (<b>d</b>) PSD amplitudes ((m s<sup>−1</sup>)<sup>2</sup> s) along the transect for the band f = 9.3000 ± 0.0156 Hz and (<b>e</b>) PSD amplitudes ((m s<sup>−1</sup>)<sup>2</sup> s) in log<sub>10</sub> scale along the transect for the band f = 9.3000 ± 0.0156 Hz indicated as (<b>e</b>) in (<b>b</b>).</p>
Full article ">Figure 2
<p>(<b>a</b>) Location B of the seismic station at the VdlS test experimental site (Valais, Switzerland). CB1 and Pra: avalanche release areas. Line and arrow: avalanche descending path. Inset: seismic station (Mark L-4C-3D sensor and Reftek–130 data acquisition system). Approximate distances: path length, 2600 m; CB1–B, 985 m. (<b>b</b>) FS seismogram and (<b>c</b>) spectrogram of a powder snow avalanche recorded at B at the VdlS experimental site (SLF, Switzerland). The horizontal axis is time. Spectrogram color scale: amplitude in (m s<sup>−1</sup>)<sup>2</sup> s with values on a log<sub>10</sub> scale. Superimposed are the exponential curves (solid black line) (Equation (7)) obtained with the corresponding Κ and <span class="html-italic">β</span> values and error margins (dashed black lines) (<a href="#geosciences-14-00294-t002" class="html-table">Table 2</a>). Enlarged image of (c,d) in <a href="#app1-geosciences-14-00294" class="html-app">Figure S3</a>.</p>
Full article ">Figure 3
<p>(<b>a</b>) Google Earth image showing the Ryggfonn experimental site (NGI, Norway). Avalanche path in yellow. In white are the approximate distances from where the avalanche seismic signals originated. Dashed line: d/d avalanches. Dotted line: d/m avalanches. The avalanche origin (shoot point (SP)) and seismometer location at TR are indicated. (<b>b</b>) Vertical component seismogram and (<b>c</b>) spectrogram of a d/m snow avalanche recorded at TR at Ryggfonn. The horizontal axis is time. Spectrogram color scale: amplitude in (m/s)<sup>2</sup> s with values on a log<sub>10</sub> scale. The exponential curves (solid white line) (Equation (7)) obtained with the corresponding Κ and <span class="html-italic">β</span> values and error margins (dashed white lines) (<a href="#geosciences-14-00294-t003" class="html-table">Table 3</a>) are superimposed.</p>
Full article ">Figure 4
<p>(<b>a</b>) Google Earth image showing the area of Laguna Beach with the location of LSS station [<a href="#B35-geosciences-14-00294" class="html-bibr">35</a>] and Landslide locations. (<b>b</b>) E-W component seismogram and (<b>c</b>) spectrogram of the Laguna Beach landslide recorded at LLS seismic station, with the inset of the first part of the avalanche. (<b>d</b>) E-W selected seismic signal (first part) (<b>e</b>) corresponding spectrogram with the exponential curve (solid black line) with the corresponding Κ and <span class="html-italic">β</span> values (Equation (7)) and the error margins superimposed (<a href="#geosciences-14-00294-t004" class="html-table">Table 4</a>). Horizontal axes are time. Spectrogram color scale: amplitude in (m s<sup>−1</sup>)<sup>2</sup> s with values on a log<sub>10</sub> scale.</p>
Full article ">Figure 5
<p>(<b>a</b>) Google Earth image showing Colima Volcano (México). In yellow, the section of the Montealegre Ravine path where the lahar seismic signals originate. The seismometer position is indicated. (<b>b</b>) Vertical component seismogram and (<b>c</b>) spectrogram of lahar 20130724 recorded at Colima Volcano (México). Horizontal axes are time. Spectrogram color scale: amplitude in (m s<sup>−1</sup>)<sup>2</sup> s with values on a log<sub>10</sub> scale. The exponential curves with the corresponding Κ and <span class="html-italic">β</span> values (Equation (7)) and the error margins (<a href="#geosciences-14-00294-t005" class="html-table">Table 5</a>) are superimposed in the spectrogram (black lines).</p>
Full article ">Figure 6
<p>(<b>a</b>) Google Earth image showing the Lattenbach debris flow catchment (Austria). In white, distance from where the debris flow seismic signals originate; geophone and infrasound sensor locations are indicated by triangles. (<b>b</b>) Vertical component seismogram, (<b>c</b>) spectrogram, (<b>d</b>) infrasound series, and (<b>e</b>) spectrogram of the 01/09/2008 debris flow in the Lattenbach catchment (Austria). Horizontal axes are time (s). Spectrogram color scale: amplitude (m s<sup>−1</sup>)<sup>2</sup> s for seismic data and in Pa<sup>2</sup> s for infrasound data with values in log<sub>10</sub> scale. The exponential curve (black line) with the corresponding Κ and <span class="html-italic">β</span> values (Equation (7)) and the error margins (black dashed lines) (<a href="#geosciences-14-00294-t006" class="html-table">Table 6</a>) are superimposed in the spectrograms.</p>
Full article ">Figure 7
<p>Flowchart describing the procedure for obtaining the template for analyzing the SON_spectrogram of the mass movement. (<b>a</b>) Theoretical SON_spectrogram generated by Equation (5). Numbers of the equations correspond to those indicated in the text. See text for the parameters used. Colors: amplitude in log<sub>10</sub> scale. (<b>b</b>) SON_spectrogram experimental data with the curve calculated with the estimated parameters K and <span class="html-italic">β</span>. The arrow lines links the steps of obtaining the theoretical and experimental SON_spectrogram.</p>
Full article ">Figure 8
<p>(<b>a</b>) SON section curves obtained from the average K and <span class="html-italic">β</span> values (Equation (7)) of <a href="#geosciences-14-00294-t007" class="html-table">Table 7</a> for the different mass movements studied. (<b>b</b>) The same curves with the origin of time at the sensor. (<b>c</b>) Inset: detail of the 50 s prior to the avalanche reaching the sensor indicated in dashed lines in (<b>b</b>). Because of the 20 s.p.s. sample rate of data, the landslide curve is plotted up to the Nyquist frequency (10 Hz).</p>
Full article ">Figure 9
<p>Example of the use of the template created with different β values (K = 1) (Equation (7)). The template is superimposed on the avalanche POW spectrogram (<a href="#geosciences-14-00294-f002" class="html-fig">Figure 2</a>). Same scales of the template and spectrogram. The template is shifted along the horizontal axis (t) to adjust the best option by eye. Black lines correspond to the curves in <a href="#geosciences-14-00294-f002" class="html-fig">Figure 2</a>. Spectrogram color scale: amplitude in (m s<sup>−1</sup>)<sup>2</sup> s with values on a log<sub>10</sub> scale.</p>
Full article ">Figure 10
<p>Determination of the parameter <span class="html-italic">β</span> of the spectrogram of the 20,120,915 lahar at the Colima Volcano (México) (<a href="#geosciences-14-00294-f005" class="html-fig">Figure 5</a>) using the template. Spectrogram color scale: amplitude in (m s<sup>−1</sup>)<sup>2</sup> s with values in log<sub>10</sub> scale. Solid lines: selected curves for the valid <span class="html-italic">β</span> value = 0.0035–0.004 s<sup>−1</sup>. Arrows indicate the direction of the sliding.</p>
Full article ">
21 pages, 18787 KiB  
Article
Snow Avalanche Susceptibility Mapping of Transportation Corridors Based on Coupled Certainty Factor and Geodetector Models
by Jie Liu, Xiliang Sun, Qiang Guo, Zhiwei Yang, Bin Wang, Senmu Yao, Haiwei Xie and Changtao Hu
Atmosphere 2024, 15(9), 1096; https://doi.org/10.3390/atmos15091096 - 9 Sep 2024
Viewed by 514
Abstract
Avalanche susceptibility assessment is a core aspect of regional avalanche early warning and risk analysis and is of great significance for disaster prevention and mitigation on proposed highways. Using sky–ground integration investigation, 83 avalanche points within the G219 Wen Quan to Horgos transportation [...] Read more.
Avalanche susceptibility assessment is a core aspect of regional avalanche early warning and risk analysis and is of great significance for disaster prevention and mitigation on proposed highways. Using sky–ground integration investigation, 83 avalanche points within the G219 Wen Quan to Horgos transportation corridor were identified, and the avalanche hazard susceptibility of the transportation corridor was partitioned using the certainty factor (CF) model and the coupled coefficient of the certainty factor–Geodetector (CF-GD) model. The CF model analysis presented nine elements of natural conditions which influence avalanche development; then, by applying the Geodetector for each of the factors, a weighting coefficient was given depending on its importance for avalanche occurrence. The results demonstrate the following: (1) According to the receiver operating characteristic (ROC) curve used to verify the accuracy, the area under the ROC curve (AUC) value for the CF-GD coupled model is 0.889, which is better than the value of 0.836 of the CF model’s evaluation accuracy, and the coupled model improves the accuracy by about 6.34% compared with the single model, indicating that the coupled model is more accurate. The results provide avalanche prevention and control recommendations for the G219 Wen Quan to Horgos transportation corridor. (2) The slope orientation, slope gradient, and mean winter temperature gradient are the main factors for avalanche development in the study area. (3) The results were validated based on the AUC values. The AUCs of the CF-GD coupled model and the CF model were 0.889 and 0.836, respectively. The accuracy of the coupled model was improved by about 6.34% compared to the single model, and the coupled CF-GD model was more accurate. The results provide avalanche control recommendations for the G219 Wen Quan to Horgos transportation corridor. Full article
Show Figures

Figure 1

Figure 1
<p>Sketch map of the study area: (<b>a</b>) map showing the location of the study area; (<b>b</b>) topographic map of the transportation corridor.</p>
Full article ">Figure 2
<p>Visit to transportation corridor to investigate avalanche distribution in 2024: (<b>a</b>) sky–ground integrated collaborative investigation; (<b>b</b>) grooved avalanche; (<b>c</b>) slope-based avalanches.</p>
Full article ">Figure 3
<p>Flowchart of the methodology.</p>
Full article ">Figure 4
<p>(<b>a</b>) Average winter snow depth; (<b>b</b>) average winter snow depth frequency ratio.</p>
Full article ">Figure 5
<p>(<b>a</b>) Average winter wind speed; (<b>b</b>) average winter wind speed frequency ratio.</p>
Full article ">Figure 6
<p>(<b>a</b>) Average winter temperature; (<b>b</b>) average winter temperature frequency ratio.</p>
Full article ">Figure 7
<p>(<b>a</b>) Surface roughness; (<b>b</b>) surface roughness frequency ratio.</p>
Full article ">Figure 8
<p>(<b>a</b>) Relief degree of land surface; (<b>b</b>) relief degree of land surface frequency ratio.</p>
Full article ">Figure 9
<p>(<b>a</b>) Surface incision; (<b>b</b>) surface incision frequency ratio.</p>
Full article ">Figure 10
<p>(<b>a</b>) Slope; (<b>b</b>) frequency ratio.</p>
Full article ">Figure 11
<p>(<b>a</b>) Aspect; (<b>b</b>) aspect frequency ratio.</p>
Full article ">Figure 12
<p>(<b>a</b>) Elevation; (<b>b</b>) elevation frequency ratio.</p>
Full article ">Figure 13
<p>Factor contribution chart.</p>
Full article ">Figure 14
<p>(<b>a</b>) Avalanche susceptibility zoning based on CF modeling; (<b>b</b>) avalanche susceptibility zoning based on a coupled CF-GD model.</p>
Full article ">Figure 15
<p>Typical avalanche map of transportation corridors.</p>
Full article ">Figure 16
<p>ROC curve analysis.</p>
Full article ">
21 pages, 9545 KiB  
Article
Universal Snow Avalanche Modeling Index Based on SAFI–Flow-R Approach in Poorly-Gauged Regions
by Uroš Durlević, Aleksandar Valjarević, Ivan Novković, Filip Vujović, Nemanja Josifov, Jelka Krušić, Blaž Komac, Tatjana Djekić, Sudhir Kumar Singh, Goran Jović, Milan Radojković and Marko Ivanović
ISPRS Int. J. Geo-Inf. 2024, 13(9), 315; https://doi.org/10.3390/ijgi13090315 - 1 Sep 2024
Viewed by 1120
Abstract
Most high-mountain regions worldwide are susceptible to snow avalanches during the winter or all year round. In this study, a Universal Snow Avalanche Modeling Index is developed, suitable for determining avalanche hazard in mountain regions. The first step in the research is the [...] Read more.
Most high-mountain regions worldwide are susceptible to snow avalanches during the winter or all year round. In this study, a Universal Snow Avalanche Modeling Index is developed, suitable for determining avalanche hazard in mountain regions. The first step in the research is the collection of data in the field and their processing in geographic information systems and remote sensing. In the period 2023–2024, avalanches were mapped in the field, and later, avalanches as points in geographic information systems (GIS) were overlapped with the dominant natural conditions in the study area. The second step involves determining the main criteria (snow cover, terrain slope, and land use) and evaluating the values to obtain the Snow Avalanche Formation Index (SAFI). Thresholds obtained through field research and the formation of avalanche inventory were used to develop the SAFI index. The index is applied with the aim of identifying locations susceptible to avalanche formation (source areas). The values used for the calculation include Normalized Difference Snow Index (NDSI > 0.6), terrain slope (20–60°) and land use (pastures, meadows). The third step presents the analysis of SAFI locations with meteorological conditions (winter precipitation and winter air temperature). The fourth step is the modeling of the propagation (simulation) of other parts of the snow avalanche in the Flow-R software 2.0. The results show that 282.9 km2 of the study area (Šar Mountains, Serbia) is susceptible to snow avalanches, with the thickness of the potentially triggered layer being 50 cm. With a 5 m thick snowpack, 299.9 km2 would be susceptible. The validation using the ROC-AUC method confirms a very high predictive power (0.94). The SAFI–Flow-R approach offers snow avalanche modeling for which no avalanche inventory is available, representing an advance for all mountain areas where historical data do not exist. The results of the study can be used for land use planning, zoning vulnerable areas, and adopting adequate environmental protection measures. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Geographical position of the Šar Mountains.</p>
Full article ">Figure 2
<p>Snow avalanches on the Šar Mountains (photo by Stanišić, M., 2023/24).</p>
Full article ">Figure 3
<p>The impact of snow avalanches on the environment of the Šar Mountains.</p>
Full article ">Figure 4
<p>Field research on the Šar Mountains.</p>
Full article ">Figure 5
<p>Maps of natural conditions and snow avalanche formation index (SAFI).</p>
Full article ">Figure 6
<p>Overlap of SAFI sites with winter precipitation and winter air temperatures.</p>
Full article ">Figure 7
<p>Combination of various datasets for the assessment of the source areas [<a href="#B51-ijgi-13-00315" class="html-bibr">51</a>].</p>
Full article ">Figure 8
<p>Flow chart with all the procedures and methods used in this research.</p>
Full article ">Figure 9
<p>The SAFI–Flow-R geospatial modeling of snow avalanches (thickness of triggered snow 0.5 m, 1 m, 3 m and 5 m).</p>
Full article ">Figure 10
<p>Settlements that are susceptible by snow avalanches.</p>
Full article ">Figure 11
<p>Validation of Snow Avalanche Modeling using ROC-AUC.</p>
Full article ">
24 pages, 32875 KiB  
Article
Integrating Sequential Backward Selection (SBS) and CatBoost for Snow Avalanche Susceptibility Mapping at Catchment Scale
by Sinem Cetinkaya and Sultan Kocaman
ISPRS Int. J. Geo-Inf. 2024, 13(9), 312; https://doi.org/10.3390/ijgi13090312 - 29 Aug 2024
Viewed by 742
Abstract
Snow avalanche susceptibility (AS) mapping is a crucial step in predicting and mitigating avalanche risks in mountainous regions. The conditioning factors used in AS modeling are diverse, and the optimal set of factors depends on the environmental and geological characteristics of the region. [...] Read more.
Snow avalanche susceptibility (AS) mapping is a crucial step in predicting and mitigating avalanche risks in mountainous regions. The conditioning factors used in AS modeling are diverse, and the optimal set of factors depends on the environmental and geological characteristics of the region. Using a sub-optimal set of input features with a data-driven machine learning (ML) method can lead to challenges like dealing with high-dimensional data, overfitting, and reduced model generalization. This study implemented a robust framework involving the Sequential Backward Selection (SBS) algorithm and a decision-tree based ML model, CatBoost, for the automatic selection of predictive variables for AS mapping. A comprehensive inventory of a large avalanche period, previously derived from satellite images, was used for the investigations in three distinct catchment areas in the Swiss Alps. The integrated SBS-CatBoost approach achieved very high classification accuracies between 94% and 97% for the three catchments. In addition, the Shapley additive explanations (SHAP) method was employed to analyze the contributions of each feature to avalanche occurrences. The proposed methodology revealed the benefits of integrating advanced feature selection algorithms with ML techniques for AS assessment. We aimed to contribute to avalanche hazard knowledge by assessing the impact of each feature in model learning. Full article
Show Figures

Figure 1

Figure 1
<p>A schematic overview of the proposed framework.</p>
Full article ">Figure 2
<p>Engelberger Aa, Meienreuss, and Göschenerreuss catchments (sub-catchments of the Reuss River Basin).</p>
Full article ">Figure 3
<p>Altitude maps of Engelberger Aa, Meienreuss, and Göschenerreuss catchments.</p>
Full article ">Figure 4
<p>LULC classification maps for Engelberger Aa, Meienreuss, and Göschenenreuss.</p>
Full article ">Figure 5
<p>Flowchart of iterative feature elimination using CatBoost with GridSearchCV.</p>
Full article ">Figure 6
<p>(<b>a</b>) The AS map of the Engelberger Aa catchment using the optimal feature set and (<b>b</b>) SHAP summary plot.</p>
Full article ">Figure 7
<p>(<b>a</b>) The AS map of the Meienreuss catchment using the optimal feature set and (<b>b</b>) SHAP summary plot.</p>
Full article ">Figure 8
<p>(<b>a</b>) The AS map of the Göschenerreuss catchment using the optimal feature set and (<b>b</b>) SHAP summary plot.</p>
Full article ">Figure A1
<p>The AS maps for Engelberger Aa, Meienreuss, and Göschenerreuss catchments, based on models with varying numbers of features (16, 15, 14, 13, and 12).</p>
Full article ">Figure A2
<p>The AS maps for Engelberger Aa, Meienreuss, and Göschenerreuss catchments, based on models with varying numbers of features (11, 10, 9, 8, and 7).</p>
Full article ">Figure A3
<p>The AS maps for Engelberger Aa, Meienreuss, and Göschenerreuss catchments, based on models with varying numbers of features (6, 5, 4, 3, and 2).</p>
Full article ">
23 pages, 8320 KiB  
Article
Validation of GPM DPR Rainfall and Drop Size Distributions Using Disdrometer Observations in the Western Mediterranean
by Eric Peinó, Joan Bech, Francesc Polls, Mireia Udina, Marco Petracca, Elisa Adirosi, Sergi Gonzalez and Brice Boudevillain
Remote Sens. 2024, 16(14), 2594; https://doi.org/10.3390/rs16142594 - 16 Jul 2024
Viewed by 1025
Abstract
Dual-frequency precipitation radar (DPR) on the Core GPM satellite provides spaceborne three-dimensional observations of precipitation fields and surface rainfall rate with quasi-global coverage. The present study evaluates the behavior of liquid precipitation intensity, radar reflectivity factor (ZKu and ZKa) and [...] Read more.
Dual-frequency precipitation radar (DPR) on the Core GPM satellite provides spaceborne three-dimensional observations of precipitation fields and surface rainfall rate with quasi-global coverage. The present study evaluates the behavior of liquid precipitation intensity, radar reflectivity factor (ZKu and ZKa) and drop size distribution (DSD) parameters (weighted mean diameter Dm and intercept parameter Nw) of the GPM DPR-derived products, version 07, from 2014 to 2023. Observations from seven Parsivel disdrometers located in different topographic zones in the Western Mediterranean are taken as ground references. Four matching techniques between satellite estimates and ground level observations were tested, and the best results were found for the so-called optimal comparison approach. Overall, GPM DPR products captured the variability of the observed DSD well at different rainfall intensities. However, overestimation of the mean Dm and underestimation of the mean Nw were observed, being much more sensitive to errors in drop diameters larger than 1.5 mm. Moreover, the lowest errors were found for radar reflectivity factor and Dm, and the highest for Nw and rainfall rate. In addition, the GPM DPR convective and stratiform classification was tested, and a substantial overestimation of stratiform cases compared to disdrometer observations were found. Full article
(This article belongs to the Special Issue Remote Sensing of Extreme Weather Events: Monitoring and Modeling)
Show Figures

Figure 1

Figure 1
<p>Digital elevation model of the region of study and the three subregions considered (mountain, plain and coast) and disdrometer sites (black dots), showing range circles of 5 km (thin red line) and 10 km (black dotted line) around each site where GPM-DPR data were collected for the present study. The lower right corner shows a map with the Köppen climate classification of the region and the disdrometer locations.</p>
Full article ">Figure 2
<p>Histograms of (<b>a</b>) Z<sub>Ku</sub>, (<b>b</b>) Z<sub>Ka</sub>, (<b>c</b>) log<sub>10</sub>(R), (<b>d</b>) D<sub>m</sub>, (<b>e</b>)10 log<sub>10</sub> (N<sub>w</sub>) and (<b>f</b>) shape parameter (µ) derived from all disdrometer and GPM DPR (DF) datasets.</p>
Full article ">Figure 3
<p>Comparison of normalized frequency distribution showing the ratio of data at each rain rate interval for disdrometer (in blue) and DPR (in red) data (top panel) and, similarly, comparison of DSD parameters D<sub>m</sub> (middle panel) and N<sub>w</sub> (bottom panel). All panels show values for the subregions plain (dashed line), coast (dotted line), mountain (semi-dashed line) and the whole region (all, thick line), for both disdrometer (DIS_) and DPR (DPR_) data.</p>
Full article ">Figure 4
<p>(<b>a</b>) Scatter density plot of raindrop size distribution measurements from all disdrometers in the D<sub>0</sub>-log(N<sub>w</sub>) space overlapped by the stratiform region (limited by turquoise dashed line) and convective region (limited by red dashed line) defined by Dolan et al. [<a href="#B48-remotesensing-16-02594" class="html-bibr">48</a>]. Disdrometer data density increase from dark to white dots, and DPR DF convective and stratiform types are indicated by red and turquoise dots, respectively. (<b>b</b>) Convective, stratiform and microphysical dominant process regions in the D<sub>0</sub>-log(N<sub>w</sub>) space according to Dolan et al. [<a href="#B48-remotesensing-16-02594" class="html-bibr">48</a>] overlapped by disdrometer (grey dots) and DF DPR (cyan dots) data.</p>
Full article ">Figure 5
<p>Scatterplots between disdrometer (x-axis) and GPM products (y-axis) considering rainfall rate (first row), Z<sub>Ka</sub> (second row), Z<sub>Ku</sub> (third row), D<sub>m</sub> (fourth row) and dBN<sub>w</sub> (fifth row) and four matching methods (point, mean 5 km, mean 10 km and optimal).</p>
Full article ">Figure 6
<p>D<sub>m</sub> vs. N<sub>w</sub> for all available disdrometer data (grey dots), GPM data of nine pixels around the disdrometers (cyan dots) and GPM data coincident with disdrometers (violet dots) under optimal method showing GPM single-frequency (<b>a</b>)- and dual-frequency (<b>b</b>)-derived estimates. The black and red dots with the error bars represent the averages and standard deviations of the disdrometer dataset and GPM 9 pixels method.</p>
Full article ">Figure 7
<p>Taylor diagram with the data obtained by the point and the optimal methods for R, Z<sub>Ka,</sub> Z<sub>Ku,</sub> D<sub>m</sub> and N<sub>w</sub> for single-frequency (<b>a</b>) and dual-frequency (<b>b</b>) GPM-derived estimates.</p>
Full article ">Figure 8
<p>(<b>a</b>) Attenuation (k) to reflectivity (Z<sub>h</sub>) ratio as a function of D<sub>m</sub> at Ku-band frequency obtained from disdrometer measurements without and with fixed shape parameter (μ = 3). (<b>b</b>) As panel (<b>a</b>) but for DFR estimated with the dual-frequency algorithm.</p>
Full article ">Figure 9
<p>Scatter density plot of log<sub>10</sub> (R) vs. D<sub>m</sub> observed by disdrometers overlayed with the relation used in the GPM DPR algorithm for ε equal to 0.2 (upper dashed line), 1.25 (solid line) and 5.0 (lower dashed line).</p>
Full article ">
21 pages, 3098 KiB  
Article
MFPANet: Multi-Scale Feature Perception and Aggregation Network for High-Resolution Snow Depth Estimation
by Liling Zhao, Junyu Chen, Muhammad Shahzad, Min Xia and Haifeng Lin
Remote Sens. 2024, 16(12), 2087; https://doi.org/10.3390/rs16122087 - 9 Jun 2024
Viewed by 827
Abstract
Accurate snow depth estimation is of significant importance, particularly for preventing avalanche disasters and predicting flood seasons. The predominant approaches for such snow depth estimation, based on deep learning methods, typically rely on passive microwave remote sensing data. However, due to the low [...] Read more.
Accurate snow depth estimation is of significant importance, particularly for preventing avalanche disasters and predicting flood seasons. The predominant approaches for such snow depth estimation, based on deep learning methods, typically rely on passive microwave remote sensing data. However, due to the low resolution of passive microwave remote sensing data, it often results in low-accuracy outcomes, posing considerable limitations in application. To further improve the accuracy of snow depth estimation, in this paper, we used active microwave remote sensing data. We fused multi-spectral optical satellite images, synthetic aperture radar (SAR) images and land cover distribution images to generate a snow remote sensing dataset (SRSD). It is a first-of-its-kind dataset that includes active microwave remote sensing images in high-latitude regions of Asia. Using these novel data, we proposed a multi-scale feature perception and aggregation neural network (MFPANet) that focuses on improving feature extraction from multi-source images. Our systematic analysis reveals that the proposed approach is not only robust but also achieves high accuracy in snow depth estimation compared to existing state-of-the-art methods, with RMSE of 0.360 and with MAE of 0.128. Finally, we selected several representative areas in our study region and applied our method to map snow depth distribution, demonstrating its broad application prospects. Full article
(This article belongs to the Special Issue Monitoring Cold-Region Water Cycles Using Remote Sensing Big Data)
Show Figures

Figure 1

Figure 1
<p>Framework of multi-scale feature perception and aggregation network.</p>
Full article ">Figure 2
<p>The structure of the residual layer in MBFE, (<b>a</b>) represents the internal structure of the residual layer without the downsampling operation, (<b>b</b>) represents the internal structure of the residual layer with the downsampling operation. <span class="html-italic">f</span> denotes the input of the residual layer, and <math display="inline"><semantics> <mrow> <msup> <mi>f</mi> <mo>′</mo> </msup> </mrow> </semantics></math> denotes the output of the residual layers.</p>
Full article ">Figure 3
<p>The structure of the MSFAA module. <span class="html-italic">f</span> denotes the input of the module, <math display="inline"><semantics> <mrow> <msup> <mi>f</mi> <mo>′</mo> </msup> </mrow> </semantics></math> denotes the output of the module.</p>
Full article ">Figure 4
<p>The structure of the HLF module.</p>
Full article ">Figure 5
<p>Location of the study area, including the Qinghai–Tibet Plateau, Xinjiang and Gansu province. The figure shows the elevation information of our study area, and the red marker is the meteorological station involved in the dataset we collected.</p>
Full article ">Figure 6
<p>Data distribution of our SRSD dataset. It is an eight-channel snow remote sensing dataset that fuses multi-source remote sensing data, including snow depth labels corresponding to 0∼42 cm.</p>
Full article ">Figure 7
<p>Multi-source data selection of this work; we randomly cropped and selected a set of data sources for display. The image on the left is a multi-spectral optical remote sensing image, the middle image is a VV single-polarization Sentinel-1 SAR image, and the image on the right displays corresponding land cover data.</p>
Full article ">Figure 8
<p>A 2D scatter diagram of measured snow depth values versus estimated by different data combinations using our method, where the green line is y = x line. The 4 channels represent the data combination of SAR + land cover, 5 channels represent the data combination of SAR + multi-spectral optical, 7 channels represent the data combination of multi-spectral optical + land cover, and our 8 channels’ data combination strategy (multi-spectral optical + SAR + Land cover) has better fitting ability by our method.</p>
Full article ">Figure 9
<p>A 2D scatter diagram of the measured snow depth values versus those estimated by different deep learning methods in the comparable studies section, which include existing snow depth estimation methods and excellent methods in image classification or image segmentation. The accuracy of snow depth estimation is also shown on the scatter plots of each method. The green line is y = x line. It is obvious that our method has better fitting ability.</p>
Full article ">Figure 10
<p>Snow depth mapping by different methods in the comparable studies section. We used purple and red matrix boxes to mark the parts with large differences in snow depth distributions in each model. We unified the snow depth scale in 0∼50 cm to observe the snow depth estimation ability of different models. According to the boxed area, we can find that the model has problems in a wide range of generalization applications when it performs similarly on our dataset; it may not be able to predict deep snow well, or it may not fit the terrain distribution well. In contrast, our proposed network can better fit terrain features and better estimate deep snow.</p>
Full article ">Figure 11
<p>Snow depth mapping by our method. (<b>a</b>,<b>b</b>) Two sets of snow depth maps show the estimation effect of our method on &gt;30 cm deep snow through local amplification of the station. (<b>c</b>,<b>d</b>) Two sets of snow depth maps show the estimation effect of our method on &lt;30 cm light snow through local amplification of the station, and the pink circle represents the location of the station. The red value denotes the measured snow depth value of the station on the day.</p>
Full article ">Figure 12
<p>Snow depth mapping by our method. (<b>a</b>,<b>b</b>) Two sets of images in this figure show the snow depth of TuoLi station area on 4 February 2015 and 27 February 2015. We can infer from the diagram that our method can perceive the changing trend in snow cover well.</p>
Full article ">
46 pages, 25578 KiB  
Review
Remote Sensing and Modeling of the Cryosphere in High Mountain Asia: A Multidisciplinary Review
by Qinghua Ye, Yuzhe Wang, Lin Liu, Linan Guo, Xueqin Zhang, Liyun Dai, Limin Zhai, Yafan Hu, Nauman Ali, Xinhui Ji, Youhua Ran, Yubao Qiu, Lijuan Shi, Tao Che, Ninglian Wang, Xin Li and Liping Zhu
Remote Sens. 2024, 16(10), 1709; https://doi.org/10.3390/rs16101709 - 11 May 2024
Cited by 1 | Viewed by 2213
Abstract
Over the past decades, the cryosphere has changed significantly in High Mountain Asia (HMA), leading to multiple natural hazards such as rock–ice avalanches, glacier collapse, debris flows, landslides, and glacial lake outburst floods (GLOFs). Monitoring cryosphere change and evaluating its hydrological effects are [...] Read more.
Over the past decades, the cryosphere has changed significantly in High Mountain Asia (HMA), leading to multiple natural hazards such as rock–ice avalanches, glacier collapse, debris flows, landslides, and glacial lake outburst floods (GLOFs). Monitoring cryosphere change and evaluating its hydrological effects are essential for studying climate change, the hydrological cycle, water resource management, and natural disaster mitigation and prevention. However, knowledge gaps, data uncertainties, and other substantial challenges limit comprehensive research in climate–cryosphere–hydrology–hazard systems. To address this, we provide an up-to-date, comprehensive, multidisciplinary review of remote sensing techniques in cryosphere studies, demonstrating primary methodologies for delineating glaciers and measuring geodetic glacier mass balance change, glacier thickness, glacier motion or ice velocity, snow extent and water equivalent, frozen ground or frozen soil, lake ice, and glacier-related hazards. The principal results and data achievements are summarized, including URL links for available products and related data platforms. We then describe the main challenges for cryosphere monitoring using satellite-based datasets. Among these challenges, the most significant limitations in accurate data inversion from remotely sensed data are attributed to the high uncertainties and inconsistent estimations due to rough terrain, the various techniques employed, data variability across the same regions (e.g., glacier mass balance change, snow depth retrieval, and the active layer thickness of frozen ground), and poor-quality optical images due to cloudy weather. The paucity of ground observations and validations with few long-term, continuous datasets also limits the utilization of satellite-based cryosphere studies and large-scale hydrological models. Lastly, we address potential breakthroughs in future studies, i.e., (1) outlining debris-covered glacier margins explicitly involving glacier areas in rough mountain shadows, (2) developing highly accurate snow depth retrieval methods by establishing a microwave emission model of snowpack in mountainous regions, (3) advancing techniques for subsurface complex freeze–thaw process observations from space, (4) filling knowledge gaps on scattering mechanisms varying with surface features (e.g., lake ice thickness and varying snow features on lake ice), and (5) improving and cross-verifying the data retrieval accuracy by combining different remote sensing techniques and physical models using machine learning methods and assimilation of multiple high-temporal-resolution datasets from multiple platforms. This comprehensive, multidisciplinary review highlights cryospheric studies incorporating spaceborne observations and hydrological models from diversified techniques/methodologies (e.g., multi-spectral optical data with thermal bands, SAR, InSAR, passive microwave, and altimetry), providing a valuable reference for what scientists have achieved in cryosphere change research and its hydrological effects on the Third Pole. Full article
(This article belongs to the Section Remote Sensing in Geology, Geomorphology and Hydrology)
Show Figures

Figure 1

Figure 1
<p>Literature on spaceborne cryosphere studies and hydrological models in HMA.</p>
Full article ">Figure 2
<p>Frequency of occurrence in the literature on spaceborne sensors for cryosphere monitoring.</p>
Full article ">Figure 3
<p>Available DEMs, surface elevation, or surface elevation difference (DH) data in HMA.</p>
Full article ">Figure 4
<p>The mean annual ground temperature (MAGT) [<a href="#B126-remotesensing-16-01709" class="html-bibr">126</a>,<a href="#B127-remotesensing-16-01709" class="html-bibr">127</a>]. (The boundary of HMA is composed of the results by Zhang [<a href="#B128-remotesensing-16-01709" class="html-bibr">128</a>], Lu [<a href="#B29-remotesensing-16-01709" class="html-bibr">29</a>], and Shean [<a href="#B59-remotesensing-16-01709" class="html-bibr">59</a>]).</p>
Full article ">Figure 5
<p>Geodetic glacier mass balance (MB in m w.e.a<sup>−1</sup>) between 2000 and 2020 in HMA with the averaged surface elevation differences by 5 km-sized hexagons from the datasets by Hugonnet et al., 2021 [<a href="#B73-remotesensing-16-01709" class="html-bibr">73</a>].</p>
Full article ">Figure 6
<p>Region—wide comparison of glacier-specific mass balance (MB, by m w.e.a<sup>−1</sup>) from five publications aggregated over three different regional boundaries in HMA (MB is marked by spots, its uncertainties are shown by the length of the bars, and different colors represents the data from the corresponding literature). (<b>a</b>) HiMAP regions [<a href="#B193-remotesensing-16-01709" class="html-bibr">193</a>]. (<b>b</b>) RGI regions [<a href="#B180-remotesensing-16-01709" class="html-bibr">180</a>]. (<b>c</b>) Regions by Kääb et al. (2015) [<a href="#B15-remotesensing-16-01709" class="html-bibr">15</a>]. (<b>d</b>) The width of the colored bars represents the periods from the five studies across HMA [<a href="#B15-remotesensing-16-01709" class="html-bibr">15</a>,<a href="#B59-remotesensing-16-01709" class="html-bibr">59</a>,<a href="#B64-remotesensing-16-01709" class="html-bibr">64</a>,<a href="#B72-remotesensing-16-01709" class="html-bibr">72</a>,<a href="#B189-remotesensing-16-01709" class="html-bibr">189</a>].</p>
Full article ">Figure 7
<p>Annual average snow water equivalent (SWE) during (<b>a</b>) 1988–2000 and (<b>b</b>) 2001–2020 (SWE is calculated from snow depth data downloaded from linkage of <a href="https://data.tpdc.ac.cn/zh-hans/data/df40346a-0202-4ed2-bb07-b65dfcda9368" target="_blank">https://data.tpdc.ac.cn/zh-hans/data/df40346a-0202-4ed2-bb07-b65dfcda9368</a> accessed on 25 February 2023).</p>
Full article ">
27 pages, 13759 KiB  
Review
Antarctic Snow Failure Mechanics: Analysis, Simulations, and Applications
by Enzhao Xiao, Shengquan Li, Ali Matin Nazar, Ronghua Zhu and Yihe Wang
Materials 2024, 17(7), 1490; https://doi.org/10.3390/ma17071490 - 25 Mar 2024
Viewed by 1205
Abstract
Snow failure is the process by which the stability of snow or snow-covered slopes is destroyed, resulting in the collapse or release of snow. Heavy snowfall, low temperatures, and volatile weather typically cause consequences in Antarctica, which can occur at different scales, from [...] Read more.
Snow failure is the process by which the stability of snow or snow-covered slopes is destroyed, resulting in the collapse or release of snow. Heavy snowfall, low temperatures, and volatile weather typically cause consequences in Antarctica, which can occur at different scales, from small, localized collapses to massive avalanches, and result in significant risk to human activities and infrastructures. Understanding snow damage is critical to assessing potential hazards associated with snow-covered terrain and implementing effective risk mitigation strategies. This review discusses the theoretical models and numerical simulation methods commonly used in Antarctic snow failure research. We focus on the various theoretical models proposed in the literature, including the fiber bundle model (FBM), discrete element model (DEM), cellular automata (CA) model, and continuous cavity-expansion penetration (CCEP) model. In addition, we overview some methods to acquire the three-dimensional solid models and the related advantages and disadvantages. Then, we discuss some critical numerical techniques used to simulate the snow failure process, such as the finite element method (FEM) and three-dimensional (3D) material point method (MPM), highlighting their features in capturing the complex behavior of snow failure. Eventually, different case studies and the experimental validation of these models and simulation methods in the context of Antarctic snow failure are presented, as well as the application of snow failure research to facility construction. This review provides a comprehensive analysis of snow properties, essential numerical simulation methods, and related applications to enhance our understanding of Antarctic snow failure, which offer valuable resources for designing and managing potential infrastructure in Antarctica. Full article
(This article belongs to the Section Mechanics of Materials)
Show Figures

Figure 1

Figure 1
<p><b>A summary of snow crystal characteristics.</b> The Nayaka diagram illustrates which snow crystal forms appear at different temperatures and supersaturation [<a href="#B38-materials-17-01490" class="html-bibr">38</a>].</p>
Full article ">Figure 2
<p>(<b>a</b>) Snow failure types’ keywords. (<b>b</b>) Authors heat picture in Snow failure research.</p>
Full article ">Figure 3
<p><b>Summary of FEM simulation method</b> [<a href="#B87-materials-17-01490" class="html-bibr">87</a>,<a href="#B89-materials-17-01490" class="html-bibr">89</a>,<a href="#B90-materials-17-01490" class="html-bibr">90</a>]. (<b>i</b>) FEM simulation flow chart and related formulate. (<b>ii</b>) Convergence analysis of Young’s modulus and density of snow sample. (<b>iii</b>) Three-dimensional representation for 3 × 3 × 3 mm<sup>3</sup> snow block, circled areas show two bonds undergoing bending deformation and histogram of stress distribution. (<b>iv</b>) Macroscopic stress and percentage of damage as function of macroscopic strain. (<b>v</b>) RVE calculations for Young’s modulus and density (right y-axis, light-grey area) vs. cube side length.</p>
Full article ">Figure 4
<p><b>A summary of DEM numerical simulation</b> [<a href="#B68-materials-17-01490" class="html-bibr">68</a>,<a href="#B70-materials-17-01490" class="html-bibr">70</a>,<a href="#B71-materials-17-01490" class="html-bibr">71</a>]. (<b>i</b>) The simulated weak snow layer (blue)—slab (gray) system and the magnified image. (<b>ii</b>) Different states of snow under pressure simulation, cohesive bonds as a function of the applied stress, and the number of broken bonds at failure as a function of the loading angle. (<b>iii</b>) The snow sample was modeled in the DEM, along with boundary conditions, and applied to the load. Individual grains are marked with varying colors for clarity. (<b>iv</b>) Three snapshots of the distribution of the local damage at different times and the final vertical profile of the damage.</p>
Full article ">Figure 5
<p><b>A summary of the FBM numerical simulation</b> [<a href="#B98-materials-17-01490" class="html-bibr">98</a>,<a href="#B99-materials-17-01490" class="html-bibr">99</a>,<a href="#B100-materials-17-01490" class="html-bibr">100</a>]. (<b>i</b>) Two models of theories of the FBM. (<b>ii</b>) Stress–strain relations of two FBM models, (a) comes from the first model in (<b>i</b>), (b) comes from the second model in (<b>ii</b>), different colors mean different sintering load. (<b>iii</b>) The ductile-to-brittle transition of the displacement-controlled FBM model. (<b>iv</b>) A snapshot of the FBM for a single load step in the form of 2D maps of the fiber bundle illustrating the sintering and load relaxation processes.</p>
Full article ">Figure 6
<p><b>Macroscopic analysis models in avalanche research.</b> (<b>a</b>) The cellular automata model simulates the flow of dense avalanches, (i) All the cylinder remains in the cell: the flow is only internal (ii) All the cylinder leaves the cell: the flow is only external. (iii) Part of the cylinder crosses the cell: there are both internal and external flows [<a href="#B111-materials-17-01490" class="html-bibr">111</a>]. (<b>b</b>) The relationship between cone penetration and cavity expansion in the CCEP model and a comparison between the bulk snow yield stress and the bulk snow effective elastic modulus obtained from the micromechanical model (MMM) and the cavity-expansion penetration model [<a href="#B112-materials-17-01490" class="html-bibr">112</a>].</p>
Full article ">Figure 7
<p><b>Macroscopic analysis model of the 3D MPM in avalanche research.</b> (<b>i</b>) Overview of the MPM algorithm. (<b>ii</b>) The three-dimensional (3D) material point method (MPM) is used to explore different avalanches on complex real terrain, and the flow characteristics of avalanches from release to deposition are comprehensively characterized. (<b>iii</b>) Four typical distribution of density ratio in avalanche [<a href="#B110-materials-17-01490" class="html-bibr">110</a>].</p>
Full article ">Figure 8
<p><b>A summary of the three-dimensional solid modeling of snow.</b> (<b>a</b>) Continuous sections and continuous micrographs are taken [<a href="#B122-materials-17-01490" class="html-bibr">122</a>]. (<b>b</b>) Examples of snow sections [<a href="#B123-materials-17-01490" class="html-bibr">123</a>]; image (A) is the partially metamorphosed snow of mid-density (−200 kg/m<sup>3</sup>), image (B) is the snow metamorphosed for 6 months at a weak temperature gradient (−1 °C/m), and image (C) is the snow metamorphosed for 100 h in an intense temperature gradient (&gt;1000 °C/m). (<b>c</b>) The SkyScan material test bench for the micro-CT system for scanning snow blocks [<a href="#B124-materials-17-01490" class="html-bibr">124</a>]. (<b>d</b>) The experimental arrangement of the micro-CT system [<a href="#B125-materials-17-01490" class="html-bibr">125</a>]. (<b>e</b>) Typical snow types presented as snow grains (left) and as surface renderings (right) [<a href="#B126-materials-17-01490" class="html-bibr">126</a>]. (<b>f</b>) Two images separated by 2 days from a time-lapse movie [<a href="#B126-materials-17-01490" class="html-bibr">126</a>], image (A) at time 577 h, image (B) at time 625 h. (<b>g</b>) Three-dimensional images of the distribution of snow grains [<a href="#B125-materials-17-01490" class="html-bibr">125</a>]. (<b>h</b>) μCT detect micro compression devices [<a href="#B127-materials-17-01490" class="html-bibr">127</a>]. (<b>i</b>) The snow grains are imaged under the microscope [<a href="#B128-materials-17-01490" class="html-bibr">128</a>]. (<b>j</b>) An overview of μCT image analysis [<a href="#B129-materials-17-01490" class="html-bibr">129</a>]; (A) is the picture of the snow vertical transect μCT and (B) is the picture of snow horizontal transect μCT. (<b>k</b>) An X-ray microtomography analysis of the isothermal densification of new snow under external mechanical stress [<a href="#B130-materials-17-01490" class="html-bibr">130</a>].</p>
Full article ">Figure 9
<p><b>A summary of the avalanche formation mechanism and release process.</b> (<b>a</b>) Different types of snow crystals [<a href="#B38-materials-17-01490" class="html-bibr">38</a>]. (<b>b</b>) Photomicrographs of metamorphic crystals frequently observed in snowpacks [<a href="#B139-materials-17-01490" class="html-bibr">139</a>]. (<b>c</b>) Snowpack crystals [<a href="#B140-materials-17-01490" class="html-bibr">140</a>]. (<b>d</b>) Wet plate avalanche [<a href="#B141-materials-17-01490" class="html-bibr">141</a>]. (<b>e</b>) The measurements of snow’s mechanical properties [<a href="#B142-materials-17-01490" class="html-bibr">142</a>]. (<b>f</b>) A schematic illustration of the applied IPCC risk concept, which is used to create risk maps and perform risk analyses [<a href="#B143-materials-17-01490" class="html-bibr">143</a>]. (<b>g</b>) Snow slab release models with pre-existing weakness [<a href="#B131-materials-17-01490" class="html-bibr">131</a>]. (<b>h</b>) Three tire–snow simulation methods. (i) F<sub>S</sub> is the shear force of snow in void (space between tread blocks), F is the frictional force between tire and snow, F<sub>µ</sub> is the digging force (edge effect generated by sipes and blocks), and F<sub>B</sub> is the braking force. (ii) Computing a continuous density field from a collection of point mass particles. (A) is the particle-based method, (B) is the way to calculate based on the size of the sample volume, and (C) is the SPH method. (iii) The DEM contact model between two contacting particles, (left) normal contact model, (right) tangential contact model [<a href="#B144-materials-17-01490" class="html-bibr">144</a>].</p>
Full article ">
10 pages, 2726 KiB  
Article
Perlite Has Similar Diffusion Properties for Oxygen and Carbon Dioxide to Snow: Implications for Avalanche Safety Equipment Testing and Breathing Studies
by Simon Walzel, Martin Rozanek and Karel Roubik
Appl. Sci. 2023, 13(23), 12569; https://doi.org/10.3390/app132312569 - 22 Nov 2023
Viewed by 970
Abstract
On average, one hundred people die each year under avalanche snow. Despite extensive global research on gas exchange in buried avalanche victims, it remains unclear how the diffusion of respiratory gases affects survival under avalanche snow. This study aims to determine how oxygen [...] Read more.
On average, one hundred people die each year under avalanche snow. Despite extensive global research on gas exchange in buried avalanche victims, it remains unclear how the diffusion of respiratory gases affects survival under avalanche snow. This study aims to determine how oxygen and carbon dioxide diffuse through snow, as well as through wet and dry perlite, which may serve as a surrogate for avalanche snow. A custom-made apparatus to study the diffusion of respiratory gases consisted of a plastic cylinder (1200 mm long, ID 300 mm) with 13 gas sampling needles evenly spaced along the axis of the cylinder filled with the tested material. Following 60 min of free diffusion, gas samples were analyzed using a vital signs monitor with a module for respiratory gas analysis (E-CAiOVX, Datex-Ohmeda, GE Healthcare, Chicago, IL, USA). A combination of 16% oxygen, 5% carbon dioxide, and 79% nitrogen was used. The rates of diffusion for both respiratory gases were comparable in snow and both forms of perlite. Oxygen propagated faster than carbon dioxide. Due to similar diffusion characteristics to snow, perlite possesses the potential to stand in as an effective substitute for soft snow for the study of respiratory dynamics, for conducting breathing experiments, and for testing avalanche safety equipment. Full article
(This article belongs to the Section Materials Science and Engineering)
Show Figures

Figure 1

Figure 1
<p>A scheme of the apparatus with the gas analyzer and gas cylinders. 1—buckle fixing the cover to the cylinder; 2—plastic mesh separating the cylinder into two parts; 3—airtight O-rings; 4—rubber plug with a metal needle; 5—three-way valve; 6—sampling point of the gas mixture in the inlet chamber; 7—gas inlet port; 8—gas outlet port; 9—pressure-reducing valve; 10—throttle valve. Dimensions are in millimeters.</p>
Full article ">Figure 2
<p>Normalized measured concentrations of O<sub>2</sub> and CO<sub>2</sub> related to the concentrations in the inlet chamber depending on the distance from the inlet chamber in S, PW, and PD. Symbol # indicates a statistically significant difference between O<sub>2</sub> and CO<sub>2</sub>.</p>
Full article ">Figure 3
<p>Measured concentrations of O<sub>2</sub> and CO<sub>2</sub> related to the concentrations in the inlet chamber depending on the distance from the inlet chamber in S, PW, and PD. Symbol # indicates a statistically significant difference between S and PD, symbol <span>$</span> indicates a statistically significant difference between PW and PD, and symbol § indicates a statistically significant difference between S and PW.</p>
Full article ">Figure 4
<p>The measured DCs for O<sub>2</sub> and CO<sub>2</sub> in snow depending on the distance from the inlet chamber.</p>
Full article ">
17 pages, 3851 KiB  
Article
Snow Avalanche Hazard Mapping Using a GIS-Based AHP Approach: A Case of Glaciers in Northern Pakistan from 2012 to 2022
by Afia Rafique, Muhammad Y. S. Dasti, Barkat Ullah, Fuad A. Awwad, Emad A. A. Ismail and Zulfiqar Ahmad Saqib
Remote Sens. 2023, 15(22), 5375; https://doi.org/10.3390/rs15225375 - 16 Nov 2023
Cited by 2 | Viewed by 2214
Abstract
Snow avalanches are a type of serious natural disaster that commonly occur in snow-covered mountains with steep terrain characteristics. Susceptibility analysis of avalanches is a pressing issue today and helps decision makers to implement appropriate avalanche risk reduction strategies. Avalanche susceptibility maps provide [...] Read more.
Snow avalanches are a type of serious natural disaster that commonly occur in snow-covered mountains with steep terrain characteristics. Susceptibility analysis of avalanches is a pressing issue today and helps decision makers to implement appropriate avalanche risk reduction strategies. Avalanche susceptibility maps provide a preliminary method for evaluating places that are likely to be vulnerable to avalanches to stop or reduce the risks of such disasters. The current study aims to identify areas that are vulnerable to avalanches (ranging from extremely high and low danger) by considering geo-morphological and geological variables and employing an Analytical Hierarchy Approach (AHP) in the GIS platform to identify potential snow avalanche zones in the Karakoram region in Northern Pakistan. The Advanced Space-borne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) was used to extract the elevation, slope, aspect, terrain roughness, and curvature of the study area. This study includes the risk identification variable of land cover (LC), which was obtained from the Landsat 8 Operational Land Imager (OLI) satellite. The obtained result showed that the approach established in this study provided a quick and reliable tool to map avalanches in the study area, and it might also work with other glacier sites in other parts of the world for snow avalanche susceptibility and risk assessments. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The occurrence of snow avalanches and damage due to snow avalanches in the study area.</p>
Full article ">Figure 2
<p>The study area map represents the location of Hispar, Gayari Sector, and Batura Glacier.</p>
Full article ">Figure 3
<p>Flowchart representation of the research methodology used.</p>
Full article ">Figure 4
<p>Thematic layers of avalanche occurrence factors at Hispar Glacier.</p>
Full article ">Figure 5
<p>Thematic layers of avalanche occurrence factors at Batura Glacier.</p>
Full article ">Figure 6
<p>Thematic layers of avalanche occurrence factors at Gayari Glacier.</p>
Full article ">Figure 7
<p>Snow avalanche susceptibility of Hispar, Gayari, and Batura Glaciers for 2012 and 2022.</p>
Full article ">
24 pages, 13409 KiB  
Article
A Study on Avalanche-Triggering Factors and Activity Characteristics in Aerxiangou, West Tianshan Mountains, China
by Jie Liu, Tianyi Zhang, Changtao Hu, Bin Wang, Zhiwei Yang, Xiliang Sun and Senmu Yao
Atmosphere 2023, 14(9), 1439; https://doi.org/10.3390/atmos14091439 - 15 Sep 2023
Cited by 3 | Viewed by 1784
Abstract
Through analyzing the triggering factors and activity characteristics of avalanches in Aerxiangou in the Western Tianshan Mountains, the formation and disaster-causing process of avalanches were studied to provide theoretical support and a scientific basis for avalanche disaster prevention. In this paper, based on [...] Read more.
Through analyzing the triggering factors and activity characteristics of avalanches in Aerxiangou in the Western Tianshan Mountains, the formation and disaster-causing process of avalanches were studied to provide theoretical support and a scientific basis for avalanche disaster prevention. In this paper, based on remote sensing interpretation and field investigation, a spatial distribution map of avalanches was established, and the induced and triggering factors in disaster-prone environments were analyzed using the certainty factor model. The degree of influence (E) of the disaster-causing factors on avalanche triggering was quantified, and the main control conditions conducive to avalanche occurrence in different periods were obtained. The RAMMS-avalanche model was used to analyze the activity characteristics at points where multiple avalanches occurred. Research results: (1) The E values of the average temperature, average snowfall, and surface roughness in February were significantly higher than those of other hazard-causing factors, reaching 1.83 and 1.71, respectively, indicating strong control. The E values of the surface cutting degree, average temperature, and average snow depth in March were all higher than 1.8, indicating that these control factors were more prominent than the other factors. In contrast, there were four hazard-causing factors with E values higher than 1.5 in April: the mean temperature, slope, surface roughness, and mean wind speed, with clear control. (2) Under the influence of the different hazard-causing factors, the types of avalanches from February–April mainly included new full-layer avalanches, surface avalanches, and full-layer wet avalanches. (3) In the RAMMS-avalanche simulation test, considering the deposition effect, compared to the previous avalanche movement path, the secondary avalanche flow accumulation area impact range changes were slight, while the movement area within the avalanche path changes was large, as were the different categories of avalanches and their different movement characteristic values. Overall, wet snow avalanches are more hazardous, and the impact force is larger. The new snow avalanches start in a short period, the sliding rate is fast, and the avalanche sliding surface (full-snow surface and face-snow) of the difference is mainly manifested in the differences in the value of the flow height. Full article
Show Figures

Figure 1

Figure 1
<p>Topographic map of the study area.</p>
Full article ">Figure 2
<p>Technology roadmap.</p>
Full article ">Figure 3
<p>Automatic weather station in Aerxiangou.</p>
Full article ">Figure 4
<p>Monthly meteorological data for Aerxiangou for the last year (January 2022–April 2023) include (<b>a</b>) Average snow depth in Aerxiangou for the past year, month by month, (<b>b</b>) Average temperature in Aerxiangou for the last year, month by month, (<b>c</b>) Average precipitation in Aerxiangou for the last year, month by month, (<b>d</b>) Average wind speed in Aerxiangou for the last year, month by month.</p>
Full article ">Figure 4 Cont.
<p>Monthly meteorological data for Aerxiangou for the last year (January 2022–April 2023) include (<b>a</b>) Average snow depth in Aerxiangou for the past year, month by month, (<b>b</b>) Average temperature in Aerxiangou for the last year, month by month, (<b>c</b>) Average precipitation in Aerxiangou for the last year, month by month, (<b>d</b>) Average wind speed in Aerxiangou for the last year, month by month.</p>
Full article ">Figure 5
<p>Hazard-causing factor classification map: (<b>a</b>) Topographic factor grading chart; (<b>b</b>) Meteorological factor classification chart (February); (<b>c</b>) Meteorological factor classification chart (March); (<b>d</b>) Meteorological factor classification chart (April).</p>
Full article ">Figure 6
<p>Degree of influence of the hazard-causing factors on avalanche.</p>
Full article ">Figure 7
<p>Map of multiple avalanche points.</p>
Full article ">Figure 8
<p>Avalanche disaster site 4#:(<b>a</b>) Disaster situation at the site; (<b>b</b>) Study area; (<b>c</b>) Cross-sectional view of the exploratory pit.</p>
Full article ">Figure 9
<p>Illustration of the snow stress: (<b>a</b>) Full-layer avalanche; (<b>b</b>) Surface avalanche.</p>
Full article ">Figure 10
<p>Characteristic values of the movement process at avalanche point 4# under the different periods: (<b>a</b>) Maximum height; (<b>b</b>) Maximum velocity; (<b>c</b>) Maximum pressure.</p>
Full article ">Figure 10 Cont.
<p>Characteristic values of the movement process at avalanche point 4# under the different periods: (<b>a</b>) Maximum height; (<b>b</b>) Maximum velocity; (<b>c</b>) Maximum pressure.</p>
Full article ">Figure 11
<p>Variation curves of the profile I–I motion eigenvalues under the different periods: (<b>a</b>) Flow height variation curve of profile I–I under the different periods; (<b>b</b>) Flow velocity curve of profile I-I flow under the different periods; (<b>c</b>) Pressure the curve of profile I–I under the different periods.</p>
Full article ">Figure 11 Cont.
<p>Variation curves of the profile I–I motion eigenvalues under the different periods: (<b>a</b>) Flow height variation curve of profile I–I under the different periods; (<b>b</b>) Flow velocity curve of profile I-I flow under the different periods; (<b>c</b>) Pressure the curve of profile I–I under the different periods.</p>
Full article ">
25 pages, 8458 KiB  
Article
The Taconnaz Rockfall (Mont-Blanc Massif, European Alps) of November 2018: A Complex and At-Risk Rockwall-Glacier-Torrent Morphodynamic Continuum
by Ludovic Ravanel, Pierre-Allain Duvillard, Laurent Astrade, Thierry Faug, Philip Deline, Johan Berthet, Maëva Cathala, Florence Magnin, Alexandre Baratier and Xavier Bodin
Appl. Sci. 2023, 13(17), 9716; https://doi.org/10.3390/app13179716 - 28 Aug 2023
Cited by 3 | Viewed by 4221
Abstract
The glacial and torrential basin of Taconnaz (Mont-Blanc massif, France) dominates the Chamonix valley. It is one of the major paths for snow avalanches in the Alps, often triggered by serac falls from the Taconnaz glacier. On 24 November 2018, the basin’s multi-risk [...] Read more.
The glacial and torrential basin of Taconnaz (Mont-Blanc massif, France) dominates the Chamonix valley. It is one of the major paths for snow avalanches in the Alps, often triggered by serac falls from the Taconnaz glacier. On 24 November 2018, the basin’s multi-risk nature was further accentuated by a new type of hazard with a rockfall triggered at c. 2700 m a.s.l. It travelled down over a distance of 1.85 km and stopped 165 m away from the construction site of a micro-hydroelectric power station. We studied the triggering conditions at the permafrost lower limit, the effects of the supra-glacial path on the flow patterns, and the fate of the scar and the deposit on torrential activity. By comparing a pre-event Structure from Motion model with a post-event LiDAR model, we estimated the volume of the scar to be 42,900 m3 (±5%). A numerical model was employed to simulate the rapid runout. It revealed the complexity of the flow, attributed to the sequestration of a part of the deposit in crevasses, the incorporation of a significant volume of ice resulting in a transition from a dry granular flow to a mud-like flow, and the presence of numerous deposit zones. Subsequent monitoring of the area after the event allowed for the documentation of the scar’s evolution, including a landslide, as well as the progressive degradation and evacuation of the deposit by the torrent without producing debris flow. The study of the triggering factors indicated glacial retreat as the probable main cause, assisted by the melting of ice lenses left by the permafrost disappearance. Finally, we present replicable methods for managing risks at the site following the event. This event improves the understanding of cascading processes that increasingly impact Alpine areas in the context of climate change. Full article
Show Figures

Figure 1

Figure 1
<p>The glacial and torrential basin of Taconnaz (Mont-Blanc massif, France). Elevations are in m a.s.l.</p>
Full article ">Figure 2
<p>The scar of the rockfall of 24 November 2018. Left: the rock wall before the event (14 July 2018; ph. M. Pététin). In yellow: rockfall source. Right: the scar after the event (25 November 2018; ph. E. Courcier).</p>
Full article ">Figure 3
<p>Path (top left: view from the Aiguilles Rouges massif; ph. C. Taillard) of the event between the scar and the front of the deposit. The deposit left in the runout zone by the dense flow is made of boulders of different sizes and shapes, with the presence of some ice blocks (inset) embedded in a cohesive matrix of fine particles. Photographs make clear the presence of an aerosol made of dust particles in suspension in the air during the flow propagation. Main images acquired by drone. Numbers: see text.</p>
Full article ">Figure 4
<p>Three-dimensional reconstruction of the collapsed rock mass.</p>
Full article ">Figure 5
<p>2018 rockfall deposit volume estimation using the comparison of the 2008 and 2018 digital terrain models (INRAE LiDAR—photogrammetric model).</p>
Full article ">Figure 6
<p>Left panel: snapshots showing the propagation of the simulated flow at <span class="html-italic">t</span><sub>0</sub> = 6.67 s (just a few seconds after release), <span class="html-italic">t</span><sub>1</sub> = 66.67 s, <span class="html-italic">t</span><sub>i+1</sub> = <span class="html-italic">t</span><sub>i</sub> + <span class="html-italic">t</span><sub>1</sub> with I = {1, 2, …, 7} and <span class="html-italic">t</span><sub>f</sub> = 566.67 s (final deposit). Right panel: Zoom on the runout zone along the torrent bed. The colour levels correspond to the spatial distributions of flow thickness with a maximum threshold fixed at 1 m for each snapshot. The observed deposit left by the dense flow is depicted in hatched area in white colour for comparison.</p>
Full article ">Figure 7
<p>Main changes during the four months after the rockfall of 24 November 2018 (ph. E. Courcier). Yellow: landslide; ochre: rockfalls. Letters: see main text.</p>
Full article ">Figure 8
<p>Evolution of the rock/ice avalanche deposit between March and October 2019. From the beginning of the monitoring period until 17 April, the entire area was snow-covered snow. Subsequently, the snow receded on the slopes, persisting only at the valley’s base, where the deposit formed. Snow avalanches may continue to occur from the gorge or the left bank slope, as seen on 9 May, when an avalanche was triggered at higher altitudes.</p>
Full article ">Figure 9
<p>Massive ice present in the scar and thermal context of the rockfall (red dot).</p>
Full article ">Figure 10
<p>(<b>Left</b>) comparison of oblique aerial photographs (ETH-Bibliothek Zürich/Dr K. Baxter, University of Dundee). (<b>Right</b>) comparison of 3D reconstructions from IGN orthophotos.</p>
Full article ">
Back to TopTop