Detection of the Minute Variations of Total Suspended Matter in Strong Tidal Waters Based on GaoFen-4 Satellite Data

: Hangzhou Bay (HZB) is the largest macro-tidal bay in China, where suspended sediment concentrations are signiﬁcantly modulated by tidal oscillations. This makes it an ideal area for the study of the impact of tide on temporal–spatial variation in suspended sediment. The GaoFen-4 (GF-4) satellite is the ﬁrst high-resolution geosynchronous orbiting satellite of China. It exhibits the unique advantages of capturing minute variations and ﬁner details of total suspended matter (TSM) due to the enhancement in spatial resolution (50 m) and observation time interval (20 s). In this study, TSM concentration of the HZB was retrieved based on the GF-4 satellite. The spatial distribution and minute variations of TSM concentration under the ebb tide from 7:30 to 7:40 a.m. on 28 August 2017, were analyzed. The results showed that the average TSM concentration inside HZB was (371.8 ± 1.8) mg/L. There was a linearly increasing trend of TSM concentration at ebb tide, with an increment of (3.96 ± 0.31) (mg/L)/min, and a more signiﬁcant increase was observed in the high TSM areas. This increase in TSM concentration was associated with both the bottom topography and tide processes. The tidal potential energy generated by the tidal range and the strong shear stress generated by the high current velocity both led to the re-suspension of the sedimentary particles, which affected the variation of TSM concentration. In addition, the inﬂuence of bottom topography changed the intensity of re-suspension and also affected the distribution of TSM concentration in HZB.


Introduction
Suspended sediments can affect the water quality, ecosystem, and bottom geomorphologic evolution in estuaries and bays [1,2]. They are an important indicator of the sediment transport and resuspension process [3,4]. The exchanges of particles such as sediment can be influenced by the tide in estuaries and bays [5]. Total suspended matter (TSM) concentration therefore provides great reference for tidal prediction, disaster warning, marine engineering construction, and underwater equipment design in strong tidal area.
The study of the inter-annual, annual, and monthly variations of TSM concentration in strong-tidal estuary, gulf, and other sea areas has been widely reported in previous literatures [6][7][8][9][10][11]. Using traditional in-situ measurement, it is difficult to detect the TSM concentration on a timely basis and on a large scale [12,13]. The satellite can provide global water leaving radiance from the sea surface, so as to provide continuous observation of TSM [11,14]. Early TSM concentration was mostly retrieved by Landsat-x_MSS [6,15], NIMBUS-7_CZCS [16], and SeaStar_SeaWiFs [8,17] satellite sensors. With the progress of optical sensor design technology, the MODIS [18,19], MERIS [9,10,20], sentinel-2A/B [11], sentinel-3_OLCI [21], and VIIRS [22] have been used widely to study the spatial and temporal change trends of TSM concentration in near shore and in large lakes. However, for highly turbid waters such as estuaries and coastal areas with significant hourly-varying characteristics, polar orbit satellites with high spatial resolution and low temporal resolution cannot meet the dynamic monitoring requirements.
In recent years, a geostationary orbit satellite has been introduced into the research. The Geostationary Ocean Color Imager (GOCI) satellite has been used to describe the diurnal variation characteristics of TSM concentration in estuaries and inshore and coastal waters [23][24][25][26]. However, the spatial resolution GOCI is 500 meters at nadir, and the image of the GOCI satellite only covers the partial sea area centered on the Korean peninsula. It is not easy to use GOCI to study the TSM distribution in small bays and coastal regions. In order to improve the low spatial resolution of GOCI image, Pan et al. [27] proposed to conduct multi-source data fusion between Landsat OLI (30 m, 16 days) and GOCI (500 m, 1 h), and the fusion results could more accurately reflect the diurnal variation characteristics of TSM. Furthermore, water conservancy, waterway, and sea reclamation projects have put forward higher requirements for the understanding of the basic law of suspended particulate transport in high-tide waters, which should have higher-frequency variations that are sometimes located in a small region. However, at present, no research has been proposed to describe higher temporal variability of TSM concentration than GOCI.
Higher temporal and spatial resolution are needed to advance the understanding of variation processes of TSM concentration under the impact of tide. GaoFen-4 (GF-4), launched in December 2015, is the first high-resolution geosynchronous orbiting satellite of China. The GF-4 satellite provides a finer resolution observation at 50 m spatial resolution. Due to the enhancement in spatial resolution, it allows the identification of the finer details of TSM in estuaries and bays. The GF-4 observes the earth surface every 20 s. This makes it possible to study high-frequency changes of TSM [28][29][30]. Furthermore, the spatial coverage of the GF-4 satellite extends to more estuaries and bays of China. In this paper, the TSM distribution and the current in the Hangzhou bay (HZB) were firstly retrieved from the GF-4 top of atmosphere (TOA) reflectance images. Then, the minute variability of TSM concentration and its mechanism was analyzed.

Study Area
The HZB, located in the northeast of Zhejiang province (29-31 • N, 121-123 • E), is the largest marginal bay along the southeast coast of China and an important transition zone between sea and land connecting the East China Sea, as shown in Figure 1. The bay is wide and shallow, with a total area of about 8500 km 2 , an average water depth of 8-10 m, a length of about 86 km, and a rapid change in width, from 100 km at the mouth of the bay to only 20 km at its narrowest point [31]. HZB is part of the mouth of the Qiantang River, to the north of which lies China's longest river, the Yangtze River. The Qiantang river and the adjacent Yangtze river in the north are the main freshwater resources entering the HZB, where the freshwater runoff is 42 and 925 km 3 /a, respectively. Due to the erosion of these two river basins, millions of tons of sediments are discharged into the bay every year. Among them, the Qiantang river receives about 7.9 × 10 6 tons of sediment each year [19], and 3.1% of the total sediment from the Yangtze river to the East China Sea is carried into the HZB by the south branch of the Yangtze river, which receives about 4.8 × 10 8 tons annually [32]. Under the influence of the complex dynamic processes, such as the south branch current of the Yangtze estuary, the runoff of the Qiantang River, and the tide of the East China Sea, the study area has typical characteristics of a strong tidal area. The tidal type of HZB is a regular semi-diurnal tide. The M 2 constituent (Table 1) is the dominant tidal component of HZB. Its seabed is an undulating and shrinking plane, and the tidal range increases rapidly from east to west. The average tidal range is 3-4 m at the mouth of the bay and 4-6 m within the bay, and the tidal peak at the mouth of the Qiantang River can reach 8.93 m. HZB is one of the most remarkable tidal bays in the world [13,32,33].
Due to strong tides and rapid water flow, a large number of suspended particles are deposited in HZB [34], which is characterized by significant periodic changes in sediment concentration [23,25,26]. Within a tidal cycle, the significant changes of TSM occur in the middle and later periods of high tide and low tide [24]. and the tidal range increases rapidly from east to west. The average tidal range is 3-4 m at the mouth of the bay and 4-6 m within the bay, and the tidal peak at the mouth of the Qiantang River can reach 8.93 m. HZB is one of the most remarkable tidal bays in the world [13,32,33].
Due to strong tides and rapid water flow, a large number of suspended particles are deposited in HZB [34], which is characterized by significant periodic changes in sediment concentration [23,25,26]. Within a tidal cycle, the significant changes of TSM occur in the middle and later periods of high tide and low tide [24].

Data Sources
The GaoFen series satellites, which were designed to collect land atmosphere and ocean high resolution imageries and to form all-weather and all-day global coverage of earth observation capability by China with 13 satellites, launched from 2013. The GF-4 satellite was launched in 2015 with the sensors of a staring camera, which has a 50 m' spatial resolution in visible or near infrared bands (bands 1-5) or 400 m' spatial resolution in the thermal infrared band (band 6). The GF-4 satellite is a geostationary satellite, and the temporal resolution is 20 s. The width of each GF-4 image is 40 km. Gf-4 data have a super high spatial and temporal resolution, which is suitable for the high frequency change study of TSM in the high turbidity sea area. Due to the limitations of transmission and data storage technology, the time interval of the data products provided is about 1 min. In this paper, we used the GF-4 satellite PMS L1A data from 7:30:00 a.m. to 7:40:00 a.m. on 28 August 2017, provided by China Centre for Resources Satellite Data and Application.
The field observation data were measured at the cruise during 2 to 13 December 2011, and 16 to 22 May 2014. The spectra of in-situ R rs were obtained by using a hand-held analytical spectral devices (ASD) spectrometer and a standard board of 30% reflectance. The in-situ TSM concentration was collected by a standard sampler and measured by the filtration gravimetric method. From each station, two parallel samples were collected and the TSM concentration was measured, respectively. The average TSM concentration of the two parallel samples was derived as the in-situ TSM concentration of the station. A blank sample was run with each sample to blank correct the data. The 71 effective TSM concentration data were finally used to calibrate the TSM inversion algorithm. The TSM was grouped with the 41 TSM data for SAI (spectral absorption index) model construction and the 30 TSM data for model validation.
The topographic map is the global topographic relief model ETOPO published by U.S. National Geophysical Data Center, which is the highest resolution known topographic relief data [34]. It includes three data with different accuracy, ETOPO1, ETOPO2, and ETOPO5, and there are two versions, the Ice Surface and the Bedrock, the former being the elevation of the ice sheet and the latter giving the elevation of the bedrock. In this paper, the Bedrock version in ETOPO1 (resolution 1 min of arc) with the highest accuracy was selected to map the bottom terrain of HZB.
Regarding the tide level data, looking up the Tide Table [35], the hourly tide level of each tide station during the collection period of GF-4 PMS images (7:00-8:00 a.m. on 28 August 2017) was obtained, and the difference value of the tide level in two adjacent hours was taken as the change of the tide level. Then, the minute change of tide level of the station was obtained after averaging.

GF-4 Remote Sensing Data Processing
The most common atmospheric correction method for ocean color remote sensing is the algorithm proposed by Gordon, which used the near-infrared (NIR) band as the dark pixel due to its good absorbability [14]. However, turbid waters in coastal regions have a strong scattering ability in the NIR band, so the hypothesis of dark pixels above is invalid [36,37]. Additionally, the algorithm using the short wave infrared (SWIR) band instead of the NIR band has been proven to be suitable for turbid waters because of their strong absorbability in the SWIR band [38][39][40]. Since the GF-4 satellite has no SWIR band, the SWIR-based dark pixel method cannot be used for its atmospheric correction directly. The Visible Infrared Imaging Radiometer Suite (VIIRS), which is on board the Suomi National Polar-Orbiting Partnership (SNPP) satellite, has three SWIR bands with band centers 1238, 1601, and 2257 nm. The quasi-synchronous VIIRS (QSV) data can be used for the atmospheric correction of other sensors in highly turbid waters [41].
The R rs of corrected reflectance was derived after correction for Rayleigh scattering, gaseous absorption, and aerosol multiple scattering effects following: where L λ is the TOA radiance measured by the sensor in band λ, d is the Sun-Earth distance at the time of the measurement showing the relative deviation from the average Sun-Earth distance d 0 equal to 1 AU, θ s is solar zenith angles, F 0 is the extraterrestrial solar irradiance for mean Sun-Earth distance d 0 , R r,λ is the reflectance of Rayleigh scattering, R a,λ is the reflectance of aerosol multiple scattering reflectance, and t v,λ is the diffuse transmittance of the atmospheric column. Among them, the Second Simulation of the Satellite Signal in the Solar Spectrum (6S) model [42] was used to estimate Rayleigh reflectance (R r,λ ). The necessary input parameters for the 6S model were image-specific geometrical conditions, atmospheric model for gaseous components, and spectral condition. The estimation of aerosol reflectance (R a,λ ) followed the QSV algorithm. First, the aerosol types and the aerosol optical thickness were estimated using the three shortwave infrared bands of VIIRS, and then the aerosol reflectance at GF-4 observing geometries according to the aerosol properties was calculated. Note that the aerosol look-up-tables (LUTs) for 80 aerosol types constructed by Ahmad et al. [43] are necessary for the QSV algorithm, but these LUTs for the GF-4 are not provided by SeaDas software (https://seadas.gsfc.nasa.gov/, accessed on 2 March 2021); thus, the aerosol reflectance and the diffuse transmittance in the GF-4 bands were estimated using the aerosol LUTs for Hyperspectral Imager for the Coastal Ocean (HICO, 400-900 nm, 128 bands) and the response function of the GF-4 bands. The spectral response function of GF-4 bands was provided by China Centre for Resources Satellite Data and Application.

The Correction of TSM Algorithm
GF-4 satellite data were used to retrieve TSM concentration. It was estimated based on multi-band remote sensing reflectance combined with spectral index with reference to the SAI algorithm suggested by Liu et al. [25]. This algorithm was based on the measured TSM concentration data in May 2013 and the GOCI images of the same day. It had been verified that it can realize the inversion of TSM concentration in HZB [25].
In this paper, the bands R rs (519), R rs (550), and R rs (770) of the GF-4 were used to replace the bands R rs (490), R rs (555), and R rs (745) of the GOCI in original algorithm (Table 2). Since the bands of the GF-4 satellite are different from the corresponding bands of the GOCI satellite in the spectral range and central wavelength, the applicability of the GF-4 bands in the SAI algorithm needs to be verified. According to the spectral range of each band of GF-4 and GOCI satellites, the spectral response function was used to conduct equivalent band processing on the in-situ R rs . The R rs of corresponding bands between the GOCI satellite and GF-4 satellite was compared, and the root mean square error (RMSE) and correlation coefficient (r) were used to evaluate the verification results. The verification results showed (Figure 2) that the band of the GOCI and GF-4 satellite had a good fitting result, that is, the bands of the GF-4 can be used for the SAI algorithm. The correlation coefficient of each corresponding band was greater than 0.95. The R rs can reach up to 0.08 1/Sr in NIR band (Figure 2c), which was consistent with the characteristics of turbidity water. The RMSE of each band increased with the increase of wavelength, that is, the uncertainty of R rs increased with the increase of wavelength, which may be related to the strong backscattering phenomenon of suspended particulate matter in the NIR band in turbidity water. Although the verification results showed that the bands of the GF-4 were very similar to the corresponding bands of GOCI, considering the seasonal characteristics of the spatial distribution of TSM concentration in HZB, only the structure of the SAI algorithm was referred, and then the two cruises in-situ data on 2 to 13 December 2011, and 16 to 22 May 2014, were used to correct the coefficients of this algorithm. The in-situ data contained not only TSM, but also the remote sensing reflectance. The in-situ TSM data were grouped with the 41 TSM data for SAI model correction and the 30 TSM data for model validation. The two sets of TSM data have similar numerical ranges, 60-2036 and 156-2238 mg/L, respectively. This large numerical range of TSM data used to develop the algorithm can reflect the response of remote sensing reflectance-derived SAI to TSM variations under different water conditions and concentration levels. These will lead to the algorithm being generally stable and usable.
The corrected SAI algorithm is shown in Figure 3a. TSM concentration was calculated by using Equation (2): The corrected SAI algorithm is shown in Figure 3a. TSM concentration was calculated by using Equation (2)   The result of model validation was shown as Figure 3b. It was shown that the C TSM calculated by SAI algorithm was consistent with the in-situ C TSM , that is, this algorithm can be used for the inversion of C TSM in HZB by the GF-4 satellite. The result of model validation was shown as Figure 3b. It was shown that the CTSM calculated by SAI algorithm was consistent with the in-situ CTSM, that is, this algorithm can be used for the inversion of CTSM in HZB by the GF-4 satellite.  The satellite derived TSM concentration within the range from Ganpu tide station to Lu Chaogang tide station (120 • 54 E-121 • 50 E) was averaged as the mean TSM concentration in HZB. Similarly, the TSM concentration within the range from Zhongjun tide station to Shengshan tide station (121 • 54 E-122 • 48 E) was averaged as the mean TSM concentration outside HZB.
The TSM difference of two GF-4 remote sensing images with adjacent time was calculated, and the mean value of the difference was calculated as the minute variations of TSM in HZB. The calculation formulas are shown as Equation (3): where ∆C TSM is the average difference of TSM concentration at two adjacent times (mg/L); C TSM t+1 and C TSM t are the TSM concentration at the adjacent times (mg/L). The current velocity was calculated by the tracer measurement method [44] based on the multi-temporal images of the GF-4 satellite. The tracer measurement method takes TSM as the tracer and calculates its moving distance under the action of water flow in a period of time, thus obtaining the velocity of current. Assuming that the time interval of image acquisition is T and the distance of the tracer TSM moving along two mutually perpendicular directions ( → x , → y ) is ∆x and ∆y, respectively, then the current velocity can be expressed as Equation (4): where the moving distance ∆x and ∆y was calculated by the maximum likelihood method [44][45][46]. The algorithm is as follows: a window of a × a was selected based on the point f 1 x i , y j of the image at the initial moment, and the correlation coefficient was calculated by traversing the image at the end moment. Taking the maximum correlation coefficient as the matching coefficient of window a × a, then the corresponding point f 10 x i + ∆x, y j + ∆y was regarded as the position after moving. The matching coefficients r a,i,j of different windows (a = 11, 13, · · · , 51) were calculated, respectively. The window a × a corresponding to the maximum matching coefficient r max was the optimal window; ∆x and ∆y corresponding to the maximum matching coefficient were the moving distance of the point f 1 x i , y j . If the maximum matching coefficient is less than 0.5, the current velocity at this point is considered invalid. The calculation formula is Equation (5): , a = 11, 13, · · · , 51, r a,i,j is the matching coefficient between f 1 x i , y j and f 10 x i + ∆x, y j + ∆y in window a × a, r max is the maximum matching coefficient, f 1 x i , y j is TSM concentration of the coordinate at the initial moment, f 10 x i + ∆x, y j + ∆y is TSM concentration after the point has moved a distance at the end moment, ∆x and ∆y are the displacement values in the horizontal direction and the vertical direction of the point, and f 1 and f 10 are the average TSM concentration values in this window.

The Spatial Distribution Characteristics of TSM in HZB
The TSM concentration inside HZB exhibited a "high-low-high" pattern from east to west, whereas it showed a slow decreasing trend at the mouth of HZB (Figure 4a). The average TSM concentration inside HZB ((371.8 ± 1.8) mg/L) was significantly higher than that at the mouth of HZB ((109.6 ± 3.0) mg/L), and there were four high TSM concentration zones and two low TSM concentration zones inside HZB at ebb tide, which was basically consistent with the results of previous studies [24,25]. Among them (Figure 4b), the high TSM concentration zone H1 was located in the area of the northeast coast of HZB and near the mouth of the bay, spanning 0.5 • latitude and showing a feathery distribution along the coast. In the center of H1, there was a relatively narrow low concentration transition zone from east to west. The average TSM concentration in the north (H1_a) was (486.6 ± 2.8) mg/L, which was higher than that the (437.3 ± 3.1) mg/L in the south (H1_b), and the high TSM concentration area was bigger than that in the south. In addition, the other three high value zones were, respectively, distributed in the neck of the bay (H3) and the south bank area (H2, H4). At the neck of HZB, TSM was distributed in a strip along the northern shore of the bay, with an average TSM concentration of (384.2 ± 2.3) mg/L. The two high-value zones near the south bank were distributed in a cluster with an average TSM concentration of (526.0 ± 3.7) and (461.0 ± 2.6) mg/L, respectively. There was a large gradient at the edge of the high value zones. The low TSM concentration was concentrated in L1 and L2 zones. L1 was located on the northwest coast near the mouth of the bay, with an average TSM concentration of (226.5 ± 1.7) mg/L. L2 was distributed between Haihuangshan tide station and Zhoushan islands, with an average TSM concentration of (158.6 ± 1.1) mg/L.
Based on the tidal process and the topography of the bottom of HZB (Figure 4a), it was found that TSM concentration had a negative correlation with the topography within the water depth of 10 m at ebb tide. The average depth of HZB is 10 meters. Among them, the shallow water areas in the bay neck, the northeast coast near the bay mouth, and the south coast are generally less than five meters, corresponding to H3, H1a, H2, and H4 of the high TSM concentration areas, respectively. In addition, there is a deep-water tidal channel with a depth of 5-10 meters in the upper reaches of the north bank and the middle part of the south bank, respectively, corresponding to L1 and L2 in the low TSM concentration zones.

Minute Scale Dynamic Change Characteristics of TSM at Ebb Tide in HZB
Based on the temporal and spatial variation of TSM concentration in minute scale ( Figures 5 and 6), it was found that TSM in HZB generally increases at ebb tide, and a more obvious increase was observed in the area with high TSM concentration. In the center of the low TSM concentration region, the concentration increased slightly. In order to quantitatively analyze the law of incremental change of TSM, linear regression analysis was performed on the average TSM concentration at ebb tide by taking the time interval as the independent variable and the TSM concentration as the dependent variable. According to the analysis results, TSM concentration showed a significant linear growth trend at ebb tide. The linear fitting formula was as Equation (6) (Figure 7a): where the intercept is TSM concentration at the initial moment, and the slope coefficient represents the increment of TSM concentration per minute of about (3.96 ± 0.31) mg/L (the increment of TSM concentration per hour was (238 ± 19) mg/L). R 2 is 0.99 and p < 0.01, indicating that the time variable could explain the 99% change in TSM concentration at ebb tide.      The in-situ TSM concentration at the four sampling stations during the 2010 cruise is shown in Figure 8. The location of the sampling station is shown in Figure 8a, and the TSM concentration variations of the four sampling stations are shown in Figure 8(b1-b4). The left axis of the graph represents tide level, and the right axis represents TSM concentration. The hourly TSM concentration increment at these four stations was calculated on average as the average increment per hour in HZB. The calculated result was (256 ± 12) mg/L, which was basically consistent with the GF-4 fitting results (the increment of TSM concentration per minute of about (3.96 ± 0.31) mg/L, that is, the increment of TSM concentration per hour was (238 ±19) mg/L).   Figure 7a is the situation inside HZB. The two transects at the neck of HZB (120 • 50 E, 121 • E) accounted for the largest proportion (43%), that is, the increment of TSM concentration in these regions directly affected the increment of TSM concentration inside HZB. In addition, the average TSM concentration in these two regions ((655.5 ± 5.7) mg/L, (507.4 ± 7.4) mg/L) was significantly higher than the total average TSM concentration inside HZB ((371.8 ± 1.8) mg/L), indicating that the variation of TSM concentration in the high TSM concentration region played a dominant role in the variation of TSM concentration inside HZB (Figure 7b). However, the increment of the five transects (121 • 10 E, 121 • 20 E, 121 • 30 E, 121 • 40 E, and 121 • 50 E) across the deep-water tidal channel was relatively small, which made a small contribution to the increase of TSM concentration at ebb tide inside HZB. The right side of the Figure 7a was the TSM concentration variation outside HZB. The increment of TSM concentration showed little change. As the ebb tide continued, the TSM concentration difference between transects of the bay mouth gradually narrowed at the late ebb tide, and the average TSM concentration increased slightly, which was less than 5 mg/L.

Discussion
The distribution and variation of TSM are influenced by many factors, such as monsoon, freshwater input, tidal processes, and bottom topography [25,47]. In this study, during the remote sensing image acquisition period, there was no typhoon or other storm in the study area. According to the characteristic that the mean wind speed and mean freshwater input are basically unchanged within a small interval, the influence of wind speed and freshwater input on TSM within 10 min was ignored.
The tidal potential energy generated by the tidal range and the strong shear stress generated by the high current velocity both lead the sediment to be disturbed and resuspended in the water, which may be the main reason for the increase of TSM concentration in HZB [5,23]. The greater the tidal range and the greater the potential energy of the tidal range, the greater the change of TSM concentration that will be caused [24]. Additionally, the higher velocity of current produces a strong shear on the seabed, which also makes the TSM concentration increase [23,48]. With the decrease of tidal range and current velocity, the re-suspension of sediment became weaker, and the TSM concentration decreased gradually. As shown in Figure 9, the variation of tide level and velocity were different at different locations in Hangzhou Bay. The change of tide levels inside HZB was mostly more than 1 cm/min, and the average change was about twice that of the mouth of the bay (Figure 9a). The maximum range of tide appears near Ganpu tide station at the bay neck (2.04 cm/min), which corresponded to the neck of HZB where TSM changed the most (mentioned in Section 3.2). The change of tide levels at the mouth of HZB was small, and the variation is mostly within the range of 0.5 to 1 cm/min. The TSM was also less variable in this water. This was consistent with the phenomenon mentioned above that the greater the tidal range, the greater the TSM concentration change. Similarly, the current velocity was the largest in the sea area near the Ganpu Tide Station and flew to the northeast and slowed down. This was also consistent with the phenomenon mentioned above that the greater the current velocity, the greater the TSM concentration change. However, in the H1, H2, and H4 area, the lower tidal range and current velocity corresponded to higher TSM concentration, and in the L1 and L2 area, the higher tidal range and current velocity corresponded to lower TSM concentration. This phenomenon was contrary to the sediment lifting effect formed by higher tidal range and current velocity. The reason may be due to the difference of bottom topography, which affected the strength of the vertical sediment lifting action [4,[49][50][51]. The tidal level changes in However, in the H1, H2, and H4 area, the lower tidal range and current velocity corresponded to higher TSM concentration, and in the L1 and L2 area, the higher tidal range and current velocity corresponded to lower TSM concentration. This phenomenon was contrary to the sediment lifting effect formed by higher tidal range and current velocity. The reason may be due to the difference of bottom topography, which affected the strength of the vertical sediment lifting action [4,[49][50][51]. The tidal level changes in Luchaogang, Haihuangshan, and Zhenhai tide stations were small (tidal range < 1.25 cm/min). However, as they were located on the shallow tidal flats, the potential energy generated by the smaller tidal range was sufficient to reach the sediment, resulting in a greater stress on the sediment and a stronger resuspension process. Relative to the other three tide stations (Zhapu, Jinshanzui, and Tanhutide stations) inside HZB, the tide level changed significantly (tidal range > 1.25 cm/min). However, they were all located above the deep-water tidal channel on the north shore of HZB, the potential energy formed by the tidal range was consumed before reaching the seafloor, and the stress formed at the bottom was small and the resuspension process caused by it was weak. In addition, the two low TSM concentration zones L1 and L2 were located at the deep-water tidal channel where the current velocity was maximum. Although a strong sediment re-suspension process may result in an increase of TSM concentration near the seabed, the concentration of the whole water was still at a low TSM level, and as the velocity slowed down, these suspended particles were more likely to settle again. The current velocity of the H1_a, H2, and H4 concentration zones was small. However, due to the shoal location, the shear stress required for particulate re-suspension was relatively small. Therefore, the small flow velocity was sufficient to suspend the deposited particulate and caused TSM to increase. Zeng et al. [25] showed that the large sandbar and tidal flat where the TSM concentration high-value zone was located in HZB had more silt, and most of it was non-viscous fine sediment. Vertical stirring was more likely to cause the resuspension of sediment particles. Xie et al. [32] found that at the bottom of the two tidal channels in the north and south of HZB, there was a layer of viscous sediment with a thickness of about 10 m, which could be re-suspended with less sediments and a higher critical stress required for resuspension. The influence of different sediment types on vertical sand lifting also explained the phenomenon of higher tidal range and current velocity corresponding to lower TSM concentration and lower tidal range and current velocity corresponding to higher TSM concentration in HZB.
In recent years, the development and utilization of new clean energy such as tidal power have been paid more and more attention and become an important means to achieve sustainable development. The clarification of TSM minute scale variation in HZB is of great significance for understanding the complex hydrodynamic processes of strong tides and high turbidity sea areas as well as the physical and biogeochemical processes of the waters. In this paper, the variation trend and influencing factors of TSM minute scale in HZB were analyzed, and the relationship between TSM concentration variation and topography, tidal level, and current velocity was described in detail. It provided a more reliable reference for tidal prediction and disaster warning in strong tidal areas, and a reference for the selection of construction technology, equipment design, quality control, and safety assurance of marine engineering in strong tidal areas and other similar environments. However, in this paper, due to the lack of synchronous in-situ data, the accuracy of the TSM inversion algorithm still needs to be improved. In the analysis of the factors affecting TSM concentration, only qualitative analysis was carried out. With the further development of high spatial and temporal resolution satellite technology, the distribution characteristics and variation rules of TSM concentration at different scales can be studied in depth. It is an important direction of our current research to understand and master the smaller scale variation characteristics of TSM concentration and its dynamic mechanism.

Conclusions
Based on the study of remote sensing of TSM concentration in HZB using GF-4 satellite data and the analysis of TSM concentration variation, the following conclusions were drawn in this paper: (1) The average TSM concentration in HZB was (371.8 ± 1.8) mg/L during the study period. There were four high TSM concentration zones and two low TSM concentration zones. The spatial distribution was affected by the bottom topography, and the high TSM zones were located on the tidal flat and bar with relatively shallow water depth, while the low TSM zones were located at the deep-water tidal channel.
(2) At ebb tide in HZB, TSM concentration showed a linear growth trend with an increase of about (3.96 ± 0.31) (mg/L)/minute. The concentration in the high TSM concentration area increased more obviously, and the variation of those high TSM concentration regions had a great influence on the TSM variation inside HZB. The concentration in the center of the low-TSM concentration regions increased slightly and contributed little to the variation of TSM.
(3) In HZB, the tidal process may be the main factor that affected TSM variation in a short time. The tidal potential energy generated by the tidal range and the strong shear stress generated by the high current velocity both lead to the re-suspension of the sedimentary particles and caused the concentration of TSM to increase. The bottom topography influenced the intensity of the tidal process and may be one of the reasons for TSM variation.