基于数据-知识驱动的高精度海底地形绘制:以南海为例

刘洋, 李三忠, 邹卓延, 索艳慧, 孙毅

海洋学研究 ›› 2024, Vol. 42 ›› Issue (3) : 142-152.

PDF(3329 KB)
PDF(3329 KB)
海洋学研究 ›› 2024, Vol. 42 ›› Issue (3) : 142-152. DOI: 10.3969/j.issn.1001-909X.2024.03.012
研究论文

基于数据-知识驱动的高精度海底地形绘制:以南海为例

作者信息 +

High-precision seafloor topographic mapping based on data-knowledge-driven: An example from the South China Sea

Author information +
文章历史 +

摘要

海底地形具有非常重要的商业、工程、军事和科学研究价值。目前,常用重力场数据反演海底地形,如自由空气重力异常和垂直重力梯度。然而,由于现有方法反演海底地形具有较强的多解性,仍然无法准确获取高精度的海底地形。该文提出了重力-密度法与随机森林结合的数据-知识驱动新方法,以重建准确的海底地形。该方法在中国南海海域进行了测试,并与重力-密度法、随机森林以及现有的SIO模型进行了对比分析。反演结果显示,数据-知识驱动提供了更好的反演性能,随机森林和重力-密度法次之,SIO模型最差。相比于重力-密度法,数据-知识驱动的平均绝对误差、平均相对误差和均方根误差分别降低了21%、25%和7%;而相比于随机森林,它们分别也降低了20%、20%和20%。此外,数据-知识驱动模型与船载测深数据具有较高的一致性,其差值大约有72%分布在±10 m范围内,占比高于其他三种模型。该结果证明了数据-知识驱动方法在海底地形反演中的可行性和有效性,有助于加快高精度海底地形的绘制。

Abstract

Seafloor topography is of considerable value in commercial, engineering, military and scientific research. Currently, gravity field data, such as free air gravity anomalies and vertical gravity gradients, are commonly used to inverse seafloor topography. However, due to the strong multi-resolution of the existing methods to inverse seafloor topography, it is still impossible to obtain accurate high-precision seafloor topography. A new data-knowledge-driven method was proposed to reconstruct accurate seafloor topography, which combines the gravity-density method with random forests. This method was applied to the South China Sea and compared with the gravity-density, random forest, and existing SIO models. The inversion results show that the data-knowledge-driven method provides better inversion performance, followed by the random forest and gravity-density methods, and the SIO model is the worst. The mean absolute error, mean relative error and root mean square error of the data-knowledge-driven are reduced by 21%, 25% and 7%, respectively, compared to those of the gravity-density method, while they are also reduced by 20%, 20% and 20%, respectively, compared to those of the random forest. In addition, the data-knowledge-driven model has a high degree of agreement with the shipboard bathymetry data, with approximately 72% of its differences distributed within ±10 m, which is higher than that of the other three models. The results demonstrate the feasibility and effectiveness of the data-knowledge-driven approach in seafloor topography inversion, which helps to accelerate the mapping of high-precision seafloor topography.

关键词

海底地形 / 机器学习 / 数据驱动 / 知识驱动 / 重力-密度法 / 随机森林 / SIO模型 / 船载测深

Key words

seafloor topography / machine learning / data-driven / knowledge-driven / gravity-density method / random forest / SIO model / shipboard bathymetry

引用本文

导出引用
刘洋, 李三忠, 邹卓延, . 基于数据-知识驱动的高精度海底地形绘制:以南海为例[J]. 海洋学研究. 2024, 42(3): 142-152 https://doi.org/10.3969/j.issn.1001-909X.2024.03.012
LIU Yang, LI Sanzhong, ZOU Zhuoyan, et al. High-precision seafloor topographic mapping based on data-knowledge-driven: An example from the South China Sea[J]. Journal of Marine Sciences. 2024, 42(3): 142-152 https://doi.org/10.3969/j.issn.1001-909X.2024.03.012
中图分类号: P714.7   

参考文献

[1]
WÖLFL A C, SNAITH H, AMIREBRAHIMI S, et al. Seafloor mapping-the challenge of a truly global ocean bathymetry[J]. Frontiers in Marine Science, 2019, 6: 283.
[2]
SANDWELL D T, SMITH W H, GILLE S, et al. Bathymetry from space: White paper in support of a high-resolution, ocean altimeter mission[J]. International Geophysics Series, 2001, 69: 1049-1062.
[3]
SMITH W. Introduction to this special issue on bathymetry from space[J]. Oceanography, 2004, 17(1): 6-7.
[4]
PEERI S, GARDNER J V, WARD L G, et al. The seafloor: A key factor in lidar bottom detection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(3): 1150-1157.
[5]
EUGENIO F, MARCELLO J, MARTIN J. High-resolution maps of bathymetry and benthic habitats in shallow-water environments using multispectral remote sensing imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(7): 3539-3549.
[6]
黄谟涛, 翟国君, 欧阳永忠, 等. 卫星测高技术应用研究回顾与展望[J]. 海洋测绘, 2004, 24(4):65-70.
HUANG M T, ZHAI G J, OU-YANG Y Z, et al. Review and prospect of application study in altimetry[J]. Hydrographic Surveying and Charting, 2004, 24(4): 65-70.
[7]
CALMANT S, BERGE-NGUYEN M, CAZENAVE A. Global seafloor topography from a least-squares inversion of altimetry-based high-resolution mean sea surface and shipboard soundings[J]. Geophysical Journal International, 2002, 151(3): 795-808.
[8]
TOZER B, SANDWELL D T, SMITH W H F, et al. Global bathymetry and topography at 15Arc Sec: SRTM15+[J]. Earth and Space Science, 2019, 6(10): 1847-1864.
[9]
AMANTE C, EAKINS B W. ETOPO1 arc-minute global relief model: Procedures, data sources and analysis[R]//NOAA Technical Memorandum NESDIS NGDC-24. NOAA, National Geophysical Data Center, 2009.
[10]
MARKS K M, SMITH W H F, SANDWELL D T. Evolution of errors in the altimetric bathymetry model used by Google Earth and GEBCO[J]. Marine Geophysical Researches, 2010, 31(3): 223-238.
[11]
SMITH W H F, SANDWELL D T. Global sea floor topography from satellite altimetry and ship depth soundings[J]. Science, 1997, 277(5334): 1956-1962.
[12]
ANDERSEN O B, KNUDSEN P. The DNSC08BAT bathy-metry developed from satellite altimetry[C]// Proceedings of the EGU-2008 Meeting: Vienna, Austria, 2008.
[13]
BECKER J J, SANDWELL D T, SMITH W H F, et al. Global bathymetry and elevation data at 30 arc seconds resolution: SRTM30_PLUS[J]. Marine Geodesy, 2009, 32(4): 355-371.
[14]
胡敏章, 张胜军, 金涛勇, 等. 新一代全球海底地形模型BAT_WHU2020[J]. 测绘学报, 2020, 49(8):939-954.
摘要
本文利用由多源卫星测高资料计算的新版全球重力异常Grav_Alti_WHU,联合船测水深资料,构建了全球75°S—70°N范围的1'×1'海底地形模型BAT_WHU2020。以船测水深、现有模型和多波束测深数据为参考,对模型精度进行了分析评价。结果表明,在中国海域及邻区(104°E—160°E,0°N—50°N),本文模型与船测水深之差值的标准差约70 m,与SIO V19.1模型精度相当,优于ETOPO1、DTU10、GEBCO_08等模型,较此前发布的BAT_VGG模型精度提高了约30%,说明本文模型构建方法可靠、数据处理准确、精度较高。在全球范围内,BAT_WHU2020模型与船测水深之差值的标准差为50~65 m,差值在±200 m范围内的比率超过95%,与SIO V19.1模型精度相当,优于ETOPO1、DTU10、GEBCO_08等模型,较BAT_VGG模型精度提高了27%~36%。以SIO V19.1模型为参考,模型之差的标准差为90~110 m,约90%格网点差值在±200 m以内,约95%格网点差值在±300 m以内,两者一致性良好。最后,讨论了地壳均衡、Parker公式高次项等对成果精度的影响,模型的真实空间分辨率,以及以多波束测深为参考的模型精度问题。分析认为,BAT_WHU2020模型空间分辨率为10~18 km,在马里亚纳海沟、麦夸里海岭地区相对精度为5%~6%。
HU M Z, ZHANG S J, JIN T Y, et al. A new generation of global bathymetry model BAT_WHU2020[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(8): 939-954.
In this paper, a 1'×1' bathymetry model BAT_WHU2020 in the range of 75°S—70°N is constructed by using the latest version of the global gravity anomaly model derived from the multi-source satellite altimetry data and the shipboard depths. The accuracy of the model is analyzed and evaluated based on ship depths, existing models and multibeam soundings. The standard deviation of the difference between the proposed model and the ship depths in China Sea and its adjacent areas (104°E—160°E, 0°N—50°N) is about 70 m, which is equivalent to the accuracy of SIO V19.1 model, superior to ETOPO1, DTU10, GEBCO_08 model, and about 30% higher than the accuracy of BAT_VGG model published before, which shows that the method in this paper is reliable, and the data processing is accurate and the accuracy is high. The standard deviation of the differences between BAT_WHU2020 model and ship depths is about 50~65 m globally, and the ratio of the difference within ±200 m is greater than 95%. It is showed that the accuracy of BAT_WHU2020 model is equivalent to SIO V19.1, better than ETOPO1, DTU10, GEBCO_08 model, and improved about 27%~36% from the BAT_VGG model. Comparing to SIO V19.1 model, the standard deviation of the model differences is about 90~110 m, about 90% of the grid differences is within 200 m, and about 95% is within 300 m. Finally, the effects of crustal isostasy and high order terms in the Parker's formula on the accuracy of the results, the accuracy of the model compared to multibeam soundings, and the spatial resolution are discussed. It indicates that the spatial resolution of BAT_WHU2020 model is about 10~18 km, and the relative accuracy is about 5%~6% around the Mariana trench and Macquarie ridge.
[15]
HSIAO Y S, KIM J W, KIM K B, et al. Bathymetry estimation using the gravity-geologic method: An investi-gation of density contrast predicted by the downward continuation method[J]. Terrestrial, Atmospheric and Oceanic Sciences, 2011, 22(3): 347.
[16]
KIM J W, VON FRESE R R B, LEE B Y, et al. Altimetry-derived gravity predictions of bathymetry by the gravity-geologic method[J]. Pure and Applied Geophysics, 2011, 168(5): 815-826.
[17]
欧阳明达. 利用海洋重力数据反演海底地形的理论与方法[D]. 郑州: 解放军信息工程大学, 2015.
OU-YANG M D. Theory and method of inversion of seabed topography using marine gravity data[D]. Zhengzhou: PLA Information Engineering University, 2015.
[18]
MCNUTT M. Compensation of oceanic topography: An application of the response function technique to the surveyor area[J]. Journal of Geophysical Research: Solid Earth, 1979, 84(B13): 7589-7598.
[19]
ARABELOS D, TZIAVOS I N. Gravity-field improvement in the Mediterranean Sea by estimating the bottom topography using collocation[J]. Journal of Geodesy, 1998, 72(3): 136-143.
[20]
KIM K B, HSIAO Y S, KIM J W, et al. Bathymetry enhancement by altimetry-derived gravity anomalies in the East Sea (Sea of Japan)[J]. Marine Geophysical Researches, 2010, 31(4): 285-298.
[21]
李倩倩, 鲍李峰. 测高重力场反演海底地形方法比较[J]. 海洋测绘, 2016, 36(5):1-4,18.
LI Q Q, BAO L F. Comparative analysis of methods for bathymetry prediction from altimeter-derived gravity anomalies[J]. Hydrographic Surveying and Charting, 2016, 36(5): 1-4, 18.
[22]
胡敏章, 李建成, 金涛勇. 应用重力地质方法反演皇帝海山的海底地形[J]. 武汉大学学报:信息科学版, 2012, 37(5):610-612,629.
HU M Z, LI J C, JIN T Y. Bathymetry inversion with gravity-geologic method in Emperor Seamount[J]. Geomatics and Information Science of Wuhan University, 2012, 37(5): 610-612, 629.
[23]
王永康, 周兴华, 唐秋华, 等. 应用重力地质法反演马里亚纳海沟地形[J]. 海洋科学进展, 2020, 38(4):708-716.
WANG Y K, ZHOU X H, TANG Q H, et al. Predicting bathymetry in Mariana trench using gravity-geologic method[J]. Advances in Marine Science, 2020, 38(4): 708-716.
[24]
欧阳明达, 孙中苗, 翟振和. 基于重力地质法的南中国海海底地形反演[J]. 地球物理学报, 2014, 57(9):2756-2765.
摘要
根据重力地质法(GGM),利用南中国海海域内63179个船测控制点水深将测高自由空间重力异常划分为长波参考场和短波残差场,并反演出了该海域112&deg;E&mdash;119&deg;E,12&deg;N&mdash;20&deg;N范围的1&rsquo;&times;1&rsquo;海底地形模型,该过程中使用的海水和海底洋壳密度差异常数1.32 g&middot;cm<sup>-3</sup>通过实测水深估计得到.利用反演得到的GGM模型对剩余的10529个检核点船测水深插值计算后与实测水深进行比较,其较差结果的均值为-1.64 m,标准差为76.95 m,相对精度为4.06%.此外,根据船测点数量、分布和海底地形的不同,选择了三个海域进行统计,结果表明:在船测控制点分布均匀的海域,GGM模型精度优于ETOPO1模型,在控制点过于分散的海域其精度会有所下降,但好于船测水深的直接格网化结果.为进一步探究检核点的较差结果中出现较大数值的成因,本文对精度较差的点位进行了单独分析,选择了两条船测航迹剖面进行了研究,并分析了检核点的水深较差、相对精度与水深和重力异常的关系,结果表明:GGM模型精度受水深和重力异常的相关性影响较小,受海底地形复杂程度影响较大,地形坡度变化平缓海域的预测精度明显高于海山地区.最后,综合GGM模型和ETOPO1模型优势,利用所有船测水深作为控制,生成了综合的海底地形模型.
OU-YANG M D, SUN Z M, ZHAI Z H. Predicting bathymetry in South China Sea using the gravity-geologic method[J]. Chinese Journal of Geophysics, 2014, 57(9): 2756-2765.
[25]
ANNAN R F, WAN X Y. Recovering bathymetry of the gulf of Guinea using altimetry-derived gravity field products combined via convolutional neural network[J]. Surveys in Geophysics, 2022, 43(5): 1541-1561.
[26]
SUN Y J, ZHENG W, LI Z W, et al. Improved the accuracy of seafloor topography from altimetry-derived gravity by the topography constraint factor weight optimization method[J]. Remote Sensing, 2021, 13(12): 2277.
[27]
AN D C, GUO J Y, LI Z, et al. Improved gravity-geologic method reliably removing the long-wavelength gravity effect of regional seafloor topography: A case of bathymetric prediction in the South China Sea[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 4211912.
[28]
CALMANT S, BAUDRY N. Modelling bathymetry by inverting satellite altimetry data: A review[J]. Marine Geophysical Researches, 1996, 18(2): 123-134.
[29]
WANG Y M. Predicting bathymetry from the earth’s gravity gradient anomalies[J]. Marine Geodesy, 2000, 23(4): 251-258.
[30]
REICHSTEIN M, CAMPS-VALLS G, STEVENS B, et al. Deep learning and process understanding for data-driven earth system science[J]. Nature, 2019, 566(7743): 195-204.
[31]
SHEN C P, APPLING A, GENTINE P, et al. Differentiable modelling to unify machine learning and physical models for geosciences[J]. Nature Reviews Earth & Environment, 2023, 4: 552-567.
[32]
YANG L, LIU M, LIU N, et al. Recovering bathymetry from satellite altimetry-derived gravity by fully connected deep neural network[J]. IEEE Geoscience and Remote Sensing Letters, 2023, 20: 1502805.
[33]
CHEN Y, ZHANG D. Integration of knowledge and data in machine learning[Z/OL]. arXiv, 2022: 2202. 10337. https://arxiv.org/abs/2202.10337.
[34]
KARPATNE A, ATLURI G, FAGHMOUS J H, et al. Theory-guided data science: A new paradigm for scientific discovery from data[J]. IEEE Transactions on Knowledge and Data Engineering, 2017, 29(10): 2318-2331.
[35]
GETTELMAN A, GEER A J, FORBES R M, et al. The future of earth system prediction: Advances in model-data fusion[J]. Science Advances, 2022, 8(14): eabn3488.
[36]
MENG C Z, SEO S, CAO D F, et al. When physics meets machine learning: A survey of physics-informed machine learning[Z/OL]. arXiv, 2022: 2203. 16797. https://arxiv.org/abs/2203.16797v1.
[37]
MORTON B, BLACKMORE G. South China Sea[J]. Marine Pollution Bulletin, 2001, 42(12): 1236-1263.
The South China Sea is poorly understood in terms of its marine biota, ecology and the human impacts upon it. What is known is most often contained in reports and workshop and conference documents that are not available to the wider scientific community. The South China Sea has an area of some 3.3 million km2 and depths range from the shallowest coastal fringe to 5377 m in the Manila Trench. It is also studded with numerous islets, atolls and reefs many of which are just awash at low tide. It is largely confined within the Tropic of Cancer and, therefore, experiences a monsoonal climate being influenced by the Southwest Monsoon in summer and the Northeast Monsoon in winter. The South China Sea is a marginal sea and, therefore, largely surrounded by land. Countries that have a major influence on and claims to the sea include China, Malaysia, the Philippines and Vietnam, although Thailand, Indonesia and Taiwan have some too. The coastal fringes of the South China Sea are home to about 270 million people that have had some of the fastest developing and most vibrant economies on the globe. Consequently, anthropogenic impacts, such as over-exploitation of resources and pollution, are anticipated to be huge although, in reality, relatively little is known about them. The Indo-West Pacific biogeographic province, at the centre of which the South China Sea lies, is probably the world's most diverse shallow-water marine area. Of three major nearshore habitat types, i.e., coral reefs, mangroves and seagrasses, 45 mangrove species out of a global total of 51, most of the currently recognised 70 coral genera and 20 of 50 known seagrass species have been recorded from the South China Sea. The island groups of the South China Sea are all disputed and sovereignty is claimed over them by a number of countries. Conflicts have in recent decades arisen over them because of perceived national rights. It is perhaps because of this that so little research has been undertaken on the South China Sea. What data are available, however, and if Hong Kong is used, as it is herein, as an indicator of what the perturbations of other regional cities upon the South China Sea are like, then it is impacted grossly and an ecological disaster has probably already, but unknowingly, happened.
[38]
WANG A M, DU Y, PENG S Q, et al. Deep water characteristics and circulation in the South China Sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 2018, 134: 55-63.
[39]
SANDWELL D T, MÜLLER R D, SMITH W H F, et al. Marine geophysics. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure[J]. Science, 2014, 346(6205): 65-67.
Gravity models are powerful tools for mapping tectonic structures, especially in the deep ocean basins where the topography remains unmapped by ships or is buried by thick sediment. We combined new radar altimeter measurements from satellites CryoSat-2 and Jason-1 with existing data to construct a global marine gravity model that is two times more accurate than previous models. We found an extinct spreading ridge in the Gulf of Mexico, a major propagating rift in the South Atlantic Ocean, abyssal hill fabric on slow-spreading ridges, and thousands of previously uncharted seamounts. These discoveries allow us to understand regional tectonic processes and highlight the importance of satellite-derived gravity models as one of the primary tools for the investigation of remote ocean basins. Copyright © 2014, American Association for the Advancement of Science.
[40]
BALMINO G, VALES N, BONVALOT S, et al. Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies[J]. Journal of Geodesy, 2012, 86(7): 499-520.
[41]
MAUS S, BARCKHAUSEN U, BERKENBOSCH H, et al. EMAG2: A 2-arc min resolution earth magnetic anomaly grid compiled from satellite, airborne, and marine magnetic measurements[J]. Geochemistry, Geophysics, Geosystems, 2009, 10(8). https://doi.org/10.1029/2009GC002471.
[42]
PARKER R L. The rapid calculation of potential anomalies[J]. Geophysical Journal International, 1973, 31(4): 447-455.
[43]
HWANG C. A bathymetric model for the South China sea from satellite altimetry and depth data[J]. Marine Geodesy, 1999, 22(1): 37-51.
[44]
KIM K B, YUN H S. Satellite-derived bathymetry prediction in shallow waters using the gravity-geologic method: A case study in the west sea of Korea[J]. KSCE Journal of Civil Engineering, 2018, 22(7): 2560-2568.
[45]
HSIAO Y S, HWANG C, CHENG Y S, et al. High-resolution depth and coastline over major atolls of South China Sea from satellite altimetry and imagery[J]. Remote Sensing of Environment, 2016, 176: 69-83.
[46]
STRYKOWSKI G, BOSCHETTI F, PAPP G. Estimation of the mass density contrasts and the 3D geometrical shape of the source bodies in the Yilgarn area, Eastern Goldfields, Western Australia[J]. Journal of Geodynamics, 2005, 39: 444-460.
[47]
BREIMAN L. Random forests[J]. Machine Learning, 2001, 45: 5-32.

基金

青岛海洋科学与技术试点国家实验室山东省专项经费项目(2022QNLM05032)
国家自然科学基金项目(42121005)
国家自然科学基金项目(91958214)
山东省自然科学基金项目(ZR2021ZD09)
111项目(B20048)

PDF(3329 KB)

Accesses

Citation

Detail

段落导航
相关文章

/