Particulate-matter emission estimates from agricultural spring-tillage operations using LIDAR and inverse modeling

Abstract. Particulate-matter (PM) emissions from a typical spring agricultural tillage sequence and a strip–till conservation tillage sequence in California’s San Joaquin Valley were estimated to calculate the emissions control efficiency (η) of the strip–till conservation management practice (CMP). Filter-based PM samplers, PM-calibrated optical particle counters (OPCs), and a PM-calibrated light detection and ranging (LIDAR) system were used to monitored upwind and downwind PM concentrations during May and June 2008. Emission rates were estimated through inverse modeling coupled with the filter and OPC measurements and through applying a mass balance to the PM concentrations derived from LIDAR data. Sampling irregularities and errors prevented the estimation of emissions from 42% of the sample periods based on filter samples. OPC and LIDAR datasets were sufficiently complete to estimate emissions and the strip–till CMP η, which were ∼90% for all size fractions in both datasets. Tillage time was also reduced by 84%. Calculated emissions for some operations were within the range of values found in published studies, while other estimates were significantly higher than literature values. The results demonstrate that both PM emissions and tillage time may be reduced by an order of magnitude through the use of a strip–till conservation tillage CMP when compared to spring tillage activities.


Introduction
As aerosols have been shown to have detrimental effects on human health and visibility, 1 many governments have set regulations on allowable ambient concentrations.In the U.S., the Environmental Protection Agency (EPA) has established National Ambient Air Quality Standards (NAAQS) for particulate matter (PM) with aerodynamic equivalent diameters ≤10 μm (PM 10 ) and PM with aerodynamic equivalent diameters ≤2.5 μm (PM 2.5 ).If an area exceeds the NAAQS, the area's air-quality governing body is required to identify the causes and to restrict anthropogenic emissions in order to reduce PM levels below the standard.
The San Joaquin Valley of California, USA, was designated as noncompliant with the PM 10 NAAQS in 1991 and given a "serious" classification in 1993. 2,3Rule 4550, one of the regulations enacted by the San Joaquin Valley Air Pollution Control District (SJVAPCD) to reduce ambient PM 10 levels in the San Joaquin Valley, required agricultural production operations to select several conservation management practices (CMPs) from a provided list, to submit their selections for SJVAPCD approval, and to implement approved CMPs.The CMPs were designed to reduce PM 10 emissions from agricultural animal and crop production activities.However, the small amount of data available in the literature concerning the emissions' reductions from the CMPs for crop production tillage activities required that the control efficiency (η) of most tillage CMPs was estimated from emissions measurements of other operations. 4The η of a CMP is the amount of particle emission reduction achieved relative to the conventional management practice.While the San Joaquin Valley Air Basin has since been classified as in attainment with the PM 10 NAAQS, its maintenance plan requires that the same strategies employed to bring it back into attainment continue to be applied.In addition, other PM 10 nonattainment areas such as Imperial Valley, California, and Phoenix, Arizona, have CMP or best management practice (BMP) rules in place for agricultural-tillage practices that are based on limited emissions measurements.
Previous agricultural-tillage PM emissions studies 5-12 have focused almost exclusively on measuring PM 10 emission rates (ERs) and emission factors (EFs) from conventional tillage operations.For this discussion, EFs are emissions based on a quantity of production (e.g., g • m −2 ) and ERs are emissions that include a time factor (e.g., g • m −2 • s −1 ).The California Air Resource Board (ARB) developed area source PM 10 emission inventory calculation methodologies for agricultural tillage and harvesting operations based on the report by Ref. 5 (see Refs. 13, and 14).References 10 and 11 are the only instances of reporting PM emissions from standard tillage operations and a CMP (strip-till conservation tillage and combined operations, respectively).Reference 6 used elastic LIDAR (light detection and ranging) data collected during tillage emissions measurements to track plume movements in the downwind vertical plane and demonstrated plume depths were greater than the elevated point sensors located downwind at 10 m above ground level (agl).The report suggested the best method for sampling fugitive dust includes a combination of elastic LIDAR and strategically placed point samplers.Reference 11 used both filter-based mass concentration point samplers and a LIDAR system to monitor tillage emissions and estimated EFs from both datasets.
The study described in Ref. 11 was specifically initiated to provide more emissions η data with respect to a CMP in Rule 4550.The focus was on a typical fall tillage operation after a row crop harvest.A companion study funded by the San Joaquin Valleywide Air Pollution Study Agency was conducted to measure η of a spring tillage CMP using the same point sensor and LIDAR methodology.Research questions which this study was designed to address included: (1) what are the magnitude, flux, and transport of PM emissions produced by agricultural practices for row crops where tillage CMPs are being implemented versus, the magnitude, flux, and transport of PM emissions produced by agricultural practices where CMPs are not being implemented?(2) What are the values of η of equipment being used to implement the "combined operations" CMP? and (3) can these CMPs for a specific crop be quantitatively compared, controlling for soil type, soil moisture, and meteorological conditions?It is important to note that the main focus of this research was to quantify η of the selected CMP, which required the emissions to be quantified, and it was not an effort to provide representative EFs for any one of the agricultural operations involved.This paper summarizes the results of the PM measurements made during the field experiment, the calculated ERs, and addresses these research questions.A full report detailing all of the sampling methodology and results is given in Ref. 15.

Site and Operation Description
This CMP η study was performed during the spring tillage sequence following the harvest of a winter wheat crop in preparation for the planting of corn.It was carried out in the San Joaquin Valley of California during May and June 2008.Two adjacent fields were used with identical crop and flood irrigation treatments over the previous several years.Both fields were cultivated in winter wheat in late 2007 and were to be planted in corn for the 2008 summer growing season.The wheat was harvested while still green for silage 4 days before the tillage processes began, resulting in standing stubble but little plant material left on the surface.The site was chosen based on producer cooperation, historically dominant northwest winds, and field layout.
The surrounding landscape was topographically flat and dominated by agricultural production, including grain and corn fields, almond orchards, grape vineyards, and commercial dairy operations.Both fields were surrounded on all sides by roads.These roads, with the exception of one, were unpaved roads used for field access by farm machinery.The paved road, which was downwind of the fields during all measurements, was heavily traveled.Railroad tracks were located to the north of this site, with two to three trains passing by per day with varying numbers of cars.USDA NRCS soil survey data list the soil in both fields as soil type 130 -Kimberlina fine sandy loam, saline-alkali. 16 The CMP selected for this study was the conservation tillage method.As described in Ref. 4, this CMP "involves using a system in which the soil is being tilled or cultivated to a lesser extent compared to a conventional system" and it is "intended to reduce primary soil disturbance operation such as plowing, disking, ripping, and chiseling."The conservation tillage CMP under study was a strip-till method which combines multiple operations to reduce the number of passes required and only disturbs the soil in strips 0.2 m wide centered every 0.8 m instead of disturbing the entire surface.Strip-till reduced both the number of passes and the tilled surface by about 75%, as well as left most of the wheat stubble still standing for ground cover.The strip-till implement used in this study was the Orthman 1-tRIPr.The cooperating farm had been using the Orthman 1-tRIPr for seedbed preparation on all of its fields for several years, with the exception of field 4, which was prepared by conventional tillage methods.
The conservation tillage CMP applied in this study consisted of three tillage activities totaling three operations across the field, excluding the building and removal of ditches.All three were monitored in separate sample periods.In comparison, the conventional tillage method as applied here had nine different tillage activities totaling 13 operations, excluding the building and removal of a ditch and field-edge borders.Monitoring of 12 of the operations occurred over nine sample periods.Note that not all operations were active throughout sample periods in which multiple operations occurred.Also, not all of the conventional tillage operations are intended to work the entire field, such as breaking down in-field borders.The term in-field borders as used here applies to low ridges of soil that separate the field into smaller areas for flood irrigation.The conventional tillage method was employed in field 4, shown in Fig. 1, and the conservation tillage CMP was used in field 5.The operations that were performed in each management practice are shown in order in Table 1, with their corresponding dates, equipment utilized, number of passes over a given area, tractor run time, total area worked, and sample period length.In cases where multiple tractors and implements were used within a sample period, they are listed in the order of use with the area worked and tractor time being summed.The lister, also called a double plow, prepares the soil for planting by creating furrows and ridges.The cultivator passes in the conventional tillage sequence function as mechanical weed control, whereas a chemical weed control (herbicide) is used in the CMP sequence.
During the first part of the lister operation, plant material not harvested along in-field borders caused clogging of the lister, decreasing effectiveness.A second tractor with a disk set was brought in to repeat the effort along the border lines to further reduce the size of residual plant material.Also note that the cultivator passes 1 and 2 and the roller pass for the first sample period on June 5 were not finished when planting in field 4 began and the second sample period was started.The combination of the cultivator and roller positions and meteorological conditions prevented significant impacts from these operations on downwind samples located near the southern end of the field.
Additionally, cultivator pass 4 was carried out the day after cultivator pass 3, but the emissions were not measured due to scheduling conflicts.It is assumed, in calculating the total PM emissions, that the ERs of both passes 3 and 4 were equal.In general, two cultivator passes are performed in sequence in opposite directions down the rows to ensure adequate weed control.
Ditches and field-edge borders were built and then broken down in both fields between May 20 and June 5 to allow for flood irrigation prior to planting.Water for irrigation was taken from the earthen holding pond of the adjacent dairy; drainage ditches on the east side of both fields returned excess water to the same holding pond.As the ditch and field-edge border construction and removal were not measured in the CMP field, the corresponding step for the conventional tillage method was not considered in the total emissions per method.Prior to any spring tillage activities, both fields 4 and 5 had in-field borders running in roughly an east/west direction.The in-field borders in field 5 were not broken down and smoothed out, but instead were used for the summer corn crop.However, in field 4 they were removed and the irrigation water moved in the furrows created by the lister.
Field personnel observed operations continually and recorded notes on tractor operation times, potential contamination issues due to traffic on surrounding dirt roads and wind-blown dust, general meteorological observations, and so on.

Instrumentation and Sample Layout
A very dominant northwest wind was found in historical data for the months of May and June from a representative meteorological monitoring station in Stratford, California, in the California Irrigation Management Information System (CIMIS).Therefore, the PM and meteorology sampling layouts were configured to measure upwind conditions to the north and west and downwind conditions to the south and east.
Meteorological measurements were made at upwind and downwind locations with the instrumentation in Table 2. Vertical temperature, humidity, and wind speed profiles were measured using two 15.3 m towers, one upwind and one downwind as shown in Fig. 1.Each tower had five humidity/temperature sensors at 1.5, 2.5, 3.9, 6.2, and 9.7 m agl and 3-cup anemometers at 2.5, 3.9, 6.2, 9.7, and 15.3 m agl.Wind direction was measured at 15.3 m using a wind vane instead of the typical 10 m due to the fact that LIDAR measurements were made at higher elevations (up to 200 m agl) and the 15.3 m measurement height was reasoned to provide a better representation  of both ground level and higher elevation wind direction than the 10-m height.Additionally, a meteorological station monitored wind speed, wind direction, temperature, relative humidity, precipitation, barometric pressure, and incoming solar radiation at 5 m agl at the air-quality trailer (AQT) location.Three pairs of three-dimensional sonic anemometers and infrared gas analyzers were deployed, one at an upwind location and one each downwind of the two fields of interest, to characterize upwind and downwind turbulences, as well as vertical fluxes of latent heat (evaporation), sensible heat, carbon dioxide, and horizontal momentum.Bulk density and soil moisture were quantified several times throughout the study, with calculations performed as described in Ref. 17. PM mass concentrations were monitored by 20 MiniVol Portable Air Samplers (Models 4.2 and 5.0, AirMetrics, Eugene, Oregon), referred to hereafter as MiniVols.They are a portable, programmable, filter-based sampler that is battery-powered and yields an integrated sample over the exposure period.Filters were exposed for the duration of each sample period (see Table 1), yielding a single mass concentration measurement per sampler per sample period.Cumulative samples of particles up to PM 2.5 or PM 10 are accomplished via an impactor plate assembly inserted just upstream of the filter; total suspended particulate-matter (TSP) may be sampled without an impactor assembly in place.The MiniVol impactor assembly is designed to operate a flow of 5.0 L • min −1 , though the flow rate is neither actively monitored nor actively controlled by the system.It is set prior to deployment by the user via a calibrated rotameter.Rotameter flow calibration was performed prior to deployment.While several studies have found that PM 2.5 and PM 10 levels reported by MiniVols are very similar to concentrations measured by federal reference method (FRM, see 40 CFR 50.6 and 50.7) monitors, the slope of the particle removal efficiency versus particle size curve of the MiniVol impactor assembly is less steep than required by FRM samplers. 18-21Therefore, PM levels reported by the MiniVols should be considered as close approximations to those that would be given by FRM samplers.Pre-and postweights for the 47 mm Teflon filters used to collect particles were quantified using a calibrated microbalance (Type MT5, Mettler-Toledo, Inc., Columbus, Ohio).Filter conditioning was carried out in accordance with guidance in 40 CFR 50 Appendix J. Sample period average mass concentrations were calculated by dividing the total mass catch (average postweight minus average preweight) by the volume of air sampled.
Particle size distribution (PSD) was measured by eight Aerosol Profilers (Model 9722, Met One Instruments, Grants Pass, Oregon), also known as optical particle counters (OPCs).The OPC sums particle counts in eight size bins over nominal 20-s sample periods.The particle diameters (d p ) for lower bin limits were 0.3, 0.5, 0.6, 1.0, 2.0, 2.5, 5.0, and 10.0 μm, with the last channel counting all particles ≥10.0 μm.The factory calibrations of signal strength versus particle size using polystyrene latex beads of known size were used due to varying atmospheric aerosol composition.The sample flows were not conditioned prior to passing through the sampling chamber during these measurements due to the dry conditions, though this is suggested in atmospheres with high relative humidity.OPC flow measurements and inter-OPC count calibrations were performed onsite and applied in postanalysis.Particle volume concentrations (V) per bin were calculated from the counts, assuming spherical particles and using the geometric mean diameter (GMD) as the representative d p .Values of V in bins up to d p ¼ k were summed to estimate the cumulative volume concentration (V k ).
The MiniVols and OPCs were deployed on towers and tripods upwind and downwind of the fields, as shown in Fig. 1.Most downwind sensors were moved between the downwind layouts, depending on the field being tilled.The AQT and the associated samplers did not move.The AQT is a 5 m × 2.5 m × 2.5 m cargo trailer used as the base of operations and equipped with tables, a refrigerator and dessicator for sample storage, and a rooftop platform for sensor deployment.Samplers were placed on tripods at 2 m agl at all locations except for those on top of the AQT at 5 m and those at the top of the towers at 9 m.Not all sample sites shown for a given setup were used in each sample period due to instrument availability limitations.However, samples were collected at a minimum of two upwind and six downwind locations during each sample period.
At most of the locations, multiple MiniVols with different impactor configurations and an OPC were collocated in order to characterize particle size and mass distributions.These data were used to calculate mass conversion factors (MCFs) for each size fraction (k), as described in detail by Ref. 21.In summary, the MCF k is calculated using PM k reported by MiniVols and V k , averaged over the MiniVol sample time, from each sample location through the following equation: where the units for PM k are μg • m −3 , V k are μm 3 • cm −3 , and MCF k are g • cm −3 .Daily average MCFs were calculated across sampling locations.
The MCF is a simplified method to account for several complex and possibly interdependent variables that affect how an aerosol mixture is measured/detected based on both optical and aerodynamic properties.It incorporates many factors, such as particle shape, density, indices of refraction different from OPC calibration aerosols, and instrument sampling efficiencies, into a single scalar value.The MCF also includes effects due to optical systems measuring particles in ambient conditions while mass concentrations are calculated based on conditioned filters.This effect may be significant in humid environments, but the effect is assumed to be negligible in warm and dry conditions such as those found during this study.The MCF k values were used to convert the OPC V k into PM k to examine concentrations on a much finer temporal scale than possible with the filter data.The Aglite LIDAR system was deployed to characterize PM concentrations in addition to point sensors.The Aglite LIDAR is a portable system using a micropulsed Nd:YAG laser with three wavelengths (λ), 355, 532, and 1064 nm.It has the capability to scan 280 deg in azimuth and from −10 deg to þ45 deg in elevation.The effective range is 500 m to 15 km with each range bin approximately 6 m in length. 22The LIDAR was placed in crosswind positions 550 m from the nearest tillage area border.It was at L1, as shown in Fig. 1, from May 17 through June 11 and at L2 for the June 18 sample period.Critical component failures prevented its use for the June 25 sample period.
The LIDAR continuously performed vertical scans on the upwind and downwind sides of the field, horizontal scans over the field, and calibration stares throughout tillage observation periods.Lines of approximate vertical scan locations are shown in Fig. 1 by the dashed lines emanating from L1 and L2; horizontal scans moved between upwind and downwind vertical scan positions at 0.75 deg from horizontal.Vertical scans started at 0.75 deg and extended up to between 15 deg and 45 deg.The maximum vertical angle varied between sample periods but was usually ≤25 deg.Images resulting from vertical scans were monitored throughout each sample period to ensure that the maximum vertical extent of the plumes was entirely captured; modifications to the maximum vertical extent were made as needed.The LIDAR beam was about 10 m agl at the closest edge of the fields at 0.75 deg in elevation.The beam was kept at or above this level due to eye safety concerns.The LIDAR system did not measure plumes below this level and, therefore, may underestimate PM flux.
A calibration stare refers to short periods (60 to 120 s) when the LIDAR beam is held adjacent to an upwind tower with collocated point sensors.Calibration stares were routinely performed throughout the sample period at 10 to 20 min intervals.In postprocessing, LIDAR return signals collected during calibration stares were calibrated to PSD and V k measurements.The process used to accomplish this is described in detail by Refs.23 and 24 and will now be summarized.The calibration process is illustrated in Fig. 2 and was carried out through the following steps: 1.The raw LIDAR signal was preprocessed to yield range (R) and background corrected return power with R. 2. Relationships between backscatter (β), extinction (α), and V k of the aerosol components were established based on OPC data from both upwind and downwind locations.The PSD of both the background and plume aerosol as a function of time was calculated, after which the α and β coefficients at the calibration range (R c ) at each LIDAR λ were calculated using Mie scattering theory applied to the PSDs.Assumptions made in these calculations were (1) all particles were spherical and (2) the bulk aerosol had a complex index of refraction equivalent to a mineral particle type (1.53-0.008i;Ref. 25).
Fig. 2 Process diagram for LIDAR PM calibration algorithm.
3. The inversion of the LIDAR data was performed using a form of Klett's solution 26 for two scatterers where α is proportional to β using the relationships found in step 2, resulting in α and β as function of R and λ [αðR; λÞ, βðR; λÞ].The backward integration Klett method was applied to R < R c and the forward integration method was used for R > R c .The backward integration method yields more stable solutions than the forward integration method and is, therefore, preferred.The calibration point R c was placed at the farthest range possible within property ownership/field of view limitations to maximize the extent of the range of interest subject to backward integration.4. The relationships from step 2 were used to convert βðR; λÞ into V k ðRÞ through a leastsquares method in the aerosol concentration retrieval step. 5. Conversion from V k ðRÞ to PM k ðRÞ was accomplished through the use of MCF k .This calibration method allows a scanning LIDAR to estimate PM k concentrations surrounding an area/source of interest at a much finer spatial scale than possible with point sensors.

Emission Calculation Methods
The PM k data from the point sensors and LIDAR were used to the estimate η of the conservation tillage CMP in this study.The point sensor PM k concentrations were coupled with an air dispersion model through inverse modeling.In typical air dispersion model applications, a source ER is supplied to a model which then calculates the resulting concentration (C sim ) at a given receptor location.Inverse modeling involves using a dispersion model and concentrations measurements around a source activity (C downwind , C upwind ) in order to estimate the observed ER (Q meas ).An initial ER (Q sim ) is supplied to the model to calculate C sim , then the following equation is used to calculate Q meas : If the model used has a proportionally linear response in C sim to changes in Q sim , the ratio ðC∕QÞ sim is a scalar value independent of the given Q sim value, i.e., there are no local maxima or minima that might influence the resulting value of Q meas .
AERMOD (American Meteorological Society and U.S. EPA Regulatory Model), a Gaussian air dispersion model that estimates C sim at a given receptor point based on meteorological conditions, source strength, and the horizontal and vertical distances of the receptor from the source, was selected to perform the inverse modeling estimation of Q meas . 27It is an air dispersion model currently recommended for regulatory modeling by the U.S. EPA.It operates in 1 h time steps, has a proportionally linear response in C sim to changes in Q sim , and assumes steady-state conditions, continuous emissions during a time step, conservation of mass, and C sim resulting from multiple sources are additive.Pollutant distribution is modeled as Gaussian in the stable boundary layer in both the horizontal and vertical directions.In the convective boundary layer, horizontal dispersion is modeled as Gaussian while vertical pollutant distribution is modeled as bi-Gaussian.The spatial resolution of C sim is controlled by the user through discrete and/or gridded receptor points.The commercially available user-interface AERMOD View from Lakes Environmental, Inc. (Waterloo, Ontario, Canada), with AERMOD version 13350, was employed.
Onsite measured wind speed, wind direction, temperature, humidity, and solar radiation were used by AERMET, the meteorological preprocessor for AERMOD, to create both surface and elevated meteorological input files.Wind, temperature, and humidity data were used from the upwind meteorology tower dataset, with wind speed from 9.7 m agl and temperature and relative humidity from 2.5 m agl.Incoming solar radiation was measured at the AQT location.Percent cloud cover was set to zero based on visual observations during the measurement periods.
The land-use classification on all sides of the site was cultivated land.Values provided to AERMET for midday albedo and Bowen ratio of 0.18 and 1.5, respectively, were suggested as average summer values under dry conditions for a fallow agricultural field in Ref. 28.The surface roughness length (z 0 ), also required by AERMET, was calculated based on wind profile measurements at the upwind meteorological tower using the following equation which relates wind speeds (u 1 , u 2 in m • s −1 ) at two heights (z 1 , z 2 in m) and was derived from the integrated logarithmic wind speed profile equation: A least sum of squares of residuals methodology was used to determine the value of z 0 that best fit measured wind speeds at the higher elevation of two paired wind speed time series over the study period.A z 0 value of 0.02 m was calculated as the arithmetic average of the values that best fit six pairings of wind speeds measured at 2.5, 3.9, 6.2, and 9.7 m agl, i.e., 2.5 and 3.9 m, 2.5 and 6.2 m, 2.5 and 9.7 m, 3.9 and 6.2 m, 3.9 and 9.7 m, and 6.2 and 9.7 m.The cup anemometer at 15.3 m malfunctioned during this deployment, rendering the data unusable for this analysis.The upwind location was selected for this analysis as the downwind tower was removed on June 12 to support another study nearby.The AQT was also removed from the southern edge of field 5 at the same time.
Tillage operations were modeled as ground level area sources with initial plume heights of 0 m and areal extents equal to the actual tilled portions of the field.Most operations covered all or most of the fields within a sample period, but some, such as the break down in-field borders operation, were intended to only work a small portion of the field surface.Tilled areas and sampler locations were measured using a hand-held GPS unit.The Q sim values for each modeled operation were based on a preliminary average ER value across all tillage operations from Ref. 11 of 8.6 μg • s −1 • m −2 per operation per pass multiplied by the number of passes over the field within a sample period.Sources were activated or deactivated in hourly time steps throughout a simulated sample period according to the tractor operation times as monitored by onsite personnel.Discrete receptors were set at each sampling location to yield C sim for inverse modeling comparisons.Uniform Cartesian receptor grids at 2 m agl and 15 m spacing between points were set from upwind sampling locations to several hundred meters downwind of the fields to visualize predicted plume movement, shape, and concentration.Hourly C sim values were averaged over the modeled sample period.Modeled plume edge effects were avoided by eliminating those locations with C sim less than 10% of the maximum C sim , adapted from suggestions by Ref. 29, from emissions calculations.
The second ER and EF calculation approach was a mass balance applied to the LIDAR PM k data.Assuming uniform background aerosol levels, average upwind concentrations were subtracted from concentrations in and around detected plumes in the downwind vertical scans.The difference was multiplied by the component of the wind perpendicular to the beam to calculate the horizontal flux of PM through the downwind vertical scanning plane.Fluxes were summed across the vertical plane, averaged over the length of the sample period, and then divided by the size of the tilled area to calculate the mean EF of PM k from the field surface.The EF was further divided by the total tractor time to calculate the mean ER of each operation.This method of calculating ERs and EFs using LIDAR is described in detail in Ref. 30.
Vertical profiles of PM mass concentration, horizontal wind speed, and wind direction are required to use the mass balance approach.The PM profile was provided by the LIDAR PM k data.Profiles of wind speed were calculated using the wind speed power law, as given by Ref. 31, where z 1 and z 2 are the lower and upper elevations (m), respectively, p is a dimensionless number that varies with atmospheric stability, and u 1 , u 2 , z 1 , and z 2 have been previously defined.
Ref. 31 lists p ≈ 0.5 for very stable conditions and p ≈ 0.15 for very unstable conditions.The horizontal wind speeds recorded at the upwind tower were used to find the values of p that best fit the time series of measured profiles up to 9.7 m agl, with imposed minimum and maximum limits of 0.1 and 0.6, respectively.Derived values of p across all sample periods ranged from 0.10 to 0.60 and averaged 0.19.These p values were then combined with the minute-averaged wind speeds in the wind speed power law to calculate the vertical profile of horizontal wind speed up to 250 m agl, though most sample periods did not require data more than 150 m agl.
Wind direction over the vertical profile was assumed to be constant.Though wind direction is known to change in a vertical profile, the influencing factors may be complex and the magnitude and direction of change highly variable.Therefore, in the absence of measured data, the assumption that wind direction did not change with increasing elevation over the 250 m profile was used.

Results and Discussion
Results of the soil analyses were almost identical between the two fields, suggesting little to no difference in the influence of soil properties on airborne PM emissions.Bulk densities averaged 1.57AE 0.05 g • cm −3 for field 4 and 1.57AE 0.08 g • cm −3 for field 5. Unless otherwise noted, error values represent one standard deviation (σ).Average soil moisture values in both fields measured immediately prior to May sample periods varied between 1.0% and 3.3%, showing very little change across operations.However, average soil moisture measured on June 5 in field 4 was 6.1% and 8.2% in field 5 on June 7.This increase was likely due to flood irrigation in both fields shortly after the May 20 sample period.A precipitation event occurred shortly after irrigation, but the quantity was not measured and the effect was assumed to be masked by the flood irrigation.
Wind conditions were favorable for the designed sampling layouts throughout the study period, as shown in Fig. 3 and Table 3.The median, minimum, and maximum values for sample period average winds were 4.0, 1.9, and 5.6 m • s −1 , respectively, for speed and 321 deg, 315 deg, and 335 deg, respectively, for direction.Sample periods were generally hot and dry, with median, minimum, and maximum sample periods average temperatures and relative humidity values of 31.4 °C, 24.7 °C, and 36.8 °C and 27%, 16%, and 40%, respectively.

Particulate-Matter Concentration Measurements
A total of 296 filter samples were collected: 116 PM 2.5 (39%), 116 PM 10 (39%), and 64 TSP (22%).Calculated PM 2.5 concentrations based on filter catch ranged from 23  Of the 296 filter samples collected, 98 (33%) did not pass quality analysis (QA) checks applied to the dataset.QA checks included visual inspection of filter surfaces, sample log inspection for noted problems, OPC time series examination for contamination (used mostly at upwind sites), concentration consistency across sampling locations, and concentration comparisons among PM 2.5 , PM 10 , and TSP at each sample location.In-depth descriptions of the QA checks are found in Ref. 15.An investigation into the cause(s) of this high rate of failure was conducted and a summary of conclusions is provided in the following paragraph.A large number of failures of near-source, downwind samples relative to the total number of downwind samples collected were found in the sample periods from May 18 to May 20.This, when combined with the results from the investigation into the large number of failures, casts doubt on the validity of the remaining near-source downwind samples from those runs.Therefore, all near-source downwind MiniVol samples for these sample periods were removed from MCF and ER calculations, rendering the upwind and far-source downwind samples that passed QA unusable for estimating PM emissions.The filter dataset used to calculate ERs and EFs totaled 131 samples (44%).Concentration ranges for this dataset were 26.7 to 149.8 μg • m −3 for PM 2.5 , 47.4 to 489.4 μg • m −3 for PM 10 , and 102.9 to 1; 896.9 μg • m −3 for TSP.The size fraction distribution of filters used to estimate emissions was nearly identical to the total sample set: 51 (39%) were PM 2.5 , 50 (38%) were PM 10 , and 30 (23%) were TSP.
Filters that did not pass QA were found to have been contaminated during one or more of the following stages: sampling, filter handling, and filter storage.Evidence of "particle bounce" was found on many PM 2.5 and PM 10 samples collected during May sample periods.Particle bounce occurs when particles that collide with the impactor plate are re-entrained in the airstream and collected on the filter downstream and result in higher reported levels than actually existed.This issue is most likely due to exposing the MiniVol samplers to dust plumes exceeding the maximum recommended exposure level and improper instrument maintenance and cleaning through the May sample periods.Corrective action in the form of inspection after each deployment and cleaning, if needed, was taken during the June sample periods; no issues associated with particle bounce were observed in the second portion of the study.Additionally, some particles were observed on top of and imbedded into the plastic annular ring around the Teflon filter material-the plastic ring is covered by the filter holder assembly during deployment.This was likely due to contamination during onsite filter storage or handling.Efforts were made to minimize this issue throughout, especially during the June sample periods.However, windblown dust did impact the handling and storage area on May 20.
The collected OPC data were used to calculate PSD, V, and V k values.Unlike the downwind MiniVol samplers, the downwind OPCs were not overwhelmed by the dust plumes from the tillage activities-the manufacturer specified range of the OPC of 0 to 3.18 × 10 8 particles m −3 was never exceeded-and thus provided usable data throughout all sample periods.Background and downwind PSD and V profiles varied throughout the study, as shown in Fig. 4. The particle concentrations in this figure were calculated as the change in number (N) per change in natural logarithm of d p (dN∕d½lnðd p Þ) where d p is the GMD per bin.The particles emitted by the tillage activities were generally large (d p > 5 μm) and, therefore, strongly dominated the volume and mass contributions of the activity to near-source atmospheric particle loadings.Three to four OPCs were in positions immediately downwind of the field under study in each sample period, with between one and four OPCs in upwind locations.
OPC time series data at upwind locations were examined for contamination from upwind activities, such as unpaved road traffic.Contamination was found in six of the 12 sample periods, with five of those occurring at the sample site adjacent to the upwind meteorological tower and immediately downwind of an unpaved road.Large, short-duration spikes indicative of contamination were removed from the upwind data in these instances to estimate the background aerosol concentration.In each instance this was performed, the estimated background levels were in very good agreement with those measured by an OPC at a different, uncontaminated upwind location.Filter samples collected at upwind locations with contamination indicated by OPC data were removed from ER estimation.
Those filter samples that passed QA, including the upwind and far-source downwind samples from May 18 to May 20, were used to estimate MCF k values if collected adjacent to an OPC.Most of the daily average MCF 10 and MCF TSP values were within the expected range of 1 to 3 g • cm −3 .However, the daily average MCF 2.5 values were much larger than expected, with individual values ranging from 3.2 to 28.2 g • cm −3 , having a mean of 14.6 AE 3.7 g • cm −3 and a median of 10.1 g • cm −3 .For comparison, the densities of pure nickel and mercury are 8.9 and 13.5 g • cm −3 , respectively.In past field campaigns, MCF 2.5 has generally been higher than MCF 10 and MCF TSP values, but these MCF 2.5 values were much higher than those seen before and account for the majority of values above 5 g • cm −3 reported in Ref. 21. Due to the nonphysically large numbers, the calculated MCF 2.5 was not used.Instead, the average soil density of 2.65 g • cm −3 given in the USDA NRCS National Soil Survey Handbook was used as a constant MCF 2.5 for all sample periods. 32Using a constant MCF 2.5 may affect the accuracy of calculated PM 2.5 concentrations, ERs, and EFs and the value of η 2.5 .However, as will be shown later, the effect on η 2.5 was assumed to be small as the values were very close to those of η 10 and η TSP .Table 4 presents the daily MCF k values used to convert V k calculated from LIDAR and OPC measurements to PM k .
The cause of the high MCF 2.5 values is unknown.No significant differences in PM 2.5 chemical composition were observed between sample periods with higher and lower MCF 2.5 values.While PSDs varied between sample periods, no trends in PSDs sufficient to explain high/low MCF 2.5 groupings were observed.Higher average MCF 2.5 values were not restricted to sample periods in which evidence of particle bounce was found (May 18 to May 20); those filters exhibiting evidence of particle bounce were removed prior to MCF calculations.MCF 10 and MCF TSP patterns tended to follow MCF 2.5 patterns, having correlation coefficients (r) of 0.64 and 0.84, respectively, but with much smaller changes in amplitude.Good negative correlations (−0.69 ≤ r ≤ −0.64) were found when comparing all MCF k values with sample duration.This means that MCF k tended to increase as sample duration decreased.
One potential explanation consistent with these relationships is contamination during filter handling and storage.If filters were equally contaminated, the greatest effect would be found on Table 4 Mass conversion factors (MCFs) used to convert optical particle measurements to mass concentrations for each sample day and averaged for the whole campaign.Error values represent the 95% confidence interval for n ≥ 3. A constant MCF 2.5 value equal to the average density of soil was used due to nonphysically high values calculated for most of the sample days. 32   those samples with the smallest mass catch, i.e., samplers with PM 2.5 impactor configurations or shorter sample times.Unfortunately, field and lab blanks were not collected to monitor for and quantify such contamination; this oversight has been corrected in subsequent studies.

Date MCF
An alternate contamination test is a comparison with independent and proximate PM measurements, though conclusions from this test are limited by comparability of sample characteristics.The closest independent monitoring site was a suburban monitoring station operated by the SJVAPCD, which reported 24-h average PM 10 concentrations on 3 days during which tillage monitoring occurred.The SJVAPCD 24-h average PM 10 and study site mean background PM 10 pairs were, respectively, 38 and 38 μg • m −3 for May 18, 34 and 47 μg • m −3 for June 5, and 38 and 163 μg • m −3 on June 11.Differences in sample period (24-h versus 1.5-h to 7.5-h and limited to daylight hours only), location setting and local sources (suburban versus rural), and instrumentation (FRM versus MiniVol) existed between the two datasets and contributed to observed differences.The differences in concentration between the two sites on May 18 and June 5 were within the range of expected values and do not support the sample contamination hypothesis.At first glance, the large difference on June 11 may be interpreted to support this hypothesis.However, all filters collected during this sample period easily passed the visual filter inspection for particle bounce and particles on and/embedded into the annular ring.If filter contamination did occur during this field study, it does not appear to have been consistent based on the comparison with proximate PM 10 measurements nor evident in the applied QA tests.Therefore, the cause of the high MCF 2.5 values is unknown and no further data exclusions were made.
Sample period-averaged OPC PM k data were ranged from 4.3 to 60.2 μg • m −3 for PM 2.5 , 41.2 to 641.1 μg • m −3 for PM 10 , and 95.3 to 3;271.9 μg • m −3 for TSP.OPC PM k values at the native 20-s averaging period had much higher ranges across all k as the plumes emitted by the roving tillage activities impacted the point samplers in short bursts, the impact of which was reduced when averaged with intervals measuring lower levels.LIDAR-derived PM k also had a high variability for the same reason.This is evident in Fig. 5, which presents PM 10 reported by all three systems at 9 m agl at the downwind tower for the June 18 sample period (14:00-16:10).The dashed black line represents the sample period average PM 10 value based on the MiniVol filter sample, the red line represents the PM 10 based on 20 s OPC data, and the blue markers represent the LIDAR-derived PM 10 for each 0.5-s signal averaging period.The higher temporal sampling frequency of the optical systems shows the timing and magnitude of individual plumes impacting the collocated group not resolved by the single MiniVol sample.The LIDAR beam was not continuously adjacent to the point sensors as it was performing upwind calibration stares and downwind vertical scans, leading to the gaps seen in the LIDAR time series.The LIDAR reported elevated PM 10 levels when collocated during plume impaction events.It also reported elevated levels when the point sensors did not-this was due to part of a plume being within the 6 m long LIDAR bin but not significantly impacting the tower.Negative PM k values were occasionally reported by the LIDAR, as seen here.These were artifacts of the optimization in the calibration procedure and are not real-negative PM concentrations are not possible.
While the scanning LIDAR was at a disadvantage compared to the point sensors in monitoring PM continuously at a given location in Fig. 5, a significant advantage of the LIDAR over the point sensors is its ability to monitor plumes over a line or area.For instance, the data shown in Fig. 6 were collected in a stare adjacent to the downwind tower (location indicated by the solid black line at constant range) over approximately 1 min of the sample period in Fig. 5.The LIDAR detected multiple, highly concentrated plumes at varying distances, only one of which appears to impact the tower.The bottom plot shows the average PM 10 concentration with range during this time.
The vertical scanning profiles used to monitor plumes emitted by the tillage activities allowed measurement of vertical and along-beam horizontal plume extents.Monitored plumes reached elevations up to 150 m agl at the downwind LIDAR scanning plane, though most remained below 100 m.Plume widths also varied.Figure 7 provides an example of two plumes captured in a single vertical scan on June 5 when two different tillage operations were performed in different areas of the field.The plume closer to the LIDAR is lower and denser than the plume farther away.
Comparisons between PM k concentrations from MiniVol, OPC, and LIDAR data at upwind and downwind locations were performed for each sample period as a check on the use of the MCF k and the LIDAR calibration procedure.Accurate estimates of PM k are necessary for accurate estimates of PM emissions.The MiniVol, as it measures PM k directly at each point, is assumed to be more representative of the actual PM k than the OPC and LIDAR.Most interinstrument comparisons revealed upwind values agreed fairly well, but greater differences were found in comparisons at downwind sample sites.Potential reasons for the similarities and differences at both upwind and downwind point sampler locations will be discussed in the following paragraphs, with an example comparison provided as well.
The similarity of PM k estimated by the OPC and LIDAR at the upwind site was expected as this is the calibration point in the LIDAR PM retrieval algorithm and both are converted to PM k through the same MCF k .The downwind sites, however, were not used as direct calibration points, and thus may exhibit greater differences between LIDAR and OPC period-averaged values.On a side note, the proximity of all upwind OPC and LIDAR PM k demonstrates the selected time interval between LIDAR calibration stares was sufficient to adequately characterize The error introduced in this inter-instrument comparison associated with using a spatially averaged MCF k was generally less than AE30% during this study based on the 95% confidence intervals (CIs) reported in Table 3.Large differences in PM 2.5 between the filter-and opticalbased methods were observed across all sample periods due to the use of a constant MCF 2.5 value that was not derived from onsite PM measurements.However, OPC and LIDAR PM 2.5 levels were close.
An example of an inter-instrument comparison is provided in Table 5 for the June 18 sample period, the period shown in Fig. 5.In this instance, average upwind OPC and LIDAR PM 10 and TSP concentrations were within 15% and 7%, respectively, of the values reported by the MiniVol.Downwind LIDAR PM 10 and TSP values were 73% and 64% of OPC levels, respectively, while the OPC PM 10 was 6% higher than the PM 10 concentration from the filter sample.The differences at the downwind location were likely caused, in large part, by the relatively small amount of time the LIDAR sampled adjacent to the point sensors.Longer stare periods at downwind locations are expected to yield better PM k accuracy and are planned for future deployments.

Estimated Emissions and Control Efficiencies
The PM k data from the LIDAR, OPCs, and MiniVols were all used to calculate PM emissions using mass balance and inverse modeling techniques.Table 6 provides the results of these calculations, as well as the summed emissions from each tillage management practice.Stated uncertainties are the 95% CIs about the average and have been provided for cases where n ≥ 3.
The summed conservation tillage sequence emissions (E ST ) consist of the following three passes: strip-till, plant, and herbicide application.The LIDAR did not detect plumes in downwind vertical scans during the herbicide application, as indicated by the no plumes observed designation.However, downwind PM samplers reported small increases in concentrations over background levels during this sample period, leading to small EFs relative to the other operations investigated.The herbicide application operation was performed by a small tractor pulling a spray applicator; as no subsurface disturbance occurred, the only active PM sources were the tractor and implement tires and the spray droplets.The summed emissions for the conventional tillage method (E CT ) include the following 13 passes in order: two break down in-field border passes, chisel, two disk passes, lister, two cultivator passes, roll, plant, fertilizer injection, and two more cultivator passes.The EFs for the break down in-field borders operation have been distributed over the entire field area (10.1 ha) instead of just the area worked (1.0 ha) to represent the emissions over the entire field in the E CT calculation.As the LIDAR system was unavailable to take measurements during the last sample period for the third cultivator pass, the emissions of the last two cultivator passes were assumed to be equal to the observed emissions for the first two cultivator passes.The inverse modeling method found that the third cultivator pass emitted from 0.6 to 3.6 times as much PM as calculated for the first and second cultivator passes, with an average of 2.1.
Some of the LIDAR-derived EFs were much higher than those calculated by inverse modeling within an operation and PM size fraction, particularly for the chisel, disk 1, and lister passes.Others were not statistically different due to overlapping CIs.E CT and E ST EFs from the LIDAR dataset were significantly higher, based on the CIs, than those calculated through inverse modeling in all cases but the PM 2.5 E ST .The LIDAR-and OPC-based PM 2.5 EFs were smaller than the MiniVol-based PM 2.5 EFs in all cases, which is likely to be related to the use of the average soil density as the MCF 2.5 in place of calculated values as discussed earlier.
While the EFs from published studies are generally not in close agreement, a high range of variability is expected from measurements made under different meteorological and soil conditions, as demonstrated in Ref. 5. Some of the PM 10 EF values calculated from this study were in agreement with those given by Refs. 5 and 10, such as the cultivate, roll, strip-till, conventional tillage planting, and conservation tillage planting passes.Other EFs were much larger than values previously reported, especially the disk 1, disk 2, chisel, and lister passes.The results from this campaign were, in general, not in as good agreement as previous studies have been.
The η values were also included in Table 6.They were calculated using the following equation based on a collection efficiency equation found in Ref. 31 and represents the PM emissions reduction of the conservation tillage compared to the conventional tillage: The strip-till conservation tillage reduced PM emissions in all size fractions by about 90%.The reduction in total tractor operation time per unit area of field was similar to PM reductions at 84% (see Table 1).The η results had a small range of <6% across methodologies and size fractions, despite large differences in summed emissions between methodologies.Limitations of the PM sensors and AERMOD contribute to limitations and uncertainties in the estimated EFs.For instance, the scanning LIDAR system did not collect data below about 10 m agl at the range of the fields due to eye-safety concerns, which results in portions of plumes not being sampled and included in EF estimates.Vice versa, the inability to locate point sensors above 10 m limits their ability to characterize the plume depth.In this case, the LIDAR and point sensors complement each other in mapping the emitted plumes, as previously demonstrated by Ref. 6.In addition, measurements at a few points may or may not represent the plume characteristics sufficiently to accurately determine the EF, particularly for a roving source such as in agricultural tillage.This was mitigated by deploying as many samplers at different sites within the downwind plume as possible.Another limitation that was identified by Ref. 11 and was also present in this analysis is AERMOD's poor simulation of elevated plumes from agricultural tillage, particularly plumes completely detached from the ground.

Conclusions
A study was conducted in California's San Joaquin Valley to estimate the PM k emissions η of a CMP relative to the conventional tillage practices.PM k concentrations resulting from a spring tillage sequence transitioning from a winter wheat silage crop to a summer corn crop were monitored.The strip-till conservation tillage CMP, consisting of three operations in three passes, was compared against the conventional tillage sequence, consisting of nine operations in 13 passes.The CMP reduced the amount of tractor operation time per unit area by 84%.Emissions were estimated through inverse modeling with point sensor PM k and through a mass balance applied to mass-calibrated LIDAR PM k .
A significant portion of the filter-based samples was rendered unusable for emissions calculations due to sampling irregularities and errors.The incompleteness of this PM k dataset prevented the calculation of total EFs per management practice and the CMP η.However, the OPC-and LIDAR-based PM k and EF datasets were sufficiently complete to calculate η values, which were all about 90%.Some of the calculated EFs were within the range found in the literature, but others were significantly higher.The total emissions per management practice and PM size fraction varied significantly, based on the 95% CIs, between the measurement and emissions estimation methodology combinations.This study demonstrated that the strip-till CMP can significantly reduce PM emissions and tractor operation time during the investigated spring tillage sequence.

Fig. 1
Fig. 1 Map of fields under study and the sample layout for each field.L1 and L2 represent light detection and ranging (LIDAR) locations used during sampling and dashed lines emanating from L1 and L2 show positions of vertical LIDAR scans.

Fig. 3
Fig.3Wind rose for the hourly averaged wind observations during the days on which samples were collected.

Fig. 5
Fig.5Time series of PM 10 concentrations as reported by the collocated OPC, MiniVol filter sampler, and LIDAR at 9 m agl on the downwind side of the tillage activity for the June 18 sample period (14:00-16:10).

Fig. 6
Fig.6(a) LIDAR-derived PM 10 in a time versus distance from the LIDAR concentration map and (b) a time series average concentration versus distance from the LIDAR graph.Data were collected in a stare past the downwind tower during the June 18 sample period.

Fig. 7
Fig. 7 Image of plumes in a vertical LIDAR scan on June 5.

Table 1
Information for each sample period regarding tillage operations, equipment used, tractor operation time, area worked and sample time.

Table 2
Manufacturer, precision, and accuracy information for deployed meteorological instrumentation.
Offset error < AE 8 cm s −1 Gain error for wind vector within 20 deg of horizontal < AE 6% of reading .2 to 3244.9 μg • m −3 ;

Table 3
Period-averaged meteorological measurements AE 1σ made at the upwind meteorological tower.Temperature, relative humidity, and wind speed were measured at 9.7 m agl and wind direction was measured at 15.3 m agl.

Table 5
Comparison of period average PM mass concentrations as reported by collocated MiniVol filter samplers and OPCs, as well as the adjacent LIDAR bin, at measurement heights of 9 m agl at upwind and downwind tower locations for the June 18 sample period.Error values provided for OPC and LIDAR PM k represent the 95% confidence interval (CI).

Table 6
Average particulate matter (PM) emission factors (EFs) and 95% confidence intervals (CI) estimated for the conventional and conservation tillage management practices.Emissions were calculated through applying a mass balance technique to mass-calibrated LIDAR data and inverse modeling with AERMOD and both filter-based particle sampler (MiniVol) and mass-calibrated OPC data.Control efficiency (η) results are also provided for the conservation tillage management practice.NPO ¼ no plumes observed.