# Polarization Lidar Detection of Agricultural Aerosol Emissions

Authors: Eduard Gregorio; Jordi Gené; Ricardo Sanz; Francesc Rocadenbosch; Patricia Chueca; Jaume Arnó; Francesc Solanelles; Joan R. Rosell-Polo
Journal: Journal of Sensors
 (2018)

---

## Abstract

Agricultural aerosol emissions can significantly impact human and animal health as well as the environment. Therefore, it is essential to adopt new sensing techniques for real-time monitoring these emissions in high temporal and spatial resolution. In recent years, light detection and ranging (lidar) technology has been used for measuring the particulate matter emitted from agricultural operations. However, conventional nonpolarized lidar systems cannot discriminate between different types of aerosols, which can lead to misinterpretation of the results. To overcome this limitation, this study applies the polarization lidar technique to monitor agricultural aerosols. A 355 nm polarization lidar system was used to measure the emissions generated during pesticide spraying operations. The results showed that depolarization ratios due to field dust (0.220–0.268) and to road dust (0.385) are clearly higher than those caused by pesticide spray drift (0.028–0.043) or by diesel exhaust (0.099), which can be used to differentiate each type of aerosol. These results support the development of new polarization lidar systems specifically designed to study the impact of agricultural activities on air quality.


---

## Body

## 1. Introduction

Agricultural and livestock farming generates considerable aerosol emissions, including particulate matter (PM) and pesticides. Agricultural activities are responsible for 5% of all particles less than 2.5μm in diameter (PM2.5) and 25% of all particles less than 10 μm (PM10) [1]. The main agricultural activities responsible for PM emissions are tilling, harvesting, agricultural fires, and emissions from fertilizer use and livestock [2]. Pesticide spray drift is another large source of pollution. Spray drift losses can be as great as 30–50% of the applied product [3]. These aerosols have a negative effect on the environment and air quality, entailing serious risks to human and animal health. For example, PM has been linked to cardiac and respiratory diseases [4]. In addition, spray drift can damage crops close to the treated area, contaminate surface water, and reach residential areas. It has been demonstrated that pesticides can travel thousands of kilometers and have been found in areas as remote as the polar regions [5].Monitoring these aerosols is extremely important due to the need to identify their sources, evaluate their effects, and propose appropriate mitigation measures. Point samplers are commonly used but unable to provide a full view of the emitted plume, and they are labor-intensive and time-consuming. Light detection and ranging (lidar) is an optical remote sensing technique that can overcome these limitations. Lidar systems can provide real-time aerosol measurements with high temporal and spatial resolution. In previous studies, nonpolarized lidar systems have been used to measure PM10 emissions from land preparation operations [6–9] to determine vertical PM concentration profiles. Several lidar studies have also been carried out to estimate the PM emissions due to livestock housing. For example, PM concentration and fluxes generated by a swine-feeding facility have been measured using a three-wavelength lidar system in combination with mass samplers and optical particle counters [10]. An alternative methodology has been recently developed to estimate the PM emission rates from animal production facilities departing from single-wavelength lidar data [11]. In another study [12], a scanning lidar system was used to obtain the backscatter coefficients of aerosols emitted by a pig farm. Lidar technology has also been applied in spray drift studies to assess the dynamics of aerially released pesticides [13–15], to study the influence of atmospheric stability over spray dispersal [16, 17], and to quantify the spray plume concentration [18]. High correlations between lidar and collector measurements of spray drift have been obtained in recent studies [19–21]. Lidar systems are also a suitable tool to validate spray-transport models [22] and dust plume dispersal models of agriculture field operations [23, 24].The development of new eye-safe scanning lidar systems specifically designed for agricultural applications [25, 26] is expected to contribute to agricultural air quality monitoring. Nevertheless, conventional lidar systems are unable to distinguish between different types of aerosols, which can lead to misinterpretation of the measurements. For example, in a terrestrial spray drift study, the return signal captured by a nonpolarized lidar system may be assumed to be due solely to the drift plume generated by the sprayer. However, the return signal could be contaminated by dust or other unexpected aerosols. Unlike conventional systems, polarization lidars can differentiate aerosols by their shape [27]. Polarization lidar systems are used to identify the presence of ice crystals in mixed-phase clouds [28], for the detection of desert dust in the troposphere [29], to study the polar stratospheric clouds [30], and to measure the contrails produced by aircraft [31], among other applications. In this work, a polarization lidar system is used to measure air emissions generated during spraying operations with an air-assisted sprayer. Although polarization lidar systems are a well-established tool for atmospheric research, to the authors’ knowledge, they have not been previously applied for the monitoring and identification of agricultural emissions. The aim of this work is to determine if the polarization lidar technique is a suitable method for unambiguously discriminating between different types of agricultural aerosol emissions.


## 2. Materials and Methods

### 2.1. Experimental Setup

The trials were conducted in two fields owned by the Institut de Recerca i Tecnologia Agroalimentàries in Gimenells (lat. 41°39′11″N, long. 0°23′28″E, elev. 259 m) located about 25 km from Lleida, Spain. Figure 1 shows the experimental setup at both study sites. Field no. 1 was an apple orchard, while field no. 2 was a flat fallow field. An air-assisted sprayer (Ilemo Arrow F-1000, Ilemo/Hardi SA, Lleida, Spain) and three nozzle types, a hollow cone (Albuz ATR Orange, Saint-Gobain, Evreux, France), an air injected low drift (Albuz TVI 80 02 Yellow), and a disc-core full cone (Teejet D3DC35, Spraying Systems Co., Wheaton, Illinois, USA), were used. Aerosol measurements were carried out using a polarization lidar system (ALS 300, Leosphere, Orsay, France). This system transmits 16 mJ laser pulses at 355 nm (tripled Nd:YAG) with a pulse repetition frequency of 20 Hz [32]. The temporal and spatial resolutions were adjusted to 1 s and 1.5 m, respectively. The range of full overlap, i.e., the distance from which the full laser beam is seen within the receiver field of view, was adjusted to 80 m. In all the trials, the lidar system was located farther than 200 m from the sprayer (Figure 1). This distance ensured that the aerosol was always measured within the full overlap area.Figure 1
Experimental setup with sprayer and lidar location: (a) field no. 1: trials T1 to T5; (b) field no. 2: trials T6 to T10.
(a)
(b)The goal of the trials conducted in field no. 1 (trials T1 to T5) was to measure the different types of aerosols emitted during terrestrial spraying operations. In trial T1, the sprayer was kept in a static position and located in alley 1 (Figure1(a)). In this trial, the tractor motor was turned on for 30 s without spraying and with the fan off, thus only diesel exhaust was emitted. Trial T2 was carried out in alley 2, and in this case, the sprayer travelled a total distance of 50 m along the alley within 30 s (dynamic trial). This trial was carried out without spraying but with the fan on, raising field dust. Trial T3 was similar to T2, but it was conducted in alley 3. T4 and T5 were dynamic trials in which spraying was carried out. Trial T4 was performed in alley 4, while trial T5 was conducted in alley 5.In trials T1 to T4, the speed of the power take-off (PTO) was 370 rpm, while in trial T5 it was increased to 540 rpm. The lidar system was kept in a static position during all trials. The laser beam was pointed horizontally and aligned parallel to the first row of trees in the apple orchard. The separation between this row and the laser beam was 5 m. The purpose of this setup was to measure the emitted aerosol that was able to travel past the treated crop and reach the neighbouring field. Five spray trials (T6 to T10) with a duration of 30 s were carried out in field no 2. As shown in Figure1(b), the sprayer position was adjusted in each trial to ensure that the spray drift always reached the measurement area.Table1 shows the experimental conditions of all the trials. For each trial, the following information is provided: experimental field, presence or absence of vegetation, alley where the sprayer was placed, operational mode of the sprayer, sprayer fan setting, spraying application, nozzle model employed, number of open nozzles, and individual flow rate.Table 1
Experimental field trials.
Trial
Location
Vegetation (alley)
Operational mode
Sprayer fan
Spraying
Nozzles
Model
Number
Flow rate (l min−1 nozzle)1
T1
Field no. 1
Yes (1)
Static
Off
No
—
—
—
T2
Field no. 1
Yes (2)
Dynamic
On
No
—
—
—
T3
Field no. 1
Yes (3)
Dynamic
On
No
—
—
—
T4
Field no. 1
Yes (4)
Dynamic
On
Yes
Albuz ATR Orange
16
1.39
T5
Field no. 1
Yes (5)
Dynamic
On
Yes
Albuz ATR Orange
16
1.39
T6
Field no. 2
No
Static
On
Yes
Albuz ATR Orange
10
1.39
T7
Field no. 2
No
Static
On
Yes
Albuz ATR Orange
5
1.39
T8
Field no. 2
No
Static
On
Yes
Teejet D3DC35
5
2.0
T9
Field no. 2
No
Static
On
Yes
Albuz TVI 80 Yellow Orange
10
1.46
T10
Field no. 2
No
Static
On
Yes
Teejet D3DC35
10
2.0
1Individual nozzle flow rate at a pressure of 1 MPa.


### 2.2. Polarization Lidar Methodology

Polarization lidar systems are based on the emission of linearly polarized light and the detection of the return power in the orthogonal and parallel planes of polarization. Range-corrected parallel and orthogonal backscatter lidar signals are given by(1)UR=KβRT2R,U⊥R=K⊥β⊥RT2R,where U and U⊥ are the range-corrected return powers in the parallel and orthogonal planes of polarization, respectively. K and K⊥ are the system constants for the two channels. β and β⊥ are the total parallel- and the total orthogonal-polarized backscatter coefficients, and T2R is the two-way atmospheric transmittance.The linear volume depolarization ratioδv is used to quantify the degree of depolarization of the backscattered light and is defined by
(2)δvR=β⊥RβR.Substituting (1) into (2), the following expressions [33] are obtained for the linear volume depolarization ratio δv and for the range-corrected total backscatter signal U:
(3)δvR=U⊥RK∗UR,(4)UR=UR+U⊥RK∗,where K∗=K⊥/K is the cross-calibration constant. Following [34], the value of K∗ depends on the optics transmissivity of the Brewster plates used to separate both polarizations and is given by
(5)K∗=U⊥RT1UR1−T11−T2+δvR,where T1 and T2 are the transmissions of the parallel polarization contributions of the two Brewster plates, estimated as T1=T2=0.92 (Dr. P. Chazette, private communication).Using vertical profile measurements with our ALS300 system and calibration of the volume depolarization in a molecular reference layer (height interval: 4000–4500 m), a cross-calibration constantK∗=20.9±0.7 was obtained. In detail, this methodology assumes δv=δm, where δm is the molecular depolarization ratio. At a wavelength of 355 nm, a value of δm=0.00395 has been considered [35]. The obtained cross-calibration constant figure is similar to the one estimated by Chazette et al. [34] for the homemade version of the ALS system.WithK∗ known, the polarization system becomes calibrated and the volume depolarization ratio δV can be readily retrieved from (3) and the measurements of U and U⊥.The volume depolarization ratio depends on the type and concentration of the aerosols present in the atmosphere as well as on the molecular contribution. Therefore, to identify different types of aerosols, it is necessary to determine the linear particle depolarization ratioδp as it only depends on particle properties. For example, large values of δp correspond to nonspherical particles or to particles with a nonhomogenous diffraction index.In vertical lidar sounding of the atmosphere, the particle depolarization ratioδp is related to the volume depolarization ratio δv as [36]
(6)δp=1+δmδvRp−1+δvδm1+δmRp−1+δv,where Rp is the backscatter ratio, given by
(7)Rp=βm+βpβm,where βm and βp are the molecular and particle backscatter coefficients, respectively.Whereas in ground-based vertical lidar sounding, subscripts “p” and “m” refer to the particle and molecular components, respectively, in the horizontal lidar sounding addressed in this work, subscripts “b” and “p” refer to the background component and to the agricultural aerosol component. The background component comprises both aerosols and molecules along the horizontal path when the tractor and the sprayer are off. Equations (6) and (7) are rewritten and adapted according to such analogy of variables:
(8)δp=1+δbδvRp−1+δvδb1+δbRp−1+δv,(9)Rp=βb+βpβb,where δb is the volume background depolarization ratio and βb is the background backscatter coefficient. δb is obtained by computing the volume depolarization ratio (3) given tractor-off measured range-corrected orthogonal and parallel power profiles, U⊥ and U, respectively, along the horizontal path, and calibration constant K∗. The volume depolarization ratio δv is measured with the tractor on. The backscatter ratio Rp is computed as the ratio of the total (background + agricultural aerosol), Utot, to the background range-corrected power, Ub,
(10)RpR≈UtotRUbR,R∈a,b,where
(11)UtotR=βb+βpRTb2RTp2R,UbR=βbTb2R,R∈a,b,where R is the range along the horizontal path and a and b are the initial and final ranges of the aerosol cloud, 0<a<R<b<Rmax, with Rmax being the maximum range. Tb2R and Tp2R are the optical transmittances due to background and particle (agricultural aerosol) extinctions, respectively. Equation (10) approximation essentially assumes that the two-way particle transmission is close to unity, Tp2R≈1, or, equivalently, that the optical depth in a,b is much smaller than unity, τP=∫abαPudu≪1, where αp is the particle extinction coefficient.


### 2.3. Lidar Data Analysis

The received lidar data was background-subtracted to remove the residual skylight and offset (noise) due to the optoelectronic receiver. Then, background-subtracted signals from polarization channels (parallel and orthogonal) were range-corrected, calibrated, and summed to obtain the total backscatter signal (4). This total signal is used to generate range-time intensity (RTI) plots of the agricultural aerosol emissions. The range-corrected background-subtracted signals from each channel were also used to calculate the volume depolarization ratio δv by applying (3). The next step is to compute the particle depolarization ratio δp (depolarization ratio of the agricultural aerosols) by applying (8).Finally, a basic error propagation study has been carried out so as to assess the impact of absolute errorsΔK∗ and ΔRp which are considered the two dominating error sources in this study, on the estimated particle depolarization ratio (8). An extensive error treatment of the depolarization lidar problem has been presented in [36] and is out of the scope of the present work. At this point, it is worth mentioning that in (8), both δv and δb are functions of K∗ via (3), which is the constitutive equation defining the volume depolarization ratio. Therefore, the absolute error on the volume depolarization can be written as
(12)Δδi=ΔK∗K∗2U⊥Ui,where subscript i=b,v stands for “background” or “volume” measurement, respectively.When (12) is combined with the basic error propagation of (8),
(13)Δδp=ΔδpΔδbΔδb+ΔδpΔδvΔδv+ΔδpΔRpΔRp,the sought-after absolute error Δδp in response to absolute errors ΔK∗ and ΔRp is obtained.


## 2.1. Experimental Setup

The trials were conducted in two fields owned by the Institut de Recerca i Tecnologia Agroalimentàries in Gimenells (lat. 41°39′11″N, long. 0°23′28″E, elev. 259 m) located about 25 km from Lleida, Spain. Figure 1 shows the experimental setup at both study sites. Field no. 1 was an apple orchard, while field no. 2 was a flat fallow field. An air-assisted sprayer (Ilemo Arrow F-1000, Ilemo/Hardi SA, Lleida, Spain) and three nozzle types, a hollow cone (Albuz ATR Orange, Saint-Gobain, Evreux, France), an air injected low drift (Albuz TVI 80 02 Yellow), and a disc-core full cone (Teejet D3DC35, Spraying Systems Co., Wheaton, Illinois, USA), were used. Aerosol measurements were carried out using a polarization lidar system (ALS 300, Leosphere, Orsay, France). This system transmits 16 mJ laser pulses at 355 nm (tripled Nd:YAG) with a pulse repetition frequency of 20 Hz [32]. The temporal and spatial resolutions were adjusted to 1 s and 1.5 m, respectively. The range of full overlap, i.e., the distance from which the full laser beam is seen within the receiver field of view, was adjusted to 80 m. In all the trials, the lidar system was located farther than 200 m from the sprayer (Figure 1). This distance ensured that the aerosol was always measured within the full overlap area.Figure 1
Experimental setup with sprayer and lidar location: (a) field no. 1: trials T1 to T5; (b) field no. 2: trials T6 to T10.
(a)
(b)The goal of the trials conducted in field no. 1 (trials T1 to T5) was to measure the different types of aerosols emitted during terrestrial spraying operations. In trial T1, the sprayer was kept in a static position and located in alley 1 (Figure1(a)). In this trial, the tractor motor was turned on for 30 s without spraying and with the fan off, thus only diesel exhaust was emitted. Trial T2 was carried out in alley 2, and in this case, the sprayer travelled a total distance of 50 m along the alley within 30 s (dynamic trial). This trial was carried out without spraying but with the fan on, raising field dust. Trial T3 was similar to T2, but it was conducted in alley 3. T4 and T5 were dynamic trials in which spraying was carried out. Trial T4 was performed in alley 4, while trial T5 was conducted in alley 5.In trials T1 to T4, the speed of the power take-off (PTO) was 370 rpm, while in trial T5 it was increased to 540 rpm. The lidar system was kept in a static position during all trials. The laser beam was pointed horizontally and aligned parallel to the first row of trees in the apple orchard. The separation between this row and the laser beam was 5 m. The purpose of this setup was to measure the emitted aerosol that was able to travel past the treated crop and reach the neighbouring field. Five spray trials (T6 to T10) with a duration of 30 s were carried out in field no 2. As shown in Figure1(b), the sprayer position was adjusted in each trial to ensure that the spray drift always reached the measurement area.Table1 shows the experimental conditions of all the trials. For each trial, the following information is provided: experimental field, presence or absence of vegetation, alley where the sprayer was placed, operational mode of the sprayer, sprayer fan setting, spraying application, nozzle model employed, number of open nozzles, and individual flow rate.Table 1
Experimental field trials.
Trial
Location
Vegetation (alley)
Operational mode
Sprayer fan
Spraying
Nozzles
Model
Number
Flow rate (l min−1 nozzle)1
T1
Field no. 1
Yes (1)
Static
Off
No
—
—
—
T2
Field no. 1
Yes (2)
Dynamic
On
No
—
—
—
T3
Field no. 1
Yes (3)
Dynamic
On
No
—
—
—
T4
Field no. 1
Yes (4)
Dynamic
On
Yes
Albuz ATR Orange
16
1.39
T5
Field no. 1
Yes (5)
Dynamic
On
Yes
Albuz ATR Orange
16
1.39
T6
Field no. 2
No
Static
On
Yes
Albuz ATR Orange
10
1.39
T7
Field no. 2
No
Static
On
Yes
Albuz ATR Orange
5
1.39
T8
Field no. 2
No
Static
On
Yes
Teejet D3DC35
5
2.0
T9
Field no. 2
No
Static
On
Yes
Albuz TVI 80 Yellow Orange
10
1.46
T10
Field no. 2
No
Static
On
Yes
Teejet D3DC35
10
2.0
1Individual nozzle flow rate at a pressure of 1 MPa.


## 2.2. Polarization Lidar Methodology

Polarization lidar systems are based on the emission of linearly polarized light and the detection of the return power in the orthogonal and parallel planes of polarization. Range-corrected parallel and orthogonal backscatter lidar signals are given by(1)UR=KβRT2R,U⊥R=K⊥β⊥RT2R,where U and U⊥ are the range-corrected return powers in the parallel and orthogonal planes of polarization, respectively. K and K⊥ are the system constants for the two channels. β and β⊥ are the total parallel- and the total orthogonal-polarized backscatter coefficients, and T2R is the two-way atmospheric transmittance.The linear volume depolarization ratioδv is used to quantify the degree of depolarization of the backscattered light and is defined by
(2)δvR=β⊥RβR.Substituting (1) into (2), the following expressions [33] are obtained for the linear volume depolarization ratio δv and for the range-corrected total backscatter signal U:
(3)δvR=U⊥RK∗UR,(4)UR=UR+U⊥RK∗,where K∗=K⊥/K is the cross-calibration constant. Following [34], the value of K∗ depends on the optics transmissivity of the Brewster plates used to separate both polarizations and is given by
(5)K∗=U⊥RT1UR1−T11−T2+δvR,where T1 and T2 are the transmissions of the parallel polarization contributions of the two Brewster plates, estimated as T1=T2=0.92 (Dr. P. Chazette, private communication).Using vertical profile measurements with our ALS300 system and calibration of the volume depolarization in a molecular reference layer (height interval: 4000–4500 m), a cross-calibration constantK∗=20.9±0.7 was obtained. In detail, this methodology assumes δv=δm, where δm is the molecular depolarization ratio. At a wavelength of 355 nm, a value of δm=0.00395 has been considered [35]. The obtained cross-calibration constant figure is similar to the one estimated by Chazette et al. [34] for the homemade version of the ALS system.WithK∗ known, the polarization system becomes calibrated and the volume depolarization ratio δV can be readily retrieved from (3) and the measurements of U and U⊥.The volume depolarization ratio depends on the type and concentration of the aerosols present in the atmosphere as well as on the molecular contribution. Therefore, to identify different types of aerosols, it is necessary to determine the linear particle depolarization ratioδp as it only depends on particle properties. For example, large values of δp correspond to nonspherical particles or to particles with a nonhomogenous diffraction index.In vertical lidar sounding of the atmosphere, the particle depolarization ratioδp is related to the volume depolarization ratio δv as [36]
(6)δp=1+δmδvRp−1+δvδm1+δmRp−1+δv,where Rp is the backscatter ratio, given by
(7)Rp=βm+βpβm,where βm and βp are the molecular and particle backscatter coefficients, respectively.Whereas in ground-based vertical lidar sounding, subscripts “p” and “m” refer to the particle and molecular components, respectively, in the horizontal lidar sounding addressed in this work, subscripts “b” and “p” refer to the background component and to the agricultural aerosol component. The background component comprises both aerosols and molecules along the horizontal path when the tractor and the sprayer are off. Equations (6) and (7) are rewritten and adapted according to such analogy of variables:
(8)δp=1+δbδvRp−1+δvδb1+δbRp−1+δv,(9)Rp=βb+βpβb,where δb is the volume background depolarization ratio and βb is the background backscatter coefficient. δb is obtained by computing the volume depolarization ratio (3) given tractor-off measured range-corrected orthogonal and parallel power profiles, U⊥ and U, respectively, along the horizontal path, and calibration constant K∗. The volume depolarization ratio δv is measured with the tractor on. The backscatter ratio Rp is computed as the ratio of the total (background + agricultural aerosol), Utot, to the background range-corrected power, Ub,
(10)RpR≈UtotRUbR,R∈a,b,where
(11)UtotR=βb+βpRTb2RTp2R,UbR=βbTb2R,R∈a,b,where R is the range along the horizontal path and a and b are the initial and final ranges of the aerosol cloud, 0<a<R<b<Rmax, with Rmax being the maximum range. Tb2R and Tp2R are the optical transmittances due to background and particle (agricultural aerosol) extinctions, respectively. Equation (10) approximation essentially assumes that the two-way particle transmission is close to unity, Tp2R≈1, or, equivalently, that the optical depth in a,b is much smaller than unity, τP=∫abαPudu≪1, where αp is the particle extinction coefficient.


## 2.3. Lidar Data Analysis

The received lidar data was background-subtracted to remove the residual skylight and offset (noise) due to the optoelectronic receiver. Then, background-subtracted signals from polarization channels (parallel and orthogonal) were range-corrected, calibrated, and summed to obtain the total backscatter signal (4). This total signal is used to generate range-time intensity (RTI) plots of the agricultural aerosol emissions. The range-corrected background-subtracted signals from each channel were also used to calculate the volume depolarization ratio δv by applying (3). The next step is to compute the particle depolarization ratio δp (depolarization ratio of the agricultural aerosols) by applying (8).Finally, a basic error propagation study has been carried out so as to assess the impact of absolute errorsΔK∗ and ΔRp which are considered the two dominating error sources in this study, on the estimated particle depolarization ratio (8). An extensive error treatment of the depolarization lidar problem has been presented in [36] and is out of the scope of the present work. At this point, it is worth mentioning that in (8), both δv and δb are functions of K∗ via (3), which is the constitutive equation defining the volume depolarization ratio. Therefore, the absolute error on the volume depolarization can be written as
(12)Δδi=ΔK∗K∗2U⊥Ui,where subscript i=b,v stands for “background” or “volume” measurement, respectively.When (12) is combined with the basic error propagation of (8),
(13)Δδp=ΔδpΔδbΔδb+ΔδpΔδvΔδv+ΔδpΔRpΔRp,the sought-after absolute error Δδp in response to absolute errors ΔK∗ and ΔRp is obtained.


## 3. Results and Discussion

### 3.1. Volume Depolarization Ratio Measurements

Figure2 presents the RTI plots of the aerosol clouds corresponding to trials T2, T5, T9, and T10. The left-hand column shows the range-corrected total backscatter power, and the right-hand column shows the volume depolarization ratio per trial. Figure 2(a) presents the dust cloud emitted during trial T2; the cloud approached the lidar position as the sprayer advanced along alley 2. Trial T5, shown in Figure 2(b), included dynamic spraying along alley 5; a mixture of spray drift and field dust was generated. In comparison with trial T2, the cloud had a greater dispersion (up to 70 m laterally) due to the greater distance between the sprayer and the laser beam.Figure 2
(Left) Range-corrected total backscatter power (arbitrary units). (Right) Volume depolarization ratio. (a) Trial T2: dynamic trial without spraying (alley 2, field no. 1). (b) Trial T5: dynamic spraying at 540 rpm (alley 5, field no. 1). (c) Trial T9: static spraying (field no. 2). (d) Trial T10: static spraying (field no. 2).
(a)
(b)
(c)
(d)Trial T9, shown in Figure2(c), was conducted in field no. 2 with the sprayer in a static position. Volume depolarization ratio had lesser values in trial T9 than in trials T2 and T5, because in trial T9 only spray drift was generated, while in trials T2 and T5 field dust was also emitted. Figure 2(d) shows trial T10, which was similar to trial T9 except in the nozzles used (Table 1) and the sprayer position (Figure 1). Interestingly, two clouds, each corresponding to a different type of aerosol, were identified. The first cloud was detected from the beginning of the trial to approximately 80 s and was due to dust ejected into the air by a passing car on a dirt road next to the field. A second cloud was detected from 85 onwards and was due to spray drift. The comparison of both clouds shows that the range-corrected total backscatter power by spray drift is an order of magnitude greater. In contrast, the greatest volume depolarization ratio is caused by the road dust.


### 3.2. Particle Depolarization Ratio Measurements

For each trial, Table2 presents the types of aerosols emitted as well as the mean value and standard deviation of the volume background depolarization ratio δb and the particle depolarization ratio δp of the aerosol as it has been defined by (8).Table 2
Type of aerosol detected in each trial, volume background depolarization ratio (δb), and particle depolarization ratio of the aerosol (δp). For each depolarization value, the mean and the standard deviation are given.
Trial
Type of aerosol
δb
δp
T1
Diesel exhaust
0.018 ± 0.0007
0.099 ± 0.0075
T2
Field dust
0.019 ± 0.0004
0.220 ± 0.0376
T3
Field dust
0.020 ± 0.0011
0.268 ± 0.0402
T4
Mixture of spray drift and field dust
0.026 ± 0.0006
0.218 ± 0.0512
T5
Mixture of spray drift and field dust
0.027 ± 0.0007
0.213 ± 0.0660
T6
Spray drift
0.013 ± 0.0003
0.035 ± 0.0080
T7
Spray drift
0.014 ± 0.0003
0.028 ± 0.0054
T8
Spray drift
0.015 ± 0.0009
0.043 ± 0.0086
T9
Spray drift
0.016 ± 0.0016
0.033 ± 0.0092
T10
Road dust/spray drift
0.025 ± 0.0008
0.385 ± 0.0568/0.034 ± 0.0082From the measurements of the particle depolarization ratio, four types of aerosols have been identified: spray drift, diesel exhaust, field dust, and road dust. Spray driftδp ranges between 0.028 and 0.043. The Lorenz-Mie theory states that spherical particles do not modify the original polarization state of the light at the backscattering direction [37]. The sphericity of the droplets was laboratory tested by means of a phase Doppler analyzer (57X10 PDA, Dantec Dynamics, Skovlunde, Denmark); 3% to 20% of the generated droplets at the nozzle output were not spherical, which explains the resulting depolarizations.Depolarization of tractor engine exhaust emissions has been determined (δp=0.099) in trial T1. To the best of our knowledge, only relative depolarization ratios of diesel exhaust have been measured in a previous study [38]. It should be noted that relative depolarizations of diesel exhaust obtained in the cited study are below the background depolarization, which does not happen in our case. As these authors point out, the composition and the depolarization ratio of exhaust emissions is greatly influenced by the combustion efficiency.Field dust is raised by the tractor as it moves along the alleys.δp of 0.220 and 0.268 have been measured in trials T2 and T3, respectively. These results are in good agreement with atmospheric lidar measurements of Saharan dust (δp=0.24–0.28 at 355 nm [36]). The high depolarization associated with the dust is explained by the irregular shape of crystalline particles that constitute it.Trials T4 and T5 are examples of spray drift tests where the cloud is contaminated by dust. In fact, resulting depolarization ratios (δp≈0.21) demonstrate that dust predominates over spraying in the monitored clouds. Polarization lidar technique opens the possibility of computing the fraction of each type of aerosol in mixture clouds as presented here. For instance, δp is currently used to separate dust from nondust in atmospheric research [39]. Its application in farming requires a laboratory characterization of δp of main agricultural aerosols.Road dust raised by the car (trial T10) presentsδp=0.385, a value similar to that obtained in previous road dust measurements (δp=0.33−0.37 at 355 nm [40]) carried out under controlled conditions. The higher δp of the road dust in comparison to field dust is attributed to a different composition of the road floor, with an important presence of gravel.Figure3 represents, by means of normalized histograms, the distribution functions of the particle depolarization ratio for all the types of aerosols analysed in this study. This graph evidences the ability to distinguish and identify each of the agricultural aerosols through δp.Figure 3
Distribution of particle depolarization ratio for spray drift (T6–T10), diesel exhaust (T1), mixture of spray drift and dust (T4 and T5), field dust (T2 and T3), and road dust (T10).As discussed in Section2.3, the two main sources of error in δp are the cross-calibration constant K∗ and the backscatter ratio Rp. K∗ can be affected by the existence of a small amount of aerosols in the calibration layer [36]. From (5), if δv was 0.005, that is, 20% higher than the molecular depolarization ratio δm, then K∗ would vary by −9%.To computeRp, ideal clouds (with no extinction) have been assumed in (10) and (11). For dust measurements, the uncertainty of Rp can be estimated by assuming an extinction coefficient of 0.15 km−1 at 355 nm [41] and a cloud depth equal to 25 m (typical value for the clouds of this study). This yields a variation of −1% on Rp.The numerical calculation of (12) and (13) for previous uncertainties on K∗ and Rp gives relative errors of about 10–12% and 17–19% on particle depolarization ratios of spray drift and dust, respectively. These are acceptable error figures, similar to those obtained in other lidar polarization experiments [33].


## 3.1. Volume Depolarization Ratio Measurements

Figure2 presents the RTI plots of the aerosol clouds corresponding to trials T2, T5, T9, and T10. The left-hand column shows the range-corrected total backscatter power, and the right-hand column shows the volume depolarization ratio per trial. Figure 2(a) presents the dust cloud emitted during trial T2; the cloud approached the lidar position as the sprayer advanced along alley 2. Trial T5, shown in Figure 2(b), included dynamic spraying along alley 5; a mixture of spray drift and field dust was generated. In comparison with trial T2, the cloud had a greater dispersion (up to 70 m laterally) due to the greater distance between the sprayer and the laser beam.Figure 2
(Left) Range-corrected total backscatter power (arbitrary units). (Right) Volume depolarization ratio. (a) Trial T2: dynamic trial without spraying (alley 2, field no. 1). (b) Trial T5: dynamic spraying at 540 rpm (alley 5, field no. 1). (c) Trial T9: static spraying (field no. 2). (d) Trial T10: static spraying (field no. 2).
(a)
(b)
(c)
(d)Trial T9, shown in Figure2(c), was conducted in field no. 2 with the sprayer in a static position. Volume depolarization ratio had lesser values in trial T9 than in trials T2 and T5, because in trial T9 only spray drift was generated, while in trials T2 and T5 field dust was also emitted. Figure 2(d) shows trial T10, which was similar to trial T9 except in the nozzles used (Table 1) and the sprayer position (Figure 1). Interestingly, two clouds, each corresponding to a different type of aerosol, were identified. The first cloud was detected from the beginning of the trial to approximately 80 s and was due to dust ejected into the air by a passing car on a dirt road next to the field. A second cloud was detected from 85 onwards and was due to spray drift. The comparison of both clouds shows that the range-corrected total backscatter power by spray drift is an order of magnitude greater. In contrast, the greatest volume depolarization ratio is caused by the road dust.


## 3.2. Particle Depolarization Ratio Measurements

For each trial, Table2 presents the types of aerosols emitted as well as the mean value and standard deviation of the volume background depolarization ratio δb and the particle depolarization ratio δp of the aerosol as it has been defined by (8).Table 2
Type of aerosol detected in each trial, volume background depolarization ratio (δb), and particle depolarization ratio of the aerosol (δp). For each depolarization value, the mean and the standard deviation are given.
Trial
Type of aerosol
δb
δp
T1
Diesel exhaust
0.018 ± 0.0007
0.099 ± 0.0075
T2
Field dust
0.019 ± 0.0004
0.220 ± 0.0376
T3
Field dust
0.020 ± 0.0011
0.268 ± 0.0402
T4
Mixture of spray drift and field dust
0.026 ± 0.0006
0.218 ± 0.0512
T5
Mixture of spray drift and field dust
0.027 ± 0.0007
0.213 ± 0.0660
T6
Spray drift
0.013 ± 0.0003
0.035 ± 0.0080
T7
Spray drift
0.014 ± 0.0003
0.028 ± 0.0054
T8
Spray drift
0.015 ± 0.0009
0.043 ± 0.0086
T9
Spray drift
0.016 ± 0.0016
0.033 ± 0.0092
T10
Road dust/spray drift
0.025 ± 0.0008
0.385 ± 0.0568/0.034 ± 0.0082From the measurements of the particle depolarization ratio, four types of aerosols have been identified: spray drift, diesel exhaust, field dust, and road dust. Spray driftδp ranges between 0.028 and 0.043. The Lorenz-Mie theory states that spherical particles do not modify the original polarization state of the light at the backscattering direction [37]. The sphericity of the droplets was laboratory tested by means of a phase Doppler analyzer (57X10 PDA, Dantec Dynamics, Skovlunde, Denmark); 3% to 20% of the generated droplets at the nozzle output were not spherical, which explains the resulting depolarizations.Depolarization of tractor engine exhaust emissions has been determined (δp=0.099) in trial T1. To the best of our knowledge, only relative depolarization ratios of diesel exhaust have been measured in a previous study [38]. It should be noted that relative depolarizations of diesel exhaust obtained in the cited study are below the background depolarization, which does not happen in our case. As these authors point out, the composition and the depolarization ratio of exhaust emissions is greatly influenced by the combustion efficiency.Field dust is raised by the tractor as it moves along the alleys.δp of 0.220 and 0.268 have been measured in trials T2 and T3, respectively. These results are in good agreement with atmospheric lidar measurements of Saharan dust (δp=0.24–0.28 at 355 nm [36]). The high depolarization associated with the dust is explained by the irregular shape of crystalline particles that constitute it.Trials T4 and T5 are examples of spray drift tests where the cloud is contaminated by dust. In fact, resulting depolarization ratios (δp≈0.21) demonstrate that dust predominates over spraying in the monitored clouds. Polarization lidar technique opens the possibility of computing the fraction of each type of aerosol in mixture clouds as presented here. For instance, δp is currently used to separate dust from nondust in atmospheric research [39]. Its application in farming requires a laboratory characterization of δp of main agricultural aerosols.Road dust raised by the car (trial T10) presentsδp=0.385, a value similar to that obtained in previous road dust measurements (δp=0.33−0.37 at 355 nm [40]) carried out under controlled conditions. The higher δp of the road dust in comparison to field dust is attributed to a different composition of the road floor, with an important presence of gravel.Figure3 represents, by means of normalized histograms, the distribution functions of the particle depolarization ratio for all the types of aerosols analysed in this study. This graph evidences the ability to distinguish and identify each of the agricultural aerosols through δp.Figure 3
Distribution of particle depolarization ratio for spray drift (T6–T10), diesel exhaust (T1), mixture of spray drift and dust (T4 and T5), field dust (T2 and T3), and road dust (T10).As discussed in Section2.3, the two main sources of error in δp are the cross-calibration constant K∗ and the backscatter ratio Rp. K∗ can be affected by the existence of a small amount of aerosols in the calibration layer [36]. From (5), if δv was 0.005, that is, 20% higher than the molecular depolarization ratio δm, then K∗ would vary by −9%.To computeRp, ideal clouds (with no extinction) have been assumed in (10) and (11). For dust measurements, the uncertainty of Rp can be estimated by assuming an extinction coefficient of 0.15 km−1 at 355 nm [41] and a cloud depth equal to 25 m (typical value for the clouds of this study). This yields a variation of −1% on Rp.The numerical calculation of (12) and (13) for previous uncertainties on K∗ and Rp gives relative errors of about 10–12% and 17–19% on particle depolarization ratios of spray drift and dust, respectively. These are acceptable error figures, similar to those obtained in other lidar polarization experiments [33].


## 4. Conclusions

The results show that the detection of agricultural aerosol emissions with the polarization lidar technique is feasible and allows the simultaneous measurement and discrimination of aerosols of different origin. Measurements of spray drift, diesel exhaust, field dust, and road dust emissions had correspondingδp values divergent enough to adequately identify them. This technique is also able to monitor the spatiotemporal evolution of PM clouds, which is a powerful method for understanding this phenomenon. The obtained results encourage further research to improve this technique by developing polarization lidar systems specifically designed for agricultural applications. Future work should also include a laboratory characterization of particle depolarization ratios due to agricultural aerosols. A robust knowledge of δp will allow to determine the fraction of each aerosol in mixed clouds usually caused by agricultural activities.



---

*Source: 1864106-2018-04-18.xml*