Water Table Elevation Research Articles

Curve Number Derivation for Watersheds Draining Two Headwater Streams in Lower Coastal Plain South Carolina, USA

AbstractThe objective of this study was to assess curve number (CN) values derived for two forested headwater catchments in the Lower Coastal Plain (LCP) of South Carolina using a three‐year period of storm event rainfall and runoff data in comparison with results obtained from CN method calculations. Derived CNs from rainfall/runoff pairs ranged from 46 to 90 for the Upper Debidue Creek (UDC) watershed and from 42 to 89 for the Watershed 80 (WS80). However, runoff generation from storm events was strongly related to water table elevation, where seasonally variable evapotranspirative wet and dry moisture conditions persist. Seasonal water table fluctuation is independent of, but can be compounded by, wet conditions that occur as a result of prior storm events, further complicating flow prediction. Runoff predictions for LCP first‐order watersheds do not compare closely to measured flow under the average moisture condition normally associated with the CN method. In this study, however, results show improvement in flow predictions using CNs adjusted for antecedent runoff conditions and based on water table position. These results indicate that adaptations of CN model parameters are required for reliable flow predictions for these LCP catchments with shallow water tables. Low gradient topography and shallow water table characteristics of LCP watersheds allow for unique hydrologic conditions that must be assessed and managed differently than higher gradient watersheds.

Accounting for exhaustive secondary data into the mapping of water table elevation

This paper presents the incorporation of a digital elevation model into the spatial prediction of water table elevation in Mazandaran province (Iran) using a range of interpolation techniques. The multivariate methods used are: linear regression (LR), cokriging (COK), kriging with an external drift (KED) and regression kriging (RK). The analysis is performed on 3 years (1987, 1997 and 2007) of water table elevation data from about 260 monitoring wells. Prediction performances of the different algorithms are compared with two univariate techniques, i.e. inverse distance weighting and ordinary kriging (OK), through cross validation and examination of the consistency of the generated maps with the natural phenomena. Significantly smaller prediction errors are obtained for four multivariate algorithms but, in particular, KED and RK outperform LR and COK for 3 years. The results show the potential for using elevation for a more precise mapping of water table elevation.

Impact of sea-level rise on saltwater intrusion length into the coastal aquifer, Partido de La Costa, Argentina

The impact to water resources of a potential 1-m rise in sea level against the low-lying coast of Partido de La Costa, Argentina was modeled using two scenarios. The first scenario was calculated under the assumption of a constant lateral flux of freshwater. A constant water-table elevation was assumed in the second scenario. Maintaining the lateral flux of freshwater from the land (the first scenario) resulted in an approximately linear increase of the inland extent of saltwater intrusion with rising sea level; saltwater penetrated landward between 25 and 40m. Meanwhile holding the water-table elevation constant (the second scenario), caused the movement of the saltwater interface to be non-linear. In this case, landward migration in excess of 200m or more might be expected. The second scenario is more likely to be the situation in Partido de La Costa. The variation of hydrogeological parameters from north to south along the barrier conspire to make the southern reaches, where both the hydraulic conductivity and aquifer thickness are greater, more sensitive to saltwater intrusion from sea-level rise than the northern part of the barrier. These findings may be applicable to similar sandy coastal aquifers in other parts of the global coastline.

Measurement and modeling of moisture content above an oscillating water table: implications for beach surface moisture dynamics

ABSTRACTThis study examined the influence of tidally‐induced oscillations of the beach water table in regulating beach surface moisture dynamics. A series of laboratory experiments were conducted to assess the influence of hysteresis and transient flow effects on surface moisture variability. The experimental apparatus utilized a column of well‐sorted fine sand partially immersed in a reservoir of water. The water level in the reservoir was raised and lowered via a diaphragm‐metering pump to simulate tidally induced fluctuations of the water table, and the moisture content profile within the column was monitored using an array of Delta‐T probes. Moisture contents at specific elevations within the column were utilized as proxies to represent various ‘surface’ elevations (relative to the high water table). Results indicate that surface moisture content behaves in a distinctly hysteretic manner. Examination of water flow scanning curves illustrated that for all surface elevations considered, higher moisture contents for a given pressure head occurred during the drying cycle than during the wetting cycle. This observation is particularly evident with shallow surface elevations (i.e. water table close to the surface) where the Haines Jump phenomenon was found to have a significant influence on moisture content dynamics. Additionally, an assessment of the accuracy of hysteretic and non‐hysteretic models to predict the measured moisture contents demonstrated that hysteretic simulations consistently provide a better representation of the observed moisture contents than non‐hysteretic simulations. A time lag was found between the respective maxima and minima in water table elevation surface moisture content. At the near surface water table positions the time lag ranged between 30 and 100 minutes, and it increased to 240 minutes (four hours) with the high water table at 60 cm below the surface. Copyright © 2013 John Wiley & Sons, Ltd.

Dynamic coastal dune spit: the impact of morphological change on dune slacks at Whiteford Burrows, South Wales, UK

Whiteford Burrows is a coastal dune spit wetland in South Wales that is susceptible to morphological change. The height of the ridge of groundwater within the sand aquifer is essentially proportional to the width of the spit. The water table elevation impacts both the frequency and duration of slack flooding events and, therefore, slack ecology. A severe late winter storm event on 17 March 1995 caused extensive erosion of the foreshore, reducing the effective width of the dune system by 4 % and the water table elevation by up to 1 m. This observed relationship allows water level elevations in the dune system to be hindcast using historical maps and air photos which record past change in dune morphology. These historical data indicate that the dunes were relatively broad in the nineteenth century and the slacks were humid and liable to regular winter flooding. The system slowly dried out towards the 1940s as the spit thinned, when subsequent widening allowed the water table to rise and once again flood slack floors in winter. Despite these changes, the alkalinity of the Whiteford Burrows dune system has inhibited organic matter accumulation and maintained conditions needed for the persistence of a diverse basiphilous vegetation assemblage in many of the slacks.

Tree Island Response to Fire and Flooding in the Short-Hydroperiod Marl Prairie Grasslands of the Florida Everglades, USA

Within the marl prairie grasslands of the Florida Everglades, USA, the combined effects of fire and flooding usually lead to very significant changes in tree island structure and composition. Depending on fire severity and post-fire hydroperiod, these effects vary spatially and temporally throughout the landscape, creating a patchy post-fire mosaic of tree islands with different successional states. Through the use of the Normalized Difference Vegetation Index (NDVI) and three predictor variables (marsh water table elevation at the time of fire, post-fire hydroperiod, and tree island size), along with logistic regression analysis, we examined the probability of tree island burning and recovering following the Mustang Corner Fire (May to June 2008) in Everglades National Park. Our data show that hydrologic conditions during and after fire, which are under varying degrees of management control, can lead to tree island contraction or loss. More specifically, the elevation of the marsh water table at the time of the fire appears to be the most important parameter determining the severity of fire in marl prairie tree islands. Furthermore, in the post-fire recovery phase, both tree island size and hydroperiod during the first year after the fire played important roles in determining the probability of tree island recovery, contraction, or loss.

Faulting and groundwater in a desert environment: constraining hydrogeology using time‐domain electromagnetic data

ABSTRACTWithin the south‐western Mojave Desert, the Joshua Basin Water District is considering applying imported water into infiltration ponds in the Joshua Tree groundwater sub‐basin in an attempt to artificially recharge the underlying aquifer. Scarce subsurface hydrogeological data are available near the proposed recharge site; therefore, time‐domain electromagnetic (TDEM) data were collected and analysed to characterize the subsurface. TDEM soundings were acquired to estimate the depth to water on either side of the Pinto Mountain Fault, a major east‐west trending strike‐slip fault that transects the proposed recharge site. While TDEM is a standard technique for groundwater investigations, special care must be taken when acquiring and interpreting TDEM data in a two‐dimensional (2D) faulted environment. A subset of the TDEM data consistent with a layered‐earth interpretation was identified through a combination of three‐dimensional (3D) forward modelling and diffusion time‐distance estimates. Inverse modelling indicates an offset in water table elevation of nearly 40 m across the fault. These findings imply that the fault acts as a low‐permeability barrier to groundwater flow in the vicinity of the proposed recharge site. Existing production wells on the south side of the fault, together with a thick unsaturated zone and permeable near‐surface deposits, suggest the southern half of the study area is suitable for artificial recharge. These results illustrate the effectiveness of targeted TDEM in support of hydrological studies in a heavily faulted desert environment where data are scarce and the cost of obtaining these data by conventional drilling techniques is prohibitive.

An investigation of monitorability issues for groundwater in the Zachs Knob syncline area, Northeast Tennessee, USA

A dye trace was conducted in 2007 to verify a site conceptual model for groundwater monitorability near an orphan landfill in Sullivan County, Tennessee. The old landfill is in the Valley and Ridge physiographic province. Multiple tracing efforts in this area dating back to 1980 were either unsuccessful or yielded ambiguous results, causing the site to be deemed un-monitorable. Available geologic information, water chemistry data, and subsurface investigation results were compiled to produce a site conceptual model for about 120 ha of land, including the old landfill area. Eosine dye was used as a mimic for potential contamination effects in order to document the monitorability of the local aquifer as part of a hydrogeologic report for a new disposal facility permit. The site exhibits characteristics of a youthful karst setting at the surface. Overflow conduits left behind from a geologically older hydrologic system were found after detailed investigation. These remnant karst features were adjusted to a higher baseflow regime and act as a constraint to the maximum water table elevations in the present hydrologic system. Eosine dye was chosen as the best surrogate leachate after background water and bench scale leachate/dye interaction tests were done. The eosine injection was conducted in a variety of ways across the site to mimic potential leachate release scenarios from the proposed liner system for a new solid waste disposal facility. Eosine was visually detected in local springs and positively detected in some domestic supply and monitoring wells. Subsurface drainage patterns indicated two main groundwater basins in the area of study. Follow-up dye traces with simultaneous injections of pyranine and sulforhodamine B were conducted in 2009 to comply with the request of State regulators. The inferred drainage patterns were confirmed with dye detections at proposed monitoring points for each of the two main groundwater basins.

Discontinuous Steady‐State Analytical Solutions of the Boussinesq Equation and Their Numerical Representation by Modflow

Known analytical solutions of groundwater flow equations are routinely used for verification of computer codes. However, these analytical solutions (e.g., the Dupuit solution for the steady-state unconfined unidirectional flow in a uniform aquifer with a flat bottom) represent smooth and continuous water table configurations, simulating which does not pose any significant problems for the numerical groundwater flow models, like MODFLOW. One of the most challenging numerical cases for MODFLOW arises from drying-rewetting problems often associated with abrupt changes in the elevations of impervious base of a thin unconfined aquifer. Numerical solutions of groundwater flow equations cannot be rigorously verified for such cases due to the lack of corresponding exact analytical solutions. Analytical solutions of the steady-state Boussinesq equation, associated with the discontinuous water table configurations over a stairway impervious base, are presented in this article. Conditions resulting in such configurations are analyzed and discussed. These solutions appear to be well suited for testing and verification of computer codes. Numerical solutions, obtained by the latest versions of MODFLOW (MODFLOW-2005 and MODFLOW-NWT), are compared with the presented discontinuous analytical solutions. It is shown that standard MODFLOW-2005 code (as well as MODFLOW-2000 and older versions) has significant convergence problems simulating such cases. The problems manifest themselves either in a total convergence failure or erroneous results. Alternatively, MODFLOW-NWT, providing a good match to the presented discontinuous analytical solutions, appears to be a more reliable and appropriate code for simulating abrupt changes in water table elevations.

Consistency of Hydrologic Relationships of a Paired Watershed Approach

Paired watershed studies are used around the world to evaluate and quantify effects of forest and water management practices on hydrology and water quality. The basic concept uses two neighboring watersheds (one as a control and another as a treatment), which are concurrently monitored during calibration (pre-treatment) and post-treatment periods. A statistically significant relationship between the control and treatment watersheds is established during calibration period such that any significant shift detected in the relationship during treatment is attributed to the treatment effects. The approach assumes that there is a consistent, quantifiable, and predictable relationship between watershed response variables. This study tests the hypothesis that the hydrologic relationships between control and treatment watersheds for daily water table elevation (WTE) and daily flow data were similar without any statistically significant difference during two different calibration (1988-1989 and 2007-2008) and treatment periods (1995-1996 and 2009), when the control and treatment watersheds were interchanged. The watersheds are two artificially drained loblolly pine forests (D1: 24.7 haand D2: 23.6 ha) located in coastal North Carolina. Results depicted significantly similar WTE regression relationships during the two calibration periods but significantly different WTE relationships during the two treatment periods with reversed control and treatment watersheds. Calibration and treatment flow relationships, and the mean treatment effects on WTE and flow, before and after treatment reversal were significantly different (α = 0.05). The study also discusses causes of differences in hydrologic relationships and treatment effects for such reversal of treatments during a 21-year span of the study on these two similar and adjacent watersheds. The observed differences in the hydrologic relationships between control and treatment watersheds before and after treatment reversal may be attributed to climate or hydrologic non-stationarity which may affect the reliability of paired watershed approach especially when the calibration periods are short.

Modelling the interaction of surface and subsurface water flow by linking TOPMODEL and MODFLOW

A coupled ground water–surface water model was used to quantify the hydrological processes in a lowland watershed. The model was constructed by coupling TOPMODEL and MODFLOW without writing new code for the integrated model. Meteorological, land use and soil type data for the Tenmile watershed, USA, were used to run TOPMODEL, which partitions precipitation into evapotranspiration, runoff and recharge. The ground water model, MODFLOW, is driven by recharge derived from TOPMODEL and provides water table elevation and baseflow. The coupled model was calibrated and validated with stream discharge and ground water head data for the period 1968 to 1972. Nash-Sutcliffe model efficiency (E) and correlation coefficient (r) were used to evaluate model behaviour. The model demonstrated good performance, with average E and r values greater than 0.7 for stream flow and ground water head during calibration and validation. Despite obvious challenges in harmonising the different spatial and temporal discretisation in surface water and ground water models, the results from this study demonstrate that the coupled model is capable of simulating hydrological processes in the study area and that with further improvements, it can be a useful tool for integrated management of surface water and ground water resources.

Winter Hydrology and NO3− Concentrations in a Forested Watershed: A Detailed Field Study in the Adirondack Mountains of New York1

Abstract: More than 85% of NO3− losses from watersheds in the northeastern United States are exported during winter months (October 1 to May 30). Interannual variability in NO3− loads to individual streams is closely related to interannual climatic variations, particularly during the winter. The objective of our study was to understand how climatic and hydrogeological factors influence NO3− dynamics in small watersheds during the winter. Physical parameters including snow depth, soil temperature, stream discharge, and water table elevation were monitored during the 2007‐2008 winter in two small catchments in the Adirondack Mountains, New York State. Snowpack persisted from mid‐December to mid‐April, insulating soils such that only two isolated instances of soil frost were observed during the study period. NO3− export during a mid‐winter rain‐on‐snowmelt event comprised between 8 and 16% of the total stream NO3− load for the four‐month winter study period. This can be compared with the NO3− exported during the final spring melt, which comprised between 38 and 45% of the total four‐month winter NO3− load. Our findings indicate that minor melt events were detectable with changes in soil temperature, streamflow, groundwater level, and snow depth. But, based on loading, these events were relatively minor contributors to winter NO3− loss. A warmer climate and fluctuating snowpack may result in more major mid‐winter melt events and greater NO3− export to surface waters.

Nitrate removal in a restored riparian groundwater system: functioning and importance of individual riparian zones

Abstract. For the design and the assessment of river restoration projects, it is important to know to what extent the elimination of reactive nitrogen (N) can be improved in the riparian groundwater. We investigated the effectiveness of different riparian zones, characterized by a riparian vegetation succession, for nitrate (NO3−) removal from infiltrating river water in a restored and a still channelized section of the river Thur, Switzerland. Functional genes of denitrification (nirS and nosZ) were relatively abundant in groundwater from willow bush and mixed forest dominated zones, where oxygen concentrations remained low compared to the main channel and other riparian zones. After flood events, a substantial decline in NO3− concentration (> 50%) was observed in the willow bush zone but not in the other riparian zones closer to the river. In addition, the characteristic enrichment of 15N and 18O in the residual NO3− pool (by up to 22‰ for δ15N and up to 12‰ for δ18O) provides qualitative evidence that the willow bush and forest zones were sites of active denitrification and, to a lesser extent, NO3− removal by plant uptake. Particularly in the willow bush zone during a period of water table elevation after a flooding event, substantial input of organic carbon into the groundwater occurred, thereby fostering post-flood denitrification activity that reduced NO3− concentration with a rate of ~21 μmol N l−1 d−1. Nitrogen removal in the forest zone was not sensitive to flood pulses, and overall NO3− removal rates were lower (~6 μmol l−1 d−1). Hence, discharge-modulated vegetation–soil–groundwater coupling was found to be a key driver for riparian NO3− removal. We estimated that, despite higher rates in the fairly constrained willow bush hot spot, total NO3− removal from the groundwater is lower than in the extended forest area. Overall, the aquifer in the restored section was more effective and removed ~20% more NO3− than the channelized section.

Hydraulic subsurface measurements and hydrodynamic modelling as indicators for groundwater flow systems in the Rotondo granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.

Changes on groundwater flow and hydrochemistry of the Gran Sasso carbonate aquifer after 2009 L’Aquila earthquake

The earthquake that struck L'Aquila on April 6 2009 (Mw 6.3) directly affected the Gran Sasso aquifer. Co-seismic and post-seismic changes in groundwater discharge and in hydrochemistry, possibly induced by the earthquake, were observed. Spot and monitoring measurements of the spring discharge, of water table level and of the main physico-chemical parameters of spring waters (T, pH, electrical conductivity, major ions and 222Rn) were thus carried out to determine the effects of the L'Aquila earthquake on groundwater at regional and local scale, to be compared with available data collected since the 1990s.Short- and mid-term effects have been observed in the groundwater flow at recharge and discharge areas. The following short term effects have been observed: i) the sudden disappearance of some springs located along the surface trace of the Paganica Fault; ii) an immediate discharge increase of the Gran Sasso highway tunnel drainages (+20%) and of other springs (+10%); iii) a progressive increase of the water table elevation (+1m) at the boundary of theGran Sasso aquifer during the following month; iv) a sudden lowering of the water table in the recharge area. Similar post-seismic effects have been recorded in the following 20 months, when spring discharge and water table remain higher than the pre-seismic ones in discharge zones.A conceptual model of the earthquake consequences on the Gran Sasso aquifer is proposed herein. The short-term hydrologic effects registered immediately after the mainshock have been caused by a pore pressure increase related to elasto-static aquifer deformation.Apart from the contribution of seasonal recharge observed in 2009-10, mid-term effects observed in the 20 months following the mainshock suggest that there was a change in groundwater hydrodynamics.Supplementary groundwater that flows toward aquifer boundaries and springs in discharge areas reflects a possible increase in hydraulic conductivity in the recharge area. This increase is probably related to fracture cleaning and/or dilatancy.Additional monitoring including hydrochemical data allows a refinement of the proposed model. The outcomes of the hydrochemical spot sampling of the pre-seismic (2001-2007), post-seismic (April 2009) and after-seismic (July and September 2009, may 2010) periods, give the following insights: i) post-seismic groundwater of Tempera spring group was more mineralised and richer in 222Rn than the pre-seismic one; ii) transient changes in pH and calcite saturation index involve the whole aquifer, moving from Tempera springs and spreading from the recharge to discharge areas, causing changes in groundwater hydrochemistry; iii) post-seismic gradual return to previous hydrochemical equilibrium.Both hydrodynamic and hydrochemical observation converge towards a non-permanent increase of the bulk hydraulic conductivity in the aquifer portion close to the Paganica Fault (recharge area and local discharge zone). This fact has caused a lowering of the water table and of the calcite saturation index in recharge areas and simultaneously an increase of water table and flow rate in discharge zones.Complete interpretation of both quantitative and hydrochemical data allows to determine the long-term consequences of this earthquake on the groundwater flow of the Gran Sasso carbonate massif.

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.

Hydraulic and thermal response of groundwater–surface water exchange to flooding in an experimental aquifer

The effects of flooding on interactions between streams and shallow aquifers are poorly studied due to the unpredictable nature of flooding events. This study examines groundwater–surface water interactions in a constructed experimental meander over several simulated flood events using hydrologic and thermal monitoring. Detailed steady-state hydraulic head measurements were taken of the stream surface and in 59 piezometer nests during baseflow, bankfull flow, and overbank flooding conditions, allowing for three-dimensional hydrologic characterization. Additionally, a lithium tracer test was conducted, and pressure transducers and thermistors were deployed to capture transient behavior. Results demonstrate extensive coupling between the stream and adjacent alluvial aquifer under all discharge conditions. Water table elevation responds rapidly to changes in stream stage and re-equilibrates to order of magnitude increases in discharge within about 1 h. Though flooding elevates the water table, steady-state hydraulic gradients within the meander are independent of stream stage. This results from the flow boundary condition imposed by the stream, which deepens in response to flooding but has a slope that is essentially independent of discharge. In addition to the stream boundary, flow within the meander is also controlled by loss of stream water to the subsurface, which directs hydraulic gradients towards the base of the meander. Finally, the temperature distribution within the meander during normal and bankfull flow mimics the water table distribution and results from advected warm stream water being progressively cooled by vertical conduction in the direction of groundwater flow. These findings suggest that areas of high reaction rates within meanders (hot spots) will be sensitive to whether a stream reach is gaining or losing water to the subsurface. Further, the location of these hot spots of is likely independent of stream stage, though other controls on reaction rates may be affected.

A geostatistical approach for groundwater head monitoring network optimisation: case of the Sfax superficial aquifer (Tunisia)

AbstractIn many developing countries, the groundwater monitoring network is randomly designed, and consequently needs to be revised and optimised to reduce operation time and cost, remove redundant piezometric data, and strengthen sparseness data zone with supplementary observation well. The geostatistical approach used in this work is based on the universal kriging variance combined with cross‐validation test. Hence, a rational interpolation of water table elevations was performed to evaluate the accuracy of data employed in the piezometric head modelling. In this study, the groundwater monitoring network of the Sfax superficial aquifer in the south east of Tunisia was optimised. The elimination of five observation wells induces the invariability in the variance of estimate due to their less contribution in interpretations of the groundwater level behaviour. However, 38 wells were added at areas of high variance of the kriging in order to ameliorate the spatial coverage of the monitoring network.

Seepage-Induced Streambank Erosion and Instability: In Situ Constant-Head Experiments

AbstractThe effects of seepage on streambank erosion and failure are less understood compared to fluvial processes, especially the linkage between surface water and groundwater mechanisms. Field data are needed to validate laboratory seepage erosion and instability conclusions and to understand how engineering tools and methods may be applied to field conditions. An innovative trench injection system was utilized to provide a constant head on a near-streambank groundwater system when filled with stream water. This research was performed on a streambank of Dry Creek, a deeply incised stream with near-vertical banks located in Mississippi. Experiments included installing a trench (2.8 m from the bank and 2 m below ground surface) and a network of tensiometers and observation wells to measure soil-water pressures and water table elevations. Bank stratigraphy consisted of a sloping, conductive loamy sand layer between cohesive streambank layers. Groundwater conditions were monitored during a series of induced...

Boussinesq equation-based model for flow in drainage layer of highway

In order to study the water flow in the drainage layer of highway under steady-state condition, one-dimensional (1D) Boussinesq equation-based model with Dupuit-Forchheimer assumption was established and the semi-analytical solutions to predict the water-table height were presented. In order to validate the model, a two-dimensional (2D) saturated flow model based on the Laplace equation was applied for the purpose of the model comparison. The water-table elevations predicted by 1D Boussinesq equation-based model and 2D Laplace equation-based model match each other well, which indicates that the horizontal flow in drainage layer is dominated. Also, it validates the 1D Boussinesq equation-based model which can be applied to predict the water-table elevation in drainage layer. Further, the analysis was conducted to examine the effect of infiltration rate, hydraulic conductivity and slope of drainage layer on the water-table elevation. The results show that water-table falls down as the ratio of Is to K decreases and the slope increases. If the aquifer becomes confined by the top of drainage layer due to the larger ratio of Is to K or smaller slope, the solution presented in this work can also be applied to approximate the water-table elevation in unconfined sub-section as well as hydraulic head in the confined sub-section.