The Role of Moisture Sources and Climatic Teleconnections in Northeastern and South-Central Iran’s Hydro-Climatology

: Iran faces climate disparities due to extreme topographic anomalies, the Caspian Sea and the Persian Gulf water bodies, inﬂuences from diverse air masses and moisture sources, and its considerable area. FLEXPART model has been utilized to determine the main marine and continental moisture sources for south-central (Shiraz box) and northeastern (Mashhad box) parts of Iran. The marine moisture sources directly inﬂuenced extreme drought and wet conditions in Shiraz and Mashhad boxes during the wet period, while no correlation was observed during the dry period. In addition to local components, extreme drought and wet conditions have also been inﬂuenced by the climatic teleconnections. Extreme drought conditions mainly occurred during the La Niña phase, while wet conditions mainly occurred during the El Niño phase. Scrutinizing the effect of marine moisture sources on the hydrology of water resources demonstrated that the moisture contribution from the Arabian Sea directly inﬂuenced the discharges of Chenar-rahdar (in the Shiraz box) and Kardeh (in the Mashhad box) rivers during the wet period. However, the Red Sea inversely correlated with the discharges of both rivers during the dry period. Hydrogeologists, hydrologists, and meteorologists can utilize the outputs of this survey to develop climatology and hydrology models in the future.


Introduction
In recent decades, numerous studies have been done to deepen our understanding about the identification of moisture sources and mechanisms controlling moisture transport in the atmosphere [1] as a major component of the water cycle. Although studies on the identification of precipitation moisture sources and their transport mechanisms are growing, few studies focus on various microscale and mesoscale parameters which influence moisture sources and their variations throughout a year. Sea surface temperature (SST), wind speed, solar radiation, as well as climatic teleconnections are among the parameters which influence moisture sources. According to the place of their origin, the major climatic teleconnections are: IOD (Indian Ocean Dipole) which occurs in the Indian Ocean; MJO (Madden-Julian Oscillation) which occurs between the Indian and Pacific Oceans; PNA (Pacific North American), ENSO (El Niño-Southern Oscillation), Southern Oscillation Index (SOI) and PDO (Pacific Decadal Oscillation) which occur in the Pacific Ocean; AMO (Atlantic Multi-decadal Oscillation), NAO (North Atlantic Oscillation), and TAV (Tropical Atlantic Variability) which occur in the Atlantic Ocean; AO (Arctic Oscillation) and AAO (Antarctica Oscillation) which occur in the polar regions and QBO (Quasi-Biennial Oscillation) which occurs in the tropical stratosphere. The influence of these climatic indices on atmospheric and hydrological parameters have been appraised in many studies in Iran [2][3][4][5][6][7][8] and elsewhere in the world [9][10][11][12][13][14][15][16][17][18].
Moisture sources are normally identified using analytical models, numerical models, and isotope techniques [19]. Isotope techniques have been widely used in hydro-climatological studies since Craig discovered the relationship between 18 O and 2 H isotopes in the water molecules of natural meteoric waters from many parts of the world [20]. On the other hand, numerical models have less often been applied compared to physical methods and are mainly limited to the studies in recent decades. Among the various recently developed numerical models, the Lagrangian dispersion model FLEXPART developed by [21,22] has a great importance. This model enables scientists to track air parcels backward and forward in time. Thus, the outputs of this model can be used to compute the gain and loss of moisture along air particle trajectories [1]. This model has been widely and successfully utilized in many projects in different parts of the world [21][22][23][24][25][26][27][28][29][30].
The moisture uptake rate from various moisture sources directly influences the amount of precipitation and also extreme drought and wet conditions in each region. In addition to climatological or meteorological drought (as the quickest response of precipitation decrease in the environment), hydrological drought has also had a great role in hydrological and hydroclimatological studies. In contrast to climatological or meteorological drought, hydrological drought needs more time to reveal itself and causes a decline in surface and groundwater levels, deteriorates water quality, and leads to numerous other economic and social obstacles. Standardized streamflow index (SSI) [31] is normally utilized to quantify extreme drought and wet conditions in surface water resources.
Iran has always been threatened by extreme drought and wet conditions due to its semi-arid and arid climate. Freshwater availability in Iran has been estimated to be around 2000 m 3 per capita yr −1 in 2000 and has been predicted to decrease to 1500 m 3 per capita per yr −1 until 2030 [32,33]. During the last few decades, numerous villages-mainly in central, eastern, and southern parts of Iran-have been evacuated due to the lack of water and millions head of livestock and hundreds of thousands of hectares of farmlands have been damaged due to far-reaching and lengthy droughts. The intense 1998-2000 drought in Iran inflicted a 3.5-billion-dollar damage and many major lakes and rivers were dried during this period across the country. Although numerous publications exist in the literature regarding spatial and temporal variations of precipitation and extreme drought and wet conditions in Iran including [3,[34][35][36][37][38][39], few studies such as [3,40] have focused on the main moisture sources of Iran's precipitations and their variations across the country.
As most of the previous investigations in Iran have studied climatological variations across the country, the objective of the present study is to identify the marine and continental moisture sources and their contribution rates in the northeastern (Mashhad box) and south-central (Shiraz box) parts of Iran. The effect of manifold marine sources on extreme drought and wet conditions and also two river discharges (Chenar-rahdar River in Shiraz box and Kardeh River in Mashhad box) has been investigated. In addition, the role of climatic teleconnections on marine moisture sources, extreme drought, and wet conditions, and also Chenar-rahdar and Kardeh river discharges have been investigated using regression methods.

Study Area
With an area of 1,648,195 km 2 in the Middle East, Iran is situated in latitudes 25 • N to 40 • N and longitudes 44 • E to 64 • E, with a semi-arid and arid climate according to Koppen's method [41]. Iran's climate has a broad spectrum ranging from the arid zones of central Iran with less than 120 mm annual precipitation to rainy jungles in the Caspian Sea coasts with more than 1000 mm average annual precipitation. This intense climate variation is due to two large mountain ranges (Alborz and Zagros) which border two extremely arid deserts in central Iran "Dasht-e Lut and Dasht-e Kavir". The Alborz mountain range prevents the flow of moisture from the Caspian Sea to the central part of Iran, while the Zagros mountain range prevents the flow of moisture from the west and southwest of Iran toward central Iran [34]. In addition, the Persian Gulf and the Oman Sea in the south, and the Caspian Sea in the north of Iran predominantly affect Iran climate and transfer considerable amount of moisture via various air masses to Iran.
Five major air masses influence Iran through the year. The Mediterranean air mass (MedT) which influences the western and southwestern parts of the country. Mediterranean cyclones contain high amounts of moisture and enhance cloudiness and precipitation mainly in the western part of Iran, while the Mediterranean frontal system and Sudan heat-flow system influence the southwestern part of Iran. The simultaneous merging of the Mediterranean frontal system and Sudan air mass cause huge precipitation events. The intensity of the Mediterranean system is weakened by the penetration of the Azores subtropical high-pressure system over Iran during April and May. The continental tropical air mass (cT) known as Sudan air mass mainly affects the southwest of Iran. The northwestern part of Iran is affected by maritime polar (mP) air mass. This air mass does not often influence Iran and brings humid cold air from northern and north parts of Europe. Continental polar (cP) air mass (known as the Siberian air mass) affects the northern and northeastern parts of Iran and its interaction with the Alborz mountain range causes huge rainfalls in Caspian Sea coast. Finally, the maritime tropical (mT) air mass rarely causes sudden intense rains in the south and southeast part of Iran in the dry period [34,42] (Figure 1). Caspian Sea in the north of Iran predominantly affect Iran climate and transfer considerable amount of moisture via various air masses to Iran. Five major air masses influence Iran through the year. The Mediterranean air mass (MedT) which influences the western and southwestern parts of the country. Mediterranean cyclones contain high amounts of moisture and enhance cloudiness and precipitation mainly in the western part of Iran, while the Mediterranean frontal system and Sudan heat-flow system influence the southwestern part of Iran. The simultaneous merging of the Mediterranean frontal system and Sudan air mass cause huge precipitation events. The intensity of the Mediterranean system is weakened by the penetration of the Azores subtropical high-pressure system over Iran during April and May. The continental tropical air mass (cT) known as Sudan air mass mainly affects the southwest of Iran. The northwestern part of Iran is affected by maritime polar (mP) air mass. This air mass does not often influence Iran and brings humid cold air from northern and north parts of Europe. Continental polar (cP) air mass (known as the Siberian air mass) affects the northern and northeastern parts of Iran and its interaction with the Alborz mountain range causes huge rainfalls in Caspian Sea coast. Finally, the maritime tropical (mT) air mass rarely causes sudden intense rains in the south and southeast part of Iran in the dry period [34,42] (Figure 1).  Extreme topographic anomalies, the Caspian Sea and the Persian Gulf water bodies, diverse air mass and moisture source influences and considerable area cause climate disparities across Iran. Two boxes, sized 1 • × 1 • , including Shiraz (A) and Mashhad (B) have been chosen in the south-central and northeastern parts of Iran, respectively ( Figure 1). The Shiraz box in the south-central part and the Mashhad box in the northeastern part of Iran have been chosen as indicators to investigate the spatial climate variations across Iran. The Shiraz box is located in a synclinal plain of the Zagros Mountain range in Fars province in the south-central part of Iran. The climate of Shiraz box is semi-arid, with the average annual precipitation of 337 mm and long-term mean annual temperature of 18 • C. Chenar-rahdar River with a wide discharge range from 0.02 to 5.28 m 3 /s −1 and the discharge average of 0.75 m 3 /s −1 (measured in Chenar-rahdar gauge station) is the main river in the Shiraz box ( Figure 1). The Mashhad box located in the northeast of Iran in Kopeh Dagh mega anticline is a semi-arid region with an average annual precipitation of 250 mm [43]. Kardeh River is one of the largest rivers in Khorasan Razavi province located in the Mashhad box with a discharge varying from 0.05 to 6.16 m 3 /s −1 and average discharge of 0.64 m 3 /s −1 (measured in Kardeh gauge station) ( Figure 1).
The monthly precipitation variations of Mashhad and Shiraz boxes ( Figure 2) confirms the climate variations across Iran. In Shiraz box, most of the precipitation occurs during winter to early spring (December to March). However, in Mashhad box, spring time precipitation (March to May) is dominant, while winter precipitation events (December to February) have a minor role. During summer (June to September), the Mashhad box demonstrates a very low amount of precipitation, while for the Shiraz box, precipitation is almost null.

Materials and Methods
To investigate the origin of moisture arriving in Shiraz and Mashhad boxes has been utilised the Lagrangian particle dispersion model (FLEXPART v9.0) developed by Stohl and James [21,22]. Respect previous versions of the model, in our experiment, FLEXPART v.9 allows considering a greater number of air particles, near 2 million, each representing a part of the whole atmospheric mass. There are diverse options for how particles are generated in the model, and what they represent (air, or a certain trace gas). FLEXPART permits the tracking of global air masses forward and/or backward in time, or it just considers the air masses residing over a target region using the atmospheric winds. Along each particle trajectory, the specific humidity (q) can be interpolated to compute the net change of the water content because of the mutual effects of evaporation (e) and precipitation (p). This calculus in Equation (1), was obtained over cycles of 10 days, considering the water vapor average residence time in the atmosphere [44,45]. In this calculation, a small error occurs for considering constant the mass (m) of every particle. Nevertheless, Stohl and James, [22] consider that the atmosphere's mass changes slightly through the addition and removal of water.
Resolving the Equation (1) for all particles (K), residing over an area A, permits to diagnose the

Materials and Methods
To investigate the origin of moisture arriving in Shiraz and Mashhad boxes has been utilised the Lagrangian particle dispersion model (FLEXPART v9.0) developed by Stohl and James [21,22]. Respect previous versions of the model, in our experiment, FLEXPART v.9 allows considering a greater number of air particles, near 2 million, each representing a part of the whole atmospheric mass. There are diverse options for how particles are generated in the model, and what they represent (air, or a certain trace gas). FLEXPART permits the tracking of global air masses forward and/or backward in time, or it just considers the air masses residing over a target region using the atmospheric winds. Along each particle trajectory, the specific humidity (q) can be interpolated to compute the net change of the water content because of the mutual effects of evaporation (e) and precipitation (p). This calculus in Equation (1), was obtained over cycles of 10 days, considering the water vapor average residence time in the atmosphere [44,45]. In this calculation, a small error occurs for considering constant the mass (m) of every particle. Nevertheless, Stohl and James, [22] consider that the atmosphere's mass changes slightly through the addition and removal of water.
Resolving the Equation (1) for all particles (K), residing over an area A, permits to diagnose the surface freshwater flux (Evaporation (E) − Precipitation (P)) according to the Equation (2) To identify the moisture sources for the Mashhad and Shiraz boxes (Figure 1), the air masses residing over each box were tracked backward in time. Where (E − P) > 0 indicate those regions where prevail the moisture uptake by air masses before reaching the Mashahd and Shiraz boxes. Conversely, those regions where the moisture loss prevails, (E − P) < 0, are considered moisture sinks. Positive values of (E − P) have only been considered, due to the main goals of this survey which tries to determine the predominant moisture sources for the Shiraz and Mashhad boxes. The analysis was performed for the wet (November-April) and dry (May-October) periods from 1981 to 2015. The length of dry and wet periods has been specified in previous studies of the precipitation in Iran [36,37] and also Iran meteorological database. To delimit the most important seasonal moisture sources according to the moisture uptake intensity, a percentile method was applied. This approach has been previously utilised for the same objective by several authors [28,46,47]. The sources are separated into marine and continental, which further allows comparing the amount of moisture uptake ((E − P) > 0) over each one.
The main advantage of FLEXPART respect Eulerian methods is that it allows to compute the budget of (e − p) along the particles trajectories [19]. Besides, it utilises turbulence and convection parametrizations, which are particularly important to simulate the moisture changes in the atmosphere [48,49]. Nevertheless, although the model has numerous advantages, it has some limitations-including its weakness to calculate the precipitation (P) and evaporation (E) severally. However, when rain falls, it normally exceeds evaporation. Thus, supposing that E and P cannot coexist in the same location at the same time, the average instantaneous rates of E and P over a longer period (e.g., years), permits to estimate average E and P [21]. Another problem discussed by these authors is that q fluctuations along individual trajectories may occur for nonphysical reasons, such as q interpolation, or trajectory errors affecting the values of (e − p). However, these errors are partially annulled when considering the many particles in an atmospheric column over the area A. Further details about the model are described by [21].
In this study, FLEXPART input's data belong to ERA-Interim reanalysis [50] at 6 h intervals (00:00, 06:00, 12:00, and 18:00 UTC) at a resolution of 1 • in latitude and longitude on 60 model levels from 0.1 to 1000 hPa has been used. The closure of the water cycle in the global atmospheric reanalysis models-ERA-Interim Re-Analysis; Modern-Era Retrospective Analysis for Research and Applications (MERRA); and Climate Forecast System Reanalysis (CFSR)-was evaluated and compared with independent observations in the period 1989-2006, by [51]. Their results revealed that generally, ERA-Interim shows a reasonable closure of the terrestrial and atmospheric water balance and a reasonable agreement with the observation datasets. When using FLEXPART with ERA interim datasets, it is expected to reduce the uncertainties of (E − P) and thus representing better the hydrological cycle.
Precipitation data for both the Mashhad and Shiraz boxes have been provided by Iran Meteorological Database from 1981 to 2014. Drought and wet conditions in the Mashhad and Shiraz boxes have been identified using the standardized precipitation index (SPI) [52] by DIP software developed in Iran's Ministry of Energy [53]. The positive SPI values demonstrate wet conditions, while negative values demonstrate drought conditions. The intensities of wet and drought conditions  (Table 1). On the other hand, the hydrological droughts in Chenar-rahdar River in the Shiraz box and Kardeh River in the Mashhad box have been investigated by the standardized streamflow index (SSI) [31]. Monthly discharges data of Chenar-rahdar (Sept. 1982 to Sept. 2013) and Kardeh (Sept. 1984 to Sept. 2008) rivers were provided by Khorasan Razavi and Fars Regional Water Authority. Other precipitation datasets utilized in this study with a resolution of 0.5 • × 0.5 • were from the Climatic Research Unit CRU3.23TS [54]. Datasets of Omega (vertical velocity) were used to diagnose the seasonal atmospheric stability over Shiraz and Mashhad boxes was available at [55]. The regression method was used to verify the effects of climatic indices, namely Bivariate ENSO (BEST) [56], North Atlantic Oscillation (NAO) [57], Southern Oscillation Index (SOI), and Quasi-Biennial Oscillation (QBO) indices (derived from [58,59]) on moisture uptake sources, extreme drought and wet conditions in the Shiraz and Mashhad boxes, and Chenar-rahdar and Kardeh river discharges.

Identification the Main Moisture Uptake Sources for the Shiraz and Mashhad Boxes and Their Contribution Rates
The dominant moisture origins have been identified for the Shiraz and Mashhad boxes (Figure 3). The p95 delimitation of maximum moisture uptake is represented by the light blue line. Thus, over these regions within the light blue line is evident that the number of trajectories that coincide is high. In the wet period, the maximum (E − P) > 0 values obtained in air masses in travel to Shiraz and Mashhad boxes occur over the boxes themselves and to the south. In the dry period, the most intense (E − P) > 0 still occur over both boxes, but it is clear according to the visual analysis of Figure 3b,d, that most important regions supplying moisture to the boxes are located to their northeast.
Those regions of great moisture uptake presented in Figure 3 have been divided considering the geographical locations. The spatial distribution of moisture uptake sources over the Shiraz and Mashhad boxes has been shown in Figure 4 and their contribution rates have been tabulated in Table 2.
During the wet period, marine moisture sources provide 52 and 40 percent of total moisture income to Shiraz and Mashhad boxes, respectively. However, during the dry period, the role of marine moisture sources highly decreases and reaches 20 and 14 percent of total moisture uptake by air masses in travel to Shiraz and Mashhad boxes, respectively [60].  During the wet period, marine moisture sources provide 52 and 40 percent of total moisture income to Shiraz and Mashhad boxes, respectively. However, during the dry period, the role of marine moisture sources highly decreases and reaches 20 and 14 percent of total moisture uptake by air masses in travel to Shiraz and Mashhad boxes, respectively [60].
Although the continental moisture sources play a dominant role over Iran mainly during the dry period, these sources do not lead to a considerable precipitation and their influences are limited to local precipitation events [34,40]. (E − P) > 0 and precipitation amount correlation in Shiraz and Mashhad boxes have also asserted this reality (Supplementary Figure S1). Because the following survey tries to determine the main moisture sources which affect Iran's precipitation throughout the Although the continental moisture sources play a dominant role over Iran mainly during the dry period, these sources do not lead to a considerable precipitation and their influences are limited to local precipitation events [34,40]. (E − P) > 0 and precipitation amount correlation in Shiraz and Mashhad boxes have also asserted this reality (Supplementary Figure S1). Because the following survey tries to determine the main moisture sources which affect Iran's precipitation throughout the year, only marine sources are considered in this study. The percentage of moisture uptake from various marine sources for Shiraz and Mashhad boxes has been tabulated in Table 2. The superior moisture uptake sources for the south-central (the Shiraz box) and northeastern (the Mashhad box) parts of the country and also for the dry and wet conditions are different. In the south of Iran, the role of moisture arriving from the Arabian Sea is dominant during the wet period, while the moisture sources located in the Mediterranean and the Arabian Seas are dominant in the north of Iran. During the dry period, the moisture uptake from the Caspian Sea and the Persian Gulf dominantly intensify in the north and south of Iran, respectively [60].  During the wet period, the roles of the Mediterranean (MedT) and Continental polar (cP) air masses in transferring moisture to the northeast of Iran (Mashhad box) are more dominant, while the continental tropical (cT) air mass transfers the dominant part of moisture to the south-central part of Iran (Shiraz box). Although moisture uptake from marine and continental sources is also high in the dry period, this moisture cannot lead to precipitation due to the influence of Azores subtropical high-pressure system over Iran causing an intense atmospheric stability over large parts of Iran and preventing precipitation [34,40,61]. Monthly precipitation distribution and ω variations over Mashhad and Shiraz boxes have been demonstrated for the dry and wet periods in Figure 5. The positive values of omega indicate descending motion of air parcels and stability in atmosphere and on the contrary negative values of Omega represent ascending motion and instability in atmosphere [62].

Studying the Role of Marine Moisture Uptake Sources on Extreme Drought and Wet Conditions in the Shiraz and Mashhad Boxes
Extreme drought and wet conditions in the Shiraz and Mashhad boxes have been studied by the calculation of 1-month, 6-month, and 12-month SPI (Supplementary Figure S2)

Studying the Role of Marine Moisture Uptake Sources on Extreme Drought and Wet Conditions in the Shiraz and Mashhad Boxes
Extreme drought and wet conditions in the Shiraz and Mashhad boxes have been studied by the calculation of 1-month, 6-month, and 12-month SPI (Supplementary Figure S2) The influence of marine moisture sources on 1-month, 6-month, and 12-month SPI has been verified by multi-regression method in Shiraz and Mashhad boxes (Supplementary Table S1). The resulted regression equations showed the highest correlation coefficient (r = 0.72) for the 6 month SPI in the Shiraz box during the wet period. The results also demonstrated that among the various moisture sources, the Arabian Sea directly influences the extreme wet and drought conditions in both the Shiraz and Mashhad boxes in the wet period p-value (<0.05). However, during the dry period, none of the moisture uptake sources influenced drought conditions p-value (>0.05) in both boxes. The direct influence of the Arabian Sea moisture source on 1-month, 6-month, and 12-month SPI for Shiraz box and 1-month SPI for Mashhad box during the wet period is due to its dominant role in providing moisture for precipitation. During the dry period, the moisture sources have no influence on SPI values due to the intense atmospheric stability caused by Azores subtropical high-pressure system over Iran which prevents precipitation. In this case, local parameters influence Iran's precipitation mainly in the Zagros region and the Caspian Sea coasts [34,61].

Studying the Role of Climatic Teleconnections on Extreme Drought and Wet Conditions in the Shiraz and Mashhad Boxes
The role of climatic teleconnections on extreme drought and wet conditions in the Shiraz and Mashhad boxes and also Chenar-rahdar and Kardeh River discharges was studied using the simple linear regression method and least square method has been used to determine the best fit line. The models with the p-value (<0.05) show meaningful correlation between climatic teleconnection indices and (1-month, 6-month, and 12-month SPI) drought conditions, while the p-value (>0.05) demonstrated insignificant correlation. Extreme drought and wet conditions in Iran are not only influenced by local climatological components but also by large-scale climatic teleconnections. The results demonstrated that among BEST, SOI, NAO, and QBO indices (independent variables) only BEST and SOI indices directly influence 1-month, 6-month, and 12-month SPI and SSI (dependent variables) in the Shiraz and Mashhad boxes ( Table 3). The regression models showed the highest correlation coefficient (r = 0.3) between BEST indices and 6-month SPI in the Shiraz box. In Table 3 positive correlation was observed between BEST and SPI and SSI values, while negative correlation was observed between SOI indices and these drought indicators. The drought conditions often occurred during the negative phase of BEST and the positive phase of SOI (La Niña phase). However, the wet conditions often occurred during the positive phase of BEST and the negative phase of SOI indices (El Niño phase).  The temporal variations of monthly BEST and SOI climatic indices with 1-month, 6-month, and 12-month SPI and SSI for Shiraz and Mashhad boxes have also been shown in Figures 6 and 7.

Studying the Effect of the Dominant Moisture Sources on the Hydrology of Water Resources in the Shiraz and Mashhad Boxes
The temporal variations of moisture uptake rate from marine sources and Chenar-rahdar River discharges in the Shiraz box and Kardeh River discharges in the Mashhad box have been shown in Figure 8. During the wet period, the discharge variations of both rivers often follow the same trend as the total moisture uptake rate (E − P) > 0 (Figure 8a,c). This is because the discharges of both Chenar-rahdar River in the Shiraz box ( Figure 8a) and Kardeh River (Figure 8c) in the Mashhad box are dominantly supplied by precipitation originating from the mentioned marine moisture sources. During the dry period, moisture sources (E − P) > 0 do not show the same trend as the discharges of Chenar-rahdar and Kardeh rivers (Figure 8b,d) respectively and cannot lead to precipitation due to the intense atmospheric stability caused by Azores subtropical high-pressure system which prevents precipitation. Chenar-rahdar and Kardeh rivers are mainly supplied by groundwater base flow and to a lower extent by infrequent precipitation events during the dry period. The temporal variations of monthly BEST and SOI climatic indices with 1-month, 6-month, and 12-month SPI and SSI for Shiraz and Mashhad boxes have also been shown in Figures 6 and 7.  The simultaneous effect of each marine moisture source (E − P) > 0 on Chenar-rahdar and Kardeh rivers in the Shiraz and Mashhad boxes, respectively, has been investigated by the multi-regression model (Supplementary Table S2). Among the various moisture origins during the wet period, only the Arabian Sea moisture source directly influenced the discharges of Chenar-rahdar and Kardeh rivers. However, during the dry period, the Red Sea has an inverse correlation with the discharges of Chenar-rahdar and Kardeh rivers. During the dry period, when the Middle East and North Africa experience higher than average weather temperature, the (E − P) >0 values of the Red Sea notably intensify due to its huge evaporation and higher SST. These conditions reinforce the Azores subtropical high-pressure system and the atmospheric stability over Iran which results in the decrease of precipitation and also leads to the decline of river discharges.

Studying the Effect of the Dominant Moisture Sources on the Hydrology of Water Resources in the Shiraz and Mashhad Boxes
The temporal variations of moisture uptake rate from marine sources and Chenar-rahdar River discharges in the Shiraz box and Kardeh River discharges in the Mashhad box have been shown in Figure 8. During the wet period, the discharge variations of both rivers often follow the same trend as the total moisture uptake rate (E − P) > 0 (Figure 8a,c). This is because the discharges of both Chenar-rahdar River in the Shiraz box ( Figure 8a) and Kardeh River (Figure 8c) in the Mashhad box are dominantly supplied by precipitation originating from the mentioned marine moisture sources. During the dry period, moisture sources (E − P) > 0 do not show the same trend as the discharges of Chenar-rahdar and Kardeh rivers (Figure 8b,d) respectively and cannot lead to precipitation due to the intense atmospheric stability caused by Azores subtropical high-pressure system which prevents precipitation. Chenar-rahdar and Kardeh rivers are mainly supplied by groundwater base flow and to a lower extent by infrequent precipitation events during the dry period.  The simultaneous effect of each marine moisture source (E − P) > 0 on Chenar-rahdar and Kardeh rivers in the Shiraz and Mashhad boxes, respectively, has been investigated by the multi-regression model (Supplementary Table S2). Among the various moisture origins during the wet period, only the Arabian Sea moisture source directly influenced the discharges of Chenar-rahdar and Kardeh rivers. However, during the dry period, the Red Sea has an inverse correlation with the discharges of Chenar-rahdar and Kardeh rivers. During the dry period, when the Middle East and North Africa experience higher than average weather temperature, the (E − P) >0 values of the Red Sea notably intensify due to its huge evaporation and higher SST. These conditions reinforce the Azores subtropical high-pressure system and the atmospheric stability over Iran which results in the decrease of precipitation and also leads to the decline of river discharges.

Conclusions
This study revealed that the Arabian Sea and the Arabian and Mediterranean seas were predominant marine moisture sources in the south-central (Shiraz box) and northeastern (Mashhad box) parts of Iran during the wet period, while moisture uptake from the Persian Gulf and the Caspian Sea were dominant in Shiraz and Mashhad boxes, respectively during the dry period. The Arabian Sea moisture source dominantly influenced drought and wet conditions and also river discharges in Shiraz and Mashhad boxes during the wet period. However, during the dry period, due to the intense atmospheric stability caused by Azores subtropical high-pressure system which prevented precipitation, no correlation was observed. The large-scale climatic teleconnections BEST and SOI also influenced drought and wet conditions in Shiraz and Mashhad boxes. Extreme wet conditions mainly occurred during the El Niño phase, while extreme drought conditions mainly occurred during the La Niña phase. Scrutinizing the role of the marine moisture sources on Chenar-rahdar and Kardeh rivers discharge demonstrated that the Arabian Sea moisture source had a direct influence on these rivers during the wet period, while no correlation was observed during the dry period. These rivers were mainly supplied by precipitation during the wet period, while groundwater base flow and, to a lower extent, infrequent precipitation supplied their discharges during the dry period.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/10/11/ 1550/s1, Table S1: The regression analysis between marine moisture uptake sources (E − P) > 0 and 1-month, 6-month, and 12-month SPI and SSI hydrological drought index in the Shiraz and Mashhad boxes; Table S2: The multi-regression analysis between marine moisture uptake sources (E − P) > 0 and the discharges of Kardeh River in the Mashhad box and Chenar-rahdar River in the Shiraz box; Figure S1: (E − P) > 0 and precipitation correlation in Shiraz (a) and Mashhad (b) boxes during the wet and dry periods; Figure S2