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

Águila, Jesús F.; Samper, Javier; Pisani, Bruno

doi:10.1007/s10040-018-1908-x

Cited by 18 publications

(4 citation statements)

References 75 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The flow model was performed with CORE 2D V5, a finite element code for transient saturated and unsaturated water flow, heat transport, and multicomponent reactive solute transport in heterogeneous and anisotropic media [29]. CORE 2D V5 has been extensively verified against analytical solutions and other reactive transport codes [30] and widely used to model groundwater flow and solute transport in real-world aquifers [31], laboratory and in situ experiments (pp. 90-91), and the long-term geochemical evolution of radioactive waste repositories.…”

Section: Computer Codementioning

confidence: 99%

Groundwater Flow Model along a Vertical Profile of the Sardas Landfill in Sabiñánigo, Huesca, Spain

Samper,

Sobral,

Pisani

et al. 2023

Water

Self Cite

View full text Add to dashboard Cite

Lindane (a harmful contaminant) was produced in Sabiñánigo (Huesca, Spain) and deposited at the Sardas landfill. This site contains a large mass of pollutants, which have an extremely large contamination potential of the Ebre River. The site has undergone numerous human interventions that have modified the natural conditions. The site exhibits complex hydrogeological patterns and has been monitored systematically for a long period of time, and a large amount of geological, hydrological, and hydrogeological data are available. Here, a 2D finite element groundwater flow model along a vertical profile heading east–west along the thalweg of the former gully is presented. The main goal is modelling groundwater flow through the landfill, the Gállego River alluvial aquifer, and its interactions with the Sabiñánigo reservoir. The numerical model confirms the prevailing conceptual hydrogeological model of the site. The main results include: (1) Groundwater flows into the landfill mainly along perimeter ditches, which do not properly drain the surface and subsurface runoff (13.84 m3/d) and from the underlying marly rock (8.84 m3/d); (2) The total landfill leachate outflow towards the alluvial floodplain underneath the front slurry wall through a shallow marl layer is equal to 17 m3/d; (3) The oscillations of the Sabiñánigo reservoir water level produce a tidal effect that results in periodic changes of the hydraulic gradient between the alluvial gravels and the reservoir; (4) Groundwater flows generally from the alluvial aquifer towards the reservoir in an average E–W direction with an average Darcy velocity equal to 5 cm/d. The flow direction, however, changes to W–E when the reservoir level rises suddenly and; (5) The hydrodynamic parameters of the alluvial silts and reservoir silting sediments are crucial in determining the influence radius of the inversion of groundwater flow direction when the reservoir level rises suddenly. Model results enhance the confidence of the conceptual model, provide the basis for detailed specific models of the landfill and the alluvial aquifer, and highlight the importance of considering the tidal effect of the reservoir level oscillations. They also provide valuable information for managing the landfill and its impact on the surrounding groundwater system.

show abstract

Section: Computer Codementioning

confidence: 99%

Groundwater Flow Model along a Vertical Profile of the Sardas Landfill in Sabiñánigo, Huesca, Spain

Samper,

Sobral,

Pisani

et al. 2023

Water

Self Cite

View full text Add to dashboard Cite

show abstract

“…This is a code for transient saturated and unsaturated water flow, heat transport, and multicomponent reactive solute transport under both local chemical equilibrium and kinetic conditions in heterogeneous and anisotropic media. CORE 2D V5 has been extensively verified against analytical solutions and other reactive transport codes [30,86,87] and is widely used to model groundwater flow and solute transport in real-world aquifers [80,88,89] and laboratory or in situ experiments [90,91], as well as to study the long-term geochemical evolution of radioactive waste repositories [92][93][94]. Moreover, the code is able to consider the porosity feedback effect caused by mineral dissolution/precipitation reactions and update the flow and transport parameters at every time step, thus enabling the investigation of the spatial and temporal evolution of fractures.…”

Section: Numerical Modelmentioning

confidence: 99%

Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media

et al. 2022

Self Cite

View full text Add to dashboard Cite

Gypsum dissolution leads to the development of karstic features within much shorter timescales than in other sedimentary rocks, potentially leading to rapid deterioration of groundwater quality and increasing the risk of catastrophes caused by subsidence. Here, we present a 2-D reactive transport model to evaluate gypsum karstification in physically and chemically heterogeneous systems. The model considers a low-permeability rock matrix composed mainly of gypsum and a discontinuity (fracture), which acts as a preferential water pathway. Several scenarios are analyzed and simulated to investigate the relevance for gypsum karstification of: (1) the dynamic update of flow and transport parameters due to porosity changes; (2) the spatial distribution of minerals in the rock matrix; (3) the time evolution of water inflows through the boundaries of the model; (4) the functions relating permeability, k, to porosity, ϕ. The average porosity of the matrix after 1000 years of simulation increases from 0.045 to 0.29 when flow, transport, and chemical parameters and the water inflows through the boundary are dynamically updated according to the porosity changes. On the contrary, the porosity of the matrix hardly changes when the porosity feedback effect is not considered, while its average increases to 0.13 if the water inflow occurs through the discontinuity. Moreover, the dissolution of small amounts of highly soluble sulfate minerals plays a major role in the development of additional fractures. The increase in hydraulic conductivity is largest for the power law with an exponent of n = 5, as well as the Kozeny-Carman and the modified Fair-atch k-ϕ relationships. The gypsum dissolution front propagates into the matrix faster when the power law with n = 2 and 3 and the Verma–Pruess k-ϕ relationships are used.

show abstract

“…The groundwater flow numerical model was performed with CORE 2D (Molinero et al 2004;Samper et al 2009Samper et al , 2012Lu et al 2011;Águila et al 2019;Naves et al 2019). The model was calibrated in two stages.…”

Section: Groundwater Flow Numerical Modelmentioning

confidence: 99%

Hydrogeology and groundwater management in a coastal granitic area with steep slopes in Galicia (Spain)

et al. 2021

Self Cite

View full text Add to dashboard Cite

Groundwater availability, management and protection are great challenges for the sustainability of groundwater resources in the scattered rural areas of the Atlantic regions of Europe where groundwater is the only option for water supply. This report presents a hydrogeological study of the coastal granitic area of Oia in northwestern Spain, which has unique geomorphological and hydrogeological features with steep slopes favoring the erosion of the weathered granite. The hydrogeological conceptual model of the study area includes: (1) the regolith layer, which is present only in the flat summit of the mountains; (2) the slope debris and the colluvial deposits, which are present in the intermediate and lowest parts of the hillside; (3) the marine terrace; and (4) the underlying fractured granite. Groundwater recharge from rainfall infiltration varies spatially due to variations in terrain slope, geology and land use. The mean annual recharge estimated with a hydrological water balance model ranges from 75 mm in the steepest zone to 135 mm in the lowest flat areas. Groundwater flows mostly through the regolith and the detrital formations, which have the largest hydraulic conductivities. Groundwater discharges in seepage areas, springs, along the main creeks and into the sea. The conceptual hydrogeological model has been implemented in a groundwater flow model, which later has been used to select the best pumping scenario. Model results show that the future water needs for domestic and tourist water supply can be safely provided with eight pumping wells with a maximum pumping rate of 700 m3/day.

show abstract

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

Cited by 18 publications

References 75 publications

Groundwater Flow Model along a Vertical Profile of the Sardas Landfill in Sabiñánigo, Huesca, Spain

Groundwater Flow Model along a Vertical Profile of the Sardas Landfill in Sabiñánigo, Huesca, Spain

Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media

Hydrogeology and groundwater management in a coastal granitic area with steep slopes in Galicia (Spain)

Contact Info

Product

Resources

About