This study presents a geostatistical simulation approach for fine-scale grain size mapping in tidal flats, which complements sparse field survey data with high-resolution optical satellite imagery and quantifies prediction uncertainty at unsampled locations. Within a multi-Gaussian regression kriging (MGRK) framework, a random forest (RF) regression model is used to estimate the trend component of grain size variability in Gaussian space. Residual components are estimated using kriging, and the trend and residual components are combined to construct conditional cumulative distribution functions for uncertainty modeling. Sequential Gaussian simulation based on the CCDFs generates alternative realizations of grain size, allowing for quantification of prediction uncertainty. The potential of this integrated approach was tested on the Baramarae tidal flat in Korea using KOMPSAT-2 imagery. Three spectral features, the green band, red band, and normalized difference water index (NDWI), explained 42.74% of the grain size variability, with NDWI identified as the most influential feature, contributing 40.8% compared with 31.7% for the red band and 27.5% for the green band. MGRK effectively captured local grain size variations, reducing the mean absolute error from 0.554 to 0.280 compared with univariate kriging based solely on field survey data, corresponding to an improvement of approximately 49.5%. The benefit of the proposed approach was validated by a reduction in prediction uncertainty, with the mean standard deviation decreasing from 0.743 in simulations based solely on field data to 0.280 in MGRK-based simulations. These findings indicate that the proposed geostatistical approach, integrating satellite-derived features, is a reliable method for fine-scale mapping of intertidal sediment grain size by providing both predictions and associated uncertainty estimates.



Keywords: tidal flat; geostatistical simulation; random forest; prediction uncertainty; satellite imagery