“…To create average yearly images, a geostatistical approach was used, which included (Soares, 2000;Paralta and Ribeiro, 2003;Guimarães et al, 2014;Quental et al, 2012;Roxo et al, 2016) For geostatistical modelling purposes the geoMS software and a parallelised version of the sequential Gaussian simulation were used (Alexandre and Almeida, 1998;Nunes and Almeida, 2010). Due to the broad dispersion of data with no preferential orientation, it was decided to compute averages for the variograms encompassing all directions together (omnidirectional variograms) (Isaaks and Srivastava, 1989;Goovaerts, 1997).…”