Hyperspectral and Lidar: Complementary Tools to Identify Benthic Features and Assess the Ecological Status of Sabellaria alveolata Reefs

Sabellaria alveolata is a sedentary gregarious tube-building species widely distributed from southwest Scotland to Morocco. This species builds what are currently considered the largest European biogenic reefs in the bay of Mont-Saint-Michel (France). As an ecosystem engineer, S. alveolata generates small to large scale topographic complexity, creating numerous spatial and trophic niches for other species to colonize. Sabellaria reefs are also under anthropogenic pressures, leading locally to massive degradation. However, stakeholders lack spatially explicit measures of reef ecological status, at adapted spatial resolution to provide key management information for this protected habitat. Traditional field surveys are extremely time-consuming and rely on expertise for visual ecological status assessment. The present study aims at using an automatic processing approach based on optical airborne data to (i) assess the potential of hyperspectral imagery to discriminate Sabellaria bioconstructions and its main ecosystem associated habitats, including different types of substrate as well as biological components and (ii) to use the combination of the hyperspectral and LiDAR signals to estimate the spatial structure of the different bioconstruction types (veneers vs. hummocks and platforms) and ecological phases (retrograding and prograding). A reef from Mont-Saint-Michel was used as a test site. We built a processing chain based on supervised classification using the Mahalanobis distance to generate an accurate distribution map (overall accuracy of 88% and a Kappa of 0.85) of 10 Sabellaria-related benthic features, including large reef developing on sand and smaller veneers encrusting rocky shore areas. Specific spectral indices were used to define the spatial distribution of the main primary producers, in particular the microphytobenthos. Joining the hyperspectral and LiDAR data led characterizing the distribution of S. alvealata’s ecological status (prograding and retrograding phases) with an overall classification accuracy and Kappa coefficient that can respectively amount to up to 93 and 0.86. In our study site, the Sabellaria reef area (between 5.52 and 6.76 ha) was dominated by retrograding phases (between 53 and 58%). Our results showed that this automatic processing chain could be relevant for the spatial characterization of other Sabellaria reef sites. Study perspectives tend toward a quantitative estimation of their ecological status index.

Sabellaria alveolata is a sedentary gregarious tube-building species widely distributed from southwest Scotland to Morocco. This species builds what are currently considered the largest European biogenic reefs in the bay of Mont-Saint-Michel (France). As an ecosystem engineer, S. alveolata generates small to large scale topographic complexity, creating numerous spatial and trophic niches for other species to colonize. Sabellaria reefs are also under anthropogenic pressures, leading locally to massive degradation. However, stakeholders lack spatially explicit measures of reef ecological status, at adapted spatial resolution to provide key management information for this protected habitat. Traditional field surveys are extremely time-consuming and rely on expertise for visual ecological status assessment. The present study aims at using an automatic processing approach based on optical airborne data to (i) assess the potential of hyperspectral imagery to discriminate Sabellaria bioconstructions and its main ecosystem associated habitats, including different types of substrate as well as biological components and (ii) to use the combination of the hyperspectral and LiDAR signals to estimate the spatial structure of the different bioconstruction types (veneers vs. hummocks and platforms) and ecological phases (retrograding and prograding). A reef from Mont-Saint-Michel was used as a test site. We built a processing chain based on supervised classification using the Mahalanobis distance to generate an accurate distribution map (overall accuracy of 88% and a Kappa of 0.85) of 10 Sabellariarelated benthic features, including large reef developing on sand and smaller veneers encrusting rocky shore areas. Specific spectral indices were used to define the spatial distribution of the main primary producers, in particular the microphytobenthos. Joining the hyperspectral and LiDAR data led characterizing the distribution of S. alvealata's ecological status (prograding and retrograding phases) with an overall classification accuracy and Kappa coefficient that can respectively amount to up to 93 and 0.86.

INTRODUCTION
Biogenic reefs are not only limited to coral reefs but also include a wide variety of species capable of building biogenic structures (or bioconstructions), including the gregarious tube-building worms of the Sabellariidae family. Over the world coastlines, large bioconstructions are found on the coasts of South America, built by the species Phragmatopoma lapidosa (Kinberg, 1866) and Sabellaria wilsoni (Lana and Gruet, 1989). On European coasts, the largest intertidal bioconstructions are built by the species Sabellaria alveolata (Linnaeus, 1767), which is presented from the north of the English coast to the south of the Moroccan coast. Sabellaria growths mainly on rocky area (Muir et al., 2016), but can also be present on soft sediments. In France, the Sainte-Anne reef (Mont-Saint-Michel bay) is known as the largest European bioconstruction. S. alveolata is an engineer species that generates small to large scale topographic complexity and creates numerous spatial and trophic niches for other species to colonize (Dubois et al., 2002(Dubois et al., , 2006Bonifazi et al., 2019). Moreover, by modifying local hydrodynamics, S. alveolata reef structures also impact the habitats of benthic assemblages located in the immediately surrounding sediments (Dubois et al., 2002;Jones et al., 2018).
S. alveolata reefs are also under direct mechanical disturbances (recreational fishing) and indirect anthropogenic pressures (shellfish cultures) that can lead to massive degradation of some reef areas. Current damages are characterized by irregular surface features with lower densities of S. alveolata tubes (Dubois et al., 2002). Gruet (1982) described an evolving S. alveolata morphotype, from "veneers, " where the tubes adhere to the substratum, to "hummocks" where the tubes radiate out from the initial settlement point before reaching "platforms" formed of extensive areas of hummocks fused together. Any one of these types can display outward signs of being in a "progradation" or "retrogradation" phases (Curd et al., 2019), either partially or totally. At European level, these formations are targeted, under the designation of Special Areas of Conservation (SACs), as a part of marine habitat (type 1170 "Reef ") to be protected by the European Union's Habitats Directive 92/43/EEC. The latter states that bioconstructions should currently be monitored in order to acquire the missing knowledge regarding their management and ultimately their preservation.
Several studies searched for indicators to characterize the ecological state of S. alveolata reefs. Desroy et al. (2011) developed a "Health Status Index" for the Mont-Saint-Michel bay bioconstructions, based on fragmentation, morphological types and epibiont coverage. However, this index is only applicable to areas where S. alveolata formations reach extensive bioconstructions, which strongly restricts its application potential. More recently, Curd et al. (2019) attempted to investigate the relationship between the structural appearance of the bioconstruction and physiological state of the tube-building polychete. In this study, physiological indicators (lipid markers and metabolic enzymes) were used as biochemical indicators of S. alveolata health status, and linked to bioconstructions types (veneers or hummocks). Biochemical analyses are useful tools to investigate worm's health status, but they are costly, time-consuming and not always connected to the reef status or phase (progradation vs. retrogradation). Other indicators were previously developed from in situ surveys. However, their use is limited to the surrounding extension of the evaluated area: analysis carried out on a larger scale (whole reef) leads to a simplification of the "landscape" and therefore potentially to wrong reef ecological assessments. Jones et al. (2018) analyzed the effects of a continuous and increasing disturbance on beta diversity indices. They also stressed the importance of spatial and temporal scales for evaluating the impact of disturbances on S. alveoata reefs. Indeed, there is a real need for effective management to develop alternative analytical tools that (i) allow for approaching metrics measured in situ to assess conservation status on a large scale and (ii) guarantee robust reproducibility of acquired information. Moreover, all these studies used destructive sampling techniques, which although still commonplace in benthic ecology, are not viable options in low-abundance habitats (Beisiegel et al., 2017).
Several studies of S. alveolata reefs have been achieved using remote sensing techniques. Acoustic tools were more specifically adopted as suitable method for Sabellaria reefs in subtidal areas (Griffin et al., 2020). The most the studies were mainly based on aerial photography (Bonnot-Courtois et al., 2005;Desroy et al., 2011) or satellite imagery (Marchand and Cazoulat, 2003;Jones, 2017) as well as UAVs (unmanned aerial vehicles) as an efficient and cost-effective tool for fine scale reef characterization (Ventura et al., 2018;Collin et al., 2019). Noernberg et al. (2010) have also exploited aerial photographs and LiDAR (Light Detection and RAnging) digital elevation model to quantify intertidal complex landform volume. A more recent work used artificial neural network modeling to map S. alveolata relative abundance (Collin et al., 2018). Such studies have been successful to investigate reef formations developing on soft-sediment, where bioconstructions sit on top of sediment are discrete objects relatively easy to identify. Yet, most of the reef formations in Europe develop over hard rocky substrata, and compete for space with other species and macroalgae (Firth et al., 2015). To the best of our knowledge, no study has tackled the problem of automatically separating reef bioconstructions from other associated benthic features in a rocky reef environment, nor has any study tried to remotely and automatically assess S. alveolata reef ecological status.
To overcome limitations in large scale S. alveolata reef characterization, we propose an approach based on the combination of two complementary techniques: (i) LiDAR a widely used tool to depict the 3D habitat structural complexity (Wedding et al., 2008;Davies and Asner, 2014;D'Urban et al., 2020) and (ii) hyperspectral imagery that allows for the characterization of bio-optical properties variability in time and space (Fyfe, 2003;Schmidt and Skidmore, 2003;Chennu et al., 2013;Bajjouk et al., 2019b). Hyperspectral imagery has been used in coastal areas to study shallow tropical coral reefs, reef-forming species with symbiotic microalgae species with primary photosynthetic pigments (Petit et al., 2017;Bajjouk et al., 2019a). Characterizing reef-forming species such as S. alveolata, a species without photosynthetic pigments and obvious spectral reflectance is a challenging task. By using hyperspectral images coupled with field observations, we also aim at better estimating the microphytobenthos (hereafter MPB) spatial distribution in sediments surrounding the bioconstructions. MPB is a key information for understanding the role of engineered habitats and species biological activities (feces and pseudofaeces) (Echappé et al., 2018).
We address here three main questions: (1) Do the bioconstructions of S. alveolata have a specific spectral signature, so that one can discriminate reefs from other associated habitats, including MPB patches (2) Does the spectral signature of S. alveolata bioconstructions allow for the identification of reef constructions affixed to intertidal rocky substrate as well as soft-sediments? (3) Can we automatically and remotely estimate the ecological health status of these bioconstructions based on hyperspectral and LiDAR signal combination? Reefs located in the bay of Mont-Saint-Michel were used as a test site. Resulting spatial patterns are discussed regarding their ecological coherence. Future developments are finely formulated for S. alveolata conservation and management perspectives.

STUDY SITE DESCRIPTION
The Bay of Mont-Saint-Michel (France) is a coastal area of the Western part of the English Channel. The study site, named Champeaux (−1.53 • , 48.73 • ), was located in the east part of this bay that hosts the largest S. alveolata reefs in Europe (Figure 1). The sedimentary environment of the bay is mainly controlled by tidal residual currents, typically characterized by an anticyclonic gyre (Noernberg et al., 2010). This site was chosen for its bioconstruction density and diversity of habitats: it includes both reefs developing on sand and bioconstructions developing on rocky shores (Lecornu et al., 2016) and stretches from the upper intertidal to the lower eulittoral zone.
In Champeaux, S. alveolata can build three-dimensional structures which present a variety of morphotypes (Gruet, 1972;Curd et al., 2019): elevated platforms, more or less coalescent ballshaped structures (hummocks) or adhering to rocks (veneers). Each morphological structure can display different phases in

Field and Laboratory Measurements
Several field surveys were carried out in 2019 across Champeaux site during the low spring tides: (i) A field spectral library was acquired on the 22nd of April.
Reflectance spectra were measured on 505 channels from 325 to 1075 nm using an Analytical Spectral Devices (ASD) field spectroradiometer that automatically averages 5 multiple records to provide a single output. The ASD was calibrated via a spectralon to provide reflectance standardized measurements. Reflectance recordings were carried out vertically, avoiding any shadow on the target and at a height allowing a homogeneous footprint. For each target, about three to five spectra were recorded to take into account intra-target variability. These measurements were taken on different homogeneous samples representing different types of benthic features: the spectral library is composed of 53 reflectance spectra of 10 different pure targets including three vegetation types (green and brown algae, MPB), S. alveolata bioconstructions developing on both rock and sediment area at different ecological phases, oyster reefs as well as substrate types (mud, sand, shells and rocks). Principal component analysis (PCA) was used as an exploratory tool to summarize and to visualize the information in the benthic features described by their spectral signatures. PCA was performed using R statistical software. (ii) A total of 221 field control points ( Figure 2B) were set on the 17th and 18th of June. Photography is associated with each spectroradiometric measurement. All the descriptive information of the target (substrates and epibionts) was also recorded for each measurement. Efforts were made to cover a wide range of ecologically diverse features of the reef providing an on-foot access. (iii) During the same summer field survey, on June 17th, samples were collected for identification of the main genera or species present in MPB communities. MPB was harvested using pieces of artificial substrates, cut in a nylon net of 100 µm mesh. Substrates consisted in small squares of 2.5 cm 2 . They were deployed in four stations located few meters apart from each other, in an area close to S. alveolata reefs and where MPB biofilms are recurrently observed at low tide. Pieces of artificial substrates were placed on wet FIGURE 2 | Data acquired on Champeaux site. Hyperspectral (A) and field data (B) were specifically acquired for the present study. lidar bathymetry data (C) were provided by SHOM, the French hydrographic office and rugosity (D) was computed from bathymetry data.
Frontiers in Marine Science | www.frontiersin.org sand 2 h before low tide. They were collected at low tide, when colonized by a biofilm of MPB cells, and immediately transferred to plastic tubes containing 4 ml of filtered seawater with a fixative (1% final concentration of Lugol's solution). Samples were brought back to the laboratory and kept at 4 • C until analysis. Before microalgal cell identification and counting, tubes were agitated in order to help for detachment of MPB cells from the substrate. Then, a sub-sample was taken, diluated 100x in filtered seawater and settled for 12 h in a Utermhöl chamber (Utermöhl, 1958). Most of supernatant was replaced by bleach (0.1% diluted in milliQ water) during one night, in order to remove a large part of organic matter residues and preserve silice frustule. Bleach was then replaced by milliQ water and microalgal cells were settled again for 12 h. Identification and counting of cells were conducted under 400 magnification, using an inverted microscope (Zeiss Axio Observer).
All field position measurements were made using a GARMIN GPS with an uncertainty of 3 m. The position recording of ground control points was done by allowing the device to capture as many satellite signals as possible and thus increase the precision of the geolocalization. Positions have then been expressed in the same projection system as that of hyperspectral images (Lambert-93 associated with the RGF93 geodetic system). Locations of ground control points were matched to images and verified by visual inspection. This last operation was performed by experts who know well the surveyed area using the ArcGIS tool.

Airborne Optical Data Acquisition
An airborne survey was carried out on June 17, 2019 as part of the BIOHERM Project, using an hyperspectral imaging sensor of the Hyspex VNIR 1600 type deployed at an altitude of 1 200 m (Figure 2A). The sensor measured radiance between 400 nm and 1 000 nm with spectral and spatial resolutions of respectively 3.7 nm and 0.5 m, providing a total of 160 spectral bands.
LiDAR data ( Figure 2C) were acquired in March 2019 as part of the Litto3D program (Litto3D, 2019). LiDAR topobathymetric data comes from a HawkEye III sensor used in dual configuration: Chiroptera (Topo-Shallow Channel) and Deep Channel. The point cloud density in topographic configuration (corresponding to the area of interest for our study) is 8 points/m 2 with a planimetric and vertical precision better than 20 cm. All data are expressed in the Lambert-93 projection system associated with the RGF93 geodetic system, and in the IGN69 altimeter system, compatible with the WGS84 world system at metric level. For the purpose of the analysis, the LiDAR point cloud was interpolated to a 1 m regular grid DEM (Digital Elevation Model). The DEM was also used to create layers describing seafloor topographic rugosity ( Figure 2D) using Arc GIS Benthic Terrain Modeler tool (NOAA and Oregon State University).

Data Pre-processing
Geometric and atmospheric corrections were performed prior to image processing (Figure 3). Orthorectification allowed generating geolocated images in the Lambert 93 system. The relative geometric precision obtained was of the order of the pixel. Atmospheric corrections were subsequently applied using the ATCOR software on each flight line in order to obtain reflectance images above the water surface. The oxygen absorption band, around 760 nm, was linearly interpolated. A post empirical correction was thereafter applied. It was carried out in two steps: (i) a linear correction based on the overlap between the lines to correct the variation in lighting conditions, (ii) a second linear correction, for each spectral band, which allowed for the correction of local atmospheric effects not or poorly taken into account by the atmospheric model. Their correction was done using measurements on a calibrated reference target installed on site. Land and water area were masked using, respectively, coastal shoreline derived from LiDAR data and Infra-Red band thresholding.

Data Processing
The imagery data processing procedure consisted of several steps. They are illustrated in the workflow diagram (Figure 3) and are detailed below.

Vegetation Index Computation
At first, spectral vegetation indices were computed in order to specifically characterize vegetation present in the study site. The principle consists in linking together certain characteristics and the measurements acquired in at least two spectral bands. Two indices were calculated using ENVI image processing software: (i) the Normalized Difference Vegetation Index (hereafter NDVI, Tucker, 1979), widely used to analyze vegetation (Cracknell, 2001;Pu et al., 2008;Kushnir et al., 2011) and (ii) the Diatom index ( Van der Wal et al., 2010;Launeau et al., 2018) as MPB sampling analysis revealed that this compartment was mainly dominated by diatoms in the study site ( Figure 4N).
To enhance discrimination of benthic features that characterize the area, topographic and spectral indices layers were merged with hyperspectral masked image, as classification input, using layer stacking processing (Figure 4). Prior to their concatenation, the hyperspectral image and the LiDAR data were georeferenced to the same projection system with a precision of 0.5 and 0.2 m respectively. They were then merged together in the same file thanks to the "Layer stacking" function implemented in the ENVI image processing software using the "Nearest Neighbor" resampling method.

Classification
The Mahalanobis distance (Richards, 1999), found to be most suitable in a majority of applications (McLachlan, 1999;Chang et al., 2004;Xiang et al., 2008;Govender et al., 2008), was used as distance function in the classification algorithm. This distance also has the advantage of taking into account the correlations of the data set (De Maesschalck et al., 2000). Indeed, hyperspectral data are measured over a large number and continuous spectral bands, each band containing potential redundant information. In addition, a Mahalanobis distance metric is scale-invariant and can adjust the geometrical distribution of data so that the distance between similar data points is small (Xing et al., 2003). This is supposed to enhance the classification performance.
To perform image classification, regions of interest (ROIs, areas of image that contain pixels of the same spectral characteristics representing the same seabed feature according to field observations) were created interactively. Ten habitat classes were considered to produce associated S. alveolata benthic habitat maps. Supervised classification clustered pixels in a dataset into classes based on previously defined training data. All pixels were classified to the closest ROI class using the Mahalanobis classification function implemented in ENVI software. For S. alveolata ecological status estimation, classification was applied in the same way but only on areas where it was previously detected.

Accuracy Assessment
The classification quality was assessed quantitatively using accuracy assessment based on comparison with field reference data none overlapping from the data used for training the classification algorithms. The kappa coefficient and the overall accuracy, the two statistics most prominently used (Foody and Mathur, 2006), were computed for classification accuracy assessment. Both the precision of localization of field data and the high spatial variability that exhibits the S. alveolata ecological status make the classification accuracy assessment of this latter very difficult to carry out. As it was stressed by Florén et al. (2015), low prediction accuracies could have several explanations but it is likely that most of them are driven by inaccurate positioning and/or the small-scale heterogeneity of the habitats as it is the case in the present study. Thus, many field data fell erroneously near the reefs when plotted on the image and thus could not be used. To overcome this issue, "hummocks" type observations were discarded because they do not constitute continuous formations and their localization lacks precision as they are occupying surfaces often lower than the image resolution (50 cm). For the S. alveolata reef status, data extraction from the classified image was done by searching the observed status within a radius of 2 m. This allowed to take into account the GPS accuracy (3 m) and thus to increase the number of kept samples for statistical evaluation. Overall accuracy and Kappa coefficient were thereafter calculated and compared for both datasets. For qualitative accuracy evaluation, we also plotted sample field data on ecological status classification map and made a comparison with the status estimated in the field.

Post-processing
As is often used in the post-processing of a thematic per-pixel classification map (e.g., Stuckens et al., 2000), a "Majority filter, " that generalizes and reduces single pixel misclassification, was applied using ArcGIS tool. This filter considers a neighborhood of pixels and assigns to the central pixel the semantic class most occurring in all pixels in the neighborhood. The kernel of the filter was the eight nearest neighbors (a three-by-three window) to the processed cell.

Spectral Reflectance Signatures
Field spectra were acquired to examine if S. alveolata bioconstructions have a sufficiently specific spectral signature leading to their identification and discrimination from the other benthic feature present in the study site. Recorded spectra are shown in Figure 5. Measures over 900 nm, less reliable due to the instrument sensitivity in this spectral region, have been dismissed.
Benthic features showed clear differences in their spectral reflectance shapes and/or their magnitudes (Figure 5). For vegetation types, reflectance values were low in the visible bands but showed a strong increase with wavelength above 673 nm reaching high levels in the entire infrared region ( Figure 5C). The spectral signature presents a significant jump between the absorption at 673 nm and the reflection in the near infrared. This jump is generally less marked on the spectral signature of the MPB population. Conversely, the absorption peak at this wavelength is stronger in the latter. The MPB also present several peaks of reflectance in the visible around 570, 602, and 646 nm when smaller specific absorption features were observed in reflectance spectra at 581 nm and 632 nm compared to macroalgae (Figure 5C). Brown algae reflectance spectra displayed lower values in the green part of the spectrum than the other vegetation types. For substrate types, spectra showed increasing reflectance with increasing wavelength in both visible and near infrared ranges (Figures 5B,D). Soft sediments displayed overall higher reflectance values than hard substrata. The highest values were recorded for shelly sediment (shells) in the infrared region. The lowest values were recorded for dark rock with an almost flat spectrum all over the range of wavelengths analyzed, without any obvious absorption feature. The spectrum corresponding to mud showed an absorption band at 673 nm similar to that of microphytobenthos. Oyster reefs exhibited greater variability in their spectral amplitude than S. alveolata reefs, especially in the infrared region, with some peaks of reflectance in the visible that were not always well discriminated. Regarding S. alveolata reefs, the spectral signatures of prograding bioconstructions showed a shape similar to that of sand with, however, a lower amplitude ( Figure 5A). With the state of degradation, their reflectance values increased and the shape of the spectrum was modified: degraded reef portions in retrograding reefs showed spectrum similar to muddy sediments, with an absorption peak at 673 nm.
According to statistical analysis, the first two dimensions of the PCA accounted for 71 and 25% of the total variation in the dataset, respectively, expressing together more than 96% of the data variability ( Figure 5F). Correlation biplots showed that the specific spectral signature of some benthic features, such as dark rock, shells, macroalgae as well as prograding S. alveolata reefs, was probably sufficient for their discrimination. However, other benthic features seemed to be more difficult to distinguish solely on the basis of their spectral signature, as for rocky substrata (Rock) with epifauna, fine sand sediment and retrograding S. alveolata reefs.

Benthtic Feature Classification and Validation
Preprocessed airborne data were classified using statistical supervised classifier based on Mahalanobis distance. As a result, dominant habitat types were automatically grouped into classes representing the pattern of Champeaux site ecosystem (Figure 6). Among the 10 benthic feature classes obtained, three were identified as soft substrates (fine sand, mud and coarse sediments), and two of them as hard substrates (rocky reefs and "other hard substrates"). Microphytobenthos was differentiated from macroalgae and displayed patches with varying levels of biomass. The pattern of class distribution also revealed the distinction between the two types of S. alveolata formations, large reefs developing on sand and adhering to rocky reefs on high shore.
The performance of the obtained classification was assessed using the confusion matrix method (Figure 7). In this approach, classified benthic habitats were compared with field observed benthic features, resulting in an overall accuracy of 88% and a Kappa of 0.85. Khi 2 test (p < 0.0001, α = 5%) corroborated the high performance of the used classifier for S. alveolata formations and its associated benthic features. Few confusions were recorded, however, mainly between S. alveolata on sand and bare sand sediments or between S. alveolata on rocky substrata and coarse sediments, but in much lower proportions than correctly classified pixels.

Remote Estimation of S. alveolatta Reef Ecological Status
In order to continuously map ecological status, supervised Mahalanobis distance based classifier was applied to each pixel of merged hyperspectral and LiDAR data. Prograding and retrograding phases were provided with 0.5 m spatial resolution. The ecological status map obtained exhibited an heterogeneous distribution with a co-presence of prograding and retrograding phases whatever the location (Figure 8). Nevertheless, prograding phase seemed to be seaward dominant, where retrograding phase was mainly present in the southeast part of the study site and closed to the shoreline.
The calculated overall accuracy and Kappa coefficient varied with the number of pixels used in the evaluation. Values of 0.77 and 0.55 were obtained, respectively, when using 22 samples, but increased to 0.93 and 0.86, respectively, taking into account a larger field dataset (113 samples). The S. alveolata area in the Champeaux site was estimated at 5.52 and 6.76 ha, respectively with and without post-processing spatial filtering (Table 1). Sabellaria reefs in this site seemed to be slightly dominated by retrograded phase (58 and 53% respectively).

DISCUSSION
By using an automatic supervised classification algorithm, we demonstrated the ability of combining hyperspectral and LiDAR signals to (i) provide a high-resolution (0.5 m) map with 10 identified habitat classes describing S. alveolata reefs and their associated benthic features; (ii) accurately discriminate two S. alveolata reef forms, large reefs developing on sand and smaller veneers encrusting rocky shore areas, as well as its main ecosystem associated habitats, and (iii) reliably and remotely estimate the spatial structure of these bioconstructions ecological prograding vs. retrograding phase. Further details are presented in the following sections.

Seabed Type Spectral Discrimination
The visual inspection of field-measured spectra acquired in the Mont-Saint-Michel bay (France) shows that most of the spectral profiles are characterized by numerous absorption bands, mainly due to the presence of pigments. Previous studies showed that differences in pigment mixtures mostly FIGURE 6 | Benthic habitat map of Champeaux site obtained from the hyperspectral airborne imagery using supervised classification based on Mahalanobis distance. explain differences observed in specific spectral signatures (Guillaumont et al., 1997;Méléder et al., 2003;Jesus et al., 2014). The absorption band observed at 673 nm is characteristic of all plants due to the ubiquitous presence of chlorophyll a, while chlorophyll c is responsible for absorption bands at 581 and 632 nm (Kuczynska et al., 2015) as displayed by brown algae and microalgae spectra. Chlorophylls a and c are dominant pigments in diatoms, haptophytes (Haptophyta or Prymnesiophyta) and to a lesser extent in some dinoflagellates (Kuczynska et al., 2015;Coupel et al., 2015). These two pigments typify microphytobenthos assemblages dominated by diatoms (Méléder et al., 2003;Jesus et al., 2014). These observations are coherent with taxonomic analyses of biofilm samples that show communities highly dominated by diatoms. Apart from chlorophylls, diatoms contain carotenoids in their pigment profile that are mainly composed of carotene, fucoxanthin and diatoxanthin under light conditions (Kuczynska et al., 2015). Carotenoids exhibit an intense absorption between 400 and 500 nm, fucoxanthin having a slightly higher absorption band (460-570 nm) than the others (Kuczynska et al., 2015). The dominance of fucoxanthin in the pigment profile of brown algae (Kirk, 1977) can also explain the main specificity of the spectral signature of this group. Lower reflectance spectra were observed for brown algae compared to other vegetation types in the blue-green part of the spectrum. This mainly results from a high concentration of fucoxanthin per cell, exhibiting high absorption in the range 460-570 nm (Kuczynska et al., 2015).
Intertidal substrates have spectral characteristics that seem to be related to both their mineral and chemical composition, as well as their structural properties (grain size and sorting parameters). Typical reflectance spectra of substrates reveal spectral shapes characterized by increasing reflectance with increasing wavelength, similar to those recorded by previous studies (Guillaumont, 1991;Méléder et al., 2003;Bajjouk et al., 2019a). Dark rocks are however an exception. Their gray black color seems to explain the very weak reflectance on all the visible and infrared spectral range. The slight absorption that appears at 673 nm is characteristic of vegetation and is observed in muddy sediment spectra, due to adhering alga cells in poor physiological state and chlorophyll a degradation products that also participate in absorption at 673 nm (Méléder et al., 2003).

Coherence of Benthic Feature Distributions
Multivariate analysis showed that although spectral richness, confusions between certain features (e.g., retrograding S. alveolata reefs vs. sand sediments) are possible. The integration of LiDAR sensor information into the hyperspectral image provides an additional characteristic on the nature of the topographic complexity that can help for improving the discrimination between benthic features. As discussed by Chang et al. (2004), materials in real world imaging could occupy a full pixel or be either mixed with other targets such as background or embedded in a single pixel as a subpixel target. Unlike SAM (Spectral Angle Mapper) and Euclidean distance that are generally used for pure pixel-based classical image processing, Mahalanobis distance does not take into account the spectral correlation among pixels. This characteristic of Mahalanobis is very useful and brings crucial information in identification or discrimination of materials in real images. The S. alveolata reef being spatially very heterogeneous in terms of abundance and structure, even at a high resolution, we therefore favor the use of Mahalanobis distance for image supervised classification.
The highly accurate overall accuracy (88% with a Kappa coefficient of 0.85) is consistent with the outputs of previous work that also successfully used this distance for image classification (Chang et al., 2004;Govender et al., 2008). Obtained results also reinforce the recommendations encouraging integration of direct observation field techniques with remotely sensed data to capture detailed habitat information at large scales, thus avoiding trade-offs between three types of scale: spatial, temporal and thematic (Lecours et al., 2015;Rhodes et al., 2015;D'Urban et al., 2020). The benthic habitat map of Champeaux obtained from the hyperspectral airborne imagery using Mahalanobis supervised classification allows the delineation of a higher number of discriminated habitats (ten) than what was identified using RGB (Red-Green-Blue) images and digital surface model, as shown for S. alveolata reefs growing atop of soft sediment exclusively (Collin et al., 2019). The most significant information that the pattern of class distribution reveals is the distinction between two types of S. alveolata formations; large reefs developing on sandy bottoms and smaller veneers encrusting rocky shore areas. This is the first time that the two types of intertidal S. alveolata reefs are automatically mapped and discriminated, revealing how relevant LiDAR information is. As stressed by Brown (2004), elevation data has the ability to enhance discrimination of intertidal species, since topographic metrics (e.g., elevation, rugosity or slope) can be computed and used as proxies for their position in relation to the tidal cycle. Visual inspection of the map shows, however, that pixels a priori automatically labeled "S. alveolata on rock" are actually located in an area (south-eastern part) without rocky substrata. This can be explained by the fact that bioconstructions in this specific area have a veneer shape but lack 3D structure, as is generally the case with S. alveolata reefs developing on sandy bottoms. Regarding reef ecological status distribution, the general orientation of the reef exhibits spatial pattern parallel to the coast. This is most likely due to the hydrodynamic conditions. Tidal currents at Champeaux are often fast and opposite to near shore direction. As a resulting effect, sand transport and probably reef development extend perpendicularly to the main current direction (Gruet, 1986). We also showed here the benefits of using specific spectral indices to improve discrimination between algal groups. By not being limited only to the commonly used NDVI-a generic index to detect vegetation-and by calculating a spectral index specific to diatoms, it was possible to separately map macroalgae and microphytobenthos populations. The presence of benthic microalgae in the sediments landward is consistent with previous works confirming the importance of S. alveolata bioconstructions in enhancing the local benthic primary production (Jones et al., 2018). Jones et al. (2018) have also demonstrated how the engineered sediment significantly affects the sedimentary characteristics of the associated reef sediment, specifically grain-size, organic matter content and ultimately microphytobenthos. 3-dimensional engineered structures act as a physical barrier: wave energy is reduced and tidal currents are slowed, increasing sedimentation between reef structures and coastline (Noernberg et al., 2010;Desroy et al., 2011;Jones et al., 2018).

Contribution of Hyperspectral and Lidar Combination to Reef Status Estimation
Estimating the ecological status of benthic habitats is a requirement of the European Union Habitat Directive, where reef-habitats are listed as priority (Council Directive 92/43/EEC). Sabellaria alveolata reefs of varying abundance are common on all European coasts and all developments toward automatic assessment of ecological status of this habitat would tremendously help large scale ecological assessments, while limiting time consuming and costly local surveys. Previous local investigation of reef ecological status based on reef types (veneer, hummock, platform) and phases (retrogradation and progradation) are costly in human resources and showed high inter-operator variability (Desroy et al., 2011). In this study, prograding and retrograding phases accurate map was provided with 0.5 m spatial resolution. An overall classification accuracy and Kappa coefficient of respectively 0.93 and 0.86 can be reached when using a validation against field data that takes into account a large sampling dataset. The sample field data plots on ecological status classification map also showed good coherence in terms of distribution pattern. This shows that the combination of LiDAR and hyperspectral reflectance variables had a robust explanatory capacity of the variability of S. alveolata reef ecological status. Indeed, LiDAR has also proven to be effective in characterizing health status of other ecosystems, namely shallow coral reefs. Using an artificial neural network classification, Collin et al. (2018) have shown that the dual combination of LiDAR surface and intensity variables had a robust explanatory power to predict the variability of coral reef states, with an overall accuracy of 0.75. We showed here that combined remote-sensing techniques can provide a high-confidence assessment of the prograding or retrograding phase of this habitat. This definitely calls for larger applications to other coastal areas.
The ecological status of S. alveolata reef was mapped continuously throughout Champeaux site. The 50 cm spatial resolution targeted in the present study outperformed the 75 m in Desroy et al. (2011). The latter authors used in situ expert-dependent estimation, while we based our method on an automatic objective processing for reef status estimation. Spatial patterns emerging from ecological status classification, also display a barrier reef on sand with a vigorous seaward front core, while the retrograded and sparse reefs are mainly located in the south-east part and next to shoreline. This seems to be consistent with previous work carried out by Collin et al. (2019) that have described a strong polarization of S. alveolata reefs with high abundance on forereef and low abundance on backreef. The forereef is subject to strong hydrodynamic and potentially higher coarse sediment resuspension, leading to increased tube-building activity, while the back reef mainly expands into more sheltered and muddier environments, less favorable to S. alveolata development (Bonnot-Courtois et al., 2008). Environmental impacts may also be exacerbated by human activities which shape the spatial pattern of reef status. Degraded reefs are mostly present in areas that are easily accessible to pedestrians. As it was observed for Mediterranean S. alveolata reef (Lisco et al., 2020), this degraded localization most likely drives the direct mechanical disturbances related to higher trampling and recreational fishing, leading to massive degradation of some reef areas (Dubois et al., 2002).

Perspectives for S. alveolata Reef Monitoring and Management
The trained classification yielded very accurate benthic feature identification. However, some pixels were wrongly classified. Processing improvements are thereby still needed to improve discrimination between mixed up classes. The strategy of field data acquisition can also be optimized to provide a sampling as balanced as possible between the different classes targeted. This can increase the efficiency of both calibration and validation steps. Another crucial criterion is the position accuracy, due to the high heterogeneity of the seafloor: an inaccuracy of only few meters often means a different substrate, a different species or a different ecological status (Florén et al., 2015). In future studies, field data should be collected and accurately positioned, by using for example DGPS (Differential Global Positioning System), to facilitate and make more robust accuracy assessment achievement. Placing benthic markers that can be identified in the imagery, as suggested by Mumby et al. (2004), is also a possibility to address this issue. Moreover, other classification algorithms, such as deep learning prediction (LeCun et al., 2015) should also be tested to examine their powerfulness compared to Mahalanobis distance based classification. This may allow for a better estimation of the areas occupied by benthic habitats, more particularly S. alveolata bioconstructions. Such a metric can be used in spatiotemporal dynamic monitoring or even integrated into ecological status indicators, as is the case for other habitats in European Directive framework (Duarte et al., 2017).
Regarding S. alveolata reef status estimation, the work reported here represents a first step which illustrates a promising approach, based on the combination of hyperspectral and LiDAR data. The review of coral reef monitoring and assessment technologies carried out by Obura et al. (2019), emphasized that satellite data combined with airborne imagery and airborne LiDAR data are now routinely used in identifying and monitoring marine protected areas as a new ways to assess benthic habitat. Quantitative estimation of S. alveolata ecological status index is an attainable goal: further work, particularly a linear mixing model approach (Bioucas-Dias et al., 2012;Bajjouk et al., 2019b), may prove to be relevant to estimate S. alveolata epibiont abundance-a key element in reef dynamics-and possibly further develop a finer metric than the prograding vs. retrograding status qualitative estimation. Integrated applications of these remote sensing techniques are still lacking to target S. alveolata ecosystem spatial ecology. As wisely pointed out by D' Urban et al. (2020), there is a strong need for policymakers to have access to high-quality environmental information, in order to make objective conservation decisions leading to ecosystem service maintenance. In addition to helping gain a better understanding of ecological processes, 3D remote sensing approaches constitute an efficient tool for recording accurate environmental key information at both large scale and high spatial resolution. These characteristics are of great interest, especially for the study of rapidly changing environments such as those associated with S. alveolata reefs. We also highlight the relevance of such information in improving communication and targeted awareness-raising through their ability to provide virtual reality system overviews.

CONCLUSION
The present study is an initial step toward designing a new alternative operational approach targeting broad-scale, highresolution, non-destructive monitoring of S. alveolata reefs. The main goal is to assess the potential of optical imagery, more specifically hyeperspectral and LiDAR data, for a highresolution (0.5 m) spatial characterization of the S. alveolata reef ecosystem.
The main benthic features showed clear differences in their field spectral reflectance shapes and/or their magnitudes. However, statistical analysis revealed the difficulty to distinguish between some types solely on the basis of their spectral signatures. The joint analysis of the LiDAR sensor information with the hyperspectral imagery provides an additional characteristic on the nature of the topographic complexity, thus enhancing the discrimination capability.
By using an automatic supervised classification, we have shown for the first time the potential of hyperspectral and Lidar signal combination (i) to accurately discriminate the two S. alveolata reef forms, large reef developing on sand and smaller veneers encrusting rocky shore areas, as well as its main ecosystem associated habitats and (ii) to reliably and remotely estimate the spatial structure of these bioconstruction prograding and retrograding dynamic phases as a proxy of ecological status. This latter is commonly required by both scientists, for better understanding of ecosystem process, and managers, for better implementation of European directives.
The study also highlighted the importance of field positioning accuracy to perform a robust statistical validation, mainly in high spatial heterogeneity as displays the S. alveolata reefs. Further work has therefore been identified to improve obtained results and to provide quantitative information on S. alveolata reef ecological status. Such information, if is provided on a large scale and at a sufficient frequency, will allow (i) improved understanding of physical and biological structural changes of these complex bioconstructions, (ii) better prediction of their spatio-temporal evolution and (iii) elaboration of appropriate measures for their preservation and socioeconomic management.

DATA AVAILABILITY STATEMENT
The online publication of data is planned at the end of the BIOHERM project in 2021. Requests to access the datasets should be directed to touria.bajjouk@ifremer.fr.

AUTHOR CONTRIBUTIONS
TB, LD, and MD designed the project strategy and performed funding acquisition. TB supervised and managed the project. TB, CJ, and SD collected field data. AD and CJ carried out the laboratory analysis. TB performed data processing and experimental validation and drafted the manuscript. SD validated the classification models. TB, CJ, SD, LD, and MD contributed to manuscript writing and editing. All authors shared the responsibility for contributing to the final version of the manuscript.

FUNDING
This work was funded by the CNES (French National Center for Space Studies) as part of the space scientific research TOSCA program to support the CHIMERE sensor and the Biodiversity mission. The study was carried out within the framework of the BIOHERM project (Ref. 4911/4500062124).

ACKNOWLEDGMENTS
We are very grateful to Yves Pastol from the SHOM (Hydrographic and Oceanographic Service of the Navy) for providing Lidar data acquired within the framework of Litto3D program. We also thank Pierre-Olivier Liabot and Simon Martin from Ifremer, for their help on field data acquisition, as well as Hytech Imaging for airborne data acquisition. We wish to warmly thank Amelia Curd, native English-speaker, for proofreading the manuscript. Finally, we would like to gratefully thank the two reviewers for their valuable comments and suggestions.