Search

Greenland is the world's largest island and over three-quarter of the entire land area is covered by the Greenland Ice Sheet (GrIS), which is the only permanent ice sheet besides Antarctica. Finally, all satellite images were mosaicked to finish the full image mosaic of Greenland. [Extracted from the article]

On 26 September 2019, a massive iceberg broke off the west side of the Amery Ice Shelf (AIS) in East Antarctica. Since 1973, the AIS calving front has steadily advanced at a rate of 1.0 km yr−1. However, the advancement rate of the central portion of the AIS increased dramatically during 2012–2015, which indicates a velocity increase prior to the calving event. Eight calving front locations from 1973 to 2018 were mapped to investigate the advancement rate of AIS over the entire observational period. Additionally, the propagation of rift A was observed unstable from 2012 to 2015. The westward propagation rate of rift A1 increased to 3.7 km yr−1 from 2015 to 2017, which was considerably faster than the other rifts near the AIS calving front. The increased advancement rate and the increasing propagation magnitude of at least one active rift appear to be precursors of this large calving event. [ABSTRACT FROM AUTHOR]

Ice Pathfinder (Code: BNU-1), launched on 12 September 2019, is the first Chinese polar observation microsatellite. Its main payload is a wide-view camera with a ground resolution of 74 m at the subsatellite point and a scanning width of 744 km. BNU-1 takes into account the balance between spatial resolution and revisit frequency, providing observations with finer spatial resolution than Terra/Aqua MODIS data and more frequent revisits than Landsat-8 OLI and Sentinel-2 MSI. It is a valuable supplement for polar observations. Geolocation is an essential step in satellite image processing. This study aims to geolocate BNU-1 images; this includes two steps. For the first step, a geometric calibration model is applied to transform the image coordinates to geographic coordinates. The images calibrated by the geometric model are the Level1A (L1A) product. Due to the inaccuracy of satellite attitude and orbit parameters, the geometric calibration model also exhibits errors, resulting in geolocation errors in the BNU-1 L1A product. Then, a geometric correction method is applied as the second step to find the control points (CPs) extracted from the BNU-1 L1A product and the corresponding MODIS images. These CPs are used to estimate and correct geolocation errors. The BNU-1 L1A product corrected by the geometric correction method is processed to the Level1B (L1B) product. Although the geometric correction method based on CPs has been widely used to correct the geolocation errors of visible remote sensing images, it is difficult to extract enough CPs from polar images due to the high reflectance of snow and ice. In this study, the geometric correction employs an image division and an image enhancement method to extract more CPs from the BNU-1 L1A products. The results indicate that the number of CPs extracted by the division and image enhancements increases by about 30% to 182%. Twenty-eight images of Antarctica and fifteen images of Arctic regions were evaluated to assess the performance of the geometric correction. The average geolocation error was reduced from 10 km to ~300 m. In general, this study presents the geolocation method, which could serve as a reference for the geolocation of other visible remote sensing images for polar observations. [ABSTRACT FROM AUTHOR]

Digital Elevation Models (DEMs) of Greenland provide the basic data for studying the Greenland ice sheet (GrIS), but little research quantitatively evaluates and compares the accuracy of various Greenland DEMs. This study uses IceBridge elevation data to evaluate the accuracies of the the Greenland Ice Map Project (GIMP)1 DEM, GIMP2 DEM, TanDEM-X, and ArcticDEM in their corresponding time ranges. This study also analyzes the impact of DEM accuracy and resolution on the accuracy of river network extraction. The results show that (1) within the time range covered by each DEM, TanDEM-X with an RMSE of 5.60 m has higher accuracy than the other DEMs in terms of absolute height accuracy, while GIMP1 has the lowest accuracy among the four Greenland DEMs, with an RMSE of 14.34 m. (2) Greenland DEMs are affected by regional errors and interannual changes. The accuracy in areas with elevations above 2000 m is higher than that in areas with elevations below 2000 m, and better accuracy is observed in the north than in the south. The stability of the ArcticDEM product is higher than those of the other three DEM products, and its RMSE standard deviation over multiple years is only 0.14 m. Therefore, the errors caused by the applications of DEMs with longer time spans are smaller. GIMP1 performs in an opposite manner, with a standard deviation of 2.39 m. (3) The river network extracted from TanDEM-X is close to the real river network digitized from remote sensing images, with an accuracy of 50.78%. The river network extracted from GIMP1 exhibits the largest errors, with an accuracy of only 8.83%. This study calculates and compares the accuracy of four Greenland DEMs and indicates that TanDEM-X has the highest accuracy, adding quantitative studies on the accuracy evaluation of various Greenland DEMs. This study also compares the results of different DEM river network extractions, verifies the impact of DEM accuracy on the accuracy of the river network extraction results, and provides an explorable direction for the hydrological analysis of Greenland as a whole. [ABSTRACT FROM AUTHOR]

Passive microwave (PMW) sensors are popularly applied to Earth observations. However, the satellite PMW radiometer data sometimes have non-negligible errors in geolocation. Coastline inflection methods (CIMs) are widely used to improve geolocation errors of PMW images. However, they commonly require accuracy satellite flight parameters, which are difficult to obtain by users. In this study, a simplified coastline inflection method (SCIM) is proposed to correct the geolocation errors without demanding for the satellite flight parameters. SCIM is applied to improve geolocation errors of FengYun-3D (FY-3D) Microwave Radiation Imager (MWRI) brightness temperature images from 2018 and 2019. It reduces the geolocation errors of MWRI images to 0.15 pixels in the along-track and cross-track direction. This means reductions of 75% and 86% of the geolocation errors, respectively. The mean brightness temperature differences between the ascending and descending MWRI images are reduced by 34%, demonstrating the improved geolocation accuracy of SCIM. The corrected images are also used to estimate Arctic sea ice concentration (SIC). By comparing with SICs retrieved from the un-corrected images, the root mean square error (RMSE) and mean absolute error (MAE) of the SICs from the corrected images are reduced from 13.7% to 10.2% and 8.9% to 6.9%, respectively. The mean correlation coefficient (R) increases from 0.91 to 0.95. All these results indicate that SCIM can reduce geolocation errors of satellite-based PMW images significantly. As SCIM is very simple and easy to be applied, it could be a useful method for users of PMW images. [ABSTRACT FROM AUTHOR]

The mountainous vegetation is important to regional sustainable development. However, the topographic effect is the main obstacle to the monitoring of mountainous vegetation using remote sensing. Aiming to retrieve the reflectance of frequently-used red–green–blue and near-infrared (NIR) wavebands of rugged mountains for vegetation mapping, we developed a new integrated topographic correction (ITC) using the SCS + C correction and the shadow-eliminated vegetation index. The ITC procedure consists of image processing, data training, and shadow correction and uses a random forest machine learning algorithm. Our study using the Landsat 8 Operational Land Imager (OLI) multi-spectral images in Fujian province, China, showed that the ITC achieved high performance in topographic correction of regional mountains and in transferability from the sunny area of a scene to the shadow area of three scenes. The ITC-corrected multi-spectral image with an NIR–red–green composite exhibited flat features with impressions of relief and topographic shadow removed. The linear regression of corrected waveband reflectance vs. the cosine of the solar incidence angle showed an inclination that nearly reached the horizontal, and the coefficient of determination decreased to 0.00~0.01. The absolute relative errors of the cast shadow and the self-shadow all dramatically decreased to the range of 0.30~6.37%. In addition, the achieved detection rate of regional vegetation coverage for the three cities of Fuzhou, Putian, and Xiamen using the ITC-corrected images was 0.92~6.14% higher than that using the surface reflectance images and showed a positive relationship with the regional topographic factors, e.g., the elevation and slope. The ITC-corrected multi-spectral images are beneficial for monitoring regional mountainous vegetation. Future improvements can focus on the use of the ITC in higher-resolution imaging. [ABSTRACT FROM AUTHOR]

Surface meltwater runoff is believed to be the main cause of the alarming mass loss in the Greenland Ice Sheet (GrIS); however, recent research has shown that a large amount of meltwater is not directly drained or refrozen but stored in the form of firn aquifers (FAs) in the interior of the GrIS. Monitoring the changes in FAs over the GrIS is of great importance to evaluate the stability and mass balance of the ice sheet. This is challenging because FAs are not visible on the surface and the direct measurements are lacking. A new method is proposed to map FAs during the 2010–2020 period by using the C-band Advanced Scatterometer (ASCAT) data based on the Random Forests classification algorithm with the aid of measurements from the NASA Operation IceBridge (OIB) program. Melt days (MD), melt intensity (MI), and winter mean backscatter (WM) parameters derived from the ASCAT data are used as the input vectors for the Random Forests classification algorithm. The accuracy of the classification model is assessed by ten-fold cross-validation, and the overall accuracy and Kappa coefficient are 97.49% and 0.72 respectively. The results show that FAs reached the maximum in 2015, and the accumulative area of FAs from 2010 to 2020 is 56,477 km2, which is 3.3% of the GrIS area. This study provides a way to investigate the long-term dynamics in FAs which have great significance for understanding the state of subsurface firn and subglacial hydrological systems. [ABSTRACT FROM AUTHOR]

In recent decades, the melting of the Greenland Ice Sheet (GrIS) has become one of the major causes of global sea-level rise. Supraglacial lakes (SGLs) are typical hydrological features produced on the surface of the GrIS during the melt seasons. The existence and evolution of SGLs play an important role in the melting process of the ice sheet surface. To understand the distribution and recent changes of SGLs in Greenland, this study developed a random forest (RF) algorithm incorporating the texture and morphological features to automatically identify SGLs based on the Google Earth Engine (GEE) platform. Sentinel-2 imagery was used to map the SGLs inventory in Greenland during the 2016–2018 melt seasons and to explore the spatial and temporal variability characteristics of SGLs. Our results show changes in SGLs from 2016 to 2018, with the total area decreasing by ~1152.22 km2 and the number increasing by 1134; SGLs are mainly distributed in western Greenland (SW, CW, NW) and northeastern Greenland (NE), where the NE region has the largest number of observed SGLs and the largest SGL was with the surface area of 16.60 km2 (2016). SGLs were found to be most active in the area with the elevation of 800–1600 m and the slope of 0–5°, and showed a phenomenon of retreating to lower elevation areas and developing to steeper slope areas. Our work provided a method for rapid inventory of SGLs. This study will help monitor the mass balance of the GrIS and predict future rapid ice loss from Greenland. [ABSTRACT FROM AUTHOR]

Iceberg D28, a giant tabular iceberg that calved from Amery Ice Shelf in September 2019, grounded off Kemp Coast, East Antarctica, from August to September of 2020. The motion of the iceberg is characterized herein by time-series images captured by synthetic aperture radar (SAR) on Sentinel-1 and the moderate resolution imaging spectroradiometer (MODIS) boarded on Terra from 6 August to 15 September 2020. The thickness of iceberg D28 was estimated by utilizing data from altimeters on Cryosat-2, Sentinel-3, and ICESat-2. By using the iceberg draft and grounding point locations inferred from its motion, the maximum water depths at grounding points were determined, varying from 221.72 ± 21.77 m to 269.42 ± 25.66 m. The largest disagreements in seabed elevation inferred from the grounded iceberg and terrain models from the Bedmap2 and BedMachine datasets were over 570 m and 350 m, respectively, indicating a more complicated submarine topography in the study area than that presented by the existing seabed terrain models. Wind and sea water velocities from reanalysis products imply that the driving force from sea water is a more dominant factor than the wind in propelling iceberg D28 during its grounding, which is consistent with previous findings on iceberg dynamics. [ABSTRACT FROM AUTHOR]

Recent studies show that many debris-covered glaciers exhibit high-magnitude differential thinning despite the presence of supraglacial debris. Existing studies have not adequately explored the variability in surface ablation, morphology, and related feedback mechanisms that incorporate topographically controlled surface irradiance and debris transport. In this study, we address these issues using a radiation-driven surface ablation model that more fully characterizes ablation dynamics by accounting for temporally-linked radiative forcing, surface geomorphological evolution and gravitational debris flux. Simulation results based on Baltoro Glacier in the central Karakoram indicate the following: 1) A debris-covered glacier can exhibit high spatial variability in surface ablation due to heterogeneous debris thickness and debris transport. A bare-ice glacier given similar initial conditions, has higher overall ablation but exhibits much less variability. 2) The topographic influence on surface ablation is non-negligible because glacier-surface topography controls irradiance and gravitational debris flux. The overall ablation on a debris-covered glacier tends to increase in response to high-frequency topographic variations due to the larger area with thin debris cover. In contrast, a bare-ice glacier exhibits decreased overall ablation in response to high-frequency topographic variations due to its high sensitivity to topographic shading. 3) Gravity-driven surface debris flux plays an important role in local debris thickness redistribution, which regulates ablation rates over the ablation season. 4) Surface ablation dynamics on a debris-covered glacier is regulated by active system couplings and feedbacks between surface morphology, melt and debris transport. Consequently, certain locations on a debris-covered glacier may be more sensitive to radiative forcing than previously thought. Simulation results suggest that nonlinear feedback responses may permit debris-covered glacier subsystems to be very sensitive to radiative forcing, thereby causing significant morphological changes to the glacier surface. • Feedbacks exist between glacier surface melt, topography and debris transport. • Numerical models are developed to study the feedbacks on glacier surface. • Gravitational process on an evolving topography redistributes debris thickness. • Ice loss on a bare-ice glacier is sensitive to topographic shading. • Complex surface topography may enhance ice loss on a debris-covered glacier. [ABSTRACT FROM AUTHOR]

Landfast ice is an important component of the Antarctic sea ice regime. It affects the Antarctic climate and ecological system. In this study, the first high-resolution, long time series of the landfast ice edge from 2006 to 2011 and 2016 to 2017 is presented. The dataset was produced based on the improved net gradient difference algorithm using 2470 SAR scenes from ENVISAT and Sentinel-1A/B as well as manual analysis of MODIS imagery to fill in SAR data gaps. The study results show that the landfast ice area in November for all studied years was approximately 49.49 ± 3.25 × 104 km2, accounting for about 3%–4% of the total Antarctic sea ice area. The maximum area was 55.70 × 104 km2 in November 2007, compared to the minimum area 44.01 × 104 km2 in 2011. The area in West Antarctica was about 40% of that in East Antarctica. The distribution of landfast ice in Antarctica has significant regional differences. The extent in the Indian Ocean sector is the maximum with a mean value of 16.49 ± 1.1 × 104 km2; however, the ratio of the landfast ice area to the sea ice area in the Pacific Ocean sector is the highest. Twenty-four landfast ice zones with groups of small, grounded icebergs were identified, most of which were located in East Antarctica, particularly along the Wilkes Land and Oates Land. Two cases are presented to illustrate how giant, grounded icebergs affected landfast ice. Results from this study are well suited to underpin the Antarctic climate or ecological system studies. • The first long time series and high-resolution circumpolar mapping of the landfast ice in Antarctica • The spatio-temporal patterns of the landfast ice in Antarctica have been revealed. • The interactions between landfast ice and giant icebergs are analyzed in two specific cases. [ABSTRACT FROM AUTHOR]

Research involving anisotropic-reflectance correction (ARC) of multispectral imagery to account for topographic effects has been ongoing for approximately 40 years. A large body of research has focused on evaluating empirical ARC methods, resulting in inconsistent results. Consequently, our research objective was to evaluate commonly used ARC methods using first-order radiation-transfer modeling to simulate ASTER multispectral imagery over Nanga Parbat, Himalaya. Specifically, we accounted for orbital dynamics, atmospheric absorption and scattering, direct- and diffuse-skylight irradiance, land cover structure, and surface biophysical variations to evaluate their effectiveness in reducing multi-scale topographic effects. Our results clearly reveal that the empirical methods we evaluated could not reasonably account for multi-scale topographic effects at Nanga Parbat. The magnitude of reflectance and the correlation structure of biophysical properties were not preserved in the topographically-corrected multispectral imagery. The CCOR and SCS+C methods were able to remove topographic effects, given the Lambertian assumption, although atmospheric correction was required, and we did not account for other primary and secondary topographic effects that are thought to significantly influence spectral variation in imagery acquired over mountains. Evaluation of structural-similarity index images revealed spatially variable results that are wavelength dependent. Collectively, our simulation and evaluation procedures strongly suggest that empirical ARC methods have significant limitations for addressing anisotropic reflectance caused by multi-scale topographic effects. Results indicate that atmospheric correction is essential, and most methods failed to adequately produce the appropriate magnitude and spatial variation of surface reflectance in corrected imagery. Results were also wavelength dependent, as topographic effects influence radiation-transfer components differently in different regions of the electromagnetic spectrum. Our results explain inconsistencies described in the literature, and indicate that numerical modeling efforts are required to better account for multi-scale topographic effects in various radiation-transfer components. [ABSTRACT FROM AUTHOR]

The temporal and spatial variability of the Antarctic coastline is a clear indicator of change in extent and mass balance of ice sheets and shelves. In this study, the Canny edge detector was utilized to automatically extract high-resolution information of the Antarctic coastline for 2005, 2010, and 2017, based on optical and microwave satellite data. In order to improve the accuracy of the extracted coastlines, we developed the Canny algorithm by automatically calculating the local low and high thresholds via the intensity histogram of each image to derive thresholds to distinguish ice sheet from water. A visual comparison between extracted coastlines and mosaics from remote sensing images shows good agreement. In addition, comparing manually extracted coastline, based on prior knowledge, the accuracy of planimetric position of automated extraction is better than two pixels of Landsat images (30 m resolution). Our study shows that the percentage of deviation (<100 m) between automatically and manually extracted coastlines in nine areas around the Antarctica is 92.32%, and the mean deviation is 38.15 m. Our results reveal that the length of coastline around Antarctica increased from 35,114 km in 2005 to 35,281 km in 2010, and again to 35,672 km in 2017. Meanwhile, the total area of the Antarctica varied slightly from 1.3618 × 107 km2 (2005) to 1.3537 × 107 km2 (2010) and 1.3657 × 107 km2 (2017). We have found that the decline of the Antarctic area between 2005 and 2010 is related to the breakup of some individual ice shelves, mainly in the Antarctic Peninsula and off East Antarctica. We present a detailed analysis of the temporal and spatial change of coastline and area change for the six ice shelves that exhibited the largest change in the last decade. The largest area change (a loss of 4836 km2) occurred at the Wilkins Ice Shelf between 2005 and 2010. [ABSTRACT FROM AUTHOR]