-
Notifications
You must be signed in to change notification settings - Fork 115
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Proposal: optional exfiltration conductance for RIV, GHB, LAK & SFR #419
Comments
Are there any publications on this that you could reference here? |
Id second this suggestion as a helpful addition that we would make use of. As for reference, the earliest that aware of would be K.R.Rushton & L.M.Tomlinson, 1979, Possible mechanisms for leakage between aquifers and rivers. Journal of Hydrology Volume 40, Issues 1–2, Pages 49-65 https://doi.org/10.1016/0022-1694(79)90087-8. Figure 2b in this paper is what I believe is being proposed (2a and 2c already possible in RIV and SFR) but may be other options also. We would use this functionality as commonly desired to have easier inflow to rivers where flows come through base and sides. Where as in ex-filtration flow is generally down through bed sediment only. Thanks |
Hi Chris, |
I can second this request. It's seems to be relatively common in groundwater modelling in the Netherlands to make a distinction between between infiltration and exfiltration. (Among the Deltares changes for iMODFLOW, an "infiltration factor" for the river package has been added.) Unfortunately, there appears to be very little scientific literature on this subject. I believe there are at least two physical processes at play for exfiltration versus infiltration:
The second process can lead to very significant resistance build-up. This is quite clear for the surface water distribution network in Dutch polders: in these cases, the surface water stage can be meters higher than the surrounding land (behind a levee). Yet nearby piezometers show no influence from the surface water, suggesting resistances in the order of thousands of days and very little infiltration. (I wonder whether in exfiltration situation some "hydraulic failure" occurs in the bed, opening the bed again.) The infiltration systems for drinkwater production in the Dutch coastal dunes is another good example, and includes actual measurements of fluxes. Pre-filtrated river water is infiltrated in (large) ponds, flows through the phreatic aquifer, and is abstracted via horizontal drains, making sure there's enough travel time to rid the water of pathogens. Though the water is pre-filtrated, a "cake" builds and infiltration fluxes measurably decrease over the years, with a large fraction of the water infiltrating via the sides of the ponds after some time. Unsurprisingly, the higher the head difference between pond and drain, the more frequent these ponds get clogged and to be emptied for maintenance. So there has been some research, but it's of a relatively practical nature, with reports in Dutch. Here are some sources from a quick search on Google scholar, by checking who cites the paper @AlastairBlack mentions: Transient Leakance and Infiltration Characteristics during Lake Bank Filtration Field evidence of a dynamic leakage coefficient for modelling river–aquifer interactions This article from my colleague Joost actually contains exactly the plot you'd want... and shows no difference between infiltration and exfiltration over a range of ~1.4 m! (Edited:) Viscosity plays something of a role, but because the measurements are for a single season, there won't be large temperature and viscosity variation. Investigating summer flow paths in a Dutch agricultural field using high To end on a practical note, for basic rivers, the behavior can be relatively easily achieved in current modflow6 by adding drainage on top of the river with a drainage elevation equal to the river stage. Maybe you could do the same with the other packages by making use of the water mover, but it seems like an awful hassle. |
FYI, for those who are trying to implement this in their models now: Zaadnoordijk, W. (2009). Simulating Piecewise-Linear Surface Water and Ground Water Interactions with MODFLOW. Ground Water. |
Hi again,
this is meant as a proposal for introducing an optional exfiltration conductance for the above mentioned boundary conditions. We often have the application case, that infiltration conductance (head/stage of boundary condition > groundwater head) is differnt from exfiltration conductance (HBound<Head) and currently there is no way in MF6 to simulate this behaviour.
My favourite solution would be a keyword in OPTIONS specifying, whether exfiltration conductance is read from a separate column in the PERIOD data.
Another simpler but still sufficient way would be to specify just an exfiltration factor in OPTIONS applied to the conductance when exfiltration takes place.
What do you think about it?
Thanks & best regards
Andreas
The text was updated successfully, but these errors were encountered: