Abstract
Seagrasses are undergoing widespread loss due to anthropogenic pressure and climate change. Since 1960, the Mediterranean seascape lost 13â50% of the areal extent of its dominant and endemic seagrass-Posidonia oceanica, which regulates its ecosystem. Many conservation and restoration projects failed due to poor site selection and lack of long-term monitoring. Here, we present a fast and efficient operational approach based on a deep-learning artificial intelligence model using Sentinel-2 data to map the spatial extent of the meadows, enabling short and long-term monitoring, and identifying the impacts of natural and human-induced stressors and changes at different timescales. We apply ACOLITE atmospheric correction to the satellite data and use the output to train the model along with the ancillary data and therefore, map the extent of the meadows. We apply noise-removing filters to enhance the map quality. We obtain 74â92% of overall accuracy, 72â91% of userâs accuracy, and 81â92% of producerâs accuracy, where high accuracies are observed at 0â25Â m depth. Our model is easily adaptable to other regions and can produce maps in in-situ data-scarce regions, providing a first-hand overview. Our approach can be a support to the Mediterranean Posidonia Network, which brings together different stakeholders such as authorities, scientists, international environmental organizations, professionals including yachting agents and marinas from the Mediterranean countries to protect all P. oceanica meadows in the Mediterranean Sea by 2030 and increase each countryâs capability to protect these meadows by providing accurate and up-to-date maps to prevent its future degradation.
Similar content being viewed by others
Introduction
Seagrasses, the submersed marine flowering plants, are undergoing substantial decline due to direct and indirect anthropogenic activities and climate change. The recent publication1 estimated the global seagrass area to date as 160,387 to 266,562 square km, which is significantly lower than the previously reported global seagrass spatial distribution ranging from 177,000 to 600,000 square km2,3. Their loss rates accelerate from 0.9% yearâ1 in the 1940s to 7% yearâ1 toward the end of the twentieth century, making them one of the most threatened ecosystems on Earth4,5. Since 1960, the Mediterranean seascape lost 13â50% of the areal extent of its most common and endemic Posidonia oceanica seagrass meadows6. If the recent heat waves in the Mediterranean Sea continue, this seagrass may disappear within 100â150 years7. Recently, efforts to protect P. oceanica have led to the creation of the Mediterranean Posidonia Network (MPN, https://medposidonianetwork.com/) including more than 50 stakeholders from 10 countries, which seeks to effectively protect 100% of P. oceanica by 2030 by introducing innovative tools and raising awareness at the local, national, and international levels about the importance of this species8. However, the lack of a multi-temporal monitoring system over large areas, and identifying the hotspots to mitigate the rapid degradation caused by anthropogenic pressures are becoming a major concern for Mediterranean coastal managers, affecting the conservation program to protect this species.
Posidonia oceanica covers 50,000 square km of coastal sandy and rocky areas consisting of 25% of the Mediterranean Sea bottom at 0â40 m depth9,10. The meadows are home to about 20% of Mediterranean sea life species and can accommodate up to 350 animal species per hectare11,12,13. They contribute to different coastal processes like sediment deposition, and attenuate currents and wave energy14. They are globally significant carbon sinks with greater organic carbon density compared to estuarine mangroves, peatlands and tropical forests15,16, thus playing a key role in climate change mitigation17. However, the meadows are susceptible to increasing heat waves and climate change18,19,20. Besides, the meadows are under substantial threat by different anthropogenic activities, i.e., 22% of construction of coastal infrastructure, 23% of water pollution, 14% of invasive species, 18% of Fishing, 5% of shipping and 18% of modifications of marine currents and hydrography21 (Fig. 1). These threats, coupled with the slow horizontal growth of P. oceanica (1â6 cm/year), underscore the critical importance of preserving the meadows in their original state22. The European Union (EU) Habitats Directive 92/43/EEC has listed P. oceanica meadows as a priority habitat. P. oceanica is also a protected species in the Natura2000 networks within the EU because of its significant influences on the surrounding biological, biogeochemical and physical processes in the littoral23.
A successful conservation programme for P. oceanica involves knowledge about the spatial extent of the meadows, understanding processes of degradation24; identifying causes of degradation25; designing and implementing effective monitoring programmes to manage, restore, or create these seagrass meadows26,27. This demands fast, cost-effective and validated methods to observe and monitor the present state, spatial pattern, and dynamics at appropriate spatial and temporal scales28,29. Multiple studies30,31,32,33,34,35,36,37,38,39 presented the potential of remote sensing data in mapping the extent of this seagrass species, complementing the traditional in-situ monitoring system. However, they were mostly for proof of concept to utilize satellite-derived data and were hardly capable of supporting the continuous monitoring system since they are computationally expensive, and require validated methods that can map across large areas with adequate accuracy and precision.
With growing archives of freely accessible earth observation data of adequate spatial and temporal resolution, along with fast-growing computational processing capabilities, we now have the potential tools to support mapping and monitoring of P. oceanica over large regions. To our knowledge, there is no use of remote sensing techniques coupled with deep-learning-based artificial intelligence, integrated within a pipeline technique that facilitates reproducibility and scalability to monitor the multi-temporal spatial extent of this seagrass species, thus supporting the conservation activity. In this study, we present an operational approach aiming to fill this gap. Our approach corresponds to a Scientific Exploitation Platform (SEP) dedicated to the extraction and processing of Earth observation (EO) products by adding one more layer on top of the existing tools to take the exploitation of the raw satellite data one step further. The proposed pipeline used Sentinel-2 data for the artificial intelligence deep-learning neural network (DLNN) to produce multi-temporal maps of P. oceanica. Here, we present the results for four major islands of the Balearic Islands, i.e. Menorca, Mallorca, Ibiza, and Formentera in Spain, and the Maltese Islands (Fig. 2). For these regions, we have trained the DLNN model and evaluated its performance under the scheme of long-term monitoring and evaluation to support conservation and ecosystem management to minimize climate change and anthropogenic impacts. This study should be of interest to MPN as well as the scientists, conservationists and/or coastal managers, who are responsible for the challenging short-term management and long-term policy decisions to protect this type of natural resources.
Results and discussion
Maps of P. oceanica
We have detected 577.29 and 74.27 square km of P. oceanica respectively in the Balearic Islands and the Maltese Islands during 2021 (Fig. 3). Our proposed method is based on pixel-by-pixel classification schemes, and the final output map has 10 m spatial resolution. We have stacked 3â9 atmospherically corrected clear Sentinel-2 level-2 images per tile (Supplementary Table S1) and applied the DLNN model with the corresponding in-situ and bathymetry data. The yearly stacked image using the median value for each pixel across the considered images reduced the potential misclassification due to non-persistent and seasonal variations in algal growth that may have similar signals. It was evident that depending on the region, our proposed approach was capable of mapping P. oceanica with 74â92% of overall accuracy, 72â91% of userâs accuracy, and 81â92% producerâs accuracy (Table 1). When a satellite pixel contained a region with sparsely distributed or unhealthy P. oceanica (according to the in-situ data), it was classified as false negative (FN) due to weak spectral signal. In contrast, when majority parts of the pixel contained healthy P. oceanica, it was classified as false positives (FP) due to strong spectral signal. However, in the north and northeast parts of Ibiza, we have observed swell noise in all the available Sentinel-2 images that caused residual noise in the yearly stacked image and resulted in misclassification. In contrast, we have noticed a different scenario in Malta. In the southeastern part of Malta, we observed the presence of spectral signals coming through the water column, which were identified by the DLNN model as P. oceanica. In this context, there could be two possible explanations, i.e., incomplete in-situ database, where the identified patches of P. oceanica were not included; consequently, it was depicted as FP during accuracy assessment; or, the presence of another seagrass species with similar signals. In case of the second one, a more detailed spectral inspection of other visible and near-infrared bands along with the RGB bands could be a solution to distinguish different habitats, which will be a topic of future studies. However, it should be noted that the in-situ data used in this study came from several years of field campaigns and might not be up-to-date in every region. As a result, in some places, the FP and FN identified in the final maps could be a true presentation of the current scenario and stressed the need for further investigation.
Our model obtained a higher userâs and producerâs accuracy at 0â25Â m depth (Table 2), which was the most representative depth limit of the seabed detection using Sentinel-2 satellite in the Mediterranean Sea37,40. As the water was transparent and Sentinel-2 could see accurately enough, the model showed higher accuracy at this depth. With increasing depth and light attenuation, there could still be some P. oceanica living on the sea bottom up to 40Â m depth41,42 but the satellite could not see them accurately enough. This explained the modelâs performance with comparatively low accuracies at that depth.
One of the advantages of our proposed method is the classification scheme. We trained the DLNN model to learn how the colour of the sea bottom with and without P. oceanica varied with depth, thus included the effect of water column implicitly. Precisely, instead of correcting the data for water column absorption, we implicitly include it in the forward modelling of the inversion problem for the numerical stability in the large depth, unlike the available literature39 that corrected water column absorption and scattering explicitly for a depth limit of 21 m. Besides, we designed the DLNN model for pixel-by-pixel classification per tile. As a result, our model was easily adaptable to other regions. We obtained accurate and reliable results while training the DLNN model with in-situ data from the target region, however, without in-situ data, the model was still capable of producing reliable maps which could provide a first-hand overview. To apprehend this, we reproduced the map of P. oceanica in the Maltese Island using the training dataset from Formentera (Fig. 4). The final output showed an overall accuracy of 64% with a userâs accuracy of 73% and a producerâs accuracy of 69%. As expected, the model performance was lower than what we observed in Table 1, however, it was still scientifically acceptable and could provide a general idea of the spatial extent of the meadows. It was evident that in both cases (Figs. 3i,j, 4) the model identified P. oceanica in the southeastern part of Malta, implying strong evidence of the presence of the seagrass meadows.
Yearly monitoring of P. oceanica and its change assessment
We applied our proposed approach to perform a time-series analysis (Fig. 5 and Supplementary Fig. S3, and Supplementary Table S2) in Formentera, as this region was highly renowned for conservation and restoration activities. While the seagrass patches mostly looked stable during 2017â2021 (â0â changes per pixel in Fig. 5b), we detected changes in the northern and northeastern parts of Formentera. Especially, in the northern part of the island at a depth ofââ¥â30 m, we detected a loss of P. oceanica meadows until 2019 and then a subsequent gain since 2020 (Supplementary Fig. S3). It should be noted that during 2018 and 2019, the quality control for stacking was challenging due to poor quality of the available satellite images, which might results in the above-mentioned change detection at the optical depth limit for Sentinel-2. However, the change assessment between 2017 and 2021 (Fig. 5a) showed that the area of the meadow in the Formentera island increased by about 8 square km and decreased by about 1 square km, with a net gain of 7 square km, which might be a result of the current conservation and restoration activities. This showed the promising capabilities of our proposed approach to support monitoring of the meadows in the short and long term.
Application in conservation and ecosystem management
In the recent past, the concept of natural recolonization (cutting and seedlings) has been introduced to restore P. oceanica meadows in the Mediterranean basin. However, most of them failed due to either choosing the wrong site where the species had never been before or a lack of long-term monitoring system due to economic constraints43. In addition, conservation and ecosystem management demand explicit and up-to-date knowledge to deal with real-time situations. Timely and reliable maps with a large spatial and temporal coverage can provide substantial ground in capacity building and making decisions. Unlike the available methodologies that were resource expensive44 and/or developed for a local scale45, our proposed method can generate maps of P. oceanica with adequate accuracy in a cost-effective manner. This can add value to the current restoration programmes by identifying suitable sites and monitoring them regularly for short-term goals and long-term evaluation. Besides, the information generated through our approach can support conservationists and ecosystem managers in making rational decisions to promote sustainability. This is particularly true for the MPN, which brings together different stakeholders such as, authorities, scientists, international environmental organizations, professionals including yachting agents and marinas from the Mediterranean countries. Being a part of the network, our approach can support the main objectives of this network, which is to protect all P. oceanica meadows in the Mediterranean Sea by 2030, and to increase each countryâs capacity to protect these meadows by providing the necessary tools to prevent its future degradation.
In this context, we present an example in Fig. 6. The figure demonstrates boat pressure over the seagrass patches in Formentera and Ibiza during JuneâAugust 2021, a peak tourist season. Each hexagon mesh presents the number of ships detected using Sentinel-1 satellite data. It was evident that the channel between Formentera and Ibiza had the highest pressure within this period, thus showing the promising capabilities to identify the hotspots of P. oceanica for conservation and sustainable coastal resource management.
This study facilitates the continuous and long-term monitoring of P. oceanica in the Mediterranean Sea in a fast and efficient manner by utilizing freely available high-resolution Sentinel-2 satellite imagery, low-resolution bathymetry data and open source free software (Python, Docker, etc.) to run the whole process automatically with minimal human interaction. Alongside the data-enriched regions, our approach can be used in data-scarce region, with environmental and other constraints that make it impossible to collect in-situ data, thus providing a first-hand overview. Besides, for precise mapping, our approach can be adaptable to high-resolution data, i.e., high-resolution bathymetry or commercial satellite data with RGB bands. Finally, our approach can be a stepping-stone for recording the long-term impact of climate change as well as different environmental drivers and stressors on the seagrass meadows following a standardized and uniform method, and ameliorate understanding of the connectivity between oceanography and seagrass ecosystem, spatial and temporal management and protection, and finally mapping other benthic ecosystems for conservation and management.
Methods
Scientific exploitation platform
To map P. oceanica in a fast, reliable and efficient way, we developed an operational pipeline (Supplementary Fig. S1) within a SEP. SEP integrated the hardware/software infrastructures that supplied the computing and storage resources needed for the exploitation and provision of the tools to manage EO datasets in a distributed environment. The platform comprised a Kubernetes cluster with many nodes hosted on-premises. The pipeline consisted of several scientific modules, such as data acquisition, data pre-processing, followed by artificial intelligence-based DLNN, and post-processing. Within the SEP, all these modules could be executed isolated or sequentially within Docker containers, encapsulating the entire environment required for the scientific workflow. This ensured scalability and the pipelineâs results could easily be reproduced with consistent outcomes. The pipeline also integrated a quality control system where a Network-file-system (NFS) volume held the raw data and the processed products, which could later be downloaded and analysed within the SEP web interface.
Data acquisition
Sentinel-2 Level-1C data were downloaded from Google Cloud Storage using their Application Programming Interface (API). This data corresponded to the top of the atmospheric reflectance (TOA) for 7 Sentinel-2 tiles covering the Balearic Islands, and 1 tile covering the Maltese Islands (Supplementary Table S1). After visual inspection, images with high cloud coverage or sun glint effects were discarded from further processing. While mapping P. oceanica, bathymetry data plays a significant role, as this species is highly depth-dependent. In this study, the freely available bathymetry data at 100Â m spatial resolution were obtained from the European Marine Observation and Data Network (EMODnet, https://www.emodnet-bathymetry.eu/). Besides, the in-situ data to develop and train the DLNN model were collected from the local government of the Balearic Islands46 (https://atlasposidonia.com/en/home/) and the Marine Database of Environment and Resources Authority (ERA) in Malta47,48. These data corresponded to single cartography in the GIS layer combining several years of field campaigns. In the case of Balearic Islands, the in-situ cartography were prepared based on the photointerpretation of aerial photographs between 0 and 5Â m depth, side scan sonar campaign between 0 and 35Â m depth, reprocessing of the information available from the LIFE Posidonia project, OCEANSNELL, Consell de Ibiza, Ecocartographies, etc. during 2008â201946,49. These data corresponded to different marine habitats following the nomenclature and coding of the Standard List of Marine Habitats of Spain (LPHME). For the Maltese Islands, the in-situ cartography corresponded to the year of 2017â2019, having the information about the habitat of P. oceanica.
Data pre-processing
Pre-processing is a crucial step while producing benthic maps using satellite data for artificial intelligence. It is even more important while using Sentinel-2 data for aquatic applications as it has low-signal-to-noise-ratio that requires a good atmospheric correction (AC) to separate the top-of-atmospheric (TOA) reflectance observed by the satellite into the signal from the atmosphere and the signal from the surface (bottom-of-atmosphere, BOA)50. In this study, we applied the ACOLITE (v2020) AC processor developed by the Royal Belgian Institute of Natural Sciences (RBINS). Based on the literature and the analysis results, we chose ACOLITE because it better reproduced the shape of the reflectance spectra than the other widely used AC processors, i.e., C2RCC (a SNAP plugin). As ACOLITE processor is image-based, it does not require any in-situ measurements. By default, ACOLITE performs AC using the âdark spectrum fittingâ (DSF) algorithm for aerosol corrections. However, this algorithm avoids severe glints for which the glint effects can still be found in the DSF-derived surface reflectance51. As sun glint was a prominent optical phenomenon in the study region, especially during spring and summertime, the optional glint correction within ACOLITE was also applied. The land and clouds were masked during the AC processing to retain the water-leaving reflectance. Therefore, the data were inspected visually. In some cases where high clouds (i.e. cirrus, cirrostratus, and cirrocumulus) and haze were left unmasked, and when these were found to degrade the region of interest near the coast, the images were discarded. Occasionally, images with high turbidity caused by river/land discharges or algae bloom were also discarded as they significantly affected the colour of the water and made it impossible to observe through the water column. These strict criteria only admitted images that were of very good quality. Therefore, to remove boats, very large waves, and other anomalies, as well as to reduce noisy patterns and artifacts in order to improve the overall quality of the image, a pixel-by-pixel median stacking was performed for the RGB bands of all the good quality atmospherically processed data (3â9 images per year). In this case, we chose median stacking as it was less affected by outliers and other anomalies compared to mean stacking, by preserving the features present in the majority of the images by taking the median value of each pixels across all the good quality images. In this way, as P. oceanica is a slow-growing species (1â6Â cm/year of horizontal growth), the chances of overestimation of the seagrass extent because of the growth was negligible. Besides, if there was a sudden destruction of the meadows due to anthropogenic activities, it would likely lead to a drastic reduction in seagrass cover, resulting in a lower median value, thus reflecting the actual change more accurately.
The EMODnet bathymetry data was reprojected and resampled to 10Â m (using the bilinear interpolation technique) to match the projection system and the spatial resolution of the Sentinel-2 data. The data was cropped after 40Â m depth to remove open and deep ocean pixels to avoid potential FP of P. oceanica, as both classes had similar spectral values37 (Supplementary Fig. S2).
The marine habitats of the in-situ data were classified into two classes, i.e. "PO" (indicating the habitat of P. oceanica) and "Non-PO" (indicating other than the habitat of P. oceanica). In-situ data corresponding to land, offshore, and shallow water with a sandy/rocky habitat or other seagrass habitats were included in the Non-PO class. The PO and Non-PO classes were labelled as 1 and 0, respectively, and rasterized using the same projection system of Sentinel-2 scenes.
Artificial intelligence-based DLNN approach
The identification of P. oceanica is a binary classification problem where pixels correspond to either PO or Non-PO. Therefore, the aim of using a DLNN model was to train a function based on the given satellite and bathymetry dataset (X) as well as target in-situ data (Y) to output a predicted class (Ŷ) of PO and Non-PO, which should match as closely as possible with the in-situ data (Y). In this study, we have randomly split the data into 80/20 ratios per Sentinel-2 tile, where 80% data were used for training and the rest 20% were used for validation of the model.
The training was done through two processes, namely, forward propagation and backpropagation to recognize patterns based on the given dataset. In forward propagation, the process moved in the forward direction through the neural network (NN) to produce a final value at the output layer, whereas at the same time, in backpropagation, the process moved from the output layer to the input layer to improve the weight (w) value generated in the forward process. Forward propagation was a value-processing step from the input layer that involved calculating a linear combination per neuron using adjustable linear coefficients and then calculating an activation function by performing a non-linear transformation for the corresponding neuron within the layers. The activation functions of all the layers except the last one were the rectified linear unit activations (ReLU), meaning that the output was set to zero if it was negative. Hence, the output coefficients from the neurons never had negative values. In the DLNN model, the last layer was a binary classifier, which means that it classified pixels as belonging to either class 1 or class 0. In this case, the activation function was a sigmoid function, which predicted the probability (Ŷ) of a pixel being class 1 or 0 based on a threshold value of 0.5.
To optimize the training process, we performed mini-batch processing, where we split the training set into smaller sets and implemented gradient descent on each batch chronologically to make the algorithm work faster. Therefore, based on the input data, each 10 m pixel of the Sentinel-2 image was classified as containing PO (class 1) or Non-PO (class 0). The result, Ŷ, was compared to Y, and the difference was quantified in a cost function. Minimizing the cost function equates to making Ŷ as similar as possible to Y. The differences between Ŷ and Y per pixel allowed the calculation of the gradient that quantified how much each coefficient in each neuron must be altered to match Ŷ to Y. This gradient was then used to update the coefficients in the neurons of all the layers through backpropagation. This was an iterative process, and therefore, at each iteration, a new prediction was made, compared to the target Y, a gradient was calculated, and the coefficients were updated. This process was carried out continuously until it produced the smallest loss value, and the output was saved in the form of a file containing the weight (w) and bias (b) values. Afterward, the derived set of coefficients from the DLNN model was used to map P. oceanica.
Post-processing and accuracy assessment
In the post-processing step, we evaluated several filters, i.e. sieve (replacing small isolated groups of pixels with the pixel value of the largest neighbour), median (ranking values of the pixels in the moving window with a specified radius and taking the middle-ranking value), bilateral (replacing the value of each pixel with a weighted average of the pixel values in the moving window), and principal component analysis (removing noise by discarding the last principal components that explain the least of variance). Among them, we chose the median filter with 3âÃâ3 pixels moving window as it significantly reduced the noise while preserving the shape of the seagrass patches. We applied the filter to the binary maps and merged the resultant maps per region to produce the regional map. Then we calculated the overall, producerâs, and userâs accuracies of the final map using the corresponding geo-referenced in-situ validation dataset. The overall accuracy expressed the ratio of the number of correctly classified pixels to the total number of given pixels regardless of the class. Producerâs accuracy expressed how often P. oceanica could be identified correctly on the classified map; whereas, the userâs accuracy expressed how often the P. oceanica class on the map would be present on the ground. Producerâs accuracy is the map accuracy from the mapmakerâs viewpoint and so it is a great statistical metric for the remote sensing scientist creating the habitat map. In contrast, the userâs accuracy is the accuracy from the map userâs viewpoint. Hence, it is more significant in a management context of a given region as it reports a quantitative probability for the actual presence of the studied habitat (in this case, P. oceanica) in the given region38. The equations for all the aforementioned classification metrics are given below:
Here, TP is the number of true positives, i.e., pixels correctly identified as PO; TN is the number of true negatives, i.e., pixels correctly identified as Non-PO; FP is the number of false positives, i.e., the number of actual Non-PO pixels classified as PO; and FN is the number of false negatives, i.e., the number of actual PO pixels classified as Non-PO.
Data availability
The necessary procedures to generate the data and reproduce the methodology have been outlined in the manuscript. The data that support the findings of this study are available from the authors on reasonable request.
References
McKenzie, L. J. et al. The global distribution of seagrass meadows. Environ. Res. Lett. 15, 074041 (2020).
Duarte, C. M. et al. Seagrass community metabolism: Assessing the carbon sink capacity of seagrass meadows. Glob. Biogeochem. Cycles 24, (2010).
Short, F. T. World Atlas of Seagrasses (University of California Press, 2003).
de los Santos, C. B. et al. Recent trend reversal for declining European seagrass meadows. Nat. Commun. 10, 3356 (2019).
Waycott, M. et al. Accelerating loss of seagrasses across the globe threatens coastal ecosystems. Proc. Natl. Acad. Sci. 106, 12377â12381 (2009).
Marbà , N., DÃaz-Almela, E. & Duarte, C. M. Mediterranean seagrass (Posidonia oceanica) loss between 1842 and 2009. Biol. Conserv. 176, 183â190 (2014).
Ibiza News. Posidonia Will Disappear From Ibiza And Formentera If Sea Temperatures Continue To Rise. Ibiza News https://www.diariodeibiza.com/ibiza/posidonia-sea-temperatures-rise/ (2023).
MPN. Mediterranean Posidonia Network. Mediterranean Posidonia Network https://medposidonianetwork.com/ (2023).
Bethoux, J. & Copin Montegut, G. Biological fixation of atmospheric nitrogen in the Mediterranean-sea. Limnol. Oceanogr. 31, 1353â1358 (1986).
MedWet. Posidonia, the lung of the Mediterranean | MedWet. https://medwet.org/2017/10/mediterranean-posidonia/ (2017).
Calvo, S., Ciraolo, G. & Loggia, G. L. Monitoring Posidonia oceanica meadows in a Mediterranean coastal lagoon (Stagnone, Italy) by means of neural network and ISODATA classification methods. Int. J. Remote Sens. 24, 2703â2716 (2003).
Francour, P. Fish assemblages of Posidonia oceanica beds at Port-Cros (France, NW Mediterranean): Assessment of composition and long-term fluctuations by visual census. Mar. Ecol. 18, 157â173 (1997).
Hemminga, M. A. & Duarte, C. M. Seagrass Ecology (Cambridge University Press, 2000).
Gacia, E., Granata, T. C. & Duarte, C. M. An approach to measurement of particle flux and sediment retention within seagrass (Posidonia oceanica) meadows. Aquat. Bot. 65, 255 (1999).
Fourqurean, J. W. et al. Seagrass ecosystems as a globally significant carbon stock. Nat. Geosci. 5, 505â509 (2012).
IUCN. Manual for the Creation of Blue Carbon Projects in Europe and the Mediterranean. (2021).
Duarte, C. M., Losada, I. J., Hendriks, I. E., Mazarrasa, I. & Marbà , N. The role of coastal plant communities for climate change mitigation and adaptation. Nat. Clim. Change 3, 961â968 (2013).
Arias-Ortiz, A. et al. A marine heatwave drives massive losses from the worldâs largest seagrass carbon stocks. Nat. Clim. Change 8, 338â344 (2018).
Jordà , G., Marbà , N. & Duarte, C. M. Mediterranean seagrass vulnerable to regional climate warming. Nat. Clim. Change 2, 821â824 (2012).
Wernberg, T. et al. An extreme climatic event alters marine ecosystem structure in a global biodiversity hotspot. Nat. Clim. Change 3, 78â82 (2013).
EEA. Grouped threats to Posidonia oceanica beds as reported by Mediterranean EU Member States under the EU Habitats DirectiveâEuropean Environment Agency. https://www.eea.europa.eu/data-and-maps/figures/grouped-threats-to-posidonia-oceanica (2015).
Marbà , N. et al. Growth and population dynamics of Posidonia oceanica on the Spanish Mediterranean coast: Elucidating seagrass decline. Mar. Ecol. Prog. Ser. 137, 203â213 (1996).
Mateo, M. A., Romeo, J., Perez, M., Littler, M. M. & Littler, D. S. Dynamics of millenary organic deposits resulting from the growth of the Mediterranean seagrass Posidonia oceanica. (1997).
Lee, K.-S., Park, S. & Kim, Y. Effects of irradiance, temperature, and nutrients on growth dynamics of seagrasses: A review. J. Exp. Mar. Biol. Ecol. 350, 144â175 (2007).
Collado-Vides, L., Caccia, V. G., Boyer, J. N. & Fourqurean, J. W. Tropical seagrass-associated macroalgae distributions and trends relative to water quality. Estuar. Coast. Shelf Sci. 73, 680â694 (2007).
Duarte, C. M. The future of seagrass meadows. Environ. Conserv. 29, 192â206 (2002).
Short, F., Carruthers, T., Dennison, W. & Waycott, M. Global seagrass distribution and diversity: A bioregional model. J. Exp. Mar. Biol. Ecol. 350, 3â20 (2007).
Kenworthy, W., Wyllie-Echeverria, S., Coles, R., Pergent, G. & Pergent-Martini, C. Seagrass conservation biology: An interdisciplinary science for protection of the seagrass biome. Seagrasses Biol. Ecol. Conserv. (2006).
Sagawa, T. et al. Using bottom surface reflectance to map coastal marine areas: a new application method for Lyzengaâs model. Int. J. Remote Sens. 31, 3051â3064 (2010).
Borfecchia, F. et al. Landsat 8 OLI satellite data for mapping of the Posidonia oceanica and benthic habitats of coastal ecosystems. Int. J. Remote Sens. 40, 1548â1575 (2019).
Cozza, R. et al. Biomonitoring of Posidonia oceanica beds by a multiscale approach. Aquat. Bot. 156, 14â24 (2019).
Dattola, L. et al. Comparison of Sentinel-2 and Landsat-8 OLI satellite images vs. high spatial resolution images (MIVIS and WorldView-2) for mapping Posidonia oceanica meadows. In Remote Sensing of the Ocean, Sea Ice, Coastal Waters, and Large Water Regions 2018 vol. 10784 252â262 (SPIE, 2018).
Fornes, A. et al. Mapping Posidonia oceanica from IKONOS. ISPRS J. Photogramm. Remote Sens. 60, 315â322 (2006).
Matarrese, R., Acquaro, M., Morea, A., Tijani, K. & Chiaradia, M. T. Applications of remote sensing techniques for mapping Posidonia oceanica meadows. In IGARSS 2008â2008 IEEE International Geoscience and Remote Sensing Symposium vol. 4 IV-906âIVâ909 (2008).
Pasqualini, V., Pergent-Martini, C., Clabaut, P. & Pergent, G. Mapping of Posidonia oceanica using aerial photographs and side scan sonar: Application off the Island of Corsica (France). Estuar. Coast. Shelf Sci. 47, 359â367 (1998).
Topouzelis, K., Makri, D., Stoupas, N., Papakonstantinou, A. & Katsanevakis, S. Seagrass mapping in Greek territorial waters using Landsat-8 satellite images. Int. J. Appl. Earth Observ. Geoinf. 67, 98â113 (2018).
Traganos, D. et al. Spatially explicit seagrass extent mapping across the entire mediterranean. Front. Mar. Sci. 9, (2022).
Traganos, D. & Reinartz, P. Mapping Mediterranean seagrasses with Sentinel-2 imagery. Mar. Pollut. Bull. 134, 197â209 (2018).
Traganos, D. & Reinartz, P. Machine learning-based retrieval of benthic reflectance and Posidonia oceanica seagrass extent using a semi-analytical inversion of Sentinel-2 satellite data. Int. J. Remote Sens. 39, 9428â9452 (2018).
Poursanidis, D., Traganos, D., Reinartz, P. & Chrysoulakis, N. On the use of Sentinel-2 for coastal habitat mapping and satellite-derived bathymetry estimation using downscaled coastal aerosol band. Int. J. Appl. Earth Observ. Geoinf. 80, 58â70 (2019).
Colantoni, P., Gallignani, P., Fresi, E. & Cinelli, F. Patterns of Posidonia oceanica (L.) Delile Beds around the Island of Ischia (Gulf of Naples) and in Adjacent Waters. Mar. Ecol. 3, 53â74 (1982).
Meinesz, A. & Laurent, R. Cartographie et état de la limite inferieure de lâherbier de Posidonia oceanica dans les Alpes-maritimes (France)âCampagne Poseïdon 1976. Bot. Mar. 21, 513â526 (1978).
Boudouresque, C.-F., Blanfuné, A., Pergent, G. & Thibaut, T. Restoration of seagrass meadows in the Mediterranean Sea: A critical review of effectiveness and ethical issues. Water 13, 1034 (2021).
Rende, S. F. et al. Ultra-high-resolution mapping of Posidonia oceanica (L.) Delile Meadows through acoustic, optical data and object-based image classification. J. Mar. Sci. Eng. 8, 647 (2020).
Hossain, M., Bujang, J. S., Zakaria, M. & Hashim, M. The application of remote sensing to seagrass ecosystems: An overview and future research prospects. Int. J. Remote Sens. 36, 61â113 (2014).
Atlas submarÃ. Atlas submarÃ: cartografÃa de Posidonia oceanica a les Illes Balear. (2020).
Dewey, S. & MacMillan, A. Environmental Survey of the National Spoil Ground off Xgħajra, Malta. A Report to Ecoserv Ltd. by Seastar Survey Ltd for a Study Commissioned by the Environment and Resources Authority. (2021).
ERA. Marine environmental monitoring: towards effective management of Maltaâs marine waters. ERA https://era.org.mt/projects/marine-environmental-monitoring-towards-effective-management-of-maltas-marine-waters-2/ (2017).
Vaquer-Sunyer, R., Barrientos, N., Conselleria De medi ambient Del govern De les illes balears & Marbà , N. Posidonia oceanica. (2021).
Vanhellmont, Q. & Ruddick, K. Acolite for Sentinel-2: Aquatic Applications of MSI Imagery. (2016).
Vanhellemont, Q. & Ruddick, K. Atmospheric correction of metre-scale optical satellite data for inland and coastal water applications. Remote Sens. Environ. 216, 586â597 (2018).
Acknowledgements
This study is a part of SIMBAD project (QSR-ESABIC-2018-001) incubated by the European Space Agency Business Incubation Centre (ESA-BIC) Madrid Region, which facilitates incentives for business projects and start-ups to use space technologies or to develop applications based on these technologies with the aim of developing new products and services unrelated to space. The authors are thankful to the support from the Parque CientÃfico de Madrid where Quasar Science Resources, S.L. is incubated. We would like to acknowledge the Ministerio para la Transición Ecológica y el Reto Demográfico (MITECO), in particular the Dirección General de Biodiversidad, Bosques y Desertificación for their financial support, as well as in the provision of necessary information for the development of the project, including very fruitful discussions which have push our developments. The first author (MC) is supported by the Industrial Doctorate Program of the Spanish Ministerio de Ciencia e Innovación (ref. DIN2020-010979/AEI/https://doi.org/10.13039/501100011033). This work is part of the first authorâs PhD within the SIMBAD project. The authors would like to thank the rest of the SIMBAD team. The authors are thankful to the European Space Agency, the European Commission, and the Copernicus programme for distributing Sentinel-2 imagery. This research would not be possible without the in-situ data kindly provided by ERA and Govierno Balear. The in-situ data for the Maltese Islands was collected within the scope of the Marine Environmental Monitoring: Towards Effective Management of Maltaâs Marine Waters project (ERA 2017). Details on the survey performed can be found in Dewey and MacMillan (2021). We would like to thank Marcial Bardolet Richter (Técnico de IBANAT. ConsellerÃa de Medi Ambient i Territori), from D.G. Biodiversidad y Espacios Naturales, ConsellerÃa de Medi Ambient i Territori de les Illas Baleares, Govern Balear, for providing the in-situ data for the Balearic Islands region.
Author information
Authors and Affiliations
Contributions
M.C., A.M., M.M., N.B., I.C., M.R. and I.d.C. contributed with ideas to design the research objectives. A.M. led the design of the architecture of the scientific modules. M.C. and A.M. developed the scientific modules for data pre-processing, the DLNN model and its accuracy assessment, whereas N.B. developed the filters and time-series analysis, with support from I.d.C. E.P. and J.B. contributed to the SEPâs development. M.M. conducted all the model runs within the Docker environment for the operational approach, with support from E.P. M.C. wrote the original manuscript with input from M.M. All authors contributed to the interpretation and discussion of results and the final manuscript. I.d.C and J.B. provided part of the financial support leading to this publication.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
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/.
About this article
Cite this article
Chowdhury, M., MartÃnez-Sansigre, A., Mole, M. et al. AI-driven remote sensing enhances Mediterranean seagrass monitoring and conservation to combat climate change and anthropogenic impacts. Sci Rep 14, 8360 (2024). https://doi.org/10.1038/s41598-024-59091-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-59091-7
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.