Next Article in Journal
Local Defogging Algorithm for the First Frame Image of Unmanned Surface Vehicles Based on a Radar-Photoelectric System
Previous Article in Journal
Target Strength Measurements of Free-Swimming Sandeel Species, Ammodytes spp., in a Large Indoor Experimental Aquarium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Retrieval of Remotely Sensed Sediment Grain Size Evolution Characteristics along the Southwest Coast of Laizhou Bay Based on Support Vector Machine Learning

1
Coastal Research Institute, Ludong University, Yantai 264025, China
2
School of Resources and Environmental Engineering, Ludong University, Yantai 264025, China
3
Key Laboratory of Ecological Restoration and Conservation of Coastal Wetlands in Universities of Shandong (Ludong University), Ludong University, Yantai 264025, China
*
Author to whom correspondence should be addressed.
Xiang Yu and Chao Zhan are co-first authors.
J. Mar. Sci. Eng. 2022, 10(7), 968; https://doi.org/10.3390/jmse10070968
Submission received: 11 May 2022 / Revised: 9 July 2022 / Accepted: 11 July 2022 / Published: 14 July 2022
(This article belongs to the Section Geological Oceanography)

Abstract

:
Grain size is the basic property of intertidal zone sediment. Grain size acts as an indicator of sedimentary processes and geomorphological evolution under human and nature interactions. The remote sensing technique provides an alternative for sediment grain-size parameter monitoring with the advantages of wide coverage and real-time surveying. This paper attempted to map the distributions of three sediment grain size contents and the mean grain size with multitemporal Landsat images along the southwestern coast of Laizhou Bay, China, from 1989 to 2015. Considering the low correlations between the measured reflectance and grain-size parameters, we used a support vector machine (SVM) to develop a nonlinear calibration model by taking several band indices as input variables. Then, the performance of the back propagation neural network (BPNN) was determined and discussed with that of the SVM. The SVM performed better than the BPNN in calibrating the four grain-size parameters based on a comparison of R2 and the root-mean-square error (RMSE). Moreover, an atmospheric correction algorithm originally proposed for case II water enabled the TM\ETM+ images to be precisely atmospherically corrected in this study. The SVM-mapped spatial-temporal grain-size variation showed a coarsening trend, which agreed with that obtained during in situ measurements in a former study. The changes in Yellow River discharge and precipitation associated with the coarsening trend were further analyzed. The yielded results showed that the coarsening trend and reduction in tidal flat area might be aggravated with overutilization. More reasonable planning would be necessary in this case.

1. Introduction

The intertidal zone is an important component of silt coastal landforms. Siltation occurs in the intertidal zone under the interaction of rivers and tides [1], and this process closely relates to hydrodynamics and sedimentology [2]. Changes in sediment grain-size parameters act as indicators of the change rate of sediment pollutants being discharged to the marine environment [3]. Similarly, these changes closely relate to nutrient cycling by affecting primary the production, survival, and growth of benthic animals, flatfish, and waders [4]. There is an urgent need for synoptic knowledge of the sediment grain-size parameter distribution in the intertidal zone [5].
Sediment grain-size parameters vary spatially and temporally in the intertidal zone [6]. In situ measurement is a general method for sediment grain-size parameter surveys. However, sampling sites following a well-defined strategy are inherently discrete and isolated. Whether spatially or temporally, it is difficult to gain an accurate and continuous change trend. In addition, the field survey is time-consuming and laborious. It is essential to develop a remotely sensed approach to periodically provide a quantitative sketch of sediment grain-size parameter distributions.
Intertidal sediment types have been classified by multispectral satellite remote sensing [7], multispectral airborne remote sensing [8], and hyperspectral airborne remote sensing [9,10]. For example, sediment types were classified into different types by supervised classification methods [11]. It was acknowledged that the validity of the training datasets dominated the supervised classification methods. The spectral features of these particular training datasets should typically work independently. Thus, user interaction is unavoidable. Unsupervised classifications with little user interaction were also adopted in producing different clusters by setting a different number of clusters and allowable dispersion around cluster centers [12,13]. The related studies focused on sediment types classification rather than grain-size parameter quantification.
Commercial satellite images with high spatial resolution have been effective for grain-size parameter retrieval, such as IKONOS [14] and WorldView-2 [15]. However, the related commercial satellite data are not free for most scholars. Sensors at low altitudes, such as aircraft and unmanned aerial vehicles (UAVs), have been validated in quantifying sediment/soil grain size by using their hyperspectral reflectance in the VNIR (visible/near-infrared) band or SWIR (shortwave infrared) band [16,17,18,19] or texture properties [20,21]. However, the narrow spatial coverage and inaccessibility of hyperspectral images decreased their utility in surveying the spatial-temporal variation in sediment grain-size parameter distributions. An alternative satellite remote sensing approach based on multispectral images rather than hyperspectral images is needed. Multispectral images were accessible and could provide periodic sketch data. However, the spectral resolution of multispectral images was inherently lower than that of hyperspectral images, resulting in low correlations between band reflectance and grain-size parameters [17]. Thus, it is difficult to develop a linear calibration model with a certain band of multispectral images. Advanced algorithms based on different band combinations were adopted. For example, multispectral satellite images have been adopted to map the sediment distribution in the intertidal areas of Wash, East England (UK) [22]. Multiple regression (MR) and spectral mixture modeling were used to quantify the amounts of each sediment type within a pixel. The yielded result was acceptable when compared with the measured sample sites. Similarly, MR calibration algorithms developed with the reflectance of TM2 and TM5 succeeded in mapping the mud content of the sediment in the Westerschelde (southwest Netherlands) with a maximum R2 of 0.72 [23]. Regression-based algorithms were developed for monitoring the sediment grain-size of intertidal flats, using information from both space-borne microwave (SAR) and optical/shortwave infrared remote sensing [24]. Moreover, principal component analysis (PCA) and wavelet neural networks (WNNs) performed well in sediment grain-size parameter retrieval along the muddy coast of the central Jiangsu Province with HJ-1 satellite images [25]. These studies demonstrated the utility of multispectral satellite images in mapping grain-size parameters.
This paper attempted to map the spatial-temporal variation in sediment grain-size parameter distributions based on Landsat images on the southwestern coast of Laizhou Bay, China, from 1989 to 2015. The main goals are described as follows: to analyze the correlations between sediment grain-size parameters and in situ reflectance measurements; to develop a nonlinear calibration model with Landsat band reflectance by using a support vector machine (SVM), which was subsequently compared with a back propagation neural network (BPNN); to retrieve sediment grain-size parameter distributions based on the validated SVM model; to analyze the spatial-temporal changes in sediment grain-size parameter distributions associated with the changes in Yellow River discharge and precipitation over 26 years. A flow diagram of the methodological steps is sketched in Figure 1.

2. Materials and Methods

2.1. Study Area

The study area, the southwestern coast of Laizhou Bay, is jointly dominated by coastal currents, Yellow River discharge, and precipitation. This area lies in the northeastern part of the Yellow River Delta, a continental weak alluvial tidal delta. Laizhou Bay sea has significant continental climate characteristics, which are greatly affected by cold waves in winter and hot waves in summer [26]. Due to climate change, the fresh water supply in the southern coast of Laizhou Bay fluctuates, and the decrease in rainfall leads to the decrease in runoff into the sea and aggravates the intensity of seawater invasion, and then leads to continuous degradation of the whole coastal wetland ecosystem. The landforms along Laizhou Bay can be divided into three categories: terrestrial landforms, buried landforms, and coastal landforms. Among the landforms, the accumulation plain is the most widely distributed area. Buried landforms are mainly buried ancient channels, which are the main channels of seawater invasion inland. The coastal landform mainly includes bedrock coast, sandy coast, and silty silt coast [27]. The intertidal zone is roughly continuous along Laizhou Bay and unvegetated. The estuary, intertidal zones, and landforms in the Yellow River Delta have frequently changed in recent decades. With the inclusion of the Yellow River Delta Efficient Ecological Economic Zone in China’s National Strategic Planning Area, urban planning, industrial structure and layout, environmental protection, and resource utilization urgently need to understand future intertidal zone sedimentary processes and predict coastal landform erosion and siltation trends.
The Yellow River Delta coast has a typical weak tide. The coastal tide is mainly dominated by the M2 tidal wave. There is no tidal point in the nearshore seawater north of pile No. 5. The average tidal range along the deltaic coast is approximately 1 m (0.6–1.3 m). Most of the shores have irregular semidiurnal tides. Only the shores near Shenxiangou have irregular full-day tides. The estuary during the flooding period is dominated by the reciprocating flow parallel to the river channel. The coasts on both sides of the Yellow River mouth are reciprocating flows of roughly parallel shorelines [28]. Depending on the observations at the northern end of the Yellow River Harbor Terminal (water depth 7 m) in 1987, the wave frequency with a wave height greater than 1.5 m was 11.8%. The normal wave was SSE-S [29]. The strong wave was NE-NNE. The maximum wave was NE with a height of 5.3 m and a cycle of 8 s. The sediment, which is transported by the Yellow River from the Loess Plateau, is mainly composed of sand, silt, and clay. A grain size of less than 0.063 mm accounts for 94.2% of the total.

2.2. Field Measurements and Laboratory Analysis

We conducted 30 field surveys to measure the surface sediment reflectance from 2 March 2013 to 19 July 2014. Sediment samples were synchronously collected from the muddy intertidal zones on the southwestern coast of Laizhou Bay. In situ reflectance measurements and sediment sample collections were conducted during ebb tide or low tide off the southwestern coast of Laizhou Bay with the surface of the intertidal zone exposed. In this paper, a total of 238 sets of intertidal zone surface sediment samples and their reflectances were obtained (Figure 2). We divided the study area into two parts, the west coast and the south coast, to delineate the spatial variation in grain-size parameters. The Xiaoqing River, which is recognized as the boundary of the abandoned southern Yellow River subdelta (ASYS) and the plain of south Laizhou Bay (PSLB), was set as the boundary of the two parts (Figure 2).

2.2.1. Sediment Grain-Size Measurement

Sediment samples distributed in the soil layers at 5–10 cm depth were gathered. Each sample was mixed with five adjacent samples and placed in a numbered polymer bag. We used a Mastersize 2000 laser analyzer to measure grain size distribution, including clay, silt, sand, mean grain size, sorting coefficient, kurtosis, and skewness [30] (Folk and Ward 1957). The principle is described as follows: The intensity of scattered light was inversely proportional to the diameter of the particle under laser beam irradiation. The intensity of scattered light was logarithmically attenuated as the particle size increased. By accepting and measuring the energy distribution of the scattered light, the particle could be obtained. The measuring range of the particle size analyzer was from 0.02 to 2000 μm. In the measurement interval, a total of 100 particle size percentage data points could be generated. The percentage of each particle size was given. The particle size distribution curve was also provided, as were the cumulative curves. Compared with other particle size analyzers, the Mastersize 2000 laser particle size analyzer had a wide range of measurements, high precision, lower sample volume, faster data acquisition, higher measurement accuracy for small particles, and an overall error of less than 1%.
We added 10 mL of 10% H2O2 to the sediment sample (0.2~0.5 g) and heated it on a hot plate to fully react to effectively remove the organic matter. Next, we added 10 mL of 10% HCl and boiled it to react and fully eliminate carbonate. The beaker was filled with distilled water and permitted to stand overnight with the supernatant removed. Ten milliliters of a sodium hexametaphosphate dispersant with a concentration of 0.05 mol/L was added and shaken for 10 min in an ultrasonic cleaner. The dispersant was later measured by using a Mastersize 2000 particle size analyzer. In this experiment, the same treatment method was applied to all samples; that is, each sample was subjected to 5 repeated measurements. Several indices, including the mean grain size (Mz), sorting coefficient (σ1), kurtosis (KG), and skewness (SKI), were obtained with Equations (1)–(4), respectively.
Mz = (ϕ16 + ϕ50 + ϕ84)/3
σ1 = (ϕ84ϕ16)/4 + (ϕ95ϕ5)/6.6
SK1 = (1/2)∗(( ϕ84 + ϕ16 − 2ϕ50)/(ϕ84ϕ16) + (ϕ95 + ϕ5 − 2ϕ50)/(ϕ95 s− ϕ5))
KG = (ϕ95ϕ5)/(2.44∗(ϕ75ϕ25))
where ϕi is a log transformation of grain diameter (D) with base two. The contents of clay, silt, and sand and the mean grain size were subsequently obtained, as shown in Figure 3. From the south coast of Laizhou Bay to the west coast, the clay content of sediment samples varied from 0 to 11.91% and showed an increasing trend (Figure 3a). The slit content of the sediment samples increased from 7.92% to 94.21% and showed a similar trend (Figure 3b). However, the sand content of sediment samples ranged from 90.95% to 1.14% and showed an obvious decreasing trend (Figure 3c). Similar to the clay content and the silt content, the mean grain size, which ranged from 1.69 φ to 6.05 φ, showed an increasing tendency (Figure 3d).

2.2.2. Measurement and Preprocessing of the Sediment Reflectance Spectrum

The in situ reflectance of sediment samples was measured by an ASD Field-Spec FR2500 Spectra-radiometer (Analytical Spectral Devices Inc., Boulder, CO, USA) with a spectral range of 350–2500 nm. A spectrum reflectance panel with 99% (Lab-sphere, Inc., North Sutton, NH, USA) served as the reference standard to adjust and optimize the spectrometer for incoming irradiation. Once the reflectance of one sediment sample was obtained, the spectrometer was optimized. The spectral resolution was 3 nm at 350–1050 nm and 10 nm at 1000–2500 nm. Each sediment sample was measured 10 times with a probe view angle of 25° and a probe distance of 15 cm on sunny and breeze days. Then, we used Savitzky–Golay smoothing to reduce the signal noise to obtain the averaged reflectance curve. The gained reflectance was subsequently matched to Landsat bands by using the spectral response function of the Landsat TM/ETM sensor [31].

2.2.3. Landsat Images and Preprocessing

In view of jet lag, eight cloud-free Landsat images with low water levels were selected (USGS, http://glovis.usgs.gov/, accessed from 13 February 1989 to 5 February 2015.) (Table 1). Digital numbers in visible bands TM1/ETM+1, TM2/ETM+2, TM3/ETM+3, NIR band TM4/ETM+4, and shortwave infrared (SWIR) bands TM5/ETM+5 and TM7/ETM+7 were transformed to apparent reflectances in view of the gains and bias coefficients of scenes, solar azimuth angle, solar spectral irradiance, and Sun–Earth distance [32]. Then, the images were georeferenced to the World Geodetic System 1984 (WGS84) by using the nearest-neighbor resampling method with the selected ground control points. A pixel size of 30 m was retained in the rectified images.
Considering complex aerosol types at land–sea junctions, we adopted an improved atmospheric correction algorithm [33] rather than FLAASH (fast line-of-sight atmospheric analysis of spectral hypercubes). MODIS Terra data were adopted here to provide essential radiometric calibration information in eliminating air molecules and aerosol scattering from the apparent reflectance of Landsat images. Although the algorithm was originally developed for atmospheric correction of case II waters [33], it precisely described the information of the specific aerosol types in tidal flats due to the sensitivity and designed spectral bands of MODIS Terra data. Thus, it was more suitable for atmospheric correction in tidal flats than FLAASH, which had four restrictive aerosol models (rural, maritime, urban, and tropospheric). The flow is shown as follows.
The scattering radiance (Lr) due to air molecules (Rayleigh scattering) was calculated using the following methodology [34].
Lr(λ) = ρr(λ) F0(λ) cos(θ0)/π
where ρr(λ), F0(λ), and θ0 denote the Rayleigh scattering reflectance, extraterrestrial solar irradiance, and solar zenith angle, respectively. F0(λ) could be obtained from the website (https://oceancolor.gsfc.nasa.gov/, accessed on 1 January 2022). θ0 was presented in the header file of Landsat images. Then, MODIS Terra data were adopted to gain ρr(λ) of the Landsat images. The Rayleigh scattering reflectance ρr(λ)∝τr(λ)∝λ−4 [35]. The conversion factor β, from the MODIS/Landsat band response-averaged Rayleigh optical thickness (τr(λ)) to that of the center wavelength, was calculated as Equation (6):
β = τr(λ)/<τr(λ)>MODIS =ρr(λ)/<ρr(λ)>MODIS
where
τr(λ) = 0.008569λ−4 (1+0.0113λ−2+0.00013λ−4)
< τ r ( λ ) > MODIS = τ r ( λ ) S ( λ ) F 0 ( λ ) d λ / S ( λ ) F 0 ( λ )
Then, β for MODIS bands were obtained. <ρr(λ)>MODIS was accessible in the MODIS Rayleigh lookup tables (LUTs) (http://seadas.gsfc.nasa.gov/, accessed on 1 January 2022). ρr(λ) could be calculated by <ρr(λ)>MODIS and β. Lr at the bands of Landsat images was derived by means of band response-averaged calculation [36] with Equation (5).
MODIS images acquired approximately synchronously with Landsat images were processed using the atmospheric correction algorithm [37] with the aerosol type provided in the aerosol LUTs. Aerosol scattering radiance (La5)) could be obtained after subtracting Lr5) from Lt5) under the assumption of black water at TM/ETM+ band 5 (λ5). Then, ρa5) was calculated by Equation (5). The relationship between the single scattering reflectance (ρas(λ)) and ρa(λ), which were accessible in the aerosol LUTs, was applied to convert ρa(λ5) to ρas(λ5). Once ρas(λ5) was obtained, ρasi) (i = 1, 2, 3, 4) was interpolated with the spectral relationship in aerosol LUTs. Similarly, ρasi) was inversely converted to ρai) using the relationship between (ρas(λ)) and ρa(λ). Then, Lai) was derived from Equation (5). Surface reflectance (Rrs) could be derived from Equation (9):
Rrs(λ) = Ls(λ)/ts(λ) ρr(λ) F0(λ) cos(θ0)
where surface radiance (Ls(λ)) and diffuse transmittance ts(λ) are given in Equations (10) and (11), respectively.
Ls(λ) = Lt5) − Lr5) − La5)
ts(λ) = exp (−(τr(λ))/2 + τoz(λ))/cos(θ0))
where τoz(λ) is a multiplication between the ozone absorption coefficient at a specific wavelength (koz(λ)) [38] and the ozone content (Uoz). Uoz is accessed from 13 February 1989 to 5 February 2015 (http://oceandata.sci.gsfc.nasa.gov/Ancillary/Meterological/ (2012), accessed on 1 January 2022). The algorithm performed well by the comparison between mean spectra of 238 sites after atmospheric correction and that of in situ measurement (Figure 4).

2.3. Sediment Grain-Size Retrieval

2.3.1. Spectral Indices

The reflectance of the Landsat image bands (b1, b2, b3, b4, b5, and b7) was related to sediment grain-size parameters. The reflectance of the Landsat image bands (b1, b3, b4, and b5) was significantly correlated with sediment grain-size parameters (r > 0.141, p < 0.05) (Figure 5a). However, the correlation coefficients were not high enough to develop stable regression models, even if the highest correlation coefficient between b4 (center wavelength at 0.83 microns) and the sand content was adopted (Figure 5b).
Spectral indices, such as SI, DI, PI, RI, and NDI, were subsequently calculated based on the Landsat image band reflectance. The indices were also correlated with sediment grain-size parameters (Table 2). By contrast, DI performed best in relating to the four grain-size parameters, with rmax values of 0.65, 0.71, 0.70, and 0.64. However, the correlation coefficients were still not high enough to promote robust linear regression models with corresponding R2 values of 0.42, 0.50, 0.49, and 0.41. Thus, we used the SVM [39] to calibrate the grain-size parameters. Seventy-five percent of all samples were randomly selected as the calibration dataset, and the remaining 25% were used as the validation dataset.

2.3.2. Estimation Models

SVM was conducted by using R, SI, DI, PI, and RI as input variables. The basic idea of SVM is described as follows: to realize some nonlinear transformation by defining an appropriate kernel function, to map the input vector to a high-dimensional space, to find the optimal classification hyperplane in high-dimensional space, and to achieve the algorithm by the inner product in high-dimensional space [40]. Once an appropriate kernel function was selected, we obtained the classification function of high-dimensional space [41].
We used the scaleForSVM normalized preprocessing function to normalize the training dataset and validation dataset. All datasets were transformed into the range between ‘0’ and ‘1’. Because the radial basis kernel function (RBF) could solve the nonlinear problem of the sample into a higher dimensional space, it could solve the nonlinear problem that the linear kernel function could not solve [42]. Thus, RBF was adopted here. The core of the SVM model was to determine the best c and g combination so that the resulting root-mean-square error was reduced. The simplest solution to this issue was to obtain c and g discrete values in a certain range. We used the grid cross-validation method to locate the best combination of c and g here (Figure 6). The SVM optimal calibration and validation models for clay, silt, sand, and mean grain size are presented in Figure 7.

3. Results

SVM was subsequently used in mapping grain-size parameters with Landsat TM/ETM+ images from 1989 to 2015. The spatial and temporal distributions of the clay, silt, sand, and mean grain size contents are shown in Figure 8, Figure 9, Figure 10 and Figure 11, respectively.

3.1. Spatial Variation in Grain-Size Parameters

In 1989, the clay content increased perpendicular to the coastline on the west coast. In contrast, the clay content showed a declining trend on the south coast. The clay content was roughly the same coastwise in 1995, 1999, 2001, and 2006. In 2009, the clay content remained roughly the same coastwise on the west coast. However, the clay content showed a declining trend perpendicular to the coastline on the south coast. The spatial variation in clay content in 2012 was similar to that in 2009, as was that in 2015. The smaller difference was that the clay content in 2015 was slightly lower than that in 2012. In 1989, the silt content showed a similar trend as the clay content on the west coast and the south coast. The silt content remained roughly the same coastwise from 1995 to 2012. In 2015, the silt content was uniform coastwise, but there was a significant decrease from that in 2012. In 1989, the sand content also showed a similar trend as the clay content on the west coast and the south coast. From 1995 to 2009, the sand content remained roughly uniform coastwise. However, these values were higher in 1995 and 2009 than in 1999, 2001, and 2006. The spatial variation in sand content increased in 2012 and especially in 2015. The spatial variation in the mean grain size in 1989, 1995, 1999, 2001, and 2012 was higher than those in 2006, 2009, and 2015. Moreover, there was an obvious increase in 2015 from 2012.

3.2. Temporal Variation in Grain Size

Mean values of the retrieved grain-size parameters were calculated and applied to describe the temporal variation in grain size during 1989 and 2015 (Table 3). The clay content decreased slightly from 5.32% in 1989 to 4.64% in 2015 and showed a decreasing trend with an average annual reduction rate of 0.005. Similarly, the silt content decreased from 59.62% in 1989 to 37.42% in 2015 and showed a clear decreasing trend with an average annual reduction rate of 0.018. The sand content increased from 48.78% in 1989 to 58.46% in 2015 and showed an increasing trend with an average annual growth rate of 0.007. Therefore, the mean grain size increased from 4.05 in 1989 to 4.10 in 2015 and showed an increasing trend with an average annual growth rate of 0.005. In general, a coarsening trend occurred in the southwestern part of Laizhou Bay over 26 years. A previous study indicated that a similar transition from viscous sand to nonviscous sand occurred here, with a decrease in the clay content and viscosity/mud ratio from 2007 to 2013 [43].
The hydrodynamics of tides on the southwestern coast of Laizhou Bay have been stable since the Yellow River estuary changed north of Qingshuigou in 1996 [44]. Although the decreasing trend of the Yellow River annual discharge accessed form 1989 to 2015 (http://www.yrcc.gov.cn/, accessed on 1 January 2022) was opposite to that of mean grain size, the negative correlation between them was not significant, indicating that grain size coarsening had no correlation with the current Yellow River discharge (Table 3, Figure 12). Similarly, although the increasing trend of precipitation was similar to that of mean grain size, the positive correlation between precipitation and mean grain size was not significant, indicating that grain size coarsening had no correlation with the current precipitation, either (Figure 12).

4. Discussion

To minimize the influence of atmospheric path radiation, we adopted an improved atmospheric correction algorithm [33] rather than FLAASH in the case of the existence of complex aerosol types at land–sea junctions. MODIS Terra data described the specific aerosol types in tidal flats due to the sensitivity and designed spectral bands. The adopted atmospheric correction algorithm for TM/ETM+ images performed well with the aid of MODIS Terra data in this study (Figure 4).
Considering the low correlations between the measured reflectance and grain-size parameters, another advanced calibration algorithm, the BPNN, was tested here for comparison with SVM. The BPNN can learn and store the mapping relationship of a large number of input–output modes [45]. The BPNN could eliminate the mathematical formula for describing and revealing the required mapping relationship beforehand. The learning method of BPNN adopted the steepest descent method to continuously adjust the weight and threshold of the network through reverse transmission so that the square sum of errors in the network was minimized [46]. Additionally, we adopted R, SI, DI, PI, and RI as input variables. The nodal number of hidden layer neurons was set between 2 and 20 based on Zhang’s method [47]. The output layer transfer function of the single hidden layer model was placed at the Purelin function. The hidden layer transfer function was set as a tansig tangent s-type function [48]. The training functions were trained functions based on the gradient descent algorithm [49]. The maximum number of iterations, learning rate, and training accuracy were set to 1000, 0.05, and 1 × 10−5, respectively. The remaining parameters were set as default values. The optimal BPNN calibration models for grain-size parameter retrieval were obtained (Figure 13, Table 4). The SVM performed better than the BPNN based on a comparison of R2 and RMSE. In addition, there were some inherent issues in the BPNN, such as how to determine local minimum points, network structure problems, overlearning, and underlearning [50,51]. The SVM transformed the problem to be solved as a problem of quadratic programming optimization. In theory, the global optimal solution was guaranteed, and the local solution and overlearning in the neural network training process were solved [52,53]. Thus, the SVM was finally adopted here in the grain-size parameter retrieval with multitemporal TM/ETM+ images over 26 years. Compared with multispectral satellite images, airborne images, such as simulated Daedalus 1268 Airborne Thematic Mapper (ATM) data, Compact Airborne Spectrographic Imager (CASI) data, and high-spatial-resolution satellite data (IKONOS), have been used in mapping grain-size parameters. Recombining the subpixel end member abundance through multivariate regression analysis significantly improved the image calibration for both sediment clay and sand content (R2 > 0.8) in the ATM images [8]. Two robust adjustment techniques (MVE and MCD multivariate M-estimators of location and scale) provided acceptable algorithms for grain-size mapping with CASI-2 data in Santander Bay, Spain [9]. Furthermore, calibration algorithms provided an estimated R2 of 0.91 for mud flat facies and 1 for sand flat facies with IKONOS data for the Hwangdo tidal flat, Cheonsu Bay, Korea. The calibration algorithms provided an estimated R2 of 0.93 for sandy loam, 0.94 for silty loam, and 0.67 for clay loam [54]. The calibration accuracies in this paper were slightly below those relying on airborne images in the aforementioned literature. However, the estimated accuracies were acceptable, with R2 values of 0.66 for clay, 0.81 for silt, 0.83 for sand, and 0.80 for mean grain size. In addition, the accessibility and cost-effectiveness of multispectral satellite images enabled the grain-size parameters to be periodically mapped with wide coverage. Notably, the effect of water on the reflectance was not considered in this paper. The influence was difficult to differentiate and eliminate from the reflectance. The reflectance of dried soil samples in the laboratory might develop an actual correlation between reflectance and grain-size parameters. However, the laboratory-acquired correlation was difficult to use to map grain-size parameters in the field. From another point of view, the water content was influenced by grain-size parameters. To some extent, there was an inverse correlation between grain size and water holding capacity [55,56]. The smaller the grain size of the soil is, the greater its water holding capacity. The effect of water content on the reflectance was inherent; however, it was affected by grain-size parameters, especially under moisture equilibrium conditions for the same sediment. A recent study had illustrated that the moisture of the same sand (a median diameter of about 250 μm) in a 1 m resolution laser image could differ from almost 0% (dry) to over 20% (fully saturated), with a strong effect on surface reflectance of a near-infrared band (1550 nm) [57]. However, the effect might be reduced here as the measured mean grain-size of water-holding sediment was less than 6.05 φ (from 2 to 66.26 μm) and the large pixel size (30 m) of the satellite imagery used here. We preferred to address this issue in further studies. Additionally, given the increasing availability of high-temporal- (almost daily) and spatial-resolution (3 m) and high-radiometric-resolution (12-bit) images from CubeSats such as PlanetScope imagery [58,59,60], the application of this imagery in mapping grain size based on spectral and textural features would be an alternative for further investigation.
The yielded results of the study showed that the area of tidal flats has been greatly reduced under the interaction of humans and nature in recent decades (Figure 8, Figure 9, Figure 10 and Figure 11). In addition, the remotely sensed coarsening trend of grain size indicated that erosion had overtaken siltation. The coarsening trend might be aggravated by the overutilization. More reasonable planning for the utilization of tidal flats is needed.

5. Conclusions

An experiment aimed at quantitatively mapping the sediment grain-size parameter distribution in the intertidal zone of the southwestern coast of Laizhou Bay is presented in this paper. A total of 238 sets of intertidal zone surface sediment samples and their reflectance spectra were measured during nearly 30 field surveys from 2 March 2013 to 19 July 2014. Then, in situ reflectance measurements were matched with the Landsat band reflectances by using the corresponding spectral response function. The generated Landsat band reflectances (b1, b3, b4, and b5) were significantly correlated with the clay, silt, sand, and mean grain size contents. However, the correlations between the original band reflectances and the sediment grain-size parameters were not high enough to develop a stable regression model, as were the correlations between several band indices (Table 2) and the sediment grain-size parameters. Then, two nonlinear fitting methods (SVM and BPNN) were compared in the sediment grain-size parameter calibration due to their efficient nonlinear fitting capability when the linear fitting did not work. However, there were still several issues to be solved in the BPNN, such as poor generalization ability and overfitting. Thanks to the SVM, the training error was minimized, and the generalization ability was maximized by setting an RBF. Moreover, in contrast with R2 and RMSE, the SVM performed better than the BPNN did.
The spatial-temporal distributions of the clay, silt, and sand contents and mean grain size over 26 years from 1989 to 2015 were finally obtained with multitemporal Landsat TM/ETM+ images based on the validated SVM model. The yielded results showed that there was an obvious decreasing trend in the contents of clay and silt on the southwestern coast of Laizhou Bay. However, there was an increasing trend in sand content and mean grain size here. The coarsening trend agrees with the conclusion obtained through in situ measurements in a former study. Although the reduction in Yellow River discharge and the increase in precipitation were not related to the coarsening of grain size significantly, the coarsening trend and the reduction in tidal flat area might be aggravated by the overutilization. More reasonable planning is required.

Author Contributions

Conceptualization, X.Y. and Q.W.; methodology, C.Z.; software, Y.L.; software, J.B.; validation, G.L., L.W. and X.L.; formal analysis, Q.W.; investigation, X.Y.; resources, C.Z.; data curation, G.L.; writing—original draft preparation, X.Y.; writing—review and editing, Y.L.; visualization, B.C.; supervision, Q.W.; project administration, Q.W.; funding acquisition, Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work is jointly supported by the National Natural Science Foundation of China (C.Z., 41901006), the Shandong University Qingchuang Science and Technology Team (C.Z., 2020KJH002), the Shandong Provincial Natural Science Foundation (X.Y., ZR2020MD082; G.L., ZR2021MD051; X.L., ZR2019PD013), the National Science Foundation of China-Shandong Unified Fund (Q.W., U1706220), Youth Innovation Team Project for Talent Introduction and Cultivation in Universities of Shandong Province (L.W., no granted number) and the National Natural Science Foundation of China (F.W., 41871163, X.L., 41901102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, Y.; Liu, Q.; Huang, R.; Pan, H.; Xu, M. Recent Evolution of Coastal Tidal Flats and the Impacts of Intensified Human Activities in the Modern Radial Sand Ridges, East China. Int. J. Environ. Res. Public Health 2020, 17, 3191. [Google Scholar] [CrossRef] [PubMed]
  2. Chiarella, D.; Longhitano, S.G.; Sabato, L.; Tropeano, M. Sedimentology and hydrodynamics of mixed (siliciclastic-bioclastic) shallow-marine deposits of Acerenza (Pliocene, Southern Apennines, Italy). Ital. J. Geosci. 2012, 131, 136–151. [Google Scholar]
  3. Woodruff, J.D.; Martini, A.P.; Elzidani, E.Z.H.; Naughton, T.J.; Kekacs, D.J.; Macdonald, D.G. Off-river waterbodies on tidal rivers: Human impact on rates of infilling and the accumulation of pollutants. Geomorphology 2013, 184, 38–50. [Google Scholar] [CrossRef]
  4. Tweedley, J.; Chambers, J.; Paice, R. Sediment Accumulation and Resuspension in the Vasse-Wonnerup Wetlands and Its Relationship to Internal Nutrient Cycling; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  5. Anthony, E.J.; Gardel, A.; Gratiot, N.; Proisy, C.; Allison, M.A.; Dolique, F.; Fromard, F. The Amazon-influenced muddy coast of South America: A review of mud-bank-shoreline interactions. Earth Sci. Rev. 2011, 103, 99–121. [Google Scholar] [CrossRef]
  6. Prodger, S.; Russell, P.; Davidson, M. Grain-size distributions on high energy sandy beaches and their relation to wave dissipation. Sedimentology 2017, 64, 1289–1302. [Google Scholar] [CrossRef] [Green Version]
  7. Ngoc, N.T.; Tue, N.T. Correlating mass physical properties with ALOS reflectance spectra for;intertidal sediments from the Ba Lat Estuary (northern Vietnam): An exploratory laboratory study. Geo-Mar. Lett. 2013, 33, 273–284. [Google Scholar] [CrossRef]
  8. Rainey, M.P.; Tyler, A.N.; Gilvear, D.J.; Bryant, R.G.; Mcdonald, P. Mapping intertidal estuarine sediment grain size distributions through airborne remote sensing. Remote Sens. Environ. 2003, 86, 480–490. [Google Scholar] [CrossRef]
  9. Castillo, E.; Pereda, R.; Luis, J.M.D.; Medina, R.; Viguri, J. Sediment grain size estimation using airborne remote sensing, field sampling, and robust statistic. Environ. Monit. Assessment 2011, 181, 431–444. [Google Scholar] [CrossRef]
  10. Deronde, B.; Kempeneers, P.; Forster, R.M. Imaging spectroscopy as a tool to study sediment characteristics on a tidal sandbank in the Westerschelde. Estuar. Coast. Shelf Sci. 2006, 69, 580–590. [Google Scholar] [CrossRef]
  11. Craven, K. Supervised classification of continental shelf sediment off western Donegal, Ireland. In Agu Fall Meeting; American Geophysical Union: Washington, DC, USA, 2017; p. OS31C-1414. [Google Scholar]
  12. Ibrahim, E.; Adam, S.; Monbaliu, J.; Maktav, D. Assessment of unsupervised classification techniques for intertidal sediments. Environ. Manag. 2009, 45, 526–540. [Google Scholar]
  13. Ibrahim, E.; Monbaliu, J. Suitability of spaceborne multispectral data for inter-tidal sediment characterization: A case study. Estuar. Coast. Shelf Sci. 2011, 92, 437–445. [Google Scholar] [CrossRef]
  14. Park, N.W.; Jang, D.H.; Chi, K.H. Integration of IKONOS imagery for geostatistical mapping of sediment grain size at Baramarae beach, Korea. Int. J. Remote Sens. 2009, 30, 5703–5724. [Google Scholar] [CrossRef]
  15. Niroumand-Jadidi, M.; Vitti, A. Grain size mapping in shallow rivers using spectral information: A lab spectroradiometry perspective. In Remote Sensing of the Ocean, Sea Ice, Coastal Waters, and Large Water Regions; SPIE: Bellingham, WA, USA, 2017. [Google Scholar]
  16. Jayson-Quashigah, P.N.; Addo, K.A.; Amisigo, B.; Wiafe, G. Assessment of short-term beach sediment change in the Volta Delta coast in Ghana using data from Unmanned Aerial Vehicles (Drone). Ocean. Coast. Manag. 2019, 182, 104952–104965. [Google Scholar] [CrossRef]
  17. Rossel, R.A.V.; Walvoort, D.J.J.; Mcbratney, A.B.; Janik, L.J.; Skjemstad, J.O. Visible, near infrared, mid infrared or combined diffuse reflectance spectroscopy for simultaneous assessment of various soil properties. Geoderma 2006, 131, 59–75. [Google Scholar] [CrossRef]
  18. Sullivan, D.G.; Shaw, J.N.; Mask, P.L.; Rickman, D.; Luvall, J.; Wersinger, J.M. Remote sensing of near-surface soil properties with the airborne terrestrial applications sensor. Proc. SPIE Int. Soc. Opt. Eng. 2004, 5544, 266–275. [Google Scholar]
  19. Sullivan, D.G.; Shaw, J.N.; Rickman, D.; Mask, P.L.; Luvall, J. Using Remote Sensing Data to Evaluate Surface Soil Properties in Alabama Ultisols. Soil Sci. 2005, 170, 954–968. [Google Scholar] [CrossRef]
  20. Bae, S.; Yu, J.; Wang, L.; Jeong, Y.; Kim, J.; Yang, D. Experimental analysis of sand grain size mapping using UAV remote sensing. Remote Sens. Lett. 2019, 10, 893–902. [Google Scholar] [CrossRef]
  21. Carbonneau, P.E.; Lane, S.N.; Bergeron, N.E. Catchment-scale mapping of surface grain size in gravel bed rivers using airborne digital imagery. Water Resour. Res. 2004, 40, WO7202. [Google Scholar] [CrossRef] [Green Version]
  22. Yates, M.G.; Jones, A.R.; Mcgrorty, S.; Goss-Custard, J.D. The Use of Satellite Imagery to Determine the Distribution of Intertidal Surface Sediments of the Wash, England. Estuar. Coast. Shelf Sci. 1993, 36, 333–344. [Google Scholar] [CrossRef]
  23. Van der Wal, D.; Herman, P.M.J. Regression-bases synergy of optical, shorwave infrared and microwave remote sensing for monitoring the grain-size of intertidal sediments. Remote Sens. Environ. 2007, 111, 89–106. [Google Scholar] [CrossRef]
  24. Kerle, N.; Skidmore, A.K.; Van der Wal, D.; Herman, P.M.J. Quantifying the particle size of intertidal sediments with satellite remote sensing in the visible light, thermal infrared and microwave spectral domain. Bull. Geol. Soc. 2015, 5, 161–172. [Google Scholar]
  25. Liu, X.X.; Zhang, D.; Han, F. Remote Sensing Study on Sediment Grain Size Distribution and its Migration Trend Analysis in Tidal Flat Based on PCA and WNN Model. Oceanol. Et Limnol. Sin. 2019, 6, 1168–1180. [Google Scholar]
  26. Li, S.N.; Wang, G.X.; Deng, W.; Hu, Y.M.; Hu, W.W. Influence of hydrology process on wetland landscape pattern: A case study in the yellow river delta. Ecol. Eng. 2009, 35, 1719–1726. [Google Scholar] [CrossRef]
  27. Cuo, L.; Zhang, Y.; Gao, Y.; Hao, Z.; Cairang, L. The impacts of climate change and land cover/use transition on the hydrology in the upper yellow river basin, China. J. Hydrol. 2013, 502, 37–52. [Google Scholar] [CrossRef]
  28. Wang, H.; Yang, Z.; Li, Y.; Guo, Z.; Sun, X.; Yan, W. Dispersal pattern of suspended sediment in the shear frontal zone off the Huanghe (Yellow River) mouth. Cont. Shelf Res. 2007, 27, 854–871. [Google Scholar] [CrossRef]
  29. Bi, N.; Yang, Z.; Wang, H.; Hu, B.; Ji, Y. Sediment dispersion pattern off the present huanghe (yellow river) subdelta and its dynamic mechanism during normal river discharge period. Estuar. Coast. Shelf Sci. 2010, 86, 352–362. [Google Scholar] [CrossRef]
  30. Folk, R.L.; Ward, W.C. Brazos River bar: A study in the significance of grain size parameters. J. Sediment. Res. 1957, 27, 3–26. [Google Scholar] [CrossRef]
  31. Huang, W.; Huang, J.; Wang, X.; Wang, F.; Shi, J. Comparability of Red/Near-Infrared Reflectance and NDVI Based on the Spectral Response Function between MODIS and 30 Other Satellite Sensors Using Rice Canopy Spectra. Sensors 2013, 13, 16023–16050. [Google Scholar] [CrossRef]
  32. Chander, G.; Markham, B. Revised Landsat-5 TM radiometric calibration procedures and postcalibration dynamic ranges. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2674–2677. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, M.; Dong, Q.; Cui, T.; Xue, C.; Zhang, S. Suspended sediment monitoring and assessment for Yellow River estuary from Landsat TM and ETM+ imagery. Remote Sens. Environ. 2014, 146, 136–147. [Google Scholar] [CrossRef]
  34. Hu, C.; Muller-Karger, F.E.; Andrefouet, S.; Carder, K.L. Atmospheric correction and cross-calibration of LANDSAT-7/ETM+ imagery over aquatic environments: A multiplatform approach using SeaWiFS/MODIS. Remote Sens. Environ. 2001, 78, 99–107. [Google Scholar] [CrossRef]
  35. Gordon, H.R.; Castaño, D.J. Aerosol analysis with the Coastal Zone Color Scanner: A simple method for including multiple scattering effects. Appl. Opt. 1989, 28, 1320–1326. [Google Scholar] [CrossRef] [PubMed]
  36. Gordon, H.R. Remote sensing of ocean color: A methodology for dealing with broad spectral bands and significant out-of-band response. Appl. Opt. 1995, 34, 8363–8374. [Google Scholar] [CrossRef]
  37. Zhang, M.W.; Tang, J.W.; Dong, Q.; Song, Q.T.; Ding, J. Retrieval of total suspended matter concentration in the Yellow and East China Seas from MODIS imagery. Remote Sens. Environ. 2010, 114, 392–403. [Google Scholar] [CrossRef]
  38. Gregg, W.W.; Carder, K.L. A Simple Spectral Solar Irradiance Model for Cloudless Maritime Atmospheres. Limnol. Oceanogr. 1990, 35, 1657–1675. [Google Scholar] [CrossRef]
  39. Tsochantaridis, I.; Hofmann, T.; Joachims, T.; Altun, Y. Support vector machine learning for interdependent and structured output spaces. In Proceedings of the International Conference on Machine Learning, Banff, AB, Canada, 4–8 July 2004. [Google Scholar]
  40. Suykens, J.A.K.; Brabanter, J.D.; Lukas, L.; Vandewalle, J. Weighted least squares support vector machines: Robustness and sparse approximation. Neurocomputing 2002, 48, 85–105. [Google Scholar] [CrossRef]
  41. Liu, H.J.; Wang, Y.N.; Lu, X.F. A method to choose kernel function and its parameters for support vector machines. In Proceedings of the International Conference on Machine Learning & Cybernetics, Guangzhou, China, 18–21 August 2005. [Google Scholar]
  42. Korenjak, M.; Taylor-Harding, B.; Binné, U.L.; Satterlee, J.S.; Stevaux, O.; Aasland, R.; White-Cooper, H.; Dyson, N.; Brehm, A. Native E2F/RBF complexes contain Myb-interacting proteins and repress transcription of developmentally controlled E2F target genes. Cell 2005, 119, 181–193. [Google Scholar] [CrossRef] [Green Version]
  43. Qing, W.; Xiaolu, W.; Xueyan, L.; Xin, W.; Chao, Z. Grain size characteristics and corsening phenomenon of inter-tidal flat surficial sediment along the abandonded southern yellow river sub-delta. Quat. Sci. 2009, 37, 353–367. [Google Scholar]
  44. Li, M.M.; Wang, Q.; Zhang, A.D.; Wang, H.Y.; Liu, Y.L.; Wang, L. Study on the evolution of the tidal morphodynamic processes in south-western Laizhou Bay based on numerical simulation in the past 50 years. Acta Oceanol. Sin. 2014, 5, 68–76. [Google Scholar] [CrossRef]
  45. Pradhan, B.; Lee, S. Regional landslide susceptibility analysis using back-propagation neural network model at Cameron Highland, Malaysia. Landslides 2010, 7, 13–30. [Google Scholar] [CrossRef]
  46. Cheng, H.L.; Park, S.C. Combination of modified BPNN algorithms and an efficient feature selection method for text categorization. Inf. Processing Manag. 2009, 45, 329–340. [Google Scholar]
  47. Zhang, Y.; Ruan, G. Bernoulli Neural Network with Weights Directly Determined and with the Number of Hidden–Layer Neurons Automatically Determined. In International Symposium on Neural Networks on Advances in Neural Networks; Springer: Berlin/Heidelberg, Germany, 2009; pp. 36–45. [Google Scholar]
  48. Garca, P.; Vidal, E. Inference of k-testable languages in the strict sense and application to syntacti.c pattern recognition. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 12, 920–925. [Google Scholar] [CrossRef]
  49. Madgwick, S.O.H.; Harrison, A.J.L.; Vaidyanathan, R. Estimation of IMU and MARG orientation using a gradient descent algorithm. In Proceedings of the IEEE International Conference on Rehabilitation Robotics, Zurich, Switzerland, 29 June–1 July 2011. [Google Scholar]
  50. Rajath, K.M.P.; Keerthi, S.R.; Aishwarya, K.M. Artificial neural networks for face recognition using PCA and BPNN. In Proceedings of the Tencon IEEE Region 10 Conference, Marina Bay Sands, Singapore, 22–26 November 2016. [Google Scholar]
  51. Wang, L.; Zeng, Y.; Chen, T. Back propagation neural network with adaptive differential evolution algorithm for time series forecasting. Expert Syst. Appl. 2015, 42, 855–863. [Google Scholar] [CrossRef]
  52. Suykens, J.A.K.; Vandewalle, J. Least Squares Support Vector Machine Classifiers. Neural Processing Lett. 1999, 9, 293–300. [Google Scholar] [CrossRef]
  53. Tong, S.; Chang, E.Y. Support vector machine active learning for image retrieval. In Proceedings of the Ninth ACM International Conference on Multimedia, Ottawa, ON, Canada, September 30–October 5 2001; pp. 107–118. [Google Scholar]
  54. Choi, J.K.; Ryu, J.H.; Lee, Y.K.; Yoo, H.R.; Han, J.W.; Chang, H.K. Quantitative estimation of intertidal sediment characteristics using remote sensing and GIS. Estuar. Coast. Shelf Sci. 2010, 88, 125–134. [Google Scholar] [CrossRef]
  55. Goldberg, D.; Gornat, B.; Rimon, D. Drip Irrigation: Principles, Design and Agricultural Practices; Elsevier: Amsterdam, The Netherlands, 1976. [Google Scholar]
  56. Mohammadi, M.H.; Vanclooster, M. Predicting the Soil Moisture Characteristic Curve from Particle Size Distribution with a Simple Conceptual Model. Vadose Zone J. 2011, 10, 594–602. [Google Scholar] [CrossRef]
  57. Smit, Y.; Ruessink, G.; Brakenhoff, L.B.; Donker, J.J.A. Measuring spatial and temporal variation in surface moisture on a coastal beach with a near-infrared terrestrial laser scanner. Aeolian Res. 2018, 31, 19–27. [Google Scholar] [CrossRef]
  58. Niroumand-Jadidi, M.; Bovolo, F.; Bruzzone, L.; Gege, P. Physics-based Bathymetry and Water Quality Retrieval Using PlanetScope Imagery: Impacts of 2020 COVID-19 Lockdown and 2019 Extreme Flood in the Venice Lagoon. Remote Sens. 2020, 12, 2381. [Google Scholar] [CrossRef]
  59. Sarah, C.; Laurence, S.; Leon, S.; Joseph, M. Tracking Dynamic Northern Surface Water Changes with High-Frequency Planet CubeSat Imagery. Remote Sens. 2017, 9, 1306–1326. [Google Scholar]
  60. Wicaksono, P.; Lazuardi, W. Assessment of PlanetScope images for benthic habitat and seagrass species mapping in a complex optically shallow water environment. Int. J. Remote Sens. 2018, 39, 5739–5765. [Google Scholar] [CrossRef]
Figure 1. A flow diagram of the methodological steps.
Figure 1. A flow diagram of the methodological steps.
Jmse 10 00968 g001
Figure 2. Sample locations along the southwestern coast of Laizhou Bay, China.
Figure 2. Sample locations along the southwestern coast of Laizhou Bay, China.
Jmse 10 00968 g002
Figure 3. The measured content of clay (a), silt (b), sand (c), and the mean grain-size (d) of the collected samples.
Figure 3. The measured content of clay (a), silt (b), sand (c), and the mean grain-size (d) of the collected samples.
Jmse 10 00968 g003
Figure 4. Mean spectra of 238 sites after atmospheric correction compared with that of in situ measurement.
Figure 4. Mean spectra of 238 sites after atmospheric correction compared with that of in situ measurement.
Jmse 10 00968 g004
Figure 5. Correlations between the band reflectances of Landsat imagery (a) and the measured grain-size parameters. The linear calibration model of sand content based on the reflectance of band 4 (b).
Figure 5. Correlations between the band reflectances of Landsat imagery (a) and the measured grain-size parameters. The linear calibration model of sand content based on the reflectance of band 4 (b).
Jmse 10 00968 g005
Figure 6. Optimal parameters of SVM model for clay (a), silt (b), sand (c), and mean grain-size (d).
Figure 6. Optimal parameters of SVM model for clay (a), silt (b), sand (c), and mean grain-size (d).
Jmse 10 00968 g006
Figure 7. Optimal parameters of SVM model for clay (a), silt (b), sand (c), and mean grain-size (d). SVM optimal calibration and validation model for the content of clay (a,b), silt (c,d), sand (e,f), and mean grain-size (g,h).
Figure 7. Optimal parameters of SVM model for clay (a), silt (b), sand (c), and mean grain-size (d). SVM optimal calibration and validation model for the content of clay (a,b), silt (c,d), sand (e,f), and mean grain-size (g,h).
Jmse 10 00968 g007aJmse 10 00968 g007b
Figure 8. Spatial and temporal distribution of clay content from 1989 to 2015.
Figure 8. Spatial and temporal distribution of clay content from 1989 to 2015.
Jmse 10 00968 g008aJmse 10 00968 g008b
Figure 9. Spatial and temporal distribution of silt content from 1989 to 2015.
Figure 9. Spatial and temporal distribution of silt content from 1989 to 2015.
Jmse 10 00968 g009aJmse 10 00968 g009b
Figure 10. Spatial and temporal distribution of sand content from 1989 to 2015.
Figure 10. Spatial and temporal distribution of sand content from 1989 to 2015.
Jmse 10 00968 g010aJmse 10 00968 g010b
Figure 11. Spatial and temporal distribution of mean grain-size from 1989 to 2015.
Figure 11. Spatial and temporal distribution of mean grain-size from 1989 to 2015.
Jmse 10 00968 g011aJmse 10 00968 g011b
Figure 12. Changes in the mean grain-size from 1989 to 2015 along with the discharge of the Yellow River and precipitation.
Figure 12. Changes in the mean grain-size from 1989 to 2015 along with the discharge of the Yellow River and precipitation.
Jmse 10 00968 g012
Figure 13. Validation results of the optimal BPNN model for clay (a), silt (b), sand (c), and mean grain-size (d).
Figure 13. Validation results of the optimal BPNN model for clay (a), silt (b), sand (c), and mean grain-size (d).
Jmse 10 00968 g013
Table 1. Specification of the Landsat imagery used in this study.
Table 1. Specification of the Landsat imagery used in this study.
DateSensorTime (UTC)Track/PathWater Level (m NAP) Tidal Stage
13 February 1989 TM02:14121/34−1.60incoming
26 March 1995TM01:51121/34−1.35outgoing
10 December 1999ETM+02:34121/34−1.77incoming
12 October 2001ETM+02:29121/34−1.80outgoing
11 November 2006ETM+02:31121/34−1.55incoming
3 November 2009ETM+02:32121/34−1.85incoming
23 August 2012ETM+02:36121/34−1.93incoming
5 February 2015ETM+02:40121/34−2.05incoming
Tidal information accessed from 13 February 1989 to 5 February 2015 is from China Oceanic Information Network (https://ocean.cnss.com.cn/, accessed on 1 January 2022).
Table 2. Correlations between the band indices and the grain-size parameters.
Table 2. Correlations between the band indices and the grain-size parameters.
Spectral IndicesExpressionsCorrelation Coefficients(r)
Mean Grain-SizeSandSiltClay
rmaxrminrmaxrminrmaxrminrmaxrmin
Roriginal reflectance0.55 *0.020.58 *0.070.58 *0.070.52 *0.07
SIRi + Rj0.48 *0.080.49 *0.090.48 *0.090.43 *0.04
DIRiRj0.65 *00.71 *0.040.70 *0.040.64 *0.01
PIRi × Rj0.43 *0.080.44 *0.080.43 *0.080.40 *0.04
RIRi/Rj0.52 *0.010.59 *00.59 *0.010.50 *0.01
NDI(RiRj)/(Ri + Rj)0.1400.15 *00.16 *00.100
* denotes a significant level (p < 0.05).
Table 3. Statistics of the SVM mapped grain-size parameters from 1989 to 2015.
Table 3. Statistics of the SVM mapped grain-size parameters from 1989 to 2015.
Clay (%)Silt (%)Sand (%)Mean Grain-Size (φ)Discharge
(109 m3)
Precipitation
(mm)
MaxMinMeanMaxMinMeanMaxMinMeanMaxMinMean
19899.651.355.3293.029.2559.6293.513.1948.786.793.254.05568.80369
19955.474.405.0376.5230.8254.1872.1820.4045.655.402.484.11567.90728
19995.474.865.0185.1530.0959.0370.8012.3538.995.862.304.1061.69373
20015.623.855.0278.4919.9855.4676.8215.2341.455.542.484.1440.89414
20065.963.185.0384.5418.8854.1380.259.4843.235.511.404.01186.70452
20096.252.624.8983.7621.7551.1275.727.8945.575.912.544.27132.90653
20126.571.614.8975.769.5152.0392.9319.6444.275.652.374.11282.50534
20155.072.894.6475.655.6037.4294.1916.8658.465.471.944.10133.60595
Table 4. Optimal node number combinations and calibration accuracies of BPNN model.
Table 4. Optimal node number combinations and calibration accuracies of BPNN model.
Grain-Size ParametersHidden Layer NumberNode NumberR2RMSE
1st2nd
clay content26110.671.95
silt content216160.78101.27
sand content211160.79120.53
mean grain-size217180.770.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, X.; Zhan, C.; Liu, Y.; Bi, J.; Li, G.; Cui, B.; Wang, L.; Liu, X.; Wang, Q. Retrieval of Remotely Sensed Sediment Grain Size Evolution Characteristics along the Southwest Coast of Laizhou Bay Based on Support Vector Machine Learning. J. Mar. Sci. Eng. 2022, 10, 968. https://doi.org/10.3390/jmse10070968

AMA Style

Yu X, Zhan C, Liu Y, Bi J, Li G, Cui B, Wang L, Liu X, Wang Q. Retrieval of Remotely Sensed Sediment Grain Size Evolution Characteristics along the Southwest Coast of Laizhou Bay Based on Support Vector Machine Learning. Journal of Marine Science and Engineering. 2022; 10(7):968. https://doi.org/10.3390/jmse10070968

Chicago/Turabian Style

Yu, Xiang, Chao Zhan, Yan Liu, Jialin Bi, Guoqing Li, Buli Cui, Longsheng Wang, Xianbin Liu, and Qing Wang. 2022. "Retrieval of Remotely Sensed Sediment Grain Size Evolution Characteristics along the Southwest Coast of Laizhou Bay Based on Support Vector Machine Learning" Journal of Marine Science and Engineering 10, no. 7: 968. https://doi.org/10.3390/jmse10070968

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop