A non-stationary geostatistical approach to multigaussian kriging for local reserve estimation