DOV-Vlaanderen / groundwater-logger-validation

Analysis on validation methods for groundwater logger data
MIT License
2 stars 2 forks source link

Drifts 02: Analysis of differences between one barometer and many others #48

Open DavorJ opened 4 years ago

DavorJ commented 4 years ago

The main idea for drift detection is to compare one barometer with some other in its vicinity, and to look at how they differ. If the mean difference after some time is large, then this might indicate a drift. For non-drifting barometers we expect white noise around the mean (which is primarily determined by the baro difference in altitude):

image

Although the two barometers are 126km apart, their differences are nicely constant over the years, which suggests that either both are drifting the same way, or that none is drifting. Here and here is an other nice one.

On the other hand, if one is drifting and the other is not, then we would get a clear slope in the mean over time like in the following image: image

Sometimes the difference patterns seem cyclic, which I do not understand. I have made a separate post for this here: #49. Anyone a clue? Anyway, there are 4428 unique informative plots based on pairwise comparisons like this. (In total there are many more -- some 200 x 200 barometer place comparisons -- but their timestamps do not always match, they do not overlap, have more than 10 % outliers based on gwloggeR::detect_outliers(), etc. So there is quite some data wrangling to be done to compare them correctly.)

Sometimes, they both seem to drift in opposite directions and thus the difference effect is reinforced:

image

And if they would drift in the same direction, we wouldn't even see a difference. I haven't found a telling example of this kind.

Anyway, pairwise comparisons are thus problematic because:

To remove these problems we compare one baro with many others. The idea is that even is one or two of the many are drifting, they would not obscure our inference about drift.

image

The dots are differences with other barometers. The red is the time-wise median. The yellow curve is a spline and mainly showing how the mean is behaving, but since the mean is heavily influenced by outliers, is not a good measure here. The median (red) is the one that should be stable.

BAOL815X is a good barometer. The red noise is quite stable. And the following one is clearly drifting:

image

This one seems to have started drifting:

image

And what about this one:

image

Don't understand what might be causing this periodic pattern.

Some problems with this approach:

ToDo