Baseline Assessment of Ecological Quality Index (EQI) of the Marine Coastal Habitats of Tonga Archipelago: Application for Management of Remote Regions in the Paciﬁc

: The loss of coral habitats and associated biodiversity have direct effects both on the physical dynamics of the coast and on natural resources, threatening the survival of local populations. Conservative actions, such as the creation of new Marine Protected Areas, are urgent measures needed to face climate change. Managers need fast and simple methods to evaluate marine habitats for planning conservation areas. Here, we present the application of an Ecological Quality Index (EQI), developed for regional-scale habitat maps of the Atlas of the Marine Coastal Habitats of the Kingdom of Tonga, by processing Copernicus Sentinel-2 imagery. Both the habitat mapping classification and the EQI application were focused on the importance of coral reef, seagrass and mangrove habitats, both as natural defense and sustenance for the local populations. Twelve main Pacific reef habitats were evaluated through a three-level EQI score assigned to six parameters: nursery ground, connectivity, species reservoir, fish attraction, biodiversity and primary production. The EQI was integrated into a developed georeferenced database associated to the QGIS software providing the ability to identify on the maps the area of interest and the associated habitats, and to quantify their ecological relevance. The EQI is proposed as a tool that can offer to stakeholders and environmental managers a simple and direct indicator of the value of the marine coastal environment. The index may be handled for management purposes of vast areas with remote and uninhabited islands.


Introduction
Marine habitat mapping is a fundamental tool used for the knowledge and management of marine coastal environments. Among the regional-scale approaches, habitat maps derived from satellite images provide the coverage of large areas in remote and uninhabited regions. This is extremely useful for the control and monitoring of vast archipelagos with scattered islands as in the Pacific Ocean. Here, governments face natural disasters, such as tsunamis, volcanic eruptions, hurricanes and floods, which have increased their intensity in recent years owing to the effects of climate change. Natural disasters may also be due to coral bleaching or to the outbreaks of alien or local species, such as the coral-eating crown-of-thorns starfish (COTS) Acanthaster sp., which has contributed to a coral cover decrease [1].
In tropical marine environments, remote sensing reaches its greatest effectiveness in the submerged environment due to the transparency of the waters. Satellite images show the extent of the coastal communities to be recognized and quantified, in some cases up to a depth of 30 m. Among the most important marine communities identified with remote sensing, particular attention is paid to coral reefs, seagrass meadows and mangrove forests. These marine communities are fundamental not only for maintaining the balance of the coasts [2] and their natural defense from storms or flooding events [3] but also for their biodiversity, which guarantees the survival of local populations whose most important resource is fishing [4]. Fishing activity, which is vital for people living on the few islands of vast archipelagos, may be locally endangered due to the illegal exploitation or excessive collection of species, such as corals, sea cucumbers, fish and mollusks [5,6]. To counter these damaging activities, Pacific communities have launched multi-year projects to safeguard their resources and encourage the development of marine spatial plans (MSPs). MSPs have, among their priorities, the establishment of marine-protected areas and the sensitization of local populations by the involvement of people in activities, such as the replanting of corals or mangroves, and the promotion of tourism activities of low impact, such as snorkeling, whale watching and diving [7][8][9][10][11].
To support these development projects, remote sensing and marine habitat mapping are extremely effective tools for providing managers with information on the entire territory of the state, including the remote and uninhabited islands and coral banks that prove to be the richer in species, and hotspots for maintaining biodiversity and pristine habitats in the Pacific [12].
Marine habitat mapping evolved in this last decade when new constellations of satellites with new sensors and higher spatial resolution made it possible to increase the collections of images begun in the 1970s with the launch of satellites dedicated to marine investigation, such as the Landsat [13]. Today, Sentinel 1, 2 and 3 have provided access to free repositories of images and dedicated software, increasing enormously the studies on marine coastal habitats. Moreover, the joint action of international bodies and private foundations have contributed to the publications of digital coral reef atlases. The Atlas of the 'Millennium Coral Reef Mapping Project' conducted in [2003][2004] was the first attempt to provide a large-scale mapping of the world's coral reefs using Landsat images [14]. The Millennium Project was designed to assess the surface of modern reefs worldwide; to do so, a classification scheme based on 126 hierarchical classes, mainly geomorphologic classes, was developed [15].
In 2013, the Khaled bin Sultan Living Oceans Foundation (KSLOF) conducted the Global Reef Expedition (GRE) to evaluate the health status of 13 coral reef locations in the world. High-resolution bathymetric maps and thematic habitat maps collecting 16 marine habitat classes were created [16,17] from WorldView-2 satellite images.
In 2019, the Allen Coral Atlas Project was published online [18]. The classification was based on coral reefs in the United Nations Environment Program (UNEP) 2010 "Coral Layer" integrated with Planet Scope satellite imagery. The Atlas included 22 geomorphic zones and six benthic habitats from the coastline up to a depth of 10 m [19,20].
Habitat maps are fundamental to know the extension of habitats in each region of interest and provide useful information to meet the global demand to face climate change and loss of biodiversity. Due to the complexity of surveying large and often inaccessible marine areas with traditional approaches such as biological censuses, habitats have been proposed as surrogates to describe the biodiversity of reef species with special attention to fish and invertebrate of commercial value [21][22][23].
Habitat maps give to stakeholders and decision makers the ability to view the extent of habitats on a map. However, for marine spatial planning, dedicated maps on the value/health status of an area should be planned. The multi-metric, indicator-based assessment approach uses parameters collected locally from water, sediment and biota that are illustrated on the maps through a score/color classification. This approach, largely adopted in Europe, is used to classify marine areas from local to regional scales [24][25][26] and to support ecosystem-based management through the existing international networks [27]. In Pacific areas, large-scale habitat mapping and multi-parametric local data are rare. To describe Tonga's special unique marine areas a qualitative approach was based on four criteria of classification: geographic explicitness, justification, information sources and legal obligations [8].
Here, we present the application of an Ecological Quality Index (EQI) to the Atlas of the Marine Coastal Habitat Mapping of the Kingdom of Tonga derived from Copernicus Sentinel-2 images. The habitat mapping allowed the quantification of the cover of each ecological habitat for the entire archipelago of Tonga (approximately 174 islands).
The EQI was evaluated for twelve main Pacific reef habitats through six ecological parameters, which were focused on the potential of reef habitats as refugia, resilience and recovery sites for reef species after a disturbance that alters the physical and biological structure of the reef. The EQI is proposed as a tool that provides to stakeholders and environmental managers a simple and direct indicator of the value of the marine environment to be handled for management purposes of vast areas with remote and uninhabited islands.

The Atlas of the Marine Coastal Habitats of the Kingdom of Tonga
The Kingdom of Tonga is located in the central South Pacific, and its archipelago spans 950 km including the Minerva Reefs. The archipelago consists of approximately 174 islands scattered across 360,000 km 2 with only 720 km 2 above sea level and only 37 inhabited islands [5,6]. The project of the Atlas of the Marine Coastal Habitat of the Kingdom of Tonga was born with the objective to use remote sensing as a basic tool for monitoring such a large area. To start a monitoring program, it is fundamental to fix a unique reference time to compare data from a chronological point of view; to meet this requirement, free satellite images available from January 2018 to December 2019 for the whole Tonga Archipelago were collected from the Sentinel-2 and Landsat free repositories on the web portals of the U.S. Geological Survey USGS-Earth Explorer [28] and the ESA-Copernicus organization [29]. Image processing followed the procedure of Immordino et al. [30]; the pre-processing steps included image radiometric, atmospheric and geometric correction, application of depth invariant bands and correction for the air-water interface through the use of ERDAS Imagine software. Two separate approaches were used for the habitat classification: the unsupervised classification was used with ENVI ® software v. 5.6.3 to classify images using bands 1-2-3-4-5-8 with bands 1 and 8 resampled from the native spatial resolution of 20 m to 10 m, and the supervised classification was used with Trimble eCognition ® software v. 10.1 to map the littoral mangrove forests. The unsupervised classification performs a subdivision of the feature space into a certain number of "spectral classes" without requiring the informative contribution of the number and identity of the classes. The fundamental hypothesis of this approach is that similar objects must correspond to similar spectral responses, and this is reflected in the space of the characteristics with values associated with the pixels of the image, which tend to be distributed in groups or clusters; the unsupervised method, therefore, finds a certain number of clusters without knowledge of the number or nature of the clusters. The criteria used for subdividing the multispectral space into "spectral classes" are proximity or similarity criteria. After the automatic assignment to the different spectral classes, the interpretation is up to the expert who, with the experience and knowledge of the characteristics of the study environment, identifies the thematic classes to be associated with the spectral classes identified with the clustering procedure [31]. Fourteen classes related to coral coastal marine habitats (terrestrial mangrove habitat excluded) have, therefore, been defined in the ISODATA algorithm of the unsupervised classification; the raster file obtained has been transformed into .shp vector format. Hence, the ENVI ® software allowed us to get a first classification of habitat classes exported in vector format as shapefiles (.shp) to be imported into QGIS freeware software for photo interpretation.
The mangrove habitat was mapped following, partially, the approach applied by Wang et al., 2018 [32]. A two-level hierarchical structure based on the spatial structure of a mangrove ecosystem and geographic object-based image analysis was utilized and modified. So, the mangrove determination in the Archipelago of Tonga was divided into two levels, namely landscape and vegetation cover type (Table 1). The object-based image analysis approach was used to implement the hierarchy of Levels 1 and 2. For this purpose, Trimble eCognition ® software was used. Creating an object-based classification mainly consists of two steps: image segmentation and object classification. Table 2 shows the algorithm (rule set) and classification processes developed for the archipelago of the Kingdom of Tonga in which some spectral indices are selected and combined.

Level 1 Water
The first level was used to discriminate vegetation and non-vegetation and produced a mask of vegetation. A chessboard segmentation with an object size of one pixel was applied in order to preserve the pixel value of coarse or moderate resolution images. Following segmentation, the objects with a Modified Normalized Difference Water Index MNDWI (1) higher than zero and a Forest Discrimination Index FDI (2) less than zero were classified as water. Moreover, due to some mudflat overlapping with water and the brightness of the mudflat being greater than water, the brightness was used as a supplement to exclude these mudflats (Brightness < 1250). The second level was used to separate mangrove from non-mangrove within the vegetation mask, and Wang et al. [32] applied the Mangrove Discrimination Index 2 (MDI2) index: By applying the MDI2 index in the few islands of the Tonga archipelago where the presence of mangroves is confirmed with the bibliography, it has been observed that, in this context, sometimes the use of the MDI2 index overestimates the extension of the mangroves. For this reason, a check with a selection of the values of Band 11 (SWIR1) was also applied together with the global basemaps on QGIS as described hereafter.
The processing allows us to obtain vector data that, in the same way as the coral unsupervised classification, was reprocessed with QGIS (version 3.4 "Madeira") through photo-interpretation tools. In this context, we used QGIS, which is an open-source, crossplatform desktop, Geographic Information System (GIS) application that supports viewing, editing and analysis of geospatial data [33]. In QGIS, all habitat classes obtained through supervised and unsupervised classification were carefully checked and corrected through QGIS tools, such as the HCMGIS plugin, that offers global basemaps from Google, Carto, ESRI and OSM Stamen that may be added as backview of the habitat mapping. Digital nautical maps, photos and videos available on the WEB were used as ancillary data to help in defining the typology and the limits of the habitat classes.
The final cartography and habitat classification formed the digital Atlas of the Marine Coastal Habitats of the Kingdom of Tonga in a shapefile (.shp) format viewable on all GIS software. Shapefile format is a digital vector storage format that stores the polygons related to each habitat class (layer). Data attributes (i.e., geographical location, cover in hectares and EQI) are linked to each layer. The digital Atlas is the representation of the shapes on a geographic georeferenced plane.

The Habitat Classification and the Ecological Quality Index (EQI)
With the aim to identify priority habitats useful for further analysis on the environmental health status, we referred mainly to the NOAA Shallow-Water Benthic Habitats Manual [34] and to the Millennium Project [14], reporting detailed illustrations and photographs of the geomorphology of the coral reef zones and of the associated biological communities of the Pacific and their ecological connections. Ancillary information, such as scientific reports, literature, web-free imagery, videos and photos, was used for interpretation and to achieve the maps of relevant habitats without field data validation [30]. To compute the Ecological Quality Index (EQI), the benthic classification and habitat mapping was adapted to the objectives of the Atlas with particular emphasis to corals, seagrasses and mangroves as priority habitats for the maintenance of biodiversity and ecosystem services [22,23].
In the two main morphological zones of the reef islands, the barrier reef and the coastal reef, the following 12 main biological habitat classes and 5 subzones were identified ( Figure 1 and Table 3):

1.
Sparse Coral/Coral Beds (SC): This class groups all colonies or beds of colonies sparse on sand or gravel bottoms around islands or in the lagoons. SC has a great importance as a fish refuge and area with coral recolonization potential after mortality events.

2.
Reef Banks (RB): These banks are found offshore. In the classification procedure, RB were considered all reefs isolated from barrier reefs or islands from a depth greater than 30 m. They can emerge at the surface with reef crest (RC) and go deeper behind the depth limit of the satellite detection. RB may be considered as submerged, coralline islands or coralline shoals. In some cases, RB may reach hundreds of meters of depth with active exchanges of larvae and nutrients from the bottom to the surface. For this reason, offshore RB may be considered as pristine areas with high biodiversity and high fish abundance that act as true marine biosphere reserves.

5.
Reef Crest (RC): It is the edge of the Barrier Reef. It is the part of the Barrier Reef more exposed to open-ocean waves. It is the most important area of the barrier for the defense of the coastline, and it is most sensitive to mortality due to low tides, elevated seawater temperatures and storms. This class is also associated to the emerging part of the Patch Reef (PR) and Reef Banks (RB).  The Ecological Quality Index (EQI) was calculated for each of the 15 habitats classes/subzones of ecological relevance as the mean of the three-level scores (1 = low value/Green, 2 = medium value/Yellow, 3 = high value/Red) assigned to the following 6 ecological parameters: nursery ground, connectivity, species reservoir, fish attraction, biodiversity and primary production ( Table 2). The ecological parameters were focused on the value of each habitat class and subclass as refugia, resilience and recovery sites for reef species (Table 2).  Table 3. Main Classes and subzones used for Tonga habitat mapping. The Ecological Quality Index (EQI) was calculated as the mean of the three-level scores assigned to six ecological parameters (nursery ground, connectivity, species reservoir, fish attraction, biodiversity, primary production).

Biologic
Biological Value (1= Low, 2 = Medium, 3 = High) Ecologi cal Nurse  Table 3. Main Classes and subzones used for Tonga habitat mapping. The Ecological Quality Index (EQI) was calculated as the mean of the three-level scores assigned to six ecological parameters (nursery ground, connectivity, species reservoir, fish attraction, biodiversity, primary production).  An additional four classes with no or scarce biological relevance were identified during the habitat classification and were included in the Atlas: Sand (S), Littoral Sand or beaches (LS), Beach Rocks (BR) and No Reefs (NR). The classification of these environments was necessary for their relationships with biological classes and for their natural value for the Marine Protected Areas. For example, sand cover in the lagoons may vary in seasons and years, and seagrasses may colonize new sand bars or disappear under them. Littoral Sand (LS) and Beach Rocks are of relevant interest from a geomorphological point of view and may be of interest both for coastal management and related tourist activities. The No Reef (NR) class is associated with rocky, basaltic reefs of volcanic island, which are boulder beaches with no value from the biological point of view.

Biological Habitat Classes
The Ecological Quality Index (EQI) was calculated for each of the 15 habitats classes/subzones of ecological relevance as the mean of the three-level scores (1 = low value/Green, 2 = medium value/Yellow, 3 = high value/Red) assigned to the following 6 ecological parameters: nursery ground, connectivity, species reservoir, fish attraction, biodiversity and primary production ( Table 2). The ecological parameters were focused on the value of each habitat class and subclass as refugia, resilience and recovery sites for reef species (Table 2).

The Atlas of the Marine Coastal Habitats of the Kingdom of Tonga
The data available for the period from 1 January 2018 to 1 December 2019 showed the greatest number of satellite images available in the Sentinel-2 repositories. Since cloud cover is observed frequently on islands in the tropics, the satellite image analysis was restricted to images with a cloud cover <10%. On the Sentinel-2 images, a further analysis was conducted on every image covering hundreds of km 2 to ensure the best, most recent and free of clouds view for each island or group of islands. Although the research focused on the period 2018-2019, for some areas, such as Hihifo-Tafai and Tofua-Kau, the only image without clouds dated back to 2017. Moreover, only one image was found for one of the two atolls of Minerva Reef. For Ata Island, the one Landsat-8 image free of clouds was utilized. With regard to mangroves, only one image dated back to 2016 for Tongatapu and one recent 2020 image for Vavau were free of clouds and were utilized. The results of the habitat classification on QGIS allowed us to visualize the habitat mapping for the Tonga Archipelago, compute the cover in hectares of the mapped habitats and estimate their ecological importance through the EQI (Figure 2).
The QGIS mapping data allowed the publication of the Atlas of the Marine Habitat of the Kingdom of Tonga [35], reporting 25 tables of the Tongan Archipelago. To obtain the Atlas, the DataPlotly plugin was used; the plugin allows the creation of interactive charts with data in shapefiles [33]. New fields were inserted into the database with color-coding of the Habitats and Ecological Quality Index that allowed a representation of tables with related pie charts. The 25 sheets composing the Atlas of the Tongan Archipelago are all the same size, 32 × 28 km, and have been included in an A3-size layout with a scale of 1:100,000. (Figure 3a).
image without clouds dated back to 2017. Moreover, only one image was found for one of the two atolls of Minerva Reef. For Ata Island, the one Landsat-8 image free of clouds was utilized. With regard to mangroves, only one image dated back to 2016 for Tongatapu and one recent 2020 image for Vavau were free of clouds and were utilized. The results of the habitat classification on QGIS allowed us to visualize the habitat mapping for the Tonga Archipelago, compute the cover in hectares of the mapped habitats and estimate their ecological importance through the EQI (Figure 2).

The Habitat Classification and the Ecological Quality Index (EQI)
A total of 104,147 hectares of the Tonga Archipelago were investigated of which 82,526 hectares (79.23%) were classified as relevant from a biological point of view. The QGIS database allowed us to compare the five main Tongan administrative regions, Ata Island and Minerva atoll (Figure 3b) on the occurrence of the 15 main biological habitats classes/subzones. The comparison showed that the Fore Reef (FR), Fringing Reef (FRI) and Reef Flat (RF) habitats were found in all of the archipelago with the highest cover percentages (Figure 4). The Sparse Coral (SC) habitat was found in Vavau, Ha'apai and Tongatapu on littoral, sandy bottoms. The Back Reef (BR) was recognized in all regions except for Eua. The Mangroves (MA) and the associated Mud (MU) habitats were limited to Vavau and Tongatapu, the regions with the greatest islands. Seagrasses (SE) were found in all regions with cover in hectares reaching the highest values in Tongatapu. Reef Banks (RB) were found only in the Vavau and Ha'apai regions.
The Ecological Quality Index per region showed that all marine areas of the Archipelago have high ecological values from 39% upward ( Figure 5). The assessment of habitats of biological importance for the whole Tonga Archipelago showed that Fore Reef, Sparse Corals and Reef Banks are the habitats with the highest sea bottom cover with values from 10,791 to 14,717 hectares, contributing to the total highest values of the EQI assigned to 71.42% of the total area covered by habitat mapping in this study. The high cover of Fore Reef and Reef Banks is related to the Tonga Ridge, which marks the eastern boundary of the large plateau connecting the Vavau region to the Ha'hapai region and the Ha'hapai region to the Tongatapu region [36]. The Sparse Corals high cover is related to large areas of shallow bottoms that are typical of all Tonga regions.   Table 2.
The Ecological Quality Index per region showed that all marine areas of the Archipelago have high ecological values from 39% upward ( Figure 5). The assessment of habitats of biological importance for the whole Tonga Archipelago showed that Fore Reef, Sparse Corals and Reef Banks are the habitats with the highest sea bottom cover with values from 10,791 to 14,717 hectares, contributing to the total highest values of the EQI assigned to 71.42% of the total area covered by habitat mapping in this study The high cover of Fore Reef and Reef Banks is related to the Tonga Ridge, which marks the eastern boundary of the large plateau connecting the Vavau region to the Ha'hapai region and the Ha'hapai region to the Tongatapu region [36]. The Sparse Corals high cover is related to large areas of shallow bottoms that are typical of all Tonga regions.

Discussion
Despite that the Sentinel-2 images have a resolution limited to 10 m for Blue, Green, Red and Near-Infrared bands, the high revisiting time of 5-10 days allows us to visualize the environmental changes due to storms or tsunamis, or to detect coral mortality events  The Ecological Quality Index per region showed that all marine areas Archipelago have high ecological values from 39% upward ( Figure 5). The assessm habitats of biological importance for the whole Tonga Archipelago showed that For Sparse Corals and Reef Banks are the habitats with the highest sea bottom cove values from 10,791 to 14,717 hectares, contributing to the total highest values of th assigned to 71.42% of the total area covered by habitat mapping in this study Th cover of Fore Reef and Reef Banks is related to the Tonga Ridge, which marks the e boundary of the large plateau connecting the Vavau region to the Ha'hapai regio the Ha'hapai region to the Tongatapu region [36]. The Sparse Corals high cover is r to large areas of shallow bottoms that are typical of all Tonga regions.

Discussion
Despite that the Sentinel-2 images have a resolution limited to 10 m for Blue, Red and Near-Infrared bands, the high revisiting time of 5-10 days allows us to vis the environmental changes due to storms or tsunamis, or to detect coral mortality due to climate change, also in remote and uninhabited areas. Hedley et al. [13] evid that the highly usable acquisition rate of Sentinel-2 makes it possible to collect s good images a year, meeting the needs for several essential coral reef scientif monitoring objectives. For example, the detection of bleaching events requires frequency of images. The revisiting time of Sentinel-2 was demonstrated to be ad

Discussion
Despite that the Sentinel-2 images have a resolution limited to 10 m for Blue, Green, Red and Near-Infrared bands, the high revisiting time of 5-10 days allows us to visualize the environmental changes due to storms or tsunamis, or to detect coral mortality events due to climate change, also in remote and uninhabited areas. Hedley et al. [13] evidenced that the highly usable acquisition rate of Sentinel-2 makes it possible to collect several good images a year, meeting the needs for several essential coral reef scientific and monitoring objectives. For example, the detection of bleaching events requires a high frequency of images. The revisiting time of Sentinel-2 was demonstrated to be adequate to follow an event that may last for 1 to 6 weeks before coral recovery or colonization by macroalgae [13]. However, to achieve detailed high-accuracy data to follow mortality events to adequate scales, an effort should be made to continue to improve the input data. To reach a monitoring of coral reefs at regional to global scales, point-based in situ observations on biodiversity, validation procedures and the use of automated procedures using artificial intelligence and machine learning approaches should be implemented [19,37]. The proposed Ecological Quality Index (EQI) was based on the most important prerequisite of reef monitoring through remote sensing: the classification system should make good ecological sense [38].
The use of Sentinel-2 images to monitor large areas of the Tonga archipelago demonstrates the possibility to perform regular surveys with satellite images. Despite the difficulty to find good images because the cloud cover is a constant disturbing factor in tropical zones, the limited time period of two years imposed at the start of the project proved to be a practical way to use satellite images in large-scale monitoring programs on coral reefs. This monitoring attribute is particularly useful for extended territories, such as Pacific coral archipelagos formed by dispersed, uninhabited islands, banks or shallow bottoms with species-rich habitats that may be a fundamental source for the recovery of areas damaged by natural and/or anthropic pressure [1,22].
However, it is important to keep in mind that the information derived from the Atlas is just a piece of a network of information. The proposed approach is useful for a first mapping of large and remote areas at a regional scale. The habitat recognition and classification without ground truth was possible thanks to the existing publications of atlases and manuals about Pacific habitats and to ancillary products as free high-resolution internet images. However, the recognition was limited to main habitat classes due to the satellite resolution that limited the reading to a maximum magnification scale of 1:100,000. Hence, information must be completed and corrected with appropriate ground truthing and detailed surveys on the areas of interest [19].
Moreover, the submerged habitats of corals and seagrasses were identified taking into consideration the operational limit of the satellite images, normally within 20-30 m depth in the best conditions of water transparency. Deeper areas, e.g., offshore zones around banks, should be investigated with other methodologies, such as multibeam or sidescan sonar, or underwater imagery to obtain observations of biodiversity [39,40]. These deep areas are of relevant ecological importance; coral species inhabiting these zones could be an important reservoir for the shallower reefs impacted by climate change, and they should be protected from overfishing and included in marine spatial planning [6,36,41].
The procedure followed in the present work using open-access, free Sentinel-2 images and freeware, open-source mapping platforms is in general agreement with the seven main attributes of the mapping framework of Lyons et al. [19]. Among these attributes, the most important is to compose a mapping framework that is easy, efficient and timely to 're-run' by the original investigators or by different individuals in order to update the outputs when new sensors or new data are available. GIS tools provide a quick and efficient visualization tool of the digital maps of the Atlas of the coastal marine habitats of Tonga and of the associated EQI. This configuration makes the Atlas accessible by everyone on laptop computers, also without an internet connection, with the advantages associated with its use in field campaigns and remote areas.
Digital maps provided by the Atlas and its database are open and available to Tongan stakeholders (scientists or administrative personnel). QGIS software gives the ability to assess the status of any area in the Tonga archipelago in terms of habitat occurrence, cover and ecological value (EQI), contributing to the evaluation of ecosystem services associated to each habitat [22,23]. The Atlas offers a baseline for monitoring environmental changes, allowing the quantification of natural or anthropic damages in term of habitat loss. The Atlas could be modified, updated and implemented with new maps on marine and terrestrial habitats, offering an efficient system for spatial planning issues and to prevent climate change. Thanks to GIS tools, past and future marine habitat maps may be uploaded and chronologically compared to the marine environment status recorded for the entire archipelago of Tonga for the years 2016-2020. Moreover, the digital Atlas could be added in the future to the WEB portals as the Pacific Environment Data Portal [42], the Pacific Data HUB [43] or the Ocean data Viewer [44], collecting data and maps on coral environments at regional or global scales.
Following the ongoing Marine Spatial Planning of Tonga, the proposed marine habitat mapping and the EQI are in accordance with the main objectives of the project [36]: to help decision makers appreciate the values of marine ecosystems and the importance of spatially planning the uses of these values; and to assist the planning processes using the accompanying data layers and raw data in Geographic Information Systems (GIS).

Conclusions
Satellite images, thanks to the innovation in image and software availability and despite some problems related to areal coverage, presence of clouds or image resolution, are consolidated tools to study and quantify the effects of climate change and human impact on earth's ecosystems. The Ecological Quality Index (EQI) proposed in this work combines the relevant biological characteristics of reefs and associated ecosystems with specific ecological classes and the ecosystem values associated with them. The digital Atlas of Marine Habitats of the Tonga Kingdom running on the freeware QGIS platform is proposed as an example of a simplified tool for government technical offices for a first working approach of the marine spatial planning of large remote areas and uninhabited and pristine islands.