Numerical Simulations of Seasonally Oscillated Groundwater Dynamics in Coastal Confined Aquifers

Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10⁻⁴ m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10⁻⁴ m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10⁻⁴ m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10⁻⁴ m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.     

Measuring a Low Horizontal Hydraulic Gradient in a High Transmissivity Aquifer

High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal-to-noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10⁻⁶ to 10⁻⁵ m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m²/d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10⁻⁶ (±0.9 × 10⁻⁶) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10⁻⁵ (±0.4 × 10⁻⁵) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10⁻⁶ to 10⁻⁵ m/m range.     

A Simple Analytical Formula for the Leakage Flux Through a Perforated Aquitard

In this methods note, we present a simple analytical formula to quantify the steady‐state leakage flux (Q) over a perforated aquitard. The flux depends on the aquitard thickness (D), the radius of the perforation (R), the hydraulic conductivity of the material inside the perforation (k_fill), the conductivities of the overlying and underlying aquifers (k₁ and k₂, respectively), and the head difference between the two aquifers (ΔH): This equation assumes an aquitard separating two homogeneous and infinite aquifers (R ? aquifer thickness) in which radial flow to and from the perforation occurs, with no other recharge or discharge boundaries near the perforation. The flux through a perforation in a hypothetical case study with D = 10 m, k₁ = 10 m/d, k₂ = 20 m/d, R = 0.072 m, and ΔH = 1 m ranges between less than 1 mL/d if the hole is backfilled with bentonite (k_fill = 10^?4 m/d), to several liters per day if the perforation is backfilled with sand from the overlying aquifer (k_fill = 10 m/d), to several m³/d if the perforation forms an open conduit (k_fill = 10⁵ m/d). The leakage fluxes calculated with this model agree well with those calculated using a numerical model (MODFLOW).     

Numerical simulation of three‐dimensional flow hydraulics in a braided channel

Results are presented from a numerical simulation of three‐dimensional flow hydraulics around a mid‐channel bar carried out using the FLUENT/UNS computational fluid dynamics (CFD) software package. FLUENT/UNS solves the three‐dimensional Reynolds‐averaged form of the Navier–Stokes equations. Turbulence closure is achieved using a RNG k–ϵ model. Simulated flow velocities are compared with measured two‐dimensional velocities (downstream and cross‐stream) obtained using an electromagnetic current meter (ECM). The results of the simulation are qualitatively consistent with the flow structures observed in the field. Quantitative comparison of the simulated and measured velocity magnitudes indicates a strong positive correlation between the two (r=0·88) and a mean difference of 0·09 m s⁻¹. Deviations between simulated and measured velocities may be identified that are both random and systematic. The former may reflect a number of factors including subgrid‐scale natural spatial variability in flow velocities associated with local bed structures and measurement uncertainty resulting from problems of ECM orientation. Model mesh configuration, roughness parameterization and inlet boundary condition uncertainty may each contribute to systematic differences between simulated and measured flow velocities. These results illustrate the potential for using CFD software to simulate flow hydraulics in natural channels with complex configurations. They also highlight the need for detailed spatially distributed datasets of three‐dimensional flow variables to establish the accuracy and applicability of CFD software. Copyright © 1999 John Wiley & Sons, Ltd.     

A comparison of hillslope drainage area estimation methods using high-resolution DEMs with implications for topographic studies of gullies

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high-resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high-resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high-resolution data. Specifically, we investigated the impact of single- or multiple-direction flow routing algorithms, DEM hydrologic-enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km² site centred on Weany Creek, a low-relief semi-arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single- or multiple-direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple-direction flow routing algorithm achieve the most realistic drainage area estimates in low-relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion.     

Sediment-transport competence of rain-impacted interrill overland flow

Laboratory experiments to determine the maximum size of sediment transported in shallow, rain-impacted flow were conducted in a recirculating flume 4·80 m long and 0·50 m wide. Rainfall intensities were varied between 51 and 138 mm h⁻¹, flow was introduced from a header tank into the flume at rates ranging from 0 to 0·64 l s⁻¹, and experiments were conducted on gradients between 3·5 and 10°. The following equation was developed: ML = (REFE)^1·6363 in which M is particle mass, L is distance moved in unit time (cm min⁻¹), RE is rainfall energy (J m⁻² s⁻¹) and FE is flow energy (J m⁻² s⁻¹). This equation can be used to predict sediment-transport competence of interrill overland flow. The equation is limited in its utility insofar as it has been developed using quartz grains and takes no account of variations in absorption of rain energy by natural ground surfaces. © 1998 John Wiley & Sons, Ltd.     

Quantitative guidance for efficient vertical flow measurements at the sediment–water interface using temperature–depth profiles

Upward discharge to surface water bodies can be quantified using analytical models based on temperature–depth (T-z) profiles. The use of sediment T-z profiles is attractive as discharge estimates can be obtained using point-in-time data that are collected inexpensively and rapidly. Previous studies have identified that T-z methods can only be applied at times of the year when there is significant difference between the streambed–water interface and deeper sediment temperatures (e.g., winter and summer). However, surface water temperatures also vary diurnally, and the influence of these variations on discharge estimates from T-z methods is poorly understood. For this study, synthetic T-z profiles were generated numerically using measured streambed interface temperature data to assess the influence of diurnal temperature variations on discharge estimation and provide insight into the suitable application of T-z methods. Results show that the time of day of data collection can have a substantial influence on vertical flux estimates using T-z methods. For low groundwater discharge fluxes (e.g., 0.1 m d⁻¹), daily transience in streambed temperatures led to relatively large errors in estimated flow magnitude and direction. For higher discharge fluxes (1.5 m d⁻¹), the influence of transient streambed temperatures on discharge estimates was strongly reduced. Discharge estimates from point-in-time T-z profiles were most accurate when the uppermost point in the T-z profile was near the bed interface daily mean (two time periods daily). Where temperature time series data are available, daily averaged T-z profiles can produce accurate discharge estimates across a wide range of discharge rates. Seasonality in shallow groundwater temperature generally had a negligible influence on vertical flow estimates. These findings can be used to plan field campaigns and provide guidance on the optimal application of T-z methods to quantify vertical groundwater discharge to surface water bodies.     

Effectiveness of post-fire salvage logging stream buffer management for hillslope erosion in the U.S. Inland Northwest Mountains

Active wildfire seasons in the western U.S. warrant the evaluation of post-fire forest management strategies. Ground-based salvage logging is often used to recover economic loss of burned timber. In unburned forests, ground-based logging often follows best management practices by leaving undisturbed areas near streams called stream buffers. However, the effectiveness of these buffers has not been tested in a post-wildfire setting. This experiment tested buffer width effectiveness with a novel field-simulated rill experiment using sediment-laden runoff (25 g/L) released over 40 min at evenly timed flow rates (50, 100 and 150 L/min) to measure surface runoff travel length and sediment concentration under unburned and high and low soil burn severity conditions at 2-, 10- and 22-month post-fire. High severity areas 2-month post-fire had rill lengths of up to 100 m. Rill length significantly decreased over time as vegetation regrowth provided ground cover. Sediment concentration and sediment dropout rate also varied significantly by soil burn severity. Sediment concentrations were 19 g/L for the highest flow 2-month post-fire and reduced to 6.9–14 g/L 10-month post-fire due to abundant vegetation recovery. The amount of sediment dropping out of the flow consistently increased over the study period with the low burn severity rate of 1.15 g L⁻¹ m⁻¹ approaching the unburned rate of 1.29 g L⁻¹ m⁻¹ by 2-year post-fire. These results suggest that an often-used standard, 15 m buffer, was sufficient to contain surface runoff and reduce sediment concentration on unburned sites, however buffers on high burn severity sites need to be eight times greater (120 m) immediately after wildfire and four times greater (60 m) 1-year post-fire. Low burn severity areas 1-year post-fire may need to be only twice the width of an unburned buffer (30 m), and 2-year post-fire these could return to unburned widths.     

Simulation Assessment of Direct Push Injection Logging for High-Resolution Aquifer Characterization

Direct push injection logging (DPIL) has become one of the most widely used approaches for obtaining vertical profiles of hydraulic conductivity (K) in environmental site investigations. Despite its widespread use, however, there has been no rigorous analysis of the underlying physical processes that take place during DPIL or how the approach would perform under different hydrogeological and operating conditions. We address these issues through a series of numerical simulations. Results show that the ratio of DPIL injection rate over pressure can be used for direct determination of K when K is >10⁻⁶ m/s. When K is <10⁻⁶ m/s and specific storage (Ss) is >10⁻³/m, the ratio becomes increasingly sensitive to Ss; in that case, additional information on Ss is needed for reliable K estimation. For unconsolidated formations of moderate K or higher, the ratio of injection rate over pressure should provide a reasonable K estimate when Ss is <10⁻³/m. Although water injection at previous depths during continuous DPIL has only a small impact on the pressure response measured at the current injection depth, probe advancement can have a significant impact when K and Ss are small. Consequently, in fine-grained materials, the advancement-generated pore water pressure increase can comprise a large portion of the measured pressure response. To diminish the impact of probe advancement in such materials, advancement speed should be kept as low as possible (e.g., 0.5 cm/s).     

Rate of sedimentation in Lake Kinneret,Israel: spatial and temporal variations

The rate of sedimentation in Lake Kinneret was measured over several years by means of sediment traps, in up to seven different locations in the lake. Gross sedimentation rates measured in the sediment traps vary from about 1·5 kg m⁻² a⁻¹ in the deepest part of the lake up to 10 kg m⁻² a⁻¹ near the mouth of the upper Jordan river. The rate of sedimentation near the Jordan's inflow is highly correlated to flow discharge in the river, while in the centre of the lake the seasonal sedimentation pattern is mainly correlated to the bloom period of Peridinium gatunense. During the bloom period of Peridinium gatunense sedimentation rates all over the lake are very similar, indicating that the Peridinium is evenly distributed in the lake. The average suspended sediment discharge of the upper Jordan river flowing into the lake is 41 000 ton a⁻¹.Copyright © 2000 John Wiley & Sons, Ltd.     

Secondary airflow and sediment transport in the lee of a reversing dune

Lee-side windspeed and sediment transport were measured over a small (1·2 m) transverse ridge in the Silver Peak dunefield, west-central Nevada, USA, using an intensive array of 25 cup anemometers and seven total flux traps. During crest-transverse and transporting flow conditions (u_0·3crest ≈ 8·4 m s⁻¹), windspeed near the surface of the lee slope averaged half (48 per cent) that of crest speeds. Dimensionless speeds in the separation zone ranged from 0·2 to 0·8 that of the outer flow (u₁₂). Along the boundary of the separation cell, windspeed increased by 10 per cent of the crest speed before separation. Equilibrium of upper and lower wake regions was not observed by the documented eight dune heights, suggesting that wake recovery may not occur over closely spaced dunes. Sediment transport measured directly on both the lee slope and interdune surfaces averaged approximately 15 per cent of crest inputs. This suggests that a significant amount (c. 70–95 per cent) of sediment transported over the crest moved as fallout. For this data set, flux was approximately proportional to the cube of the near-surface windspeed (u_0·3) and in general there was an order of magnitude difference between flux measured at the crest and that measured within the separation zone. Transport direction in the separation zone was acutely oblique to the incident direction owing to secondary flow deflection. Beyond the interdune, transport direction progressed from oblique to crest-transverse. This indicates that an appreciable amount of sediment may move laterally along the lee slope and interdune corridor under crest-transverse flows. Regarding the grain size and sorting properties of transported sediment, there was no significant difference in mean grain size over the dune, although in general particles were finer and more poorly sorted in the lee. Copyright © 1999 John Wiley & Sons, Ltd.     

Sedimentation in arctic proglacial lakes: Mittivakkat Glacier,south‐east Greenland

Several sediment cores were collected from two proglacial lakes in the vicinity of Mittivakkat Glacier, south‐east Greenland, in order to determine sedimentation rates, estimate sediment yields and identify the dominant sources of the lacustrine sediment. The presence of varves in the ice‐dammed Icefall Lake enabled sedimentation rates to be estimated using a combination of X‐ray photography and down‐core variations in ¹³⁷Cs activity. Sedimentation rates for individual cores ranged between 0·52 and 1·06 g cm⁻² year⁻¹, and the average sedimentation rate was estimated to be 0·79 g cm⁻² year⁻¹. Despite considerable down‐core variability in annual sedimentation rates, there is no significant trend over the period 1970 to 1994. After correcting for autochthonous organic matter content and trap efficiency, the mean fine‐grained minerogenic sediment yield from the 3·8 km² basin contributing to the lake was estimated to be 327 t km⁻² year⁻¹. Cores were also collected from the topset beds of two small deltas in Icefall Lake. The deposition of coarse‐grained sediment on the delta surface was estimated to total in excess of 15 cm over the last c. 40 years. In the larger Lake Kuutuaq, which is located about 5 km from the glacier front and for which the glacier represents a smaller proportion of the contributing catchment, sedimentation rates determined for six cores collected from the centre of the lake, based on their ¹³⁷Cs depth profiles, were estimated to range between 0·05 and 0·11 g cm⁻² year⁻¹, and the average was 0·08 g cm⁻² year⁻¹. The longer‐term (c. 100–150 years) average sedimentation rate for one of the cores, estimated from its unsupported ²¹⁰Pb profile, was 0·10–0·13 g cm⁻² year⁻¹, suggesting that sedimentation rates in this lake have been essentially constant over the last c. 100–150 years. The average fine‐grained sediment yield from the 32·4 km² catchment contributing to the lake was estimated to be 13 t km⁻² year⁻¹. The ¹³⁷Cs depth profiles for cores collected from the topset beds of the delta of Lake Kuutuaq indicate that in excess of 27 cm of coarse‐grained sediment had accumulated on the delta surface over the last approximately 40 years. Caesium‐137 concentrations associated with the most recently deposited (uppermost) fine‐grained sediment in both Icefall Lake and Lake Kuutuaq were similar to those measured in fine‐grained sediment collected from steep slopes in the immediate proglacial zone, suggesting that this material, rather than contemporary glacial debris, is the most likely source of the sediment deposited in the lakes. This finding is confirmed by the ¹³⁷Cs concentrations associated with suspended sediment collected from the Mittivakkat stream, which are very similar to those for proglacial material. Copyright © 2000 John Wiley & Sons, Ltd.     

Flow over a rounded backward-facing step,using a <Emphasis Type="Italic">z</Emphasis>-coordinate model and a <Emphasis Type="Italic">σ</Emphasis>-coordinate model

Homogeneous, nonrotating flow over a backward-facing rounded step is simulated using the 2D vertical version of two general circulation models, a z-coordinate model—the Massachusetts Institute of Technology general circulation model (MITgcm)—and a σ-coordinate model—the Bergen Ocean Model (BOM). The backward-facing step is a well-known testcase since it is geometrically simple but still embodies important flow characteristics such as separation point, reattachment length, and recirculation of the flow. The study compares the core of the two models and uses constant eddy viscosities and diffusivities. The Reynolds numbers ranges from 2·10² to 2·10⁶. The results correspond with previously published results having a relatively stationary separation point and a fluctuating reattachment length due to downslope propagating eddies released from the reattachment zone for Reynolds numbers higher than or equal to 2 · 10⁴. For Reynolds number within the laminar regime, the flow is stationary. The discrepancies between the models increase by enhancing Reynolds numbers. The σ-coordinate model experiences a reduction in eddy sizes with increasing resolution and Reynolds numbers in correspondence with published experiments, while the size of the eddies are independent of the Reynolds number using the MITgcm. Due to mixing generated by the staircase topography, the z-coordinate model gives a better convergence of the separation point and reattachment length compared with the BOM; however, this conclusion might change with the inclusion of a relevant turbulence scheme.     

The impact of flow discharge on the hydraulic characteristics of headcut erosion processes in the gully region of the Loess Plateau

Headcut erosion is associated with major hydraulic changes induced by the gully head of concentrated flow. However, the variation in the hydraulic characteristics of the headcut erosion process is still not clear in the gully region of the Loess Plateau. A series of rainfall combined scouring experiments (flow discharges ranging from 3.6 to 7.2 m³ hr⁻¹, with 0.8 mm min⁻¹ rainfall intensity) were conducted on experimental plots to clarify the variation in the hydraulic parameters induced by gully head and erosion processes under different flow discharges. The results showed that concentrated flows in the catchment area and gully bed were turbulent (Reynolds number ranging from 1,876 to 6,693) and transformed between supercritical and subcritical (Froude number ranging from 0.96 to 3.73). The hydraulic parameters, such as the flow velocity, Reynolds number, shear stress, stream power, Darcy–Weisbach friction factor, and unit stream power in the catchment area were 0.45–0.59 m s⁻¹, 2086–6693, 1.96–5.33 Pa, 0.89–2.86 W m⁻², 0.08–0.16, and 0.023–0.031 m s⁻¹, respectively. When the concentrated flows dropped from the gully head, the hydraulic parameters in the gully bed decreased by 3.39–26.07%, 1.49–29.99%, 65.19–67.14%, 67.25–74.96%, 28.53–61.31%, and 67.82–77.14%, respectively, which contributed to the flow energy consumption at the gully head. As flow discharge increased, Reynolds number, shear stress, and stream power increased, while flow velocity, Froude number, unit stream power, and Darcy–Weisbach friction factor did not. The flow energy consumption at the gully head was 9.66–10.13, 13.25–13.74, 15.68–16.41, and 19.28–20.25 J s⁻¹, respectively, under different flow discharges and accounted for 60.58–68.50% of the flow energy consumption of the experimental plots. Generally, the sediment discharges increased rapidly at the initial stage, then increased slowly, and finally reached a steady state condition, which showed a significant declining logarithmic trend with experimental duration (P<.01) and increased with increasing flow discharge. Accordingly, the flow energy consumption was significantly correlated with the sediment yield. These findings could improve our understanding of the hydraulic properties and flow energy characteristics of headcut erosion.     

Quantification of the hydraulic diffusivity of a bentonite-sand mixture using the water head decrease measured upon sudden flow interruption

This paper presents a new modelling approach to quantify the hydraulic diffusivity of low-permeability unconsolidated porous media under confined saturated-flow conditions in the laboratory. The derived analytical solution for the transient variation of the hydraulic head after flow interruption was applied to experimental data obtained from continuous measurements of the water pressure at two locations in the soil column. Three soil samples made of a mixture of natural bentonite (at different mass fractions) and medium sand were studied during a series of stepwise constant flow rates of water. The numerical results well fit the experimentally measured decrease of the dimensionless hydraulic head. The study shows that the increase of the mass fraction of bentonite in the soil sample from 10 to 30% is accompanied by a strong decrease of the hydraulic diffusivity from 2.4 × 10⁻² to 1.1 × 10⁻³ m² s⁻¹, which is clearly due to the decrease of the hydraulic conductivity of the soil sample. The specific storages obtained for each of the three samples are in the same order of magnitude and seem to decrease with the increase of mass fraction of bentonite. However, they clearly reflect the predominant portion of the compressibility of the porous medium compared with that of water. Compared with reported literature values for clayey soils, the specific storage values in this study are slightly higher, varying within the range of 2 × 10⁻³ to 8.1 × 10⁻³ m⁻¹._. The experimental results also give insight into the limitations of the modelling approach. In the case of low-permeability soils (K < 2 × 10⁻⁶ ms⁻¹) and steady-flow conditions with low Reynolds numbers, for example, Re < 0.003, it is recommended to choose a time duration for flow interruption between subsequent flow rate steps of longer than 5 s. For high-permeability porous media, to increase the precision of the quantified hydraulic diffusivity, it might be useful to select a measuring frequency significantly higher than 1 Hz.     

Potentials and limitations of modelling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil

Root water uptake patterns are often studied with simulation models of the unsaturated soil water flow, as they are difficult to measure directly. Calibration of these models is not straightforward and causes uncertainties in simulated uptake distributions. In this paper we study how uncertainties in the calibration of the SWIF model affect uncertainty intervals in simulated uptake patterns of an Austrian pine stand (Pinus nigra var. nigra) on a sandy soil. After calibrating and validating SWIF with a large data set of more than 125 000 measured soil water contents over a three year period, uncertainty ranges in simulated soil water dynamics and root water uptake distributions were estimated with a Monte Carlo analysis. In general, uncertainties in root uptake patterns were small (typically <2 10⁻⁴ m³ m⁻³ day⁻¹) and were higher for trees with a shallow rooting system (0·8 m) than for trees with a deep rooting system (2·5 m). Uncertainties arose mainly from uncertainties in simulated soil water fluxes and from variations in the reduction of uptake during periods of drought. Uncertainties in soil water contents were far higher (typically 0·01 m³ m⁻³) than uncertainties in uptake, illustrating that uncertainties in uptake parameters and those in the distribution of water uptake hardly affect the modelling of soil water dynamics. Root water uptake models should therefore be validated against measured uptake distributions, which can be determined on sandy soils during dry periods with a high water use when soil fluxes are negligible to uptake. Copyright © 2000 John Wiley & Sons, Ltd.     

Rill erosion processes on a constantly saturated slope

Rill erosion processes on saturated soil slopes are important for understanding erosion hydrodynamics and determining the parameters of rill erosion models. Saturated soil slopes were innovatively created to investigate the rill erosion processes. Rill erosion processes on saturated soil slopes were modelled by using the sediment concentrations determined by sediment transport capacities (STCs) measurement and the sediment concentrations at different rill lengths. Laboratory experiments were performed under varying slope gradients (5°, 10°, 15°, and 20°) and unit-width flow rates (0.33, 0.67, and 1.33 × 10⁻³ m³ s⁻¹ m⁻¹) to measure sediment concentrations at different rill lengths (1, 2, 4, and 8 m) on saturated soil slopes. The measured sediment concentrations along saturated rills ranged from 134.54 to 1,064.47 kg/m³, and also increased exponentially with rill length similar to non-saturated rills. The model of the rill erosion process in non-saturated soil rills was applicable to that in saturated soil rills. However, the sediment concentration of the rill flow increased much faster, with the increase in rill length, to considerably higher levels at STCs. The saturated soil rills produced 120–560% more sediments than the non-saturated ones. Moreover, the former eroded remarkably faster in the beginning section of the rills, as compared with that on the non-saturated soil slopes. This dataset serves as the basis for determining the erosion parameters in the process-based erosion models on saturated soil slopes.     

A stream tracer technique employing ionic tracers and specific conductance data applied to the Maimai catchment,New Zealand

The stream tracer technique and transient storage models (TSMs) have become common tools in stream solute and hyporheic exchange studies. The expense and logistics associated with water sample collection and analysis often results in limited temporal resolution of stream tracer breakthrough curves (BTCs). Samples are often collected without a priori or real‐time knowledge of BTC information, which can result in poor sample coverage of the critical shoulder (initial rise) and tail (post‐steady state fall) of the BTC. We illustrate the use of specific conductance (SC) measurements as a surrogate for conservative dissolved tracer (Br⁻) samples. The advantages of collecting SC data for use in the TSM are (1) cost, (2) ease of data collection, and (3) well‐defined breakthrough curves, which strengthen TSM parameter optimization. This method is based on developing an ion concentration (IC)–SC relationship from limited discrete tracer solute samples. SC data can be collected on a more frequent basis at no additional analysis cost. TSM simulations can then be run for the conservative tracer data derived from SC breakthrough curves and the IC–SC relationship. This technique was tested in a 120 m reach of stream (2–60 m subreaches) in the Maimai M15 catchment, New Zealand during baseflow recession. Dissolved LiBr was injected for 12·92 h, with Br⁻ as the conservative ion of interest. Four TSM simulations using the OTIS model are optimized using UCODE to fit (1) Br⁻ data derived from the Br⁻–SC relationship (n = 1307 observations at each of two stream sampling sites), (2) all stream Br⁻ data collected (n = 58 in upper reach, n = 60 in lower reach), (3) half of the stream Br⁻ data collected, and (4) 20 stream Br⁻ samples from each site. No two simulations resulted in the same optimal parameter values. Results suggest that the greater the frequency of observations, the greater the confidence in estimated parameter values. Br⁻–SC simulations resulted in the best overall model fits to the data, with the lowest calculated error variance of 6·37, narrowest 95% parameter estimate confidence intervals, and the highest correlation coefficient of 0·99 942, among the four simulations. This is largely due to the improved representation of the shoulder and tail of the BTC with this method. The IC–SC correlation method is robust in situations in which (1) changes in background SC data can be accounted for, and (2) the data used to define the IC–SC relationship are representative of the range of data collected. This method provides more efficient sample analysis, improved data resolution, and improved model results compared to the alternative stream tracer data gathering methods presented. Additionally, we describe a new parameterization of the cross‐sectional area of the stream during flow recession, as a function of discharge, based on a stream hydraulic geometry relationship. This variant of the OTIS model provides a more realistic representation of stream dynamics during unsteady discharge. Copyright © 2005 John Wiley & Sons, Ltd.     

Chloride source delineation in an urban-agricultural watershed: Deicing agents versus agricultural contributions

Analyses (n = 525) of chloride (Cl⁻), bromide (Br⁻), nitrate as nitrogen (NO₃-N), sodium (Na⁺), calcium (Ca²⁺) and potassium (K⁺) in stream water, tile-drain water and groundwater were conducted in an urban-agricultural watershed (10% urban/impervious, 87% agriculture) to explore potential differences in the signature of Cl⁻ originating from an urban source as compared with an agricultural source. Only during winter recharge events did measured Cl⁻ concentrations exceed the 230 mg/L chronic threshold. At base flow, nearly all surface water and tile water samples had Cl⁻ concentrations above the calculated background threshold of 18 mg/L. Mann–Whitney U tests revealed ratios of Cl⁻ to Br⁻ (p = .045), to NO₃-N (p < .0001), to Ca²⁺ (p < .0001), and to Na⁺ (p < .0001) to be significantly different between urban and agricultural waters. While Cl⁻ ratios indicate that road salt was the dominant source of Cl⁻ in the watershed, potassium chloride fertilizer contributed as an important secondary source. Deicing in watersheds where urban land use is minimal had a profound impact on Cl⁻ dynamics; however, agricultural practices contributed Cl⁻ year-round, elevating stream base flow Cl⁻ concentrations above the background level.     

A comparison of the effectiveness of the roots of two grass species in reducing soil erosion on alluvial fans in south-east China

The effects of root systems on soil detachment by overland flow are closely related to vegetation types. The objective of this study was to quantify the effects of two gramineous roots (Paspalum mandiocanum with shallow roots and Pennisetum giganteum with deep roots) on soil detachment capacity, rill erodibility, and critical shear stress on alluvial fans of benggang in south-east China. A 4-m-long and 0.12-m-wide flume was used. Slope steepness ranged from 9% to 27%, and unit flow discharge ranged from 1.39 × 10⁻³ to 4.19 × 10⁻³ m² s⁻¹. The mean detachment capacities of P. mandiocanum and P. giganteum lands were 18% and 38% lower than that of bare land, respectively, and the effects of root on reducing soil detachment were mainly reflected in the 0- to 5-cm soil layer. The most important factors in characterizing soil detachment capacity were root length density and soil cohesion, and soil detachment capacity of the two grass lands could be estimated using flow shear stress, soil cohesion, and root length density (NSE = 0.90). With the increase in soil depth, rill erodibility increased, whereas shear stress decreased. The mean rill erodibilities of P. mandiocanum and P. giganteum lands were 81% and 61% as much as that of bare land, respectively. Additionally, rill erodibilities of the two grass lands could be estimated as an exponential function by root length density and soil cohesion (NSE = 0.88). The mean critical shear stress of P. mandiocanum and P. giganteum lands was 1.29 and 1.39 times that of bare land, respectively, and it could be estimated with a linear function by root length density (NSE = 0.76). This study demonstrated that planting of the two grasses P. mandiocanum and P. giganteum could effectively reduce soil detachment and enhance soil resistance to erosion on alluvial fans, with the deep roots of P. giganteum being more effective than the shallow roots of P. mandiocanum. The results are helpful for understanding the influencing mechanism of root systems on soil detachment process.