# Spatiotemporal Variation of Particulate Fallout Instances in Sfax City, Southern Tunisia: Influence of Sources and Meteorology

Authors: Moez Bahloul; Iness Chabbi; Ali Sdiri; Ridha Amdouni; Khaled Medhioub; Chafai Azri
Journal: Advances in Meteorology
 (2015)

---

## Abstract

Particles deposition in the main industrial zone of Sfax City (southern Tunisia) was studied by weekly monitoring particulate fallout instances at twenty sites from November 11, 2012, to April 15, 2013. Very high fluctuation in those particle fluxes, ranging from 0.376 to 9.915 g/m2, was clearly observed. Spatiotemporal distribution of the deposited particulate fluxes and the exposure of each site to the main industrial plumes (i.e., phosphate treatment plant “SIAPE,” soap industry “SIOS-ZITEX,” and lead secondary melting industry “FP Sfax Sud”) indicated the concomitant effects of surrounding industrial sources. In addition, the highest particulate deposition seemed to be associated with predominant strong cyclonic situations. Those deposition rates exceeded the levels recorded in the case of strong stabilities, considered as responsible for pollutant accumulation.


---

## Body

## 1. Introduction

Air pollution usually results from diverse factors including mining activities, metallurgy, chemical industries, road and air traffic, household waste incineration, and industrial wastes [1, 2]. The composition of various gaseous and particulate pollutants that are emitted to the atmosphere depended on the origin, the emission conditions (i.e., flow of distribution near sources), and the exposure to the industrial plumes. In addition to the above-mentioned factors, weather conditions (e.g., wind speed and direction, temperature inversion, anticyclones, and cyclones) may effectively contribute to air pollution. It is well known that pollutants distribution in the atmosphere depended on dispersal phenomena that are, in turn, directly related to meteorological parameters (i.e., daylight hours, wind, temperature and humidity of the air, atmospheric stability, and precipitations [3]). During the atmospheric transport, the composition of original pollutants can be modified by mixtures of substances and aggregates from different origins, with various sizes and chemical compositions depending on the sources, the weather, and the environment [4]. In fact, chemical reactions between pollutants and atmospheric compounds may take place to produce secondary pollutants. In these conditions, the origin of atmospheric pollution becomes quite ambiguous [5]. Therefore, it may be more convenient to study not only the atmospheric composition but also the fallout settling. Within this context, the present work has been undertaken to evaluate the atmospheric particulate fallout instances near the main industrial zone of Sfax City (southern Tunisia). To achieve this goal, solar saltern was chosen as receiving environment because of (1) the homogenous support of the deposits (aquatic surface) and (2) the exposition to the main sources of pollution.


## 2. Material and Methods

Spatial and temporal variation of the deposition fluxes above Sfax solar saltern was recorded by monitoring the weekly rate deposition of particulates using collectors installed in twenty sites (i.e., PI2, TS2, TS35, TS36, TS38, S6, TS3, TS16, TS20, TS25, PM2, R2, R3, PI1, TS7, TS10, TS12, TS13, TS32, and TS42, Figure1). These site abbreviations were conserved as predetermined by the technical staff of the solar saltern. Study samples were collected during the period extending from November 11, 2012, to April 15, 2013, during which Sfax area received limited rainfall storms and low Saharan dusts input [22]. Sampling procedure was carried out according to the NF standard 43-007, 2008, by using DIEM plates. “Air Quality—Ambient Air—Determination of the Mass of Dry Atmospheric Depositions—Sampling on Deposit Plates—Preparation and Treatment” clearly described the placement of plates, coated with petroleum jelly, to analyze and subsequently treat the collected samples. The particulate phase dry deposition flux was measured using two smooth deposition plates (22.2 × 7.5 cm2 for each one) that were pointed at 1.5-meter height into a metal holder. The dimensions of each greased strip placed on the deposition plate were 10 × 5 cm2. Two greased strips mounted on deposition plates with a total collection area of 100 cm2 were used for each deposition sample. This type of deposition plate was used successfully as a surrogate surface by others to directly measure particulate dry deposition [23, 24].Figure 1
Map of the study area and sampling location in Sfax, Tunisia.For the purpose of the present study, deposition plates and strips were washed with alcohol and dichloromethane to remove all the impurities. Finally, strips were dried at 105°C, adequately numbered, and stored in separate containers for subsequent use. Those strips are made of aluminum (5 × 10 cm2) and covered with petroleum jelly. According to the NF standard 43-007, 2008, strips have to be placed. Afterward, GPS coordinates of all sampling points were registered for further mapping.After the sampling period of 7 days, strips were collected and submitted to the following laboratory treatments to collect dust samples:(i)
Connect the Erlenmeyer flask with screw to the pump.
(ii)
Weigh the virgin filter in a precision balance (three weightings).
(iii)
Place the filter on top of the Buckner funnel; place the top part and the spring that allows the joining of both parts.
(iv)
Using a pipette, place a small amount of dichloromethane on the surface of the strip containing dust.
(v)
Scrape with a spatula.
(vi)
When the strip is completely cleared, proceed to the filtration connecting the pump.
(vii)
Remove the filter from the Buckner funnel and weigh the sample with a precision balance (three weightings).Nine campaigns of measure were carried out by covering the selected period on a weekly basis. The weekly fluxes of atmospheric fallout instances were expressed in g/m2. Meteorological data such as air temperature, humidity, atmospheric pressure, wind velocity, and dominant wind direction was obtained from the meteorological station of Sfax Airport, the closest station to the sampling sites.


## 3. Climatic Characteristics of Sfax City

Sfax, a southern city of Tunisia (latitude 34°43′N; longitude 10°46′E), borders the Mediterranean Sea (Figure1). Its semiarid climate is therefore influenced by marine exposure generating a well ventilated area by low to moderate wind, rarely exceeding 5 m/s [22]. However, the influence of both continental and maritime air masses may lead to hot and dry summer with relatively cold and wet winter. Average temperature was estimated to be 12.3°C for winter season and 24.9°C for summer. It is to be mentioned that hot season is much longer than cold one, lasting more than five months per year (May—September) with an average temperature >24.2°C. As for rainfall, the region received about 217 mm/year with high seasonal and annual fluctuations.


## 4. Industrial Activities

Sfax is considered as the most industrial city of Tunisia. The main industrial activities included phosphate treatment (“SIAPE” plant), soap manufacturing (“SIOS-ZITEX”), and lead secondary melting industry (“FP Sfax Sud”). Most of those activities are located to the southern edge of the city; they threatened the environment because of their high emission of sulfur oxides (SOx) and particulate matter that largely exceeded the Tunisian standards [21]. SIAPE factory released about 4.5 tons/day of particulate matter with high amounts of sulphate, phosphorous compounds, and SOx, largely exceeding the permissible emission standards (8 and 20 times for particulate matter and SOx, resp.). For instance, soap industry discharged about 5.2 tons per day of particulate matter (>10 times the permissible emission standards value). Concerning the main source of particulate matter emissions, one can adopt the following sequence: SIOS-ZITEX (46%) > SIAPE (39%) > lead melting industry (15%) (Figure 2). In this regard, Azri et al. [25] stated that SIAPE and the lead melting industry were the main sources of heavy metals, as confirmed by the high enrichment factor of particulate matter with regard to heavy metals. They demonstrated that SIAPE was the main source of heavy metals (Zn, Ni, Cd, and Cu) with more than 9 g of Zn per ton of treated raw materials (equivalent to 8 kg/day of Zn emitted in the atmosphere); the emission factors of Ni, Cd, and Cu ranged between 2 and 3.2 g/ton. As for lead foundry, the emission of Pb was evaluated to be 68 kg/ton (equivalent to 200 kg/day of Pb released in the atmosphere). In contrast, the emission factors of the selected metals (i.e., Pb, Zn, Ni, Cd, and Cu) by the SIOS-ZITEX plant were insignificant.Figure 2
Annual fluxes of particulate emissions from the main pollutant sources [21].


## 5. Results and Discussions

Atmospheric particulate fallout instances as well as their spatiotemporal distribution were studied in detail to find out their most significant impact on the southern urban zone of Sfax area, southern Tunisia. Depending on current pollutants sources and weather conditions, the obtained results are described below.


### 5.1. Distribution of the Atmospheric Particulate Fallout Instances

Atmospheric particulate fallout instances recorded in the southern area of Sfax City showed highly variable values; the weekly flow values ranged from 0.376 to 9.915 g/m2 (Table 1) with an average of 1.750 g/m2. This value was lower than the permissible values fixed by the French (AFNOR) and German (TA-LUFT) standards (i.e., 7 and 2.45 g/m2 for AFNOR and TA-LUFT, resp.). The maximum weekly flux largely exceeded the relevant standards (1.5 to 4 times) in several sites. Furthermore, the registered fluxes were significantly different from those reported in other studies (Table 2) [6–20]. Based on their behavior (i.e., trends and amplitudes), spatial distribution of particulate fallout fluxes can be classified into three main groups (Figure 3):(i)
The first group included TS3, TS16, TS20, TS25, PM2, R2, R3, and PI1 study sites with weekly fluxes ranging between 0.376 and 3.76 g/m2. It presented about 42% of the studied sites with somewhat similar trends.
(ii)
The second group represented about 26% of the studied sites (i.e., PI2, TS35, TS36, TS38, and S6) with higher weekly fluxes (about 9.915 g/m2) and high variable trends.
(iii)
The third group covered 32% of the studied sites (i.e., TS7, TS10, TS2, TS12, TS13, TS32, and TS42) characterized by moderate fluxes (5.827 g/m2) combined with moderately variable trends.Table 1
Variation of the weekly fluxes (g/m2) of dry particulate fallout during the period of study.
Sites
Campaigns
Weekly average fluxes (g/m2)
C1
C2
C3
C4
C5
C6
C7
C8
C9
TS3
1.128
2.444
0.940
2.256
2.256
0.376
1.128
0.376
0.752
1.295
TS16
0.799
2.632
0.752
1.128
0.752
0.376
0.752
0.564
0.376
0.903
TS20
0.564
3.008
0.940
1.316
1.880
1.504
0.564
0.752
0.940
1.274
TS25
0.470
2.256
1.128
1.504
3.760
0.940
0.376
0.752
0.752
1.326
PM2
0.376
1.504
1.034
1.410
1.880
1.222
1.410
0.752
1.128
1.190
R2
0.752
1.128
1.010
0.376
3.759
0.376
0.376
0.940
0.376
1.010
R3
0.752
2.256
0.987
0.893
1.880
0.799
2.257
0.752
0.752
1.258
PI1
3.383
0.376
1.445
0.634
0.752
1.504
5.639
1.128
0.376
1.693
PI2
1.880
5.146
1.880
0.505
2.256
0.940
2.256
8.271
1.880
2.779
TS35
1.880
9.915
1.662
0.570
1.504
1.034
3.947
4.699
1.128
2.927
TS38
1.880
8.717
1.771
1.601
0.940
1.128
3.102
6.485
1.316
2.993
S6
1.316
7.519
1.717
2.632
0.376
1.081
3.524
5.592
1.504
2.807
TS2
2.162
0.376
0.376
1.504
0.752
5.639
1.128
1.125
1.316
1.597
TS7
2.232
2.820
1.880
0.752
1.504
3.383
0.376
0.376
0.752
1.564
TS10
2.303
5.263
0.376
2.256
1.128
1.128
1.504
0.750
1.504
1.801
TS12
2.444
2.632
0.752
0.752
0.728
0.376
0.752
7.519
1.128
1.146
TS13
1.880
3.806
3.008
0.940
0.705
0.887
0.376
1.498
1.504
1.622
TS32
0.752
5.827
1.730
0.376
0.752
1.504
3.419
1.128
0.752
1.804
TS36
0.752
4.135
1.744
1.504
0.564
1.292
3.313
3.360
1.128
1.977
TS42
3.007
4.981
2.369
1.128
0.658
1.398
1.897
2.244
0.752
2.048
C1: from 11 to 17 November 2012.
C2: from 25 November to 1 December 2012.
C3: from 2 to 8 December 2012.
C4: from 30 December 2012 to 5 January 2013.
C5: from 27 January to 2 February 2013.
C6: from 3 to 9 February 2013.
C7: from 24 February to 2 March 2013.
C8: from 12 to 18 March 2013.
C9: from 9 to 15 April 2013.Table 2
Fluxes of particulate deposits measured at different conditions from several areas over the world.
Sites
Weekly dry particulate deposits measured at different conditions (g/m2)
Authors
Southern urban Sfax
0.376–9.915
This work
Urban Constantine (Algeria)
7.735
[6]
Urban Dakar
3.99
[7]
Urban New Zealand
5.131–37.681
[8]
Urban British Columbia (Canada)
2.541
[9]
Urban Australia
13.223
[10]
Urban China
4.669–17.892
[11]
Urban Miami, Florida,USA
0.273
[12]
Urban Sapporo, Japan
1.211
[13]
La Mède (France)Urban zone
0.166
[14]
Kozani (western Greece)Urban zone
1.33
[15]
Vegoritis (western Greece)Rural zone
0.77
[15]
Petrana(western Greece)Rural zone
0.791
[15]
Mallorca(Balearic Islands, Spain) Rural zone
0.119
[16]
Tafira, Gran Canaria, Canary Islands (Spain)Rural zone
0.175
[17]
Montiers-sur-Saulx (France)Rural zone
0.031
[18]
Marienau, Forbach (France)Industrial zone
1.435
[19]
Quan-Xin (Taiwan)Industrial zone
2.192
[20]Figure 3
Spatial variation of deposits particles fluxes above the solar saltern.Cluster analysis was carried out using the ITCF statistical software package [26] to classify variables (i.e., sites of collection) based on the average linkage method. Pearson correlation coefficients were also used instead of the Euclidean distance that may lead to erroneous conclusions [27]. The establishment of dendrogram related to the study sites allowed deciphering the same main groups previously identified (Figure 4).Figure 4
Classification of studied sites.Figure5 illustrated the average flow of each site, showing a significant spatial variability of deposits. The linear regression between particulate deposit fluxes and their respective average deviations (Figure 6) shows that the higher the flux, the higher its standard deviation. The exposure to the industrial plumes (i.e., SIOS-ZITEX, FP Sfax Sud, and SIAPE) of each deposit site indicated that the highest fluxes values were recorded in PI2, TS35, TS36, TS38, and S6 study sites; they are strongly associated with the particulate fallout instances of SIAPE (Figure 7). It is worth noting that the exposure frequency was calculated through the occurrence of industrial plumes crossing each site downstream dominant winds. The methodology was based, first, on the emplacement of the center of the study period dominant wind rose at each point representing the industrial sources presented in Figure 1, second, on the computation of the frequency (%) of dominant wind directions joining each industrial source and study sites receiving maxima of deposit fallout instances, and, finally, on the establishment of the industrial plume trajectories (the plume axe was chosen as superimposed to the dominant wind direction) simulated by screen 3 software [28] for identifying the other exposed sites. Therefore, they are the most threatened by the aspect of the particulate deposits (quantitative and qualitative aspects). Further examination of the studied sites indicated that “SIAPE” was the main source of metals released in the atmosphere of the Sfax urban zone [25].Figure 5
Variation of weekly average fluxes (and their standard deviations) of deposits particles.Figure 6
Linear regression between the weekly averages fluxes of deposits and their standard deviations.Figure 7
Spatial variation of weekly average fluxes of deposits and frequencies of exposures to the industrial plumes.Temporal variation of particulate deposit fluxes demonstrated a weekly average value of 0.376 to 9.915 g/m2, showing large fluctuations (Figure 8) within the same site and between sites. The increase of particulate deposit fluxes can be attributed to the combined effect of both neighboring sources and the airflow properties which has been for a long time governed by meteorological conditions.Figure 8
Temporal variation of deposits particles fluxes above the solar saltern.


### 5.2. Influence of the Regional Meteorology and the Effect of the Particular Conditions

Figure8 showed that temporal particulate fallout instances are highly variable. Based on particulate deposits, two main periods (i.e., C2 and C8) can be easily recognized for their high fluxes usually associated with the predominant strong cyclonic episodes. Those fluxes values were much higher than those recorded under stable conditions, supposed to be favourable for the accumulation of various pollutants [29]. During C2 and C8 periods, low surface pressures over Sfax (associated with relatively higher wind speeds) were predominant (Figures 9 and 10); unstable weather conditions can be adapted as possible explanation (Figure 11). The meteorological factor in favor of high particles deposition fluxes in these periods can be found, on the one hand, in the highest frequencies of industrial plumes, especially of SIAPE (23% of total observations) crossing the sites (PI2, TS35, TS36, TS38, and S6) (Figure 7), associated with the relative increase of wind speeds reaching 8 m/s. On the other hand, these high fluxes can be accentuated by emission conditions and the existence of obstacles. It was proven in the case of plumes emitted at low altitudes that the wind increases the deposition of particles instead of insuring their dispersion, creating hence risk zones of locally increasing pollution levels [30]. The presence of various obstacles enhanced particles deposition speed [31–33]. Previous works at urban Sfax carried out by JICA [34] proved that the actual heights of SIAPE chimneys (between 30 and 70 m) were lower than those necessary for the good diffusion of emissions (equivalent height ≥ 100 m). Furthermore, the existence of phosphogypsum deposit (resulting from processing plants of phosphates) having 20 m height and an area of many hectares, located on the side of Sfax solar saltern at 200 m of the SIAPE, would constitute an obstacle against the emitted effluents (especially by southwest wind directions). Studies carried out by Azri et al. [21] proved that, under extremely unstable atmospheric conditions associated with relatively high wind velocities reaching 8 m/s, very high amounts of various gaseous and particulate pollutants were recorded at both upstream and downstream phosphogypsum deposit, especially at exposed sites to SIAPE plume. Concentrations of such pollutants can be multiplied by a factor of 10. This result was also proved elsewhere [35]. Furthermore, these relatively high wind velocities may enhance earth dust transportation by exceeding the threshold velocity (7 m/s [35, 36]).Figure 9
Surface synoptic maps over Sfax on (a) 26/11/2012, (b) 27/11/2012, (c) 28/11/2012, (d) 29/11/2012, (e) 30/11/2012, and (f) 01/12/2012.
(a)
(b)
(c)
(d)
(e)
(f)Figure 10
Surface synoptic maps over Sfax on (a) 13/03/2013, (b) 14/03/2013, (c) 15/03/2013, (d) 16/03/2013, and (e) 17/03/2013.
(a)
(b)
(c)
(d)
(e)Figure 11
Variations of atmospheric surface pressure and wind speed during the study period.Succinct analysis of the synoptic maps revealed that, since November 26, 2012 (i.e., C2 period), Tunisia has been under the influence of many successive surface depressions that took place over the Mediterranean Sea (Figure9). These depressions were associated with the Mediterranean front that separated the European cold air masses from the Saharan hot ones. A low-pressure zone over Europe (at 500 hPa) was characterized by strong winds. As for the temperature, it clearly decreased during the subsequent days (Figure 12). On November 26, 2012, the rapidly evolving surface depression moved towards the east of Italy. From November 27 to December 1, 2012, many successive surface depressions moved towards the east of the Mediterranean Sea.Figure 12
Temperature evolution at 850 hPa over Sfax in C2 period.On March 13, 2013 (i.e., C8 period), surface depression was centred over Spain. It was persistent over western Mediterranean Sea until March 16, 2013. Then, it moved on March 17 toward the central and eastern Mediterranean Sea. On March 14 and 15, 2013, the meteorological situation can be described as a vertical superposition over western Spain of surface depressions and low geopotentials. The evolution of temperature at 850 hPa (as seen in C2 period) clearly witnessed a decline in the successive days (Figure13). This particular period, associated with relatively high wind speeds (8 m/s), was characterized, over Sfax, by a cloudy sky and low temperatures (a decrease of 3 to 6°C for maximum registered temperatures at 14:00 LT). These environmental conditions favored the dilution phenomenon that will lease to the decrease of the air constituent concentrations, except for ozone concentrations enriched by the landing stratospheric air masses [37, 38]. However, significant increase of particulate deposition was observed over PI2, TS35, TS36, TS38, and S6 study sites, as seen in C2 situation. This can be explained by high frequencies of industrial plumes especially of SIAPE (23% of total observations) and the increase of surface wind velocities (about 8 m/s).Figure 13
Temperature evolution at 850 hPa over Sfax in C8 period.On March 16, 2013, the surface depression and the low geopotential area, under the dominance of western and northwestern winds, moved to the east. Since March 18 2013, it moved towards Eastern Europe while Tunisia returned to stable weather.Results presented above were refined by a factorial analysis of correspondences applied both to atmospheric particulate fallout instances collected at the study sites (d
(
P
I
2
), d
(
T
S
2
), d
(
T
S
35
), d
(
T
S
36
), d
(
T
S
38
), d
(
S
6
), d
(
T
S
3
), d
(
T
S
16
), d
(
T
S
20
), d
(
T
S
25
), d
(
P
M
2
), d
(
R
2
), d
(
R
3
), d
(
P
I
1
), d
(
T
S
7
), d
(
T
S
10
), d
(
T
S
12
), d
(
T
S
13
), d
(
T
S
32
), and d
(
T
S
42
)) through the nine campaigns and to binary numbers related to registered classes for atmospheric surface pressure, wind velocity, wind sectors, air temperature, and humidity. Three classes P1–P3 related to the atmospheric surface pressure were selected (P1 ∈ [1,002; 1,014 hPa]; P2 ∈ [1,014; 1,020 hPa]; P3 ∈ [1,020; 1,030 hPa]). For the wind velocity, three classes V1–V3 were chosen (V1 ∈ [1; 3 m/s]; V2 ∈ [3; 5 m/s]; V3 ∈ [5; 8 m/s]). Four wind sectors S1–S4 were also selected (S1 ∈ [0; 90°]; S2 ∈ [90; 180°]; S3 ∈ [180; 270°]; S4 ∈ [270; 360°]). For the temperature, two classes T1-T2 were chosen (T1 ∈ [11; 17.5°C]; T2 ∈ [17.5; 24°C]). For the relative humidity, four classes have been selected (RH1 ∈ [34; 43%]; RH2 ∈ [43; 52%]; RH3 ∈ [52; 61%]; RH4 ∈ [61; 70%]).The projection over the 1 × 2 factorial plane (presenting the maximum of inertia) of all selected particulate deposit fluxes and meteorological variables shows distinct data groups (Figure14):(i)
Group G1 was characterized by the concomitant effect of the predominant strong cyclonic episodes and SIAPE effluent fallout instances distinguished by an unstable atmosphere (P1) associated with relatively high velocities (V3) favourable to the accentuation deposit fallout instances in PI2, TS35, TS36, TS38, and S6 study sites. The association of these component parameters is well pronounced under the predominance of the western wind’s sector (S3) which drained the industrial plumes, especially of SIAPE. This group was more pronounced in C2 and C8 campaigns.
(ii)
Group G2 was representative of the periods characterized by the effect of the steady to very steady atmosphere. It was pronounced in most campaigns (except for C2 and C8). During these campaigns, atmospheric surface pressure was higher than 1,014 hPa (P2 and P3), implying a stable to very stable atmosphere. Wind speeds are relatively lower and reach values oscillating between 1 and 5 m/s (V1 and V2). The deposit fallout instances (d
(
T
S
3
), d
(
T
S
16
), d
(
T
S
20
), d
(
T
S
25
), d
(
P
M
2
), d
(
R
2
), d
(
R
3
), d
(
P
I
1
), d
(
T
S
7
), d
(
T
S
10
), d
(
T
S
2
), d
(
T
S
12
), d
(
T
S
13
), d
(
T
S
32
), and d
(
T
S
42
)) in the remaining study sites (compared to those of the first group) were shown to be influenced by the dominance of the eastern, southern, and western wind’s sector (S1, S2, and S4). The temperature and the humidity were shown without significant effect.Figure 14
Projection of variables in the 1 × 2 factorial plane (representing 57.03% of the total variance) during the study period (threshold of significance = 0.5 forp
<
0.05 and n
=
9).Significant adverse effects of the local industrial activities threatened the environment in Sfax area, especially the southern edge of the city where several activities are being performed. Therefore, it is necessary to undertake a sustainable management policy that would meet the national environmental standards and regulations. More in-depth studies can provide evidence of the challenges related to environmental sustainability and the implementation of active measures to reduce the adverse effects of particulates emission. It is therefore recommended to take the following considerations into account:(i)
Ensuring adequate treatment of the industrial emissions, especially for SIAPE factory.
(ii)
Strengthening the control network of both dry and wet deposits in Sfax urban zones and its suburbs to master pollutants distribution.
(iii)
Taking all the measures required for the protection of the environment.


## 5.1. Distribution of the Atmospheric Particulate Fallout Instances

Atmospheric particulate fallout instances recorded in the southern area of Sfax City showed highly variable values; the weekly flow values ranged from 0.376 to 9.915 g/m2 (Table 1) with an average of 1.750 g/m2. This value was lower than the permissible values fixed by the French (AFNOR) and German (TA-LUFT) standards (i.e., 7 and 2.45 g/m2 for AFNOR and TA-LUFT, resp.). The maximum weekly flux largely exceeded the relevant standards (1.5 to 4 times) in several sites. Furthermore, the registered fluxes were significantly different from those reported in other studies (Table 2) [6–20]. Based on their behavior (i.e., trends and amplitudes), spatial distribution of particulate fallout fluxes can be classified into three main groups (Figure 3):(i)
The first group included TS3, TS16, TS20, TS25, PM2, R2, R3, and PI1 study sites with weekly fluxes ranging between 0.376 and 3.76 g/m2. It presented about 42% of the studied sites with somewhat similar trends.
(ii)
The second group represented about 26% of the studied sites (i.e., PI2, TS35, TS36, TS38, and S6) with higher weekly fluxes (about 9.915 g/m2) and high variable trends.
(iii)
The third group covered 32% of the studied sites (i.e., TS7, TS10, TS2, TS12, TS13, TS32, and TS42) characterized by moderate fluxes (5.827 g/m2) combined with moderately variable trends.Table 1
Variation of the weekly fluxes (g/m2) of dry particulate fallout during the period of study.
Sites
Campaigns
Weekly average fluxes (g/m2)
C1
C2
C3
C4
C5
C6
C7
C8
C9
TS3
1.128
2.444
0.940
2.256
2.256
0.376
1.128
0.376
0.752
1.295
TS16
0.799
2.632
0.752
1.128
0.752
0.376
0.752
0.564
0.376
0.903
TS20
0.564
3.008
0.940
1.316
1.880
1.504
0.564
0.752
0.940
1.274
TS25
0.470
2.256
1.128
1.504
3.760
0.940
0.376
0.752
0.752
1.326
PM2
0.376
1.504
1.034
1.410
1.880
1.222
1.410
0.752
1.128
1.190
R2
0.752
1.128
1.010
0.376
3.759
0.376
0.376
0.940
0.376
1.010
R3
0.752
2.256
0.987
0.893
1.880
0.799
2.257
0.752
0.752
1.258
PI1
3.383
0.376
1.445
0.634
0.752
1.504
5.639
1.128
0.376
1.693
PI2
1.880
5.146
1.880
0.505
2.256
0.940
2.256
8.271
1.880
2.779
TS35
1.880
9.915
1.662
0.570
1.504
1.034
3.947
4.699
1.128
2.927
TS38
1.880
8.717
1.771
1.601
0.940
1.128
3.102
6.485
1.316
2.993
S6
1.316
7.519
1.717
2.632
0.376
1.081
3.524
5.592
1.504
2.807
TS2
2.162
0.376
0.376
1.504
0.752
5.639
1.128
1.125
1.316
1.597
TS7
2.232
2.820
1.880
0.752
1.504
3.383
0.376
0.376
0.752
1.564
TS10
2.303
5.263
0.376
2.256
1.128
1.128
1.504
0.750
1.504
1.801
TS12
2.444
2.632
0.752
0.752
0.728
0.376
0.752
7.519
1.128
1.146
TS13
1.880
3.806
3.008
0.940
0.705
0.887
0.376
1.498
1.504
1.622
TS32
0.752
5.827
1.730
0.376
0.752
1.504
3.419
1.128
0.752
1.804
TS36
0.752
4.135
1.744
1.504
0.564
1.292
3.313
3.360
1.128
1.977
TS42
3.007
4.981
2.369
1.128
0.658
1.398
1.897
2.244
0.752
2.048
C1: from 11 to 17 November 2012.
C2: from 25 November to 1 December 2012.
C3: from 2 to 8 December 2012.
C4: from 30 December 2012 to 5 January 2013.
C5: from 27 January to 2 February 2013.
C6: from 3 to 9 February 2013.
C7: from 24 February to 2 March 2013.
C8: from 12 to 18 March 2013.
C9: from 9 to 15 April 2013.Table 2
Fluxes of particulate deposits measured at different conditions from several areas over the world.
Sites
Weekly dry particulate deposits measured at different conditions (g/m2)
Authors
Southern urban Sfax
0.376–9.915
This work
Urban Constantine (Algeria)
7.735
[6]
Urban Dakar
3.99
[7]
Urban New Zealand
5.131–37.681
[8]
Urban British Columbia (Canada)
2.541
[9]
Urban Australia
13.223
[10]
Urban China
4.669–17.892
[11]
Urban Miami, Florida,USA
0.273
[12]
Urban Sapporo, Japan
1.211
[13]
La Mède (France)Urban zone
0.166
[14]
Kozani (western Greece)Urban zone
1.33
[15]
Vegoritis (western Greece)Rural zone
0.77
[15]
Petrana(western Greece)Rural zone
0.791
[15]
Mallorca(Balearic Islands, Spain) Rural zone
0.119
[16]
Tafira, Gran Canaria, Canary Islands (Spain)Rural zone
0.175
[17]
Montiers-sur-Saulx (France)Rural zone
0.031
[18]
Marienau, Forbach (France)Industrial zone
1.435
[19]
Quan-Xin (Taiwan)Industrial zone
2.192
[20]Figure 3
Spatial variation of deposits particles fluxes above the solar saltern.Cluster analysis was carried out using the ITCF statistical software package [26] to classify variables (i.e., sites of collection) based on the average linkage method. Pearson correlation coefficients were also used instead of the Euclidean distance that may lead to erroneous conclusions [27]. The establishment of dendrogram related to the study sites allowed deciphering the same main groups previously identified (Figure 4).Figure 4
Classification of studied sites.Figure5 illustrated the average flow of each site, showing a significant spatial variability of deposits. The linear regression between particulate deposit fluxes and their respective average deviations (Figure 6) shows that the higher the flux, the higher its standard deviation. The exposure to the industrial plumes (i.e., SIOS-ZITEX, FP Sfax Sud, and SIAPE) of each deposit site indicated that the highest fluxes values were recorded in PI2, TS35, TS36, TS38, and S6 study sites; they are strongly associated with the particulate fallout instances of SIAPE (Figure 7). It is worth noting that the exposure frequency was calculated through the occurrence of industrial plumes crossing each site downstream dominant winds. The methodology was based, first, on the emplacement of the center of the study period dominant wind rose at each point representing the industrial sources presented in Figure 1, second, on the computation of the frequency (%) of dominant wind directions joining each industrial source and study sites receiving maxima of deposit fallout instances, and, finally, on the establishment of the industrial plume trajectories (the plume axe was chosen as superimposed to the dominant wind direction) simulated by screen 3 software [28] for identifying the other exposed sites. Therefore, they are the most threatened by the aspect of the particulate deposits (quantitative and qualitative aspects). Further examination of the studied sites indicated that “SIAPE” was the main source of metals released in the atmosphere of the Sfax urban zone [25].Figure 5
Variation of weekly average fluxes (and their standard deviations) of deposits particles.Figure 6
Linear regression between the weekly averages fluxes of deposits and their standard deviations.Figure 7
Spatial variation of weekly average fluxes of deposits and frequencies of exposures to the industrial plumes.Temporal variation of particulate deposit fluxes demonstrated a weekly average value of 0.376 to 9.915 g/m2, showing large fluctuations (Figure 8) within the same site and between sites. The increase of particulate deposit fluxes can be attributed to the combined effect of both neighboring sources and the airflow properties which has been for a long time governed by meteorological conditions.Figure 8
Temporal variation of deposits particles fluxes above the solar saltern.


## 5.2. Influence of the Regional Meteorology and the Effect of the Particular Conditions

Figure8 showed that temporal particulate fallout instances are highly variable. Based on particulate deposits, two main periods (i.e., C2 and C8) can be easily recognized for their high fluxes usually associated with the predominant strong cyclonic episodes. Those fluxes values were much higher than those recorded under stable conditions, supposed to be favourable for the accumulation of various pollutants [29]. During C2 and C8 periods, low surface pressures over Sfax (associated with relatively higher wind speeds) were predominant (Figures 9 and 10); unstable weather conditions can be adapted as possible explanation (Figure 11). The meteorological factor in favor of high particles deposition fluxes in these periods can be found, on the one hand, in the highest frequencies of industrial plumes, especially of SIAPE (23% of total observations) crossing the sites (PI2, TS35, TS36, TS38, and S6) (Figure 7), associated with the relative increase of wind speeds reaching 8 m/s. On the other hand, these high fluxes can be accentuated by emission conditions and the existence of obstacles. It was proven in the case of plumes emitted at low altitudes that the wind increases the deposition of particles instead of insuring their dispersion, creating hence risk zones of locally increasing pollution levels [30]. The presence of various obstacles enhanced particles deposition speed [31–33]. Previous works at urban Sfax carried out by JICA [34] proved that the actual heights of SIAPE chimneys (between 30 and 70 m) were lower than those necessary for the good diffusion of emissions (equivalent height ≥ 100 m). Furthermore, the existence of phosphogypsum deposit (resulting from processing plants of phosphates) having 20 m height and an area of many hectares, located on the side of Sfax solar saltern at 200 m of the SIAPE, would constitute an obstacle against the emitted effluents (especially by southwest wind directions). Studies carried out by Azri et al. [21] proved that, under extremely unstable atmospheric conditions associated with relatively high wind velocities reaching 8 m/s, very high amounts of various gaseous and particulate pollutants were recorded at both upstream and downstream phosphogypsum deposit, especially at exposed sites to SIAPE plume. Concentrations of such pollutants can be multiplied by a factor of 10. This result was also proved elsewhere [35]. Furthermore, these relatively high wind velocities may enhance earth dust transportation by exceeding the threshold velocity (7 m/s [35, 36]).Figure 9
Surface synoptic maps over Sfax on (a) 26/11/2012, (b) 27/11/2012, (c) 28/11/2012, (d) 29/11/2012, (e) 30/11/2012, and (f) 01/12/2012.
(a)
(b)
(c)
(d)
(e)
(f)Figure 10
Surface synoptic maps over Sfax on (a) 13/03/2013, (b) 14/03/2013, (c) 15/03/2013, (d) 16/03/2013, and (e) 17/03/2013.
(a)
(b)
(c)
(d)
(e)Figure 11
Variations of atmospheric surface pressure and wind speed during the study period.Succinct analysis of the synoptic maps revealed that, since November 26, 2012 (i.e., C2 period), Tunisia has been under the influence of many successive surface depressions that took place over the Mediterranean Sea (Figure9). These depressions were associated with the Mediterranean front that separated the European cold air masses from the Saharan hot ones. A low-pressure zone over Europe (at 500 hPa) was characterized by strong winds. As for the temperature, it clearly decreased during the subsequent days (Figure 12). On November 26, 2012, the rapidly evolving surface depression moved towards the east of Italy. From November 27 to December 1, 2012, many successive surface depressions moved towards the east of the Mediterranean Sea.Figure 12
Temperature evolution at 850 hPa over Sfax in C2 period.On March 13, 2013 (i.e., C8 period), surface depression was centred over Spain. It was persistent over western Mediterranean Sea until March 16, 2013. Then, it moved on March 17 toward the central and eastern Mediterranean Sea. On March 14 and 15, 2013, the meteorological situation can be described as a vertical superposition over western Spain of surface depressions and low geopotentials. The evolution of temperature at 850 hPa (as seen in C2 period) clearly witnessed a decline in the successive days (Figure13). This particular period, associated with relatively high wind speeds (8 m/s), was characterized, over Sfax, by a cloudy sky and low temperatures (a decrease of 3 to 6°C for maximum registered temperatures at 14:00 LT). These environmental conditions favored the dilution phenomenon that will lease to the decrease of the air constituent concentrations, except for ozone concentrations enriched by the landing stratospheric air masses [37, 38]. However, significant increase of particulate deposition was observed over PI2, TS35, TS36, TS38, and S6 study sites, as seen in C2 situation. This can be explained by high frequencies of industrial plumes especially of SIAPE (23% of total observations) and the increase of surface wind velocities (about 8 m/s).Figure 13
Temperature evolution at 850 hPa over Sfax in C8 period.On March 16, 2013, the surface depression and the low geopotential area, under the dominance of western and northwestern winds, moved to the east. Since March 18 2013, it moved towards Eastern Europe while Tunisia returned to stable weather.Results presented above were refined by a factorial analysis of correspondences applied both to atmospheric particulate fallout instances collected at the study sites (d
(
P
I
2
), d
(
T
S
2
), d
(
T
S
35
), d
(
T
S
36
), d
(
T
S
38
), d
(
S
6
), d
(
T
S
3
), d
(
T
S
16
), d
(
T
S
20
), d
(
T
S
25
), d
(
P
M
2
), d
(
R
2
), d
(
R
3
), d
(
P
I
1
), d
(
T
S
7
), d
(
T
S
10
), d
(
T
S
12
), d
(
T
S
13
), d
(
T
S
32
), and d
(
T
S
42
)) through the nine campaigns and to binary numbers related to registered classes for atmospheric surface pressure, wind velocity, wind sectors, air temperature, and humidity. Three classes P1–P3 related to the atmospheric surface pressure were selected (P1 ∈ [1,002; 1,014 hPa]; P2 ∈ [1,014; 1,020 hPa]; P3 ∈ [1,020; 1,030 hPa]). For the wind velocity, three classes V1–V3 were chosen (V1 ∈ [1; 3 m/s]; V2 ∈ [3; 5 m/s]; V3 ∈ [5; 8 m/s]). Four wind sectors S1–S4 were also selected (S1 ∈ [0; 90°]; S2 ∈ [90; 180°]; S3 ∈ [180; 270°]; S4 ∈ [270; 360°]). For the temperature, two classes T1-T2 were chosen (T1 ∈ [11; 17.5°C]; T2 ∈ [17.5; 24°C]). For the relative humidity, four classes have been selected (RH1 ∈ [34; 43%]; RH2 ∈ [43; 52%]; RH3 ∈ [52; 61%]; RH4 ∈ [61; 70%]).The projection over the 1 × 2 factorial plane (presenting the maximum of inertia) of all selected particulate deposit fluxes and meteorological variables shows distinct data groups (Figure14):(i)
Group G1 was characterized by the concomitant effect of the predominant strong cyclonic episodes and SIAPE effluent fallout instances distinguished by an unstable atmosphere (P1) associated with relatively high velocities (V3) favourable to the accentuation deposit fallout instances in PI2, TS35, TS36, TS38, and S6 study sites. The association of these component parameters is well pronounced under the predominance of the western wind’s sector (S3) which drained the industrial plumes, especially of SIAPE. This group was more pronounced in C2 and C8 campaigns.
(ii)
Group G2 was representative of the periods characterized by the effect of the steady to very steady atmosphere. It was pronounced in most campaigns (except for C2 and C8). During these campaigns, atmospheric surface pressure was higher than 1,014 hPa (P2 and P3), implying a stable to very stable atmosphere. Wind speeds are relatively lower and reach values oscillating between 1 and 5 m/s (V1 and V2). The deposit fallout instances (d
(
T
S
3
), d
(
T
S
16
), d
(
T
S
20
), d
(
T
S
25
), d
(
P
M
2
), d
(
R
2
), d
(
R
3
), d
(
P
I
1
), d
(
T
S
7
), d
(
T
S
10
), d
(
T
S
2
), d
(
T
S
12
), d
(
T
S
13
), d
(
T
S
32
), and d
(
T
S
42
)) in the remaining study sites (compared to those of the first group) were shown to be influenced by the dominance of the eastern, southern, and western wind’s sector (S1, S2, and S4). The temperature and the humidity were shown without significant effect.Figure 14
Projection of variables in the 1 × 2 factorial plane (representing 57.03% of the total variance) during the study period (threshold of significance = 0.5 forp
<
0.05 and n
=
9).Significant adverse effects of the local industrial activities threatened the environment in Sfax area, especially the southern edge of the city where several activities are being performed. Therefore, it is necessary to undertake a sustainable management policy that would meet the national environmental standards and regulations. More in-depth studies can provide evidence of the challenges related to environmental sustainability and the implementation of active measures to reduce the adverse effects of particulates emission. It is therefore recommended to take the following considerations into account:(i)
Ensuring adequate treatment of the industrial emissions, especially for SIAPE factory.
(ii)
Strengthening the control network of both dry and wet deposits in Sfax urban zones and its suburbs to master pollutants distribution.
(iii)
Taking all the measures required for the protection of the environment.


## 6. Conclusions

The spatial and temporal evolution of the particulate fallout instances in southern urban area of Sfax was studied to find out its environmental impact on the neighboring areas. It was found that particulate fluxes showed significantly variable trends. Based on their behavior (trends and amplitudes), spatial distribution of particulate fallout fluxes can be classified into three main patterns:(i)
The first pattern, covering 42% of sites (TS3, TS16, TS20, TS25, PM2, R2, R3, and PI1), was characterized by fluxes ranging between 0.376 and 3.76 g/m2.
(ii)
The second pattern concerning 26% of sites (PI2, TS35, TS36, TS38, and S6) recorded high fluxes, reaching 9.915 g/m2.
(iii)
The third pattern was an intermediate group characterized by moderate fluxes reaching 5.827 g/m2; it represented 32% of sites (TS7, TS10, TS2, TS12, TS13, TS32, and TS42).Such distribution patterns seemed to be governed by the combined effects between surrounding industrial sources, exposure to industrial plumes, and local airflow characteristics.Under relatively predominant strong cyclonic situations, the increase in fluxes exceeded the levels recorded under the conditions of relatively strong stabilities that usually enhance pollutants accumulation. The meteorological factor in favour of high fluxes of particles deposition can be found in the highest frequencies of industrial plumes, associated with the relative increase of wind speeds and the effect of the phosphogypsum deposit constituting an obstacle against the emitted effluents.



---

*Source: 471396-2015-08-04.xml*