Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.
more »
« less
A new depression‐dominated delineation ( D ‐cubed) method for improved watershed modelling
Prior to hydrologic modelling, topographic features of a surface are derived, and the surface is divided into sub‐basins. Surface delineation can be described as a procedure, which leads to the quantitative rendition of surface topography. Different approaches have been developed for surface delineation, but most of them may not be applicable to depression‐dominated surfaces. The main objective of this study is to introduce a new depression‐dominated delineation (D‐cubed) method and highlight its unique features by applying it to different topographic surfaces. The D‐cubed method accounts for the hierarchical relationships of depressions and channels by introducing the concept of channel‐based unit (CBU) and its connection with the concept of puddle‐based unit (PBU). This new delineation method implements a set of new algorithms to determine flow directions and accumulations for puddle‐related flats. The D‐cubed method creates a unique cascaded channel‐puddle drainage system based on the channel segmentation algorithm. To demonstrate the capabilities of the D‐cubed method, a small laboratory‐scale surface and 2 natural surfaces in North Dakota were delineated. The results indicated that the new method delineated different surfaces with and without the presence of depressional areas. Stepwise changes in depression storage and ponding area were observed for the 3 selected surfaces. These stepwise changes highlighted the dynamic filling, spilling, and merging processes of depressions, which need to be considered in hydrologic modelling for depression‐dominated areas. Comparisons between the D‐cubed method and other methods emphasized the potential consequences of use of artificial channels through the flats created by the depression‐filling process in the traditional approaches. In contrast, in the D‐cubed method, sub‐basins were further divided into a number of smaller CBUs and PBUs, creating a channel‐puddle drainage network. The testing of the D‐cubed method also demonstrated its applicability to a wide range of digital elevation model resolutions. Consideration of CBUs, PBUs, and their connection provides the opportunity to incorporate the D‐cubed method into different hydrologic models and improve their simulation of topography‐controlled runoff processes, especially for depression‐dominated areas.
more » « less- PAR ID:
- 10039757
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 31
- Issue:
- 19
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 3364-3378
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Calculating flow routing across a landscape is a routine process in geomorphology, hydrology, planetary science, and soil and water conservation. Flow-routing calculations often require a preprocessing step to remove depressions from a DEM to create a “flow-routing surface” that can host a continuous, integrated drainage network. However, real landscapes contain natural depressions that trap water. These are an important part of the hydrologic system and should be represented in flow-routing surfaces. Historically, depressions (or “pits”) in DEMs have been viewed as data errors, but the rapid expansion of high-resolution, high-precision DEM coverage increases the likelihood that depressions are real-world features. To address this long-standing problem of emerging significance, we developed FlowFill, an algorithm that routes a prescribed amount of runoff across the surface in order to flood depressions if enough water is available. This mass-conserving approach typically floods smaller depressions and those in wet areas, integrating drainage across them, while permitting internal drainage and disruptions to hydrologic connectivity. We present results from two sample study areas to which we apply a range of uniform initial runoff depths and report the resulting filled and unfilled depressions, the drainage network structure, and the required compute time. For the reach- to watershed-scale examples that we ran, FlowFill compute times ranged from approximately 1 to 30 min, with compute times per cell of 0.0001 to 0.006 s.more » « less
-
more » « less
Abstract Delineating accurate flowlines using digital elevation models is a critical step for overland flow modeling. However, extracting surface flowlines from high‐resolution digital elevation models (HRDEMs) can be biased, partly due to the absence of information on the locations of anthropogenic drainage structures (ADS) such as bridges and culverts. Without the ADS, the roads may act as “digital dams” that prevent accurate delineation of flowlines. However, it is unclear what variables for terrain‐based hydrologic modeling can be used to mitigate the effect of “digital dams.” This study assessed the impacts of ADS locations, spatial resolution, depression processing methods, and flow direction algorithms on hydrologic connectivity in an agrarian landscape of Nebraska. The assessment was conducted based on the offset distances between modeled drainage crossings and actual ADS on the road. Results suggested that: (a) stream burning in combination with the D8 or D‐Infinity flow direction algorithm is the best option for modeling surface flowlines from HRDEMs in an agrarian landscape; (b) increasing the HRDEM resolution was found significant for facilitating accurate drainage crossing near ADS locations; and (c) D8 and D‐Infinity flow direction algorithms resulted in similar patterns of drainage crossing at ADS locations. This research is expected to result in improved parameter settings for HRDEMs‐based hydrologic modeling.
-
null (Ed.)Abstract. Depressions – inwardly draining regions – are common to many landscapes. When there is sufficient moisture, depressions take the form of lakes and wetlands; otherwise, they may be dry. Hydrological flow models used in geomorphology, hydrology, planetary science, soil and water conservation, and other fields often eliminate depressions through filling or breaching; however, this can produce unrealistic results. Models that retain depressions, on the other hand, are often undesirably expensive to run. In previous work we began to address this by developing a depression hierarchy data structure to capture the full topographic complexity of depressions in a region. Here, we extend this work by presenting the Fill–Spill–Merge algorithm that utilizes our depression hierarchy data structure to rapidly process and distribute runoff. Runoff fills depressions, which then overflow and spill into their neighbors. If both a depression and its neighbor fill, they merge. We provide a detailed explanation of the algorithm and results from two sample study areas. In these case studies, the algorithm runs 90–2600 times faster (with a reduction in compute time of 2000–63 000 times) than the commonly used Jacobi iteration and produces a more accurate output. Complete, well-commented, open-source code with 97 % test coverage is available on GitHub and Zenodo.more » « less
-
more » « less
Abstract Thousands of small wetland depression features (cypress domes) dot the low‐relief karst of Big Cypress National Preserve (BICY) in South Florida, USA. We hypothesized that these wetland depressions are organized in a regular pattern, which is atypical of wetlandscapes elsewhere. Regular patterning implies the existence of coupled feedbacks operating at different spatial scales, with local wetland depression expansion (facilitation via karst dissolution) limited by competition among adjacent depressions for finite water resources (inhibition). We sought to test the hypothesis that wetlands in BICY exhibit regular patterning, and to quantify pattern properties to evaluate competing genesis mechanisms. We tested four predictions about landscape structure and geometry using high‐resolution Light Detection and Ranging elevation data from six 2.25‐km2domains across BICY. Specifically, we predicted (1) feature overdispersion resulting from competition between adjacent basins; (2) truncated wetland area distributions due to growth inhibition feedbacks; (3) periodicity in surface elevation indicating a characteristic pattern wavelength; and (4) elevation bimodality indicating distinct upland and wetland states. All four predictions were strongly supported. Depressions were significantly overdispersed and efficiently fill the landscape, generating hexagonal patterning. Wetland areas followed truncated power law scaling, indicating incremental constraints on basin expansion, in contrast to depression areas elsewhere. Variogram and radial spectrum analyses revealed clear periodicity (~150‐ to 250‐m wavelength) in surface elevations. Finally, surface elevations were consistently bimodal with elevation divergence of 10 to 40 cm. Regular patterning of wetland depressions across BICY is clear, implying long‐term biogeomorphic control on landform structure in this karst landscape.
