National Elevation Data Research Articles

Integration of heterogeneous terrain data into Discrete Global Grid Systems

ABSTRACT The Canadian Digital Elevation Model (CDEM) and the High-Resolution Digital Elevation Model (HRDEM) released by Natural Resources Canada are primary terrain data sources in Canada. Due to their different coverage, datums, resolutions, and accuracies, a standardized framework for national elevation data across various scales is required. This study provides new insights into the adoption of Discrete Global Grid Systems (DGGS) to facilitate the integration of multi-source terrain data at various granularities. In particular, the Icosahedral Snyder Equal Area Aperture 3 Hexagonal Grid (ISEA3H) was employed, and quantization, integration, and aggregation were conducted on this framework. To demonstrate the modeling process, an experiment was undertaken for two areas in Ontario, taking advantage of parallel computing which was beneficial from the discreteness of DGGS cells. The accuracy of the modeled elevations was estimated by referring to the ground-surveyed values and was included in the spatially referenced metadata as an indicator of data quality. This research can serve as a guide for future development of a national elevation service, providing consistent, multi-resolution elevations and avoiding complex, duplicated pre-processing at the user’s end. Future investigation into an operational integration platform to support real-world decision-making, as well as the DGGS-powered geospatial datacube, is recommended.

Automated extraction of hydrographically corrected contours for the conterminous United States: the US Geological Survey US Topo product

ABSTRACTThe US Topo is a new generation of digital topographic maps delivered by the US Geological Survey (USGS). These maps include contours in the traditional 7.5-min quadrangle format. The process for producing digital elevation contours has evolved over several years, but automated production of contours for the US Topo product began in 2010. This process, which is quite complex yet fairly elegant, automatically chooses the proper USGS quadrangle, captures the corresponding 1/3 as grid points from the national elevation data set (3D Elevation Program), and adjusts elevation data to better fit water features from the National Hydrography Dataset. After additional processing, such as identifying and tagging depressions, constructing proper contours across double-line streams, and omitting contours from water bodies, contours are automatically produced for the quadrangle. The resulting contours are then compared to the contours on the original (legacy) topographic map sheets, or to the 10-m contours from the original map sheets. Where the elevation surface used to generate the contours has been derived from the previously published contours for a quadrangle, the generated contours match the legacy contours quite well. Where newer elevation sources, such as lidar, originate the elevation surface, generated contours may vary significantly from the previous cartographically produced contours due to more accurate representations of the surface, and less reliance on cartographic interpretation.

Regional rates of young US forest growth estimated from annual Landsat disturbance history and IKONOS stereo imagery

Forests of the Contiguous United States (CONUS) have been found to be a large contributor to the global atmospheric carbon sink. The magnitude and nature of this sink is still uncertain and recent studies have sought to define the dynamics that control its strength and longevity. The Landsat series of satellites has been a vital resource to understand the long-term changes in land cover that can impact ecosystem function and terrestrial carbonstock. We combine annual Landsat forest disturbance history from 1985 to 2011 with single date IKONOS stereoimagery to estimate the change in young forest canopy height and above ground live dry biomass accumulation for selected sites in the CONUS. Our approach follows an approximately linear growth rate following clearing over short intervals and does not estimate the distinct non-linear growth rate over longer intervals.We produced canopy height models by differencing digital surface models estimated from IKONOS stereo pairs with national elevation data (NED). Correlations between height and biomass were established independently using airborne LiDAR, and then applied to the IKONOS-estimated canopy height models. Graphing current biomass against time since disturbance provided biomass accumulation rates. For 20 study sites distributed across five regions of the CONUS, 19 showed statistically significant recovery trends (p is less than 0.001) with canopy growth from 0.26 m yr-1to 0.73 m yr-1. Above ground live dry biomass (AGB) density accumulation ranged from 1.31 t/ha yr-1 to 12.47 t/ha yr-1. Mean forest AGB accumulationwas 6.31 t/ha yr-1 among all sites with significant growth trends. We evaluated the accuracy of our estimates by comparing to field estimated site index curves of growth, airborne LiDAR data, and independent model predictions of C accumulation. Growth estimates found with this approach are consistent with site index curves and total biomass estimates fall within the range of field estimates. This is aviable approach to estimate forest biomass accumulation in regions with clear-cut harvest disturbances.

Extracting Ground-Level DEM From SRTM DEM in Forest Environments Based on Mathematical Morphology

The demands for bare-earth or ground-level digital elevation models (DEMs) are significant for various environmental and ecological studies. As one of the most widely used global DEMs, the Shuttle Radar Topography Mission (SRTM) DEM is available to the public. However, the SRTM DEM is not a bare-earth DEM because it includes man-made structures and vegetation. The objective of this paper is to develop a mathematical morphology-based approach to generate ground-level DEM (GLD) from the SRTM DEM in forest environments. The proposed algorithm is implemented as follows. First, an initial GLD is derived from the SRTM DEM by using morphological operations with a single structuring element. Second, homogeneous forest patches are generated by applying watershed segmentation to pseudocanopy height (PCH) that is obtained by subtracting the initial GLD from the SRTM DEM. Based on segmented patches, a refined GLD is derived by using morphological operations with adaptive structuring elements of different sizes. Third, PCH is updated with the refined GLD and then resegmented into forest patches, from which a final GLD is obtained by subtracting the updated mean PCH from the SRTM DEM. Finally, bare-earth DEMs from the National Elevation Data Set (NED) are used to validate the extracted GLD. The results show that the root mean square error of the final GLD compared with the NED is significantly reduced for two study sites. This type of GLD would be applicable to large-scale environmental studies where accurate topographical information is not available.

Deciphering the Precision of Stereo IKONOS Canopy Height Models for US Forests with G-LiHT Airborne LiDAR

Few studies have evaluated the precision of IKONOS stereo data for measuring forest canopy height. The high cost of airborne light detection and ranging (LiDAR) data collection for large area studies and the present lack of a spaceborne instrument lead to the need to explore other low cost options. The US Government currently has access to a large archive of commercial high-resolution imagery, which could be quite valuable to forest structure studies. At 1 m resolution, we here compared canopy height models (CHMs) and height data derived from Goddard’s airborne LiDAR Hyper-spectral and Thermal Imager (G-LiHT) with three types of IKONOS stereo derived digital surface models (DSMs) that estimate CHMs by subtracting National Elevation Data (NED) digital terrain models (DTMs). We found the following in three different forested regions of the US after excluding heterogeneous and disturbed forest samples: (1) G-LiHT DTMs were highly correlated with NED DTMs with R2 > 0.98 and root mean square errors (RMSEs) < 2.96 m; (2) when using one visually identifiable ground control point (GCP) from NED, G-LiHT DSMs and IKONOS DSMs had R2 > 0.84 and RMSEs of 2.7 to 4.1 m; and (3) one GCP CHMs for two study sites had R2 > 0.7 and RMSEs of 2.6 to 3 m where data were collected less than four years apart. Our results suggest that IKONOS stereo data are a useful LiDAR alternative where high-quality DTMs are available.

Evaluation of data sources for mapping glacial meltwater features

The meltwater system of disintegrating ice sheets provides an important source of information for the reconstruction of ice-retreat patterns during deglaciation. Recent method development in glacial geomorphology, using satellite imagery and digital elevation models (DEMs) for glacial landform mapping, has predominantly been focused on the identification of lineation and other large-scale accumulation features. Landforms created by meltwater have often been neglected in these efforts. Meltwater features such as channels, deltas and fossil shorelines were traditionally mapped using stereo interpretation of aerial photographs. However, during the transition into the digital era, driven by a wish to cover large areas more economically, meltwater features were lost in most mapping surveys. We have evaluated different sets of satellite images and DEMs for their suitability to map glacial meltwater features (lateral meltwater channels, eskers, deltas, ice-dammed lake drainage channels and fossil shorelines) in comparison with the traditional mapping from aerial photographs. Several sets of satellite images and DEMs were employed to map the landform record of three reference areas, located in northwestern Scotland, northeastern Finland and western Sweden. The employed satellite imagery consisted of Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Satellite Pour l'Observation de la Terre (SPOT) 5 and Indian Remote Sensing (IRS) 1C, and the DEMs used were from NEXTMap Britain, Panorama, National elevation data set of Sweden and National Land Survey of Finland. ASTER images yielded better results than the panchromatic band of Landsat 7 ETM+ in all three regions, despite the same spatial resolution of the data. In agreement with previous studies, this study shows that DEMs display accumulation features such as eskers suitably well. Satellite images are shown to be insufficiently detailed for the interpretation of smaller features such as meltwater channels. Hence, satellite imagery and DEMs of intermediate resolution contain meltwater system information only at a general level that allows for the identification of landforms of medium to large sizes. It is therefore pertinent that data with an appropriate spatial and spectral resolution are accessed to fulfil the need of a particular mapping effort. Stereo interpretation of aerial photographs continues to be an advisable method for local meltwater system reconstructions; alternatively, it can be replaced by mapping from high-resolution DEMs such as NEXTMap Britain. For regional to sub-continental reconstructions, the use of ASTER satellite imagery is recommended, because it provides both spectral and spatial resolutions suitable for the identification of meltwater features on a medium to large scale.

Regional distribution of forest height and biomass from multisensor data fusion

Elevation data acquired from radar interferometry at C‐band from SRTM are used in data fusion techniques to estimate regional scale forest height and aboveground live biomass (AGLB) over the state of Maine. Two fusion techniques have been developed to perform post‐processing and parameter estimations from four data sets: 1 arc sec National Elevation Data (NED), SRTM derived elevation (30 m), Landsat Enhanced Thematic Mapper (ETM) bands (30 m), derived vegetation index (VI) and NLCD2001 land cover map. The first fusion algorithm corrects for missing or erroneous NED data using an iterative interpolation approach and produces distribution of scattering phase centers from SRTM‐NED in three dominant forest types of evergreen conifers, deciduous, and mixed stands. The second fusion technique integrates the USDA Forest Service, Forest Inventory and Analysis (FIA) ground‐based plot data to develop an algorithm to transform the scattering phase centers into mean forest height and aboveground biomass. Height estimates over evergreen (R2 = 0.86, P < 0.001; RMSE = 1.1 m) and mixed forests (R2 = 0.93, P < 0.001, RMSE = 0.8 m) produced the best results. Estimates over deciduous forests were less accurate because of the winter acquisition of SRTM data and loss of scattering phase center from tree‐surface interaction. We used two methods to estimate AGLB; algorithms based on direct estimation from the scattering phase center produced higher precision (R2 = 0.79, RMSE = 25 Mg/ha) than those estimated from forest height (R2 = 0.25, RMSE = 66 Mg/ha). We discuss sources of uncertainty and implications of the results in the context of mapping regional and continental scale forest biomass distribution.

Statistical fusion of lidar, InSAR, and optical remote sensing data for forest stand height characterization: A regional‐scale method based on LVIS, SRTM, Landsat ETM+, and ancillary data sets

A method is presented to characterize forest stand heights in a 110,000 km2 region in the eastern United States surrounding the Chesapeake Bay area, driven by a statistical fusion model solely based on remote sensing data. The predicted map was tested against ground survey data from the Forest Inventory and Analysis (FIA) plot network. Input data to the model were 2003 medium footprint lidar data from the Laser Vegetation Imaging Sensor (LVIS) sensor, interferometric radar data from the 2000 Shuttle Radar Topography Mission (SRTM), 1999–2001 Landsat ETM+ data, and ancillary data sets of land cover and canopy density developed for the 2001 National Land Cover Database. In the presented approach, the interferometric synthetic aperture radar (InSAR), optical, and ancillary data sets were masked to the forested areas of the study region and used to segment the raster data stack. The generated image objects closely represented quasi‐homogenous forest stands. For a small region in the study area covered by an LVIS acquisition, LVIS lidar data were then used within the established segments to extract lidar‐based mean forest stand heights. Subsequently these LVIS mean stand heights were used as the response variable to the statistical prediction model (randomForest) which had segment‐based metrics like mean InSAR height (derived from SRTM minus ground digital elevation model data from the National Elevation Data set), mean optical reflectance (derived from Landsat ETM+ Tassled Cap Data), and ancillary metrics as predictive variables. The model developed over the area where LVIS data were available was then applied to map the entire study region. Independent validation of the model was performed in two ways. First, splitting of the model data stack into training and independent testing populations, i.e., testing on LVIS data. This test was deemed to describe the model performance within the LVIS swath. Second, predicted heights were compared to plot height metrics derived from FIA data in the entire study region, thus testing the validity of the model across the larger study area. Results, which are somewhat tampered by the time disconnect between the various data collections, showed the validity and usefulness of this approach. Independent LVIS testing resulted in a correlation coefficient r = 0.83 with an RMSE of 3.0 m (9% error), independent FIA data tested with r = 0.71 with an RMSE of 4.4 m (13% error).

Comparative analysis of SRTM–NED vegetation canopy height to LIDAR‐derived vegetation canopy metrics

Vegetation canopy heights derived from the SRTM 30 m grid DEM minus USGS National Elevation Data (NED) DTM were compared to three vegetation metrics derived from a medium footprint LIDAR data (LVIS) for the US Sierra Nevada forest in California. Generally the SRTM minus NED was found to underestimate the vegetation canopy height. Comparing the SRTM–NED‐derived heights as a function of the canopy percentile height (shape/vertical structure) derived from LVIS, the SRTM SAR signal was found to penetrate, on average, into about 44% of the canopy and 85% after adjustment of the data. On the canopy type analysis, it was found that the SRTM phase scattering centres occurred at 60% for red fir, 53% for Sierra mixed conifer, 50% for ponderosa pine and 50% for montane hardwood‐conifer. Whereas analysing the residual errors of the SRTM–NED minus the LVIS‐derived canopy height as a function of LVIS canopy height and cover it was observed that the residuals generally increase with increasing canopy height and cover. Likewise, the behaviour of the RMSE as a function of canopy height and cover was observed to initially increase with canopy height and cover but saturates at 50 m canopy height and 60% canopy cover. On the other hand, the behaviour of the correlation coefficient as a function of canopy height and cover was found to be high at lower canopy height (<15 m) and cover (<20%) and decrease rapidly making a depression at medium canopy heights (>15 m and <50 m) and cover (>20% and <50%). It then increases with increasing canopy height and cover yielding a plateau at canopies higher than 50 m and cover above 70%.

Two-dimensional, high-resolution modeling of urban dam-break flooding: A case study of Baldwin Hills, California

Modeling of dam-break flooding in an urban residential area in southern California is presented. Modeling is performed using BreZo, an unstructured grid, Godunov-type, finite volume model that solves the shallow-water equations. The model uses terrain data from a 1.5 m Light Detection and Ranging (LiDAR) Digital Terrain Model (DTM) and contour data depicting the reservoir and breach geometry. A spatially distributed Manning coefficient based on a landcover classification derived from digital orthophotos and vector data (e.g., parcel outlines) is also used, and the interception of flow by storm drains is modeled with sink terms in the 2D continuity equation. The model is validated with flood extent and stream flow measurements, and a sensitivity analysis is completed to identify the necessary level of data and model complexity for accuracy purposes. Results show street depressions in the land surface should be resolved by the computational mesh for flood extent and stream flow accuracy. A ca. 5 m resolution mesh that spans streets by approximately 3 cells achieves a good balance between accuracy and computational effort. Results also show that heterogeneous resistance is important for stream flow accuracy, and the interception of overland flow by storm sewers is important for flood extent accuracy. The sensitivity of predictions to several additional factors such as the reservoir level, breach geometry and DTM source (LiDAR, National Elevation Data, Shuttle Radar Topography Mission Data) is also reported.

Evaluation of on-line DEMs for flood inundation modeling

Recent and highly accurate topographic data should be used for flood inundation modeling, but this is not always feasible given time and budget constraints so the utility of several on-line digital elevation models (DEMs) is examined with a set of steady and unsteady test problems. DEMs are used to parameterize a 2D hydrodynamic flood simulation algorithm and predictions are compared with published flood maps and observed flood conditions. DEMs based on airborne light detection and ranging (LiDAR) are preferred because of horizontal resolution, vertical accuracy (∼0.1 m) and the ability to separate bare-earth from built structures and vegetation. DEMs based on airborne interferometric synthetic aperture radar (IfSAR) have good horizontal resolution but gridded elevations reflect built structures and vegetation and therefore further processing may be required to permit flood modeling. IfSAR and shuttle radar topography mission (SRTM) DEMs suffer from radar speckle, or noise, so flood plains may appear with non-physical relief and predicted flood zones may include non-physical pools. DEMs based on national elevation data (NED) are remarkably smooth in comparison to IfSAR and SRTM but using NED, flood predictions overestimate flood extent in comparison to all other DEMs including LiDAR, the most accurate. This study highlights utility in SRTM as a global source of terrain data for flood modeling.

Geographic Information Systems/National Elevation Data Route Mileage Verification

Maintaining the nation’s roads and highways is a critical national infrastructure activity. To do so requires that accurate inventories of roadway features as well as the roadway network itself be maintained. This paper addresses the accuracy of road network mileages, that is, the lengths of various road segments and routes. A study was conducted using Geographic Information Systems to determine the accuracy of using national elevation data (NED) to calculate these lengths. The advantage of doing so is that, if successful, costly and time consuming field studies could be eliminated or reduced. This savings could positively impact all State Departments of Transportation. The scope of the study includes interstate roads as well as all other road types including secondary roads. The primary contribution of the paper is its finding that using NED data provides significant advantages over field studies and provides a very acceptable accuracy level. It also identifies areas of questionable accuracy and indicates why. Additionally, the paper provides some insight into secondary road data which is an often neglected area.

Characterization and Quantification of Data Voids in the Shuttle Radar Topography Mission Data

The goal of this study was to characterize and quantify the occurrence of data voids in data from the Shuttle Radar Topography Mission (SRTM) for the conterminous United States. For this purpose, SRTM data and corresponding data from the national elevation data were downloaded in 21 samples spatially organized to cover the main topography of the U.S. Void locations in SRTM data were compared to terrain attributes and subsequently the area of individual data voids to the same attributes. It was found that data voids amounted to 0.3% of the total dataset. Data voids were found in all topographic settings but more often in slopes steeper than approximately 20/spl deg/ that face south (170/spl deg/), and also in flat areas such as lakes and rivers. It was also found that more than 50% of all data voids were composed of connected pixels in groups less than six pixels. The largest data voids could be attributed to water bodies, while the rest could be explained by terrain-radar interaction characteristics.