Application of the water table fluctuation method to characterize groundwater recharge in the Pampa plain,Argentina

Abstract

The water table fluctuation (WTF) method is based on accepting that rises of a water table are due to recharge water reaching the groundwater. To apply the method, an estimate of the specific yield of the zone of fluctuation of the groundwater level is required. In this paper, a method for estimation of the specific yield (S_y) is proposed; it consists of a graphical procedure which relates rises in groundwater level to the precipitation from which they originated. The method presents more reliable S_y values as the number of events measured increases. Eighteen years of daily measurements were analysed to obtain a S_y value of 0.09, which was used to apply the WTF method. The obtained recharge values show consistency with values calculated by other authors for the same region.

Editor D. Koutsoyiannis

Citation Varni, M., Comas, R., Weinzettel, P., and Dietrich, S., 2013. Application of water table fluctuation method to characterize the groundwater recharge in the Pampa plain, Argentina. Hydrological Sciences Journal, 58 (7), 1445–1455.     

Effect of water table drawdown on peatland dissolved organic carbon export and dynamics

Peatlands play an important role in the global carbon cycle, and loss of dissolved organic carbon (DOC) has been shown to be important for peatland carbon budgets. The objective of this study was to determine how net production and export of DOC from a northern peatland may be affected by disturbance such as drainage and climate change. The study was conducted at a poor fen containing several pool–ridge complexes: (1) control site–no water table manipulation; (2) experimental site–monitored for one season in a natural state and then subjected to a water table drawdown for 3 years; (3) drained site–subjected to a water table drawdown 9 years prior to monitoring. The DOC concentration was measured in pore water along a microtopographic gradient at each site (hummock, lawn and hollow), in standing water in pools, and in discharge from the experimental and drained sites. The initial water table drawdown released ～3 g of carbon per square metre in the form of DOC, providing a large pulse of DOC to downstream ecosystems. This value, however, represents only 1–9% of ecosystem respiration at this site. Seasonal losses of DOC following drainage were 8–11 g of carbon per square metre, representing ～17% of the total carbon exchange at the experimental study site. Immediately following water table drawdown, DOC concentrations were elevated in pore water and open water pools. In subsequent seasons, DOC concentration in the pool declined, but remained higher than the control site even 11 years after water‐table drawdown. This suggests continued elevated net DOC production under lower water table conditions likely related to an increase in vegetation biomass and larger water table fluctuations at the experimental and drained sites. However, the increase in concentration was limited to initially wet microforms (lawns and hollows) reflecting differences in vegetation community changes, water table and soil subsidence along the microtopographic gradient. Copyright © 2008 John Wiley & Sons, Ltd and Her Majesty the Queen in right of Canada.     

A hydrographic method to identify the groundwater net recharge,barometric effect and evapotranspiration from a complicated semidiurnal water table fluctuation

A hydrographic method was proposed to separate out the hourly scaled groundwater level changes caused by net recharge, barometric effects and evapotranspiration from a semidiurnal water table fluctuation. A characteristic midnight time, with a turning point of the barometric pressure change and high relative air humidity, which meant that neither the barometric effect nor groundwater evapotranspiration occurred, was proposed for determining the net recharge rate r_net. Then, the barometric efficiency f_bar was estimated using the other time period without evapotranspiration, and the evapotranspiration rate r_ETG was finally obtained using the remainder of the water level changes. A case example illustrated that estimation of f_bar using the proposed method was more accurate than that using the traditional error analysis method, which may result in a significant underestimation under the condition of the present water level changes. Additionally, the abnormal semidiurnal fluctuations, more specifically, two step-down fluctuations, which may be a common pattern when the groundwater level is controlled by net recharge, barometric effects and evapotranspiration, can be well understood using the three components separated out. The results also showed that nighttime groundwater evapotranspiration, accounting for an average of 23% of that during the daytime, cannot be ignored.     

One‐dimensional expression to calculate specific yield for shallow groundwater systems with microrelief

Although the importance to account for microrelief in the calculation of specific yields for shallow groundwater systems is well recognized, the microrelief influence is often treated very simplified, which can cause considerable errors. We provide a general one‐dimensional expression that correctly represents the effect of a microrelief on the total specific yield that is composed of the soil and surface specific yield. The one‐dimensional expression can be applied for different soil hydraulic parameterizations and soil surface elevation frequency distributions. Applying different van Genuchten parameters and a simple linear microrelief model, we demonstrate that the specific yield is influenced by the microrelief not only when surface storage directly contributes to specific yield by (partial) inundation but also when water levels are lower than the minimum surface elevation. Compared with a simplified representation of the soil specific yield, in which a mean soil surface is assumed for the calculation of soil specific yield, the correct representation can lead to lower as well as higher soil specific yields depending on the specific interaction of the soil water retention characteristics and the microrelief. The new equation can be used to obtain more accurate evapotranspiration estimates from water level fluctuations and to account for the effect of microtopographic subgrid variability on simulated water levels of spatially distributed hydrological models. Copyright © 2015 The Authors Hydrological Processes Published by John Wiley & Sons Ltd     

Groundwater recharge estimation in arid hardrock‐alluvium aquifers using combined water‐table fluctuation and groundwater balance approaches

This paper proposes an approach to estimate groundwater recharge using an optimization‐based water‐table fluctuation method combined with a groundwater balance model in an arid hardrock‐alluvium region, located at the Oman–United Arab Emirates border. We introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long‐term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm³, respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12‐year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.     

Drivers of peatland water table dynamics in the central Andes,Bolivia and Peru

Cushion plant dominated peatlands are key ecosystems in tropical alpine regions of the Andes in South America. The cushion plants have formed peat bodies over thousands of years that fill many valley bottoms, and the forage produced by the plants is critical for native and nonnative domesticated mammals. The sources and flow paths of water supporting these peatlands remain largely unknown. Some studies have suggested that glacier meltwater streams support some peatlands, and that the ongoing loss of glaciers and their meltwaters could lead to the loss or diminishment of peatlands. We analysed the hydrologic regime of 10 peatlands in four mountain regions of Bolivia and Peru using groundwater monitoring. Groundwater levels in peatlands were relatively stable and within 20 cm of the ground surface during the rainy season, and many sites had water tables 40–90 cm below the ground surface in the dry season. Topographic and groundwater elevations in the peatlands demonstrated that the water source of all 10 peatlands was hillslope groundwater flowing from lateral moraines, talus, colluvium, or bedrock aquifers into the peatlands. There was little to no input from streams, whether derived from glacier melt or other sources, and glacier melt could not have recharged the hillslope aquifers supporting peatlands. We measured the stable water isotopes in water samples taken during different seasons, distributed throughout the catchments, and the values are consistent with this interpretation. Our findings indicate that peatlands in the study region are recharged by hillslope groundwater discharge rather than stream water and may not be as vulnerable to glacial decline as other studies have indicated. However, both glaciers and peatlands are susceptible to changing thermal and precipitation regimes that could affect the persistence of peatlands.     

Moisture dynamics and hydrophysical properties of a transplanted acrotelm on a cutover peatland

The natural carbon storage function of peatland ecosystems can be severely affected by the abandonment of peat extraction, influencing peatland drainage, leading to large and persistent sources of atmospheric CO₂. Moreover, these cutover peatlands have a low and variable water table position and high tension at the surface, creating harsh ecohydrological conditions for vegetation re‐establishment, particularly peat forming Sphagnum moss. Standard restoration techniques aim to restore the peatland to a carbon accumulating system through various water management techniques to improve hydrological conditions and by reintroducing Sphagnum at the surface. However, restoring the hydrology of peatlands can be expensive due to the cost of implementing the various restoration techniques. This study examines a peat extraction‐restoration technique where the acrotelm is preserved and replaced directly on the cutover peat surface. An experimental peatland adopting this acrotelm transplant technique had both a high water table and peat moisture conditions providing sufficient water at the surface for Sphagnum moss. Average water table conditions were higher at the experimental site (?8·4 ± 4·2 cm) compared to an adjacent natural site (?12·7 ± 6·0 cm) suggesting adequate moisture conditions at the restored surface. However, the experimental site experienced high variability in volumetric moisture content (VMC) in the capitula zone (upper 2 cm) where large diurnal changes in VMC (～30%) were observed, suggesting possible disturbance to the peat matrix structure during the extraction‐restoration process. However, soil–water retention analysis and physical peat properties (porosity and bulk density) suggest that no significant differences existed between the natural and experimental sites. Any structural changes within the peat matrix were therefore minimal. Moreover, low soil‐water tensions were maintained well above the laboratory measured critical Sphagnum threshold of 33% (?100 mb) VMC, further indicating favourable conditions for Sphagnum moss survival and growth. Copyright © 2007 John Wiley & Sons, Ltd.     

Water uptake by saltcedar (Tamarix ramosissima) in a desert riparian forest: responses to intra‐annual water table fluctuation

There is considerable interest in naturalizing flow regime on managed rivers to slow the spread of saltcedar (Tamarix ramosissima) invasion in southwestern USA or to preserve riparian forests dominated by saltcedar and other species in northwestern China. However, little is known about the responses of established saltcedar in water sources to frequent intra‐annual fluctuation of water table resulting from this new, more dynamic flow regime. This study investigates how saltcedar at a riparian site in the middle reaches of the Heihe River, northwest China, responds in water sources use to intra‐annual water table fluctuations. Stable oxygen isotope was employed to determine accurate depth at which saltcedar obtains its water supply, and soil moisture monitoring was used to determine sources of plant‐available soil water. We found that the primary zone of water uptake by saltcedar were stable at 25–60 cm depth, but the water sources used by saltcedar switched between groundwater and soil moisture with the water table fluctuations. Saltcedar derived its water from groundwater when water table was at depth less than 60 cm but switched to soil moisture at 25–60 cm depth when water table declined. It is supposed that the well‐developed clay layer at 60–80 cm depth constrained lateral roots of saltcedar to the soil layers above 60 cm, while the fine‐textured soils at this site, which were periodically resaturated by rising groundwater before the stored soil moisture had become depleted, provided an important water reservoir for saltcedar when groundwater dropped below the primary zone of fine roots. The root distribution of saltcedar may also be related to local groundwater history. The quick decline in water table in the early 1980s when the riparian saltcedar had established may strand its roots in the shallow unsaturated zone. We suggested that raising the water table periodically instead of maintaining it invariably above the rooting depth could sustain desired facultative phreatophytes while maximizing water deliveries. Copyright © 2015 John Wiley & Sons, Ltd.     

A new analytical solution of tidal water table fluctuations in a coastal unconfined aquifer

This paper presents a new perturbation solution of the non-linear Boussinesq equation for one-dimensional tidal groundwater flow in a coastal unconfined aquifer. Built upon the work of Parlange et al. [Parlange, J.-Y., Stagnitti, F., Starr, J.L., Braddock, R.D., 1984. Free-surface flow in porous media and periodic solution of the shallow-flow approximation, J. Hydrol., 70, 251–263], the solution adopts a new perturbation parameter that is by definition less than unit, and thus is applicable to a wider range of physical conditions within the constraint of the Boussinesq approximation. This approach avoids a secular term in the third-order perturbation equation of Parlange et al. (1984), enabling the derivation of the third- and higher-order solutions. In comparison with a numerical (“exact”) solution, the new perturbation solution is shown to be slightly more accurate than that of Parlange et al. (1984) with the second-order approximation. The obtained third-order solution exhibits considerable improvement in accuracy. In relatively simple analytical forms, the present perturbation solution will help to understand better the non-linear characteristics of tidal water table fluctuations in as modeled by the non-linear Boussinesq equation coastal unconfined aquifers.     

基于Forward-Modeling方法的黑河流域水储量变化特征研究

黑河流域陆地水储量变化对流域下游等周边区域水资源的合理利用以及经济和社会发展等有着重要的意义.本文利用2003年1月至2013年12月的GRACE RL05数据反演了黑河流域陆地水储量长时间序列的变化,并针对重力场模型和数据处理中产生的信号泄漏问题,采用Forward-Modeling方法进行了改正并恢复泄漏信号;将GRACE获得的泄漏信号恢复前后的黑河流域水储量变化结果与全球水文模型GLDAS和CPC进行比较分析,结果表明泄漏信号改正后的结果与水文模型结果的时间序列相关性均有明显提高,从其空间分布结果可以看出Forward-Modeling方法有效地恢复初始信号、增强被湮没的信号,泄漏信号误差减小;通过分析黑河流域水储量变化的长时间序列结果,发现其具有明显的阶段性变化特征,即2003—2006年呈明显下降趋势,约为-0.86 cm·a-1,在2007—2010年趋于平衡状态,而2011—2013年则呈现缓慢上升趋势约为0.14 cm·a-1;联合GRACE数据和GLDAS数据反演了黑河流域地下水储量变化,并与全球降雨数据GPCC进行了比较分析,两者相关性可达到0.88以上.     

Quantification of the hydraulic diffusivity of a bentonite-sand mixture using the water head decrease measured upon sudden flow interruption

This paper presents a new modelling approach to quantify the hydraulic diffusivity of low-permeability unconsolidated porous media under confined saturated-flow conditions in the laboratory. The derived analytical solution for the transient variation of the hydraulic head after flow interruption was applied to experimental data obtained from continuous measurements of the water pressure at two locations in the soil column. Three soil samples made of a mixture of natural bentonite (at different mass fractions) and medium sand were studied during a series of stepwise constant flow rates of water. The numerical results well fit the experimentally measured decrease of the dimensionless hydraulic head. The study shows that the increase of the mass fraction of bentonite in the soil sample from 10 to 30% is accompanied by a strong decrease of the hydraulic diffusivity from 2.4 × 10⁻² to 1.1 × 10⁻³ m² s⁻¹, which is clearly due to the decrease of the hydraulic conductivity of the soil sample. The specific storages obtained for each of the three samples are in the same order of magnitude and seem to decrease with the increase of mass fraction of bentonite. However, they clearly reflect the predominant portion of the compressibility of the porous medium compared with that of water. Compared with reported literature values for clayey soils, the specific storage values in this study are slightly higher, varying within the range of 2 × 10⁻³ to 8.1 × 10⁻³ m⁻¹._. The experimental results also give insight into the limitations of the modelling approach. In the case of low-permeability soils (K < 2 × 10⁻⁶ ms⁻¹) and steady-flow conditions with low Reynolds numbers, for example, Re < 0.003, it is recommended to choose a time duration for flow interruption between subsequent flow rate steps of longer than 5 s. For high-permeability porous media, to increase the precision of the quantified hydraulic diffusivity, it might be useful to select a measuring frequency significantly higher than 1 Hz.     

Discharge and water table predictions using a generalised TOPMODEL formulation

A simple, generalized saturated zone formulation is presented in this paper to relax the assumption of an exponential function originally made in TOPMODEL. This saturated zone model is based on the concept of a ‘discharge:relative storage’ (QΔS) function which is derived empirically, using recession curve analysis, and may be of arbitrary form. The generalized formulation is applied to the Seternbekken MINIFELT catchment in Norway, where detailed distributed water table data have been measured. These water table data are used to suggest an empirical, power law modification of the topographic a/tan β index. Results for the simulation through time of discharges and water table depths at a few locations show that the generalized saturated zone formulation is as efficient a simulator of the observed data as a conventional TOPMODEL, but requires one parameter less to be calibrated. The simulation of detailed water table distributions is only approximate in both cases. The modified power law index shows only a small improvement but provides a basis for a discussion of possible sources of error in the TOPMODEL assumptions for this site. © 1997 John Wiley & Sons, Ltd.     

A new analytical expression for ultimate specific yield and shallow groundwater drainage

The ultimate specific yield, i.e. the specific yield under steady‐state condition, is generally dependent on the water table depth in the shallow groundwater systems. This paper presents a new expression to describe this relationship. The expression is based on the relationship between soil water content–soil depths below ground surface near the saturation zone, i.e. water content profile, which can be fitted by the van Genuchten model. Because this model is highly nonlinear, a Taylor series expansion was used to solve the integral equation for specific yield in our study. The new expression was verified by a drainage experiment, and the results demonstrated that the new expression is reliable for eolian sand when the range of water table fluctuation is less than 30 cm. The expression is easily applicable because knowledge of the water content profile near the saturation zone is the only requirement. Compared with the expression obtained from the Brooks–Corey model, the new expression can be used for a variety of aquifer and soil media. The improved accuracy of the specific yield provides a better estimate of discharge rates in shallow groundwater systems with water table fluctuations. Copyright © 2014 John Wiley & Sons, Ltd.     

一种通过地形指数提取流域自由水蓄水容量空间分布的方法

对于分布式水文模型而言,如何获得参数的空间分布是模型应用的重点和难点问题.本文将分水源参数中的敏感参数——自由水蓄水容量为研究对象.建立地形指数与自由水蓄水容量的函数关系,以此提取流域内的自由水蓄水容量空间分布.最后利用本方法提取了陕西省陈河流域的自由水蓄水容量空间分布,并将之作为栅格型新安江模型的参数进行洪水模拟演算.应用结果表明本文提出的方法得到了理想的模拟结果.该方法以物理规律为基础能较为准确地计算出流域内自由水蓄水容量的空间分布,为分布式模型的发展奠定了坚实的基础.     

The westerly fluctuation and water vapor transport over the Qilian-Heihe valley

The westerly fluctuation and the atmospheric water vapor transport over the Qilian-Heihe valley are analyzed and the results show that, in the water vapor transport stream field from Jun to September, this valley is in the westerly stream and the water vapor comes from westerlies water transport via the Black Sea and the Caspian Sea. The net water vapor transport is less net import and different from most areas of the northwest China. The interannual changes in water vapor transport over the valley arise from the westerly fluctuation, and have a positive relationship to the interannual changes in westerly wind speed. The cold air actions from the Mongol low pressure are the primary system that controls the westerly water vapor transport. Its action chain is that, the Mongol low pressure is strengthened → the circulation meridionality will be increased → the cold air will move southwards → the westerly will be stronger → the wind convergence of direction and speed will be stronger → the water vapor convergence transport will be increased → the local water vapor content will be increased. The interannual changes in atmospheric water vapor transport over the valley rely mainly on the convergence transport, but the effect of advection transport is less. The interannual changes of strong or weak westerly affect mainly the convergence transport, and then make the atmospheric water vapor net transport increase or decrease over the Qilian-Heihe valley.     

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the contrary, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending on parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). This effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance. Published 2020. This article is a U.S. Government work and is in the public domain in the USA     

Assessment of the impact of pools on the water isotopic signature of a boreal patterned peatland

Recent studies have shown that boreal peatlands exhibit considerable chemical variability but without clear spatial pattern. This chemical heterogeneity illustrates the complex hydrological behaviour of peatlands, particularly patterned fen. Isotopic, chemical and physical tracers were used to describe the hydrological behaviour of a small boreal headwater catchment (13 ha) during the snow‐free period with a special emphasis on the downstream patterned fen. Results showed that shallow pools were mixed every day during the summer, particularly during nights or discharge periods. Despite large water storage capacities in pools, which should induce large buffer effect, hydrological behaviour of patterned fen is more similar to a piston flow process. This is probably because of the division of the fen into successive small cascading streamflow reservoirs. The consequences were a rapid change of the chemical signature throughout the fen, particularly upstream. A spatial pattern was observed downstream in early summer. The isotopic signature passed from an upstream depleted and homogeneous signature to a progressively enriched downstream signature. However, this pattern was not identified during the wetter period (late summer), probably because the discharge, which dominated the water budget, decreased the surface water residence time and flushed a large proportion of stored surface water. We developed for this patterned fen a conceptual model of the surface flow to explain these particular mixing effects and the implications on the dynamics of the chemical signature. To further our understanding of similar boreal headwater catchments, future work should include the development of a multiple mixed‐reservoir model. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental study on the aseismic capacity of a wooden house using shaking table

The aseismic capacity of a typical Korean wooden house built using traditional construction methods is quantitatively estimated. Tenon joints were used in wooden frames. Two 1:4 scale models were tested for rock and soil foundation conditions. Scaled real earthquake time histories were input for the tests. The natural frequency and modal damping ratio of the wooden house in the elastic range were 1·66 Hz and 7 per cent, respectively. The Peak Ground Acceleration (PGA) at the collapse of the house at the soil site was 0·25g, whereas PGA for moderate damage at the rock site was 0·6g. A significant reduction in acceleration response and increase in displacement response was observed for rock and soil foundation conditions, respectively. The wooden house studied is much more vulnerable at soil sites than at rock sites due to the rich low-frequency contents of the input motion and the flexible characteristics of the wooden house. Non-linear dynamic analyses using the modified Double-Target model were compared with test results. The modified Double-Target model appropriately simulates the non-linear inelastic behaviour of a wooden house with tenon joints. Copyright © 1999 John Wiley & Sons, Ltd.     

Modelling of benzene distribution in the subsurface of an abandoned gas plant site after a long term of groundwater table fluctuation

The non‐aqueous phase liquid simulator was used to model and interpret the occurrence of a thin benzene‐contaminated soil layer 9.0 m below the groundwater table in an abandoned gas plant site. The simulator was first evaluated in column tests under similar conditions to the contaminated site. Saturation–capillary pressure (S–P) relationships were extended from the laboratory scale of the column tests to the field scale of the subsurface at the abandoned site. Dynamic boundary conditions were established in order to prevent the model from generating excessive vertical velocities. The modelled benzene layer formation process agreed well with the in situ observations. With falling and then rising of the water table, benzene release from the surface migrated downward and then upward and distributed itself below and above the water table. Biochemical degradation of benzene made the distribution discontinuous in the subsurface. These two factors resulted in the thin benzene‐contaminated layer below the groundwater table. Copyright © 2012 John Wiley & Sons, Ltd.     

高邮车逻井水位大幅波动异常分析

2016—2019年高邮车逻井水位出现多次大幅波动变化。结合周边自然环境和农业活动的调查工作,对异常成因和影响水位变化的因素进行了跟踪分析。研究认为,该井水位出现的大幅波动变化异常是高邮湖水位短期大幅波动产生的荷载影响所致,与水利部门对高邮湖水位调控间关系密切,而非地震地球物理异常。