ICESat Research Articles

Exploring the Relationship between Forest Canopy Height and Canopy Density from Spaceborne LiDAR Observations

Forest structure is a useful proxy for carbon stocks, ecosystem function and species diversity, but it is not well characterised globally. However, Earth observing sensors, operating in various modes, can provide information on different components of forests enabling improved understanding of their structure and variations thereof. The Ice, Cloud and Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS), providing LiDAR footprints from 2003 to 2009 with close to global coverage, can be used to capture elements of forest structure. Here, we evaluate a simple allometric model that relates global forest canopy height (RH100) and canopy density measurements to explain spatial patterns of forest structural properties. The GLA14 data product (version 34) was applied across subdivisions of the World Wildlife Federation ecoregions and their statistical properties were investigated. The allometric model was found to correspond to the ICESat GLAS metrics (median mean squared error, MSE: 0.028; inter-quartile range of MSE: 0.022–0.035). The relationship between canopy height and density was found to vary across biomes, realms and ecoregions, with denser forest regions displaying a greater increase in canopy density values with canopy height, compared to sparser or temperate forests. Furthermore, the single parameter of the allometric model corresponded with the maximum canopy density and maximum height values across the globe. The combination of the single parameter of the allometric model, maximum canopy density and maximum canopy height values have potential application in frameworks that target the retrieval of above-ground biomass and can inform on both species and niche diversity, highlighting areas for conservation, and potentially enabling the characterisation of biophysical drivers of forest structure.

TanDEM-X PolarDEM 90 m of Antarctica: generation and error characterization

Abstract. We present the generation and validation of an updated version of the TanDEM-X digital elevation model (DEM) of Antarctica: the TanDEM-X PolarDEM 90 m of Antarctica. Improvements compared to the global TanDEM-X DEM version comprise filling gaps with newer bistatic synthetic aperture radar (SAR) acquisitions of the TerraSAR-X and TanDEM-X satellites, interpolation of smaller voids, smoothing of noisy areas, and replacement of frozen or open sea areas with geoid undulations. For the latter, a new semi-automatic editing approach allowed for the delineation of the coastline from DEM and amplitude data. Finally, the DEM was transformed into the cartographic Antarctic Polar Stereographic projection with a homogeneous metric spacing in northing and easting of 90 m. As X-band SAR penetrates the snow and ice pack by several meters, a new concept for absolute height adjustment was set up that relies on areas with stable penetration conditions and on ICESat (Ice, Cloud, and land Elevation Satellite) elevations. After DEM generation and editing, a sophisticated height error characterization of the whole Antarctic continent with ICESat data was carried out, and a validation over blue ice achieved a mean vertical height error of just −0.3 m ± 2.5 m standard deviation. The filled and edited Antarctic TanDEM-X PolarDEM 90 m is outstanding due to its accuracy, homogeneity, and coverage completeness. It is freely available for scientific purposes and provides a high-resolution data set as basis for polar research, such as ice velocity, mass balance estimation, or orthorectification.

Ground-based corner cube retroreflector optimal design for the geolocation validation of satellite photon counting lidar.

The corner cube retroreflector (CCR) can be applied in the on-orbit geolocation validation of satellite photon counting lidar, which provides insight into the fidelity of laser footprint geolocation on the ground. A novel parameter optimization method, to the best of our knowledge, is proposed to enhance the performance of the CCR. Based on the velocity aberration effect of the satellite and the far-field diffraction pattern (FFDP) of the CCR, a mathematical model with respect to the received energy and design parameters of the CCR, including the aperture, the dihedral angle errors (DAE), and the curvature radius of the front face (CROF), is derived. We can achieve the optimal design parameters of the CCR through searching for the locations of the maximal and uniform FFDP on the circumference of the aberration position. We resolve the optimal aperture, the DAE, and the CROF of the CCR as 14.3 mm, 2.4 in., and 6000 m for the Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) lidar. The received energy derived from the derived parameters has a significant improvement of around 26 times compared with the CCR parameters employed in the ICESat-2. The results demonstrate that the proposed algorithm is effective, and the regression models of the optimal aperture and DAE for different satellite altitudes are conveniently employed in the design of the CCR for the geolocation validation of satellite photon counting lidar.

Satellite-based estimation of roughness lengths and displacement heights for wind resource modelling

Abstract. Wind turbines in northern Europe are frequently placed in forests, which sets new wind resource modelling requirements. Accurate mapping of the land surface can be challenging at forested sites due to sudden transitions between patches with very different aerodynamic properties, e.g. tall trees, clearings, and lakes. Tree growth and deforestation can lead to temporal changes of the forest. Global or pan-European land cover data sets fail to resolve these forest properties, aerial lidar campaigns are costly and infrequent, and manual digitization is labour-intensive and subjective. Here, we investigate the potential of using satellite observations to characterize the land surface in connection with wind energy flow modelling using the Wind Atlas Analysis and Application Program (WAsP). Collocated maps of the land cover, tree height, and leaf area index (LAI) are generated based on observations from the Sentinel-1 and Sentinel-2 missions combined with the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2). Three different forest canopy models are applied to convert these maps to roughness lengths and displacement heights. We introduce new functionalities for WAsP, which can process detailed land cover maps containing both roughness lengths and displacement heights. Validation is carried out through cross-prediction analyses at eight well-instrumented sites in various landscapes where measurements at one mast are used to predict wind resources at another nearby mast. The use of novel satellite-based input maps in combination with a canopy model leads to lower cross-prediction errors of the wind power density (rms = 10.9 %–11.2 %) than using standard global or pan-European land cover data sets for land surface parameterization (rms = 14.2 %–19.7 %). Differences in the cross-predictions resulting from the three different canopy models are minor. The satellite-based maps show cross-prediction errors close to those obtained from aerial lidar scans and manually digitized maps. The results demonstrate the value of using detailed satellite-based land cover maps for micro-scale flow modelling.

A Method of Extracting High-Accuracy Elevation Control Points from ICESat-2 Altimetry Data

Due to its high ranging accuracy, spaceborne laser altimetry technology can improve the accuracy of satellite stereo mapping without ground control points. In the past, full-waveform ICE, CLOUD, and Land Elevation Satellite (ICESat) laser altimeter data have been used as one of the main data sources for global elevation control. As a second-generation satellite, ICESat-2 is equipped with an altimeter using photon counting mode. This can further improve the application capability for stereo mapping because of the six laser beams with high along-track repetition frequency, which can provide more detailed ground contour descriptions. Previous studies have addressed how to extract high-accuracy elevation control points from ICESat data. However, these methods cannot be directly applied to ICESat-2 data because of the different modes of the laser altimeters. Therefore, in this paper, we propose a method using comprehensive evaluation labels that can extract high-accuracy elevation control points that meet the different level elevation accuracy requirements for large scale mapping from the ICESat-2 land-vegetation along-track product. The method was verified using two airborne lidar data sets. In flat, hilly, and mountainous areas, by using our method to extract the terrain elevation, the root-mean-square error of elevation control points decrease from 1.249–2.094 m, 2.237–3.225 m, and 2.791–4.822 m to 0.262–0.429 m, 0.484–0.596 m, and 0.611–1.003 m, respectively. The results show that the extraction elevations meet the required accuracy for large scale mapping.

Accurate Refraction Correction—Assisted Bathymetric Inversion Using ICESat-2 and Multispectral Data

Shallow-water depth information is essential for ship navigation and fishery farming. However, the accurate acquisition of shallow-water depth has been a challenge for marine mapping. Combining Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) bathymetry data with multispectral data, satellite-derived bathymetry is a promising solution through which to obtain bathymetric information quickly and accurately. This study proposes a photon refraction correction method considering sea-surface undulations to address errors in the underwater photons obtained by the ICESat-2. First, the instantaneous sea surface and beam emission angle are integrated to determine the sea-surface incidence angle. Next, the distance of photon propagation in water is determined using sea-surface undulation and Snell’s law. Finally, position correction is performed through geometric relationships. The corrected photons were combined with the multispectral data for bathymetric inversion, and a bathymetric map of the Yongle Atoll area was obtained. A bathymetric chart was created using the corrected photons and the multispectral data in the Yongle Atoll. Comparing the results of different refraction correction methods with the data measured shows that the refraction correction method proposed in this paper can effectively correct bathymetry errors: the root mean square error is 1.48 m and the R2 is 0.86.

Retrieving building height in urban areas using ICESat-2 photon-counting LiDAR data

Building height is a pivotal factor in studying urban form and understanding the impacts of the vertical characteristics of urban areas on the environment. The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) employ photon-counting LiDAR to collect Earth’s surface elevation data globally. However, method at ground and building photon classification and subsequently building height extraction from photon-counting LiDAR data is still missing. This study aimed to develop a methodological framework to retrieve building height from ICESat-2 data under various data acquisition conditions stratified by day/night condition, strong/weak beam. First, the noise removal algorithm adopted by ATL03 product is utilized to filter out noise photons and retain signal photons in the raw dataset. Second, the random sample consensus (RANSAC) algorithm and statistical method are adopted to classify the signal photons into ground and building photons according to their characteristics and spatial neighborhood relationship. The building height is then computed by the elevation of building and ground photons. Finally, estimated heights are evaluated using the reference building height derived from terrestrial laser scanning (TLS) data. The results indicate that the proposed methodological framework can effectively identify building and ground photons under different building types. Strong consistence between estimated and reference building height is obtained by quantifying the estimation precision of building height under different data acquisition conditions, root mean square error (RMSE) ranging between 0.35 m and 0.45 m, indicating that the ATL03 data can be utilized to derive the building height in urban areas for all acquisition times and beam intensity. Further analysis demonstrates that strong/weak beam significantly influences the building height estimation compared to day/night condition. Overall, this study provides a method for estimating building height in urban areas using ICESat-2 data, and the findings approve the strong capability of ICESat-2 data in estimating building height.

The Prospect of Global Coral Reef Bathymetry by Combining Ice, Cloud, and Land Elevation Satellite-2 Altimetry With Multispectral Satellite Imagery

Empirical methods for estimating shallow-water bathymetry using passive multispectral satellite imagery are robust and globally applicable, in theory, but they require copious local measurements of water depth for algorithm calibration. Such calibration data have historically been unavailable for most locations, but NASA’s Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2), a satellite-based LiDAR, might hold unique promise to fill this critical data gap. Although ICESat-2 was not designed as a marine altimeter, its ATLAS sensor consists of six green (532 nm) lasers that can penetrate a water surface and return photons reflected by the seabed, thereby generating bathymetric profiles. Utilizing TCarta’s NSF SBIR-funded Space-Based Laser Bathymetry Extraction Tool and ICESat-2’s ATL03 geolocated photon data product, we have compared ICESat-2 bathymetric retrievals with a portfolio of soundings acquired in situ using a vessel-mounted single-beam echosounder. This analysis demonstrated very high correlation (R2 = 0.96) between the field and space-based bathymetry data. The comparisons were made at multiple Caribbean and Pacific coral reef sites over water depths ranging from 1 to 20 m. Results suggest that ICESat-2 could be an effective approach for calibrating and validating empirical and radiative transfer methods, alike, for estimating shallow-water bathymetry from remote sensing imagery, thereby enabling the immediate potential for shallow-water bathymetric mapping of Earth’s reefs.

Deriving Antarctic Sea‐Ice Thickness From Satellite Altimetry and Estimating Consistency for NASA's ICESat/ICESat‐2 Missions

AbstractWith the launch of the Ice, Cloud and Land Elevation Satellite‐2 (ICESat‐2), densely measuring the Antarctic ocean all year long, monthly sea ice thickness changes can be inspected. However, the consistency between ICESat‐2 and its predecessor ICESat remains unknown if a long‐term change of the Antarctic sea ice thickness is to be examined. Due to the complex physical structure of Antarctic sea ice, an improved One‐Layer Method (OLMi) for sea‐ice thickness retrieval, with an ∼0.3 m uncertainty, is proposed in this study. Using this method, the monthly sea ice variation in each sea sector of Antarctica from ICESat‐2 is inspected, revealing detailed characteristics of Antarctic sea ice over sector scale such as bi‐modal distributions. The different retrieving methods for the total freeboard from ICESat are found to be ∼0.036 m lower than that from ICESat‐2. The circum‐Antarctic sea ice thickness shows a decreasing tendency (2003–2020), especially in summer with −0.014 m/yr.

Estimation of thermodynamic and dynamic contributions to sea ice growth in the Central Arctic using ICESat-2 and MOSAiC SIMBA buoy data

The fine spatial resolution of the ICESat-2 (IS2) satellite altimeter allows monitoring the evolution of sea ice thickness with detailed dynamic information (e.g. ridges and leads). In this study, we first assess the ability of IS2 to estimate thermodynamic ice growth and dynamic thickening during the ice-growing season in the central Arctic Ocean. As an indicator of the thermodynamic ice growth, we use 10 thermistor string-based sea ice mass balance array (SIMBA) buoys deployed at a scale of ~50 km from the Icebreaker Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. We collect IS2 data within 20 km buffer distance from the individual buoys, and calculate the mode, median, and mean of the IS2-derived ice thickness. The IS2 modal thickness shows the least bias (−0.169 m) with the buoy ice thickness, representing level ice thickness. In addition, the increasing rate of the IS2 modal thickness is close to the thermodynamic ice growth with a small bias of −0.054 cm/day. However, the increasing rates of the IS2 median and mean thickness are greater than the thermodynamic ice growth by about 0.114 cm/day and 0.198 cm/day, respectively, because they also include ice growth caused by thickness redistribution during dynamic deformation. The dynamic contributions may account for 26.1 ± 10.3% and 34.4 ± 10.1% of the total increase of the IS2 median and mean thickness, respectively. Within a ~ 50 km radius area from the MOSAiC Central Observatory, IS2 measurements exhibit that the ridge fraction increased from <2% in November to ~4% in March (~0.029%/day of average increasing rate) and ridge height increased about 0.047 cm/day during the same period. However, lead formation does not show significant contributions to the dynamic ice thickening because leads are temporary features lasting only 2–3 days. Although there are considerable uncertainties in IS2 ice thickness estimation and IS2-buoy thickness comparison, this study emphasizes the importance of combining measurements by IS2 and SIMBA buoys to explain the regional sea ice mass balance with separating the thermodynamic and dynamic contributions.

Toward geomorphic mapping of reef habitat by laser altimeter on ICESat-2 Satellite

The main objectives of Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) are to quantify the change of ice sheets, sea ice, and their behavior. It uses precise laser measurements of surface elevation to achieve the goals. Using 532 nm laser pulse, ICESat-2 has been providing other data of Earth System components as well such as vegetation canopy height, land height, and nearshore bathymetry. This study aims were to investigate the potential use of laser photon data of ICESat-2 for geomorphic mapping and reef rugosity indication. About eight tracks of data photon of Kapoposang island and four tracks of Kodingarenglompo were investigated to indicate the existence of coral reef and its geomorphic zone. The results show that the photon’s depth for Kapoposang island is easier to interpret than Kadingarenglompo due to the clarity of the water column. The deepest depth photon achieves for Kapoposang reef was up to 17 meters while at Kodingarenglompo just about 6 meters. The depth photon variability at reef area is significantly different from sand and submerged vegetation area. The depth dan slope variability are easy to detect using this lidar altimetry which in turn very helpful on classifying the geomorphic zone of reef habitat. Hence these results may lead to study more about the use of ICESat-2 data for reef habitat mapping.

Assessing the agreement of ICESat-2 terrain and canopy height with airborne lidar over US ecozones

Despite its critical importance to carbon storage modeling, forest vertical structure remains poorly characterized over large areas. Canopy height estimates from current satellite missions such as ICESat-2 (Ice, Cloud, and land Elevation Satellite-2) offer promise to close this knowledge gap, but their validation is critically important to inform their measurement uncertainties and scientific utility. Using existing airborne laser scanning (ALS) data, the agreement of a variety of terrain and aboveground canopy height metrics including summary height statistics and percentiles, from ICESat-2’ Land, Water and Vegetation Elevation product (ATL08) product was assessed in 12 sites across six major biomes in the United States. The agreement between ATL08 and ALS heights was assessed using the mean bias (Bias, ATL08 – ALS), the mean absolute error (MAE) and their percent equivalents, percent bias (pBias) and percent MAE (pMAE), respectively. In general, the agreement between ATL08 and ALS terrain heights was high (Bias 0.18 m, pBias 0.1%) while canopy heights showed lower agreement (Bias −1.71 m, pBias −15.9%). Analyses by biome, time of acquisition and beam strength of the ICESat-2 photon data also showed generally higher agreement for ATL08 terrain than canopy heights. Analyses also showed the performance of ATL08 heights varied with canopy cover with ATL08 terrain heights showing the best agreement when canopy cover was between 40 and 70% while the best performance for ATL08 canopy heights was observed when canopy cover was greater than 80%. This observation, coupled with analyses by biome, indicate that ATL08 canopy heights are more suitable in relatively dense canopy environments such as conifer and broadleaf forests than relatively sparse environments such a temperate grassland and Savannas. Higher level canopy height percentiles (95th and 98th) showed higher agreement (mean Bias −12.5%) with ALS heights than lower percentiles (minimum, 25th, mean pBias ~39.2%). These findings indicate that ATL08 canopy heights show more promise for routine canopy height characterization using the 95th and 98% percentiles but is limited in characterizing intermediate vertical structure. The observed performance differences between ATL08 terrain and canopy heights are attributed to differences in photon sampling rates over terrain and canopy surfaces which, compounded with background noise in ICESat-2 photon data, led to different effectiveness for ATL08 processing routines in filtering terrain and off-terrain points. This assessment of the impact of a variety of factors provides the vegetation community with an understanding of the capabilities and limitations of height estimates from the ICESat-2 ATL08 product.

Planetary Boundary Layer Height Estimates From ICESat-2 and CATS Backscatter Measurements

The lowest layer of the atmosphere in which all human activity occurs is called the Planetary Boundary Layer (PBL). All physical interactions with the surface, such as heat and moisture transport, pollution dispersion and transport happen in this relatively shallow layer. The ability to understand and model the complex interactions that occur in the PBL is very important to air quality, weather prediction and climate modeling. A fundamental and physically important property of the PBL is its thickness or height. This work presents two methods to obtain global PBL height using satellite lidar data from the Ice, Cloud and land Elevation Satellite-2 (ICESat-2) and the Cloud-Aerosol Transport System (CATS). The first method is a straightforward backscatter threshold technique and the second is a machine learning approach known as a Convolutional Neural Network. The PBL height retrievals from the two methods are compared with each other and with PBL height from the NASA GEOS MERRA-2 reanalysis. The lidar-retrieved PBL heights have a high degree of spatial correlation with the model heights but are generally higher over ocean (∼400 m) and over northern hemisphere high latitude regions (∼1,000 m). Over mid-latitude and tropical land areas, the satellite estimated PBL heights agree well with model mid-day estimates. This work demonstrates the feasibility of using satellite lidar backscatter measurements to obtain global PBL height estimates, as well as determining seasonal and regional variability of PBL height.

Penetration of interferometric radar signals in Antarctic snow

Abstract. Synthetic aperture radar interferometry (InSAR) is an efficient technique for mapping the surface elevation and its temporal change over glaciers and ice sheets. However, due to the penetration of the SAR signal into snow and ice, the apparent elevation in uncorrected InSAR digital elevation models (DEMs) is displaced versus the actual surface. We studied relations between interferometric radar signals and physical snow properties and tested procedures for correcting the elevation bias. The work is based on satellite and in situ data over Union Glacier in the Ellsworth Mountains, West Antarctica, including interferometric data of the TanDEM-X mission, topographic data from optical satellite sensors and field measurements on snow structure, and stratigraphy undertaken in December 2016. The study area comprises ice-free surfaces, bare ice, dry snow and firn with a variety of structural features related to local differences in wind exposure and snow accumulation. Time series of laser measurements of NASA's Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 show steady-state surface topography. For area-wide elevation reference we use the Reference Elevation Model of Antarctica (REMA). The different elevation data are vertically co-registered on a blue ice area that is not affected by radar signal penetration. Backscatter simulations with a multilayer radiative transfer model show large variations for scattering of individual snow layers, but the vertical backscatter distribution can be approximated by an exponential function representing uniform absorption and scattering properties. We obtain estimates of the elevation bias by inverting the interferometric volume correlation coefficient (coherence), applying a uniform volume model for describing the vertical loss function. Whereas the mean values of the computed elevation bias and the elevation difference between the TanDEM-X DEMs and the REMA show good agreement, a trend towards overestimation of penetration is evident for heavily wind-exposed areas with low accumulation and towards underestimation for areas with higher accumulation rates. In both cases deviations from the uniform volume structure are the main reason. In the first case the dense sequence of horizontal structures related to internal wind crust, ice layers and density stratification causes increased scattering in near-surface layers. In the second case the small grain size of the top snow layers causes a downward shift in the scattering phase centre.

Automated Water Level Monitoring at the Continental Scale from ICESat-2 Photons

Of the approximately 6700 lakes and reservoirs larger than 1 km2 in the Contiguous United States (CONUS), only ~430 (~6%) are actively gaged by the United States Geological Survey (USGS) or their partners and are available for download through the National Water Information System database. Remote sensing analysis provides a means to fill in these data gaps in order to glean a better understanding of the spatiotemporal water level changes across the CONUS. This study takes advantage of two-plus years of NASA’s ICESat-2 (IS-2) ATLAS photon data (ATL03 products) in order to derive water level changes for ~6200 overlapping lakes and reservoirs (&gt;1 km2) in the CONUS. Interactive visualizations of large spatial datasets are becoming more commonplace as data volumes for new Earth observing sensors have markedly increased in recent years. We present such a visualization created from an automated cluster computing workflow that utilizes tens of billions of ATLAS photons which derives water level changes for all of the overlapping lakes and reservoirs in the CONUS. Furthermore, users of this interactive website can download segmented and clustered IS-2 ATL03 photons for each individual waterbody so that they may run their own analysis. We examine ~19,000 IS-2 derived water level changes that are spatially and temporally coincident with water level changes from USGS gages and find high agreement with our results as compared to the in situ gage data. The mean squared error (MSE) and the mean absolute error (MAE) between these two products are 1 cm and 6 cm, respectively.

Evaluation of a Statistical Approach for Extracting Shallow Water Bathymetry Signals from ICESat-2 ATL03 Photon Data

In this study we present and validate a simple empirical method to obtain bathymetry profiles using the geolocated photon data from the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) mission, which was launched by NASA in September 2018. The satellite carries the Advanced Topographic Laser Altimeter System (ATLAS), which is a lidar that can detect single photons and calculate their bounce point positions. ATLAS uses a green laser, causing some of the photons to penetrate the air–water interface. Under the right conditions and in shallow waters (&lt;40 m), these photons are reflected back to ATLAS after interaction with the ocean bottom. Using ICESat-2 data from four different overflights above the Heron Reef, Australia, a comparison with SDB data showed a median absolute deviation of approximately 18 cm and Root Mean Square Errors (RMSEs) down to 28 cm. Crossovers between two different overflights above the Heron Reef showed a median absolute difference of 13 cm. For an area north-west of Sisimiut, Greenland, the comparison was done with multibeam echo sounding data, with RMSEs down to between 35 cm, and correspondingly showed median absolute deviations between 33 and 49 cm. The proposed method works well under good conditions with clear waters such as in the Great Barrier Reef; however, for more difficult areas a more advanced machine learning technique should be investigated in order to find an automated method that can distinguish between bathymetry and other signals and noise.

Satellite and UAV-based remote sensing for assessing the flooding risk from Tibetan lake expansion and optimizing the village relocation site

Climate change in recent decades led to the remarkable expansions for most lakes in endorheic basins of the Tibetan Plateau (TP). Enlarged lake inundation areas may pose adverse effects and potential threats on the local human living environment, especially for high-risk villages adjacent to rapidly expanding lakes. Taking a rapidly expanding lake, Angzi Co in the central TP as a study case, we investigated the flooding risk of lake growth on the local living environment and proposed an optimized solution of village relocation selection on the basis of satellite and unmanned aerial vehicle (UAV) remote sensing. The detection of spatiotemporal variations of Angzi Co using optical and altimetric satellite observations revealed a significant area and water level increase by 81.28 km2 and 5.78 m, respectively, from 2000 to 2020. We also assessed the vertical accuracy of multi-source digital elevation model (DEM) products using Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) altimetry data and further examined the flooding risk and potential influences of lake expansion on adjacent settlements (Guozha Village). Results indicated that UAV-DEM achieves excellent advantages in depicting details of lake shoreline variations and simulating potential submergence regions, followed by Advanced Land Observing Satellite World 3D DEM (AW3D DEM). Moreover, assuming that Angzi Co maintains the water level at a growth rate of 0.29 m/a (the average change rate during 2000–2020), the village will be submerged in approximate 10 years based on our assessment. Furthermore, we designed an optimal relocation site southwest of Guozha Village and approximately 3 km away based on the GIS-MVDA method and field investigations. An initial remote sensing-based approach for assessing the flooding risk from dramatic lake expansions in the TP and optimizing the village relocation site was proposed in this study to provide an essential scientific reference for formulating risk mitigation solutions under future climate change scenarios.

Digital terrain model elevation corrections using space-based imagery and ICESat-2 laser altimetry

The most prevalent global surface models are derived from space-based technologies. The spatial and temporal coverage of radar and imaging systems have significant advantage over other sensors and platforms. However, these systems are often are challenged in certain environmental conditions to produce accurate elevations as compared to what might achieved with laser altimetry. Radar derived DEM (Digital Elevation Model) elevation accuracy is often less in vegetated regions, as the wavelengths associated with radar mapping missions do not fully penetrate the canopy. Elevation errors are also substantial over dynamic topography. Correction models using airborne lidar reference datasets can be effective for localized studies to improve the DEM but often the data is not readily available or seasonally/temporally irrelevant. This work presents an automated method for correcting digital terrain models derived from the Shuttle Radar Topographic Mission (SRTM) elevation products using NASA's newest Earth observing laser altimeter, the Ice, Cloud and Land Elevation Satellite-2 (ICESat-2). ICESat-2 elevations, in concert with Landsat 8 (global imagery), create a model based correction strategy for the SRTM derived elevations using geographically correlated canopy cover and surface slope information. The results are validated at the study site using high-resolution, high fidelity airborne lidar datasets as a reference surface. The application of the correction model on the radar measurements results in nearly 50% improvement in elevation accuracy for this region. Additionally this established proof of concept provides a starting point for further research in how this method, with ICESat-2 data, can be extended to other environmental regions, other radar or image derived elevation products and inform future techniques for similar application at the global scale.

ICESat‐2 Atmospheric Channel Description, Data Processing and First Results

AbstractThe Advanced Topographic Laser Altimeter System (ATLAS) was launched aboard the Ice Cloud and land‐Elevation Satellite‐2 (ICESat‐2) satellite in September 2018. ATLAS is a single wavelength (532 nm) lidar system designed to acquire high resolution measurements of the earth's surface while also obtaining atmospheric backscatter from molecules, clouds, and aerosols. Because ATLAS is optimized for altimetry, the atmospheric data acquired is unique in many respects and requires non‐standard analysis techniques. For example, the high repetition rate laser limits the vertical extent of the profiles to just 14 km and causes atmospheric scattering from above 15 km to be added to the scattering in the lower 0–14 km profile. In addition, the limited vertical range of the acquired profiles renders it difficult to compute the magnitude of the solar background and hinders the application of standard calibration techniques. Despite these limitations, methods have been developed to successfully produce data products that have value to the atmospheric community for cloud and aerosol research and are currently available at the National Snow and Ice Data Center (NSIDC). In this paper we describe the ICESat‐2 atmospheric channel and the methods used to process the ATLAS raw photon count data to obtain calibrated backscatter and higher level products such as layer heights and type, blowing snow, column optical depth and apparent surface reflectance.

Quantifying Surface-Height Change Over a Periglacial Environment With ICESat-2 Laser Altimetry.

We use Ice, Cloud, and land Elevation Satellite 2 (ICESat‐2) laser altimetry crossovers and repeat tracks collected over the North Slope of Alaska to estimate ground surface‐height change due to the seasonal freezing and thawing of the active layer. We compare these measurements to a time series of surface deformation from Sentinel‐1 interferometric synthetic aperture radar (InSAR) and demonstrate agreement between these independent observations of surface deformation at broad spatial scales. We observe a relationship between ICESat‐2‐derived surface subsidence/uplift and changes in normalized accumulated degree days, which is consistent with the thermodynamically driven seasonal freezing and thawing of the active layer. Integrating ICESat‐2 crossover estimates of surface‐height change yields an annual time series of surface‐height change that is sensitive to changes in snow cover during spring and thawing of the active layer throughout spring and summer. Furthermore, this time series exhibits temporal correlation with independent reanalysis datasets of temperature and snow cover, as well as an InSAR‐derived time series. ICESat‐2‐derived surface‐height change estimates can be significantly affected by short length‐scale topographic gradients and changes in snow cover and snow depth. We discuss optimal strategies of post‐processing ICESat‐2 data for permafrost applications, as well as the future potential of joint ICESat‐2 and InSAR investigations of permafrost surface‐dynamics.