Dynamics of Land Use in the Basin of Ancient Lagoons (Toho, Todougba, Ahouangan, Dati, Djonou) in South Benin from 1990 to 2020 and Their Current Vulnerability to Pollution

The management of living resources and space requires continuous assess-ments in order to adjust deviations as necessary for sustainable use. It is in this context that the present study was conducted and aims on the one hand to assess the dynamics of land use in the basin of the five-finger lagoons and on the other hand to assess the degree of their vulnerability to organic pollution. Thus, the satellite images were collected in free access on the site http://earthexplorer.usgs.gov/ of the USGS (US Geological Survey) of NASA. Image correction was performed with ENVI 5.1 software and classification was performed with ArcGIS 10.8 software. The degree of vulnerability was evaluated through the weighting of different parameters (slope, land use, type of pass, runoff, and density of the hydrobiological network) by the mul-ti-criteria method of Saaty (Analytic Hierarchy Process Methods) and the application of mapping. It is found under a very high degree of vulnerability (73% - 100%) and do not augur well for the living resources. The other lagoons, notably the Dathi and Djonou lagoons, are also under stress from the pressures of surrounding activities. It is therefore important to assess the level of pollution of the Toho and Todougba lagoons through pollution indices in order to make decisions that will allow their sustainable use, especially for aquaculture.


Introduction
Wetlands are ecosystems that provide various economic, social and ecological functions directly or indirectly for humanity. They also ensure the physical-chemical and biological regulation of the physical environment, the conservation of biological diversity and constitute a very productive environment on the planet (Akotossode et al., 2018;Nfotabong-Atheull et al., 2013;Onana et al., 2014). However, population growth and its corollaries (land use, agriculture, fish farming, fisheries, need for food security and others) coupled with the effects of climate change (climate disruption and others) constitute permanent stress factors increasingly growing on wetlands (Rao et al., 2007). As a result, the current state of a wetland is highly dependent on how it has been used over the years.
In Benin, the basin of the complex of ancient lagoons still called the five-finger lagoon (including the Toho and Todougba lagoons) is the object of a growing attraction for various activities, such as fish farming in cages, agriculture, tourism (for its ornithological richness) with the development of hotels, places of rest and recreation, and crafts. These activities are likely to create modifications in the functioning of the basin in this case of the aquatic ecosystems which are the final receptacles of the whole drainage network. These modifications can lead to losses of biodiversity, fragmentation of biological habitats (Benayas et al., 2009;Dauda et al., 2014;Yoboue et al., 2018), of pollution (Degefu et al., 2011;Djihouessi et al., 2019) and other impacts that affect the biocenosis. These susceptibilities may be exacerbated by the extent of greenhouse gas production in wetlands (Deemer et al., 2016).
The Toho and Todougba lagoons which constitute the major water body of this basin (Hagemeijer et al., 1997) through the potentialities they offer to fish farming and fishing, are experiencing the beginning of organic pollution (Achoh et al., 2021). Thus, the more degradation pressure there is on the basin, the more organic matter will be exchanged between the different lagoons of the complex, which will be polluted and become unsuitable for aquatic life.
The objective of the present work is to evaluate the effect of anthropogenic pressure on the five-finger lagoon basin over the last thirty years (1990-2020) based on remote sensing and to deduce the current vulnerability to pollution of these lagoons. Journal of Geoscience and Environment Protection

Study Area
The Ancient Lagoon Basin, also known as the Five Finger Lagoon, is located in the heart of Ramsar Site 1017. It is located between 6˚23'34.89''N and 2˚10'24.97''W. It is bounded to the north by the commune of Tori-Bossito and Kpomassè, to the south by the commune of Cotonou, to the east by the commune of Abomey-Calavi and to the west by the arrondissement of Pahou (commune of Ouidah). The basin has five lagoons: Toho, Todougba, Ahouangan, Dati and Djonou. In addition to these lagoons, the basin has a vast network of swamps, artificial lakes, flood plains and vegetation consisting of palm groves (Elaeis guineensis), fallow fields, savannah and classified forest (Hagemeijer et al., 1997). Several activities are carried out in the area, including agriculture, fishing, aquaculture, livestock, market gardening, handicrafts and hotel and restaurant activities.
The Toho and Todougba lagoons constitute the major ecosystems of this basin and are therefore very attractive for cage fish farming (Aïzonou et al., 2019) and other activities, particularly tourism. The communication between the different arms of the aquatic ecosystems of the basin conditions the exchange of dissolved matter. The descriptions of the study environment are presented in

Data Collection
The images used in this study were extracted by the Global Visualization Viewer on the site https://earthexplorer.usgs.gov/ of the USGS (US Geological Survey) of NASA. They were obtained on different satellites (Landsat 5 TM, Landsat 7 TM, Landsat ETM +, Landsat OLI-TIRS) whose specifications are summarized in Table 1.

Land Use Classification
Using ArcGIS 10.8 software, the maximum likelihood classification algorithm was used to classify the four satellite images. The land cover classification was performed using the supervised classification method using ENVI 5.1 software.
Nine ( the value "1" and testify that the classification performed strongly reflects the realities of the field. The rate of change (Tc) of land cover in a given period was calculated by applying the formula used by Toyi et al. (2013).

Assessment of Vulnerability to Pollution of Lagoons
The vulnerability of a surface water to pollution depends on several parameters including slope, soil type, land use, runoff, density of the watershed network (Ake et al., 2020;Ake et al., 2019;Anoh et al., 2012Anoh et al., , 2018Macary et al., 2010).
Based on different works done, the vulnerability assessment was done in three steps namely: identification, classification and assignment of scores to the different parameters; determination of the weighting coefficient by the multi-criteria method of Saaty (1977) (Analytic Hierarchy Process Methods) and application of mapping for vulnerability assessment.

Identification, Classification and Assignment of Notes to the Different Parameters
Slope: The slope map was made after processing a map of the Digital Elevation Model (DEM) of the basin with ArcGIS 10.8. Then, the classification was made with the module "3D analysis tool" of the same software and obtained six (6) classes ranging from 0% to 6.99671% and presented in Figure 3.
Soil type: The soil map of the study area was obtained from the soil profile of Benin. The area has three main soil categories (Figure 4), namely -Highly depleted ferralitic soils -Ferralitic soils that are strongly rejuvenated -Ferruginous soils Annual runoff: Runoff was determined using the Thornthwaite method. It was determined by multiplying the runoff coefficient by the average annual rainfall for the last ten years. Let R be the annual runoff (mm), α runoff coefficient (%) and Pa the annual rainfall (mm).   The runoff coefficient for Benin, the former Dahomey is 69% (Roose, 1973;Roose & Lelong, 1976). The average annual rainfall obtained from the Water Directorate is 1200 mm·an −1 . Thus, the calculated annual runoff is: R = 828 mm. This factor has only one class with a rating of 3.
Density of the hydrographic network: The density of drainage was obtained by making the map from the "line-density" tool of the ArcGis software. The density varies between −6 and 44 km. The mapping of the parameters was followed by a reclassification of the different values of these five parameters based on the work of Schoen and Codvelle (2001) and allowed to represent the classification maps of the said parameters from the tool "Reclassify" of the module "Spatial Analyst" under ArcGis before their weighting. The classification resulted in six (6) Figure 5.

Scoring of the Different Parameters Classes
The assignment of scores to the different classes of parameters was inspired by the work of Schoen & Codvelle (2001) on the vulnerability of surface waters. Thus, the scores assigned vary from 1 to 4, taking into account the importance of the different parameters of vulnerability to pollution of the ecosystems studied. Table 2 shows the scores assigned to the different classes of vulnerability parameters.

Weighting of the Different Classes of Parameters by Saaty's Multicriteria Method
The multicriteria method AHP (Analytic Hierarchy Process) of (Saaty, 1977) was used for the weighting of the parameters and is summarized in two steps, namely the realization of the binary combinations and the development of the priority matrices.
Realization of the binary combinations: For the binary combination, the parameters were compared two by two according to their level of importance. These comparisons lead to the constitution of reciprocal square matrices. According to this comparison, when two parameters have the same importance, the method (Saaty, 1977) requires that the value of "1" be assigned to both parameters. When one parameter is more important than the other, it takes a higher value between 1 and 10; the other takes the inverse of this value. This method produces standardized weights that sum to "1". The matrix resulting from the pairwise comparison of the different parameters is summarized in Table 3.
Determination of the weighting coefficients: The weighting coefficient (WC) of a parameter corresponds to its level of involvement in the vulnerability to pollution of the aquatic ecosystem studied. It is determined in two steps: the determination of the Eigenvector noted (Vp) and the calculation of the weighting coefficient (Cp).
The formula for determining the eigenvector is: With n = the total number of parameters compared and a 1 , a 2 , a 3 , …, a n the different values attributed to the parameters during the comparison.
As for the weighting coefficient (Cp), it is calculated by the following formula: With RI the random index and CI the consistency index.
The values of the random index are a function of the number of vulnerability parameters compared (Saaty, 1977). Considering the number of parameters compared in this study which is 5, the value of RI is 1.12 (Saaty, 1977).   The CR Consistency ratio is less than 10% (0.23%), the assigned judgments are good to better express the vulnerability of the basin to pollution.

Vulnerability Assessment
The vulnerability index (Iv) is obtained from the following formula: From the vulnerability index, the classification map was made and obtained six (06) classes between 0.23 and 21.80. The conversion of the classes obtained into percentages was done using the method used by Jourda et al. (2007). This conversion allows a better understanding of the degree of vulnerability. The With Ii = index to be identified, Imax = the maximum index (21.80) and Imin = the minimum index (0.23)

Results
The diachronic study of land use in the basin of ancient lagoons (Toho, Todougba, Ahouangan, Dati, Djonou) carried out over the last thirty years (1990-2020) revealed significant variations in time and space ( Figure 6). Indeed, nine (09) classes were defined from the classification carried out by 10-year intervals. This classification proved to be very significant according to the Kappa index, which is of the order of 0.9 for all the treatments carried out, testifying to the convergence between the field data and the satellite data.
From the analysis of the data obtained (Table 5 and Table 6), it appears that over the last thirty years , the area of dense forests, open forests, savannahs and plantations have significantly decreased by 98.93%, 100%, 99.60% and 77.10% respectively. On the other hand, the other classes, notably bare soil, settlements, fields and fallow land, swamps and water bodies have recorded a significant increase in their surface area. Thus, the soils are completely bare (Tc = 1264%), the settlements are more erected (Tc = 364%). The water bodies and swamps have extended their bed further with respective rates of change of 23.33% and 54.18%.
From the above, we can see that the surface area of water bodies and swamps has increased over the last thirty years, as well as that of bare land and settlements, while the surface area of forests, savannahs and plantations has decreased.
Index and degree of vulnerability of the lagoons of the basin to pollution The different classes obtained from the vulnerability index map vary from [0.23 -4.89] to [18.50 -21.80] with 0.23 and 21.80 being the minimum and maximum values respectively. From the analysis of the vulnerability index map (Figure 7) and Table 7, it appears that, with the exception of the coastal zone, high degrees of vulnerability (73% to 100%) are observed in the first two arms of the complex, notably the Toho and Todougba lagoons, and in a few surrounding pockets. The Ahouangan lagoon, which is the third arm, also recorded pockets of high vulnerability. For the Dathi and Djonou lagoons, high degrees of vulnerability are observed in the vicinity of the two water bodies.

Discussion
Land use reveals variations in time and space over a given area. In the basin of the ancient lagoons, there has been a sharp increase in the area of bare land at the expense of forests, savannahs and plantations. Settlements too have seen a Journal of Geoscience and Environment Protection   sharp increase; this is a response to the continuous search for living space by humans especially in the current context of rampant demography (Akotossode et al., 2018). Although the area of water bodies has seen a sawtooth variation over the last thirty years in the basin, the factor to which they are more vulnerable is pollution. Indeed, the impacts of the destruction of vegetation and the occupation of space by urban areas are felt in terms of organic pollution of aquatic ecosystems in the basin. This is shown in the work of Achoh et al. (2021) which reveals the pollution of the Toho and Todougba lagoons (two major lagoons in the basin) through the evaluation of indices. Thus, the more pressure there is on the occupation of space by humans, the greater the production of organic matter and the more polluted the ecosystems of the basin will be. The evaluation of the vulnerability of the lagoons to organic pollution was carried out in order to detect areas of high pressure on the aquatic ecosystems of the basin. The results obtained on the degree of vulnerability, particularly on the first two arms of the lagoons (Toho and Todougba lagoons), show that the two ecosystems are highly vulnerable to organic pollution (73% to 100%). This very high degree of vulnerability is related to the strong pressure exerted on these two lagoons through various activities, including mainly cage aquaculture, which provides more than 50% of the total quantity of fish farmed and marketed in the Republic of Benin (DPH/MAEP, 2019). Pollution indicators were revealed in this session of the basin by Capo-Chichi et al. (2018) who had noticed an abun-dance of macroinvertebrates associated with organic-rich environments. Yoboue et al. (2018) also reported that cage aquaculture is a source of organic pollution. Therefore, the development of cage aquaculture on Toho and Todougba lagoons explains the high degree of vulnerability observed. For the other arms, the high vulnerability in the vicinity is related to the high land occupation especially in 2020 where the percentage of land occupation by settlements between 1990 and 2020 is 364.05% while forests and savannahs tend to disappear completely.

Conclusion
The land use dynamics of the Five Fingers Lagoon Basin revealed that forests and savannahs have drastically reduced in area at the expense of settlements and other factors between 1990 and 2020. With population growth and the vital needs of humans, this dynamic has a significant effect on the degree of vulnerability to pollution of the lagoons in the basin. Thus, the Toho and Todougba lagoons, which are under the influence of cage fish farming, recognized as a major source of organic pollution, are subject to a very high degree of vulnerability to pollution. It is necessary to evaluate the pollution indices in these two ecosystems of the basin which are heavily exploited in order to propose measures for a sustainable use. Finally, we suggest that an authorization for the installation of fish farmers on the lagoon should be established and made mandatory in order to better control the flow of organic matter due to fish farming. It is then necessary to determine with rigor, the carrying capacity of the two lagoons in order to have a look on the intensity of production to better control the flow of organic matter related to the fish feeding. It will still be useful to recommend a single type of feed to be used throughout the ecosystem in order to avoid disparities in the contribution of polluting factors, such as phosphorus and nitrogen. It is also advisable to pay attention to the space occupation flow in order to curb as much as possible the collateral anthropic impact.