High spatial resolution photo mosaicking for the monitoring of coralligenous reefs

Coralligenous reefs are characterized by large bathymetric and spatial distribution, as well as heterogeneity; in shallow environments, they develop mainly on vertical and sub-vertical rocky walls. Mainly diver-based techniques are carried out to gain detailed information on such habitats. Here, we propose a non-destructive and multi-purpose photo mosaicking method to study and monitor coralligenous reefs developing on vertical walls. High-pixel resolution images using three different commercial cameras were acquired on a 10 m2 reef, to compare the effectiveness of photomosaic method to the traditional photoquadrats technique in quantifying the coralligenous assemblage. Results showed very high spatial resolution and accuracy among the photomosaic acquired with different cameras and no significant differences with photoquadrats in assessing the assemblage composition. Despite the large difference in costs of each recording apparatus, little differences emerged from the assemblage characterization: through the analysis of the three photomosaics twelve taxa/morphological categories covered 97–99% of the sampled surface. Photo mosaicking represents a low-cost method that minimizes the time spent underwater by divers and capable of providing new opportunities for further studies on shallow coralligenous reefs.


Introduction
Photography represents an efficient and powerful way to disseminate scientific thoughts and findings both inside and outside the scientific community (Wilder 2009). It also constitutes an objective tool for measuring and describing natural phenomena, reducing human error often affecting in situ observations. In the last thirty years, the number of studies using photography and video as methodologies for data collection has drastically increased in marine science; species, communities, and habitats can be advantageously investigated through these techniques for different scopes (Bicknell et al. 2016;Durden et al. 2016;Bayley and Mogg 2019).
Benthic habitat mapping is a fundamental step for managing and preserving coastal zones, as well as examining unknown environments (Park et al. 2011). Among the underwater domains, there is a strong correlation between the spatial distribution of marine organisms and geomorphological characteristics and features of seabed, especially on hard substrates (e.g., rocky and biogenic environments). In fact, biogenic habitats usually show more diverse geological and biological characteristics compared to sandy nearshore environments (Meadows et al. 2012).
Coralligenous can be considered one of the most important and characteristic ecosystems of the Mediterranean Sea due to biodiversity, productivity, recreational and commercial values. The term coralligenous refers to  (Martin et al. 2014;Ingrosso et al. 2018). They are produced primarily by the accumulation of encrusting algae growing in dim light conditions, and secondarily by bio-constructor animals (bryozoans, serpulids, and corals) which help to develop and consolidate the bioconstruction (Ballesteros 2006). Coralligenous platforms and reefs are the two main morphologies described in the literature: The former are built over horizontal sedimentary substrates, the latter develop mainly on vertical rocky substrates located in shallower waters (Ballesteros 2006;Bracchi et al. 2015; SPA/RAC 2019). European directives stress the importance of assessing the ecological quality of coralligenous habitats; as for platforms monitoring is commonly carried out through the use of Remotely Operated Vehicles (ROVs) (Gennaro et al. 2020). Producing high-quality habitat mapping of reefs is not a trivial task because of their irregular bathymetric distribution and morphology which may result in operational constraints during mapping campaigns and sampling surveys. This aspect is especially true when coralligenous reefs are found on vertical walls and overhangs. In fact, in these complex habitats acoustic (multibeam and side-scan sonar) data derived from ship-based research surveys, ROVs, and Autonomous Underwater Vehicles (AUVs) cannot always provide exhaustive information at community level. In addition, even if the aforementioned remote techniques have been extensively employed to successfully map subhorizontal and deep bioconstructions (Brown et al. 2011;Bonacorsi et al. 2012;Bracchi et al. 2016), these vehicles have high operational costs, and the lack of availability and impracticality of deploying AUVs and ROVs leads many ecologists to rely on images taken by inexpensive and readily available hand-held cameras . In fact, on vertical biogenic outcrops, often located in shallower environments, diver-based photographic assessments remain a valuable method for community studies and environmental quality estimation. To date, visual assessment and video-photography are the most diffused sampling methods on coralligenous reefs: They increase the number of samples per survey, facilitate sampling over large spatial and temporal scales, and provide permanent records of the communities under study (Kipson et al. 2011;Zapata-Ramírez et al. 2013;Doxa et al. 2016;Tribot et al. 2016;Cerrano et al. 2019;Valisano et al. 2019). Conversely, the trade-off between making an adequate sampling effort and keeping acceptable safety measures is one of the major constraints associated with divers-based underwater sampling on coralligenous reefs (Parravicini et al. 2010). Undisturbed coralligenous reefs are characterized by high small spatial scale patchiness which implies a considerable minimum sampling effort to characterize the structure of the assemblage at a given site (Kipson et al. 2011;Piazzi et al. 2004Piazzi et al. , 2016Piazzi et al. , 2017Cecchi et al. 2014;Casas-Güell et al. 2016). In this framework, the application of a technique able to provide a detailed view of benthic assemblages with a relative time and cost efficiency will acquire a key role during the acquisition of ecological data and monitoring of such complex environments.
In the past decade, underwater photography techniques and diving equipment have been considerably improved by modern technological advances, leading to large development and spread of digital cameras coupled with a notable improvement in the realm of Computer Vision (e.g., Scale Invariant Feature Transform -SIFT-and Structure from Motion -SfM-algorithms), resulting in new powerful systems for processing underwater images. In particular, visually presenting information of a scene with an unbroken view of the whole region surrounding the observer by using panoramic images can provide to marine ecologists a useful tool to view and navigate underwater landscapes. Furthermore, it contributes to overcoming the light attenuation and backscatter limit imposed underwater. Most of the applications can be divided into two main categories: (1) those attempting to create a consistent continuous image (i.e., a two-dimensional photomosaic), possibly at the expense of minor local distortions; and (2) those aiming at the accurate recovery of three-dimensional information to extract quantitative data about the seafloor, principally based on photogrammetric approaches, which have a substantially higher level of complexity and costs. Two-dimensional photomosaic consists of aligning and stitching several images together with overlapping fields of view (mosaicking), in order to create a single composite high-resolution image (Escartín et al. 2008;Rende et al. 2015). Image stitching algorithms built on feature-based techniques, consists of three main components: image feature extraction (i.e. the process to establish correspondences between points, lines, edges, or other geometric entities to detect corresponding feature points in an image pair), registration (i.e. the process of aligning the images captured from different viewpoints), and blending (i.e. the process applied across the stitch to make both the seams less apparent and to smooth the transition between images). During image registration, image pairs are compared to find the translations that can be used for the alignments of images. After image registration, images are blended to form a single image (Adel et al. 2014). Many feature detectors have been proposed over the years such as SIFT (Lowe 2004), SURF (Bay et al. 2008), FAST (Rosten and Drummond 2006), and ORB (Kulkarni et al. 2013). The well-known SIFT (Scale Invariant Feature Transform) is very robust and exhibits invariance to image translation, rotation, scaling transformations, and good robustness to light changes, noise, and affine transformation (Xing and Huang 2010;Zhao et al. 2019). Such features make SIFT a widely used extraction method for feature points from the past 15 years.
Underwater photomosaics derived from open source panorama software such as MicrosoftÓ Image Composite Editor (ICE) coupled with Side Scan Sonar acoustic data have proved to be a non-destructive and cost-effective method for ground-truthing in marine habitats mapping (Pergent et al. 2017). Moreover, a panoramic view of the seabed can be also applied to quantify seagrass distribution and abundance improving micro-cartographic analysis (Rende et al. 2015) as well as to accurately map tropical reefs (Judilla et al. 2012) and archeological sites (Vasilijevic et al. 2015).
In this paper, a cost-effective sampling technique based on photomosaics obtained through the image stitching approach for studying coralligenous reefs is introduced and evaluated. Although photogrammetric surveys were carried out on coralligenous reefs to study Axinella cannabina (Esper, 1794), Paramuricea clavata (Risso, 1826), and Corallium rubrum (Linnaeus, 1758) facies (Feral et al. 2014;Palma et al. 2018;Royer et al. 2018), neither tests nor applications of photomosaic technique are currently reported.
Data collected by photomosaics and the most common photoquadrats technique were employed to assess the composition of coralligenous assemblages, and the results were compared to check for potential discrepancies between the two methods. Potential effects on mosaics related to the use of different commercial cameras, representative of different market segments, were also assessed. Finally, the effectiveness of photomosaics, practical and theoretical advantages, opportunities, limitations, and challenges for future researches are discussed.

Study area
The study was carried out along the east coast of Giglio Island (central Tyrrhenian Sea, Italy), inside the restricted access area defined after the 2012 Costa Concordia shipwreck and during its salvage operations (Fig. 1). Any touristic or recreational activity (i.e., diving, sailing, and swimming) was forbidden in the wreckage area to avoid any disturbance during the wreck removal phases and to maximize safety in the case of an accident. A brief description of the morphological features and the evolution of the study area is given by Avio et al. (2017) All the pictures used in this study were acquired along the northern vertical wall of the Gabbianara ridge, between 19 and 25 m depth, where thick coralligenous reefs are present (50 m long 9 15 m height approximately). The bioconstruction was impacted by the presence of the wreck and its salvaging operations that caused fine sediment dispersion and the partial loss of benthic organisms (Casoli et al. 2017(Casoli et al. , 2020Penna et al. 2018). We decided to test here this new sampling technique in virtue of applying in future to assess the recolonization and the resilience processes on the aforementioned coralligenous reefs.

Data collection
Several objects of known dimensions were arranged on the reef as references to scale the scene consisting of four ground control points (GPCs, 20 9 18 cm) and one 20 cm ruler.
Samplings were carried out by scientific scuba divers during late spring of 2019, according to STAR procedure (STAndaRdized coralligenous evaluation procedure; Piazzi et al. 2019b), a standardized protocol for coralligneous reefs monitoring. The divers swam at a low speed of 10-15 m/min, describing parallel and regularly spaced transect to allow the images overlay, following a boustrophodonic pattern (''mow the lawn'' method; Henderson et al. 2013;Burns et al. 2015). Pictures were acquired with a time-lapse of 1 s at a regular distance of about 50 cm from the vertical wall, with pure nadiral orientation. Acquisition parameters have been chosen following the best combination described by Marre et al. (2019). The portion of the reefs interested by the investigation covered a surface of 10.21 m 2 .
Image acquisition for photomosaics construction was carried out using three different cameras which differed in sensors size/pixel resolution and cost: Canon EOS 5D, Sony a7-III, and GoPro Hero 5 (Table 1). Canon EOS 5D and Sony a7-III settings during image acquisition were kept the same: f/8, shutter speed 1/100 s., and ISO 250. GoPro Hero 5 was used in automatic mode. For all three cameras, the white balance was set in automatic mode and eventually corrected later during image post processing.
All images were stitched into a single, high-resolution image mosaic using Agisoft PhotoScan Professional Edition V. 1.4.0 (Agisoft 2012). It performs 3D reconstruction of objects from images and employs the whole photogrammetric workflow. At the first stage, PhotoScan detects points in the source photos which are stable under viewpoint and lighting variations and generates a descriptor for each point based on its local neighborhood. These descriptors are used later to detect correspondences across the photos. The tie point extraction and detection employ an improved version of SIFT (Chiabrando et al. 2015;   Murtiyoso et al. 2018). Furthermore, the dense matching method used may be a variant of the Semi-Global Matching (SGM) algorithm. We chose this software because of a user-friendly interface with several simplified functionalities and parameters, which is advantageous for many nonexpert users. It should be noted that we did not focus our work directly on 3D output [e.g., dense point clouds, triangular mesh Models, and Digital Terrain Models (DTMs)], but we used the 2D orthomosaics for digital camera comparison and analysis of the assemblage. Ten images of a calibration checkboard were acquired underwater for each camera with the same target to lens distance used in the following data acquisition. Subsequently, these images were used for the calibration of optics in Agisoft Lens, a tool implemented by Agisoft package which allowed the estimation of camera calibration parameters, including the non-linear distortion coefficients. These parameters were then inputted into Photoscan as.xml file to speed up pre-calibration steps aimed at removing lens distortion.
To assess the accuracy of the photomosaic and compare it with the commonly used method of photoquadrats (Parravicini et al. 2009;Kipson et al. 2011;Deter et al. 2012;Trygonis and Sini 2012;Piazzi et al. 2017Piazzi et al. , 2019b, five photographic samples of 0.2 m 2 (50 9 40 cm, for a total of 1 m 2 ) were randomly acquired in the same portion of the reef by Canon EOS 5D (Supplementary materials, Figure S1): The camera was arranged within a PVC frame to minimize parallax error. The distance from the substrate, as well as lighting configuration, was the same used for photomosaic image acquisition. The comparison among methods was made with Canon EOS 5D exclusively, to avoid any difference due to camera and light apparatus.
Photomosaics in.tiff format were imported in ESRI ArcMap 10.2.2, where georeferencing and scaling routine were carried out as follows: two points placed 10 cm apart, acquired by real-time kinematic (RTK) GPS out of the water in real coordinates (WGS84, UTM fuse 32 N) were associated with two points in the photomosaic with the same linear distance (10 cm), which was measured in the picture by the 20 cm ruler placed in the scene. This allowed a reference for the scaling routine rather than a correct georeferentiation of the photomosaics: The pixel to distance (pixel/cm) conversion was carried out through the Georeferencing tools in ArcMap. Besides, five photoquadrats (50 9 40 cm) were successively over imposed on the Canon photomosaic selecting conspicuous points where control points were added through the specific tool for georeferencing raster images available in ArcGIS. To avoid excessive distortion of raster photoquadrats, we selected a first-order polynomial transformation to transform the raster dataset to the map coordinates. The total error is computed by taking the root mean square (RMS) sum of all the residuals to compute the RMS error. This value describes how consistent the transformation is between the different control points, so we tried to keep this value as close as possible to zero, adding no less than 3 control points. The five photoquadrats used to test the accuracy of the photomosaic were haphazardly distributed on the studied reefs and covered 1 m 2 of surface.

Polygon data digitization
Sessile organisms were identified to the lowest taxonomic level possible allowed by the images. As for conspicuous organisms, the species level was achieved, whereas similar taxa that could not be identified through photographs were aggregated into taxonomical/morphological categories. To assess the composition of coralligenous assemblage in both photomosaics and photoquadrats, polygons were manually outlined by the freehand drawing tool in ArcGIS10.2.2, and then assigned by the user to a specific taxon or morphological category. This method is in use in different image analysis software (Trygonis and Sini 2012). The percentage of substrate surface covered by a given taxa or category (percentage cover, PC) was automatically calculated in ESRI ArcMap 10.2.2 by using Calculate area tool, included in the Spatial statistic toolbox.

Accuracy assessment of photomosaic
Once the polygon data digitization was complete, the main taxa/morphological categories were identified on both photomosaic and ground-truth data (i.e., photoquadrats) acquired by Canon EOS 5-D. We assessed the accuracy of the taxonomic classification by comparing percentage cover derived from the 5 quadrats cut out of the photomosaic image with manually classified photoquadrats. A confusion matrix was calculated to evaluate the accuracy of the classification made on photomosaic image including: the overall accuracy, which is defined as the proportion of all correctly classified objects and the total sample size; and the Kappa index of agreement (KIA), which is defined as the agreement of the classification results with the corresponding reference data. The proposed categories for the assessment of the classification performance measured by the Kappa value were the following: poor (B 0.1), slight (0.1-0.20), fair (0.21-0.40), moderate (0.41-0.60), substantial (0.61-0.80), and excellent (C 0.81) (Sim and Wright 2005). Differences in the structure of assemblages were assessed through a one-way multivariate permutational analysis of variance (PERMANOVA). PC data were not transformed in order to stress the importance of the abundance of taxa/morphological categories in determining the differences among Photomosaic and Photoquadrats (Methods, as fixed factors).

Cameras comparison
The number of taxa/morphological categories (N) and PC were, respectively, counted and measured in each photomosaic to compare the three cameras used in the study. The assessments were carried out separately in each photomosaic by a single observer, to avoid any analyst effect. Furthermore, the number of patches (NP) was counted in each photomosaic, where a patch was defined as a surface covered by specific taxa/morphological category well distinguished from the surrounding organisms. Taxa/morphological categories with no distinctive and recognizable borders (i.e., largely diffused crustose algae and sediment) were excluded from this analysis. The three parameters were assessed by ESRI ArcMap 10.2.2. Then, the three photomosaics were divided into 9 sub-samples, and the differences in terms of NP among cameras were tested through a one-way multivariate PERMANOVA. Furthermore, pairwise tests were used to discriminate among cameras. A similarity percentage-species contribution analysis (SIMPER) was carried out to identify the taxa/morphological categories that mostly contributed to the dissimilarity among cameras.

Results
Analyses were carried out on a portion of coralligenous reef measuring 10.21 m 2 surface area. The number of pictures acquired by each camera and used to build the three photomosaics ranged from 110 to 120. The three mosaics had a very high resolution, 0.36 mm/pixel, 0.21 mm/pixel, and 0.50 mm/pixel, for Canon EOS 5D, Sony a7-III, and GoPro Hero 5, respectively (Supplementary materials, Figure S2). The features of the cameras utilized during the study are summarized in Table 1.

Method comparison
Accuracy assessment carried out on the mosaic derived from Canon EOS-5D imagery showed high fidelity in class identification when compared to the five ground-truth photoquadrats with average overall accuracy and KIA value of 76.49% (± 10.76 SD) and 0.71 (± 0.15 SD), respectively. Major errors were due both to the largest (e.g., Peyssonnelia spp.) and most variable (Pseudochlorodesmis furcellata, mucilage, Reteporella grimaldii) taxa/categories which often obscured smaller categories, such as serpulids or erect ascidians (Halocinthia papillosa), leading to category-cover misclassification. No significant differences in the composition of the assemblage were found among photomosaic and photoquadrats methods ( Table 2). The analysis of the five replicates acquired by photomosaic and photoquadrats methods allowed the identification of 27 taxa/morphological categories of the total amount of 33 taxa/morphological categories identified on the whole study area: percentage cover (PC) measured in the five photomosaics and photoquadrats replicates showed an overall similar pattern (Fig. 2).

Assemblage assessment through photomosaics
Thirty-three taxa/morphological categories were identified and mapped by Canon photomosaic, followed by Sony (N = 31) and GoPro (N = 27) photomosaics.
The assemblages assessed through the three photomosaics were overall similar (Fig. 3). Twelve taxa accounted for the 97.2%, 97.7%, and 99.0% of percentage cover (PC), in Canon, Sony, and GoPro, respectively, whereas twentyone taxa showed PC less than 0.6% when present (supplementary materials, Tab. S1). Algae, bare surfaces, and sediment characterized the coralligenous reef of the studied site. The most common algal taxa were: Peyssonnelia spp., P. furcellata, mucilage, Palmophyllum crassum, and Flabellia petiolata. As for benthic invertebrates, exclusively R. grimaldii and Eunicella cavolini showed PC higher than 1%. The other taxa were quite rare or scarcely reported in the study area. In particular, small erect bryozoans, Ceramiales (erect red algae), articulated Corallinales, and the bryozoan Margaretta cereoides were not reported in GoPro photomosaic, whereas solitary corals and boring sponges were absent in both Sony and GoPro photomosaics.
The number of patches (NP) significantly differed among cameras (Table 3, Fig. 4); pairwise test highlighted that these differences were between Canon and GoPro cameras. Serpulids, encrusting bryozoans, F. petiolata, R. grimaldii, P. crassum, mucilage, and encrusting sponges were the taxa/morphological categories responsible for the dissimilarity between cameras (SIMPER analyses, Supplementary Materials Table SXX). Considering that PC values described similar assemblages, small patches were mostly responsible for the differences here described.

Discussion
Image mosaicking represents an multi-purpose and lowcost technique to survey and sample coralligenous reefs. It offers new perspectives to study the diversity, ecological status, dynamic processes, growth, and competition phenomena occurring in one of the most important Mediterranean bioconstructions. The combination of large-scale images (tens of m 2 ) and high-resolution benthic information (sub-centimeter scale) is the main advantage of photomosaic technique. This perfectly copes with coralligenous reefs features, such as small scale spatial heterogeneity due to the patchy distribution of species (Ferdeghini et al. 2000;Piazzi et al. 2004), high diversity, minimum sampling effort (Kipson et al. 2011), presence of both large colonial and small organisms (Casoli et al. 2016b;Sini et al. 2019). Larger size of sampling unit increases advantages in using digital photography as sampling method by reducing bias errors and producing similar results to those obtained by visual assessment (VA) (Parravicini et al. 2009). Furthermore, permanent high-quality records of the seafloor allow the application/replication of measures through different data extraction techniques, according to researchers' needs (i.e., accuracy, resolution, time). In the present study, we carried out percentage cover estimation by image analysis (manual drawn), but also other data extraction techniques could be performed, such as VA, frequency of occurrence, and point-intercept transects (Van Rein et al. 2012). The performance of photomosaic and photoquadrat techniques was not significantly different. The high precision, highlighted by overall accuracy and KIA, reveals that photomosaic is consistent with the photographic-based methods commonly used to study coralligenous reefs. Therefore, sub-sampling the large areas covered by photomosaics can be useful for the application of the most common indexes assessing the ecological quality of coralligenous habitats (Deter et al. 2012;Gatti et al. 2015;Ferrigno et al. 2017;Montefalcone et al. 2017;Piazzi et al. 2017;Sartoretto et al. 2017). Another feature of this method is the possibility to identify changes at the landscape level, repeating observations on the same portion of the seabed. Avoiding the deploying of several permanent markers or plots as for the photoquadrats method (Garrabou et al. 2002;García-Gómez et al. 2020) may allow minimizing the disturbance in community dynamic studies. The advantages and disadvantages of the two methods were summarized in Table 4. Both the methods represent cost-effective sampling strategies, with costs varying according to the camera model and light apparatus used. The comparison among three cameras showed high precision in the description of the coralligenous assemblage, thus supporting the suitability of low-cost cameras, as well. Nevertheless, differences were reported for low abundant or scattered taxa/categories, exclusively. The GoPro Hero 5 showed lower accuracy than Canon EOS 5D for both composition and abundance of taxa: camera features and light apparatus were the main drivers of the reported differences. Blurry pictures could have influenced the identification of smaller organisms. On the other hand, differences among Sony a7-III and Canon EOS 5D may be attributable to light, even though results were not consistent with the camera resolution (24.2 Mpixel and 12.8 Mpixel, respectively): strobes provide more power and higher angle of coverage. It produced an increase in PC of holes in both GoPro and Sony (8.92% and 8.16%, respectively) if compared to Canon (7.63%), and the resulting classification of darkening areas as cavities instead of benthic organisms. The use of an appropriate illumination system is pivotal to avoid loss of information of smaller or cryptic species. The slight differences among camera models can also be explained in terms of camera orientation, which is one of the most important parameters affecting reprojection errors. According to Marre et al. (2019), a pure nadiral orientation was used during image acquisition given the best results in their methodological study. In the case of more complex environments such as coralligenous reefs, details in holes, crevices, and overhang can be lost without acquiring oblique images. Nevertheless, the camera orientation was set according to the aim of the study that was to create a two-dimensional photomosaic of the coralligenous wall and not a three-dimensional photogrammetric model of horizontal surfaces using SfM routine (Enochs et al. 2015; Palma et al. 2018Palma et al. , 2019Royer et al. 2018;Marre et al. 2019Marre et al. , 2020aVentura et al. 2021). One of the main limitations of the proposed method is the impossibility of absolutely georeferencing of points on the vertical seafloor by using classical ground control points (GCPs), in which real coordinates are usually estimated for the three axes (x, y, and z) with a high degree of accuracy (in the order of mm with professional-grade GNSS systems during aerial campaigns). The use of marker buoys, which would have allowed the position of the sampling area in the real space and the collection of the positions of reference points from the surface, cannot be achieved due to currents and depth which often characterize infralittoral sciaphilous assemblages on coralligenous outcrops. Despite the high cost, the georeferencing of points on the vertical seafloor can be carried out through forward-looking Multi Beam Echo Sounder (MBES) from ROVs which would allow a traditional photogrammetric approach: such method was efficiently applied in deep waters to achieve terrain descriptors and map species distribution on full vertical environments (Huvenne et al. 2011;Robert et al. 2017). The multi-purpose features of the photo mosaicking technique and the possibility to observe the same surface over time are particularly well-adapted to assess the recolonization and dynamic process occurring on coralligenous reefs characterized by bare surfaces. The operation aimed at stabilizing the bottom where the wreck of Costa Concordia caused fine sediment release in both 2013 and 2014 (Casoli et al. 2017), which partially buried the basal portion of the studied bioconstruction. From January 2015 to April 2018, the seafloor of the restricted area was cleaned by the disturbances derived from the wreck removal (Remediation phase): A total of 25,000 t of sediment were removed. As a consequence of the sediment removal, dead portions of coralligenous reefs were revealed (Fig. 5), likewise the sediment-exposed terrace described after exceptional stormy events (Teixidó et al. 2013;Betti et al. 2020). Such bare surfaces represent unique substrates to study the natural recovery of coralligenous reefs: The application of photomosaic technique will help to comprehend the community succession patterns.
Photomosaic represents a safe technique up to 40 m depth: On average, 8 ± 3 min were necessary to cover areas of approximately 10 m 2 characterized by high resolution, as described in this study. Such average was timed in a total of 9 preparatory data acquisition, included the three photomosaics used in the present work. Photomosaic surveys can be efficiently applied to cover large surfaces up to 40 m depth while remaining within the safety curve ( Fig. 6; Casoli, unpublished data). The recent progress in technical dive, as well as the use of Close Circuit Underwater Breathing Apparatus (CCUBA; i.e., rebreather), if coupled with photomosaic might contribute to further increase the size of the sampling area. In light of increasing divers' safety, the application of AUVs or ROVs technologies might represent a step forward, including georeferenced pictures, no depth-related constraints, and considerably increasing the investigated area. Besides species distribution in deep or inaccessible environments, they allow accurate reconstruction of seafloor features on vertical structures (Robert et al. 2017). On the contrary, limitations in the use of AUVs or ROVs systems for shallow coastal habitats studies are represented by their High-resolution benthic information 4 4 Large-scale images 4 7 Application of different data extraction techniques 4 7 From species to landscape investigations 4 7 Georeferentiation 4(horizontal seafloors) 7 Application of ecological indexes 4 4 Image analysis effort 7 4 Fig. 5 The dead part of concretion at the bottom of coralligenous reefs revealed after the sediment removal. The boundary between living bioconstruction and bare surfaces indicated the burial level challenging costs and maneuverability issues; seabed complexity and presence of vertical walls might impact the outcome other than the recovery of the AUVs. The time required for image processing and alignment was about 1 h for each photomosaic, according to time reported by Palma et al. (2018) and Marre et al. (2019). On the contrary, the most time-consuming step was image analysis which required [ 6 h. It depended on the accuracy and resolution of the process itself, but it could drastically vary depending upon data extraction techniques. The implementation of an automated benthic species identification for coralligenous reefs will drastically reduce such time-consuming analyses, allowing to minimize limits related to this tool (Chirayath and Li 2019;Williams et al. 2019;Marre et al. 2020b).
Technology advances support the researchers to easily and efficiently recognize ecological processes and changes in the marine environment. This study outlined an accurate, reproducible, and cost-effective method to study and monitor coralligenous communities at different depths reached by SCUBA divers. The multi-purpose features of photomosaic allow different investigations, from species to seascape level. In particular, this method could be applied to assess the effects of both climate anomalies and human disturbances effects on protected habitats, such as coralligenous reefs, through the definition of fixed monitoring stations. The photomosaic method improves the opportunity of collection of field data, not excluding standard sampling procedures or the application of indexes aimed to assess the ecological quality of the investigated communities.
Conflict interest All the authors, as scientists of the Sapienza University of Rome and CIBM, were commissioned by the Titan-Micoperi group and Costa Crociere S.p.A. to assess the environmental risks of engineering operations and to monitor the environmental conditions in the impacted areas in order to protect marine habitats during the salvage of the Concordia wreck. The authors' work was controlled and continuously shared with an Italian governmental body (The 'Osservatorio'), which fulfilled the public function of verifying the performance of the monitoring.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.