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 (102)

Search Parameters:
Keywords = snow avalanche

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
22 pages, 336 KiB  
Review
What Is New in Altitude- and Cold-Related Illnesses of Travel: Appraisal and Summary of the Updated Guidelines from the Wilderness Medical Society
by Arghavan Omidi, Gregory D. Hawley, Dylan Kain, Farah Jazuli, Milca Meconnen, Mark Polemidiotis, Nam Phuong Do, Olamide Egbewumi and Andrea K. Boggild
Int. J. Environ. Res. Public Health 2025, 22(2), 284; https://doi.org/10.3390/ijerph22020284 - 14 Feb 2025
Viewed by 238
Abstract
Wilderness medicine is a rapidly evolving field and has benefited from expanded research efforts. Moreover, with an escalating occurrence of severe and cataclysmic global climatologic events, human illness arising from interaction with wilderness and recreational environments warrants increasing consideration. Within the last decade, [...] Read more.
Wilderness medicine is a rapidly evolving field and has benefited from expanded research efforts. Moreover, with an escalating occurrence of severe and cataclysmic global climatologic events, human illness arising from interaction with wilderness and recreational environments warrants increasing consideration. Within the last decade, the Wilderness Medical Society (WMS) has aggregated research findings and created guidelines on prevention measures and therapeutic options for acute altitude illness, frostbite injuries, and avalanche and non-avalanche snow burials. As new research emerges, some guidelines have been updated to reflect the most current and sound scientific conclusions. In this review, we have synthesized the evidence-based guidelines and have reviewed the quality of the guidelines according to the Appraisal of Guidelines for Research and Evaluation (AGREE) II framework. Further research efforts can expand the scope of evidence-based practice in travel medicine and ideally standardize the implementation of recommendations within both pre-travel and post-travel medical practices. Full article
(This article belongs to the Special Issue The Impact of Extreme Weather and Climate on Human Health)
19 pages, 9749 KiB  
Article
Numerical Simulation of Debris Flow Behavior over a Series of Groundsills
by Chyan-Deng Jan, Yi-Chao Zeng and Litan Dey
Water 2025, 17(3), 293; https://doi.org/10.3390/w17030293 - 21 Jan 2025
Viewed by 494
Abstract
Debris flows propagating in natural environments often encounter irregular terrain features, such as bottom roughness and man-made structures like groundsills, which significantly influence their behavior and dynamics. In practice, groundsills are commonly used as debris flow mitigation structures. This study examines the effects [...] Read more.
Debris flows propagating in natural environments often encounter irregular terrain features, such as bottom roughness and man-made structures like groundsills, which significantly influence their behavior and dynamics. In practice, groundsills are commonly used as debris flow mitigation structures. This study examines the effects of a beam-type groundsill array on the flow behavior of sediment mixtures in an inclined channel using numerical simulations. The sediment mixtures, modeled as Bingham fluids, were tested as they flowed over groundsill arrays with varying densities, characterized by the spacing-to-height ratio (d/h) ranging from 2 to 10. The findings indicate that interaction with the groundsills produces a hydraulic jump-like flow, reaching a height approximately 2.2 times the approach flow depth across different array densities. High-density arrays (d/h4) substantially hindered flow propagation, reducing front velocities but leading to sediment buildup upstream of the groundsills. Conversely, low-density arrays (d/h>4) facilitated smoother flow with higher velocities. These insights into the relationship between array density, flow behavior, and sediment trapping provide valuable guidance for optimizing groundsill array designs to effectively reduce the mobility of gravity-driven flows of non-Newtonian fluids (such as snow avalanches, debris, lava, or mudflows) and mitigate the associated risks. Full article
(This article belongs to the Section Hydraulics and Hydrodynamics)
Show Figures

Figure 1

Figure 1
<p>A series of groundsills constructed along the debris flow gully in Hualien County, Taiwan.</p>
Full article ">Figure 2
<p>Particle size distribution curve of the clay–silt materials used as sediment in the rheometer tests.</p>
Full article ">Figure 3
<p>(<b>a</b>) Schematic diagram illustrating the rheological measurement setup using a rotational rheometer [<a href="#B21-water-17-00293" class="html-bibr">21</a>]. (<b>b</b>) Rheological characteristics of the sediment mixture with <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.30</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Schematic diagrams of the experimental setup for simulating sediment mixture flow over a beam-type array of groundsills in an inclined channel: (<b>a</b>) side view and (<b>b</b>) top view at a fixed channel slope, <math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>15</mn> <mo>°</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Illustrations of beam-type groundsill arrays with (<b>a</b>) high array density (closely spaced groundsills), (<b>b</b>) low array density (widely spaced groundsills) in 2D and 3D, and (<b>c</b>) schematic diagram for the relation between <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> </mrow> </semantics></math> and channel slope, <math display="inline"><semantics> <mrow> <mi>θ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Configuration of the numerical model domain for simulating flow over an array of groundsills at a fixed slope: (<b>a</b>) boundary conditions (S: symmetry and W: wall conditions) and (<b>b</b>) computational mesh.</p>
Full article ">Figure 7
<p>Visualization of the instantaneous velocity magnitude (m/s) of sediment mixtures propagating over a low array density (specifically, <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>) of groundsills at different simulation times.</p>
Full article ">Figure 8
<p>Temporal evolution of the front position of sediment mixtures for varying groundsill array densities.</p>
Full article ">Figure 9
<p>Contour plots for the fraction of sediment mixture at simulation times of 0.7, 1.0, and 1.5 s for flow over groundsill arrays with densities of (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>, (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>, (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>8</mn> </mrow> </semantics></math>, and (<b>e</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 9 Cont.
<p>Contour plots for the fraction of sediment mixture at simulation times of 0.7, 1.0, and 1.5 s for flow over groundsill arrays with densities of (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>, (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>, (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>8</mn> </mrow> </semantics></math>, and (<b>e</b>) <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Contour plots for fraction of sediment mixture at simulation time <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>30</mn> <mo> </mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> for flow over groundsill arrays with densities of <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mn>6</mn> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mn>10</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 11
<p>Streamwise profiles of (<b>a</b>) flow depth and (<b>b</b>) depth-averaged velocity of sediment mixtures propagating over groundsill arrays at <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1.5</mn> <mo> </mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>, comparing cases with and without groundsills for different array densities.</p>
Full article ">Figure 12
<p>Schematic representation of the formation of hydraulic jump-like flows at the first groundsill for array density <math display="inline"><semantics> <mrow> <mi>d</mi> <mo>/</mo> <mi>h</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>.</p>
Full article ">
19 pages, 10750 KiB  
Article
Snow Avalanche Hazards and Avalanche-Prone Area Mapping in Tibet
by Duo Chu, Linshan Liu, Zhaofeng Wang, Yong Nie and Yili Zhang
Geosciences 2024, 14(12), 353; https://doi.org/10.3390/geosciences14120353 - 18 Dec 2024
Viewed by 636
Abstract
Snow avalanche is one of the major natural hazards in the mountain region, yet it has received less attention compared to other mountain hazards, such as landslides, floods, and droughts. After a comprehensive overview of snow avalanche hazards in Tibet area, the spatial [...] Read more.
Snow avalanche is one of the major natural hazards in the mountain region, yet it has received less attention compared to other mountain hazards, such as landslides, floods, and droughts. After a comprehensive overview of snow avalanche hazards in Tibet area, the spatial distribution and main driving factors of snow avalanche hazards in the high mountain region in Tibet were presented in the study first. Snow avalanche-prone areas in Tibet were then mapped based on the snow cover distribution and DEM data and were validated against in situ observations. Results show that there are the highest frequencies of avalanche occurrences in the southeastern Nyainqentanglha Mountains and the southern slope of the Himalayas. In the interior of plateau, avalanche development is constrained due to less precipitation and much flatter terrain. The perennially snow avalanche-prone areas in Tibet account for 1.6% of the total area of the plateau, while it reaches 2.9% and 4.9% of the total area of Tibet in winter and spring, respectively. Snow avalanche hazards and fatalities appear to be increasing trends under global climate warming due to more human activities at higher altitudes. In addition to the continuous implementation of engineering prevention and control measures in pivotal regions in southeastern Tibet, such as in the Sichuan–Tibet highway and railway sections, enhancing monitoring, early warning, and forecasting services are crucial to prevent and mitigate avalanche hazards in the Tibetan high mountain regions, which has significant implications for other global high mountain areas. Full article
(This article belongs to the Section Natural Hazards)
Show Figures

Figure 1

Figure 1
<p>Study area.</p>
Full article ">Figure 2
<p>Annual mean SCF in Tibet.</p>
Full article ">Figure 3
<p>Mean SCF in winter (<b>a</b>) and spring (<b>b</b>) in Tibet.</p>
Full article ">Figure 4
<p>Perennial snow avalanche-prone areas in Tibet.</p>
Full article ">Figure 5
<p>Winter snow avalanche-prone areas in Tibet.</p>
Full article ">Figure 6
<p>Spring snow avalanche-prone areas in Tibet.</p>
Full article ">Figure 7
<p>Field investigation on snow cover and snow avalanches in the Parlung Zangbo and Sangchu River basins. (<b>a</b>) A typical channeled snow avalanche; (<b>b</b>) snow avalanche bridge; (<b>c</b>) five channeled snow avalanches; (<b>d</b>,<b>e</b>) the destruction to forests on the mountain slope by snow avalanches.</p>
Full article ">Figure 8
<p>Snow avalanche deposits at near the Langqiu village in the Sentinel-2 image (left). (<b>a</b>) Avalanche deposit in the location 1 in the Sentinel-2 image; (<b>b</b>) avalanche deposit in the location 2 in the Sentinel-2 image.</p>
Full article ">Figure 9
<p>Snow avalanche deposits at Galongla section from Zhamo to Metok highway in the Sentinel-2 image (left). (<b>a</b>) Avalanche deposit in the location 1 in the Sentinel-2 image; (<b>b</b>) avalanche deposit in the location 2 in the Sentinel-2 image.</p>
Full article ">Figure 10
<p>Perennial snow avalanche-prone area in the Parlung Zangbo and Sangchu River basins in southeastern Tibet.</p>
Full article ">Figure 11
<p>Snow avalanche-prone areas in winter in the Parlung Zangbo and Sangchu River basins in southeastern Tibet.</p>
Full article ">Figure 12
<p>Snow avalanche-prone areas in spring in the Parlung Zangbo and Sangchu River basins in southeastern Tibet.</p>
Full article ">
20 pages, 98934 KiB  
Article
Automated Snow Avalanche Monitoring and Alert System Using Distributed Acoustic Sensing in Norway
by Antoine Turquet, Andreas Wuestefeld, Guro K. Svendsen, Finn Kåre Nyhammer, Espen Lauvlund Nilsen, Andreas Per-Ola Persson and Vetle Refsum
GeoHazards 2024, 5(4), 1326-1345; https://doi.org/10.3390/geohazards5040063 - 17 Dec 2024
Viewed by 1011
Abstract
Avalanches present substantial hazard risk in mountainous regions, particularly when avalanches obstruct roads, either hitting vehicles directly or leaving traffic exposed to subsequent avalanches during cycles. Traditional detection methods often are designed to cover only a limited section of a road stretch, hampering [...] Read more.
Avalanches present substantial hazard risk in mountainous regions, particularly when avalanches obstruct roads, either hitting vehicles directly or leaving traffic exposed to subsequent avalanches during cycles. Traditional detection methods often are designed to cover only a limited section of a road stretch, hampering effective risk management. This research introduces a novel approach using Distributed Acoustic Sensing (DAS) for avalanche detection. The monitoring site in Northern Norway is known to be frequently impacted by avalanches. Between 2022–2024, we continuously monitored the road for avalanches blocking the traffic. The automated alert system identifies avalanches affecting the road and estimates accumulated snow. The system provides continuous, real-time monitoring with competitive sensitivity and accuracy over large areas (up to 170 km) and for multiple sites on parallel. DAS powered alert system can work unaffected by visual barriers or adverse weather conditions. The system successfully identified 10 road-impacting avalanches (100% detection rate). Our results via DAS align with the previous works and indicate that low frequency part of the signal (<20 Hz) is crucial for detection and size estimation of avalanche events. Alternative fiber installation methods are evaluated for optimal sensitivity to avalanches. Consequently, this study demonstrates its durability and lower maintenance requirements, especially when compared to the high setup costs and coverage limitations of radar systems, or the weather and lighting vulnerabilities of cameras. Furthermore the system can detect vehicles on the road as important supplemental information for search and rescue operations, and thus the authorities can be alerted, thereby playing a vital role in urgent rescue efforts. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Map showing Existing Cable (blue) and New Cable extension (orange). Photos taken during the installation are attached from (i) the cabinet at the northern end of the monitoring system, (ii) the vehicle warning system end of the north section, (iii) the vehicle warning system beginning of the south section (iv) an example photo of microtrenching. (<b>b</b>) Map of Norway and the region surrounding the avalanche monitoring zone. Important places are marked, including Holmbuktura, the location of the installation. (<b>c</b>) A cross section sketch showing the details of microtrench cable installation (iv). Direct buried new cable is installed at 15 cm depth and plastic tube covered installation is done at 20 cm depth from the surface.</p>
Full article ">Figure 2
<p>Aerial overview of the Holmbuktura region detailing characteristic avalanche paths and the avalanche monitoring setup. The image on the left (<b>a</b>) shows a comprehensive view of the valley with shaded areas for avalanche zones in north and south. The paths (1–5) along the slope show 5 characteristic avalanche paths, delineating the primary areas of avalanche activity. The cyan line represents the trajectory of the sensor cable installation, placed to capture both the dynamics of avalanches and the road traffic activity. The plot on the right (<b>b</b>) shows the altitude evolution along 5 selected paths, giving the impression of the topography of the region. Image © 2024 Google Earth, Image Landsat/Copernicus, Image © 2024 Maxar Technologies, Image © 2024 CNES/Airbus.</p>
Full article ">Figure 3
<p>Simplified flowchart of the automated avalanche detection and monitoring system. Data is continuously collected and processed through edge computing in two separate modules: (1) vehicle detection and (2) avalanche detection, which operate independently to avoid interference. Detected avalanches and vehicles are then transferred to a central repository and messaging module. This module evaluates risk levels, checks for stranded or at-risk vehicles, and prepares necessary visualizations and alerts. If the risk level exceeds a predefined threshold, the system sends alerts, including plots and messages, via SCADA message system and email using 4G communication.</p>
Full article ">Figure 4
<p>Examples of signals recorded during monitoring with the DAS system in Holmbuktura are shown. The strain rate waterfall plot (Z) highlights features of different events: (<b>a</b>) avalanche activity in the north, (<b>b</b>) avalanche activity in the south, (<b>c</b>) a passenger car, and (<b>d</b>) a snowplow.</p>
Full article ">Figure 5
<p>The power spectral density (PSD) was computed from signals recorded during monitoring with the DAS system in Holmbuktura. The signals represent distinct events, specifically: (<b>a</b>) avalanche activity in the north, (<b>b</b>) avalanche activity in the south, (<b>c</b>) a passenger car, and (<b>d</b>) a snowplow.</p>
Full article ">Figure 6
<p>Most energetic traces from all avalanches are presented as raw signals. In (<b>a</b>) avalanche signals are presented and marked with “Zone N” and “Zone S” showing where the avalanches happened. Event 0 is an avalanche which stopped right before the road it is presented for comparison. Corresponding mean frequency of the 200 s trace is computed and marked on the end of trace. In (<b>b</b>) we present the spectrogram of all avalanches. (<b>c</b>,<b>d</b>) we compare the power spectra of north avalanches and south avalanches respectively. The associated log-averaged power spectra are also plotted.</p>
Full article ">Figure 7
<p>Co-located direct buried “D” (red) and piped loopback cable “P” (blue) traces from avalanches only hitting the southern section are presented (<b>a</b>). Corresponding mean frequency of the entire trace is computed and marked on the trace as well. On right we compare the power spectra from direct buried cable (<b>b</b>) and piped buried (<b>c</b>). The associated log-average power spectra are also plotted.</p>
Full article ">Figure 8
<p>Detailed analysis of the most energetic trace from Event 5. The avalanche signal is analyzed using 20 s sliding time window to investigate avalanche dynamics. In (<b>a</b>), the normalized signal is shown in the time domain; (<b>b</b>) presents the mean frequency of the 20 s time window sliding every 1 s; and (<b>c</b>) displays the power spectra of selected time windows. The colored boxes in (<b>a</b>) indicate time windows, which are highlighted with markers in the mean frequency plot (<b>b</b>) and in the power spectra plot (<b>c</b>) in corresponding colors.</p>
Full article ">Figure A1
<p>Comparison of avalanche dates with historical data of environmental variables. Temperature, snow depth, rain and wind speed from the region covering October 2022 to May 2024 is obtained from OpenMeteo [<a href="#B69-geohazards-05-00063" class="html-bibr">69</a>] is presented. We have plotted the 200 h moving averaged data to visualize long term trends.</p>
Full article ">
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 927
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 937
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
Cited by 2 | Viewed by 1531
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 816
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
Cited by 1 | Viewed by 883
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
Cited by 1 | Viewed by 1661
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 1127
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 1599
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
Cited by 1 | Viewed by 1263
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 ">
44 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 3244
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
Cited by 2 | Viewed by 1493
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 ">
Back to TopTop