<p>I'm interested in automating spatial interpolation of groundwater level measurements. I'm wondering if anyone knows a method for automatically taking into account things like dry wells, where I can't get a depth to water measurement, but where I do know that the depth to water is greater than the depth of the screened interval of the well.</p>

<p>The interpolation would be constrained at that location so that the water table would be below the bottom of the well. In addition to saying things like "<em>the value at $(x_1,y_1) = z_1$</em>" I could also specify "<em>the value at $(x_2, y_2) &lt; z_2$</em>" or "<em>at $(x_3, y_3) &gt; z_3$</em>".</p>

<p>I'd like to know if anyone knows any interpolation algorithms that can do this, and if there are any tools out there to implement them.</p>

CC BY-SA 3.0
2016-02-05T17:25:09.263Z
It isn't clear what you need help with, are you looking for code that is already developed to manage this, or are you asking for help understanding what the algorithm should include? Thinking in terms of elevations will help, one measures depth to water, but the data that gets utilized is elevation, so dry wells will have a groundwater elevation of <z2 as you talk about it in your question, with z2 being the bottom of the screen interval.
blaydes
CC BY-SA 3.0
2016-02-05T17:45:45.687Z
Right - I do convert to altitude before interpolating. If there is code for this already, of course that would be nice, but any ideas on how to integrate that kind of information into an interpolation algorithm would be welcome, and I could work on the code myself if there aren't already any tools out there.
bklag
CC BY-SA 3.0
2016-02-11T02:31:58.977Z
Automated mapping of the water table is a difficult task. Probably a bigger constraint that should be applied is that it can't rise above the land surface. If you have water level trends in the wells that have gone dry, they can be used with climate data to estimate the current water level. If you can estimate the range of probable level with it truncated by the well depth then you could probably try some kind of stochastic solution.
haresfur
CC BY-SA 3.0
2016-02-11T19:52:24.363Z
You're right, keeping it below land surface where there are no surface-water bodies is one of the big challenges to automating it. The thing is, if I was doing it manually I would be making certain assumptions (like adjusting things in areas where the initial interpolation went over land surface to get the water table below) anyway. I guess I just need a systematic way to apply those assumptions that I could program in as pre- and/or post-processing steps, instead of trying to find an interpolator that has that kind of thing as part of the algorithm.
bklag
CC BY-SA 3.0
2016-03-09T18:20:04.527Z
Are you still working on this? I have some suggestions that could help direct your decisions, unfortunately in the end you will need to make a decision about how you correct for any apparent errors that appear and present a justification for what you decide.
blaydes
<geology><hydrology><water><groundwater><hydrogeology>
How to integrate dry wells into groundwater level mapping automatically (inequality constraints on the interpolation)?
bklag
# How to integrate dry wells into groundwater level mapping automatically (inequality constraints on the interpolation)?

(Asked by: username_0 on 2016-02-04)

I'm interested in automating spatial interpolation of groundwater level measurements. I'm wondering if anyone knows a method for automatically taking into account things like dry wells, where I can't get a depth to water measurement, but where I do know that the depth to water is greater than the depth of the screened interval of the well.

The interpolation would be constrained at that location so that the water table would be below the bottom of the well. In addition to saying things like "*the value at $(x\_1,y\_1) = z\_1$*" I could also specify "*the value at $(x\_2, y\_2) < z\_2$*" or "*at $(x\_3, y\_3) > z\_3$*".

I'd like to know if anyone knows any interpolation algorithms that can do this, and if there are any tools out there to implement them.

**Tags:** geology, hydrology, water, groundwater, hydrogeology

> **Comments:**
> **username_1**: It isn't clear what you need help with, are you looking for code that is already developed to manage this, or are you asking for help understanding what the algorithm should include? Thinking in terms of elevations will help, one measures depth to water, but the data that gets utilized is elevation, so dry wells will have a groundwater elevation of <z2 as you talk about it in your question, with z2 being the bottom of the screen interval.  
> **username_0**: Right - I do convert to altitude before interpolating. If there is code for this already, of course that would be nice, but any ideas on how to integrate that kind of information into an interpolation algorithm would be welcome, and I could work on the code myself if there aren't already any tools out there.  
> **username_2**: Automated mapping of the water table is a difficult task. Probably a bigger constraint that should be applied is that it can't rise above the land surface. If you have water level trends in the wells that have gone dry, they can be used with climate data to estimate the current water level. If you can estimate the range of probable level with it truncated by the well depth then you could probably try some kind of stochastic solution.  
> **username_0**: You're right, keeping it below land surface where there are no surface-water bodies is one of the big challenges to automating it. The thing is, if I was doing it manually I would be making certain assumptions (like adjusting things in areas where the initial interpolation went over land surface to get the water table below) anyway. I guess I just need a systematic way to apply those assumptions that I could program in as pre- and/or post-processing steps, instead of trying to find an interpolator that has that kind of thing as part of the algorithm.  
> **username_1**: Are you still working on this? I have some suggestions that could help direct your decisions, unfortunately in the end you will need to make a decision about how you correct for any apparent errors that appear and present a justification for what you decide.  

---

## 0 answers

