The use of GIS and remote sensing for the assessment of waterlogging in the dryland irrigated catchments of Farafra Oasis,Egypt

This paper investigates the interplay of the hydrogeological characteristics, soil properties and recent land reclamation projects on the distribution of waterlogging and salinization within the Farafra Oasis. The multi‐temporal remote sensing data and field observations show that new reclaimed areas have been recently cultivated in distant areas from the old agricultural land. These new cultivations have developed widespread waterlogging, seepage channels and soil salinization. Analyses of the Shuttle Radar Topography Mission digital elevation model (DEM) showed that both old and new agricultural areas are located within same closed drainage basin. The fluvial channels of these catchments, which were developed during wet climatic pluvial, have largely been obliterated by the prevailing aridity and often buried under aeolian deposits. However, the new cultivations have been developed on the fingertips of these fluvial channels, while the old fields occupy the low‐level playas. The soil of the new cultivated areas is mainly lithic with a high calcium carbonate content, thus limiting the downward percolation of excess irrigation water and therefore developing perched water table and seepage through the palaeo‐channels. The automatically extracted drainage networks from DEM resemble fluvial patterns and coincide with the seepage channels slowly heading toward old cultivation. The inactive alluvial channels and landforms have to be considered when planning for new cultivation in dryland catchments to better control waterlogging and salinization hazard. It is highly recommended that newly developed seepage channels have to be detected and intercepted before reaching old agriculture areas. Therefore, the ‘dry‐drainage’ concept can be implemented as the seepage water can be conveyed into nearby playas reserved for evaporation. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling subsurface flow conditions in a salinized catchment in south‐western Australia,with a view to improving management practices

Finite element modelling of the saturated–unsaturated surface–subsurface flow mechanisms operative in a small salinized catchment in south‐western Australia was used to help define the flow system and explain the causes of waterlogging and salinization there. Data available at the site from a previous study were used to obtain a first approximation to the flow system. Altering the properties of some of the strata gave a closer calibration. It was found that the modelled saturated hydraulic conductivity of the B horizon in the duplex soil zone needed to be at least an order of magnitude lower than that measured in order to reproduce the perching conditions observed in the field. Also, the model indicated the influence of a doleritic dyke, whose presence was confirmed by field measurement. Our analysis showed that there were two main flow systems operating in the hillslope. The first, and most dominant, was the recharge occurring through the upslope gradational soil zone and percolating down to both the deeply weathered regolith and the basal aquifer. The second flow system is an unsaturated flow system operating in the high permeability A horizon in the downslope duplex soil zone. The first system is primarily responsible for the saline seepage zone in the valley bottom. The second contributes to the waterlogging and perching occurring upslope of the seepage zone. Vertical flow through the higher permeability B horizon in the gradational soil zone in the upper slopes is a major contributor of recharge. Recharge by flow through macropores occurs where, but only where, perched aquifers develop and allow the macropores to be activated. Areas with perched aquifers occurred in downslope locations and near a doleritic dyke located upslope. Thus, the area where macropore recharge occurred was not large. The recharge rate required to maintain the piezometric levels at present values is only about 30 mm/yr (about 5% of the annual rainfall). The piezometric levels under the upper part of the catchment varied greatly with only small changes in recharge rate. A 50% reduction in recharge rate had the effect of reducing the length of the seepage zone at the end of winter by 40%. Changes in recharge rate had little effect on the extent of the perched aquifer at the end of winter. Deep‐rooted perennial forages, shrubs or trees on the gradational soil zone in the upper part of the catchment and on the zones upslope of geological barriers to flow would be required to reduce the recharge and to allow for rehabilitation of the saline valley floor. Waterlogging associated with the perched water table in the bottom part of the catchment would be best addressed by tree plantations located just upslope of the salinized zone in the valley floor. Copyright © 1999 John Wiley & Sons, Ltd.     

Predicting Vertical LNAPL Distribution in the Subsurface under the Fluctuating Water Table Effect

The present study proposes a methodology for predicting the vertical light nonaqueous-phase liquids (LNAPLs) distribution within an aquifer by considering the influence of water table fluctuations. The LNAPL distribution is predicted by combining (1) information on air/LNAPL and LNAPL/water interface elevations with (2) the initial elevation of the water table without LNAPL effect. Data used in the present study were collected during groundwater monitoring undertaken over a period of 4 months at a LNAPL-impacted observation well. In this study, the water table fluctuations raised the free LNAPL in the subsurface to an elevation of 206.63 m, while the lowest elevation was 205.70 m, forming a thickness of 0.93 m of LNAPL-impacted soil. Results show that the apparent LNAPL thickness in the observation well is found to be three times greater than the actual free LNAPL thickness in soil; a finding that agrees with previous studies reporting that apparent LNAPL thickness in observation wells typically exceeds the free LNAPL thickness within soil by a factor estimated to range between 2 and 10. The present study provides insights concerning the transient variation of LNAPL distribution within the subsurface and highlights the capability of the proposed methodology to mathematically predict the actual LNAPL thickness in the subsurface, without the need to conduct laborious field tests. Practitioners can use the proposed methodology to determine by how much the water table should be lowered, through pumping, to isolate the LNAPL-impacted soil within the unsaturated zone, which can then be subjected to in situ vadose zone remedial treatment.     

Bridging the gap between numerical solutions of travel time distributions and analytical storage selection functions

Recent advancements in analytical solutions to quantify water and solute travel time distributions (TTDs) and the related StorAge Selection (SAS) functions synthesize catchment complexity into a simplified, lumped representation. Although these analytical approaches are efficient in application, they require rarely available long‐term and high‐frequency hydrochemical data for parameter estimation. Alternatively, integrated hydrologic models coupled to Lagrangian particle‐tracking approaches can directly simulate age under different catchment geometries and complexity, but at a greater computational expense. Here, we bridge the two approaches, using a physically based model to explore the uncertainty in the estimation of the SAS function shape. In particular, we study the influence of subsurface heterogeneity, interactions between distinct flow domains (i.e., the vadose zone and saturated groundwater), diversity of flow pathways, and recharge rate on the shape of TTDs and the SAS functions. We use an integrated hydrology model, ParFlow, linked with a particle‐tracking model, SLIM, to compute transient residence times (or ages) at every cell in the domain, facilitating a direct characterization of the SAS function. Steady‐state results reveal that the SAS function shape shows a wide range of variation with respect to the variability in the structure of subsurface heterogeneity. Ensembles of spatially correlated realizations of hydraulic conductivity indicate that the SAS functions in the saturated groundwater have an overall weak tendency toward sampling younger ages, whereas the vadose zone gives a strong preference for older ages. We further show that the influence of recharge rate on the TTD is tightly dependent on the variability of subsurface hydraulic conductivity.     

The role of ground water in arid/semiarid ecosystems, Northwest China

Ground water plays an important role in water supply and the ecology of arid to semiarid areas such as Northwest China, where the landscape is fragile due to frequent drought in the past few decades. This paper discusses the role of ground water in these ecosystems, including the effect of condensation water and water table depth on the growth of plants and degree of soil salinity. The paper also discusses the controlling process for land desertification and soil salinization in Northwest China. Water table depth is a key factor controlling the water balance, ground water flow, and salt transport in the vadose zone. The suitable water table depth for vegetation growth, which can prevent land desertification and soil salinization, is within a range of 2 to 4 m; the optimal depth is approximately 3 m. As examples, changes in ecosystems owing to water resources development in Tarim and Manas basins, Xinjiang, China, are discussed.     

Seasonal shifts in the solute ion ratios of vadose zone rock moisture from the Eel River Critical Zone Observatory

One of the greatest challenges in critical zone studies is to document the moisture dynamics, water flux,and solute chemistry of the unsaturated, fractured and weathered bedrock that lies between the soil and groundwater table. The central impediment to quantifying this component of the subsurface is the difficulty associated with direct observations. Here, we report solute chemistry as a function of depth collected over a full year across the shale-derived vadose zone of the Eel River Critical Zone Observatory using a set of novel sub-horizontal wellbores,referred to as the vadose zone monitoring system. The results of this first geochemical glimpse into the deep vadose zone indicate a dynamic temporal and depth-resolved structure. Major cation concentrations reflect seasonal changes in precipitation and water saturation, and normalized ratios span the full range of values reported for the world's largest rivers.     

A Three-Dimensional Analytical Model to Aid in Selecting Monitoring Locations in the Vadose Zone

Monitoring of the vadose zone is a potentially complex, time-consuming, and expensive problem. The location of monitoring points and selection of monitoring instruments can be optimized by using computer models. Numerical models developed for this purpose in the past have often been expensive and difficult to use. This paper describes a fast, three-dimensional, approximate analytical solution to the moisture content in the unsaturated zone. An analytical solution is available for steady-state drainage, whereas an approximate analytical solution is available for the transient case. The model will handle an arbitrary distribution of fluid sources, as well as vertical and horizontal impermeable boundaries.
The model may be applied to predict the incursion of fluid from accidental leakage or infiltration over large areas from unlined ponds and land treatment sites. The model is quite useful as an aid in designing monitoring or premonitoring programs near hazardous waste sites. Examples are presented to demonstrate the model's utility in estimating the maximum spread of a contaminant, the extent to which the fluid may spread with depth, the regions of high and low capillary pressure, and the non-linear behavior of the saturation when drainage from several sources in considered. These results are useful for the placement of monitoring locations and the selection of appropriate instruments, and as a tool in working with regulatory agencies to design monitoring programs. A glimpse of the future is necessary for today's planning. Computer models are some of the most useful crystal balls we have available.     

Importance of the vadose zone in analyses of unconfined aquifer tests

Analytical models commonly used to interpret unconfined aquifer tests have been based on upper-boundary (water table) conditions that do not adequately address effects of time-varying drainage from the vadose zone. As a result, measured and simulated drawdown data may not agree and hydraulic parameters may be inaccurately estimated. A 72-hour aquifer test conducted in Cape Cod, Massachusetts, in a slightly heterogeneous, coarse-grained, glacial outwash deposit was found to be a good candidate for testing models with different upper-boundary conditions. In general, under the commonly invoked assumption of instantaneous drainage, measured and simulated drawdowns were found to agree with one another only at late time and early time. In the intermediate-time range, because of delayed drainage, measured drawdowns always exceeded simulated values, most noticeably in piezometers located near the water table. To reduce these discrepancies, an analytical model was developed that can fully account for time-varying drainage given that the aquifer is not strongly heterogeneous. The approach is flexible as the model, which makes use of empirical relations, does not constrain drainage to follow any particular functional relation. By this approach, measured and simulated drawdowns agree over the complete time range, and the estimated parameters are consistent with prior studies and with what is known about the aquifer geometry, stratigraphy, and composition. By properly accounting for vadose zone drainage, it was found that realistic estimates of all hydraulic parameters, including specific yield, could be obtained with or without the use of late-time data.     

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.     

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.     

A Mathematical View of Water Table Fluctuations in a Shallow Aquifer in Brazil

Detailed monitoring of the groundwater table can provide important data about both short‐ and long‐term aquifer processes, including information useful for estimating recharge and facilitating groundwater modeling and remediation efforts. In this paper, we presents results of 4 years (2002 to 2005) of monitoring groundwater water levels in the Rio Claro Aquifer using observation wells drilled at the Rio Claro campus of São Paulo State University in Brazil. The data were used to follow natural periodic fluctuations in the water table, specifically those resulting from earth tides and seasonal recharge cycles. Statistical analyses included methods of time‐series analysis using Fourier analysis, cross‐correlation, and R/S analysis. Relationships could be established between rainfall and well recovery, as well as the persistence and degree of autocorrelation of the water table variations. We further used numerical solutions of the Richards equation to obtain estimates of the recharge rate and seasonable groundwater fluctuations. Seasonable soil moisture transit times through the vadose zone obtained with the numerical solution were very close to those obtained with the cross‐correlation analysis. We also employed a little‐used deep drainage boundary condition to obtain estimates of seasonable water table fluctuations, which were found to be consistent with observed transient groundwater levels during the period of study.     

Evidence for climatic and hillslope‐aspect controls on vadose zone hydrology and implications for saprolite weathering

Through the delivery of water in snowmelt, climate should govern the rate and extent of saprolite formation in snow‐dominated mountain watersheds, yet the mechanisms by which water flows deeply into regolith are largely unexplored. In this study we link rainfall, snow depth, and water content data from both soil and shallow saprolite to document vadose zone dynamics in two montane catchments over 2 years. Measurements of snow pack thickness and soil moisture reveal strong contrasts between north‐ and south‐facing slopes in both the timing of meltwater delivery and the duration of significant soil wetting in the shallow vadose zone. Despite similar magnitudes of snowmelt recharge, north‐facing slopes have higher sustained soil moisture compared to south‐facing slopes. To help interpret these observations, we use a 2D numerical model of vadose zone dynamics to calculate the expected space–time moisture patterns on an idealized hillslope under two wetting scenarios: a single sustained recharge pulse versus a set of short pulses. The model predicts that the duration of the recharge event exerts a stronger control on the depth and residence time of water in the upper unsaturated zone than the magnitude of the recharge event. Model calculations also imply that water should move more slowly through the subsurface and downward water flux should be substantially reduced when water is applied in several pulses rather than in one sustained event. The results suggest that thicker soil and more deeply weathered rock on north‐facing slopes may reflect greater water supply to the deep subsurface. Copyright © 2015 John Wiley & Sons, Ltd.     

The role of capillary forces in steady-state and transient seepage flows

Steady and transient seepage problems were investigated numerically using a dimensionless formulation for water flow in variably saturated, two-dimensional anisotropic and homogenous porous media. The dimensionless formulation combines the aspect ratio and the anisotropy ratio in one dimensionless parameter, M, and accounts for capillarity effects using the ratio of the domain height to the height of capillary fringe as another dimensionless parameter, alpha. Two domain geometries were considered: rectangular and trapezoidal. It was found that M and alpha play major roles in the development and the height of the seepage face. The height of the steady-state seepage face increased with the decreasing value of alpha. The effects of capillarity on the development of the transient seepage face were investigated by lowering the open water level on one side of the domain at various uniform velocities while the other side was kept at a constant value. It was found for most scenarios that decoupling between the water table and the open water level first occurred for domains with large alpha. In other words, the seepage face heights of these systems were larger than those with a small alpha. These results are the opposite of the steady-state results and are due to the fact that drainage of the pores was the major mechanism controlling the drop of the transient water table, especially at earlier times. A criterion for the decoupling of the water table from the falling open water level (i.e., the formation of the transient seepage face) was developed by Dracos (1965). It states that decoupling does not occur as long as the falling speed, VF, of the open water level is less than Vd = K(o) sin2 B/phi where K(o) is the saturated hydraulic conductivity of the soil, phi is the porosity, and beta is the angle of the exit face with the horizontal (i.e., beta = pi/2 implies a vertical face). Our investigation revealed that decoupling occurs for smaller falling velocities than Vd, and is caused by, among others, the fact that the water table is not tangent to the exit face but rather intersects it at a nonzero angle.     

Prediction of regional-scale groundwater recharge and nitrate storage in the vadose zone: A comparison between a global model and a regional model

Extensive nitrogen loads at the soil surface exceed plant uptake and soil biochemical capacity, and therefore lead to nitrogen accumulation in the deep vadose zone. Studies have shown that stored nitrogen in the vadose zone can eventually reach the water table and affect the quality of groundwater resources. Recently, global scale models have been implemented to quantify nitrate storage and nitrate travel time in the vadose zone. These global models are simplistic and relatively easy to implement and therefore facilitate analysis of the considered transport processes at a regional scale with no further requirements. However, the suitability of applying these models at a regional scale has not been tested. Here, we evaluate, for the first time, the performance and utility of global scale models at the regional scale. Applied to the Loess Plateau of China, we compare estimates of groundwater recharge and nitrate storage derived from global scale models with results from a regional scale approach utilizing the Richards and advection-dispersion equations. The estimated nitrate storage was compared to nitrate observations collected in the deep vadose zone (>50 m) at five sites across the Loess Plateau. Although both models predict similar spatial patterns of nitrate storage, the recharge fluxes were three times smaller and the nitrate storage was two times higher compared with the regional model. The results suggest that global scale models are a potentially useful screening tool, but require refinement for local scale applications.     

A diagnosis of sub-surface water table dynamics in low hydraulic conductivity soils in the sugar cane fields of Pongola,South Africa

Water and land are the two natural resources restraining crop production in South Africa. With the increasing demand for food, emphasis has shifted from the sole reliance on rain fed crop production, to irrigation. The deterioration in irrigation water quality from surface water sources is, however, posing a big challenge to the sustainability of irrigated crop production. This is because more water is required for leaching, resulting in shallow water tables in agricultural lands. The installation of well designed subsurface drainage systems alone is not enough; the provision of timely maintenance is also necessary. In this study, the extent and severity of problems as a consequence of shallow water tables and their possible causes were investigated at three sugarcane fields in Pongola, South Africa, having low hydraulic conductivity soils. Also investigated were soil salinity levels and the temporal variation in the salinity of the irrigation water. A water table map of a 32 ha sugarcane field was generated, using observed water table depth (WTD) data from 36 piezometers monitored from September 2011 to February 2012. Out of the total 32 ha under cultivation, 12% was found to be affected by shallow WTDs of less than the 1.0 m design WTD. The inability of natural drainage to cope with subsurface drainage needs and the poor maintenance of subsurface drainage systems contributed to the shallow water tables in the area. Furthermore, the currently adopted drainage design criteria also proved unsatisfactory with mean observed water table depth and drainage discharge (DD) of 20% and 50%, respectively, less than their respective design levels. The salinity of the irrigation water was, on average, 32% higher than threshold tolerance level of sugarcane. The root zone soil salinity levels at the three study sites were greater than the 1.7 dS m⁻¹ threshold for sugar cane. The subsurface drainage design criteria adopted at the site needs to be revisited by ensuring that the slope of the land is taken into consideration in the drainage design in addition to adhering to a recommended maintenance schedule.     

Real-time forecasting of water table depth and soil moisture profiles

We present a method for real-time forecasting of water table depth and soil moisture profiles. The method combines a simple form of data-assimilation with a moving window calibration of a deterministic model describing flow in the unsaturated zone and local as well as regional drainage. The local drainage level is calibrated on-line using a moving window calibration. Assigning more weight to the last available measurements then yields a form of model adaptation that is in between on-line calibration and data-assimilation (i.e. a simplified form of Newtonian nudging). Five-day hydrological forecasts are performed based on 5-day weather forecasts, while on-line observations of phreatic level and soil moisture content are assimilated on a daily basis. Advantages of the proposed method are that it improves the real-time forecasts compared to off-line calibration and ordinary moving window calibration and that it yields physically consistent soil moisture profiles.     

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.     

Recharge,upflux and water table response for shallow water table conditions in southwest Florida

A disproportionate increase or decrease in water table in response to minor water input or drainage is observed in shallow water table conditions inside drainage lysimeters. This increase happens because the capillary fringe of the shallow water table reaches up to or near the surface (Wieringermeer effect). The correlations between water table level changes and rainfall, seepage irrigation, drip irrigation, and drainage were analysed. Correlations with rainfall, seepage irrigation, and drainage were high (R² ranged from 0·46 to 0·97). Drip irrigation had low correlations due to the low rates of application (R² ranged from 0·26 to 0·44). Conventional methods of calculating recharge, such as multiplying the specific yield with the water table fluctuations, cannot be used for Wieringermeer effect situations. A method using water balance data and soil moisture at different depths in the lysimeters was developed to estimate recharge and upflux. The recharge results were used to develop the apparent specific yield S_ya, which could be used to calculate consequent recharge events from water table fluctuation data. Combining the water table fluctuation relationships developed with the S_ya value will allow the prediction of recharge from rainfall and irrigation events without the need for soil moisture equipment. Copyright © 2005 John Wiley & Sons, Ltd.     

Tide-induced airflow in a two-layered coastal land with atmospheric pressure fluctuations

Tide-induced airflow is commonly seen in coastal lands and affects ground stability especially with a less permeable pavement on the ground surface. A tide-induced airflow model in a two-layered unsaturated zone consisting of a highly permeable layer underneath a less permeable layer was established by Li and Jiao [Li HL, JJ Jiao. One-dimensional airflow in unsaturated zone induced by periodic water table fluctuation. Water Resour Res 2005;41:W04007. doi:10.1029/2004WR003916] to describe the one-dimensional airflow with constant atmospheric pressure at the ground surface. In this study, we expand the Li and Jiao model by considering the realistic atmospheric pressure fluctuations and the initial condition. A new transient solution to the airflow model is developed for an initial boundary value problem (IBVP). The transient solution can be used not only to calculate the subsurface air pressure at a future time with a known initial condition, but also to evaluate the asymptotic air pressure variations when time becomes long. The amplitude ratio and phase lag of the subsurface air pressure relative to the tide-induced hydraulic head variations inside the unconfined aquifer below the unsaturated zone are investigated. The results reveal that effect on the subsurface pressure due to changes of atmospheric pressure amplitude depends on the configurations of air resistance in the less permeable layer and the air-filled porosity difference in the two layers. The introduction of atmospheric pressure fluctuations into the airflow model leads to insignificant influence on water table level. A field application of the new solution at Hong Kong International Airport in Hong Kong, China is demonstrated. It indicates that the new transient solution can be conveniently used to evaluate the subsurface air pressure with discrete atmospheric pressure data at the ground surface.     

Vadose Zone Monitoring: An Early Warning System

Currently, vadose zone monitoring is required under the Resource Conservation and Recovery Act (RCRA) only at land treatment facilities. Contaminant leak detection through ground water monitoring is very important, but it is considered to be after the fact. Remedial action costs can be reduced considerably by monitoring the vadose zone for compounds that exhibit high rates of movement. Volatile organic compounds (VOCs) exhibit this property and are present at many municipal landfills, recycling facilities, and treatment storage and disposal facilities (TSDFs). Through the authors'personal experience, it has been noted that gaseous phase transport of VOCs through the vadose zone is at least an order of magnitude greater than advective transport of VOCs in ground water. Therefore, VOCs in soil gas are an effective early warning leak detection parameter. Downward movement of leachate can be intercepted by porous cup lysimeters. Attenuation in the vadose zone slows the apparent movement of contaminants; however, it is only a matter of time before leachate reaches the water table. The authors believe that soil-gas and pore-water monitoring should and eventually will be required at all RCRA sites. If vadose zone monitoring becomes an additional requirement under RCRA, both the facility owner and the taxpayer will benefit. During the interim, facility owners can benefit by employing vadose zone monitoring techniques coupled with either qualitative or quantitative chemical analyses.