Automated, objective texture segmentation of multibeam echosounder data - Seafloor survey and substrate maps from James Island to Ozette Lake, Washington Outer Coast

Without knowledge of basic seafloor characteristics, the ability to address any number of critical marine and/or coastal management issues is diminished. For example, management and conservation of essential fish habitat (EFH), a requirement mandated by federally guided fishery management plans (FMPs), requires among other things a description of habitats for federally managed species. Although the list of attributes important to habitat are numerous, the ability to efficiently and effectively describe many, and especially at the scales required, does not exist with the tools currently available. However, several characteristics of seafloor morphology are readily obtainable at multiple scales and can serve as useful descriptors of habitat. Recent advancements in acoustic technology, such as multibeam echosounding (MBES), can provide remote indication of surficial sediment properties such as texture, hardness, or roughness, and further permit highly detailed renderings of seafloor morphology. With acoustic-based surveys providing a relatively efficient method for data acquisition, there exists a need for efficient and reproducible automated segmentation routines to process the data. Using MBES data collected by the Olympic Coast National Marine Sanctuary (OCNMS), and through a contracted seafloor survey, we expanded on the techniques of Cutter et al. (2003) to describe an objective repeatable process that uses parameterized local Fourier histogram (LFH) texture features to automate segmentation of surficial sediments from acoustic imagery using a maximum likelihood decision rule. Sonar signatures and classification performance were evaluated using video imagery obtained from a towed camera sled. Segmented raster images were converted to polygon features and attributed using a hierarchical deep-water marine benthic classification scheme (Greene et al. 1999) for use in a geographical information system (GIS). (PDF contains 41 pages.)

Benthic habitats of Buck Island Reef National Monument

NOAA’s Center for Coastal Monitoring and Assessment’s Biogeography Branch has mapped and characterized large portions of the coral reef ecosystems inside the U.S. coastal and territorial waters, including the U.S. Caribbean. The complementary protocols used in these efforts have enabled scientists and managers to quantitatively compare different marine ecosystems in tropical U.S. waters. The Biogeography Branch used these same general protocols to generate three seamless habitat maps of the Bank/Shelf (i.e., from 0 ≤50 meters) and the Bank/Shelf Escarpment (i.e., from 50 ≤1,000 meters and from 1,000 ≤ 1,830 meters) inside Buck Island Reef National Monument (BIRNM). While this mapping effort marks the fourth time that the shallow-water habitats of BIRNM have been mapped, it is the first time habitats deeper than 30 meters (m) have been characterized. Consequently, this habitat map provides information on the distribution of mesophotic and deep-water coral reef ecosystems and serves as a spatial baseline for monitoring change in the Monument. A benthic habitat map was developed for approximately 74.3 square kilometers or 98% of the BIRNM using a combination of semi-automated and manual classification methods. The remaining 2% was not mapped due to lack of imagery in the western part of the Monument at depths ranging from 1,000 to 1,400 meters. Habitats were interpreted from orthophotographs, LiDAR (Light Detection and Ranging) imagery and four different types of MBES (Multibeam Echosounder) imagery. Three minimum mapping units (MMUs) (100, 1,000 and 5,000 square meters) were used because of the wide range of depths present in the Monument. The majority of the area that was characterized was deeper than 30 m on the Bank/Shelf Escarpment. This escarpment area was dominated by uncolonized sand which transitioned to mud as depth increased. Bedrock was exposed in some areas of the escarpment, where steep slopes prevented sediment deposition. Mesophotic corals were seen in the underwater video, but were too sparsely distributed to be reliably mapped from the source imagery. Habitats on the Bank/Shelf were much more variable than those seen on the Bank/Shelf Escarpment. The majority of this shelf area was comprised of coral reef and hardbottom habitat dominated by various forms of turf, fleshy, coralline or filamentous algae. Even though algae was the dominant biological cover type, nearly a quarter (24.3%) of the Monument’s Bank/Shelf benthos hosted a cover of 10%-<50% live coral. In total, 198 unique combinations of habitat classes describing the geography, geology and biology of the sea-floor were identified from the three types of imagery listed above. No thematic accuracy assessment was conducted for areas deeper than about 50 meters, most of which was located in the Bank/Shelf Escarpment. The thematic accuracy of classes in waters shallower than approximately 50 meters ranged from 81.4% to 94.4%. These thematic accuracies are similar to those reported for other NOAA benthic habitat mapping efforts in St. John (>80%), the Main Eight Hawaiian Islands (>84.0%) and the Republic of Palau (>80.0%). These digital maps products can be used with confidence by scientists and resource managers for a multitude of different applications, including structuring monitoring programs, supporting management decisions, and establishing and managing marine conservation areas. The final deliverables for this project, including the benthic habitat maps, source imagery and in situ field data, are available to the public on a NOAA Biogeography Branch website (http://ccma.nos.noaa.gov/ecosystems/coralreef/stcroix.aspx) and through an interactive, web-based map application (http://ccma.nos.noaa.gov/explorer/biomapper/biomapper.html?id=BUIS). This report documents the process and methods used to create the shallow to deep-water benthic habitat maps for BIRNM. Chapter 1 provides a short introduction to BIRNM, including its history, marine life and ongoing research activities. Chapter 2 describes the benthic habitat classification scheme used to partition the different habitats into ecologically relevant groups. Chapter 3 explains the steps required to create a benthic habitat map using a combination of semi-automated and visual classification techniques. Chapter 4 details the steps used in the accuracy assessment and reports on the thematic accuracy of the final shallow-water map. Chapter 5 summarizes the type and abundance of each habitat class found inside BIRNM, how these habitats compare to past habitat maps and outlines how these new habitat maps may be used to inform future management activities.

Hydro-acoustic remote sensing of benthic biological communities on the shallow South East Australian continental shelf

Information regarding the composition and extent of benthic habitats on the South East Australian continental shelf is limited. In this habitat mapping study, multibeam echosounder (MBES) data are integrated with precisely geo-referenced video ground-truth data to quantify benthic biotic communities at Cape Nelson, Victoria, Australia. Using an automated decision tree classification approach, 5 representative biotic groups defined from video analysis were related to hydro-acoustically derived variables in the Cape Nelson survey area. Using a combination of multibeam bathymetry, backscatter and derivative products produced highest overall accuracy (87%) and kappa statistic (0.83). This study demonstrates that decision tree classifiers are capable of integrating variable data types for mapping distributions of benthic biological assemblages, which are important in maintaining biodiversity and other system services in the marine environment.

Comparison of automated classification techniques for predicting benthic biological communities using hydroacoustics and video observations

The effective management of our marine ecosystems requires the capability to identify, characterise and predict the distribution of benthic biological communities within the overall seascape architecture. The rapid expansion of seabed mapping studies has seen an increase in the application of automated classification techniques to efficiently map benthic habitats, and the need of techniques to assess confidence of model outputs. We use towed video observations and 11 seafloor complexity variables derived from multibeam echosounder (MBES) bathymetry and backscatter to predict the distribution of 8 dominant benthic biological communities in a 54 km² site, off the central coast of Victoria, Australia. The same training and evaluation datasets were used to compare the accuracies of a Maximum Likelihood Classifier (MLC) and two new generation decision tree methods, QUEST (Quick Unbiased Efficient Statistical Tree) and CRUISE (Classification Rule with Unbiased Interaction Selection and Estimation), for predicting dominant biological communities. The QUEST classifier produced significantly better results than CRUISE and MLC model runs, with an overall accuracy of 80% (Kappa 0.75). We found that the level of accuracy with the size of training set varies for different algorithms. The QUEST results generally increased in a linear fashion, CRUISE performed well with smaller training data sets, and MLC performed least favourably overall, generating anomalous results with changes to training size. We also demonstrate how predicted habitat maps can provide insights into habitat spatial complexity on the continental shelf. Significant variation between patch-size and habitat types and significant correlations between patch size and depth were also observed.

Integrating time-series hydroacoustics and video observations for detecting changes in benthic habitats

The ability to quantify change in marine benthic habitats must be considered a key goal of marine habitat mapping activities. Changes in distribution of distinct suites of benthic biological species may occur as a result of natural or human induced processes and these processes may operate at a range of temporal and spatial scales. It is important to understand natural small scale inter-annual patterns of change in order to separate these signals from potential patterns of longer term change. Work to describe these processes of change from an acoustic remote sensing stand point has thus far been limited due to the relatively recent availability of full coverage swath acoustic datasets and cost pressures associated with multiple surveys of the same area. This paper describes the use of landscape transition analysis as a means to differentiate seemingly random patterns of habitat change from systematic signals of habitat transition at a shallow (10–50 m depth) 18 km2 study area on the temperate Australian continental shelf between the years 2006 and 2007. Supervised classifications for each year were accomplished using independently collected high resolution (3 m cell-size) multibeam echosounder (MBES) and video-derived reference data. Of the 4 representative biotic classes considered, signals of directional systematic changes were observed to occur between a shallow kelp dominated class, a deep sessile invertebrate dominated class and a mixed class of kelp and sessile invertebrates. These signals of change are interpreted as inter-annual variation in the density and depth related extent of canopy forming kelp species at the site, a phenomenon reported in smaller scale temporal studies of the same species. The methods applied in this study provide a detailed analysis of the various components of the traditional change detection cross tabulation matrix allowing identification of the strongest signals of systematic habitat transitions across broad geographical regions. Identifying clear patterns of habitat change is an important first step in linking these patterns to the processes that drive them.

Evaluation of four supervised learning methods for benthic habitat mapping using backscatter from multi-beam sonar

An understanding of the distribution and extent of marine habitats is essential for the implementation of ecosystem-based management strategies. Historically this had been difficult in marine environments until the advancement of acoustic sensors. This study demonstrates the applicability of supervised learning techniques for benthic habitat characterization using angular backscatter response data. With the advancement of multibeam echo-sounder (MBES) technology, full coverage datasets of physical structure over vast regions of the seafloor are now achievable. Supervised learning methods typically applied to terrestrial remote sensing provide a cost-effective approach for habitat characterization in marine systems. However the comparison of the relative performance of different classifiers using acoustic data is limited. Characterization of acoustic backscatter data from MBES using four different supervised learning methods to generate benthic habitat maps is presented. Maximum Likelihood Classifier (MLC), Quick, Unbiased, Efficient Statistical Tree (QUEST), Random Forest (RF) and Support Vector Machine (SVM) were evaluated to classify angular backscatter response into habitat classes using training data acquired from underwater video observations. Results for biota classifications indicated that SVM and RF produced the highest accuracies, followed by QUEST and MLC, respectively. The most important backscatter data were from the moderate incidence angles between 30° and 50°. This study presents initial results for understanding how acoustic backscatter from MBES can be optimized for the characterization of marine benthic biological habitats. © 2012 by the authors.

Quantitative experimental comparison of single-beam, sidescan, and multibeam benthic habitat maps

Map comparison is a relatively uncommon practice in acoustic seabed classification to date, contrary to the field of land remote sensing, where it has been developed extensively over recent decades. The aim here is to illustrate the benefits of map comparison in the underwater realm with a case study of three maps independently describing the seabed habitats of the Te Matuku Marine Reserve (Hauraki Gulf, New Zealand). The maps are obtained from a QTC View classification of a single-beam echosounder (SBES) dataset, manual segmentation of a sidescan sonar (SSS) mosaic, and automatic classification of a backscatter dataset from a multibeam echosounder (MBES). The maps are compared using pixel-to-pixel similarity measures derived from the literature in land remote sensing. All measures agree in presenting the MBES and SSS maps as the most similar, and the SBES and SSS maps as the least similar. The results are discussed with reference to the potential of MBES backscatter as an alternative to SSS mosaic for imagery segmentation and to the potential of joint SBES–SSS survey for improved habitat mapping. Other applications of map-similarity measures in acoustic classification of the seabed are suggested.

Comparison of a self-processed EM3000 multibeam echosounder dataset with a QTC view habitat mapping and a sidescan sonar imagery, Tamaki Strait, New Zealand

A methodology for automatically processing the data files from an EM3000 multibeam echosounder (Kongsberg Maritime, 300 kHz) is presented. Written in MatLab, it includes data extraction, bathymetry processing, computation of seafloor local slope, and a simple correction of the backscatter along-track banding effect. The success of the latter is dependent on operational restrictions, which are also detailed. This processing is applied to a dataset acquired in 2007 in the Tamaki Strait, New Zealand. The resulting maps are compared with a habitat classification obtained with the acoustic ground-discrimination software QTC View linked to a 200-kHz single-beam echosounder and to the imagery from a 100-kHz sidescan sonar survey, both performed in 2002. The multibeam backscatter map was found to be very similar to the sidescan imagery, quite correlated to the QTC View map on one site but mainly uncorrelated on another site. Hypotheses to explain these results are formulated and discussed. The maps and the comparison to prior surveys are used to draw conclusions on the quality of the code for further research on multibeam benthic habitat mapping.

Integrating multibeam backscatter angular response, mosaic and bathymetry data for benthic habitat mapping

Multibeam echosounders (MBES) are increasingly becoming the tool of choice for marine habitat mapping applications. In turn, the rapid expansion of habitat mapping studies has resulted in a need for automated classification techniques to efficiently map benthic habitats, assess confidence in model outputs, and evaluate the importance of variables driving the patterns observed. The benthic habitat characterisation process often involves the analysis of MBES bathymetry, backscatter mosaic or angular response with observation data providing ground truth. However, studies that make use of the full range of MBES outputs within a single classification process are limited. We present an approach that integrates backscatter angular response with MBES bathymetry, backscatter mosaic and their derivatives in a classification process using a Random Forests (RF) machine-learning algorithm to predict the distribution of benthic biological habitats. This approach includes a method of deriving statistical features from backscatter angular response curves created from MBES data collated within homogeneous regions of a backscatter mosaic. Using the RF algorithm we assess the relative importance of each variable in order to optimise the classification process and simplify models applied. The results showed that the inclusion of the angular response features in the classification process improved the accuracy of the final habitat maps from 88.5% to 93.6%. The RF algorithm identified bathymetry and the angular response mean as the two most important predictors. However, the highest classification rates were only obtained after incorporating additional features derived from bathymetry and the backscatter mosaic. The angular response features were found to be more important to the classification process compared to the backscatter mosaic features. This analysis indicates that integrating angular response information with bathymetry and the backscatter mosaic, along with their derivatives, constitutes an important improvement for studying the distribution of benthic habitats, which is necessary for effective marine spatial planning and resource management.

Accounting for uncertainty in volumes of seabed change measured with repeat multibeam sonar surveys

Seafloors of unconsolidated sediment are highly dynamic features; eroding or accumulating under the action of tides, waves and currents. Assessing which areas of the seafloor experienced change and measuring the corresponding volumes involved provide insights into these important active sedimentation processes. Computing the difference between Digital Elevation Models (DEMs) obtained from repeat Multibeam Echosounders (MBES) surveys has become a common technique to identify these areas, but the uncertainty in these datasets considerably affects the estimation of the volumes displaced. The two main techniques used to take into account uncertainty in volume estimations are the limitation of calculations to areas experiencing a change in depth beyond a chosen threshold, and the computation of volumetric confidence intervals. However, these techniques are still in their infancy and, as a result, are often crude, seldom used or poorly understood. In this article, we explored a number of possible methodological advances to address this issue, including: (1) using the uncertainty information provided by the MBES data processing algorithm CUBE, (2) adapting fluvial geomorphology techniques for volume calculations using spatially variable thresholds and (3) volumetric histograms. The nearshore seabed off Warrnambool harbour - located in the highly energetic southwest Victorian coast, Australia - was used as a test site. Four consecutive MBES surveys were carried out over a four-months period. The difference between consecutive DEMs revealed an area near the beach experiencing large sediment transfers - mostly erosion - and an area of reef experiencing increasing deposition from the advance of a nearby sediment sheet. The volumes of sediment displaced in these two areas were calculated using the techniques described above, both traditionally and using the suggested improvements. We compared the results and discussed the applicability of the new methodological improvements. We found that the spatially variable uncertainty derived from the CUBE algorithm provided the best results (i.e. smaller confidence intervals), but that similar results can be obtained using as a fixed uncertainty value derived from a reference area under a number of operational conditions.

Mappatura degli habitat e degli oggetti antropici della laguna di Venezia

L’utilizzo del Multibeam Echo sounder (MBES) in ambienti di transizione poco profondi, con condizioni ambientali complesse come la laguna di Venezia, è ancora in fase di studio e i dati biologici e sedimentologici inerenti ai canali della laguna di Venezia sono attualmente scarsi e datati in letteratura. Questo studio ha lo scopo di mappare gli habitat e gli oggetti antropici di un canale della laguna di Venezia in un intervallo di profondità tra 0.3 e 20 m (Canale San Felice) analizzando i dati batimetrici e di riflettività (backscatter) acquisiti da ISMAR-Venezia nell’ambito del progetto RITMARE. A tale scopo il fondale del canale San Felice (Venezia) è stato caratterizzato dal punto di vista geomorfologico, sedimentologico e biologico; descrivendo anche l’eventuale presenza di oggetti antropici. L’ecoscandaglio utilizzato è il Kongsberg EM2040 Dual-Compact Multibeam in grado di emettere 800 beam (400 per trasduttore) ad una frequenza massima di 400kHZ e ci ha consentito di ricavare ottimi risultati, nonostante le particolari caratteristiche degli ambienti lagunari. I dati acquisiti sono stati processati tramite il software CARIS Hydrographic information processing system (Hips) & Sips, attraverso cui è possibile applicare le correzioni di marea e velocità del suono e migliorare la qualità dei dati grezzi ricavati da MBES. I dati sono stati quindi convertiti in ESRI Grid, formato compatibile con il software ArcGIS 10.2.1 (2013) che abbiamo impiegato per le interpretazioni e per la produzione delle mappe. Tecniche di ground-truthing, basate su riprese video e prelievi di sedimento (benna Van Veen 7l), sono state utilizzate per validare il backscatter, dimostrandosi molto efficaci e soddisfacenti per poter descrivere i fondali dal punto di vista biologico e del substrato e quindi degli habitat del canale lagunare. Tutte le informazioni raccolte durante questo studio sono state organizzate all’interno di un geodatabase, realizzato per i dati relativi alla laguna di Venezia.

Bathymetry of Lake Constance – a High-Resolution Survey in a Large, Deep Lake

Within the project “Tiefenschärfe – hochauflösende Vermes- sung Bodensee” a high-resolution seamless terrain model is created using airborne topobathymetric laserscanning and multibeam echosounder (MBES) techniques. The project visu- alizes the enormous wealth of features of underwater land- scapes of lakes. The combination of hydroacoustic (multibeam echosounder) and laser-optic (topobathymetric laserscan- ning) methods was used for the first time in a freshwater body of this size. Opportunities, limitations and restrictions of these high-resolution methods are presented.

Penetration velocity and deceleration of penetrometer stations N1 to N97 during cruise Senckenberg_11_2008

In-situ geotechnical measurements of surface sediments were carried out along large subaqueous dunes in the Knudedyb tidal inlet channel in the Danish Wadden Sea using a small free-falling penetrometer. Vertical profiles showed a typical stratification pattern with a resolution of ~1 cm depicting a thin surface layer of low sediment strength and a stiffer substratum below (quasi-static bearing capacity equivalent: 1-3 kPa in the top layer, 20-140 kPa in the underlying sediment; thickness of the top layer ca. 5-8 cm). Observed variations in the thickness and strength of the surface layer during a tidal cycle were compared to mean current velocities (measured using an acoustic Doppler current profiler, ADCP), high-resolution bathymetry (based on multibeam echo sounding, MBES) and qualitative estimates of suspended sediment distributions in the water column (estimated from ADCP backscatter intensity). The results revealed an ebb dominance in sediment remobilization, and a general accretion of the bed towards low water. A loose top layer occurred throughout the tidal cycle, likely influenced by bedload transport and small events of suspended sediment resettlement (thickness: 6 +-2 cm). Furthermore, this layer showed a significant increase in thickness (e.g. from 8 cm to 16 cm) related to periods of overall deposition. These findings imply that dynamic penetrometers can conveniently serve to (1) quantify potentially mobile sediments by determining the thickness of a loose sediment surface layer, (2) unravel sediment strength development in potentially mobile sediments and (3) identify sediment accumulation. Such data are an important complement and add a new geotechnical perspective during investigations of sediment remobilization processes in highly dynamic coastal environments.