Next Article in Journal
MDANet: A High-Resolution City Change Detection Network Based on Difference and Attention Mechanisms under Multi-Scale Feature Fusion
Next Article in Special Issue
Comparative Evaluation of Semi-Empirical Approaches to Retrieve Satellite-Derived Chlorophyll-a Concentrations from Nearshore and Offshore Waters of a Large Lake (Lake Ontario)
Previous Article in Journal
Automatic Building Roof Plane Extraction in Urban Environments for 3D City Modelling Using Remote Sensing Data
Previous Article in Special Issue
Characterizing Water Composition with an Autonomous Robotic Team Employing Comprehensive In Situ Sensing, Hyperspectral Imaging, Machine Learning, and Conformal Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Retrieval of Total Suspended Matter Concentration Based on the Iterative Analysis of Multiple Equations: A Case Study of a Lake Taihu Image from the First Sustainable Development Goals Science Satellite’s Multispectral Imager for Inshore

1
College of Geoscience and Surveying Engineering, China University of Mining and Technology–Beijing, Beijing 100083, China
2
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
3
Beijing Institute of Space Mechanics and Electricity, Beijing 100083, China
4
College of Resources Environment and Tourism, Capital Normal University, Beijing 100048, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(8), 1385; https://doi.org/10.3390/rs16081385
Submission received: 10 March 2024 / Revised: 5 April 2024 / Accepted: 11 April 2024 / Published: 14 April 2024
(This article belongs to the Special Issue Remote Sensing Band Ratios for the Assessment of Water Quality)

Abstract

:
Inland waters consist of multiple concentrations of constituents, and solving the interference problem of chlorophyll-a and colored dissolved organic matter (CDOM) can help to accurately invert total suspended matter concentration ( C tsm ). In this study, according to the characteristics of the Multispectral Imager for Inshore (MII) equipped with the first Sustainable Development Goals Science Satellite (SDGSAT-1), an iterative inversion model was established based on the iterative analysis of multiple linear regression to estimate C tsm . The Hydrolight radiative transfer model was used to simulate the radiative transfer process of Lake Taihu, and it analyzed the effect of three component concentrations on remote sensing reflectance. The characteristic band combinations B 6 / 3 and B 6 / 5 for multiple linear regression were determined using the correlation of the three component concentrations with different bands and band combinations. By combining the two multiple linear regression models, a complete closed iterative inversion model for solving C tsm was formed, which was successfully verified by using the modeling data (R2 = 0.97, RMSE = 4.89 g/m3, MAPE = 11.48%) and the SDGSAT-1 MII image verification data (R2 = 0.87, RMSE = 3.92 g/m3, MAPE = 8.13%). And it was compared with iterative inversion models constructed based on other combinations of feature bands and other published models. Remote sensing monitoring C tsm was carried out using SDGSAT-1 MII images of Lake Taihu in 2022–2023. This study can serve as a technical reference for the SDGSAT-1 satellite in terms of remote sensing monitoring of C tsm , as well as monitoring and improving the water environment.

Graphical Abstract

1. Introduction

Inland lakes are essential freshwater storage and supply sources [1] and play an indispensable role in satisfying human drinking water needs, in providing agricultural irrigation and industrial water, etc. However, with the rapid development of climate change, industrialization, and urbanization, the water quality and ecological environment of inland lakes are increasingly threatened [2]. In response to this challenge, water quality monitoring is seen as an effective means to protect freshwater resources. Among them, total suspended matter is an important indicator for measuring water quality and the level of pollution in water quality monitoring. It mainly consists of solid substances suspended in water such as insoluble organic matter, inorganic matter, sediment, microorganisms, and clay, which have a direct impact on light transmission in water, aquatic vegetation growth, and the primary productivity of water bodies, as well as being closely related to the ecological health of water environment [3,4]. Therefore, monitoring total suspended matter concentration ( C tsm ) in water bodies is of great significance for understanding the current situation and development trend of water quality, in addition to managing water pollution [5].
Satellite remote sensing technology has been widely used in the field of water environment monitoring in recent years, due to its capacity to acquire important information about water on a regular, timely, and synchronous basis [6]. To improve the effectiveness of remote sensing monitoring, high spatial resolution is typically required [7]. The first Sustainable Development Goals Science Satellite (SDGSAT-1) was developed by the International Research Center of Big Data for Sustainable Development Goals (CBAS). The satellite is equipped with the Multispectral Imager for Inshore (MII), which has a spatial resolution of 10 m, a width of up to 300 km, and seven different spectral bands, making it highly effective in water environment monitoring. It has unique advantages in the detection of C tsm [8], which provides strong support for remote sensing monitoring of C tsm .
The optical properties of inland water bodies are quite complex [9], and there are a variety of optically active components in water bodies, among which total suspended matter, chlorophyll-a, and colored dissolved organic matter (CDOM) are three major optically active substances that affect the optical properties of water bodies [10], and the absorption and scattering of these components jointly affect the reflectivity of water bodies. At present, scholars at home and abroad have developed a variety of remote sensing methods for retrieving suspended matter concentration, mainly including analytical models, empirical models, and semi-empirical/semi-analytical models [11,12,13]. Among them, the analytical model is based on the theory of radiative transfer in water bodies and uses the interaction between apparent optical properties and inherent optical properties to estimate the concentration of suspended matter [14,15]. Despite its obvious physical significance, the construction process of the algorithm is complicated, and its application is usually challenging [16,17]. The empirical model often adopts the statistical regression method to construct a functional relationship between the concentration of suspended matter and the reflectance of a single band or a combination of bands [18,19], which is straightforward and easy to use, but it is relatively dependent on modeling data, and its applicability is greatly limited [20,21]. The semi-empirical/semi-analytical model not only combines the absorption and scattering process of light but also uses the empirical relationship to describe the relationship between C tsm and the reflectivity [22,23]. In contrast, it is a simplified form of the analytical model and has more theoretical basis than the empirical model, so it is widely used [24].
Lake Taihu, as a typical shallow inland eutrophication lake, has an average depth of less than 2 m [25], and its reflectance is a mixed reflectance composed of several components [26]. The current mainstream inversion models of C tsm do not consider the interference contribution of chlorophyll-a and CDOM deeply, and the directly established relationship model between water reflectance and C tsm is not enough to accurately reflect the information of C tsm . In fact, when C tsm in inland water is monitored by remote sensing, the common influence of other substances in the water body must be taken into account.
The Hydrolight radiative transfer model adopts the absorption and scattering characteristics of each component of the water body to simulate the radiative transfer process of light in the water body. It can also simulate the spectral contribution of a single component, which is useful for analyzing the complex interactions between different components of the water body [27]. Water reflectance can be regarded as a multivariate linear combination of several components in water [28], and the Hydrolight model can be used to determine the relationship between the concentration of a single component and its contribution to reflectance. Therefore, the extraction of C tsm can be viewed as a multivariate linear issue.
Based on the above analysis, this paper constructs an iterative inversion method of C tsm based on multiple linear regression and iterative analysis based on the spectrum setting of SDGSAT-1 MII and the Hydrolight radiation transfer simulation experiment in combination with the measured data in Lake Taihu. Multiple linear regression was used to establish the relationship between the water reflectance and the contribution of a single component, and these component contribution relations were replaced into a completely closed iterative inversion model for solving C tsm through joint analysis, to solve the problem of separating the reflectance contribution of other substances in the inversion process of C tsm . This study is expected to bring into play the application potential of the SDGSAT-1 satellite in water quality parameter inversion and provide an important technical reference for improving the water quality of Lake Taihu and controlling the water environment.

2. Materials and Methods

2.1. Study Area

Lake Taihu (30°56′~31°34′N, 119°54′~120°36′E, Figure 1), with a water area of 2338.1 km2, is the third largest freshwater lake in China [29], surrounded by many central cities of the Yangtze River Delta, such as Shanghai, Suzhou, Wuxi, Changzhou, Huzhou, and Jiaxing. It occupies an important position in the process of high-quality development and construction of the surrounding economy and the Yangtze River Delta.
With the economic development and urbanization of the surrounding basin, the water environment of Lake Taihu has been gradually damaged, and the problem of water pollution has become increasingly prominent. Among them, total suspended matter is an important index that intensifies the eutrophication of Lake Taihu and causes water pollution. Due to the shallow water characteristics of Lake Taihu, it is easy to be disturbed by wind and waves, resulting in the suspension of sediment and in perennial turbidity in the lake and significant suspended matter characteristics [30]. The increase in C tsm further induces the flow and release of other substances in the lake, promotes the proliferation of algae in the lake, intensifies the eutrophication process of the lake, and ultimately leads to the deterioration of water quality and water pollution. Therefore, monitoring C tsm in Lake Taihu is of great significance for improving the water environment and controlling water pollution.

2.2. Field Measurement Data

This study used two measured data sets: the in situ data set from Lake Taihu collected in August 2013 and October 2015 for model construction, and the data collection gathered synchronously with the satellite transit on 27 July 2022, for model verification. In Figure 1, the dots with different colors represented the distribution of the actual sampling points at the three sampling times. The following is a detailed description of the two-part data set:
(1) In August 2013 and October 2015, field data of a total of 54 valid samples were collected in Lake Taihu. The collection process mainly includes the field collection of water spectral data and water samples and the subsequent laboratory analysis of the collected water samples to obtain the concentration data of different components of the water body and the inherent optical data.
A spectrophotometer was used to detect chlorophyll-a concentration ( C chla ). C tsm was determined using the calcination weighing method. When measuring the spectral absorption of chlorophyll-a and suspended matter, the water sample was first filtered by a GF/F glass fiber filter with a diameter of 25 mm and then measured by an ultraviolet/visible dual-beam integrating sphere spectrophotometer [31]. For the measurement of the scattering characteristics of total suspended matter and chlorophyll-a, the original data of the measurement should be corrected and then calculated [32]. CDOM concentration ( C cdom ) is usually stated as its absorption coefficient at 440 nm, which was measured by an ultraviolet/visible light spectrophotometer and then converted into an absorption coefficient [33]. The CDOM absorption coefficient is calculated according to the following formula:
a c d o m ( λ ) = 2.303 D ( λ ) L c
where a cdom ( λ ) is the absorption coefficient of CDOM. D(λ) is the absorbance of the water sample at wavelength λ, and Lc is the optical range of the cuvette.
(2) To validate the application effect of the inversion model on satellite images, the measured C tsm data were determined using water samples collected in Lake Taihu on 27 July 2022 and consisted of 16 sample points. The field sampling date was 27 July 2022, from 7–11 a.m., and the transit of the SDGSAT-1 satellite over Lake Taihu was on 27 July 2022, at 8:51 a.m. Therefore, the synchronization between the field sampling data and the satellite transit time on 27 July 2022 was achieved [34]. The latitude and longitude information of the measured sample points was used to spatially match with the image pixels using geolocalization.
The purpose of water spectral measurement is to obtain water remote sensing reflectance R rs ( λ ) . The water surface spectral data are measured by ASD spectrometer through the above-water method [35]. In the case of avoiding direct solar reflection and ignoring external influences such as solar flares and white caps, the formula for calculating the above-water remote sensing reflectance is as follows:
R r s ( λ ) = ρ p ( λ ) ( L s w ( λ ) r s L s k y ( λ ) ) π L p ( λ )
where R rs ( λ ) is the reflectance of the water body; λ is the wavelength; and L sw ( λ ) , L sky ( λ ) , and L p ( λ ) are the signal values measured by the spectrometer facing water, sky, and reference plate, respectively. rs is the reflectance of the air–water boundary facing the skylight, which can be 0.022 for a calm water surface, and 0.025 for a wind speed of about 5 m/s, and can be considered between 0.026 and 0.028 for wind speeds of about 10 m/s. ρ p ( λ ) is the reflectivity of the reference plate.

2.3. Satellite Remote Sensing Data

This study obtained SDGSAT-1 MII images through the data open system provided by the International Research Center for Big Data for Sustainable Development Goals and selected SDGSAT-1 MII L4 images of Lake Taihu from February 2022 to November 2023 for remote sensing monitoring C tsm in Lake Taihu.
Before using remote sensing image data, a series of pretreatment work should be performed, including radiometric calibration, atmospheric correction, water extraction, and so on. The first is radiation calibration, which obtains the calibration gain values of each band of the image from the MII image meta-file. The specific parameters are shown in Table 1. The center wavelength of each channel is the wavelength corresponding to the maximum spectral response of the remote sensor. The formula for converting the observed value DN from satellite load to apparent radiance L is as follows:
L = D N × g a i n
where DN is the observed value of the satellite load, L is the apparent radiance, and gain is the gain value of the band.
As new remote sensing data, the SDGSAT-1 MII image is less researched at home and abroad, and there is no mature atmospheric correction algorithm specially applied to this image. FLAASH model is the most commonly used atmospheric correction model based on the MODTRAN radiative transfer algorithm [36], which applies to different types of sensors, and it can eliminate most of the atmospheric effects effectively. After that, the FLAASH module was used for atmospheric correction of the image to obtain the corrected reflectivity, and then, the output reflectivity image was divided by π to approximate the remote sensing reflectivity of the water [34,37]. Finally, normalized differential water body index (NDWI) was used to extract the mask from Lake Taihu [38].

2.4. Hydrolight Radiation Transfer Simulation

2.4.1. Model Method and Parameter Setting

The Hydrolight model is a tool developed using the radiation transmission theory of water bodies. This software can simulate the light transmission process in a wide range of water environments using the absorption and scattering properties of each component of water bodies. In this paper, the CASE 2 IOPS model of the Hydrolight model is used to simulate the remote sensing reflectance of Lake Taihu. The inherent optical parameters, water quality parameters, environmental conditions, and the spectral range of the absorption model and the scattering model under this model are set as follows:
  • Absorption model
The absorption coefficient a of the water body can be expanded as the sum of all absorption components [39], which is composed of four absorption contributions of pure water, chlorophyll-a, total suspended matter, and CDOM, and each absorption component can be expressed as the product of the inherent optical quantity per unit and the concentration of related components. Therefore, it can be expressed by the following formula:
a λ = a w λ + a c h l a * λ C c h l a + a c d o m * λ C c d o m + a t s m * λ C t s m
where a(λ) is the spectral total absorption. a w ( λ ) , a chla * ( λ ) , and a tsm * ( λ ) represent the absorption coefficient of pure water at wavelength λ, the specific absorption coefficient of chlorophyll-a, and the specific absorption coefficient of total suspended matter, respectively, and the values in the literature [40] are used for a w ( λ ) , and the mean values of the specific absorption coefficients of all sample points are used for a chla * ( λ ) and a tsm * ( λ ) , and the three-parameter curves are shown in Figure 2. C chla and C tsm that used the concentration data of 54 measured sample points were used. C cdom could be expressed as the absorption coefficient of CDOM at a wavelength of 440 nm [41], using 54 measured sample points as the input data of C cdom . a cdom * ( λ ) is the specific absorption coefficient of CDOM at wavelength λ, which conforms to the exponential attenuation model with the absorption coefficient of CDOM at a wavelength of 440 nm, and the specific relationship is as follows:
a c d o m * λ = a c d o m 440 exp S g λ 440
where a cdom ( 440 ) is the absorption coefficient of CDOM at a wavelength of 440 nm. S g is the exponential function slope of the CDOM absorption spectrum, and the calculated average S g is about 0.015. The mean data of a cdom * ( λ ) obtained for 54 sets of sample points are shown in Figure 2.
  • Backscattering model
The backscattering coefficient of the water body is composed of the backscattering contributed by pure water, chlorophyll-a, and total suspended matter [26], which can be expressed as follows
b b λ = b ˜ w b w λ + b ˜ c h l a b c h l a * λ C c h l a + b ˜ t s m b t s m * λ C t s m
where b b ( λ ) is the backscattering coefficients. b w ( λ ) , b chla * ( λ ) , and b tsm * ( λ ) are the scattering coefficients of pure water, specific scattering coefficients of chlorophyll-a, and specific scattering coefficients of total suspended matter at wavelength λ. The values of the model in the literature [42] are used for b w ( λ ) , and the mean values of the specific scattering coefficients of all sample points are used for b chla * ( λ ) and b tsm * ( λ ) . The three-parameter curves are shown in Figure 3. b ~ w , b ~ chla , and b ~ tsm are the backscattering proportions of pure water, chlorophyll-a, and total suspended matter, and the empirical values in reference [43] are set as 0.5, 0.005, and 0.028, respectively.
  • External environmental condition
The RADTRAN model was selected for the sky radiation transmission model, the default values of the weather parameters were adopted, and the refraction index of the water body was set to 1.34.
  • The range of the simulated spectrum
The band range of simulated remote sensing reflectance was 400 nm to 800 nm, and the simulation interval was 1 nm.

2.4.2. Influence of Single-Component Concentration Change on Reflectance

Based on the measured concentration ranges of three components, total suspended matter, chlorophyll-a, and CDOM, the Hydrolight model was used to simulate the influence of single-component concentration changes on remote sensing reflectance. In the simulation process, other input parameters remained unchanged except for the concentration of three components. First, a uniform minimum initial concentration was set. C tsm was 10 g/m3, C chla was 0.2 mg/m3, and a cdom ( 440 ) was 0.2 m−1. Subsequently, only the concentration of one component was changed, while the other two components remained unchanged. The change settings of C tsm , C chla , and a cdom ( 440 ) are shown in Table 2, and the remote sensing reflectance simulation data sets under different component concentrations were generated.
The relative change rate is an indicator used to measure the relative change degree of variables, and its calculation formula is usually the division of the change value by the previous value [44]. Therefore, in order to evaluate the influence degree of total suspended matter, chlorophyll-a, and CDOM on the reflectance under different concentrations, the index of relative change rate R rs Δ ( λ ) is introduced to quantify the relative change in the remote sensing reflectance of water body under different conditions. The specific calculation formula is shown as follows:
R r s Δ λ = R r s E λ R r s e λ R r s e λ 100 %
where R rs Δ ( λ ) represents the relative change rate, R rs E ( λ ) represents the remote sensing reflectance of the water body at wavelength λ after the concentration change, and R rs e ( λ ) represents the remote sensing reflectance of the water body at wavelength λ before the concentration change. The high change rate indicates that the remote sensing reflectance is easily affected by the concentration, while the low change rate indicates that the remote sensing reflectance is not easily affected by the concentration.
In order to have a clear understanding of the effect of single-component concentration changes on band reflectance, simulated R rs ( λ ) and R rs Δ ( λ ) were used to analyze the role of single-component concentration changes on the remote sensing reflectance of the center wavelengths of the sensor’s 2nd–6th bands.

2.5. Sensor Channel Convolution and Multi-Band Combination Enhancement

2.5.1. Sensor Channel Convolution

Since the energy actually received by each band of the satellite sensor is the sum of the energy received at each wavelength within the wavelength range of the band, it is necessary to use the spectral response function of the SDGSAT-1 MII sensor (Figure 4) for the convolution of the spectral reflectance simulated by Hydrolight and convert it into the remote sensing reflectance of the corresponding channel of the satellite. The calculation formula is as follows:
B i = λ 1 λ 2 R r s λ f s λ d λ λ 1 λ 2 f s λ d λ i = 2 , 3 , 4 , , 6
where B i is the equivalent remote sensing reflectance of i channel of the satellite; R rs ( λ ) is the reflectance of the simulated spectrum with an interval of 1 nm; and the spectral range of Hydrolight simulation basically covers the SDGSAT-1 MII channel 2–6. λ 1 and λ 2 are wavelengths at both ends of channel i of the satellite; f s ( λ ) is the spectral response function of the SDGSAT-1 MII sensor.

2.5.2. Multi-Band Combination Enhancement

Compared with the original single-band data, the combination of multi-band mathematical operations helps to eliminate the interference between various bands to a certain extent. Band ratio operation and band difference operation are both effective spectral processing methods, which can highlight the correlation between water quality parameters and water reflectivity [45].
Therefore, ratio reflectance combination and difference reflectance combination were obtained by processing the ratio of and difference in the convolved SDGSAT-1 MII channel reflectance data according to the following formula:
B i / j = B i B j ( i , j = 2 , 3 , 4 , , 6 ,   a n d   i j )
B i - j = B i B j ( i , j = 2 , 3 , 4 , , 6 ,   a n d   i j )
where B i / j is the ratio reflectance combination, and B i - j is the difference reflectance combination.

2.6. Correlation Analysis

Correlation analysis can quantitatively measure the intensity of linear relationship between variables [46], which provides a basis for the in-depth analysis of the relationship between the component concentration and reflectivity and for subsequent modeling and interpretation. Correlation analysis was conducted between the original band and multi-band combination and the measured C chla , C tsm , and a cdom ( 440 ) in this study.
In the correlation analysis, in addition to the correlation coefficient, a statistic is also obtained, that is, p. The t-test is usually used to determine whether the correlation between the variables is significant, and the level of significance is used to obtain the statistic p to assess whether the correlation between the variables is significant or not. The correlation coefficient, statistic t, and p are calculated by the following formula:
r = n = 1 N x n x ¯ y n y ¯ n = 1 N x n x ¯ 2 n = 1 N y n y ¯ 2
t = r N 2 1 r 2
p = 2 1 t t x d x
where r is the correlation coefficient, and n is the sample sequence number; N is the total sample size; x represents different original bands and combinations of bands; and y represents C chla , C tsm , and a cdom ( 440 ) . When p > 0.05, there was no significant correlation between the variables. When p < 0.01, it means that the correlation between variables at 0.01 level is significant, and the significance is indicated by **; when p < 0.05, it means that the correlation between variables at the level of 0.05 is significant, and the significance is indicated by *.

2.7. Iterative Inversion Method of Ctsm

2.7.1. Iterative Inversion Method Based on Multiple Linear Regression Analysis

Remote sensing reflectance of the water body is the mixed reflectance of multiple components of the water body, and the mixed reflectance R can be regarded as the linear sum of the contribution of C tsm and the contribution of other components concentration, as shown in the following formula:
R = ( R ) t s m + ( R ) Other   components
where R can be the reflectance in the form of different bands and different combinations of bands, viewed as consisting of multiple single-component contributions to reflectance. ( R ) tsm represents the single-component remote sensing reflectance contributed by C tsm , and ( R ) other   components represents the component remote sensing reflectance contributed by the concentration of other components except C tsm .
In order to reduce parameter dependence, the contribution of a single component is regarded as a function of its concentration. For example, the relationship between the single-component reflectance contributed by C tsm and C tsm is given by the following equation:
( R ) t s m = h C t s m
where h ( C tsm ) represents the functional relationship between the single-component remote sensing reflectance contributed by C tsm and C tsm .
The Hydrolight model can simulate the remote sensing reflectance of a single component using a single-component concentration, and the relationship between the remote sensing reflectance of a single component and the concentration of a single component can be obtained. By establishing the multiple linear regression relationship between the mixed remote sensing reflectance and the remote sensing reflectance contributed by a single component, C tsm is taken as the unknown in the multiple regression relationship to solve C tsm by using different bands or combinations of bands. Therefore, in theory, it is feasible to invert C tsm based on multiple linear regression.
The purpose of this study is to reduce the influence of other component concentrations on the inversion of C tsm . The most critical thing is to select the appropriate band or band combination as the characteristic band or the characteristic band combination to establish the multiple linear regression relationship between the characteristic band or the characteristic band combination and the single-component remote sensing reflectance, according to the analysis results of the influence of single-component concentration change on the reflectance and the correlation analysis results. This step is a prerequisite for retrieving C tsm .
On the premise that the characteristic band or combination of bands is determined as the mixed reflectance, based on the relationship between the concentration of one component and the remote sensing reflectance of one component, and the multiple linear regression relationship between the mixed reflectance and the reflectance of one component, a closed relationship can be formed with C tsm as an iterative variable by solving C tsm in order to unite the multiple regression model, and an iterative functional relationship can be formed, as shown in the following equation:
C t s m m + 1 = φ C t s m m m = 0 ,   1 ,   2 ,   ,   M
where C tsm ( m ) represents the input value of the model, C tsm ( m + 1 ) represents the output value of the model, and m represents the calculation times of the model to solve C tsm , that is, the number of iterations.
By iterating the function relation, we can derive the next value from the previous value one-by-one and control the cycle process in the form of an infinite approximation solution until the circulation is stopped, when C tsm remains unchanged. The final output result C tsm ( m + 1 ) represents the inversion result of C tsm .

2.7.2. Evaluation Indicators

Three indexes, namely, determination coefficient ( R 2 ) , root mean square error (RMSE), and mean absolute percentage error (MAPE), are used to evaluate the effect of the inversion model. Their specific calculation formulas are as follows:
R 2 = 1 h = 1 H y h y h 2 h = 1 H y ¯ y h 2
RMSE = 1 H h = 1 H y h y h 2
MAPE = 100 % H h = 1 H | y h y h y h |
where y h represents the retrieved value of C tsm at the hth sample point; y h represents the measured value of C tsm . at the hth sample point. H is the number of sample points.

3. Results

3.1. Validation of Remote Sensing Reflectance Simulations

According to the method described in Section 2.4.1, 54 sets of simulated R rs ( λ ) of Lake Taihu in the wavelength range of 400–800 nm were obtained and compared with the measured R rs ( λ ) . The comparison results are shown in Figure A1. The simulated R rs ( λ ) and the measured R rs ( λ ) presented relatively similar curve characteristics, and the spectral trend was relatively consistent. Figure 5a shows the fitting results of the simulated R rs ( λ ) to the measured R rs ( λ ) in the range of 400–800 nm, and there is an obvious linear relationship between the simulated R rs ( λ ) and the measured R rs ( λ ) (R2 = 0.906), and the simulation error is small (RMSE = 0.0031 sr−1). Figure 5b shows the fitting results of the simulated R rs ( λ ) to the measured R rs ( λ ) at the center wavelengths of the satellite bands 2–6, and these scatters are uniformly distributed along the 1:1 line. Error between simulated R rs ( λ ) and measured R rs ( λ ) at the center wavelength of the satellite band are shown in Table 3, with 0.826–0.906 for R2, and RMSE in the range of 0.0025–0.0038 sr−1.
Figure 6 shows the measured data of C chla , C tsm , and a cdom ( 440 ) of 54 sample points, as well as the correlation coefficients between simulated R rs ( λ ) and measured R rs ( λ ) of the 54 sample points. Measured C tsm ranged from 15–145 g/m3, measured C chla ranged from 0.27–133 mg/m3, and measured a cdom ( 440 ) ranged from 0.22–1.33 m−1. The statistical results showed that the correlation coefficients between simulated R rs ( λ ) and measured R rs ( λ ) of the 54 groups of sample points were above 0.93, and the proportion of samples with correlation coefficients greater than 0.95 reached 83.3%.
The above results show that the simulated R rs ( λ ) could be used as reliable basic data for correlation theoretical analysis.

3.2. Analysis of Influence Characteristics of Remote Sensing Reflectance

Based on the measured single-component concentration range, an experiment on the effect of single-component concentration changes on R rs ( λ ) was designed, and the simulated R rs ( λ ) data sets were obtained under the conditions that C tsm varied from 10 to 150 g/m3, C chla varied from 0.2 to 140 mg/m3, and a cdom ( 440 ) varied from 0.2 to 1.3 m−1. Using these simulated R rs ( λ ) data sets, the changes in R rs ( λ ) and R rs Δ ( λ ) at different concentrations were analyzed, as well as the effect of a single-component concentration on the center wavelength of the sensor’s bands 2–6, and the results are shown in Figure 7 and Figure 8.
In the range of 400–800 nm, R rs ( λ ) and R rs Δ ( λ ) are shown in Figure 7a,d under different C tsm . As C tsm increased, R rs ( λ ) showed a gradual increase trend, and R rs Δ ( λ ) showed a gradual decrease trend. When C tsm was low, R rs Δ ( λ ) was more obvious, and the maximum could cause a 90% relative change. With the increase in concentration, R rs Δ ( λ ) gradually decreased and finally tended to be stable. R rs ( λ ) and R rs Δ ( λ ) are shown in Figure 7b,e under different C chla . Before 700 nm, the increase in C chla led to a decrease in R rs ( λ ) , and R rs Δ ( λ ) decreased. After 700 nm, with the increase in C chla , R rs ( λ ) showed an upward trend, but R rs Δ ( λ ) gradually decreased, resulting in a positive change degree of R rs ( λ ) up to 12%. R rs ( λ ) and R rs Δ ( λ ) under different a cdom ( 440 ) are shown in Figure 7c,f. The influence of a cdom ( 440 ) on R rs ( λ ) was weak on the whole. With the increasing wavelength, R rs ( λ ) gradually decreased, and R rs Δ ( λ ) gradually decreased. After 700 nm, R rs ( λ ) is basically no longer affected by the change in a cdom ( 440 ) .
On the center wavelength of the sensor’s bands 2–6, R rs ( λ ) and R rs Δ ( λ ) at the center wavelength of the five bands are shown in Figure 8a,d under different C tsm . With the increase in C tsm , R rs ( λ ) at the center wavelength of the five bands showed an increasing trend at the same time, and R rs ( 773 ) showed a higher degree of change. R rs ( λ ) and R rs Δ ( λ ) at the center wavelength of the five bands are shown in Figure 8b,e under different C chla . As C chla increased, the trends and extent of changes in R rs ( λ ) at the center wavelengths of the four bands except 773 nm were more consistent. Under different a cdom ( 440 ) , R rs ( λ ) and R rs Δ ( λ ) at the center wavelength of the five bands are shown in Figure 8c,f. R rs ( 773 ) and R rs Δ ( 773 ) remained essentially unchanged with an increase in a cdom ( 440 ) .

3.3. Correlation between Remote Sensing Reflectance and Concentration of Three Components

The results of the correlation analysis between B i , B i / j , and B i - j and C tsm , C chla , and C cdom , respectively, are shown in Figure A2, from which the following can be inferred:
(1) The results of correlation analysis between B i and the different concentrations are shown in Figure A2a. C tsm showed a positive correlation with all five bands, with the strongest significant positive correlation with B 6 (r = 0.985 (p < 0.01)); C chla had the strongest correlation with B 3 (r = −0.465 (p < 0.01)); C cdom showed a negative correlation with all five bands, and with B 2 with the strongest correlation (r = −0.649 (p < 0.01)).
(2) In the correlation analysis of B i / j with the different concentrations (Figure A2b), 55% of B i / j had a significant correlation with C tsm above 0.7; 15% of B i / j had a significant correlation with C chla above 0.7; and C cdom had a correlation above 0.7 with just two ratio combinations of bands.
(3) The results of the correlation analysis between B i - j and concentrations are shown in Figure A2c. The number of B i - j with significant correlation with C tsm higher than 0.7 accounted for 50%; the number of B i - j with significant correlation with C chla higher than 0.7 accounted for 40%; and the correlation between C cdom and B i - j was generally low, among which the strongest correlation was with B 6 - 4 (r = −0.51 (p < 0.01)).
(4) A comparative analysis of the correlations between concentrations and different bands and combinations of bands showed that three bands ( B 2 , B 3 , B 4 ), as well as six combinations of ratios ( B 2 / 3 , B 2 / 6 , B 3 / 2 , B 3 / 6 , B 4 / 6 , B 5 / 6 ) and eight combinations of differences ( B 2 - 4 , B 2 - 6 , B 3 - 6 , B 4 - 2 , B 4 - 6 , B 6 - 2 , B 6 - 3 , B 6 - 4 ), were significantly correlated with the three component concentrations at the same time, indicating that it is inappropriate to rely on the bands and combinations of bands related to C tsm alone for the inversion of C tsm , and it is not possible to avoid the effects of the other components.
(5) With a further comprehensive comparison of the correlation characteristics of the three component concentrations, as well as by combining the results of the characterization of the effects of the three component concentrations on reflectance in Section 3.2, it can be concluded that C cdom has a weak effect on reflectance, while there is a more obvious correlation between C chla and C tsm and remote sensing reflectance, and the two are the dominant factors affecting remote sensing reflectance.
Based on the above conclusions, in order to eliminate the other component contributions from the remote sensing reflectance and then invert C tsm independently, the basic idea of constructing C tsm model in this study includes the following:
(1) Firstly, the band or combination of bands that did not correlate with C cdom but showed a correlation with C tsm and C chla were selected, and the results are shown in Figure 9, which were regarded as the combined contribution of the two components of C tsm and C chla , so that the influence of C cdom was no longer taken into account.
(2) Next, the bands or band combinations that showed a significantly strong correlation (r > 0.70 (p < 0.01)) with both C chla and C tsm were selected as factors, and multiple linear regression and iterative resolution methods were used to remove the signals of C chla from the remote sensing reflectance, and thus extract C tsm .
Finally, B 6 / 3 and B 6 / 5 were identified as two band combination characterization factors for subsequent multiple linear regression and iterative parsing studies.

3.4. Multivariate Iterative Inversion Model Construction and Evaluation

3.4.1. Multivariate Iterative Inversion Model Construction

  • Reflectance simulation of single-component contribution
B 6 / 3 and B 6 / 5 were regarded as the mixed reflectance of both C tsm and C chla . In order to accurately invert C tsm , it is necessary to simulate the contribution of C tsm and C chla , respectively, followed by contribution elimination through multivariate iteration.
Based on the measured concentration data, Hydrolight was used to simulate the single-component contributions to obtain the single-component reflectance of the contribution from C tsm and the single-component reflectance of the contribution from C chla (Figure 10 and Figure 11), respectively, and the relationship that exists between the component concentrations and the respective contributions to the reflectance is shown in the following equations:
( B 6 / 3 ) t s m = 0.0066 C t s m + 0.0207   R 2 = 0.996
( B 6 / 5 ) t s m = 0.0028 C t s m + 0.1391   R 2 = 0.995
( B 6 / 3 ) c h l a = 0.0041 C c h l a + 0.0065   R 2 = 0.999
( B 6 / 5 ) c h l a = 0.0054 C c h l a + 0.1552   R 2 = 0.997
where ( B 6 / 3 ) tsm and ( B 6 / 5 ) tsm are single-component reflectances contributed by C tsm , and ( B 6 / 3 ) chla and ( B 6 / 5 ) chla are single-component reflectances contributed by C chla .
The results showed that there was a significant linear relationship between the component concentrations and their respective contributions to remote sensing reflectance, with R2 being greater than 0.99 in all cases.
  • Multiple linear regression analysis
According to Equation (14) in Section 2.7.1, B 6 / 3 and B 6 / 5 are equivalent to the mixed reflectance R, B 6 / 3 can be regarded as the linear sum of the reflectance contributed by each of C tsm and C chla , i.e., a linear combination of ( B 6 / 3 ) tsm and ( B 6 / 3 ) chla , and B 6 / 5 can be thought of as the linear sum of the reflectance contributed by each of C tsm and C chla , i.e., a linear combination of ( B 6 / 5 ) tsm and ( B 6 / 5 ) chla . Multiple linear regression analysis between the mixed reflectance R and the single-component reflectance using Equation (14). The objective of this study is to solve C tsm by combining Equation (15) from Section 2.7.1 with a multiple linear regression analysis, and therefore, a single-component-contributing reflectance is used as the dependent variable in performing the multiple linear regression, and a characteristic band combination and another single component contributing reflectance are used as independent variables.
Considering the correlation between B 6 / 3 and C tsm (r = 0.942, p < 0.01, Figure 9) was higher than that with C chla (r = 0.709, p < 0.01, Figure 9), a multiple regression was established with ( B 6 / 3 ) tsm as the dependent variable, and B 6 / 3 and ( B 6 / 3 ) chla as the two independent variables. Meanwhile, B 6 / 5 correlation with C chla (r = 0.912, p < 0.01, Figure 9) was greater than the correlation with C tsm (r = 0.779, p < 0.01, Figure 9), and therefore, ( B 6 / 5 ) chla was used as the dependent variable and a multiple regression relationship was established with B 6 / 5 and ( B 6 / 5 ) tsm as the two independent variables.
B 6 / 3 and B 6 / 5 were obtained by ratio combination based on satellite channel data after the convolution of 54 sets of simulated spectra, ( B 6 / 3 ) tsm and ( B 6 / 5 ) tsm simulated from C tsm of the 54 sample sites (Figure 10), and ( B 6 / 3 ) chla and ( B 6 / 5 ) chla simulated from C chla of the 54 sample sites (Figure 11), which correspond to the decomposition of the component contributions to B 6 / 3 and B 6 / 5 .
Based on the above analysis and on the principle of Ordinary Least Squares (OLSs), the multivariate linear regression equations were established by Python, and the two regression equations obtained were expressed as follows.
( B 6 / 3 ) t s m = 1.07305 B 6 / 3 0.9504 ( B 6 / 3 ) c h l a 0.06868 R 2 = 0.977
B 6 / 5 c h l a = 1.05341 B 6 / 5 0.89225 ( B 6 / 5 ) t s m + 0.09336 R 2 = 0.991
The corresponding R2 of the two regression equations were 0.977 and 0.991, indicating that the models were well fitted. The residual is the difference between the predicted value of the dependent variable and the cause variable obtained after multiple linear regression. The residuals of the two models were located between −0.045–0.059 sr−1 and −0.015–0.034 sr−1, respectively, with smaller residual values, which indicated that the multiple linear regression model was able to predict the dependent variable more accurately.
  • Construction of the iterative inversion model
Based on the multiple linear regression analysis, the two multiple regression models (Equations (24) and (25)) were combined with the relationship between the concentration and the respective contributing reflectance, thus forming a complete closed iterative inversion model, which is shown in the following equation:
C t s m ( m + 1 ) = φ C t s m m = 162.58333 ( B 6 / 3 ) ( 115.17283 ( B 6 / 5 ) 0.27315 C t s m ( m ) 5.85233 ) m = 0 , 1 , 2 , , M
The input parameter settings of the iterative inversion model in different application scenarios are shown in Table 4. In the iterative inversion model, B 6 / 3 , B 6 / 5 , and C tsm ( 0 ) were used as the input parameters of the model, and this iterative functional relationship was convergent and able to calculate C tsm iteratively. When the successive output values gradually converged to be equal, i.e., the difference between the output value and the input value gradually tended to be zero, which indicated that the algorithm was approaching the solution until the output value was kept unchanged, the iterative process was stopped, and the last guessed value C tsm ( m + 1 ) was used as the final inversion result of C tsm .

3.4.2. Analysis of the Initial Value Setting of Ctsm(0)

Firstly, the model was applied to the modeling data, and the calculation method was adopted by numerical iteration one-by-one. The input data of the model included B 6 / 3 , B 6 / 5 , and C tsm ( 0 ) . And C tsm ( m + 1 ) was initialized to 0 as the output value. In the iterative process, both C tsm ( m ) and C tsm ( m + 1 ) were updated automatically according to the iterative formula (Equation (26)).
In order to analyze the impact of the initial value C tsm ( 0 ) on the model, the C tsm ( 0 ) of each group was set as 1 g/m3, 10 g/m3, 50 g/m3, 100 g/m3, 200 g/m3, and 500 g/m3, respectively. Figure 12 shows the change in the difference between the 54 sets of output values and input values during iteration, using 54 groups of curves of different colors. Under different values of C tsm ( 0 ) , the difference between the output value C tsm ( m + 1 ) and the input value C tsm ( m ) decreased gradually with the iteration. When the output value and the input value were consistent, the iteration process stopped, and each group of values was iterated 29 times. The model achieved the same output results under different C tsm ( 0 ) settings, and the number and time of calculation did not change with the specified initial value of C tsm ( 0 ) . The results showed that the iterative inversion model had no dependence on the initial value of C tsm ( 0 ) , and the initial value of C tsm ( 0 ) could be set arbitrarily.
In addition, the relative error (RE) was used to quantify the convergence accuracy in the iterative process, which was defined as the difference between the output value of the current iteration and the final output value divided by the final output value, and was expressed as a percentage. When C tsm ( 0 ) = 1 g/m3, RE changes in the iteration process are shown in Figure 13. With the progress of iteration, RE decreased gradually. When the number of iterations was five, RE reached −0.15%. When the number of iterations was 10, RE decreased to −0.0002%. When the number of iterations reached 29, RE dropped to 0, indicating that the model had successfully converged to the final solution and reached a stable output state.
By analyzing the measured C tsm and the retrieved C tsm ( m + 1 ) of the iterative inversion model, the accuracy of the iterative inversion model was evaluated (Figure 14). There was a highly significant linear relationship between the measured value of C tsm and the inversion value, R2 was 0.97, and the difference between the two was small, and the obtained RMSE and MAPE were 4.89 g/m3 and 11.48%, respectively, indicating that the model had good accuracy.

3.4.3. Valuation of Ctsm Estimation Methods for SDGSAT-1 MII Image

The iterative inversion model constructed in this paper was applied to the pre-processed SDGSAT-1 MII image of 27 July 2022, and verified by using the ground and satellite synchronization data by adopting the pixel-by-pixel iterative calculation method. The comparison results between the measured C tsm and the retrieved C tsm are shown in Figure 15. They were uniformly distributed along the 1:1 line, and their R2 was 0.87, showing extremely high linear correlation and consistency, and the error level was relatively low, with a RMSE of 3.92 g/m3 and a MAPE of 8.13%. The verification results showed that the iterative inversion model proposed in this paper showed a good application effect and good application accuracy on the SDGSAT-1 MII image.

4. Discussion

4.1. Comparison between Different Characteristic Band Combinations

The aim of this paper is to eliminate the other component contributions from the reflectance and thus invert the suspension concentration independently, and based on the results of the analyses in Section 3.2 and Section 3.3, it is concluded that B 6 , B 6 / 2 , B 6 / 3 , B 6 / 4 , and B 6 / 5 do not correlate with C cdom but rather correlate with C tsm and C chla . Thus, it is possible to consider B 6 , B 6 / 2 , B 6 / 3 , B 6 / 4 , and B 6 / 5 as C tsm and C chla as the combined contribution of the two components (Figure 9); thus, the effect of C cdom is no longer considered. Furthermore, to remove the signal of C chla present in reflectance, B 6 / 3 and B 6 / 5 , which showed significant strong correlation (r > 0.70 (p < 0.01)) with both C tsm and C chla , were selected as factors for multiple linear regression and iterative resolution.
In order to evaluate the advantage of constructing an iterative inversion model for the inversion of C tsm using B 6 / 3 and B 6 / 5 as characteristic band combinations, from the band and band combinations that do not correlate with C cdom but correlate with C tsm and C chla (Figure 9), other band and band combinations were selected as new characteristic band combinations to construct new iterative inversion models for comparative studies.
Considering that only B 6 / 5 had a greater correlation with C chla than C tsm , B 6 / 5 was left unchanged. Among the bands and band combinations with higher correlation with suspended matter concentration than C chla , the correlations of B 6 , B 6 / 2 , B 6 / 3 , and B 6 / 4 with C tsm were all greater than 0.9, but there was a difference in the correlation with C chla , where the correlation of B 6 / 3 was greater than 0.7 and the correlation of B 6 with C chla was only 0.471, and the correlations of B 6 / 2 , B 6 / 4 , and C chla were greater than 0.65, and B 6 , B 6 / 2 and B 6 / 4 were used as the characteristic band combinations for constructing the iterative inversion model together with B 6 / 5 , respectively.
Three iterative inversion models based on B 6 and B 6 / 5 , B 6 / 2 and B 6 / 5 , and B 6 / 4 and B 6 / 5 were formed by repeating the previous processing process. The iterative inversion of C tsm applied to the modeling data are shown in Figure 16 by the three newly constructed iterative inversion models, the 54 sets of curves with different colors represent the changes in the difference between the 54 sets of output values and input values during the iteration process. When models were applied to the modeled data, C tsm ( 0 ) was uniformly set as 1 g/m3, and the three iterative inversion models carried out 29 iterations for each set of values. Figure 17 shows the fitting of the three iterative inversion models between the final output iterative inversion values and the measured values of each group, and the results are relatively close. Subsequently, the three newly constructed iterative inversion models were applied to the pre-processed SDGSAT-1 MII image of 27 July 2022, and verified using the ground and satellite synchronization data. The results are shown in Figure 18, where a few points deviate from the 1:1 line to a large extent.
By comparing the modeling errors of the four iterative inversion models and the errors applied to the SDGSAT-1 MII image (Table 5), the comparison results showed that the modeling accuracy and application effect of the three new iterative inversion models were not as good as the iterative inversion models based on B 6 / 3 and B 6 / 5 , which indicated the iterative inversion model based on B 6 / 3 and B 6 / 5 showed higher accuracy and a better ability of retrieving C tsm .

4.2. Comparison with Other Inversion Models

Referring to the inversion models of C tsm proposed by other researchers in recent years, the measured data used in this study were applied to the corresponding models to obtain the adapted models (Table 6), which were applied for the estimation of C tsm in Lake Taihu, and the results of the fit between measured C tsm and retrieved C tsm are shown in Figure 19. Among these adapted models, the Hou_2017 model [18] had the highest modeled accuracy in Lake Taihu, but its error was still large compared with the iterative inversion model proposed in this paper when models were applied to the modeled data.
These models were applied to the pre-processed SDGSAT-1 MII image of 27 July 2022, using B 4 as R rs ( 555 ) , and B 5 as R rs ( 645 ) and R rs ( 660 ) , and validating them using the synchronized ground-satellite data (Figure 20). The RMSE and MAPE between measured C tsm and inverted C tsm are listed in Table 7. The results showed that the Shi_2015 model [47] had the smallest error in the imagery among the adapted models. However, all of these adapted models performed poorly in terms of imagery applicability when compared to the iterative inversion model proposed in this paper.

4.3. Temporal and Spatial Distribution Characteristics of Ctsm in Lake Taihu

The iterative inversion model based on B 6 / 3 and B 6 / 5 is applied to the SDGSAT-1 MII images available in 2022–2023 for the observation of C tsm in Lake Taihu, and the spatial distribution results are shown in Figure 21. There were some seasonal differences in the spatial distribution of C tsm in Lake Taihu, with low concentrations mainly concentrated in the center of the lake. In May, June, July, and November, C tsm in the center of the lake was low and lower than that in the surrounding waters. The reason for this was that the vegetation growing in the water could effectively prevent suspended sediment from being resuspended, thus slowing down the increase in C tsm ; in January and February, the aquatic vegetation was in the sinking stage, coupled with the frequent monsoon winds, and the suspended sediment was resuspended again under the action of the wind; thus, C tsm in this period showed a phenomenon of accumulation in the central area of the lake.
In order to analyze the temporal changes in C tsm in Lake Taihu, the inversion results of C tsm were divided into six concentration intervals, including 0–15 g/m3, 15–30 g/m3, 30–45 g/m3, 45–60 g/m3, 60–75 g/m3, and > 75 g/m3, and the proportion of image pixels of C tsm of each interval in each quarter was counted as a percentage of the overall effective image pixels. The effective pixels did not include the algal bloom pixels and the pixels covered by clouds and fog. Figure 22 shows the results of the proportion of C tsm pixels in each concentration interval in different seasons in Lake Taihu. In general, C tsm was higher in spring, summer, and winter, while it was lower in fall. Due to the small amount of available SDGSAT-1 MII images in 2022 and 2023, detailed seasonal changes will need to be determined through the long-term monitoring of data.

5. Conclusions

Considering that the composition and optical properties of inland water are affected by the concentration of different components in water, in order to overcome the interference of C chla and C cdom in the inversion study of C tsm , this study established an iterative inversion model for estimating C tsm based on multiple linear regression based on the characteristics of SDGSAT-1 MII sensor. The Hydrolight model was used to simulate the radiative transfer process of Lake Taihu, and the influence characteristics and correlation analysis of C tsm , C chla , and C cdom on remote sensing reflectance were carried out. On this basis, the characteristic band combinations for multiple linear regression analysis were determined. The iterative inversion model of solving C tsm is formed by combining two multivariate linear regression models. The iterative inversion model was validated using modeling data and image data, and it was compared with the iterative inversion model built with other band combinations as characteristic band combinations, as well as with other published models. Furthermore, the spatiotemporal distribution characteristics of C tsm in Lake Taihu were investigated using SDGSAT-1 MII images taken in 2022–2023. The primary conclusions of this study are as follows:
(1) The simulation of the Hydrolight model produced 54 groups of 400–800 nm simulated remote sensing reflectance of Lake Taihu, and validation revealed that the simulated spectra had a high similarity with the shape of the measured spectral curves, and the correlation between the simulated spectra and the measured spectra of each group was over 0.93, and the coefficient of determination of the two fitted together was R2 of 0.906, which provided the necessary data for the accurate inverse inversion of C tsm .
(2) Through the characterization of the influence of remote sensing reflectance and the analysis of the correlation between different bands and band combinations and the concentrations of the three components of total suspended matter, chlorophyll-a, and CDOM, the band combinations that did not correlate with C cdom and showed a significant correlation with C tsm and C chla (r > 0.7 (p < 0.01)), B 6 / 3 , and B 6 / 5 were selected as the characteristic band combinations, which provided a theoretical basis for the role of removing chlorophyll-a and CDOM in the inverse study of C tsm .
(3) The two multiple linear regression models established with B 6 / 3 and B 6 / 5 as the combination of the characteristic bands fit well, with R2 of 0.977 and 0.991, respectively; combining the relationship between the concentration of a single component and its contribution to the remote sensing reflectance, the two multiple regression models were joined to form a complete closed iterative inversion model. And applying the iterative inversion model to the modeling data and the image data, the R2 was 0.97 and 0.87, the RMSE was 4.89 g/m3 and 3.92 g/m3, and the MAPE was 11.48% and 8.13%, respectively, with strong consistency between the retrieved C tsm and the measured C tsm .
(4) Compared with the iterative inversion model based on B 6 and B 6 / 5 , the iterative inversion model based on B 6 / 2 and B 6 / 5 , the iterative inversion model based on B 6 / 4 and B 6 / 5 , and other published models, the iterative inversion model established in this paper based on B 6 / 3 and B 6 / 5 had higher accuracy and was more suitable for the SDGSAT-1 MII images of Lake Taihu, and C tsm of Lake Taihu in 2022–2023 remote sensing monitoring results showed that there were differences in the spatial distribution of C tsm , with low concentrations mainly concentrated in the center region of the lake, and lower C tsm in autumn than in the other three seasons.
Overall, the iterative inversion model developed for SDGSAT-1 MII in this study had good applicability in the Lake Taihu area. Due to the limited number of measured C tsm synchronized with the time of SDGSAT-1 satellite transit over Lake Taihu, the collection time of field sample points failed to cover four seasons, and it is a challenge to fully validate the model using multiple images. In future research, we endeavor to collect more measured C tsm synchronized with the SDGSAT-1 satellite transit time of Lake Taihu in different seasons to explore and validate the applicability of the inversion model on SDGSAT-1 satellite more comprehensively, and to provide a theoretical basis and technical support for remote sensing monitoring and improvement of the water environment in Lake Taihu and other inland water.

Author Contributions

Conceptualization, X.H. and J.L.; methodology, X.H., J.L. and Y.S. (Yuan Sun); validation, Y.S. (Yuan Sun), Y.B. and Y.S. (Yonghua Sun); formal analysis, X.H. and J.L.; resources, J.L.; writing—original draft preparation, X.H.; writing—review and editing, J.L., Y.S. (Yuan Sun), Y.B., Y.S. (Yonghua Sun) and X.C.; visualization, X.H.; supervision, J.L., X.C. and Y.Y.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Civil Aerospace Technology Advance Research Project of China and the National Natural Science Foundation of China (Grant No. 41971391).

Data Availability Statement

The original contributions presented in the study are included in the article material, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Figure A1. Comparison between simulated R rs ( λ ) and measured R rs ( λ ) . (a1a8), (b1b8), (c1c8), (d1d8), (e1e8), (f1f8), and (g1g6) are the comparisons of 54 groups of R rs ( λ ) .
Figure A1. Comparison between simulated R rs ( λ ) and measured R rs ( λ ) . (a1a8), (b1b8), (c1c8), (d1d8), (e1e8), (f1f8), and (g1g6) are the comparisons of 54 groups of R rs ( λ ) .
Remotesensing 16 01385 g0a1
Figure A2. Correlation analysis results: (a) correlation between B i and the concentrations of the three components; (b) correlation between B i / j and the concentrations of the three components; and (c) correlation between B i   - j and the concentrations of the three components.
Figure A2. Correlation analysis results: (a) correlation between B i and the concentrations of the three components; (b) correlation between B i / j and the concentrations of the three components; and (c) correlation between B i   - j and the concentrations of the three components.
Remotesensing 16 01385 g0a2

References

  1. Sun, M.; Zhang, L.; Yang, R.; Li, X.; Zhao, J.; Liu, Q. Water resource dynamics and protection strategies for inland lakes: A case study of Hongjiannao Lake. J. Environ. Manag. 2024, 355, 120462. [Google Scholar] [CrossRef]
  2. He, Y.; Lu, Z.; Wang, W.; Zhang, D.; Zhang, Y.; Qin, B.; Shi, K.; Yang, X. Water clarity mapping of global lakes using a novel hybrid deep-learning-based recurrent model with Landsat OLI images. Water Res. 2022, 215, 118241. [Google Scholar] [CrossRef] [PubMed]
  3. Soomets, T.; Uudeberg, K.; Jakovels, D.; Brauns, A.; Zagars, M.; Kutser, T. Validation and comparison of water quality products in baltic lakes using sentinel-2 msi and sentinel-3 OLCI data. Sensors 2020, 20, 742. [Google Scholar] [CrossRef]
  4. Zhai, M.; Zhou, X.; Tao, Z.; Lv, T.; Zhang, H.; Li, R.; Huang, Y. Retrieve of total suspended matter in typical lakes in China based on broad bandwidth satellite data: Random forest model with Forel-Ule Index. Front. Environ. Sci. 2023, 11, 1132346. [Google Scholar] [CrossRef]
  5. Li, N.; Zhang, Y.; Shi, K.; Zhang, Y.; Sun, X.; Wang, W.; Huang, X. Monitoring water transparency, total suspended matter and the beam attenuation coefficient in inland water using innovative ground-based proximal sensing technology. J. Environ. Manag. 2022, 306, 114477. [Google Scholar] [CrossRef] [PubMed]
  6. Tian, L.; Wai, O.W.; Chen, X.; Li, W.; Li, J.; Li, W.; Zhang, H. Retrieval of total suspended matter concentration from Gaofen-1 Wide Field Imager (WFI) multispectral imagery with the assistance of Terra MODIS in turbid water–case in Deep Bay. Int. J. Remote Sens. 2016, 37, 3400–3413. [Google Scholar] [CrossRef]
  7. Xiong, Y.; Ran, Y.; Zhao, S.; Zhao, H.; Tian, Q. Remotely assessing and monitoring coastal and inland water quality in China: Progress, challenges and outlook. Crit. Rev. Environ. Sci. Technol. 2020, 50, 1266–1302. [Google Scholar] [CrossRef]
  8. Ge, K.; Liu, J.; Wang, F.; Chen, B.; Hu, Y. A Cloud Detection Method Based on Spectral and Gradient Features for SDGSAT-1 Multispectral Images. Remote Sens. 2022, 15, 24. [Google Scholar] [CrossRef]
  9. Shen, M.; Luo, J.; Cao, Z.; Xue, K.; Qi, T.; Ma, J.; Liu, D.; Song, K.; Feng, L.; Duan, H. Random forest: An optimal chlorophyll-a algorithm for optically complex inland water suffering atmospheric correction uncertainties. J. Hydrol. 2022, 615, 128685. [Google Scholar] [CrossRef]
  10. Silveira Kupssinskü, L.; Thomassim Guimarães, T.; Menezes de Souza, E.; Zanotta, D.C.; Roberto Veronez, M.; Gonzaga, L., Jr.; Mauad, F.F. A method for chlorophyll-a and suspended solids prediction through remote sensing and machine learning. Sensors 2020, 20, 2125. [Google Scholar] [CrossRef]
  11. Liu, J.; Liu, J.; He, X.; Pan, D.; Bai, Y.; Zhu, F.; Chen, T.; Wang, Y. Diurnal dynamics and seasonal variations of total suspended particulate matter in highly turbid hangzhou bay waters based on the geostationary ocean color imager. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2018, 11, 2170–2180. [Google Scholar] [CrossRef]
  12. Wen, Z.; Wang, Q.; Liu, G.; Jacinthe, P.-A.; Wang, X.; Lyu, L.; Tao, H.; Ma, Y.; Duan, H.; Shang, Y. Remote sensing of total suspended matter concentration in lakes across China using Landsat images and Google Earth Engine. ISPRS-J. Photogramm. Remote Sens. 2022, 187, 61–78. [Google Scholar] [CrossRef]
  13. Tang, Y.; Pan, Y.; Zhang, L.; Yi, H.; Gu, Y.; Sun, W. Efficient Monitoring of Total Suspended Matter in Urban Water Based on UAV Multi-spectral Images. Water Resour. Manag. 2023, 37, 2143–2160. [Google Scholar] [CrossRef]
  14. Wang, C.; Wang, D.; Yang, J.; Fu, S.; Li, D. Suspended Sediment within Estuaries and along Coasts: A Review of Spatial and Temporal Variations based on Remote Sensing. J. Coast. Res. 2020, 36, 1323–1331. [Google Scholar] [CrossRef]
  15. Xu, J.; Zhang, B.; Song, K.; Wang, Z.; Duan, H.; Chen, M.; Yang, F.; Li, F. Bio-Optical model of total suspended matter based on reflectance in the near infrared wave band for Case-II waters. Spectrosc. Spectr. Anal. 2008, 28, 2273–2277. [Google Scholar]
  16. Werdell, P.J.; McKinna, L.I.; Boss, E.; Ackleson, S.G.; Craig, S.E.; Gregg, W.W.; Lee, Z.; Maritorena, S.; Roesler, C.S.; Rousseaux, C.S. An overview of approaches and challenges for retrieving marine inherent optical properties from ocean color remote sensing. Prog. Oceanogr. 2018, 160, 186–212. [Google Scholar] [CrossRef]
  17. Chen, S.; Han, L.; Chen, X.; Li, D.; Sun, L.; Li, Y. Estimating wide range Total Suspended Solids concentrations from MODIS 250-m imageries: An improved method. ISPRS-J. Photogramm. Remote Sens. 2015, 99, 58–69. [Google Scholar] [CrossRef]
  18. Hou, X.; Feng, L.; Duan, H.; Chen, X.; Sun, D.; Shi, K. Fifteen-year monitoring of the turbidity dynamics in large lakes and reservoirs in the middle and lower basin of the Yangtze River, China. Remote Sens. Environ. 2017, 190, 107–121. [Google Scholar] [CrossRef]
  19. Li, W.; Yang, Q.; Ma, Y.; Yang, Y.; Song, K.; Zhang, J.; Wen, Z.; Liu, G. Remote Sensing Estimation of Long-Term Total Suspended Matter Concentration from Landsat across Lake Qinghai. Water 2022, 14, 2498. [Google Scholar] [CrossRef]
  20. Wang, C.; Chen, S.; Li, D.; Wang, D.; Liu, W.; Yang, J. A Landsat-based model for retrieving total suspended solids concentration of estuaries and coasts in China. Geosci. Model Dev. 2017, 10, 4347–4365. [Google Scholar] [CrossRef]
  21. Adhikari, A.; Menon, H.B. A Bio-optical Numerical Approach for Remote Retrieval of Total Suspended Matter from Turbid Waters. J. Indian Soc. Remote Sens. 2022, 50, 1773–1786. [Google Scholar] [CrossRef]
  22. Pahlevan, N.; Sarkar, S.; Franz, B.; Balasubramanian, S.V.; He, J. Sentinel-2 MultiSpectral Instrument (MSI) data processing for aquatic science applications: Demonstrations and validations. Remote Sens. Environ. 2017, 201, 47–56. [Google Scholar] [CrossRef]
  23. Jiang, D.; Matsushita, B.; Pahlevan, N.; Gurlin, D.; Lehmann, M.K.; Fichot, C.G.; Schalles, J.; Loisel, H.; Binding, C.; Zhang, Y. Remotely estimating total suspended solids concentration in clear to extremely turbid waters using a novel semi-analytical method. Remote Sens. Environ. 2021, 258, 112386. [Google Scholar] [CrossRef]
  24. Liu, Z.; Chen, J.; Cui, T.; Tang, J.; Wang, L. A combined semi-analytical algorithm for retrieving total suspended sediment concentration from multiple missions: A case study of the China Eastern Coastal Zone. Int. J. Remote Sens. 2021, 42, 8004–8033. [Google Scholar] [CrossRef]
  25. Kang, L.; Zhu, G.; Zhu, M.; Xu, H.; Zou, W.; Xiao, M.; Zhang, Y.; Qin, B. Bloom-induced internal release controlling phosphorus dynamics in large shallow eutrophic Lake Taihu, China. Environ. Res. 2023, 231, 116251. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, W.; Matsushita, B.; Chen, J.; Fukushima, T. A relaxed matrix inversion method for retrieving water constituent concentrations in case II waters: The case of Lake Kasumigaura, Japan. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3381–3392. [Google Scholar] [CrossRef]
  27. Liu, C.-C.; Miller, R.L. Spectrum matching method for estimating the chlorophyll-a concentration, CDOM ratio, and backscatter fraction from remote sensing of ocean color. Can. J. Remote Sens. 2008, 34, 343–355. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Ma, R.; Duan, H.; Loiselle, S.; Xu, J. A spectral decomposition algorithm for estimating chlorophyll-a concentrations in Lake Taihu, China. Remote Sens. 2014, 6, 5090–5106. [Google Scholar] [CrossRef]
  29. Qin, B.; Deng, J.; Shi, K.; Wang, J.; Brookes, J.; Zhou, J.; Zhang, Y.; Zhu, G.; Paerl, H.W.; Wu, L. Extreme climate anomalies enhancing cyanobacterial blooms in eutrophic Lake Taihu, China. Water Resour. Res. 2021, 57, e2020WR029371. [Google Scholar] [CrossRef]
  30. Ma, R.; Duan, H.; Liu, Q.; Loiselle, S.A. Approximate bottom contribution to remote sensing reflectance in Taihu Lake, China. J. Gt. Lakes Res. 2011, 37, 18–25. [Google Scholar] [CrossRef]
  31. Clementson, L.A.; Parslow, J.S.; Turnbull, A.R.; McKenzie, D.C.; Rathbone, C.E. Optical properties of waters in the Australasian sector of the Southern Ocean. J. Geophys. Res. Oceans 2001, 106, 31611–31625. [Google Scholar] [CrossRef]
  32. Stramska, M.; Stramski, D.; Mitchell, B.G.; Mobley, C.D. Estimation of the absorption and backscattering coefficients from inߚwater radiometric measurements. Limnol. Oceanogr. 2000, 45, 628–641. [Google Scholar] [CrossRef]
  33. Kim, J.; Jang, W.; Kim, J.H.; Lee, J.; Cho, K.H.; Lee, Y.-G.; Chon, K.; Park, S.; Pyo, J.; Park, Y. Application of airborne hyperspectral imagery to retrieve spatiotemporal CDOM distribution using machine learning in a reservoir. Int. J. Appl. Earth Obs. Geoinf. 2022, 114, 103053. [Google Scholar] [CrossRef]
  34. Zhang, F.; Li, J.; Wang, C.; Wang, S.; Wang, Z.; Zhang, B. Multitype inland water atmospheric correction and water quality estimation based on HY-1C CZI images. Natl. Remote Sens. Bull. 2023, 27, 79–91. [Google Scholar] [CrossRef]
  35. Tang, J.-W.; Tian, G.-L.; Wang, X.-Y.; Wang, X.-M.; Song, Q.-J. The methods of water spectra measurement and analysis I: Above-Water method. J. Remote Sens. 2021, 1, 37–44. [Google Scholar]
  36. Pancorbo, J.; Lamb, B.T.; Quemada, M.; Hively, W.D.; Gonzalez-Fernandez, I.; Molina, I. Sentinel-2 and WorldView-3 atmospheric correction and signal normalization based on ground-truth spectroradiometric measurements. ISPRS-J. Photogramm. Remote Sens. 2021, 173, 166–180. [Google Scholar] [CrossRef]
  37. Wang, R.; Shen, Q.; Peng, H.; Yao, Y.; Li, J.; Wang, M.; Shi, J.; Xu, W. Study on the applicability of multi-source high-resolution satellite images for monitoring black and odorous water body. Natl. Remote Sens. Bull. 2022, 26, 179–192. [Google Scholar] [CrossRef]
  38. Zhao, Z.; Yang, J.; Wang, M.; Chen, J.; Sun, C.; Song, N.; Wang, J.; Feng, S. The PCA-NDWI Urban Water Extraction Model Based on Hyperspectral Remote Sensing. Water 2024, 16, 963. [Google Scholar] [CrossRef]
  39. Werdell, P.J.; Franz, B.A.; Bailey, S.W.; Feldman, G.C.; Boss, E.; Brando, V.E.; Dowell, M.; Hirata, T.; Lavender, S.J.; Lee, Z. Generalized ocean color inversion model for retrieving marine inherent optical properties. Appl. Optics. 2013, 52, 2019–2037. [Google Scholar] [CrossRef]
  40. Pope, R.M.; Fry, E.S. Absorption spectrum (380–700 nm) of pure water. II. Integrating cavity measurements. Appl. Opt. 1997, 36, 8710–8723. [Google Scholar] [CrossRef]
  41. Mohd-Shazali, S.M.; Madihah, J.-S.; Ali, N.; Cheng-Ann, C.; Brewin, R.J.; Idris, M.S.; Noir, P.P. Dynamics of absorption properties of CDOM and its composition in Likas estuary, North Borneo, Malaysia. Oceanologia 2022, 64, 583–594. [Google Scholar] [CrossRef]
  42. Smith, R.C.; Baker, K.S. Optical properties of the clearest natural waters (200–800 nm). Appl. Opt. 1981, 20, 177–184. [Google Scholar] [CrossRef] [PubMed]
  43. Li, Y.; Wang, Q.; Huang, J.; Lyu, H.; Wei, Y. Optical Characteristics and Remote Sensing of Water Color in Lake Taihu; Science Press: Beijing, China, 2010; pp. 129–137. [Google Scholar]
  44. Miller, J.D.; Thode, A.E. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. Environ. 2007, 109, 66–80. [Google Scholar] [CrossRef]
  45. Su, H.; Lu, X.; Chen, Z.; Zhang, H.; Lu, W.; Wu, W. Estimating coastal chlorophyll-a concentration from time-series OLCI data based on machine learning. Remote Sens. 2021, 13, 576. [Google Scholar] [CrossRef]
  46. Wang, J.-H.; Li, C.; Xu, Y.-P.; Li, S.-Y.; Du, J.-S.; Han, Y.-P.; Hu, H.-Y. Identifying major contributors to algal blooms in Lake Dianchi by analyzing river-lake water quality correlations in the watershed. J. Clean. Prod. 2021, 315, 128144. [Google Scholar] [CrossRef]
  47. Shi, K.; Zhang, Y.; Zhu, G.; Liu, X.; Zhou, Y.; Xu, H.; Qin, B.; Liu, G.; Li, Y. Long-term remote monitoring of total suspended matter concentration in Lake Taihu using 250 m MODIS-Aqua data. Remote Sens. 2015, 164, 43–56. [Google Scholar]
  48. Li, J.; Chen, X.; Tian, L.; Huang, J.; Feng, L. Improved capabilities of the Chinese high-resolution remote sensing satellite GF-1 for monitoring suspended particulate matter (SPM) in inland waters: Radiometric and spatial considerations. ISPRS-J. Photogramm. Remote Sens. 2015, 106, 145–156. [Google Scholar] [CrossRef]
Figure 1. Study area and measured sample point distribution.
Figure 1. Study area and measured sample point distribution.
Remotesensing 16 01385 g001
Figure 2. Absorption spectrum.
Figure 2. Absorption spectrum.
Remotesensing 16 01385 g002
Figure 3. Scattering spectrum.
Figure 3. Scattering spectrum.
Remotesensing 16 01385 g003
Figure 4. SDGSAT-1 MII spectral response function.
Figure 4. SDGSAT-1 MII spectral response function.
Remotesensing 16 01385 g004
Figure 5. The simulated R rs ( λ ) is fitted with the measured R rs ( λ ) : (a) in the 400–800 nm and (b) at band center wavelength.
Figure 5. The simulated R rs ( λ ) is fitted with the measured R rs ( λ ) : (a) in the 400–800 nm and (b) at band center wavelength.
Remotesensing 16 01385 g005
Figure 6. The measured concentration data of 54 sample points and the correlation coefficient between the simulated R rs ( λ ) and the measured R rs ( λ ) for 54 sample points.
Figure 6. The measured concentration data of 54 sample points and the correlation coefficient between the simulated R rs ( λ ) and the measured R rs ( λ ) for 54 sample points.
Remotesensing 16 01385 g006
Figure 7. R rs ( λ ) and R rs Δ ( λ ) under single-component concentration change. (ac) are R rs ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) ; (df) are R rs Δ ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) .
Figure 7. R rs ( λ ) and R rs Δ ( λ ) under single-component concentration change. (ac) are R rs ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) ; (df) are R rs Δ ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) .
Remotesensing 16 01385 g007
Figure 8. R rs ( λ ) and R rs Δ ( λ ) under single-component concentration change, λ = 448 nm, 509 nm, 569 nm, 668 nm, and 773 nm. (ac) are R rs ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) ; (df) are R rs Δ ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) .
Figure 8. R rs ( λ ) and R rs Δ ( λ ) under single-component concentration change, λ = 448 nm, 509 nm, 569 nm, 668 nm, and 773 nm. (ac) are R rs ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) ; (df) are R rs Δ ( λ ) obtained under different C tsm , C chla , and a cdom ( 440 ) .
Remotesensing 16 01385 g008
Figure 9. Bands and band combinations that are not related to C cdom but to C tsm and C chla .
Figure 9. Bands and band combinations that are not related to C cdom but to C tsm and C chla .
Remotesensing 16 01385 g009
Figure 10. The relationship between C tsm and the single-component remote sensing reflectance contributed by C tsm : (a) ( B 6 / 3 ) tsm and (b) ( B 6 / 5 ) tsm .
Figure 10. The relationship between C tsm and the single-component remote sensing reflectance contributed by C tsm : (a) ( B 6 / 3 ) tsm and (b) ( B 6 / 5 ) tsm .
Remotesensing 16 01385 g010
Figure 11. The relationship between C chla and the single-component remote sensing reflectance contributed by C tsm : (a) ( B 6 / 3 ) chla and (b) ( B 6 / 5 ) chla .
Figure 11. The relationship between C chla and the single-component remote sensing reflectance contributed by C tsm : (a) ( B 6 / 3 ) chla and (b) ( B 6 / 5 ) chla .
Remotesensing 16 01385 g011
Figure 12. The change in the difference between each set of output value C tsm ( m + 1 ) and input value C tsm ( m ) during iteration.: (a) C tsm ( 0 ) = 1 g/m3; (b) C tsm ( 0 ) = 10 g/m3; (c) C tsm ( 0 ) = 50 g/m3; (d) C tsm ( 0 ) = 100 g/m3; (e) C tsm ( 0 ) = 200 g/m3; and (f) C tsm ( 0 ) = 500 g/m3.
Figure 12. The change in the difference between each set of output value C tsm ( m + 1 ) and input value C tsm ( m ) during iteration.: (a) C tsm ( 0 ) = 1 g/m3; (b) C tsm ( 0 ) = 10 g/m3; (c) C tsm ( 0 ) = 50 g/m3; (d) C tsm ( 0 ) = 100 g/m3; (e) C tsm ( 0 ) = 200 g/m3; and (f) C tsm ( 0 ) = 500 g/m3.
Remotesensing 16 01385 g012
Figure 13. The change in RE in the process of iterative convergence.
Figure 13. The change in RE in the process of iterative convergence.
Remotesensing 16 01385 g013
Figure 14. The measured C tsm are fitted to the retrieved C tsm of the iterative model.
Figure 14. The measured C tsm are fitted to the retrieved C tsm of the iterative model.
Remotesensing 16 01385 g014
Figure 15. The measured C tsm is fitted to the retrieved C tsm from the SDGSAT-1 MII image, based on the iterative model.
Figure 15. The measured C tsm is fitted to the retrieved C tsm from the SDGSAT-1 MII image, based on the iterative model.
Remotesensing 16 01385 g015
Figure 16. Change in the difference between the output value and the input value during the iteration process: (a) based on B 6 and B 6 / 5 ; (b) based on B 6 / 2 and B 6 / 5 ; and (c) based on B 6 / 4 and B 6 / 5 .
Figure 16. Change in the difference between the output value and the input value during the iteration process: (a) based on B 6 and B 6 / 5 ; (b) based on B 6 / 2 and B 6 / 5 ; and (c) based on B 6 / 4 and B 6 / 5 .
Remotesensing 16 01385 g016
Figure 17. Fitting of measured C tsm and retrieved C tsm when models are applied to the modeled data.
Figure 17. Fitting of measured C tsm and retrieved C tsm when models are applied to the modeled data.
Remotesensing 16 01385 g017
Figure 18. Fitting of measured C tsm and retrieved C tsm when models are applied to the SDGSAT-1 MII image of 27 July 2022.
Figure 18. Fitting of measured C tsm and retrieved C tsm when models are applied to the SDGSAT-1 MII image of 27 July 2022.
Remotesensing 16 01385 g018
Figure 19. Fitting of measured C tsm and retrieved C tsm when the adapted models are applied to the modeled data [18,47,48].
Figure 19. Fitting of measured C tsm and retrieved C tsm when the adapted models are applied to the modeled data [18,47,48].
Remotesensing 16 01385 g019
Figure 20. Fitting of measured C tsm and retrieved C tsm when the adapted models are applied to the SDGSAT-1 MII image of 27 July 2022 [18,47,48].
Figure 20. Fitting of measured C tsm and retrieved C tsm when the adapted models are applied to the SDGSAT-1 MII image of 27 July 2022 [18,47,48].
Remotesensing 16 01385 g020
Figure 21. Spatial distribution of C tsm in Lake Taihu, 2022–2023.
Figure 21. Spatial distribution of C tsm in Lake Taihu, 2022–2023.
Remotesensing 16 01385 g021
Figure 22. Quarterly statistics of the proportion of pixels in each range of C tsm .
Figure 22. Quarterly statistics of the proportion of pixels in each range of C tsm .
Remotesensing 16 01385 g022
Table 1. Band center wavelength and band gain of SDGSAT-1 MII.
Table 1. Band center wavelength and band gain of SDGSAT-1 MII.
BandCenter Wavelength (μm)Gain
By December 2022After December 2022
B10.4060.0515601330.052084676
B20.4480.0362413530.038928845
B30.5090.0233168350.025864978
B40.5690.0158496660.017501881
B50.6680.0160963810.016499392
B60.7730.0197190390.021554446
B70.8480.0138114580.015360482
Table 2. Experimental design for characteristic analysis of influence of remote sensing reflectance.
Table 2. Experimental design for characteristic analysis of influence of remote sensing reflectance.
VariableInitial ConcentrationConcentration Change Settings
C tsm (g/m3)1020, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150
C chla (mg/m3)0.21, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140
a cdom ( 440 ) (m−1)0.20.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4
Table 3. Error between simulated R rs ( λ ) and measured R rs ( λ ) at the center wavelength of the band.
Table 3. Error between simulated R rs ( λ ) and measured R rs ( λ ) at the center wavelength of the band.
R rs ( λ ) R2RMSE (1/sr)
R rs ( 448 ) 0.8260.0025
R rs ( 509 ) 0.8370.0031
R rs ( 569 ) 0.8240.0032
R rs ( 668 ) 0.8620.0038
R rs ( 773 ) 0.9060.0029
Table 4. Model input parameter setting.
Table 4. Model input parameter setting.
ClarificationModel Application Scenarios
One-Dimensional Reflectance ValuesSDGSAT-1 MII Images
C tsm ( 0 ) a set of one-dimensional valuesa remote sensing reflectance image
B 6 / 3
B 6 / 5
calculation methoditerative on a value-by-value basisiterative image-by-image pixel
Table 5. Error statistics for four iterative inversion models.
Table 5. Error statistics for four iterative inversion models.
Characteristic Band CombinationsThe Modeled DataThe SDGSAT-1 MII Image
R2RMSE (g/m3)MAPE
(%)
R2RMSE (g/m3)MAPE
(%)
B 6   and   B 6 / 5 0.918.0516.930.828.4831.06
B 6 / 3   and   B 6 / 5 0.974.8911.480.873.928.13
B 6 / 2   and   B 6 / 5 0.956.0314.050.848.6730.47
B 6 / 4   and   B 6 / 5 0.965.4912.580.866.6323.06
Table 6. The model and accuracy of C tsm inversion model after adaption.
Table 6. The model and accuracy of C tsm inversion model after adaption.
Model NameTypologyFormulasEvaluation Indicators
Shi_2015 [47]original modelCtsm = 9.65 e 58.81 R rs ( 645 ) R2 = 0.70; RMSE = 14.3; MAPE = 23
adapted modelCtsm = 8.98 e 53.766 R rs ( 645 ) R2 = 0.66; RMSE = 16.96; MAPE = 27.01
Li_2015 [48]original modelCtsm = −20.7 + 2.8 e 61.9 R rs ( 660 )  R2 = 0.95; RMSE = 14.7
adapted modelCtsm = 12.05 e 53.177 R rs ( 660 ) R2 = 0.64; RMSE = 17.59; MAPE = 29.45
Hou_2017 [18]original modelCtsm = 132.83 ( R rs ( 645 ) / R rs ( 555 ) ) 2 − 52.618 ( R rs ( 645 ) / R rs ( 555 ) )  R2 = 0.88; RMSE = 34.2; MRE = 33.79
adapted modelCtsm = 1.21 e 4.11 ( R rs ( 645 ) / R rs ( 555 ) ) R2 = 0.82; RMSE = 12.44; MAPE = 16.97
Table 7. Error statistics when models are applied to the SDGSAT-1 MII image of 27 July 2022.
Table 7. Error statistics when models are applied to the SDGSAT-1 MII image of 27 July 2022.
ModelRMSE (g/m3)MAPE (%)
in this paper3.928.13
Shi_2015 [47]10.9137.9
Li_2015 [48]18.0057.62
Hou_2017 [18]15.2444.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, X.; Li, J.; Sun, Y.; Bao, Y.; Sun, Y.; Chen, X.; Yan, Y. Retrieval of Total Suspended Matter Concentration Based on the Iterative Analysis of Multiple Equations: A Case Study of a Lake Taihu Image from the First Sustainable Development Goals Science Satellite’s Multispectral Imager for Inshore. Remote Sens. 2024, 16, 1385. https://doi.org/10.3390/rs16081385

AMA Style

Hu X, Li J, Sun Y, Bao Y, Sun Y, Chen X, Yan Y. Retrieval of Total Suspended Matter Concentration Based on the Iterative Analysis of Multiple Equations: A Case Study of a Lake Taihu Image from the First Sustainable Development Goals Science Satellite’s Multispectral Imager for Inshore. Remote Sensing. 2024; 16(8):1385. https://doi.org/10.3390/rs16081385

Chicago/Turabian Style

Hu, Xueke, Jiaguo Li, Yuan Sun, Yunfei Bao, Yonghua Sun, Xingfeng Chen, and Yueguan Yan. 2024. "Retrieval of Total Suspended Matter Concentration Based on the Iterative Analysis of Multiple Equations: A Case Study of a Lake Taihu Image from the First Sustainable Development Goals Science Satellite’s Multispectral Imager for Inshore" Remote Sensing 16, no. 8: 1385. https://doi.org/10.3390/rs16081385

APA Style

Hu, X., Li, J., Sun, Y., Bao, Y., Sun, Y., Chen, X., & Yan, Y. (2024). Retrieval of Total Suspended Matter Concentration Based on the Iterative Analysis of Multiple Equations: A Case Study of a Lake Taihu Image from the First Sustainable Development Goals Science Satellite’s Multispectral Imager for Inshore. Remote Sensing, 16(8), 1385. https://doi.org/10.3390/rs16081385

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