Chemistry of streams draining grassland and forest catchments at Plynlimon,mid-Wales

Abstract

The chemistry of streamwater, bulk precipitation, throughfall and soil waters has been studied for three years in two plantation forest and two moorland catchments in mid-Wales. Na and CI are the major ions in streamwater reflecting the maritime influence on atmospheric inputs. In all streams, baseflow is characterised by high pH waters enriched in Ca, Mg, Si and HCO³. Differences in baseflow chemistry between streams reflect the varying extent of calcite and base metal sulphide mineralization within the catchments. Except for K, mean stream solute concentrations are higher in the unmineralized and mineralized forest catchments compared with their respective grassland counterparts. In the forest streams, storm flow concentrations of H⁺ are approximately 1.5 times and Al four times higher than in the moorland streams. Annual catchment losses of Na, Cl, SO⁴, NO³, Al and Si are greatest in the forest streams. In both grassland and forest systems, variations in stream chemistry be explained by mixing waters from different parts of the catchment, although NO³ concentrations may additionally be controlled by N transformations occurring between soils and streams. Differences in stream chemistry and solute budgets between forest and moorland catchments are related to greater atmospheric scavenging by the trees and changes in catchment hydrology consequent on afforestation. Mineral veins within the catchment bedrock can significantly modify the stream chemical response to afforestation.     

Temporal variability of solute transport under vadose zone conditions

The heterogeneity of the solute flux field in the horizontal plane at the field scale has been documented in several field studies. On the other hand, little information is available on the persistence of certain solute transport scenarios over consecutive infiltration cycles. This study was initiated to analyse the recurrence of solute leaching behaviour as estimated in two soil column tests emphasizing the preferential flow phenomenon. Twenty-four small-sized soil samples were subjected to two consecutive unsaturated steady-state flow leaching experiments with bromide as tracer. Observed breakthrough curves (BTCs) were analysed by the method of moments and by the advection–dispersion equation (ADE) to classify solute behaviour. Frequency distributions of the parameters indicating the solute velocity were heavily skewed or bimodal, reflecting the broad variability of the leaching scenarios, including some with pronounced preferential solute breakthrough. Exclusion of the preferential flow columns from our calculations revealed an average amount of 37% of immobile water. The large-scale BTCs derived from assembling the individual concentration courses of each run showed similar features, such as an early bromide breakthrough. However, two distinct apices, viz. one preferential and one matrix, were observed only in the first run, whereas the concentration decrease between the peaks was missing from the second run. A change in soil structure with continuous leaching was presumed to modify the interplay of the various flow domains, thereby altering the spreading of the BTCs. Correlation analysis between parameters of both tests suggests that preferential transport conditions are likely to occur at the same locations in the field over several infiltration cycles, whereas the ‘classical’ or expected matrix flow is time variant and therefore seems to be hardly predictable. © 1998 John Wiley & Sons, Ltd.     

Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.     

Ground-Water Chemistry of Dibdiba Formation, North Kuwait

Chemical studies have been carried out on a number of water wells from the Dibdiba Formation northeast of Kuwait. Water salinity of this formation ranges between 3,300 mg/l to 7,000 mg/l, increasing with depth. The water entrapped in Dibdiba Formation is mainly sodium chloride type which can be differentiated into three different groups according to the ranges of concentration of the main cations and anions. These groups are: (3331113) sodium chloride water type in which Cl > Na, group (3333113) sodium chloride water type in which Na > Cl. In both groups the sequence of dominant cations is Na > Ca > Mg. Group (3333111) sodium chloride water type has Na > Cl and the sequence of dominant cations is Na > Mg > Ca. Chemical ratios of Ca/Mg, Na/Cl, and C1/HCO₃ were worked out for Dibdiba ground water. The ratio of Ca/Mg Dibdiba Formation ranges from 0.4 to 8.58, the ratio of Na/Cl ranges between 0.98 to 1.33, and the ratio of C1/HCO₃ is 232. A plot of chemical analysis on a trilinear diagram shows that Dibdiba Formation ground-water chemical properties are dominated by alkalies (Na > Ca) and strong acid (Cl > SO₄). Water chemistry may reflect the history of the flow path, indicating regional flow as shown by increasing Na⁺, Cl^-, SO^–₄ and where Ca⁺, Mg⁺ achieve equilibrium.     

Tracing solute infiltration using a combined method of dye tracer test and electrical resistivity tomography in an undisturbed forest soil profile

An accurate prediction of solute infiltration in a soil profile is important in the area of environmental science, groundwater and civil engineering. We examined the infiltration pattern and monitored the infiltration process using a combined method of dye tracer test and electrical resistivity tomography (ERT) in an undisturbed field soil (1 m × 1 m). A homogeneous matrix flow was observed in the surface soil (A horizon), but a preferential flow along macropores and residual rock structure was the dominant infiltration pattern in the subsurface soil. Saturated interflow along the slopping boundaries of A and C1 horizons and of an upper sandy layer and a lower thin clay layer in the C horizon was also observed. The result of ERT showed that matrix flow started first in A horizon and then the infiltration was followed by the preferential flows along the sloping interfaces and macropores. The ERT did not show as much detail as the dye‐stained image for the preferential flow. However, the area with the higher staining density where preferential flow was dominant showed a relatively lower electrical resistivity. The result of this study indicates that ERT can be applied for the monitoring of solute transportation in the vadose zone. Copyright © 2010 John Wiley & Sons, Ltd.     

Tracer‐based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope

Runoff from boreal hillslopes is often affected by distinct soil boundaries, including the frozen boundary and the organic‐mineral boundary (OMB), where highly porous and hydraulically conductive organic material overlies fine‐grained mineral soils. Viewed from the surface, ground cover appears as a patchwork on sub‐meter scales, with thick, moss mats interspersed with lichen‐covered, silty soils with gravel inclusions. We conducted a decameter‐scale subsurface tracer test on a boreal forest hillslope in interior Alaska to quantify locations and mechanisms of transport and storage in these soils, focusing on the OMB. A sodium bromide tracer was added as a slug addition to a pit and sampled at 40 down‐gradient wells, screened primarily at the OMB and within a 7 × 12 m well field. We maintained an elevated head in the injection pit for 8.5 hr to simulate a storm. Tracer breakthrough velocities ranged from <0.12 to 0.93 m hr^?1, with the highest velocities in lichen‐covered soils. After 12 hr and cessation of the elevated head, the tracer coalesced and was only detected in thick mosses at a trough in the OMB. By 24 hr, approximately 17% of the tracer mass could be accounted for. The majority of the mass loss occurred between 4 and 12 hr, while the tracer was in contact with lichen‐covered soils, which is consistent with tracer transport into deeper flow paths via preferential flow through discrete gravelly areas. Slow breakthroughs suggest that storage and exchange also occurred in shallow soils, likely related to saturation and drainage in fine‐grained mineral soils caused by the elevated hydraulic head. These findings highlight the complex nature of storage and transmission of water and solutes from boreal hillslopes to streams and are particularly relevant given rapid changes to boreal environments related to climate change, thawing permafrost and increasing fire severity.     

Comparing dispersivities and soil chloride concentrations of turfgrass-covered undisturbed and disturbed soil columns

Research relating to soil leaching properties under turfgrass conditions has often been conducted on disturbed soils where macropore structure has been destroyed. The objective of this study was to compare the solute movement characteristics of undisturbed and disturbed soil columns covered with turfgrass. Dispersivities and chloride (Cl) breakthrough curves of undisturbed and disturbed soils were investigated. Soil columns were excavated into three sections after testing, for which the mean bulk density was 1.33 Mg M⁻³ for the undisturbed columns and 1.16 Mg m⁻³ for the disturbed columns. The dispersivity for the undisturbed columns was over three times greater than for the disturbed columns. Chloride concentration found in Layer 1 (0–6.7 cm), Layer 2 (6.7–13.4 cm), and Layer 3 (13.4–20.0 cm) were 2.8, 5.3, and 4.8 times higher, respectively, for the disturbed soils than for the undisturbed. Applying conclusions from solute movement studies using repacked columns covered with turfgrass to actual undisturbed field conditions could lead to errors in interpretation because of the effect of macropores.     

Variability of unsaturated Bromide fluxes as measured through a layered volcanic vadose zone in New Zealand

The measured drainage fluxes through a layered volcanic vadose zone exhibited high spatial variability as a consequence of heterogeneous flow conditions. The drainage flux variability was quantified using automated equilibrium tension lysimeters, installed in close‐proximity and resulted in high variability in the Br masses recovered from a conservative tracer experiment. The primary cause of the heterogeneous flow was attributed to textural changes occurring at the interface between volcanic layers, resulting in development of funnel‐flow patterns, and further enhanced by the existence of hydrophobic conditions. The Br recoveries in individual automated equilibrium tension lysimeters were used to determine the corresponding variable sizes of the surface areas contributing drainage to the lysimeters. The tracer experiment confirmed the existence of unsaturated lateral transport occurring at the sloping interface of the coarse Taupo Ignimbrite material with the silty Palaeosol layer at approximately 4.2 m depth. This study demonstrates that measurements of both flux and solute concentrations at multiple locations are essential when heterogeneous flow is suspected to be present, to be able to determine reliable estimates of contaminant leaching through the vadose zone at the plot scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of antecedent soil moisture on hydraulic conductivity in a series of texture‐contrast soils

Antecedent soil moisture significantly influenced the hydraulic conductivity of the A1, A2e and B21 horizons in a series of strong texture‐contrast soils. Tension infiltration at six supply potentials demonstrated that in the A1 horizon, hydraulic conductivity was significantly lower in the ‘wet’ treatment than in the ‘dry’ treatment. However in the A2e horizon, micropore and mesopore hydraulic conductivity was lower in the ‘dry’ treatment than the ‘wet’ treatment, which was attributed to the precipitation of soluble amorphous silica. In the B21 horizon, desiccation of vertic clays resulted in the formation of shrinkage cracks which significantly increased near‐saturated hydraulic conductivity and prevented the development of subsurface lateral flow in the ‘dry’ treatment. In the ‘wet’ treatment, the difference between the hydraulic conductivity of the A1 and B21 horizons was reduced; however, lateral flow still occurred in the A1 horizon due to difficulty displacing existing soil water further down the soil profile. Results demonstrate the need to account for temporal variation in soil porosity and hydraulic conductivity in soil‐water model conceptualisation and parameterisation. Copyright © 2011 John Wiley & Sons, Ltd.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite-element analysis of the transport of water and solutes in title-drained soils

The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems.     

A Geochemical Approach to Determine Sources and Movement of Saline Groundwater in a Coastal Aquifer

Geochemical evaluation of the sources and movement of saline groundwater in coastal aquifers can aid in the initial mapping of the subsurface when geological information is unavailable. Chloride concentrations of groundwater in a coastal aquifer near San Diego, California, range from about 57 to 39,400 mg/L. On the basis of relative proportions of major‐ions, the chemical composition is classified as Na‐Ca‐Cl‐SO₄, Na‐Cl, or Na‐Ca‐Cl type water. δ²H and δ¹⁸O values range from ?47.7‰ to ?12.8‰ and from ?7.0‰ to ?1.2‰, respectively. The isotopically depleted groundwater occurs in the deeper part of the coastal aquifer, and the isotopically enriched groundwater occurs in zones of sea water intrusion. ⁸⁷Sr/⁸⁶Sr ratios range from about 0.7050 to 0.7090, and differ between shallower and deeper flow paths in the coastal aquifer. ³H and ¹⁴C analyses indicate that most of the groundwater was recharged many thousands of years ago. The analysis of multiple chemical and isotopic tracers indicates that the sources and movement of saline groundwater in the San Diego coastal aquifer are dominated by: (1) recharge of local precipitation in relatively shallow parts of the flow system; (2) regional flow of recharge of higher‐elevation precipitation along deep flow paths that freshen a previously saline aquifer; and (3) intrusion of sea water that entered the aquifer primarily during premodern times. Two northwest‐to‐southeast trending sections show the spatial distribution of the different geochemical groups and suggest the subsurface in the coastal aquifer can be separated into two predominant hydrostratigraphic layers.     

Water transport and tracer mixing in volcanic ash soils at a tropical hillslope: A wet layered sloping sponge

Andosol soils formed in volcanic ash provide key hydrological services in montane environments. To unravel the subsurface water transport and tracer mixing in these soils we conducted a detailed characterization of soil properties and analyzed a 3-year data set of sub-hourly hydrometric and weekly stable isotope data collected at three locations along a steep hillslope. A weakly developed (52–61 cm depth), highly organic andic (Ah) horizon overlaying a mineral (C) horizon was identified, both showing relatively similar properties and subsurface flow dynamics along the hillslope. Soil moisture observations in the Ah horizon showed a fast responding (few hours) “rooted” layer to a depth of 15 cm, overlying a “perched” layer that remained near saturated year-round. The formation of the latter results from the high organic matter (33–42%) and clay (29–31%) content of the Ah horizon and an abrupt hydraulic conductivity reduction in this layer with respect to the rooted layer above. Isotopic signatures revealed that water resides within this soil horizon for short periods, both at the rooted (2 weeks) and perched (4 weeks) layer. A fast soil moisture reaction during rainfall events was also observed in the C horizon, with response times similar to those in the rooted layer. These results indicate that despite the perched layer, which helps sustain the water storage of the soil, a fast vertical mobilization of water through the entire soil profile occurs during rainfall events. The latter being the result of the fast transmissivity of hydraulic potentials through the porous matrix of the Andosols, as evidenced by the exponential shape of the water retention curves of the subsequent horizons. These findings demonstrate that the hydrological behavior of volcanic ash soils resembles that of a “layered sponge,” in which vertical flow paths dominate.     

Groundwater dating with 3H and SF6 in relation to mixing patterns,transport modelling and hydrochemistry

A comprehensive study of a sandy aquifer of deltaic origin in southern Poland included water chemistry, isotopes, dissolved trace gases and transport modelling. Tritium, sulphur hexafluoride (SF₆) and freons (F‐11, F‐12 and F‐113) showed the presence of modern waters in the recharge areas and shallow confined parts of the aquifer. The presence of older Holocene waters farther from the recharge areas was indicated by lack of ³H, SF₆ contents ≤0·02 fmol l⁻¹ and relatively low ¹⁴C values. The discharge from the system is by upward seepage in the valley of a major river. Pre‐Holocene waters of a cooler climate, identified on the basis of δ¹⁸O, δ²H, ¹⁴C, Ne and Ar data, were found in some distant wells. Concentrations of N₂, Ne and Ar determined by gas chromatography were used for calculating the noble gas temperatures, air excess needed for correction of SF₆, and nitrogen content released by denitrification process. The time series of ³H content available for some wells supplied quantitative information on age distributions and the total mean ages of flow through the unsaturated and saturated zones. The derived ³H age distributions turned out to be very wide, with mean values in the range of about 30 to 160 years. For each well with determined ³H age, the SF₆ data showed either a lower age range or the possibility of a lower age as expected due to shorter travel times of SF₆ through the unsaturated zone, which most probably also resulted in different types of age distributions of these tracers. Freons appeared to be of little use for individual age determinations. A quantitative estimation of two‐component mixing from SF₆ ³H relations is not possible unless the travel time of ³H through the unsaturated zone is comparable to that of SF₆. The ratio of integrals of the response function over the age range with tracer and the whole response function yields the ratio of water with tracer to total flow of water. That ratio is a tracer‐dependent function of time. Transport modelling of SF₆ tracer done with MT3D code yielded initially large discrepancies between calculated and measured tracer concentrations. Some discrepancies remained even after calibration of the transport model with SF₆. Simulation of tritium contents with a calibrated transport model yielded reasonable agreement with measured contents in some wells and indicated a need for further investigations, particularly in the eastern part of the aquifer. The existence of distinct hydrochemical zones is consistent with the tracer data; young waters with measurable ³H and SF₆ contents are aerobic and of HCO₃ Ca or HCO₃ SO₄ Ca types. Slightly elevated Na and Cl contents, as well as the highest concentrations of SO₄ and NO₃ within this zone are due to anthropogenic influences. Anaerobic conditions prevail in the far field, under the confining cover, where pre‐bomb era Holocene waters dominate. In that zone, dissolved oxygen, NO₃ and U contents are reduced, and Fe, Mn and NH₄ contents increase. In the third zone, early Holocene and glacial waters occur. They are of HCO₃ Ca Na or HCO₃ Na types, with TDS values higher than 1 g l⁻¹ and Na content higher than 200 mg l⁻¹, due to either small admixtures of ascending or diffusing older water or freshening of marine sediments, a process that is probably occurring till the present time. Copyright © 2005 John Wiley & Sons, Ltd.     

Runoff‐generating processes in hydromorphic soils on a plot scale: free gravity‐driven versus pressure‐controlled flow

Draining soil water is an important runoff generator. This study aims to describe runoff‐generating processes on a plot scale (1 m²) in hydromorphic soils with different initial soil water contents. We irrigated 16 hydromorphic soils in the northern Pre‐Alps in Switzerland and recorded the variations in water content with time domain reflectometry (TDR) at five different depths per plot. Sprinkling was repeated three times at approximately 23‐h intervals and lasted for 1 h with a volume flux density of 70 mm h^?1. The comparison between the measured water content of the drainages with two physically based models revealed which of the flow processes dominated during water recessions. We distinguished between vertical drainage, lateral outflow and infiltration without drainage. Approximately 45% of all recorded time series of soil water content did not drain within approximately 20 h after the end of irrigation, about 25% drained laterally and 10% of the outflow was vertical. The drainage of the remaining 20% was the result of both lateral and vertical water flow (≈12%), or was not interpretable with the approaches applied (≈8%). Vertical flow was only observed in layers without any or with just a few hydromorphic features. Lateral draining horizons had approximately half the storage capacity and amplitude of water recession of those with vertical flow. Vertical flow was only observed in the upper soil horizons. Thus, vertical flow transmitted water to layers with lateral outflow and did not delay runoff by deep percolation. Increasing initial soil moisture increased the percentage of water content recordings according to a lateral outflow slightly, while vertical flow was less frequent. Copyright © 2010 John Wiley & Sons, Ltd.     

FLOW SEPARATION IN UNDISTURBED SOIL USING MULTIPLE ANIONIC TRACERS. PART 3: UNSTEADY CORE-SCALE INFILTRATION EXPERIMENTS

Solute transport through structured, undisturbed soil has been studied in transient, unsaturated experiments using columns from grass and woodland sites on the Lancaster University campus. Three anionic tracers have been used, bromide (Br⁻) and two fluorinated organic acids (pentraflurobenzoic acid and 2,6-diflurobenzoic acid). The process of displacement of stored water from undisturbed columns was investigated using successive inputs of different tracers under similar antecedent conditions. The results indicated that initial breakthrough was rapid, with a relative concentration of 0.8 being reached between 0.4 and 0.5 pore volumes of discharge. It was found that there was an apparent continued discharge of ‘old’ water, stored in the column before any additions of tracer, even after the addition of a total of 4.9 and 5.4 pore volumes of water for the grass and woodland columns, respectively. The implications of the results of these tracer studies for modelling solute transport in structured soils are considered.     

Hydrogeochemical balance of forest umbrisol profiles (‘Sierra de Gata’, central western Spain)

Three techniques for obtaining soil water solutions (gravitational and matrical waters extracted using both in situ tension lysimeters and in vitro pressure chambers) and their later chemical analysis were performed in order to know the evolution of the soil‐solution composition when water moves down through the soil, from the Ah soil horizon to the BwC‐ or C‐horizons of forest soils located in western Spain. Additionally, ion concentrations and water volumes of input waters to soil (canopy washout) and exported waters (drainage solutions from C‐horizons) were determined to establish the net balance of solutes in order to determine the rates of leaching or retention of ions. A generalized process of sorption or retention of most components (even Cl^?) was observed, from the soil surface to the C‐horizon, in both gravitational and matrical waters, with H₄SiO₄, Mn²⁺, Na⁺, and SO₄^2? being the net exported components from the soil through the groundwater. These results enhance the role of the recycling effect in these forest soils. The net percentages of elements retained in these forest soils, considering the inputs and the outputs balance, were 68% K⁺, 85% Ca²⁺, 58% Mg²⁺, 7% Al³⁺, 5% Fe³⁺, 34% Zn²⁺, 57% Cl^?, and 20% NO₃^?, and about 75% of dissolved organic carbon was mineralized. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

PRINCIPES NOUVEAUX POUR DISTINGUER CE QUI EST Dû A L'ÉCOULEMENT SOUTERRAIN SUR L'HYDROGRAPHE DES FLEUVES

Abstract

Throughflow has been measured from three soil horizons on a 12 slope with impermeable, bedrock. Storm flow comes from the 10–45 cm horizon and is controlled by the upslope extent of saturated conditions. Base flow comes from the 45–75 cm horizon and is supplied by slow unsaturated flow from the whole soil mass to a small constant zone of saturation.

Differences between input and output stream hydrographs over 270 metres of channel are attributed to throughflow and correlate well with measured values providing a basis for separating throughflow components from the stream hydrograph. Observed stream flows contain no overland flow or ground water flow components. The main basin flood peak is not generated within this control section of channel but is produced in the headwater zone (0.1 km²) by the faster runoff characteristics of the soils in that area and by topographic factors which lead to rapid channel extension.