Swath bathymetry with GLORIA

For many years, GLORIA has been producing sonar images of the deep ocean floor. In the mid-1980's, the SeaMARC II system came to prominence producing depth values as well as sonar images. The basic method compares the phases of the signals returning from the seafloor to two rows of transducers. The phase differences are converted into angles of arrival and together with the arrival times converted into range and depth values. This capability has now been added to the GLORIA system. The fact that GLORIA uses a 2s FM pulse means the backscattered reverberation can come from a strip of seafloor up to 1.5 km wide. To accommodate this, overlapping complex FFT's are used to produce a time-frequency matrix for the returning signals. In this matrix, a constant range feature appears as a diagonal. Phases are then calculated using a least-mean-squares estimate along diagonals. The main source of error and bias is due to surface reflection, and this is taken into account. The GLORIA swath bathymetry system was tested on two cruises and it was possible to produce depth contours with a good level of confidence. The total swath width was over eight water depths and would have been greater with a more favorable velocity profile. Comparison with other bathymetry data (such as multibeam systems) showed excellent correlation, having a standard deviation of only 4% of total water depth     

Swath Bathymetry Processing of GLORI-B and SeaBeam 2000

A new bathymetry processing software package has been developed to postprocess new GLORI-B swath bathymetry data using preexisting techniques. GLORI-B bathymetry is calculated using an interferometry (phase delay) method using the modified GLORIA towfish which has parallel rows of transducers on both sides. We describe four types of artifacts observed during the first use of this new system during Legs 5 and 6 of the Gloria Expedition which surveyed the fastest spreading segment of the global seafloor spreading system and the broad chain of volcanoes near Easter Island. These artifacts include cross-track bias, along-track bias, a 'dropped edge' effect, and random noise. We describe and illustrate how we minimize these artifacts. We merge the SeaBeam 2000 bathymetry data with the GLORI-B bathymetry data to produce a final bathymetric mosaic which covers about 243,400 km² and shows a different style of diffuse widely spread volcanism not previously observed along hotspot chains. The data are used in several studies describing seamount morphology, elastic thickness of the lithosphere, tectonic and geochemical evolution of the area, and mantle flow from a hotspot to a superfast seafloor spreading center.     

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.     

Correction for depth biases to shallow water multibeam bathymetric data

Vertical errors often present in multibeam swath bathymetric data. They are mainly sourced by sound refraction, internal wave disturbance, imperfect tide correction, transducer mounting, long period heave, static draft change, dynamic squat and dynamic motion residuals, etc. Although they can be partly removed or reduced by specific algorithms, the synthesized depth biases are unavoidable and sometimes have an important influence on high precise utilization of the final bathymetric data. In order to confidently identify the decimeter-level changes in seabed morphology by MBES, we must remove or weaken depth biases and improve the precision of multibeam bathymetry further. The fixed-interval profiles that are perpendicular to the vessel track are generated to adjust depth biases between swaths. We present a kind of postprocessing method to minimize the depth biases by the histogram of cumulative depth biases. The datum line in each profile can be obtained by the maximum value of histogram. The corrections of depth biases can be calculated according to the datum line. And then the quality of final bathymetry can be improved by the corrections. The method is verified by a field test.     

Merged GLORIA sidescan and Hydrosweep pseudo-sidescan: Processing and creation of digital mosaics

We have replaced the usual band of poor-quality data in the near-nadir region of our GLORIA long-range sidescan-sonar imagery with a shaded-relief image constructed from swath bathymetry data (collected simultaneously with GLORIA) which completely cover the nadir area. We have developed a technique to enhance these pseudo-sidescan images in order to mimic the neighbouring GLORIA backscatter intensities. As a result, the enhanced images greatly facilitate the geologic interpretation of the adjacent GLORIA data, and geologic features evident in the GLORIA data may be correlated with greater confidence across track. Features interpreted from the pseudo-sidescan may be extrapolated from the near-nadir region out into the GLORIA range where they may nt have been recognized otherwise, and therefore the pseudo-sidescan can be used to ground-truth GLORIA interpretations. Creation of digital sidescan mosaics utilized an approach not previously used for GLORIA data. Pixels were correctly placed in cartographic space and the time required to complete a final mosaic was significantly reduced. Computer software for digital mapping and mosaic creation is incorporated into the newly-developed Woods Hole Image Processing System (WHIPS) which can process both low- and high-frequency sidescan, and can interchange data with the Mini Image Processing System (MIPS) most commonly used for GLORIA processing. These techniques are tested by creating digital mosaics of merged GLORIA sidescan and Hydrosweep pseudo-sidescan data from the vicinity of the Juan Fernandez microplate along the East Pacific Rise (EPR).     

Postprocessing and corrections of bathymetry derived from sidescansonar systems: application with SeaMARC II

A procedure for postprocessing bathymetry data provided by a phase-measuring sidescan sonar system is presented. The data were collected with the SeaMARC II system, and are generally characterized by a high level of noise and uneven spatial sampling. Before any spatial filtering is applied, data are selected to remove most of the obvious artifacts and to retain instantaneous depth profiles whose slant ranges increase monotonically from a central location to the edges of the swath. An extrapolation scheme, patterned after a potential field, is proposed to fill gaps in the coverage or to extend the bathymetric swath to that of the corresponding sidescan image when regridding the data to a rectangular frame. To fill the near nadir gap typically found in these data, a specific interpolation methodology is developed that takes into account the slant range of the first bottom return as received by the sidescan sonar itself or by a shipboard echo-sounder. Spatial low-pass filtering is applied through convolutions with parabolic windows whose width is proportional to the footprint of the acoustic beam along track and roughly 1/8 of the swath width across track. Mismatches of contour lines between adjacent tracks are reduced through a statistical method design to correct systematic profile errors     

Autonomous bathymetry survey system

This paper describes the Autonomous Bathymetry Survey System (AutoSurvey), a system that provides automation of swath sonar bathymetric surveys. This system enables faster surveying of an area through environmentally adaptive techniques while ensuring adequate coverage and data quality. AutoSurvey assesses data quality and coverage in real time and generates next-trackline waypoints based on actual system performance. The need for real-time performance assessment is discussed. A primary factor considered is the effect of the environment on swath bathymetry system performance, which is difficult to predict a priori. The system's features, design, and implementation are discussed in this paper. Simulation and sea trial results are presented, as well as an analysis of the system's ability to reduce survey time     

The effect of bathymetric filtering on nearshore process model results

Nearshore wave and flow model results are shown to exhibit a strong sensitivity to the resolution of the input bathymetry. In this analysis, bathymetric resolution was varied by applying smoothing filters to high-resolution survey data to produce a number of bathymetric grid surfaces. We demonstrate that the sensitivity of model-predicted wave height and flow to variations in bathymetric resolution had different characteristics. Wave height predictions were most sensitive to resolution of cross-shore variability associated with the structure of nearshore sandbars. Flow predictions were most sensitive to the resolution of intermediate scale alongshore variability associated with the prominent sandbar rhythmicity. Flow sensitivity increased in cases where a sandbar was closer to shore and shallower. Perhaps the most surprising implication of these results is that the interpolation and smoothing of bathymetric data could be optimized differently for the wave and flow models. We show that errors between observed and modeled flow and wave heights are well predicted by comparing model simulation results using progressively filtered bathymetry to results from the highest resolution simulation. The damage done by over smoothing or inadequate sampling can therefore be estimated using model simulations. We conclude that the ability to quantify prediction errors will be useful for supporting future data assimilation efforts that require this information.     

Extensional faulting and segmentation of the Mid-Atlantic Ridge north of the Kane Fracture Zone (24° 00′ N to 24° 40′ N)

The combination of multi-beam echo-sounder swath bathymetry and high-resolution deep-towed sidescan sonar provides a powerful database from which to examine mid-ocean ridge processes. We have used such a database, gathered from the Mid-Atlantic Ridge north of the Kane Fracture Zone (the MARNOK area), to examine the relationship between tectonic, volcanic, and bathymetric segmentation. We have identified structural domains, with different fault distributions, and neovolcanic segments that are distinct from the 2nd or 3rd order bathymetric segmentation.From their mutual relationships, a model is proposed for the magmatic accretion of oceanic crust at slow spreading ridges that relates the local melt supply to the tectonic style. We suggest that these are mutually interactive, and determine whether volcanic extrusion along the ridge is continuous and slow, or episodic and rapid.     

利用相邻测线重叠区域进行多波束测深横摇运动残差改正

针对多波束测深条带边缘波束易受到姿态和声速等多种误差影响、相对中央波束数据质量较低的问题,本文提出一种利用相邻测线重叠区域对多波束测深数据边缘波束进行横摇运动残差改正的模型,提高边缘波束测深数据的质量。使用沿航向的测深点匹配插值模型,完成中央波束测深点与边缘波束测深点的匹配,得到边缘波束测深误差值;使用横摇运动残差改正模型,实现顾及姿态角的条件下补偿波束入射角。计算实例表明：本文模型能够较为准确地提取边缘波束测深误差值,改正后的海底地形削弱了误差导致的上下起伏,有效地减少了影响边缘波束的多种误差,具有实际的工程应用价值。     

A new method for weakening the combined effect of residual errors on multibeam bathymetric data

Multibeam bathymetric system (MBS) has been widely applied in the marine surveying for providing high-resolution seabed topography. However, some factors degrade the precision of bathymetry, including the sound velocity, the vessel attitude, the misalignment angle of the transducer and so on. Although these factors have been corrected strictly in bathymetric data processing, the final bathymetric result is still affected by their residual errors. In deep water, the result usually cannot meet the requirements of high-precision seabed topography. The combined effect of these residual errors is systematic, and it’s difficult to separate and weaken the effect using traditional single-error correction methods. Therefore, the paper puts forward a new method for weakening the effect of residual errors based on the frequency-spectrum characteristics of seabed topography and multibeam bathymetric data. Four steps, namely the separation of the low-frequency and the high-frequency part of bathymetric data, the reconstruction of the trend of actual seabed topography, the merging of the actual trend and the extracted microtopography, and the accuracy evaluation, are involved in the method. Experiment results prove that the proposed method could weaken the combined effect of residual errors on multibeam bathymetric data and efficiently improve the accuracy of the final post-processing results. We suggest that the method should be widely applied to MBS data processing in deep water.     

The filtering and compressing of outer beams to multibeam bathymetric data

Some errors and noises are often present in multibeam swath bathymetric data. Echo detection error (EDE) is one of the main errors. It causes the depth error to become bigger in outer beams and looks like sound refraction. But depth errors due to EDEs have a trumpet-shaped appearance, instead of a curved appearance that is caused by the sound refraction errors. EDEs, including systematic acoustic signal detection errors and internal noises, cannot be removed during the correction of sound refraction. It causes depth inconsistencies between adjacent swaths and degrades precision of outer beams. Sometimes, the bathymetric errors caused by EDEs do not even meet the requirements of IHO (International Hydrographic Organization). Therefore, a post-processing method is presented to minimize the EDEs by filtering outliers and compressing outer beams of multibeam bathymetric data. The outliers caused by internal noises are removed by an automatic filter algorithm first. Then the outer beams are compressed to reduce systematic acoustic signal detection errors according to their depths, the calculated depth line and standard deviations (SDs). The automatic filter process is important for calculating the depth line. The selection of inner beams to calculate the average SD of beam depths is crucial to achieving compressing goals. The quality of final bathymetric data in outer beams can be improved by these steps. The method is verified by a field test.     

机载激光测深精度外部检核方法

给出了单波束测深的原理,分析了单波束观测数据预处理模型,提出了采用单波束测深成果检核机载激光数据质量的技术方法,并以我国自行研制的机载激光测深系统为例,给出了该系统在某海区试验数据的外部检核结果.针对两种测深手段之间明显存在系统性偏差的问题,提出了以单波束测深成果为控制,对机载激光测深数据系统偏差进行校正和补偿的处理方...     

Acoustic imagery using a multibeam bathymetric system

Historically, measurement and collection of deep‐ocean acoustic imagery are accomplished by towed sidescan systems. Recently, work has been performed to extract acoustic imagery from current hull‐mounted wide‐swath bathymetric sonars with minimal hardware modification. Past work of deriving acoustic imagery from swath sonars has been performed primarily with SeaBeam's sixteen 2^2/3 ° preformed beams. The Navy is investigating the feasibility of extracting an acoustic image from the Sonar Array Survey Systems (SASS), a high‐resolution (1^o beams) wide‐fan (90°) bathymetric system. Due to the large data volume (approximately 1 MB per ping), SASS normally discards the raw acoustic returns once bathymetry is calculated. In early 1991 the Naval Air Development Center (NADC) installed the hardware on board the USNS Maury to capture and record the raw acoustic signal (inphase and quadrature) from the SASS's 144 hydrophones for later inversion to a backscatter image. Preliminary qualitative mosaics of the sidescan images show promising results and warrant further development.     

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%.     

Signal processing strategies for a bathymetric sidescan sonar

A dominant source of errors in swath bathymetry is acoustic interference. In 1989 the author published an analysis of these errors and predicted depth accuracies for a system which reduced their effect by averaging. This present paper shows how a considerable improvement in performance may be obtained by a variety of signal processing strategies that include the use of several widely spaced receivers and the elimination of the most unsatisfactory measurements before averaging. Simulations show how impressive sea bed profiles can be produced with a single ping, even at low signal-to-interference ratios     

Quantitative backscatter measurements with a long-range side-scansonar

It is shown that useful relative backscatter strengths can be calculated from GLORIA long-range side-scan sonar data using a simple acoustic model. The calculation was performed on GLORIA side-scan sonar data collected during 1987 in the southern Indian Ocean. GEOSECS hydrographic information was used to access the effects of refraction (ray bending and aspherical spreading signal losses). Sea Beam bathymetry was used to correct the effective insonified area and compute the grazing angle. A major difficulty in performing this calculation over the terrain chosen (mid-ocean ridge topography) was one of adjusting navigation so that small features in Sea Beam and GLORIA data matched. Preliminary results show a 10-dB falloff in backscatter strength with decreasing grazing angle (10°-40°) at 6.5 kHz over what must presumably be a rough surface (extruded basalts and breccias)     

The tectonic fabric of the Barbados Ridge accretionary complex

A map showing near-surface structure and other surficial geological features on the Barbados Ridge accretionary complex has been constructed from GLORIA long-range side-scan sonographs, detailed bathymetric surveys using Seabeam, and single and multichannel seismic reflection sections. The map shows major variations in the size, orientation and style of structures across the accretionary complex that are principally related to the increasing thickness of the complex and the sediment drape upon it, and variations along the complex that derive from the northward thinning of sediment accreted from the ocean floor and the effects of ridges in the oceanic basement. The development of crossfaulting is mainly associated with basement ridges. Mud diapirism is abundant in the southern part of the complex and is characteristic of that part of the complex into which turbidites from the Orinoco submarine fan are accreted. Slumping on a variety of scales is another feature of the complex.The map is presented at a scale of 1:1063415 (4 inches per degree of longitude) and also shows bathymetry at 200 m contour intervals.     

Digital preprocessing techniques for GLORIA II sonar images

Image processing techniques are discussed that correct distortions in GLORIA II side scan sonar imagery including water column offset, slant-range distortion, multiple returns, aspect ratio, speckle noise, striping, and cross-track power drop-off. The software operates within NASA's ELAS image processing system and is applied to the original 12-bit GLORIA II data. Procedures are discussed for generating large scale mosaics and three-dimensional overlays with sea floor bathymetry. The results are shown in four sonographs acquired off the southern coast of California.     

Interference fringes on GLORIA side-scan sonar images from the Bering Sea and their implications

GLORIA side-scan sonographs from the Bering Sea Basin show a complex pattern of interference fringes sub-parallel to the ship's track. Surveys along the same trackline made in 1986 and 1987 show nearly identical patterns. It is concluded from this that the interference patterns are caused by features in the shallow subsurface rather than in the water column. The fringes are interpreted as a thin-layer interference effect that occurs when some of the sound reaching the seafloor passes through it and is reflected off a subsurface layer. The backscattered sound interferes (constructively or desctructively) with the reflected sound. Constructive/destructive interference occurs when the difference in the length of the two soundpaths is a whole/half multiple of GLORIA's 25 cm wavelength. Thus as range from the ship increases, sound moves in and out of phase causing bands of greater and lesser intensity on the GLORIA sonograph. Fluctuations (or wiggles) of the fringes on the GLORIA sonographs relate to changes in layer thickness. In principle, a simple three dimensional image of the subsurface layer may be obtained using GLORIA and bathymetric data from adjacent (parallel) ship's tracks. These patterns have also been identified in images from two other systems; SeaMARC II (12 kHz) long-range sonar, and TOBI (30 kHz) deep-towed sonar. In these, and other cases world-wide, the fringes do not appear with the same persistence as those seen in the Bering Sea.