The Gavkhouni Wetland Dryness and Its Impact on Air Temperature Variability in the Eastern Part of the Zayandeh-Rud River Basin, Iran

: The Gavkhouni wetland provides many environmental and economic beneﬁts for the central region of Iran. In recent decades, it has completely dried up several times with substantial impacts on local ecosystems and climate. Remote sensing-based Land Surface Temperature (LST), Normalized Difference Water Index (NDWI) and Normalized Difference Vegetation Index (NDVI) in combination with in-situ data were used to investigate the trend of the Gavkhouni wetland dryness and the associated impact on the variability of local air temperature. The results indicate that the wetland has increasingly experienced drier conditions since the year 2000. The wetland was almost completely dry in 2009, 2011, 2015 and 2017. In addition, the results show that Gavkhouni wetland dryness has a signiﬁcant impact on local climate, increasing the mean seasonal air temperature by ~1.6 ◦ C and ~1 ◦ C in spring and summer, respectively. Overall, this study shows that remote sensing imagery is a valuable source for monitoring dryness and air temperature variations in the region. Moreover, the results provide a basis for effective water allocation decisions to maintain the hydrological and ecological functionality of the Gavkhouni wetland. Considering that many factors such as latitude, cloud cover, and the direction of prevailing winds affect land surface and air temperatures, it is suggested to use a numerical climate model to improve a regional understanding of the effects of wetland dryness on the surrounding climate.


Introduction
Wetlands, as aquatic ecosystems, provide a variety of benefits to wildlife and humans in many different ways [1,2]. Some of the services provided by wetlands include reducing greenhouse gasses, regulating climate [3], purifying water, storing pollutants in their soils, and vegetation maintaining water supply [4], providing breeding habitat for various species of fauna, and flora, and biodiversity function [5]. Wetlands are one of the vital micro-components of climate as they have their own specific moisture cycles with aquatic plants [6]. Changing land uses, such as forests, agricultural lands, built-up areas, and land cover, including soil, snow, or water, can affect climate [7].
In recent decades, wetlands have been heavily impacted by anthropogenic activities and climate change [5], lack of an integrated management program for wetlands and the rivers that feed them, droughts, and other factors cause them to dry up.
Land Surface Temperature (LST) as an essential key parameter of the land environment [8], and surface LST in micro, local, regional and global scales [9] is affected by changes in land surface characteristics such as vegetation cover and land use/cover [10]. Therefore, a decrease in the wetland area or reduction of the amount of water in the wetland may be associated with changes in LST. Hence, the study of changes in LST, which plays an important role in the interaction between the surface and the atmosphere [11], is important, especially for aquatic ecosystems in arid and semi-arid regions, which face various natural phenomena such as climate change, drought and also anthropogenic influences. Remote sensing-based LST is a useful tool for studying land surface changes in spatial and temporal terms. Despite the limited ground data, remote sensing techniques are used significantly to understand the changes in land surface properties and also to identify the behavior of climate variables due to the provision of extensive information at high resolution [12]. LST is used in various studies for different scientific purposes, e.g., vegetation and land use/cover changes [13][14][15][16][17][18][19][20][21][22][23], the urban heat island [24][25][26], drought monitoring [27][28][29][30][31], frost risk [32], land degradation [33,34], hydrological purposes [35,36].
Among the various remote sensing indices, NDWI (Normalized Difference Water Index) is one of the most water practical indices used in many studies to track and predict wetland hydrological changes [37][38][39][40]. In addition, NDVI (Normalized Difference Vegetation Index) has been used as a vegetation indicator to investigate changes in vegetation conditions [33,41].
The Gavkhouni wetland located in the eastern part of the Zayandeh-Rud River Basin is fed by the Zayandeh-Rud River. In recent years, the wetland has dried up due to rainfall deficits, severe water shortages, human intervention and manipulation [42], and the transfer of chemical and physical pollutants, which have affected the conditions of the surrounding area and reduced the wetland's services (https://www.iranaqua.ir, accessed on 1 August 2021). Due to the lack of in-situ meteorological data to study the air and land surface temperature changes, remote sensing techniques are used to understand the overall state of change in the wetland's region. A few studies, especially using remote sensing data, were conducted for this wetland, such as Mirahsani et al. [43], who focused on drought analysis using SPI (The Standardized Precipitation Index) and VCI (The Vegetation Condition Index) indices, and Shiran et al. [44], who showed the relationship between NDVI and LST and land properties.
A comprehensive understanding of wetland ecosystem services, especially from the ecological point of view, requires knowledge of the interaction between wetland and climatic parameters of the region. Hence, research is needed on the impacts of the Gavkhouni wetland dryness on regional climate variability to make appropriate river and wetland management decisions in this region. Therefore, the main objective of the current study is quantifying the drying out of the Gavkhouni wetland and its impact on air temperature variability over the study area.

Study Area
The main feature of the Zayandeh-Rud River Basin, which is located in the most central and arid region of Iran, is the Zayandeh-Rud River (Figure 1a), which originates from the Zagros Mountains in the west at an altitude of 4500 m (32 • 22 N, 50 • 04 E) and flows 350 km to the east to end in the shallow place called Gavkhouni wetland at 1470 m [45]. The study area shown in Figure 1a is located in the eastern part of the Zayandeh-Rud River Basin and focuses on the Gavkhouni wetland and its surrounding area (52 • 39 E, 32 • 25 N). The Gavkhouni wetland, which was included in the international Ramsar Convention as an international wetland in 1975 [42] has an area of 476 km 2 and is situated downstream of the Zayandeh-Rud River. According to the reports of the Iran Meteorological Organization, the maximum depth of stored water in the wetland is about 80 to 150 cm.  The climate of the Zayandeh-Rud River Basin has a dramatic range. The average rainfall is ~1500 mm and is comprised of snow falling in the west of the basin, and ~110 mm of rain in the eastern parts, including the Gavkhouni wetland [42].
The Zayandeh-Rud River Basin is one of the most populous and industrialized watersheds in the country of Iran [45]. In the past several decades, the water and vegetation of the basin have been highly abundant, resulting in the identification of more than 910 plant species, 390 genera, and 83 families, including native plants, with genetic, medicinal, ornamental, edible, and industrial values. These plants also serve as a natural filtration system (https://www.iranaqua.ir, accessed on 1 August 2021). In contrast, the Zayandeh-Rud River has recently been dried up for several years and has turned the wetland into a salt pan, which has negatively affected the climate, environment, life of related organisms and the natural ecosystem of it and the surrounding region. Especially changes in the biological and environmental conditions of the Gavkhouni wetland as the most critical wetland in this arid region, have many different consequences, such as a drastic decline of flora and fauna in the wetland, which has been observed in recent years. Figure 2 qualitatively shows the environmental conditions of the wetland. The climate of the Zayandeh-Rud River Basin has a dramatic range. The average rainfall is~1500 mm and is comprised of snow falling in the west of the basin, and~110 mm of rain in the eastern parts, including the Gavkhouni wetland [42].
The Zayandeh-Rud River Basin is one of the most populous and industrialized watersheds in the country of Iran [45]. In the past several decades, the water and vegetation of the basin have been highly abundant, resulting in the identification of more than 910 plant species, 390 genera, and 83 families, including native plants, with genetic, medicinal, ornamental, edible, and industrial values. These plants also serve as a natural filtration system (https://www.iranaqua.ir, accessed on 1 August 2021). In contrast, the Zayandeh-Rud River has recently been dried up for several years and has turned the wetland into a salt pan, which has negatively affected the climate, environment, life of related organisms and the natural ecosystem of it and the surrounding region. Especially changes in the biological and environmental conditions of the Gavkhouni wetland as the most critical wetland in this arid region, have many different consequences, such as a drastic decline of flora and fauna in the wetland, which has been observed in recent years. Figure 2 qualitatively shows the environmental conditions of the wetland.

Trends of NDWI and NDVI for Presence of Water and Vegetation
The remote sensing data used, and the processing applied for trend analysis were based on the Climate Engine Web Application (http://climateengine.org, accessed on 1 August 2021).
For the trend analysis, the values of two spectral indices (NDWI for the presence of water and NDVI for vegetation) in row 38 and path 163 (row and column of the Landsat image tile) were taken from time-series images of Landsat 4 and Landsat 5 (Thematic Mapper; TM), Landsat 7 (Enhanced Thematic Mapper; ETM+), and Landsat 8 (Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS); OLI/TIRS) at a spatial resolution of 30 m from January 1984 to August 2020.
Initially, the values of the indices were derived, and then the time series trends (linear method), long-term averages for the period 1984-2020, and monthly and seasonal averages were subtracted from the time series to show the overall trends in the drying of wetland and the decrease in vegetation cover in its surrounding area.
As shown in Figure 1b, four points were selected near the Gavkhouni wetland and within the wetland (points W1-W2 and V1-V2 were selected for NDWI and NDVI, respectively). In addition, the monthly time series of river discharge at the Varzaneh hydrometric station (latitude = 32.4236°, longitude = 52.6633° and elevation = 1495 m) as the latest station of the Zayandeh-Rud River before the entrance of the Gavkhouni wetland in the period of 1970 to 2019 were obtained from the Regional Water Company of Esfahan (https://www.esrw.ir, accessed on 1 August 2021).

Calculation of LST
The data and the processes used to calculate LST were based on the Earth Explorer application (EE) and an open-source software of remote sensing (QGIS). EE is a web application that allows users to retrieve Earth imagery across all available geospatial data types published by the U.S. Geological Survey (USGS) (https://earthexplorer.usgs.gov, accessed on 1 August 2021). According to the overall results of the wetland drying trends (Section 2.2.1), 16 days in four seasons were selected in the last decade, i.e., 2009, 2011, 2015, and 2017 (Table 1), when the wetland was almost dry. In this section, 16 Landsat images were collected from Landsat 5 (TM) and Landsat 8 (OLI/TIRS) with a spatial resolution of 30 m and cloud cover of less than 10% to investigate the spatio-temporal variability of LST and to show the climate conditions of the studied area. LST variable was calculated according to Equation (1) [46] as follows:

Trends of NDWI and NDVI for Presence of Water and Vegetation
The remote sensing data used, and the processing applied for trend analysis were based on the Climate Engine Web Application (http://climateengine.org, accessed on 1 August 2021).
For the trend analysis, the values of two spectral indices (NDWI for the presence of water and NDVI for vegetation) in row 38 and path 163 (row and column of the Landsat image tile) were taken from time-series images of Landsat 4 and Landsat 5 (Thematic Mapper; TM), Landsat 7 (Enhanced Thematic Mapper; ETM+), and Landsat 8 (Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS); OLI/TIRS) at a spatial resolution of 30 m from January 1984 to August 2020.
Initially, the values of the indices were derived, and then the time series trends (linear method), long-term averages for the period 1984-2020, and monthly and seasonal averages were subtracted from the time series to show the overall trends in the drying of wetland and the decrease in vegetation cover in its surrounding area.
As shown in Figure 1b, four points were selected near the Gavkhouni wetland and within the wetland (points W1-W2 and V1-V2 were selected for NDWI and NDVI, respectively). In addition, the monthly time series of river discharge at the Varzaneh hydrometric station (latitude = 32.4236 • , longitude = 52.6633 • and elevation = 1495 m) as the latest station of the Zayandeh-Rud River before the entrance of the Gavkhouni wetland in the period of 1970 to 2019 were obtained from the Regional Water Company of Esfahan (https://www.esrw.ir, accessed on 1 August 2021).

Calculation of LST
The data and the processes used to calculate LST were based on the Earth Explorer application (EE) and an open-source software of remote sensing (QGIS). EE is a web application that allows users to retrieve Earth imagery across all available geospatial data types published by the U.S. Geological Survey (USGS) (https://earthexplorer.usgs.gov, accessed on 1 August 2021). According to the overall results of the wetland drying trends (Section 2.2.1), 16 days in four seasons were selected in the last decade, i.e., 2009, 2011, 2015, and 2017 (Table 1), when the wetland was almost dry. In this section, 16 Landsat images were collected from Landsat 5 (TM) and Landsat 8 (OLI/TIRS) with a spatial resolution of 30 m and cloud cover of less than 10% to investigate the spatio-temporal variability of LST and to show the climate conditions of the studied area. LST variable was calculated according to Equation (1) [46] as follows: where T s is the LST in degree Celsius ( • C), TB is at-sensor brightness temperature (degree Kelvin, • K, Equation (2)), λ is the wavelength of emitted radiance for which the peak response and the average of the limiting wavelength [47] will be used, ε λ is the land surface emissivity calculated from NDVI and Pv as the percent of vegetation calculated using NDVI, NDVImin, and NDVImax, ρ = h × c/σ, where h is Plank's constant (6.626 × 10 −34 J s), and c is the velocity of light (2.998 × 10 8 m/s), and σ is Stefan-Boltzmann's constant (1.38 × 10 −23 J/K) [48]: In the above Equation (2), K 1 and K 2 are the band-specific thermal conversion constants obtained from metadata in W/(m 2 sr), and L λ is spectral at-sensor radiance in W/(m 2 sr µm) [13].

Calculation of Air Temperature
Air temperature (Tair_calc) variability in the region was derived by a linear regression model between the observed air temperature (Tair_obs) at two available meteorological stations (Varzaneh and Hasanabad, Figure 1b) and the remotely sensed LST product at these locations. The coefficient of determination (R 2 ) and root mean square error (RMSE) were used to evaluate the accuracy of the model developed for estimating air temperature based on Landsat LST products for the studied area. The equations of these metrics are as follows as Equations (3) and (4) where X 0 is the observed data, and X s is the simulated data.

Trend of Observed Temperatures
To understand the effects of the dryness of the Gavkhouni wetland on the temperature of its surrounding geographical areas, according to the maps of LST involving 16 cases, two summer cases were selected in which the wetland was dry and almost wet, respectively. According to the reports of the Iran Meteorological Organization, the prevailing wind direction in Varzaneh is west in winter, spring, and autumn, and east in summer. Therefore, two cases were considered in summertime to obtain an overview of the effect of the wetland dryness on its surroundings, especially in the west toward the northwest parts. Next, the observed air temperatures of the meteorological stations (Varzaneh and Hasanabad) at 06:00 and 12:00 GMT were compared with the long-term seasonal average of observed air temperatures (Tair_obs_m) ( Table 2 and Figure 1b). Then, the trends of the observational mean monthly minimum (Tmin_m) and maximum air temperature values (Tmax_m) were calculated using the non-parametric Mann-Kendall test [49]. The observed monthly air temperature data were collected for the whole of Iran in 40 meteorological stations (including 38 synoptic stations and two combined synoptic-climatology stations), for Isfahan province (four stations including Isfahan, Kashan, Naein, and Varzaneh) and for Varzaneh station, as three areas from 1970 to 2019. Then the data was divided into 1970-1999 and 2000-2019 as two periods before and after 2000 ((based on the river discharge (m 3 /s) at Varzaneh station, and the trends of the NDWI and NDVI indices (Section 2.2.1). Finally, the mean value of Sen's slope (the slope in Mann-Kendall trend test), for all 12 months for three areas was derived and compared.

Comparison of Observational Air Temperature in Dryness and Wetness Conditions of Wetland
To identify the difference in air temperature adjacent to the wetland area in two different wet and dry conditions of the wetland, groups of the observational air temperatures in dry and wet conditions were used. To find out the dry and wet conditions, the monthly values of river discharge (m 3 /s) at the hydrometric station closest to the Gavkhouni wetland from 1970 to 2020, and the results of the presence of water and vegetation cover (Section 2.2.1) were considered. For this purpose, all available Landsat imagery from 1982 to 2020 was also used to determine the years in which the wetland was full of water or was almost dry. Next, the mean monthly air temperatures of these two situations were compared using the Student's t-test (t test). To eliminate the role of global warming during recent decades (2000s and 2010s), the values of the global warming rate for each season and both the minimum and maximum temperatures were derived using the results of the temperature increase rate or Sen's slope estimator (Section 2.2.4). In other words, to find out the effect of the wetland dryness on air temperature variations in the region, the global warming rate for all data was removed, and the remaining values were used in a paired two-sample t test. Figure 3 shows the time series analysis of NDWI and NDVI values at the selected points (W1; in the river and near the entrance of the wetland (Figure 3a), W2; in the middle of the wetland (Figure 3b), V1 and V2; near the river and the entrance of the wetland, respectively (Figure 3c,d)) between January 1984 and August. In Figure 3, a significant downward trend can be seen at all points, but the long-term average value of NWDI was negative and close to zero, corresponding to the non-water area. While the long-term value of NDVI was positive and close to 0.3, corresponding to shrublands and grassland. In addition, the highest values of NDWI (~0.3, corresponding to the water feature) and NDVI (≥0.4) were observed in the 1990s and 1980s (especially in 1988(especially in , 1989(especially in , 1993(especially in , and 1994(especially in for NDWI, and 1997(especially in , 1993(especially in , and 1994   warming rate for all data was removed, and the remaining values were used in a paired two-sample t test. Figure 3 shows the time series analysis of NDWI and NDVI values at the selected points (W1; in the river and near the entrance of the wetland (Figure 3a), W2; in the middle of the wetland (Figure 3b), V1 and V2; near the river and the entrance of the wetland, respectively (Figure 3c,d)) between January 1984 and August. In Figure 3, a significant downward trend can be seen at all points, but the long-term average value of NWDI was negative and close to zero, corresponding to the non-water area. While the long-term value of NDVI was positive and close to 0.3, corresponding to shrublands and grassland. In addition, the highest values of NDWI (~0.3, corresponding to the water feature) and NDVI (≥0.4) were observed in the 1990s and 1980s (especially in 1988(especially in , 1989(especially in , 1993(especially in , and 1994(especially in for NDWI, and 1997(especially in , 1993(especially in , and 1994   In the most recent years (especially from 2008/2009 to 2018), the values of NDWI and NDVI were below their long-term averages. Figure 4 illustrates the seasonal and monthly averages of the NDWI and NDVI indices over the long term at four points. Figure 4a shows that the highest values of the seasonal average of NDWI at point W1 were in autumn, followed by summer, winter, and spring. At point W2 the highest values were observed in winter, spring, summer, and autumn because the Zayandeh-Rud River was dry for a long time and no water flowed to fill the Gavkhouni wetland. While the highest values of NDVI at both points V1 and V2 were in summer, followed by autumn, spring, and winter. In the most recent years (especially from 2008/2009 to 2018), the values of NDWI and NDVI were below their long-term averages. Figure 4 illustrates the seasonal and monthly averages of the NDWI and NDVI indices over the long term at four points. Figure 4a shows that the highest values of the seasonal average of NDWI at point W1 were in autumn, followed by summer, winter, and spring. At point W2 the highest values were observed in winter, spring, summer, and autumn because the Zayandeh-Rud River was dry for a long time and no water flowed to fill the Gavkhouni wetland. While the highest values of NDVI at both points V1 and V2 were in summer, followed by autumn, spring, and winter.   Figure 4b shows the monthly average values of the NDWI and NDVI indices over the long-term at four selected points, where the NDWI at point W1 had high values in March, June, and July, while it was observed at point W2 in September and October. The high values of NDVI at both points V1 and V2 were related to June to November (summer and autumn months). As can be seen in Figure 4b, the changes at point W2 were similar to those at points V1 and V2, which means that the wetland (W2) was dry during the summer months (negative values of NDWI). Since the Zayandeh-Rud River (which fed the wetland) was dry, the vegetation covers near the river and the wetland decreased almost with a time lag. Figure 5 shows the monthly discharge values at Varzaneh as the nearest hydrometric station to the Gavkhouni wetland, from 1970 to 2019. As shown in Figure 5a, the highest discharge values were recorded in spring (especially May, and then April) and winter months (especially December). The highest discharge values were related to 1976, then 1993, 1987, 1988, and 1980 with maximum values of 107, 90.5, 82.5, 64, and 63 (m 3 /s), respectively. Figure 5a,b shows that most years before 2000 were wet, while the wetland condition after that year was dry.   Figure 4b shows the monthly average values of the NDWI and NDVI indices over the long-term at four selected points, where the NDWI at point W1 had high values in March, June, and July, while it was observed at point W2 in September and October. The high values of NDVI at both points V1 and V2 were related to June to November (summer and autumn months). As can be seen in Figure 4b, the changes at point W2 were similar to those at points V1 and V2, which means that the wetland (W2) was dry during the summer months (negative values of NDWI). Since the Zayandeh-Rud River (which fed the wetland) was dry, the vegetation covers near the river and the wetland decreased almost with a time lag. Figure 5 shows the monthly discharge values at Varzaneh as the nearest hydrometric station to the Gavkhouni wetland, from 1970 to 2019. As shown in Figure 5a, the highest discharge values were recorded in spring (especially May, and then April) and winter months (especially December). The highest discharge values were related to 1976, then 1993, 1987, 1988, and 1980 with maximum values of 107, 90.5, 82.5, 64, and 63 (m 3 /s), respectively. Figure 5a,b shows that most years before 2000 were wet, while the wetland condition after that year was dry.   Figure 4b shows the monthly average values of the NDWI and NDVI indices over the long-term at four selected points, where the NDWI at point W1 had high values in March, June, and July, while it was observed at point W2 in September and October. The high values of NDVI at both points V1 and V2 were related to June to November (summer and autumn months). As can be seen in Figure 4b, the changes at point W2 were similar to those at points V1 and V2, which means that the wetland (W2) was dry during the summer months (negative values of NDWI). Since the Zayandeh-Rud River (which fed the wetland) was dry, the vegetation covers near the river and the wetland decreased almost with a time lag. Figure 5 shows the monthly discharge values at Varzaneh as the nearest hydrometric station to the Gavkhouni wetland, from 1970 to 2019. As shown in Figure 5a, the highest discharge values were recorded in spring (especially May, and then April) and winter months (especially December). The highest discharge values were related to 1976, then 1993, 1987, 1988, and 1980 with maximum values of 107, 90.5, 82.5, 64, and 63 (m 3 /s), respectively. Figure 5a,b shows that most years before 2000 were wet, while the wetland condition after that year was dry.

LST and Tair_calc Spatiotemporal Variations
To investigate the spatial and temporal variations of land surface and air temperatures with a spatial resolution of 30 m over the Gavkhouni wetland and its surrounding area, as Figures 6 (as an example for summer season) and S1-S3 demonstrate the Landsat LST ( • C) at 06:48-7:03 GMT (10:18-10:33 local time) on 16 selected days in summer, winter, spring, and autumn, respectively. In all seasons and on almost all selected days, LST in the wetland was lower than in the surrounding area, also at Varzaneh and Hasanabad stations. In the study area, during winter, spring, summer, and autumn, the areas with the highest values of LST were most prevalent in the outer areas near the wetland corresponding to the sand dunes (as an example shown in Figure 6(e1-1)) and black rock (Figure 6(e1-2)). Some parts of the Gavkhouni wetland, mainly after the entrance, and narrow parts in the east and southeast of the wetland were cooler than other parts, while some parts, especially in the middle and south of the wetland, had a higher temperature class. The minimum values of LST were also observed outside the wetland, for example, some parts in the west of the study area (white rocks in the mountainous part ( Figure 6(e1-3))) and some narrow parts near Varzaneh (in the northwest of the wetland) and Hasanabad towns (west) (built-up area ( Figure 6(e1-4))) which had high albedo and ultimately decreased the LST values.

LST and Tair_calc Spatiotemporal Variations
To investigate the spatial and temporal variations of land surface and air temperatures with a spatial resolution of 30 m over the Gavkhouni wetland and its surrounding area, as per Figures 3-5, some cases were selected in 2009, 2011, 2015 and 2017, that approximately corresponded to the dry condition of the wetland and the available and usable Landsat imagery. The maps LST and the subsequently calculated air temperature (Tair_calc) are presented and discussed below.
Figures 6 (as an example for summer season) and S1-S3 demonstrate the Landsat LST (°C) at 06:48-7:03 GMT (10:18-10:33 local time) on 16 selected days in summer, winter, spring, and autumn, respectively. In all seasons and on almost all selected days, LST in the wetland was lower than in the surrounding area, also at Varzaneh and Hasanabad stations. In the study area, during winter, spring, summer, and autumn, the areas with the highest values of LST were most prevalent in the outer areas near the wetland corresponding to the sand dunes (as an example shown in Figure 6(e1-1)) and black rock (Figure 6(e1-2)). Some parts of the Gavkhouni wetland, mainly after the entrance, and narrow parts in the east and southeast of the wetland were cooler than other parts, while some parts, especially in the middle and south of the wetland, had a higher temperature class.
The minimum values of LST were also observed outside the wetland, for example, some parts in the west of the study area (white rocks in the mountainous part ( Figure 6(e1-3))) and some narrow parts near Varzaneh (in the northwest of the wetland) and Hasanabad towns (west) (built-up area ( Figure 6(e1-4))) which had high albedo and ultimately decreased the LST values. As shown in Figure S4, the calculated air temperature (Tair_calc) was obtained using the linear regression between LST and observed air temperature (Tair_obs) for the selected 16 days at two stations (Varzaneh and Hasanabad) which showed the strong correlation (r = 0.97, r 2 = 0.94) and RMSE = 3 • C. The difference between LST and Tair_obs in 16 days showed higher values in Varzaneh than the Hasanabad station.
At these two stations, and in all seasons, the values of LST were higher than those of Tair_calc and Tair_obs. Moreover, the values of Tair_calc and Tair_obs were close in all seasons except winter, and the RMSE was about 3 • C at two stations.
As shown in Figure S5, the values of LST at Hasanabad were lower than that of Varzaneh station on most of the selected days and months. Figures 6 and S1-S3 also show that in this study area, as with LST (Figures 6 and S1-S3), the lowest and highest Tair_calc values were found in almost all the seasons inside and outside the wetland, respectively.

Two Summer Cases
For this purpose, according to Table 2 and Figures 6 and S1-S3, the days 27 June 2015 and 28 July 2009 were selected ( Figure 6(a1,c1)) where the maximum and minimum values of LST occurred at 06:00 GMT for the entire wetland and for the entire study area in summer, respectively. These cases also had a wet and dry winter, respectively. Table 2 and Figure S5 show that in summer at 06:00 GMT, when the wetland was not completely dry (27 June 2015), the Tair_obs value of Hasanabad was cooler than its Tair_obs_m (30 • C) and Varzaneh with differences of 0.7 • C and 2.1 • C, respectively. On these days, there was no wind in Varzaneh. When the wetland was almost completely dry (28 July 2009), the Tair_obs value of Hasanabad (32.1 • C) was warmer than that of Tair_obs_m and Varzaneh stations with a difference of 2.1 • C.
In summer at noontime, when the wetland was not completely dry (27 June 2015), the Tair_obs value of Hasanabad was slightly warmer than the Tair_obs_m value (34.5 • C) with a difference of 0.1 • C and also cooler than Varzaneh (35.8 • C). While the wetland was completely dry (28 July 2009), the Tair_obs value of Hasanabad was 2.5 • C higher than the Tair_obs_m value and lower than Varzaneh (38.2 • C), and the air temperature of Varzaneh was 2.1 • C warmer than the Tair_obs_m value. On 28 July 2009 and 27 June 2015, at Varzaneh, there were southeasterly and southerly winds at noon (130 and 170 degrees in dry and wet conditions, respectively). Due to the small difference in altitude between Varzaneh and Hasanabad, the wind directions of Hasanabad were considered the same as those of Varzaneh station.  Figure 7a,b in all three areas, the slope of Tmax_m was positive in half of the months (December, January, February, March, May, and October), while the slope of Tmin_m was positive in all months, and the slope was more significant in Varzaneh than in the province and Iran in most months, i.e., February and April to December with a relatively significant difference (0.02 to 0.06 • C/year). At Varzaneh station, the maximum slope values of Tmax_m were observed in January (0.07), and then in February (0.06), the maximum slope values of Tmin_m were related to October (0.1), September and July (0.09) (autumn and summer months). In the years before 2000, the rates of change of Tmax_m in all three areas were positive only in December, January, and February (generally winter months) and negative in most months. Furthermore, the slope of this temperature was more negative in most months in Isfahan province than in Varzaneh and the country. The slope values of Tmin_m were positive in all three areas in six months (August, September, October, November, December, and January) and negative in spring and almost summer months in the province of Isfahan (Figure 7c,d). At Varzaneh station, the slope values of Tmin_m in the cold months of the year such as October and December were also higher than those in the province and the country by the difference (0.06, 0.04) and (0.03 and 0.01), respectively (Figure 7c,d), while the slope values in December, January, and November (during the winter) were negative in the post-2000 period (Figure 7e,f).

Mann-Kendall Trend Test Results
According to Figure 7e,f, in the 20 years after 2000, the slope values of Tmax_m were positive up to half of the months in three areas (generally summer and winter months: December, January, February, June, July, and September), and the rates of Tmin_m changes in three areas, in the five months of March, May, June, September, and October were positive, with the magnitude of the slope in Varzaneh being larger than in Isfahan province and Iran only in June (0.02 and 0.06, respectively). Furthermore, the greatest magnitude difference in slope values of Tmax_m between Varzaneh with the province In the years before 2000, the rates of change of Tmax_m in all three areas were positive only in December, January, and February (generally winter months) and negative in most months. Furthermore, the slope of this temperature was more negative in most months in Isfahan province than in Varzaneh and the country. The slope values of Tmin_m were positive in all three areas in six months (August, September, October, November, December, and January) and negative in spring and almost summer months in the province of Isfahan (Figure 7c,d). At Varzaneh station, the slope values of Tmin_m in the cold months of the year such as October and December were also higher than those in the province and the country by the difference (0.06, 0.04) and (0.03 and 0.01), respectively (Figure 7c,d), while the slope values in December, January, and November (during the winter) were negative in the post-2000 period (Figure 7e,f).
According to Figure 7e,f, in the 20 years after 2000, the slope values of Tmax_m were positive up to half of the months in three areas (generally summer and winter months: December, January, February, June, July, and September), and the rates of Tmin_m changes in three areas, in the five months of March, May, June, September, and October were positive, with the magnitude of the slope in Varzaneh being larger than in Isfahan province and Iran only in June (0.02 and 0.06, respectively). Furthermore, the greatest magnitude difference in slope values of Tmax_m between Varzaneh with the province and the country were associated with April (0.03) and March (0.04)-December (0.03), respectively. The greatest magnitude difference in slope values of Tmin_m of Varzaneh with the province and the country were associated with June with values of 0.02 and 0.05, respectively.

t Test Results
For performing the t test between two groups, the years 1976, 1977, 1978, 1979, 1980, 1987, 1988, 1989, 1993, and 1994 were selected for the wet group, and 2000, 2001, 2002, 2009, 2010, 2011, 2015, 2016, 2017, and 2018 were chosen for the dry group, that showed the Gavkhouni wetland was wet before the year 2000 and has dried up in the recent decades. Table 3 shows the seasonal mean and t-values of Tmin_m and Tmax_m obtained from the two-sample t test at Varzaneh station. The difference between the seasonal mean value of Tmax_m when the wetland was dry and when it had water was significant only in spring (March, April, and May) with a value of 1.4 • C (p < 0.05 in 2-tailed significance). However, when the wetland was dry, the seasonal mean Tmin_m was significantly increased in most seasons, i.e., winter (December, January, and February), spring (March, April, May), and summer (June, July, and August) with values of 2.64 • C, 1.61 • C, and 1 • C, respectively.

Discussion
This study used observed data of temperature and river discharge in combination with remote sensing-based indices and LST to reveal the trend of the Gavkhouni wetland dryness and its impact on the temperature of the surrounding areas. The results show that the drying of the wetland over the last two decades has led to a significant increase in air temperatures in the surrounding areas, especially in the west and northwest in spring and summer seasons.
Since 2000, the smallest wetland extent was observed in 2009, followed by 2011, 2017 and 2015. During all seasons (Figures 6 and S1-S3), LST and air temperature in the Gavkhouni wetland were lower than for the surrounding area, which can be explained by the high albedo of the salt deposits which are exposed during dry periods. In general, water has a low albedo and thus absorbs more incident visible and near-infrared solar radiance Hence, the albedo increases due to the dryness of the wetland and the salt deposited there, leading to a decrease of LST. Therefore, the wetland cannot balance its ambient temperatures. The values of LST and Tair_calc in the selected years show that the year 2009 was one of the driest years not only in Iran but also in the entire Middle East [50]. The results also indicated a good correlation between LST and Tair_obs (r 2 = 094, r = 0.97) at the stations near the Gavkhouni wetland ( Figures S4 and S5). Jafari and Hasheminasab [34] and Shiran et al. [44] also showed that LST variable is a good tool for monitoring arid environments. In addition, the RMSE was by an average of 3 • C, but in general Tair_calc values were consistent with changes in LST for most days and most parts of the studied area. The long-term seasonal average of observed air temperature (Tair_obs_m) at Varzaneh and Hasanabad stations at 06:00 and 12:00 GMT (Table 2) showed that at Hasanabad at 06:00 Tair_obs_m was higher than Varzaneh in all seasons except summer, while at noon it was cooler than Varzaneh station in all seasons. The results of the Mann-Kendall (Sens's slope estimator) from 1970 to 2019, in general, showed that the slope values of Tmax_m changes before the year 2000 were often negative, while after the year 2000 were positive in half months. The slope values of Tmin_m changes in Varzaneh in the period before the year 2000 were higher in the cold months of the year than in the province and the country, while in the last 20 years were lower than the province and the country, which is related to the dryness of the wetland and its failure to play a role in regulating the surrounding temperature. In addition, in the warm months of the year, in the 30 years before and 20 years after 2000, the slope of Tmin_m changes in Varzaneh station was less and more, respectively, than both the province and country. In the 50 years, the maximum of the slope values of Tmin_m changes between the three areas were observed in Varzaneh station and months of October (0.1), September (0.09), and July (0.09), also in this period, the maximum difference of the slope of Tmin_m between Varzaneh and the country (0.05 • C/year) was related to July (month of summer). On 28 July 2009, when the wetland was in dry conditions, the Tair_obs values at Hasanabad and Varzaneh stations at noontime had a difference of 2.5 • C and 2.1 • C, respectively, from the long-term average values. To confirm that the wetland dryness affects the climate in the surrounding areas, a t test analysis was used to compare the region's temperature for the dry and wet wetland years. Considering the significant results of the t test for the average of the minimum air temperature in Varzaneh in most seasons, especially in spring and summer (Table 3), and also the Sen's slope results, it can be said that in recent decades with the dryness of the Gavkhouni wetland, the average minimum air temperature in spring and summer has increased by about 1.6 • C and 1 • C, respectively. Therefore, from the analysis of the increase in air temperature at Varzaneh and Hasanabad stations on 28 July 2009, It can be concluded that the increase in temperature on this day at noontime was caused not only by global warming but also by the dryness of the Gavkhouni wetland.
In general, the Gavkhouni wetland has been experienced dryness conditions mostly since the 2000s and 2010s decades. According to the NDVI and NDWI results, the negative trend from 1984 to 2020 indicates that the Zayandeh-Rud River has experienced dryness during this period that dried up the Gavkhouni wetland which was a permanent wetland, and the vegetation near the wetland and its entrance decreased. These decreasing trends show that the spectral indices were almost consistent with the other evidence and observed data such as the time series of discharge values at Varzaneh station (Figure 5a,b) which clearly showed the wetland dryness. Besides the drought occurrences, resulted from climate change impact, the human activities in the basin have influenced the drying of the Gavkhouni wetland. Abou Zaki et al. [51] also reported that even in arid and semi-arid climates, water scarcity is mainly due to anthropogenic reasons rather than climate. As presented in other studies, the Gavkhouni wetland dryness is related to population growth and bad management of water resources, such as upstream diversions and unreasonable transmissions of water resources in the basin [42], expansion of agricultural land and irrigation patterns causing hydrological droughts [51], dam construction in upstream of the Zayandeh-Rud River [34]. Related to the remote sensing indices, their accuracy has been confirmed in other studies, such as Eckert et al. [33] and Shiran et al. [44] which indicated that the analysis of NDVI time series can be useful and appropriate to understand the changes in vegetation cover, and also in the study of Jiang et al. [41] which showed the NDVI can accurately indicate the vegetation condition. Abou Zaki et al. [51] also reported that NDWI is a reliable tool for studying the area of surface water bodies.

Conclusions
This study was conducted to determine the effects of the Gavkhouni wetland dryness on temperature variations in the surrounding areas. In the last two decades, in addition to global warming impact (such as droughts caused by higher temperatures and evapotranspiration), wetland dryness caused by human activities has led to an increase in air and land surface temperatures in spring and summer, especially in the western and northwestern parts. From 1970 to 2019, the increase in air temperatures in the region was higher than in the Isfahan province and Iran in general. During the warm months, especially in July, this rate was about 0.02 to 0.06 • C/year higher than the average for Isfahan province and the country. From the current study, it can be concluded that monitoring changes in surface and air temperatures, as well as studying the drying trend of water bodies and their relationship with climate variability in surrounding areas, can be used to develop an integrated plan for managing wetlands and the river that feeds them. The overall results of this study provide a basis for making effective decisions about allocating specific amounts of water to the Gavkhouni wetland to maintain its climate and environmental functions. Considering that many factors such as latitude, cloud cover and direction of prevailing winds affect land surface and air temperatures, it is proposed to use a numerical climate model to improve regional understanding of the impact of wetland dryness on the surrounding climate.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.