Antarctic Snow Failure Mechanics: Analysis, Simulations, and Applications

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.

Concentration and thickness of sea ice in the Weddell Sea from SSM/I passive microwave radiometer data

This study evaluated feasibility statistically and analyzed, during the freezing period, the relationship between brightness temperature (Tb) data of the 37V polarisation and the GR3719 (Gradient Ratio 37V and 19V) obtained by Special Sensor Microwave/Imager from F11 and F13 satellites with sea ice thickness (SIT) data obtained in the Weddell Sea through Antarctic Sea Ice Processes and Climate program. The multiple linear regression (MLR) was applied at 1,520 points, with 70% of these points being randomly separated to generate the MLR and 30% to carry out the validation. To perform the temporal mapping, the MLR was applied only to pixels with sea ice concentration (SIC) ≥ 90%, obtained through the fraction image calculated from the spectral linear mixing model (SLMM) using the Tb in the channels and polarizations 19H, 19V and 37V. The results of the SLMM validation process for estimating the SIC were σ = 10.5%, RMSE = 11.0%, and bias = -2.8%, and the SIT based on the MLR, the results were R² = 0.57, RMSE = 0.268 m, and bias = 0.103 m. In the SIT mapping, we highlight the trend of thickness reduction on the east coast of the Antarctic Peninsula during the period 1992-2009.

An Improved Algorithm for the Retrieval of the Antarctic Sea Ice Freeboard and Thickness from ICESat-2 Altimeter Data

ICESat-2 altimeter data could be used to estimate sea ice freeboard and thickness values with a higher measuring accuracy than that achievable with data provided by previous altimeter satellites. This study developed an improved algorithm considering variable lead proportions based on the lowest point method to derive the sea surface height for the retrieval of Antarctic sea ice freeboard and thickness values from ICESat-2 ATL-07 data. We first collocated ICESat-2 tracks to corresponding Sentinel-1 SAR images and calculated lead (seawater) proportions along each track to estimate the sea surface height in the Antarctic Ocean. Then, the Antarctic sea ice freeboard and thickness were estimated based on a local sea surface height reference and a static equilibrium equation. Finally, we assessed the accuracy of our improved algorithm and ICESat-2 data product in the retrieval of the Antarctic sea ice thickness by comparing the calculated values to ship-based observational sea ice thickness values acquired during the 35th Chinese Antarctic Research Expedition (CHINARE-35). The results indicate that the Antarctic sea ice freeboard estimated with the improved lowest point method was slightly larger than that estimated with the ICESat-2 data product algorithm. The root mean squared error (RMSE) of the improved lowest point method was 35 cm with the CHINARE-35 measured sea ice thickness, which was smaller than that determined with the ICESat-2 data product algorithm (65 cm). Our improved algorithm could provide more accurate data on the Antarctic sea ice freeboard and thickness, thus supporting Antarctic sea ice monitoring and the evaluation of its change under global effects.

A Model-Based Climatology of Low-Level Jets in the Weddell Sea Region of the Antarctic

Low-level jets (LLJs) are climatological features in polar regions. It is well known that katabatic winds over the slopes of the Antarctic ice sheet are associated with strong LLJs. Barrier winds occurring, e.g., along the Antarctic Peninsula may also show LLJ structures. A few observational studies show that LLJs occur over sea ice regions. We present a model-based climatology of the wind field, of low-level inversions and of LLJs in the Weddell Sea region of the Antarctic for the period 2002–2016. The sensitivity of the LLJ detection on the selection of the wind speed maximum is investigated. The common criterion of an anomaly of at least 2 m/s is extended to a relative criterion of wind speed decrease above and below the LLJ. The frequencies of LLJs are sensitive to the choice of the relative criterion, i.e., if the value for the relative decrease exceeds 15%. The LLJs are evaluated with respect to the frequency distributions of height, speed, directional shear and stability for different regions. LLJs are most frequent in the katabatic wind regime over the ice sheet and in barrier wind regions. During winter, katabatic LLJs occur with frequencies of more than 70% in many areas. Katabatic LLJs show a narrow range of heights (mostly below 200 m) and speeds (typically 10–20 m/s), while LLJs over the sea ice cover a broad range of speeds and heights. LLJs are associated with surface inversions or low-level lifted inversions. LLJs in the katabatic wind and barrier wind regions can last several days during winter. The duration of LLJs is sensitive to the LLJ definition criteria. We propose to use only the absolute criterion for model studies.

Estimation of sea-ice thickness and volume in the Sea of Okhotsk based on ICESat data

Sea-ice thickness in the Sea of Okhotsk is estimated for 2004-2008 from ICESat derived freeboard under the assumption of hydrostatic balance. Total ice thickness including snow depth (h _tot ) averaged over 2004-2008 is 95 cm. The interannual variability of h _tot is large; from 77.5 cm (2008) to 110.4 cm (2005). The mode of h _tot varies from 50-60 cm (2007 and 2008) to 70-80 cm (2005). Ice thickness derived from ICESat data is validated from a comparison with that observed by Electromagnetic Induction Instrument (EM) aboard the icebreaker Soya near Hokkaido, Japan. Annual maps of h _tot reveal that the spatial distribution of h _tot is similar every year. Ice volume of 6.3 × 10¹¹ m³ is estimated from the ICESat derived h _tot and AMSR-E derived ice concentration. A comparison with ice area demonstrates that the ice volume cannot always be represented by the area solely, despite the fact that the area has been used as a proxy of the volume in the Sea of Okhotsk. The ice volume roughly corresponds to that of annual ice production in the major coastal polynyas estimated based on heat budget calculations. This also supports the validity of the estimation of sea-ice thickness and volume using ICESat data.

Theoretical surface type classifier based on a waveform model of a satellite laser altimeter and its performance in the north of Greenland

Current land-cover classification methods using ICESat/GLAS's (Ice, Cloud, and land Elevation Satellite/Geoscience Laser Altimeter System) datasets are based on empirical thresholds or machine learning by training multiple GLAS parameters, e.g., the reflectivity and elevation of the target and width, amplitude, kurtosis, and skewness of the return waveform. A theoretical classifier is derived based on a waveform model of an actual laser altimeter illuminating the sea surface. With given system parameters and the sea surface wind corresponding to the location of a laser footprint (the wind can be calculated by using the National Centers for Environmental Prediction dataset), a precise theoretical waveform can be generated as a reference. Compared with the measured waveform, a weighted total difference, which is very sensitive to small-scale sea ice within the laser footprint, can be calculated to classify the GLAS measured data as open water. In the north of Greenland, after discarding the saturated GLAS data, the new theoretical classifier performed better [overall accuracy (OA)=95.62%, Kappa coefficient=0.8959] compared to the classical support vector machine (SVM) classifier (OA=90.44%, Kappa=0.7901), but the SVM classifier showed a better result for the user's accuracy of sea ice. Benefiting from the synergies of the theoretical and SVM classifiers, the integrated theoretical and SVM classifier achieved excellent accuracy (OA=98.21%, Kappa=0.9588). In the future, the new ICESat-2 photon counting laser altimeter will also construct a "waveform" (elevation distribution) by selecting the photon cloud, and thus, this new analytical method will be potentially useful for detecting open water in the Arctic.