RTM Gravity Forward-Modeling Using Topography/Bathymetry Data to Improve High-Degree Global Geopotential Models in the Coastal Zone

We apply the residual terrain modeling (RTM) technique for gravity forward-modeling to successfully improve high-resolution global gravity fields at short spatial scales in coastal zones. The RTM scheme is combined with the concept of rock-equivalent topography, allowing to use a single uniform constant mass-density in the RTM forward-modeling, both at land and sea. SRTM30_PLUS bathymetry is merged with higher-resolution SRTM V4.1 land topography, and expanded into spherical harmonics to degree 2160, yielding a new and consistent high-degree RTM reference surface. The forward-modeling performance is demonstrated in coastal zones of Greece and Canada using ground-truth vertical deflections, gravity from land and shipborne gravimetry, and geoid heights from GPS/leveling, with improvements originating from bathymetry clearly identified. We demonstrate that the SRTM30_PLUS bathymetry carries information on gravity field structures at spatial scales less than 5 arc minutes, which can be used to augment EGM2008 in (rugged) coastal zones, both over land and marine areas. This may be of value (i) to partially reduce the signal omission error in EGM2008/GOCE-based height transfer in areas devoid of dense gravity data, (ii) to fill the gap between land gravity and shipborne gravity along rugged coastlines, and (iii) for the development of next-generation altimetric gravity fields.     

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.     

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σ).     

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σ).     

The global distribution of seamounts based on 30 arc seconds bathymetry data

Seamounts and knolls are ‘undersea mountains’, the former rising more than 1000 m from the seafloor. These features provide important habitats for aquatic predators, demersal deep-sea fish and benthic invertebrates. However most seamounts have not been surveyed and their numbers and locations are not well known. Previous efforts to locate and quantify seamounts have used relatively coarse bathymetry grids. Here we use global bathymetric data at 30 arc-sec resolution to identify seamounts and knolls. We identify 33,452 seamounts and 138,412 knolls, representing the largest global set of identified seamounts and knolls to date. We compare estimated seamount numbers, locations, and depths with validation sets of seamount data from New Zealand and Azores. This comparison indicates the method we apply finds 94% of seamounts, but may overestimate seamount numbers along ridges and in areas where faulting and seafloor spreading creates highly complex topography. The seamounts and knolls identified herein are significantly geographically biased towards areas surveyed with ship-based soundings. As only 6.5% of the ocean floor has been surveyed with soundings it is likely that new seamounts will be uncovered as surveying improves. Seamount habitats constitute approximately 4.7% of the ocean floor, whilst knolls cover 16.3%. Regional distribution of these features is examined, and we find a disproportionate number of productive knolls, with a summit depth of <1.5 km, located in the Southern Ocean. Less than 2% of seamounts are within marine protected areas and the majority of these are located within exclusive economic zones with few on the High Seas. The database of seamounts and knolls resulting from this study will be a useful resource for researchers and conservation planners.     

Predicting Global Seafloor Topography Using Multi-Source Data

A new one-minute global seafloor topography model was derived from vertical gravity gradient anomalies (VGG), altimetric gravity anomalies, and ship soundings. Ship soundings are used to constrain seafloor topography at wavelengths longer than 200 km and to calibrate the topography to VGG (or gravity) ratios at short wavelengths area by area. VGG ratios are used to predict seafloor topography for wavelength bands of 100–200 km and to suppress the effect of crust isostasy. Gravity anomalies are used to recover seafloor topography at wavelengths shorter than 100 km. The data processing procedure is described in detail in this paper. The accuracy of the model is evaluated using ship soundings and existing models, including General Bathymetric Charts of the Oceans (GEBCO), DTU10, ETOPO1, and SIO V15.1. The results show that, in the discussed regions, the accuracy of the model is better than ETOPO1, GEBCO, and DTU10. Additionally, the model is comparable with V15.1, which is generally believed to have the highest accuracy. In the north-central Pacific Ocean, the accuracy of the model increased by approximately 29.5% compared with the V15.1 model. This indicates that a more accurate seafloor topography model can be formed by combining gravity anomalies, VGG, and ship soundings.     

A new bathymetric model for the central Fram Strait

Based on data from R/V Polarstern multibeam sonar surveys between 1984 and 1997 high resolution bathymetry has been generated for the central Fram Strait. The area insonified covers approx. 36,500 km² between 78–80°N and 0–7.5°E allowing the creation of a Digital Terrain Model (DTM) with 100 m grid spacing. The DTM was utilized for contouring and generation of a new series of bathymetric charts (AWI Bathymetric Charts of the Fram Strait, AWI BCFS) at a scale of 1:100,000. The paper starts with a brief introduction to the regional setting of the study area comprising information on the local links between bathymetry, sea ice transport and water mass exchange. The bathymetric feature names used in this article and how they were chosen is outlined. Next, the input data and processing applied are described. Thereafter the newly created grid and contour data are put into context with existing data sets. Finally the main bathymetric features of the area are characterized and the generated data products available for public disposal are specified.     

Swath mapping of the seafloor and its application to deep-bottom fisheries in New Caledonia

A bottom longline fishery operated in the EEZ of New Caledonia from 1988 to 1991. Fishing focused on five seamounts with summits at depths ranging from 500 to 750 m. The target species was alfonsino, Beryx splendens. As the soundings available from marine charts were not detailed enough, the fishing masters had to make their own charts in order to set the gear in the right location. A series of 11 scientific cruises devoted to a research program on alfonsino was launched in late 1991. During the first scientific cruises, several days were spent mapping the seabed to improve the existing knowledge of the topography of three seamounts by coupling the EDO echosounder depth measurements to the GPS positions. As this procedure is slow, it is applicable to limited areas or, if a wider grid is used, to large zones to locate major structures such as ridges and seamounts. The emergence of the multibeam echosounder has greatly improved seabed mapping performance. This tool seems to fit the requirement of exploratory deep-bottom fishing particularly well, as it covers large areas while providing details of the bottom. The EM 12 was used to obtain bathymetry and imagery of the SE portion of the EEZ of New Caledonia, covering an approximate area of 70,000 km² in two weeks. For depths ranging from 500 to 4,500 m, the results were impressive. They confirmed the known major features and provided greater topographical detail, revealing fine unknown structures. They also provided information on the type of substratum, information which might influence the fishing strategy. Finally, they made it possible to obtain an accurate estimate of the exploitable area which in turn led to new stock estimates.     

ONWET: A Simple Integrated Tool for Beach Morphology and Wave Dynamics Analysis

The analysis of beach profile combined with wave data provides a clear understanding of beach topography and physical behavior of the coast. Numerous commercial and noncommercial computer programs and add-on tools are available to analyze beach profile and wave data. But, most of them are specifically designed for particular coastal application. Moreover, searching of open source software for diverse coastal application is a tedious process. To make it easily accessible, we have developed an integrated Matlab GUI for beach morphology and wave analysis. The tool “ONWET (OceaN WavE Tool)” mainly includes wave processes, beach morphology, spatiotemporal variation in beach profile, wave refraction, and sediment transport rate. The tool has automated inbuilt functions and is designed to accommodate importing/exporting of data (results both in data and graphical format). It includes user friendly options and can be freely download from our institutional website. The tool is very useful for researchers, scientists, engineers, and environmentalists who are working in the field of coastal studies and beach nourishment.     

多波束数据的海底数字地形模型构建

提出大批量多格式原始多波束数据的DTM构建方法，以满足大区域、大比例尺海底地形制图的需要。在多种格式原始多波束数据接口和系统内部标准数据结构的基础上，通过对数据文件、数据种类和数据运算量的有效组织管理，实施边读入、边权重配赋的网格插值，分析了高斯、指数和距离平方反比权重函数的适用性，及最小值、最大值和平均值的实用性。在权重配赋网格插值基础上，提出分形fBM和张力样条配合使用的方程式插值方法，保证DTM数据的有效外延和地形分辨。整套算法效率高，并能有效保证DTM的精度，对存储在外部介质的数据遍历一次即可完成网格插值。     

多光谱遥感水深反演的数据匹配方法研究

在多光谱遥感水深反演中,数据匹配是建立统计相关模型的基础。采用两种多光谱遥感影像数据,研究了像元找水深点的顺序匹配法和水深点找像元的顺序匹配、二分匹配和像元坐标法等数据匹配方法,并开展相应的实验验证。结果表明,4种匹配方法一定程度上能满足当前多光谱遥感水深反演的数据匹配需要,其中,水深点找像元的像元坐标法耗时最短,匹配效率最高。     

A generic swath‐mapping data format

There is a pressing need for standardization of data derived from bathy‐metric swath‐mapping systems. Currently several dozen multibeam and sidescan sonar data formats exist within the oceanographic community, and more can be expected as new systems are developed. Without some standardization of swath‐mapping data formats, the capability for use and integration of data from different systems will be severely compromised.

This paper presents a strategy for organizing swath bathymetry data in a logical modular fashion that will allow data from all current swath bathymetric sonar systems to be stored and accessed in a common fashion. We have chosen the approach of defining compact efficient modules for each logically independent portion of a data record and storing it in a manner that is portable between diverse computer architectures and operating systems. This approach is extensible to accommodate new types of data. Although specifically developed for swath bathymetry, this format is also capable of supporting digital sidescan data and other types of swath data.     

An Evaluation of Publicly Available Global Bathymetry Grids

We evaluate the strengths and weaknesses of six publicly available global bathymetry grids: DBDB2 (Digital Bathymetric Data Base; an ongoing project of the Naval Research Laboratory), ETOPO2 (Earth Topography; National Geophysical Data Center, 2001, ETOPO2 Global 2’ Elevations [CD-ROM]. Boulder, Colorado, USA: U.S. Department of Commerce, National Oceanic and Atmospheric Administration), GEBCO (General Bathymetric Charts of the Oceans; British Oceanographic Data Centre, 2003, Centenary Edition of the GEBCO Digital Atlas [CD-ROM] Published on behalf of the Intergovernmental Oceanographic Commission and the International Hydrographic Organization Liverpool, UK), GINA (Geographic Information Network of Alaska; Lindquist et al., 2004), Smith and Sandwell (1997), and S2004 (Smith, unpublished). The Smith and Sandwell grid, derived from satellite altimetry and ship data combined, provides high resolution mapping of the seafloor, even in remote regions. DBDB2, ETOPO2, GINA, and S2004 merge additional datasets with the Smith and Sandwell grid; but moving from a pixel to grid registration attenuates short wavelengths (<20 km) in the ETOPO2 and DBDB2 solutions. Short wavelengths in the GINA grid are also attenuated, but the cause is not known. ETOPO2 anomalies are offset to the northeast, due to a misregistration in both latitude and longitude. The GEBCO grid is interpolated from 500 m contours that were digitized from paper charts at 1:10 million scale, so it is artificially smooth; yet new efforts have captured additional information from shallow water contours on navigational charts. The S2004 grid merges the Smith and Sandwell grid with GEBCO over shallow depths and polar regions, and so is intended to capture the best of both products. Our evaluation makes the choice of which bathymetry grid to use a more informed one. The U.S. government right to retain a non-exclusive royalty-free license in and to any copyright is acknowledged.     

测线间非采样点水深值的平滑效应研究

为了表述和评价水深原图及海图水深数据对海底地形的表征能力,引入了水深插值平滑效应新概念,分析了航海图海道测量及海图编辑综合的水深数据选取特性,给出了测线间水深插值平滑效应的评价指标和方法,并进行了实例统计计算,计算结果表明了方法的可行性和有效性。     

Robust Automatic Reduction of Multibeam Bathymetric Data Based on M-estimators

Multibeam echosounders have commonly been employed for a wide range of applications including offshore survey, navigation, hydrogeology, and oceanography. Because the tremendous volume of the bathymetric data is demanding for some purposes and requires significant storage space, the data reduction plays a prominent role in practice. Additionally, the multibeam soundings are inevitably contaminated with sporadic outliers, and as such, the data cleaning can be challenging especially in shallow waters. We present a speedily robust method for reliably reducing the volume of the bathymetric data within grid cells. In this respect, robust M-estimators are recursively applied to the data in a patch-wise manner to alleviate the undesirable effects of the outlying observations. Accordingly, the reduced bathymetry is automatically made unaffected by the possible outliers once their equivalent weights have been downweighted. The performance of the presented method has been demonstrated by synthetic datasets and an experimental dataset collected by an ATLAS FS 20/100 kHz shallow-water multibeam echosounder in the offshore waters of Kish wharf. The reliability, efficiency, and capability of the proposed method have been verified, which makes it quite possible to meet the IHO requirements for special-order seafloor mapping.     

Numerical simulation of a low-lying barrier island's morphological response to Hurricane Katrina

Tropical cyclones that enter or form in the Gulf of Mexico generate storm surge and large waves that impact low-lying coastlines along the Gulf Coast. The Chandeleur Islands, located 161 km east of New Orleans, Louisiana, have endured numerous hurricanes that have passed nearby. Hurricane Katrina (landfall near Waveland MS, 29 Aug 2005) caused dramatic changes to the island elevation and shape. In this paper the predictability of hurricane-induced barrier island erosion and accretion is evaluated using a coupled hydrodynamic and morphodynamic model known as XBeach. Pre- and post-storm island topography was surveyed with an airborne lidar system. Numerical simulations utilized realistic surge and wave conditions determined from larger-scale hydrodynamic models. Simulations included model sensitivity tests with varying grid size and temporal resolutions. Model-predicted bathymetry/topography and post-storm survey data both showed similar patterns of island erosion, such as increased dissection by channels. However, the model under predicted the magnitude of erosion. Potential causes for under prediction include (1) errors in the initial conditions (the initial bathymetry/topography was measured three years prior to Katrina), (2) errors in the forcing conditions (a result of our omission of storms prior to Katrina and/or errors in Katrina storm conditions), and/or (3) physical processes that were omitted from the model (e.g., inclusion of sediment variations and bio-physical processes).     

基于回波时间的多波束测深与声纳数据匹配方法

以多波束精确的水深数据为参照源,采用原始回波时间对多波束测深数据与其同源声纳数据进行匹配,从而获得高精度和高分辨率的海底影像数据,并避免了传统声纳图像处理过程中斜距改正所带来的几何形变。匹配结果采用光照图输出,并与三维水深图、原始声纳图像和CARIS处理后的声纳图像进行比较分析。该方法有效地提高了多波束数据的利用率,增强了对海底地形的探测分辨率。     

基于CUBE算法的多波束测深数据自动处理研究

对CUBE算法自动处理多波束测深数据的模型建立、格网节点的多重估计和最优估值选取准则进行了详细介绍,深入分析了多重估计的实用性,并通过实测数据对该算法进行实现.利用了抗差Kalman滤波改进CUBE算法.通过模拟数据对改进的CUBE算法进行实验,验证了算法改进的必要性.     

海底地形格网数据精度指标要求分析

海底地形不易进行实地探测及核实,海底地形格网数据产品精度指标的确定不具备实地核实的条件,为了确定海底地形格网数据精度指标,通过大量的数据实验,分析得出单位面积内,一定比例的极值水深差值的平均值与海底地形之间有一定的线性关系,从而提出了通过单位面积水深插值来划分海底地形的一种量化方式,并依据实验数据,分析误差规律,给出了海底地形格网数据中误差的精度指标及方法。     

Processing and analysis of Simrad multibeam sonar data

The common approach to analysing data collected with multibeam and sidescan sonars is to visually interpret charts of contoured bathymetry and mosaics of seabed images. However, some of the information content is lost by processing the data into charts because this involves some averaging; the analysis might uncover more information if done on the data at an earlier stage in the processing. Motivated by this potential, I have created a software system which can be used to analyse data collected with Simrad EM1000 (shallow water) and EM12 (deep water) multibeam sonars, as well as to generate bathymetry contour charts and backscatter mosaics. The system includes data preprocessing, such as navigation filtering, depth filtering (removal of outlying values), and amplitude mapping using the multibeam bathymetry to correctly position image pixels across the swath. The data attributes that can be analysed include the orientation and slope of the seafloor, and the mean signal strength for each sounding. To determine bathymetry attributes such as slope, the soundings across a number of beams and across a series of pings are grouped and a least-squares plane fitted to them. Bathymetric curvature is obtained by detrending the grouped data using the least-squares plane and fitting a paraboloid to the residuals. The magnitudes and signs of the paraboloid's coefficients reveal depressions and hills and their orientations. Furthermore, the seafloor geology can be classified using a simple combination of these attributes. For example, flat-lying sediments can be classified where the backscatter, slope and curvature fall below specified values.