基于EGM模型的多波束无验潮测深技术应用

多波束水深测量中受潮汐因素的影响,测量垂直基准是变化的,具有瞬时性。传统多波束测量,需在测区内设立一个或多个验潮站进行同步水位观测,最终将水深归算到深度基准面上。针对多波束水深测量中垂直基准转换的复杂性问题,文中基于地球重力场模型,结合测区内实测的GNSS/水准数据,通过插值算法建立了测区范围内似大地水准面精化模型,构建了多波束无验潮水深测量的垂直基准转换模型。通过实例表明,该方法有效地消除了潮汐、动态吃水及涌浪等因素影响,直接获取深度基准面的水深值,提高工作效率,可满足近岸多波束水深测量的工作需求。     

Derivation of Bathymetry from High-resolution Optical Satellite Imagery and USV Sounding Data

Remote sensing bathymetry inversion can quickly obtain water depth data of large areas, but this process relies on a large number of in-situ depth data points. USV-based (Unmanned Surface Vehicle) technique can obtain the bathymetry data of shallow water where ordinary ships are inaccessible, but this technique is inefficient and generally only data along survey line can be collected. The combination of USV and high-resolution remote sensing provides a new solution for water depth surveying and mapping around an island. This paper focuses on the key techniques, using USV sounding data and GeoEye-1 multispectral remote sensing images covering the region of Wuzhizhou island in the experiment. The results show that the MAE (Mean Absolute Error) of USV sounding is 0.25 m, while the MRE (Mean Relative Error) is 1.41%, and the MRE of remote sensing bathymetry aided by USV sounding can be controlled within 20%. Errors are mainly from areas shallower than 5 m, and are also affected by the USV sounding position accuracy. It shows that it is feasible to combine the USV sounding and high-resolution remote sensing bathymetry, and this technique has broad application prospects in the field of bathymetry in large shallow areas.     

Predicting Bathymetry from the Earth's Gravity Gradient Anomalies

Historically, prediction of ocean floor depth, or bathymetry, has been based on the isostatic modeling and linearized relationships between gravity anomalies and bathymetry. The need for isostatic modeling limits the application of the resulting bathymetry predictions as constraints in geophysical models. An alternative technique making use of the Earth's vertical gravity gradient for predicting bathymetry is explored in this paper. This technique is based on the fact that the observed gravity gradient anomalies result primarily from local mass concentrations on the ocean floor, and that mass compensation by the oceanic crust has an insignificant effect on the gravity gradients, and can be neglected. The resulting bathymetry prediction therefore is independent of isostatic modeling assumptions, allowing it to be used as a constraint on models of lithospheric compensation and for other geodetic and geophysical applications.     

Predicting Bathymetry from the Earth's Gravity Gradient Anomalies

Historically, prediction of ocean floor depth, or bathymetry, has been based on the isostatic modeling and linearized relationships between gravity anomalies and bathymetry. The need for isostatic modeling limits the application of the resulting bathymetry predictions as constraints in geophysical models. An alternative technique making use of the Earth's vertical gravity gradient for predicting bathymetry is explored in this paper. This technique is based on the fact that the observed gravity gradient anomalies result primarily from local mass concentrations on the ocean floor, and that mass compensation by the oceanic crust has an insignificant effect on the gravity gradients, and can be neglected. The resulting bathymetry prediction therefore is independent of isostatic modeling assumptions, allowing it to be used as a constraint on models of lithospheric compensation and for other geodetic and geophysical applications.     

ICESat-2激光卫星与光学遥感影像融合水深测量

针对浅海测深的数据特点和应用需求,以我国南海甘泉岛为例,研究了利用ICESat-2(Ice,Cloud,and Land Elevation Satellite-2)激光卫星数据和光学遥感影像开展主被动融合水深测量的方法。首先通过信号点提取、水面/水底识别、水下点折射改正等步骤处理ICESat-2数据,获得水深值,随后以激光点作为控制,计算光学水深反演模型参数,最后由点及面地获取大范围高精度水深。实验表明,甘泉岛区域主被动融合测深中误差优于1.30m,基于激光卫星数据的主被动融合测深方法能够为浅水水深测量提供新手段。     

Influence of bathymetric fluctuations on coastal storm surge

Natural events constantly alter nearshore bathymetric properties. Hurricanes particularly affect bathymetry as they pass over a body of water. To compute an accurate forecast or recreate a hurricane's effects through hindcasting techniques, an operational bathymetry data set must be known in advance. However, obtaining and maintaining current and accurate bathymetric data can be costly and difficult to manage. In this paper we examine the extent to which variations in nearshore bathymetry affect the storm surge at the coast. A common question for wave and surge modeling is, “how good is the bathymetric data?” If we can allow for a range of fluctuations in the bathymetry without significantly adjusting the results of the surge predictions, we can potentially save months of field work and millions of dollars. A one-dimensional (1D) analytical solution for waves and water level is developed for initial testing. In the 1D case we find that as long as the amplitudes of the bathymetric fluctuations are less than 60% of the original depth, the surge at the coast is within ± 10% of the surge generated on the initial bottom slope. If the fluctuation produces a hole, a deepening of the local bathymetry, within 80% of the local water depth, the coastal storm surge calculated is still within 10% of the unperturbed value computed for bottom slopes shallower than 1:20. In addition, we find there is an optimum distance offshore for each sloped profile that corresponds to a depth between 25 and 40 m, beyond which the effects of bathymetric fluctuations begin to decrease. A coupled 2D modeling system is implemented to test our hypothesis along a realistic coastline. After selecting three study sites, we vary the bathymetry at the selected locations by ± 20%. Consistent with the 1D tests, the storm surge at the shoreline varies by less than 5%.     

Bathymetry Data Correction Using Global Optimization Method

The bathymetry data of marine bodies have been collected over a century, and the collected data have a wide range of resolution and accuracy. Acquisition of bathymetry data is very costly and time-consuming. One can use the old, low-quality bathymetry data to fill the gap in high-quality, recently acquired bathymetry data after correcting the old data to improve its quality so that it is comparable to the high-quality data. The old data correction can be treated as a nonlinear inverse problem. Simulated annealing (SA) global optimization method was used here in solving this problem. The two sets of data that were used are project survey (PS) and Vietnamese Navy Chart (VNC) data. The PS data were collected in 2000 in an offshore survey from the Vietnam coast in the South China Sea (SCS). The VNC data were obtained by digitizing VNC that was published in 1981. Inverse distance weighted (IDW) interpolation method was used for forward modeling. Weperformed the SA algorithm run starting at a high "temperature," then lowering the "temperature" gradually up to the "critical temperature" and then staying there for the rest of the run. The best model chosen by the algorithm showed an improvement of 63% from the original model. We then constructed a digital bathymetry model (DBM) of the study area with the combined corrected VNC and the PS data.     

Bathymetry Data Correction Using Global Optimization Method

The bathymetry data of marine bodies have been collected over a century, and the collected data have a wide range of resolution and accuracy. Acquisition of bathymetry data is very costly and time-consuming. One can use the old, low-quality bathymetry data to fill the gap in high-quality, recently acquired bathymetry data after correcting the old data to improve its quality so that it is comparable to the high-quality data. The old data correction can be treated as a nonlinear inverse problem. Simulated annealing (SA) global optimization method was used here in solving this problem. The two sets of data that were used are project survey (PS) and Vietnamese Navy Chart (VNC) data. The PS data were collected in 2000 in an offshore survey from the Vietnam coast in the South China Sea (SCS). The VNC data were obtained by digitizing VNC that was published in 1981. Inverse distance weighted (IDW) interpolation method was used for forward modeling. Weperformed the SA algorithm run starting at a high "temperature," then lowering the "temperature" gradually up to the "critical temperature" and then staying there for the rest of the run. The best model chosen by the algorithm showed an improvement of 63% from the original model. We then constructed a digital bathymetry model (DBM) of the study area with the combined corrected VNC and the PS data.     

Gas-hydrate stability thickness map along the Indian continental margin

The gas-hydrate stability thickness (GHST) map along the Indian continental margin is prepared from available bathymetry, sea-bottom temperature and geothermal gradient data. The bottom-simulating reflector (BSR) often marks the base of gas-hydrate stability zone. The prior information about the stability thickness in a particular area will help in identifying BSR on seismic data. The map is also useful to the exploration scientists to set a depth window within which proxies for gas-hydrate can be looked into. A GHST map was initially prepared in 1998 based on the-then available data. A lot of new data has been generated by various organizations under the Indian National Gas Hydrate Programs for the advancement of exploration and exploitation activities. By incorporating the new data from the published and available documents, we have modified the GHST map along the Indian margin. Besides filling the data gap, the new map shows the gas-hydrate stability zone in the Andaman offshore. In addition, we show maps of sea-bottom temperature, sediment thickness, geothermal gradient and heat flow to provide a bird’s eye view of these parameters along the continental margin of India.     

机载LiDAR测深在海洋测绘中应用的可行性分析

针对海洋测绘中遇到的环境恶劣、任务繁重等情况,分析了机载Li DAR测深应用到海洋测绘中的可行性。介绍机载Li DAR测深原理、优势、典型设备,针对系统最大穿透深度、最浅探测深度、测点密度及测深精度等关键技术指标进行分析。最终得出结论,机载Li DAR测深符合海洋测绘的需求,应用到海洋测绘中是切实可行的。     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

Abstract

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^? ≤ ? ≤ 35^? and 15^? ≤ λ ≤ 25^?. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8–20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).     

Automatic Sounding Generalization in Nautical Chart Considering Bathymetry Complexity Variations

Sounding acts as the main feature in a digital nautical chart as it describes the concerned marine topography for the safety of navigation. Unlike the geometry-oriented selection of point feature, the generalization of soundings for chart compiling is expected to be context-oriented, which means bathymetry complexity variations across the study region should be preserved in the sounding selection process. However, such variations are not explicitly accessible to automated systems. This paper proposes an approach that effectively analyzes and measures bathymetry complexity from sounding data, with a focus on topography variations among different regions. The presented approach first divides the exploring region into several subregions, by adopting techniques of computational geometry and graph theory. Then, the approach quantitatively measures the bathymetry complexity of the subregions from grid-based digital terrain model. Finally, a composite bathymetry complexity index integrating aspects of steepness and depth variation is developed to guide the operation of sounding selection in different subregions. Generally, when seafloor is rugged with steep slopes, the number of soundings is high. While in flatter areas, a smaller amount of soundings is retained. The potential of our approach is demonstrated by an application to a real data set.     

两种测高重力异常反演海底地形方法比较

针对基于测高重力异常反演海底地形理论众多、选取标准无法确定的情况,利用中国南海海域内的测高重力异常和船测水深数据研究比较了重力地质法(GGM)和SmithSandwell (SAS)法两种精度高、计算速度相对较快的海底地形反演理论。其中,GGM方法的密度差异常数Δρ由向下延拓技术确定为2.15 g·cm-3,SAS方法采用移去-恢复技术得到反演波段内重力异常和水深数据。结果表明:测线分布条件一定时,水深多在-1 000 m左右或反演区域岛礁、海山等复杂海底地形较多时选取SAS方法,水深主要在-3 000 m以深的区域或海底地形复杂程度不高时选取GGM方法则能获取更好的效果,其效果最优处与船测水深在检核点处的差值最优平均值能达-0.61 m,标准差可达14.67 m。     

基于SPOT-6遥感影像的近海水深反演

水深反演对于浅海地形调查、海岸带保护和开发具有重要的意义。本文选取南海东锣岛海域为研究区,采用SPOT-6 6 m高分辨率数据,基于SPOT-6多光谱数据和水深实测值分别建立蓝-绿、蓝-红等6个波段比值模型。结果显示,绿-红波段的比值模型精度最高,R²值达到0.706 4,这对于反演热带海洋地区水深具有一定借鉴意义。基于该模型的反演结果对不同水深范围内的平均相对误差进行比较,结果表明:0~5 m水深范围反演误差较高,主要由于该范围水体情况较复杂,5～10 m 水深范围内平均相对误差最小为13.62%,而随着水深的增加,反演的误差增大,分析误差的主要来源是海水中的悬浮颗粒物、黄色物质、叶绿素质量浓度等方面的影响。     

High-resolution seafloor mapping using the DSL-120 sonar system: Quantitative assessment of sidescan and phase-bathymetry data from the Lucky Strike segment of the Mid-Atlantic Ridge

High-resolution, side-looking sonar data collected near the seafloor (100 m altitude) provide important structural and topographic information for defining the geological history and current tectonic framework of seafloor terrains. DSL-120 kHz sonar data collected in the rift valley of the Lucky Strike segment of the Mid-Atlantic Ridge near 37° N provide the ability to quantitatively assess the effective resolution limits of both the sidescan imagery and the computed phase-bathymetry of this sonar system. While the theoretical, vertical and horizontal pixel resolutions of the DSL-120 system are <1 m, statistical analysis of DSL-120 sonar data collected from the Lucky Strike segment indicates that the effective spatial resolution of features is 1–2 m for sidescan imagery and 4 m for phase-bathymetry in the seafloor terrain of the Mid-Atlantic Ridge rift valley. Comparison of multibeam bathymetry data collected at the sea-surface with deep-tow DSL-120 bathymetry indicates that depth differences are on the order of the resolution of the multibeam system (10–30 m). Much of this residual can be accounted for by navigational mismatches and the higher resolving ability of the DSL-120 data, which has a bathymetric footprint on the seafloor that is 20 times smaller than that of hull-mounted multibeam at these seafloor depths (2000 m). Comparison of DSL-120 bathymetry with itself on crossing lines indicates that residual depth values are ±20 m, with much of that variation being accounted for by navigational errors. A DSL-120 survey conducted in 1998 on the Juan de Fuca Ridge with better navigation and less complex seafloor terrain had residual depth values half those of the Lucky Strike survey. The quality of the bathymetry data varies as a function of position within the swath, with poorer data directly beneath the tow vehicle and also towards the swath edges.Variations in sidescan amplitude observed across the rift valley and on Lucky Strike Seamount correlate well with changes in seafloor roughness caused by transitions from sedimented seafloor to bare rock outcrops. Distinct changes in sonar backscatter amplitude were also observed between areas covered with hydrothermal pavement that grade into lava flows and the collapsed surface of the lava lake in the summit depression of Lucky Strike Seamount. Small features on the seafloor, including volcanic constructional features (e.g., small cones, haystacks, fissures and collapse features) and hydrothermal vent chimneys or mounds taller than 2 m and greater than 9 m² in surface area, can easily be resolved and mapped using this system. These features at Lucky Strike have been confirmed visually using the submersible Alvin, the remotely operated vehicle Jason, and the towed optical/acoustic mapping system Argo II.     

基于水深分段选择因子的多光谱影像反演水深

岛礁周边海底地形精密测量是海洋测绘的难点问题之一,利用多光谱卫星能精确测定岛礁周边水深地形,是传统岛礁水深测量有效的补充方式之一。本文以蜈支洲岛周围水域为研究区域,提出了基于BP神经网络水深分段选取反演因子的方法,利用WorldView-2数据和多波束实测数据建立多模型,整合预测的各段水深获得海底地形。实验表明,本文方法能充分利用水体的光谱特性,与单一BP神经网络模型结果相比,MRE降低了6.4%,RMSE降低了1.2 m,是一种可行的多光谱水深反演的方法     

海底地形测量成果的质量检核评估(三):自容式压力验潮仪数据的精密处理

水位控制是海底地形测量资料处理的重要技术环节之一,潮汐观测数据质量控制则是提高水位改正精度的关键。针对当前自容式压力验潮仪使用中经常会遇到的潮汐观测数据缺失、波浪扰动影响等技术难题以及海水密度精密改正、基准面传递确定等技术需求,开展了相关理论方法研究并开发了较实用的潮汐数据处理模块。结果表明,经多个环节的精密处理后潮汐观测数据精度提高明显,适用于海底地形测量水位改正及其他海洋工程应用。     

Three-Dimensional Geographically Weighted Inverse Regression (3GWR) Model for Satellite Derived Bathymetry Using Sentinel-2 Observations

Current trends of development of satellite derived bathymetry (SDB) models rely on applying calibration techniques including analytical approaches, neuro-fuzzy systems, regression optimization and others. In most of the cases, the SDB models are calibrated and verified for test sites, that provide favorable conditions for the remote derivation of bathymetry such as high water clarity, homogenous bottom type, low amount of sediment in the water and other factors. In this paper, a novel 3-dimensional geographical weighted regression (3GWR) SDB technique is presented, it binds together methods already presented in other studies, the geographically weighted local regression (GWR) model, with depth dependent inverse optimization. The proposed SDB model was calibrated and verified on a relatively difficult test site of the South Baltic near-shore areas with the use of multispectral observations acquired by a recently launched Sentinel-2 satellite observation system. By conducted experiments, it was shown that the proposed SDB model is capable of obtaining satisfactory results of RMSE ranging from 0.88 to 1.23[m] depending on the observation and can derive bathymetry for depths up to 12m. It was also shown, that the proposed approach may be used operationally, for instance, in the continuous assessment of temporal bathymetry changes, for areas important in the context of ensuring local maritime safety.     

一种新型无人船海洋水深测量技术

在海洋测深中,由于波浪和潮汐的影响,调查船或无人船所测量的瞬时水深不能直接作为海图水深。本文提出了一种新型的无人船海洋水深测量技术,以评估搭载RTK和单波束测深仪的无人船用于海洋水深测量的潜能。首先,使用无人船所搭载RTK的厘米级精度高程数据,通过低通滤波剔除波浪信息,而获得海平面高程。然后,基于潮汐表和无人船海平面高程,构建了一种参考椭球面和海图的两个基准面之差的获取方法;在常规的海洋调查中,该基准面差通常需要由长期的验潮获得。最后,利用海图基准和无人船测量的瞬时水深的转换关系,计算出海图水深。在海南省蜈支洲岛周边海域,利用自研发的无人船“USBV”开展了相关海上实验,以验证所提出的技术方法。实验结果验证了该无人船海洋水深测量技术。     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^ˆ ≤ ϕ ≤ 35^ˆ and 15^ˆ ≤ λ ≤ 25^ˆ. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8-20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).