Computation of calibration parameters for multibeam echo sounders using the least squares method

A new method for field calibration of multibeam echo sounders is presented. The method applies least-squares adjustment computation of static offsets in parameters such as roll, pitch, heading and the horizontal shift of the transducer relative to the reference point of the positioning system. All the parameters are derived simultaneously. An investigation demonstrates the precision and the flexibility of the method. Flexibility is here related to the design of calibration lines, i.e., their number, their direction and length. The strength of the method is its ability to utilize all the depth measurements in an optimal way. The variance of the parameters can be derived directly from least-squares adjustment, but limited knowledge of the stochastic properties of the observations can lead to optimistic precision estimates. An alternative method to compute quality estimates is proposed and demonstrated.     

混浊水域多波束反向散射强度传播损失改正

为了改善现有的反向散射强度传播损失改正方法在混浊水域适应性较弱的问题,提出了混浊水域声传播损失计算模型。首先利用不同深度的海洋环境参数构建声波吸收系数剖面,然后基于声速,沿波束传播路径,对每个波束分层计算传播损失。实验分析表明,本方法传播损失改正效果最佳,其Spearman等级相关系数绝对值仅为0.04,远小于传统模型和TVG改正,大大减弱了混浊水域情况下多波束回波强度与传播距离的相关性,有效改善了多波束声呐图像的质量。     

Analysis of seafloor backscatter strength dependence on the survey azimuth using multibeam echosounder data

The sediment backscatter strength measured by multibeam echosounders is a key feature for seafloor mapping either qualitative (image mosaics) or quantitative (extraction of classifying features). An important phenomenon, often underestimated, is the dependence of the backscatter level on the azimuth angle imposed by the survey line directions: strong level differences at varying azimuth can be observed in case of organized roughness of the seabed, usually caused by tide currents over sandy sediments. This paper presents a number of experimental results obtained from shallow-water cruises using a 300-kHz multibeam echosounder and specially dedicated to the study of this azimuthal effect, with a specific configuration of the survey strategy involving a systematic coverage of reference areas following “compass rose” patterns. The results show for some areas a very strong dependence of the backscatter level, up to about 10-dB differences at intermediate oblique angles, although the presence of these ripples cannot be observed directly—neither from the bathymetry data nor from the sonar image, due to the insufficient resolution capability of the sonar. An elementary modeling of backscattering from rippled interfaces explains and comforts these observations. The consequences of this backscatter dependence upon survey azimuth on the current strategies of backscatter data acquisition and exploitation are discussed.     

两种多波束反向散射强度数据归一化方法的建立与分析

多波束反向散射强度数据应用广泛,但由于受到角度响应的影响,导致生成的多波束声呐图像质量偏低,且现有角度响应改正方法在复杂海底底质环境下适应性较差。为此本文对散射强度进行分析,给出了两种多波束反向散射强度数据归一化方法,分别为基于高斯拟合以及角度响应的散射强度改正方法,前者主要是基于散射强度的变化规律进行改正,而后者则是基于声波的散射机理进行改正。实验结果表明两种方法较传统改正方法精度均有约30%的提升,并且角度响应方法较高斯拟合方法改正精度更高,但计算效率有所下降。以上实验验证了两种方法的有效性,实现了散射强度数据的归一化,提升了多波束声呐图像的质量。     

Processing of high-frequency multibeam echo sounder data for seafloor characterization

Processing simultaneous bathymetry and backscatter data, multibeam echosounders (MBESs) show promising abilities for remote seafloor characterization. High-frequency MBESs provide a good horizontal resolution, making it possible to distinguish fine details at the water-seafloor interface. However, in order to accurately measure the seafloor influence on the backscattered energy, the recorded sonar data must first be processed and cleared of various artifacts generated by the sonar system itself. Such a preprocessing correction procedure along with the assessment of its validity limits is presented and applied to a 95-kHz MBES (Simrad EM 1000) data set. Beam pattern effects, uneven array sensitivities, and inaccurate normalization of the ensonified area are removed to make possible further quantitative analysis of the corrected backscatter images. Unlike low-frequency data where the average backscattered energy proves to be the only relevant feature for discriminating the nature of the seafloor, high-frequency MBES backscatter images exhibit visible texture patterns. This additional information involves different statistical distributions of the backscattered amplitudes obtained from various seafloor types. Non-Rayleigh statistics such as K-distributions are shown to fit correctly the skewed distributions of experimental high-frequency data. Apart from the effect of the seafloor micro-roughness, a statistical model makes clear a correlation between the amplitude statistical distributions and the signal incidence angle made available by MBES bathymetric abilities. Moreover, the model enhances the effect of the first derivative of the seafloor backscattering strength upon statistical distributions near the nadir and at high incidence angles. The whole correction and analysis process is finally applied to a Simrad EM 1000 data set.     

Toward a standard line for use in multibeam echo sounder calibration

A procedure is suggested in which a relative calibration for the intensity output of a multibeam echo sounder (MBES) can be performed. This procedure identifies a common survey line (i.e., a standard line), over which acoustic backscatter from the seafloor is collected with multiple MBES systems or by the same system multiple times. A location on the standard line which exhibits temporal stability in its seafloor backscatter response is used to bring the intensity output of the multiple MBES systems to a common reference. This relative calibration procedure has utility for MBES users wishing to generate an aggregate seafloor backscatter mosaic using multiple systems, revisiting an area to detect changes in substrate type, and comparing substrate types in the same general area but with different systems or different system settings. The calibration procedure is demonstrated using three different MBES systems over 3 different years in New Castle, NH, USA.     

An uncertainty model for deep ocean single beam and multibeam echo sounder data

Comparing single beam and multibeam echo sounder data where surveys overlap we find that: 95% of multibeam measurements are repeatable to within 0.47% of depth; older single beam data can be at least as accurate as multibeam; single beam and multibeam profiles show excellent agreement at full-wavelengths longer than 4 km; archival sounding errors are not Gaussian; 95% of archival soundings in the northwest Atlantic are accurate to within 1.6% of depth; the 95th percentile error is about five times greater in pre-1969 data than in post-1968 data; many of the largest errors are located over large seafloor slopes, where small navigation errors can lead to large depth errors. Our uncertainty model has the form σ ² = a ² + (bz)² + (cs)², where 2σ is approximately the 95th percentile error, z is the depth, s is the slope, and a, b, c are constants we determine separately for pre-1969 and post-1968 data.     

Control of the repeatability of high frequency multibeam echosounder backscatter by using natural reference areas

The increased use of backscatter measurements in time series for environmental monitoring necessitates the comparability of individual results. With the current lack of pre-calibrated multibeam echosounder systems for absolute backscatter measurement, a pragmatic solution is the use of natural reference areas for ensuring regular assessment of the backscatter measurement repeatability. This method mainly relies on the assumption of a sufficiently stable reference area regarding its backscatter signature. The aptitude of a natural area to provide a stable and uniform backscatter response must be carefully considered and demonstrated by a sufficiently long time-series of measurements. Furthermore, this approach requires a strict control of the acquisition and processing parameters. If all these conditions are met, stability check and relative calibration of a system are possible by comparison with the averaged backscatter values for the area. Based on a common multibeam echosounder and sampling campaign completed by available bathymetric and backscatter time series, the suitability as a backscatter reference area of three different candidates was evaluated. Two among them, Carré Renard and Kwinte, prove to be excellent choices, while the third one, Western Solent, lacks sufficient data over time, but remains a valuable candidate. The case studies and the available backscatter data on these areas prove the applicability of this method. The expansion of the number of commonly used reference areas and the growth of the number of multibeam echosounder controlled thereon could greatly contribute to the further development of quantitative applications based on multibeam echosounder backscatter measurements.     

Validation of automated supervised segmentation of multibeam backscatter data from the Chatham Rise,New Zealand

Using automated supervised segmentation of multibeam backscatter data to delineate seafloor substrates is a relatively novel technique. Low-frequency multibeam echosounders (MBES), such as the 12-kHz EM120, present particular difficulties since the signal can penetrate several metres into the seafloor, depending on substrate type. We present a case study illustrating how a non-targeted dataset may be used to derive information from multibeam backscatter data regarding distribution of substrate types. The results allow us to assess limitations associated with low frequency MBES where sub-bottom layering is present, and test the accuracy of automated supervised segmentation performed using SonarScope® software. This is done through comparison of predicted and observed substrate from backscatter facies-derived classes and substrate data, reinforced using quantitative statistical analysis based on a confusion matrix. We use sediment samples, video transects and sub-bottom profiles acquired on the Chatham Rise, east of New Zealand. Inferences on the substrate types are made using the Generic Seafloor Acoustic Backscatter (GSAB) model, and the extents of the backscatter classes are delineated by automated supervised segmentation. Correlating substrate data to backscatter classes revealed that backscatter amplitude may correspond to lithologies up to 4 m below the seafloor. Our results emphasise several issues related to substrate characterisation using backscatter classification, primarily because the GSAB model does not only relate to grain size and roughness properties of substrate, but also accounts for other parameters that influence backscatter. Better understanding these limitations allows us to derive first-order interpretations of sediment properties from automated supervised segmentation.     

Comparisons of sea surface height variability observed by pressure-recording inverted echo sounders and satellite altimetry in the Kuroshio Extension

Satellite-measured along-track and gridded sea surface height (SSH) anomaly products from AVISO are compared with in situ SSH anomaly measurements from an array of 43 pressure-recording inverted echo sounders (PIESs) in the Kuroshio Extension. PIESs measure bottom pressure (P _bot) and round-trip acoustic travel time from the sea floor to the sea surface (τ). The P _bot and τ measurements are used to estimate, respectively, the mass-loading and steric height variations in SSH anomaly. All comparisons are made after accurate removal of tidal components from all data. Overall good correlations are found between along-track and PIES-derived SSH anomalies with mean correlation coefficient of 0.97. Comparisons between the two measurements reveal that the mass-loading component estimated from P _bot is relatively small in this geographical region. It improves regression coefficients about 5?% and decreases mean root-mean-squared (rms) differences from 7.8 to 6.4?cm. The AVISO up-to-date gridded product, which merges all available satellite measurements of Jason-1, Envisat, Geosat Follow-On, and TOPEX/Poseidon interlaced, shows better correlations and smaller rms differences than the AVISO reference gridded product, which merges only Jason-1 and Envisat. Especially, the up-to-date gridded product reveals 6.8?cm rms improvement on average at sites away from Jason-1 ground tracks. Gridded products exhibit low correlation (0.75–0.9) with PIES-derived SSH in a subregion where the SSH fluctuations have relatively high energy at periods shorter than 20?days.     

Multisource multibeam backscatter data: developing a strategy for the production of benthic habitat maps using semi-automated seafloor classification methods

The establishment of multibeam echosounders (MBES) as a mainstream tool in ocean mapping has facilitated integrative approaches towards nautical charting, benthic habitat mapping, and seafloor geotechnical surveys. The bathymetric and backscatter information generated by MBES enables marine scientists to present highly accurate bathymetric data with a spatial resolution closely matching that of terrestrial mapping, and can generate customized thematic seafloor maps to meet multiple ocean management needs. However, when a variety of MBES systems are used, the creation of objective habitat maps can be hindered by the lack of backscatter calibration, due for example, to system-specific settings, yielding relative rather than absolute values. Here, we describe an approach using object-based image analysis to combine 4 non-overlapping and uncalibrated (backscatter) MBES coverages to form a seamless habitat map on St. Anns Bank (Atlantic Canada), a marine protected area hosting a diversity of benthic habitats. The benthoscape map was produced by analysing each coverage independently with supervised classification (k-nearest neighbor) of image-objects based on a common suite of 7 benthoscapes (determined with 4214 ground-truthing photographs at 61 stations, and characterized with backscatter, bathymetry, and bathymetric position index). Manual re-classification based on uncertainty in membership values to individual classes—especially at the boundaries between coverages—was used to build the final benthoscape map. Given the costs and scarcity of MBES surveys in offshore marine ecosystems—particularly in large ecosystems in need of adequate conservation strategies, such as in Canadian waters—developing approaches to synthesize multiple datasets to meet management needs is warranted.     

Morphology of pockmarks along the western continental margin of India: Employing multibeam bathymetry and backscatter data

This study addresses the morphology of pockmarks along the western continental margin of India using multibeam bathymetry and backscatter data. Here, for the first time we have utilized the application of ArcGIS (Geographical Information System) for understanding the morphology of pockmarks from the western continental margin of India. The pockmarks observed in water depths of 145–330 m are circular, elliptical or elongated in plan-view, with an average length and width of 157 (±72) m and 83 (±19) m respectively. The average pockmark relief and perimeter are 1.9 (±0.9) m and 412 (±181) m, respectively. The pockmarks have average areas and volumes of 10 759 m² and 15 315 m³ respectively. Spatial separation that coincides with 210 m isobath divides the pockmarks into two groups with differing distributions and morphologies. These pockmarks originated from seepages of biogenic or thermogenic gas or from pore fluids from deeper sources, migrated vertically along the faults. Besides a possible structural control, the pockmark morphologies are also affected by bottom currents and/or by submarine slumping. The average acoustic backscatter strength from pockmark centre is higher (−35 dB) than the average backscatter of the total area (−40 dB), which suggests their possible linkage to the precipitation of diagenitic minerals from biodegradation of seepage material.     

Evaluation of image-based multibeam sonar backscatter classification for benthic habitat discrimination and mapping at Stanton Banks,UK

In recent years, efforts have increased to develop quantitative, computer-directed methods for segmentation of multibeam (MBES) backscatter data. This study utilises MBES backscatter data acquired at Stanton Banks (UK) and subsequently processed through the QTC-Multiview software environment in a bid to evaluate the program's ability to perform unsupervised classification. Statistical comparison with ground-truth data (grab, stills and video) enabled cross validation of acoustic segmentation and biological assemblages observed at the site. 132 unspecified variables were extracted from user-specified rectangular patches of the backscatter image, reduced to three vectors by PCA, then clustered and classified by the software. Multivariate analyses of ground-truth data were conducted on 75 stills images and 51 grab samples. Video footage coincident with the stills was divided into 30 s segments and coded by dominant substrate and species. Cross tabulation determined the interrelationship between software classifications, multivariate analysis of the biological assemblages and coded video segments. Multiview optimally identified 19 classes using the automated clustering engine. These were revised to 6 habitats a posteriori, using combined analysis of ground-truth data and Multiview data products. These habitats broadly correspond to major physiographic provinces within the region. Multivariate statistical analysis reveals low levels of assemblage similarity (<35%) for samples occurring within Multiview classes, irrespective of the mode of acquisition. Coded video data is more spatially appropriate than the other methods of ground-truthing investigated, although it is less well suited to the extraction of truly quantitative data. Multivariate analysis indicates assemblages within physiographically distinct Multiview classes have a low degree of biological similarity, supporting the notion that abiotic proxies may be contraindicative of benthic assemblage variations. QTC-Multiview performs well as a mechanism for computer-assisted segmentation of MBES backscatter imagery into acoustic provinces; however a degree of caution is required prior to ascribing ecological significance to these classifications.     

Seafloor change detection using multibeam echosounder backscatter: case study on the Belgian part of the North Sea

To characterize seafloor substrate type, seabed mapping and particularly multibeam echosounding are increasingly used. Yet, the utilisation of repetitive MBES-borne backscatter surveys to monitor the environmental status of the seafloor remains limited. Often methodological frameworks are missing, and should comprise of a suite of change detection procedures, similarly to those developed in the terrestrial sciences. In this study, pre-, ensemble and post-classification approaches were tested on an eight km² study site within a Habitat Directive Area in the Belgian part of the North Sea. In this area, gravel beds with epifaunal assemblages were observed. Flourishing of the fauna is constrained by overtopping with sand or increased turbidity levels, which could result from anthropogenic activities. Monitoring of the gravel to sand ratio was hence put forward as an indicator of good environmental status. Seven acoustic surveys were undertaken from 2004 to 2015. The methods allowed quantifying temporal trends and patterns of change of the main substrate classes identified in the study area; namely fine to medium homogenous sand, medium sand with bioclastic detritus and medium to coarse sand with gravel. Results indicated that by considering the entire study area and the entire time series, the gravel to sand ratio fluctuated, but was overall stable. Nonetheless, when only the biodiversity hotspots were considered, net losses and a gradual trend, indicative of potential smothering, was captured by ensemble and post-classification approaches respectively. Additionally, a two-dimensional morphological analysis, based on the bathymetric data, suggested a loss of profile complexity from 2004 to 2015. Causal relationships with natural and anthropogenic stressors are yet to be established. The methodologies presented and discussed are repeatable and can be applied to broad-scale geographical extents given that broad-scale time series datasets become available.     

Combining pixel and object based image analysis of ultra-high resolution multibeam bathymetry and backscatter for habitat mapping in shallow marine waters

Habitat mapping data are increasingly being recognised for their importance in underpinning marine spatial planning. The ability to collect ultra-high resolution (cm) multibeam echosounder (MBES) data in shallow waters has facilitated understanding of the fine-scale distribution of benthic habitats in these areas that are often prone to human disturbance. Developing quantitative and objective approaches to integrate MBES data with ground observations for predictive modelling is essential for ensuring repeatability and providing confidence measures for habitat mapping products. Whilst supervised classification approaches are becoming more common, users are often faced with a decision whether to implement a pixel based (PB) or an object based (OB) image analysis approach, with often limited understanding of the potential influence of that decision on final map products and relative importance of data inputs to patterns observed. In this study, we apply an ensemble learning approach capable of integrating PB and OB Image Analysis from ultra-high resolution MBES bathymetry and backscatter data for mapping benthic habitats in Refuge Cove, a temperate coastal embayment in south-east Australia. We demonstrate the relative importance of PB and OB seafloor derivatives for the five broad benthic habitats that dominate the site. We found that OB and PB approaches performed well with differences in classification accuracy but not discernible statistically. However, a model incorporating elements of both approaches proved to be significantly more accurate than OB or PB methods alone and demonstrate the benefits of using MBES bathymetry and backscatter combined for class discrimination.     

基于多波束背向散射强度信号的海底表层沉积物粒度分类研究——以澳洲Joseph Bonaparte湾为例

近海海底地形探测与沉积物精确分类对涉海工程建设、生物栖息地反演以及海底资源勘查与开发具有重要的现实意义。以澳洲Joseph Bonaparte湾为例,利用多波束测深技术获取了该海湾约880 km2水域的水深数据与背向散射强度信号,结合同步采集的54个海底表层沉积物样品,通过随机决策树模型对该海域海底表层沉积物进行了分类研究。结果表明:（1）利用随机决策树模型分析该海域沉积物类型与背向散射强度的关系时,当模型内部参数设置:树的总数为200,最小分裂节点为2,每棵树的最大分裂级数为5时,可提高预测准确率;（2）该参数设置下,利用13°和37°入射角的背向散射强度预测该海域沉积物类型时,准确率最高,其值为83.3%,且在研究海域,砂质砾和砾质砂分布在背向散射强度较强的深槽或海沟等地区,而砾质泥质砂和含砾泥质砂主要分布在背向散射强度较弱的浅水海域。分析还发现,当水深数据作为预测海底表层沉积物类型的特征变量时,有可能降低最终预测结果的准确率。     

日本海西南海域现场观测和卫星高度计获取的海面高度距平的比较研究

压力传感逆式回声仪(pressure-sensor-equipped inverted echo sounders,PIES)可以用来测量海底压力和声波从海底到海面的传播时间。海底压力和声波传播时间分别被用来估计水体质量变化(正压)和比容变化(斜压)对海面高度距平的贡献。对由PIES在日本海西南海域现场观测数据得到的海面高度距平(PIES SLA)与卫星高度计海面高度距平(Sat SLA)进行了比较研究。利用相关分析法,对PIES SLA和沿轨T/P卫星、沿轨ERS-2卫星测得的海面高度距平(TP SLA、ERS-2 SLA)进行了比较;对PIES SLA和AVISO网格化海面高度距平进行了比较,估计可能的误差来源,并分析PIES SLA正压部分和斜压部分对SLA的贡献。比较发现,PIES SLA和Sat SLA的相关系数较高,且均方根误差较小,并且对特定区域和特定站点产生误差可能的原因进行了进一步的探讨。通过研究,有以下结论:(1)相对于湾流和黑潮地区,这一区域正压部分对海面高度的贡献相对较大;(2)如果再考虑斜压变化对海面高度的贡献,PIES SLA和Sat SLA相关系数会有所提升;(3)在高能区PIES SLA和Sat SLA相关系数较高,符合得相对比较好。总的来说,在日本海地区,PIES SLA和Sat SLA相关系数较高,具有较高的一致性,能为我国海洋二号(HY-2)等卫星高度计的校验提供一种可靠的方式。该研究对于PIES的研发和设计以及对于PIES的布放位置的选择都有一定的借鉴意义。     

Analysis of simrad EM12 multibeam bathymetry and acoustic backscatter data for seafloor mapping,exemplified at the Mid-Atlantic Ridge at 45° N

The EM12 multibeam echosounder can record acoustic backscatter information as well as high resolution bathymetry. The dataset presented, from the axis of the Mid-Atlantic Ridge at 45° N, was the first EM12 survey of a mid-ocean ridge. This paper presents methods for utilising the backscatter information. Data processing enables the production of a mosaic of acoustic backscatter, and visualisation techniques are investigated to provide initial qualitative views of the combined backscatter and bathymetry datasets. The co-registration of the backscatter and bathymetry data enables quantitative analysis of their relationships. Various sites of different geological type have been selected and their angular acoustic backscattering relationships estimated, including the effect on backscatter of incidence angle, its regional variability with bottom type and the influence of bottom slope. Incidence angles and bottom type are shown to affect backscatter to a similar degree, while slopes appear to contribute little. The geometry of hull-mounted systems, such as the EM12, is significantly different from that of conventional sidescan sonars, such as GLORIA, and the backscatter images from the two types differ in various respects. Because of the wide variations in incidence angle that are common with hull-mounted systems, and the importance of incidence angle in determining backscatter strength, it is vital to consider the effect of incidence angle during interpretation.     

Calibration techniques and sampling resolution requirements for groundtruthing multibeam acoustic backscatter (EM3000) and QTC VIEW™ classification technology

Both acoustic and sediment surveys were carried out in the Broughton Archipelago, British Columbia, in order to map a former aquaculture site and calibrate acoustic surveys with georeferenced sediment properties. The acoustic surveys included EM3000 Multibeam (including backscatter) and QTC VIEW™ (Series IV) technologies, while the geotechnical survey entailed Van Veen grab sampling of surface sediments and associated analyses. The two acoustic technologies were consistent in their ability to identify distinct regions of seafloor characterized by rock outcrops, consolidated substrates, or gel-mud depositional fields. Both multibeam backscatter data and QTC VIEW™ number-coded classifications were extracted across a range of circular areas located at each georeferenced sampling station (radii: 2, 3, 4, 5, 8, 12, 16, 20 m). Statistical correlations were observed between backscatter and certain geotechnical properties, such as sediment porosity, sediment grain size fractions (<2 μm, silt content), and particulate sulfur concentration. The areal resolution of backscatter extraction was explored in terms of determining a sensitive calibration technique between backscatter and sediment properties. In general the highest r² values between backscatter and sediment variables were observed across extraction radii between 8 and 20 m. Such groundtruthing techniques could be used to interpolate seafloor characteristics between sampling stations and provide a steering tool for sampling designs associated with benthic monitoring programs.     

Modelling the distribution of hard seabed using calibrated multibeam acoustic backscatter data in a tropical,macrotidal embayment: Darwin Harbour,Australia

Spatial information on the distribution of seabed substrate types in high use coastal areas is essential to support their effective management and environmental monitoring. For Darwin Harbour, a rapidly developing port in northern Australia, the distribution of hard substrate is poorly documented but known to influence the location and composition of important benthic biological communities (corals, sponges). In this study, we use angular backscatter response curves to model the distribution of hard seabed in the subtidal areas of Darwin Harbour. The angular backscatter response curve data were extracted from multibeam sonar data and analysed against backscatter intensity for sites observed from seabed video to be representative of “hard” seabed. Data from these sites were consolidated into an “average curve”, which became a reference curve that was in turn compared to all other angular backscatter response curves using the Kolmogorov–Smirnov goodness-of-fit. The output was used to generate interpolated spatial predictions of the probability of hard seabed (p-hard) and derived hard seabed parameters for the mapped area of Darwin Harbour. The results agree well with the ground truth data with an overall classification accuracy of 75% and an area under curve measure of 0.79, and with modelled bed shear stress for the Harbour. Limitations of this technique are discussed with attention to discrepancies between the video and acoustic results, such as in areas where sediment forms a veneer over hard substrate.