How to generate the probability map of a variable under a threshold using gstat?

Hi All,

I want to produce a probability map of a soil variable under a
threshold. I know I can make it in ArcGIS, but not sure using gstat. Any
tips are highly appreciated.

Regards,

Yong Li