Indicator Kriging in ArcGis Pro (2.6)

Question

Dear all,

I have a set of boring data that I want to use to interpolate to find the hydraulic conductivity of the aquifer. Particularly, I want to look at the semi-impermeable layer that separates the 1st from the 2nd aquifer. I divided my 3D dataset in layers, so that I can work in 2D (which should be OK in a sedimentary aquifer). I have transformed my data into a binary set where 0=sand and 1=clay.

It seems to me this should be a relatively simple 2d indicator kriging exercise. I only find an esri course on EBK online, not on indicator kriging. Thus, I am not sure how to play with the parameters in the geostatistical wizard and I don't know how to interpret for example the semivaiogram and the indicator prediction for a binary dataset.

Does anyone have tips for documentation/tutorials on this? Or does anyone have tips for my work?

Thanks a lot,

Suzanne

SuzanneFaber · Answer

Hi Eric, thanks, I'll have a look at where I get to with EBK and let you know!

EricKrause · Answer

It may not work with your data, but the idea I had was to Empirical Bayesian Kriging with a Probability map output.  If you had a continuous variable with a threshold indicating the cutoff between sand and clay, you could just interpolate the continuous variable with EBK and calculate the probability that each new location exceeds that threshold using the predicted value and the standard error.  The Probability output type does this calculation automatically.

The biggest problem with Indicator Kriging is that when you reclassify from a continuous variable to a binary variable, you lose lots of information.  This is why Indicator Kriging usually doesn't perform well.  But if you only have indicator variables to start out and didn't create them with a reclassification, indicator kriging is about the best you can do.

And just to add to the complexity, there's another type of kriging called Probability Kriging.  It takes a continuous variable and a threshold.  It then reclassifies the values to 0 and 1 and performs indicator kriging.  However, it also uses the raw values of the continuous variable as a cokriging variable.  This essentially tries to bring back in some of the information that was lost during the reclassification.  In my experience, just using EBK with the Probability output type almost always yields best results (again, if you have a continuous variable), but there are a bunch of options available.

SuzanneFaber · Answer

Dear Eric,

thanks a lot for your explanation. What you explain about indicator kriging makes sense and I see that in this case, the indicator kriging did not supply me with a very trustable probability map!

The data that I have, have mainly come from visual analysis of soil samples. So basically, if clay was seen/felt, I recorded clay, if sand was seen I recorded sand. I also have some data from cone penetration tests (measuring resistance, with also gives me a soil classification more than an actual hydraulic conductivity K value). I also used publicly available data on the local aquifers' conductivities. In the set of K values I created from this data, clay has a value of 0.005 m/d and sand ranged from (a little bit of) fine sand (3 m/d), through to sand (15 m/s), all the way to course (50 m/d) and very course (80m/d) sand with pebbles.

I am not sure how I could improve my data in order to get better kriging results.. but curious to hear what you think!

Suzanne