Estimating the role of a ditch network in groundwater recharge in a Mediterranean catchment using a water balance approach

Water balance variables were monitored in a farmed Mediterranean catchment characterized by a dense ditch network to allow for the separate estimation of the diffuse and concentrated recharge terms during flood events. The 27 ha central part of the catchment was equipped with (i) rain gauges, (ii) ditch gauge stations, (iii) piezometers, (iv) neutron probes, and (v) an eddy covariance mast including a 3D sonic anemometer and a fast hygrometer. The water balance was calculated for two autumnal rain and flood events. We also estimated the uncertainty of this approach with Monte Carlo simulations. Results show, that although ditch area represents only 6% of the total study area, concentrated recharge appeared to be the main source of groundwater recharge. Indeed, it was 40–50% of the total groundwater recharge for autumnal events, which are the major annual recharge events. This indicate that both, concentrated and diffuse recharge should be taken into account in any hydrological modeling approach for Mediterranean catchments. This also means that, since they collect overland flow that is often largely contaminated by chemicals, ditches may be a place where groundwater contamination is likely to occur. The uncertainty analysis indicates that recharge estimates based on water balance exhibit large uncertainty ranges. Nevertheless, Monte Carlo simulations showed that concentrated recharge was higher than expected based on their area.     

The impact of ditch blocking on the hydrological functioning of blanket peatlands

Ditch blocking in blanket peatlands is common as part of peatland restoration. The effects of ditch blocking on flow regimes and nearby water tables were examined in a field trial. After an initial 6‐month monitoring period, eight ditches had peat dams installed 10 m apart along their entire length (dammed), four of these ditches were also partially infilled through bank reprofiling (reprofiled). Four ditches were left open with no dams or reprofiling (open). These 12 ditches and the surrounding peat were monitored for 4 more years. An initial five‐fold reduction in discharge occurred in the dammed and the reprofiled ditches with the displaced water being diverted to overland flow and pathways away from the ditches. However, there was a gradual change over time in ditch flow regime in subsequent years, with the overall volume of water leaving the dammed and the reprofiled ditches increasing per unit of rainfall to around twice that which occurred in the first year after blocking. Hence, monitoring for greater than one year is important for understanding hydrological impacts of peatland restoration. Overland flow and flow in the upper ~4 cm of peat was common and occurred in the inter‐ditch areas for over half of the time after ditch blocking. There was strong evidence that topographic boundaries of small ditch catchments, despite being defined using a high‐resolution Light Detection And Ranging‐based terrain model, were not always equivalent to actual catchment areas. Hence, caution is needed when upscaling area‐based fluxes, such as aquatic carbon fluxes, from smaller scale studies including those using ditches and small streams. The effect of ditch blocking on local water tables was spatially highly variable but small overall (time‐weighted mean effect <2 cm). Practitioners seeking to raise water tables through peatland restoration should first be informed either by prior measurement of water tables or by spatial modelling to show whether the peatland already has shallow water tables or whether there are locations that could potentially undergo large water‐table recoveries.     

A simple overland flow calculation method for distributed groundwater recharge models

A method to improve the calculation of overland flow in distributed groundwater recharge models is presented and applied to two sub‐catchments in the Thames Basin, UK. Recharge calculation studies tend to simulate the runoff flow component of river flow in a simplistic way, often as a fraction of rainfall over a particular period. The method outlined in this study intends to improve the calculation of groundwater recharge estimates in distributed recharge models but does not present an alternative to complex overland flow simulators. This method uses seasonally varying coefficients to calculate runoff for specified hydrogeological classes or runoff zones, which are used to model baseflow index variations across the basin. It employs a transfer function model to represent catchment storage. Monte Carlo simulation was applied to refine the runoff values. Decoupling the runoff zones between the two sub‐catchments produces a better match between the simulated and observed values; however, the difference between observed runoff and the simulated output indicates other factors, such as landuse and topographical characteristics that affect the generation of runoff flow, need to be taken into account when classifying runoff zones. British Geological Survey © NERC 2011. Hydrological Processes © 2011 John Wiley & Sons, Ltd.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Effect of bedrock flow on catchment rainfall‐runoff characteristics and the water balance in forested catchments in Tanzawa Mountains,Japan

In this paper, we examined the role of bedrock groundwater discharge and recharge on the water balance and runoff characteristics in forested headwater catchments. Using rigorous observations of catchment precipitation, discharge and streamwater chemistry, we quantified net bedrock flow rates and contributions to streamwater runoff and the water balance in three forested catchments (second‐order to third‐order catchments) underlain by uniform bedrock in Japan. We found that annual rainfall in 2010 was 3130 mm. In the same period, annual discharge in the three catchments varied from 1800 to 3900 mm/year. Annual net bedrock flow rates estimated by the chloride mass balance method at each catchment ranged from ?1600 to 700 mm/year. The net bedrock flow rates were substantially different in the second‐order and third‐order catchments. During baseflow, discharge from the three catchments was significantly different; conversely, peak flows during large storm events and direct runoff ratios were not significantly different. These results suggest that differences in baseflow discharge rates, which are affected by bedrock flow and intercatchment groundwater transfer, result in the differences in water balance among the catchments. This study also suggests that in these second‐order to third‐order catchments, the drainage area during baseflow varies because of differences between the bedrock drainage area and surface drainage area, but that the effective drainage area during storm flow approaches the surface drainage area. Copyright © 2012 John Wiley & Sons, Ltd.     

Contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments

We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 John Wiley & Sons, Ltd.     

Evaluating hydrologic performance of bioretention cells in shallow groundwater

Bioretention cells, which are generally effective in controlling surface runoff and recharging groundwater, have been widely adopted as low impact development practices. However, shallow groundwater has limited their implementation in some locations due to the potential problems of a reduction in surface runoff control, groundwater pollution, and continuous groundwater drainage through the underdrain. Many guidelines have established minimum requirements for the groundwater depth below bioretention cells, but they may not be optimized for certain environmental conditions and bioretention cell designs. This study made use of a variably saturated flow model to examine the hydrologic performance of a single bioretention cell in shallow groundwater with event‐based simulations, considering a wide range of initial groundwater depths, media and in situ soil types, surface runoff loads, and underdrain sizes. Performance indicators (e.g., runoff reduction, time for infiltrated water to reach the bioretention cell bottom and the groundwater table, and height and dissipation time of groundwater mound) were evaluated to examine the processes of runoff generation, the formation and dissipation of groundwater mounds, and the bioretention cell's performance in a shallow groundwater environment. The most influential factors were the initial groundwater depth, the hydraulic conductivity of the media soil, and the rainfall runoff load. With a deeper initial groundwater table, infiltrated water took longer to reach the bioretention cell bottom and groundwater table. Groundwater mounds, however, took longer to dissipate even though they were smaller. The groundwater quality can be better protected if relatively less‐permeable soil types (e.g., sandy loam) are used as the media, although it may compromise the performance in runoff quantity control. However, only very high surface runoff loads would cause concerns regarding a reduction in runoff quantity control and possible groundwater contamination due to the shallow groundwater. A distance of 1.5–3 m between the bioretention cell bottom and the groundwater table is generally sufficient. The results of this study could help to guide the planning and design of bioretention cells in areas of shallow groundwater.     

Field data and flow system response in clay (vertisol) shale terrain,north central Texas,USA

The water budget in clay shale terrain is controlled by a complex interaction between the vertisol soil layer, the underlying fractured rock, land use, topography, and seasonal trends in rainfall and evapotranspiration. Rainfall, runoff, lateral flow, soil moisture, and groundwater levels were monitored over an annual recharge cycle. Four phases of soil–aquifer response were noted over the study period: (1) dry‐season cracking of soils; (2) runoff initiation, lateral flow and aquifer recharge; (3) crack closure and down‐slope movement of subsurface water, with surface seepage; (4) a drying phase. Surface flow predominated within the watershed (25% of rainfall), but lateral flow through the soil zone continued for most of the year and contributed 11% of stream flow through surface seepage. Actual flow through the fractured shale makes up a small fraction of the water budget but does appear to influence surface seepage by its effect on valley‐bottom storage. When the valley soil storage is full, lateral flow exits onto the valley‐bottom surface as seasonal seeps. Well response varied with depth and hillslope position. FLOWTUBE model results and regional recharge estimates are consistent with an aquifer recharge of 1·6% of annual precipitation calculated from well heights and specific yield of the shale aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigating hydrological regimes and processes in a set of catchments with temporary waters in Mediterranean Europe

Abstract

Seven catchments of diverse size in Mediterranean Europe were investigated in order to understand the main aspects of their hydrological functioning. The methods included the analysis of daily and monthly precipitation, monthly potential evapotranspiration rates, flow duration curves, rainfall—runoff relationships and catchment internal data for the smaller and more instrumented catchments. The results showed that the catchments were less “dry” than initially considered. Only one of them was really semi-arid throughout the year. All the remaining catchments showed wet seasons when precipitation exceeded potential evapotrans-piration, allowing aquifer recharge, “wet” runoff generation mechanisms and relevant baseflow contribution. Nevertheless, local infiltration excess (Hortonian) overland flow was inferred during summer storms in some catchments and urban overland flow in some others. The roles of karstic groundwater, human disturbance and low winter temperatures were identified as having an important impact on the hydrological regime in some of the catchments.     

Hydrological impact of roadside ditches in an agricultural watershed in Central New York: implications for non‐point source pollutant transport

Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint‐source (NPS) pollution by acting as hydrological links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field‐ and watershed‐scale hydrology, and pollutant transport in a 38 km² agricultural watershed in south‐central NY. Salient findings include the following: (i) 94% of road ditches discharged to natural streams, effectively doubling the drainage density; (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches; and (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. II: modelling the effects on groundwater flow

Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.     

Stormflow generation involving pipe flow in a zero‐order basin of Peninsular Malaysia

Hydrological responses in a zero‐order basin (ZOB), a portion of whose discharge emerged via preferential flow through soil pipes, were examined over a 2‐year period in Peninsular Malaysia to elucidate primary stormflow generation processes. Silicon (Si) and specific conductance (EC) in various runoff components were also measured to identify their sources. ZOB flow response was dependent on antecedent precipitation amount; runoff increased linearly with precipitation during events >20 mm in relatively wet antecedent moisture conditions. Runoff derived from direct precipitation falling onto saturated areas accounted for <0·2% of total ZOB flow volume during the study period, indicating the predominance of subsurface pathways in ZOB flow. ZOB flow (high EC and low Si) was distinct from perennial baseflow via bedrock seepage (low EC and high Si) 5 m downstream of the ZOB outlet. Pipe flow responded quickly to ZOB flow rate and was characterized by a threshold flow capacity unique to each pipe. Piezometric data and pipe flow records demonstrated that pipes located deeper in the soil initiated first, followed by those at shallower depths; initiation of pipe flow corresponded to shallow groundwater rise above the saprolite‐soil interface. Chemical signatures of pipe flow were similar to each other and to the ZOB flow, suggesting that the sources were well‐mixed soil‐derived shallow groundwater. Based upon the volume of pipe flow during storms, the combined contribution of the pipes monitored accounted for 48% of total ZOB flow during the study period. Our results suggest that shallow groundwater, possibly facilitated by preferential flow accreted above the saprolite–soil interface, provides dominant stormflow, and that soil pipes play an important role in the rapid delivery of solute‐rich water to the stream system. Copyright © 2006 John Wiley & Sons, Ltd.     

Impacts of forest conversion and agriculture practices on water pathways in Southern Brazil

Land‐use/cover change (LUCC), and more specifically deforestation and multidecadal agriculture, is one of the various controlling factors of water fluxes at the hillslope or catchment scale. We investigated the impact of LUCC on water pathways and stream stormflow generation processes in a subtropical region in southern Brazil. We monitored, sampled and analysed stream water, pore water, subsurface water, and rainwater for dissolved silicon concentration (DSi) and ¹⁸O/¹⁶O (δ¹⁸O) signature to identify contributing sources to the streamflow under forest and under agriculture. Both forested and agricultural catchments were highly responsive to rainfall events in terms of discharge and shallow groundwater level. DSi versus δ¹⁸O scatter plots indicated that for both land‐use types, two run‐off components contributed to the stream discharge. The presence of a dense macropore network, combined with the presence of a compact and impeding B‐horizon, led to rapid subsurface flow in the forested catchment. In the agricultural catchment, the rapid response to rainfall was mostly due to surface run‐off. A 2‐component isotopic hydrograph separation indicated a larger contribution of rainfall water to run‐off during rainfall event in the agricultural catchments. We attributed this higher contribution to a decrease in topsoil hydraulic conductivity associated with agricultural practices. The chemical signature of the old water component in the forested catchment was very similar to that of the shallow groundwater and the pore soil water: It is therefore likely that the shallow groundwater was the main source of old water. This is not the case in the agricultural catchments where the old water component had a much higher DSi concentration than the shallow groundwater and the soil pore water. As the agricultural catchments were larger, this may to some extent simply be a scale effect. However, the higher water yields under agriculture and the high DSi concentration observed in the old water under agriculture suggest a significant contribution of deep groundwater to catchment run‐off under agriculture, suggesting that LUCC may have significant effects on weathering rates and patterns.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Hillslope hydrology and wetland response of two small zero‐order boreal catchments on the Precambrian Shield

Two Precambrian Shield zero‐order catchments were monitored from January 2003 to July 2004 to characterize their hydrological and biogeochemical characteristics prior to a forest management experiment. Hydrometric observations were used to examine temporal trends in hillslope‐wetland connectivity and the hillslope runoff processes that control wetland event response. The hillslope groundwater flux from the longer transect (E1) was continuous throughout the study period. Groundwater fluxes from a shorter and steeper hillslope (E0) were intermittent during the study period. Large depression storage elements (termed micro‐basins) located on the upper hillslope of the E1 catchment appeared to be at least partly responsible for the observed rapid wetland runoff responses. These micro‐basins were hydrologically connected to a downslope wetland by a subsurface channel of glacial cobbles that functioned as a macropore channel during episodic runoff events. The runoff response from the hilltop micro‐basins is controlled by antecedent water table position and water is quickly piped to the wetland fringe through the cobble channel during high water table conditions. During periods of low water table position, seepage along the bedrock–soil interface from the hilltop micro‐basin and other hillslopes maintained hillslope–wetland connectivity. The micro‐basins create a dynamic variable source‐area runoff system where the contributing area expands downslope during episodic runoff events. The micro‐basins occupied 30% of the E1 catchment and are a common feature on the Precambrian Shield. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table

Abstract

Many of the hydrological and ecological functions of alluvial flood plains within watersheds depend on the water flow exchanges between the vadoze soil zone and the shallow groundwater. The water balance of the soil in the flood plain is investigated, in order to evaluate the main hydrological processes that underlie the temporal dynamics of soil moisture and groundwater levels. The soil moisture and the groundwater level in the flood plain were monitored continuously for a three-year period. These data were integrated with the results derived from applying a physically-based numerical model which simulated the variably-saturated vertical water flow in the soil. The analysis indicated that the simultaneous processes of lateral groundwater flow and the vertical recharge from the unsaturated zone caused the observed water table fluctuations. The importance of these flows in determining the rises in the water table varied, depending on soil moisture and groundwater depth before precipitation. The monitoring period included two hydrological years (September 2009–September 2011). About 13% of the precipitation vertically recharged the groundwater in the first year and about 50% in the second. The difference in the two recharge coefficients was in part due to the lower groundwater levels in the recharge season of the first hydrological year, compared to those observed in the second. In the latter year, the shallow groundwater increased the soil moisture in the unsaturated zone due to capillary rise, and so the mean hydraulic conductivity of the unsaturated soil was high. This moisture state of soil favoured a more efficient conversion of infiltrated precipitation into vertical groundwater recharge. The results show that groundwater dynamics in the flood plain are an important source of temporal variability in soil moisture and vertical recharge processes, and this variability must be properly taken into account when the water balance is investigated in shallow groundwater environments.

Citation Pirastru, M. and Niedda, M., 2013. Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table. Hydrological Sciences Journal, 58 (4), 898–911.     

Nitrate fluctuations at the water table: implications for recharge processes and solute transport in the Chalk aquifer

This study investigates fluctuations in nitrate concentration at the water table to improve understanding of unsaturated zone processes in the Chalk aquifer. Sampling was conducted using a novel multi‐level sampler during periods of water table rise over 5 years at a vertical resolution of 0.05 m. Nitrate concentration increased as the water table seasonally recovered, with similar inter‐annual trends with depth. The rising water table activated horizontal fractures facilitating the delivery of water elevated by up to 10 mg/l of nitrate with respect to the adjacent groundwater below. These fractures are considered to activate via piston displacement of water from the adjoining matrix. Hydrograph analysis identified 16 events which significantly perturbed the water table within 24–48 h of rainfall. Consistent nitrate concentrations indicate recharge through persistent fracture flow from the surface was not generally the primary driver of the rapid water table response during these events. Instead, the response was attributed to the piston displacement of porewater immediately above the water table. However, a single event in November 2012 delivered relatively dilute recharge indicating rapid persistent fracture flow following rainfall was possible to a depth of 14–15 m. Decreases in porewater nitrate concentration around fracture horizons and the dilution of many groundwater samples with respect to porewaters indicate a fresher source of water at depth. This was considered most likely to be a result of near surface water bypassing the matrix because of widespread mineralization on fracture surfaces, which retard water and solute exchange. Therefore, persistent fracture flow maybe considered a frequent process, operating independently of the matrix, and is not necessarily event driven. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal variability of groundwater dynamics in a sub‐Mediterranean mountain catchment

The temporal and spatial dynamics of groundwater was investigated in a small catchment in the Spanish Pyrenees, which was extensively used for agriculture in the past. Analysis of the water table fluctuations at five locations over a 6‐year period demonstrated that the groundwater dynamics had a marked seasonal cycle involving a wetting‐up period that commenced with the first autumn rainfall events, a saturation period during winter and spring and a drying‐down period from the end of spring until the end of the summer. The length of the saturation period showed great interannual variability, which was mainly influenced by the rainfall and evapotranspiration characteristics. There was marked spatial variability in the water table, especially during the wetting‐up period, which could be related to differences in slope and drainage area, geomorphology, soil properties and local topography. Areas contributing to runoff generation were identified within the catchment by field mapping of moisture conditions. Areas contributing to infiltration excess runoff were correlated with former cultivated fields affected by severe sheetwash erosion. Areas contributing to saturation excess runoff were characterized by a marked spatial dynamics associated with catchment wetness conditions. The saturation spatial pattern, which was partially related to the topographic index, was very patchy throughout the catchment, suggesting the influence of other factors associated with past agricultural activities, including changes in local topography and soil properties. The relationship between water table levels and stream flow was weak, especially during the wetting‐up period, suggesting little connection between ground water and the hydrological response, at least at some locations. The results suggest that in drier and human‐disturbed environments, such as sub‐Mediterranean mountains, saturation patterns cannot be represented only by the general topography of the catchment. They also suggest that groundwater storage and runoff is not a succession of steady‐state flow conditions, as assumed in most hydrological models. Copyright © 2013 John Wiley & Sons, Ltd.     

Field‐Scale Relationships Among Soil Properties and Shallow Groundwater Quality

It is important to understand the link between land surface/soil properties and shallow groundwater quality. To that end, soil properties and near‐water‐table groundwater chemistry of a shallow, unconfined aquifer were measured on a 100‐m grid on a 64‐ha irrigated field in southeastern North Dakota. Soil properties and hydrochemistry were compared via multivariate analysis that included product‐moment correlations and factor analysis/principal component analysis. Topographic low areas where the water table was in close proximity to the soil surface generally had higher apparent electrical conductivity (EC_a) and higher percent silt and clay than higher positions on the landscape. The majority of the groundwater was characterized by Ca‐ and Mg‐HCO₃ type water and was associated with topographic high areas with lower EC_a and net groundwater recharge. Small topographic depressions were areas of higher EC_a (net groundwater discharge) where salts that precipitated via evapotranspiration and evaporative discharge dissolved and leached to the groundwater during short‐term depression‐focused recharge events. At this site, groundwater quality and soil EC_a were related to surface topography. High‐resolution topography and EC_a measurements are necessary to characterize the land surface/soil properties and surficial groundwater quality at the field‐scale and to delineate areas where the shallow groundwater is most susceptible to contamination.