Topography is a fundamental input to hydrologic models critical
for generating realistic streamflow networks as well as infiltration and
groundwater flow. Although there exist several national topographic datasets
for the United States, they may not be compatible with gridded models that
require hydrologically consistent digital elevation models (DEMs). Here, we
present a national topographic dataset developed to support gridded
hydrologic simulations at 1 km and 250 m spatial resolution over the contiguous
United States. The workflow is described step by step in two parts: (a) DEM
processing using a Priority Flood algorithm to ensure hydrologically
consistent drainage networks and (b) slope calculation and smoothing to
improve drainage performance. The accuracy of the derived stream network is
evaluated by comparing the derived drainage area to drainage areas reported
by the national stream gage network. The slope smoothing steps are evaluated
using the runoff simulations with an integrated hydrologic model. Our DEM
product started from the National Water Model DEM to ensure our final
datasets will be as consistent as possible with this existing national
framework. Our analysis shows that the additional processing we provide
improves the consistency of simulated drainage areas and the runoff
simulations that simulate gridded overland flow (as opposed to a network
routing scheme). The workflow uses an open-source R package, and all output
datasets and processing scripts are available and fully documented. All of
the output datasets and scripts for processing are published through CyVerse
at 250 m and 1 km resolution. The DOI link for the dataset is
Topography is one of the most important inputs to hydrologic simulations; it defines watershed boundaries and shapes river networks. Surface flow travel times and runoff characteristics are very sensitive to hillslope characteristics (D'Odorico and Rigon, 2003; Freer et al., 2002; Frei et al., 2010; Gupta and Mesa, 1988). In addition to shaping surface flow networks, groundwater fluxes and residence times are also strongly driven by topographic gradients (e.g., Condon and Maxwell, 2015). High-resolution elevation data are not difficult to find; for example the National Elevation Dataset provides a 30 m resolution digital elevation model (DEM) across the US (Gesch et al., 2002). It is well established that topographic datasets require processing to be suitable for hydrologic simulation because flow networks and slopes can be sensitive to noise in the DEM and can be affected by the resolution and spatial gridding of a DEM (Habtezion et al., 2016; Sørensen and Seibert, 2007; Thompson et al., 2001; Vaze et al., 2010; Wolock and McCabe, 2000; Wu et al., 2008; Zhang and Montgomery, 1994). This processing generally consists of some combination of steps to (1) remove erroneous local minima that impeded flow (Kenny et al., 2008; Lindsay, 2016a, b); (2) lower or “burn” in the drainage network to ensure that flow occurs along identified stream segments (Lindsay, 2016b; Woodrow et al., 2016), and (3) smooth the DEM to remove noise introduced by the sampling resolution (Gallant, 2011; Lindsay et al., 2019).
Although there are many tools and methods available to process topography for hydrologic applications, processing is generally site specific, and we lack a national topographic dataset for the US that is designed specifically for gridded physical hydrology simulations. At the national scale there are four main topographic datasets for the US. First, as previously noted the National Elevation Dataset (NED), provides 30 m national DEM primarily derived from the U.S. Geological Survey (USGS) 10 and 30 m DEMs (Gesch et al., 2002). This provides a high-resolution national topographic dataset that includes hydrologic information such as flow directions and accumulation; however, it is not directly processed for DEM-based hydrologic simulations which may require more smoothing.
Second, the National Hydrography Dataset Plus (NHDPlus) dataset provides a collection of geospatial data, including various key features of stream network including elevation, drainage area and watershed boundary, in multiple resolutions, i.e., NHDPlus High Resolution in 10 m and NHD Medium Resolution in 30 m. It is derived from the 10 m USGS elevation data and the national complete Watershed Boundary Dataset (Viger et al., 2016). Because of its public availability, high resolution and large spatial coverage, the dataset has been applied in many hydrologic studies considering runoff, river routing and flood inundation (David et al., 2011; Garousi-Nejad et al., 2019; Smith et al., 2013; Zhang et al., 2018). The NHDPlus dataset provides the most complete hydrologic mapping of the US; however, spatial inconsistences and noise in these datasets make it difficult to use these data directly for gridded hydrologic modeling. The NHDPlus stream network is derived from various topographic map sources, leading to spatial inconsistencies and inaccuracies in river network (Moore et al., 2019). Spatial discrepancies have also been found between NHDPlus and the local higher-resolution light detection and ranging data (lidar)-derived stream network (Samu, 2012). It was found that the NHDPlus stream network is not well aligned with DEM-used height above nearest drainage (HAND), which can lead to discrepancies in hydraulic properties for NHDPlus reaches (Garousi-Nejad et al., 2019).
Third, the National Water Model (NWM) is a national hydrologic modeling framework that simulates terrestrial hydrology at 250 m spatial resolution. This model has its own DEM and stream mask that are derived from the 30 m DEM in NED and stream network in NHDPlus (Gochis et al., 2018). The NWM DEM was developed using the TauDEM topographic processing tool (Tarboton, 2005) to generate a connected drainage network (Gochis et al., 2018). Although the NWM DEM does include hydrologically based topographic processing, the NWM configuration that this dataset was assembled for uses a network routing approach for streamflow (as opposed to gridded streamflow simulations directly applied to the DEM) (Gochis et al., 2018; Johnson et al., 2019). As a result, stream gradients are calculated on a reach-by-reach basis using an abstracted stream network rather than resulting directly from the gridded DEM. For hydrologic models that rely on partial differential equations (PDEs) to simulate gridded streamflow, additional smoothing and DEM processing is required to improve simulation performance. These models require a topographic dataset which represents these stream reaches accurately within the DEM itself. Also, in PDE-based simulations fluxes occur across grid cell faces, so improved performance can be achieved by directly aligning the DEM processing with grid cell faces.
Finally, the Hydrological data and maps based on SHuttle Elevation
Derivatives at multiple Scales (HydroSHEDS) provide hydrographic
information for global-scale applications. It offers a set of hydrological
datasets, including river networks, watershed boundaries, drainage
directions and flow accumulations at
Here we present a topographic dataset of the contiguous United States (CONUS) designed for PDE-based hydrologic modeling applications that simulate surface and subsurface flows with a single, consistent topographic input. The workflow is discussed step by step in the following sections in two parts: (a) DEM processing by Priority Flood algorithm and (b) slope calculation and smoothing. Several cases are implemented for both parts with comparative analysis to illustrate the improvements of additional procedures. We compared the drainage area with measurements and applied slope into a hydrologic model to evaluate the surface flow simulation.
This section covers the detailed information about input datasets including DEM and stream network maps (Sect. 2.1), topographic processing methodology for DEM adjustment and slope calculation (Sect. 2.2), and evaluation metrics used to assess the accuracy of the drainage network and drainage performance (Sect. 2.3).
Our approach requires three primary inputs: (1) a gridded DEM of topography, (2) a mask of known stream locations, and (3) a mask of lakes and terminal points for endorheic basins (shown in Fig. 1). All analysis covers CONUS plus areas in Canada and Mexico that drain to CONUS. This spatial coverage is consistent with the NOAA National Water Model (NWM) (Sampson and Gochis, 2018). All analysis is completed and available at 1 km and 250 m spatial resolution, but for the purposes of discussion we focus on the 1 km analysis throughout the results section. Our input datasets are primarily derived from the NWM version 1.2. The NWM inputs are derived from publicly available national datasets and have already been processed for hydrologic consistency within their modeling framework, which provides an ideal starting point for this work.
Note that we intentionally start from NWM datasets as opposed to a raw DEM product in order to ensure that the resulting products are as consistent as possible with this growing framework of national hydrologic datasets. Our processing only adjusts the NWM processing where it is needed to improve performance of gridded simulations. This approach is intended to minimize DEM-based discrepancies between modeling approaches and to facilitate future model comparisons and coupling between different modeling approaches.
We start from a 250 m DEM that was developed for the NWM V1.2, which was developed from the 30 m resolution National Elevation Dataset (Gesch et al., 2002) and processed using the TauDEM (Tarboton, 2005) to ensure drainage in a D8 routing scheme (i.e., with flow allowed out of every grid cell in eight directions). The NED dataset provides elevation of bare ground without considering the tall vegetation. We start from this DEM directly for 250 m analysis, and for the 1 km analysis we upscale the NWM DEM by selecting the minimum elevation within every 1 km cell (Condon and Maxwell, 2019). The derived stream network is significantly affected by the upscaling approach (Moretti and Orlandini, 2018). The DEM is upscaled based on the minimum as opposed to the average because this is appropriate in identifying the lowest points and capturing hydrologic features such as streams.
The stream network mask is derived from a rasterized version of the NWM stream vector subset based on Strahler stream order thresholds (Horton, 1945; Strahler, 1957) (note that here too we start from the NWM network to ensure consistency between our final product and existing frameworks). The stream mask is used to guide the drainage patterns in the DEM for the topographic processing in this study. By experimenting with different stream densities, we found that the stream networks consisting of fifth-order and higher streams provide the best guidance in 1 km resolution processing and a denser network consisting of third-order and above streams is ideal for 250 m resolution processing. Figure 1 shows the fifth-order stream network that was used for the 1 km DEM processing. For both the third- and fifth-order stream masks additional manual adjustments were made to improve the resulting drainage network, which is discussed in more detail in Sect. 3.1.
The final input is a mask of lakes and terminal points for endorheic basins. Our topographic processing ensures that all grid cells in the domain drain to either an ocean cell, internal lake or a terminal point within endorheic basins, referred to as “sinks”. The lakes and sinks presented in our domain mask provide end points for all internally draining (i.e., not draining to the ocean) cells. We selected 13 major lakes from the NWM lake shapefile, which was developed based on NHDPlus. Only the major terminal lakes which are critical for watershed delineation were included in this topographic processing. The extents of these lakes are rasterized as shown in Fig. 1 with the borders treated as outlets, and the interior lake cells are excluded from topographic processing.
For the more arid endorheic basins, such as are commonly found in the Great Basin of the western US, intermittent rivers drain to much smaller water bodies which may not be perennially inundated. For these basins, we designate single cells at the river outlet as sink cells which were added to the domain manually based on the fifth-order stream mask. A total of 262 sinks were added to the domain as shown in Fig. 1. The significance of these sinks for topographic processing is discussed further in Sect. 3.1.
Map of the DEM dataset, stream network, lakes and sinks used for processing the 1 km hydrographic dataset.
While some hydrologic processing has already been implemented in NWM DEM as
described in the previous section, additional processing is needed for the
purpose of our processing (1) to ensure a fully connected drainage network
for fluxes occurring across cell faces and (2) to improve runoff performance
for gridded overland flow simulations. In our processing, we require
drainage be accomplished in D4 (i.e., only in the primary north, south, east,
west directions); this is useful for partial differential equation (PDE)-based hydrologic simulations where fluxes occur across cell faces (although
it should be noted that each grid cell has both an
Figure 2 summarizes the topographic processing
workflow for the final dataset we present. All of the processing steps were
completed using the PriorityFlow V1.0 toolbox
(Condon and Maxwell, 2019), which is available
on GitHub ( We first derive a new stream network to be used for the stream elevation
smoothing, using a drainage area threshold (stream_th). The flow direction files
resulting from the Priority Flood processing can be used to calculate the
drainage area for every grid cell in the domain. From this, a stream network
of any density can be derived based on the stream_th. This allows us to apply
smoothing to a denser stream network than was used in the earlier Priority
Flood processing script. The smaller the stream_th is, the denser the identified
stream network will be, and the more cells will be smoothed. This stream
network can then be divided into segments and associated sub-watersheds.
Here we use a drainage area threshold of 100 km Stream segments are smoothed by adjusting the elevations along the segment
to maintain a constant slope from the starting to the ending point of each
segment. After these adjustments are made to the stream reaches an
additional elevation processing pass must be made to ensure that all of the
bank cells along a stream reach still drain to the stream. Here we apply an
epsilon threshold (bank_epsilon) to require that adjacent cells are at least
bank_epsilon greater than their downstream neighbor. The bank adjustment step traverses
up a stream segment checking that all cells draining to every stream cell
(based on flow direction) are at least bank_epsilon (0.1 m in this dataset) higher.
Anywhere this is not the case, neighboring cells are filled accordingly, and
processing continues up the bank until it reaches a point where all upstream
neighboring cells are fully draining.
Workflow diagram outlining the major processing steps used and the outputs generated; input parameters are listed in italics and described in more detail in Table 1.
Parameters for the topographic processing and values in this dataset.
As noted in Sect. 2.1, there are already multiple
national elevation datasets that are well established and have been
previously validated. Our goal is to start from these datasets to develop a
national topographic dataset that is as consistent with these products as
possible but is designed for gridded overland flow simulation. Because we
are starting from established elevation models (DEM for NWM V1.2), our
evaluation is focused on the hydrologic improvements to the dataset.
Specifically, we evaluate (1) the extent to which the resulting drainage
network matches observations and (2) the drainage characteristics of the
resulting dataset. Details on the approach for evaluating each of these
metrics follow.
The model is configured with a single shallow layer (0.1 m) with an
impervious surface (permeability in 10
We explored a wide range of parameter combinations in our testing to obtain a reliable topographic dataset for hydrologic modeling. Table 1 outlines all of the parameters that were used for the final topographic processing reflected in the datasets published here. For the purposes of this discussion, we focus on the processing steps that had the largest impact on the accuracy of the resulting drainage network spatially (Sect. 3.1) and the runoff dynamics (Sect. 3.2). The goal here is not to present a sensitivity study of parameters in topographic processing. Rather this discussion is provided to illustrate the key topographic processing steps that can improve hydrologic performance. This section demonstrates the improvements of the published dataset over existing topographic products, specifically for hydrologic modeling applications. Furthermore, the workflow described above and the results presented here can help guide others developing DEMs for different domains or resolutions.
In this section we explore the impact of incorporating a prior stream
network to guide the topographic processing. Figure 3 shows the difference in drainage area between the processed DEM and
reported stream gage area when the processing is run without an initial
stream mask included in the PriorityFlow adjustment step.
Figure 3a shows the percent difference in drainage
area for all 7542 gages used for comparison, while
Fig. 3b shows only 1195 gages with drainage area
over 5000 km
Drainage area comparison between processed with no initial stream
network and USGS gage measurements.
To address this, we use the fifth-order stream mask to guide the topographic processing of the 1 km domain (third order for the 250 m domain). Figure 4 shows the stream gage performance maps by adding fifth-order stream mask. Adding the stream mask to processing increases the number of gages with area matches within 10 % from 5066 in the original processing without any stream mask provided to 5189 (880 to 937 for the large stream gages). Furthermore, it decreases the number of gages with area differences greater than 20 % to 204 overall and 158 among the large gages. This is a clear improvement over the processing without a stream mask. However, there are still many locations with poor drainage area matches, including large gages.
Drainage area comparison between processed with initial fifth-order stream network and USGS gage measurements.
Therefore, as a final step we applied a series of manual fixes to the stream network. Three primary types of fixes were applied: (1) adding sinks to the domain to increase the number of terminal points for internally draining basins, (2) manually modifying the stream mask to reflect higher-resolution information, (3) removing erroneous stream outlets that were created where the headwater cells in the stream mask intersects the boundary of the domain. This last adjustment corrected some of the anomalously high drainage areas in the norther portion of the domain (blue box in Fig. 4). For discussion purposes we will focus on the first two adjustments that impacted a much larger portion of the stream network.
Sinks were added primarily in the Great Basin portion of the domain (the red
box in Fig. 4), where endorheic basins are common.
This arid portion of the domain is characterized by ephemeral streams (i.e.,
streams that do not flow year-round) that are often poorly mapped.
Figure 5 shows a portion of the Great Basin before
and after sinks were added. Figure 5a shows the
input topography and cells with drainage area over 10
Panels
By far the most common manual update that was made in the processing step was to manually update the location of the stream mask. This was accomplished by visually comparing the fifth-order stream network to the higher-order streams and NHD streamlines and extending the fifth-order stream mask as needed to reflect this additional information in headwater regions. Note that simply applying a denser stream mask nationally (for example using the third-order stream mask for the 1 km domain processing) is not a viable solution here. When translated to the 1 km grid these denser stream masks result in too many grid cells being classified as stream cells and a decreased quality of the stream network overall. Therefore, we work primarily from the fifth-order stream mask but incorporate higher-resolution information as needed in an iterative process comparing our resulting drainage network to the stream gage areas. Figure 6 gives an example of a location at the northern portion of the domain where the fifth-order stream mask was expanded to improve performance of stream network derivation. Figure 6a shows the fifth-order stream network in black lines and third-order stream network in red lines. Figure 6b shows the linear drainage network that results from the initial processing with the fifth-order stream. In this flat region, it is difficult to obtain the correct drainage path with only the topographic information. By adding the third-order stream network, the drainage path presents a more reasonable pattern and an improved drainage area matching in Fig. 6c. A similar process was completed for the 250 m resolution starting from the third-order streams and extending using the second-order streams.
Panels
Figure 7 shows the final drainage area evaluation
after all the adjustments to the stream network are incorporated. Overall
this demonstrates a strong improvement in the drainage area network relative
to both the case with no stream network at all and the original fifth-order stream mask. Nearly all (90.0 %) of the large drainage area gages
(from 880 with no stream mask to 1075 gages with the final stream mask) now
have area agreements within 10 % and 74.2 % of the gages overall (from
5066 with no stream mask to 5595 gages with the final stream mask).
Figure 7c summarizes the performance across the
three test cases grouped by drainage area. As can be seen here, the final
channel network improves performance across all drainage area categories.
Performance remains most variable in the smallest drainage basins (
The 250 m domain is not plotted here but is included in the processed
datasets. Similar analysis of the 250 m DEM shows that, as would be expected,
a higher-resolution DEM provides better information for drainage network
derivation. There are 6048 gages that have area agreements within 10 % using
the unmodified third-order stream mask directly from NWM. For reference
this is 982 more gages than the 1 km DEM with the unmodified stream mask.
Similar improvements have been found by applying all the processes to 250 m
resolution. Here too the processing does not significantly improve
performance for gages less than 100 km
Drainage area comparison between processed with final-order stream
network and USGS gage measurements.
In addition to the drainage network location we evaluate how the topographic
processing influences the runoff characteristics of the domain. As described
in Sect. 2.3, this behavior is evaluated using
runoff tests and assessing anomalously high ponding depths. Note that in
this section we evaluate how our processing reduces ponding locations. This
is not meant to imply that we are trying to get rid of all ponding in the
domain. Rather we are looking for locations where the DEM resolution is
leading to what we expect to be non-physically realistic ponding –
specifically, discontinuities along the drainage network or anomalous
locations along the hillslope. Here we consider the impact of smoothing
along the stream (step 4 in the topographic processing), the flat fix step
applied to the rest of the domain (step 5c in the topographic processing)
and removing the secondary slope along the stream cells (step 5b in the
topographic processing). For reference we compare four cases listed below
with progressively more processing steps applied.
Note that in all of these cases processing steps 1–3 in the topographic
processing remain the same; therefore, the flow directions and location of
the drainage network remain unchanged. The smoothing steps evaluated here do
not alter the primary flow directions; they simply adjust the gradients
along flow paths.
The ponding depth from runoff tests from four slope cases.
Ponding depths from a small domain located in Colorado at hour 20 of the
runoff test are displayed in Fig. 8 to illustrate
the local impact of smoothing. The domain is 16 254 km
As would be expected, the resulting stream network shape is the same, and the largest ponding depths occur along the main stem of the drainage network in all cases. What is notable in this figure though are the spatial differences in how this ponding occurs. In the no smoothing case there are discontinuities along the main stem of the stream as well as and many localized high ponding points across the domain (Fig. 8a). Applying the stream smoothing in the three subsequent cases smooths these discontinuities along the main drainage network where smoothing was applied. An example of this improvements is shown in the blue box of the domain in Fig. 8a and b and the magnified figure in the bottom-right part. Adding the flat fix step decreases the number of isolated ponding points outside the stream network (an example is shown in the red boxes in Fig. 8b and c). Note that even with this smoothing step applied there are still isolated ponding locations. This is acceptable as all cells still have the ability to drain even if it is slow and because our goal is not to achieve uniformly fast drainage everywhere. Rather the flat fix step is intended only to address those locations where large anomalies in the processed DEM result in significant ponding which may not be physically realistic. Finally, the “removed secondary slope along stream” step does not change the smoothing properties of the domain; it simply removes the secondary slope in the stream to force water in the stream drain only in the primary direction, which increases the drainage speed of stream cells.
To further illustrate the impact of the stream smoothing step, Figure 9a plots the elevation along the main stem shown in Fig. 9c before and after smoothing is applied. Segments are plotted in different colors in Fig. 9c and separated by dashed lines in Fig. 9a. As can be seen here, in both cases the stream path is monotonically decreasing. However before stream smoothing is applied, there are many large steps along the stream path. With smoothing constant slopes are applied along each stream segment by adjusting the DEM along the stream. Note that the same elevation smoothing along the stream is carried out in the three smoothed test cases. Figure 9b compares the ponding depth along the stream for all four cases. Comparing the no smooth to the add stream smooth cases, it can be seen that the stream smoothing step has the largest impact on the ponding depths in the stream. Figure 9 illustrates that it significantly reduces the ponding depth variability from cell to cell. Adding the flat fix increases the ponding depth especially shown at the downstream. This is because the flat fix step increases the drainage speed of the domain as a whole, resulting in larger ponding depth in the stream. Finally, removing the secondary slope along the stream results in the smoothest pressure variation along the stream network by draining all water only to the primary direction.
The small domain shown in Figs. 8 and 9 illustrates the major behaviors seen across the domain in response to the smoothing steps. Figure 10 summarizes the results across the CONUS simulation. The background of Fig. 10 shows the ponding depth from the runoff test with the final slope at hour 20, and the barplots summarize the percentage decrease in cells with ponding greater than 0.1 m outside the main stream network after 20 h of the runoff simulation relative to the baseline, no smooth case. The stream cells are excluded from this analysis as larger ponding depth is expected in the stream network. Results are summarized by major watershed outlined in black on the figure. Still, smoothing the stream increases the drainage rate across the domain so the add stream smooth case does significantly decrease the number of ponding points outside the domain. The add flat fix case impacts a small number of cells in most basins, which is also to be expected given that it is targeting isolated discontinuities in the DEM. The flat fix step has the largest relative impact in steeper domains such as the upper Colorado due to the large variability of slopes between neighbor cells here. Removing the secondary slope has the largest impact in flatter portions of the domain such as the Mississippi Embayment.
The background map is the ponding depth from runoff test at hour 20 with the final slope; barplots are the percentage decrease of anomalies' ponding cells (ponding depth over 0.1 m) outside the main stream network relative to the baseline (no smooth case) in HUC2s. The scale of the barplot in each HUC is the same, with the grey bar shown in the right of figure.
Demonstration of output datasets over CONUS:
Summary of gridded outputs included in the dataset.
Figure 11 shows the final elevation map, slope and
drainage area that result from the topographic processing. These maps plot
the 1 km outputs, but all results are also available at 250 m resolution.
Elevations and slopes are the primary inputs for the hydrologic modeling,
but additional datasets defining the drainage networks and watersheds are
also generated through the processing step (as shown by the final outputs in
Fig. 2). While there are other national datasets
that derive drainage networks and watersheds, we include these outputs here
too because they are on the same grid and derived through the same
processing steps as the DEM and therefore provide spatial information which
is perfectly aligned with the final processed topography. The dataset
includes a gridded map of elevation, slope in the
Figure 12 summarizes the differences between the final DEM of this study and the original DEM, including both the fraction of grid cells where the elevation changes and how much elevations were adjusted by. DEM adjustments are split into two categories according to the processing step where they occurred: (1) the Priority Flood adjustment step where elevations were adjusted to ensure complete drainage and (2) the subsequent smoothing steps. Results are summarized by regional hydrologic unit codes (HUC2s). The contiguous US is divided into hydrologic units identified by a unique HUC. The first level of classification divides the contiguous US into 18 major geographic regions, known as HUC2s. Overall the Priority Flood algorithm changed the elevation in 11.6 % of the grid cells, while the stream smoothing steps impacted 17.4 % of the cells. In both cases the median slope adjustment magnitude is small (1.74 m for Priority Flood and 3.12 m for smoothing). The Priority Flood algorithm frequently results in larger elevation adjustments than the stream smoothing step. Note also in Fig. 12b that the Priority Flood algorithm only increases elevations from the original to the processed DEMs, whereas the stream smoothing step can both increase and decrease elevations. Spatially the largest fraction of grid cells is adjusted in the flatter basins, for example the Mississippi Embayment (HUC 8), while the largest elevation adjustments occur in the steeper western portions of the domain such as Colorado and California (HUCs 13–18).
All of the topographic processing was completed with the PriorityFlow R
library, which is available on GitHub (
This study presents a topographic dataset developed for physically based hydrologic simulations. We combine a Priority Flood topographic processing algorithm with a priori stream network mask and USGS stream gage drainage areas to develop DEM and flow direction rasters which will closely match observed drainage networks. We also apply a series of slope smoothing steps along stream reaches and globally to improve surface runoff performance. These smoothing steps are unique from other global topographic smoothing approaches because we directly consider flow direction and stream locations to provide different smoothing along the stream reaches. The resulting DEM is designed to capture hydrologic features and improve runoff simulations for PDE-based hydrologic models which rely on gridded topography for overland flow simulations. This performance is evaluated nationally using reported drainage areas from the USGS stream gage network and runoff simulations. All outputs from the processing are available across the contiguous US at both 1 km and 250 m resolution. The processing workflow developed here uses the open-source R package PriorityFlow and is fully documented with the published datasets so that others can modify processing for different resolutions or specific subdomains as desired.
JZ and LEC were in charge of processing the dataset including processing the DEM and conducting the runoff simulations. JZ prepared the datasets, wrote the manuscript, plotted the figures and did the analysis. LEC provided the PriorityFlow R package, documented the dataset, and helped with the analysis and figure plotting. HT helped in evaluating the runoff simulations and reviewed the manuscript and documentation. RMM provided overall supervision in processing the datasets and reviewed the manuscript.
The authors declare that they have no conflict of interest.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work was supported by the U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research, and Biological and Environmental Sciences IDEAS project, the Sustainable Systems Scientific Focus Area under award number DE-AC02-05CH11231; and the US National Science Foundation Office of Advanced Cyberinfrastructure under award number OAC-1835855.
This research has been supported by the U.S. Department of Energy Office of Science (grant no. DE-AC02-05CH11231) and the US National Science Foundation Office of Advanced Cyberinfrastructure (grant no. OAC-1835855).
This paper was edited by Ge Peng and reviewed by four anonymous referees.