Skip to main content
Log in

Parametric and numerical analysis of the estimation of groundwater recharge from water-table fluctuations in heterogeneous unconfined aquifers

Analyse paramétrique et numérique de l’évaluation de la recharge des eaux souterraines à partir des fluctuations de nappe dans des aquifères libres hétérogènes

Análisis paramétrico y numérico de la estimación de la recarga subterránea a partir de las fluctuaciones de la superficie freática en acuíferos libres heterogéneos

非均质非承压含水层地下水位波动地下水补给估算的参数和数值分析

Análise paramétrica e numérica da estimativa de recarga das águas subterrâneas a partir de flutuações do nível freático em aquíferos livres heterogêneos

  • Paper
  • Published:
Hydrogeology Journal Aims and scope Submit manuscript

Abstract

The groundwater recharge produced by discrete precipitation events in unconfined aquifers is often estimated from water-table fluctuations (WTFs) recorded in shallow wells. This recharge estimate is prone to uncertainties when recharge is not instantaneous, when there is groundwater drainage, and when there are other processes producing WTFs. A numerical analysis of these uncertainties is presented, which accounts for noninstantaneous recharge and the changes in the stage of a river connected to the unconfined aquifer. This analysis is performed for idealized synthetic unconfined aquifers with one-dimensional (1-D) and 2-D numerical flow models which account for the anisotropy and the spatial heterogeneity of the hydraulic conductivity (K). The logarithm of K is assumed to be a Gaussian random field with a spherical semivariogram. Groundwater recharge may be grossly underestimated with the WTF data when the recharge is not instantaneous. Estimation errors are especially important near the river. On the other hand, the recharge may be largely overestimated when the river stage rises simultaneously during the recharge episode. The errors increase with the variance of the Ln K value and depend on the main direction of anisotropy and the spatial connectivity of the most permeable areas near the river. The errors are large along the most permeable zones connected to the river in the main direction of anisotropy. The recharge estimation errors are largest when the main direction of anisotropy is perpendicular to the river and are smallest when the main direction of anisotropy is parallel to the river.

Résumé

La recharge des eaux souterraines produite par des événements discrets de précipitation dans des aquifères libres est souvent estimée à partir des fluctuations de niveau de nappe (WTFs) enregistrées dans les puits peu profonds. Cette évaluation de recharge est encline aux incertitudes lorsque la recharge n’est pas instantanée, lorsqu’il y a drainage d’eaux souterraines, et lorsqu’il y a d’autres processus produisant des fluctuations de nappe. Une analyse numérique de ces incertitudes est présentée, qui prend en compte une recharge non-instantanée et les changements de niveau d’un cours d’eau en connexion avec l’aquifère libre. Cette analyse est réalisée pour un aquifère libre synthétique idéalisé, à l’aide de modèles numériques unidimensionnels (1-D) et 2-D d’écoulement prenant en compte l’anisotropie et l’hétérogénéité spatiale de la conductivité hydraulique K. On pose l’hypothèse que le logarithme de K est un champ aléatoire gaussien avec un semi-variogramme sphérique. La recharge des eaux souterraines peut être largement sous-estimée avec les données de niveau de nappe quand la recharge n’est pas instantanée. Les erreurs d’évaluation sont particulièrement importantes à proximité du cours d’eau. D’autre part, la recharge peut-être largement surestimée quand le niveau du cours d’eau monte simultanément pendant l’épisode de recharge. Les erreurs augmentent avec la variance de la valeur du Ln K et dépendent de la direction principale de l’anisotropie et de la connectivité spatiale des zones les plus perméables à proximité du cours d’eau. Les erreurs sont grandes le long des zones les plus perméables en relation avec le cours d’eau, dans la direction principale de l’anisotropie. Les erreurs d’évaluation de la recharge sont les plus grandes quand la direction principale de l’anisotropie est perpendiculaire au cours d’eau et les plus petites quand la direction principale de l’anisotropie est parallèle au cours d’eau.

Resumen

La recarga subterránea que se produce en acuíferos libres en episodios discretos de precipitación se suele estimar a partir de las fluctuaciones de la superficie freática (FSF) medidas en pozos someros. La recarga que se obtiene con este método tiene incertidumbres cuando la recarga no es instantánea, existe drenaje subterráneo e intervienen otros procesos que producen fluctuaciones de la superficie freática. En este artículo se presenta un análisis numérico de estas incertidumbres que tiene en cuenta la existencia de una recarga no instantánea y los cambios en el nivel de un río conectado al acuífero libre. El análisis se realiza para acuíferos libres, sintéticos e idealizados mediante modelos numéricos en una (1-D) y dos (2-D) dimensiones que tienen en cuenta la anisotropía y la heterogeneidad espacial de la conductividad hidráulica, K. Se supone que el logaritmo de K es una función aleatoria Gaussiana con un semivariograma esférico. La recarga subterránea puede ser ampliamente subestimada con el método FSF cuando la recarga no es instantánea. Los errores de estimación son especialmente importantes cerca del río. Por otro lado, la recarga puede ser sobreestimada en gran medida si el nivel del río asciende durante el episodio de recarga. Los errores de estimación aumentan con la varianza del Ln K y dependen de la dirección principal de anisotropía y de la conectividad espacial de las zonas más permeables cerca del río. Los errores son grandes a lo largo de las zonas más permeables conectadas al río a lo largo de la dirección principal de anisotropía. Los errores en la estimación de la recarga son máximos cuando la dirección principal de anisotropía es perpendicular al río y son mínimos cuando la dirección principal de anisotropía es paralela al río.

摘要

非承压含水层中离散降水事件产生的地下水补给往往是根据浅井记录的地下水位波动(WTF)来估算的。当补给不是瞬时的,当有地下水排水时,以及当有其他过程产生地下水位波动时,这种补给估计很容易产生不确定性。对这些不确定因素进行了数值分析,解释了非瞬时补给和与无承压含水层相连的河流阶段的变化。本文采用一维(1-D)和二维数值流模型,考虑了导水率的各向异性和空间异质性,对理想的合成无侧限含水层进行了数值分析。K的对数是一个具有球面半变异函数的高斯随机场。当补给不是瞬时,WTF数据可能会严重低估地下水补给。估计误差在河流附近尤为重要。另一方面,在补给过程中,当河段同时上升时,补给可能被大大高估。误差随LnK值的变化而增大,主要取决于各向异性的主方向和河流附近最易渗透的区域的空间连通性。沿主要各向异性方向与河流相连的最渗透带的误差较大。当主各向异性方向垂直于河流时,补给量估计误差最大,当主各向异性方向与河流平行时,计算误差最小。

Resumo

A recarga das águas subterrâneas produzida por eventos pontuais de precipitação em aquíferos livres é frequentemente estimada a partir das flutuações do nível freático (water table fluctuations - WTFs) registadas em poços pouco profundos. Essa estimativa de recarga está sujeita a incertezas quando a recarga não é instantânea, quando há drenagem de águas subterrâneas e quando há outros processos que produzem flutuações no nível freático. É apresentada uma análise numérica dessas incertezas considerando a recarga não instantânea e as mudanças no nível de um rio conectado ao aquífero livre. Esta análise é realizada para aquíferos livres sintéticos idealizados com modelos numéricos de fluxo unidimensional (1-D) e 2-D, que representam a anisotropia e a heterogeneidade espacial da condutividade hidráulica, K. Presume-se que o logaritmo de K seja um campo aleatório gaussiano com um semivariograma esférico. A recarga das águas subterrâneas pode ser subestimada grosseiramente com os dados da WTF quando a recarga não é instantânea. Erros de estimativa são especialmente importantes perto do rio. Por outro lado, a recarga pode ser superestimada quando o nível do rio sobe simultaneamente durante o episódio de recarga. Os erros aumentam com a variância do valor de Ln K e dependem da direção principal da anisotropia e da conectividade espacial das áreas mais permeáveis próximas ao rio. Os erros são grandes ao longo das zonas mais permeáveis conectadas ao rio na direção principal da anisotropia. Os erros de estimativa de recarga são maiores quando a direção principal da anisotropia é perpendicular ao rio e são menores quando a direção principal da anisotropia é paralela ao rio.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12

Similar content being viewed by others

References

Download references

Acknowledgements

The comments and corrections of the two anonymous reviewers are greatly acknowledged.

Funding

The research leading to this work has received funding from ENRESA, the Spanish Ministry of Economy and Competitiveness (Project CGL2016-78281), the FEDER funds and the Galician Regional Government (Ref: ED431C 2017/67 from “Consolidación e estruturación de unidades de investigación competitivas”). The first author had a contract from the FPI Program of the Spanish Ministry of Economy and Competitiveness.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Javier Samper.

Electronic supplementary material

ESM 1

(PDF 2931 kb)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Águila, J.F., Samper, J. & Pisani, B. Parametric and numerical analysis of the estimation of groundwater recharge from water-table fluctuations in heterogeneous unconfined aquifers. Hydrogeol J 27, 1309–1328 (2019). https://doi.org/10.1007/s10040-018-1908-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10040-018-1908-x

Keywords

Navigation