1. Introduction
Mineral dust aerosols are one of the main factors affecting the Earth’s climate system, air quality, biogeochemical cycles (Martin, 1991; Archer et al., 2000; Kok et al., 2018; Middleton, 2019), atmospheric chemistry (Dentener et al., 1996; Dickerson et al., 1997; Martin et al., 2003), and ecosystem health (Prospero, 1999). Southwest Asia’s arid deserts are among the nine main dust sources worldwide (Prospero et al., 2002; Goudie et al., 2021) and contribute approximately 20% of the total dust emissions globally (Tanaka and Chiba, 2006). The dust emission rate from these deserts is approximately equal to the dust emission rate from the Bodele depression and central Saharan dust sources in Africa (Hamidi et al., 2017). The Rub Al Khali desert along the Saudi Arabia-Oman border, shown in Figure 1a, has been identified as the most prolific dust source on the Arabian Peninsula (AP) according to the Total Ozone Mapping Spectrometer (TOMS) (Goudie and Middleton, 2001; Giles, 2005), the MODIS instrument (Butt et al., 2017; Ali et al., 2017), and ground-based measurements (Notaro et al., 2013; Prakash et al., 2015).
Studies of sources of dust in Southwest Asia have classified these sources into two main groups: local sources and remote sources (Notaro et al., 2013). Local sources include the Rub Al Khali, the Ad Dahna, the An Nafud, and the Syrian and Iraqi deserts, while remote sources include the Sahara (Washington et al., 2003; Osetinsky and Alpert, 2006; Hannachiet al., 2011; Awad and Mashat, 2014; Yu et al., 2013; Beegum et al., 2018; Mashat et al., 2021) (Fig. 1a). Specifically, among local sources, the alluvial plain of the Tigris and Euphrates Rivers in Iraq is recognized as the main dust source in the Middle East region (Sissakian et al., 2013; Zoljoodi et al., 2013; Namdari et al., 2018; Mohammadpour et al., 2021).
Using backward trajectories, Notaro et al. (2013) identified the Rub Al Khali desert as a major local dust source for western Saudi Arabia and the Sahara Desert as a major remote source (Mashat et al., 2021). For northern and eastern Saudi Arabia, the Syrian and Iraqi Deserts were identified as the primary dust sources.
Dust storms are prominent over Saudi Arabia during February-June, with a peak along the southern coast of the Red Sea occurring in mid-winter, a peak across northern Saudi Arabia around the An Nafud Desert occurring in spring, and a peak over eastern Saudi Arabia around the Ad Dahna Desert occurring in early summer (Notaro et al., 2013).
The main atmospheric dust systems influencing the Middle East region are frontal dust systems with strong winds associated with intense baroclinicity. These systems form in non-summer seasons because of the presence of high pressure over the southeasternAP area and low pressure over the eastern Mediterranean region or due to the Siberian High over central Asia and low pressure over the eastern Mediterranean Sea (Middleton, 1986; Hermida et al., 2018; Francis et al., 2019; Hamidi, 2019; Mohammadpour et al., 2021, 2022). During February-April, the Mediterranean storm track is active, with passing cyclones and associated cold fronts carrying Saharan dust to Saudi Arabian stations along the northern coast of the Red Sea. Studies have shown that the most severe frontal dust storms occur during March and April (Hamidi, 2019) or they are generally prominent from February to June (Notaro et al., 2013; Yassin et al., 2018).
According to Abdi et al. (2011), dust storms are primarily triggered by dynamic lifting in the cool season in response to cold fronts and their associated mid-latitude troughs, or by diurnal vertical mixing in the warm season in response to solar heating. In turn, the radiative heating of dust aerosols can influence synoptic weather patterns (Islam and Almazroui et al., 2012; Liu et al., 2014; Zhao et al., 2015), e.g., increasing the Saudi Arabian heat low (Mohalfi et al., 1998; Francis et al., 2020).
Al-Jumaily and Ibrahim (2013) analyzed several dust storms in Iraq (northern AP) and found that their main driver is the passage of low-pressure systems over Iran, which carry cool air from the northern region toward a warmer region or warmer air over eastern Syria and Iraq (northern AP), following a temperature gradient (a front) formed or a releasing mechanism of dust formed along the Syrian Desert and over the eastern AP (Al-Dahna Desert) (Rao et al., 2001; Nabavi et al., 2017; Hamidi, 2019; Kaskaoutis et al., 2019; Hamzeh et al., 2021).
The concentration of atmospheric dust is tightly correlated with wind velocity (Kutiel and Furman, 2003; Kalenderski et al., 2013; Jish et al., 2015; Gandham et al., 2020). Strong surface cyclones can also stir up dust clouds or lift dust into the atmosphere along convergence zones, in monsoon regions, between cold air masses, and in association with cyclones and tropical anticyclonic air masses.
Atmospheric analysis of the synoptic patterns associated with or even favoring dust outbreaks over the AP is important for issuing forecasts and warning the public in cases of severe storms, and it could provide useful information about the atmospheric characteristics of dust events (Houssos et al., 2015; Hermida et al., 2018; Mashat et al., 2018).
Previous studies of dust systems influencing the AP (e.g., Goudi and Middleton, 2006; Hamidi et al., 2013; El Kenawy et al., 2014; Mashat et al., 2018; Rashki et al., 2019) did not differentiate between the atmospheric dust systems influencing different climate regions on the AP, but indicated that a combination of more than one of the atmospheric systems could influence dust systems. Therefore, the purpose of this study was to detect the specific atmospheric regimes influencing the individual regions of the AP (northern, western, and eastern regions) and to determine the spatial distribution of atmospheric dust systems.
This paper is organized as follows. In section 2 the data and methodology are described. Section 3 presents the statistical results of the surface measurements and the synoptic features of the different regions. The final two sections present the discussion and conclusions of the study.
2. Data and methodology
In this study, two sets of data were used: ground-based dust measurements from 27 surface stations in Saudi Arabia and 6-h meteorological reanalysis data from the National Centers for Environmental Prediction (NCEP/NCAR) from 1978 to 2008, with a horizontal spatial resolution of 2.5º × 2.5º (Kalnay et al., 1996; Kistler et al., 2001).
Daily surface observations of dust from 1978 to 2008 were retrieved from meteorological stations in Saudi Arabia (Fig. 1b). The surface station observations of dust and sand have been classified into two categories: the dust type, including suspended dust, and blowing dust, which have World Meteorological Organization (WMO) report codes (WW codes) 06 and 07, respectively; and the dust storm type, including dust storm observations with WW codes 09 and 30-32, and severe dust storm observations with WW codes 33-35 (WMO, 2005).
According to the WMO classification (WMO, 2005), Rashki et al. (2017), and Mashat et al. (2016), the characteristics of these types are as follows: (1) suspended dust (WW 06), widespread dust suspended in the atmosphere and not raised at or near the stations; (2) blowing dust (WW 07), dust or sand raised at the time of observation; (3) dust storms (WW 09 and 30-32), large dust or sand quantities raised by strong winds with visibility < 1000 m but > 200 m; and (4) severe dust storms (WW 33-35), very strong wind lifting a large number of soil particles, with visibility < 200 m.
The synoptic evaluation of the composite maps for each region was performed based on the sea level pressure (SLP); the wind components at 850, 500, and 250 hPa; the vertical velocity at 850 hPa; the geopotential heights at 850 and 500 hPa; and the static stability in the 1000- to 500-hPa layer. The formula of the static stability used depends on Eq. (1):
−(/θ) × (/∂) (1)
where T is temperature and θ the potential temperature at pressure (P) levels of 1000 and 500 hPa (Gates, 1961).
To describe the spatial characteristics of the dust system, the AP can be divided into three dust regions based on previous dust studies (Middleton, 1986; Notaro et al., 2015; Awad and Mashat, 2016; Mashat et al., 2017, 2018, 2019) and not into five subregions as in Alobaidi et al. (2017) based on wet season climatological features. The northern region of dust (NRD) contains 10 stations (red stations in Fig. 1b), the western region of dust (WRD) contains 12 stations (blue stations in Fig. 1b), and the eastern region of dust (ERD) contains five stations (black stations in Fig. 1b), a list of detailed information for the stations used is in Table SI in the supplementary material.
This study statistically assessed the two observed dust types from two perspectives: spatially and temporally. In the first analysis, the spatial distribution of the two types at each station was represented by the frequency of each type. The size of the observation sample (circles for dust and triangles for dust storms; Fig. 1b) is proportional to the ratio of occurrences with respect to the relative observations at the other stations, i.e., the largest symbol size corresponds to the highest frequency, and the smallest symbol size corresponds to the lowest frequency. The temporal distributions, at annual and monthly scales, were calculated as the total number of observations over the entire area, i.e., as the number of observations per day from all stations.
In the synoptic study, to select the dust storm events associated with typical atmospheric systems in the entire AP and for each climate region, the following conditions were strictly followed:
Widespread dust events affected most of the AP, rather than only one or two climate regions, where these dust events affected at least one station in each of the three climate regions.
In each climate region, a dust event was defined in cases with at least one observation of a dust storm at the stations within the specific climate region but not observed at any station in other climate regions. For example, the cases in the northern region were those observed only at stations in the northern region and not observed at any of the stations belonging to the western or eastern regions.
For the case study, the second condition was used, but at least three stations had to have observed the dust storm.
This study defined the synoptic characteristics of AP dust events and the regional dust storm characteristics via the following steps:
Comparing the characteristics of composite events for different regions to determine how much the atmospheric dust characteristics differ among the different regions and to determine the specific characteristics that distinguish each region from the others.
Determining the common atmospheric characteristics associated with dust storm events influencing the AP.
3. Results
The spatial distribution of ground-based dust events (Fig. 1b) shows that most of the dust is concentrated in the northern and eastern regions at the Tabuk, Riyadh and Dhahran stations, in agreement with the surface observation results from Notaro et al. (2013), the aerosol optical depth (AOD) results from Ali et al. (2017), and the results from TOMS AI (Prospero et al., 2002; Washington et al., 2003). A moderate number of dust observations are present at the northeastern stations, i.e., Turaif, Arar, Hail, and Hafrbaten, and at the western stations, especially Wadi-Aldawaser and Najran, which are close to the desert (Mashat et al., 2018). Additionally, the stations on the Red Sea or close to the Red Sea, such as Alwejh, Makkah, Yenbo, and Abha, have few dust observations (Notaro et al., 2013; Alobaidi et al., 2017).
Additionally, the spatial distribution of dust storm observations (Fig. 1b) shows that most of these are concentrated on the eastern boundary of the AP at the Turaif, Arar, Aljouf, Rafha, Alqaisumah, and Alhasa stations, as well as some other stations such as Tabuk, Hafrbaten, Riyadh(s), Gassim and Wadi-Aldawaser. In contrast, most stations close to the Red Sea, such as Alwejh, Makkah, Albaha, Almadinah, Jeddah, Altaif, Khamis-Mushait and Abha, have few observations of dust storms. Generally, the distribution of dust storms appeared similar to the annual distribution presented by Hermida et al. (2018) or the distribution of AOD (Ali et al., 2017).
The annual distribution of dust in Fig. 2a shows two distinct periods with different frequencies: from 1978 to 1994 and from 1995 to 2008. The first period has an annual average value of 1011.88 observations, corresponding to 37.48 observations per station. The annual average of observations in the second period is 1066, corresponding to 39.48 observations per station.
The annual distribution of dust storms in Fig. 2b shows three distinct periods. The first one, from 1978 to 1984, has an annual average of 52.43 observations or approximately two observations per station. The second period, from 1985 to 1994, has an annual average of 123.56 observations or approximately 4.6 observations per station. The third period, from 1995 to 2008, has an annual average of 31.14 observations or approximately 1.2 observations per station.
While the above results indicate that the period from 1985 to 1994 had a large number of dust phenomena, the period from 1978 to 1984 had the smallest number of dust phenomena observations. Furthermore, although the period from 1995 to 2005 had a large number of dust events, it had the smallest number of dust storm observations.
The monthly distribution (Table I) of dust shows that the average is approximately 326.4 350.1, and 359.8 observations/season in March, April, and May, respectively, indicating that the monthly average increases from cold months (March) to warm months (May), in agreement with the findings of Notaro et al. (2013) and Rezazadeh et al. (2013).
Generally, the annual distribution of spring months (Fig. 3) shows a wavy shape with a peak in the early 1990s and lows in the early 1980s and late 1990s.
In detail, although the monthly pattern is similar to the annual distribution of dust (Fig. 2a), the distribution has high variability in May (Fig. 3e) and low variability in March (Fig. 3a). In addition, the distribution indicated two trends: a rapidly increasing rate from 1978 to 1994 and a slowly decreasing rate from 1995 to 2004.
In March, from the annual distribution of dust storms (Fig. 3b), two periods of different trend rates were distinguished: an increasing rate period from 1978 to 1994, with an annual average of 30.18 cases or 1.1 cases for each station/season; and a low decreasing rate period from 1995 to 2008, with an annual average of 8.1 cases or 0.3 cases for each station/season.
In April (Fig. 3d), two distinct periods were distinguished: the first from 1978 to 2000, with an annual average of 26.29 cases or 0.97 cases for each station; and the second from 2001 to 2008, with an annual average of 8 cases or 0.3 cases for each station.
In May (Fig. 3f), except for the four years from 1991 to 1994, which had a large number of cases with an annual average of 71.25 or 2.64 for each station, two periods were distinguished. The first one, from 1978 to 1990, had a low increasing rate, with an annual average of 23.46 cases or 0.87 cases for each station. The second period, from 1995 to 2008, had a very low decreasing rate, with an annual average of 9.57 cases or 0.33 cases for each station.
From the previous discussion, the monthly distribution appears similar to the first three clusters from Hermida et al. (2018), while the maximum number of dust cases was observed in May, not in March, as in Notaro et al. (2013).
3.1 Climate synoptic characteristics of widespread cases
In this subsection, the widespread cases (WSC), i.e., those cases influencing all three climate regions of the AP, were selected to represent the general dust storm events that affected the AP.
3.1.1 Mean sea level pressure and maximum wind at 250 hPa
Generally, widespread dust cases on the AP were influenced by two low-pressure centers (Fig. 4a): the northern center was located over Iraq and northeastern Saudi Arabia, with a lowest pressure value of 1006 hPa; and the southern center was located over the southwestern AP, with a lowest pressure value of 1005 hPa. Additionally, the previous low-pressure system was surrounded by either high-pressures or ridges, forming high-pressure gradient areas surrounding the AP.
The Azores high, with a high-pressure center of 1016 hPa over the southern Mediterranean and Libya locally known as Libyan high pressure (Washington and Todd 2005), formed a western pressure gradient area over the northwestern AP. The Arabian Sea high-pressure center, with a pressure value of1010 hPa, interacted with the low-pressure system and formed the southern AP pressure gradient area. The ridge of the Siberian high, with a pressure value of 1008 hPa, over Iran and the Arabian Gulf (AG) (or Persian Gulf) interacted with the low-pressure system and formed the eastern pressure gradient area.
The maximum wind at 250 hPa (Fig. 4a), with a maximum value of 82 knots, has a trough shape over the northern Red Sea and an exit region over Iran. The maximum wind was observed to be a factor associated with dust storms over the region (Francis et al., 2019).
3.1.2 Geopotential height and vertical velocity at 850 hPa
At a pressure level of 850 hPa (Fig. 4b), the northern cyclone located over Turkey, with a value of 1450 gpm, exhibited a deep trough over the AP, and in conjunction with the subtropical anticyclone over Africa, with a geopotential height of 1530 gpm, it formed an area of geopotential gradient over the western AP. The trough of the northern cyclone, in conjunction with the anticyclonic center with 1500 gpm over the Arabian Sea and southeastern AP, formed an area of geopotential gradient over the eastern AP and AG. The wind speed over the AP at this pressure level increased from south to north until reaching the northern AP, and then it decreased.
The distribution of omega at a pressure level of850 hPa (Fig. 4b) shows that the vertical velocity followed the surface pressure systems, where upward motion accompanied low pressure, and downward motion accompanied the high-pressure centers or ridges. The highest upward motions (-1.0 and -1.2 Pa s-1)were located over the southwestern AP and Iraq/Iran, respectively, while the highest downward motion (+1.0 Pa s-1) was located over the Arabian Sea.
3.1.3 Geopotential height at 500 hPa and the static stability layer from 1000 to 500 hPa
At a pressure level of 500 hPa (Fig. 4c), the trough ofthe northern cyclone, with a geopotential valueof 5560 gpm, extended southward over the Red Sea and eastern Africa, while the subtropical anticyclone over Africa, with a geopotential value of 5880 gpm, formed a ridge over the eastern AP and AG. The wind speed at this level increased over the AP from south to north, ranging from 5 to 45 knots, and decreased, between 5 and 10 knots, from west to east.
The spatial distribution of the static stability (Fig. 4c) shows a low-stability area located over the AP, with a value of 0.4 deg db-1 and surrounded by high-stability areas. The high-stability areas surrounding the AP had values of 0.85, 1.1 and 0.85 deg db-1over the Mediterranean Sea, the Arabian Sea and the northern region, respectively. Furthermore, two high-stability centers of 0.65 and 0.55 deg db-1 were formed over the Red Sea and AG, respectively. The interaction between the low stability area overthe AP and the surrounding high-stability areas formed areas of a stability gradient around the AP, with a weak gradient area over the eastern AP.
3.2 Climate synoptic characteristics of the classified regions
In the following sub-section it can be seen that for each climate region the composite of the dust event satisfied the identifying conditions for the climate regions.
3.2.1 Mean sea level pressure and maximum wind at 250 hPa
Analysis of the atmospheric systems associated with the dust cases influencing the NRD (Fig. 5a) shows that this region was influenced by a low-pressure system, with a value of 1008 hPa and a deep trough over the AP, i.e., 2 hPa lower and located north of the corresponding system in WSC. The Azores high is associated with a high-pressure center with a value of 1017 hPa over northern Libya (Libyan high pressure), i.e., 1 hPa higher than the corresponding system in WSC, and a ridge over the eastern Mediterranean region. The Siberian high, with a value of 1019 hPa, is associated with a ridge over the AG. The Sudan low, with a value of 1005 hPa, is associated with the RST over the western Red Sea and extends northward to southern Egypt, i.e., extending northward farther than the WSC. The southern AP low pressure, with a value of 1007 hPa (less deep than the WSC), has two troughs, a weak one over the western AP and a strong one over the eastern AP, and is linked with the northern low-pressure system with a value of 1008 hPa, i.e., 2 hPa higher than the corresponding system in the WSC. This atmospheric regime is known as the frontal dust system (Hamidi et al., 2013).
The southern low-pressure system forms a strong pressure gradient area over the eastern Mediterranean region with the Azores ridge, and the north-south low-pressure system forms a pressure gradient with the Siberian ridge over Iraq and western Iran. These pressure gradients have been considered a characteristic of dust events because they increase the wind speed over a large area (Zoljoodi et al., 2013; Al-Yahyai and Charabi, 2014; Hermida et al., 2018).
The maximum wind at 250 hPa (Fig. 5a) shows that the core with the highest wind speed is located over the northern AP, with a value of 77 knots, i.e., 5 knots less than that of the WSC, and it forms a trough-shaped structure over the eastern Mediterranean region.
Analysis of the atmospheric systems influencing the WRD (Fig. 5b) shows that the Azores high is associated with a high-pressure center, with a value of 1016 hPa over Libya and a ridge over the eastern Mediterranean region. The Siberian high, with a value of 1019 hPa, has a pronounced ridge over the AG; i.e., its ridgeline is located eastward of that of the WSC. The Sudan low, with a value of 1005 hPa, is associated with a pronounced RST over the western Red Sea and extends northward to southern Egypt, i.e., extending farther northward than that of the WSC. The southern AP low pressure, with a value of 1007 hPa, has a pronounced trough over the western AP, i.e., deeper than that of the WSC, and it forms a gradient area with the Azores ridge over the eastern Red Sea, northern AP and Levant. Both the Siberian and Azores ridges form a wave-shaped pressure gradient area with a southern AP low-pressure system over the northern AP.
The core of the maximum wind of 250 hPa (as shown in Fig. 5b), is located over the eastern AP and northern AG and appears south of the composite results for the NRD. Additionally, an elongated core with a value of 72 knots, i.e., 5 knots slower than that of the WSC, extends over the Saharan region to central Egypt. The maximum wind forms a pronounced trough-shaped structure over the eastern Mediterranean and Red seas.
In the analysis of dust in the ERD (Fig. 5c), the Azores high is associated with a high-pressure center over southern Italy, with a value of 1016 hPa, and a pronounced ridge over the eastern Mediterranean and Red seas. The Siberian high, with a value of 1018 hPa, is associated with a ridge over the AG. The Sudan low, with a value of 1005 hPa, is associated with an RST over the western Red Sea, and the southern AP low pressure, with a value of 1006 hPa (i.e., less deep than the WSC by 1 hPa), has two cores over the southwestern and southeastern AP (rather than only one core, as in the WSC), representing one of the main dust systems that influence the region (Hamidi et al., 2014). Additionally, the southern low pressure has two troughs over the eastern and western AP and, with the Siberian and Azores systems, forms a pronounced wave-shaped pressure gradient over the AP, which, as noted by Hamidi et al. (2013), might transport dust from the Tigris and Euphrates area.
The maximum wind at 250 hPa (Fig. 5c) has a core over the northern AP with a value of 72 knots, and it is located to the south of its locations in the previously discussed situations. This southward shift of the maximum wind was explained by Perrone (1981) and Francis et al. (2019) as the southward progression of the polar jet in the AP, which enhanced the shamal wind. The maximum wind forms a trough shape over the eastern AP. The atmospheric systems associated with the dust events of this region are similar to clusters 1 and 2 (Hermida et al., 2018), and the relevant results suggest that the remote source of dust in these systems is the Syrian and Iraqi deserts, while the main source is the Rub’al Khali Desert.
3.2.2 Geopotential height and vertical velocity at 850 hPa
Analysis of the atmospheric systems affecting the dust conditions in the NRD (Fig. 6a) shows that the northern cyclone, defined by the 1460 gpm contour, has a pronounced trough over the northern AP, in agreement with the April cases (Beegum et al., 2018). Obviously, the NRD northern cyclone is less deep; its trough weakens southward and appears to the west of the corresponding trough in the WSC. The subtropical anticyclone has two branches. One branch is located over Africa and has a value of1530 gpm and a ridge extending eastward to the western AP. The other branch is located over the AP and has a value of 1510 gpm and a ridge extending over the northern AP to Iraq, i.e., farther north than that of the WSC. The cyclonic and anticyclonic systems form two areas with geopotential gradients: one over the eastern Mediterranean and one over the northern AP and AG. The distribution of wind values shows that the wind decreases from west to east and from south to north.
The vertical velocity (Fig. 6a) forms a series of areas with upward and downward motions around the AP, including downward motion over the Red Sea, upward motion over the western AP and Levant, and downward motion over the AG. The highest values of downward motion are 0.4 and 1.0 Pa s-1 over the Red Sea and AG, respectively. The highest values of upward motion are -1.0 and -1.2 Pa s-1 over the southern AP and Levant, respectively.
The composite of the dust cases in the WRD (Fig. 6b)shows that the African branch of the subtropical anticyclone, with a value of 1530 gpm, shrank westward to eastern Africa and that the AP branch, with a value of 1510 gpm, shrank toward the eastern AP. The northern cyclone trough, with a value of 1490 gpm, flattened over the northern AP, while the tropical cyclone trough extended northward to the central AP. Obviously, the WRD northern cyclone weakened with depth, while the tropical cyclone deepened in comparison with the corresponding system in the WSC. The wind in this situation increased from west to east but was weaker than that in the NRD. The high wind in the NRD could be explained by frontal systems (Hamidi, 2019)
Although the general features of vertical velocity (Fig. 6b) are similar to those for the dust composite of the NRD, the magnitude of the vertical velocity changed to -1.2 and -0.8 Pa s-1 over the AP and Levant, respectively. The downward motion became 0.6 and 1.0 Pa s-1 over the Red Sea and AG, respectively. These changes show that the vertical velocity in the dust composite for the WRD intensified over the western regions, the Red Sea, and the western AP and AG, while it decreased over the Levant.
The composite of the dust cases influencing the ERD (Fig. 6c) shows that a subtropical anticyclone over Africa, with a value of 1530 gpm, extends eastward to the western AP, i.e., it is 30 gpm higher than the corresponding system in the WSC. A subtropical anticyclonic center, with a value of 1510 gpm, covers the eastern part of the AP. A tropical cyclone, with a value of 1490 gpm, is located over the tropical region southwest of the Red Sea. A northern cyclone, with a value of 1470 gpm, has a weak trough over the Levant, i.e., less deep than that in the WSC. In addition, a pronounced cyclonic center over Iran, with a value of 1470 gpm, forms an area of geopotential gradient over the AG and western Iran with the Arabian anticyclonic center. The wind speed in this composite decreases from west to east.
In this composite (Fig. 6c), the features of the vertical velocity are similar to those in the other composites, but the magnitude and area of vertical velocity are different. The downward motion over the Red Sea, with a value of 0.6 Pa s-1, extends slightly eastward over the AP. The downward motion over the AG, with a value of 0.8 Pa s-1, is smaller. The upward motion over the AP has a value of 1.0 Pa s-1 and is slightly larger in response to similar mechanisms discussed previously.
3.2.3 Geopotential height at 500 hPa and the static stability layer from 1000 to 500 hPa
The northern cyclone in the NRD composite (details in Fig. S1a in the supplementary material), has a trough over the eastern Mediterranean and Egypt, i.e., to the west of that in the WSC. Additionally, a belt of low stability is observed over the Sahara, the AP and Iran, surrounded to the north and south by belts of high stability.
In the WRD (Fig. S1b) the trough of the northern cyclone is oriented more eastward and appears relatively flatter than those in the NRD and WSC. Additionally, the areas of low stability over the northern and southern AP shrank and expanded, respectively.
Based on the composite of dust events in the ERD (Fig. S1c), the trough of the northern cyclone shifted eastward compared with that in the above composites. Additionally, a new weak trough appeared over the eastern AP. The stability distribution of the ERD is similar to that in the composite for the other regions but increased over the Red Sea and decreased over the AG.
3.3 Case study
For each case study of the climate regions, the composite of the dust events was defined for dust cases for which at least three stations observed dust storm events within the specific climate region, but there were no dust storm events observed at any station in other climate regions.
3.3.1 Mean sea level pressure and maximum wind at 250 hPa
The synoptic features of the NRD cases (Fig. 7a) show that the northern cyclone over the Levant, with a value of 1006 hPa, is associated with a northwest-southeast trough over the northeastern AP. The southern cyclone, with a value of 1007 hPa, is associated with a weak trough over the southwestern AP. The Sudan low, with a value of 1005 hPa, is associated with an RST extending northward into southern Egypt. The Azores high is associated with a high-pressure center over northern Libya, which has a value of 1018 hPa and a ridge over the eastern Mediterranean region. The Siberian high, with a value of 1019 hPa, is associated with a ridge over western Iran. In conjunction with the northern cyclone over the Levant, these high-pressure or ridge systems form two areas with high-pressure gradients over the eastern Mediterranean and over Iraq. The maximum wind (Fig. 7a), with a value of 84 knots, i.e., greater than WSC by 2 knots, forms a trough shape over the northern AP.
In the WRD cases (Fig. 7b), the high-pressure center over Libya, with a value of 1016 hPa, extends over northern Libya and Egypt, and a southwest-northeast oriented ridge extends over the eastern Mediterranean. Additionally, the high pressure center forms a pronounced ridge over the Red Sea, reaching south to approximately 18º N. The Siberian high is associated with a high pressure center, with a value of 1011 hPa, over the northern AG region. The southern cyclone, with a value of 1007 hPa, has a pronounced trough over the western AP. The Sudan low, with a value of 1005 hPa, is associated with an RST extending northward to northern Sudan. With the trough of southern low pressure, the high-pressure ridge over the Red Sea and the high-pressure center over northern AG form two areas with pressure gradients over the western AP and the central AP. The maximum wind (Fig. 7b) has a core with a value of 77 knots over the northern AP, i.e., 5 knots less than the WSC, and it forms a weak trough shape over the northern Red Sea.
The atmospheric systems that characterize the ERD cases (Fig. 7c) show that the high-pressure center over the eastern Mediterranean, with a value of 1015 hPa, is associated with a weak pressure gradient over the eastern Mediterranean and a pronounced ridge over the Red Sea, extending southward toapproximately 17º N. The Siberian high, with a value of 1017 hPa, has a weak, flat ridge over the eastern AP. In contrast, the Sudan low, with a value of1007 hPa, has a wide trough over Sudan. The southern low pressure, with a value of 1005 hPa, has a core over the southeastern AP and forms two troughs over the AP. A weak trough is located over the eastern region, and a relatively strong trough is located over the western region, although the eastern trough has a relatively high-pressure gradient area extending north-south over the eastern AP. The maximum wind (Fig. 7c) has a core with a value of 85 knots that extends over northern Africa and the northern AP, i.e., greater than WSC by 2 knots, and it forms a trough shape over the AP.
3.3.2 Geopotential height and vertical velocity at 850 hPa
In the NRD cases (Fig. 8a), the northern cyclone, with a value of 1440 gpm, has a pronounced trough over the northern AP. The tropical cyclone, with a value of 1490 gpm, has a trough over Sudan. The subtropical anticyclone over Africa, with a value of 1530 gpm, has a ridge extending eastward to the eastern Red Sea. The subtropical anticyclone center over the eastern AP, with a value of 1510 gpm, appears to relate to the African anticyclone and forms an anticyclonic, bell-shaped structure over the southern region. The above systems interacted with each other, and the northern cyclone trough formed a distinct shape over the eastern Mediterranean, Levant and northern AP. The wind over the AP decreased from west to east and from north to south. Additionally, a pronounced anticyclone wind pattern formed over the eastern AP.
The vertical velocity in the NRD cases (Fig. 8a) shows a series of downward and upward motions from west to east of the AP. The downward motion over eastern Africa and the Red Sea are extensions of downward-motion cores over the Mediterranean and Arabian Sea, with values of 0.6 and 1.0 Pa s-1, respectively. Most of the AP, from the west, influenced the upward motion with a value of -1.0 Pa s-1, while the northern AP and Levant have a high upward motion of -1.2 Pa s-1. In contrast, the AG has a high downward motion of 1.0 Pa s-1. This distribution of vertical velocity forms a high-gradient vertical velocity area over the AP, with a relatively large scale over the Levant and northern AP.
In the WRD cases (Fig. 8b), the African anticyclone shrank westward, and the AP anticyclone shrank until it was over the eastern AP. The northern cyclone trough shrank northward and became flat. The tropical cyclone extended northward over the entire western AP. The wind over the AP decreased from south to north over the western AP, and a pronounced anticyclone wind pattern formed over the eastern AP. Also, the core of high vertical velocity shifted southward over the central Red Sea, where it exhibited a downward-motion value of0.8 Pa s-1 and undertook an upward motion over the western AP, where it exhibited a value of -1.2 Pa s-1.The downward motion over the AG expanded. This distribution formed a high-gradient area of vertical velocity over the AP, especially over the western AP.
In the ERD cases (Fig. 8c), the African anticyclone, with a value of 1530 gpm, extended eastward over the western AP. A new core of the northern cyclone, with a value of 1440 gpm, located over the Caspian Sea formed two troughs: a weak trough over the northern AP and a strong trough over Iran. The tropical cyclone, with a value of 1490 gpm, extended northward over the central AP to the east of the WRD. These systems formed an area of geopotential gradient over AG and Iran and became stronger northwestward over the eastern AP. It is also observed that the core of the maximum upward motion shifted southward, with a value of -1.1 Pa s-1 over the AP, and a high-value area, with a value of -1.1 Pa s-1, was also observed over Turkey. In addition, the values of downward motion over the Red Sea and AG were 1.0 Pa s-1 and 0.8 Pa s-1, respectively. Therefore, the downward motion increased over the Red Sea and decreased over the AG.
3.3.3 Geopotential height at 500 hPa and the static stability layer at 1000-500 hPa
The northern cyclone associated with the NRD cases (see details in suppl Fig. S2a) has a deep trough over the eastern Mediterranean and eastern Africa. The subtropical anticyclone has pronounced ridges over western Africa and the AP. Additionally, the stability distribution shows a low stability over the AP, Iran and Pakistan. The low-stability areas over the AP interacted with the surrounding high-stability areas forms areas with stability gradients around the AP.
For the WRD cases (Fig. S2b), in comparison with that of the NRD cases, the trough of the northern cyclone is shifted eastward, while the subtropical anticyclone is shifted eastward and is relatively flatter. Additionally, the African stability decreased, and the stability generated two new centers of low stability over the northern and southern AP.
For the ERD cases (Fig. S2c), in comparison with that of the WRD cases, the trough of the northern cyclone is shifted eastward. The Arabian ridge is flat, while the African ridge becomes pronounced. Additionally, the low-stability area over the AP is shifted to the southern region and the high-stability areas extend over the Mediterranean and the Red Sea.
4. Discussion
The results of this study were based on a statistical study of surface observations of dust and dust storm types for selected cases affecting individual regions, which influenced more than two regions, and on meteorological data from the NCEP/NCAR reanalysis data set to determine the synoptic characteristics of the regional composites and widespread composite.
The spatial distributions of dust and dust storms demonstrated that most of the area north of 24º N is strongly affected by dust phenomena, in agreement with the findings of Kutiel and Furman (2003). Additionally, the Dawasir station, which is close to the desert, is also strongly affected (Notaro et al., 2013).
The spring monthly average of dust (dust storm) events indicated that the dust (dust storm) event frequency increased with warming from month to month, in agreement with the findings of Notaro (2015). The annual distribution of dust (dust storms) showed a wavy shape, with a high, increasing rate in the first study period until 1994 and a low, increasing (decreasing) rate in the second study period until 2004.
Although the monthly annual distribution of dust (dust storms) stations could be divided into two periods with different rates, the distribution in March exhibited pronounced changes, while the distributions in April and May exhibited small changes.
The synoptic climatology of the different regions indicated that the pronounced atmospheric systems in the NRD are the northern low-pressure system and its connection with the southern low-pressure system over the AP. In addition, the different atmospheric systems that interacted and produced areas of pressure gradient around the northern AP were considered a factor in the transport of dust from Africa by Awad and Mashat (2014) or a factor in the transport of dust from the Syrian and Iraqi Deserts to the northern AP by Hamidi et al. (2013) and Notaro et al. (2013). Moreover, the pressure gradient lines were oriented more eastward than those in the WSC, and the distributions of atmospheric systems are similar to those of frontal dust systems (Hamidi et al., 2013; Prakash et al., 2015; Yu et al., 2015; Almazroui and Awad, 2016; Nabavi et al., 2016; Beegum et al., 2018; Hermida et al., 2018), Mediterranean storm tracks (Notaro etal., 2013;), and Saharan cyclone tracks (Hannachi et al.,2011); therefore, they are called frontal systems.
For the WRD composite, the most pronounced atmospheric system was the RST. In addition, the southern low pressure formed a pronounced trough over the western AP, with a wavy pressure gradient over the northern AP, in agreement with the results of Saaroni et al. (1998), Mashat and Awad (2015), Prakash et al. (2015), and Mashat et al. (2018). Other notable features included the pronounced Siberian ridge over the AG and the maximum wind core located over the eastern AP in a position farther north than that in the NRD composite. Therefore, this system is called the RST system.
The atmospheric dust systems in the ERD composite are characterized by a southern low-pressure system that forms wavy isobars and a pronounced trough over the eastern AP. The trough and core of the maximum wind are located over the eastern and southern AP, respectively, unlike in the other situations. Additionally, the southern low-pressure system deepens to 1 hPa, and the Azores high-pressure center is located farther north than in the WRD situation. This appears similar to the third shamal system described by Hamidi et al. (2013); additionally, it could be considered a thermal low system (Mohalfi et al., 1998; Johnson, 2003; Das et al., 2015; Jin and Sun, 2016) and therefore it is called the thermal low system.
At a pressure level of 850 hPa, most of the AP in the NRD, or the frontal system, was affected by the subtropical anticyclone, and the northern AP was affected by the trough of the northern cyclone. In the WRD composite or the RST system, the subtropical anticyclone and northern cyclone retreated from the AP, while the tropical cyclone trough extended northward over the western AP. However, in the ERD composite or the thermal low system, the African branch of the subtropical anticyclone extended eastward to the western AP. The northern cyclone had a trough over the northern AP, and the cyclonic center over Iran formed a geopotential gradient area with the Arabian branch of the subtropical anticyclone over the AP and AG. This situation was explained by Hamidi et al. (2014) to be a factor responsible for emitting dust from local deserts or transporting dust from remote deserts. The interaction between the atmospheric systems oriented the geopotential gradient over the AG more in a north-south direction than in the other composites.
Although the features of the vertical velocity at 850 hPa are generally similar between the different dust composites, the highest values and their positions differ between the different composites. The centers of the core of the strongest vertical velocity (downward and upward) varied between the different dust situations. The cores were in the northern, central, and southern AP in the frontal system, the RST system, and the thermal low system, respectively. The distribution of the downward and upward motion forms two vertical cells, with the upward branch and two downward branches in the surrounding areas. Generally, the axis of the two vertical circles, with the axis passing between the downward-motion areas, rotated anticlockwise as the dust region shifted from the frontal system to the RST system and then to the thermal low system.
The atmospheric systems at 500 hPa indicated that the Middle East was influenced by the northern cyclone trough in the dust cases, but this trough shifted eastward and became less pronounced as the dust region shifted from the frontal system to the RST system and then to the thermal low system. In addition, the Arabian subtropical ridge became flat and shifted eastward as the dust region shifted anticlockwise from the frontal system to the RST system and then to the thermal low system. The wind speed at the 500-hPa level reached its highest value when the dust was affected by the thermal low system and reached its lowest value when the dust was affected by the RST system.
The relationship between the stability over the Mediterranean, Red Sea and AG plays an important role in the influence of dust storms on the area. This study shows that, when the stability over the Mediterranean decreased, dust storms influenced the northern region, but when the stability over the Mediterranean and AG increased, dust storms influenced the western region. Furthermore, dust storms influenced the eastern region when the stability increased over the Mediterranean and Red Sea regions. Additionally, the results indicated that the formation of dust events requires not only low stability (Alharbi et al., 2013) but also an interaction between regions of low and high stability.
The case study clearly differentiates among the synoptic features of the different regional dust patterns and confirms that low-pressure systems are the main systems associated with the frontal system, and the RST and its interactions with the Azores and Siberian ridges are the main synoptic features associated with the RST system. The deeper southern AP low-pressure system and its interaction with the Azores system and the high-pressure center over the northeastern AP are the main synoptic features associated with the thermal low system.
In addition, although the maximum wind forms a trough shape in all cases, the extended position of the core follows the dust region. The core is located over the eastern AP in the RST system, over the northern AP in the frontal system and over Africa and the northern AP in the thermal low system.
At the upper levels, the distribution of the atmospheric systems associated with dust storms in the different AP regions indicated that the northern cyclone trough enhanced the surface low pressure in the frontal system and that the extended tropical cyclone trough over the western region enhanced the RST at the surface in the RST system. However, the southern low pressure in the thermal low system was accompanied by a ridge between the trough over the AG and the trough over the central AP, enhancing the northwesterly wind over the eastern AP.
5. Conclusion
The spatial distribution of dust and dust storms over the AP demonstrated that most of the area north of24º N is strongly affected by dust phenomena, as is the Wadi-Aldawaser station, which is close to the desert. The temporal distribution demonstrated two periods for dust (dust storms), from 1978 to 1995 and from 1995 to 2008, with a high increasing rate and a low increasing (or decreasing) rate, respectively.
Synoptically, the spatial distribution of the atmospheric dust systems can be divided into three different systems with different characteristics, each mainly influencing the climate regions of the AP.
The northern system (or frontal system) consists of a surface low-pressure system, high-level troughs and a jet stream at the head of the troughs. This type of system forms an integrated system in which the maximum wind enhances the mid-level upward motion and increases the lifting force of the dust blown from local deserts, such as the Ad Dahna and the Syrian and Iraq deserts.
In the western system (or RST system), the RST and western trough of the southern low-pressure system were connected and deepened in the atmosphere via the northern extension of the upper tropical cyclone, and the trough-shaped maximum wind and upper-level trough were located over the western region (Shalaby et al., 2015; Mashat et al., 2019).
In the eastern system (or thermal low system), the southern low-pressure system formed a pronounced trough over the AG, accompanied by a subtropical anticyclonic center in the upper layers, which enhanced the thermal low (El Kenawy et al., 2014; El Kenawy and McCabe, 2016) or increased the surface temperature (Babu et al., 2016; Hasanean and Almazroui, 2016). Additionally, the pressure and geopotential gradients were oriented north-south, which might explain the formation of southerly shamal winds (Shao et al., 2011; Yu et al., 2016; Rashki et al., 2019) or high wind speeds (El-Sabh and Murty, 1989; Beegum et al., 2018; Hermida et al., 2018).
Widespread cases and previous studies (e.g., Goudi and Middleton, 2006; Hamidi et al., 2013; Mashat et al., 2018) have indicated that a combination of more than one of these systems could influence dust systems over the AP.
Generally, all cases are accompanied by pressure and geopotential gradients as factors involved in the generation of the wind responsible for emitting dust from the sources and transporting it to the affected regions (El-Sabh and Murty 1989; Hamidi et al., 2014, Beegum et al., 2018; Hermida et al., 2018). The shape and core of the maximum wind at250 hPa follow the main atmospheric systems accompanying dust storm events in the different regions. The southern low-pressure system is the common atmospheric feature accompanying all dust systems. If this low pressure is located over Rub Al Khali or over another dominant dust area in the AP, local dust-source contributions are expected to influence the AP in all dust systems (Prospero et al., 2002; Giles, 2005; Notaro et al., 2013).
The atmospheric mechanisms influencing the dust and the expected dust sources in the Sahara, Syria, Iraq, and northern and eastern AP are considered. The dust systems over different regions can be classified into two mechanisms: dust-generating systems (frontal and thermal low systems) or dust-transporting systems (RST systems).