Experimental Assessment of Dust Emissions on Compacted Soils Degraded by Traffic

: Haul traffic on earthworks runways during construction works is an important factor of dust emission. Compacted soils surface become progressively degraded as the number of wheels passing increases. Fine particles are then segregated from the soil surface and lifted when the shear stress generated by the flow above the surface increases, leading to the worsening of air quality and reduction of visibility. Laboratory tests were performed to assess dust emissions on traffic degraded soils. Mixtures of kaolin clay and sand were compacted using a laboratory roller compactor and were degraded using a vehicle simulator. Models describing the evolutions of soil degradation and Particle Size Distribution (PSD) during traffic were established. Then, the velocity profiles above each soil sample were obtained in a wind tunnel. The experimental results were analyzed to determine the Reynolds shear stresses generated by the turbulence of the flow. PSD, degradation and stresses were implemented in the Convective Turbulent Dust Emission (CTDE) model to estimate the dust emission flux of the soils for several passes of the wheel. A comparison between results from the model and field measurements underlines that turbulence is not the main contributor to dust emissions when a vehicle is in motion.


Introduction
Vehicle traffic on compacted soils is one of the main contributors to particulate emission into the atmosphere. Excluding wind erosion, traffic on unpaved road accounts for almost 30% of dust emissions in the United States [1]. Circulation of haul trucks on earthworks runways generates also an important amount of dust. This is ranked as the second most common risk for many professionals of the construction industry [2]. Indeed, particles raised into suspension have both health and environmental impacts [3,4]. Moreover, a dust plume can reduce visibility [5,6], resulting in traffic hazards [7,8]. Soil watering is a good way to limit particle lift [9]. Nevertheless, it involves the use of large volumes of water. A better understanding of the unpaved roads dust emission mechanisms would help optimize soil watering.
The United States Environmental Protection Agency (USEPA) developed a model (called AP-42) to estimate PM10 (particle diameters less than or equal to 10 µm) emissions generated by vehicle traffic on unpaved roadways [10]. However, this model is purely empirical and shows significant discrepancies with in-situ measurements on rural unpaved roads [11] as well as on earthworks haul roads [12]. Studies have shown the effect of vehicle characteristics [11,13] and soil properties [14] on dust emission. These approaches, based on in-situ measurements, were mainly empirical. Therefore, the appearance of the particle by soil degradation and the associated lift are two mechanisms that are not yet sufficiently well-modeled.
To date, the most comprehensive analytical models of dust generation were developed in the context of wind erosion studies. When considering the soil-atmosphere interaction, many particles lift models are based on a threshold wind velocity friction [15−17]. However, this mechanism does not consider the turbulence. Therefore, it is relatively far from the process of particle entrainment by vehicle activity [14], since turbulence is often considered as the major factor influencing particles lift [11,18]. In order to quantify the dust emissions related to vehicle traffic, it, therefore, seems more relevant to use a turbulence-based model. In the present study, the model developed by Klose and Shao [19] is used as a basis to assess dust emissions.
Particles subjected to lift are generated by the abrasion of the soil surface due to repeated passes of vehicles [20]. However, based on the literature review, the process of soil surface degradation by tires has not yet been studied.
The preliminary aim of this study is to characterize experimentally the compacted soil degradation by traffic. The processes of particle generation and particle size evolution are studied. Three clay mixtures with different sand content were compacted using a laboratory roller compactor. Compacted soil samples were then degraded by a wheel in a traffic simulator. Stresses experienced by earthworks runways during truck circulation were reproduced. Models to quantify the evolution of soil degradation and Particle Size Distribution (PSD) during traffic were established.
The second objective of this study is to examine the interaction of the degraded soils with the atmosphere. Flow velocities were measured using LASER Doppler Velocimetry (LDV) at various locations above the soil samples in a wind tunnel. Reynolds shear stresses were determined near the soil surface to assess the lift forces experienced by the dust particles.
Degradation and PSD evolutions, as well as the particle cohesion and lift forces, were implemented into the Klose and Shao's model [19], in order to quantify the dust emission of the traffic degraded soils.
In section 2, the convective turbulent dust emission model and the related assumptions are detailed. In section 3, experimental facilities and measurement techniques are presented. In section 4, the results are detailed and discussed. In the last section, conclusions and perspectives are presented.

Description of the Klose and Shao's Model
Convective Turbulent Dust Emission (CTDE) model establishes that particle lift initiate when interparticle cohesive forces become smaller than lifting forces generated by turbulence [19]. According to a force balance, the emission flux φd (kg.m −2 .s −1 ) of particles having a diameter d can be expressed as Equation (1): where Nd is the particle number concentration per unit of volume (m −3 ), Tp is the particle response time (s), f is the lifting force expressed in newtons (N), fi is the interparticle cohesive force (N) and δvc is the thickness of the viscous sublayer area where the particles are sheared. δvc is given by Equation (2): with ν being the air kinematic viscosity (15.6 × 10 −6 m 2 .s −1 at 25 °C) and u * the friction velocity (m.s −1 ).
According to [21], Tp is given by Equation (3): with being the particle density (kg.m −3 ) and the air density (1.184 kg.m −3 at 25 °C). The force exerted on the particle is given by Equation (4): where is the Reynolds shear stress, that is, the instantaneous vertical flux of horizontal momentum given by Equation (5): In Equation (5), u', v' and w' are, according to the Reynold decomposition, the instantaneous fluctuations of the flow velocity (m.s −1 ) in the longitudinal, vertical and spanwise directions respectively. These fluctuations depend on time, so the Reynolds shear stresses are not constant. Thus, is a stochastic quantity that obeys a probability distribution p(τ) depending on the Probability Distribution Functions (PDF) of u', v' and w'.
The cohesive force (fi) mainly depends on the Van der Walls interactions, electrostatic forces as well as capillary and chemical binding forces [21]. These interactions are affected by many parameters including, but not limited to, the particle size, the particle shape, the mineral composition and the surface roughness. Therefore, cohesive forces are difficult to estimate. It is more convenient to treat these forces as stochastic variables following a probability distribution p(fi). Finally, according to [19], the total convective dust emission flux can be given by Equation (6): where p(d) is the PSD of the soil and p(f) is the probability distribution of f. This model has shown a good ability to predict dust emissions during atmospheric convective turbulent events in Taklimakan Desert (China), Horqin Sandy Land (China) and in the Murray-Darling River Basin (Australia). Coefficients of determination (r²) up to 0.71 were obtained when comparing model predictions with observations both at a site-based scale [22,23] and a regional scale [24].
In the present study, the originality of the approach is to use this model to quantify dust emissions on traffic degraded soils. Input parameters to be implemented in the model are: PSD, the quantity of particles and probability distributions of cohesive and lifting forces. The following section details the determination of these parameters.

Particle Size Distribution (PSD) and Quantity of Particles Subjected to Lift
The PSD of soil can be obtained by sedimentation [25], sieving [26] or by techniques such as LASER diffraction [27]. The latter technique was used in the present study. In the context of dust emission by vehicle-induced traffic, the PSD of soil particles that get lifted is different than the PSD of the initial soil. Indeed, due to the soil moisture content and the compaction load, compacted soil is composed of particles bonded together to form aggregates. These aggregates are pulled out of the compacted soil surface under the shearing action of a wheel rolling over the surface [28]. The PSD of these aggregates is comparatively coarser than that of the initial soil. By increasing the number of wheel passing, finer particles are formed by aggregate crumbling. This leads to an evolution of PSD over time. In the same way, the number of particles segregated from the surface evolves as the soil degradation by traffic increases. The experimental approach for the study of PSD and degradation evolutions during traffic is detailed in section 3.1.

Cohesive Forces
The stochastic behavior of interparticle cohesive forces was dealt with in the study of Zimon [29]. The experiments were performed with glass particles submitted to a centrifugal acceleration to induce their flight. By gradually increasing the centrifuge force and measuring the proportion of particles lifted, the statistical distribution of cohesive forces was estimated. According to the results obtained, the cohesive force can be described using a log-normal distribution defined by Equation (7): where is the median value of the cohesive force (fi) and its geometric standard deviation (N). The log-normal distribution for cohesive forces is also proposed by other researchers for other types of particles [30][31][32]. Dealing with Zimon's data, Shao [21] proposed an empirical law to estimate and as a function of particle size (Equation (8)): ( ) = 10 (10 × e .
where d is expressed in µm. Figure 1 shows the distributions of fi for particles having diameters of 1 µm, 2 µm, 3.5 µm, 7.5 µm and 20 µm. The geometric standard deviation for the log-normal distribution is inversely proportional to the particle diameter. The distribution functions p(fi) cover a wide range for fine particles, demonstrating the need to consider the stochastic behavior of cohesive forces when examining dust emission.

Lifting Forces
Lifting forces due to the wind were determined through wind tunnel experiments by measuring velocities of the air above traffic-degraded soils (section 3.2). However, for practical reason, only the longitudinal component u' was measured. In order to assess the Reynolds shear stresses, assumptions about vertical and transverse fluctuations were made and discussed below.
For simplification in the study, it is assumed that turbulence is isotropic (u' = v' = w'). In reality, the presence of a mean flow orientated in horizontal direction introduce anisotropies in the turbulent velocity, the longitudinal fluctuation u' being greater than other fluctuations. Thus, an idealised turbulent state is considered that overestimates the Reynolds shear stresses (conservative hypothesis). It is worth noting that this assumption is often used in numerical models [33−35]. Reynolds shear stresses are therefore expressed in a simpler way Equation (10): where τ is the product of two variables, both having a distribution function pu'. Thus, the Reynolds shear stresses distribution function p(τ) is a product distribution [36] (p. 160), given by Equation (11): where x is the set of values that can be taken by u'.
In the next section, the experimental facilities and measurements techniques are described.

Compaction
In civil engineering, soils are compacted to improve their bearing capacity [37]. A good quality of compaction is particularly important for soils submitted to stresses induced by vehicle traffic, such as unpaved roads or runways on earthworks site. In this study, soils were compacted using the roller compactor of the University Gustave Eiffel ( Figure 2). This device allows reproducing the stress tensor rotation in the soil during compaction, which is similar to field conditions [38]. Soil samples were prepared with mixtures of two materials: Kaolin clay (defined as K ) and Hostun HN38 sand (defined as S ). Three different mixtures were studied: a pure clay sample S0K100 (100% of clay) and two different sand-clay mixtures S50K50 (50% of clay and 50% of sand) and S75K25 (25% of clay and 75% of sand). These three soil mixtures cannot fully represent natural soils but are representative of the fine soils encountered usually on construction sites, and thus chosen herein. According to the French guide of road earthworks [39], S0K100 is classified A3 (medium-plastic clay), S50K50 is classified A2 (fine clayed sand) and S75K25 is classified A1 (fine sand). The PSD of the samples was measured by LASER diffraction (section 3.1.2). These soils were compacted in 500 × 180 × 28.5 mm 3 rectangular samples. The compacting process consisted of pouring the required quantity of soil to be compacted into the container placed inside the base. A displacement table gradually lifted the container while a smooth wheel moved over the surface until the soil was compacted to a height of 28.5 mm in the container.
Five samples (1 of S0K100, 1 of S50K50 and 3 of S75K25) were compacted at the Proctor optimum, which is the optimum moisture content for maximum soil density after compaction [40]. The compaction parameters are detailed in Table 1. After compaction, the soils were air-dried for 24 h and then degraded with a tire in a traffic simulator. All the tests were carried out during summer under the same weather conditions (approximately 25 °C during compaction and degradation).

Soil Degradation by Traffic Simulation
After compaction, the soil samples were degraded in the traffic simulator VECTRA at the University Gustave Eiffel. This apparatus (Figure 3a) simulates the motion of a wheel equipped with a 400 mm diameter tire. The movement is generated according to the following procedure: the continuous rotation of a gear motor unit is transformed into a reciprocating movement by an oscillating bearing carrying a telescopic arm. This arm transmits its movement to a mobile carriage that is mounted on a rail. The wheel is fixed to this mobile carriage and moves longitudinally along the axis of the rail by reciprocating passes over a distance of 0.410 ± 0.005 m with a frequency of 1.0 ± 0.1 Hz (Figure 3b). This corresponds to an average speed of approximately 4 km/h, which is the speed limit of the apparatus. Such a low speed induces a longer contact between the wheel and the ground. As a consequence, the normal stress distribution is greater in the soil samples than it would have been in the field [41], but the frictional stress is lower. This apparatus is generally used to study the rutting of bituminous materials with a smooth-tire. The present tests were devoted to the study of the degradation of soil with sculpted tires. This is the novelty and the originality of the present study. Two tires (Mitas FL-08 and Continental IC 10, named M and C, respectively) were used during the tests (Figure 3c and d). They were inflated to a pressure of 600 kPa as specified by the manufacturer. The contact action of the tire on the ground was achieved by means of an actuator located under the sample. The pressure applied by the actuator was 320 kPa (minimum value of the device), which corresponds to a tire/ground contact pressure of 898 kPa for the Mitas FL-08 tire and 753 kPa for the Continental IC10 tire. These contact pressures did not consider the penetration of the tire into the ground. They were assessed by dividing the contact force by the contact area with a plate. The simulated contact pressures were between 2 and 4 times higher than the non-penetrating contact pressures of Michelin XADN tires, which are used on many articulated dump trucks and varying between 200 and 400 kPa for conventional loads [42]. The tests carried out reproduced, therefore, the movement of rubber tires on the soil with three main scale effects: high stresses, low driving speeds and no engine torque.
Samples I-V were subjected to a total number of 20,000 wheel passes while sample IV underwent 10,000 passes (Table 1). These numbers of passes were chosen to simulate the similar soil stress that can be expected on a construction site over a long period of time (one month approximately). This was defined according to Midwest Research Institute measurements at four earthmoving sites in the United States [43], assuming between 100 and 200 truck passes per day (i.e., between 2000 and 4000 passes per month). Considering three wheel passes per truck (an articulated dump truck has three wheels on each side), it leads to about 10 000 wheel passes per month. The Midwest Research Institute's estimates were carried out on construction sites during periods of low activity [43]. It is assumed that the soil on a construction site during a period of high activity can be subjected to at least twice as much truck traffic (i.e., 20,000 wheel passes).
For each sample, the traffic simulator was stopped after 20, 60, 120, 200, 300, 400, 500, 600, 800, 1000, 2000, 4000 and 10,000 wheel passes. At each stop of the device, the particles detached from the soil were collected by lightly sweeping the surface with a fine brush and subsequently weighed. The evolution of the mass of the removed particles allows the quantification of soil degradation by traffic cycles. The particles detached after 20,000 wheel passes for samples I and II were placed in a LASER diffractometer to measure their PSD. Similarly, the particle sizes after 2000, 10,000 and 20,000 passes on sample III were also measured. The corresponding PSD was measured according to ISO 13,320 [27] and are presented in Section 4.1.2.

Soil-atmosphere Interaction: Wind Tunnel Experiments
After compaction and traffic degradation, samples I, III, IV and V were placed in a wind tunnel and the airflow above the soil was investigated. For this purpose, the boundary layer was characterized. The tests were carried out in the subsonic wind tunnel at ESTACA West Campus (Figure 4). The test section had a squared cross-section of 0.3 × 0.3 m² and was 1 m long. The airflow was generated by a 3kW engine. The mean turbulence intensity in the test section outside the boundary layer was low (<1%). Before wind tunnel experiments, the soil samples were sprayed with a lacquer (Struers TM ) to prevent the boundary layer from being disrupted by dust emissions. Roughness measurements with a rotating LASER profilometer before and after lacquer application showed that this did not alter the surface roughness.
Velocity measurements were recorded using a 2D LASER velocimeter provided by DANTEC® Dynamics. This system was mounted on a 3D displacement table. Two upstream flow velocities were considered (U∞ = 8 m.s −1 et U∞ = 16 m.s −1 ). They were defined to fit with typical truck speeds on construction sites (between 30 and 60 km.h −1 ). The coordinate system was such that x corresponded to the horizontal axis (positive downstream), y and z being the vertical (positive upward) and spanwise (positive from the left to the right looking from the entrance of the test section) directions, respectively. The origin O was taken at the center of the leading edge of the sample, on the channel centreline. The vertical profiles of the velocity were measured at the centreline (z = 0) of the test section at four given positions depicted by points 1 to 4 in Figure 4. The first point was used as a reference on the PVC floor of the wind tunnel (x = −0.03 m) while the three other points were located on the soil sample (x = 0.05, 0.15 and 0.27 m). In Figure 4, the flow is from the left to the right.  (Table 2). From a theoretical point of view, the boundary layer profile developing above the soil can be approximated using Equation (12) where is the Van Kármán constant ( ≈ 0.41) et B a constant determined by experiments. The wind friction velocity u* was determined using the method developed by Djenedi et al. [44]. It can be applied regardless of the surface roughness and the Reynolds number of the flow. By a trialand-error method, the value of u* was found by matching each vertical velocity profile with the normalized velocity defect form given by Equation (13): where δ is the boundary layer thickness (m) and h is a universal function that applies to all turbulent boundary layers [44] and defined by Equation (14): where the coefficients pi (I = 1 The determination of the flow regime was carried out using the shape factor H Equation (15): where δ* is the displacement thickness (m) and θ is the momentum thickness (m), given by Equation (16) and Equation (17), respectively: * = (1 − ) The flow is turbulent if H < 1.4 and laminar if H > 2.6 [45]. From flow velocity measurements, the statistical distributions of their corresponding fluctuations were determined. The statistical distribution of Reynolds shear stresses were determined according to Equation (11). Then, for each velocity profile, the fluctuations closest to the soil were considered to apply the dust emission model (Equation 6). This process is detailed in Section 4.3.1.

Detachment of Particles from the Soil Surface
Surface degradation was quantified using the D parameter which corresponds to the cumulative mass of detached particles in relation to the ground surface covered by the tire during its movement (0.04 m² for each tire in the present study). Figure 5 shows the evolution of D as a function of the number of wheel passes for each of the samples.
where N is the number of wheel passes and ci (I = 1, 2, …, 5) are coefficients determined from the experiments and detailed in Table 3. Table 3. Parameters of the degradation model (Equation (18)).  Figure 5 shows a good agreement between Equation (18) and the experimental data. For the S75K25 soil (samples III, IV and V), the fitting curve was plotted considering an average of the three samples. Indeed, for this soil, a good reproducibility of the experimental data on samples was obtained with the same test conditions (samples III and IV) as well as for sample V which was degraded with a different tire ( Figure 5).

Soil
Overall, the behavior was the same for the three soils: a regular increase of the surface degradation with the number of wheel passes was observed, up to a certain threshold (about 200 passes) where an inflection appeared. During the first traffic cycles, the compacted soil surface seems to resist the shear generated by the tire, but gradually the tire tread pattern starts appearing on the surface. Visually, after about ten passes, the first aggregates segregated and crumbled, leading to a significant increase in soil degradation. After a few hundred cycles, the soil mainly became plasticized: fewer particles were pulled out as rutting (structural deformation) became the main mode of degradation.
The difference between the three soils relies on the quantification of degradation: after 100 wheel passes, degradation was worth several tens, several hundreds and several thousands of grams of particles per unit area for S50K50, S75K25 and S0K100 respectively. S50K50 soil had a well-graded PSD with a good balance between sand and clay. This gave a microstructure to the S50K50 sample that was not as porous as the S75K25 sample [28]. The imbalance between the amount of sand and clay in the S75K25 soil resulted in the presence of more macropores in this material. Indeed, sand particles give rise to large pores between them and S75K25 did not contain enough clay to fill those pores. Macropores are areas of weakness in the granular matrix. Their presence increases the risk of aggregate detachment. The S0K100 soil had a PSD consisting only of clay particles with low porosity. However, this soil had a shear strength which was almost 20% lower than S50K50 soil and 30% lower than S75K25 soil [38]. Thus, the S0K100 soil surface was much less resistant to the stresses generated by wheel traffic. A field test would be required to validate the degradations observed in the laboratory tests.
After studying the appearance of particles on the soil surface during traffic, the PSD of these particles is analyzed in the following section. Figure 6 shows the initial PSD of the three soils and the PSD of the loose particles at 20,000 wheel passes on the S0K100 (Figure 6a) and S50K50 (Figure 6b) samples. Figure 6c shows the same curves, as well as the PSD of the loose particles at 2000 and 10,000 wheel passes on the S75K25 sample. Overall, the detached particles had a lower percentage of fine particles (<10 µm) and a higher percentage of coarse particles (>10 µm) compared to the original soil. This phenomenon was particularly marked for both soils containing sand compared to clay. This was explained by the fact that the sand and clay particles agglomerate to form coarser particles (aggregates) that were torn off by the tire. Figure 6c shows that this phenomenon of particle aggregation decreased as the number of wheel passes increased. There was, therefore, a progressive crumbling of the aggregates with the friction generated by the tire. The gap between the PSD curves that correspond to the degraded state and to the initial state was assessed using the empirical function "difference" (Equation 19). This function was obtained by considering the soil properties and the number of wheel passes. For more details regarding the methodology adopted, refer to the Appendix provided at the end of the present paper.

PSD of Particles Segregated from the Soils during Traffic Degradation
where %clay is the mass percentage of particles smaller than 2 µm and %sand is the mass percentage of particles between 20 µm and 2000 µm [46], dmin and dmax represent the minimum and maximum soil particle diameter (µm) respectively. If the PSD of the initial soil p(d)initial soil is known, it is possible to predict the granulometry of the particles detached after N wheel passes, that is p(d)N passes (Equation 20): (20) Figure 7 compares the data provided in Figure 6 with the prediction results (dashed lines). The correlation between the model and experimental results was examined by studying the difference between the area under the model curve (p(d)model) and the area under the experimental curve (p(d)exp), according to Equation (21): Except for 20 µm < d < 60 µm, where the deviation between the theoretical and experimental curves were obvious for S50K50 and S75K25samples, the average deviations were 6% for S0K100 and 30% for S50K50. For S75K25, the average deviations between the model and the experimental data were 15%, 32% and 47% for 20,000 passes, 10,000 passes and 2000 passes, respectively. Overall, the model was in good agreement with the experimental results, which indicates that it was able to reproduce satisfactorily the mechanisms of degradation for the three soils.  Table 4 presents the main characteristics of the boundary layer determined from the velocity measurements (parameters H, u * and δvs calculated from Equations (15), (13) and (2), respectively). Considering H, all profiles corresponded to either a fully turbulent flow or a transitional regime. Apart from a few situations (at positions x = 0.05 m and 0.27 m on samples IV and V), the thicknesses of the viscous sublayer were similar in terms of magnitude for the different samples considering the same inflow velocity. For the sake of simplicity, for each velocity, the mean over all values was chosen to be implemented into the dust emission estimation model. That is, δvs = 300µm for U∞ = 8 m/s and δvs = 135 µm for U∞ = 16 m/s. It is worthwhile to note that a variation of 50% of this parameter in the model led to a variation of the dust flow of less than 2%.

Estimation of Dust Emissions from Studied Soils
Wind tunnel measurements were used to determine the velocity fluctuations above the samples and to derive Reynolds shear stresses generated by the turbulence developing in the vicinity of the surface. Figure 8 shows an example of these estimated Reynolds stresses for sample IV and U∞ = 8 m/s. Figure 8a shows the streamwise turbulence intensity above the sample. Figure 8b shows the PDF of the velocity fluctuations corresponding to the point depicted by a green cross (x = 0.27 m; y = 0.0001 m) in Figure 8a. The velocity fluctuations appeared to follow a normal distribution defined by Equation (22): where µ et σ are the median value and the geometric standard deviation of the distribution. It was then possible to determine the PDF of τ (Figure 8c) using Equation (11). Figure 9 shows the corresponding PDF of τ for the 4 samples and for the position x = 0.27 m. The analyses were performed at this position as it was the farthest from the wind tunnel entrance. The boundary layer was supposed to be the most developed. In addition, according to Table 4, all velocity profiles at this location were considered as turbulent. The results depicted in Figure 9 appear to be relatively similar from one sample to another. Then, the mean values were considered and plotted (dotted line in Figure 9). Based on these average values, p(f) was determined and implemented in the model given by Equation (6).   Figure 10 shows that estimated dust emissions for clay (S0K100) were 2 to 3 orders of magnitude larger than those for sandy soils (S50K50 and S75K25). This was due to the high degradation capacity of clay. Moreover, this soil contained a significant proportion of PM10. For the three soils, the turbulence generated for the largest velocity (U∞ = 16 m/s) led to a dust flux about 30 times greater than for the smallest velocity (U∞ = 8 m/s). This value is found to be high when compared with the results from Etyemezian et al. [47] who carried out in-situ measurements of dust concentrations related to vehicle traffic on unpaved roads. Their work showed that dust emissions were approximately correlated with vehicle speed to the power 3 and therefore a doubling of the speed results in an eightfold increase in emissions. This shows that turbulence is not the main contributor to dust emissions when a vehicle is in motion. The shear between the tires and the ground must have a major influence, which is not considered in the tests presented herein.

Comparison with Field Data
Gillies et al. [11] carried out measurements for the assessment of dust emissions generated by vehicle traffic on an unpaved road. They established emission factors (EF), expressed in kilograms of PM10 emitted per kilometer travelled by the vehicle (kg.vkt −1 , kilogram per vehicle kilometer travelled), a unit used in many other studies [11,13,14,47]. The results of these measurements were related to US Army truck circulation and are presented in Table 5. Emission factors per unit area (EF in kg.m −2 ) are presented for vehicle speeds of 30 km.h −1 (≈ 8 m.s −1 ) and 60 km.h −1 (≈ 16 m.s −1 ) in order to make a comparison with the present experimental results. Typically, the time required for the dust plume to dissipate in the wake of a truck is a few tens of seconds [48,49]. For the emission factors in Table 5, this corresponds to dust emission fluxes between 0.0001 and 0.001 kg.m −2 .s −1 . These fluxes are in the same orders of magnitude as those estimated for S50K50 and S75K25 soils at 1000 and 10,000 wheel passes (see Figure 10). Nevertheless, the measurements of Gillies et al. [11] were carried out on a soil with a silt content (particle diameter < 75 µm) less than 7%, whereas S50K50 and S75K25 soils had silt contents of 44% and 73%, respectively. Thus, it is reasonable to think that vehicle traffic on these soils would have led to higher dust emissions than those estimated in Figure 10. Indeed, on this figure, dust emissions were generated by turbulence corresponding to an airflow having a speed of about 30 km.h −1 and 60 km.h −1 . When a truck moves at these speeds, turbulence is generated in the wake of the vehicle, which was not being considered during the wind tunnel experiments.

Conclusions and Perspectives
In this paper, the degradation of compacted soils submitted to traffic was studied. Three soils were compacted by roller compaction in the laboratory: a pure clay (S0K100) and two sand-clay mixtures with a clay content of 50% and 25% (S50K50 and S75K25). The soils were then stressed using a traffic simulator. Models were established to quantify the surface degradation as well as the evolution of the PSD of the detached particles as a function of the number of wheel passes. Due to the low wheel speed in the traffic simulator (4 km/h), field measurements would be required to validate the models established in the laboratory. Moreover, the study is based on mixtures of pure clay and sand. This is a first step in the development of soil degradation models. Future work will have to be carried out with soil samples collected from construction sites.
The soil/atmosphere interaction was studied in a wind tunnel. Inflow velocities of 8 m.s −1 and 16 m.s −1 were considered above the traffic degraded soils. The study was focused on the characterization of the turbulence intensities and Reynolds shear stresses generated in the near ground turbulent flow.
The obtained results of soil degradation and flow turbulence were implemented in the CTDE model to estimate dust emissions on the three soils. On the one way, the results pointed out that emissions for clay were several orders of magnitude higher than those for sandy soils. On the other way, turbulence developing for an incoming flow velocity of 16 m.s −1 led to dust emissions 30 times higher than those estimated at 8 m.s −1 .
The results derived from the model were compared with dust emissions from truck traffic on unpaved roads at similar speeds. The comparisons highlight the importance of taking into account the turbulence generated in the wake of trucks. This point will need to be addressed in future works. Similarly, environmental factors and shear generated by the tires are needed to be taken into account for a better estimate of emissions. The present study is a first step toward the development of an analytical model of dust emission generated by truck traffic. The use of this model should leadto a better management of the water resources necessary for dust abatement on earthmoving worksites.

Appendix: Model for the determination of the PSD of particles detached from a soil by traffic
From the Figure 6, the differences between the PSD of the detached particles and that of the initial soil were plotted for each soil sample. The solid curves in Figure A1 were obtained. Figure A1. Deviation from the initial soil particle size distribution at 20,000 wheel passes on the sample (a) S0K100, (b) S50K50 and (c) S75K25. For the latter sample, the deviations at 10,000 passes (d) and 2,000 passes (e) are also presented. The dashed curves represent the model of Equation (31).
For each sample, a sinusoidal shape with a variable amplitude is observed. This evolution was approximated using a function according to the following steps: 1. A function was constructed allowing the description of the sinusoidal variation. This function looked like: where A is a function depending on the type of soil and the number of wheel passes.
2. The scales in Figure A1 were semi-logarithmic, so the function took the following form:

= × ln ( ) × sin [ln ( )]
3. The curves in Figure A1 were defined between dmin et dmax and can be approximated as having a period equal to dmax-dmin, which gave the following function: 4. On the definition domain, the curves were first negative and then positive. This was the inverse of the behaviour of the function given by Equation (25). A final modification was therefore necessary: 5. It was therefore a question of determining the amplitude function A which was written as: with fsoil a function depending on the type of soil and a function depending on the number of wheel passes N.
6. The function was determined from the curves in Figures A1-c, A1-d and A1-e. Assuming that fsoi l= 1, the value of corresponded to the amplitude A of Equation (26). Then, was determined by trial-and-error in order to minimize the average weighted deviation between experimental and theoretical curves. The weighted deviation was defined by: A weighted deviation criterion was chosen in order to minimize the difference between experimental and theoretical curves for the amplitude peaks, which were the most important parts of the curve to model. According to that, it appeared that the best approximation for the function was a second-degree function: ( ) = 9,72.10 − 4,92.10 + 1,09 7. The fsoil function was determined using the curves in Figures A1-a, A1-b et A1-c. These three curves corresponded to the case where g (N = 20,000) = 0.49. The amplitude of the difference between 20,000 passes and initial soil was low for S0K100 (clay), medium for S75K25 and high for S50K50. Thus, it was considered that this amplitude depended on the product of the percentage of clay by the percentage of sand in the soil (%clay × %sand). According to this definition and based on the particle size distributions of the three soils, S0K100 had 20% clay and 7% sand ( %clay × %sand = 140 ), S50K50 had 13% clay and 36,6% sand (% × % = 475,8) and S75K25 had 6,7% clay and 64% sand ( % × % = 428,8). The amplitude function of Equation (26) to approximate the curves in Figures A1-a As with the estimation of the function, the fsoil function was assessed by trial-and-error to minimize the weighted deviation (Equation (28)). The best approximation for fsoil was: Finally, the function for modelling the difference between the particle size distribution of the loose particles and the initial soil was defined by: The comparisons between the model estimates and the experimental results were shown in Figure A1. Relatively good agreements were found.