Coastal elevation data are essential for a wide variety of applications, such as coastal management, flood modelling, and adaptation planning. Low-lying coastal areas (found below 10 m +Mean Sea Level (MSL)) are at risk of future extreme water levels, subsidence and changing extreme weather patterns. However, current freely available elevation datasets are not sufficiently accurate to model these risks. We present DeltaDTM, a global coastal Digital Terrain Model (DTM) available in the public domain, with a horizontal spatial resolution of 1 arcsecond (∼30 m) and a vertical mean absolute error (MAE) of 0.45 m overall. DeltaDTM corrects CopernicusDEM with spaceborne lidar from the ICESat-2 and GEDI missions. Specifically, we correct the elevation bias in CopernicusDEM, apply filters to remove non-terrain cells, and fill the gaps using interpolation. Notably, our classification approach produces more accurate results than regression methods recently used by others to correct DEMs, that achieve an overall MAE of 0.72 m at best. We conclude that DeltaDTM will be a valuable resource for coastal flood impact modelling and other applications.@en

Digital Elevation Models (DEMs) are a necessity for modelling many large-scale environmental processes. In this study, we investigate the potential of data from two spaceborne lidar altimetry missions, ICESat-2 and GEDI—with respect to their vertical accuracies and planimetric data collection patterns—as sources for rasterisation towards creating global DEMs. We validate the terrain measurements of both missions against airborne lidar datasets over three areas in the Netherlands, Switzerland, and New Zealand and differentiate them using land-cover classes. For our experiments, we use five years of ICESat-2 ATL03 data and four years of GEDI L2A data for a total of 252 million measurements. The datasets are filtered using parameter flags provided by the higher-level products ICESat-2 ATL08 and GEDI L3A. For all areas and land-cover classes combined, ICESat-2 achieves a bias of −0.11 m, an MAE of 0.43 m, and an RMSE of 0.93 m. From our experiments, we find that GEDI is less accurate, with a bias of 0.09 m, an MAE of 0.98 m, and an RMSE of 2.96 m. Measurements in open land-cover classes, such as “Cropland” and “Grassland”, result in the best accuracy for both missions. We also find that the slope of the terrain has a major influence on vertical accuracy, more so for GEDI than ICESat-2 because of its larger horizontal geolocation error. In contrast, we find little effect of either beam power or background solar radiation, nor do we find noticeable seasonal effects on accuracy. Furthermore, we investigate the spatial coverage of ICESat-2 and GEDI by deriving a DEM at different horizontal resolutions and latitudes. GEDI has higher spatial coverage than ICESat-2 at lower latitudes due to its beam pattern and lower inclination angle, and a derived DEM can achieve a resolution of 500 m. ICESat-2 only reaches a DEM resolution of 700 m at the equator, but it increases to almost 200 m at higher latitudes. When combined, a 500 m resolution lidar-based DEM can be achieved globally. Our results indicate that both ICESat-2 and GEDI enable accurate terrain measurements anywhere in the world. Especially in data-poor areas—such as the tropics—this has potential for new applications and insights.@en

River ice is a major contributor to flood risk in cold regions due to the physical impediment of flow caused by ice jamming. Although a variety of classifiers have been developed to distinguish ice types using HH or VV intensity of SAR data, mostly based on data from RADARSAT-1 and -2, these classifiers still experience problems with breakup classification, because meltwater development causes overlap in co-polarization backscatter intensities of open water and sheet ice pixels. In this study, we develop a Random Forest classifier based on multiple features of Sentinel-1 data for three main classes generally present during breakup: rubble ice, sheet ice and open water, in a case study over the Athabasca River in Canada. For each ice stage, intensity of the VV and VH backscatter, pseudo-polarimetric decomposition parameters and Grey Level Co-occurrence Matrix texture features were computed for 70 verified sample areas. Several classifiers were developed, based on i) solely intensity features or on ii) a combination of intensity, pseudo-polarimetric and texture features and each classifier was evaluated based on Recursive Feature Elimination with Cross-Validation and pair-wise correlation of the studied features. Results show improved classifier performance when including GLCM mean of VV intensity, and VH intensity features instead of the conventional classifier based solely on intensity. This highlights the complementary nature of texture and intensity for the classification of breaking river ice. GLCM mean incorporates spatial patterns of the co-polarized intensity and sensitivity to context, while VH intensity introduces cross-polarized surface and volume scattering signals and is less sensitive to wind than the commonly used co-polarized intensity. We conclude that the proposed method based on the combination of texture and intensity features is suitable for and performs well in physically complex situations such as breakup, which are hard to classify otherwise. This method has a high potential for classifying river ice operationally, also for data from other SAR missions. Since it is a generic approach, it also has potential to classify river ice along other rivers globally.@en