How processing digital elevation models can affect simulated water budgets

For regional models, the shallow water table surface is often used as a source/sink boundary condition, as model grid scale precludes simulation of the water table aquifer. This approach is appropriate when the water table surface is relatively stationary. Since water table surface maps are not readily available, the elevation of the water table used in model cells is estimated via a two-step process. First, a regression equation is developed using existing land and water table elevations from wells in the area. This equation is then used to predict the water table surface for each model cell using land surface elevation available from digital elevation models (DEM). Two methods of processing DEM for estimating the land surface for each cell are commonly used (value nearest the cell centroid or mean value in the cell). This article demonstrates how these two methods of DEM processing can affect the simulated water budget. For the example presented, approximately 20% more total flow through the aquifer system is simulated if the centroid value rather than the mean value is used. This is due to the one-third greater average ground water gradients associated with the centroid value than the mean value. The results will vary depending on the particular model area topography and cell size. The use of the mean DEM value in each model cell will result in a more conservative water budget and is more appropriate because the model cell water table value should be representative of the entire cell area, not the centroid of the model cell.     

Vadose zone-attenuated artificial recharge for input to a ground water model

Accurate representation of artificial recharge is requisite to calibration of a ground water model of an unconfined aquifer for a semiarid or arid site with a vadose zone that imparts significant attenuation of liquid transmission and substantial anthropogenic liquid discharges. Under such circumstances, artificial recharge occurs in response to liquid disposal to the vadose zone in areas that are small relative to the ground water model domain. Natural recharge, in contrast, is spatially variable and occurs over the entire upper boundary of a typical unconfined ground water model. An improved technique for partitioning artificial recharge from simulated total recharge for inclusion in a ground water model is presented. The improved technique is applied using data from the semiarid Hanford Site. From 1944 until the late 1980s, when Hanford's mission was the production of nuclear materials, the quantities of liquid discharged from production facilities to the ground vastly exceeded natural recharge. Nearly all hydraulic head data available for use in calibrating a ground water model at this site were collected during this period or later, when the aquifer was under the diminishing influence of the massive water disposals. The vadose zone is typically 80 to 90 m thick at the Central Plateau where most production facilities were located at this semiarid site, and its attenuation of liquid transmission to the aquifer can be significant. The new technique is shown to improve the representation of artificial recharge and thereby contribute to improvement in the calibration of a site-wide ground water model.     

Delineating ground water recharge from leaking irrigation canals using water chemistry and isotopes

Across the Great Plains irrigation canals are used to transport water to cropland. Many of these canals are unlined, and leakage from them has been the focus of an ongoing legal, economic, and philosophical debate as to whether this lost water should be considered waste or be viewed as a beneficial and reasonable use since it contributes to regional ground water recharge. While historically there has been much speculation about the impact of canal leakage on local ground water, actual data are scarce. This study was launched to investigate the impact of leakage from the Interstate Canal, in the western panhandle of Nebraska, on the hydrology and water quality of the local aquifer using water chemistry and environmental isotopes. Numerous monitoring wells were installed in and around a small wetland area adjacent to the canal, and ground water levels were monitored from June 1992 until January 1995. Using the water level data, the seepage loss from the canal was estimated. In addition, the canal, the monitoring wells, and several nearby stock and irrigation wells were sampled for inorganic and environmental isotope analysis to assess water quality changes, and to determine the extent of recharge resulting from canal leakage. The results of water level monitoring within study wells indicates a rise in local ground water levels occurs seasonally as a result of leakage during periods when the canal is filled. This rise redirects local ground water flow and provides water to nearby wetland ecosystems during the summer months. Chemical and isotopic results were used to delineate canal, surface, and ground water and indicate that leaking canal water recharges both the surface alluvial aquifer and upper portions of the underlying Brule Aquifer. The results of this study indicate that lining the Interstate Canal could lower ground water levels adjacent to the canal, and could adversely impact the local aquifer.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

Heat flow and ground water movement in the Bohemian cretaceous basin (Czechoslovakia)

Sub-surface temperature fields may be considerably affected by active ground water systems, thereby seriously hampering the interpretation of heat flow data. Quantitative evaluation of the convective component of heat transfer is thus very important in cases such as large sedimentary basins with vast underground water circulation. We propose in this study a simple model of horizontal aquifer. This model was used to examine the effect of the lateral convection on the surface heat flow near the recharge zone of basinal margins. The perturbation of the heat flow field above the aquifer was calculated for various aquifer geometry and various flow velocities and the regional scale dependence of the perturbation on the hydraulic properties of the aquifer was demonstrated. The model was applied to the Bohemian Cretaceous Basin and it was shown that within a few kilometres from the recharge zones the observed surface heat flow may be underestimated by up to several tens of percent. The procedure was tested in two locations in this area, in an attempt to make hydrogeological corrections to the measured heat flow values in several boreholes.     

Prediction of the effect of water supplying from Shirindare dam on the Bojnourd aquifer using MODFLOW2000

In recent years, the water demand has been increasing considerably in Bojnourd, capital of Northern Khorasan province in NE of Iran, and the extracted water from Bojnourd alluvial aquifer, with an area of 65 km², is not sufficient for residents. The required water is going to be supplied from Shirindare dam, located out of Bojnourd aquifer’s catchment area; therefore, the groundwater levels will rise in some parts of the aquifer, due to the return flow of supplied water, which will cause serious problems for the city. In this paper, the groundwater flow system of Bojnourd aquifer has been numerically simulated using MODFLOW code in GMS interface. The model, primarily, was calibrated for a steady state condition for the mean values of oneyear period (Sep. 2009 to Sep. 2010) which has a steady condition with low stresses on the aquifer. Then the model was run/calibrated for transient conditions for a two year period (Sep. 2007 to Sep. 2009). After determining the hydraulic properties of the aquifer and confirming their validity, different management scenarios, were applied to the model. Results reveal that groundwater levels in the urban area will rise by over 3 m, by infiltrating 40% of supplied water from the Shirindare dam into the aquifer. To manage the rising water levels, two different management scenarios were applied to the aquifer model. In doing so and with proper management of aquifer exploitation during critical situations, not only will the groundwater level drop; also the city of Bojnourd can develop urban landscaping by constructing sports/cultural camping area using the extra pumped water.     

Modelling coastal ground‐ and surface‐water interactions using an integrated approach

Beach water table fluctuations have an impact on the transport of beach sediments and the exchange of solute and mass between coastal aquifer and nearby water bodies. Details are given of the refinement of a dynamically integrated ground‐ and surface‐water model, and its application to study ground‐ and surface‐water interactions in coastal regions. The depth‐integrated shallow‐water equations are used to represent the surface‐water flow, and the extended Darcy's equation is used to represent the groundwater flow, with a hydrostatic pressure distribution being assumed to apply for both these two types of flows. At the intertidal region, the model has two layers, with the surface‐water layer being located on the top of the groundwater layer. The governing equations for these two types of flows are discretized in a similar manner and they are combined to give one set of linear algebraic equations that can be solved efficiently. The model is used to predict water level distributions across sloping beaches, where the water table in the aquifer may or may not decouple from the free water surface. Five cases are used to test the model for simulating beach water table fluctuations induced by tides, with the model predictions being compared with existing analytical solutions and laboratory and field data published in the literature. The numerical model results show that the integrated model is capable of simulating the combined ground‐ and surface‐water flows in coastal areas. Detailed analysis is undertaken to investigate the capability of the model. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of Nitrate-Nitrogen Movement Near High-Capacity Irrigation Wells

AGalerkin finite-element model coupled with a particle tracking routine was developed to analyze the flow and transport dynamics near a high-capacity irrigation well. The model was used to compute the head distribution around the pumping well, to determine the area of influence, and to define ground water flowlines during short-term pumping periods typical of those used to collect water quality samples from high-capacity wells. In addition to hypothetical example results, the model was used to qualitatively analyze data obtained from pump-and-sample experiments conducted in an unconfined alluvial aquifer within the Platte River valley of south-central Nebraska where nitrate-nitrogen (NO₃-N) contamination is prevalent.
Simulation results of both the hypothetical and field cases suggest that short-term pumping events, impact a limited volume of aquifer. The area of influence and flowlines are affected by aquifer anisotropy, pumping rate, and well construction characteristics). Ground water above or below the screened intervals does not enter a partially penetrating well in anisotropic aquifers. In aquifers where NO₃-N concentration varies vertically and horizontally, waler quality samples from an irrigation, or other high-capacity, well provide only limited information about ground water contamination. A numerical model is thus recommended for calculating the area of influence and determining flowlines around high-capacity wells so that information derived from water quality samples collected at the wellhead can be better interpreted.     

Influence of aquifer properties on phytoremediation effectiveness

Recent research has shown that planting deep-rooted trees, such as poplar, can take up and degrade important ground water pollutants such as trichloroethylene (TCE) as they transpire water from the capillary fringe of shallow contaminated aquifers. The effect of hydrogeologic factors on the minimum plantation area needed to prevent downgradient migration of contaminated ground water is not well known. Accordingly, the objective of this research was to identify the hydrogeologic parameters that control phytoremediation effectiveness. We used a numerical ground water flow model to evaluate the effect that natural variations in hydrogeologic parameters and growing season duration have on the minimum plantation area required for capture. We found that the plantation area that was needed to completely capture a ground water contamination plume was directly proportional to aquifer horizontal hydraulic conductivity, saturated thickness, and ground water gradient. The plantation area needed for capture increased nonlinearly with increasing plume width, aquifer anisotropy, and decreasing growing season duration. The plantation area needed for capture was generally insensitive to aquifer-specific yield and storativity. Steady-state simulations can be used to predict the plantation area needed for capture in many applications. A particularly important finding of this work is that evapotranspiration fluxes through plantations appropriately sized to contain the plume substantially exceeded the ground water flux through the plume itself.     

Fluctuation of groundwater in an urban coastal city of India: a GIS‐based approach

Groundwater is the most important and valuable natural resources especially in coastal urban environment where surface water is insufficient to satisfy the water requirement. Puri city is located on the east coast of India where groundwater is the only source available to meet city water supply. As the city is situated on the sandy aquifer, quality of groundwater is deteriorating because of anthropogenic activities, lack of sewerage system, etc. The objective of the study was to assess the groundwater fluctuation during post‐monsoon and summer with respect to hydrogeological conditions, topography, and groundwater consumption pattern of the city. For this assessment and analysis, Geographic Information System (GIS) was used to visualize topography of the area through digital elevation model (DEM) and distribution of groundwater contours spatially and temporally. The probable areas prone to contamination were identified based on aquifer property and depths to water table below ground. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of water use on arsenic release to well water in a confined aquifer

Field-based experiments were designed to investigate the release of naturally occurring, low to moderate (< 50 microg/L) arsenic concentrations to well water in a confined sandstone aquifer in northeastern Wisconsin. Geologic, geochemical, and hydrogeologic data collected from a 115 m2 site demonstrate that arsenic concentrations in ground water are heterogeneous at the scale of the field site, and that the distribution of arsenic in ground water correlates to solid-phase arsenic in aquifer materials. Arsenic concentrations in a test well varied from 1.8 to 22 microg/L during experiments conducted under no, low, and high pumping rates. The quality of ground water consumed from wells under typical domestic water use patterns differs from that of ground water in the aquifer because of reactions that occur within the well. Redox conditions in the well can change rapidly in response to ground water withdrawals. The well borehole is an environment conducive to microbiological growth, and biogeochemical reactions also affect borehole chemistry. While oxidation of sulfide minerals appears to release arsenic to ground water in zones within the aquifer, reduction of arsenic-bearing iron (hydr)oxides is a likely mechanism of arsenic release to water having a long residence time in the well borehole.     

Delineation of regional arid karstic aquifers: an integrative data approach

This research integrates data procedures for the delineation of regional ground water flow systems in arid karstic basins with sparse hydrogeologic data using surface topography data, geologic mapping, permeability data, chloride concentrations of ground water and precipitation, and measured discharge data. This integrative data analysis framework can be applied to evaluate arid karstic aquifer systems globally. The accurate delineation of ground water recharge areas in developing aquifer systems with sparse hydrogeologic data is essential for their effective long-term development and management. We illustrate the use of this approach in the Cuatrociénegas Basin (CCB) of Mexico. Aquifers are characterized using geographic information systems for ground water catchment delineation, an analytical model for interbasin flow evaluation, a chloride balance approach for recharge estimation, and a water budget for mapping contributing catchments over a large region. The test study area includes the CCB of Coahuila, Mexico, a UNESCO World Biosphere Reserve containing more than 500 springs that support ground water-dependent ecosystems with more than 70 endemic organisms and irrigated agriculture. We define recharge areas that contribute local and regional ground water discharge to springs and the regional flow system. Results show that the regional aquifer system follows a topographic gradient that during past pluvial periods may have linked the Río Nazas and the Río Aguanaval of the Sierra Madre Occidental to the Río Grande via the CCB and other large, currently dry, upgradient lakes.     

Delineating a recharge area for a spring using numerical modeling, Monte Carlo techniques, and geochemical investigation

Recharge areas of spring systems can be hard to identify, but they can be critically important for protection of a spring resource. A recharge area for a spring complex in southern Wisconsin was delineated using a variety of complementary techniques. A telescopic mesh refinement (TMR) model was constructed from an existing regional-scale ground water flow model. This TMR model was formally optimized using parameter estimation techniques; the optimized "best fit" to measured heads and fluxes was obtained by using a horizontal hydraulic conductivity 200% larger than the original regional model for the upper bedrock aquifer and 80% smaller for the lower bedrock aquifer. The uncertainty in hydraulic conductivity was formally considered using a stochastic Monte Carlo approach. Two-hundred model runs used uniformly distributed, randomly sampled, horizontal hydraulic conductivity values within the range given by the TMR optimized values and the previously constructed regional model. A probability distribution of particles captured by the spring, or a "probabilistic capture zone," was calculated from the realistic Monte Carlo results (136 runs of 200). In addition to portions of the local surface watershed, the capture zone encompassed areas outside of the watershed--demonstrating that the ground watershed and surface watershed do not coincide. Analysis of water collected from the site identified relatively large contrasts in chemistry, even for springs within 15 m of one another. The differences showed a distinct gradation from Ordovician-carbonate-dominated water in western spring vents to Cambrian-sandstone-influenced water in eastern spring vents. The difference in chemistry was attributed to distinctive bedrock geology as demonstrated by overlaying the capture zone derived from numerical modeling over a bedrock geology map for the area. This finding gives additional confidence to the capture zone calculated by modeling.     

A subregional-scale method to assess aquifer vulnerability to pesticides

A method to predict aquifer vulnerability to pesticide contamination at the subregional scale was developed. The assessment method was designed to incorporate relevant hydrologic and pesticide-transport information and to use generally available data. The method assumes steady-state advection of pesticides in the vadose zone, including sorption and biological decay. The solution is presented as a vulnerability index (VI) that increases as the aquifer vulnerability increases. The hydrologic input data for the VI model are obtained from the soil survey geographic database. Pesticides were grouped into three leachability classes using a leachability ratio (half-life divided by organic carbon partition coefficient). Pesticide transformation is assumed to occur in the surface layer. The influence of vertical transport in the remainder of the vadose zone has been incorporated by applying a multiplying factor to the VI that varies with depth to ground water. Hydraulic conductivity is used as a surrogate for soil-water velocity for practical purposes. The performance of the VI model is evaluated using ground water data from Weld County, Colorado. The model is demonstrated to be successful at predicting relative vulnerability, defined as the magnitude of pesticide concentration and percent of wells in a unit that exhibit a pesticide detection. Areas of low, medium, and high vulnerability are assigned. Results indicate that the vulnerability classifications are most dependent on the leachability ratio, hydraulic conductivity, and organic carbon content.     

A critical analysis of the cumulative rainfall departure concept

Evaluation of trends in time-series, such as precipitation or ground water levels, is an essential element in many hydrologic evaluations, including water resource studies and planning efforts. The cumulative rainfall departure (CRD) from normal rainfall is a concept sometimes utilized to evaluate the temporal correlation of rainfall with surface water or ground water levels. Permutations of the concept have been used to estimate recharge or aquifer storativity, and in attempts to explain declining ground water levels. The cumulative departure concept has hydrologic meaning in the short term, as a generalized evaluation of either meager or abundant rainfall, and when utilized in connection with a detailed water budget analysis can be used in a predictive fashion. However, the concept can be misapplied if extended over lengthy periods. Misapplication occurs because of several factors including the separation of the mean and median in nonnormal distributions, how the choice of beginning and end points of the data can affect the results, the lack of consideration that above-average rainfall can reset the hydrologic system without mathematically eliminating the accumulated deficit, and the lack of support for the necessary inference that rainfall events and hydrologic levels widely separated in time are linked. Standard statistical techniques are available to reliably determine trends and can provide rigorous statistical measures of the significance of conclusions. Misuse of the CRD concept can lead to erroneous and unsupported conclusions regarding hydrologic relationships and can potentially result in misguided water resource decision-making.     

Neural networks to simulate regional ground water levels affected by human activities

In arid regions, human activities like agriculture and industry often require large ground water extractions. Under these circumstances, appropriate ground water management policies are essential for preventing aquifer overdraft, and thereby protecting critical ecologic and economic objectives. Identification of such policies requires accurate simulation capability of the ground water system in response to hydrological, meteorological, and human factors. In this research, artificial neural networks (ANNs) were developed and applied to investigate the effects of these factors on ground water levels in the Minqin oasis, located in the lower reach of Shiyang River Basin, in Northwest China. Using data spanning 1980 through 1997, two ANNs were developed to model and simulate dynamic ground water levels for the two subregions of Xinhe and Xiqu. The ANN models achieved high predictive accuracy, validating to 0.37 m or less mean absolute error. Sensitivity analyses were conducted with the models demonstrating that agricultural ground water extraction for irrigation is the predominant factor responsible for declining ground water levels exacerbated by a reduction in regional surface water inflows. ANN simulations indicate that it is necessary to reduce the size of the irrigation area to mitigate ground water level declines in the oasis. Unlike previous research, this study demonstrates that ANN modeling can capture important temporally and spatially distributed human factors like agricultural practices and water extraction patterns on a regional basin (or subbasin) scale, providing both high-accuracy prediction capability and enhanced understanding of the critical factors influencing regional ground water conditions.     

Use of Computers in the Field

Performing long-term pump tests, single-well aquifer tests and long-term ground water monitoring were, in the past, considered costly and labor-intensive tasks. However, because of recent advances in portable data acquisition systems, these field tests can be performed in a relatively cost-effective and efficient manner. Certain data acquisition systems (DAS) have unique options available. Recently developed and commonly used computers are the Envirolab®, the Paroscientific®, and the In-Situ®. The use of implementing DAS in ground water field studies has proved to be successful when compared to the many disadvantages of performing those tasks manually. Successful acquisition of aquifer characterization data using computerized methods is highly dependent on the implementation of proper field programming and operational procedures. In order to provide high qualify and meaningful data, the field operator must be familiar with proper equipment implementation and programming instructions. Once familiar with these procedures, the use of DAS in the field provides a cost-effective and efficient tool for aquifer testing.     

Modeling axisymmetric flow and transport

Unmodified versions of common computer programs such as MODFLOW, MT3DMS, and SEAWAT that use Cartesian geometry can accurately simulate axially symmetric ground water flow and solute transport. Axisymmetric flow and transport are simulated by adjusting several input parameters to account for the increase in flow area with radial distance from the injection or extraction well. Logarithmic weighting of interblock transmissivity, a standard option in MODFLOW, can be used for axisymmetric models to represent the linear change in hydraulic conductance within a single finite-difference cell. Results from three test problems (ground water extraction, an aquifer push-pull test, and upconing of saline water into an extraction well) show good agreement with analytical solutions or with results from other numerical models designed specifically to simulate the axisymmetric geometry. Axisymmetric models are not commonly used but can offer an efficient alternative to full three-dimensional models, provided the assumption of axial symmetry can be justified. For the upconing problem, the axisymmetric model was more than 1000 times faster than an equivalent three-dimensional model. Computational gains with the axisymmetric models may be useful for quickly determining appropriate levels of grid resolution for three-dimensional models and for estimating aquifer parameters from field tests.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Delineating and quantifying ground water discharge zones using streambed temperatures

Streambed temperature mapping, hydraulic testing using minipiezometers, and geochemical analyses of interstitial water of the streambed were used to delineate the pattern of ground water discharge in a sandy streambed and to develop a flux-based conceptual model for ground water/surface water interactions. A new and simple empirical method was used to relate fluxes obtained from minipiezometer data to streambed temperatures. The relationship allowed flux to be calculated at locations where only streambed temperature measurements were made. Slug testing and potentiomanometer measurements at 34 piezometers indicated ground water discharge ranged from 0.03 to 446 L/m2/day (and possibly as high as 7060 L/m2/day) along a 60 m long by 11 to 14 m wide reach of river. Complex but similar plan-view patterns of flux were calculated for both summer and winter using hundreds of streambed temperatures measured on a 1 by 2 m grid. The reach was dominated by ground water discharge and 5% to 7% of the area accounted for approximately 20% to 24% of the total discharge. < 12% of the total area consisted of recharge zones or no-discharge zones. A conceptual model for ground water/surface water interactions consisting of five different behaviors was developed based on the magnitude and direction of flux across the surface of the streambed. The behaviors include short-circuit discharge (e.g., high-flow springs), high discharge (e.g., preferential flowpaths), low to moderate discharge, no discharge (e.g., horizontal hyporheic or ground water flow), and recharge. Geological variations at depth played a key role in determining which type of flow behavior occurred in the streambed.