More Like this

1. Computing water flow through complex landscapes – Part 3: Fill–Spill–Merge: flow routing in depression hierarchies

   https://doi.org/10.5194/esurf-9-105-2021  

   Barnes, Richard; Callaghan, Kerry L.; Wickert, Andrew D. (January 2021, Earth Surface Dynamics)

   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
2. Computing water flow through complex landscapes – Part 1: Incorporating depressions in flow routing using FlowFill

   https://doi.org/10.5194/esurf-7-737-2019  

   Callaghan, Kerry L.; Wickert, Andrew D. (January 2019, Earth Surface Dynamics)

   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
3. The Role of the Critical Zone Structure on the Hydrology and Pool Patterning of Boreal Peatlands

   Comas, X; Slater, L.; Reeve, A. (October 2021, FastTIMES)

   Northern peatlands are unique ecosystems typically located in boreal latitudes that sustain unique habitats and species. They are a critical component of the global carbon cycle, acting both as a carbon reservoir (via peat accumulation after sequestration of carbon dioxide from the atmosphere) and as a carbon source to the atmosphere (by releasing methane and carbon dioxide, two greenhouse gases produced by microbes that thrive in the conditions found in peatlands). Although the ecology of peatlands has been well studied for many decades, the geological controls on peatland development and groundwater flow patterns are not completely understood. In Maine (USA), peatlands began forming about 10,000 years ago following the retreat of the ice sheets at the end of the last ice age. They formed in depressions, often starting as lakes or wetlands, within the landscape carved by glaciers and draped with sediments. Almost two decades of research in several peatlands in Maine suggest that glacial landforms buried beneath peatlands may play a key role in regulating both the hydrology (including the presence of open water pools at the surface) and release of methane gasses from peatlands. Subsurface images using an array of hydrogeophysical methods (including ground-penetrating radar, GPR) constrained with direct coring has revealed the presence of buried esker complexes beneath (or close to) surface pools. Hydrological measurements further suggest that the presence of these permeable esker deposits (mainly gravel and sand) may enhance the connection of peatland water to underlying groundwater and help sustain these pools. Whereas geological maps show that the presence of esker systems in Maine are widespread, satellite-derived digital elevation models (DEMs) reveal that they are often proximal to peatland boundaries, further suggesting that they may commonly extend below the peat formation and control hydrology and carbon dynamics in more peatlands than previously thought. Better understanding of how the critical zone influences coupled water and carbon cycling in northern peatlands may improve the understanding of their contribution to radiative forcing of climate as the climate warms. 

   more » « less
4. MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

   https://doi.org/10.57931/1904379  

   Prusevich, Alexander; Lammers, Richard; Glidden, Stanley (January 2022, MultiSector Dynamics - Living, Intuitive, Value-adding, Environment)

   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
5. Characterization of Flowpath Using Geochemistry and 87Sr/86Sr Isotope Ratios in the Yalahau Region, Yucatan Peninsula, Mexico

   https://doi.org/10.3390/w12092587  

   McKay, Jessica; Lenczewski, Melissa; Leal-Bautista, Rosa Maria (September 2020, Water)

   null (Ed.)

   The Yalahau region, located in the northeastern portion of the Yucatán Peninsula, hosts a series of elongated depressions trending north/south in the direction of Isla Holbox, identified as the Holbox Fracture Zone. Previous studies have explored the geomorphology and various hydrologic characteristics of the Yucatán Peninsula; however, there is a paucity of data concerning the interior region where the fractures are located. Strontium isotope ratios and major ion geochemistry data of the surface water and groundwater of this region serve as a hydrogeochemical fingerprint, aiding in constraining the hydrological boundaries, determining flow paths, and characterizing hydrogeochemical processes that impact the composition of the groundwater within the region. 87Sr/86Sr isotope ratios indicate a different signature than the surrounding bedrock Sr ratio, suggesting that the flow throughout the Yalahau region is moving through channels faster than that of much of the Yucatán. Through major ion geochemistry and 87Sr/86Sr isotope ratios, we were able to delineate at least two flow paths within the Yalahau region and identify a point of saline intrusion at least 35 km from the coast. Gaining an understanding of the hydrogeochemistry and water flow regions is crucial in determining the impact of various activities (e.g., extensive tourism, drinking water withdrawal, wastewater discharge/injection) that occur within the Yucatán Peninsula. 

   more » « less