Crevasses are formed by glacier movement and the stresses within glacier ice. Knowledge of the crevasses’ distribution is critical for understanding the glacier and ice shelf stability. In this study, we propose an automated crevasse extraction framework based on Sentinel-1 SAR imagery and an improved U-Net network. The spatial distribution of crevasses on Antarctic ice shelves in 2020 was mapped with a spatial resolution of ~40 m, and the characteristics of crevasses on the Nickerson Ice Shelf, Jelbart Ice Shelf, Amery Ice Shelf, Thwaites Glacier, and Shackleton Ice Shelf were analyzed. The results indicated the extraction accuracy of our method was 84.2% and the F1 score was 72.5%. Compared with previous published studies, the identification of the crevasse areas had good visual consistency. However, in some scenes, the recall rate was relatively lower due to the quality of the SAR image, terrain surrounding the crevasses, and observation geometry. The crevasses on different ice shelves had different characteristics in terms of length, density, type, and spatial pattern, implying the different stress structures of ice shelves. The Thwaites Glacier and the Nickerson Ice Shelf in the West Antarctica Ice Sheet (WAIS) had shorter ice crevasses, whereas the lengths of ice crevasses on the Jelbart Ice Shelf and the Amery Ice Shelf in the East Antarctica Ice Sheet (EAIS) were relatively long. Nevertheless, there are more closely spaced crevasses on the ice shelf in WAIS compared to that in the EAIS. For the distribution of crevasse types, the Nickerson Ice Shelf and the Shackleton Ice Shelf had various forms of crevasses. There were mainly transverse crevasses developed on the Jelbart Ice Shelf and the Amery Ice Shelf. This study provides a helpful reference and guidance for automated crevasse extraction. The method proposed by this study manifests great application potential and the efficacy of producing a time-series crevasse data set with higher spatial resolution and larger coverage. In the future, more Sentinel-1 SAR imagery will be applied and the effect of temporal and spatial variations in crevasses on the stability of ice shelves will be investigated, which will contribute to project the ice shelf stability and explore the sea level rise implications of recent and future cryosphere changes.

With global warming, supraglacial lakes play an important role in ice sheet stability and climate change. They are not only the main factors affecting mass balance and sea-level rise but also the key units of surface runoff storage and mass loss. To automatically map the spatiotemporal distribution of supraglacial lakes in Greenland, this paper proposes an attention-based U-Net model with Sentinel-1 SAR imagery. The extraction results show that compared with the traditional network, this method obtains a higher validation coefficient, with an F1 score of 0.971, and it is spatiotemporally transferable, able to realize the extraction of supraglacial lakes in complex areas without ignoring small lakes. In addition, we conducted a case study in the Jakobshavn region and found that the supraglacial lake area peaked in advance between spring and summer due to extreme melting events from 2017 to 2021. Meanwhile, the supraglacial lakes near the 79°N Glacier tended to expand inland during the melting season.

Synthetic aperture radar tomography (TomoSAR) has been proven to be a useful way to reconstruct vertical structure over forest areas with P-band images, on account of its three-dimensional imaging ability. In the case of a small number of non-uniformly distributed acquisitions, compressive sensing (CS) is generally adopted in TomoSAR. However, the performance of CS depends on the selected hyperparameter, which is closely related to the noise of a pixel. In this paper, to overcome this limitation, we propose a sparse iterative covariance-based estimation (SPICE) approach based on the wavelet and orthogonal sparse basis (W&O-SPICE) for application over forest areas. SPICE is a sparse spectral estimation method that achieves a high vertical resolution, and takes account of the noise adaptively for each resolution cell. Thus, it does not require the user to select a hyperparameter. Furthermore, the used sparse basis not only ensures the sparsity of the forest canopy scattering contribution, but it can also keep the original sparse information of the ground contribution. The proposed method was tested in simulated experiments and the results demonstrated that W&O-SPICE can successfully reconstruct the vertical structure of a forest. Moreover, three P-band fully polarimetric airborne SAR images with non-uniformly distributed baselines were applied to reconstruct the vertical structure of a tropical forest in Mabounie, Gabon. The underlying topography and forest height were estimated, and the root-mean-square errors (RMSEs) were 6.40 m and 4.50 m with respect to the LiDAR digital terrain model (DTM) and canopy height model (CHM), respectively. In addition, W&O-SPICE showed a better performance than W&O-CS, beamforming, Capon, and the iterative adaptive approach (IAA).

In polar regions, vegetation is especially sensitive to climate dynamics and thus can be used as an indicator of the global and regional environmental change. However, in Antarctica, there is very little information on vegetation distribution and growth status. To fill this gap, we evaluated the ability of both linear and nonlinear spectral mixture analysis (SMA) models, including a group of newly developed modified Nascimento’s models for Antarctic vegetated areas (MNM-AVs), in estimating the abundance of major Antarctic vegetation types, i.e., mosses and lichens. The study was conducted using WorldView-2 satellite data and field measurements over the Fildes Peninsula and its surroundings, which are representative vegetated areas in Antarctica. In MNM-AVs, we introduced secondary scattering components for vegetation and its background to account for the sparsity of vegetation cover and reassigned their coefficients. The new models achieved improved performances, among which MNM-AV3 achieved the lowest error for mosses (lichens) abundance estimation with RMSE = 0.202 (0.213). Compared with MNM-AVs, the linear model performed particularly poor for lichens (RMSE = 0.322), which is in contrast to the case of mosses (RMSE = 0.212), demonstrating that spectral signals of lichens are more prone to mix with their backgrounds. Abundance maps of mosses and lichens, as well as a map of moss health status for the entire study area, were then obtained based on MNM-AV3 with around 80% overall accuracy. Moss areas account for 0.7695 km² in Fildes and 0.3259 km² in Ardley Island; unhealthy mosses amounted to 40% (49%) of the area in the summer of 2018 (2019), indicating considerable environmental stress.

Forest height is an essential input parameter for forest biomass estimation, ecological modeling, and the carbon cycle. Tomographic synthetic aperture radar (TomoSAR), as a three-dimensional imaging technique, has already been successfully used in forest areas to retrieve the forest height. The nonparametric iterative adaptive approach (IAA) has been recently introduced in TomoSAR, achieving a good compromise between high resolution and computing efficiency. However, the performance of the IAA algorithm is significantly degraded in the case of a small tomographic aperture. To overcome this shortcoming, this paper proposes the robust IAA (RIAA) algorithm for SAR tomography. The proposed approach follows the framework of the IAA algorithm, but also considers the noise term in the covariance matrix estimation. By doing so, the condition number of the covariance matrix can be prevented from being too large, improving the robustness of the forest height estimation with the IAA algorithm. A set of simulated experiments was carried out, and the results validated the superiority of the RIAA estimator in the case of a small tomographic aperture. Moreover, a number of fully polarimetric L-band airborne tomographic SAR images acquired from the ESA BioSAR 2008 campaign over the Krycklan Catchment, Northern Sweden, were collected for test purposes. The results showed that the RIAA algorithm performed better in reconstructing the vertical structure of the forest than the IAA algorithm in areas with a small tomographic aperture. Finally, the forest height was estimated by both the RIAA and IAA TomoSAR methods, and the estimation accuracy of the RIAA algorithm was 2.01 m, which is more accurate than the IAA algorithm with 3.25 m.