the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
HHU24SWDSCS: A shallow-water depth model over island areas in South China Sea retrieved from Satellite-derived bathymetry
Abstract. Accurate shallow-water depth information for island areas is crucial for maritime safety, resource exploration, ecological conservation, and offshore economic activity. Traditional approaches like shipborne sounding and airborne bathymetric light detection and ranging (LiDAR) surveys are expensive, time-consuming, and are limited in politically sensitive regions. Moreover, satellite altimetry-predicted depths exhibit large errors over shallow waters. In contrast, satellite-derived bathymetry (SDB), estimated from multispectral imagery, provides a rapid, open source, and cost-effective technique to fully characterize the bathymetry of a region. Given the scarcity of in-situ water-depth data for the South China Sea (SCS), a shallow-water depth model, HHU24SWDSCS, was developed by integrating 1298 Ice, Cloud, and land Elevation Satellite (ICESat-2) tracks with 70 Sentinel-2 multispectral images. The model covers >120 islands and reefs in the SCS, with a resolution of 10 m. Validation against independent ICESat-2 depth data produced a root mean square error for the model of 0.81–1.35 m (<5 % of the maximum depth), with an average coefficient of determination of 0.91. Validation against independent airborne LiDAR bathymetry data revealed an accuracy of 1.01 m for the Lingyang Reef. Further comparisons with existing bathymetry models revealed the superior performance of the model. While the existing bathymetry models exhibit errors up to tens of meters or larger for island regions, and should therefore be used with caution, the HHU24SWDSCS model exhibited good accuracy in shallow waters across the SCS. This model thus provides a reference for mapping shallow-water depth close to islands and provides fundamental support for research in oceanography, geodesy, and other disciplines.
- Preprint
(12948 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 08 Jan 2025)
-
RC1: 'Comment on essd-2024-443', Anonymous Referee #1, 09 Dec 2024
reply
The paper uses ICESat data for a training data set for a linear model (multiple linear regression) for SDB from Sentinel-2. The paper shows good results for estimating SDB for several areas in the South China sea. Demonstrations of methods for applying ICESat data for tuning, and for execution are a useful addition to the SDB field.
The methods section has insufficient information to reproduce the application to the satellite. That is necessary, given this is a methods paper. The most noteworthy is how the apparently multiple h_0 and h_i coefficients developed for each ICESat trackline were applied around an island. This should come right after line 340. This is non-trivial, there are 3 h_i coefficients (3 bands).
Also, the limitations of some of the statistical validation need to be identified.
Methods:
The Sentinel-2 reflectance source and atmospheric correction are not explained.
Was one image used for each location?
The selection of areas for the regression is not clear. Perhaps because text is split up between lines 170 and 320-335.
Line 170 What does this mean? “We used the GEBCO_2023 model to identify and remove deep-water effects (>100 m) in SDB estimation”. Weren’t the NDWI and ICESat depths used to do this?
Line 320-335. This is not clear. Dividing the track into segments “based on water depth variation trend (from ascend to descend)”. Divide how? This is critical to how the correlation coefficients will be determined. The red band will disappear much sooner than green or blue. “Ascend to descend” should be changed, they are actions so it doesn’t make sense. “Shallow to deep”?
After that section, how is SDB determined for the whole island? H_0 and h_i were determined for each track “segment”. Then what? Were they interpolated or averaged?
And were the correlation coefficients determined for all locations on the track within the shallow water mask. The shallow water mask was determined by the intersection of the NDWI and ICESat?. And how was ICESat screened, line 174?
The split of data was “80% training and 20% validation”. What does this mean? Was this random, were non-overlapping ICESat transects left out of training? If not, and a random split was used, the validation is not independent. It fails to consider spatial autocorrelation (there are a lot of papers on this topic), which would bias in favor of the results. The study “validation” does not need to be redone, but this problem needs to be clearly identified, and text calling it a validation should be changed. Perhaps saying that “Model consistency was evaluated. “
Figure 12 and Table 3 do provide one independent validation, as the lidar was not used for training.
On statistisics. In spite of the popularity of Rsquared (R2) as a validation metric, it is both a poor error metric and it is redundant to RMSE (and so unnecessary). And R2 cannot be compared for samples with different ranges (variance in X). Many statisticians have reported this; King 1986 (https://www.jstor.org/stable/2111095) is a good example. There are several descriptions of the problem on the web (R2 is the fit of the line against the variance in the data, so a wider range of data will have a higher R2). Figure 9 shows the problem. Compare 9e to 9b. Occurring to R2, 9e (0.938) outperforms 9b (0.878). However, 9e has twice the error, 1.631 vs 0.802 for 9b. R2 does not provide useful information. Why? 9e has twice the range of depths, so the squared variance is much greater. It’s ok to leave the R2 in the figures, because there are people who are desperate to see it, but leave any comparisons of R2 out of the text. Remove R2 reference from 362-374, 410-425, 529. This problem should be stated at line 341: e.g., “R2 is actually redundant with RMSE. However, R2 also varies with data range, so unlike RMSE, R2 values cannot be meaningfully compared between different samples. R2 values are included because they are familiar.“
Line 335. “The effects of deep-water areas were then removed to minimize the influence of bottom reflection on SDB estimation”. What effects were removed from what? Does this text belong before line 312? (“average deep-water reflectance”).
Figure 12 and 13 captions are not clear, please include the letters in the caption. It would be even better to label each column.
Citation: https://doi.org/10.5194/essd-2024-443-RC1
Data sets
HHU24SWDSCS: A shallow-water depth model over island areas in South China Sea retrieved from Satellite-derived bathymetry Yihao Wu et al. https://doi.org/10.5281/zenodo.13852568
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
208 | 48 | 6 | 262 | 4 | 4 |
- HTML: 208
- PDF: 48
- XML: 6
- Total: 262
- BibTeX: 4
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1