Hello all, 
My question is about how can I model the effect of groundwater abstraction on river flow, in fact, how can I make Fluid-Transfer BC become the calculated variable. 
I have an unconfined, phreatic type aquifer. It is highly connected to the river and recharge comes mostly from it. In a future scenario I have some pumping wells and I want to know how much water can I extract without draining the river. 
It would be very helpful if someone could provide me a paper with a similar problem. 
Thanks, 
Idoia.