Trends in the Recent Evolution of Coastal Lagoons and Lakes in Galicia (NW Iberian Peninsula)

: Coastal lagoons are habitats of great environmental value. However, they are currently subject to major threats, particularly due to increasing sea levels. This study aims to identify changes—both natural and induced by anthropic activity—and their impact on the recent evolution of three di ﬀ erent types of coastal lagoons in Galicia (Louro, Vix á n, and Xuño). The application of information obtained through laser imaging detection and ranging (LiDAR) techniques suggests that the outer limits of the three lagoon systems have not experienced any relevant changes in the last 60 years (i.e., no occupation of the lagoon area has been identiﬁed). However, the internal conﬁguration of these wetland areas has experienced some alterations. A generalized increase in the area occupied by macrophytic communities ( Phragmites australis , Scirpus maritimus , Juncus maritimus , etc.) has been observed. Image interpretation by geographic information systems (GIS) and ﬁeld surveys suggest that the area currently occupied by macrophytes experienced a 7% to 63% increase at the expense of the free water body. This loss of ﬂooded area is consistent with the increase in sedimentation rates associated with the convergence of several causes, such as the abandonment of traditional macrophyte biomass harvesting and agricultural activities around the lagoons, the expansion of riparian forests, and sediment contributions by erosion due to recurrent forest ﬁres within the drainage basins of each lagoon. Finally, water and sediment composition suggest that, of the three studied lagoons, two of them (Louro and Vix á n) are included within the deﬁnition of “coastal lagoons” (habitat code 1150) by the Habitats Directive (Directive 92 / 43 / EEC), while the Xuño lagoon should be considered a “natural eutrophic lake” (habitat code 3150). Authors' Contributions: Conceptualization, A.P.A., P.M., X.L.O.P and P.F.S.; methodology, A.G.P and P.F.S.; formal analysis, P.M., P.F.S, A.G.P and X.L.O.P.; investigation, A.P.A., P.F.S., P.M. and X.L.O.P; writing— original draft preparation, P.F.S., X.L.O.P and A.G.P.; writing—review and editing, P.F.S., A.G.P, X.L.O.P.


Introduction
Coastal wetlands are fragile and currently highly endangered ecosystems. For centuries, coastline anthropization has led to the alteration or destruction of many of these habitats. According to [1], around 50% of the world's wetlands have been lost since 1900. Although the rate of wetland loss has decelerated in Europe and has been low in North America since the 1980s, it remains high in Asia due to the ongoing, fast-paced, large-scale transformation of natural coastal and inland wetlands. In Europe, the progressive destruction of wetlands can be traced back at least to Roman times [1,2], while it has been taking place since the 17th century in North America and southern Africa [3,4], and for at least 2000 years in China [5]. Additionally, global change, including the associated increase in sea

Environmental Setting
The three lagoons are located in the Atlantic coast of Galicia (NW Iberian Peninsula). This coastline features an important number of lagoons of great environmental interest, all of which are included in the Natura 2000 network. The geological substrate consists of granites and metamorphic rocks (mainly schist). This region has a maritime Mediterranean climate, annual average temperature is 14.6-16.2 • C, and average rainfall ranges from 1244 to 1828 mm/year ( Figure A1). Seventy percent of this rainfall is concentrated in the period between October and March, while intense drought periods occur in the summer [12].
Topographically speaking, the Louro lagoon, with an average elevation of 1.84 m above sea level (masl) calculated using laser imaging detection and ranging (LiDAR) data, and the Vixán lagoon (3.14 masl) are located between two rocky protrusions connected by a sandbar closing the direct outlet of small freshwater tributaries to the sea (Figure 1). Both of them can be classified within the coastal lagoon typology, showing a triangular (∆) shape with its apex pointing landward and its base along the seaward shore, which corresponds to a sandbar forming a beach, a well-developed dune system, and a lagoon behind it that is separated from the sea by narrow barriers of land [13]. A variety of coastal water bodies exists, ranging from areas totally enclosed by land to areas primarily exposed to the sea. The enclosed features, which are primarily shielded from the sea, are called coastal ponds and lakes. Bodies of water with outlets to the sea are called coastal bays or lagoons, depending on their shapes. Oertel in [13,14] called these features "barrier lagoons" because of their association with coastal barriers. In general, coastal lakes and lagoons form secondary coastlines behind the main ocean coastline. The inner shores of these secondary coastlines are shielded from direct exposure to ocean waves, the connection between these lagoons and the sea occurs through a channel (Figure 1), and the main difference between them lies in its functioning. In the Louro lagoon, the channel is still active; therefore, in winter, concurrently with the period of maximum water availability in the system, water overflows to the ocean, while occasional overwash or tidal pumping can lead to ephemeral periods of seawater entrance. Conversely, the channel in the Vixán lagoon is permanently closed; only prior to 1992, when the Corrubedo Dune Complex and Carregal and Vixán Lagoons Natural Park was designated, locals used to open it to drain nearby croplands. Nowadays, this activity is completely prohibited.
The Xuño lagoon, with a 5.28 m elevation asl, is smaller in size and is located on the limit of the dune system associated with the Areas Longas beach (Porto do Son, A Coruña). It can be classified within the type known as cat's eye ponds due to its small size and lenticular shape; alternatively, it could also be classified as a coastal back-barrier perched lake [15]. Its elevation is higher than in the case of the other lagoons and, therefore, seawater never reaches it. This lagoon has a small channel through which freshwater overflows into another channel connecting another coastal lagoon of larger size (San Pedro de Muro lagoon) with the sea. coastal water bodies exists, ranging from areas totally enclosed by land to areas primarily exposed to the sea. The enclosed features, which are primarily shielded from the sea, are called coastal ponds and lakes. Bodies of water with outlets to the sea are called coastal bays or lagoons, depending on their shapes. Oertel in [13,14] called these features "barrier lagoons" because of their association with coastal barriers. In general, coastal lakes and lagoons form secondary coastlines behind the main ocean coastline. The inner shores of these secondary coastlines are shielded from direct exposure to ocean waves, the connection between these lagoons and the sea occurs through a channel (Figure 1), and the main difference between them lies in its functioning. In the Louro lagoon, the channel is still active; therefore, in winter, concurrently with the period of maximum water availability in the system, water overflows to the ocean, while occasional overwash or tidal pumping can lead to ephemeral periods of seawater entrance. Conversely, the channel in the Vixán lagoon is permanently closed; only prior to 1992, when the Corrubedo Dune Complex and Carregal and Vixán Lagoons Natural Park was designated, locals used to open it to drain nearby croplands. Nowadays, this activity is completely prohibited. The Xuño lagoon, with a 5.28 m elevation asl, is smaller in size and is located on the limit of the dune system associated with the Areas Longas beach (Porto do Son, A Coruña). It can be classified within the type known as cat's eye ponds due to its small size and lenticular shape; alternatively, it could also be classified as a coastal back-barrier perched lake [15]. Its elevation is higher than in the case of the other lagoons and, therefore, seawater never reaches it. This lagoon has a small channel through which freshwater overflows into another channel connecting another coastal lagoon of larger size (San Pedro de Muro lagoon) with the sea.

Lagoon Delimitation and Zoning
The functional area of each lagoon was determined using information obtained using LiDAR (laser imaging detection and ranging) by the Spanish National Geographical Institute (Instituto Geográfico Nacional de España, IGN). These data allowed for the generation of a digital terrain model (DTM) with 1 m spatial resolution. From this file, the functional area of each lagoon was delimited, considering the elevation of these areas, as well as the slopes associated with their borders.
Zoning of each one was performed by photointerpretation of available aerial images and orthophotographs from the years 1956,1984,1989,2001,2004,2008,2010,2014, and 2017, except for the Louro lagoon, for which the 2001 aerial image was not available. Photograms of the first two years consist of black and white images; images for the first year were obtained by the American flight, B series (1956)(1957), while those for the second year were obtained from the national flight (1980)(1981)(1982)(1983)(1984)(1985)(1986). Photograms for 1989 correspond to the coastal flight (1989)(1990)(1991), while 2001 images come from the five-year flight (1999)(2000)(2001)(2002)(2003). From 2004 onwards, color images were obtained from the National Aerial Orthophotography Program (Plan Nacional de Ortofotografía Aérea, PNOA) of the National Geographical Institute. The images used in this study represent the total available number of high-resolution images of the study areas for the analyzed period. In this case, we decided not to select satellite images because their spatial resolution would not allow for their correct comparison with aerial images and orthophotographs.
Each one of the units was delimited using ArcGIS 10.5 at 1:1000 scale for all the dates, using hectares (ha) as the measurement unit to analyze variations in area coverage. The delimitation process in aerial images was done manually using different remote sensing techniques, such as color, texture, and context. It is worth noting that uncertainty is greater for older dates, due to the lower resolution of the images and the greater complexity of black and white photointerpretation [16].
To determine the temporal changes that occurred in each lagoon, the following units were delimited: 1.
Free water: area permanently covered by a water sheet.

2.
Phragmites: area colonized by the community defined by Phragmites australis reed beds, characterized by the highest freshwater influence.

3.
Juncus/Scirpus: area with a higher seawater influence, occupied by Juncus maritimus, Juncus acutus, Scirpus maritimus, and S. lacustris, which often appear interspersed in the study areas.

4.
Sand or dune area: sedimentary area, either forming dunes or intertidal flats.
The selection of the number of categories, as well as their definition, was based on the analysis of all the available images for each lagoon and the field surveys of the areas on different dates. However, this delimitation refers to the dominant species within each section and does not exclude the presence of other species sharing the same space. From this spatial information, the area and proportion covered by each category were estimated for each study year. These data allow for analyzing variations that have occurred over time and explanations as to their possible causes. Finally, seasonal variations were also analyzed by comparing images taken in late summer (September/October), corresponding to the period with lower water availability, with images taken in late spring (June).

Geochemical Analysis of Water and Sediments
Lagoon water and sediments were sampled on different months during the 2009-2011 period. Water samples were taken in Teflon jars that had been previously treated with HCl (5%) and then rinsed several times with Milli-Q water. In the laboratory, samples were filtered through a 0.45 µm pore size.
Redox potential, pH, and electrical conductivity were determined in situ using a multiparameter probe (Hanna instruments). Salinity was determined using a Hanna refractometer. Concentrations of nutrients (NH 4 + , NO 3− , and PO 4 3− ) and chlorides were determined by ion chromatography using a Dionex 4500i system (Dionex Inc., Sunnyvale, CA, USA). Sediment samples were taken using 150 cm long PVC (polyvinyl chloride) tubes with 80 mm inner diameter. In the laboratory, sediment samples were cut in 3-to 5-cm thick segments. Sediment samples were used to determine the influence (both current and past) of seawater on each lagoon, using the degree of pyritization (concentration of pyritic S) in sediments relative to organic C content as proxy [17]. The organic C/pyritic S ratio is a proxy that allows distinguishing anoxic freshwater sediments (with low concentrations of pyrite, FeS 2 ) from anoxic marine sediments (with high pyrite contents) [17][18][19].
Organic C content was determined using a Leco autoanalyzer (after removal of carbonates with HCl 1M), and pyritic S content was determined by sequential Fe extraction, which allows for extracting Fe associated with pyrite (FeS 2 ) (see [20,21] for additional details). From pyritic Fe content, pyritic S content was determined stoichiometrically.

Delimitation and Evolution of the Lagoons in the Last 60 Years
The evolution of the occupation of the different habitats identified within the Louro lagoon is shown in Figure 2. The Louro lagoon is divided into two subsystems. The first one is located in the distal section of the lagoon and constitutes the permanently flooded area. The construction of a dyke before the 1950s contributed to this flooding. The second subsystem is constituted by a sand flat behind the dunes, which remains flooded from early autumn, when it starts to fill with the first rains, until late July. In late September, the water-covered area of the Louro lagoon is at its minimum ( Figure 2b).
More specifically, it is worth noting that no changes have been observed regarding the area and configuration of the lake basin; however, strong interannual oscillations have been observed in terms of area covered by free water. The photointerpretation analysis shows that the degree of flooding of the lagoon varies substantially depending on the season and the water availability. Thus, the minimum free water body was found in September 2004, with 4.4 ha, while the maximum was found in 2017, with 15.99 ha, representing a 40% variation relative to the total area of the lagoon. These variations can be clearly observed in Figure 2e, which represents free water bodies for 2004 and 2014, two years that diverged widely in terms of water availability.
The area covered by the different plant communities also showed important interannual variations. Progressive colonization is observed in the area behind the dunes and in the sand flat, which was reduced by 63% between 1956 and 2017. The area covered by Phragmites australis was virtually stable (10-13%), while the area occupied by the vegetation community dominated by Scirpus maritimus and Juncus maritimus experienced a marked increase since 1956, going from 2.19 ha to over 4 ha in the first decade of the 21st century.
Likewise, the area dominated by Paspalum vaginatum greatly increased from occupying 6% of the lagoon in 1956 to 56% in 2010. In the Xuño lagoon, the permanently flooded area represented a minimum portion of the lagoon (1%-3%) during the 1956-2008 period ( Figure 3). From 2008 onwards, this area was reduced to less than 1%. This trend continued between 2014 and 2017, when the water sheet area decreased from 0.008 ha to 0.004 ha. These values clearly show the advanced siltation state of the lagoon. As for the area colonized by the Phragmites australis community, its distribution varies throughout the analyzed period, reaching its maximum extension (61%) in 2017 and its minimum extension (28%) in 1989. As in the case of the aforementioned formation, the vegetation community dominated by Juncus maritimus and Scirpus maritimus varied irregularly during the 1956-2017 period, reaching its maximum extension in 2001 (41%) and decreasing to its minimum in 2014 (31%).
Likewise, aquatic vegetation also followed an irregular pattern during this period. The maximum area covered by aquatic vegetation was 31% (0.66 ha) in 2008, while its minimum was 3% (0.07 ha) in 2017. Moreover, this lagoon showed strong seasonality in terms of degree of flooding, In the Xuño lagoon, the permanently flooded area represented a minimum portion of the lagoon (1-3%) during the 1956-2008 period (Figure 3). From 2008 onwards, this area was reduced to less than 1%. This trend continued between 2014 and 2017, when the water sheet area decreased from 0.008 ha to 0.004 ha. These values clearly show the advanced siltation state of the lagoon. As for the area colonized by the Phragmites australis community, its distribution varies throughout the analyzed period, reaching its maximum extension (61%) in 2017 and its minimum extension (28%) in 1989. As in the case of the aforementioned formation, the vegetation community dominated by Juncus maritimus and Scirpus maritimus varied irregularly during the 1956-2017 period, reaching its maximum extension in 2001 (41%) and decreasing to its minimum in 2014 (31%).
Likewise, aquatic vegetation also followed an irregular pattern during this period. The maximum area covered by aquatic vegetation was 31% (0.66 ha) in 2008, while its minimum was 3% (0.07 ha) in 2017. Moreover, this lagoon showed strong seasonality in terms of degree of flooding, ranging from complete flooding of the lagoon basin in the winter to the lack of free water bodies by the end of the summer (Figure 3b). ranging from complete flooding of the lagoon basin in the winter to the lack of free water bodies by the end of the summer (Figure 3b). In the Vixán lagoon, the free water body experienced a generalized decrease during the 1956-2017 period (Figure 4). In 1984, the estimated area covered by free water was 3.81 ha (24%), while in 2017 it was 3.30 ha (21%). This means a loss of more than 0.50 ha of free water body.
Like the aforementioned lagoons, the Vixán lagoon also experiences a marked seasonal variability, which can represent a 30%-40% decrease in the free water body depending on the intensity of the summer drought, during which water remains only in the most distal area of the lagoon.
The dune and sand area also experience strong variations, going from a 12% occupation in 1989 to their almost complete disappearance in the last study years. The area covered by the dune environment was increasingly occupied by the community composed of Scirpus maritimus and several species in the genus Juncus, whose area increased from 24% in 1956 to 31% in 1989. From that year onwards, its proportion has remained stable at around 30%. Finally, an expansion of the Phragmites australis In the Vixán lagoon, the free water body experienced a generalized decrease during the 1956-2017 period (Figure 4). In 1984, the estimated area covered by free water was 3.81 ha (24%), while in 2017 it was 3.30 ha (21%). This means a loss of more than 0.50 ha of free water body.
Like the aforementioned lagoons, the Vixán lagoon also experiences a marked seasonal variability, which can represent a 30-40% decrease in the free water body depending on the intensity of the summer drought, during which water remains only in the most distal area of the lagoon.
The dune and sand area also experience strong variations, going from a 12% occupation in 1989 to their almost complete disappearance in the last study years. The area covered by the dune environment was increasingly occupied by the community composed of Scirpus maritimus and several species in the genus Juncus, whose area increased from 24% in 1956 to 31% in 1989. From that year onwards, its proportion has remained stable at around 30%. Finally, an expansion of the Phragmites australis community has also been detected during the 2010-2017 period, reaching 53% of the total area, while in 1989 it represented only 41%.

Physicochemical Characterization and Composition of Water
Water pH ranged from slightly acidic values in the Xuño lagoon (6.1-7.3) to alkaline values in the Louro lagoon (7.8-9.8; Table 1). The extremely high values observed in the Louro lagoon are related to the increase in photosynthetic activity during the summer. Redox potential ranged between 213-485 mV, with both maximum and minimum values corresponding to the Xuño lagoon. The lowest values were found in late summer, when the free water volume is very small and organic load is high [22].
The main difference in water composition was observed for parameters related to ion concentration. Thus, salinity, electrical conductivity, and chlorides showed their highest values in the

Physicochemical Characterization and Composition of Water
Water pH ranged from slightly acidic values in the Xuño lagoon (6.1-7.3) to alkaline values in the Louro lagoon (7.8-9.8; Table 1). The extremely high values observed in the Louro lagoon are related to the increase in photosynthetic activity during the summer. Redox potential ranged between 213-485 mV, with both maximum and minimum values corresponding to the Xuño lagoon. The lowest values were found in late summer, when the free water volume is very small and organic load is high [22].
The main difference in water composition was observed for parameters related to ion concentration. Thus, salinity, electrical conductivity, and chlorides showed their highest values in the Louro lagoon, followed by Vixán and Xuño. The low salinity (<1% ), electrical conductivity (0.27-0.88 dS m −1 ), and chloride concentration (47-228 mg L −1 ) of water in the Xuño lagoon is evidence to its continental origin, while the high conductivity and salinity of the water in the Vixán and Louro lagoons are characteristic of brackish waters and indicate a clear marine influence.
Nutrient concentration in the water (NH 4 + , NO 3 − , and PO 4 3− ) showed abnormally high values in the three lagoons (Table 1), but particularly in the Vixán and Xuño lagoons in late summer. Finally, it is worth highlighting the marked seasonality shown by all of them as a consequence of a major decrease in water volume at the end of the summer period.

Organic C/Pyritic S Ratio in Sediments
Sediments from the three lagoons are characterized by pH values that are close to neutral (pH 6.0-7.5) and a redox potential that is characteristic of strongly reduced sedimentary environments (Eh <100 mV) (data not shown [22]). Figure 5 shows the organic C/pyritic S ratio for surface and deep sediments (0-100 cm). The results show extremely low contents of pyritic S in sediments from the Xuño lagoon, which suggest that this lagoon system has always been a freshwater lagoon. Conversely, the high pyrite content observed in the Vixán and Louro lagoons suggest that they have been clearly receiving some marine influence since their formation. Figure 5. Pyritic S to organic C ratio in sediments. Solid lines represent the typical organic C/pyritic S ratio for normal modern marine sediments (C/S ratio = 2.8 ± 0.8, data from [18,23]), and the dashed line represents the C/S ratio for freshwater sediments, data from [17].

Discussion
Despite the current environmental protection of lagoons and coastal lakes, these systems are progressing towards their collapse. The main observed cause is clogging, which leads to the disappearance of the free water body and to a progressive loss of biodiversity (mainly of aquatic birds [24]). Future actions must be taken in accordance with the origin and dynamic of each lagoon system in particular.
In relation to the possibility of improving these analyses, the images used in this paper were taken with irregular periodicity. Unfortunately, for this reason, variations in the area covered by the lagoons could not be analyzed in more detail. It would be useful to obtain images in different seasons to appropriately test seasonal variations. For this purpose, UAVs (unmanned aerial vehicles) could be a highly suitable tool.

Types of Lagoons According to their Location and Geochemistry (Coastal Lagoons and Natural Eutrophic Lakes)
The Habitats Directive (more formally known as Council Directive 92/43/EEC on the conservation of natural habitats and of wild fauna and flora) establishes the types of habitats of community interest that must be protected by the European Union member states. Among them are coastal lagoons (habitat code 1150). According to the Interpretation Manual of European Union Habitats (2003), "lagoons are expanses of shallow coastal salt-water, of varying salinity and water volume, wholly or partially separated from the sea by sand banks or shingle, or, less frequently, by rocks. Salinity may vary from brackish water to hypersalinity depending on rainfall, evaporation and through the addition of fresh seawater from storms, temporary flooding of the sea in winter or tidal exchange". The results regarding water and sediment composition in the three lagoons indicate that Louro and Vixán clearly fulfill the requirement of being a saline or brackish environment (Table 1; Figure 5). However, the so-called Xuño "lagoon" is, and has always been, a freshwater lake with no direct connection with the sea. Therefore, its current classification within habitat 1150 ("Coastal lagoons") by the environmental authorities of the regional government (Xunta de Galicia) is incorrect, and it should be assigned to the category "Natural eutrophic lakes with Magnopotamion or Figure 5. Pyritic S to organic C ratio in sediments. Solid lines represent the typical organic C/pyritic S ratio for normal modern marine sediments (C/S ratio = 2.8 ± 0.8, data from [18,23]), and the dashed line represents the C/S ratio for freshwater sediments, data from [17].

Discussion
Despite the current environmental protection of lagoons and coastal lakes, these systems are progressing towards their collapse. The main observed cause is clogging, which leads to the disappearance of the free water body and to a progressive loss of biodiversity (mainly of aquatic birds [24]). Future actions must be taken in accordance with the origin and dynamic of each lagoon system in particular.
In relation to the possibility of improving these analyses, the images used in this paper were taken with irregular periodicity. Unfortunately, for this reason, variations in the area covered by the lagoons could not be analyzed in more detail. It would be useful to obtain images in different seasons to appropriately test seasonal variations. For this purpose, UAVs (unmanned aerial vehicles) could be a highly suitable tool.

Types of Lagoons According to Their Location and Geochemistry (Coastal Lagoons and Natural Eutrophic Lakes)
The Habitats Directive (more formally known as Council Directive 92/43/EEC on the conservation of natural habitats and of wild fauna and flora) establishes the types of habitats of community interest that must be protected by the European Union member states. Among them are coastal lagoons (habitat code 1150). According to the Interpretation Manual of European Union Habitats (2003), "lagoons are expanses of shallow coastal salt-water, of varying salinity and water volume, wholly or partially separated from the sea by sand banks or shingle, or, less frequently, by rocks. Salinity may vary from brackish water to hypersalinity depending on rainfall, evaporation and through the addition of fresh seawater from storms, temporary flooding of the sea in winter or tidal exchange". The results regarding water and sediment composition in the three lagoons indicate that Louro and Vixán clearly fulfill the requirement of being a saline or brackish environment (Table 1; Figure 5). However, the so-called Xuño "lagoon" is, and has always been, a freshwater lake with no direct connection with the sea. Therefore, its current classification within habitat 1150 ("Coastal lagoons") by the environmental authorities of the regional government (Xunta de Galicia) is incorrect, and it should be assigned to the category "Natural eutrophic lakes with Magnopotamion or Hydrocharition-type vegetation" (habitat code 3150). However, this is not an isolated case. Other lagoons in the Atlantic coast of the Iberian Peninsula, such as the Doniños and Traba lakes (both included in the Spanish Natura 2000 network), were also classified as coastal lagoons, but have recently been considered as coastal back-barrier perched lakes [15,24]. This typology is defined as freshwater bodies that are elevated over sea-level and are not directly subjected to the inflow of seawater, a definition that fits the characteristics of the Xuño lake.
Ultimately, lagoon classification is not always easy and clear, and neither is their categorization as listed in Annex I of the Habitats Directive, which mixed multidisciplinary information and used various environmental proxies in the majority of cases (see e.g., [15,24,25]).

Evolution and Environmental Implications of the Reduction in the Free Water Body and of Changes in the Areas Colonized by Vegetation
The studied lagoons and lake represent different scenarios in terms of hydrologic dynamics and vegetation-covered area. In all of them, an increase in macrophyte coverage at the expense of the free water body was observed. The three lagoons are experiencing a progressive siltation and eutrophication process, which alters their natural evolution, and which could eventually affect the natural values inherent to these lagoons. It is worth highlighting the loss of habitat quality to host stable populations of aquatic birds (e.g., diving Anatidae, Fulicula atra, Tachybaptus ruficollis) and marsh birds (especially Emberiza schoeniclus and Circus aeruginosus) inhabiting these lagoons.
Siltation and eutrophication are caused by several factors. The first one of them is the abandonment of traditional practices such as the annual harvesting of aerial biomass produced by macrophytes (mainly Pragmites australis, Scirpus maritimus, and Juncus maritimus) for cattle bedding. Systematic macrophyte harvesting reduced siltation and eutrophication of the lagoons [26]. These species are excellent nutrient extractors [27,28]. Other elements worth mentioning are changes in land use around the lagoons which have taken place in recent decades. During the 1960s and 1970s, areas surrounding the lagoons were consistently used for crop growing throughout the year (Figure 3). However, since the last decades of the 20th century, this activity was progressively abandoned, leading to the rapid regeneration of riparian forests composed mainly of common alder (Alnus glutinosa) and willow (Salix atrocinerea). The expansion of tree-covered areas around the lagoons led to a decrease in water availability due to the large water volume extracted from the ground by these species due to evapotranspiration [29]. Finally, siltation has also been compounded by sediment contributions due to recurrent forest fires within the drainage basin of the lagoons, particularly in the Louro and Vixán areas.
The evolution of coastal lagoons does not represent a problem in the Galician context. In recent years, numerous studies have analyzed variations and their implications in the characteristics of lagoons, focusing on issues such as eutrophication. Clear examples of study cases are the analyses performed in the south of Spain [30][31][32], although diverse authors have studied other areas in the Iberian Peninsula [33] and the United Kingdom [34].

Conclusions
The following main conclusions can be drawn from this study:

•
The three lagoons are experiencing a progressive siltation and eutrophication process, which decreases water quality and substantially deteriorates their environmental function. The progressively increasing coverage by macrophytic communities at the expense of the area covered by free water body and quality is compromising the lagoons' capacity to host stable populations of aquatic and marsh bird species.

•
An adequate management of these spaces requires, first and foremost, a thorough understanding of their origin and dynamics. Knowledge of these aspects is currently incomplete, which has led to their incorrect classification in terms of habitat type, according to directive 92/43/CEE.

•
The allocation of lagoons and lakes to the categories listed in Annex I of the Habitats Directive is very complex and requires a multidisciplinary approach.