Observed Water Table Research Articles

A Stochastic Modelling Technique for Predicting Groundwater Table Fluctuations with Time Series Analysis

This paper demonstrates a methodology for modeling the groundwater table fluctuations observed in a shallow unconfined aquifer of Bangladesh. For this purpose, Kushtia district of Bangladesh has been selected as a case study area. A groundwater table monitoring well is selected in each upazilla (sub-district) under the study area. The time series water table observations collected on a weekly basis during the period from 1999 to 2006 from each well is used for experiments. Empirical result demonstrates that observed data of groundwater level exhibits cyclic patterns and shows an annual periodicity. Box and Jenkins univariate stochastic models widely known as ARIMA, are applied to simulate groundwater table fluctuations in all monitoring wells under consideration in the study area. The predicted data using the best models are compared to the observed time series and the resulting accuracy is checked using different error parameters. The simulation revealed that ARIMA models generate reasonably accurate forecasts in terms of numerical accuracy. The results also showed that the predicted data represent the actual data very well for each monitoring well.

Evidence for the onset and persistence with depth of preferential flow in unsaturated fractured porous media

Two distinct types of fracture flow can occur in unsaturated fractured porous media: non-preferential flow, whereby the fractures and matrix wet up in equilibrium, and preferential flow, whereby water in the fractures bypasses the matrix. This has important implications for how infiltration, recharge, groundwater flooding and contaminant transport are modelled. Taking the UK Chalk as a case study, we explore evidence for the occurrence of unsaturated preferential fracture flow, considering separately its initiation in the near surface, and persistence at 20–30 m depth. We postulate a link between the apparent hysteretic response of the Chalk soil moisture characteristic, which was observed and simulated, and the initiation of preferential recharge. Focusing on observed water table responses to an extreme rainfall event on 20th July 2007, we show by inverse modelling of recharge that preferential flow persisted to a depth of 20 m below ground level, but not to 30 m, resulting in markedly different water table responses at two boreholes. These findings both lend support to our conceptualisation. However, since the observations from which preferential flow can be inferred are few and are indirect, quantification of the controls on the onset and depth persistence of preferential flow remains a significant challenge.

A combined field and modeling study of groundwater flow in a tidal marsh

Abstract. Bald mud beaches were found among the mangrove marshes in Dongzhaigang National Nature Reserve, Hainan, China. To investigate the possible reasons for this phenomenon, the intertidal zones of a mangrove transect and a bald beach transect with similar topography and tidal actions were selected for comparison study. Along both transects, observed water table variations were significant in the high and low intertidal zones and negligible in the middle intertidal zones. Despite the same tidal actions and above-mentioned similarities, observed groundwater salinity was significantly smaller along the mangrove transect (average 23.0 ppt) than along the bald beach transect (average 28.5 ppt). These observations invite one hypothesis: the hydraulic structure of tidal marsh and freshwater availability may be the main hydrogeological factors critical to mangrove development. Two-dimensional numerical simulations corroborated the speculation and gave results in line with the observed water table. The two transects investigated were found to have a mud-sand two-layered structure: a surface zone of low-permeability mud and an underlying high-permeability zone that outcrops at the high and low tide lines. The freshwater recharge from inland is considerable along the mangrove transect but negligible along the bald beach transect. The high-permeability zone may provide opportunity for the plants in the mangrove marsh to uptake freshwater and oxygen through their roots extending downward into the high-permeability zone, which may help limit the buildup of salt in the root zone caused by evapotranspiration and enhance salt removal, which may further increase the production of marsh grasses and influence their spatial distribution. The bald beach is most probably due to the lack of enough freshwater for generating a brackish beach soil condition essential to mangrove growth. It is also indicated that seawater infiltrated the high-permeability zone through its outcrop near the high intertidal zone, and discharged from the tidal river bank in the vicinity of the low tide line. These processes thereby formed a tide-induced seawater-groundwater circulation, which likely provided considerable contribution to the total submarine groundwater discharge (SGD). Finally, implications and uncertainties behind this study were summarized for future examinations.

Groundwater Flow in a Tidally Influenced Gravel Beach in Prince William Sound, Alaska

The hydraulics of a gravel beach on Eleanor Island in Prince William Sound (PWS), Alaska, was investigated by measurements of the water table and by numerical simulations using the software SUTRA. Sediment samples were obtained from six wells at the approximate three depths: 0.30, 0.60, and 0.90 m, and they were sieved in the lab to generate the grain size distribution, and subsequently the saturated hydraulic conductivity, K0. The results indicated that the beach consisted of two layers: an upper layer with a high K0 and a lower layer with a low K0. In the landward portion of the intertidal zone, K0 was around 10 � 2 m∕s and 10 � 5 m∕s for the upper and lower layers, respectively. In the seaward portion of the intertidal zone, the K0 values were 10 � 4 m∕s and 10 � 6 m∕s in the upper and lower layers, respectively. The agreement of the model with the observed water table data was very good. A decrease in the slope of the beach in the seaward direction caused water to flow out from the lower layer and propagate seaward within the upper layer. This ensured that the water table seaward of the exit point was replenished and remained above the interface of the layers. The results suggest that the porosity of the beach decreases from 0.33 at the surface to around 0.05 at 1.0 m deep. The absence of oil on this beach is probably because the beach water table remained above the layers' interface during low tides. Using the calibrated model, it was found that the flow from land into the beach is around 6 m 3 ∕day (per meter parallel to the shoreline) and that the flow out of the lower layer varied with time reaching a maximum of 13 m 3 ∕ðday·mÞ at low tide. The beach exchanges up to 7 m 3 ∕ðday·mÞ with the sea. DOI: 10.1061/(ASCE)HE.1943-5584.0000454. © 2012 American Society of Civil

Water subsidies from mountains to deserts: their role in sustaining groundwater-fed oases in a sandy landscape

In arid regions throughout the world, shallow phreatic aquifers feed natural oases of much higher productivity than would be expected solely from local rainfall. In South America, the presence of well-developed Prosopis flexuosa woodlands in the Monte Desert region east of the Andes has puzzled scientists for decades. Today these woodlands provide crucial subsistence to local populations, including descendants of the indigenous Huarpes. We explore the vulnerability and importance of phreatic groundwater for the productivity of the region, comparing the contributions of local rainfall to that of remote mountain recharge that is increasingly being diverted for irrigated agriculture before it reaches the desert. We combined deep soil coring, plant measurements, direct water-table observations, and stable-isotopic analyses (2H and 18O) of meteoric, surface, and ground waters at three study sites across the region, comparing woodland stands, bare dunes, and surrounding shrublands. The isotopic composition of phreatic groundwaters (delta2H: -137 per thousand +/- 5 per thousand) closely matched the signature of water brought to the region by the Mendoza River (-137 per thousand +/- 6 per thousand), suggestin that mountain-river infiltration rather than in situ rainfall deep drainage (-39 per thousand +/- 19 per thousand) was the dominant mechanism of recharge. Similarly, chloride mass balances determined from deep soil profiles (> 6 m) suggested very low recharge rates. Vegetation in woodland ecosystems, where significant groundwater discharge losses, likely >100 mm/yr occurred, relied on regionally derived groundwater located from 6.5 to 9.5 m underground. At these locations, daily water-table fluctuations of 10 mm, and stable-isotopic measurements of plant water, indicated groundwater uptake rates of 200-300 mm/yr. Regional scaling suggests that groundwater evapotranspiration reaches 18-42 mm/yr across the landscape, accounting for 7 17% of the Mendoza River flow regionally. Our study highlights the reliance of ecosystem productivity in natural oases on Andean snowmelt, which is increasingly being diverted to one of the largest irrigated regions of the continent. Understanding the ecohydrological coupling of mountain and desert ecosystems here and elsewhere should help managers balance production agriculture and conservation of unique woodland ecosystems and the rural communities that rely on them.

Numerical modeling of shallow water table behavior with Lisse effect

Air entrapment is an important consideration in environments with shallow water tables and sandy soil, like the condition of highly conductive sandy soils and flat topography in Florida, USA. It causes water table rises in soils, which are significantly faster and higher than those in soils without air entrapment. Two numerical models, Integrated Hydrologic Model (IHM) and HYDRUS-1D (a single-phase, one-dimensional Richards’ equation model) were tested at an area of west central Florida to help further understanding the shallow water table behavior during a long term air entrapment. This investigation employed field data with two modeling approaches to quantify the variation of air pressurization values. It was found that the air pressurization effect was responsible at time up to 40 cm of water table rise being recorded by the observation well for these two models. The values of air pressurization calculated from IHM and HYDRUS-1D match the previously published values. Results also indicated that the two numerical models did not consider air entrapment effect (as the predictive parameters remain uncertain) and thus results of depth to water table from these models did not compare to the observations for these selected periods. Incorporating air entrapment in prediction models is critical to reproduce shallow water table observations.

Hydrodynamics in a gravel beach and its impact on the Exxon Valdez oil

This paper investigated the interaction of groundwater and seawater in a tidally influenced gravel beach. Field observations of water table, pore water salinity were performed. The two‐dimensional finite element model MARUN was used to simulate observed water table and salinity. Based on field observations and model calibrations, a two‐layered beach structure was identified which is characterized by a high‐permeability surface layer underlain by a low‐permeability lower layer. The salt wedge seaward of the low tide line was almost invariant in comparison with the strong fluctuations of the salinity plume in the surface layer of the intertidal zone. The presence of the two layers prevented the presence of a freshwater discharge “tube” between the upper saline plume and salt wedge. This is in contrast with the previous works where freshwater discharge tube was observed. The tide‐induced submarine groundwater discharge (SGD) was estimated at 9 m3 d−1 m−1, a large value that is probably due to the large tidal range of ∼4.8 m and the very permeable surface layer. The freshwater‐seawater dynamics revealed here may provide new insights into the complexity, intensity, and time scales of mixing between fresh groundwater and seawater in tidal beaches. The simulated water table of the beach was higher than the interface between the surface layer and the lower layer, which prevented Exxon Valdez oil from penetrating into the lower layer in 1989.

Hydrodynamic factors affecting the persistence of the Exxon Valdez oil in a shallow bedrock beach

We report a field study and numerical modeling of multicomponent flow in a tidal gravel beach in Knight Island, Prince William Sound, Alaska, where oil from the 1989 Exxon Valdez oil spill persisted. Field measurements of water table, salinity, and tracer (lithium) concentration were obtained for around a week during the summer of 2008. The numerical model MARUN was used to simulate the field observations. On the basis of field experiments and numerical simulations, the beach was identified to have a two‐layered hydraulic structure: a high‐permeability surface layer underlain by a low‐permeability lower layer. The hydraulic conductivity was found to be 5 × 10−2 m s−1 for the surface layer and 7 × 10−6 m s−1 for the lower layer. The simulations reproduced the observed water table, salinity, and lithium concentrations accurately. The small flow entering the beach from the land side resulted in a beach water table dropping below the interface of the two layers. This seems to be the major reason for the presence of oil in the lower layer. The exchange flow between the beach and the sea due to tidal influence was ∼2.12 m3 d−1 m−1. The patterns of inflow and outflow rates showed that the maximum seawater‐groundwater exchange occurred in the middle to high intertidal zone, which explains the persistence of oil in the lower intertidal zone. To explore bioremediation of the beach with nutrient amendment, a numerical simulation of nutrient application on the beach surface was conducted, where the applied nutrient concentration was 5,000 mg L−1. The results showed that the nutrient concentration remaining in oiled areas after a week was larger than 50 mg L−1, which is larger than that needed for maximum microbial growth (2–10 mg L−1). This implies that the bioremediation via nutrient application on the beach surface could be adopted if nutrients were the only limiting factor.

Bancos elevados y tubería de aireación para la producción de yautía (Xanthosoma spp.) en un Ultisol de la zona montañosa con pobre drenaje y nivel freático temporero.

La yautía (Xanthosoma spp.) es parte importante de la dieta puertorriqueña, con alto valor en el mercado. Sin embargo, la producción local ha disminuido dramáticamente en décadas recientes, atribuible en parte a una condición de raíz conocida como "mal seco". Dicha condición se acentúa en suelos pobremente drenados, particularmente si ocurren períodos alternos de sequía y exceso de agua. En un suelo Ultisol muy pobremente drenado de la zona montañosa, se evaluaron los rendimientos de yautía bajo tres sistemas de manejo: 1) siembra convencional en suelo sin bancos ele vados; 2) siembra en bancos elevados; y 3) siembra en bancos ele vados con tubería para aireación en la base de los bancos. Todos los tratamientos incluyeron riego por goteo para evitar estrés hídrico. Se instalaron tubos perforados de PVC para observación de niveles freáticos. Los rendimientos de cormelos mercadeables fueron 5.38, 19.94 y 17.61 t/ha para los tratamientos 1, 2 y 3, respectivamente. Los rendimientos de los dos tratamientos con bancos (con o sin tubería de aireación) fuer on iguales estadísticamente (p = 0.05), pero ambos tratamientos produjeron rendimientos estadísticamente mayores que en el tratamiento sin bancos. El nivel freático promedio permaneció a menos de 50 cm desde la superficie durante la mayor parte del desarrollo del cultivo, y en varias fechas de muestreo se encontró a menos de 30 cm desde la superficie. Una inspección manual del suelo al momento de cosecha indicó que el suelo se mantuvo en condición suelta y friable aun 10 meses luego de la labranza inicial. Los resultados de este estudio indican que en suelos pobremente drenados con buena estructura (baja impedancia mecánica) se pueden obtener altos rendimientos de yautía sembrando en bancos elevados bajo riego por goteo,sin necesidad de instalar tubería de aireación. Se documenta también la importancia de medir los niveles freáticos en suelos pobremente drenados de la zona montañosa.

Modelling the Effects of Groundwater-Based Urban Supply in Low-Permeability Aquifers: Application to the Madrid Aquifer, Spain

The European Union Water Framework Directive establishes the obligation to all Member States to develop and implement catchment-scale management plans. These aim at reaching a good status of all water bodies in Member States by 2015 or, at the latest, by 2027. Numerical models provide a suitable approach to evaluate the possibilities of achieving this goal by enabling users to deal with complex hydrological interrelations in a dynamic manner. This paper presents a modelling-based approach to examine the past and future effects of groundwater-based urban supply in Metropolitan Madrid, Spain. Monthly scale model calibration is based on 32 years of water table observations. Care is taken to address the effect of intensive pumping on groundwater levels and streamflows, as well as the role of groundwater resources in meeting urban demands during dry spells. The paper concludes by reflecting on the implications and probabilities of meeting the aforementioned deadlines, as well as on the sustainability of past and future pumping trends. Policy considerations aside, modelling results suggest that groundwater extractions could be augmented considerably without causing a significant decrease in streamflows.

A simple hydrologic framework for simulating wetlands in climate and earth system models

Wetlands are ecosystems of important functions in the earth’s climate system. Through relatively high evapotranspiration, they affect surface water and energy exchange with the atmosphere directly influencing the physical climate. Through CH4, CO2 and N2O fluxes, they regulate the biogeochemical cycles, indirectly influencing the physical climate. However, current models do not explicitly include the water table, present under all large and stable wetlands; model wetlands are identified as flat land with wet soil resulting from precipitation events. That is, the wetlands are only ‘wetted’ from above but not from below by the high water table. Furthermore, without the knowledge of the water table position, estimates of CH4 and other gases (e.g., CO2 and N2O) are poorly constrained. We present a simple hydrologic framework for simulating wetlands based on water table depth. A synthesis of hydrologic controls on wetlands highlights the key role that groundwater plays. It directly feeds wetlands, supports surface-water fed wetlands by maintaining a saturated substrate, and links land drainage to sea level by impeding drainage in lowlands. Forced by routine climate model output (precipitation–evapotranspiration-surface runoff), land topography, and sea level, we simulate the present-day water table in North America at the 1 km scale. We validate the simulation with water table observations and compare regions of shallow water table to mapped wetlands. Our results show that the framework captures the salient features of wetland distribution and extent at regional and continental scales, a direct result of large-scale groundwater convergence that nourishes the lowlands even in arid climates. The low requirement of forcing and computation make the framework easy to adopt in climate and earth system models for simulating wetland responses to climate and sea level change for the present, paleo reconstructions, and future projections.

Design and risk assessment tool for vegetative treatment areas receiving agricultural wastewater: Preliminary results

Vegetative treatment areas (VTAs) are commonly being used as an alternative method of agricultural process wastewater treatment. However, it is also apparent that to completely prevent discharge of pollutants to the surrounding environment, settling of particulates and bound constituents from overland flow through VTAs is not sufficient. For effective remediation of dissolved agricultural pollutants, VTAs must infiltrate incoming wastewater. A simple water balance model for predicting VTA soil saturation and surface discharge in landscapes characterized by sloping terrain and a shallow restrictive layer is presented and discussed. The model accounts for the cumulative effect of successive rainfall events and wastewater input on soil moisture status and depth to water table. Nash-Sutcliffe efficiencies ranged from 0.65 to 0.81 for modeled and observed water table elevations after calibration of saturated hydraulic conductivity. Precipitation data from relatively low, average, and high annual rainfall years were used with soil, site, and contributing area data from an example VTA for simulations and comparisons. Model sensitivity to VTA width and contributing area (i.e. barnyard, feedlot, silage bunker, etc.) curve number was also investigated. Results of this analysis indicate that VTAs should be located on steeper slopes with deeper, more-permeable soils, which effectively lowers the shallow water table. In sloping landscapes (>2%), this model provides practitioners an easy-to-use VTA design and/or risk assessment tool that is more hydrological process-based than current methods.

Detection of sub-surface waterlogging using Terra-1 MODIS data

Waterlogging due to rising ground watertable, being a sub-surface phenomenon, is not amenable to detection by optical remote sensing. Microwave and thermal sensor data have, however, shown some promise in the detection of sub-surface waterlogging. The present study was taken up to evaluate the potential of near-IR, short-wave IR (SWIR) and thermal-IR data from Moderate Resolution Imaging Spectrometer (MODIS) aboard Terra-1 acquired during day-and-night time postmonsoon data for detection of sub-surface waterlogging. The approach involves retrieval of day-and-night land surface temperature (LST), generation of normalized difference of channel-2 and 6 (ND26); 2 and 7 (ND27); ground truth collection involving concurrent ground water table observations to satellite date of pass, thresholding of normalized differences (NDs) and correlating the NDs with depth of ground water table. Amongst various spectral indices, day and night-time LST difference (DLST) and night-time LST have been found to correlate well with the incidence of waterlogging (water table depth < 2m), followed by normalized difference of band-2 (841–876 nm) and band-7 (2105–2155 nm). The sensitivity of threshold limits for these indices was maximum for DLST followed by ND26 and ND27. Poor accuracy of detecting sub-surface waterlogging with thermal bands during day time is attributed to the non-corresponding of the time of Terra MODIS data acquisitions with thermal maxima of the terrain. Though the ND27 gave better accuracy to detect subsurface waterlogging, it is very sensitive to threshold limits.

Approach Using Active Groundwater Storage for Hydrologic Model Calibration in West-Central Florida

Hydrologic model calibration is always a challenging and tedious process especially for the calibration of complex models, which includes continuous hydrograph models, requires sophisticated calibration methods. The Hydrologic Simulation Program-FORTRAN (HSPF) is one of the popular and powerful time variable hydrologic models. However, in order to improve the assessment of hydrologic activities in shallow ground water settings, the model needs to be reliably calibrated for ground water contribution. Little guidance is provided in the literature concerning the manner of this contribution. In fact, the most common calibration of HSPF uses subjective parameter fitting and focuses on the attainment of statistical goodness of fit of runoff fluxes and water levels, ignoring ground water components. The goal of this research is using a different approach to calibrate HSPF with observed water table records. In this study, HSPF is applied on a small area in west-central Florida and calibrated by comparing active ground water storage to well elevation records in range land and forested land covers. The Nash-Sutcliffe efficiency and correlation coefficient computed using observed and simulated daily flows are 0.91 and 0.96 at Peace River, respectively, also with good fair results for other stations in the model domain. The study shows that improved calibration of the model can be achieved if active ground water storage and well records are compared for timing and magnitude of fluctuations.

Relative Elevation Topographic Surface Modelling of a Large Alluvial River Floodplain and Applications for the Study and Management of Riparian Landscapes

This paper presents a novel and useful spatial modelling technique to create a topographic surface that estimates a floodplain's elevation relative to the average low-flow water surface elevation of a river channel. This model was applied to a 121 km study area of the middle Sacramento River, California, USA, where it was tested as a surrogate for observed water table depth and an observed 3.3 year recurrence interval flood inundation surface using independent data sets. The modelled relative elevation topographic surface correlated significantly (p < 0.005) to observed well water depths suggesting that the modelled surface reflected a reasonable approximation of vertical distance to the water table. Results from a flood inundation pattern analysis indicated an overall accuracy of 79% for correctly predicting inundated and non-inundated zones. The model was then used to measure relative channel bank height and the distribution of riparian plant communities to examine landscape ecological relationships.

Contribution to irrigation from shallow water table under field conditions

The mathematical model SWBACROS was applied to estimate the contribution of a shallow groundwater to the water needs of a maize crop. The model was applied with the top and boundary conditions defined by the observed irrigation/rainfall events and the observed water table depth. The simulated water contents of the top zone were very close to the observed values. Furthermore the model was applied with an assumed free drainage bottom boundary condition. The difference of the computed water content profiles under the presence and absence of the water table gave a very good estimate of the capillary rise. It was found that under the specific field conditions about 3.6 mm/day of the water in the root zone originated from the shallow water table, which amounts to about 18% of the water, which was transpired by the maize crop.

Incorporating water table dynamics in climate modeling: 1. Water table observations and equilibrium water table simulations

Soil moisture is a key participant in land‐atmosphere interactions and an important determinant of terrestrial climate. In regions where the water table is shallow, soil moisture is coupled to the water table. This paper is the first of a two‐part study to quantify this coupling and explore its implications in the context of climate modeling. We examine the observed water table depth in the lower 48 states of the United States in search of salient spatial and temporal features that are relevant to climate dynamics. As a means to interpolate and synthesize the scattered observations, we use a simple two‐dimensional groundwater flow model to construct an equilibrium water table as a result of long‐term climatic and geologic forcing. Model simulations suggest that the water table depth exhibits spatial organization at watershed, regional, and continental scales, which may have implications for the spatial organization of soil moisture at similar scales. The observations suggest that water table depth varies at diurnal, event, seasonal, and interannual scales, which may have implications for soil moisture memory at these scales.

Simulation of water table level and peat temperatures in boreal peatlands

Hydrology controls the physical, chemical and biological processes in peatlands and hence could be the most important process regulating their function, development and characteristic biogeochemistry. Models describing hydrological processes and soil heat exchange phenomena are an important tool in understanding the peatland biogeochemical cycles of C and N. Presented in this paper is a peat soil climate model that uses weather data, mire site and peat characteristics as input data. The model is related to the heterogeneous peat characteristics (e.g. bulk density, degree of humification, remains of Sphagnum, Carex, Eriophorum and woody tissues), which in turn vary with depth. Evapotranspiration is related to weather parameters and tree stand characteristics. The model simulations were validated against the observed data collected during 1993 and 1994 growing seasons at a fen and bog sites at the Lakkasuo mire complex located in central Finland. The observed and simulated day-to-day variations in water table levels and soil temperatures during both seasons appeared overall to be in phase. As the model lacks the data needed to characterize the peat matrix hysteresis, the model was found lacking in response to wetting and drying cycles. The mean differences between simulated and observed water table levels during 1994 were −0.3 ± 1.3 and −0.5 ± 2.3 cm for the fen and bog sites, respectively. The model was found to be sensitive to mire surface characteristics and evapotranspiration, particularly for the bog site with an unsaturated zone deeper than at the fen. The absolute mean differences between the simulated and measured peat temperatures from 5 to 150 cm were less than 1.0 °C with a maximal standard deviation of 1.6 °C. The model deviations for the upper layers showed larger variations compared to deeper layers, implying a greater accuracy in defining the lower boundary of the thermal regime within the peat column.

Evaluation of DRAINMOD for predicting water table heights in irrigated fields at the Jordan Valley

A detailed field experiment was carried out in the Jordan Valley, south of Lake Kinneret, Israel for evaluation of the water management model DRAINMOD. This field was chosen to represent the local agro-climate conditions of that zone. Banana crop was grown and was irrigated daily with about 3200 mm/year and 0.5 leaching fraction. Subsurface drainage system with 2.5 m drain depth and 160 m drain spacing existed in the field. The water table depth was measured with about 100 piezometers, in which most of them were observed weekly, and four were continuosly recording piezometers. Five identical drainage plots were selected, out of 10 existing, as replicates for the evaluation of DRAINMOD. Deviations in a range of 0.3–1.7 m between observed water table depth and that simulated by DRAINMOD were found in four out of the five replicates. A reasonable agreement was found only in one drainage plot out of the five tested. These findings contradict the world wide convention that DRAINMOD simulation is in a good agreement with observed field data. An additional study was therefore conducted to explore the reasons for these large deviations. Three reasons were suggested: (i) a strong side effect by the Jordan River, which flows some 350 m west to the test field; a very steep 4.6% gradient was found toward the Jordan River; (ii) presence of sandy permeable layers below the depth of the drains which magnifies the boundary condition effect of the Jordan River; (iii) a very significant component of deep and lateral seepage (more than 50% of the yearly irrigation plus rainfall). A combination of these three reasons was suggested as an explanation to the apparent large disagreement. It was therefore recommended not to use DRAINMOD or similar vertical flow models for simulation of water table depths in irrigated fields with subsurface drain pipe systems in the Jordan Valley.

Evapotranspiration of Two Vegetation Covers in a Shallow Water Table Environment

A method is introduced to estimate evapotranspiration (ET) in shallow water table environments. The method involves measuring the diurnal fluctuations in total soil moisture above the water table to estimate (i) the net lateral and vertical subsurface flux in the aquifer and (ii) evapotranspiration from the vegetation cover. In a hillslope discharge zone, the net lateral subsurface flux was calculated from the recovery rate of soil moisture between midnight and 0400 h. Evapotranspiration was then estimated from a daily water balance in a soil column that included the water table. The method was tested on two vegetation covers, a pasture in a groundwater recharge area, and a riparian zone with woody vegetation in a groundwater discharge area. A moisture probe carrying eight sensors was used in each area to estimate the total soil moisture in a sandy soil environment. The observed water table fluctuated between land surface and a depth of 1.2 m during the study period, allowing observation and estimation of the total soil moisture in a soil column that included the water table. The results of this investigation support another hypothesis that, in humid, shallow water table environments, ET demand may be supported by adjacent ecosystems. This method provided reasonable results for the two landscapes investigated and was able to capture the variability of evapotranspiration in heterogeneous vegetation covers. It provided a relatively inexpensive alternative to characterize ET within regionally heterogeneous but microhomogenous landscapes. Though tested for coarse‐textured soil, the method involving soil moisture monitoring can be easily adapted to other soil types with shallow water table. Another advantage of using this method is that ET can be successfully estimated without detailed knowledge of soil hydraulic properties, subsurface flow patterns, or vegetation characteristics.