Impact of an Extreme Flood Event on Streambank Retreat: Cedar River,Nebraska, USA

The 2010 dam breach and consequent anomalous flood event on the Cedar River in Nebraska, USA provided an opportunity to study the following objectives: (1) evaluate the impact of an extreme flood event on streambank retreat along a 45 km stretch relative to the average annual retreat; (2) quantify the changes in streambank retreat for each km segment downstream of the breach; and (3) examine the influence of riparian vegetation and radius of curvature on meander bank erosion rate. During the hydrologic event, discharge peaked at nearly three times greater than the next highest recorded rate and equated to a return period of 2,000 years. Aerial images and ArcGIS were utilized to calculate the average annual streambank retreat for each year during the preflood (2006–2010), flood (2010), and postflood (2010–2016) periods. The 2010 flood period had a significantly higher average annual streambank retreat of 2,820 m²/km/yr than the preflood and postflood periods, which, respectively, measured 576 and 384 m²/km/yr. From 2006 to 2016, 29% of all streambank erosion was from this one extreme flood event, thus demonstrating the impact that one extreme flood event can have on streambank retreat and the geomorphology of a stream system.     

Spatial Relationships of Levees and Wetland Systems within Floodplains of the Wabash Basin,USA

Given the unique biogeochemical, physical, and hydrologic services provided by floodplain wetlands, proper management of river systems should include an understanding of how floodplain modifications influence wetland ecosystems. The construction of levees can reduce river–floodplain connectivity, yet it is unclear how levees affect wetlands within floodplains, let alone the cumulative impacts within an entire watershed. This paper explores spatial relationships between levee and floodplain wetland systems in the Wabash Basin, United States. We used a hydrogeomorphic floodplain delineation technique to map floodplain extents and identify wetlands that may be hydrologically connected to river networks. We then spatially examined the relationship between levee presence, wetland area, and other river network attributes within discrete subbasins. Our results show that cumulative wetland area is relatively constant in subbasins that contain levees, regardless of maximum stream order within the subbasin. In subbasins that do not contain levees, cumulative wetland area increases with maximum stream order. However, we found that wetland distributions around levees can be complex, and further studies on the influence of levees on wetland habitat may need to consider finer resolution spatial scales.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

CHANNEL MORPHOLOGY EVOLUTION AND OVERBANK FLOW IN THE GEORGIA PIEDMONT1

ABSTRACT: Historic changes in stream channel morphology were investigated in the Georgia Piedmont to better understand the hydrologic processes and functioning of the region's riverine systems. USGS gaging station data and channel geomorphology data were collected from thirty study sites in the Upper Oconee River Basin for flood frequency analysis. Historic and modern (i.e., present-day) channel capacity discharge (i.e., overbank flow) was calculated using Manning's equation and historic channel cross-section records. The recurrence interval for overbank flow was estimated for each site from flood frequency data. Results indicate that channel expansion has occurred throughout the basin, especially in upper reaches. Recurrence intervals for modern overbank events were variable and generally high ranging from < 2 to > 500 years for first to third order streams. They were less variable and lower for fourth and fifth order streams, ranging from < 2 to 3 years. Potential depositional thresholds were identified that exemplify the complex response of sediment distribution patterns throughout the basin. Results indicate overbank flows occur less frequently now than they once did due to historic accelerated sedimentation and subsequent channel expansion. One application of these findings is that these basin processes are likely applicable across the region and may impact the hydrologic functioning of associated Piedmont riverine wetlands that depend on flooding regimes.     

Capturing LiDAR‐Derived Hydrologic Spatial Parameters to Evaluate Playa Wetlands

The digital elevation model data from traditional stereo photogrammetric methods are inadequate in providing accurate vertical parameters to feed hydrologic models for low‐lying, extremely flat areas. High‐resolution light detection and ranging (LiDAR) data provide the robust capability of capturing small variations in low‐relief playa wetlands. The Rainwater Basin in south‐central Nebraska includes a complex of seasonally shallow playa wetlands that attract millions of migratory waterfowl every spring and fall. This research focuses on the development of a procedure with applicable protocols to produce LiDAR‐derived three‐dimensional wetland maps and to extract the critical surface parameters (i.e., watershed boundaries, flow direction, flow accumulation, and drainage lines) for playa wetlands. The topo‐hydrologic conditions of playa wetlands were evaluated at the watershed level. The results show that in the Rainwater Basin, 70.7% of the historic hydric soil footprints identified in the Soil Survey Geographic (SSURGO) database were not functioning as topographically depressional wetlands. This finding was confirmed by a recent five‐year Annual Habit Survey showing that 69.8% of the historic hydric soil footprints did not function during the spring migratory bird seasons between 2004 and 2009. The majority of playa wetlands' topographic conditions have been substantially changed and the SSURGO data cannot fully reflect current topographic reality in the Rainwater Basin.     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

Identification of Putative Geographically Isolated Wetlands of the Conterminous United States

Geographically isolated wetlands (GIWs) are wetlands completely surrounded by uplands. While common throughout the United States (U.S.), there have heretofore been no nationally available, spatially explicit estimates of GIW extent, complicating efforts to understand the myriad biogeochemical, hydrological, and habitat functions of GIWs and hampering conservation and management efforts at local, state, and national scales. We used a 10‐m geospatial buffer as a proxy for hydrological or ecological connectivity of National Wetlands Inventory palustrine and lacustrine wetland systems to nationally mapped and available stream, river, and lake data. We identified over 8.3 million putative GIWs across the conterminous U.S., encompassing nearly 6.5 million hectares of wetland resources (average size 0.79 ± 4.81 ha, median size 0.19 ha). Putative GIWs thus represent approximately 16% of the freshwater wetlands of the conterminous U.S. The water regime for the majority of the putative GIWs was temporarily or seasonally flooded, suggesting a vulnerability to ditching or hydrologic abstraction, sedimentation, or alterations in precipitation patterns. Additional analytical applications of this readily available geospatially explicit mapping product (e.g., hydrological modeling, amphibian metapopulation, or landscape ecological analyses) will improve our understanding of the abundance and extent, effect, connectivity, and relative importance of GIWs to other aquatic systems of the conterminous U.S.     

Changing the Course of Rivers in an Asian City: Linking Landscapes to Human Benefits through Iterative Modeling and Design

Concerns over water scarcity, climate change, and environmental health risks have prompted some Asian cities to invest in river rehabilitation, but deciding on the end goals of rehabilitation is a complex undertaking. We propose a multidisciplinary framework linking riparian landscape change to human well‐being, providing information relevant to decision makers, in a format that facilitates stakeholder involvement. We illustrate this through a case study of the densely settled, environmentally degraded, and flood prone Ciliwung River flowing through metropolitan Jakarta, Indonesia. Our methodology attempts to respond to this complexity through an iterative approach, strongly based on conceptualization and mathematical modeling. Nested hydrologic, hydrodynamic, and water quality models provide outputs at catchment‐, corridor‐, and localized site‐scales. Advanced 3‐D landscape modeling is used for procedural design and precise visualization of proposed changes and their impacts, as predicted by the mathematical models. Finally, participatory planning and design methods allow us to obtain critical stakeholder feedback in shaping a socially acceptable approach. Our framework aims at demonstrating that a change in paradigm in river rehabilitation is possible, and providing future scenarios that balance concerns over flooding, water quality, and ecology, with the realities of a rapidly growing megacity.     

Evaluating the Applicability of a Two‐dimensional Flow Model of a Highly Heterogeneous Domain to Flow and Environmental Management

Two‐dimensional simulation of highly heterogeneous domains, especially those with disparate length scales, roughness conditions, and geometries, often leads to challenges such as long computation times and numerical instability. Simulation of challenging domains is often needed to guide flood management and environmental regulation agencies in operation and potential domain modifications. This work evaluates the ability of a two‐dimensional unsteady hydrodynamic model to represent long‐duration transient flows over a domain with highly heterogeneous roughness, geometric characteristics, and length scales through bed roughness representation. The domain includes 13 km of Cache Creek and the 14.5 km² Cache Creek Settling Basin, which traps both sediment and mercury. Calibration under different bed roughness methods, validation, and modeling results of bathymetric modification scenarios are presented. The modeling approach's performance supports its application as a tool for management of similar domains, such as settling basins, leveed floodplains, and reservoirs. Accurate representation of flow dynamics can also inform environmental management that involves transport of sediments, nutrients, and heavy metals. This study found that a two‐dimensional unsteady flow model can accurately represent long‐duration transient flow in a large settling basin with highly heterogeneous characteristics without parsing of the domain or flow events simulated.     

Dynamic floodplain vegetation model development for the Kootenai River, USA

The Kootenai River floodplain in Idaho, USA, is nearly disconnected from its main channel due to levee construction and the operation of Libby Dam since 1972. The decreases in flood frequency and magnitude combined with the river modification have changed the physical processes and the dynamics of floodplain vegetation. This research describes the concept, methodologies and simulated results of the rule-based dynamic floodplain vegetation model "CASiMiR-vegetation" that is used to simulate the effect of hydrological alteration on vegetation dynamics. The vegetation dynamics are simulated based on existing theory but adapted to observed field data on the Kootenai River. The model simulates the changing vegetation patterns on an annual basis from an initial condition based on spatially distributed physical parameters such as shear stress, flood duration and height-over-base flow level. The model was calibrated and the robustness of the model was analyzed. The hydrodynamic (HD) models were used to simulate relevant physical processes representing historic, pre-dam, and post-dam conditions from different representative hydrographs. The general concept of the vegetation model is that a vegetation community will be recycled if the magnitude of a relevant physical parameter is greater than the threshold value for specific vegetation; otherwise, succession will take place toward maturation stage. The overall accuracy and agreement Kappa between simulated and field observed maps were low considering individual vegetation types in both calibration and validation areas. Overall accuracy (42% and 58%) and agreement between maps (0.18 and 0.27) increased notably when individual vegetation types were merged into vegetation phases in both calibration and validation areas, respectively. The area balance approach was used to analyze the proportion of area occupied by different vegetation phases in the simulated and observed map. The result showed the impact of the river modification and hydrological alteration on the floodplain vegetation. The spatially distributed vegetation model developed in this study is a step forward in modeling riparian vegetation succession and can be used for operational loss assessment, and river and floodplain restoration projects.     

Hydrological,Physical, and Chemical Functions and Connectivity of Non‐Floodplain Wetlands to Downstream Waters: A Review

We reviewed the scientific literature on non‐floodplain wetlands (NFWs), freshwater wetlands typically located distal to riparian and floodplain systems, to determine hydrological, physical, and chemical functioning and stream and river network connectivity. We assayed the literature for source, sink, lag, and transformation functions, as well as factors affecting connectivity. We determined NFWs are important landscape components, hydrologically, physically, and chemically affecting downstream aquatic systems. NFWs are hydrologic and chemical sources for other waters, hydrologically connecting across long distances and contributing compounds such as methylated mercury and dissolved organic matter. NFWs reduced flood peaks and maintained baseflows in stream and river networks through hydrologic lag and sink functions, and sequestered or assimilated substantial nutrient inputs through chemical sink and transformative functions. Landscape‐scale connectivity of NFWs affects water and material fluxes to downstream river networks, substantially modifying the characteristics and function of downstream waters. Many factors determine the effects of NFW hydrological, physical, and chemical functions on downstream systems, and additional research quantifying these factors and impacts is warranted. We conclude NFWs are hydrologically, chemically, and physically interconnected with stream and river networks though this connectivity varies in frequency, duration, magnitude, and timing.     

Streamflow,Sediment Transport,and Geomorphic Change during the 2011 Flood on the Missouri River Near Bismarck–Mandan,ND

Geomorphic change from extreme events in large managed rivers has implications for river management. A steady‐state, quasi‐three‐dimensional hydrodynamic model was applied to a 29‐km reach of the Missouri River using 2011 flood data. Model results for an extreme flow (500‐year recurrence interval [RI]) and an elevated managed flow (75‐year RI) were used to assess sediment mobility through examination of the spatial distribution of boundary or bed shear stress (τ_b) and longitudinal patterns of average τ_b, velocity, and kurtosis of τ_b. Kurtosis of τ_b was used as an indicator of planform channel complexity and can be applied to other river systems. From differences in longitudinal patterns of sediment mobility for the two flows we can infer: (1) under extreme flow, the channel behaves as a single‐thread channel controlled primarily by flow, which enhances the meander pattern; (2) under elevated managed flows, the channel behaves as multithread channel controlled by the interaction of flow with bed and channel topography, resulting in a more complex channel; and (3) for both flows, the model reach lacks a consistent pattern of deposition or erosion, which indicates migration of areas of erosion and deposition within the reach. Despite caveats and limitations, the analysis provides useful information about geomorphic change under extreme flow and potential implications for river management. Although a 500‐year RI is rare, extreme hydrologic events such as this are predicted to increase in frequency.     

Nitrogen Subsidies from Hillslope Alder Stands to Streamside Wetlands and Headwater Streams,Kenai Peninsula,Alaska

We examined nitrogen transport and wetland primary production along hydrologic flow paths that link nitrogen‐fixing alder (Alnus spp.) stands to downslope wetlands and streams in the Kenai Lowlands, Alaska. We expected that nitrate concentrations in surface water and groundwater would be higher on flow paths below alder. We further expected that nitrate concentrations would be higher in surface water and groundwater at the base of short flow paths with alder and that streamside wetlands at the base of alder‐near flow paths would be less nitrogen limited than wetlands at the base of long flow paths with alder. Our results showed that groundwater nitrate‐N concentrations were significantly higher at alder‐near sites than at no‐alder sites, but did not differ significantly between alder‐far sites and no‐alder sites or between alder‐far sites and alder‐near sites. A survey of ¹⁵N stable isotope signatures in soils and foliage in alder‐near and no‐alder flow paths indicated the alder‐derived nitrogen evident in soils below alder is quickly integrated downslope. Additionally, there was a significant difference in the relative increase in plant biomass after nitrogen fertilization, with the greatest increase occurring in the no‐alder sites. This study demonstrates that streamside wetlands and streams are connected to the surrounding landscapes through hydrologic flow paths, and flow paths with alder stands are potential “hot spots” for nitrogen subsidies at the hillslope scale.     

Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

Does Wetland Location Matter When Managing Wetlands for Watershed‐Scale Flood and Drought Resilience?

Wetland protection and restoration strategies that are designed to promote hydrologic resilience do not incorporate the location of wetlands relative to the main stream network. This is primarily attributed to the lack of knowledge on the effects of wetland location on wetland hydrologic function (e.g., flood and drought mitigation). Here, we combined a watershed‐scale, surface–subsurface, fully distributed, physically based hydrologic model with historical, existing, and lost (drained) wetland maps in the Nose Creek watershed in the Prairie Pothole Region of North America to (1) estimate the hydrologic functions of lost wetlands and (2) estimate the hydrologic functions of wetlands located at different distances from the main stream network. Modeling results showed wetland loss altered streamflow, decreasing baseflow and increasing stream peakflow during the period of the precipitation events that led to major flooding in the watershed and downstream cities. In addition, we found that wetlands closer to the main stream network played a disproportionately important role in attenuating peakflow, while wetland location was not important for regulating baseflow. The findings of this study provide information for watershed managers that can help to prioritize wetland restoration efforts for flood or drought risk mitigation.     

MANAGEMENT OF FLOOD CONTROL SUMPS AND POLLUTANT TRANSPORT1

ABSTRACT: Levee sump systems are used by many riverine communities for temporary storage of urban wet weather flows. The hydrologic performance and transport of stormwater pollutants in sump systems, however, have not been systematically studied. The objective of this paper is to present a case study to demonstrate development and application of a procedure for assessing the hydraulic performance of flood control sumps in an urban watershed. Two sumps of highly variable physical and hydraulic characteristics were selected for analysis. A hydrologic modeling package was used to estimate the flow hydrograph for each outfall as part of the flow balance for the sump. To validate these results, a water balance was used to estimate the total runoff using sump operational data. The hydrologic model calculations provide a satisfactory estimate of the total runoff and its time‐distribution to the sump. The model was then used to estimate pollutant loads to the sump and to the river. Although flow of stormwater through a sump system is regulated solely by flood‐control requirements, these sumps may function as sedimentation basins that provide purification of stormwater. A sample calculation of removals of several conventional pollutants in the target sumps using a mass balance approach is presented.     

Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

Evaluation of SWAT Impoundment Modeling Methods in Water and Sediment Simulations

Worldwide studies show 80%–90% of all sediments eroded from watersheds is trapped within river networks such as reservoirs, ponds, and wetlands. To represent the impact of impoundments on sediment routing in watershed modeling, Soil and Water Assessment Tool (SWAT) developers recommend to model reservoirs, ponds, and wetlands using impoundment tools (ITs). This study evaluates performance of SWAT ITs in the modeling of a small, agricultural watershed dominated by lakes and wetlands. The study demonstrates how to incorporate impoundments into the SWAT model, and discusses and evaluates involved parameters. The study then recommends an appropriate calibration sequence, i.e., landscape parameters calibration, followed by pond/wetlands calibration, then channel parameter calibrations, and lastly, reservoir parameter calibration. Results of this study demonstrate not following SWAT recommendation regarding modeling water land use as an impoundment depreciates SWAT performance, and may lead to misplaced calibration efforts and model over‐calibration. Further, the chosen method to model impoundments’ outflow significantly impacts sediment loads in the watershed, while streamflow simulation is not very sensitive. This study also allowed calculation of mass accumulation rates in modeled impoundments where the annual mass accumulation rate in wetlands (2.3 T/ha/yr) was 39% higher than mass accumulation rate in reservoirs (1.4 T/ha/yr).     

Clogging of an Effluent Dominated Semiarid River: A Conceptual Model of Stream‐Aquifer Interactions1

Abstract: Water managers in arid and semiarid regions increasingly view treated wastewater (effluent) as an important water resource. Artificial recharge basins allow effluent to seep into the ground relieving stressed aquifers, however these basins frequently clog due to physical, chemical, and biological processes. Likewise effluent is increasingly used to maintain perennial base flow for dry streambeds, however, little is known about the impact of effluent on streambed hydraulic conductivity and stream‐aquifer interactions. We address this issue by investigating: if a clogging layer forms, how the formation of a clogging layer alters stream‐aquifer connections, and what hydrologic factors control the formation and removal of clogging layers. We focused on the Upper Santa Cruz River, Arizona where effluent from the Nogales International Waste Water Treatment Plant sustains perennial flow. Monthly sampling, along a 30 km river reach, was done with two foci: physical streambed transformations and water source identification using chemical composition. Historical dataset were included to provide a larger context for the work. Results show that localized clogging occurs in the Upper Santa Cruz River. The clogging layers perch the stream and shallow streambed causing desaturation below the streambed. With these results, a conceptual model of clogging is established in the context of a semiarid hydrologic cycle: formation during the hot premonsoon months when flow is nearly constant and removal by large flood flows (>10 m³/s) during the monsoon season. However, if the intensity of flooding during the semiarid hydrologic cycle is lessened, the dependent riparian area can experience a die off. This conceptual model leads us to the conclusion that effluent dominated riparian systems are inherently unstable due to the clogging process. Further understanding of this process could lead to improved ecosystem restoration and management.