Search

Studying glacier changes is essential for understanding global glacier mass balance, patterns in climate change, and the potential consequences they may have. The glaciers of Geladandong Mountain serve as a crucial water source for the Yangtze River and its downstream regions. Any changes in the glacier at the river source will directly impact the ecological environment and water recycling process downstream. This study presents an approach for monitoring Geladandong Mountain glacier elevation changes using ICESat-2 data combined with a Binned Multilevel Regression Fitting Method (BMRFM). We analysed multi-year ICESat-2 data, considering various topographic features, to understand glaciers' annual and seasonal elevation changes. Our research reveals significant spatial heterogeneity in glacier elevation changes, influenced by altitude, slope, and aspect factors. It demonstrates the effectiveness of ICESat-2 in glacier monitoring, offering insights into the impacts of climate change on glacier dynamics. Moreover, the results indicate a predominant trend of ongoing melting at lower altitudes, with a small amount of accretion at higher elevations. Over the last two decades, the overall elevation change rate is −0.23 ± 0.12 metres per year. Seasonal analysis shows substantial glacier thickening from April to August, primarily due to increased precipitation offsetting elevated temperatures. This study offers a robust method for monitoring glaciers using satellite data, thereby advancing the precision and reliability of glacier research. [ABSTRACT FROM AUTHOR]

Satellite-based estimation of forest variables including forest biomass relies on model-based approaches since forest biomass cannot be directly measured from space. Such models require ground reference data to adapt to the local forest structure and acquired satellite data. For wide-area mapping, such reference data are too sparse to train the biomass retrieval model and approaches for calibrating that are independent from training data are sought. In this study, we compare the performance of one such calibration approach with the traditional regression modelling using reference measurements. The performance was evaluated at four sites representative of the major forest biomes in Europe focusing on growing stock volume (GSV) prediction from time series of C-band Sentinel-1 and Advanced Land Observing Satellite Phased Array L-band Synthetic Aperture Radar (ALOS-2 PALSAR-2) backscatter measurements. The retrieval model was based on a Water Cloud Model (WCM) and integrated two forest structural functions. The WCM trained with plot inventory GSV values or calibrated with the aid of auxiliary data products correctly reproduced the trend between SAR backscatter and GSV measurements across all sites. The WCM-predicted backscatter was within the range of measurements for a given GSV level with average model residuals being smaller than the range of the observations. The accuracy of the GSV estimated with the calibrated WCM was close to the accuracy obtained with the trained WCM. The difference in terms of root mean square error (RMSE) was less than 5% units. This study demonstrates that it is possible to predict biomass without providing reference measurements for model training provided that the modelling scheme is physically based and the calibration is well set and understood. [ABSTRACT FROM AUTHOR]

The ICESat-2 and SWOT satellite earth observation missions have provided highly accurate water surface slope (WSS) observations in global rivers for the first time. While water surface slope is expected to remain constant in time for approximately uniform flow conditions, we observe time varying water surface slope in many river reaches around the globe in the ICESat-2 record. Here, we investigate the causes of time variability of WSSs using simplified river hydraulic models based on the theory of steady, gradually varied flow. We identify bed slope or cross section shape changes, river confluences, flood waves, and backwater effects from lakes, reservoirs, or the ocean as the main non-uniform hydraulic situations in natural rivers that cause time changes of WSSs. We illustrate these phenomena at selected river sites around the world, using ICESat-2 data and river discharge estimates. The analysis shows that WSS observations from space can provide new insights into river hydraulics and can enable the estimation of river discharge from combined observations of water surface elevation and WSSs at sites with complex hydraulic characteristics. [ABSTRACT FROM AUTHOR]

The sea ice extent in the Weddell Sea exhibited a positive trend from the start of satellite observations in 1978 until 2016 but has shown a decreasing trend since then. This study analyzes seasonal and interannual variations in sea ice thickness using ICESat-2 laser altimetry data over the Weddell Sea from 2019 to 2022. Sea ice thickness was calculated from ICESat-2's ATL10 freeboard product using the Improved Buoyancy Equation. Seasonal variability in ice thickness, characterized by an increase from February to September, is more pronounced in the eastern Weddell sector, while interannual variability is more evident in the western Weddell sector. The results were compared with field data obtained between 2019 and 2022, showing a general agreement in ice thickness distributions around predominantly level ice. A decreasing trend in sea ice thickness was observed when compared to measurements from 2003 to 2017. Notably, the spring of 2021 and summer of 2022 saw significant decreases in Sea Ice Extent (SIE). Although the overall mean sea ice thickness remained unchanged, the northwestern Weddell region experienced a noticeable decrease in ice thickness. [ABSTRACT FROM AUTHOR]

Mapping forest canopy height is critical for climate modeling and forest management, and tropical forests present unique challenges for remote sensing due to their dense vegetation and complex structure. The advent of ICESat-2 and GEDI, two advanced lidar datasets, offers new opportunities for improving canopy height estimation. In this study, we used footprint-level canopy height products from ICESat-2 and GEDI, combined with features extracted from Landsat-8, PALSAR-2, and FABDEM products. The AutoGluon stacking ensemble learning algorithm was employed to construct inversion models, generating 30 m resolution continuous canopy height maps for the tropical forests of Puerto Rico. Accuracy validation was performed using the high-resolution G-LiHT airborne lidar products. Results show that tropical forest canopy height inversion remains challenging, with all models yielding relative root mean square errors (rRMSE) exceeding 0.30. The stacking ensemble model outperformed all base learners, and the GEDI-based map had slightly higher accuracy than the ICESat-2-based map, with RMSE values of 4.81 and 4.99 m, respectively. Both models showed systematic biases, but the GEDI-based model exhibited less underestimation for taller canopies, making it more suitable for biomass estimation. The proposed approach can be applied to other forest ecosystems, enabling fine-resolution canopy height mapping and enhancing forest conservation efforts. [ABSTRACT FROM AUTHOR]

The empirical formula (EF) method, which do not rely on topographic data, stands as the prevailing technique for estimating lake water storage variation (LWSV). However, for smaller lakes, the sporadic monitoring frequency of satellite altimetry fails to adequately support this method, presenting a challenge in accurately gauging LWSV. Using Lake Chahannur, a lake in China with an area smaller than 50 km2, as a case study, seven schemes based on the EF method and the Area-Volume-Height (A-V-H) curve method were designed to estimate the LWSV of this undersized lake. The efficacy and precision of each scheme were evaluated against field-measured elevations. Findings reveal that due to the limited satellite altimetry monitoring, both the EF method and the H-driven A-V-H curve schemes struggle to provide consistent and comprehensive estimations. In the A-driven A-V-H curve schemes, terrain data from SRTM DEM suffers from mask processing and substantial errors, with the former posing challenges for shrinking lakes and the latter significantly compromising estimation accuracy. While field-measured elevations boast high precision, the interpolation process leads to terrain maps lacking in detail, with site density becoming a crucial factor influencing the accuracy of LWSV estimation. The combination of terrain reconstruction and A-driven pattern emerges as the most promising, boasting high accuracy, rich detail, and significantly reduced reliance on satellite altimetry monitoring, making it particularly suitable for small lakes. Chahannur's bottom elevation ranges between 1271.71 and 1273.44 m, and the lake shows a downward trend in water volume from 1991 to 2020, with fluctuations totaling approximately 35 million m3. This study serves as a vital addition to the field of LWSV estimation, potentially broadening the scope of estimation from large-scale lakes to a wider array of global surface water bodies. [ABSTRACT FROM AUTHOR]

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) has great potential for development due to its advantages of the use of multiple beams, low energy consumption, high repetition frequency, and high measurement sensitivity. However, the weak photon signal emitted by the photon counting lidar is susceptible to the background noise caused by the sun and the atmosphere, which can seriously affect the processing and application of laser data. This paper proposes an improved DBSCAN clustering algorithm for denoising single photon laser point clouds in mountainous areas. Firstly, a grouping method based on elevation and distance statistics is proposed to reduce the influence of terrain undulations on denoising accuracy. Finally, an automatic radius search method is put forward to determine clustering radius of each group, automatically find the optimal radius, and improve the existing DBSCAN clustering method. The method proposed in this paper is compared with the classical DBSCAN algorithm. The results show that the proposed algorithm significantly improves denoising accuracy in mountainous areas and effectively filters out most background noise. [ABSTRACT FROM AUTHOR]

In September 2018, NASA launched the ICESat-2 satellite into a 500 km high Earth orbit. It carried a truly unique lidar system, i.e., the Advanced Topographic Laser Altimeter System or ATLAS. The ATLAS lidar is capable of detecting single photons reflected from a wide variety of terrain (land, ice, tree leaves, and underlying terrain) and even performing bathymetric measurements due to its green wavelength. The system uses a single 5-watt, Q-switched laser producing a 10 kHz train of sub-nanosecond pulses, each containing 500 microjoules of energy. The beam is then split into three "strong" and three "weak" beamlets, with the "strong" beamlets containing four times the power of the "weak" beamlets in order to satisfy a wide range of Earth science goals. Thus, ATLAS is capable of making up to 60,000 surface measurements per second compared to the 40 measurements per second made by its predecessor multiphoton instrument, the Geoscience Laser Altimeter System (GLAS) on ICESat-1, which was terminated after several years of operation in 2009. Low deadtime timing electronics are combined with highly effective noise filtering algorithms to extract the spatially correlated surface photons from the solar and/or electronic background noise. The present paper describes how the ATLAS system evolved from a series of unique and seemingly unconnected personal experiences of the author in the fields of satellite laser ranging, optical antennas and space communications, Q-switched laser theory, and airborne single photon lidars. [ABSTRACT FROM AUTHOR]

The intense changes in glaciers in the southeastern Tibetan Plateau (SETP) have essential impacts on regional water resource management. In order to study the seasonal fluctuations of glaciers in this region and their relationship with climate change, we focus on the Yigong Zangbo River Basin in the SETP, extract the annual and seasonal variations of glaciers in the basin during 2018–2023, and analyze their spatio-temporal characteristics through the seasonal-trend decomposition using the LOESS (STL) method. Finally, combining the Extreme Gradient Boosting (XGBoost) model and the Shapley additive explanations (SHAP) model, we assess the comprehensive impact of meteorological factors such as temperature and snowfall on glacier changes. The results indicate that glaciers in the Yigong Zangbo River Basin experienced remarkable mass loss during 2018–2023, with an average annual melting rate of −0.83 ± 0.12 m w.e.∙yr−1. The glacier mass exhibits marked seasonal fluctuations, with increases in January–March (JFM) and April–June (AMJ) and noticeable melting in July–September (JAS) and October–December (OND). The changes over these four periods are 2.12 ± 0.04 m w.e., 0.93 ± 0.15 m w.e., −1.58 ± 0.19 m w.e., and −1.32 ± 0.17 m w.e., respectively. Temperature has been identified as the primary meteorological driver of glacier changes in the study area, surpassing the impact of snowfall. This study uses advanced altimetry data and meteorological data to monitor and analyze glacier changes, which provides valuable data for cryosphere research and also validates a set of replicable research methods, which provides support for future research in related fields. [ABSTRACT FROM AUTHOR]

The photon-counting laser altimeter on the Ice, Cloud and Land Elevation 2 (ICESat-2) satellite offers a solution for tracking the dynamic water level variations in medium and large inland lakes. We utilize the monthly ATL13 global inland water data of ICESat-2 satellite from 2018 to 2020 to estimate and analyze the water level change in Poyang Lake. The measured data from Hukou, Xingzi and Kangshan hydrological stations are used for verification and error correction, and the water level and rainfall data of each station are combined to analyze the dynamic variation of Poyang Lake water level and reveal the underlying drivers. The results show that, the annual water level of Poyang Lake varied sharply with obvious seasonal variations and an overall upward trend;the high water level period was from June to October, which peaked from July to September. The linear correlation coefficient of water levels between ICESat-2 and measured data is above 0. 846, rising to 0. 974 after error correction. The Root Mean Square Error (RMSE) is 1. 660 m, 1. 073 m, and 0. 836 m for Hukou, Xingzi, and Kangshan stations, respectively;error correction and recalculation can decrease the RMSE to 0. 663 m, 0. 659 m, and 0. 440 m for Hukou, Xingzi, and Kangshan stations, respectively, enhancing the measurement accuracy by nearly one meter. The variation of water level in Poyang Lake is highly correlated with the change of rainfall, the reduced precipitation during periods from January to February and October to December corresponds to the declining water level in dry season, while the increased rainfall from March to October corresponds to the water level rise in wet season, and the precipitation concentration period from July to September aligns with the peak of water level in Poyang Lake. [ABSTRACT FROM AUTHOR]

Objectives: To achieve the accurate monitoring of sea ice conditions in the key area of Arctic passage, high resolution sea ice thickness is essential. However, sea ice thickness with kilometer-scale resolution cannot meet the requirement at present. Methods: Considering the relationship between sea ice thickness and image albedo, this paper tried to establish a regression model between the ICESat-2 ATL10 sea ice freeboard product and the high resolution Sentinel-2 optical image albedo by taking advantage of the intensive altimetry data of ICESat-2 along the orbit to obtain dense sea ice freeboard. Based on high resolution sea ice thickness calculated from dense freeboard with the hydrostatic equilibrium model and combined with high resolution sea ice concentration derived from Sentinel-2 image, we classified ice areas into different navigational categories according to the ice-break capability of ship and compared optimal route design at multiple spatial scales. Results: Regression models established from the ATL10 sea ice freeboard and the Sentinel-2 image albedo have good fitting accuracy. The R² of the regression model is higher than 0.5, the average deviation is less than 0.05 m, and the root mean square error of the accuracy validation is less than 0.2 m. High resolution sea ice parameters derived by the proposed method can describe the distribution of fine leads between floating ice which are difficult to be captured by low resolution sea ice parameters. Conclusions: High resolution sea ice parameters can describe the details of sea ice conditions more accurately. Therefore, the proposed high resolution inversion method of sea ice parameters can improve the capability of navigation planning for ships in Arctic passage, further improving the navigation safety. [ABSTRACT FROM AUTHOR]

The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) is an active spaceborne remote sensing system that utilizes photon-counting LiDAR to capture highly detailed information about under-vegetation terrain and forest structure over vast spatial regions. It facilitates the accurate retrieval of terrain elevation and canopy height information, critical for assessing the global carbon budget and understanding the role of forests in climate change mitigation. However, challenges arise from the characteristics of the ICESat-2 photon-counting LiDAR data, such as their linear distribution, extensive spatial coverage, and substantial residual noise. These challenges hinder the performances of the state-of-the-art methods when applied on ICESat-2 data for extracting ground or top of canopy, while they perform well on airborne LiDAR that is featured with planar distribution, small coverage, and high signal-to-noise ratio. Consequently, this study proposes a novel algorithm termed Adaptive Linear Cloth Simulation Filtering (ALCSF), for the automated extraction of ground and top-of-canopy photons from ICESat-2 signal photons. The ALCSF algorithm innovatively introduces a cloth strip model as a reference to accommodate the distribution characteristics of ICESat-2 photons. Additionally, it employs a terrain-adaptive strategy to adjust the rigidity of the cloth strip by utilizing terrain slope information, thus making ALCSF applicable to large-scale areas with significant topographical changes. Furthermore, the proposed ALCSF addresses noise interference by simultaneously considering the movability of particles of the cloth strip model and the photon distribution during iterative adjustments of the cloth strip. The performance of the ALCSF is evaluated by comparing it with the ICESat-2 Land–Vegetation Along-Track Products (ATL08) across twelve datasets that encompass various times of day and scenes. In the results, the ALCSF exhibits notable improvements over ATL08 products, effectively reducing the root mean square error (RMSE) of ground elevation by 21.8% and canopy height by 25.8%, with superior performance in preserving terrain details. This highlights the significance of ALCSF as a valuable tool for enhancing the accuracy of ICESat-2 land and vegetation products, ultimately contributing to the estimation of the global carbon budget in future studies. [ABSTRACT FROM AUTHOR]

Antarctica has an active subglacial hydrological system, with interconnected subglacial lakes fed by subglacial meltwater. Subglacial hydrology can influence basal sliding, inject freshwater into the sub‐ice‐shelf cavity, and impact sediment transport and deposition which can affect the stability of grounding lines (GLs). We used satellite altimetry data from the ICESat, ICESat‐2, and CryoSat‐2 missions to document the second recorded drainage of Engelhardt Subglacial Lake (SLE), which began in July 2021 and discharged more than 2.3 km3 of subglacial water into the Ross Ice Shelf cavity. We used differential synthetic aperture radar interferometry from RADARSAT‐2 and TerraSAR‐X alongside ICESat‐2 repeat‐track laser altimetry (RTLA) and REMA digital elevation model strips to detect 2–13 km of GL retreat since the previous drainage event in 2003–06. Combining these satellite observations, we evaluated the mechanism triggering SLE drainage, the cause of the observed GL retreat, and the interplay between subglacial hydrology and GL dynamics. We find that: (a) SLE drainage was initiated by influx from a newly identified upstream lake; (b) the observed GL retreat is mainly driven by the continued retreat of Engelhardt Ice Ridge and long‐term dynamic thinning that caused a grounded ice plain to reach flotation; and (c) SLE drainage and GL retreat were largely independent. We also discuss the possible origins and influence of a 27 km grounded promontory found to protrude seaward from the GL. Our observations demonstrate the importance of high‐resolution satellite data for improving the process‐based understanding of dynamic and complex regions around the Antarctic Ice Sheet margins. Plain Language Summary: Large volumes of water flow beneath the Antarctic Ice Sheet through an interconnected network of rivers and lakes. This water system impacts slipperiness at the base of the ice, affecting how fast it moves. It also delivers freshwater into the ocean, directly contributing to sea‐level rise and increasing melt beneath the floating ice shelves. In this study, we use satellite data to track the 2021–24 drainage of Engelhardt Subglacial Lake in West Antarctica. This lake is located close to the Ross Ice Shelf grounding line, the point at the edge of the ice sheet where the ice first lifts off the bedrock and starts to float on the ocean. This region of the grounding line has retreated by up to 13 km since the last lake drainage in 2003–06. Here, we investigate what caused the drainage, the reasons for the grounding line retreat, and whether the two processes are connected. We also report the growth of a grounded promontory extending 27 km out to sea, which may be evidence of a former ice stream moraine or melt channel. These findings help to improve our limited understanding of the relationship between subglacial hydrology and grounding line dynamics in Antarctica. Key Points: Satellite altimetry detects second recorded drainage of Engelhardt Subglacial Lake in July 2021, triggered by an influx from an upstream lakeThe grounding line retreated by 2–13 km since the last drainage in 2003, linked to the retreat of Engelhardt Ice Ridge and ice plain ungroundingSatellite observations suggest that the 2021–24 lake drainage and grounding line retreat were largely independent processes [ABSTRACT FROM AUTHOR]

Abstract The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) has great potential for development due to its advantages of the use of multiple beams, low energy consumption, high repetition frequency, and high measurement sensitivity. However, the weak photon signal emitted by the photon counting lidar is susceptible to the background noise caused by the sun and the atmosphere, which can seriously affect the processing and application of laser data. This paper proposes an improved DBSCAN clustering algorithm for denoising single photon laser point clouds in mountainous areas. Firstly, a grouping method based on elevation and distance statistics is proposed to reduce the influence of terrain undulations on denoising accuracy. Finally, an automatic radius search method is put forward to determine clustering radius of each group, automatically find the optimal radius, and improve the existing DBSCAN clustering method. The method proposed in this paper is compared with the classical DBSCAN algorithm. The results show that the proposed algorithm significantly improves denoising accuracy in mountainous areas and effectively filters out most background noise. Graphical Abstract

ABSTRACTSpaceborne photon-counting light detection and ranging (LiDAR) have been extensively applied in shallow-water bathymetry. The density of underwater topographic photons (UTP) varies and is discontinuous due to sunlight noise, beam intensity, and seabed reflectivity, which differ from the land photon distribution due to the attenuation of water. Therefore, a general method for extracting and grading UTP is still lacking. We propose an active contour method combined with a variable convolution kernel method to calculate the photon range by considering the energy contributions of adjacent photons. Adaptive parameters under different observation conditions were determined to obtain the optimal convolution kernel using a kernel ridge regression model. This implies that the number of photons contained in the buffer zone was largest after the extracted UTP was fitted to a curve. Quantitative and qualitative verifications proved that the method performed well under different conditions. The photons obtained by the energy functional and the curve obtained by the fitting method were then used to grade the photons. Finally, an online developed UTP dataset and extraction framework were proposed to provide an applicable method for current and subsequent spaceborne photon-counting LiDAR.

Abstract Peatlands in the Canadian boreal forest are being negatively impacted by anthropogenic climate change, the effects of which are expected to worsen. Peatland types and sub‐classes vary in their ecohydrological characteristics and are expected to have different responses to climate change. Large‐scale modelling frameworks such as the Canadian Model for Peatlands, the Canadian Fire Behaviour Prediction System and the Canadian Land Data Assimilation System require peatland maps including information on sub‐types and vegetation as critical inputs. Additionally, peatland class and vegetation height are critical variables for wildlife habitat management and are related to the carbon cycle and wildfire fuel loading. This research aimed to create a map of peatland sub‐classes (bog, poor fen, rich fen permafrost peat complex) for the Canadian boreal forest and create an inventory of peatland vegetation height characteristics using ICESat‐2. A three‐stage hierarchical classification framework was developed to map peatland sub‐classes within the Canadian boreal forest circa 2020. Training and validation data consisted of peatland locations derived from various sources (field data, aerial photo interpretation, measurements documented in literature). A combination of multispectral data, L‐band SAR backscatter and C‐Band interferometric SAR coherence, forest structure and ancillary variables was used as model predictors. Ancillary data were used to mask agricultural areas and urban regions and account for regions that may exhibit permafrost. In the first stage of the classification, wetlands, uplands and water were classified with 86.5% accuracy. In the second stage, within the wetland areas only, peatland and mineral wetlands were differentiated with 93.3% accuracy. In the third stage, constrained to only the peatland areas, bogs, rich fens, poor fens and permafrost peat complexes were classified with 71.5% accuracy. Then, ICESat‐2 ATL08 spaceborne lidar data were used to describe regional variations in peatland vegetation height characteristics and regional and class‐wise variations based on a boreal forest wide sample. This research introduced a comprehensive large‐scale peatland sub‐class mapping framework for the Canadian boreal forest, presenting the first moderate resolution map of its kind.

The leaf area index (LAI) is a critical variable for forest ecosystem processes. Passive optical and active LiDAR remote sensing have been used to retrieve LAI. LiDAR data have good penetration to provide vertical structure distribution and deliver the ability to estimate forest LAI, such as the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2). Segment size and beam type are important for ICESat-2 LAI estimation, as they affect the amount of signal photons returned. However, the current ICESat-2 LAI estimation only covered a limited number of sites, and the performance of LAI estimation with different segment sizes has not been clearly compared. Moreover, ICESat-2 LAIs derived from strong and weak beams lack a comparative analysis. This study derived and evaluated LAI from ICESat-2 data over the National Ecological Observatory Network (NEON) sites in North America. The LAI estimated from ICESat-2 for different segment sizes (20, 100, and 200 m) and beam types (strong beam and weak beam) were compared with those from the airborne laser scanning (ALS) and the Copernicus Global Land Service (CGLS). The results show that the LAI derived from strong beams performs better than that of weak beams because more photon signals are received. The LAI estimated from the strong beam at the 200 m segment size shows the highest consistency with those from the ALS data (R = 0.67). Weak beams also present the potential to estimate LAI and have moderate agreement with ALS (R = 0.52). The ICESat-2 LAI shows moderate consistency with ALS for most forest types, except for the evergreen forest. The ICESat-2 LAI shows satisfactory agreement with the CGLS 300 m LAI product (R = 0.67, RMSE = 1.94) and presents a higher upper boundary. Overall, the ICESat-2 can characterize canopy structural parameters and provides the ability to estimate LAI, which may promote the LAI product generated from the photon-counting LiDAR. [ABSTRACT FROM AUTHOR]

There now exists two decades of basin-wide coverage of Arctic sea ice from three dedicated polar-orbiting altimetry missions (ICESat, CryoSat-2, and ICESat-2) launched by NASA and ESA. Here, we review our retrieval approaches and discuss the composite record of Arctic ice thickness (2003–2023) after appending two more years (2022–2023) to our earlier records. The present availability of five years of snow depth estimates—from differencing lidar (ICESat-2) and radar (CryoSat-2) freeboards—have benefited from the concurrent operation of two altimetry missions. Broadly, the dramatic volume loss (5500 km3) and Arctic-wide thinning (0.6 m) captured by ICESat (2003–2009), primarily due to the decline in old ice coverage between 2003 and 2007, has slowed. In the central Arctic, away from the coasts, the CryoSat-2 and shorter ICESat-2 records show near-negligible thickness trends since 2007, where the winter and fall ice thicknesses now hover around 2 m and 1.3 m, from a peak of 3.6 m and 2.7 m in 1980. Ice volume production has doubled between the fall and winter with the faster-growing seasonal ice cover occupying more than half of the Arctic Ocean at the end of summer. Seasonal ice behavior dominates the Arctic Sea ice's interannual thickness and volume signatures. [ABSTRACT FROM AUTHOR]

Forest Canopy Height (FCH) is a crucial parameter that offers valuable insights into forest structure. Spaceborne LiDAR missions provide accurate FCH measurements, but a significant challenge is their point-based measurements lacking spatial continuity. This study integrated ICESat-2's ATL08-derived FCH values with multi-temporal and multi-source remote sensing (RS) datasets to generate continuous FCH maps for northern forests in Iran. Sentinel-1/2, ALOS-2 PALSAR-2, and FABDEM datasets were prepared in Google Earth Engine (GEE) for FCH mapping, each possessing unique spatial and geometrical characteristics that differ from those of the ATL08 product. Given the importance of accurately representing the geometrical characteristics of the ATL08 segments in modeling FCH, a novel Weighted Kernel (WK) approach was proposed in this paper. The WK approach could better represent the RS datasets within the ATL08 ground segments compared to other commonly used resampling approaches. The correlation between all RS data features improved by approximately 6% compared to previously employed approaches, indicating that the RS data features derived after convolving the WK approach are more predictive of FCH values. Furthermore, the WK approach demonstrated superior performance among machine learning models, with random forests outperforming other models, achieving a coefficient of determination (R2) of 0.71, root mean square error (RMSE) of 4.92 m, and mean absolute percentage error (MAPE) of 29.95%. Furthermore, in contrast to previous studies using only summer datasets, this study included spring and autumn data from Sentinel-1/2, resulting in a 6% increase in R2 and a 0.5-m decrease in RMSE. The proposed methodology filled the research gaps and improved the accuracy of FCH estimations. [ABSTRACT FROM AUTHOR]

The huge volume of data from the Ice, Cloud and land Elevation Satellite-2 (ICESat-2), designed for mapping polar ice, sea ice, and continental vegetation, requires a highly automated data analysis and reliable terrain classification. In particular, we have developed a method to identify 4 distinct terrain categories in observed terrain, namely ocean, land, sea ice, and ice sheets. This study performed the following efforts: first, the spatial distribution characteristics for each of the 4 categories within individual ICESat-2 "major frames" along the orbit were extracted; second, these features were fed into Classification and Regression Tree (CART) and Random Forest (RF) for training; and lastly, post-processing enhancement was used to improve the classification results. Based on the 76,891 major frame samples (10,764,740 m along track) acquired via various ICESat-2 datasets, the accuracy of the two model were calculated using ten-fold cross-validation. The results indicate that the RF algorithm obtained higher classification accuracy (average accuracy [AA] = 0.9353, overall accuracy [OA] = 0.9342, and Cohen's Kappa coefficient [kappa] = 0.9122) when compared with the CART algorithm (AA = 0.9066, OA = 0.9057, and kappa = 0.8743). Overall, our approach can effectively reduce the workload of human field investigation or visual inspection of altimetry data, improve the accuracy for Earth surface classification, and add to the variety of ways to obtain global surface information from ICESat-2 data. [ABSTRACT FROM AUTHOR]

Building height is a crucial indicator when studying urban environments and human activities, necessitating accurate, large-scale, and fine-resolution calculations. However, mainstream machine learning-based methods for inferring building heights face numerous challenges, including limited sample data and slow update frequencies. Alternatively, satellite laser altimetry technology offers a reliable means of calculating building heights with high precision. Here, we initially calculated building heights along satellite orbits based on building-rooftop contour vector datasets and ICESat-2 ATL03 photon data from 2019 to 2022. By integrating multi-source passive remote sensing observation data, we used the inferred building height results as reference data to train a random forest model, regressing building heights at a 10 m scale. Compared with ground-measured heights, building height samples constructed from ICESat-2 photon data outperformed methods that indirectly infer building heights using total building floor number. Moreover, the simulated building heights strongly correlated with actual observations at a single-city scale. Finally, using several years of inferred results, we analyzed building height changes in Tianjin from 2019 to 2022. Combined with the random forest model, the proposed model enables large-scale, high-precision inference of building heights with frequent updates, which has significant implications for global dynamic observation of urban three-dimensional features. [ABSTRACT FROM AUTHOR]

Peatlands in the Canadian boreal forest are being negatively impacted by anthropogenic climate change, the effects of which are expected to worsen. Peatland types and sub‐classes vary in their ecohydrological characteristics and are expected to have different responses to climate change. Large‐scale modelling frameworks such as the Canadian Model for Peatlands, the Canadian Fire Behaviour Prediction System and the Canadian Land Data Assimilation System require peatland maps including information on sub‐types and vegetation as critical inputs. Additionally, peatland class and vegetation height are critical variables for wildlife habitat management and are related to the carbon cycle and wildfire fuel loading. This research aimed to create a map of peatland sub‐classes (bog, poor fen, rich fen permafrost peat complex) for the Canadian boreal forest and create an inventory of peatland vegetation height characteristics using ICESat‐2. A three‐stage hierarchical classification framework was developed to map peatland sub‐classes within the Canadian boreal forest circa 2020. Training and validation data consisted of peatland locations derived from various sources (field data, aerial photo interpretation, measurements documented in literature). A combination of multispectral data, L‐band SAR backscatter and C‐Band interferometric SAR coherence, forest structure and ancillary variables was used as model predictors. Ancillary data were used to mask agricultural areas and urban regions and account for regions that may exhibit permafrost. In the first stage of the classification, wetlands, uplands and water were classified with 86.5% accuracy. In the second stage, within the wetland areas only, peatland and mineral wetlands were differentiated with 93.3% accuracy. In the third stage, constrained to only the peatland areas, bogs, rich fens, poor fens and permafrost peat complexes were classified with 71.5% accuracy. Then, ICESat‐2 ATL08 spaceborne lidar data were used to describe regional variations in peatland vegetation height characteristics and regional and class‐wise variations based on a boreal forest wide sample. This research introduced a comprehensive large‐scale peatland sub‐class mapping framework for the Canadian boreal forest, presenting the first moderate resolution map of its kind. [ABSTRACT FROM AUTHOR]

Most work to date on satellite-derived bathymetry (SDB) depth change estimates water depth at individual times t1 and t2 using two separate models and then differences the model estimates. An alternative approach is explored in this study: a multi-temporal Sentinel-2 image is created by "stacking" the bands of the times t1 and t2 images, geographically coincident reference data for times t1 and t2 allow for "true" depth change to be calculated for the pixels of the multi-temporal image, and this information is used to fit a single model that estimates depth change directly rather than indirectly as in the model-differencing approach. The multi-temporal image approach reduced the depth change RMSE by about 30%. The machine learning modelling method (categorical boosting) outperformed linear regression. Overfitting of models was limited even for the CatBoost models having the maximum number of variables examined. The visible Sentinel-2 spectral bands contributed most to the model predictions. Though the multi-temporal stacked image approach produced clearly superior depth change estimates compared to the conventional approach, it is limited only to those areas for which geographically coincident multi-temporal reference/"true" depth data exist. [ABSTRACT FROM AUTHOR]

Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) has excellent potential for obtaining water depth information around islands and reefs. Combining the density-based spatial clustering of applications with noise algorithm (DBSCAN) and multiscale quadtree analysis, we propose a new photon-counting lidar denoising method to discard the large amount of noise in ICESat-2 data. First, the kernel density estimation (KDE) is used to preprocess the point cloud data, and a threshold is set to remove the noise photons on the sea surface. Next, the DBSCAN algorithm is used to preliminarily remove underwater noise photons. Then, the quadtree segmentation and Otsu algorithm are used for fine denoising to extract accurate bottom signal photons. Based on ICESat-2 pho-ton-counting data from six typical islands and reefs worldwide, the proposed method outperforms other algorithms in terms of denoising effect. Compared to in situ data, the determination coefficient (R2) reaches 94.59%, and the root mean square error (RMSE) is 1.01 m. The proposed method can extract accurate underwater terrain information, laying a foundation for offshore bathymetry. [ABSTRACT FROM AUTHOR]

Satellite-derived bathymetry (SDB) is an effective means of obtaining global shallow water depths. However, the effect of inherent optical properties (IOPs) on the accuracy of SDB under different water quality conditions has not been clearly clarified. To enhance the accuracy of machine learning SDB models, this study aims to assess the performance improvement of integrating the quasi-analytical algorithm (QAA)-derived IOPs using the Sentinel-2 and ICESat-2 datasets. In different water quality experiments, the results indicate that four SDB models (the Gaussian process regression, neural networks, random forests, and support vector regression) incorporating QAA-IOP parameters equal to or outperform those solely based on the remote sensing reflectance (Rrs) datasets, especially in turbid waters. By analyzing information gains in SDB, the most effective inputs are identified and prioritized under different water qualities. The SDB method incorporating QAA-IOP can achieve an accuracy of 0.85 m, 0.48 m, and 0.74 m in three areas (Wenchang, Laizhou Bay, and the Qilian Islands) with different water quality. Also, we find that incorporating an excessive number of redundant bands into machine learning models not only increases the demand of computing resources but also leads to worse accuracy in SDB. In conclusion, the integration of QAA-IOPs offers promising improvements in obtaining bathymetry and the optimal feature selection should be carefully considered in diverse aquatic environments. [ABSTRACT FROM AUTHOR]

NASA's Ice, Cloud, and land Elevation Satellite (ICESat-2), designed for surface altimetry, plays a pivotal role in providing precise ice sheet elevation measurements. While its primary focus is altimetry, ICESat-2 also offers valuable atmospheric data. Current conventional processing methods for producing atmospheric data products encounter challenges, particularly in conditions with low signal or high background noise. The thresholding technique traditionally used for atmospheric feature detection in lidar data uses a threshold value to accept signals while rejecting noise, which may result in signal loss or false detection in the presence of excessive noise. Traditional approaches for improving feature detection, such as averaging, lead to a trade-off between detection resolution and accuracy. In addition, the discrimination of cloud from aerosol in the identified features is difficult given ICESat-2's single wavelength and lack of depolarization measurement capability. To address these challenges, we demonstrate atmospheric feature detection and cloud–aerosol discrimination using deep learning-based semantic segmentation by a convolutional neural network (CNN). The key findings from our research are the effectiveness of a deep learning model for feature detection and cloud–aerosol classification in ICESat-2 atmospheric data and the model's surprising capability to detect complex atmospheric features at a finer resolution than is currently possible with traditional processing techniques. We identify several examples where the traditional feature detection and cloud–aerosol discrimination algorithms struggle, like in scenarios with several layers of vertically stacked clouds, or in the presence of clouds embedded within aerosol, and demonstrate the ability of the CNN model to detect such features, resolving the boundaries between adjacent layers and detecting clouds hidden within aerosol layers at a fine resolution. [ABSTRACT FROM AUTHOR]

The ICESat‐2 (Ice, Cloud and Land Elevation Satellite‐2) photon‐counting laser altimeter technology required the design and development of very sophisticated onboard algorithms to collect, store and downlink the observations. These algorithms utilize both software and hardware solutions for meeting data volume requirements and optimizing the science achievable via ICESat‐2 measurements. Careful planning and dedicated development were accomplished during the pre‐launch phase of the mission in preparation for the 2018 launch. Once on‐orbit all of the systems and subsystems were evaluated for performance, including the receiver algorithms, to ensure compliance with mission standards and satisfy the mission science objectives. As the mission has progressed and the instrument performance and data volumes were better understood, there have been several opportunities to enhance ICESat‐2's contributions to Earth observation science initiated by NASA and the ICESat‐2 science community. We highlight multiple updates to the flight receiver algorithms, the onboard software for signal processing, that have extended ICESat‐2's data capabilities and allowed for advanced science applications beyond the original mission objectives. Plain Language Summary: NASA launched its second Earth observing laser altimeter in 2018 with mission objectives of collecting observations in support of Earth science as a window into climate change impacts on our planet. Pre‐launch studies focused on specific instrument settings and on‐board data processing to support the mission objectives without violating data volume constraints. Once the instrument was on‐orbit and operational, evaluation of the algorithms for success in signal detection, signal finding, and signal telemetry was undertaken. In response to the evaluation, updates have been made to optimize the data provided by the mission. Key Points: Since the 2018 ICESat‐2 (Ice, Cloud and Land Elevation Satellite‐2) launch multiple updates have been made to the onboard flight receiver algorithm parameters to improve access and utility of the data for a multi‐disciplinary science communityThe adjustments have been primarily in the in the vertical telemetry band settings and have facilitated enhanced observations of blowing snow and increased detection opportunities of bathymetry in nearshore environmentsThe parameter changes have been made to mitigate data losses in certain situations and to advance science applications outside of the primary science objectives of the ICESat‐2 mission [ABSTRACT FROM AUTHOR]

Accurate forest structural parameters (such as forest height and canopy cover) support forest carbon monitoring, sustainable forest management, and the implementation of silvicultural practices. The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2), which is a spaceborne Light Detection and Ranging (LiDAR) satellite, offers significant potential for acquiring precise and extensive information on forest structural parameters. However, the ICESat-2 ATL08 product is significantly influenced by the geographical environment and forest characteristics, maintaining considerable potential for enhancing the accuracy of forest height estimation. Meanwhile, it does not focus on providing canopy cover data. To acquire accurate forest structural parameters, the Terrain Signal Neural Network (TSNN) framework was proposed, integrating Computer Vision (CV), Ordering Points to Identify the Clustering Structure (OPTICS), and deep learning. It encompassed an advanced approach for detecting terrain vegetation signals and constructing deep learning models for estimating forest structural parameters using ICESat-2 ATL03 raw data. First, the ATL03 footprints were visualized as Profile Raster Images of Footprints (PRIF), implementing image binarization through adaptive thresholding and median filtering denoising to detect the terrain. Second, the rough denoising buffers were created based on the terrain, combining with the OPTICS clustering and Gaussian denoising algorithms to recognize the terrain vegetation signal footprints. Finally, deep learning models (convolutional neural network (CNN), ResNet50, and EfficientNetB3) were constructed, training standardized PRIF to estimate forest structural parameters (including forest height and canopy cover). The results indicated that the TSNN achieved high accuracy in terrain detection (coefficient of determination (R2) = 0.97) and terrain vegetation signal recognition (F-score = 0.72). The EfficientNetB3 model achieved the highest accuracy in forest height estimation (R2 = 0.88, relative Root Mean Squared Error (rRMSE) = 13.5%), while the CNN model achieved the highest accuracy in canopy cover estimation (R2 = 0.80, rRMSE = 18.5%). Our results have significantly enhanced the accuracy of acquiring ICESat-2 forest structural parameters, while also proposing an original approach combining CV and deep learning for utilizing spaceborne LiDAR data. [ABSTRACT FROM AUTHOR]

The Ice, Cloud, and Land Elevation Satellite‐2 (ICESat‐2) mission has collected surface elevation measurements for over 5 years. ICESat‐2 carries an instrument that emits laser light at 532 nm, and ice and snow absorb weakly at this wavelength. Previous modeling studies found that melting snow could induce significant bias to altimetry signals, but there is no formal assessment on ICESat‐2 acquisitions during the melting season. We performed two case studies over the Greenland Ice Sheet to quantify bias in ICESat‐2 signals over snow: one to validate Airborne Topographic Mapper (ATM) data against Next Generation Airborne Visible/Infrared Imaging Spectrometer (AVIRIS‐NG) grain sizes, and a second to estimate ICESat‐2 bias relative to ATM. We used snow optical grain sizes derived from ATM and AVIRIS‐NG to attribute altimetry bias to snowpack properties. For the first case study, the mean and standard deviation of optical grain sizes were 340 ± 65 µm (AVIRIS‐NG) and 670 ± 420 µm (ATM). A mean altimetry bias of 4.81 ± 1.76 cm was found for ATM, with larger biases linked to increases in grain size. In the second case study, we found a mean grain size of 910 ± 381 µm and biases of 6.42 ± 1.77 cm (ICESat‐2) and 9.82 ± 0.97 cm (ATM). The grain sizes and densities needed to recreate biases with a model are uncommon in nature, so we propose that additional surface attributes must be considered to characterize ICESat‐2 bias over snow. The altimetry biases are within the accuracy requirements of the ICESat‐2 mission, but we cannot rule out more significant errors over coarse‐grained snow. Plain Language Summary: The Ice, Cloud, and Land Elevation Satellite‐2 (ICESat‐2) mission has been used to measure changes in land ice, vegetation cover, and sea ice, and there is growing interest to use ICESat‐2 for snow science. ICESat‐2 uses a green laser to estimate the elevation of the Earth's surface, and recent work has shown that snow can introduce errors in green laser measurements such as ICESat‐2. In this study, we used ICESat‐2 and airborne data to (a) identify errors in the ICESat‐2 data and (b) link the errors to snow properties over the Greenland Ice Sheet. We found that ICESat‐2 errors are linked to two snow properties: the size of individual snow grains and the density of the snow. Two retrieval methods were used to estimate snow grain size, though the methods show significant differences when large snow grains are observed. Snow grain size and density are insufficient to fully explain ICESat‐2 biases, so we propose that other factors also play a role. The ICESat‐2 errors are within the accuracy requirements of ICESat‐2, but more significant errors may be possible during the Northern Hemisphere melting season. Key Points: Very large snow grains and densities are needed to explain centimeter‐level differences between ICESat‐2 and airborne lidarBiases from snow grain size and other factors may be influenced by different height measurement methods between ICESat‐2 and airborne lidarSnow grain size retrievals strongly differ between airborne lidar and AVIRIS‐NG over regions with coarse‐grained snow [ABSTRACT FROM AUTHOR]

Forecasting flood and drought events requires accurate modeling tools. Hydraulic river models are based on estimates of riverbed geometry which are traditionally collected in situ. The novel Ice, Cloud and Land Elevation Satellite 2 [ICESat‐2] lidar altimetry mission with 6 simultaneous high‐resolution laser beams provides the opportunity to define river cross‐section geometries as well as observe water surface elevation [WSE] and water surface slope spatially resolved along the river chainage. This paper describes a method to utilize terrain altimetry and water surface slope estimates to define complete river geometries from ICESat‐2 data products, using the diffusive wave approximation to calculate depth in the submerged section not penetrated by the lidar. Exemplifying the method, cross‐sections are defined for a stretch of the Mekong River. Hydrodynamic model results of the stretch are compared with ICESat‐2 WSE estimates and in situ gauging station time series. Insights in river characteristics from satellite imagery and the ICESat‐2 slope estimates allow for fine‐tuning of the cross‐sections using spatially varying Manning numbers. The final model achieves a root mean square error against the ICESat‐2 WSE of 0.676 m and average Kling‐Gupta Efficiency against gauging station time series of 0.880. The method is limited by the diffusive wave approximation resulting in inaccurate cross‐section estimates in sections with supercritical flow or significant acceleration. Errors can be identified from ICESat‐2 WSE estimates and reduced with additional cross‐sections. Combined with hydrological models, the method will allow for cross‐section definition without in situ data. Plain Language Summary: The depth and width of the river channel are important factors when seeking to predict floods. To predict water level in a river, computer models for flood forecasting must be informed with river channel geometry at cross‐sections along the river. Such cross‐sections are traditionally measured in the field which is expensive and time consuming. This paper describes a method to estimate river cross‐sectional shape from land and water height measurements from the satellite ICESat‐2. When the satellite path crosses the river, the precise laser instrument onboard outlines the surface to indicate the river channel shape, but it does not penetrate the water. We found that it was possible to estimate the shape of the submerged part of the cross‐section using the water surface slope obtained from the satellite's unique instrument. Evaluating the method on a stretch of the Mekong River, we were able to model water level along the river with accuracy sufficient for flood forecasting. In some sections, where the flow speed changes quickly, a hydraulic model with the defined cross‐sections does not reproduce observed water levels. Using satellite crossings from days with steadier flow reduces the error. The method allows for building hydraulic river models without cross‐section surveys. Key Points: Land elevation and water surface slope estimates from the ICESat‐2 lidar altimetry mission are combined to define river cross‐sectionsRMSE of 0.676 m and KGE of 0.88 is achieved for a hydrodynamic model of a stretch of the Mekong River using the derived cross‐sectionsSections with significant flow acceleration lead to inaccurate cross‐sections, but errors are reduced with additional ICESat‐2 crossings [ABSTRACT FROM AUTHOR]

Copyright of Journal of Remote Sensing is the property of Editorial Office of Journal of Remote Sensing & Science Publishing Co. and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Accurate bathymetric data is crucial for various aspects such as marine resource exploitation and marine ecological conservation. Currently, satellite-derived bathymetry (SDB) based on empirical and physical models has been widely utilized in constructing underwater terrain in shallow seas. However, the application of such SDB models is limited in remote island reef areas lacking in-situ measurement data. To overcome this issue, the manuscript proposes an unconstrained SDB optimization method without in-situ measurement data, utilizing satellite multispectral imagery (Geoeye-1) and spaceborne LiDAR data (ICESat-2). By classifying the seafloor substrate in coral reef areas into sandy and coral, based on the depth invariant index (DII), we employ an adaptive logarithmic ratio model for unconstrained SDB. The ICESat-2 LiDAR data are then used to correct the SDB results, achieving bathymetry optimization in the coral reef area of the Xisha Islands. Additionally, the proposed method is applied to Yuanzhi Island of the Xisha Islands, and the accuracy of the bathymetric results is evaluated against ALB (Airborne LiDAR Bathymetry) data. The findings demonstrate that compared to conventional methods, our method can improve the accuracy of SDB results with good adaptability. In the Yuanzhi Island area, the proposed method yields SDB results with an R2 of 0.93, an MAE (Mean Absolute Error) of 0.94, and an RMSE (Root Mean Square Error) of 1.12 m, compared to ALB data. The average error is less than 10 % of the maximum depth, essentially meeting the requirements of the International Hydrographic Organization (IHO) standards for depth measurement error when depth is

Supervised classification methods can distinguish between noise and signal in ice, cloud, and land elevation satellite-2 (ICESat-2) data across various feature perspectives and autonomously optimize parameters. Nevertheless, model generalization remains a significant limitation for practical applications. This study focuses on developing a universal denoising model for ICESat-2 using machine learning algorithms and analyzing its spatial transferability under various forest and terrain conditions. A photon-denoising feature parameter system is developed based on the analysis of the three-dimensional distribution of photons in forested regions. This system reduces the parameters dependent on absolute physical quantities and increases those that are less influenced by terrain and forest features to enhance the model’s transferability. Subsequently, automated machine learning algorithms (AutoML) are used for model selection and parameter optimization across six non-parametric regression models. We evaluate the accuracies of the local, global, and transfer models in estimating canopy height across four representative forested areas in China. Results show that the algorithm can effectively distinguish between signal and noise photons. The estimated canopy heights from signal photons are highly consistent with heights obtained using airborne laser scanning (ALS), exhibiting a Pearson correlation coefficient (r) of 0.89, root mean square errors (RMSE) of 3.75 m, relative root mean square error (rRMSE) of 0.27, relative bias (rBias) of −0.11 and mean Bias of −1.45 m. Notably, the accuracy of canopy height estimation by the global model has increased by an average of 21 % compared to ICESat-2 land-vegetation along-track products (ATL08). Furthermore, the model exhibits significant spatial transfer capabilities, with the accuracies of the transfer model exceeding those of ATL08 by margins ranging from 4 % to 41 %. This study marks a significant advancement in photon-denoising methodologies, providing a robust and transferable solution for large-scale environmental data analysis.

Our understanding of the impact of climate change on water availability and natural hazards in high-mountain regions is limited due to the spatial and temporal scarcity of ground observations of precipitation and snow. Freely available, satellite-based information about the snowpack is currently mainly limited to indirect measurements of snow-covered area or very coarse-scale snow water equivalent (SWE), but only for flat areas in lowlands without vegetation cover. Novel space-based laser altimeters, such as ICESat-2, have the potential to provide high-resolution snow depth data in worldwide mountain regions where no ground observations exist. However, these space-based laser altimeters come with spatial gaps between ground tracks, obtained without repetition at a give location. To overcome these drawbacks, here, we present a combined probabilistic data assimilation and deep learning approach to reconstruct spatio-temporal SWE from observations of snow depth along ground tracks, imitating ICESat-2 tracks in view of a potential future global application.Our approach is based on assimilating SWE and snow cover information in a degree-day model with an iterative ensemble smoother (IES) which allows temporally reconstructing SWE along hypothetical ground tracks separated by 3 km. As input, the degree-day model uses daily precipitation and downscaled air temperature from the ERA5 reanalysis. A feedforward neural network (FNN) is then used for spatial propagation of the daily mean and standard deviation of the updated SWE ensemble members obtained from the IES. The combined IES-FNN approach provides uncertainty-aware spatio-temporally continuous estimates of SWE.We tested our approach in the alpine Dischma valley (Switzerland) using high-resolution snow depth maps obtained from photogrammetric techniques mounted on airplanes and unmanned aerial system observations. Our results show that the IES-FNN model provides reliable estimates at a resolution of approximately 100 m. Even assimilating only one SWE observation during the year (combined with satellite-based melt-out date estimates) produces satisfying results when evaluating the IES-FNN SWE reconstructions on independent dates and smaller (

The accuracy of digital elevation models (DEMs) is crucial for practical applications in complex terrain. In this study, various correction models were employed to evaluate and correct the 90 m SRTM3 DEM and TanDEM-X DEM data in the Qinghai-Tibet Plateau region. A simple and universal three-dimensional data co-registration method was utilized for horizontal alignment. The study demonstrated a significant reduction in horizontal discrepancies between datasets after data co-registration. By comparing the performance of five different correction methods, it was found that the RF model exhibited excellent accuracy and robustness, making it particularly suitable for applications that require high precision. The XGB and BPNN models offer a good balance between accuracy and time efficiency, making them suitable for research scenarios that are less sensitive to precision requirements. In contrast, the MLR and IDW methods performed poorly and are not recommended for high-precision applications. Furthermore, the corrected DEMs showed improved resistance to the influence of high slopes, enhancing their applicability in complex terrain areas. This study provides an important reference for the correction of high-quality DEMs in the Qinghai-Tibet Plateau region and offers valuable insights for DEM research in similar areas.

The ICESat-2 satellite equipped with a new photon-counting laser altimeter has received much attention as a source of accurate elevation observations. However, in this research field, there is a lack of an open-source high-accuracy elevation control point dataset with the specific quality requirements at a global scale. To this end, using ICESat-2 altimeter data as the main data source, we constructed and organized a dataset as a useful supplement for this research field. The dataset was generated by a methodology based on detection environment evaluation, photon spatial analysis, and the redundant observation statistics. The dataset includes more than 600 million elevation control points and covers the global land areas, except for Greenland and Antarctica. The dataset has been validated by multiple digital elevation models (DEMs) from around the world (sourced from airborne LiDAR data). The results show that the dataset has high-accuracy elevation control points. The overall root-mean-square error (RMSE) of the original elevations of ICESat-2 is about 1.384–4.820 m, but the overall RMSE of the elevation control points in the new dataset is about 0.279–0.642 m. Moreover, the results obtained in this study show that the dataset is suitable for application within high vegetation cover areas.

ABSTRACTTo tackle the challenge of denoising spaceborne photon-counting laser altimeter point clouds with uneven noise density, this study proposes a denoising method based on adaptive parameter density clustering, which utilizes numerical simulations to achieve self-adaptation of key parameters (neighborhood radius [Formula: see text] and minimum number of points [Formula: see text]). First, taking the directional adaptive ellipse DBSCAN (DAE-DBSCAN) as an example, photons with different background photon count rates ([Formula: see text]) are used to traverse [Formula: see text] and [Formula: see text] to calculate their optimal values ([Formula: see text] and [Formula: see text] with the highest denoising accuracy). Then, a mathematical prediction model of [Formula: see text], [Formula: see text] and [Formula: see text] was established. The actual background photon count rates were introduced into the key parameter prediction model to obtain the optimal [Formula: see text] and [Formula: see text]. Finally, a denoising experiment was conducted using the simulated photons and the ATLAS data. The results show that the proposed method had higher accuracy than the constant parameter denoising method, with an [Formula: see text] >0.95. Even for photons of complex mountainous terrain with a high background photon count rate, the denoising accuracy was still higher than 0.9. The proposed method improves the denoising accuracy of photons with different noise densities by adapting density clustering parameters.

ABSTRACTOpen source Global Digital Elevation Models (GDEMs) serve as an important base for studies in geosciences. However, these models contain vertical errors due to various reasons. In this study, data from two Satellite LiDAR altimetry systems, GEDI and ICESat-2, were used to improve the vertical accuracy of GDEMs. Three different machine learning methods, namely an Artificial Neural Network (ANN), Extreme Gradient Boosting (XGBoost), and a Convolutional Neural Network (CNN), were employed to improve existing DEM data with satellite LiDAR data. The methodology was tested in five areas with varying characteristics. Ground control data were selected from high accuracy DEMs generated from Airborne LiDAR and GNSS data. The use of ANN method improved the vertical accuracy of SRTM data from 6.45 to 3.72 m in Test area-4. Similarly, the CNN method demonstrated an improvement in the vertical accuracy of bare ground SRTM data increasing from 3.4 to 0.6 m in Test area-4. In Test area-5, the ANN method improved the vertical accuracy of SRTM data with slopes between 30 and 60%, increasing from 3.8 to 0.5 m. Notably, the results underscore the successful improvement of GDEMs across all test areas.

The forest structure of Canada's northern unmanaged boreal forests are not well characterized despite their critical importance for ecosystem services such as carbon storage and caribou habitat, as well as their vulnerability to climate change and increasing disturbance rates, especially wildfire. Spaceborne lidar observations acquired from the Advanced Topographic Laser Altimeter System (ATLAS) instrument onboard the Ice, Cloud and Elevation Satellite-2 (ICESat-2), enables monitoring of forest structure in these remote northern ecosystems and characterization of long-term post-fire recovery; however, ICESat-2 is a sampling instrument, acquiring data transects rather than wall-to-wall data coverage. By imputing ICESat-2 estimates of canopy height using annual time-series (1984–2021) of spectral indices, change metrics, and land cover metrics from Landsat surface reflectance composites, we derived annual, spatially explicit, wall-to-wall estimates of canopy height for the period 1984–2021 across a 19.6 Mha area of boreal forests in northwestern Ontario and a 630,000 ha area of managed boreal forest further south that had coincident airborne laser scanning (ALS) data for validation. The accuracy of the imputation of derived canopy height estimates in the northwestern site was assessed using a reserved set of ICESat-2 observations (r = 0.76; RMSD = 3.04 m, pRMSD = 33.9%, MD = −0.10 m, and pMD = −1.0%). Using coincident ALS data, the accuracy of the imputed canopy heights in the southern study site were also assessed (r = 0.81; RMSD = 3.65 m, pRMSD = 25.87%), MD = −0.70 m, pMD = 4.96%). Examination of height dynamics post fire indicated that canopy height decreased 8–10 years post fire and recovered between 60% and 85% of pre-fire canopy height within 30 years post fire. The approach presented could readily be extended to similar northern boreal forest areas to provide broad-scale characterizations of fire impacts on forest structure and subsequent recovery. [ABSTRACT FROM AUTHOR]

Inland water level and its dynamics are key components in the global water cycle and land surface hydrology, significantly influencing climate variability and water resource management. Satellite observations, in particular altimetry missions, provide inland water level time series for nearly three decades. Space-based remote sensing is regarded as a cost-effective technique that provides measurements of global coverage and homogeneous accuracy in contrast to in-situ sensors. The advent of Open-Loop Tracking Command (OLTC), and Synthetic Aperture Radar (SAR) mode strengthened the use of altimetry missions for inland water level monitoring. However, it is still very challenging to obtain accurate measurements of water level over narrow rivers and small lakes. This scoping systematic literature review summarizes and disseminates the research findings, highlights major results, and presents the limitations regarding inland water level monitoring from satellite observations between 2018 and 2022. Following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guideline and through a double screening process, 48 scientific publications were selected meeting the eligibility criteria. To summarize the achievements of the previous 5 years, we present fundamental statistical results of the publications, such as the annual number of publications, scientific journals, keywords, and study regions per continent and type of inland water body. Also, publications associated with specific satellite missions were analyzed. The findings show that Sentinel-3 is the dominant satellite mission, while the ICESat-2 laser altimetry mission has exhibited a high growth trend. Furthermore, publications including radar altimetry missions were charted based on the retracking algorithms, presenting the novel and improved methods of the last five years. Moreover, this review confirms that there is a lack of research on the collaboration of altimetry data with machine learning techniques. [ABSTRACT FROM AUTHOR]

TanDEM-X bistatic interferometric system successfully generated a high-precision, high-resolution global digital elevation model (DEM). However, in forested areas, two core problems make it difficult to obtain sub-canopy topography: (1) the penetrability of short-wave signals is limited, and the DEM obtained in dense forest areas contains a significant forest signal, that is, the scattering phase center (SPC) height; and (2) the single-baseline and single-polarization TanDEM-X interferometric synthetic aperture radar (InSAR) data cannot provide sufficient observations to make the existing physical model reversible for estimating the real surface phase, whereas the introduction of optical data makes it difficult to ensure data synchronization and availability of cloud-free data. To overcome these problems in accurately estimating sub-canopy topography from TanDEM-X InSAR data, this study proposes a practical method of sub-canopy topography estimation based on the following innovations: (1) An orthogonal polynomial model was established using TanDEM-X interferometric coherence and slope to estimate the SPC height. Interferometric coherence records forest height and dielectric property information from an InSAR perspective and has spatiotemporal consistency with the InSAR-derived DEM. (2) Introduce Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) data to provide more observational information and automatically screen ICESat-2 control points with similar forest and slope conditions in the local area to suppress forest spatial heterogeneity. (3) A weighted least squares criterion was used to solve this model to estimate the SPC height. The results were validated at four test sites using high-precision airborne light detection and ranging (LiDAR) data as a reference. Compared to the InSAR-derived DEM, the accuracy of the sub-canopy topography was improved by nearly 60%, on average. Furthermore, we investigated the necessity of local modeling, confirming the potential of the proposed method for estimating sub-canopy topography by relying only on TanDEM-X and ICESat-2 data. [ABSTRACT FROM AUTHOR]

The highly accurate data of topography and bathymetry are fundamental to ecological studies and policy decisions for coastal zones. Currently, the automatic extraction and classification of signal photons in coastal zones is a challenging problem, especially the surface type classification without auxiliary data. The lack of classification information limits large-scale bathymetric applications of ICESat-2 (Ice, Cloud, and Land Elevation Satellite-2). In this study, we propose a photon extraction–classification method to process geolocated photons in coastal areas from the ICESat-2 ATL03 product. The basic idea is to extract the signal photons using an adaptive photon clustering algorithm, and the extracted signal photons are classified based on the accumulated histogram and triangular grid. We also generate the bottom profile using the weighted interpolation. In four typical coastal areas (artificial coast, natural coast, island, and reefs), the extraction accuracy of a signal photons exceeds 0.90, and the Kappa coefficients of four surface types exceed 0.75. This method independently extracts and classifies signal photons without relying on auxiliary data, which can greatly improve the efficiency of obtaining bathymetric points in all kinds of coastal areas and provide technical support for other coastal studies using ICESat-2 data. [ABSTRACT FROM AUTHOR]

In the summer of 2020, ESA changed the orbit of CryoSat‐2 to align periodically with NASA's ICESat‐2 mission, a campaign known as CRYO2ICE, which allows for near‐coincident CryoSat‐2 and ICESat‐2 observations in space and time over the Arctic until summer 2022, where the CRYO2ICE Antarctic campaign was initiated. This study investigates the Arctic CRYO2ICE radar and laser freeboards acquired by CryoSat‐2 and ICESat‐2, respectively, during the winter seasons of 2020–2022, and derives snow depths from their differences along the orbits. Along‐track snow depth observations can provide high‐resolution snow depth distributions which are vital for air‐ice‐ocean heat and momentum transfer, understanding light transmission, and snow‐ice‐interactions. Generally, ICESat‐2 is backscattered at a surface above the elevation of the CryoSat‐2 signal. CRYO2ICE snow depths are thinner than the daily model‐ or passive‐microwave‐based snow depth composites used for comparison, with differences being most pronounced in the Atlantic and Pacific Arctic. Satellite‐derived and model‐based snow estimates show similar seasonal accumulation over first‐year ice, but CRYO2ICE has limited seasonal accumulation over multi‐year ice which is linked to a slow increase in ICESat‐2, and to some extent CryoSat‐2, freeboards. We present a first estimation of spaceborne along‐track snow depth estimates with average uncertainty of 10–11 ± 2–3 cm for 7‐km segments, with random and systematic contributions of 7 and 4 cm. These observations show the potential for along‐track dual‐frequency observations of snow depth from the future Copernicus mission CRISTAL; but they also highlight uncertainties in radar penetration and the correlation length scales of snow topography that still require further research. Plain Language Summary: Estimates of snow depth on sea ice are currently either outdated or limited in resolution, thus a need to derive high‐resolution snow depth on sea ice is crucial. Former studies have computed snow depth on sea ice using the difference in penetration between radar (Ku‐band) and laser (or Ka‐band radar) altimeters from different satellite missions, assuming the Ku‐band radar and laser/Ka‐band are reflected at the bottom and top of the snow pack, respectively. Those studies have resulted in monthly composites of snow depth due to a limited overlap of individual tracks between the satellite missions. Since the CRYO2ICE (CryoSat‐2/ICESat‐2 Resonance) campaign was initiated in July 2020, we have for the first time the possibility of investigating near‐coincident observations of spaceborne radar and laser altimeters. These results are important to initiate relevant discussions on snow depth retrieval in preparation for the future dual‐frequency altimeter mission, CRISTAL (Copernicus Polar Ice and Snow Topography Altimeter), expected to launch in 2028. Key Points: CRYO2ICE (C2I) orbits are investigated over sea ice for the first time and snow depth is computed along the orbitsC2I snow depth has an uncertainty of 10–11 cm along 7 km segments with 5–95% confidence interval of 7–16 cm for both winter seasonsAlong‐track snow depth is compared with modeled and passive‐microwave estimates and C2I is thinner for the Pacific and Atlantic Arctic [ABSTRACT FROM AUTHOR]

High-quality height reference data are embedded in the accuracy verification processes of most remote sensing terrain applications. The Ice, Cloud, and Land elevation Satellite 2 (ICESat-2)/ATL08 terrain product has shown promising results for estimating ground heights, but it has not been fully evaluated. Hence, this study aims to assess and enhance the accuracy of the ATL08 terrain product as a height reference for the newest versions of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), and TanDEM-X (TDX) DEMs over vegetated mountainous areas. We used uncertainty-based filtering method for the ATL08 strong and weak beams to enhance their accuracy. Then, the results were evaluated against a reference airborne LiDAR digital terrain model (DTM), by selecting 10,000 points over the entire area and comparing the accuracy of ASTER, SRTM, and TDX DEMs assessed by the LiDAR DTM to the accuracy of the ASTER, SRTM, and TDX DEMs assessed by the ATL08 strong beams, weak beams, and all beams. We also detected the impact of the terrain aspect, slope, and land cover types on the accuracy of the ATL08 terrain elevations and their relationship with height errors and uncertainty. Our findings show the accuracy of the ATL08 strong beams was enhanced by 43.91%; while the weak beams accuracy was enhanced by 74.05%. Furthermore, slope strongly influenced ATL08 height errors and height uncertainty; especially on the weak beams. The errors induced by the slope significantly decreased when the uncertainty levels were reduced to <20 m. The evaluations of ASTER, SRTM, and TDX DEMs by ATL08 strong and weak beams are close to those assessed by LiDAR DTM points within 0.6 m for the strong beams. These findings indicate that ATL08 strong beams can be used as a height reference over vegetated mountainous regions. [ABSTRACT FROM AUTHOR]

The new-era space-borne LiDAR systems, viz. ICESat-2 and GEDI, offer new possibilities for mapping terrain and canopy heights through geophysical data products, viz ATL08 (ICESat-2) and L2A (GEDI). Additionally, GEDI provides an above-ground biomass density (AGBD) product (L4A) derived through parametric calibration of L2A metrics. Detailed comparisons of these data products among different forest types and other influencing factors (viz. acquisition parameters and ground conditions) are essential for continued improvement and broader use. In this study, we perform a detailed accuracy assessment of these data products over tropical dry deciduous forests in the Central Indian region during leaf-off and leaf-on seasons with the reference airborne LiDAR data. The GEDI L4A (AGBD) product is validated against the reference AGBD map derived from the field AGBD estimates and canopy height model from airborne LiDAR. Our results suggest that regardless of leaf condition, strong or power beams during nights from both systems (GEDI and ICESat-2) are highly capable of retrieving terrain height with RMSE 2.8–3.8 m and low bias of − 0.2 m to + 1.4 m. Nevertheless, GEDI canopy height retrievals were strongly linked to leaf availability with significant underestimations (bias > − 5.8 m) during the leaf-off season against a bias of ± 1 m during the leaf-on season. In contrast, strong night beams from ICESat-2 were found to retrieve canopy heights accurately with a bias of − 0.4 m and RMSE of 3.5 m during the leaf-off season. The GEDI-AGBD estimates were found to be substantially mismatched with the reference AGBD map with an overall RMSE of 46%. The atmospheric conditions and topographic slope strongly influenced the accuracy of both systems. [ABSTRACT FROM AUTHOR]

The southeastern Tibetan Plateau (SETP), which hosts the most extensive marine glaciers on the Tibetan Plateau (TP), exhibits enhanced sensitivity to climatic fluctuations. Under global warming, persistent glacier mass depletion within the SETP poses a risk to water resource security and sustainability in adjacent nations and regions. This study deployed a high-precision ICESat-2 satellite altimetry technique to evaluate SETP glacier thickness changes from 2018 to 2022. Our results show that the average change rate in glacier thickness in the SETP is −0.91 ± 0.18 m/yr, and the corresponding glacier mass change is −7.61 ± 1.52 Gt/yr. In the SETP, the glacier mass loss obtained via ICESat-2 data is larger than the mass change in total land water storage observed by the Gravity Recovery and Climate Experiment follow-on satellite (GRACE-FO), −5.13 ± 2.55 Gt/yr, which underscores the changes occurring in other land water components, including snow (−0.44 ± 0.09 Gt/yr), lakes (−0.06 ± 0.02 Gt/yr), soil moisture (1.88 ± 1.83 Gt/yr), and groundwater (1.45 ± 0.70 Gt/yr), with a closure error of −0.35 Gt/yr. This demonstrates that this dramatic glacier mass loss is the main reason for the decrease in total land water storage in the SETP. Generally, there are decreasing trends in solid water storage (glacier and snow) against stable or increasing trends in liquid water storage (lakes, soil moisture, and groundwater) in the SETP. This persistent decrease in solid water is linked to the enhanced melting induced by rising temperatures. Given the decreasing trend in summer precipitation, the surge in liquid water in the SETP should be principally ascribed to the increased melting of solid water. [ABSTRACT FROM AUTHOR]

Satellite data have exhibited a great potential to track lake water-level changes. By integrating multi-source satellite datasets (i.e., Ice, Cloud and land Elevation Satellite-2 (ICESat-2) lidar datasets, Landsat imagery, and Global Surface Water Dataset (GSWD)), we proposed a method to track annual water-level changes of Lake Mead, USA over 1984–2018. Specifically, we first extracted ICESat-2 profiles as elevation information. Second, we used Landsat imagery to extract annual lake boundaries and calculate lake areas. Third, we constructed the relationship between elevations and inundation frequencies provided by the GSWD for ICESat-2 points, and applied it to extract annual water levels by using inundation frequency of Landsat derived lake boundaries. Finally, we created the relationship between lake water level and lake areas via a cubic polynomial model to refine the annual lake water levels. Our water level results agreed well with the in-situ data (R2 = 0.99, p < 0.01, RMSE = 0.64 m). Our results showed that the ICESat-2 data performed well in providing accurate elevation information around the shore of lake. Using satellite datasets, our proposed method has the potential to derive water levels for water bodies worldwide in the future. [ABSTRACT FROM AUTHOR]