Understanding Mn-nodule distribution and evaluation of related deep-sea mining impacts using AUV-based hydroacoustic and optical data

In this study, shipand autonomous underwater vehicle (AUV)-based multibeam data from the German ferromanganese-nodule (Mn-nodule) license area in the Clarion–Clipperton Zone (CCZ; eastern Pacific) are linked to ground-truth data from optical imaging. Photographs obtained by an AUV enable semi-quantitative assessments of nodule coverage at a spatial resolution in the range of meters. Together with high-resolution AUV bathymetry, this revealed a correlation of small-scale terrain variations (< 5 m horizontally, < 1 m vertically) with nodule coverage. In the presented data set, increased nodule coverage could be correlated with slopes> 1.8 and concave terrain. On a more regional scale, factors such as the geological setting (existence of horst and graben structures, sediment thickness, outcropping basement) and influence of bottom currents seem to play an essential role for the spatial variation of nodule coverage and the related hard substrate habitat. AUV imagery was also successfully employed to map the distribution of resettled sediment following a disturbance and sediment cloud generation during a sampling deployment of an epibenthic sledge. Data from before and after the “disturbance” allow a direct assessment of the impact. Automated image processing analyzed the nodule coverage at the seafloor, revealing nodule blanketing by resettling of suspended sediment within 16 h after the disturbance. The visually detectable impact was spatially limited to a maximum of 100 m distance from the disturbance track, downstream of the bottom water current. A correlation with high-resolution AUV bathymetry reveals that the blanketing pattern varies in extent by tens of meters, strictly following the bathymetry, even in areas of only slightly undulating seafloor (<1 m vertical change). These results highlight the importance of detailed terrain knowledge when engaging in resource assessment studies for nodule abundance estimates and defining mineable areas. At the same time, it shows the importance of high-resolution mapping for detailed benthic habitat studies that show a heterogeneity at scales of 10 to 100 m. Terrain knowledge is also needed to determine the scale of the impact by seafloor sediment blanketing during mining operations.


Nodule abundance estimation in relation to benthic structures
The deep ocean is an area of economic interest due to its potential reserve of metal resources. Before deep-sea mining can be conducted, a better understanding is required of the ecological role of the deep sea as the largest habitat on Earth. One focus lies on impacts of ferromanganese-nodule (Mn-nodule) mining which recently has been studied in international projects like MIDAS (FP7 project 603 418) and Mining Impact (JPI Oceans project). Mn nodules form a hard substrate for sessile fauna (Purser et al., 2016;Vanreusel et al., 2016) and their removal is expected to impact respective fauna, but mobile fauna is impacted as well (Bluhm et al., 1995). Quantifying Mn-nodule occurrence and understanding distribution patterns on the seafloor are thus required for Published by Copernicus Publications on behalf of the European Geosciences Union.
environmental baseline studies as well as ecological and impact assessments. Several studies correlate bathymetry and nodule occurrence, revealing a complex/non-coherent interrelation which mainly depends on the considered spatial scale. Most studies have focused on nodule occurrence variability between very different terrain settings such as seamounts, valleys, plains and undulating terrain (Halbach, 1988;Pattan and Kodagali, 1988;Skornyakova and Murdmaa, 1992;Sharma and Kodagali, 1993;Park et al., 1997;Jung et al., 2001;Kim et al., 2012). The detected Mn-nodule variability has been associated with sediment deposition properties, e.g., assuming increased accumulation of sediment in flat or depression areas compared to sloping seafloor (Frazer and Fisk, 1981;Widmann et al., 2014). Local sediment accumulation influences the dominant formation type of the Mn nodules (diagenetic vs. hydrogenetic), their size and metal concentration (Jung et al., 2001;Kim et al., 2012;Mewes et al., 2014;Widmann et al., 2014), but detailed small-scale investigations (1 to 100 m scale) are not commonly done. Okazaki and Tsune (2013) propose to perform such investigations to improve our knowledge about Mn-nodule formation processes and the affecting parameters. Moreover, the substrate changes considered in this study provide relevant information for estimating size and heterogeneity of local-scale habitats.
The study analyzes ship-based bathymetric data for largescale background information together with autonomous underwater vehicle (AUV)-obtained high-resolution multibeam (MB) and optical data to reveal detailed nodule coverage patterns within a 12 km 2 area. The resolution of the acoustical AUV data enables the identification of vertical morphological undulations of less than 1 m with a lateral resolution of 3-5 m. Results from optical and hydroacoustic data were analyzed to assess correlations between Mn-nodule coverage and small-scale morphology, and to extrapolate Mn-nodule occurrences to a wider area.
An equivalent approach was applied for an environmental impact study on sediment blanketing during a simulated "mining operation". Mn-nodule mining will affect the seafloor and benthic fauna in several ways. A removal of the uppermost sediment layer (5-20 cm) will cause habitat loss for sessile fauna which depends on nodules as hard substrate ( Fig. 1; Vanreusel et al., 2016) and for organisms living in the uppermost "fluffy" sediment layer. The suspended sediment plume can clog filter organs of suspension feeders and the redepositioning of suspended sediment will bury sessile organisms. These physical impacts could further be accompanied by biogeochemical disequilibria in the water column and the sediment surface that will impact the local environment on short and long timescales (Shirayama and Fukushima, 1997;Kotlinski and Stoyanova, 1998;Sharma et al., 2001;Thiel and Forschungsverbund Tiefsee-Umweltschutz, 2001). The release of toxic substances or metals might also add to the impact. Together, these effects can increase mortality with unknown short-term and cumulative effects (Markussen,Figure 1. Sessile benthic organisms depending on manganese nodules as a hard substrate habitat. Images are from the German claim area in the Clarion-Clipperton Zone (photos: remotely operated vehicle (ROV) Kiel 6000, GEOMAR Helmholtz Centre for Ocean Research Kiel, Germany). cal observations using video and photo material. These data finally contributed to the development of sediment plume distribution models (Lavelle et al., 1981;Jankowski et al., 1996) and led to interpolated blanketing maps Yamazaki et al., 1997). Observations from these largerscale (as compared to this study) BIEs -OMI, OMA, OMCO (1978), DISCOL (1989), BIE-II (1993), JET (1994), IOM-BIE (1995), INDEX (1997) and MMAJ (1997) -indicated different distribution distances of the created bottom plume ranging from several tens or hundreds of meters Trueblood and Ozturgut, 1997;Sharma, 2001) up to several kilometers away from the disturbances (Burns, 1980;Lavelle et al., 1981;Jankowski et al., 1996;Yamazaki and Kajitani, 1999;Yamada and Yamazaki, 1998).
The use of different disturbance gear for different duration per BIE leads to inconsistent interpretations (Jones, 2000). Unfortunately, definitions for "a sediment plume" differ with regards to minimum particle size and amount of particles, and thus parameters are applied differently in model approaches. While some models calculate the distribution of the resuspended material until all sediment particles have settled from the water column (Andrew Dale, SAMS, personal communication, 2016), others define certain particle concentration thresholds (Burns, 1980). A concise plume definition could be based on thresholds (e.g., with regard to amount of particles, shape and size distribution) that benthic organisms could tolerate on short but also longer cumulative timescales. Defining such thresholds requires in situ experimental data of the reaction of benthic fauna to sediment plumes in the area of potential mining. Detailed studies do not exist but are essential for estimating the ecological consequences of deepsea mining.
Plume model results are based on several assumptions to include parameters describing the environment. Particle sizes and settling velocities are key factors in modeling plume distribution distances (Jankowski et al., 1996) and uncertainties can lead to misinterpretations. Our current understanding of the behavior of resuspended particles in the deep sea is based on laboratory experiments. Often, these experiments struggle to correctly determine settling velocities of flocculating particles and/or they rely on specific deep-sea sediments and might not account for the correct environmental parameters. Nevertheless, modeling the distribution of a plume induced by Mn-nodule mining is the only way to predict the possibly impacted area. Meaningful models need to incorporate all environmental aspects and need to operate at the highest possible resolution.
The study presented here focuses on an area within which the resettled sediment was visually observable in deep-sea photographs of the seafloor. Two AUV photo surveys over the same area were conducted before and after the deployment of an epibenthic sledge (EBS), that created a sediment plume. The two data sets are directly compared to determine the scale of the visible disturbance. The correlation of the photo data with AUV-obtained bathymetry data reveals the influence of the local terrain variability on the sediment blanketing pattern and thus the sediment plume spreading.

Study area description
The study area is part of the Prospective Area 1 (PA1) within the German license area in the eastern Clarion-Clipperton Zone (CCZ). The wider area (Fig. 2) has a mean water depth of 4240 m with abundant isolated or chains of seamounts of various heights. Mostly N-S-trending parallel graben and horst structures originating from the East Pacific Rise (Rühlemann et al., 2011, Fig. 2) can be seen as well. The PA1 itself is a plateau-like area, elevated approximately 150 m above the surrounding terrain and classified as a "mineable plateau" in Fig. 2d. Small-scale undulations (< 100 m) on the plateau are smooth with slopes of < 10 • in wide parts of the area. A stronger relief is associated with seamounts or ridge structures (Fig. 2). To the west, the plateau is bounded by a deep graben structure, whereas towards the east and the south the terrain slopes down to depressions with irregularly scattered seamounts of different sizes (red areas in Fig. 2d).
High-resolution studies using data from several AUV deployments were carried out within the mineable plateau (black square in Fig. 2). This plateau is characterized by slopes of less than 3 • and the area is considered suitable for Mn-nodule collector systems . The AUV study area is located west of a ridge that follows the characteristic N-S lineation and is only slightly structured, showing a smoothly undulating terrain (Fig. 3). Towards the north, the AUV study area is bounded by a slightly elevated area, which dips towards the south, leading into a wider depression. Smaller basins/depressions are found throughout the entire area. Slopes are generally < 10 • and those exceeding this value are found only in association with a group of seamounts in the southwest (Fig. 3). The zoomed-in view in Fig. 3 is based on newly acquired EM122 data (Greinert, 2016) showing a terrain more variable than initially assumed from the data set shown in Fig. 2 (data from an EM120, 2 • by 2 • beam angle system). The more variable data also point towards a further terrain differentiation and potentially a less homogenous Mn-nodule coverage within the area. The acquired AUV multibeam data that were processed to a horizontal resolution of 5 m provide the small-scale morphology. The resulting bathymetric maps were used for correlating nodule coverage and sediment blanketing patterns of the disturbance experiment.  Table A2 (sl indicates slope). The terms "mineable" and "unmineable" are defined by the slope threshold (here, mineable indicates where the slope ≤ 3 • , and unmineable indicates where the slope > 3 • ; this is a rather conservative threshold, as current discussions mention 7 • as more realistic). Black squares mark the study area (center at 117 • 1 W 11 • 51 N) shown in Fig. 3. tem (12 kHz, 0.5 • along-and 1 • across-track beam angle; 55 m cell size). The swath angle was set to 120 • and the survey speed was about 8 kn. The REMUS 6000 type AUV Abyss was deployed for the high-resolution mapping and photo surveys (Linke and Lackschewitz, 2016). A RE-SON Seabat 7125 MBES system was used (200 kHz, 2 • along-track and 1 • across-track beam angle). The AUV-based multibeam surveys were conducted at an altitude of 80 m above the seafloor. Bathymetric maps were produced with the software packages GMT 5.2 (Wessel et al., 2013) and Ar-cGIS 10.2. The MBES data were analyzed in ArcGIS 10.2 as a floating point raster (see the Appendix for cell size of different regions; Table A1). Data were projected as Universal Transverse Mercator (UTM) coordinate system (zone 11N) to enable spatial analysis. Bathymetric first-order derivatives (slope, aspect) and second-order derivatives (bathymetric position index -BPI, vector ruggedness measure -VRM, total curvature, plan curvature, profile curvature) were calculated for each region and subregion using tools in the System for Automated Geoscientific Analyses Geographic Information System (SAGA GIS) and ArcGIS (spatial analyst toolbox and the "benthic terrain modeler" add-on toolbox; Wright et al., 2012, see the Appendix and Table A4 for further details). The AUV camera system DeepSurveyCam (Kwasnitschka et al., 2016) was used for visual seafloor inspection by two photo surveys. Photos were taken from 7 to 9 m altitude at a mean speed of 3 kn, gaining more than 50 000 usable photos. Two sub-areas (A1 and A2), located approximately 5 km apart, were photographed extensively. One part of both AUV surveys in sub-area A2 followed exactly the same track before and approximately 16 h after the deployment of an epibenthic sledge Schoening, 2017a, b). Photos were automatically analyzed for the Mn-nodule coverage on the seafloor (percent coverage) and nodule size distribution (different quantiles in cm 2 ; see the Appendix for details). The Compact Morphology Nodule Delineation (CoMoNoD) algorithm was used for this task . This nodule delineation method is based on color differences between the nodules and the sediment background. Derived quantitative values are georeferenced and can thus be jointly analyzed with the AUVobtained bathymetry. This allows understanding correlations between nodule abundance and the decameter-scale morphological changes in the AUV studied area.
As a side product of benthic sampling using a B-EBS type sledge (Brenke, 2005), a sediment plume was created. The sledge itself has a length of 360 cm and a width of 120 cm with a weight of approximately 420 kg in water. It creates a pressure of about 13 g cm −2 onto the sediment surface (see Fig. A3 in the Appendix). The sledge was towed during station SO239_024-EBS across sub-area A2 (Fig. 4a) from west to east at approximately 0.5 m s −1 leaving a track of approxi-mately 20 cm in depth (Martínez Arbizu and Haeckel, 2015). The resettlement of the sediment plume was visually mapped using the camera system and automated image analysis mentioned above.

AUV-based bathymetry and overview of Mn-nodule coverage
The water depth within the AUV-mapped area ranges from 4110 to 4143 m ( Fig. 4a) with 93 % of the area showing slopes < 3 • (Fig. 4b). Steeper slopes between 7 and 10 • occur locally in the east and towards the north as well as in randomly distributed pit structures, which occur throughout the area (Fig. 4). These pit structures take up approximately 10 % of the area and attain sizes from several tens of meters to 150 m in diameter with a maximum depth of 4 m. They occur exclusively within larger depressions (approximately 50 % of the area) as visualized by the BPI map (Fig. 4c). The two AUV photo surveys provide visual data from within the high-resolution MBES map (black lines, Fig. 4a) covering a depth range from 4134 to 4114 m. The first AUV camera survey (SO239_019_Abyss168) provided data over two extensive sub-areas in the west (A1) and east (A2), as well as one survey line connecting both areas (Figs. 4a and 5). The second survey remapped parts of sub-area A2 (SO239_028_Abyss169) and additionally mapped a similarly sized area further south (Martínez Arbizu and Haeckel, 2015). Based on the automated image analyses, the majority of the seafloor shows nodule coverage values between 8 and 17 % (Fig. 5a). Values below and above this range (< 1 %) have been neglected, since they are caused by unusual objects (like EBS tracks or organisms) in the images.
In the following examinations, the threshold between "low" and "high" Mn-nodule coverage is set at 12.5 %, which is the analyzed mean coverage value of the considered range. In the eastern A2 sub-area, a greater proportion of higher coverage values (13-16 %) can be observed. A positive correlation was found between Mn-nodule coverage and median size of the nodules (Fig. 6a). The correlation decreases with increasing nodule size, indicating a compensation of the size by a decreasing number of occurrences. This is consistent with findings of former studies, e.g., by Okazaki and Tsune (2013).

Broad-scale variability (less detailed, correlation
with ship-based bathymetric data, resolution 100-1000 m) Parts of dive SO239_019_AUV2 run across the entire working area providing data from different terrains that can be linked to the ship-based bathymetric information. The correlation between photo analysis and this less resolving bathymetry indicates a trend of decreasing nodule coverage at elevations and steeper sloping areas (Fig. 5). Video data acquired during previous cruises provide similar observations (Kuhn, 2015). The distribution pattern seen in the imagery also points towards small-scale Mn-nodule coverage variability which is possibly related to minor topographic changes in meter to sub-meter scale. As only one track covers the central region of the working area, clear correlations between Mn-nodule occurrence and large-scale ship-based bathymetry are difficult to assess. Finding clear correlations is further complicated by the uncertainty of the AUV navigation (up to 30 m), which prevents a precise georeferencing of the photos between sub-areas A1 and A2. More robust visual reference data could be provided by conducting a sparse mesh survey across the entire area or by a contiguous photo mosaic across different terrains.

Local-scale variability (more detailed, correlation with AUV-based bathymetric data, resolution 1-100 m)
The assessment of small-scale Mn-nodule coverage heterogeneity was based on the western A1 and eastern A2 sub-areas; here, overlapping photo mosaics and AUV-based bathymetric data in meter resolution exist (Fig. 4a). Sub-area A2 (700 m×500 m, 0.35 km 2 ) is bound to the east by a 5-7 m high "ridge" with a relatively steep slope (slope A2E, 3-7 • ). The western part of this area (A2W) shows only minor morphological variation and a total relief of approximately 2 m (Fig. 4). Despite the rather small relief changes, variations in Mn-nodule coverage can be observed (Figs. 5 and 7). Figure 7a illustrates the detailed bathymetry of the studied area with red dots indicating lower Mn-nodule coverage (≤ 12.5 %) as indicted by image inspection. Those areas with a BPI50 > 0, slopes ≤ 1.8 • and plan curvature values > −0.02 radians m −1 were found to show the best correlation with the lower Mn-nodule coverage in sub-area A2 ( Fig. 6d, f, h). A NW-SE-oriented, elongated patch in the central part of A2 that corresponds to a flat-topped (slope ≤ 1.5 • ), slightly convex-shaped elevated structure (< 1 m above the surrounding terrain) shows a low Mn-nodule coverage. A higher Mn-nodule coverage instead occurs at steeper slopes (> 1.5 • ) and in morphological depressions (negative BPI values, negative plan and total curvature), indicating a sediment depositional environment. Two distinct depression structures (Pit 1 and Pit 2 structures in Fig. 7), both approximately 60-80 m in diameter and 1-2 m deep, show a different pattern; here, the visible Mn-nodule coverage is significantly lower (0-8 %). The almost-spherical pit structures are bound by slopes of > 2 • and thus produce slightly increased vector ruggedness measurement values > 1.0 × 10 −4 (VRMs; Fig. A4a) and the lowest observed BPI-values (Figs. 4c and A4b). Similar structures are observed throughout the entire study area (blue shaded areas in Fig. 7 and circular features seen in the slope map of Fig. 4b).
Based on additional ROV and benthic camera surveys, it is assumed that these pit structures exhibit very few to no Mn nodules (Peukert, 2016).
No further correlation between Mn-nodule coverage and bathymetric derivatives was found, and no relation to absolute water depth could be observed (Fig. 6b). However, Fig. 6b shows a significantly lower variability of Mn-nodule coverage for water depths shallower than 4019 m (only approximately 4 % variability, compared to 6-7 % variability in deeper areas); these areas correspond to steeper slopes (approximately 2 to 3.5 • ) associated with the eastern bounding elevation (slope A2E, Fig. 4a) of sub-area A2. Along this west-facing slope the Mn-nodule coverage clearly increases with increasing depth. Areas featuring low slope values show higher variability in Mn-nodule coverage (Fig. 6b).
A lower Mn-nodule coverage (< 12.5 %) is predicted for the green areas marked in Fig. 7a when using the BPI50, slope and plan curvature classification of the A2 sub-area (Fig. 7b). Although the resulting area does not match completely with the areas of low coverage derived from the photo analyses (red dots/shades in Fig. 7), it represents the best correlation that could be achieved. Based on this result, a Mnnodule coverage of < 12.5 % can be expected in 39 % of the study area (green shades in Fig. 7) and is likely to be very low or zero in at least 1 % of the area (blue shaded parts).
In sub-area A1 (230×600 m, 0.138 km 2 ), no correlation is observed between the photo-analyzed Mn-nodule coverage (red) and the seafloor classification of A2 (green; Fig. 7c). In addition, scatterplots ( Fig. 6c-h) show different dependen- cies between Mn-nodule coverage to BPI, slope and plan curvature between A1 and A2. In both areas, though, the coverage attains more uniform values towards steeper slopes (Fig. 7e, f). A stronger correlation is shown in A2 and an inverse correlation in A1. In both areas, the highest variability but also the lowest values of Mn-nodule coverage occur in generally flat areas (curvature values around 0, low slope values; Fig. 7e, g).
Comparing the terrain statistics of areas A1 and A2 (Fig. 8) reveals differences in their bathymetric settings, which might cause this discrepancy. Sub-area A2 mainly slopes towards west to southwest, as indicated by the aspect distribution and to a lesser degree in "opposite" northeast to east directions caused by the general pattern of N-Sstriking graben and horst structures. In comparison, the main slope direction in A1 is towards a southerly direction. The slope distribution in A1 indicates a dominance of slopes up to 1 • . This is reflected by the large area of flat seafloor as determined by the AUV-BPI440 value distribution (Fig. 8). In A2, slopes are steeper, the terrain is more variable and a larger number of depressions are observed compared to A1. The VRM shows similar values for both areas (Table A6).
These differences in bathymetric derivative values point at a lower terrain variability in A1, confirmed by the more consistent depth values in A1 relative to A2. Considering the generally lower Mn-nodule coverage within A1 (Fig. 5a) it is concluded that lower Mn-nodule coverage correlates with lower terrain variability and lower slope values. This generalized observation is consistent with findings for A2. Although a direct one-to-one relationship valid in both sub-areas could not be derived, the general trend indicates higher Mn-nodule coverage with more variable terrain, along smooth slopes and in concave-shaped terrain (depressions).

Sediment plume resettling
To evaluate sediment plume resettling, results of the automated image-based Mn-nodule detection before the EBS disturbance (SO239_019_Abyss168 with 6061 usable photos) and after the EBS disturbance (SO239_028_Abyss169 with 10 783 usable photos) were compared (Fig. 9a). Areas with the lowest analyzed coverage were associated with sediment blanketing that covers the Mn nodules completely (here de- Figure 6. Scatterplots indicating relationships between Mn-nodule percent coverage (%) and eight other nodule and terrain values: median nodule size in cm 2 (a), depth within the A2 sub-area (b, color-coded by slope), BPI50 (c, d), slope (e, f) and plan curvature (g, h). Charts (c, e, g) relate to sub-area A1; charts (b, d, f, h) to sub-area A2. The sub-areas show different correlations. fined as < 8 % Mn-nodule coverage; green shaded areas in Fig. 9).
The AUV tracks of the photo surveys run perpendicular to the EBS track. A strong sediment blanketing can be observed close to the disturbance track ( Figs. 9 and 10). The photo mosaic shows a sharp transition between low (no) and higher Mn-nodule coverage north of the EBS track over a distance of 5 to 20 m (Fig. 10). South of the track the transition from complete sediment blanketing to areas without visible sediment cover is gradually fading out with increasing distance to the EBS track; a slight sediment blanketing of Mn nodules can be observed up to 70 m away from the EBS track (Figs. 9 and 10). This pattern indicates a southward-directed bottom current, which is confirmed by acoustic Doppler current profiler (ADCP)-based current measurements (station SO239_005; see Fig. A5 in the Appendix). An upwardlooking 300 kHz ADCP (15 min ensembles, 2 m bin sizes) was positioned 500 m SE of the EBS track at the time of the EBS deployment (Martínez Arbizu and Haeckel, 2015).
The combination of AUV-obtained bathymetry and imagery reveals a distinct blanketing pattern depending on the small-scale morphology (Fig. 9). In section A2W, where the total relief is only 1-2 m, the resuspended material was distributed 20 to 30 m towards the north of the EBS track and 40 to 50 m towards the south of it. West of A2W, within Pit  Table A5. the seafloor slope increases to maximum 6 • and the water depth decreases by several meters (Fig. 9). This morphological change causes the sediment plume to cover the seafloor only up to 6 m towards the north of the EBS track. The greatest distance at which sediment has been deposited away from the EBS track occurs within Pit 1 (Fig. 9). Here, the least visible Mn-nodule coverage extends up to 70 m south of the EBS track, significantly further than in the flat A2W section and slope area A2E. ters) over an extensive area of 0.49 km 2 that is completely photomapped. This highly detailed insight is of importance for a spatially detailed evaluation of the small-scale habitat distribution and potentially allows a better resource assessment. However, potential uncertainties for the absolute numbers of Mn-nodule size and coverage are explained here. Photographs only provide information of the sediment surface and thus will not be able to detect buried/sedimentcovered Mn nodules (Sharma and Kodagali, 1993;Sharma et al., 2010Sharma et al., , 2013, resulting in a potential underestimation of the absolute Mn-nodule abundance (Kuhn and Rathke, 2017). For absolutely accurate resource assessments and verification of the presented results, detailed sampling based on this study would need to follow. With respect to mapping Mn nodules as hard grounds for sessile fauna, photographs give a realistic quantitative representation of size/coverage and spatial nodule abundance changes. Nevertheless, the automated image analysis CoMoNoD has some uncertainties . In general, the nodule identification employed here is based on contrast differences between the nodules (dark) and the sediment (bright). Mn nodules on the seafloor could be located too close to each other to be correctly separated by the applied algorithm and, depending on the image quality and the contrast thresholds, quantitative coverage and size distribution values can be inaccurate (Sharma et al., 2010;Schoening et al., 2012Schoening et al., , 2016Schoening et al., , 2017Tsune and Okazaki, 2014) (Fig. A6). Turbidity in the water, backscatter from particles, the water properties and the altitude of the camera (AUV) impact image quality (Edwards et al., 2003;Kwasnitschka et al., 2016). Nevertheless, optical imaging provides information on relative changes in seafloor Mn-nodule coverage and nodule sizes. This makes the presented technique an effective monitoring tool for habitat and environmental impact assessments that investigate the resettling of suspended sediment. We detected a distinct trend of higher coverage correlating with larger nodules (Fig. 7a) that could be a result of imperfect segmentation if the nodule density is too high.

Correlation between bathymetry and nodule occurrence
4.2.1 Regional-scale correlation between ship-based bathymetry and Mn-nodule coverage/size In general, properties such as sedimentation rate (Frazer and Fisk, 1981;Mewes et al., 2014), type and thickness of the sediment (Frazer and Fisk, 1981;Jeong et al., 1994) are believed to determine Mn-nodule growth (von Stackelberg and Beiersdorf, 1991); for sediment deposition environments the interplay between bottom currents and bathymetry plays an important role (Halbach, 1988). The depositional properties vary on a regional scale, considering large geomorphological terrain types, but are also impacted on a local scale of only a few kilometers and even less (Craig, 1979;Frazer and Fisk, 1981;Sharma and Kodagali, 1993;Mewes et al., 2014). Varying considerations of scale and regional differences in nodule exposure between different oceans across different studies have thus led to partly contradicting statements of the relationship between the Mn-nodule coverage/size and bathymetric settings.
Several investigations report small Mn nodules and low coverage in depressions and plains which are considered as sediment accumulation sites, in contrast to seamounts, slopes and crests (Pattan and Kodagali, 1988;Sharma and Kodagali, 1993). Other studies discussed comparatively larger diagenetic Mn nodules in plains which are also considered as sediment accumulation areas. More abundant but smaller hydrogenetic Mn nodules have been observed in more rugged terrain (Skornyakova and Murdmaa, 1992;Kim et al., 2012;Widmann et al., 2014). Such terrains are interpreted to increase current velocities and turbulence caused by channel effects reducing sediment accumulation. Mewes et al. (2014) present a correlation between Mnnodule size and sedimentation rate, where large nodules correlate with a smaller amount of clay fraction in the sediments that they interpreted to be caused by stronger bottom currents/lower sedimentation rate. A similar observation is presented by Skornyakova and Murdmaa (1992), who state that diagenetic/large Mn-nodule formation is linked to a periodical redistribution of the surface sediment layer.
With respect to the large scale of the ship-based bathymetry in Fig. 2, the working area of this study is located in a sediment-accumulating flat terrain with smooth bathymetry, characterized by the occurrence of medium to large (> 4 cm) Mn nodules . However, a more detailed view allows the identification of terrain variability on a scale of several tens to hundreds of meters that enables a more detailed assessment of the associated Mn-nodule coverage variability (Fig. 5). Larger nodules/higher coverage values occur in depressions and at sloping seafloor when compared to broad-scale bathymetry. Larger nodule sizes could be the result of stronger bottom currents preventing/reducing the deposition of sediment on nodules and/or favoring nodule growth. For another area in the German claim, box core (BC) samples taken by the Federal Institute of Geosciences and Natural Resources (BGR Hannover, Germany), revealed larger (diagenetic) nodules in a very broad-scale flat terrain. This area has been classified by (Widmann et al., 2014) as an area of sediment accumulation and is compared to a rougher, supposedly sediment "winnowing" area, with many smaller nodules formed hydrogenetically. The interpretation of sediment accumulating and winnowing areas is based on broad-scale shipbased bathymetry of much coarser resolution compared to this study.
According to the study by Skornyakova and Murdmaa (1992), the data presented here indicate lower sedimentation rates associated with stronger bottom currents in the depressions supporting the growth of larger Mn nodules. Increased bottom currents within the depressions could possibly be induced by convergent channeling or turbulence of bottom currents, which contradicts the assumption of lower current strength and therefore higher sedimentation within depressions.

Local-scale correlation between AUV-based bathymetry and Mn-nodule coverage/size
Variability in Mn-nodule coverage within several tens of meters or less can be correlated with AUV-based bathymetry.
In sub-area A2, patches of low Mn-nodule coverage correlate with low bathymetric elevations even when the relief differs by less than 1 m. The strongest correlation between low Mn-nodule coverage was determined with slightly convex-shaped elevated structures (surfaces < 1 • slope, positive plan curvature and positive BPI values). These parameters most likely define local-scale sedimentation environment affecting the local balance between sediment accumulation and erosion. The presented data show that favorable nodule growth/occurrence conditions coincide with gentle sloping sites and low relief depressions, where sediment is assumed to accumulate slowly. Within sub-area A2, a smaller variability of Mn-nodule coverage can be observed in correlation with slope A2E towards the east. This is in agreement with observations by Sharma and Kodagali (1993), who also observed more uniform nodule coverage in sloping areas. The authors point out that this could be a result of a larger exposure of the Mn nodules rather than absolute difference, since they discovered discrepancies between direct sampling and results of photo analyses.
Rather special for the presented data set are the pronounced pit structures, observed throughout the AUVmapped area with very little to no Mn nodules observed at the sediment surface. This is in contradiction to the wider depressions, where a higher Mn-nodule coverage was observed. The existence of such pronounced depressions most likely leads to a reduction of bottom current velocities resulting in a higher sediment deposition of suspended sediment and potentially even sediment slumping from the sides. This could result in sedimentation rates too high for Mn-nodule formation (Halbach, 1988;Mewes et al., 2014) or a simple cover/burial of previously formed and still-existing Mn nodules below the sediment surface. The formation process of the pits is unclear, but could be karst structures , which are younger than the Mn-nodule forma-tion which would point towards Mn-nodule burial within the pits. At the same time, a higher sedimentation rate in a low current regime would also mean a higher accumulation of clay size particles, which are proposed to not be favorable for nodule growth (Mewes et al., 2014). Another possibility could be that these pit structures are pockmarks, formed by pore water release (Harrington, 1985;Hovland and Judd, 1988) with a significant change in local pore water geochemical properties and eventually warmer temperature that prevented Mn-nodule formation in the past. Unfortunately, the pit structures could not be sampled in more detail and it is unknown whether Mn nodules exist at all or if different geochemical conditions are present within the pits. Similar, but larger, structures exist in the disturbance and recolonization experiment (DISCOL) area (Greinert, 2015), showing very similar geochemical conditions as other Mn-covered areas in both highly detailed sediment surface analyses as well as deeper sediment cores.

Comparison between sub-areas A1 and A2
When comparing the relationships between the bathymetric derivatives and the Mn-nodule coverage, it becomes evident that correlations visible in A2 cannot be seen in A1 (Fig. 6a, c), where areas of lower nodule coverage could not be matched with distinct terrain types. This result points towards additional parameters that influence Mn-nodule occurrence. Geochemical processes could be involved that drive the Mn-nodule formation; these in turn depend on the sediment properties (composition, sedimentation rate, porosity, etc.). Bottom currents could additionally influence the sedimentation rate and affect the geochemical processes in the benthic boundary layer and Mn-nodule surface. Local differences in the hydrodynamic regime near the bottom seem likely, as the bathymetric derivatives vary between the two sub-areas. Sub-area A2 is bound towards east and north by elevated terrain (7 to 10 m higher) which could have a focusing effect on bottom currents eventually causing a more erosional environment. In contrast, sub-area A1 is unbound by elevated terrain within 2 km distance. This might cause a stronger influence of seasonally changing bottom currents, preventing a clearer correlation of Mn-nodule coverage with the seafloor morphology.

Broad-vs. small-scale correlation
The observations made on a broad scale (several hundreds of meters; grid cell size of 55 m) show that high Mn-nodule coverage correlates with depressions (Fig. 5), which is consistent with observations on a smaller scale (scale of tens of meters; grid cell size of 5.5 m) for sub-area A2 (Fig. 7). Outside of A2 decreasing Mn-nodule coverage correlates with steeper sloping areas, which is contradicting to observations on a small scale, where the lowest Mn-nodule coverage correlates with extremely low slope angles of less than 1.8 • . This contradicting finding highlights that simple and generalized correlations between Mn-nodule occurrence and bathymetric but also geochemical properties in the sediment might not be possible on the regional scale (10 to 1000 km) but on the local scale (100 m-10 km). This is because the formation parameters also change on such local scales which are not possible to accurately predict using ship-based multibeam data, "sparse" box coring (distances of few kilometers) and limited information about current regimes.

Sediment plume resettling
The approach of conducting a photo mosaic survey before and after a seafloor disturbance proved successful for detecting sediment blanketing visually, offering the possibility to accurately map the area of strongest plume impact. This area is characterized by the sediment plume transport direction and resettling of the majority of the sediment. Very fine particles within the sediment plume might be dispersed much further; more detailed biological studies need to evaluate which sediment concentrations and grain sizes will impact benthic organisms on long timescales (cumulative effects) outside the visually clearly detectable impacted areas.
The thickness of the resettled sediments could not be determined from the AUV-based images or ROV-based video footage during the cruises. Video observations from other, similar areas point towards a sediment cover on millimeter or sub-millimeter range that can still be detected in images (e.g., in laboratory experiments; Yamazaki et al., 1997).

Morphology influence on sediment transport
The extent of the visible sediment blanketing, that varies over several tens of meters, can be related to a focusing of the sediment plume settling or the prevention of it through smallscaled morphological changes in form of barriers (steeper slopes facing against the current) or the opening of plume transportation pathways (sloping terrain with the current). Varying terrain in general will modify the current regime near the bottom and thus the settling properties of the sediment plume; it might also enhance the interactions between the particles due to increased turbulence that might stimulate increased flocculation and thus scavenging of very small particles that otherwise would be much further distributed. The shorter transport of sediment in north-and southward directions from the EBS track along slope A2E implies that the transportation of the suspension load follows the slope downhill. In subsection A2W, where the terrain is very smooth (the relief changes by 1 to 2 m), a dependency of the sediment blanketing extent to structures of the undulating seafloor could still be observed. At the western end of subsection A2W, the east-facing slopes act as barrier for an undisturbed migration of the sediment plume with the bottom current towards the south. The spreading of the sediment blanketing is wider in the east of subsection A2W where the seafloor is al-most horizontal, before slightly dipping towards the east and into the Pit 1 structure. The slopes considered show angles of less than 2 • and the morphological variability is sometimes less than 1 m. More distinct features, like Pit 1 (Fig. 9), cause a more variable sediment plume dispersal. The sediment blanketing within this 2 m deep feature does not exceed the southward edge of the depression. The resuspended sediment seems trapped within this feature with possible additional suspension load coming from the neighboring eastward slope.

Estimation of plume height
In a first approach, we estimated the plume height generated by the EBS by considering the extent of the observed sediment blanketing and measured bottom current velocities at the time of the EBS deployment (31 mm s −1 ; measured by ADCP). Former models from the CCZ reported settling velocities of particles in a sediment plume in the range of 0.1 to 1 mm s −1 derived from visual and experimental data (Lavelle et al., 1981;Oebius et al., 2001). Preliminary results of particle size analysis from a comparable site within PA1 indicate a median grain size of 29 µm (Benjamin Gillard, Jacobs University Bremen, Germany, personal communication, 2016). Following Stokes' law and disregarding aggregation of the particles, the determined median particle size for the area would translate to sinking velocities of approximately 1 mm s −1 . Assuming an average dispersal width of 30-50 m downstream, as indicated for the A2W subsection ( Fig. 9), this would require a plume height of approximately 0.96 to 1.6 m. Aggregation processes leading to larger particle sizes are likely to occur, which, due to increased friction, would sink slower than similar-sized Stokes particles, but that would scavenge a substantial amount of very small particles (Laurenz Thomsen, Jacobs University Bremen, Germany, personal communication, 2016). As part of studies in the south Pacific DISCOL area, lander-based ADCP backscatter measurements detected a passing-by sediment plume induced by a similar EBS experiment as discussed here. These data indicate a plume height between 1.5 and 2 m (Greinert, 2015).

Implications for possible mining scenario
It can be assumed that, due to the higher turbulence caused by the deployment of an industrial collector system and the continuous release of suspended material into the water column during mining, the dynamic behavior of the sediment plume could be altered and adjusted in such a way that the suspended sediment is resettling in the fastest possible way, keeping the dispersion to a minimum. Determining the dynamic behavior of the plume under different collectordispersion scenarios by monitoring in situ and under real mining conditions is thus essential to improve our understanding and model capacity with regards to the near-and far-field plume distribution and finally to evaluate ecological short-and long-term impacts.
These ecological impacts can be significantly confined to a small area by reducing the height of the sediment plume, increasing the settling velocity and aggregation of particles (scavenging the very fine sediment fraction). Vertical discharge of sediment after its separation from the Mn nodules should be avoided; instead a horizontal discharge close to the bottom (< 10 m from the bottom; below stable stratification above the well-mixed bottom boundary layer) with a velocity as slow as possible (speed of the collector) should be aimed for. One first implementation of this concept was the setup of the Deep-Sea Sediment Resuspension System (DSSRS) disturber (Brockett and Richards, 1994) deployed in a few large-scale benthic impact experiments (BIE-II, JET, IOM-BIE, INDEX).
As indicated by our results, a low-height sediment plume will be trapped in small depressions. Thus, detailed knowledge of the local morphology on small scales is a prerequisite to correctly determine the area and thickness of resettling sediment. This is also relevant in planning adjacent mining tracks from a miner's point of view, since strong sediment blanketing might bury adjacent nodules to be mined. According to our results, this impact will be highest in sediment accumulation sites, but even on flat areas with slopes of less than 3 • the distribution of the sediment plume and the resulting sediment blanketing distance will vary on a range of several tens of meters. In areas with steeper slopes (e.g., 10 • ), the sediment blanketing distance can be even wider.
In our very small-scale experiment, the EBS created a local impact with clearly visible sediment blanketing within 100 m downstream off the track. This localized impact is also the result of only partial resuspension of the surface sediments that was directly caused by the EBS (1.2 m in width). Observations of EBS tracks during another experiment revealed that a larger part of the sediment is compressed by the EBS and pushed aside with only a smaller (unknown) fraction being suspended (Boetius, 2015). It can be speculated that resedimentation of particles outside the visible blanketing area is minor, will happen over longer time and thus might not have a significant effect on the benthic organisms and the ecosystem (short-and long-term cumulative impacts on specific fauna still need to be determined).
The actual scenario of disturbance will be different during real-case mining during which the top 10-20 cm of the sediment are removed, then filtered for nodules and then discharged at the seafloor. One single track will be about 17 m wide as, e.g., planned in a German concept , whereas the track width of the EBS was only 1.2 m. As not only one track will be mined, but the collector system will operate constantly in a lawn-mowing pattern of long tracks scraping off the seafloor surface, the entire mined area will see a strong impact (Jankowski and Zielke, 2001). Considering local topography, bottom currents, optimizing particle settling (You, 2004) for fast and effective flocculation by the collector, and the cleaning of the sediment plume from the water column by settling phytodetritus from plankton blooms (increased flocculation), the size of the impacted area and the impact itself caused by the sediment blanketing outside the mined area might be rather small (< 10 km) and controllable. For a final validation, an experimental setup closer to the expected mining conditions is needed (Sharma, 2011); the presented study shows that we have the understanding, tools and the methodologies at hand to perform monitoring studies needed for such a realistic deep-sea mining experiment.

Conclusions
We conclude that, for both of our study topics, the Mn-nodule distribution to terrain comparison as well as the redeposition of sediments indicate that Mn-nodule coverage and sediment blanketing vary measurably on a very small scale (several tens to hundreds of meters), even if the seafloor terrain changes are minor (less than 1 m vertical change). This supports the second conclusion that spatial scale needs to be considered when discussing possible parameters that influence Mn-nodule coverage as such, and that relations found in one region most likely cannot be generalized to other regions and across different scales. Confirming previous studies, our data also show no simple relationship between Mn-nodule coverage and the seafloor morphology even when working on the same spatial scale. It needs to be realized that a complex interrelation between morphological characteristics and local environmental conditions (physical, chemical, sedimentological) influences the visually detectable Mn-nodule coverage at the seafloor surface. Reasons for this are that (1) variable amounts of Mn nodules have formed under different geochemical, bottom current or sedimentological conditions in different places, (2) Mn nodules might have dissolved in certain areas because of changing geochemical conditions or mechanical erosion, (3) Mn nodules were buried by sediment, whereas the sediment deposition pattern is influenced by the seafloor terrain and its interplay with bottom currents, and (4) the existence of Mn nodules (abundance, size, total coverage) itself influences sediment erosion, making denser covered areas with large nodules more resistant against sediment erosion.
With respect to the sediment plume study, it became obvious that a visible blanketing occurs in a limited distance (here < 100 m) away from the disturbance track and that the blanketing pattern strongly depends on bottom current direction, strength, small-scale bathymetry and initial plume height. From these observations, it can be concluded that each sediment plume disposal via an exhaust/diffusor of the collector should occur horizontally as close to the bottom as possible, rather than on top of the vehicle, blowing the sediment particles high into the water column, aiming for a finer dispersal over larger areas. Our studies also highlight that the per-formed disturbance experiment cannot be scaled up to a real mining scenario and that more detailed studies are required to understand and quantify the cumulative impact of unsettled particles on filter-feeding organisms beyond the clearly visual blanketing area.
On the technical side, the study showed that we have the needed tools and techniques at hand to map the seafloor for Mn-nodule resource assessments and a better understanding of Mn-nodule distribution, as well as for assessing mining impacts visually. It became clear that without such highresolution techniques valid assessments cannot be carried out. Areas that appeared suitable for mining (slopes ≤ 3 • ) in ship-based bathymetric data showed steeper relief (slopes > 3 • ) in higher-resolution AUV-based data. For an environmentally friendly deep-sea mining, such high-resolving maps are a prerequisite to accurately define areas that need to be protected and maneuver mining infrastructure around them considering the actual bottom currents and sediment settling areas during the mining.
Code and data availability. Source code for the automated nodule delineation is available in Pangaea (Schoening, 2017c). The data used in this work are available at Pangaea. This includes MBES data (Greinert, 2016), optical imagery  and image-derived nodule coverage (Schoening, 2017a, b).

A1 Calculation of the bathymetric derivatives
The slope was calculated using the algorithm included in the spatial analyst toolbox (Burrough, 1986) of ArcGIS (Table A4 output in degrees). Curvature as a second-order derivative of the bathymetry represents the slope of the slope. It has also been determined with the ArcGIS spatial analyst toolbox. For each cell, a fourth-order polynomial is fit to a surface composed of a 3 × 3 cell window. From this surface, the tool calculates the coefficients (Table A3), which are set into relation with the elevation values for every cell (Zevenbergen and Thorne, 1987). The two maximum slopedependent curvature values of the plan curvature, perpendicular to the maximum slope, and the profile curvature parallel to the maximum slope direction were calculated. The plan curvature defines flow convergence (concave surface, values < 0) and divergence (convex values > 0). Profile curvature affects the acceleration of the flow with values < 0 indicating a concave-shaped surface and values > 0 indicating convexshaped surface. In addition, the total curvature of a surface has been calculated which is also > 0 when convex shaped or < 0 when concave shaped. Curvature values of 0 are indicative of flat surfaces.
For calculating aspect, BPI and terrain ruggedness, the Ar-cGIS benthic terrain modeler (BTM) add-in (Wright et al., 2012) and its incorporated algorithms were used. The classification performed by the BTM is based on manually set properties (Tables A1 and A2) of the derivatives' slope, BPI (fine scale and broad scale) and water depth. This simple classification process provides sufficient information to distinguish different terrain settings of the study area (Figs. 2 and 3). The classification of the AUV-mapped study area (Fig. 6, Table A5) was performed to reveal areas of lower Mn-nodule coverage and is based on the derivatives BPI, slope and plan curvature, since these morphological parameters showed the best correlation with Mn-nodule coverage (Fig. 7). For the BPI calculation of the AUV-mapped area, the algorithm used within the BTM was modified (without integer rounding) to preserve the small-scaled features (Wilson et al., 2007).
The aspect is defined as the inclination direction of the maximum rate of change in depth from each cell to its neighbors, the slope inclination (Burrough, 1986). The algorithm calculates an aspect value for each cell of a raster and incorporates the respective adjacent cells in both horizontal directions from the center cell (dz/dx and dz/dy) ( Table A4).
The BPI describes the relative topographic variability of a central grid cell to a circular annulus with an inner and outer radius; both are manually defined (Table A1) (Weiss, 2001;Wright et al., 2012). For classification, a broad BPI (BBPI; large radii) and a fine-scale BPI (FBPI, small radii) are calculated and standardized. Positive values indicate that the central grid cell is elevated with respect to the mean annulus height values; negative values indicate depressions. The BPI is usually subscribed with the applied scale factor (grid cell size multiplied by outer annulus radius).
The terrain ruggedness was calculated for the AUV bathymetric data set (Fig. A4) using the algorithm for the vector terrain measurement (VRM) of the BTM (Sappington et al., 2007). It incorporates slope and aspect heterogeneity of the terrain and is defined as the magnitude of a resultant normalized vector from the decomposed x, y, z components of the cells and their slope and aspect, normalized to the number of cells in the neighborhood (Sappington et al., 2007). The terrain ruggedness is a unitless measure, ranging from 0 (flat) to 1 (most rugged). All derivative results were displayed and evaluated in ArcGIS 10.2.

A2 Interpretation of Mn-nodule size results
Considering a potential error in correctly detecting nodules by the CoMoNoD algorithm, the application of quantiles of the size distribution allows a more robust interpretation of the data. It is suggested not to use size values of the smallest and largest 1 % of the quantile calculation due to the abovementioned error source. The graph in Fig. A1 illustrates the quantiles of the calculated sizes of two images, which clearly differ from each other. The graph correctly displays a size difference between both images, indicating a greater number of larger nodules for image 29 302. This shows that applying CoMoNoD to calculate nodule sizes is reasonable. In this case, the best differentiation exists for the 50-75 % quantile. Towards larger and smaller size values the two curves approach each other which points towards the detection of similar -non-nodule -features in both images. Therefore, the median size values are considered to best represent the Mnnodule size distribution differences between images/areas. Without ground-truth data from sampling, computed size values should not be used as absolute values for resource assessment. However, they can be used to quantify nodule size distributions within seafloor areas and hence to compare variations in nodule distribution and coverage. Figure A1. Two example images (bottom) which clearly differ in nodule size and coverage. The graph shows the size distribution as calculated by the CoMoNoD algorithm. The most significant difference is observed in the 75 % quantile. Figure A2. Scatterplots indicating the relation between Mn-nodule coverage and total curvature (a, b) and profile curvature (c, d) in the A1 and A2 sub-areas. Only in area A2 weak correlations could be observed.      (2001)  Curvature Spatial analyst K = ((δ 2 Z/δS 2 )/(1 + δZ/δS) 2 ) 3/2 Zevenbergen and Thorne (1987) Plan curvature Spatial analyst 2 × (Dsin 2 θ + Ecos 2 θ − F sin θ cos θ) Zevenbergen and Thorne (1987) Profile curvature Spatial analyst −2 × (Dcos 2 θ + Esin 2 θ − F sin θ cos θ) Zevenbergen and Thorne (1987)