Investigation of Aeolian Dust Deposition Rates in Different Climate Zones of Southwestern Iran

: Dust and atmospheric particles have been described in southwestern Iran primarily in terms of load, concentration and transport. The passive deposition, however, has been discussed inadequately. Therefore, the relationships between different climate zones in southwestern Iran and dust deposition rates were quantiﬁed between 2014 and 2017 using both space- (second modern-era retrospective analysis for research and applications, version 2 reanalysis model) and ground-based (eolian ground deposition rate) tools. In addition, the surface meteorological records, including the wind patterns favoring the occurrence of dust events, were examined. A hot desert climate (BWh), hot semi-arid climate (BSh), and temperate hot and dry summer climate (Csa) were identiﬁed as the three dominant climate regions in the study area, exhibiting the highest average dust deposition rates. In this study, correlations between the most relevant climate patterns and deposition rate weather parameters were found to describe a region’s deposition rate when a dust event occurred. Based on these results, the BSh and Csa regions were found to be associated with the seasonal cycle of dust events in March, April, and May, revealing that in the long run meteorological conditions were responsible for the varying dust deposition rates. Relatively, precipitation and temperature were the two major factors inﬂuencing dust deposition rates, not wind speed. Moreover, the peak seasonal deposition rates in the spring and summer were 8.40 t km − 2 month − 1 , 6.06 t km − 2 month − 1 , and 3.30 t km − 2 month − 1 for the BWh, BSh, and Csa climate regions, respectively. However, each of these climate types was directly related to the speciﬁc quantity of the dust deposition rates. Overall, the highest dust deposition rates were detected over the years studied were 100.80 t km − 2 year − 1 , 79.27 t km − 2 year − 1 , and 39.60 t km − 2 year − 1 for BWh, BSh, and Csa, respectively. of spring > summer, the BWh region results did not ﬁt this rule. Rather, its dominant deposition rates occurred in the order of summer > spring for 2014 and 2017. These results show that the lowest dust deposition rates occurred during the fall and winter, while the highest occurred in the the and ground (GDR) conducted for a time corresponding velocity near the surface. an agreement between m m s observed


Introduction
Over hundreds of years, dislocated eolian dust materials create positive and negative feedback between dust particles in the atmosphere and climate, which can influence terrestrial ecosystems [1]. Fertilized dust provides essential nutrients for plant growth, causing the revegetation of desert surfaces [2,3]. However, eolian dust also increases soil salinity, reduces air quality, and impairs human health e.g., [4][5][6][7].
Globally, eolian dust transport cycles and dust loads are closely connected to climate change e.g., [8,9] and have reciprocal effects. On the one hand, variations in meteorological episodes influence dust concentration [10,11], while on the other hand, the amount of eolian dust suspended in the atmosphere influences the climate, e.g., [12].
In particular, particles can travel over long distances before they are removed from the atmosphere by dry and/or wet deposition [13]. Thus, a particle's longevity in the atmosphere alters the radiative balance by scattering and absorbing incoming solar radiation [14] or changing the optical properties of clouds, e.g., [15]. Changing the long-term optical properties may push the climate system toward warmer and wetter conditions or amplify sudden changes in the climate, e.g., [16][17][18]. A multi-model climate study showed that climate change is correlated with the negative impacts of eolian dust in the atmosphere, which is pronounced by a decrease in the mean precipitation over land in warmer climates, thus increasing the lifetime of aerosols with a lower wet deposition rate, e.g., [19]. In both scenarios, airborne dust is slowly removed from the atmosphere to the Earth's surface at a certain dry deposition rate.
According to [20], the annual atmospheric dust load originating in the Sahara Desert, the Arabian Peninsula, Central Asia, Northern China, and Australia is estimated to be more than 200 million tons. This high concentration affects both the air quality [21,22] and global climate, e.g., [17,23]. Ref. [23] suggested that major sources of dust in the Middle East can be identified by analyzing the surface meteorological records at each station as follows: DE = DIS + BD + DS + SDS. (1) From 1998 to 2003, the synoptic records of land surface observations were classified into four definitions based on the maximum number of dust events (DE) under the following conditions [24]: • Dust in suspension (DIS), which is identified as widespread dust particles in suspension that influence visibility, limiting it to approximately 10 km.

•
Blowing dust (BD) is defined as a raised dust particle that reduces visibility to 1-10 km. • Dust storms (DS) occur when a large amount of dust particles is raised, reducing visibility to between 200 m and 1000 m. • Severe dust storms (SDS) represent large quantities of lifting dust particles that reduce visibility to less than 200 m.
In the Middle East, which is one of the regions most affected by dust, the Arabian Peninsula ranks in the global top five of the most significant dust source regions, e.g., [23,24]. Ref. [25] examined five years of data, from 2009 to 2014, finding that the southwestern regions of Iran were most affected by dust, characterized by annual mean concentrations of more than 300 µg/m 3 of suspended atmospheric particles. [26], defined the same value as 372 µg/m 3 in the south (Ahwaz) and approximately 229 µg/m 3 in the west (Kermanshah).
Similarly, the deposition rate is a comparatively essential point for identifying eolian activity and environmental quality, e.g., [27,28]. In contrast to atmospheric contributions, few experiments have measured the deposition rate from ground observations (Table 1).  [6] more, the influence of climate zones on ground deposition rates has hardly been utilized in this region. Therefore, this study provides new and valuable insights into the spatial and temporal distributions of dust deposition correlated to climate factors in southwestern Iran.
The main objectives of this study are: (1) to investigate the spatial and temporal variability of dust deposition rates, (2) to determine the seasonal and spatial variability of the dust deposition rate in relation to the climate zones, and (3) to identify the correlation between the most relevant climate patterns and weather parameters to the deposition rate.

Study Area
The study area is located between 35°30′00″ N, 45°30′00″ E and 30°30′00″ N, 49°30′00″ E in southwestern Iran, covering an area of approximately 117,000 km 2 , as shown in Figure 1 [44]. The eastern part of this area is dominated by the Zagros mountain range (running from the northwest to the south). Because of this geographic setting, the study area occupies a broad latitudinal range of climates. In general, the climate of the study area is characterized as semi-dry hot and occasionally humid, which is especially common in the south, while the region's winters can be cold and dry [45]. The area has great potential for agricultural development because of the permanent river flows that contribute to the area's fertility. The observed seasonal dust events indicate that southwestern Iran is strongly affected by annual dust events from the large deserts in Iraq, Saudi Arabia, and Syria, especially during the hot season, e.g., [46][47][48]. Thus, southwestern Iran is recognized as a global hotspot of eolian dust.  [39]. Top; climate map customized for study area reproduced courtesy of [49]. G01 to G10 shows the distribution of dust collectors. Down: Map of the study area, including latitude, longitude, elevation, and the main tectonic elements of the region. Reproduced courtesy of the U.S. Geological Survey data explorer [50].  Figure 1 shows the distribution of ground dust collectors (G01-G10) within the Köppen climate classification system in the study area [49]. The study area was classified on the basis of standard climate zone properties of the Köppen climate system, which has also been used in previous research [45,51]. In the Köppen climate classification scheme, a climate is defined with three letters, the first of which represents a world climate at a 1 km resolution, of which there are five different groups [49]: A: tropical, B: arid (dry), C: warm (temperate), D: continental (cold), and E: polar. While most of these groups are classified based on the temperature of a region, group B emphasizes the aridity of a region. The second letter in the Köppen classification is the seasonal precipitation type, represented by steppe (S) desert (W) rain force without dry season (f), monsoon (m), savanna dry winter (w), and dry summer (s). Finally, the third letter indicates the level of heat as cold (k), hot (h), and, in general, hot summer (a), warm summer (b), cold summer (c), and very cold winter (d). Herein, the climate of the study area was divided into three main climate zones ( Table 2). Climate of study area was classified as hot desert climate (BWh), hot semi-arid climate (BSh), and temperate hot and dry summer climate (Csa) using the Köppen climate classification scheme [45], wherein each letter represents one of five main climate group at a 1 km resolution. Note that "T" is temperature, "P" is precipitation, and " C" denotes Celsius. P is paired with subscripts of "ann" annually, "th" threshold , "wmin, wmax" winter minimum and maximum, respectively, and "smin, smax" summer minimum and maximum, respectively.

Climates of Study Area
As shown in Figure 1, the study area was categorized into three main climate groups (B, C, and D). The dust event history and availability of consistent monthly data collection (access to data for approximately four years) are the main issues during the research period. Samplers were constructed and placed in the area to enable more exposure to the frequent dust events. According to the definition of the study area setting, the high-altitude areas of the Dsa and Dsb climate zones (belonging to the D group) did not contribute to the data collection process, and studies on climatic regions belonging to the D group were excluded.
The hot desert climate (BWh) in the southern region of the northern coastal plains of the Persian Gulf covers the locations of G08, G09, and G10 between latitudes of 30 N and 31 N. A hot semi-arid climate zone (BSh) in the western areas is indicated over the locations of G01, G02, and G07, which are located between 34 N and 32 N. The temperate hot summer zone (Csa) is located between 33 N and 34 N from the extreme western regions of G04 and G05 to G06.
According to the recent land use and land cover (LULC) pattern and in accordance with the Land Use/Cover Area framework of the Statistical Survey (LUCAS), the study area has five major classifications: artificial (6.5%), industrial (0.05%), wet land (4.5%), vegetation (27%), and bare land (61%). As the deposition rate may have been provided using passive samplers (gravimetric method of sampling), the samples should not be influenced by artificial or natural obstacles. Therefore, the samplers were placed at unhindered areas in a radius of 1 km and 2 m high in order to minimize obstacle impacts. Dust particles were collected before settling on the surface. Detailed information on the LULC of each station is provided by [52].

Methods and Data
To identify the impact of the climate on the deposition rates in southwestern Iran, the ground-based deposition rate was recorded and classified into three climate zones in the study area. Wind, precipitation, and temperature are important for determining the synoptic conditions responsible for the dust deposition rate. Thus, the maximum, minimum, and average values for all the climate zone-based data derived from the operating meteorological stations were calculated by the I.R. Iran Meteorological Organization. The data used include the monthly average ground deposition rate (GDR), wind speed and direction (northward and eastward winds), precipitation, and temperature during the study period of 2014-2017. To evaluate the correlation between the deposition rates and climate zones, the results of ten dust sampling sites were reclassified using climate gradients. The rates distinguished with seasonal variation responses to the regional climate provided using the second modern-era retrospective analysis for research and applications (MERRA-II) model output had reasonable accuracy. By evaluating the relationship between high temporal variation (model output) and seasonal rate (from GDR), the spatiotemporal pattern of the dust deposition rate over different climate zones revealed the probable uncertainty of the deposition rate.

Ground Survey Setup
Both gravimetric and directional dust samplers (total of 20) were deployed to satisfy the requirements of ASTM international D1356-05 which is formaly known as American Society for Testing and Material [53], and began operation in 2014. Each sampler was mounted 2 m above the ground surface and kept away from possibly anthropogenic impact (traffic, human activities, and commotion). The dust was collected passively on a single piece of paper (waterproof) with a diameter of 20 cm. The directional samplers collected dust from the wind direction to the plate, while gravimetric samplers collected particles returning to the surface by gravity.
As summarized in Table 3, among the 10 ground stations, G03 and G06 appeared on the boundary of two climate zones (BSh and Csa). According to the map scale and climate resolution of 1 km, the stations belong to Csa, which encompasses a range of topographic and meteorological conditions. The detailed meteorological data and GDR datasets from the stations in southwestern Iran were analyzed to assess the seasonality of the dust deposition rates. Dust samples are labeled from G01 to G10. The climate ranges are classified according to Köppen (1936) [49].
Atmosphere 2021, 12, 229 6 of 23 As shown in Figure 2, the monthly dust collection was weighted and classified according to the climate zone. Data from the weighted GDRs were used to investigate the spatiotemporal distribution of dust events. Dust samples are labeled from G01 to G10. The climate ranges are classified according to Köppen (1936) [49].
As shown in Figure 2, the monthly dust collection was weighted and classified according to the climate zone. Data from the weighted GDRs were used to investigate the spatiotemporal distribution of dust events.

MERRA-II Application
As of late, several new reanalysis applications, including those from the National Center for Environmental Prediction, North American Regional Reanalysis, was released by NASA [54]. In particular, in 1980, the second modern-era retrospective analysis for research and applications (MERRA-II) reanalysis tool was introduced. It was called M2TMNXSLV, and replaced the original MERRA reanalysis tool with an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system with a numerical model. For three main variables, including aerosols, atmospheric dynamic measurements, and dust, meteorological observations can be obtained. The MERRA-II has a spatial resolution of 0.5° × 0.025° at 550 nm, and hourly to monthly temporal data are available. As these data were integrated from a variety of satellite systems with numerical models in order to produce a temporally and spatially consistent composite of the climate [54,55], MERRA-II is ideal for investigating climate variability. Given

MERRA-II Application
As of late, several new reanalysis applications, including those from the National Center for Environmental Prediction, North American Regional Reanalysis, was released by NASA [54]. In particular, in 1980, the second modern-era retrospective analysis for research and applications (MERRA-II) reanalysis tool was introduced. It was called M2TMNXSLV, and replaced the original MERRA reanalysis tool with an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system with a numerical model. For three main variables, including aerosols, atmospheric dynamic measurements, and dust, meteorological observations can be obtained. The MERRA-II has a spatial resolution of 0.5 • × 0.025 • at 550 nm, and hourly to monthly temporal data are available. As these data were integrated from a variety of satellite systems with numerical models in order to produce a temporally and spatially consistent composite of the climate [54,55], MERRA-II is ideal for investigating climate variability. Given this dynamic, the MERRA-II dataset using GEOS-5 (version 5.12.4) can also function extensively with an expanding array of application areas, such as ecological forecasting, air quality, wind flow patterns, meteorological observations, and dust, as part of its routine reports [55][56][57]. Because various parameters of MERRA-II were evaluated, validating its results, the dry deposition values obtained using MERRA-II with latitude and longitude of the climate zone were obtained for each climate area using the NASA official website [58]. Further, the MERRA-II results were used to evaluate the correlation between the dry deposition rates and GDR during the period 2014-2017.

Results
This section presents the results of the weather pattern analysis for the study area, as well as the ground-based and MERRA-II deposition rates. After a detailed description of the meteorological and dust deposition rate data in Sections 3.1 and 3.2, their connection was analyzed for further discussion in Section 4.  Figure 3 illustrates the average interannual monthly temperatures of the three climate zones (BWh, BSh, and Csa). The data compare the temperatures in each climate zone over time.

Results
This section presents the results of the weather pattern analysis for the study area, as well as the ground-based and MERRA-II deposition rates. After a detailed description of the meteorological and dust deposition rate data in Sections 3.1 and 3.2, their connection was analyzed for further discussion in Section 4. Figure 3 illustrates the average interannual monthly temperatures of the three climate zones (BWh, BSh, and Csa). The data compare the temperatures in each climate zone over time. The fluctuating trend with the highest standard deviation (Std) was recorded for the temperate-hot-summer zone (13 °C < Csa < 14 °C). Meanwhile, the steppe hot zone temperature Std remained between 12 °C and 13 °C. Moreover, the desert-hot zone had a smaller Std than that of the other two zones with an Std range of 11 °C < BWh < 13 °C.

Meteorological Setting
As shown in Figure 4, the precipitation and average temperatures of the three climate zones were calculated from 2014 to 2017. The seasonal distribution of the precipitation data showed that the highest amounts of rainfall were detected during the winter and spring, while the summer months received very little precipitation. Moreover, the average The fluctuating trend with the highest standard deviation (Std) was recorded for the temperate-hot-summer zone (13 • C < Csa < 14 • C). Meanwhile, the steppe hot zone temperature Std remained between 12 • C and 13 • C. Moreover, the desert-hot zone had a smaller Std than that of the other two zones with an Std range of 11 • C < BWh < 13 • C.
As shown in Figure 4, the precipitation and average temperatures of the three climate zones were calculated from 2014 to 2017. The seasonal distribution of the precipitation data showed that the highest amounts of rainfall were detected during the winter and spring, while the summer months received very little precipitation. Moreover, the average annual precipitation in the hot desert climate (BWh) was between 107 mm in 2015 and 228 mm in 2014. That of the hot steppe area (BSh) was considerably higher, typically around 540 mm year −1 , with the exception of 2015, which had average precipitation of only 228 mm. The warm temperate climate zone (Csa) followed the same trend with belowaverage precipitation in 2015 (381 mm), while the precipitation in 2014 reached 667 mm. This temporal pattern was also detected by NOAA researchers [59], who identified 2015 as the warmest year on record for Asia and South America, with a global temperature deviation of 0.90 • C above average for the period 1901-2000. New monthly records were set in each month of 2015, with the exceptions of January and April In addition, the summer of 2015 was the driest season in central Europe from 1901 to 2015 [60]. Conversely, in monsoon-influenced regions, more extreme precipitation evens related to the temperature anomaly were recorded [61]. Recently, the global total precipitation from January 2015 to December 2018 was compared with the data from 1951 to 2010 [62] revealing that the region east of the Persian Gulf (including the study area) was, among other regions (North America, northern South America, Europe, northeast Australia, and southern and south-west Africa), characterized by a significant negative precipitation anomaly during recent years. This indicates a negative correlation between summer precipitation and air temperature, although it is non-linear. of 0.90 °C above average for the period 1901-2000. New monthly records were set in each month of 2015, with the exceptions of January and April In addition, the summer of 2015 was the driest season in central Europe from 1901 to 2015 [60]. Conversely, in monsooninfluenced regions, more extreme precipitation evens related to the temperature anomaly were recorded [61]. Recently, the global total precipitation from January 2015 to December 2018 was compared with the data from 1951 to 2010 [62] revealing that the region east of the Persian Gulf (including the study area) was, among other regions (North America, northern South America, Europe, northeast Australia, and southern and south-west Africa), characterized by a significant negative precipitation anomaly during recent years. This indicates a negative correlation between summer precipitation and air temperature, although it is non-linear. Wind characteristics are also highly important synoptic parameters as they heavily influence the temporal distribution and intensity of dust events and dust deposition rates. Herein, clear wind data distinctions were found between the three climate zones. Figure 5 shows the frequency of the average monthly wind speeds over all four years of observation. The overall mean wind speeds for BWh, BSh, and Csa were 6.09 m s −1 , 4.55 m s −1 , and 4.67 m s −1 , respectively. Note that the data distribution is skewed toward lower wind speeds, suggesting that the most frequent wind speeds are below the mean values. The Stds of the wind speed variation are representative of the unsteady wind speeds in BWh, BSh, and Csa, which are 1.09 m s −1 , 0.90 m s −1 , and 0.81 m s −1 , respectively. Wind characteristics are also highly important synoptic parameters as they heavily influence the temporal distribution and intensity of dust events and dust deposition rates. Herein, clear wind data distinctions were found between the three climate zones. Figure 5 shows the frequency of the average monthly wind speeds over all four years of observation. The overall mean wind speeds for BWh, BSh, and Csa were 6.09 m s −1 , 4.55 m s −1 , and 4.67 m s −1 , respectively. Note that the data distribution is skewed toward lower wind speeds, suggesting that the most frequent wind speeds are below the mean values. The Stds of the wind speed variation are representative of the unsteady wind speeds in BWh, BSh, and Csa, which are 1.09 m s −1 , 0.90 m s −1 , and 0.81 m s −1 , respectively. Westward winds in Figure 6 showed a similar seasonal variation pattern in the BWh and BSh climate zones, with maximum wind speeds of 3.5 m s −1 and 1.8 m s −1 , respectively, during the summer and minimum speeds during the winter and spring. Meanwhile, the Csa climate was characterized by maximum eastward wind speeds of 1.2 m s −1 in the summer and minima in the spring and winter. Further, the northward wind speed over the Csa region exhibits bimodal seasonal variation with two maxima, one during the spring during the summer and minimum speeds during the winter and spring. Meanwhile, the Csa climate was characterized by maximum eastward wind speeds of 1.2 m s −1 in the summer and minima in the spring and winter. Further, the northward wind speed over the Csa region exhibits bimodal seasonal variation with two maxima, one during the spring and another during the winter, while the minima were recorded in fall and summer. Conversely, regions with southward wind speeds revealed a complex temporal dynamic with double peaks during the spring and summer of 2015 and 2016, and during the summer and fall of 2014. During summer in 2017, however, a unimodal seasonal variation with one maximum was observed. Westward winds in Figure 6 showed a similar seasonal variation pattern in the BWh and BSh climate zones, with maximum wind speeds of 3.5 m s −1 and 1.8 m s −1 , respectively, during the summer and minimum speeds during the winter and spring. Meanwhile, the Csa climate was characterized by maximum eastward wind speeds of 1.2 m s −1 in the summer and minima in the spring and winter. Further, the northward wind speed over the Csa region exhibits bimodal seasonal variation with two maxima, one during the spring and another during the winter, while the minima were recorded in fall and summer. Conversely, regions with southward wind speeds revealed a complex temporal dynamic with double peaks during the spring and summer of 2015 and 2016, and during the summer and fall of 2014. During summer in 2017, however, a unimodal seasonal variation with one maximum was observed. To further illustrate the seasonal differences within and between the three climate zones, seasonal wind roses were designed using the frequency of different wind speeds as a percentage of the total wind ( Figure 7). The predominant seasonal wind directions in the BWh are along the northwest axis in every season, with average frequencies of 54% (winter), 52% (spring), 50% (summer), and 34% (fall). During summer and fall, southern To further illustrate the seasonal differences within and between the three climate zones, seasonal wind roses were designed using the frequency of different wind speeds as a percentage of the total wind ( Figure 7). The predominant seasonal wind directions in the BWh are along the northwest axis in every season, with average frequencies of 54% (winter), 52% (spring), 50% (summer), and 34% (fall). During summer and fall, southern and southeastern wind directions became more dominant, reaching up to 25%. The average seasonal wind speed in the BWh was the lowest in the fall (2.0 m s −1 ) and winter (2.3 m s −1 ) and the highest in spring (2.9 m s −1 ) and summer (3.0 m s −1 ).
In all seasons, the dominant wind direction in the Csa climate zone was southeast, with average frequencies of 40%, 40%, 39%, and 35% in the winter, spring, summer, and fall, respectively. Calms and events with no discernable wind direction occurred around the year, with frequencies of 18.5% (winter), 18.5% (spring), 9.7% (summer), and 31.8% (fall). Further, the Csa wind rose plots illustrated that the average seasonal wind speed was the lowest during the fall (1.8 m s −1 ) and highest (2.5 m s −1 ) during the spring and summer, while winter had an average of 2.3 m s −1 .
The predominant seasonal wind direction over the BSh was also along the southeast axis, with frequencies of 52% and 40% in the fall and winter, respectively. However, the most frequent single wind directions were detected during the spring (east: 22.5%; west: 20%) and summer (west-northwest: 50%; northwest: 10%). The average wind speeds were 1.9 m s −1 , 2.0 m s −1 , 2.3 m s −1 , and 3.4 m s −1 during the fall, winter, spring, and summer, respectively. Calms and events with no discernable wind direction were recorded in every season at frequencies of 22.2% (winter), 26.4% (spring), 10% (summer), and 28.9% (fall). summer, while winter had an average of 2.3 m s −1 .
The predominant seasonal wind direction over the BSh was also along the southeast axis, with frequencies of 52% and 40% in the fall and winter, respectively. However, the most frequent single wind directions were detected during the spring (east: 22.5%; west: 20%) and summer (west-northwest: 50%; northwest: 10%). The average wind speeds were 1.9 m s −1 , 2.0 m s −1 , 2.3 m s −1 , and 3.4 m s −1 during the fall, winter, spring, and summer, respectively. Calms and events with no discernable wind direction were recorded in every season at frequencies of 22.2% (winter), 26.4% (spring), 10% (summer), and 28.9% (fall).

Ground Dust Deposition Rate
The GDR was determined for each station used for meteorological measurements during the same study period. The temporal variability of the GDR is shown in Figure 8. The average dust deposition rates were 0.35 mg cm −2 month −1 , 0.64 mg cm −2 month −1 , and 0.86 mg cm −2 month −1 for the Csa, BSh, and BWh climate zones, respectively, during the period 2014-2017. Further, the Csa, BSh, and BWh zones had dispersion values of 0.23 mg Figure 7. Seasonal wind-roses for BWh, Csa, and BSh climate zones, illustrating the percentage of time that winds blow from a particular direction at a certain speed. The wind rose maps contain centroid circles, wherein the small circle in the center is the calm value in percent of total wind speed, which occurs when speeds are less than 2.0 m s −1 . The centroid circles (track) relative to the eight wind sectors represent wind severity and directions. Different color arrows signify different wind speeds, as shown in the legend.

Ground Dust Deposition Rate
The GDR was determined for each station used for meteorological measurements during the same study period. The temporal variability of the GDR is shown in Figure 8. The average dust deposition rates were 0.35 mg cm −2 month −1 , 0.64 mg cm −2 month −1 , and 0.86 mg cm −2 month −1 for the Csa, BSh, and BWh climate zones, respectively, during the period 2014-2017. Further, the Csa, BSh, and BWh zones had dispersion values of 0.23 mg cm −2 month −1 , 0.34 mg cm −2 month −1 , and 0.36 mg cm −2 month −1 , respectively.

MERRA-II Deposition Rate
The MERRA-II model application data were downloaded and prepared according to the study area's coordinate system in order to determine the time series of the monthly deposition rate against seasonal variation. Figure 9 shows the monthly and annual average deposition rates. The maximum deposition rates occurred during the spring and summer, whereas the lowest values were observed in the fall and winter. In particular, the average dust deposition rate in BWh fluctuated between 4.86 t km −2 year −1 and 5.61 t km −2 year −1 , while in the steppe hot climate (BSh), the year-to-year deposition rate varied on average from 2.96 t km −2 year −1 to 3.78 t km −2 year −1 . The temperate hot summer zone (Csa) was characterized by an average deposition rate of less than 2.05 t km −2 year −1 , or less than 50% that of the hot desert zone. An important feature of the seasonal variation patterns for each climate zone is that their trends have the same seasonality, but at different deposition rate levels.
G06 from Csa. The GDRs are weighted monthly in mg cm −2 from data collected onto a plate with a diameter of 20 cm. Dotted lines show average deposition rates (black dotted line) and standard deviation (brown dotted line). Diagonal lines is created for each climate region (BWh, BSh, and CSa) by a linear trend line revealing the overall direction of std and average deposition rates in blue and yellow respectively.

MERRA-II Deposition Rate
The MERRA-II model application data were downloaded and prepared according to the study area's coordinate system in order to determine the time series of the monthly deposition rate against seasonal variation. Figure 9 shows the monthly and annual average deposition rates. The maximum deposition rates occurred during the spring and summer, whereas the lowest values were observed in the fall and winter. In particular, the average dust deposition rate in BWh fluctuated between 4.86 t km −2 year −1 and 5.61 t km −2 year −1 , while in the steppe hot climate (BSh), the year-to-year deposition rate varied on average from 2.96 t km −2 year −1 to 3.78 t km −2 year −1 . The temperate hot summer zone (Csa) was characterized by an average deposition rate of less than 2.05 t km −2 year −1 , or less than 50% that of the hot desert zone. An important feature of the seasonal variation patterns for each climate zone is that their trends have the same seasonality, but at different deposition rate levels. The recorded trends of the GDR and MERRA II-output agree well with the seasonal dust distribution. Both datasets show the same summer maximum deposition rate factor, which is consistent with the climate pattern variation. Therefore, climate variability is closely related to deposition rate [10,11] and strongly governed by meteorological factors, e.g., [9,63]. Moreover, large-scale precipitation over land during cold seasons (fall and winter) resulted in lower deposition rates. The recorded trends of the GDR and MERRA II-output agree well with the seasonal dust distribution. Both datasets show the same summer maximum deposition rate factor, which is consistent with the climate pattern variation. Therefore, climate variability is closely related to deposition rate [10,11] and strongly governed by meteorological factors, e.g., [9,63]. Moreover, large-scale precipitation over land during cold seasons (fall and winter) resulted in lower deposition rates.

Mapping Seasonal GDR and Role of Climate
As shown in Figure 10, during the summer, the maximum percentages of the accumulated deposition rates were 35%, 28%, and 30% of the annual deposition rate in the BWh, BSh, and Csa, respectively. Further, in the spring, peak monthly dust deposition rates were 28%, 32%, and 31% in the BWh, BSh, and Csa regions, respectively.
This seasonal contribution rate is in accordance with the findings of previous studies conducted in eastern Iran, western Afghanistan [64,65], and Khuzestan [66]. Ref. [67] also obtained similar results for the Gansu Province and Gobi Desert (China), wherein the highest deposition rates were observed during the spring (32.51%) and the lowest values occurred in the fall. Although the prominent dust deposition rates occurred in the order of spring > summer, the BWh region results did not fit this rule. Rather, its dominant deposition rates occurred in the order of summer > spring for 2014 and 2017.

Mapping Seasonal GDR and Role of Climate
As shown in Figure 10, during the summer, the maximum percentages of the accu mulated deposition rates were 35%, 28%, and 30% of the annual deposition rate in th BWh, BSh, and Csa, respectively. Further, in the spring, peak monthly dust depositio rates were 28%, 32%, and 31% in the BWh, BSh, and Csa regions, respectively. This seasonal contribution rate is in accordance with the findings of previous studie conducted in eastern Iran, western Afghanistan [64,65], and Khuzestan [66]. Ref. [67] als obtained similar results for the Gansu Province and Gobi Desert (China), wherein th highest deposition rates were observed during the spring (32.51%) and the lowest value occurred in the fall. Although the prominent dust deposition rates occurred in the orde of spring > summer, the BWh region results did not fit this rule. Rather, its dominant dep osition rates occurred in the order of summer > spring for 2014 and 2017.
A study in Shuwaikh city [68] showed a similar pattern of monthly deposition rate in the summer (July and August 2009). Conversely, in western Taiwan [69] and north western Spain [70], observations showed unusually high winter deposition rates .

Typical Climate Terms Versus GDR
During the smpling period, the correlation coefficients obtained between the mete orological parameters, such as wind direction, precipitation, and relative temperatur and the GDRs were determined to be the main issues affecting the GDR, as shown in Fig  ures 11-14. Because climates are formed during long-term weather parameters, the influ ence of these meteorological parameters on the GDR during the study periodwas consid ered. Ultimately, the correlation values of 62%, −28%, and 36% were provide findings o [71], who showed the major effects of decreasing rainfall on dust deposition rates. A study in Shuwaikh city [68] showed a similar pattern of monthly deposition rates in the summer (July and August 2009). Conversely, in western Taiwan [69] and northwestern Spain [70], observations showed unusually high winter deposition rates.

Typical Climate Terms Versus GDR
During the smpling period, the correlation coefficients obtained between the meteorological parameters, such as wind direction, precipitation, and relative temperature, and the GDRs were determined to be the main issues affecting the GDR, as shown in Figures  11-14. Because climates are formed during long-term weather parameters, the influence of these meteorological parameters on the GDR during the study periodwas considered. Ultimately, the correlation values of 62%, −28%, and 36% were provide findings of [71], who showed the major effects of decreasing rainfall on dust deposition rates.
The correlation plot ( Figure 11) shows that positive and statistically significant correlations were detected between deposition rate and temperature. The discrimination of the relative variation in the GDR of eolian dust was further evident in the observed trends of the maximum and minimum temperatures and precipitation over the zonal surface. The correlation coefficient and relative Root Mean Square Error (RMSE) between the dust deposition and maximum temperatures were 62% and 9.48, respectively. Despite the variation in regional temperature, statistically significant positive trends were identified. Similar results were published by [6], who found that the dust deposition rate was considerably correlated (53%) with the average monthly air temperature. Furthermore, the results presented herein are in accordance with the results of [72], who found that a rise in the study area temperature region is associated with decreasing precipitation. This observation was expected, as the BWh and Csa regions have less evaporative cooling because of their lengthy sunshine duration. As illustrated in Figure 12, a statistically significant negative correlation exists between the monthly deposition rate and monthly precipitation (R = −0.28%, p < 0.01), which occurred because the dry deposition rate was minimized during the same period as precipitation [73]. This phenomenon may relate to the dry and wet deposition rates when returning to the surface via rain or gravity. The correlation plot ( Figure 11) shows that positive and statistically significant correlations were detected between deposition rate and temperature. The discrimination of the relative variation in the GDR of eolian dust was further evident in the observed trends of the maximum and minimum temperatures and precipitation over the zonal surface. The correlation coefficient and relative Root Mean Square Error (RMSE) between the dust deposition and maximum temperatures were 62% and 9.48, respectively. Despite the variation in regional temperature, statistically significant positive trends were identified. Similar results were published by [6], who found that the dust deposition rate was considerably correlated (53%) with the average monthly air temperature. Furthermore, the results presented herein are in accordance with the results of [72], who found that a rise in the study area temperature region is associated with decreasing precipitation. This observation was expected, as the BWh and Csa regions have less evaporative cooling because of their lengthy sunshine duration. As illustrated in Figure 12, a statistically significant negative correlation exists between the monthly deposition rate and monthly precipitation (R = −0.28%, p < 0.01), which occurred because the dry deposition rate was minimized during the same period as precipitation [73]. This phenomenon may relate to the dry and wet deposition rates when returning to the surface via rain or gravity. The correlation analysis suggests that a strong negative correlation exists between the deposition rate and precipitation over the climate zones, with a correlation coefficient of −0.28 and an RMSE of 10.58. In general, higher precipitation is correlated with lower deposition values, and vice versa. Thus, this relationship is logical as the BWh region is associated with less precipitation than the Csa and BSh regions. The mean annual precipitation varied by more than 420 mm in high-altitude (45.90° E 33.00° N) regions to less than 120 mm in arid regions (47.60° E 29.97° N). These results show that the lowest dust deposition rates occurred during the fall and winter, while the highest occurred in the  In addition to the horizontal components of the local vector wind velocities, the vertical wind and pressure velocities were measured and compared [82]. Moving from the earth's surface produced a vertical velocity with units of length/s that was near positive. In this scenario, the wind velocity variation significantly affects the dust transport process by lifting the particles [83,84] and degrading the air quality [22]. Moreover, the vertical velocity considerably affected the lifting force, which is primarily related to the differential velocities between particles and vertical flow [85].   In addition to the horizontal components of the local vector wind velocities, the vertical wind and pressure velocities were measured and compared [82]. Moving from the earth's surface produced a vertical velocity with units of length/s that was near positive. In this scenario, the wind velocity variation significantly affects the dust transport process by lifting the particles [83,84] and degrading the air quality [22]. Moreover, the vertical velocity considerably affected the lifting force, which is primarily related to the differential velocities between particles and vertical flow [85].  The correlation analysis suggests that a strong negative correlation exists between the deposition rate and precipitation over the climate zones, with a correlation coefficient of −0.28 and an RMSE of 10.58. In general, higher precipitation is correlated with lower deposition values, and vice versa. Thus, this relationship is logical as the BWh region is associated with less precipitation than the Csa and BSh regions. The mean annual precipitation varied by more than 420 mm in high-altitude (45.90 • E 33.00 • N) regions to less than 120 mm in arid regions (47.60 • E 29.97 • N). These results show that the lowest dust deposition rates occurred during the fall and winter, while the highest occurred in the spring and summer. This is consistent with the recent contribution of [74], who stated that precipitation could extract particles from the air. These findings are also consistent with previous studies by [73] and [75], which discuss dry deposition. Table 4 illustrates the average seasonal deposition rate in the three climate regions as related to the wind velocity. The results can be interpreted as follows: 1. BWh corresponds to eastward wind speeds (≤3.0 m s −1 ), with maximum deposition rates of 34.5 t km −2 and 6.2 t km −2 , as determined by the GDR and model output, respectively. 2. The BSh zone is related to westward wind speeds of ≤2.5 m s −1 , and has average deposition rates of 24.9 t km −2 and 3.3 t km −2 , as observed via GDR and the model output, respectively. 3. Finally, the Csa climate zone corresponds to eastward and westward wind speeds of ≤3.4 m s −1 , which are characterized by a lower dust portion than that of the BWh and BSh winds. Thus, the wind velocity correlated well with the observed seasonal deposition rates. Although wind speed and direction may play important roles in replenishing dust particles and positively influencing atmospheric concentration, settling dust particles by gravity do not require wind. A comparison between the model output and ground deposition rate (GDR) was conducted for a given time and the corresponding wind velocity near the surface. These data indicate an agreement between the model output and GDR ranges. A consistency index of wind speed (2.0 m s −1 < v < 3.5 m s −1 ) was observed in the westward wind direction, and a minimum peak percentage of calm events was found for each climate zone.
The soil entrainment thresholds for different soil texture classifications have been determined [76,77], and reviewed by [78]. However, [76] classified wind velocities to determine the terms in which the particles would be set in motion. Conversely, [78] evaluated the threshold wind velocity dynamics to identify the characteristics of threshold dynamics relative to sediment transportation. Thus, the threshold velocity can vary widely, both spatially and temporally, because the soil and aerodynamic properties of the surface vary [79].
Previously, dust intensity has been calculated using a combination of visibility and wind speed [80]. Ref. [67] suggested that the total suspended particles (TSP) concentration increased with increasing wind speed. Meanwhile, a study on Northeast Asian Dust found that wind velocities of 5 m s −1 to 10 m s −1 impacted the TSP (<500 µg m −3 ) and decreased visibility to 2000 m [81].
The results presented herein are consistent with those of [80], revealing that significant deposition occurs as a result of gravity when the wind speed is weak (Figure 13).
In addition to the horizontal components of the local vector wind velocities, the vertical wind and pressure velocities were measured and compared [82]. Moving from the earth's surface produced a vertical velocity with units of length/s that was near positive.
In this scenario, the wind velocity variation significantly affects the dust transport process by lifting the particles [83,84] and degrading the air quality [22]. Moreover, the vertical velocity considerably affected the lifting force, which is primarily related to the differential velocities between particles and vertical flow [85].
Thus, as expected, negative and statistically significant correlations were found between the deposition rate and vertical wind pressure velocity (−42%, p < 0.01). Lower values of vertical velocity pressure led to increased particle deposition (Figure 14). In accordance with the correlation plot, the dust-in-suspension and deposition rate events do not require any wind speed. In the BWh region, the vertical pressure velocity was relatively weak, reaching less than 0.01 Pa s −1 . Moreover, under the minimum vertical velocity, dust-in-suspension occurs, which is in accordance with previous research findings [86,87].
In this study, the typical weather patterns for each study area relative to the deposition rate were also examined, as summarized in Table 5, wherein the minimum, maximum, average, and standard values of temperature, precipitation, and GDR in the BWh, BSh, and Csa regions are listed for 2014-2017. Specifically, the correlation between weather patterns and the main climate zone parameter [88] indicated that high temperatures, precipitation, and wind direction are the most influential factors in the classification zone. In the BWh zone, the high dust deposition rates in the spring range from 16 t km −2 to 21 t km −2 with a mean value of 8.4 t km −2 month −1 . Meanwhile, the mean deposition rate of the BSh zone was 6.6 t km −2 , which is also considerably high. The lowest amount of dust deposition occurred in the Csa zone, at a value of 3.5 t km −2 .
Moreover, the climate pattern contributions and the classification of the three zones based on temperature and precipitation were evaluated. On average, the annual minimum precipitation and maximum temperature in the BWh region occurred at latitudes between 30 N and 32 N. Meanwhile, the same average temperature value (55 • C) was observed for the BSh zone, with a monthly precipitation value between 32 N and 33 N of 33.30 mm. In addition, the Csa climate boundary from 33 N to 34.5 N revealed a maximum monthly precipitation of 42 mm and an average temperature of 18.2 • C. These results could be considered to be positive consequences of high deposition rates, as higher quantities were deposited in arid and desert zones [89]. In accordance with [90] and [91], increasing air temperature led to increased evaporation and decreased precipitation, thus drying the surface. Higher temperatures also led to droughts in the subtropical dry zones of warmer climates, e.g., [72,92].

Significance of Relationship between MERRA-II and GDR
To evaluate the GDR results, a comparison between the MERRA-II and GDR over the research area during a period of high seasonal dust deposition rates was conducted. Although these physical air-to-surface processes are complex and depend on the dust concentration, some results show that the higher the dust concentration, the higher the dust deposition [87,93]. However, small changes in the statistical factors can cause changes in the deposition rates, as evidenced by the observed deposition rate trends found in the ground observations of the climate zonal surfaces. Thus, the eolian gradient between the MERRA-II and GDR results was considered ( Figure 15). In 2014, this gradient revealed a positive and statistically significant correlation with a relatively low RMSE (r = 67%, RMSE = 0.091) at the p < 0.01 probability level. In 2015, there was also a significant correlation (r = 65%), and a low RMSE (0.01) at the p < 0.01 probability level. Although a similar significance correlation (r = 66%, RMSE = 0.102) was obtained for 2016, different positive and statistically significant correlations were detected in 2017 from the MERRA-II and GDR data (r = 74%, RMSE = 0.08) at the p < 0.05 probability level.
Atmosphere 2021, 12, x FOR PEER REVIEW 20 of 25 positive and statistically significant correlation with a relatively low RMSE (r = 67%, RMSE = 0.091) at the p < 0.01 probability level. In 2015, there was also a significant correlation (r = 65%), and a low RMSE (0.01) at the p < 0.01 probability level. Although a similar significance correlation (r = 66%, RMSE = 0.102) was obtained for 2016, different positive and statistically significant correlations were detected in 2017 from the MERRA-II and GDR data (r = 74%, RMSE = 0.08) at the p < 0.05 probability level. A comparison of the deposition rates in the study area climate zones indicated that dust concentrations in the atmosphere represent the dust deposition rate factors [73,87]. In accordance with [75], the deposition rate is affected by climate factors. For example, the high eolian deposition rate is caused by the influencing factors of subtropical dry zones. The significance of the varying deposition rates among the climate zones can be evaluated in the context of climate regions. The results presented herein show that deposition rates vary significantly and are closely related to climate variability, which is consistent with A comparison of the deposition rates in the study area climate zones indicated that dust concentrations in the atmosphere represent the dust deposition rate factors [73,87]. In accordance with [75], the deposition rate is affected by climate factors. For example, the high eolian deposition rate is caused by the influencing factors of subtropical dry zones.
The significance of the varying deposition rates among the climate zones can be evaluated in the context of climate regions. The results presented herein show that deposition rates vary significantly and are closely related to climate variability, which is consistent with the findings of [7,23,75], confirming that the BWh region had the highest deposition rate at 8.4 t km 2 month −1 . Meanwhile, the BSh and Csa regions had lower deposition rates of 6.60 t km −2 month −1 and 3.30 t km −2 month −1 , respectively. In addition to these findings, the highest dust deposition rates followed the order of BWh > BSh > Csa during the research period.

Conclusions and Remarks
This study focused on the enhancement of the dust deposition rates of three climate zones, BWh, BSh, and Csa, in southwestern Iran. To better understand the role of eolian dust deposition rates relationship with the climatic region and its impact on air quality, significant efforts were made to increase both space-and ground-based observations. Herein, the GDRs and dust event frequencies, including the mass and weight compositions, were measured. Experiments on the links between dust deposition and climate patterns were quantitatively assessed for the GDR and MERRA-II model. By comparing these results, a high correlation coefficient was revealed for deposition rates during the research period. This strong correlation was calculated to be 0.67, 0.65, 0.66, 0.74 at a p < 0.01 significance level from 2014 to 2017, respectively. Then, the effect of the climate on the deposition rate is quantified, revealing the annual mean deposition rates to be 100.80 t km −2 year −1 , 79.27, and 39.60, for the BWh, BSh, and Csa zones, respectively. The seasonal variability of the deposition rates revealed that the highest rates occurred in the warm season and the lowest occurred in the cold season. The total annual dust deposition rates, however, were negatively related to the total annual precipitation, while they had a positive correlation with the annual temperature. Therefore, the atmospheric particle lifetime is controlled by precipitation and temperature, two of the main factors that characterize climate zones, which had a major influence on the dust deposition rates in the study area.
Thus, the results herein regarding the processes responsible for the dust deposition rates, and the climate pattern feedback provide insights concerning new sources, deposition rates, and their associated climate effects. Developing new sources are formed when a high volume of atmospheric particles is deposited by gravity over a land surface with no obstacles.
Although short-term changes in atmospheric conditions may influence the deposition rate (especially during simultaneous dust events), long-term atmospheric behavior corresponds to a large deposition rate. As eolian deposition rates are sensitive to climate zones, additional GDR observations on climate regimes should be performed to obtain precise information on dust plumes. Analyzing the seasonal and spatial variability of the dust activity in different climate zones is important not only for managing the impact of eolian dust but also for developing proper LULC strategies that can mitigate future dust impacts. This study might also contribute to the estimation of future dust deposition, as long as the GDR continues to be the main factor representing deposition trends. In the future, studies on the prediction of climate change and dust deposition rates should be combined.