Quantifying Coral Reef Composition of Recreational Diving Sites: A Structure from Motion Approach at Seascape Scale

Recreational diving is known to have both direct and indirect impacts on coral habitats. Direct impacts include increasing sedimentation, breaks and diseases that lead to a decrease in the richness and abundances of hard corals. Indirect impacts include urban development, land management and sewage disposal. The ecological effects of scuba diving on the spatial composition metrics of reef benthic communities are less well studied, and they have not been investigated at seascape scale. In this study, we combine orthomosaics derived from Structure from Motion (SfM) photogrammetry and data-mining techniques to study the spatial composition of reef benthic communities of recreational diving sites at seascape scale (>25 m 2 ). The study focuses on the case study area of Ponta do Ouro Partial Marine Reserve (Mozambique). Results showed that scuba-diving resistant taxa (i.e., sponges and algae) were abundant at small (>850 m 2 ) and highly dived sites (>3000 dives yr − 1 ), characterized by low diversity and density, and big organisms with complex shapes. Fragile taxa (i.e., Acropora spp.) were abundant at low (365 dives yr − 1 ) and moderately dived sites (1000–3000 dives yr − 1 ) where the greater depth and wider coral reef surfaces attenuate the abrasive effect of waves and re-suspended sediments. Highest taxa diversity and density, and lowest abundance of resistant taxa were recorded at large (>2000 m 2 ) and rarely dived sites. This study highlights the potential applications for a photogrammetric approach to support monitoring programs at Ponta do Ouro Partial Marine Reserve (Mozambique), and provides some insight to understand the influence of scuba diving on benthic communities.


Introduction
Coral reefs across the world are facing multiple threats with anthropogenic activities playing a key role in the type and magnitude of the impacts. Among local anthropogenic threats, scuba diving directly affects the composition of benthic reef communities by increasing sedimentation, breaks and diseases [1,2]. Previous studies on coral reef composition have focused on high-resolution taxonomic the complexity of organism shapes across a seascape and their fractal dimensions [34]. Therefore, investigating the spatial composition of reefs using SfM could greatly contribute to our understanding of the reef functioning by helping to assess the pressures that resulted in the dimensions, shapes, abundance richness observed at a seascape scale.
Therefore, this study investigated the spatial composition of benthic reef communities at a seascape scale (>75 m 2 ). We used high-resolution orthomosaics generated with SfM photogrammetry to examine the differences in composition across recreational diving sites. We hypothesized that (i) seascape-based metrics would differ across sites as in Table 1, and (ii) taxonomical and functional benthic groups would differ across sites, with greater abundance of resistant organisms to sites with high levels of diving activity and to small and shallow reefs. Table 1. Seascape metrics considered to be dependent on anthropogenic pressures and environmental conditions. "Patch" refers to each of the single organisms within a digitized polygon (Section 2.2.3). "*" indicates seascape metrics calculated on the virtual random quadrats using FRAGSTATS v4 and their description [34].  Indicates the probability that two entities, (i.e., pixels) taken at random from the same seascape, belong to different patch types.
Large values indicate high probability that two pixels are from different patch types [33]. Indicates the complexity of the organisms across an area. e.g., if small and large patches have similar and simple geometric shapes, the index will be small, indicating that as the patch area increases, the patch perimeter increases too but by small increment [34].

Study Area and Experimental Design
The study was conducted in the Ponta do Ouro Partial Marine Reserve (PPMR) (Mozambique). The protected area includes some of the southernmost coral reefs of the African continent (Figure 1), and is characterized by low profile outcrops originating from the fossilization of sand dunes that runs parallel to the coastline, and are covered by a non-accretive "coral carpet" and high numbers of endemic species [35][36][37][38][39]. High energy hydrodynamic conditions, characterizing the North-South coastline between Ilha da Inhaca and Ponta do Ouro, shape local topography and sediment transport [37,39,40]. As a result, studies have reported negligible impacts of climate change on Maputaland reefs compared to other reefs across the tropics, because of local hydrodynamic conditions (i.e., Agulhas current, local upwellings, and a swell-generated inshore counterflow) preventing extreme water heating [41][42][43]. Communities are often dominated by soft corals that are tolerant to strong wave energy and sediment re-suspension [19]. The PPMR is an international destination for scuba divers and nature tourism since the end of the civil war in the 1990s when 30k to 40k dives year −1 were recorded the 1995 [40], 80k-90k dives year −1 were recorded in 1998 [44] with 42.5k-62k dives year −1 were counted in 2001-2002 [45]. More recently, an average of 26k (± 3k) dives year −1 were reported by the PPMR (Table 2). In this study, seven dive sites were selected within a range of average depths comprised between −12 m to −25 m below mean sea level (BMSL), and areas between 425 m 2 and 3500 m 2 . The number of divers year −1 ranged from none to >11k reported in the literature for 2001-2002 and 2011-2016 ( Figure 1, Table 3).
High-resolution data were collected at each site in one single dive between April and June 2015 covering the entire coral reef area or surveying those areas commonly visited by divers at the larger study sites. Dive sites were mapped as in [23] using a modified Diver Propulsion Vehicle (DPV) equipped with action cameras (GoPro Hero3 Silver-Woodman Labs, Inc., San Mateo, CA, USA), a tablet (Asus Google Nexus 7 tablet-Google Inc., Mountain View, CA, USA), an illumination system and one GPS buoy (Garmin Ltd., Lenexa, KS, USA) ( Figure 2). A detailed description of the vehicle in Palma et al. [23]. The camera was set to time-lapse mode, recording nadir images at 1 Hz frequency with a focal length equivalent to 21 mm. The position of each recorded image was estimated by coupling the camera recording time and the GPS clock [46]. The pitch, yaw and roll data of the sampling device were recorded logging the inertial navigation data obtained from the gyroscope and the accelerometer of the tablet and by using the UBICA Underwater Position System (UUPS) [47]. The geographic position of the sampling was recorded via a towed buoy with an integrated Etrex10 GPS (Garmin Ltd., Lenexa, KS, USA). PixGPS (BR Software, Asker, Norway) was used to tag the recorded images with the GPS positioning and the rotation data. The coordinates for each of the frames were used to georeference (scale, translate and rotate) the imagery into the coordinate system defined by the World Geodetic System (WGS84) and to minimize geometric distortions [48]. The Wide Area Augmentation System (WS) mode at 1 Hz frequency was used to record the position of the buoy. The diver maintained an average swimming speed of 0.75 m s −1 and an average distance to the sea bottom of 2.7 m, with the range being between 1.5 m and 3 m and with each frame covering 3 m 2 ( Figure 2a). The sampling path was adapted to the seascape characteristics to maximize overlapping images, by carrying out both parallel and perpendicular transects . Table 2. Dives per year and site. "*" data provided by [49]." **" data provided by the management authority of the Ponta do Ouro Partial Marine Reserve (Mozambique), 18    For each site, high-resolution orthomosaics (pixel size < 0.01 m) were obtained by processing the geo-tagged images of the underwater surveys as in [23]. The Images were aligned to generate spare point cloud setting Photoscan version 1.4.1 (Agisoft LLC., St. Petersburg, Russia) [23] to high accuracy, generic pair selection, limiting the key points identification to 50,000 common features and unlimiting the number of the outcoming tie points cloud. Meshes were calculated setting the software surface type generation to height field, high face count and enabling interpolations. Adaptive orthomosaics texturing was used to export tiled (4096 pixels × 4096 pixels) georeferenced orthomosaics using the projected reference system Universal Transverse Mercator (UTM) fuse 36 Southern Hemisphere, defined by the World Geodetic System (WGS84).
First, the reef area belonging to solid substrate was digitalized using ArcGIS 10.1 software (Redlands, CA, USA). Then, to study the composition of the coral reefs at seascape scale, virtual random quadrats of 5 m × 5 m were randomly positioned three times over each ortho-mosaic ( Figure 2b). The quadrat size 5 m × 5 m has proved to be representative of the coral reefs at the PPMR [23], and allowed to sample a minimum of 80% of hard substrates in all quadrat replicates and reef sites. The centroids of the virtual quadrats were distributed random in the reef sites using the ArcGIS tool "Create Random Points" (Redlands, CA, USA). We discarded the quadrats with more than 20% of their area outside the digitalized reef area [23]. Within each quadrat, the benthic organisms were identified in taxonomic and morphological classes, digitalized using ArcGIS 10.1 software (Redlands, CA, USA), and finally rasterized to obtain categorical maps (Table 4, Figure 2c1) [23].
The classes were summarized in functional groups (FGs, [19,50]) and for each class, abundance and the summary statistics of surface and perimeter were calculated and used for multivariate analysis (Figure 2(c2,c3)) (Section 2.2.2). FGs classes are based on a series of taxonomic, morphological, physiological, and behavioral characteristics reflecting adaptations to physical disturbance. "Resistant" to physical impacts organisms included sponges (SP), algal patches (TA) and "other invertebrates" (OI) (e.g., gastropods, tunicates, echinoderms, hexacorals Stichodactyla spp., Atriolum sp.). These taxa show high regenerative physiological resistance to physical damages because of their skeletal structure and have been often reported in relation to ecological shifts in coral reefs communities [2,[51][52][53]. The "moderately resistant" organisms to physical impacts included soft corals (SCD, SCP, SCC, SCR, SCF) (e.g., Sarcophyton sp., Sinularia sp., Lobophyton spp.), massive corals (CM) (e.g., Platygyra spp., Montastrea spp., Galaxea spp., Favites spp., Favia spp. and Turbinaria sp.), folios corals (CF) (e.g., Montipora spp., Echinopora spp.), and encrusting corals (CE) (e.g., Porites spp.). Soft corals are often described as sensitive organisms to sediment burying and wave actions [54] but their flexibility of the colonies, their low profile suggest they are resistant animals to hydrodynamic condition [19]. Massive, encrusting and folios corals show resistance to physical impacts due to a stout or flat morphological growth. However, they can be dislodged or more commonly scratched by accidental impact with divers equipment [2,13,55]. Finally, the "fragile" "group included Acropora spp. (AC), branched corals (CB) (e.g., Pocillopora spp.), and tabular corals (CTA) (i.e., Turbinaria sp.). Acropora spp. and branched-coral taxa have tree-shaped carbonate skeletons that increase the structural complexity of reefs but, as erected forms, they are also the most exposed to breaks and fragmentation [2,11,12,[55][56][57]. Among the fragile group, tabular corals and free-living fungiid colonies are the most subject to breaks, for example during severe weather conditions entire colonies might be tipped over [58].  [23]; (b) processed orthoimage of the dive site Blacks with the example of three virtual quadrats (5 m × 5 m) randomly positioned across the reef area; (c) GIS-based analysis of virtual sampling (i, iv), digitalization and taxonomy analysis of the benthic organisms (ii, iii, v), data rasterization (vi); (d) summary of spatial data processing of (1) clustering analysis and (2) seascape index calculation; (e) interpretation of the resulting clusters accordingly to (1) scuba divers pressures and (2) environmental pressures. Table 4. Taxonomic and morphological classes identified from the orthomosaics for each site. The table also reports the classes of resistance obtained from the analysis of the quadrats.

Resistant-to-Physical-Impact Categories Functional Groups (FGs) Classes Code
Fragile Acropora spp. Data on FGs (i.e., diversity, surfaces and perimeter, average, maximal, minimal, standard deviation, count), and reef sites (i.e., depth, extension) were analyzed using a principal component analysis (PCA- [59]) and the distance-based clustering X-Means [60] to look for similarities in community composition between diving sites and possible diving impacts (Figure 2d,e). Data of number of dives per reef site were not included in the analysis. The approach, which involves a series of steps performed sequentially to the harmonized data was applied at seascape scale and carried out using the DataMiner system of the D4Science e-Infrastructure (Italy) [61,62].
First, a PCA was performed to reduce the dimensionality of the dataset and to explore the patterns in data variability among diving sites [59,63]. Secondly, X-Means was applied to the PCA-transformed vectors associated with the largest variance, to identify the variables indicating similarities of benthic composition between sites [63]. In this study, we searched for several clusters between 1 and 10, and applied a tolerance threshold of 25% over the centroids of the clusters that represent a summary of the characteristics of the cluster [63]. The clusters were evaluated with respect to expert opinion. Two independent experts evaluated the clusters by assigning the clusters' results for each diving site. A confusion matrix was used to assess the agreement between the experts' assessments and the X-Means clusters as the percentage of matching assignments (absolute percentage of agreement). The Cohen's Kappa [64] was calculated to estimate the agreement between the experts and the model compared to purely random assignments [65]. To look at patterns of taxonomic and functional benthic composition over the clusters we used the categories in Table 4.

Estimation of Seascape Metrics
The composition of the reefs was investigated by calculating key seascape metrics such as diversity, aggregation and shape on the categorical maps using the FRAGSTATS v4 software (Amherst, MA, USA) [34]. The set of metrics was selected based on their relevance to seascape ecology composition. FRAGSTATS v4 computes a variety of metrics for categorical maps and has been used in a wide range of disciplines including seascape ecology [23,[66][67][68][69]. In this study, the categorical maps were processed at patch (the single individual digitized in a polygon), class (taxonomical and morphological FGs, Table 4), and seascape level (denotes all the patches of different classes within the quadrats) by considering the disposition, the size and the shape of the patches within the quadrats. Patch density, Simpson Diversity Index, Patch Richness, Shannon's diversity Index, Perimeter-Area Fractal Dimension, and Mean patch fractal dimension at seascape level, were compared between reef sites by applying a One-Way ANOVA test using the aov function in the R "stats" v3.5.1 package [70].
To address the ANOVA assumptions, particularly the assumption of heterogeneity in the residuals, Patch Density, Patch Richness, Simpson Diversity Index, Perimeter-Area Fractal Dimension were log-transformed, but this was not necessary for the rest of seascape metrics ( Table 1). Validation of underlying assumptions of normality and homoscedasticity of tests residuals was applied following Zuur et al. [71] (shapiro.test function from the R "stats" v3.5.1 package [72] and leveneTest in R "car" v3.0-2 package [73]) and subsequently, post-hoc Tukey tests were applied to compare which specific treatments differ significantly using the TukeyHSD function from the R (New Zealand) "stats" v3.5.1 package [72].

Seascape Composition of Diving Sites
Thousands of images were processed to generate the orthomosaics at each diving sites. 2102 images were aligned over 2186 for Doodle, 1165 images were aligned over 1190 for Blacks, 2108 images were aligned over 2165 for Steps, 1758 images were aligned over 1780 for Maverick, 2444 images were aligned over 2475 for Paradise Ledge, 1490 images were aligned over 1651 for Texas, and 3662 images were aligned over 3394 for Techobanine. From the analysis and digitalization of the orthomosaics, 11 functional groups were identified (Table 4). Doodles and Blacks showed the higher frequencies of sponges and algal patches among sites and lower abundances of organisms with <0.5 individual m −2 (Table 5). Steps presented an average density of 2.37 individual m −2 and was characterized by soft digitate and crested corals (SCD, SCC) and several classes of invertebrates (OI) including acidians and anemons. Branched coral (i.e., Pocillopora sp.) and Acropora were also present at site (Table 5). At Maverick and Paradise Ledge, organism densities averaged between 4.4 individual m −2 and 4.85 individual m −2 and soft corals (SCD, SCC, SCF, SCR) dominate the seascape. Acropora (ACS, ACD) and non-Acropora branched corals with folios and massive hard corals were also present ( Table 5). Texas has the highest average density with 6.47 individual m −2 mainly represented by branched corals (CB), soft corals (SCD, SCC, SCR, SCF) and encrusting sponges (SPE). Tabular corals (CTU), Acroporas (ACS) and alges (TA) were also present (Table 5). Finally, Techobanine recorded 3.56 individual m −2 mainly represented by branched corals (ACC, CB, ACS), massive and tabular corals (CM, CTU) and soft corals (SCD, SCC, SCF) ( Table 5).

Evaluation of Clusters and Agreement with the Experts
From the analysis, the reef was classified in three clusters, labelled, semantically described (Table 6), and evaluated with respect to expert opinion (Table 7). Clusters evaluation showed moderate to a substantial agreement with expert opinion, indicating reliable semantic interpretations of the vectors identified in the clusters variations ( Table 7). The marginal agreement is recorded between the two experts due to Expert 2 who classified most of the sites as moderately impacted.
Cluster number 1 groups very impacted sites (Doodles and Blacks) and contains 23 vectors. It indicates the whole set of benthic taxa that are very resistant and resilient to physical pressures (i.e., algae, sponges), with low density and diversity, and generally more complex morphological shapes per organism.
Moderately impacted cluster (Cluster 2) contains 33 vectors and three sites (Steps, Techobanine and Texas). It represents the set of taxa that are low resistant but highly resilient to pressures, i.e., branched hard coral. They are sensitive to breaks but they have high growth rates [74] and show higher density and diversity, lower morphological complexity per organism and higher fractal complexity at seascape scale than cluster 1 (Figures 3 and 4).
Finally, the low impacted cluster (cluster 3) describes 32 vectors and two sites (Paradise ledge and Maverick). It refers to taxa that are moderately resistant and low resilient to physical pressures (i.e., soft corals that have low growth rate and low recruitment rate but highly resistant to high energy hydrodynamic conditions [19]), show the highest level of taxa diversity and density, and a high level of fractal complexity at seascape scale. Table 6. X-Means clusters identified in the multivariate analysis, labelled and semantically described.

Clusters
Label Description  Very impacted sites (cluster 1) showed the greatest patch surfaces, and similar distribution to moderately impacted sites (cluster 2, Figures 3a and 4). Doodles and Blacks (cluster 1) were also characterized by patches with larger perimeters and higher variability in the first and third quartile than Steps, Techobanine, Texas (cluster 2), Paradise ledge and Maverick (cluster 3) (Figures 3b and 4). In fact, sites in clusters 2 and 3 were characterized by patches with smaller perimeters and lower variabilities (Figures 3b and 4). The low impacted sites (cluster 3) showed the highest abundance of benthic organisms and high variability among the two sites (Figures 3c and 4). Conversely, in very and moderately impacted sites, organism abundances have overall similar distributions among sites (cluster 1 and 2, Figures 3a,c and 4). Across sites, the resistant group accounted for >85% of the patches surface, the 80% of their perimeter and more than >60% in abundances in the very impacted cluster (Figure 5a,b). Moderately resistant organisms constituted only the 10% of the patches surface at very impacted cluster, while contributed >30% in the moderately and more than 60% in the low impacted clusters (Figure 5a). The moderately resistant group had a similar pattern also for patches perimeter and abundances (Figure 5b,c).
Among sites, the resistant group had the highest contribution in percentage to mean patch areas and abundances to Blacks (95%), Doodle (>75%), and Steps (>50%) (Figures 6 and 7). The moderately fragile group mostly accounted in percentage to mean patch areas and abundances to Paradise Ledge (50%), Maverick (70%), and Texas (60%) (Figures 7 and 8). Finally, Paradise Ledge and Techobanine showed the highest contribution in percentage to mean patch areas (45% in Paradise Ledge, 70% in Techobanine) of the fragile organisms group (Figures 7 and 8).

Seascape Metrics
Most of the seascape metrics were significantly different across reefs (i.e., patch richness, patch density, Shannon's diversity index and perimeter-area fractal dimension), except for the Simpson Diversity Index and the mean patch fractal dimension index. These indices indicated a similar average complexity of organisms across the seascapes. (Figure 9). Patch density, patch richness and Shannon's Diversity Index lowered at very impacted sites (cluster 1, Blacks and Doodle) (Figure 9). Patch density showed significant differences (F(6,14) = 23.04, p < 0.01) between sites with larger seascape complexity at Paradise Ledge (52 ± 1.8 organisms) and lower at Black (6 ± 1.0 patches) and Doodles (6 ± 1.5 patches) (Figure 9). Patch richness exhibited significant differences in the number of classes of organisms between sites (F(6,14) = 37.40, p < 0.01), with five classes at Black 5 (±0.4) and Doodle 6 (±0.4) and over 10 in the other diving sites (Figure 9). SDHI index indicates that Maverick and Techobanine have higher number of classes of organisms, equally distributed within the sites, than the other sites ( Figure 9). Finally, Doodle has the highest morphological complexity of the patches (2.23 ± 0.54) associated (perimeter-area fractal dimension index F(6,14) = 3.54, p = 0.02) (Figure 9).

Discussion
This research assessed the spatial composition of benthic reef communities at seascape scale and within recreational diving sites using SfM orthomosaics. Coral reefs exhibited different spatial compositions of benthic communities, with patterns emerging between diving sites. In this section, we discuss the variation of taxonomic and functional descriptors across diving sites (Section 4.1), and the potential applications of the proposed approach (Section 4.2).

Linking Spatial Benthic Composition to Recreational Diving Sites
Differences between sites were recorded by FGs and their metrics (i.e., surface, perimeter, abundance) in the clusters, indicating that scuba diving and local environmental conditions could influence the resultant spatial community structures and that specific taxonomic and functional descriptors provide further insight into the factors affecting coral reefs. The seascape indices did not show significant variability within sites suggesting that the selected sampling units and the number of replicates were suitable for describing the communities of the selected coral reefs. The resistant-to-physical-impact categories provided further insights into seascape composition, and partially support our second hypothesis that the contribution of resistant organism tend to occur in small, shallow and highly dived coral reef sites (Figure 9).
Dimensional factors, such as patch surface and perimeter are indicators of benthic assemblage composition and spatial complexity. The larger the organism, the higher is its planar complexity as a result of intra-and inter-species competition within the seascape [75]. The seascape composition factors (i.e., functional groups and abundance), are indicators of benthic taxa diversity and organisms density. Lower values in taxa diversity are related to more homogeneous seascapes, whereas higher values of organism densities suggested lower human impacts (i.e., low number of dives per year) and environmental pressures (i.e., deeper and larger coral reefs).
Low densities of organisms with large and complex shapes were recorded in Doodle and Blacks which were grouped in the very impacted cluster (cluster 1). These sites are characterized by lower taxa diversity and by the highest number of dives per year and are the closest sites to the launching point on Ponta do Ouro beach. Doodle has been highly frequented since 1995 while Blacks has been dived since 2011 with moderately intensities (Table 3). Both sites are subjected to strong environmental pressures due the number of dives per year in relation to small surface areas (450 m 2 to 850 m 2 ) at moderately shallow average depth (−17 m to −19 m). The waves and the abrasive action of re-suspended sediment limit the settlement and growth of resistant-to-physical-impact organisms such as sponges and algae. Previous studies have demonstrated that the high-amplitude swells originating in the Southern Ocean [35,55,76] promote suspension and deposition of sediment, mechanisms of abrasions that limit at the reef-sediments interface the recruitments of benthic larvae, also for resistant-to-physical pressure taxa [19]. Doodles and Blacks also presented the highest frequencies of organisms resistant to physical impacts (i.e., sponges, algal patches). The few fragile organisms identified in cluster 1 sites are tabular corals which have a limited three-dimensional complexity shape and are potentially more resistant to the physical impacts of scuba diving than Acroporas and similar branched corals.
Higher FGs diversity and a high presence of Acropora colonies were recorded for the moderately impacted cluster (cluster 2). This cluster included three sites very distant from each other ( Figure 1): Techobanine is at 19.6 km, Steps at 3.6 km and Texas at 9.6 km from the Ponta do Ouro launching point. Steps is characterized by the lowest organisms' density while Texas has the highest. Both sites present high variability which suggests that there are spatial differences in the benthic community composition (patch richness and patch density) over the sampled seascapes ( Figure 9). Texas also showed smaller average organism surface and higher fractal dimensions than the other two sites, suggesting greater variability in the shapes and dimensions of organisms that could be related to the habitat recovering from extreme environmental events (i.e., cyclones and storms) or to intense scuba diving of the particular site in the past ( Table 2). Techobanine showed the highest frequency of hard branched corals and colonies of Acroporas, probably favored by the lower average depth of the site (−12 m) and by the very low number of dives per years recorded in the last decade. The low density and large dimensions of the organisms is explained by the higher growth rate of these hard corals [77] which mitigates the impacts caused by the high energy hydrodynamic conditions of the area.
Among cluster 2, Steps is the only site with high number of scuba dives since 1995 (>3000 dives year −1 ), while Texas has less than 300 dives year −1 in the last decade and Techobanine is annually rarely visited (Tables 2 and 3). Steps is smaller than Techobanine and Texas, but it develops along a rocky edge. There are small canyons and cracks in the rocky substrate which provide sheltered areas against waves action and scuba diving direct impacts, resulting in a more heterogeneous community with moderately fragile and resistant organisms [78]. Cluster 3 showed higher level of taxa diversity and organism density probably as a result of the low scuba-diving pressure and the more benign environmental conditions. The reefs are characterized by several sizes of corals, small perimeter variability and low fractality within the landscape. Paradise ledge and Maverick are the deepest (−18 m to −20 m) and largest sites, with an average of 365 dive year −1 at Paradise Ledge in the last 12 years. Maverick is a rarely dived site close to Paradise Ledge.

Implications for Reef Management
This study proposes a new transferable approach based on SfM to study the spatial composition of coral communities and contributes towards the identification of reefs where diving activities are impacting on the community structure.
Those responsible for managing reefs are increasingly considering both the revenue derived from scuba diving as well as the impacts on the habitats in their management plans as it is one of the most profitable activities in marine parks [79]. Global policies are needed to face the challenges of global warming and consequently massive bleaching events on reefs [80]. However, more local actions are needed to preserve and manage the reefs in terms of the impacts caused by human activities. There is therefore a clear need to provide marine managers with cost effective and easily deployable approaches that can be applied to coral reefs.
Our method is the first to provide relevant, cost effective information to marine managers that will enable them to study the changes in composition and functionality across reefs at PPMR. It will provide an objective basis for conservation actions and long-term data collection programs.
Recently, the village of Ponta do Ouro has been connected via a new road to Maputo. The new road has reduced the travel time to less than two hours by car or buses while previously it was a 7 h journey in an all-terrain vehicle. Managers of the PPMR are expecting greater numbers of visitors with consequent negative impacts on the local natural habitats including the coral reefs (personal communication from the Director of the Ponta do Ouro Partial Marine Reserve). It will be important for the PPMR managers to establish and document a baseline condition survey of all the potentially affected reefs. This can be carried out quickly and cost effectively with the method we have developed. The seascape approach will also facilitate the development of a prioritized action plan. It will require the collection of quantitative, georeferenced, and representative samples, including compositional data to provide a functional perspective. The SfM-based approach applied in this way will leads to a better understanding of the short, medium and long-term implications of benthic changes to management actions.
Managers of the PPMR can benefit from the combined SfM-approach and statistical classification at any stage of the analysis. First looking at the clusters that describe: (i) taxonomic and morphologic composition of the benthic communities, (ii) seascape spatial composition and organism shape complexity, (iii) classes of fragility to physical impacts. Second by looking at the final matrices which embed the summary of the related ecological data and provide insights into the spatial variability of the coral reef benthic communities in the protected area. From the above information, managers of the PPMR can prioritize conservation measures at each site.
The information collected in this study has indicated the need for further data collection and analysis if the overall management of the reefs is to be effective and to ensure that the methodology developed is able to inform the most appropriate management approaches for the reefs. This includes: establishing a baseline before there is an increase in people visiting the area; identifying the magnitude of the various risks to the extent and condition of the reefs; distinguishing between the impacts caused by the various factors for example natural variability by season and by year, human impacts from diving, environmental pressures from natural events such as storms, climate change; determining the relative impact of the different factors for each location. This information would then enable appropriate risk management approaches to be developed. These could include: restricting the numbers of dives allowed at certain locations and during particular environmental conditions; allowing recovery periods for areas where no diving is allowed; access for beginners being restricted to particular areas; clear instructions and educational programs for divers and those running diving schools so that they are aware of and avoid certain activities known to cause damage [81]. To this end, information collected from scuba diving can support marine conservation in the context of marine citizen science [82]. Trained divers involved in monitoring programs (i.e., Reef Check) have reduced cost, time, equipment and simplified logistics [83,84]. In this direction, recreational diving represents a potential to continue monitoring and surveillance programs [82], and could be an opportunity to promote photogrammetry at seascape scale to ensure that the risk mitigation approaches are effective. The monitoring approach we have developed can play an essential role within the overall reef management approach. The monitoring approach needs to include: determining the baseline conditions; identifying and distinguishing between the risks to the condition of the reefs from all sources; introducing appropriate reef management actions; periodic monitoring to ensure that the management measures are being effective and modifying them as required; long-term monitoring to assess natural variability and the impact of underlying long-term environmental changes and pressures.

Conclusions
This study has reported the application of SfM-based approach to marine habitats and its advantages. SfM demonstrated to be a flexible, affordable and reliable technology for rapid surveys, providing informative data to support spatial analysis of ecological communities and to investigate the effects of human stressors. Our results showed: • substantial agreement with expert opinions; • differences in the benthic composition between diving sites in terms of: (i) organism diversity, surface, perimeter, abundance, (ii) functional groups categories (e.g., resistant, moderately resistant, fragile); • Resistant-to-physical-impact categories (i.e., sponges and algae) were most abundant at small diving sites and highly visited by divers. These sites were characterized by big organisms and with complex shapes, low taxa diversity and density.

•
Fragile-to-physical-impact categories (i.e., Acropora spp.) were most abundant at diving sites low or moderately visited by divers. These sites are characterized by complex geomorphology and moderate hydrodynamic conditions. • The highest taxa diversity and density, and the lowest abundance of resistant-to-physical-impact FGs were recorded at large and rarely dived sites. The sites exhibited different seascape metrics (i.e., patch density, patch richness, Shannon Diversity Index, Perimeter-Area Fractal Index), with general patterns emerging in terms of responses to diving pressure and environmental conditions.
To conclude, this study suggests that having effective coral ecosystem management systems in place at PPMR will be essential if these important marine habitats are to survive and flourish. This is especially the case with the potential significant increase in human impacts arising from the increasing accessibility of the area to visitors. By focusing on strengthening the resilience of marine ecosystems at a seascape scale, marine management can optimize goods and services as new environmental challenges arise [85,86].
Author Contributions: M.P. conceived, designed and analyzed the datasets. C.M. contributed to the data analysis and writing of the manuscript. M.R.C. and P.L. contributed to editing and revising the manuscript. U.P. and J.F. assessed the clusters. G.C. and C.C. provided their statistical and ecological perspectives and revised the manuscript. All authors read and approved the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.