Interactive Changes in Climatic and Hydrological Droughts, Water Quality, and Land Use/Cover of Tajan Watershed, Northern Iran

: In response to novel and complex uncertainties, the present research is conducted to characterize the most signi ﬁ cant indicators of watershed health including drought, water quality, and vegetation for the Tajan watershed, Mazandaran, Iran. The Standardized Precipitation Index (SPI) and Stream ﬂ ow Drought Index (SDI) are, respectively, used to quantify the meteorological and hydrological droughts in the present (1993–2020) and future (2023–2050) employing optimistic RCP2.6 and pessimistic RCP8.5 scenarios. To concoct discharge data for the future, IHACRES v1.0 software is used with a Nash–Sutcli ﬀ e coe ﬃ cient (NSE) of 0.48 and a coe ﬃ cient of determination (R 2 ) of 0.58. Maps of land use and Normalized Di ﬀ erence Vegetation Index (NDVI) are also prepared using Landsat images. Subsequently, the surface water quality is assessed using AqQA v1.1.0 software. The results show the di ﬀ erence in the severity of future meteorological droughts in di ﬀ erent stations. In addition, the predominance of non-drought (SDI ≥ 0) or mild drought ( − 1 ≤ SDI < 0) is indicated for future hydrology. The land use changes show a decrease in rangeland ( − 5.47%) and an increase in residential land (9.17%). The water quality analysis also indicates an increase in carbonate ions in the watershed outlet. Communicating the relationships between study indicators, which is a big gap in the current watershed management approach, avoids future failures and catastrophes.


Introduction
The proper functioning of a watershed influences the appropriate and principled access to ecosystem services, characteristics, functions, or ecological processes that contribute to human well-being [1,2].Assessing the watershed's ability to perform expected functions depends on overall watershed health.However, in recent decades, various natural and human-induced driving forces such as climate change, drought, flood, deforestation, population growth, overgrazing, erosion, and sedimentation have significantly impacted the health of watersheds [2,3].
Climate change and human activities, specifically changes in land use/cover, are widely recognized as the most significant factors impacting watershed health [4,5].The destruction of natural resources and changes in land use, along with industrial growth and development, has led to an increase in greenhouse gas production [6].Due to the rise in greenhouse gases, the climate in various regions is changing, leading to global warming.According to the Intergovernmental Panel on Climate Change (IPCC), the Earth's average surface temperature has increased by approximately 0.6 °C in the last century.It is projected that this increase will range from 1.4 to 5.7 °C by the year 2100.It is anticipated that these changes in temperature and precipitation will have a significant impact on the planet's ecosystems [7].
Iran has not been immune to these changes.In recent years, it has experienced numerous natural hazards, such as floods, droughts, and pests, due to changes in climatic variables [8][9][10][11].The frequency and severity of these hazards are on the rise.Climate plays a significant role in various natural hazards, including drought, which can have a profound impact on the economy and people's lives.This is particularly true in sectors such as agriculture, animal husbandry, and industry, where it can cause damage and destruction in different regions.Iran is currently facing one of its biggest climatic challenges in recent years: drought [12,13].
"A drought is a period of abnormally dry weather characterized by a prolonged deficiency of precipitation below a certain threshold over a large area and a period longer than a month" [14].Drought is also classified according to the affected domain as meteorological, agricultural, hydrological, and socioeconomic.Each individual component, as well as the overall total, has an impact on various aspects of the watershed.For instance, climatic changes have a significant impact on the temporal and spatial distribution of precipitation in various regions, leading to numerous problems.One of the most significant differences between drought and other natural disasters is its gradual process and slow speed.It may take several years or even longer for the effects to become evident, making it different from other hazards [15].It is a sluggish, risky phenomenon, so the mental impacts on living organisms from drought are more delicate and prolonged than other disasters [16].However, the applicability of different time scales in the assessment of drought conditions is important since each can capture something different, which is highly relevant to, for example, agriculture or hydrology.Recent research has highlighted that droughts are increasing over shorter time scales, but their overall duration is shortening, indicating that there is a sustained short-term climatological trend in drought [17,18].Longer time scales for drought identification show a reduced occurrence of drought, but each occurrence is prolonged [17,18].
Land use change, as another dynamic component of an ecosystem, has been affected by human activities, and along with the phenomenon of climate change, it has increased the severity of risks and environmental disturbances [19][20][21][22][23].As a result of the increase in population in the past years and the increase in demand for resources, food, and energy, a large area of wilderness forest and rangelands was subjected to excessive pressure and destruction, and the result of these changes is clearly visible in the occurrence of today's dangers [23][24][25][26].The most important effects of land use change that are currently occurring as human threats and reducing the health of the ecosystem are the phenomena of dust, increase in runoff and flooding, decrease in oxygen and increase in CO2, decrease in biodiversity, etc. [24].Furthermore, the impact of climate change and land use/ cover causes changes in other watershed components such as runoff and water quality [5,24,[27][28][29][30][31][32][33].The land use and vegetation in an area can directly impact the water quality and overall health of the watershed [34][35][36][37][38][39].Therefore, as it is known, the various components in a watershed are related to each other in a complex way, and changes in each of these components cause changes in other components.
The literature also depicts the interactive impact of watershed health indicators in different conditions.For example, Mishra et al. [25] studied the potential effects of extreme weather-related disasters (e.g., drought, floods, and wildfires) on various water quality in the Carolinas, the United States of America.The results showed that flooding leads to high levels of pollution-related bacteria, especially in the river environment.The impact of drought on water quality indicators in different land use environments can be different, which highlights the dominant role of watershed characteristics.In addition, in Thailand, Khadka et al. [40] predicted the effects of climate and land use changes on future droughts (2021-2050) in one of the main branches of the Mekong River.Both climate and land use changes had significant effects on drought-sensitive agriculture.In addition, in the Metropolitan region of Chile, drought indices were used to describe meteorological and hydrological droughts for the period 1985 to 2015.To understand the relationship between drought and water quality, the correlation between daily discharge and surface water quality observations was evaluated.A significant negative relationship (p<0.05) was observed between discharge and electrical conductivity and major ions in most of the region.Hydrological stations located in irrigation areas exceed the allowed threshold [41].
Assessing the ecosystem conditions affected by natural and human-induced hazards, such as drought and land use changes, is particularly important at the watershed level and over different time periods to facilitate effective decision-making and management.The Tajan watershed, a crucial ecosystem in Mazandaran province, faces destruction and decline due to various factors [42].These factors include severe land use changes driven by increasing land prices, deforestation, reduced water resources caused by climate change, increased damages from floods and droughts, changes in surface water quality due to land use changes, and the use of agricultural pesticides and industrial effluents [43].Therefore, the study of interactive changes in the most impact-causing driving forces of watershed health degradation in this watershed is highly crucial, but these driving forces have not been considered yet.This research also tries to fill the gaps related to the prediction of climatic and hydrological drought simultaneously in different scenarios and its relationship with land use and the greenness index.Our findings suggest various scenarios to address and adapt to these changes to mitigate destruction, and the promotion of watershed sustainability could be proposed.

Study Area
This research was conducted in the Tajan watershed (53°18′ to 53°05′ E longitude and 36°29′ to 36°09′ N latitude), which covers an area of 3980 Km 2 (Figure 1).This watershed, located upstream of Sari and Neka cities in Mazandaran province, is one of the largest watersheds in Northern Iran.It is bordered by the Alborz Mountains to the south and the Caspian Sea to the north.The watershed is extensive, spanning across six cities: Neka, Ghaemshahr, Sari, Behshahr, Damghan, and Savadkuh.The average annual temperature of the watershed is approximately 15 °C, having a hot and humid climate.The region receives an average rainfall of 832 mm annually [42].The altitude ranges from −2 m to 3718 m asl, with the outlet closed above the Caspian Sea inlet.The majority of the studied area consists of forest lands, with agricultural lands limited to river banks or sloping areas in certain regions.
The primary occupation of local communities, particularly in the upstream villages, is predominantly agriculture and animal husbandry.Due to favorable weather conditions, the tourism industry has been somewhat active in recent years, leading to changes in land use in various areas of this watershed.In addition to the well-known Tajan river, this area is also home to other major branches such as Chahardange, Dodange, Lajim, and Zarmrud.These rivers flow through different villages before eventually reaching the Caspian Sea.This watershed is selected for the present study because it contains numerous ecosystems and important water sources and rivers that provide drinking water and support agriculture for local and urban communities downstream.

Data Sources and Processing
The essential climatic data (1993-2020) needed for climate modeling, forecasting, and meteorological drought investigation, including precipitation, maximum and minimum temperature, and sunshine hours or radiation, are collected from the Mazandaran Meteorology General Office (https://www.mazmet.ir/en/,accessed on 5 April 2022).To predict future climate patterns, climate data are collected from three synoptic and rain gauge stations (i.e., Kordkhil, Soleyman tangeh, and Rig Cheshmeh), located in or near the study area (Figure 1).These stations were selected based on their suitable distribution and consistent time period.
Lars-WG v6.0 software is utilized to forecast future climate conditions (2023-2050) based on two climate scenarios: optimistic (RCP2.6)and pessimistic (RCP8.5)[44].DrinC v1.7 software is also used to determine drought thresholds based on the SPI standard table, using current and future data.
Additionally, existing hydrometric stations in the region with more complete data are used to investigate the trend of flow discharge changes.The data from the hydrometric stations are obtained from the Mazandaran Regional Water Authority (https://www.mzrw.ir/?l=EN, accessed on 5 April 2022).The time series of the hydrological data spans from 1997 to 2017 for the base period and from 2021 to 2040 for the future period.The IHACRES hydrological model was utilized to forecast future discharge.The trend of future discharge data is predicted using long-term data on precipitation and temperature, along with observed discharge.The hydrological drought at the outlet station (Kordkhil) is then calculated using the SDI in DrinC v1.7 software [15].After analyzing the stream flow data and using the standard table for the SDI, the thresholds for hydrological drought were determined.
AqQA v1.1.0software was used to analyze the trends and characteristics of surface water quality at two stations of Kordkhil and Parvich abad (Figure 1), over 20 years.Finally, the trend of land use/cover changes and their impact on parameters such as water quality and runoff were analyzed using Landsat satellite images .

Climatic Drought Calculation
The Standardized Precipitation Index (SPI) is a powerful tool for analyzing rainfall data [45].This index assigns a numerical value to rainfall, allowing for comparisons between areas with different climates.It also establishes a quantitative index for assessing current drought conditions.The SPI quantifies the lack of precipitation across various time scales.The main feature of this index is its flexibility in measuring droughts across various time scales, as droughts can vary greatly in terms of duration.It is important to detect and monitor these risks using various time scales.The SPI for each region is calculated using long-term rainfall statistics for various time periods.The longterm statistic is fitted to the gamma probability function, which is then used to calculate the cumulative probability of rainfall for a specific station, month, and time scale.When calculating the SPI, the first step is to select a suitable probability distribution function that can accurately represent a long series of observed rainfall data.The Gamma distribution has been widely used in various studies as an effective model for describing rainfall data.The cumulative probability of the observed rainfall data is calculated.The cumulative probability is transformed using the inverse normal Gaussian function with a mean of 0 and a variance of 1 to calculate the SPI for 3, 6, and 12 months [8,12].The gamma distribution is defined by the probability density function (Equation ( 1)).

𝑓(𝑥, 𝑎
where  is the shape and β shows the scale parameters, and these scale parameters can be obtained by maximum likelihood estimation.The standard classification of this index is presented in Table 1.

Hydrological Drought Calculation
The Streamflow Drought Index (SDI) is used to investigate hydrological drought.This index is similar to the SPI in terms of calculations, but it is recommended for use on a monthly and seasonal time scale.The SDI, developed by Nalbantis [47], is based on the standard monthly flow (Equation ( 2)).
This equation assumes that in the time series of monthly river flows, i represents the hydrological year and k represents the corresponding month.Vk and Sk represent the average and standard deviation, respectively, of the total flow volume for the base period k.Vik refers to the cumulative volume of river flow during the base period.The standard classification of SDI is presented in Table 2. [47,48].

Condition Classification SDI Value
Extreme drought SDI < −2

IHACRES Hydrological Model
IHACRES is a conceptual rainfall-runoff model that was initially developed by Jakeman et al. [49].The model has been successfully applied to watersheds with various sizes, ranging from 490 m 2 in China to 10,000 Km 2 in England, and with time intervals ranging from 6 min to 1 month.To implement this model, data inputs such as precipitation, temperature, and flow discharge are required.Temperature and precipitation are used to simulate flow discharge, which is then used to calibrate and evaluate the accuracy of the model [50].
IHACRES comprises two nonlinear reduction modules and a linear hydrograph module.To implement this model, rainfall and temperature are first converted to effective rainfall in each time step using the nonlinear module.Then, the single hydrograph's linear module converts it to surface runoff within the same time step [51].This model consists of six parameters.Three parameters, c (moisture storage capacity of the watershed), tw (the time it takes for the watershed to dry), and f (watershed temperature adjustment factor), are related to the nonlinear loss component of the model.The other three parameters, tq, ts (the time it takes for fast and slow flow to decrease, respectively), and Vs (the volume of slow flow that contributes to river flow), are related to the linear function component of the model [52].In this research, the default value of key parameters of the model was used.The discharge in the selected hydrometric station is predicted using the IHACRES model.This prediction is based on climate forecasts, including precipitation and temperature data, as well as the current discharge.
For the simulation, the data of daily precipitation, average daily temperature, and daily discharge of the statistical period (1997-2017) of Kordkhil station were entered into IHACRES v1.0 software, and a time period of 2 years (730 days) was determined for model warm-up.In addition, the period of 1997-2008 was chosen as the calibration period, and the period of 2009-2017 was chosen as the validation period.
To simulate the daily discharge in the Tajan watershed from 2021 to 2040, we used precipitation, temperature, and discharge data from the Kordkhil station.This station was chosen as the main and outlet station for the watershed due to its suitable data.The data needed to run the rainfall-runoff model and predict future discharge includes future precipitation, future temperature, and the base discharge of the desired station.The precipitation and temperature are modeled in previous stages for the future time.The model was calibrated for the period 1997-2012.The period was determined through a trial-and-error process, and the period with the most successful simulation was ultimately chosen.The model was evaluated in the remaining period.The IHACRES model execution diagram is shown in Figure 2.

Surface Water Quality Assessment
The surface water quality is examined by analyzing data from two hydrometric stations: Kordkhil station (downstream) and Parvich abad station (upstream) (Figure 1).We used a suitable time period (2001-2021) to investigate the trend of quality changes.The data analyzed in this research to determine water quality include electrical conductivity (EC), total dissolved solids (TDS), chloride (Cl − ), sodium (Na + ), potassium (K + ), sodium absorption ratio (SAR), and bicarbonate.The substances measured in this study include bicarbonate (HCO3), calcium (Ca 2+ ), magnesium (Mg 2+ ), sulfate (SO4), and acidity (pH + ).To assess the surface runoff potential for drinking water and agriculture, various charts such as Schuler (for drinking water) and Wilcox (for agricultural water) are utilized in AqQA v1.1.0software.Moreover, to check the trend of ions in the past years, the Mann-Kendall test [53,54] and ordinary-least-square (OLS) regression slope [55] are used.

Trend Analysis and Predicting Land Use/Cover Change
Available sensors and satellite images were utilized for this study, ensuring compatibility with other data.Land use data are prepared using Landsat satellite images based on the TM, ETM+, and OLI sensors for three span times: 1991, 2010, and 2020.First, pre-processing was conducted on satellite images, and field surveys were conducted in June and July 2019 to obtain training samples.The land use classes in the study area are monitored using the maximum likelihood algorithm, which is a conventional method known for its high accuracy [24,29,56].
To investigate the impact of drought on vegetation in the studied watershed, the Normalized Difference Vegetation Index (NDVI) time series [57] from Landsat satellite images was analyzed [58] over a 20-year period (2001-2020).The LCM method was utilized in TerrSet 19.0.8 software to forecast future land use conditions.To achieve this, this study initially examined land use changes over the study period.Subsequently, various scenarios were evaluated to determine the most accurate implementation submodel and independent variables.These variables, including the Digital Elevation Model (DEM), slope, distance from residential areas, waterways, roads, and agricultural lands, are found to have a strong correlation with land use changes in the study area.The transformation of each land use and the prediction of land use changes are also modeled using a multi-layer perceptron function and Markov chain.Kappa coefficients and overall accuracy were used to assess the accuracy of the land use maps that are created.TerrSet 19.0.8 software and ENVI 5.5.2 software [56] were utilized for all the mentioned steps.The findings indicate that the level of drought has varied between −3.3 (indicating very severe drought) and 2.4 (indicating extremely humid conditions) in recent years.Moreover, in Figure 4 and Table 3, the forecast of the future drought condition for a period of 27 years (2023-2050) under the influence of two optimistic (RCP2.6)and pessimistic (RCP8.5)scenarios is shown.(5).Severe drought (6).Extreme drought (7).

Results of Climatic Drought Assessment
The SPI for Rig Cheshmeh station in the pessimistic scenario range from 2.76 to −2.11.Its lowest value is expected to occur in 2029-2030 during a very severe drought, while the highest value is expected in 2039-2040 (extremely humid).In addition, the index values for the Rig Cheshmeh station in the optimistic scenario vary from 2.77 to −2.07.Its lowest value will occur in 2029-2030 (very severe drought) and the highest value will occur in 2039-2040 (extremely humid).Therefore, it can be expected that, in the water year 2029-2030, there is a possibility of a very severe drought in this region.As a result, it is crucial to take the necessary measures to address this potential climatic challenge.
The SPI values for Soleyman tangeh station in the pessimistic scenario range from 1.89 to −2.16.The lowest value (−2.16) was obtained in 2049-2050 during a very severe drought, while the highest value (1.89) was obtained in 2041-2042 during very wet conditions.In addition, the values of this index for Soleyman tangeh station in the optimistic scenario vary from 1.94 to −2.18.The lowest value occurs in 2049-2050, indicating a very severe drought, while the highest value is observed in 2041-2042, indicating relatively wet conditions.Therefore, it can be expected that in the water year 2049-2050, there is a possibility of a severe climatic drought in this region.In addition, the results show no significant difference in the future drought index between the two scenarios.
The SPI values for Kordkhil station in the pessimistic scenario vary from 1.58 to −2.85.The lowest value (−2.85) is expected in 2032-2033, indicating a very severe drought.On the other hand, the highest value (1.58) is projected in 2029-2030, indicating a very humid condition.In the optimistic scenario, the Kordkhil station will experience its lowest value of this index in 2032-2033 (−2.92, indicating a very severe drought) and its highest value in 2029-2030 (1.59, indicating a very wet period).Therefore, the results of Kordkhil station are somewhat different compared to the other two stations that are located at a higher altitude, and this difference is probably due to the proximity of this station to the Caspian Sea and the climate of that region being affected by the humidity of the coast.It can also be expected that there will be a very severe climatic drought in the water year 2032-2033.Therefore, the results of this research regarding the increase in drought periods compared to the base period are in accordance with the findings of [43].

Results of Hydrological Drought Assessment
Figure 5 illustrates the drought trend based on observation data from two hydrometric stations using the SDI.The data span a 30-year period (1990-2021) and is analyzed in 12-month increments.In this study, two hydrometric stations were used to accurately measure the amount of discharge in both the upstream and downstream areas of the watershed.Therefore, the Parvich abad station (upstream), the Kordkhil station (downstream), and the outlet of the watershed are selected.Table 4 shows the results of the rainfall-runoff model simulation.Herein, both values are measured in mm y −1 .The coefficient values obtained in the two stages of evaluation and recalibration were 0.48 and 0.53, respectively, which are consistent with the findings of Tramblay et al. [59].This model is very sensitive to the time matching of rainfall and discharge, and due to the lack of comprehensive data in the study area, it was not possible to implement it for different stations with similar time periods.In addition, the data obtained from the prediction of the Lares-WG climate model as well as the observational data of discharge have uncertainty, which has caused a slight deviation in the model results.* Q represents the portion of precipitation that contributed to the formation of the watershed river throughout the measurement and evaluation period, while P represents the total precipitation during that period.
Figure 7 shows the results of analyzing the hydrological drought trend at the Kordkhil station.The analysis was conducted using the SDI over a 20-year period, with time intervals of 3, 6, and 12 months.The analysis considered both the RCP8.5 and RCP2.6 scenarios.The figures clearly indicate that the most significant drought changes in the future period, in both scenarios, will occur every three months.The magnitude of these changes will range from −1.5 (moderate drought) to 3.5 (non-drought).The results indicate that climate change has not had a significant impact on discharge.The likely cause of this is heavy rainfall events resulting from climate change, which leads to an increase in flash floods in the river.The research findings only show the effect of climate change on discharge, and it is possible that other factors like changes in land use may also have an influence.The results of this research are somewhat consistent with the findings of previous research [4] regarding the limited impact of climate change on runoff.Xue et al. [4] discovered that human activities, such as alterations in land use, have had a greater influence on runoff than climate change.Zhang et al. [33]   The hydrological drought situation and index values for the predicted period (2021-2040) of Kordkhil station are presented in Table 5.The analysis was conducted using a time step of 12 months and two scenarios, RCP2.6 and RCP8.5, with the SDI.The values of this index for Kordkhil station in the pessimistic scenario range from 3.15 to −1.49.The lowest value, indicating a moderate drought, is expected to occur in 2034-2035 (−1.49), while the highest value, indicating no drought, is projected for 2028-2029 (3.15).The index values for Kordkhil station in the optimistic scenario range from 3.18 to −1.62.The lowest value, indicating severe drought, is projected to occur in 2034-2035 (−1.62), while the highest value, indicating no drought, is expected in 2028-2029 (3.18).Since this station is situated at the outlet of the watershed and receives most of the runoff from the upstream watershed, the river's discharge is typically high and less influenced by changes.Additionally, the region's proximity to the Caspian Sea has led to increased rainfall, resulting in higher runoff.These findings are consistent with the research conducted by Kazemzadeh and Malekian [19].Based on the available results, there is a high likelihood of severe hydrological drought in both scenarios in 2034-2035.This drought is likely caused by changes in precipitation during this period.Therefore, it is crucial to take necessary measures in the year before, during, and after to minimize the impact on water farmers along the Tajan river.Prior notifications should be given to encourage appropriate farming practices during those years.The investigation of the hydrological drought index (SDI) compared to the meteorological drought index (SPI) in the future indicates their compatibility in most years.There is a notable increase in precipitation from 2028 to 2030 (Figure 8).There is also the possibility of a wet year between 2039 and 2041.The most severe drought is projected to occur between 2034 and 2036.Therefore, in the future, we can expect alternating periods of drought and wet years.It is important to take measures to adapt and mitigate the associated risks.The climate fluctuations result in heavy rainfall, floods, and increased stream flow during certain times of the year.The amount of flow in an area is not solely determined by precipitation; factors such as land use and cover also significantly influence runoff [15,48].

Results of Surface Water Quality Assessment
The results of changes in surface water quality in the Tajan watershed at two stations, Parvich abad (upstream) and Kordkhil (downstream), is reported in Table 6.The physical characteristics of both stations indicated the presence of calcium carbonate in the water, suggesting the existence of limestone formations in the upstream areas.The water has a high hardness due to the presence of calcium carbonate, with a concentration of 273 mg l −1 at the outlet.
The water salinity is higher at the outlet station of the watershed.The electrical conductivity at the Kordkhil station is about 770 µmohs cm −1 , indicating a high risk of salinity.The changes in the sodium absorption ratio (SAR) at both stations indicate an increase in water salinization at the outlet station (Kordkhil).This raises concerns about potential soil degradation if the downstream lands are irrigated with this water.Land use/cover changes have the great impact on surface water quality.The expansion of agricultural lands resulting from the use of chemical fertilizers, along with soil erosion caused by agriculture can have an impact on water quality (Moradzadeh et al., 2024).The results also indicate an increase in the dissolution of solids at the downstream station.The changes indicate a decline in water quality at the downstream station, likely due to the convergence of multiple rivers at that location.The conventional Wilcox diagram was used to assess the water quality for agricultural purposes.The water quality at Kordkhil station (left) is divided into two ranges, C2S1 and C3S1 (Figure 9).Given the low sodium content indicated by S1, this water is suitable for all types of soil and poses a low risk of soil degradation.C2 also indicates water with moderate salinity, which can be suitable for agricultural irrigation with proper leaching.C3 indicates that the water has high salinity, making it unsuitable for irrigating agricultural lands, particularly those with inadequate drainage.This poses a risk of land salinization.The water quality results from the Parvich abad (right) station indicate that it is slightly better than the Kordkhil station.Therefore, the water quality at this station falls within the C2S1 class range.Therefore, this water can be safely used for agricultural irrigation.The higher quality of this station is attributed to its higher altitude and lower solute dissolution in the water.
The drinking water quality assessment of the Tajan watershed in two stations, Kordkhil and Parvich abad, is depicted in Figure 10 and Table 7.The results illustrate significant fluctuations in water quality at the Kordkhil station over the years, with the most notable decline occurring in recent years (2017)(2018)(2019)(2020)(2021).This indicates a high concentration of ions in the water of the area under study.The concentration range of ions varies between 0.007 and 0.25 meq kg −1 .The concentration range of ions in Parvich abad station (0.003-0.15 meq kg −1 ) is lower than that of Kordkhil station.Additionally, there is no specific trend in the water quality changes at this station.A decline in water quality can be observed in 2017-2018.As evident from Tables 8 and 9, the trend of water ions at Kordkhil station over the past 20 years (2001-2021) has shown a significance level of 0.05 in all parameters, except for magnesium (Mg).This decrease in the Mg parameter suggests a decline in water quality during this period.
Because of its location at the outlet of the watershed, Kordkhil station receives water from various rivers in the upstream areas of the watershed.These rivers have different land use and formations, which results in lower water quality compared to the higher stations.Most of the land use changes and human activities are concentrated in the downstream areas of the watershed, directly impacting the water quality of this station.To compare the water quality downstream and upstream of the watershed, we also analyzed the trend of water quality at the Parvich abad station, which is located upstream of the Shahid Rajaei Dam.The results of this analysis showed that, unlike the Kordkhil station, the trend of ions in this station was decreasing in all parameters except for calcium (Ca) and pH.These results can also be attributed to the reduced influence of human activities on the water entering this station.The trend diagram of water quality parameters in Kordkhil station, as an example, is also shown in Figure 11.In addition, Figure 12 shows the trend of some water quality parameters with two types of climatic and hydrological drought.The results demonstrate the different behavior of water quality parameters and drought indices in the study area.

Results of Land Use/Cover Changes Assessment
The study of land use changes revealed that the dominant land uses are forest, agriculture, rangeland, garden, and residential (Figure 13).The results of land use changes, based on the continuity scenario, indicate that if future changes follow the trends of the past 29 years, the area of forest lands will decrease by 19.77 km 2 , which is equivalent to a 2.02% decrease in the current forest land area over the next 20 years.In this scenario, the rangeland area may decrease by −5.74%.Residential, agricultural, and garden lands are projected to increase by 0.45, 6.82, and 9.17, respectively, during this period (Figure 12, Table 10).
If the current trend of land use change in the studied area continues, approximately 213.21 km 2 of rangelands and 64.84 km 2 of forests will be converted into agricultural lands by 2040 (Table 10).The least amount of change occurs when garden lands are converted into residential lands.The expansion of agricultural and residential areas, often located along riverbanks, has contributed to the increased use of chemical fertilizers and the discharge of domestic sewage into rivers, leading to environmental pollution.An increase in sediment in surface water and a decrease in water and soil quality are additional consequences of expanding agricultural and residential land, which is also demonstrated by previous research [20][21][22]34].Therefore, land use management in a region can significantly impact the health of the watershed.The results of investigating the impact of drought on NDVI in the studied watershed over a 20-year period (2000-2020) is shown in Figures 14 and 15.The results indicate a recent increase in vegetation cover, likely due to the expansion of agricultural lands or gardens in the study area.The level of agricultural and garden lands has been increasing in recent years.Despite a decrease in precipitation in recent years, the amount of greenness has not significantly decreased, except in 2021-2022 (Figures 14 and 15).

NDVI
The changes in the NDVI compared to the meteorological index of SPI from 2000 to 2021 in various stations clearly shows that the NDVI does not align with the SPI, particularly in recent years (Figure 16 and Table 11).It can be concluded that the drought had minimal impact on the greenness index of the Tajan watershed (Figure 16 and Table 11).This is likely because the watershed has ample surface water resources.These results contradict the findings of Zhang et al. [33] regarding the adverse impact of drought on vegetation.Considering that the NDVI only measures vegetation, it is possible that the expansion of agricultural and garden lands during this period has led to an increase in the NDVI.Therefore, this index should be considered in conjunction with other factors when assessing the impact of drought and the overall health of the watershed.In this regard, Ji and Peters [57] reported a significant effect of seasonality on the NDVI-SPI relationship.The high correlations were obtained only in the middle of the growing season, which reflects the sensitivity to water availability in this time.

Conclusions
The analysis of both climate scenarios (RCP2.6 and RCP8.5) indicates the potential for severe droughts in certain years in the future, despite the possibility of wet climate conditions in other years.In addition, the hydrological drought in the Tajan river has exhibited fluctuations, with some years experiencing a decrease in flow and others showing an increase.The results of land use analysis indicate an increase in agricultural, garden, and residential land, while forests and rangelands have decreased.Therefore, the reduction in permanent forest and rangeland cover may lead to increased runoff and flooding, resulting in higher flow rates during certain times of the year.Therefore, a drought in a region experiencing simultaneous changes in climate and land use can yield varying outcomes.The changes in NDVI did not comply with the SPI.Changes in a watershed have a significant impact on surface water quality, making it an important health factor to consider.Given the observed increase in agricultural land, particularly in rainfed areas, in the studied region, it is anticipated that without proper measures, water quality will decline in the future.The changes mentioned above pose a threat to the health of watersheds from the standpoint of the number and intensity of floods, erosion and sedimentation, decreased water penetration into aquifers, increased pollutants, and decreased water resources.

Figure 2 .
Figure 2. Diagram of the implementation of the hydrological model of IHACRES.

Figure 3
Figure 3 displays the trend of climatic drought based on observational data from various stations between 1993 and 2020.Changes in meteorological drought conditions over a 27-year period indicate climate fluctuations, which are signs of climate change in the region.Based on the obtained results, it cannot be definitively stated that the next 30 years will be consistently dry or wet, although it is necessary to use adaptation approaches to mitigate the long-term effects of drought.

Figure 3 .
Figure 3. Climatic drought trends in the basic observational data of different stations.

Figure 4 .
Figure 4.The change trend of SPI in different time steps and different scenarios of (A) Soleyman tangeh, (B) Rig Cheshmeh, and (C) Kordkhil stations.

Figure 5 .
Figure 5. Hydrological drought trend of observational data of selected stations.

Figure 6
Figure6shows the comparison graph of the simulated and observed discharge of Kordkhil station.The IHACRES model accurately simulates the discharge, and the two graphs show a high level of agreement, as depicted in Figure6.
Figure7shows the results of analyzing the hydrological drought trend at the Kordkhil station.The analysis was conducted using the SDI over a 20-year period, with time intervals of 3, 6, and 12 months.The analysis considered both the RCP8.5 and RCP2.6 scenarios.The figures clearly indicate that the most significant drought changes in the future period, in both scenarios, will occur every three months.The magnitude of these changes will range from −1.5 (moderate drought) to 3.5 (non-drought).The results indicate that climate change has not had a significant impact on discharge.The likely cause of this is heavy rainfall events resulting from climate change, which leads to an increase in flash floods in the river.The research findings only show the effect of climate change on discharge, and it is possible that other factors like changes in land use may also have an influence.The results of this research are somewhat consistent with the findings of previous research[4] regarding the limited impact of climate change on runoff.Xue et al.[4] discovered that human activities, such as alterations in land use, have had a greater influence on runoff than climate change.Zhang et al.[33] also demonstrated the significant influence of climate change on runoff in comparison to land use.The results of the present study also indicate that the discharge of the Tajan watershed is unlikely to change significantly in the future due to climate change.However, there is a possibility that other factors such as land use change or water transfer may decrease the discharge of the Tajan river.

Figure 8 .
Figure 8.The trend of future changes in SPI relative to SDI in a 12-month time step based on two scenarios RCP2.6 and RCP8.5.

Figure 9 .
Figure 9. Agricultural water quality in Kordkhil (top) and Parvich abad (down) stations presented in Wilcox diagrams.

Figure 11 .
Figure 11.The trend of water quality parameters of Kordkhil station in the period 2001-2021 (Red and green lines respectively indicating the Sen's slope and Mann-Kendall trend).

Figure 12 .
Figure 12.Simultaneous trends of some water quality parameters and drought.

Figure 13 .
Figure 13.The trend of land use changes in the Tajan watershed in the basic and future periods.

Figure 15 .
Figure 15.NDVI map of Tajan watershed in the period of 2000-2020.

Figure 16 .
Figure 16.NDVI and SPI changes trend with a 12-month time step in different stations.

Table 3 .
Meteorological drought situation of Tajan watershed stations in the future using two climate scenarios.

Table 4 .
The results of calibration and evaluation of the model at the Kordkhil station of the Tajan watershed in the period of 1997-2017.

Table 5 .
Hydrological drought condition of Kordkhil station in the future under two different scenarios.

Table 6 .
Physical characteristics of water quality in two studied stations.

Table 8 .
Trend analysis of water quality parameters of Kordkhil station in the period 2001-2021 at the confidence level of 95%.
These arrows indicate the upward and downward trend.

Table 9 .
Trend analysis of water quality parameters of Parvich abad station in the period 2001-2021 at the confidence level of 95%.
* These arrows indicate the upward and downward trend.

Table 10 .
The difference in land use area in the two periods of 2020 and 2040.

Table 11 .
Statistical analysis of changes in SPI and NDVI.