Analysis of the Temporal Changes of Inland Ramsar Sites in Turkey Using Google Earth Engine

: Ramsar Convention (RC) is the first of modern intergovernmental agreement on the conscious use and conservation of natural resources. It provides a platform for contracting parties working together to develop the best available data, advice, and policy recommendations to increase awareness of the benefits of wetlands in nature and society. Turkey became a party of the RC in 1994, and in the years 1994 to 2013, 14 wetlands that reached the Ramsar criteria were recognized as Ramsar sites (RS). With this study, all inland RS in Turkey from 1985 to 2020 were examined, and changes in the water surface areas were evaluated on the GEE cloud computing platform using Landsat satellite images and the NDWI index. The closest meteorological station data to each RS were evaluated and associated with the surface area changes. The reasons for the changes in these areas, besides the meteorological effects, have been scrutinized using management plans and pub-lications. As a result, inland wetlands decreased at different rates from 1985 to 2020, with a total loss of 31.38% and 21571.0 ha for the spring months. Since the designation dates of RS, the total amount of water surface area reduction was 27.35 %, constituting 17,758.90 ha.


Introduction
Wetlands are essential components of the environment in terms of preserving the natural equilibrium and biodiversity with their variety of functions, including drinking water provision, sediment and nutrient storage, and flood and climate change control, as well as being used for natural remediation, fishing, feeding groundwater, and cultural purposes. They are also the most critical global genetic reserves, inhabiting 40% of all species and 12% of animal species [1]. Throughout all these roles and drives, wetlands contribute remarkably to regional and national economies. Scientific studies have indicated that 64% of the world's wetlands have disappeared since the 1900s. The leading causes for the degradation and loss of wetlands are changes in land-use; significantly increased agriculture; water diversion through dams, dikes, and canalization, particularly in river valleys and coastal areas; air and water pollution; and disposal of excess nutrients [2]. Efforts have been made to improve and protect the remaining wetlands, and strategies have been developed to appreciate the importance of wetlands. The Ramsar Convention is the first of modern intergovernmental agreement on the conscious use and conservation of natural resources. It was developed in the 1960s as a response to the increasing concerns on the conversion and destruction of wetlands, as well as for the impacts on people and biodiversity, especially on water birds [3]. It was accepted in 1971, and put into force in 1975. The United Nations Educational, Scientific, and Cultural Organization (UNESCO) acts as the legal depositary of the convention, which provides the single most global framework for intergovernmental cooperation on wetland issues and a platform of 171 contracting parties working together to develop the best available data, advice, and policy

1
A wetland should be considered internationally important if it contains a representative, rare, or unique example of a natural or near-natural wetland type found within the appropriate biogeographic region.
2 A wetland should be considered internationally important if it supports vulnerable, endangered, or critically endangered species or threatened ecological communities.

3
A wetland should be considered internationally important if it supports populations of plant and/or animal species important for maintaining the biological diversity of a particular biogeographic region.

4
A wetland should be considered internationally important if it supports plant and/or animal species at a critical stage in their life cycles, or provides refuge during adverse conditions. Specific criteria are based on water birds.

5
A wetland should be considered internationally important if it regularly supports 20,000 or more water birds. 6 A wetland should be considered internationally important if it regularly supports 1% of the individuals in a population of one species or subspecies of water bird.

7
A wetland should be considered internationally important if it supports a significant proportion of indigenous fish subspecies, species, or families; life-history stages; species interactions; and/or populations that are representative of wetland benefits and/or values, and thereby contributes to global biological diversity.

8
A wetland should be considered internationally important if it is an important source of food for fish, spawning grounds, nurseries, and/or migration paths on which fish stocks, either within the wetland or elsewhere, de-pend.

9
A wetland should be considered internationally important if it regularly supports 1% of the individuals in a population of one species or subspecies of wetland-dependent nonavian animal species.
Turkey's natural wetland area covers almost 2.3 % of the country's surface area, with 17,269 natural wetlands managed by the General Directorate of Nature Conservation and National Parks (GDNCNP) under the Ministry of Agriculture and Forestry (MoAF). Nine hundred twenty-one of these wetlands have a size greater than 8 ha, constituting 99.48% of the overall wetland areas [6]. Turkey become a party of the Ramsar Convention in 1994, and 14 wetlands have been recognized as Ramsar sites in the years 1994 to 2013. Turkey implemented the Wetlands Protection Regulation in 2002 to fulfill the obligations. The Ramsar sites of Turkey are shown in Figure 1 and are listed in Table 2 with the site numbers, protection areas, wetland types, coordinates, designation dates, and the criteria it meets [5,7].  In Turkey, almost half of the wetlands have lost their ecological and economic characteristics as a result of drying, filling, and interference in water systems in the last 40 years [8]. The most significant cause of this loss was the deliberate drying activities in the 1950s in order to combat the malaria disease. The main factors threatening wetlands in Turkey nowadays are unplanned water infrastructure projects like dams built on rivers feeding wetlands, as well as interbasin water transfers, excessive water use in agriculture, pollution, illegal hunting and fishing, other large-scale investments such as highways and mines, and management problems [9]. The factors that threaten the well-being of Turkey's Ramsar sites compiled from the existing management plans and reports are listed in Table  3. Wetland Management Plans (MP) are prepared under the coordination of the GDNCNP so as to protect wetlands, identify the problems they face, develop solutions, and ensure rational use by protecting ecological functioning and biological diversity. Wetland management plans contain all technical, administrative, and hardware structures necessary to manage the areas.
Twelve of the Ramsar sites in Turkey have management plans, and the implementation and updating times of these management plans are given in Table 3, together with the different protection statuses of these sites. Factors threatening the Ramsar sites in Turkey are given in Table 4, with the items of the threats taken from Table 3 for each Ramsar Site. It is critical to acquire correct information in order to identify wetlands of various scales (global, regional, or national), monitor them, and formulate accurate management plans. Inventory, assessment, and monitoring are vital components of effective wetland management, and provide essential data and information to support management decisions [9]. It is confirmed by many studies that remote sensing (RS) technology and data are significant and essential for determining, analyzing, and monitoring the change in wetlands through time. Since the first Earth-observing satellite (Landsat1) launched in 1972 [23], optical sensors have proved their potential to monitor large-scale land use/landcover (LULC) change on the Earth's surface. Satellite sensors of varying spatial, temporal, and spectral resolutions have been used to extract and analyze information regarding surface water from past to present. Landsat satellites (MSS, TM, ETM +, and OLI) that continuously provide medium resolution images are among the most widely used optical sensors in environmental research in the last five decades. The use of Landsat satellite images has an important place in numerous studies where the water surface areas of wetlands are extracted or the temporal changes are determined [24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42][43].
The ability of cloud computing programs has been significantly improved in recent years, and they have shown great application potential in large-scale land cover mapping [44]. Cloud computing platforms such as Google Earth Engine (GEE) provide the infrastructure to access and process large amounts of regularly updated Earth Observation data rapidly in a systematic and reproducible manner [45]. The bulk of the collection comprises RS imagery, including Landsat archive and Sentinel 1 and 2 data. GEE facilitates the use of all images archived on the platform for global, national, and regional land cover mapping and land cover change monitoring. There are studies in which the water surface areas from a regional to global scale are extracted in the GEE cloud platform. In 2016, a dataset of water body changes based on the GEE platform for global large-scale, long series was published by Pekel [46]. They used 3 million Landsat images and quantified the changes on the global surface water over the past 32 years at a 30 m resolution. All of the available Landsat images (7534 scenes) were used for the Hetao Plain in China for determining the long-term changes of the surface water area. Normalized Difference Water Index (NDWI) and Normalized Vegetation Index (NDVI) were applied to map the opensurface water from 1989 to 2019 in the GEE [47]. In the study conducted at the Tarim River Basin of China, 56284 Landsat scenes were used to determine the available surface area and generate a 30 m annual water frequency map from 1992 to 2019 [48]. Moreover, the maximal and minimal water extent in 1990, 2000, 2010, and 2017 in the Middle Yangtze River Basin in China were calculated on the GEE platform by processing 2343 scenes of Landsat images [49]. Additionally, 75593 scenes of Landsat images were processed to investigate the long-term changes of open-surface water bodies in the Yangtze River Basin in China from 1984 to 2018 [50].
Studies showing temporal changes in the Ramsar Sites of Turkey are quite limited [51][52][53][54][55][56][57][58][59][60][61][62][63][64][65][66][67][68][69][70], and a study has not yet been conducted at a national scale. With this study, all of the Ramsar sites of Turkey were given with their various features, the Inland Ramsar sites from 1985 to 2020 were examined, and the results were compared. The study aims to determine the change in the Ramsar sites over the past 35 years, as well as the effects of being designated a Ramsar site by detecting the changes since the date that they were designated as a Ramsar site. Therefore, 1985 and the designation years were considered as the reference dates. The change in water surface areas from 1985 to 2020 was evaluated using a GEE cloud computing platform, considering the most available water supply in spring months (March, April, and May) and the least available in the autumn period (September and October). The nearby meteorological station data (temperature, precipitation, and evaporation) for each Ramsar site were evaluated and associated with the water surface area changes. The reasons for areal changes, besides the meteorological effects, were scrutinized by using the available management plans and publications.

Materials
This study was carried out on the GEE cloud platform. Landsat 5 TM and Landsat 8 OLI satellite images were used, and the water surface areas between 1985 and 2020 were determined for the spring and autumn periods each year for the inland Ramsar sites. In Turkey, it was observed that wetlands had the highest surface area in the spring period and the lowest in autumn. Therefore, images from March, April, and May were used for representing the water-rich period, whereas images from September and October were used to represent the water-poor period. Seasonal composite images were produced, and water surface areas were extracted using these images. The characteristics and number of images used are given in Table 5. The closest meteorological stations to the Ramsar sites are numbered and given in Figure 1. Precipitation, evaporation, and temperature parameters for 1985-2020 were obtained from the Turkish State Meteorological Service. Meteorological stations' annual average temperature values are given in Figure 2 with their meteorological station numbers' associated Ramsar sites. Although there are great differences between the seasons and regions in Turkey, the long-term annual average total precipitation for 1981-2010 was 646mm. While there was a drought period in 2006, 2007, and 2008, positive anomalies were observed in precipitation between 2009 and 2012. In particular, 2009 was recorded as the wettest year with 804 mm precipitation, and 2008 was the driest year with 506 mm. The average temperature, which was 12.7 ° C between 1970-1978, increased to 13.8 ° C between 2006 and 2012 [73].
The temperature data for 13 meteorology stations were examined, an upward trend was observed for all of them, and the average value of the 13 stations was 13.4 °C in 1985 and 15.3 °C in 2020.
The average precipitation and evaporation values, and evaporation/precipitation rates for the years 198-2020 are compiled in Table 6. When the Ramsar sites were evaluated according to their water surface areas, it can be seen that they had different water surface area sizes. The surface areas can be divided into three groups, those with a water surface area less than 100 ha, areas in the range of 1000-10,000 ha, and large areas with 10,000 ha and above. There are three Ramsar sites under 100 ha, namely Meke Maar, Kızoren Obrouk, and Lake Kuyucuk. There are eight Ramsar areas between 1000 ha and 10,000 ha, and there are three with 10,000 ha and above, namely Lake Burdur, Lake Kus (Manyas), and Lake Uluabat. As marine/coastal Ramsar Sites interact with the sea and have different characteristics compared with inland Ramsar sites, they will be discussed in a different study, and the inland Ramsar Sites were analyzed in this study.
Digital elevation models (DEM) were created using Google Earth Pro elevation data to indicate the topography and altitude of the area from the sea level in the inland wetlands, and the water surface areas were shown on these models.

Methods
There are several techniques defining surface water extraction using Landsat data described in the literature. Water classification methods for optical imagery could be categorized into four basic types, namely thematic classification, linear unmixing, singleband thresholding, and spectral water indices [74]. Spectral indices are commonly preferred because of their ease of use and being computationally less time-consuming than others. Many indices have been developed specifically to exploit the unique spectral signature of water compared with other land cover types. NDWI was introduced in 1996 to delineate open water features using the green (Band 2) and near-infrared (Band 4) of Landsat TM [75]. A threshold value of zero was proposed for extracting surface water using the raw digital number of Landsat, where all positive NDWI values would be classified as water and negative values as non-water. Many studies comparing the other indices show that NDWI gives better results than the other indices [76][77][78][79]; therefore, NDWI was used to extract the water surface areas in this study. NDWI is expressed as follows [75]: where Green is a green band such as TM Band 2 (OLI, Band 3), and NIR is a near-infrared band such as TM Band 4 (OLI, Band 5).

Results
In Table 7, the water surface areas of inland Ramsar sites in spring are given for 1985, the designation years (DY), and 2020. The water surface areas of the Ramsar Sites in autumn are given in Table 8, for 1985, DY, and 2020. DY and 1985 were taken as the reference years and the changes in surface water areas were presented in ha (Table 7) and percentages ( Figure 3) by comparing the surface area values from 1985 to DY, and DY to 2020.
When examining the water surface areas changes and percent changes for the spring months, we considered the following: • The total water surface area in the inland Ramsar (Table 8) and percent changes (Figure 4) for the autumn months, we considered the following: The total water surface area in all the inland Ramsar sites was 52,520.65 ha in 1985, and decreased to 40,787.25 ha in 2020. The change was 22.34 % with a loss of 11,733.4 ha. When the changes in inland wetlands from 1985 to 2020 were examined, only the area of Lake Kus (Manyas) was increased by 12.59% and all of the others decreased. In Lake Kus (Manyas), Meke Maar, and Lake Kuyucuk, the water completely disappeared.
• When the inland wetlands from 1985 to DY were examined, there was an increase of 16.10% (208.2.8 ha) in Lake Kus (Manyas), 1.57% (179.28 ha) in Lake Uluabat, 0.17% (2.15 ha) in Nemrut Caldera, and 2140.05% (84.96 ha) in Lake Kuyucuk. There was a decrease in the other sites. • The change from DY to 2020 was 19.91% with a loss of 10,137.04 ha, there were spatial losses at different rates, and Lake Burdur ranked first with an area of 5142.85 ha. The areas in the different water surface area groups that underwent the most changes, Meke Maar, Lake Seyfe, and Lake Burdur, were evaluated together with the annual water surface area changes and annual meteorological data, as shown in Figures 5,  6, and 7, respectively, between 1985-2020. When Tables 7 and 8 and Figures 3 and 4 are examined, it can be seen that the most dramatic percent change occurred in Lake Kuyucuk. Nemrut Caldera, on the other hand, was the area where the least change was seen in the spring and autumn months. The annual surface area changes of these two were examined between 1985 to 2020 with the annual water surface area changes and the meteorological parameters (Figure 8 and Figure 9). The remaining four Ramsar sites were similarly examined (Figures 10-13).     When the smaller Ramsar sites below 100 ha were examined, Meke Maar decreased the most from DY to 2020 in both spring and autumn. While there was a water surface area of 2.17 ha with a decrease of 95.72% from the spring months of 2005 to 2020, it was completely dry in autumn 2020 (Figures 3 and 4).
Meke Maar, a crater lake formed by filling the islets in the middle of an extinct volcano crater [80], had the most significant change within the Ramsar sites. Its annual water surface area change for both the spring and autumn between 1985-2020, and the corresponding meteorological data representing the same terms, are given in Figure 5 (a) and (b), respectively. The water surface areas belonging to 1985, DY (2005), and 2020 are given on the produced DEM in Figure 5 (c). A photo from the Meke Maar is also added in Figure  5  The amount of water present decreased tremendously since 2005, as seen in Figure 5 (a). There was no water left in the autumn period from 2013. Considering the annual evaporation and precipitation values of the nearby meteorological station, the evaporation was approximately four times that of the precipitation (Table 6), as the averages of 35 years reflect that precipitation was 293.13 mm, whereas evaporation was 1192.77 mm. The average precipitation was less than half of Turkey's average (646mm). The lack of water during the autumn months was basically due to the high evaporation. Even though there were no evaporation data available from 2012 to 2020, temperature values showed that there would be no expected change in evaporation. Similarly, there was also no change in the precipitation values since 2012. In this case, it can be highlighted that meteorological parameters were not the primary reason for water loss in the area. Primary reason for this change was the excessive withdrawal of groundwater for irrigational purposes [8].
When Ramsar sites in the range of 1000-10,000 ha were examined, the greatest decrease from DY to 2020 was seen in Lake Seyfe. While the decrease was 71.11 % from DY (1994) to the spring months of 2020, it was observed that the area was completely dried in autumn. From 1985 to 2020, the total reduction rate was 78.6% with 6086.52 ha. Lake Seyfe's annual water surface area change for both the periods of spring and autumn between 1985-2020, and the corresponding meteorological data from 1985 to 2020, are given in Figure 6 (a) and (b), respectively. The water surface areas belonging of 1985, designation year (1994), and 2020 are given in the produced DEM in Figure 5 (c). A view from the lake is added in Figure 6  The amount of water decreased tremendously since 1988 considering the 1985-2020 surface areas of Lake Seyfe (Figure 6(a)). When meteorological data were examined, according to the averages of 35 years, precipitation was 386.26 mm, evaporation was 1422.52 mm, and the amount of evaporation was approximately 3.6 times more than the precipitation. There was a decrease in water, especially in the autumn months, since 2001, as well as disappearances in some years related to low water levels and high evaporation rates. It was stated in the management plan [20] that the water level dropped to 60-70 cm in the summer season, most of which turned into salty swamp and even dried up due to the low rainfall of the region, and ended up with drying the streams feeding the lake and the high evaporation.
When the Ramsar areas of Turkey with an area of above 10,000 ha were examined, Burdur Lake was observed to have the largest water surface area in 1985. However, it has lost 8214.93 ha area with a 39.7% decrease from 1985 to 2020. In Lake Burdur, the rate of change from DY (1994) to 2020 was 29.9%. It is one of the deepest lakes of the country, bearing a tectonic structure with a maximum depth of 110 m and an average depth of 40 m [19]. Lake Burdur's annual water surface area change for both spring and autumn for 1985-2020, and the meteorological station data from 1985 to 2020 are given in Figures 7(a)  and (b), respectively.
The water surface areas belonging to 1985, DY (1994), and 2020 are given on the produced DEM in Figure 7(c), and a view of the lake is in Figure 7(d). Studies carried out in Burdur Lake have revealed that the lake water is salty and arsenic in proportions, and cannot be used as domestic and agricultural irrigation water [19]. Therefore, it has no agricultural or domestic use. When the graph showing the annual water surface areas in Figure 7 is examined, a decrease every year in the spring and autumn months can be seen. When the meteorological data were examined (Figure 7 (b)), there was no noticeable change in the annual average precipitation trend, but there was an increase in temperature and evaporation over time. Looking at the averages of 35 years, precipitation was 418.7mm, while evaporation was 1245.44 mm, indicating that evaporation was about three times more than precipitation. In the Burdur Lake management plan [19], it was stated that the dams built on the rivers supplying water to the lake feeding were interrupted and that it was only fed by surface waters coming from precipitation, and the lake was rapidly shrinking as a result of the low precipitation and high evaporation rates.
When Figure 3 and Figure 4 are examined, the biggest change can be seen from DY to 2020 in the spring and autumn months in Lake Kuyucuk, while there was almost no change in Nemrut Caldera.
The annual water surface area changes for both spring and autumn of Lake Kuyucuk between 1985-2020, as well as the Kars meteorological station data from 1985 to 2020, are given in Figure 8(a) and (b), respectively. The water surface areas of 1985, to the designation year (2009), to year 2020 are given on the produced DEM in Figure 8 (c). A photo from the lake is added in Figure 9 (d). Lake Kuyucuk is a small freshwater lake located in the middle of a slightly wavy open land, fed by small springs, streams, and precipitations; its deepest point is 13m [21]. When the annual water surface areas of Lake Kuyucuk were examined, water was seen in the spring months, while the amount of water decreased in some years in autumn, while in some years there was no water left. It has been reported that floods occurred in the region in some years because of excessive rainfall, and the lake area was filled; however, the effective volume of the lake decreased because of the transfer of sediment to the lake through overflows [21]. The reason for the excessive amount of water in the area in DY 2009 was because it was the year with the highest rainfall between 1985-2020 (Figure 8(b)). The construction of a dam (Kars) nearby Lake Kuyucuk started in 2010, and the water retention in this dam lake started in July 2018 [85]. From this date onwards, there was no water left in Kuyucuk Lake in autumn (Figure 8(a)).
The Ramsar site where the water surface area changed the least was the Nemrut Caldera, as shown in Figure 9 (a). There was a decrease of −1.3% from 1985 to 2020 ( Figure  3). Its annual water surface area changed for both spring and autumn between 1985-2020, and the representative meteorological station data from 1985 to 2020 are given in Figure 9 (a) and (b), respectively. The water surface areas belonging to 1985, DY (2013), and 2020 are given on the produced DEM in Figure 9 (c). A photo from the Nemrut Caldera is added in Figure 9 (d). Nemrut Mountain is known as one of the largest craters of the world, with a width of about 10 km. The Crater Lake located on the mountain is the second-largest caldera lake in the world [87]. The western half of the base is covered with a lake. There are five lakes at the summit, two of which are continuous, and three are seasonal. The largest of the Nemrut lakes is the Nemrut Lake in the shape of a half-moon. The average depth of this lake is around 100 m. As the lake is far from settlements, there are no structures and facilities nearby that may cause pollution [11]. When the graph showing the annual water surface areas in Figure 9a is examined, it can be seen that there was a slight decrease every year in both the spring and autumn months. When the meteorological data were examined (Figure 9(b)), there was no noticeable change in the annual total precipitation trend, but there was an increase in temperature and evaporation over time. Over the average of 35 years, precipitation was 401.38 mm, while evaporation was 1072.99 mm, indicating that evaporation was about 2.67 times more than precipitation. Therefore, the increased evaporation trend can be considered a factor in reducing the water surface area. However, although the caldera is a closed basin and there is no outflow from the rivers, there is a noticeable decrease in the water level and this should be investigated [87].
Lake Kus (Manyas) is the second-largest Ramsar site in terms of water surface area after Burdur Lake. The lake is a broad and shallow freshwater lake with an average depth of 3 m. Its water is always turbid because it contains colloidal clay. Lake Manyas, which is very rich in plankton and bottom creatures, has allowed for the development of wildlife. Lake Manyas is the first area to be recognized as a "bird paradise" in Turkey, so it pioneered the recognition of birds, wetlands, and nature. Lake Manyas's annual water sur-face area change for both spring and autumn between 1985-2020, as well as the representative meteorological station data from 1985 to 2020, are given in Figures 10 (a) and (b), respectively. The water surface areas belonging to 1985, DY (1994), and 2020 are given in the produced DEM in Figure 10 (c). A photo from the Lake Kus (Manyas) is added in Figure 10  When the graph showing the annual water surface areas in Figure 10 (a) is examined, it can be seen that there was a decrease every year in the spring, and the trend did not change much in the autumn months. When the meteorological data were examined (Figure 10 (b)), an increase in the annual total precipitation trend and temperature was found. The average precipitation (35 years) was 704.16 mm and the average evaporation (27 years) was 1141.87 mm, indicating that evaporation was about 1.62 times more than precipitation. While the lake had a water surface area of 15,269.35 ha in 1985, it decreased to 14,919.41 ha with a 2.3% loss in 2020. The main reason for the high water level difference in the spring and autumn months is the withdrawal of water from the lake for agricultural activities [10].
Lake Uluabat is the third-largest Ramsar Site in Turkey; it is located on the bird migration route entering Anatolia from the northwest, and is the lake with the most extensive lotus beds in Turkey. The deepest part of the lake is 6 m, and it is a turbid, eutrophic freshwater lake. The water resources at and around the lake, as well as the small streams that reach the lake during rainy periods, keep the water level in balance [89]. Lake Uluabat annual water surface area change for both spring and autumn between 1985-2020, as well as the representative meteorological station data from 1985 to 2020, are given in Figure 11 (a) and (b), respectively. The water surface areas belonging to 1985, DY (1998), and 2020 are given on the produced DEM in Figure 11 (c). A photo from Lake Uluabat is added in Figure 11(d) [90]. When the graph showing the annual water surface areas in Figure 11 (a) is examined, it can be seen that there was a decreasing trend in the spring, and the trend did not change much in the autumn months. When the meteorological data were examined (Figure 11 (b)), there was an increase in the annual total precipitation trend and annual average temperature. A decreasing trend was observed in evaporation. The average for precipitation (35 years) was 707.13 mm, and the average for evaporation (32 years) was 1149.62 mm, indicating that the evaporation was about 1.63 times more than precipitation. While the lake had a water surface area of 13,479.7 ha in 1985, it decreased to 11,784.49 ha with a 12.6% loss in 2020. The area covered by the lake is shrinking because of the drainage waters coming from the agricultural areas and the sediment load carried by the streams. Environmental risks continue in the lake, whose water quality has deteriorated due to pollution, and protection cannot be provided to sustain the biodiversity in the lake [10,14].
Sultan Marshes is a rare ecosystem where fresh and saltwater ecosystems coexist and form the feeding, sheltering, and hatching area of 301 bird species, including endangered species. Its annual water surface area changes for both spring and autumn between 1985-2020, as well as the representative meteorological station data from 1985 to 2020, are given in Figures 12 (a) and (b). The water surface areas belonging to 1985, DY (1994), and 2020 are given on the produced DEM in Figure 12 (c). A photo from the Sultan Marshes is added in Figure 12  When the graph showing the annual water surface areas in Figure 12 (a) is examined, it can be seen that there was a decrease every year in both the spring and autumn months. In general, the water surface area became small in the autumn months and even decreased to almost non-existence in some years. Furthermore, when the meteorological data were examined (Figure 12 (b)), there were decreases in the annual total evaporation trend. On the other hand, the annual average temperature trend increased as in other meteorological stations, and a stable precipitation trend was observed below the average for Turkey. After the all-year meteorological data were examined, precipitation was 386.26 mm while evaporation was 1422.52 mm, indicating that evaporation was about 3.68 times more than precipitation.
The increase in population, the expansion of agricultural areas as much as possible, and the inclusion of irrigated agricultural products in field agriculture significantly deteriorated the ecosystem. Water is supplied to the fields from the dams on the springs feeding the Sultan Marshes, and some of it is recycled to the lakes through drainage channels. However, with the increase in the water required for agricultural activities, sufficient water cannot be sent to the lakes, which is the main reason for the drying up of Sultan Marshes [10,17].
The Kızoren obrouk, which has the smallest surface area among the Ramsar Sites, is one of about twenty obrouks in the Konya Closed Basin, with a width of 300 m and a depth of 145 m. As in other obrouks, the water color, which is dark blue and green in the first months of summer, turns into a clear indigo color as the summer progresses. The Kızoren obrouk annual water surface area change for both spring and autumn periods between 1985-2020, as well as the representative meteorological station data from 1985 to 2020, are given in Figures 13 (a) and (b), respectively. The water surface areas belonging to 1985, DY (2006), and 2020 are given on the produced DEM in Figure 13 (c). A photo from the obrouk is added in Figure 13  Because of their formation, the obrouks have the same level of water as the groundwater. The water level decreased significantly in Kızoren Obrouk, as in other lakes in the Konya Basin, where the groundwater level drops to a depth of 45 m. In the past, the level of obrouk water was at a level that a pump could draw. Decreased groundwater due to the intensive agricultural activities carried out in the Konya Closed Basin also reduced the level of the obrouk water, making it impossible to draw water from the obrouks [10]. The water level, which was only 1-2 m from the surface in 1988, is now 28 m below the surface level [93]. When the graph showing the annual water surface areas in Figure 13 (a) is examined, it can be seen that there was a decrease every year both in the spring and autumn months. In addition, the graph shows that there is more water surface area in autumn than in spring. The reason for this is the decrease in the groundwater level in the spring months due to agricultural irrigation. Because of the narrowing structure of the obrouk, the water surface area decreases as the water level decreases.

Discussion
The water surface area change in Turkey's inland Ramsar sites was investigated over a 35-year period. In general, there was a decrease in all Ramsar sites in various ratios from 1985 to 2020. As can be summarized in Table 4, all Ramsar Sites are faced with numerous threats. After reviewing the related articles, reports, and management plans, improper water use can be accounted as the main reason underlying the decrease in the water surface area, except for Burdur Lake, where water cannot be used for domestic and agricultural usage. Human interventions without considering the balance of protection threaten the existence of these areas; therefore, the ecosystem services they offer are decreasing day by day.
When the meteorological parameters related to the changes in the Ramsar areas were examined, it was seen that the decrease in water surface areas in many sites was also related to meteorological events. During the 35-year period, although there was no extreme change in general, there was an increase in the amount of evaporation at some stations, and no significant change was observed in the precipitation trend. In general, the amount of reduction was higher in areas with high evaporation/precipitation rates. When the temperature data for 13 meteorological stations were examined, an upward trend was observed for all of them, and the average value increase was 1.9 °C in 35 years.
The Turkish State Meteorological Service produces regional climate projections using global model outputs such as RCP (Representative Concentration Pathway), RegCM4 (Regional Climate Model-4), and the downscaling method [73]. Model outputs were produced by selecting the period of 1971-2000 as the reference period, and 2013-2099 for the projections. These studies indicate that temperatures in Turkey will increase by 0.5 to 3.0 °C until 2050, and 0.5 to 4.0 °C until 2100 with the optimistic scenario (RCP 4.5). According to the pessimistic scenario (RCP 8.5), the expected increase will be higher; by 0.9 to 3.5 °C until 2050, and by 0.9 to 6.3 °C until 2100. In terms of precipitation amounts, as stated in the optimistic scenario, positive anomalies are expected throughout the country until 2040-2050. After these years, decreases are expected in the average precipitation; as the number of evaporation increases with the increase in temperature, the danger of drying increases in shallow wetlands. Although a two-dimensional image analysis was performed for the surface area changes in this study, the water volume for the maintenance and stability of wetlands, which is planned in the continuation of this study, should be taken into account.
As the change in the water cycle that will occur due to climate change will affect water resources, it is essential to protect/manage wetlands in a rational way in order not to destroy them. Regarding the protection/management of these areas, many areas have more than one protection status (as shown in Table 4). In almost all wetlands of Turkey, management plans were prepared gradually, and protection zones were determined with the obligations arising from the Ramsar Convention and national laws. Management plans of twelve Ramsar Sites have been completed (Table 4), and plans for Kızoren Obrouk and Meke Maar are under preparation. However, with this study, it has been revealed that there were significant aerial losses in the inland wetlands from the dates they were declared Ramsar sites until 2020. In addition, when the implementation dates of management plans in Table 4 are compared with the water surface area graphs of the sites, it can be seen that the water surface areas decreased after the management plans were revealed and put into implementation. In the presence of these management plans and multiple protection statuses, the loss seen is thought-provoking. This results from the fact that these plans are not effective enough, they are not taken into account by society, and the process is not well managed by the authorities. Consequently, these areas have not been adequately protected, and management plans are not fully implemented.

Conclusions
It is crucial to adopt the rational use principle for protecting wetlands and to prepare the management plans of the areas accordingly. The management plan should be understood as the fundamental law that organizes all conservation and utilization activities to be carried out within the wetlands of concern. To maintain the protection-utilization balance, all stakeholders benefiting from the wetland should adopt the area where they live; continue their economic activities; and act with awareness of the economic, social, and ecological importance of such areas. For this reason, effective communication and cooperation between relevant institutions and organizations, local governments, and non-governmental organizations at every stage of the preparation and implementation of wetland management plans is necessary so that disputes are resolved, and compromise and coordination are provided. Many of the Ramsar sites are affiliated with different ministries or directorates, as they have more than one protection status simultaneously, resulting in a confusion of authority in the management of the sites. This issue should be taken into account and resolved first.
Another critical issue is that the developments should be monitored and reviewed concurrently with the implemented wetland management plans. Monitoring shows the progress wetland management plans have made towards achieving their goals, making it easier to identify whether everything is in order and whether adaptation is required for facing the dynamic conditions. Unfortunately, there has not yet been an established mechanism for monitoring wetland management plans in Turkey. As seen in this study, remote sensing tools and methods are of great importance in monitoring wetlands and determining changes, so they should be used effectively by relevant institutions and organizations for sustainable management and conservation.