|
|
||||||||
a Agriculture and Agri-Food Canada, P.O. Box 3000, Lethbridge, AB, Canada T1J 4B1
b Department of Earth and Atmospheric Sciences, University of Alberta, Edmonton, AB, Canada T6G 2E3
c Agricultural Research Service, United States Department of Agriculture, Watkinsville, GA 30677
* Corresponding author (mcginn{at}agr.gc.ca)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: (C/F)sim, theoretical ratio of concentration to emission (Farea/Fpoint), ratio emissions based on area to that based on point configure of source GC, gas chromatograph h, obstacle height LS, Lagrangian stochastic xs, source separation distance
| INTRODUCTION |
|---|
|
|
|---|
A few studies have measured CH4 emissions using close approximations to real farms. Kinsman et al. (1995) calculated CH4 emitted from a barn housing dairy cows, by monitoring inflow and exhaust air concentration and airflows. However, this technique, aside from being technologically challenging, cannot be applied to the majority of more modern "open" dairy barns or beef feedlots where the airflow is passive. A more appropriate approach for open facilities is to release a tracer gas at a known rate (Marik and Levin, 1996; Kaharabata et al., 2000). This assumes the relationship of the tracer flux to the tracer concentration in the downwind plume is similar for that of the target gas (where the concentration is also measured). Although this ratiometric technique is theoretically straightforward, it relies on matching the tracer and target gas sources that can be difficult to achieve (because the true emission configuration is unknown) and requires extensive setup time.
A more flexible technique for quantifying emissions is to model the dispersion of a target gas from the source (e.g., Kaharabata et al., 2000; Flesch et al., 2004), so that a downwind concentration of gas can establish the emission rate. This "inverse-dispersion" technique has the potential advantage of simplicity, as it requires only a single gas concentration measurement and basic wind information. However, most of the applications of this technique have been for emission sources over a uniform surface (e.g., Wilson et al., 1982). For more complex settings, like a farm with buildings and trees, further evaluation of the inverse-dispersion approach is needed.
The objective of our study was to examine the application of an inverse-dispersion technique that could be easily deployed to determine farm-scale CH4 emissions (or any other trace gas). Such a technique would allow the quick evaluation of mitigation practices of entire farms and provide a means of understanding the variability that exists between different farm types.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Measurements
Three open-path lasers (GasFinder2; Boreal Laser, Spruce Grove, AB, Canada), tuned for detecting average CH4 concentration along a path (1 m above the ground), were situated 75 to 100 m away from the nearest farm building (Fig. 1). Lasers were aligned to enclose the north, east, and south perimeters of the farm; an array of retroreflectors terminated the laser path. When wind direction was westerly, the laser along the north perimeter was turned to a second array of retroreflectors located 200 m to the north of the laser (Fig. 1). In this manner, with minimum effort, there was always one laser recording concentration of the target gas upwind of the farm (background concentration) and one recording downwind (plume) concentration. The average concentration along the laser path was scanned approximately every 5 s and averaged every 15 min. Laser correction factors were determined repeatedly throughout the experiment. When possible we adjusted each laser to match our gas chromatograph (GC) measured CH4 concentrations (see below), if the wind direction indicated both the laser and GC samples measured the same concentration (i.e., background). These corrections varied from 0.83 to 0.97. These corrections were needed to remove bias instrument error that was evident when the three lasers were set side-by-side (2-m separation) in the field with no obvious upwind source. This side-by-side comparison resulted in average CH4 concentrations of 1488 (±93), 1560 (±99), and 1593 (±112) µg m3, which were significantly different from each other (P < 0.05). Accuracy analysis (Allen and Raktoe, 1981) was used to decompose the differences. On average, this analysis attributed 80% to bias (offset), 13% to random, and 7% to regression error. Assuming nothing changed for the on-farm study, the corrections to the lasers using the GC data removed the largest (bias) error. According to the laser specifications, the accuracy of the laser is 661.2 µg m3 per m (at 20° and 101 KPa), or 3.3 and 1.8 µg m3 over path distances of 200 to 360 m used in our study, respectively. However, our experience suggests that when working with lasers with different histories, and at the low concentrations associated with on-farm CH4 plumes, it is necessary to revisit the between-laser precision for each new study. We estimate in our case that maximum differences between lasers from the side-by-side comparison were 79 µg m3, but during our on-farm study these differences were most likely closer to 20 µg m3 after correcting for bias using the GC data as a common reference.
The average alongwind horizontal velocity (U); average wind direction (dir); the standard deviation of alongwind, acrosswind, and vertical velocity fluctuations (
u,
v, and
w, respectively); friction velocity (u*); surface roughness (zo); and the MoninObukhov stability length (L) were calculated from data collected at 10-Hz sampling interval (CR23X datalogger; Campbell Scientific, Logan, UT) using a three-dimensional sonic anemometer (CSAT3; Campbell Scientific). These statistics were calculated following Flesch et al. (2004), using the same 15-min averaging period as for gas concentrations. The anemometer was located at the laser locations in the southeast corner of the farm (Fig. 1). At this location, the influence of the farm buildings was deemed to be far enough away to have minimal impact on the sonic data.
Additional line-averaged concentrations of CH4 and sulfur hexafluoride (SF6) were measured downwind of the farm by drawing air through regulated intakes (using flow-holder valves) every 20 m along a long polyethylene tube (0.6-cm i.d.) at a 1-m height. On 18 November the sample line was lengthened from 100 m (5 intakes) to 200 m (10 intakes) long. A 12-vdc pump (TD3LS7; Brailsford and Company, Rye, NY) was connected to one end of the sample line. The flow rate through each intake was adjusted to be equal and approximately 2 L min1. A 2-L polyethylene bottle was connected to the exhaust of the pump to provide a 1-min time filter. A 20-mL air sample was taken every 5 min using a syringe and septum port. The content of the syringe was immediately discharged into an evacuated vial (6.8 mL volume) through a double septum cap (follows from Rochette and Bertrand, 2003), and the time and date recorded on the vial. Every 30 min, the background concentration upwind of the farm was sampled using three air samples drawn at 1.5-m height with a 20-mL syringe (after purging the syringe three times between sampling). At the end of each day, the vials were taken to the laboratory for analysis of CH4 and SF6 concentrations on a GC the following day (within 24 h of sampling). The over pressured vial provided a sufficiently large air volume to pass through the electron capture detector (ECD) column in the GC.
Before the start of the experiment, two lines of high-density polyethylene tubing (6-mm i.d.) were laid out to allow independent release of SF6 and CH4 from 10 release points throughout the barn and pen areas (Fig. 1). The release lines where connected to mass-flow controllers and cylinder tanks of ultrapure SF6 and CH4. The gas flow at each release point was adjusted using flow hold valves to be one-tenth of the gas flow set by mass-flow controller, where the mass-flow controller was 0.2 L min1 for SF6 (GFCS-011806; Aalborg Instruments and Controls, Orangeburg, NY) and 54.2 L min1 for CH4 (GFCS-011246; Aalborg Instruments and Controls). For each gas, three release points (20 m from each end and at the middle) were located 1 m above the peak of the roof. There were also three release points on the east side of the barn along the ventilation openings (aligned with roof locations). Four release points (2-m height) were spaced along the south perimeter of the open pen that housed the pregnant, nonlactating cows. The mass-flow controller units were checked in the laboratory by releasing gas over a long period of time and measuring the mass loss of the cylinder.
Emission Calculations
We used the inverse-dispersion technique outlined by Flesch et al. (2004). A Lagrangian stochastic (LS) model predicts the theoretical ratio of the targeted gas concentration to the emission rate (C/F)sim. An emission rate (FLS; g s1) is then determined after measuring the target gas concentration in the plume (Cobs; g m3):
![]() | [1] |
The commercially available software product WindTrax (Thunderbeach Scientific, 2005) was used in this calculation. A global positioning system was used to mark the farm structures, our gas release points (10 points), and the laser paths around the farm. WindTrax determined (C/F)sim in one of two modes, depending on how we modeled the emission source. Treated as 10 point sources (the actual configuration), WindTrax calculated 5000 trajectories downwind of each emission point (forward mode). Alternatively, when we treated the farm as two surface area sources, WindTrax calculated 10 000 trajectories upwind of the laser path (backward mode). These sets of trajectories define (C/F)sim.
The WindTrax model assumes idealized wind conditions appropriate for simple undisturbed terrain (using MoninObukhov wind relationships). However, in our farm setting the buildings and trees created wind complexity. Another difficulty in calculating (C/F)sim was that we had to assume an emission source configuration (i.e., location and relative proportion of emissions). While we knew the "synthetic" SF6 and CH4 release configuration, we did not know the details of the "natural" CH4 emissions; this is needed to separate the synthetic from the total emissions. Considering these complications, Flesch et al. (2005) concluded, "...with careful selection of a measurement location, an idealized inverse-dispersion technique can give the emissions from a farm complex within an error of ±10%...." They list three suggestions for dealing with wind and source complexity, and making accurate measurements:
Our site meets the first condition, but not the last two. At times our concentration measurements were less than 10 h and 1 xs from the farm. Although the measurement location was closer than recommended, it was chosen to ensure a detectable concentration rise for the CH4 lasers.
Sampling Periods
Our analysis scheme was based on 15-min sampling periods ending at 15, 30, 45, and 60 min past the hour (with corresponding wind measurements). The CH4 gas was turned on for 35 min, starting 5 min before a 15-min period, to allow the downwind concentrations to reach equilibrium before the measurement began. At the end of the 35 min the CH4 gas was shut off for 25 min. This cycle was repeated each hour (Fig. 2
). This onoff release pattern for CH4 was adopted to separate the farm CH4 emission from the known synthetic emission. After gas release stopped at the end of a 15-min period, we waited for concentrations to fall back to an equilibrium level (typically 5 min) before calculating concentration. From this pattern we created two "gas-on" and two "gas-off" observations each hour. Assuming the natural CH4 emissions were constant during the hour, the difference between the gas-on and gas-off emissions was taken as our estimate of the synthetic gas release rate (and should equal the actual 54.2 L min1 release rate). Because there were no natural SF6 emissions from the farm, it was not necessary to cycle between gas-on and gas-off conditions. Therefore, the SF6 gas was left on for the complete measurement period.
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
Gas Recovery Using the Actual Source Configuration
In the following discussion we judge the success of our inverse-dispersion calculations in terms of recovery rate, which expresses the inferred emission rate as a percentage of the known release rate. When we treat our farm emissions as coming from 10 point sources, using the actual outlet locations (including the height of release) in WindTrax, the average recovery rate for CH4 was 86% (Table 1). More impressive was the 91, 105, and 89% recovery on three of the four days. The cause of the larger error of 59% recovery rate on 1 December was not identified. In terms of wind conditions, there was nothing dramatically different on this day. On average the SF6 results were better, with a perfect recovery of 100% (±29%; Table 1). The better recovery for SF6 compared with CH4 was expected, since we did not have to deal with confounding natural emissions.
|
The perfect average recovery of SF6 hides a large day-to-day variability, from 148% on 9 November to 74% on 1 December. This large variability (compared with CH4) may be partially due to the short sampling line. While the CH4 concentration path length was between 250 and 360 m, and positioned so that the downwind laser usually captured the full width of the farm emission plume, the SF6 line was only 100 to 200 m long, which was not long enough to capture the full plume. This increases the uncertainty in our calculations, as the SF6 concentration along the sampling line would be sensitive to lateral dispersion from the farm (whereas if the sample line spans the plume it is insensitive to lateral dispersion). Likewise, our calculated FLS would be similarly sensitive to the modeled dispersion. This is an important consideration, as modeling horizontal dispersion is more error prone than modeling vertical dispersion (Flesch and Wilson, 2005).
Inadequate sampling may also cause some of the SF6 uncertainty. From Fig. 3 we can see the strong temporal variability in CH4 concentration, which also exists for SF6. Given this variability, our SF6 sampling frequency of once every 5 min may be insufficient to give an accurate 15-min average (to correspond with the 15-min wind statistics).
Gas Recovery Using a Generic Source Configuration
In theory, accurate inverse-dispersion calculations require that we model the actual emission source configuration. In many real-world situations we do not have the privilege of knowing these details. For example, we do not know the proportion of natural CH4 emissions coming from the barn and the dry stock pens. Flesch et al. (2005) discussed a solution to this problem: locate the concentration measurements where the modeled (C/F)sim will be insensitive to the source configuration details. Generally, this sensitivity is reduced as we move downwind. Our dataset allows us to consider the problem of unknown source configuration, and how it influences inverse-dispersion calculations.
As a hypothetical exercise, we applied the inverse-dispersion technique to our farm, but now assume a generic emission configuration. From a distance we might broadly identify two emission sources: the milking barn and the dry stock pens (i.e., where the cattle are). In WindTrax we thus treat the farm as two surface sources, corresponding to the barn and pens, and assume each has the same real emission rate. Although this does not match either the natural or synthetic source configuration, it has obvious similarities in terms of the relative spread of release, and is a rational choice given limited information. The question is how much this assumption affects the emission calculations.
We recalculated the CH4 and SF6 recoveries using WindTrax with the generic configuration (Table 2). Not surprisingly, there was a decline in accuracy. We found the average CH4 recovery was reduced from 86 to 66%, and SF6 recovery from 100 to 77%. From Flesch et al. (2005) we conclude that this undesired sensitivity is the result of concentration being measured too close to the farm.
|
|
| CONCLUSIONS |
|---|
|
|
|---|
Several factors were important to our successful application. First, the farm studied was relatively isolated on the landscape. This meant a clearly defined ambient wind regime, and locating our anemometer in this regime was straightforward. It also meant that we did not have other nearby emission sources (for CH4) to confound our concentration signals. A second factor was that we used relatively long line-average concentration measurements. The advantage of this should not be ignored. As argued by Flesch and Wilson (2005), utilizing crosswind average concentration reduces one of the major sources of uncertainty in the inverse-dispersion technique, specifically inaccurate modeling of lateral dispersion. Another factor in our success was ignoring periods known to be problematic for inverse-dispersion calculations (i.e., low winds, extreme stabilities, and measurements of downwind concentration that do not sample enough of the farm plume). One difficulty is the tracking of trajectories when wind directions are not consistent over the averaging time, common under low winds and very stable periods. Shortening the averaging time may avoid this particular problem in the inverse-dispersion calculation.
A fourth factor in our success was having a concentration measurement location that allowed us to ignore wind complexity. Even though we were as close as 75 m to the farm, this was evidently far enough (and much closer than recommended by Flesch et al., 2005). But to achieve accurate results at this location, it was necessary to use the actual source configuration in the dispersion calculations (i.e., 10 outlets). In many real-world situations these details will not be known, requiring assumptions about the source configuration. This will mean locating the measurement location further from the farm, keeping in mind that the location is a compromise between being far enough from the source to ensure the plume is uniform, yet close enough to ensure the target gas concentration is detectable above background concentration.
The selection of a proper concentration measurement location is probably the foremost issue in applying the inverse-dispersion technique. Concentration must be sufficiently high (relative to background) for accurate detection with the measuring devices. This creates a tension between wanting to be close to better detect concentration rise, and wanting to be far to minimize dispersion modeling errors caused by wind and source complexity. Our results highlight how the accuracy of the modeled source configuration can influence this distance decision. We conclude that the more one knows the source geometry and their relative strengths (e.g., the relative location, size, and type of animals, barn construction), and incorporating these to give an accurate source configuration in the dispersion model, the greater the likelihood of a successful application.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. M. McGinn, K. A. Beauchemin, T. K. Flesch, and T. Coates Performance of a Dispersion Model to Estimate Methane Loss from Cattle in Pens J. Environ. Qual., July 23, 2009; 38(5): 1796 - 1802. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. A. Harper, T. K. Flesch, J. M. Powell, W. K. Coblentz, W. E. Jokela, and N. P. Martin Ammonia emissions from dairy production in Wisconsin J Dairy Sci, May 1, 2009; 92(5): 2326 - 2337. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. McGinn, T. K. Flesch, B. P. Crenna, K. A. Beauchemin, and T. Coates Quantifying Ammonia Emissions from a Cattle Feedlot using a Dispersion Model J. Environ. Qual., October 16, 2007; 36(6): 1585 - 1590. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Soil Science Society of America Journal | Journal of Plant Registrations | The Plant Genome | |||