Sediment samples were collected in the summer of 2017 over the entire study area at 500 m spacing and in the channels. Grain size analysis of the collected sediment samples was conducted using a Malvern Mastersizer 2000© laser beam grain sizer. The median grain size (d50 in µm) at each sample location was then interpolated over the entire study area. Interpolation was made with the SAGA-GIS software (version 7.9.0). According to the cross-validation error, the best method for the D50mm interpolation was the Modified Quadratic Shepard. The 10-fold validation provided an R² of 0.93, an NMRSE of 24.5, an RMSE of 83.9, an MRE of 7039 with the fit set to “node”, the quadratic neighbours and weighting neighbours set to 50 and the spatial resolution was set to 10 m. The resultant interpolation map was then categorized following the nomenclature of Blott and Pye (2001) provided in the file style_sediment_map.qml.