EVALUATION OF THE PRECISION OF PHASE-MEASURING BATHYMETRIC SIDE SCAN SONAR RELATIVE TO MULTIBEAM ECHOSOUNDERS

Arthur Ayres Neto

Lieutenant Commander Christopher Florentino

Lieutenant Commander Vitor Bravo Pimentel

Published

Last Updated

Abstract
The International Hydrographic Organization (IHO) and the Directorate of Hydrography and Navigation (DHN) consider that where the under-keel clearance is critical, bathymetric surveys should provide full seafloor coverage. Multibeam Echosounders (MBES) and Phase-Measuring Bathymetric Side Scan (PMBS) are capable sensors to achieve this coverage, however, they have different technologies. The MBES is widely accepted by the international community and considered the standard system in surveys that require the highest IHO accuracy. The aim of this paper is to assess whether the PMBS can achieve comparable results with the MBES, focusing on Nautical Chart production. Despite “noisy” data and some effects, this research suggests that PMBS can be used to meet IHO requirements, under some special conditions.

1. Introduction

An important component of a hydrographic survey is the bathymetric measurements. Currently, Multibeam Echosounder (MBES) systems are the most used sensor for morphological investigation and seafloor mapping. It fulfills the International Hydrographic Organization specifications (IHO, 2008) and the Brazilian Hydrographic Service requirements, due to its swath and low uncertainty.

However, the Phase-Measuring Bathymetric Sidescan (PMBS) has been presented as an alternative technology to the MBES. Some studies were presented by Gostnell (2005), Gostnell, Yoos and Brodet (2006), Madricarlo, Foglini and Tonielli (2011), Dodd (2013), Brisson, Wolfe and Staley (2014), Ai, Armstrong and Fleury (2015), Jerram and Schmidt (2015a), Ma, Xu and Xu (2016) and Bongiovanni and Schmidt (2016), where issues regarding data acquisition, processing methodologies and potential uses of this equipment were discussed.

Even after positive results, the PMBS technology remains under discussion. The equipment is applied by some hydrographic services (e.g. USACE – USA) and not fully adhered to by others (NOAA – USA). In Canada, regulations are provided to be followed, but no specific equipment is mentioned to achieve them. In Brazil, the present standard procedures do not define the acceptance of any equipment, but several requirements are usually achieved by MBES (DHN, 2017). The use of PMBS is, therefore, very incipient in the country, owing to both the existing world controversy and the lack of national studies with methodical hydrographic rigor focusing on regional scenarios (oceans, rivers, lakes and seas).

This work adds to previous efforts in the analysis of PMBS performance, seeking to improve consistency to the evaluation of the PMBS suitability for Special Order hydrographic surveys. To accomplish this, the IHO and the Brazilian Hydrographic Service specifications (defined by the Directorate of Hydrography and Navigation) were strictly observed.

2. Theoretical concepts

MBES are sonars capable of providing bathymetric data by crossing the acoustic waves through a transmitting and a receiving sensor, which are arranged orthogonally, in a technique known as mills cross (Bjørnø, 2011; L-3 COMMUNICATIONS SEABEAM INSTRUMENTS, 2000; Simmons et al., 2017). In this arrangement, the beams are formed at the intersections between the transmitting and the receiving wavefronts (L-3 COMMUNICATIONS SEABEAM INSTRUMENTS, 2000). These beams, after the two-way travel time, provide the seafloor signal detection that is used to calculate the depth, which is done by amplitude or phase detection (IHO, 2005).

The PMBS is a sensor capable of simultaneously providing backscatter and bathymetric information (Jerram and Schmidt, 2015a). It consists of a transmitter and more than one receiver spaced at a distance multiple of the signal wavelength “λ” (KONGSBERG MARITIME, 2014; Lurton, 2000; Wilby, 1999). Georeferenced depth information is provided by measuring the echo’s phase difference at the receiver transducers (IHO, 2005).

In a quick comparison, the PMBS has a lower cost (Ai and Parent, 2011); has a greater swath capability in shallow water (Brisson; Wolfe; Staley, 2014; Gostnell, 2005; Jerram and Schmidt, 2015b); provides a faster survey with better visualization at the edges of the swath (Brisson; Wolfe; Staley , 2014); has a higher data density (USACE, 2013), leading to a larger across-track resolution (1-10 mm) (Jerram and Schmidt, 2015a; Jerram and Schmidt, 2015b; Ma; Xu; Xu, 2016); has the same backscatter view of a sidescan sonar (Brisson; Wolfe; Staley, 2014; Jerram and Schmidt, 2015b); and uses the principle of interferometry, which is also used in MBES phase detection of the outer beams.

Despite these characteristics, the PMBS is not yet fully accepted by the international hydrographic community (Brazil included). This is due to a number of factors such as low detection capacity in nadir (Dodd, 2013, Ma, Xu, Xu, 2016); very noisy data (Dodd, 2013; Jerram and Schmidt, 2015b); depth ambiguity (layover); large volume of data, which requires better computing capability (Jerram and Schmidt, 2015b; Ma; Xu; Xu, 2016); and strong vertical uncertainty (Brisson; Wolfe; Staley, 2014; Jerram and Schmidt, 2015b; NOAA, 2017).

3. Methodology

3.1 TESTS AND AREAS

In order to evaluate the acceptance of PMBS by the Brazilian Hydrographic Service, tests were developed to assess certain critical capacities, while comparing the results with a simultaneously acquired MBES survey. To perform this experiment, a vessel-pole-mounted PMBS and vessel-organic hull-mounted MBES were used.

Five test sites were chosen within the Guanabara Bay (Figure 1). The Guanabara Bay has an area nearly 384km2, with a central channel approximately 30 meters deep. The seafloor has an average depth is 5.7 meters and consists of quartz sand at the entrance to the Atlantic Ocean, associated to the effects of waves and tides. The remainder of the seafloor, is mostly mud, resulting from Holocene transgression and river sedimentation (Kjerfve et al., 1997).

Figure 1 – Location of the test areas and the tide gauge stations. Source: Adapted by authors.

Test Area 1 aims to verify the ability of PBMS to precisely detect and delineate bathymetric features. The area is about 13 meters deep and approximately 0.375km2 wide (750m x 500m), with numerous rock outcrops with varied sizes and shapes.

Test Area 2 is designed to check the ability of PBMS to resolve slopes with shallow to steep gradients. The chosen area covers around 0.350km2 (550m x 700m) in the vicinity of Villegagnon Island, with varied gradients and depths ranging from 1.5m to 39m. In this way, several slope gradients (approximately 4 to 10 degrees) were covered, providing data for different analyses.

Test Area 3 aims to observe the maximum distance the PMBS can detect objects that are pronouncedly off the seafloor and estimate their height. The detection of objects with substantial vertical relief is a fundamental factor for safety of navigation. For this purpose, an area of 0.1Km2 (250m x 400m) and with an average depth of 15 meters was selected. The area comprises 5 pillars of the Presidente Costa e Silva Bridge which rise above the water. This allows simulation of targets that the survey vessels cannot approach due to the risks involved.

Test Area 4 seeks to detect the shoalest depths of pronounced vertical features that do not reach the water surface. Thus, a wreck area of 0.330Km2 (600m x 550m) was chosen, with an average depth around the object of 30 meters, with a minimum known depth of 20.8 meters.

Test Area 5 aims to verify the efficiency of the PMBS compared to MBES. The survey line spacing was chosen according to interferometric parameters. It is known from the literature that the PMBS has a maximum swath width of approximately 10-12 times the water depth, while the MBES has 5-6 times the water depth (Jerram and Schmidt, 2015b). The purpose was to verify the useful coverage efficiency of the PMBS. For this test, we choose an essentially flat area of approximately 1.4km2 (1,750m x 800m) with depths varying from 3 to 20 meters.

Test Areas 3 and 4 were also used for the resolution tests; verifying bathymetric resolution and data density.

3.2 EQUIPMENTS

The survey was performed on August 02nd and 03rd of 2017. The equipment used were:

A blister-mounted Kongsberg EM2040 multibeam echosounder, operating at 400KHz central frequency; INS Kongsberg Seapath 300 with MRU5 and 2.582m spaced heading antenna; GPS positioning with Fugro DGNSS 3710 correction; AML Micro X sound velocity sensor; and Kongsberg SIS navigation and acquisition software.

A pole-mounted Kongsberg Geoacoustics Geoswath 4R phase-measuring bathymetric sidescan, operating at 250KHz central frequency; INS Applanix POS MV WAVE MASTER 2 and 1.9m spaced heading antenna; GPS positioning with RTK Marinestar Fugro correction; Valeport Mini SVS 25mm sound velocity sensor; and Geoswath Plus acquisition software.

The Sound Velocity (SV) Profile information was acquired using the AML Minos X sensor at 10cm resolution. SV data was processed into the SIS software and uploaded into the data during acquisition in MBES and during the post-processing in PMBS. The profiles were acquired before commencing the survey in each area and in every time the SV varied more than 2 m/s at the face -of- transductor sound velocity sensor.

The tide information was collected at tide gauge stations of Boqueirão, Fiscal Island and the Ponta da Armação Naval Complex (CNPA, in Portuguese) (Figure 1).

Throughout the survey, the weather was fine, with low wind speeds (0-3 on the Beaufort scale).

3.3 PROCESSING DATA

MBES data were processed in CARIS HIPS & SIPS 10.2. The PMBS data were processed in the GEOSWATH PLUS software. Both data were then exported to CARIS HIPS & SIPS 10.2 for the CUBE surface generation at different resolutions.

3.4 QUALITY CONTROL

The quality control was performed by measuring the precision and accuracy of each surface node.

3.4.1 Precision

The precision was verified by two methods: 1 – calculating the CUBE hypothesis TPU (Calder and Mayer, 2013); 2 – calculating the sample uncertainty of each node.

In the first method, if 95% of the nodes had a TPU within the recommended range for a given IHO Order, the surface would be accepted for that Order. In the second method, the node sample uncertainty was calculated from the standard deviation multiplied by 1.96 (in order to have 95% confidence, considering a Normal distribution) of all the considered valid depth information that is contained within a node (Pereira, 2016). In practical terms, it represents the dispersion of the data and is equivalent to the first one. If 95% of the nodes show sample uncertainty within the range recommended for a particular IHO Order, the surface will be considered in that Order.

While the first method takes into account the uncertainties that give rise to depth information, the second is based on the coherence of the bathymetric data.

3.4.2 Accuracy

The accuracy of PMBS surfaces was evaluated in three ways: (1) observing the dispersion of the MBES data (here considered as the reference surface) relative to the PMBS depth surface (here called layer DepthPMBS); (2) by calculating the depth difference surface between the MBES (DepthMBES) and the PMBS (DepthPMBS) surfaces; (3) by calculating the surfaces that contain the difference of the shallowest depth information, hereinafter referred to as Shoal layers (ShoalMBES and ShoalPMBS).

The multibeam data dispersion from the DepthPMBS layer was done using the “QC Report” tool of the CARIS HIPS & SIPS software (CARIS, 2017). The program checks the statistics of some standard lines (MBES lines) against the surface that will be evaluated (DepthPMBS). Based on the surface depth values, another two above-and-below parallels surfaces are virtually created, ranging the TVU IHO Orders. If, the percentage of pings of each line that lies within this range is higher than or equal to 95%, the surface can be considered representative of the order it was evaluated against (Pereira, 2016).

The difference surface is calculated by subtracting the values of the MBES and PMBS depth layer (DepthMBES – DepthPMBS). Such procedure provides a qualitative comparison opportunity.

The layer containing the difference of the shallowest depth is very useful for hydrographic services. In Brazil, to ensure navigation safety, the surface layer used to create the nautical chart is the one that contains the shoalest depths found at each node. Although this surface may not be the most accurate or suitable when it comes to the real seafloor identification, it is useful to guarantee there will be no obstacles shallower than those selected and reported during the data processing.

4. Results

Table 1 gives an integrated perspective of the results. It is possible to observe that the analysed criteria from all areas have very similar statistical results.

Table 1 – Statistical Results for all surveyed areas. Organized by the authors.

Regarding the accuracy of each test site, it is possible to observe that the mean hypothesis uncertainty calculated by the CUBE algorithm was lower than the one recommended for Special Order for both methods, so that all surface nodes satisfy this requirement. Regarding the node standard deviation, the MBES showed a considerably smaller deviation (Figure 2b) than the PMBS (Figure 2c)., meaning that the valid data from the second sensor (PMBS) are more scattered. Due to the small dispersion of the node standard deviation, the MBES data was classified as Special Order, while the PMBS data can only be classified as Order 1A. Another aspect to be observed is that the average density of the pings per node was 2.4 to 5.6 times greater for the PMBS than that of the MBES.

When the accuracy is analyzed, it is observed the PMBS surface reached Special Order. The difference between the depth surfaces (Depth layers) ranged from 0 to 0.1m, with the MBES surface being shallower. This indicates that the two sensors estimated the depths in a remarkably similar way. The situation changes when we consider the difference between the Shoal layers. In this case, the difference varies between 0.5m and 0.7m proving again the dispersion of the PMBS data. If the data considered valid from the PMBS had cohesion like the MBES, the range of depth difference between the Shoal layers would be equivalent to the differences observed for the Depth layers.

Figure 2 a) Rock outcrops in Test Area 1. The yellow line marks the profiles shown in the frames bellow; b) a MBES data profile; c) a PMBS data profile; d) a PMBS data profile with the reject data

4.1 AREA 1 AND 2 – LAYOVER EFFECT STUDY CASE

A layover effect occurs when two different echoes arrive at the same time at the transducer. As there are no sectoral angles of reception, as in the multibeam, the sensor processes the two signals as noise. Sometimes, this ambiguity also happens when the sonar beam reaches the top of a tall feature before it reaches the base of it and is more critical in nadir and steep gradient regions. From Figure 3, one can better visualize how this effect occurs in sloping features. In this figure, the seafloor is represented in blue, an imaginary seafloor with higher slope represented by the dotted green line, the sensor by the yellow dot and the wavefront in orange. It is possible to notice that, in the case of the blue seafloor whose slope is smaller than the opening angle θ, the “echo B” will have traveled a greater distance than the “echo A”. If the green seafloor is considered, we have the opposite situation and “path B” would be shorter than “path A”, arriving at the sensor at the same time as some other echo reflected by the seafloor, creating ambiguity and, therefore, noise. Thus, it can be concluded that if the distance travelled “B” ≤ “A”, there will be information ambiguity, and that if the slope of some seabed feature is greater than the complement of the grazing angle α, the layover effect of will occur.

The PMBS useful swath is 10 to 12 times the water depth, so there are incident angles ranging from approximately 90° to 5.7° to 4.3°. Thus, it is known that features with slopes approximately greater than 85.7° will always be subject to layover effect.

Figure 3 – Layover Effect principle

However, adjacent survey lines tend to reduce the layover effect if the gradient is not so abrupt. In a 200% coverage situation, the terrain will be ensonified at least by two swaths that will have points with different aperture angles, so that while one of the swaths may show layover effect, the other may not. When the seafloor or target has a slope greater than the maximum incident angle (around 5°), this effect cannot be remedied. Therefore, this effect was not seen in Area 2, where the seafloor slope is very low, but can be seen in the rock outcrops of Area 1 where the slopes of the rock edges are larger than the PMBS maximum transmission opening angle. An interesting example is the rock shown in Figure 4a (inside the red circle). On the softer gradient side, there is no layover effect, while in the most abrupt side it was evident.

It can be observed in the Shoal Depth surfaces difference (Figure 4a) that the shape of the rocks does not get very well defined in the PBMS data compared to the MBES data (Figure 4b), resulting in depth misestimation. The largest ones lie exactly at the steep edge of the rocks (Figure 4a). Because these areas are very steep, they are susceptible to the layover effect. In Test Area 2, due to the softer gradients (less than 15 degrees of slope), the differences between the Depth layers are not so noticeable.

Figure 4 a) Shoal Difference Surface Layer with the rock’s edges highlighted; and b) MBES Depth Layer showing the analyzed rocky field in the Test Area 1.

Probably the layover occurred in both areas, but due to the overlapping of the survey lines, it was not noticed in Area 2. In Area 1 (Figure 4b), despite the overlapping of different swaths, the slope of the stone edges is beyond the grazing angle, a situation in which the layover will always appear, independent of overlaps.

4.2 AREA 3

Area 3 stands out for the existence of the bridge pillars (Figures 1, 5a and 5b). These features simulate targets with a vertical projection that reaches the water surface, making it impossible for the survey vessel to approach to the abutments.

The detection of these objects by the PMBS is more pronounced than in the MBES, as can be seen through the of comparative bathymetric profiles (Figure 5c) between the surfaces generated by both sensors. The Shoal layer in the PMBS shows that this sensor could image and quantify the height of bridge pillars almost to the water surface, although they were not totally identified in the Depth layer. This is the result of different acquisition technology and data processing technique. Automatic processing was not able to detect the abutments, requiring manual operator intervention. However, once this parameter is adjusted, it is observed that the PMBS was able to detect the targets to approximately 0.5m below the sea surface. The MBES Shoal layer depicts the targets to approximately 10m water depth, while the Depth layer recorded only a small elevation.

Due to the proximity between the pillars of the bridge, another situation that can be simulated is the survey of a narrow channel or river, with no space to run adjacent lines. In the lines with no swath overlapping, the nadir gap is visible (extreme left and right of Figure 5b). In these cases, the area immediately below the sonar has a lack of bathymetric information, which is a disadvantage of the PMBS system relative to MBES.

It is also noticeable the PMBS swath is larger. This makes it easier to survey in the vicinity of targets or features that pose risks to the vessel and in places very close to river banks or to the coast. In these cases, the PMBS imaging will cover a larger area thus posing less risk to the equipment and survey team.

Figure 5 – a) The MBES Depth layer surface; b) The PMBS Depth layer surface. In the a) and b) frames, the depths are ranging from 1m – dark red – until 20m – violet. The red dashed line indicates the profile in frame c; and c) A comparative frame showing the profiles. There are six different profiles: MBES Depth layer (black), MBES Shoal layer (green), PMBS Depth layer (red), PMBS Shoal layer (blue), Depth difference layer (purple) and Shoal difference layer (yellow).

4.3 AREA 4

This area is characterized by the presence of a shipwreck (Figures 1, 6a and 6b). The multibeam shoal layer is slightly shallower than the corresponding PMBS shoal layer (Figure 6c), reflecting the amount of noise captured by the PMBS system during acquisition. Although the target has been fully ensonified, the true depth information was neither properly quantified nor transformed into a real bathymetric information. The masking of true depth information is quite common on such situations when the data is rejected as noise due to inconsistencies in the depth determination in a way similar to the layover effect.

In flat or low gradient areas over the wreck, on both layers (Depth and Shoal), the bathymetry acquired with PMBS is shallower than that acquired with MBES, as observed in the other areas.

Figure 6 – The stern of the wreck. a) The MBES Depth layer surface; b) Shoal Difference layer surface. The red line indicates the profile in frame c; and c) A comparative frame showing six different bathymetric profiles: MBES Depth layer (black), MBES Shoal layer (green), PMBS Depth layer (red), PMBS Shoal layer (blue), Depth difference layer (purple) and Shoal difference layer (yellow).

4.4 AREA 5

The efficiency test has the objective of comparing the time to survey a given area or the size of a surveyed area in a given time. In this work, we used the second approach.

It is known the PMBS and the MBES swath may be different for a given water depth. Figure 7 shows that the PMBS was able to practically achieve full coverage using the same swath width. For a better analysis of the results, the area was divided into 3 sectors: from 25.2m to 15m (Figure 8a), from 15m to 5m (Figure 8b), from 5m to 2m (Figure 8c). Figure 8a shows the

PMBS was able to easily cover 100% of the area with no need to change the swath. Meanwhile, it is clear the MBES would need more survey lines and, therefore, more time, to totally cover the same area. Figure 8b shows that MBES leaves huge gaps between the survey lines, whereas the PMBS shows only small nadir gap. In Figure 8c the gaps between the MBES lines are much more pronounced, while the gap between PMBS lines is practically non-existent.

Figure 7 – In the left, the MBES depth layer surface and in the right, the PMBS depth layer surface. The colours are related to depths ranging. From 25.2m until 15m in blue, from 15m until 5m in green and from 5m until 2m in purple.
Figure 8 – a) Test Area 5 field, where the depths vary from 25.2 to 15m; b) Test Area 5 field, where the depths vary from 15 to 5m; and c) Test Area 5 field, where the depths vary from 15 to 5m. In frames a), b) and c), the red line indicates the where the profile was taken. In the profiles graphic, there are two different lines: MBES Depth layer (black) and PMBS Depth layer (red).

While the MBES had a constant swath of approximately 3 times the water depth for the whole area, the PMBS achieved swaths of approximately 3.4 times the water depth in areas 20m deep, 6.7 times the water depth in areas 9m deep and 12 times the water depth in areas 3.5m deep. Therefore, the test showed that PMBS is more efficient when the water depth is 15 meters or shallower. When the whole area is considered, the PMBS efficiency was 37.4% higher than the MBES.

4.5 RESOLUTION TEST

For the resolution test, it is interesting that the areas have noticeable bathymetric features. Thus, we will only address different resolutions in the Areas 1 and 4.

Table 2 correlates resolution and ping density per node for both areas. In a first analysis, we used a quantitative approach to evaluate the surfaces considering that, according to Calder and Mayer (2003), in the CUBE interpolator, the minimum desirable number of pings to compose a cell is 11. According to this criterion, resolutions of 0.1m on both sensors and 0.2m on MBES are discouraged.

Table 2 – Data density per node for different resolutions. Organized by the authors.

Nevertheless, surfaces must also have their resolutions assessed qualitatively through visual analysis by observing how the morphological features behave in each image shown in Figure 9. As the resolution improves, there is an increase in the number of blanked areas and noise. In contrast, it is possible to have a better notion of the seafloor features, as well their details and texture.

Another aspect to be explored in the image visualization is the difference in the target’s size. Bigger targets (such as seen in Area 4) tend to be noisier, masking some target details (Figure 9, Test Area 4, PMBS frames). In general, despite the greater number of pings per node the PMBS data provided greater statistical robustness, the data is more scattered, generating noise. If we degrade the resolution to suppress the noise, details are also be smoothed.

Achieving a trade-off between the size of blanked areas, amount of noise and the visualization of target details is what will provide the optimal resolution for each case. Factors such as target size, depth, survey purpose should also be considered when defining the resolution to be used during the survey. Although the PMBS can produce a much larger number of pings per node than MBES, it presented a much noisier data, making it difficult to work with resolutions better than 30cm.

Figure 9 – Different depth layers resolutions for Test Areas 1 and 4. The rocky images are from Test Area 1 and the wreck images are from Test Area 4. The resolutions are in the right side.

5. Discussion

The sensors were evaluated for their relative precision in different conditions. The assessment was done in two different ways: with respect to the hypothesis and the node sample uncertainties. The first was meant to verify how the TVU was propagated, while the second was focused on the evaluation of the data dispersion. The results showed that MBES data was able to achieve Special Order on both criteria. The PMBS hypothesis uncertainty also matched Special Order since all the information that supported the bathymetric information had low uncertainty. However, its data showed high dispersion level, making the node sample uncertainty not able to reach Special Order standard, but only Order 1A. To better understand this issue, Figure 10 depict the MBES and PMBS hypothesis and node sample uncertainty for Area 1.

Figure 10 – MBES and PMBS surface assessed for Special Order. Left column is the result for hypothesis uncertainty. Right column shows the results for node sample uncertainty. The nodes that achieved Special Order in each criteria are in green. The ones that did not are in red.

It was expected to see the nodes related to a very steep seafloor in red in all situations and for both systems, as seen for the existing rock outcrops in this area. The relevant fact here is that PMBS surface achieve Special Order in the hypothesis uncertainty criteria, but did not for the node sample uncertainty, even under flat seafloor condition. This shows that not only is the PMBS data dispersion is to high, as demonstrated by Brisson; Wolfe; Staley (2014); Jerram and Schmidt (2015b) and NOAA (2017), but also how careful one must be when evaluating a bathymetric surface. Figure 11 shows, on the other hand, that PMBS data matched Order 1A of the IHO standard, indicating that these systems may be used if the demands for the survey are not so rigorous (Table 1).

Figure 11 – PMBS surface assessed for Order 1a in the node sample uncertainty criteria. The nodes that achieved Order 1a are in green. The ones that did not are in red.

The accuracy was analyzed in 3 ways (Table 1). First, using CARIS HIPS & SIPS QC Report tool, it is observed that all areas surveyed with the PMBS were within the Special Order. The second way was to compare the Depth layers surface generated with MBES and PMBS data. This showed that the depth difference ranged between 0 and 10cm, with the PMBS layer always representing shallower depths. The third way was by comparing the layers representing the shallowest depths (shoal layers), since this is the preferential data exported to the nautical charts. In this case the depth difference ranged among 0.5m to 0.7m for all study areas. This results from the data dispersion level, which is high for the PMBS Shoal layer. The problem here is that exporting the PMBS Shoal layer to the nautical chart would decrease the risks to navigation safety, but would have an impact by decreasing the navigable draft of a channel

Another important aspect is the amount of noise in the raw data. The whole PMBS data shows a large amount of noise (see Figure 2d), as already demonstrated in Dodd (2013) and Jerram and Schmidt (2015b). This noise can be generated by several factors: ambient noise, the vessel and the equipment itself (Lurton, 2000); multipath (mainly in shallow and targeted locations) (Denbigh, 1989); shifting footprint effect, baseline decorrelation effect (Lurton, 2000); and layover effect (Ărăcin and Calin, 2000; Saucan et al., 2015; Woock and Frey, 2010). In our research the layover effect was mainly observed in the edges of very steep morphological features, while noise is mostly widespread in the data

In this work, it was considered that a good resolution for the PMBS would be 0.3m. Despite delivering a higher data density than the multibeam (USACE, 2013; Jerram and Schmidt, 2015a; Jerram and Schmidt, 2015b; Ma; Xu; Xu, 2016), the PMBS has a very large amount of spurious data, which became more evident in surfaces with higher resolutions. This fact limits the resolution improvement, restricting the PMBS to the MBES surface resolution level. The resolutions of 0.1m and 0.2m did not provide the minimum amount of pings (Table 2) to ensure the node’s statistical robustness according to Calder and Mayer (2003). At the same time, these resolution levels enhance the noise in the data, giving the surface a very “rough texture”. On the other hand, surface resolution of 0.4m or higher decrease the ability of the data to show the targets details resulting in a “out-of-focus” surface (Figure 9).

The PMBS efficiency was greater than the MBES especially at depths shallower than 15 meters. This is essentially due to its larger swaths in shallower regions (Brisson; Wolfe; Staley, 2014; Gostnell, (2005); Jerram and Schmidt, (2015b)). When tested in a flat area, the nadir gaps observed in the PMBS data, makes the need for 100% overlapping mandatory as stated by Dodd (2013) and Ma, Xu, Xu (2016).

Due to the amount of bad data in the PMBS (Figure 2a), manual processing is practically impossible. Automatic and statistical processing, despite being strongly recommended, can still exhibit flaws. We suggest that data processing shall be done statistically, to have most of the noisy data rapidly cleaned, and manually, to further clean the data and better delineate the bathymetric features. The development of a better automatic technique for PMBS data cleaning is a challenge that must be addressed in order to reduce the processing time needed to generate the final bathymetric surface.

Thus, based on the results of the present study and considering the PMBS characteristics presented above, it is possible to conclude that PMBS is mostly recommended for restricted areas, with relative low gradients, shallower than 15 meters, which have been already mapped, since small features or structures edges can be masked by the noisiness or layover effect.

In an area where the dangers to navigation are already known, the PMBS can be used, since the uncertainties related to the target’s shape can be diluted in the sounding selection or absorbed by the scale of the nautical chart. As an example, in this survey, the edge’s uncertainty (when compared to the MBES surface) of the rock outcrops seen in Area 1 was around 1m. Considering the graphic reduction value in Brazil is 0.28mm, the scales that can address such uncertainty would be the ones smaller than 1:3,571 (1/0.00028). If the edge uncertainty cannot be properly estimated, one can use the THU (Total Horizontal Uncertainty) as the value for graphic reduction. Considering the maximum recommended THU for the order 1A (5m + 5% of the depth), the scale would be even more restrictive. In this case scales approximately ≤ 1: 25000 would be accepted. However, this parameter must be based on the analyst’s experience and the purpose of the survey. Table 3 summarizes all the conclusions of the present work.

Table 3 – PMBS surveying area suggestion. Organized by the authors.

It should also be noted that countries that accept bathymetric sidescan sonars have more specific standards which are more restrictive in their criteria and requirements than the IHO S-44 itself. Despite accepting PMBS surveys, the United States (NOAA, 2017), Canada (CHS, 2013) and New Zealand (LINZ, 2010) are examples of countries that are in the situation above. Brazilian hydrographic standards, especially NORMAM-25 (Rev. 2), mostly replicate the specifications and the Orders recommended in S-44 (5th Edition) (IHO, 2008), being more lenient than the standards of the aforementioned countries. The adoption of PMBS systems for Special Order surveys in Brazil must be followed by more systematic tests with methodical hydrographic rigor.

6. Conclusion

A hydrographic survey is a complex process of obtaining information about the environment. The bathymetric measurement, as an integral part of any hydrographic survey, stands out in its complexity, due to the different variables involved and the high level of precision and accuracy required. In turn, end users are continuously demanding more reliable, accurate and less expensive surveys, since the financial and administrative consequences are high.

In this context, the PMBS arises as an option to the existing traditional MBES method. It intends to provide greater efficiency and cost reduction. However, a closer analysis reveals that one technology is not capable to simply replace the other, at least for now.

The PMBS technology has advantages and disadvantages compared to the MBES. We concluded that it is possible to operate and perform bathymetric surveys using PMBS systems, but under less restrictive conditions.

There are other factors that may drive to different results: changes in the bottom type and operating frequency may affect the bottom detection level; survey platform – using PMBS systems on an AUV would deliver very different results, especially due to the its acoustic positioning; rougher sea state and weather conditions demand more sophisticated attitude sensors affecting depth accuracy. It is important that users and equipment technicians be aware of the PMBS limitations and potentialities, exploiting its advantages to the maximum and minimizing its disadvantages. Some drawbacks, like the amount of noise and consequently the dispersion of valid data, can be mitigated by processing (filters or other algorithms) and remain an issue for future investigations.

7. References

– Ai, Y and Parent, M. (2017). Performance Analysis Using an Uncertainty Model for Hydro- Chart 5000. Sea Technology. Arlington, United States, viewed 12 November 2017, https://geo-matching.com/uploads/default/m/i/migrationvwmfdr.pdf.

– Ai, Y.; Armstrong, S. and Fleury, D. (2015). Evaluation of the Klein HydroChart 3500 interferometric bathymetry sonar for NOAA sea floor mapping. In: OCEANS 2015 – MTS/IEEE Washington, USA.

– Ărăcin, A. and Calin, A. (2013). Interferometric Bathymetry – principles and utility. Recent Advances in Geodesy and Geomatics Engineering. Proceedings of the 1st European Conference of Geodesy & Geomatics Engineering (GENG ’13), Antalya, Turkey: WSEAS Press. p. 198–205.

– Bjørnø, L. (2011). Developments in sonar and array technologies. In: 2011 IEEE Symposium on Underwater Technology (UT’11) and Workshop on Scientific Use of Submarine Cables and Related Technologies (SSC’11), Tokyo, Japan.

– Bongiovanni, C. and Schmidt, V. (2016). Evaluating Outside Source Interferometric Data for Chart Updates. In: Canadian Hydrographic Conference (CHC-2016), Halifax, Canada.

– Brisson, L. N.; Wolfe, D. and Staley, M. P. S. M. (2014). Interferometric swath bathymetry for large scale shallow water hydrographic surveys. In: Canadian Hydrographic Conference (CHC-2014), St. John’s, Canada.

– CHS (2013). Canadian Hydrographic Service STANDARDS FOR HYDROGRAPHIC SURVEYS. 2nd Edition. 2013. 17 p.

– Calder, B. R. and Mayer, L. A. (2003). Automatic processing of high-rate, high-density multibeam echosounder data. Geochemistry, Geophysics, Geosystems: G3, Washington, USA: American Geophysical Union & Geochemical Society, v. 4, n. 6, 2003. pp. 1-22.

– CARIS. (2017). CARIS HIPS and SIPS 10.2 Reference Guide, 494p.

– Denbigh, P. N. (1989). Swath Bathymetry: Principles of Operation and an Analysis of Errors. IEEE Journal of Oceanic Engineering, v. 14, n. 4, p. 289–298.

– DHN. (2017). NORMAM-25 rev.2., viewed 20 May 2020, https://www.marinha.mil.br/dhn/?q=pt-br/node/266.

– Dodd, D. (2013). Uncertainty Evaluation of the EdgeTech 4600 Swath Bathymetry System. In: US Hydrographic Conference (U.S.Hydro-2013), New Orleans, United States.

– Gostnell, C. (2005). Efficacy of an interferometric sonar for hydrographic surveying: Do interferometers warrant an in-depth examination? – A Preliminary Study. In: US Hydrographic Conference (U.S.Hydro-2005), San Diego, United States.

– Gostnell, C.; Yoos, J. and Brodet, S. (2006). NOAA test and evaluation of interferometric sonar technology. In: Canadian Hydrographic Conference (CHC-2006), Halifax, Canada.

– Jerram, K. and Schmidt, V. (2015a). Object Detection with Phase-Measuring Bathymetric Sidescan Sonar Depth Data. (White Paper – CCOM), viewed 04 January 2018, http://sandy.ccom.unh.edu/publications/library/Object_detection_with_PMBS.pdf.

– Jerram, K. and Schmidt, V. (2015b). Storm Response Surveying with Phase- Measuring Bathymetric Sidescan Sonar. (White Paper – CCOM), viewed 04 January 2018, <http://sandy.ccom.unh.edu/publications/library/Response_surveying_with_PMBS.pdf>.

– Kjerfve, B. ; Ribeiro C. H. ; Dias G. T. M. ; Filippo A. M. and Quaresma V. S. (1997). Oceanographic Characteristics of an Impacted Coastal Bay: Baía de Guanabara, Rio de Janeiro, Brazil. Continental Shelf Research, Great Britain: Elsevier Science Ltd, v. 17, n. 13, p. 1609–1643.

– KONGSBERG MARITIME. (2014). GeoSwath Plus wide swath bathymetry system – Instruction Manual. 639p.

– L-3 COMMUNICATIONS SEABEAM INSTRUMENTS. (2000). Multibeam Sonar Theory of Operation. East Walpole, United States.

– LINZ. (2010). Contract Specifications for Hydrographic Surveys, Version 1.2, 110p.

– Lurton, X. (2000). Swath Bathymetry Using Phase Difference: Theoretical Analysis of Acoustical Measurement Precision. IEEE Journal of Oceanic Engineering, v. 25, n. 3, p. 351–363, viewed 23 January 2018: <https://ieeexplore.ieee.org/document/855385/>.

– Ma, K.; Xu, W. and Xu, J. (2016). The Comparison between Traditional and Interferometric Multibeam Systems. In: 6th International Conference on Sensor Network and Computer Engineering (ICSNCE 2016), Xi’an, China.

– Madricardo, F.; Foglini, F. and Tonielli, R. (2011). Bathymetric Surveys in Super-Shallow Water: Case Study of The Lagoon of Venice, Italy. In: 4th International Conference & Exhibition on Underwater Acoustic Measurements: Technologies and Results. Kos Island, Greece.

– NATIONAL OCEANIC AND ATMOSPHERIC ADMINISTRATION (NOAA). (2017). Hydrographic Surveys Specifications and Deliverables, viewed 16 December 2017, http://www.nauticalcharts.noaa.gov/hsd/specs/specs.htm.

– IHO. (2008). IHO Standards for Hydrographic Surveys (S-44). 5th ed. Monaco.

– IHO. (2005). Manual on Hydrography (C-13). 1st ed. corrections to 2011. Monaco.

– Pereira, D. L. S. (2016). Incertezas e Controle de Qualidade na Medição de Profundidades por Meio de Ecobatímetros Multifeixe. Anais Hidrográficos, Niterói: DHN, v. 73. p. 53-65.

– Saucan, A. A.; Sintes C.; Chonavel T. and Caillec J-M L. (2015). Model-Based Adaptive 3D Sonar Reconstruction in Reverberating Environments. IEEE Transactions on Image Processing, IEEE Signal Processing Society, v. 24, n. 10, p. 2928–2940.

– Simmons, S. M. ; Parson, D. R. ; Best, J. L. ; Oberg, K, A. ; Czuba, J, A. and Keevil, G.M. (2017). An evaluation of the use of a multibeam echo-sounder for observations of suspended sediment. Applied Acoustics, Great Britain: Elsevier Science Ltd, 2017. v. 126, p. 81–90.

– UNITED STATES ARMY CORPS OF ENGINEERS (USACE). (2013). EM 1110-2-1003 -Engineering and Design Hydrographic Surveying. Washington, Un ited States, 700p.

– Wilby, A. D. (1999). The Advantages, Challenges and Practical Implementation of an Interferometric Swath Bathymetry System. In: OCEANS’99 MTS/IEEE, Seatle, United States.

– Woock, P. and Frey, C. (2010). Deep-sea AUV navigation using side-scan sonar images and SLAM. In: OCEANS’10 IEEE, Sidney, Australia, viewed 16May 2018, https://ieeexplore.ieee.org/document/5603528/.

8. Authors Biography

Lieutenant Commander Vitor Bravo Pimentel, Head of the Operation Department of the Brazilian Antartic Program Polar Vessel “Almirante Maximiano”. B.Sc. in Naval Science (Naval College, 2006), he completed the Hydrographic Course IHO/FIG/ICA Category A (Directorated of Hydrography and Navigation – DHN, 2009) and a M.Sc. in Ocean and Earth Dynamics (Fluminense Federal University – UFF, 2018). He is also a Ph.D. student of Hydrography at UFF with the researching theme “Nautical Charts Update Through Phase Measurement Bathymetric Sonar”. Email: vitor.pimentel@marinha.mil.br

Lieutenant Commander Christopher Florentino, Head of Cartographic Products Section at the Directorate of Hydrography and Navigation (DHN). B.Sc. in Cartographic Engineer (Faculty of Science and Technology at São Paulo State University, 2005) joined the Brazilian Navy (Corps of Engineers) in 2006, and earned a M.Sc. in Ocean and Earth Dynamics (Fluminense Federal University – UFF, 2018). As a current Ph.D. student of Hydrography at UFF, he researches the “Application of multiple analyses focus on quality assessment of national bathymetric surfaces in order to support the Marine Spatial Data Infrastructures (MSDI) initiatives”. Email: christopher@marinha.mil.br

Arthur Ayres Neto, Associate Professor in Marine Geology and Geophysics with emphasis on morphological and geological characterization of the seafloor with acoustic techniques at the Fluminense Federal University (UFF) since 2006. B.Sc in Geology (Federal University of Rio de Janeiro – UFRJ, 1988); M.Sc. in Marine Geology and Geophysics (UFF, 1993); Ph.D in Marine Geophysics (Kiel University, Germany, 1999). Worked for 19 years on the offshore industry with hydrography in the development of subsea engineering projects. Email: aayres@id.uff.br