Polarization Lidar Detection of Agricultural Aerosol Emissions

Research Group in AgroICT & Precision Agriculture, Department of Agricultural and Forest Engineering, Universitat de Lleida (UdL)-Agrotecnio Center, Lleida, Spain Remote Sensing Laboratory, Department of Signal Theory and Communications, Universitat Politècnica de Catalunya (UPC)/IEEC, Barcelona, Spain Centro de Agroingeniería, Instituto Valenciano de Investigaciones Agrarias (IVIA), Moncada, Valencia, Spain Department of Agriculture, Livestock, Fisheries and Food, Generalitat de Catalunya, Lleida, Spain


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 (PM 2.5 ) and 25% of all particles less than 10 μm (PM 10 ) [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 PM 10 emissions from land preparation operations [6][7][8][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 threewavelength 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][14][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][20][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 wellestablished 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.

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

2
Journal of Sensors ensured that the aerosol was always measured within the full overlap area.
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 (Figure 1(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 Figure 1(b), the sprayer position was adjusted in each trial to ensure that the spray drift always reached the measurement area.
Table 1 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.
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 T 2 R 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 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: 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 where T 1 and T 2 are the transmissions of the parallel polarization contributions of the two Brewster plates, estimated as T 1 = T 2 = 0 92 (Dr.P. Chazette, private communication).

Journal of Sensors
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 constant K * = 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.
With K * 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] where R p is the backscatter ratio, given by 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: 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 tractoroff 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 R p is computed as the ratio of the total (background + agricultural aerosol), U tot , to the background range-corrected power, U b , where 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 < R max , with R max being the maximum range.T 2 b R and T 2 p R 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, T 2 p R ≈ 1, or, equivalently, that the optical depth in a, b is much smaller than unity, τ P = b a α P u du ≪ 1, where α p is the particle extinction coefficient.

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 rangecorrected 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 ΔR p 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 where subscript i = b, v stands for "background" or "volume" measurement, respectively.When ( 12) is combined with the basic error propagation of (8), Journal of Sensors the sought-after absolute error Δδ p in response to absolute errors ΔK * and ΔR p is obtained.

Results and Discussion
3.1.Volume Depolarization Ratio Measurements.Figure 2 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 righthand 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.Trial T9, shown in Figure 2(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 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.

Particle Depolarization Ratio Measurements.
For each trial, Table 2 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).
From 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.
Figure 3 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 .
As discussed in Section 2.3, the two main sources of error in δ p are the cross-calibration constant K * and the backscatter ratio R p .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 compute R p , ideal clouds (with no extinction) have been assumed in (10) and (11).For dust measurements, the uncertainty of R p 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 R p .
The numerical calculation of ( 12) and ( 13) for previous uncertainties on K * and R p 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].

Conclusions
The results show that the detection of agricultural aerosol emissions with the polarization lidar technique is feasible

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
1Individual nozzle flow rate at a pressure of 1 MPa.

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.