Application of Statistical Techniques to Study Stable Isotopes ( 18 O and 2 H) Characteristics of Precipitation in Iran (Southwest Asia)

: Various climatic and geographic parameters inﬂuence precipitation in Iran, which makes the interpretation of stable isotope signatures in precipitation very complicated. Thus, precipitation sampling stations for stable isotope analyses in Iran have been classiﬁed by cluster analysis (CA) into 10 clusters, based on their stable isotope characteristics. The classiﬁcation of stations by CA also has a close correlation with the Koppen climatic zones across Iran. Finally, the stations in each cluster were plotted on the GMWL and EMMWL. This study shows that classifying precipitation sampling stations can simplify the interpretation of stable isotopes in the precipitation of regions with complicated climatic systems.


Introduction
Iran is a semi-arid and arid country in the Middle East that has historically faced water shortage crises. Although the annual precipitation in Iran is low (341 mm) [1], this amount of precipitation is also unevenly distributed across the country. In the Caspian Sea coastal area, mainly its western parts, the amount of annual precipitation is more than 1800 mm, while the precipitation notably decreases to less than 100 mm in the central part of Iran [1]. The dominant moisture sources and air masses causing precipitation also show large variations across Iran. During the wet and cold period (November to April), cP (continental polar), cT (continental tropical), mP (maritime polar), and MedT air masses influence Iran. However, during the dry and hot period (May to October), a mT (maritime tropical) air mass influences this country [1][2][3]. The cP air mass mainly transfers the moisture from the Caspian Sea, the Mediterranean Sea, and to the lesser extent the Black Sea to Iran. The cT air mass mainly transfers the moisture of the Arabian Sea, the Persian Gulf, the Red Sea, and the Oman Sea to Iran. The mP air mass transfers the moisture of the Black Sea and high latitude water bodies such as the North Atlantic Ocean to Iran. Finally, the mT air mass transfers the moisture of the Oman Sea and the Indian Ocean to Iran [4][5][6][7][8][9].
The integrated moisture fluxes over Iran from 1981 to 2015 are shown in Figure 1. The pattern of moisture fluxes over Iran during the dry period is different from that of the cold and wet period. During the dry period, strong moisture fluxes are observed over the Arabian Sea and the Indian Ocean, causing monsoons in India, which sometimes influence the southeast of Iran, Figure 1a. However, in the cold and wet period, moisture fluxes are observed over the Black Sea, the Mediterranean Sea, the Red Sea, and the Persian Gulf, Figure 1b. The large variations observed in the precipitation amount across Iran have a significant influence on the climate of this country. According to the Koppen classification of climate zones, Iran has been classified into various zones, including BWh (arid, desert, and hot), BWk (arid, desert, and cold), BSh (arid, steppe, and hot), BSk (arid, steppe, and cold), Csa (temperate, dry summer, and hot summer), Csb (temperate, dry summer, and warm summer), Cfa (temperate, no dry season, and hot summer), Dsa (continental, dry summer, and hot summer), Dsb (continental, dry summer, and warm summer), Dsc (cold, dry summer, and cold summer), Dfb (continental, no dry season, and warm summer), and Dfc (continental, no dry season, and cold summer) [11,12].
Due to Iran's low average annual precipitation amount, significant spatial variations of precipitation amount across the country, and the various air masses influence Iran consideration of all aspects of precipitation has a dominant role. Therefore, precipitation characteristics across Iran should be studied using accurate and reliable methods. The stable isotopes technique can present an accurate and comprehensive view of Iran's precipitation characteristics.
The application of stable isotopes in water resource studies has been utilized since Craig, in 1961 [13], found that 18 O/ 16 O has a very close and strong correlation with 2 H/H variations in fresh water molecules [14]. The stable isotopes technique presents crucial information regarding the climate condition of the moisture's origin, as well as the climate condition in precipitation sampling sites. There have been many studies regarding the application of stable isotopes in precipitation in Iran. Most of these studies were local small scale investigations [6,7,9,[15][16][17][18][19] or large scale studies that covered all of Iran [3,20,21]. The aim of this research was to perform a comprehensive study on the stable isotope characteristics of precipitation in Iran using statistical techniques.

Materials and Methods
In this study, stable isotopes ( 18 O and 2 H) and d-excess were studied in 34 stations across Iran. Stable isotope data are presented in delta notation (δ), which is the relative deviation of the sample from the standard (Vienna standard mean ocean water (VSMOW)) by Equation (1): The large variations observed in the precipitation amount across Iran have a significant influence on the climate of this country. According to the Koppen classification of climate zones, Iran has been classified into various zones, including BWh (arid, desert, and hot), BWk (arid, desert, and cold), BSh (arid, steppe, and hot), BSk (arid, steppe, and cold), Csa (temperate, dry summer, and hot summer), Csb (temperate, dry summer, and warm summer), Cfa (temperate, no dry season, and hot summer), Dsa (continental, dry summer, and hot summer), Dsb (continental, dry summer, and warm summer), Dsc (cold, dry summer, and cold summer), Dfb (continental, no dry season, and warm summer), and Dfc (continental, no dry season, and cold summer) [11,12].
Due to Iran's low average annual precipitation amount, significant spatial variations of precipitation amount across the country, and the various air masses influence Iran consideration of all aspects of precipitation has a dominant role. Therefore, precipitation characteristics across Iran should be studied using accurate and reliable methods. The stable isotopes technique can present an accurate and comprehensive view of Iran's precipitation characteristics.
The application of stable isotopes in water resource studies has been utilized since Craig, in 1961 [13], found that 18 O/ 16 O has a very close and strong correlation with 2 H/H variations in fresh water molecules [14]. The stable isotopes technique presents crucial information regarding the climate condition of the moisture's origin, as well as the climate condition in precipitation sampling sites. There have been many studies regarding the application of stable isotopes in precipitation in Iran. Most of these studies were local small scale investigations [6,7,9,[15][16][17][18][19] or large scale studies that covered all of Iran [3,20,21]. The aim of this research was to perform a comprehensive study on the stable isotope characteristics of precipitation in Iran using statistical techniques.

Materials and Methods
In this study, stable isotopes ( 18 O and 2 H) and d-excess were studied in 34 stations across Iran. Stable isotope data are presented in delta notation (δ), which is the relative deviation of the sample from the standard (Vienna standard mean ocean water (VSMOW)) by Equation (1): The analytical standard uncertainties for precipitation samples were ±0.1‰ and ±1‰ for δ 18 O and δ 2 H, respectively. Stable isotopes data have been gathered from PhD and MSc theses, as well as scientific papers.
As the number of studied stations across Iran is large, these stations were classified based on their stable isotope characteristics using cluster analysis (CA). Then, based on station classification by CA analysis and Koppen climatic zones across Iran, the stable isotope characteristics of precipitation in Iran were studied.

Results and Discussion
In this study, stable isotope signatures in 34 sampling stations across Iran were studied, as shown in Table 1. Studying the variations of stable isotopes in precipitation across Iran is difficult due to the various climate zones that exist in this country. Therefore, CA was used to classify the studied stations based on their stable isotope characteristics. CA classified the studied stations into 10 groups, see Table 2. Studying the spatial distribution of clusters across Iran and comparing it with the Koppen climatic zones map presented valuable results, Figure 2.
A significant number of the studied stations are in the first cluster, which are located in the arid, steppe, and cold region (BSK) climatic zone. These stations are located in the southern and southwestern parts of the Zagros region. Precipitation in these stations occurs by moisture supply from the Persian Gulf, the Red Sea, and the Arabian Sea via cT and MedT air masses. The average stable isotope contents of δ 18 O and δ 2 H in the stations located in this cluster are −5.34‰ and −23.51‰, respectively. Furthermore, the d-excess also shows high values of 19.30‰.
In the second cluster, the studied stations are located in the central part of Iran and classified as BSk and BWk climate zones. These stations receive a lower amount of annual precipitation compared to the other stations in Iran. In this cluster, the role of moisture originating from low latitude water bodies such as the Persian Gulf, the Red Sea, and the Arabian Sea is significant. The very high d-excess values in the stations located in this cluster (22.75‰) also confirm this. However, the δ 18 O and δ 2 H in the stations located in this cluster show depleted values of −6.95‰ and −33.30‰, respectively. Iran and comparing it with the Koppen climatic zones map presented valuable results, Figure 2.  A significant number of the studied stations are in the first cluster, which are located in the arid, steppe, and cold region (BSK) climatic zone. These stations are located in the southern and southwestern parts of the Zagros region. Precipitation in these stations occurs by moisture supply from the Persian Gulf, the Red Sea, and the Arabian Sea via cT and MedT air masses. The average stable isotope contents of δ 18 O and δ 2 H in the stations located in this cluster are −5.34‰ and −23.51‰, respectively. Furthermore, the d-excess also shows high values of 19.30‰.
In the second cluster, the studied stations are located in the central part of Iran and classified as BSk and BWk climate zones. These stations receive a lower amount of annual precipitation compared to the other stations in Iran. In this cluster, the role of moisture The stations in the third cluster are also classified in the BSk and BWk climate zones. In addition to the cT and MedT air masses, which influence precipitation in this region, cP air mass also has a dominant role in the stations located in this cluster. The cP air mass significantly transfers the moisture of water bodies with a low sea surface temperature (SST), such as the Caspian Sea, to these stations. The average δ 18 O and δ 2 H values in precipitation in the stations in this cluster are −6.90‰ and −42.03‰, respectively. The low d-excess values (13.40‰) in the studied stations in this cluster also confirm the role of moisture originating from low SST water bodies, such as the Caspian Sea.
The fourth cluster consists of the Shahroud, Gorgan, and Rasht stations. The Rasht station is located in the Cfa, Gorgan is located in BSk, and Shahroud is located in BWk climate zones according to the Koppen classification. Rasht, in the western part of the Caspian Sea central area, receives the highest amount of precipitation compared to the other stations across Iran, while Shahroud, located in an arid zone, receives a very low amount of annual precipitation.
In the fourth cluster, the role of the cP air mass and moisture originating from the Caspian Sea is dominant in providing moisture for precipitation. The very low d-excess values (7.70‰) in the precipitation of the stations in this cluster also confirm its origin in the moisture supply from the Caspian Sea. The average δ 18  Rafsanjan is the only station in the seventh cluster that is located in the central part of Iran in the BWh climate zone. The central part of Iran receives a very low amount of precipitation compared to the other parts of Iran. Two large deserts, "Dashteh-Lut" and "Dashteh-Kavir", also exist in this part of Iran. Rafsanjan receives moisture from low latitude water bodies such as the Persian Gulf, the Red Sea, and the Arabian Sea, characterized by high SST. The average δ 18 O and δ 2 H values in this station are −0.5‰ and 12.90‰, respectively. However, the d-excess in the precipitation events in this station shows high values (16.90‰), which confirms their origin from low latitude water bodies.
In the eighth cluster, the Hamadan and Zarivar stations are located in the western part of Iran, in the Zagros region. These stations are mainly under the influence of the MedT air mass and the Mediterranean Sea, with high SST. The average δ 18 O and δ 2 H values in this cluster of stations are −8.45‰ and −53.40‰, respectively. The average d-excess in these stations shows a high value of 14.10‰.
In the ninth cluster is only the Hashtgerd station. The role of the cP air mass and moisture originated from the Caspian Sea is dominant in providing moisture for the precipitation over this station. This is the reason for the very low d-excess value (2.50‰) in this station. The average δ 18 O and δ 2 H values in the station in this cluster are −6.50‰ and −49.70‰, respectively. Although the Hashtgerd station is located beside third cluster stations, it has been classified in a separate cluster (ninth cluster). This is because Hashtgerd is in the boundary zone, between climatic zones BWk and BSk.
In the last cluster (tenth cluster) is only the Sabalan station. This station is located in a very cold zone with a high amount of annual precipitation. This cluster shows the most depleted isotope values compare to other stations, Table 2. It is mainly under the influence of the cP, MedT, and mP air masses and moisture originating from the Mediterranean Sea, Black Sea, and to a lesser extent the Caspian Sea. The d-excess in the precipitation events of this station show high values (17.10‰).

Studying Stable Isotope Signatures of Precipitation in Iran Based on the Koppen Climatic Zones
The stable isotopes in precipitation events across Iran are under the influence of various local and regional parameters. Plotting the average δ 18 O and δ 2 H values in the studied stations based on their clusters on the global meteoric water lines (GMWL) and Eastern Mediterranean meteoric water line (EMMWL) demonstrate valuable results, Figure 3.
tions, it has been classified in a separate cluster (ninth cluster). This is because Hashtgerd is in the boundary zone, between climatic zones BWk and BSk.
In the last cluster (tenth cluster) is only the Sabalan station. This station is located in a very cold zone with a high amount of annual precipitation. This cluster shows the most depleted isotope values compare to other stations, Table 2. It is mainly under the influence of the cP, MedT, and mP air masses and moisture originating from the Mediterranean Sea, Black Sea, and to a lesser extent the Caspian Sea. The d-excess in the precipitation events of this station show high values (17.10‰).

Studying Stable Isotope Signatures of Precipitation in Iran Based on the Koppen Climatic Zones
The stable isotopes in precipitation events across Iran are under the influence of various local and regional parameters. Plotting the average δ 18 O and δ 2 H values in the studied stations based on their clusters on the global meteoric water lines (GMWL) and Eastern Mediterranean meteoric water line (EMMWL) demonstrate valuable results, Figure 3. The stations located in the Zagros mountain range in the first, second, sixth, and seventh clusters mainly plot on the EMMWL. This is because the precipitation events in the stations located in this part of Iran are mainly provided by the moisture originating from the Mediterranean Sea, as well as from other water bodies such as the Persian Gulf, the Arabian Sea, and the Red Sea. The Shirin bahar station in the sixth cluster shows significant deviation from both GMWL and EMMWL, due to the huge secondary evaporation effect. Although the stations located in the fifth cluster are also located in the Zagros region, a very mild deviation can also be observed from EMMWL. This is due to the secondary evaporation of precipitation events in the stations that exist in this cluster.
On the other hand, stations located in the northern part of Iran, including the ninth, third, and fourth clusters, plot on or in the vicinity of the GMWL. There are only two exceptions, the Rasht and Gorgan stations in the fourth cluster. The Rasht station plots on the EMMWL, which is due to the effect of local events, such as the numerous lagoons that influence stable isotope contents of precipitation in this area. The Gorgan station deviates The stations located in the Zagros mountain range in the first, second, sixth, and seventh clusters mainly plot on the EMMWL. This is because the precipitation events in the stations located in this part of Iran are mainly provided by the moisture originating from the Mediterranean Sea, as well as from other water bodies such as the Persian Gulf, the Arabian Sea, and the Red Sea. The Shirin bahar station in the sixth cluster shows significant deviation from both GMWL and EMMWL, due to the huge secondary evaporation effect. Although the stations located in the fifth cluster are also located in the Zagros region, a very mild deviation can also be observed from EMMWL. This is due to the secondary evaporation of precipitation events in the stations that exist in this cluster.
On the other hand, stations located in the northern part of Iran, including the ninth, third, and fourth clusters, plot on or in the vicinity of the GMWL. There are only two exceptions, the Rasht and Gorgan stations in the fourth cluster. The Rasht station plots on the EMMWL, which is due to the effect of local events, such as the numerous lagoons that influence stable isotope contents of precipitation in this area. The Gorgan station deviates from both the GMWL and EMMWL due to the evaporation effect. Finally, the Sabalan station that is located at high elevation in the Sablan mountain hills in the northern part of Iran demonstrates high d-excess values and plots in the vicinity of the EMMWL. This is because this region is dominantly under the influence of moisture transported from the Mediterranean Sea.

Conclusions
This study shows that stable isotope signatures in precipitation across Iran are significantly under the influence of the climatic conditions of the sampling stations, as well as precipitation moisture source regions. Therefore, the precipitation characteristics show spatial variations across the country. Using cluster analysis (CA), the studied stations were classified into 10 clusters according to the stable isotope characteristics of precipitation. Stations located in each cluster show some differences in δ 18 O and δ 2 H values and d-excess. Furthermore, the spatial variations of stable isotope signatures in precipitation have been linked to the climatic zones where the studied stations are located. The method used in this study can also be applied in other parts of the world to classify precipitation sampling stations based on their stable isotope characteristics.