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.

Seafloor mapping data, Victoria

The collection includes three major data types

- Seafloor structure information: Collation of seafloor structure information (bathymetry and softness hardness) collected using multibeam sonar systems as part of the Victorian Marine Habitat Mapping Project and bathymetric light detection and ranging data (LiDAR) collected as part of the future coats program. The geographic is Victorian State waters (~1400km2 multibeam sonar & ~4000 km2 for LiDAR as of 2012

- Video information: ~800 linear kilometres of precisely georeferenced towed video data from depths of 10-100m in Victorian coastal waters. 300 baited stereo video drops in Hopkins and Discovery Bay site locations. All video has been classified to the lowest taxonomic resolution possible with the video systems employed and converted to spatial data layers in the GIS environment.

- Habitat maps: Predicted habitat models integrating seafloor structure and video information developed using classification algorithms and remote sensing technologies.

Corner inlet and Nooramunga habitat mapping project

Deakin University and the University of Tasmania were commissioned by Parks Victoria (PV) to create two updated habitat maps for areas within the Corner Inlet and Nooramunga Marine and Coastal Park and Ramsar area. The team obtained a ground-truth data set using in situ video and still photographs. This dataset was used to develop and assess predictive models of benthic marine habitat distributions incorporating data from both ALOS (Advanced Land Observation Satellite) imagery atmospherically corrected by CSIRO and LiDAR (Light Detection and Ranging) bathymetry. This report describes the results of the mapping effort as well as the methodology used to produce these habitat maps.

Overall accuracies of habitat classifications were good, returning overall accuracies >73 % and kappa values > 0.62 for both study localities. Habitats predicted with highest accuracies included Zosteraceae in Nooramunga (91 %), reef in Corner Inlet (80 %), and bare sediment (no-visible macrobiota/no-visible seagrass classes; both > 76 %). The majority of classification errors were due to the misclassification of areas of sparse seagrass as bare sediment. For the Corner Inlet study locality the no-visible macrobiota (10,698 ha), Posidonia (4,608 ha) and Zosteraceae (4,229 ha) habitat classes covered the most area. In Nooramunga no-visible seagrass (5,538 ha), Zosteraceae (4,060 ha) and wet saltmarsh (1,562 ha) habitat classes were most dominant.

In addition to the commissioned work preliminary change detection analyses were undertaken as part of this project. These analyses indicated shifts in habitat extents in both study localities since the late 1990s/2000. In particular, a post-classification analysis highlighted that there were considerable increases in seagrass habitat (primarily Zosteraceae) throughout the littoral zones and river/creek mouths of both study localities. Further, the numerous channel systems remained stable and were free of seagrass at both times. A substantial net loss of Posidonia in the Corner Inlet locality is likely but requires further investigation due to potential misclassifications between habitats in both the 1998 map (Roob et al. 1998) and the current mapping. While the unsupervised Independent Components Analysis (ICA) change detection technique indicated some changes in habitat extent and distribution, considerable areas of habitat change observed in the post-classification approach are questionable, and may reflect misclassifications rather than real change. A particular example of this is an apparent large decrease in Zosteraceae and increase in Posidonia being related to the classification of Posidonia beds as Zosteraceae in the 1998 mapping. Despite this, we believe that changes indicated by both the ICA and post-classification approaches have a high likelihood of being ‘actual’ change. A pattern of gains and losses of Zosteraceae in the region north of Stockyard channel is an example of this. Further analyses and refinements of approaches in change detection analyses such as would improve confidence in the location and extent of habitat changes over this time period.

This work has been successful in providing new baseline maps using a repeatable method meaning that any future changes in intertidal and shallow water marine habitats may be assessed in a consistent way with quantitative error assessments. In wider use, these maps should also allow improved conservation planning, advance fisheries and catchment management, and progress infrastructure planning to limit impacts on the Inlet environment.

Assessing the use of milk fatty acids to infer the diet of the Australian sea lion (Neophoca cinerea)

Information on the diet of threatened species is important in devising appropriate management plans to ensure their conservation. The Australian sea lion (Neophoca cinerea) is Australia’s only endemic and globally one of the least numerous pinniped species. However, dietary information is currently limited because of the difficulty in using traditional methods (identification of prey hard parts from scats, regurgitates and stomach samples) to reliably provide dietary information. The present study assessed the use of fatty acid (FA) analysis to infer diet using milk samples collected from 11 satellite tracked Australian sea lions from Olive Island, South Australia. Satellite tracking revealed that females foraged in two distinct regions; ‘inshore’ regions characterised by shallow bathymetry (10.7 ± 4.8 m) and ‘offshore’ regions characterised by comparatively deep bathymetry (60.5 ± 13.4 m). Milk FA analysis indicated significant differences in the FA composition between females that foraged inshore compared with those that foraged offshore. The greatest differences in relative levels of individual FAs between the inshore and offshore groups were for 22 : 6n-3 (6.5 ± 1.2% compared with 16.5 ± 1.9% respectively), 20 : 4n-6 (6.1 ± 0.7 compared with 2.5 ± 0.7 respectively) and 22 : 4n-6 (2.4 ± 0.2% compared with 0.8 ± 0.2% respectively). Using discriminant scores, crustacean, cephalopod, fish and shark-dominated diets were differentiated. The discriminant scores from Australian sea lions that foraged inshore indicated a mixed fish and shark diet, whereas discriminant scores from Australian sea lions that foraged offshore indicated a fish-dominated diet, although results must be interpreted with caution due to the assumptions associated with the prey FA dataset. FA analysis in combination with satellite tracking proved to be a powerful tool for assessing broad-scale spatial dietary patterns.

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.

Corner Inlet and Nooramunga Habitat Mapping Project – Part 2 : High resolution mapping and change detection analysis using 2011 Rapid Eye Satellite Image compared to 2009 ALOS Satellite Image

Building on a habitat mapping project completed in 2011, Deakin University was commissioned by Parks Victoria (PV) to apply the same methodology and ground-truth data to a second, more recent and higher resolution satellite image to create habitat maps for areas within the Corner Inlet and Nooramunga Marine and Coastal Park and Ramsar area. A ground-truth data set using in situ video and still photographs was used to develop and assess predictive models of benthic marine habitat distributions incorporating data from both RapidEye satellite imagery (corrected for atmospheric and water column effects by CSIRO) and LiDAR (Light Detection and Ranging) bathymetry. This report describes the results of the mapping effort as well as the methodology used to produce these habitat maps.

Overall accuracies of habitat classifications were good, with error rates similar to or better than the earlier classification (>73 % and kappa values > 0.58 for both study localities). The RapidEye classification failed to accurately detect Pyura and reef habitat classes at the Corner Inlet locality, possibly due to differences in spectral frequencies. For comparison, these categories were combined into a ‘non-seagrass’ category, similar to the one used at the Nooramunga locality in the original classification. Habitats predicted with highest accuracies differed from the earlier classification and were Posidonia in Corner Inlet (89%), and bare sediment (no-visible seagrass class) in Nooramunga (90%). In the Corner Inlet locality reef and Pyura habitat categories were not distinguishable in the repeated classification and so were combined with bare sediments. The majority of remaining classification errors were due to the misclassification of Zosteraceae as bare sediment and vice versa. Dominant habitats were the same as those from the 2011 classification with some differences in extent. For the Corner Inlet study locality the no-visible seagrass category remained the most extensive (9059 ha), followed by Posidonia (5,513 ha) and Zosteraceae (5,504 ha). In Nooramunga no-visible seagrass (6,294 ha), Zosteraceae (3,122 ha) and wet saltmarsh (1,562 ha) habitat classes were most dominant.

Change detection analyses between the 2009 and 2011 imagery were undertaken as part of this project, following the analyses presented in Monk et al. (2011) and incorporating error estimates from both classifications. These analyses indicated some shifts in classification between Posidonia and Zosteraceae as well as a general reduction in the area of Zosteraceae. Issues with classification of mixed beds were apparent, particularly in the main Posidonia bed at Nooramunga where a mosaic of Zosteraceae and Posidonia was seen that was not evident in the ALOS classification. Results of a reanalysis of the 1998-2009 change detection illustrating effects of binning of mixed beds is also provided as an appendix.

This work has been successful in providing baseline maps at an improved level of detail using a repeatable method meaning that any future changes in intertidal and shallow water marine habitats may be assessed in a consistent way with quantitative error assessments. In wider use, these maps should also allow improved conservation planning, advance fisheries and catchment management, and progress infrastructure planning to limit impacts on the Inlet environment.

Yaringa and French Island Marine National Park habitat mapping

Deakin University and the Department of Primary Industries were commissioned by ParksVictoria (PV) to create two updated habitat maps for Yaringa and French Island MarineNational Parks. The team obtained a ground-truth data set using in situ video and still photographs. This dataset was used to develop and assess predictive models of benthic marine habitat distributions incorporating data from World-View-2 imagery atmospherically corrected by CSIRO and LiDAR (Light Detection and Ranging) bathymetry. In addition, the team applied an unsupervised classification approach to an aerial photograph to assess the differences between the two remote sensors. This report describes the results of the mapping as well as the methodology used to produce these habitat maps.This study has provided mapping of intertidal and subtidal habitats of Yaringa and FrenchIsland MNPs at a 2 m resolution with fair to good accuracies (Kappa 0.40-0.75). These were combined with mangrove and saltmarsh habitats recently mapped by Boon et al. (2011) to provide compete-coverage habitat maps of Yaringa and French Island MNPs.The mapping showed that Yaringa MNP was dominated by mangroves, wet saltmarsh and dense Zostereaceae, covering 33%, 29% and 19%, respectively. Similarly, intertidalvegetation and subtidal vegetation (dominated by Zosteraceae) covered 26% and 25% ofFrench Island MNP. However, as a result of turbidity and missing satellite imagery 27% ofFrench Island MNP remains unmapped.The coupling of WV-2 and LiDAR reduced potential artefacts (e.g. sun glint causing whiteand black pixels known as the “salt and pepper effect”). The satellite classification appeared to provide better results than the aerial photography classification. However, since there is a two-year difference between the capture of the aerial photography and the collection of the ground-truth data this comparison is potentially temporally confounded. It must also be noted that there are differences in costs of the data,the spatial resolution between the two datasets (i.e. WV-2 = 2 m and the Aerial = 0.5 m) and the amount spectral information contained in the data (i.e. WV-2 = 8 bands and the aerial = 4 bands), which may ultimately determine its utility for a particular project.The spatial assessment using FRAGSTATs of habitat patches within Yaringa MNP provides a viable and cost effect way to assess habitat condition (i.e. shape, size and arrangement).This spatial assessment determined that dense Zosteraceae and NVSG habitat classeswere generally larger in patch size and continuity than the medium/sparse Zosteraceaehabitat. The application spatial techniques to time-series mapping may provide a way toremotely monitor the change in the spatial characteristics of marine habitats.This work was successful in providing new baseline habitat maps using a repeatable method meaning that any future changes in intertidal and shallow water marine habitats may be assessed in a consistent way with quantitative error assessments. In wider use, these maps should also allow improved conservation planning, fisheries and catchment management, and contribute toward infrastructure planning to limit impacts on Western Port.

Habitat classification of temperate marine macroalgal communities using bathymetric LiDAR

Here, we evaluated the potential of using bathymetric Light Detection and Ranging (LiDAR) to characterise shallow water (<30 m) benthic habitats of high energy subtidal coastal environments. Habitat classification, quantifying benthic substrata and macroalgal communities, was achieved in this study with the application of LiDAR and underwater video groundtruth data using automated classification techniques. Bathymetry and reflectance datasets were used to produce secondary terrain derivative surfaces (e.g., rugosity, aspect) that were assumed to influence benthic patterns observed. An automated decision tree classification approach using the Quick Unbiased Efficient Statistical Tree (QUEST) was applied to produce substrata, biological and canopy structure habitat maps of the study area. Error assessment indicated that habitat maps produced were primarily accurate (>70%), with varying results for the classification of individual habitat classes; for instance, producer accuracy for mixed brown algae and sediment substrata, was 74% and 93%, respectively. LiDAR was also successful for differentiating canopy structure of macroalgae communities (i.e., canopy structure classification), such as canopy forming kelp versus erect fine branching algae. In conclusion, habitat characterisation using bathymetric LiDAR provides a unique potential to collect baseline information about biological assemblages and, hence, potential reef connectivity over large areas beyond the range of direct observation. This research contributes a new perspective for assessing the structure of subtidal coastal ecosystems, providing a novel tool for the research and management of such highly dynamic marine environments. © 2014 by the authors; licensee MDPI, Basel, Switzerland.

Assessment of habitat representation across a network of marine protected areas with implications for the spatial design of monitoring

Networks of marine protected areas (MPAs) are being adopted globally to protect ecosystems and supplement fisheries management. The state of California recently implemented a coast-wide network of MPAs, a statewide seafloor mapping program, and ecological characterizations of species and ecosystems targeted for protection by the network. The main goals of this study were to use these data to evaluate how well seafloor features, as proxies for habitats, are represented and replicated across an MPA network and how well ecological surveys representatively sampled fish habitats inside MPAs and adjacent reference sites. Seafloor data were classified into broad substrate categories (rock and sediment) and finer scale geomorphic classifications standard to marine classification schemes using surface analyses (slope, ruggedness, etc.) done on the digital elevation model derived from multibeam bathymetry data. These classifications were then used to evaluate the representation and replication of seafloor structure within the MPAs and across the ecological surveys. Both the broad substrate categories and the finer scale geomorphic features were proportionately represented for many of the classes with deviations of 1-6% and 0-7%, respectively. Within MPAs, however, representation of seafloor features differed markedly from original estimates, with differences ranging up to 28%. Seafloor structure in the biological monitoring design had mismatches between sampling in the MPAs and their corresponding reference sites and some seafloor structure classes were missed entirely. The geomorphic variables derived from multibeam bathymetry data for these analyses are known determinants of the distribution and abundance of marine species and for coastal marine biodiversity. Thus, analyses like those performed in this study can be a valuable initial method of evaluating and predicting the conservation value of MPAs across a regional network.

Exploring spatiotemporal trends in commercial fishing effort of an abalone fishing zone: a GIS-based hotspot model

Assessing patterns of fisheries activity at a scale related to resource exploitation has received particular attention in recent times. However, acquiring data about the distribution and spatiotemporal allocation of catch and fishing effort in small scale benthic fisheries remains challenging. Here, we used GIS-based spatio-statistical models to investigate the footprint of commercial diving events on blacklip abalone (Haliotis rubra) stocks along the south-west coast of Victoria, Australia from 2008 to 2011. Using abalone catch data matched with GPS location we found catch per unit of fishing effort (CPUE) was not uniformly spatially and temporally distributed across the study area. Spatial autocorrelation and hotspot analysis revealed significant spatiotemporal clusters of CPUE (with distance thresholds of 100's of meters) among years, indicating the presence of CPUE hotspots focused on specific reefs. Cumulative hotspot maps indicated that certain reef complexes were consistently targeted across years but with varying intensity, however often a relatively small proportion of the full reef extent was targeted. Integrating CPUE with remotely-sensed light detection and ranging (LiDAR) derived bathymetry data using generalized additive mixed model corroborated that fishing pressure primarily coincided with shallow, rugose and complex components of reef structures. This study demonstrates that a geospatial approach is efficient in detecting patterns and trends in commercial fishing effort and its association with seafloor characteristics.

Forests of the sea: predictive habitat modelling to assess the abundance of canopy forming kelp forests on temperate reefs

 Large brown seaweeds (kelps) form forests in temperate and boreal marine systems that serve as foundations to the structure and dynamics of communities. Mapping the distributions of these species is important to understanding the ecology of coastal environments, managing marine ecosystems (e.g., spatial planning), predicting consequences of climate change and the potential for carbon production. We demonstrate how combining seafloor mapping technologies (LiDAR and multibeam bathymetry) and models of wave energy to map the distribution and relative abundance of seaweed forests of Ecklonia radiata can provide complete coverage over hundreds of square kilometers. Using generalized linear mixed models (GLMMs), we associated observations of E. radiata abundance from video transects with environmental variables. These relationships were then used to predict the distribution of E. radiata across our 756.1km2 study area off the coast of Victoria, Australia. A reserved dataset was used to test the accuracy of these predictions. We found that the abundance distribution of E. radiata is strongly associated with depth, presence of rocky reef, curvature of the reef topography, and wave exposure. In addition, the GLMM methodology allowed us to adequately account for spatial autocorrelation in our sampling methods. The predictive distribution map created from the best GLMM predicted the abundance of E. radiata with an accuracy of 72%. The combination of LiDAR and multibeam bathymetry allowed us to model and predict E. radiata abundance distribution across its entire depth range for this study area. Using methods like those presented in this study, we can map the distribution of macroalgae species, which will give insight into ecological communities, biodiversity distribution, carbon uptake, and potential sequestration.

Predictive mapping of abalone fishing grounds using remotely-sensed LiDAR and commercial catch data

Defining the geographic extent of suitable fishing grounds at a scale relevant to resource exploitation for commercial benthic species can be problematic. Bathymetric light detection and ranging (LiDAR) systems provide an opportunity to enhance ecosystem-based fisheries management strategies for coastally distributed benthic fisheries. In this study we define the spatial extent of suitable fishing grounds for the blacklip abalone (Haliotis rubra) along 200 linear kilometers of coastal waters for the first time, demonstrating the potential for integration of remotely-sensed data with commercial catch information. Variables representing seafloor structure, generated from airborne bathymetric LiDAR were combined with spatially-explicit fishing event data, to characterize the geographic footprint of the western Victorian abalone fishery, in south-east Australia. A MaxEnt modeling approach determined that bathymetry, rugosity and complexity were the three most important predictors in defining suitable fishing grounds (AUC = 0.89). Suitable fishing grounds predicted by the model showed a good relationship with catch statistics within each sub-zone of the fishery, suggesting that model outputs may be a useful surrogate for potential catch.

Spatial abundance patterns and recruitment of a virus-affected commercial mollusc fishery

 Infectious pathogens figure prominently among those factors threatening marine wildlife. Mass mortality events caused by pathogens can fundamentally alter the structure of wild fish stocks and depress recruitment rates and yield. In the most severe instances, this can precipitate stock collapses resulting in dramatic economic losses to once valuable commercial fisheries. An outbreak of a herpes-like virus among commercially fished abalone populations in the south-west fishery of Victoria, Australia, during 2006-2007, has been associated with high mortality rates among all cohorts. Long-term records from fishery-independent surveys of blacklip abalone Haliotis rubra (Leach) enabled abundance from pre- and post-viral periods to be analysed to estimate stock density and biomass. The spatial distribution of abundance in relation to physical habitat variables derived from high-resolution bathymetric LiDAR data was investigated. Significant differences were observed in both measures between pre- and post-viral periods. Although there was some limited evidence of gradual stock improvement in recent years, disease-affected reefs have remained below productivity rates prior to the disease outbreak suggesting a reduction in larval availability or settlement success. This was corroborated by trends in sublegal sized blacklip abalone abundance that has yet to show substantial recovery post-disease. Abundance data were modelled as a function of habitat variables using a generalised additive model (GAM) and indicated that high abundance was associated with complex reef structures of coastal waters (<15 m). This study highlights the importance of long-term surveys to understand abalone recovery following mass mortality and the links between stock abundance and seafloor variability.

Nearshore-offshore-basin species diversity and body size variation patterns in Late Permian (Changhsingian) brachiopods

 B.V. Body size is a fundamental and defining character of an organism, and its variation in space and time is generally considered to be a function of its biology and interactions with its living environment. A great deal of body size related ecological and evolutionary research has been undertaken, mostly in relation to extant animals. Among the many body size-related hypotheses proposed and tested, the size-bathymetry relationship is probably the least studied. In this study, we compiled a global body size dataset of Changhsingian (Late Permian, ca. 254. Ma-252. Ma) brachiopod species from low-latitude areas (30°S-30°N) and analyzed their species diversity and body size distribution patterns in relation to the nearshore-offshore-basin bathymetric gradient. The dataset contained 1768 brachiopod specimens in 435 species referred to 159 genera and 9 orders, from 135 occurrences (localities) of 18 different palaeogeographic regions. Treating the whole of the Changhsingian Stage as a single time slice, we divided the nearshore-offshore-basin bathymetric gradient into three broad depth-related environments: nearshore, offshore and basinal environments, and compared how the species diversity and body size varied along this large-scale bathymetric gradient.Here, we report an array of complex patterns. First, we found a clear overall inverse correlation between species diversity and water depth along the nearshore-offshore-basin gradient, with most species concentrating in the nearshore environment. Second, when the median sizes of all low-latitude brachiopod species from the three environments were compared, we found that there was no significant size difference between the nearshore and offshore environments, suggesting that neither the wave base nor the hydrostatic pressure exerts a critical influence on the body size of brachiopods. On the other hand, the median sizes of brachiopods from the nearshore environment and, to a lesser extent, the offshore environment were found to be significantly larger than that of basinal brachiopods. This trend of significant size reduction in basinal brachiopods mirrors the relative low species diversity in the basinal environment, and neither can be easily explained by the tendency of decreasing food availability towards deeper sea environments. Rather, both trends are consistent with the hypothesis of an expanding Oxygen Minimum Zone (OMZ) in the bathyal (slope to deepsea) environments, where hypoxic to anoxic conditions are considered to have severely restricted the diversification of benthos and favored the relative proliferation of small-sized brachiopods. Finally, a significant difference was also found between eurybathic and stenobathic species in their body size response to the nearshore-offshore-basin gradient, in that eurybathic species (species found in all three environments) did not tend to change their body size significantly according to depth, whereas stenobathic forms (species restricted to a single environment) exhibit a decline in body size towards the basinal environment. This pattern is interpreted to suggest that bathymetrically more tolerant species are less sensitive to depth control with respect to their body size change dynamics, in contrast to stenobathic species which tend to grow larger in shallower water depths.

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.