Deep-sea polymetallic nodule mining is in the exploration phase at present with some groups proposing a move towards extraction within years1. Management of this industry requires evidence of the long-term effects on deep-sea ecosystems2, but the ability of seafloor ecosystems to recover from impacts over decadal scales is poorly understood3. Here we show that, four decades after a test mining experiment that removed nodules, the biological impacts in many groups of organisms are persistent, although populations of several organisms, including sediment macrofauna, mobile deposit feeders and even large-sized sessile fauna, have begun to re-establish despite persistent physical changes at the seafloor. We also reveal that areas affected by plumes from this small-scale test have limited detectable residual sedimentation impacts with some biological assemblages similar in abundance compared to control areas after 44 years. Although some aspects of the modern collector design may cause reduced physical impact compared to this test mining experiment, our results show that mining impacts in the abyssal ocean will be persistent over at least decadal timeframes and communities will remain altered in directly disturbed areas, despite some recolonization. The long-term effects seen in our study provide critical data for effective management of mining activities, if they occur, including minimizing direct impacts and setting aside an effective network of protected areas4,5.

Mineral prospectivity mapping constitutes an efficient tool for delineating areas of highest interest to guide future exploration. Multiple knowledge-driven approaches have been applied for the creation of prospectivity maps for deep-sea ferromanganese (Fe-Mn) crusts over the last decades. The results of a data-driven approach making use of an extensive data collection exercise on occurrences of Fe-Mn crusts in the World Ocean and recent increase in global marine datasets are presented. A Random Forest machine learning algorithm is applied, and results compared with previously established expert-driven maps. Optimal predictive conditions for the algorithm are observed for (i) a forest size superior to a hundred trees, (ii) a training dataset larger than 10%, and (iii) a number of predictors to be used as nodes superior to two. The confusion matrix and out-of-bag errors on the remaining unused data highlight excellent predictive capabilities of the trained model with a prediction accuracy for Fe-Mn crusts of 87.2% and 98.2% for non-crusts locations, with a Kohen’s K index of 0.84, validating its application for prediction at the World scale. The slope of the seafloor, sediment thickness, sediment type, biological productivity, and abyssal mountain constitute the five strongest explanatory variables in predicting the occurrence of Fe-Mn crusts. Most ‘hand-drawn’ knowledge-driven prospective areas are also considered prospective by the random forest algorithm with notable exceptions along the coast of the American continent. However, poor correlation is observed with knowledge-driven GIS-based criterion mapping as the Random Forest considers un-prospective most target areas from the GIS approach. Overall, the Random Forest prediction performs better in predicting a high chance of Fe-Mn crust occurrence in ISA licensed area than the GIS approach, which constitutes an external validation of the predictive quality of the random forest model.

The mesophotic domain is a poorly explored part of the oceans, notably in the Mediterranean Sea. Benthic communities in these depths are not well documented and as such are under higher risk from anthropogenic impacts. Hard substrate habitats in this depth window are not common and are a key ecotope. The Malta Plateau in the central Mediterranean, which is characterized by low sedimentation rates, offers a potentially unexplored domain for these features. Bathymetric and backscatter data offshore of the eastern coast of the island of Malta were used to identify > 1,700 small structures in mesophotic depths. These structures were verified to be biogenic mounds by dives. The mounds extend from several meters to tens of meters in diameter and occur in two main depth windows −40 to 83 meters below present sea level (mbpsl) and 83–120 mbpsl—each formed probably in a different stage during the last glacial cycle. The mounds are composed of interlocking bioconstruction by encrusting organisms and are colonized by sponges and various cold water corals (most of which are protected; e.g., Madrepora oculata). This unique and important habitat is currently under grave threat by human activity, most immediately by trawling and anchoring activity.

Accurate knowledge of the spatial extent of seagrass habitats is essential for monitoring and management purposes given their ecological and economic significance. Extent data are typically presented in binary (presence/absence) or arbitrary, semi-quantitative density bands derived from low-resolution satellite imagery, which cannot resolve fine-scale features and intra-habitat variability. Recent advances in consumer-grade unoccupied aerial vehicles (UAVs) have advanced our ability to survey large areas at higher resolution and at lower cost. This has improved the accessibility of mapping technologies to developing coastal nations, where a large proportion of the world’s seagrass habitats are found. Here, we present the application of UAV-gathered imagery to determine seagrass habitat extent and percent of canopy cover. Four contrasting sites were surveyed in the Turneffe Atoll Marine Reserve, Belize, and seagrass canopy cover was ground truthed from in situ quadrats. Orthomosaic images were created for each site from the UAV-gathered imagery. Three modelling techniques were tested to extrapolate the findings from quadrats to spatial information, producing binary (random forest) and canopy cover (random forest regression and beta regression) habitat maps. The most robust model (random forest regression) had an average absolute error of 6.8–11.9% (SE of 8.2–14), building upon previous attempts at mapping seagrass density from satellite imagery, which achieved errors between 15–20% approximately. The resulting maps exhibited great intra-habitat heterogeneity and different levels of patchiness, which were attributed to site energetics and, possibly, species composition. The extra information in the canopy cover maps provides greater detail and information for key management decisions and provides the basis for future spatial studies and monitoring programmes.

Currently, several satellite-derived bathymetry (SDB) workflows are based on a variety of satellite imagery which are analyzed by empirical or analytical methods. The latest availability of PRISMA hyperspectral data provides a new opportunity for testing their application in shallow water bathymetry mapping. Here we utilize two Level-2 PRISMA scenes from the Caribbean Sea capturing seafloor areas with diverse benthic features and we analyze them using the shallow water analytical models provided by the water-color simulator (WASI) software. The presented study examines the influence of spatial resolution and end-member spectra on the SDB output. Consequently, in one study area we apply inversion using additional reference spectra and in the other study area we exploit the PRISMA panchromatic band for producing a pan-sharpened, hyperspectral cube for bathymetry inversion. The results show good correlation with reference bathymetry data (sonar and admiralty chart) suggesting that PRISMA imagery has a clear potential in optical bathymetry studies. The use of appropriate end-member spectra assists in enhancing the accuracy of SDB, and pan-sharpened PRISMA imagery assists in improving the results when detailed bathymetry is required. PRISMA imagery can be effectively analyzed with open-source software WASI-2D and thus contribute new bathymetric data to regional-scale seafloor mapping projects.

Carbon capture and storage is a key mitigation strategy proposed for keeping the global temperature rise below 1.5 degrees C. Offshore storage can provide up to 13% of the global CO2 reduction required to achieve the Intergovernmental Panel on Climate Change goals. The public must be assured that potential leakages from storage reservoirs can be detected and that therefore the CO2 is safely contained. We conducted a controlled release of 675 kg CO2 within sediments at 120 m water depth, to simulate a leak and test novel detection, quantification and attribution approaches. We show that even at a very low release rate (6 kg day(-1)), CO2 can be detected within sediments and in the water column. Alongside detection we show the fluxes of both dissolved and gaseous CO2 can be quantified. The CO2 source was verified using natural and added tracers. The experiment demonstrates that existing technologies and techniques can detect, attribute and quantify any escape of CO2 from sub-seabed reservoirs as required for public assurance, regulatory oversight and emissions trading schemes.

A huge amount of seabed acoustic reflectivity data has been acquired from the east to the west side of the southern Adriatic Sea (Mediterranean Sea) in the last 18 years by CNR-ISMAR. These data have been used for geological, biological and habitat mapping purposes, but a single and consistent interpretation of them has never been carried out. Here, we aimed at coherently interpreting acoustic data images of the seafloor to produce a benthic habitat map of the southern Adriatic Sea showing the spatial distribution of substrates and biological communities within the basin. The methodology here applied consists of a semi-automated classification of acoustic reflectivity, bathymetry and bathymetric derivatives images through object-based image analysis (OBIA) performed by using the ArcGIS tool RSOBIA (Remote Sensing OBIA). This unsupervised image segmentation was carried out on each cruise dataset separately, then classified and validated through comparison with bottom samples, images, and prior knowledge of the study areas.

Information on the spatial distribution of habitats and vulnerable species is important for conservation planning. In particular, detailed knowledge on connectivity of marine ecosystems in relation to depth and seafloor characteristics is crucial for any proposed conservation and management actions. Yet, the bulk of the seafloor remains under-sampled, unstudied and unmapped, thereby limiting our understanding of connections between shallow and deep-water communities. Recent studies on mesophotic coral ecosystems (MCEs) have highlighted the Western Indian Ocean as a particularly understudied marine region. Here we utilise an autonomous underwater vehicle (AUV) to collect in-situ temperature, oxygen concentration, bathymetry, acoustic backscatter and photographic data on benthic communities from shallow (<30 m) and mesophotic (30-150 m) depths at selected sites in the Greater Pemba Channel, Tanzania. Further, we use generalised additive models (GAMs) to determine useful predictors of substratum (hard and sand) and benthic community type (coral, turf algae, fleshy algae, fish). Our results revealed the presence of a complex seafloor characterised by pockmarks, steep slopes, submarine walls, and large boulders. Photographs confirmed the presence of MCE composed of corals, algae and fishes on the eastern margins of the Pemba Channel. The GAMs on the presence and absence of benthic community explained 35%-91% of the deviance in fish and fleshy algae assemblages, respectively. Key predictors of the distribution of hard substrata and the coral reef communities were depth, showing the upper boundary of MCEs present at 30-40 m, and seafloor slope that showed more occurrences on steep slopes. The upper 100 m of water column had stable temperatures (25-26 degrees C) and oxygen concentrations (220-235 mu mol/l). We noted the presence of submarine walls, steeply inclined bedrock, which appeared to support a highly bio-diverse community that may be worthy of particular conservation measures. Our results also highlight the capability of using marine robotics, particularly autonomous vehicles, to fill the knowledge gap for areas not readily accessible by divers or with surface vessels, and their potential application for the initial survey and subsequent monitoring of Marine Protected Areas.

The ocean floor, its species and habitats are under pressure from various human activities. Marine spatial planning and nature conservation aim to address these threats but require sufficiently detailed and accurate maps of the distribution of seabed substrates and habitats. Benthic habitat mapping has markedly evolved as a discipline over the last decade, but important challenges remain. To test the adequacy of current data products and classification approaches, we carried out a comparative study based on a common dataset of multibeam echosounder bathymetry and backscatter data, supplemented with groundtruth observations. The task was to predict the spatial distribution of five substrate classes (coarse sediments, mixed sediments, mud, sand, and rock) in a highly heterogeneous area of the south-western continental shelf of the United Kingdom. Five different supervised classification methods were employed, and their accuracy estimated with a set of samples that were withheld. We found that all methods achieved overall accuracies of around 50%. Errors of commission and omission were acceptable for rocky substrates, but high for all sediment types. We predominantly attribute the low map accuracy regardless of mapping approach to inadequacies of the selected classification system, which is required to fit gradually changing substrate types into a rigid scheme, low discriminatory power of the available predictors, and high spatial complexity of the site relative to the positioning accuracy of the groundtruth equipment. Some of these issues might be alleviated by creating an ensemble map that aggregates the individual outputs into one map showing the modal substrate class and its associated confidence or by adopting a quantitative approach that models the spatial distribution of sediment fractions. We conclude that further incremental improvements to the collection, processing and analysis of remote sensing and sample data are required to improve map accuracy. To assess the progress in benthic habitat mapping we propose the creation of benchmark datasets.

In this paper we present the seabed maps of the shallow-water areas of Lampedusa and Linosa, belonging to the Pelagie Islands Marine Protected Area. Two surveys were carried out ("Lampedusa 2015" and "Linosa 2016") to collect bathymetric and acoustic backscatter data through the use of a Reson SeaBat 7125 high-resolution multibeam system. Ground-truth data, in the form of grab samples and diver video-observations, were also collected during both surveys. Sediment samples were analyzed for grain size, while video images were analyzed and described revealing the acoustic seabed and other bio-physical characteristics. A map of seabed classification, including sediment types and seagrass distribution, was produced using the tool Remote Sensing Object Based Image Analysis (RSOBIA) by integrating information derived from backscatter data and bathy-morphological features, validated by ground-truth data. This allows to create a first seabed maps (i.e. benthoscape classification), of Lampedusa and Linosa, at scale 1:20 000 and 1:32 000, respectively, that will be checked and implemented through further surveys. The results point out a very rich and largely variable marine ecosystem on the seabed surrounding the two islands, with the occurrence of priority habitats, and will be of support for a more comprehensive maritime spatial planning of the Marine Protected Area.

The flanks of the ultraslow-spreading Mid-Cayman Spreading Center (MCSC) are characterized by domal massifs or oceanic core complexes (OCCs). The most prominent of these, Mount Dent, comprises lower-crustal and upper-mantle lithologies and hosts the Von Damm vent field similar to 12km west of the axial deep. Here, presented autonomous underwater vehicle-derived swath sonar (multibeam) mapping and deep-towed side-scan sonar imagery lead to our interpretation that: (i) slip along the OCC-bounding detachment fault is ceasing, (ii) the termination zone, where detachment fault meets the hanging wall, is disintegrating, (iii) the domed surface of the OCC is cut by steep north-south extensional faulting, and (iv) the breakaway zone is cut by outward facing faults. The Von Damm vent field and dispersed pockmarks on the OCC's south flank further suggest that hydrothermal fluid flow is pervasive within the faulted OCC. On the axial floor of the MCSC, bright acoustic backscatter and multibeam bathymetry reveal: (v) a volcanic detachment hanging wall, (vi) a major fault rifting the southern flank of Mount Dent, and (vii) a young axial volcanic ridge intersecting its northern flank. These observations are described by a conceptual model wherein detachment faulting and OCC exhumation are ceasing during an increase in magmatic intrusion, brittle deformation, and hydrothermal circulation within the OCC. Together, this high-resolution view of the MCSC provides an instructive example of how OCCs, formed within an overall melt-starved ultraslow spreading center, can undergo magmatism, hydrothermal activity, and faulting in much the same way as expected in magmatically more robust slow-spreading centers elsewhere. Plain Language Summary Mid-ocean ridges create the seafloor of the ocean basins, and the processes involved range from the eruption of volcanoes to the generation of submarine hydrothermal vents. A small spreading center in the Caribbean Seathe Mid-Cayman Spreading Centerspreads at such a slow rate that instead of a steady eruption of volcanic material, tectonic faults pull the deeper parts of the Earth to the surface forming Mount Everest-scale mountains. In 2010 and 2013 we collected data from the most prominent massif in the Cayman Trough using cutting-edge mapping tools such as an autonomous underwater vehicle, or AUV. Our data show that the massif is complex internally, breaking up by tectonic forces, intruded by magma, and eroded at the surface. All of these processes are related to the flow of fluids through the crust that sustain the Von Damm vent field, a globally distinct vent field in terms of its biology and chemistry. Mount Dent is likely similar to many other such complexes around the world, including other deep, slow-spreading centers.

Carbonate escarpments are submarine limestone and dolomite cliffs that have been documented in numerous sites around the world. Their geomorphic evolution is poorly understood due to difficulties in assessing escarpment outcrops and the limited resolution achieved by geophysical techniques across their steep topographies. The geomorphic evolution of carbonate escarpments in the Mediterranean Sea has been influenced by the Messinian salinity crisis (MSC). During the MSC (5.97–5.33Ma), the Mediterranean Sea became a saline basin due to a temporary restriction of the Atlantic-Mediterranean seaway, resulting in the deposition of more than one million cubic kilometres of salt. The extent and relative chronology of the evaporative drawdown phases associated to the MSC remain poorly constrained. In this paper we combine geophysical and sedimentological data from the central Mediterranean Sea to reconstruct the geomorphic evolution of the Malta Escarpment and infer the extent and timing of evaporative drawdown in the eastern Mediterranean Sea during the MSC. We propose that, during a MSC base-level fall, fluvial erosion formed a dense network of canyons across the Malta Escarpment whilst coastal erosion developed extensive palaeoshorelines and shore platforms. The drivers of geomorphic evolution of the Malta Escarpment after the MSC include: (i) canyon erosion by submarine gravity flows, with the most recent activity taking place <2600cal. years BP; (ii) deposition by bottom currents across the entire depth range of the Malta Escarpment; (iii) tectonic deformation in the southern Malta Escarpment in association with a wrench zone; (iv) widespread, small-scale sedimentary slope failures preconditioned by oversteepening and loss of support due to canyon erosion, and triggered by earthquakes. We carry out an isostatic restoration of the palaeoshorelines and shore platforms on the northern Malta Escarpment to infer an evaporative drawdown of 1800–2000m in the eastern Mediterranean Sea during the MSC. We interpret the occurrence of pre-evaporite sedimentary lobes in the western Ionian Basin as suggesting that either evaporative drawdown and canyon formation predominantly occurred before salt deposition, or that only the latest salt deposition at the basin margin occurred after the formation of the sedimentary lobes.

27 The flanks of the ultraslow-spreading Mid-Cayman Spreading Center (MCSC) are characterized by 28 domal massifs, or oceanic core complexes (OCCs). The most prominent of these, Mt. Dent, comprises 29 lower-crustal and upper-mantle lithologies and hosts the Von Damm vent field (VDVF) ~12 km west of the 30 axial deep. Here-presented AUV-derived swath sonar (multibeam) mapping and deep-towed side-scan sonar 31 imagery lead to our interpretation that: (i) slip along the OCC-bounding detachment fault is ceasing, (ii) the 32 termination zone, where detachment fault meets the hanging wall, is disintegrating, (iii) the domed surface 33 of the OCC is cut by steep north-south extensional faulting, and (iv) the breakaway zone is cut by outward34 facing faults. The VDVF and dispersed pockmarks on the OCC’s south flank further suggest that 35 hydrothermal fluid flow is pervasive within the faulted OCC. On the axial floor of the MCSC, bright 36 acoustic backscatter and multibeam bathymetry reveal: (v) a volcanic detachment hanging wall, (vi) a major 37 fault rifting the southern flank of Mt. Dent, and (vii) a young axial volcanic ridge intersecting its northern 38 flank. These observations are described by a conceptual model wherein detachment faulting and OCC 39 exhumation are ceasing during an increase in magmatic intrusion, brittle deformation, and hydrothermal 40 circulation within the OCC. Together, this high-resolution view of the MCSC provides an instructive 41 example of how OCCs formed within an overall melt-starved ultraslow spreading center can undergo 42 magmatism, hydrothermal activity, and faulting in much the same way as expected in magmatically more 43 robust slow-spreading centers elsewhere. 44 45

The decarbonisation of our energy supply is reliant on new technologies that are raw material intensive and will require a significant increase in the production of metals to sustain them. Ferromanganese (FeMn) crusts are seafloor precipitates, enriched in metals such as cobalt and tellurium, both of which have a predicted future demand above current production rates. In this study, we investigate the texture and composition of FeMn crusts on Tropic Seamount, a typical Atlantic guyot off the coast of western Africa, as a basis for assessing the future mineral resource potential of Atlantic Seamounts. The majority of the summit is flat and covered by FeMn crusts with average thicknesses of 3–4 cm. The crusts are characterized by two dominant textures consisting of either massive pillared growth or more chaotic, cuspate sections of FeMn oxides, with an increased proportion of detrital and organic material. The Fe, Mn, and Co contents in the FeMn oxide layers are not affected by texture. However, detrital material and bioclasts can form about 50% of cuspate areas, and the dilution effect of this entrained material considerably reduces the Fe, Mn, and Co concentrations if the bulk samples are analyzed. Whilst Tropic Seamount meets many of the prerequisites for a crust mining area, the thickness of the crusts and their average metal composition means extraction is unlikely to be viable in the near future. The ability to exploit more difficult terrains or multiple, closely spaced edifices would make economic feasibility more likely.

The most significant breakthroughs in science are often made as a result of technological developments and innovation. A new capacity to gather more data, measure more precisely or make entirely new observations generally leads to new insights and fundamental understanding. The future of ocean research and exploration therefore lies in robotics: marine robotic systems can be deployed at depths and in environments that are out of direct reach for humans, they can work around the clock, and they can be autonomous, freeing up time and money for other activities. To advance the field of submarine geomorphology, the two types of robots that currently make the biggest difference are Remotely Operated Vehicles (ROVs) and Autonomous Underwater Vehicles (AUVs). Other autonomous or robotic systems are available for marine research (e.g. gliders, autonomous surface vehicles, benthic crawlers etc.), but their application for geomorphological studies is less extensive. This chapter gives an overview of the main characteristics of ROVs and AUVs, their advantages and disadvantages, and their main applications for geomorphological research. In comparison to the other geomorphological methods discussed in this book, however, it has to be made clear that ROVs and AUVs are not so much methods per se, instead they are platforms from which existing and new approaches can be applied.

The Chain Fracture Zone is a 300‐km‐long transform fault that offsets the Mid‐Atlantic Ridge. We analyzed new multibeam bathymetry, backscatter, gravity, and magnetic data with 100% multibeam bathymetric data over the active transform valley and adjacent spreading segments as part of the Passive Imaging of the Lithosphere Asthenosphere Boundary (PI‐LAB) Experiment. Analyses of these data sets allow us to determine the history and mode of crustal formation and the tectonic evolution of the transform system and adjacent ridges over the past 20 Myr. We model the total field magnetic anomaly to determine the age of the crust along the northern ridge segment to better establish the timing of the variations in the seafloor fabric and the tectonic‐magmatic history of the region. Within the active transform fault zone, we observe four distinct positive flower structures with several en échelon fault scarps visible in the backscatter data. We find up to −10 mGal residual Mantle Bouguer Anomaly in the region of the largest positive flower structure within the transform zone suggesting crustal thickening relative to the crustal thinning typically observed in fracture zones in the Atlantic. The extensional/compressional features observed in the Chain Transform are less pronounced than those observed further north in the Vema, St. Paul, and Romanche and may be due to local ridge segment adjustments.

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

Submerged flanks of volcanic islands are prone to hazards including submarine landslides that may trigger damaging tsunamis and sediment-laden seafloor flows (called ‘turbidity currents’). These hazards can break seafloor infrastructure which is critical for global communications and energy transmission. Small Island Developing States are particularly vulnerable to these hazards due to their remote and isolated nature, small size, high population densities and weak economies. Despite their vulnerability, few detailed offshore surveys exist for such islands, resulting in a geohazard ‘blindspot’, particularly in the South Pacific. Understanding how these hazards are triggered is important; however, pin-pointing specific triggers is challenging as most studies have been unable to link continuously between onshore and offshore environments, and focus primarily on large-scale eruptions with sudden production of massive volumes of sediment. Here we focus on an area with a high volcanic sediment supply, which is also similar to locations (volcanic or not) where human-induced vegetation change over-supplies sediments to coastal margins. We address these issues by integrating the first detailed (2 m x 2 m) bathymetry data acquired from Tanna Island, Vanuatu with a combination of terrestrial remote sensing data, onshore and offshore sediment sampling, and documented historical events. Mount Yasur on Tanna has experienced low-magnitude Strombolian activity for at least the last 600 years. We find clear evidence for submarine landslides and turbidity currents, yet none of the identified triggers are related to major volcanic eruptions, in contrast to conclusions from several previous studies. Instead we find that cascades of non-volcanic events (including outburst floods with discharges of >1000 m3/s, and tropical cyclones), that may be separated by decades, are more important for preconditioning and triggering of landslides and turbidity currents in oversupplied sedimentary regimes such as at Tanna. We conclude with a general model for how submarine landslides and turbidity currents are triggered at volcanic and other heavily eroding mountainous islands. Our model highlights the often-ignored importance of outburst floods, non-linear responses to land-use and climatic changes, and the complex interactions between a range of coastal and tectonic processes that may overshadow volcanic regimes.

Multibeam sonar systems now routinely record seafloor backscatter data, which are processed into backscatter mosaics and angular responses, both of which can assist in identifying seafloor types and morphology. Those data products are obtained from the multibeam sonar raw data files through a sequence of data processing stages that follows a basic plan, but the implementation of which varies greatly between sonar systems and software. In this article, we provide a comprehensive review of this backscatter data processing chain, with a focus on the variability in the possible implementation of each processing stage. Our objective for undertaking this task is twofold: (1) to provide an overview of backscatter data processing for the consideration of the general user and (2) to provide suggestions to multibeam sonar manufacturers, software providers and the operators of these systems and software for eventually reducing the lack of control, uncertainty and variability associated with current data processing implementations and the resulting backscatter data products. One such suggestion is the adoption of a nomenclature for increasingly refined levels of processing, akin to the nomenclature adopted for satellite remote-sensing data deliverables.