Multivariate geostatistics for the predictive modelling of the surficial sand distribution in shelf seas

Multivariate geostatistics have been used to obtain a detailed and high-quality map of the median grain-size distribution of the sand fraction at the Belgian Continental Shelf. Sandbanks and swales are the dominant geomorphological features and impose a high-spatial seafloor variability. Interpolation over complex seafloors is difficult and as such various models were investigated. In this paper, linear regression and ordinary kriging (OK) were used and compared with kriging with an external drift (KED) that makes use of secondary information to assist in the interpolation. KED proved to be the best technique since a linear correlation was found between the median grain-size and the bathymetry. The resulting map is more realistic and separates clearly the sediment distribution over the sandbanks from the swales. Both techniques were also compared with a simple linear regression of the median grain-size against the bathymetry. An independent validation showed that the linear regression yielded the largest average prediction error (almost twice as large as with KED). most static sedimentological maps, our approach allows for defining grain-size classes that can be adapted according to the needs of various applications. These relate mainly to the mapping of soft substrata habitats and of the most suitable aggregates for extraction. This information is highly valuable in a marine spatial planning context.