Quantification of vertical water fluxes in the vadose zone using particle‐size distribution and pedology‐based approaches to model soil heterogeneities

One‐dimensional flow simulations were conducted at four locations of the shallow alluvial aquifer of the upper Rhine River (at the Erstein polder) to quantify the time‐dependent moisture distribution, the water flux and the water volume infiltrated in the unsaturated zone as a function of soil heterogeneities during a five‐day‐long flooding event. Three methods of estimating the hydraulic parameters of soil in the vadose zone were tested. They are based on the following: (1) experimental data, (2) soil particle‐size distribution and (3) pedology information on soils. Water fluxes calculated from modelling approaches 2 and 3 were compared with those of the experiment‐based values and the effect of these differences on the arrival time and velocity of water at the water table were analysed. Major differences in water fluxes were found among the methods of estimating the hydrodynamic parameters. At the Terrace location, the groundwater recharge predicted using soil data from methods 1 and 2 are approximately 4500 and 2400 mm, respectively. Flow simulations using soil data and the experiment‐based method show the highest velocities of infiltrating water at the soil surface and largest volume of groundwater infiltration but result in the lowest centres of the moisture content mass. The results obtained using soil data based on the pedological method are similar to those calculated using soil parameters based on the particle‐size distribution of extracted soil samples. Water pressure profiles calculated on Terrace and Channel location, 3 and 7 days after the inundation event agreed reasonably well with those observed when using hydrodynamic parameters from the experiment‐based method. However, the flow model using the pedology‐based parameters largely underestimates the time needed to achieve hydrostatic conditions of the soil water profile once water flooding at the soil surface stops. This can be mainly attributed to the low values of estimated van Genuchten parameter α. Copyright © 2012 John Wiley & Sons, Ltd.     

Artificial recharge efficiency assessment by soil water balance and modelling approaches in a multi-layered vadose zone in a dry region

To assess recharge through floodwater spreading, three wells, approx. 30 m deep, were dug in a 35-year-old basin in southern Iran. Hydraulic parameters of the layers were measured. One well was equipped with pre-calibrated time domain reflectometry (TDR) sensors. The soil moisture was measured continuously before and after events. Rainfall, ponding depth and the duration of the flooding events were also measured. Recharge was assessed by the soil water balance method, and by calibrated (inverse solution) HYDRUS-1D. The results show that the 15 wetting front was interrupted at a layer with fine soil accumulation over a coarse layer at the depth of approx. 4 m. This seemed to occur due to fingering flow. Estimation of recharge by the soil water balance and modelling approaches showed a downward water flux of 55 and 57% of impounded floodwater, respectively.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Application of VLEACH to Vadose Zone Transport of VOCs at an Arizona Superfund Site

Cleanup standards for volatile organic compounds in thick vadose zones can be based on indirect risk (transport to ground water) when contamination is below depths of significant direct risk. At one Arizona Superfund site, a one-dimensional vadose zone transport model (VLE-ACH) was used to estimate the continued transport of VOCs from the vadose zone to ground water. VLEACH is a relatively simple and readily available model that proved useful for estimating indirect risk from VOCs in the vadose zone at this site. The estimates of total soil concentrations used as initial conditions for VLF.ACH incorporated a variety of data from the site. Soil gas concentrations were found to be more useful than soil matrix data for estimating total soil concentrations at this arid-zone site. A simple mixing cell model was used with the VLEACH-derived mass loading estimates from the vadose zone over time to estimate the resulting changes in ground water concentrations. For this site, the results of the linked VLEACH/mixing cell simulations indicate it is likely that the federal MCI. for TCE will be exceeded in underlying ground water if remedial action on I he vadose zone is not pursued.     

Hydrological cycle and water balance estimates for the megadune–lake region of the Badain Jaran Desert,China

The hydrology and water balance of megadunes and lakes have been investigated in the Badain Jaran Desert of China. Field observations and analyses of sand layer water content, field capacity, secondary salt content, and grain size reveal 3 types of important natural phenomenon: (a) vegetation bands on the leeward slope of the megadunes reflect the hydrological regime within the sandy vadose zone; (b) seepage, wet sand deposits, and secondary salt deposits indicate the pattern of water movement within the sandy vadose zone; (c) zones of groundwater seeps and descending springs around the lakes reflect the influence of the local topography on the hydrological regime of the megadunes. The seepage exposed on the sloping surface of the megadunes and gravity water contained within the sand layer confirm the occurrence of preferential flow within the vadose zone of the megadunes. Alternating layers of coarse and fine sand create the conditions for the formation of preferential flows. The preferential flows promote movement of water within the sand layer water that leads to deep penetration of water within the megadunes and ultimately to the recharging of groundwater and lake water. Our results indicate that a positive water balance promotes recharge of the megadunes, which depends on the high permeability of the megadune material, the shallow depth of the surface sand layer affected by evaporation, the occurrence of rainfall events exceeding 15 mm, and the sparse vegetation cover. Water balance estimates indicate that the annual water storage of the megadunes is about 7.5 mm, accounting for only 8% of annual precipitation; however, the shallow groundwater per unit area under the megadunes receives only 3.6% of annual precipitation, but it is still able to maintain a dynamic balance of the lake water. From a water budget perspective, the annual water storage in the megadunes is sufficient to serve as a recharge source for lake water, thereby enabling the long‐term persistence of the lakes. Overall, our findings demonstrate that precipitation is a significant component of the hydrological cycle in arid deserts.     

Evaluating soil water routing approaches in watershed-scale,ecohydrologic modelling

Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.     

A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.     

Physical controls on septic leachate movement in the vadose zone at the hillslope scale,Putnam County,New York,USA

The fate and transport of contaminants in the vicinity of septic fields remains poorly understood in many hydrogeomorphological environments. We report hydrometric data from an intensive hillslope‐scale experiment conducted between 29 August and 11 November 1998 at a residential leach field in New York State. The objective of our study was to characterize water flux within the vadose zone, understand the physical controls on the flux, and predict how this ultimately will affect subsurface water quality. Soil‐water flux was calculated using matric potential measurements from a network of 25 tensiometer nests, each nest consisting of three tensiometers installed to depths of 10, 50 and 130 cm. Unsaturated hydraulic conductivity curves were derived at each depth from field‐determined time‐domain reflectometry–tensiometry moisture‐release curves and borehole permeametry measurements. Flownets indicated that a strong upward flux of soil water occurred between rainstorms. Following the onset of (typically convective) rainfall, low near‐surface matric potentials were rapidly converted to near‐saturated and saturated conditions, promoting steep vertical gradients through the near‐surface horizons of the hillslope. Lateral hydraulic gradients were typically 10 times smaller than the vertical gradients. Resultant flow vectors showed that the flux was predominantly vertical through the vadose zone, and that the flux response to precipitation was short‐lived. The flux response was controlled primarily by the shape of the unsaturated hydraulic conductivity curves, which indicated a rapid loss of conductivity below saturation. Thus, soil water had a very high residence time in the vadose zone. The absence of rapid wetting at 130 cm and the delayed and small phreatic zone response to rainfall indicated that water movement through macropores did not occur on this hillslope. These results are consistent with a Cl tracing experiment, which demonstrated that the tracer was retained in the vadose zone for several months after injection to the system. Copyright © 2002 John Wiley & Sons, Ltd.     

High-characterization of Vadose Zone Dynamics in Limestone Underground Quarries by Time Domain Reflectometry

The aim of this paper is to study the vadose zone dynamics during a hydrological cycle. The application of the Time Domain Reflectometry (TDR) method to determine the water content of porous rock has been widely investigated. More than 657 point measurements of rock water content observed during a hydrological cycle and distributed among three abandoned underground quarries in Gironde, France show a permanently undersaturated limestone. Two periods of maximum water content correspond to two occurring effective precipitations. The dephasing and amplitude attenuation of the hydraulic wave with the depth can be modelled and explained by the physical properties of the porous medium. Maps of the spatial distribution of water content show the heterogeneity of water flow in a vadose zone.     

Extinction depth and evapotranspiration from ground water under selected land covers

In many landscapes, vegetation extracts water from both the unsaturated and the saturated zones. The partitioning of evapotranspiration (ET) into vadose zone evapotranspiration and ground water evapotranspiration (GWET) is complex because it depends on land cover and subsurface characteristics. Traditionally, the GWET fraction is assumed to decay with increasing depth to the water table (DTWT), attaining a value of 0 at what is termed the extinction depth. A simple assumption of linear decay with depth is often used but has never been rigorously examined using unsaturated-saturated flow simulations. Furthermore, it is not well understood how to relate extinction depths to characteristics of land cover and soil texture. In this work, variable saturation flow theory is used to simulate GWET for three land covers and a range of soil properties under drying soil conditions. For a water table within half a meter of the land surface, nearly all ET is extracted from ground water due to the close hydraulic connection between the unsaturated and the saturated zones. For deep-rooted vegetation, the decoupling of ground water and vadose zone was found to begin at water table depths between 30 and 100 cm, depending on the soil texture. The decline of ET with DTWT is better simulated by an exponential decay function than the commonly used linear decay. A comparison with field data is consistent with the findings of this study. Tables are provided to vary the extinction depth for heterogeneous landscapes with different vegetation cover and soil properties.     

Hydraulic model tests for propagation of flow and sediment in floods due to breaking of a natural landslide dam during a mountainous torrent

During mountain torrents, large-magnitude floods may result from heavy rainfall and cause the breakage of landslide dams naturally formed by heavy rainfall, earthquakes, and so on. The characteristics of longitudinal spreading of clear water discharge and changes in flow depth must be clarified because the changes in peak depth have not yet been examined in steep-slope torrents and because there are few data on spreading of flash floods and related sedimentation in mountainous torrents. In the present study, experimental data were collected through hydraulic model tests over a rigid bed, and the spreading of water, fine sediment, bed load, and large boulders due to flooding are discussed assuming that flash flooding/debris flows occur in the upstream reach. The effects of changes in flow width, such as expansions and contractions in the flow width, as well as changes in meandering channels, sediment transportation, and spreading flow depth resulting from bores are examined using flume data for a steep-slope torrent. The data obtained in the present study reveal that fine sediment components are transported to the downstream reach if large-magnitude floods occur and that the spreading rate and peak lags of the fine sediment and water level indicate the occurrence of a flood in the upstream reach.     

Coastal aquifer response to extreme storm events in Emilia‐Romagna,Italy

With global warming and sea level rise, many coastal systems will experience increased levels of inundation and storm flooding, especially along sandy lowland coastal areas, such as the Northern Adriatic coast (Italy). Understanding how extreme events may directly affect groundwater hydrology in shallow unconfined coastal aquifers is important to assess coastal vulnerability and quantify freshwater resources. This study investigates shallow coastal aquifer response to storm events. The transitory and permanent effects of storm waves are evaluated through the real time monitoring of groundwater and soil parameters, in order to characterize both the saturated and unsaturated portions of the coastal aquifer of Ravenna and Ferrara (southern Po Delta, Italy). Results highlight a general increase in hydraulic head and soil moisture, along with a decrease in groundwater salinity and pore water salinity due to rainfall infiltration during the 2 days storm event. The only exceptions are represented by the observation wells in proximity to the coastline (within 100 m), which recorded a temporary increase in soil and water salinity caused by the exceptional high waves, which persist on top of the dune crest during the storm event. This generates a saline plume that infiltrates through the vadose zone down to the saturated portion of the aquifer causing a temporary disappearance of the freshwater lens generally present, although limited in size, below the coastal dunes. Despite the high hydraulic conductivity, the aquifer system does not quickly recover the pre‐storm equilibrium and the storm effects are evident in groundwater and soil parameters after 10 days past the storm overwash recess.     

Constraining a Flow Model with Field Measurements to Assess Water Transit Time Through a Vadose Zone

The modeling of thick vadose zones is particularly challenging because of difficulties in collecting a variety of measured sediment properties, which are required for parameterizing the model. Some models rely on synthetic data, whereas others are simplified by running as homogeneous sediment domains and relying on a single set of sediment properties. Few studies have simulated flow processes through a thick vadose zone using real and comprehensive data sets comprising multiple measurements. Here, we develop a flow model for a 7-m-thick vadose zone. This model, combining the numerical codes CTRAN/W with SEEP/W, includes the measured sediment hydraulic properties of the investigated vadose zone and incorporates the actual climate and subsurface conditions of the study site (precipitations, water-table elevations, and stable isotope data). The model is calibrated by fitting the simulated and measured vertical profiles of water content. Our flow model calculates a transit time of 1 year for the travel of water through the 7-m vadose zone; this estimate matches stable isotope-based results obtained previously for this site. A homogeneous sediment domain flow model, which considers only a single set of sediment properties, produces a transit time that is approximately half the duration of that of the heterogeneous flow model. This difference highlights the importance of assuming heterogeneous material within models of thick vadose zones and testifies to the advantage gained when using real sediment hydraulic properties to parametrize a flow model.     

Soil water dynamics and water balance on a tropical coral island

Understanding soil water dynamics and the water balance of tropical coral islands is important for the utilization and management of their limited freshwater resources, which is only from rainfall. However, there is a significant knowledge gap in the influence of soil water on the water cycle of coral islands. Soil water dynamics and the water balance of Zhaoshu Island, Xisha Archipelago were thus investigated using soil moisture measurements and the Hydrus-1D model from October 2018 to September 2019. Over the study period, vegetation transpiration, soil evaporation, groundwater recharge and storage in the vadose zone were approximately 196, 330, 365 and 20 mm, occupying 22%, 36%, 40% and 2% of annual rainfall total (911 mm), respectively. For the wet season (from May to October) these values became 75, 202, 455 and 40 mm, occupying 10%, 26% and 59% and 5% of the seasonal rainfall total (772 mm), respectively. During the dry season (from November to April), a dry soil layer between 40 and 120 cm depth of the soil profile was identified that prevented water exchange between the upper soil layers and the groundwater resulting in the development of deep roots so that vegetation could extract groundwater to supplement their water requirements. Vegetation not only consumes all dry season rainfall (140 mm) but extracts water deeply from groundwater (90 mm) as well as from the vadose layer (20 mm). As such, the vegetation appears to be groundwater-dependent ecosystems. The research results aid us to better understand the process of water dynamics on coral islands and to protect coral island ecosystems.     

Seepage response along an alluvial valley in a semi‐arid catchment in north‐central California

Over a period of 12 months, soil moisture content and potential was monitored in an annual‐grass‐dominated 20 ha catchment in order to determine flow paths leading to exfiltration at the catchment outlet. Water was found to enter the catchment valley either through flow originating in the slopes or through surface infiltration during rainfall events. Although subsurface flow from the slopes to the catchment outlet occurred throughout the year, surface recharge was restricted to a few events during the wet season. In the deeper saturated profile of the valley, flow was directed upwards along the valley edges and gradually became horizontal towards the central axis of the valley. During the peak of the rainfall season, horizontal flow close to the catchment outlet intercepted the gradually sloping surface, resulting in exfiltration. Plants influenced the hydrology of the catchment by removing moisture from the root zone during spring and early summer, resulting in evapotranspiration losses from the vadose zone. Heterogeneities within the valley soil were evident as variable‐permeability layers that resulted in a seasonally confined water table within the valley. This investigation shows that the vadose zone plays an important role in redistributing surface recharge and emphasizes the importance of accounting for effective moisture in low‐yielding catchments with ephemeral surface runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

A scale‐invariant property of the water retention curve in weakly heterogeneous vadose zones

Abrupt changes of hydraulic properties in a vadose zone are modelled within a stochastic framework, which regards the saturated conductivity and parameters related to the distribution of soil pores as stationary, log‐normally distributed, random space functions. As a consequence, flow variables become random fields, and we aim at deriving an effective Richards equation. To obtain the latter, we adopt a perturbation expansion truncated at the first order (weakly heterogeneous media), which leads to the effective hydraulic conductivity and water retention curves. Overall, the effective properties are scale dependent. However, within the proposed framework, we demonstrate that the inflection point of the laboratory scale retention curve is not affected by the heterogeneity of the vadose zone. Finally, to illustrate the quantitative implications of our results, we consider a monitoring experiment at field scale, and we show how our approach leads to an effective water retention curve, which differs significantly from that which would be obtained without accounting for the above scale‐invariance property.     

Analysis of Water Saturation, NAPL Content, Degradation Half-Life, and Lower Boundary Conditions on VOC Transport Modeling: Implications for Closure of Soil Venting Systems

Simulations using a one-dimensional, analytical, vadose zone, solute-transport screening code (VFLUX) were conducted to assess the effect of water saturation, NAPL saturation, degradation half-life, and boundary conditions at the vadose zone/ground water interface on model output. At high initial soil concentrations, model output was significantly affected by input parameters and lower boundary conditions yet still resulted in consistent decision-making to initiate or continue venting application. At lower soil concentrations, however, typical of what is observed after prolonged venting application, differences in model input and selection of lower boundary conditions resulted in inconsistent decision-making. Specifically, under conditions of low water saturation, use of a first-type, time-dependent lower boundary condition indicated that the primary direction of mass flux was from ground water to the vadose zone, suggesting little benefit from continued venting application. Use of a finite, zero-gradient lower boundary condition, though, indicated continued mass flux from the vadose zone to ground water, suggesting a continued need for venting application. In this situation, sensitivity analysis of input parameters, selection of boundary conditions, and consideration of overall objectives in vadose zone modeling become critical in regulatory decision-making.     

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.     

Simulating soil‐water movement under a hedgerow surrounding a bottomland reveals the importance of transpiration in water balance

The objective of this study was to quantify components of the water balance related to root‐water uptake in the soil below a hedgerow. At this local scale, a two‐dimensional (2D) flow domain in the x–z plane 6 m long and 1·55 m deep was considered. An attempt was made to estimate transpiration using a simulation model. The SWMS‐2D model was modified and used to simulate temporally and spatially heterogeneous boundary conditions. A function with a variable spatial distribution of root‐water uptake was considered, and model calibration was performed by adjusting this root‐water uptake distribution. Observed data from a previous field study were compared against model predictions. During the validation step, satisfactory agreement was obtained, as the difference between observed and modelled pressure head values was less than 50 cm for 80% of the study data. Hedge transpiration capacity is a significant component of soil‐water balance in the summer, when predicted transpiration reaches about 5·6 mm day^?1. One of the most important findings is that hedge transpiration is nearly twice that of a forest canopy. In addition, soil‐water content is significantly different whether downslope or upslope depending on the root‐water uptake. The high transpiration rate was mainly due to the presence of a shallow water table below the hedgerow trees. Soil‐water content was not a limiting factor for transpiration in this context, as it could be in one with a much deeper water table. Hedgerow tree transpiration exerts a strong impact not only on water content within the vadose zone but also on the water‐table profile along the transect. Results obtained at the local scale reveal that the global impact of hedges at the catchment scale has been underestimated in the past. Transpiration rate exerts a major influence on water balance at both the seasonal and annual scales for watersheds with a dense network of hedgerows. Copyright © 2007 John Wiley & Sons, Ltd.     

An analytical solution for predicting the transient seepage from a subsurface drainage system

Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems.