Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geomorphology”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hydrobiogeochemical Controls on the Delivery of Dissolved Organic Matter to Boreal Headwater Streams

Understanding controls on the delivery of dissolved organic matter (DOM) from terrestrial to aquatic systems is key for constraining carbon balances and fluxes in boreal headwater catchments. The largest export of DOM generally occurs during intense periods of water delivery to the landscape (e.g., rain events); however, the timing and magnitude of DOM is complicated by geomorphological, hydrometeorological, and biogeochemical variability. To better assess mechanisms controlling the delivery of DOM to boreal streams, DOM was investigated in two morphologically distinct catchments of an experimental forest watershed in western Newfoundland, one dominated by low relief wetlands and the other by steep hillslopes. The DOM was compared during two fall storms of similar magnitude but contrasting moisture conditions during the autumn transition. Variable concentrations and optical character highlighted antecedent conditions are important for the delivery of DOM to boreal streams. Concentration-discharge relationships with stream hydrograph separation analysis suggest preferential flowpaths through shallow mineral horizons as a key pathway for delivery of DOM in the hillslope dominated catchment compared to rapid input from near-stream wetlands in the low relief catchment. Loss of DOM via preferential pathways suggests an additional mechanism for surface soil carbon loss not directly informed by lateral flow and soil hydrology, which has important implications for our conceptualization and modeling of carbon fluxes in boreal systems. These results suggest watershed scale responses of DOM in boreal headwater catchments are complex, and better linkages between catchment scale hydrology and landscape hydrobiogeochemistry are required to constrain terrestrial to aquatic carbon fluxes in boreal landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Reduced Sediment Settling in Turbulent Flows Due To Basset History and Virtual Mass Effects

The behavior of suspended particles in turbulent flows is a recalcitrant problem spanning wide-ranging fields including geomorphology, hydrology, and dispersion of particulate matter in the atmosphere. One key mechanism underlying particle suspension is the difference between particle settling velocity (w s ) in turbulence and its still water counterpart (w so ). This difference is explored here for a range of particle-to-fluid densities (1 – 10) and particle diameter to Kolmogorov micro-eddy sizes (0.1 – 10). Conventional models of particle fluxes that equate w s to w so result in eddy diffusivities and turbulent Schmidt numbers contradictory to laboratory experiments. Incorporating virtual mass and Basset history forces resolves these inconsistencies, providing clarity as to why w s /w so is sub-unity for the aforementioned conditions. The proposed formulation can be imminently used to model particle settling in turbulence, especially when sediment distribution outcomes over extended time scales far surpassing turbulence time scales are sought.

Geology↗
DOE OSTI2024

A Tale of Two Catchments: Causality Analysis and Isotope Systematics Reveal Mountainous Watershed Traits That Regulate the Retention and Release of Nitrogen

Abstract Mountainous watersheds are characterized by variability in functional traits, including vegetation, topography, geology, and geomorphology, which determine nitrogen (N) retention, and release. Coal Creek and East River are two contrasting catchments within the Upper Colorado River Basin that differ markedly in total nitrate (NO 3 − ) export. The East River has a diverse vegetation cover, and sinuous floodplains, and is underlain by N‐rich marine shale. At 0.21 ± 0.14 kg ha −1 yr −1 , the East River exports ∼3.5 times more NO 3 − relative to the conifer‐dominated Coal Creek (0.06 ± 0.02 kg ha −1 yr −1 ). While this can partly be explained by the larger size of the East River, the distinct watershed traits of these two catchments imply different mechanisms controlling the aggregate N‐export signal. A causality analysis shows physical and biogenic processes were critical in determining NO 3 − export from the East River catchment. Stable isotope ratios of NO 3 − (δ 15 N NO3 and δ 18 O NO3 ) show the East River catchment is a strong hotspot for biogeochemical processing of NO 3 − at the hillslope soil‐saprolite. By contrast, the conifer‐dominated Coal Creek retained nearly all atmospherically deposited NO 3 − , and its export was controlled by catchment hydrological traits (i.e., snowmelt periods and water table depth). The conservative N‐cycle within Coal Creek is likely due to the abundance of conifer trees, and smaller riparian regions, retaining more NO 3 − overall and reduced processing prior to export. This study highlights the value of integrating isotope systematics to link watershed functional traits to mechanisms of watershed element retention and release.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Topography Controls Variability in Circumpolar Permafrost Thaw Pond Expansion

Abstract One of the most conspicuous signals of climate change in high‐latitude tundra is the expansion of ice wedge thermokarst pools. These small but abundant water features form rapidly in depressions caused by the melting of ice wedges (i.e., meter‐scale bodies of ice embedded within the top of the permafrost). Pool expansion impacts subsequent thaw rates through a series of complex positive and negative feedbacks which play out over timescales of decades and may accelerate carbon release from the underlying sediments. Although many local observations of ice wedge thermokarst pool expansion have been documented, analyses at continental to pan‐Arctic scales have been rare, hindering efforts to project how strongly this process may impact the global carbon cycle. Here we present one of the most geographically extensive and temporally dense records yet compiled of recent pool expansion, in which changes to pool area from 2008 to 2020 were quantified through satellite‐image analysis at 27 survey areas (measuring 10–35 km 2 each, or 400 km 2 in total) dispersed throughout the circumpolar tundra. The results revealed instances of rapid expansion at 44% (15%) of survey areas. Considered alone, the extent of departures from historical mean air temperatures did not account for between site variation in rates of change to pool area. Pool growth was most clearly associated with upland (i.e., hilly) terrain and elevated silt content at soil depths greater than one meter. These findings suggest that, at short time scales, pedologic and geomorphologic conditions may exert greater control on pool dynamics in the warming Arctic than spatial variability in the rate of air temperature increases.

Abolt, C. J.↗
DOE OSTI2021

Warm-season net CO 2 uptake outweighs cold-season emissions over Alaskan North Slope tundra under current and RCP8.5 climate

Abstract Arctic warming has increased vegetation growth and soil respiration during recent decades. The rate of Arctic warming will likely amplify over the 21st century. Previous studies have revealed that the most severe Arctic warming occurred during the cold season (September to May). The cold-season warming has posited significant CO 2 emissions to the atmosphere via respiration, possibly offsetting warm-season (June to August) net CO 2 uptake. However, prevailing Earth system land models poorly represent cold-season CO 2 emissions, making estimates of Arctic tundra annual CO 2 budgets highly uncertain. Here, we demonstrate that an improved version of the energy exascale Earth system model (E3SM) land model (ELMv1-ECA) captures the large amount of cold-season CO 2 emissions over Alaskan Arctic tundra as reported by two independent, observationally-constrained datasets. We found that the recent seven-decades warming trend of cold-season soil temperature is three times that of the warm-season. The climate sensitivity of warm-season net CO 2 uptake, however, is threefold higher than for the cold-season net CO 2 loss, mainly due to stronger plant resilience than microbial resilience to hydroclimatic extremes. Consequently, the modeled warm-season net CO 2 uptake has a larger positive trend (0.74 ± 0.14 gC m −2 yr −1 ) than that of cold-season CO 2 emissions (0.64 ± 0.11 gC m −2 yr −1 ) from 1950 to 2017, supported by enhanced plant nutrient uptake and increased light- and water-use efficiency. With continued warming and elevated CO 2 concentrations under the representative concentration pathway (RCP) 8.5 scenario, the increasing rate of warm-season net CO 2 uptake is more than twice the rate of cold-season emissions (1.33 ± 0.32 gC m −2 yr −1 vs 0.50 ± 0.12 gC m −2 yr −1 ), making the modeled Alaskan Arctic tundra ecosystem a net CO 2 sink by 2100. However, other geomorphological and ecological disturbances (e.g. abrupt permafrost thaw, thermokarst development, landscape-scale hydrological changes, wildfire, and insects) that are not considered here might alter our conclusion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Surface parameters and bedrock properties covary across a mountainous watershed: Insights from machine learning and geophysics

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

58 GEOSCIENCES↗
DOE OSTI2022

The dynamic atmospheric and aeolian environment of Jezero crater, Mars

Despite the importance of sand and dust to Mars geomorphology, weather, and exploration, the processes that move sand and that raise dust to maintain Mars’ ubiquitous dust haze and to produce dust storms have not been well quantified in situ, with missions lacking either the necessary sensors or a sufficiently active aeolian environment. Perseverance rover’s novel environmental sensors and Jezero crater’s dusty environment remedy this. In Perseverance’s first 216 sols, four convective vortices raised dust locally, while, on average, four passed the rover daily, over 25% of which were significantly dusty (“dust devils”). More rarely, dust lifting by nonvortex wind gusts was produced by daytime convection cells advected over the crater by strong regional daytime upslope winds, which also control aeolian surface features. One such event covered 10 times more area than the largest dust devil, suggesting that dust devils and wind gusts could raise equal amounts of dust under nonstorm conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Coastal typologies and surface and subsurface characteristics of the Alaskan Beaufort Sea Coast

This dataset was generated to classify the Alaskan Beaufort Sea Coast (ABSC) into a set of distinct coastal typologies, to understand the surface and subsurface characteristics and variability of the ABSC, and to quantify relationships between these characteristics and historical rates of shoreline change. This geospatial dataset contains two csv files of points along the ABSC at a 50 m spacing, one for points sheltered by a barrier island and one for points exposed to the open ocean. Each point has a lat/lon location, and we have attributed to each point average values for elevation, historical long-term shoreline change rates, shoreline orientation, landcover, mean annual ground temperature, geomorphic unit, lithology, geology, ecological landscape unit, maximum thaw settlement potential, massive ice content, and segregated ice content. Each point is also assigned to one of 16 coastal typologies, determined by a hierarchical clustering algorithm on the elevation, shoreline change, orientation, and ground temperature data. There are 9 sheltered typologies and 7 exposed typologies, identified by an integer label in the last column of each csv file. The other two csv files contain the integer IDs and classes for the landcover and geomorphology datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

East River Pore Water FTICR-MS Data associated with "Geomorphic and Hydrogeochemical Controls on Fine-Scale Hyporheic Microbiome Assembly and Function"

Pore water samples were collected from distributed locations around Meander A in East River (Crested Butte, CO, USA) during the summer of 2018. This dataset consists of the characterization of dissolved organic matter using 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS; in .xml format) analyzed through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory).This study aimed to understand organic matter assemblages differ within a streambed with differing geomorphology and the interplay between streambed organic matter and microbiome function. Any published work that utilizes data presented in this dataset should cite the dataset using the DOI number. If using the dataset, please acknowledge the U.S. Department of Energy (DOE) Office of Science (BER) Subsurface Biogeochemical Research (SBR) program, the Environmental Molecular Sciences Laboratory (EMSL), and Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS). The publication DOI will be added to this abstract when available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

rabpro: global watershed boundaries, river elevation profiles, and catchment statistics

River and Basin Profiler (rabpro) is a Python package to delineate watersheds, extract river flowlines and elevation profiles, and compute watershed statistics for any location on the Earth’s surface. As fundamental hydrologically-relevant units of surface area, watersheds are areas of land that drain via aboveground pathways to the same location, or outlet. Delineations of watershed boundaries are typically performed on digital elevation models (DEMs) that represent surface elevations as gridded rasters. Depending on the resolution of the DEM and the size of the watershed, delineation may be very computationally expensive. With this in mind, we designed rabpro to provide user-friendly workflows to manage the complexity and computational expense of watershed calculations given an arbitrary coordinate pair. In addition to basic watershed delineation, rabpro will extract the elevation profile for a watershed’s mainchannel flowline. This enables the computation of river slope, which is a critical parameter in many hydrologic and geomorphologic models. Finally, rabpro provides a user-friendly wrapper around Google Earth Engine’s (GEE) Python API to enable cloud-computing of zonal watershed statistics and/or time-varying forcing data from hundreds of available datasets. Altogether, rabpro provides the ability to automate or semi-automate complex watershed analysis workflows across broad spatial extents.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Biogeochemical Processes Across Aquatic Interfaces

The aquatic interfaces exposing terrestrial soils to oxic-anoxic regime shifts represent biogeochemical “hotspots” that are extremely sensitive to climate and environmental change. However, processes and interaction across theses aquatic interfaces are poorly understood and underrepresented in current Earth system models. In this project, we aim to develop predictive understanding of the feedbacks between microbial systems and geochemical environments that determine emergent ecosystem behaviors and resilience in response to disturbances. We use experimental, mechanistic modeling and meta-analysis tools to elucidate interactions among soil, water, geomorphology and microbiology that regulate the molecular transformations and fluxes of carbon, nutrients, and redox-sensitive compounds across aquatic interfaces.

58 GEOSCIENCES↗
DOE OSTI2022

Restoration Hydro: A Watershed Approach to Standard Modular New Hydropower

The objectives of FOA DE- FOA-0001836- “Standard Modular Hydropower” included designing a standardized, modular, and environmentally compatible hydropower schematic for implementation in greenfield sites that generate up to 10 MW of capacity. Utilizing funds competitively awarded under DOE’s Water Power Technologies Office, the Natel Energy team developed a concept for modular new stream reach (NSR) hydropower that incorporates multi species upstream and downstream fish passage, improved river channel connectivity, and recreational modules. The in-stream design of the collective modules minimized site specific design and maximized the opportunities for modularity. Financial data was also presented using actual costs from regional suppliers, with figures provided in 2022 dollars. While the project team did not address potential permitting process improvements, the site selection criteria did consider established barriers to hydropower development such as tribal and preserved lands, interconnection proximity, and endangered species to exclude or deprioritize. The project’s design schematic met the objectives of the FOA, and presented a unique solution that targets alluvial pockets as natural features for sustainable development. Natel’s concept also incorporated the company's fish-safe Restoration Hydro Turbine for safe downstream passage, while featuring a rock arch that integrates fish passage, water, recreation, and grade control modules (including sediment). Alignment with the Department of Energy Office of Energy Efficiency and Renewable Energy (EERE) “Innovative Design Concepts for Standard Modular Hydropower and Pumped-Storage Hydropower” Program: According to the Hydropower Vision (DOE, 2016), approximately 16 GW of hydropower growth is possible with the development of technology solutions that balance efficiency, economics, and environmental sustainability. The desired outcome of the SMH program is transformational innovation specifically in the site identification, conceptual, and detailed design phases of technology development lifecycles (DOE, 2018). In developing the SMH design schematic, the team aimed to address the opportunities outlined in the Vision through an inverted design philosophy; rather than singularly prioritizing efficiency and power production, the team focused on integrating hydropower with restoration of degraded streams to optimal ecosystem function and provision of exceptional recreation value as design criteria. To achieve this, Restoration Hydro incorporates the principles of nature-based engineering (WWAP, 2018) and biomimicry (Biomimicry NL.) to strategically deploy complementary combinations of permanent, semi-permanent, and ephemeral low-head structures - such as natural and engineered log jams - that harness geomorphological and hydrological processes at the landscape-scale. Primary applications of Restoration Hydro include: 1) restoration of degraded watersheds’ natural ecological function and enhancement of hydrological connectivity; and 2) creation of associated co-benefits to hydro production, including increased groundwater recharge, improved sediment transport and management, improved water security and water quality. Restoration Hydro projects build upon proven watershed restoration engineering techniques by integrating hydropower turbines into low-head structures using innovative and evolving civil works concepts that facilitate fish and sediment passage, and in some cases create additional revenue-generating recreational opportunities. Powering low-head structures creates a directly monetizable layer of economic value in the form of flexible, reliable, renewable energy on top of the already high-value water, environmental and recreational benefits of watershed and river restoration. The approach aims to create a virtuous, self-reinforcing cycle whereby Restoration Hydro projects support the scaling of ecosystem restoration activities, creating a water-energy-carbon multiplier effect that, through the principles of adaptive change management: 1) improves the resilience of landscapes and downstream population centers for changing hydrological cycles; 2) creates a reliable energy resource that facilitates the integration of intermittent renewable power sources into grids; and 3) supports climate change mitigation through grid decarbonization and enhanced ecosystem carbon capture and retention.

13 HYDRO ENERGY↗
DOE OSTI2026

Reservoir Sediment Management and Monitoring Database

Overview This dataset compiles dam sediment management and monitoring information from surveys, case studies, and journal articles. Additionally, features described by the National Inventory of Dams (i.e., presence of sluice gates) are included to indicate known infrastructure features that may address sediment releases. The location and description of records from downstream monitoring gages are catalogued in order to help with tracking conditions over time (e.g., before and after management actions, as operations change, etc.). The data help address national scale understanding of challenges and solutions related to the accumulation of sediment behind a dam as well as downstream passage. Sediment trapping causes problems as it reduces storage capacity, disrupts dam and reservoir function, impedes access for recreation, alters water quality/habitat conditions, and contributes to riverbank and coastal erosion within the reservoir. Data compilation from a variety of sources is a first step towards assessing system-wide efficacy of management solutions. This dataset was developed under the Water Power Technologies Office funded effort which began as a Seedling on Reservoir Sedimentation Data, and was supported by the Reservoir Sedimentation Modeling Framework and Data Analysis project. These projects have addressed challenges in describing sediment transport, trapping, and management at dams throughout the US. Methodology An outer join on dams/reservoirs with surveys and survey reports (documented in the RESSED database, USBR or USACE databases, project websites, etc.) with the National Inventory of Dams, based on the NIDID to determine dams with documented management and/or sluice gates. Additional dams with documented management activity were identified through review of technical articles from the past 25 years in Journal of Hydrology, Journal of Water Resources Planning and Management, Geomorphology, Journal of Hydraulic Engineering, Water, Journal of Cleaner Production, International Journal of Sediment Research, Nature Scientific Reports, Earth Surface Processes and Landforms, and Environmental Science and Pollution Research. Individual records were created for each survey or management activity documented. To evaluate downstream sediment monitoring records, the nhdPlusTools and dataRetrieval packages in R were used to find gages within 10km of each dam in the management database. Length of record and location of matched gages were retrieved for those parameters relevant to sediment concentration or total sediment discharge.

Hansen, Carly [ORNL] (ORCID:0000000193280838)↗
DOE OSTI2021

Disastrous Flash Floods Triggered by Moderate to Minor Rainfall Events. Recent Cases in Coastal Benguela (Angola)

The present work focuses on two recent flash floods in coastal Benguela (Angola), both triggered by moderate rainfall but which had disastrous consequences for local populations (namely 71 deaths in 2015 and 17 in 2019). The research involved a regional survey to establish the effects of these floods combined with a geomorphological and socio-economic analysis of the most affected areas to understand the main forcing factors. The two flash floods produced major damage in restricted sectors within very small coastal catchments (<16 km2). The prevalence of fine-grained sedimentary rocks, relatively steep hills, thin soil cover, and vegetation scarcity are natural factors that promote surface runoff. However, socio-economic conditions are most likely the main reasons of flood damage. Namely, rapid population growth with poor planning and making use of low-quality construction materials, the high waste yields that are not properly managed and the absence of flood risk awareness. In the small valleys around the fast-growing cities of coastal Benguela, hazardous flash floods occur recurrently, even after moderate precipitation. Most affected areas are determined by local conditions that compromise drainage at the time of the rainfall event, being very difficult to predict.

Dinis, Pedro A. (ORCID:0000000175587369)↗
DOE OSTI2020

Longitudinal, Lateral, Vertical, and Temporal Thermal Heterogeneity in a Large Impounded River: Implications for Cold-Water Refuges

Dam operations can affect mixing of the water column, thereby influencing thermal heterogeneity spatially and temporally. This occurs by restricting or eliminating connectivity in longitudinal, lateral, vertical, and temporal dimensions. We examined thermal heterogeneity across space and time and identified potential cold-water refuges for salmonids in a large impounded river in inland northwestern USA. To describe these patterns, we used thermal infrared (TIR) imagery, in situ thermographs, and high-resolution, 3-D hydraulic mapping. We explained the median water temperature and probability of occurrence of cool-water areas using generalized additive models (GAMs) at reach and subcatchment scales, and we evaluated potential cold-water refuge occurrence in relation to these patterns. We demonstrated that (1) lateral contributions from tributaries dominated thermal heterogeneity, (2) thermal variability at confluences was approximately an order of magnitude greater than of the main stem, (3) potential cold-water refuges were mostly found at confluences, and (4) the probability of occurrence of cool areas and median water temperature were associated with channel geomorphology and distance from dam. These findings highlight the importance of using multiple approaches to describe thermal heterogeneity in large, impounded rivers and the need to incorporate these types of rivers in the understanding of thermal riverscapes because of their limited representation in the literature.

Mejia, Francine H.↗
DOE OSTI2021

Early Holocene cold snaps and their expression in the moraine record of the eastern European Alps

Glaciers preserve climate variations in their geological and geomorphological records, which makes them prime candidates for climate reconstructions. Investigating the glacier–climate system over the past millennia is particularly relevant first because the amplitude and frequency of natural climate variability during the Holocene provides the climatic context against which modern, human-induced climate change must be assessed. Second, the transition from the last glacial to the current interglacial promises important insights into the climate system during warming, which is of particular interest with respect to ongoing climate change. Evidence of stable ice margin positions that record cooling during the past 12 kyr are preserved in two glaciated valleys of the Silvretta Massif in the eastern European Alps, the Jamtal (JAM) and the Laraintal (LAR). We mapped and dated moraines in these catchments including historical ridges using beryllium-10 surface exposure dating (10Be SED) techniques and correlate resulting moraine formation intervals with climate proxy records to evaluate the spatial and temporal scale of these cold phases. The new geochronologies indicate the formation of moraines during the early Holocene (EH), ca. 11.0 ± 0.7 ka (n = 19). Boulder ages along historical moraines (n = 6) suggest at least two glacier advances during the Little Ice Age (LIA; ca. 1250–1850 CE) around 1300 CE and in the second half of the 18th century. An earlier advance to the same position may have occurred around 500 CE. The Jamtal and Laraintal moraine chronologies provide evidence that millennial-scale EH warming was superimposed by centennial-scale cooling. The timing of EH moraine formation coincides with brief temperature drops identified in local and regional paleoproxy records, most prominently with the Preboreal Oscillation (PBO) and is consistent with moraine deposition in other catchments in the European Alps and in the Arctic region. This consistency points to cooling beyond the local scale and therefore a regional or even hemispheric climate driver. Freshwater input sourced from the Laurentide Ice Sheet (LIS), which changed circulation patterns in the North Atlantic, is a plausible explanation for EH cooling and moraine formation in the Nordic region and in Europe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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

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.

58 GEOSCIENCES↗
DOE OSTI2020

Simulator for Hydrologic Unstructured Domains (SHUD v1.0): numerical modeling of watershed hydrology with the finite volume method

Abstract. Hydrologic modeling is an essential strategy for understanding and predicting natural flows, particularly where observations are lacking in either space or time or where complex terrain leads to a disconnect in the characteristic time and space scales of overland and groundwater flow. However, significant difficulties remain for the development of efficient and extensible modeling systems that operate robustly across complex regions. This paper introduces the Simulator for Hydrologic Unstructured Domains (SHUD), an integrated, multiprocess, multiscale, flexible-time-step model, in which hydrologic processes are fully coupled using the finite volume method. SHUD integrates overland flow, snow accumulation/melt, evapotranspiration, subsurface flow, groundwater flow, and river routing, thus allowing physical processes in general watersheds to be realistically captured. SHUD incorporates one-dimensional unsaturated flow, two-dimensional groundwater flow, and a fully connected river channel network with hillslopes supporting overland flow and baseflow. The paper introduces the design of SHUD, from the conceptual and mathematical description of hydrologic processes in a watershed to the model's computational structures. To demonstrate and validate the model performance, we employ three hydrologic experiments: the V-catchment experiment, Vauclin's experiment, and a model study of the Cache Creek Watershed in northern California. Ongoing applications of the SHUD model include hydrologic analyses of hillslope to regional scales (1 m2 to 106 km2), water resource and stormwater management, and interdisciplinary research for questions in limnology, agriculture, geochemistry, geomorphology, water quality, ecology, climate and land-use change. The strength of SHUD is its flexibility as a scientific and resource evaluation tool where modeling and simulation are required.

58 GEOSCIENCES↗