The gas emission rates of ammonia (NH3) and methane (CH4) from an
artificial source covering a surface area of 254 m2 were determined by
inverse dispersion modeling (IDM) from point-sampling and line-integrated
concentration measurements with closed- and open-path analyzers. Eight
controlled release experiments were conducted with different release rates
ranging from 3.8±0.21 to 17.4±0.4 mg s-1 and from 30.7±1.4 to 142.8±2.9 mg s-1 for NH3 and CH4,
respectively. The distance between the source and concentration measurement
positions ranged from 15 to 60 m. Our study consisted of more than 200
fluxes averaged over intervals of 10 or 15 min. The different releases
cover a range of different climate conditions: cold (<5∘C), temperate (<13∘C), and warm (<18∘C). As the average of all releases with all instrument types, the CH4 recovery rate QbLS/Q was 0.95±0.08 (n=19).
There was much more variation in the recovery of NH3, with an average
of 0.66±0.15 (n=10) for all the releases with the
line-integrated system. However, with an improved sampling line placed close to the source an average recovery rate of 0.82±0.05 (n=3) was
obtained for NH3. Under comparable conditions, the recovery rate
obtained with an open-path analyzer was 0.91±0.07 (n=3). The
effects of measurement distance, physical properties of the sampling line,
and deposition are discussed.
Introduction
The global agricultural system is currently facing one of its biggest humanitarian challenges: feeding the world's rising population while
preserving the environment and climate for future generations (FAO, 2017). The agricultural sector is a major contributor to global greenhouse gas (GHG) emissions (15 %) and ammonia (NH3) emissions (64 %) (OECD and FAO, 2019), leading to air pollution, climate change, deforestation, and loss of biodiversity (Aneja et al., 2009).
The European Union has established a reduction target for 2030 to reduce the
GHG emissions by at least 55 % (EEA, 2019), compared to
1990, and NH3 emissions by 19 % (NEC Directive 2016/2284/EU, 2016), compared to 2005. Agriculture must contribute to GHG emission reductions, and valid estimates of GHG emissions are important for national inventory regulation strategies and for selecting efficient mitigation techniques.
Choosing the appropriate methodology to quantify gaseous emissions can be a
challenge. In particular agricultural sources are challenging, as the sources
often are small and inhomogeneous, exhibit non-steady emissions over time
(e.g., NH3 emissions after slurry application; Hafner, 2018), and
are influenced by other sources in close vicinity. Most of the methodologies
have restrictions on the measurement location and/or the source and involve
complex instrumentation setup (e.g., fast-response analyzers, measurements
at multiple heights). The micrometeorological mass balance (MMB) method
(Desjardins et al., 2004) requires measuring concentration at multiple positions several meters above the ground, which is a challenge for obtaining high time resolution, and it ignores the horizontal turbulent transport (Hu et al., 2014). The tracer flux ratio method (TRM), which has also been used to
measure agricultural emissions (Vechi et al., 2022; Fredenslund et al., 2019; Delre et al., 2018), is a relatively labor-intensive and costly method typically with short intense measurement periods. In the case of dynamic emissions, this is not sufficient for resolving the temporal variations in emissions over days or weeks.
The inverse dispersion method (IDM) based on backward Lagrangian Stochastic
(bLS) dispersion modeling (e.g., Flesch et al., 2004, 1995) has been widely used for the assessment of NH3 and methane (CH4) emissions from many agricultural sources: dairy housing (Bühler et al., 2021; VanderZaag et al., 2014; Harper et al., 2009), cattle feedlot (McGinn et al., 2019; Todd et al., 2011; Van Haarlem et al., 2008; Flesch et al., 2007; McGinn et al., 2007), application of liquid animal manure (Kamp et al., 2021; Carozzi et al., 2013; Sintermann et al., 2011; Sanz et al., 2010), grazed pasture (McGinn et al., 2011; Voglmeier et al., 2018), rice field (Yang et al., 2019), lagoon (Ro et al., 2014; Wilson et al., 2001), composting stockpiles (Sommer et al., 2004), agricultural biodigester (Baldé et
al., 2016b; Flesch et al., 2011), farm (Flesch et al., 2005), and stored liquid manure (Lemes et al., 2022; Baldé et al., 2016a; Grant et al., 2015; McGinn et al., 2008).
IDM has been tested in controlled release experiments with different
conditions: ground-level source without obstacles (Flesch
et al., 2014; McBain and Desjardins, 2005; Flesch et al., 2004), ground-level source surrounded by a fence (Flesch et al., 2005; McBain and Desjardins, 2005), elevated source (Gao et al., 2008; McBain and Desjardins, 2005), and multiple emission sources (Hu et al., 2016; Ro et al., 2011; Gao et al., 2008) and to quantify the effect of NH3 deposition (Häni
et al., 2018).
IDM is a function of the geometry and location of source and downwind
concentration sensor (including height for the sensor) and the turbulence
characteristics in the surface layer. The statistical properties of the flow
in the atmospheric surface layer for the IDM are defined by the friction
velocity (u∗), roughness length (z0), Obukhov length (L),
and wind direction (Flesch et al., 2004). Emissions are derived from concentration measurements up- and downwind of the source, which could be determined with point or line-integrated measurements from closed- or open-path analyzers. IDM assumes an ideal atmospheric surface layer, which means (i) a horizontally homogeneous and flat surface, (ii) homogeneity and quasi-stationarity with respect to the turbulence characteristics, and (iii) spatially uniform emissions from a confined source (Flesch et al., 2004). Therefore, there should not be any obstacles (e.g., trees, buildings) in close vicinity to the source to fulfill the required IDM assumptions. Additionally, IDM has the limitation that there should not be any other sources of the same gas species that affects up- and downwind concentration differently. The IDM is simple, flexible (Harper et al., 2011), and robust even in non-ideal conditions and has a reported accuracy of 100±10 % when it is properly used (e.g., place of instruments, filtering
criteria) (Harper et al., 2010). Moreover, IDM is a direct measurement method that does not alter the physical properties of the source, and it is applicable for both small and large emissions of any shape of sources (Flesch et al., 2004) as opposed to indirect enclosure methods (e.g., chamber measurements).
Concentration measurements are mostly done with an open-path optical system
(e.g., Baldé et al., 2018; Bühler et al., 2021), because long path lengths (>50 m) enable a higher emission plume coverage and avoids
internal surfaces (e.g., tubes, pumps) where NH3 can adsorb (Shah et al., 2006; Vaittinen et al., 2014). However, open path has a limitation
on low concentration measurements (<10 ppb for CH4 and NH3) (Bai et al., 2022) and requires complex calibrations to reduce the uncertainty of the measurements (Häni et al., 2021; DeBruyn et al., 2020). In addition, it requires intensive labor to move and optically align the instruments to different positions depending on the predominant wind direction. Commercially available open-path analyzers exhibit limitations with respect to acceptable detection limits (Häni et al., 2021). Closed-path analyzers have rarely been used together with the IDM (Ro et al., 2011) due to its limitation caused by the adsorption of NH3 in the system. In addition, closed-path analyzers have only been used for point measurements, which challenges the ability to catch the emission plume and makes it sensitive to wind direction accuracy.
Data filtering is needed to ensure accuracy of the IDM, which is related to
the meteorological conditions (e.g., wind speed, atmospheric stability) and
wind direction. The quality criteria for filtering are based on the
atmospheric conditions in a measurement interval to ensure the assumptions
of the model are adequately met, which also lower the uncertainty of the
resulting data. Different criteria have been used in previous studies: Flesch et al. (2005) recommend to remove data where u*<0.15 m s-1, L<10 m, and z0>1 m, whereas McBain and Desjardins (2005) recommend u*<0.19 m s-1, L≤
3 m, and z0>1 m. Flesch et al. (2014) suggest the filtering criteria for the night of u*<0.05 m s-1 and the gradient between measured and MO-calculated temperature (ΔΔTthres)=0.05 K. Bühler et al. (2021) removed data where u*<0.05 m s-1, L<2 m, z0>0.1 m, standard deviation of the horizontal wind components (u,v) divided by u*(σu,v/u*)>4.5, and Kolmogorov constant (C0) >10.
This study aimed to assess the applicability and performance of a
closed-path analyzer used with a sampling system that allows for line-integrated concentration measurements used with the IDM for determining
emission rates of CH4 and NH3. This novel measuring system will
allow for measuring emissions from sources with low emission rates and will
have good flexibility for moving it around the source depending on the wind
direction in order to increase the probability of catching the emission
plume. This novel method is assessed by eight controlled releases of
CH4 and NH3 combined with up- and downwind measurements in
different positions using point and line-averaged concentration provided with
closed- and open-path analyzers. The use of CH4 and NH3 and open-
and closed-path systems to measure concentration will give us an opportunity
to (i) test the system of the line-averaged concentration measurement with a
closed-path analyzer and (ii) evaluate potential loss of NH3 downwind
from the source by deposition and/or gas-to-particle conversion, processes
that will not occur for inert CH4. This controlled-release study is the
first to compare the performances of open-path and line-integrated
closed-path systems for measuring emissions of NH3 and CH4.
Material and methodsSite descriptions
From November 2019 to March 2022, eight controlled release experiments were
performed at different grassland sites under varying conditions (see
Table 1). Five releases (I-DK to IV-DK and VIII-DK)
took place at Aarhus University campus in Viborg, Denmark, on two different fields
(56∘29′34.5′′ N, 9∘34′28.3′′ E and 56∘29′36.4′′ N, 9∘34′15.9′′ E). Three releases (V-CH to VII-CH) were performed at Bern University of Applied Sciences, Switzerland (46∘59′35.1′′ N, 7∘27′43.1′′ E). At all sites, the terrain was horizontally flat, and the height of the canopy varied between 15 and 25 cm for the different experiments. Obstacles upwind of the artificial source were more than 100 m away in all experiments. There were no significant sources near the experiment sites.
Date, gas cylinder description, ammonia and methane
release rate (RR), source and canopy height, downwind distance from source
to instruments, type of system attached to the cavity ring-down spectroscopy
(CRDS), and instrumentation of each controlled-release experiment (CRE).
G2103, G2202-i, G4301, and G2508 are different CRDS models, and GF corresponds to
GasFinder and MD to miniDOAS. Times are GMT+1.
CREDateGas cylinder NH3 RRCH4 RRSourceCanopyDistance fromSystemInstrumentsheightheightsource edgewithCRDSContent(bar)Total(mg s-1)(mg s-1)(cm)(cm)(m)I-DK29 November 20195 % NH3 and6214.6±0.3–02050Point 40∘2 G2103I-DK11:50–12:5095 % N2±2 %aII-DK29 November 201999 % CH4 and621–30.7±1.402050Point 40∘G2201-i andII-DK14:00–14:301 % N2±2 %aG4301III-DK12 October 20205 % NH3 and6213.8±0.2168.7±3.702535 to 6016 m line at 40∘G2103, G4301,III-DK11:45–15:1595 % CH4±2 %a(Line 1)and G2508IV-DK20 July 202110 % NH3 and62217.4±0.4142.8±2.9501815 to 3012 m line at 40∘G2103, G4301,IV-DK10:30–16:0090 % CH4±2 %a(Line 2)and G2508V-CH9 October 202110 % NH3 and27215.2±0.3128.9±2.701515 to 3016 m line at 40∘G4301, G2508,V-CH10:00–12:1090 % CH4±2 %bto 60(Line 1)3 GF, and 3 MDVI-CH9 October 202110 % NH3 and27413.2±0.3111.8±2.201515 to 3016 m line at 40∘G4301, G2508,VI-CH14:20–16:5090 % CH4±2 %bto 60(Line 1)3 GF, and 3 MDVII-CH9 October 202127413.2±0.3111.8±2.2501517:20–17:5010 % NH3 and15 to 3016 m line at 40∘G4301, G2508,11 October 202190 % CH4±2 %bto 60(Line 1)3 GF and 3 MD15:10–16:20VIII-DK22 April 202210 % NH3 and62214.5±0.3118.9±2.85071512 m line at 40∘3 G250812:30–15:0090 % CH4±2 %a(Line 2) and12 m line at 80∘with heatedinlets (Line 3)
a Air Liquide, Horsens, Denmark. b Carbagas, Bern, Switzerland.
Instrumentation
In this study, different models of cavity ring-down spectroscopy (CRDS)
analyzers from Picarro (Picarro Inc., Santa Clara, CA, USA) were used to
measure up- and downwind NH3 and CH4 concentration
(Table 1). Model G2201-i and model G4301 measure CH4 concentration, G2103 measures NH3 concentration, and G2509 measures CH4 and NH3 simultaneously. The CRDS is a closed-path analyzer with continuous absorption that measures concentrations at approximately 0.5 Hz. The CRDS analyzer consists of a laser and an optical cavity chamber with highly reflective mirrors, which give an effective path length of several kilometers. The light is absorbed in the cavity, and the decay of light intensity is called the ring-down time, which is directly related to the concentration of the specific compound. It has been frequently used to study agricultural emissions (e.g., Kamp et al., 2021; Pedersen et al., 2020; Kamp et al., 2019; Sintermann et al., 2011). Calibration of these CRDS instruments is conducted using a certified standard gas and a dilution system with NH3-free air before each release. Mass flow controllers (Bronkhorst EL-FLOW, Ruurlo, the Netherlands) were used to obtain a range of desired concentrations in all calibrations. The standard gas contained 10±0.31 ppm NH3 (Air Liquide, Horsens, Denmark) and 100±2 ppm CH4 (Air Liquide, Horsens, Denmark). Calibration showed high linearity of the instruments with R2=1. The CRDS instruments used pairwise for upwind and downwind measurements agreed within ±5 %, and no correction of the instruments was therefore done; see Figs. S5 and S6 in the Supplement.
In experiments V-CH to VII-CH, the downwind CH4 concentration was
measured with three GasFinder3 analyzers (GF3, Boreal Laser Inc., Edmonton,
Canada) and the downwind NH3 concentration with three miniDOAS
instruments (Sintermann et al., 2016). The GF3 analyzer is an open-path tunable diode laser device that measures line-integrated CH4 concentrations over path lengths of 5 to 500 m (i.e., single path length between sensor and retroreflector) with a temporal resolution of 0.3 to 1 Hz. The retroreflectors used in the experiments were
equipped with seven corner cubes, suitable for path lengths around 50 m. The
GasFinder devices have been widely used to measure emissions from different
types of agricultural sources with the IDM (Bühler et al., 2021; McGinn et al., 2019; VanderZaag et al., 2014; Harper et al., 2010; Flesch et al., 2007). The performance of the GF3 instruments is discussed in detail by Häni et al. (2021). The GF3 instruments were intercompared with the calibrated CRDS instrument by measuring ambient concentrations over at least 1 d and corrected accordingly. The applied corrections were c=cGF+0.14, c=cGF+0.22, and c=cGF+0.12 for the GF3 placed 15, 30, and 60 m from the source, respectively.
The miniDOAS instrument is an open-path device that measures NH3, NO,
and SO2 in the UV region between 190 and 230 nm based on the
differential optical absorption spectroscopy (DOAS; Platt and
Stutz, 2008) technique. It provides path-averaged concentrations for path
lengths between 15 and 50 m, with around 10 to 20 scans per second
averaged over 1 min. Ammonia emissions from agricultural sources (Kamp et al., 2021; Kupper et al., 2021; Voglmeier et al., 2018) and from an artificial source (Häni et al., 2018) were measured with miniDOAS analyzers. Further details on the instrument is given in Sintermann et al. (2016). The miniDOAS instruments were intercompared with the calibrated CRDS instrument by measuring ambient concentrations over at least 1 d. The miniDOAS offset concentration from the reference period 8 October 2021 from 15:30 to 17:00 (the time zone for all instances in the text is GMT+1) was added (3.2 µgm-3).
Gas release from an artificial source
The artificial source area had a gas distributor unit at the center, and
eight 1/4 in. polytetrafluoroethylene (PTFE) tubes leave the distributor to get a circular shape of the source area. Each tube contained three critical
orifices (100 µm diameter, stainless steel, Lenox Laser, USA) in
series with a 3 m distance between them. In total, the 24 orifices covered a
circular area of 254 m2.
Gas was released from a gas cylinder, and the flow was controlled with a mass
flow controller (in Denmark: Bronkhorst EL-FLOW; Ruurlo, the Netherlands; in
Switzerland: red-y smart controller, Vögtlin Instruments GmbH, Aesch,
Switzerland). The source height, the content of the gas cylinders, and the
release rate for each experiment are given in Table 1.
Setup
In the upwind position of all the experiments and in the downwind position
of the I-DK and II-DK experiments, the CRDS measured the concentration from a
single point 1.5 m above ground through a polytetrafluoroethylene (PTFE)
tube that was insulated and heated to approximately 40∘C. In the rest of the experiments, the CRDS measured downwind concentration from a sampling line system of PTFE tubes (Lines 1 and 2) or polyvinylidene fluoride (PVDF) tubes (Line 3). All tubes were insulated and heated (40 ∘C or
80 ∘C). The difference between the point and line-integrated system
is the number of positions where the gas sample is taken from. The point
system had only one inlet, while the line-integrated had several. Three
different versions of the line-integrated system (line) were built and used
during this research. In the III-DK, V-DK, VI-CH, and VII-CH experiments, the
sampling line system consisted of a 16 m tube with nine inlets, 2 m between
each inlet (Line 1). In the IV-DK and VIII-DK experiments, the sampling lines
were 12 m long with seven inlets, 2 m between each inlet (Lines 2 and 3). The inlets are made of critical orifices (0.25 mm i.d. for I-DK to VII-CH and 0.5 mm ID for VIII-DK polyetheretherketone (PEEK)) that guarantee
uniform flow through each inlet (Lines 1 to 3). In the VIII-DK
experiment, the sampling line system including the inlets was heated to
80 ∘C (Line 3).
Figure 1 shows the position of the source area relative to the sampling position, and the arrow indicates the wind direction during the experiments. The downwind concentrations were measured in one, two, or three distance (Table 1). In the V-CH, VI-CH, and VII-CH, downwind concentrations were measured at the same time at 15, 30, and 60 m distance from the edge of the source with multiple GF3 and miniDOAS instruments; one CRDS instrument was placed 15 m downwind (Fig. 1). The distance between the reflector and
the laser and detector of the GF3 and miniDOAS at the downwind position parallel
to the CRDS sampling line was also 16 m. For the other two downwind
positions the path lengths were 15 and 50 m, respectively. The height of
the measurement paths of all the open-path instruments was between 1.2 and
1.5 m. The background concentration of NH3 was stable with no sources
in close vicinity, thus in the three experiments, the average concentration
of each instrument 10 min before the release of each experiment was used as
the NH3 upwind concentration for the miniDOAS and the CRDS instruments.
In the V-CH, VI-CH, and VII-CH experiments, the measured NH3 background
concentration was 2.7 and 4.1 mg m-3 for the miniDOAS and 2.1 and 4.8 mg m-3 for the CRDS, respectively. The background concentration for V-CH and VI-CH was the same, since they were carried out on the same day (see Table 1).
Position of the orifices of the artificial source, ultrasonic
anemometer (sonic), and the concentration analyzer used in the eight
controlled-release experiments (CRE) of this study. Three types of analyzers
have been used: cavity ring-down spectrometer (CRDS), GasFinder (GF), and
miniDOAS (MD). The arrow indicates the wind direction during each
experiment.
In Denmark, the three wind components were measured at 16 Hz with a 3D
ultrasonic anemometer (WindMaster, Gill, Hampshire, UK) at a 1.5 and 1.7 m
height. In addition to concentration and wind, air temperature and atmospheric pressure were also measured. In Switzerland, the wind components
were measured at 20 Hz with a 3D ultrasonic anemometer (WindMaster, Gill,
Hampshire, UK) at a 2 m height. Air temperature and atmospheric pressure were
obtained from a meteorological station nearby the experiment site.
A global positioning system (in Denmark: GPS Trimble R10; Sunnyvale,
California, USA; in Switzerland: GPS Trimble Pro 6; Sunnyvale, California,
USA) was used to record the position of all instruments and the individual
critical orifices of the source.
Inverse dispersion method
The measured gas emission rates (Q) from the artificial source were
calculated in 15 min (experiments conducted in Denmark) or 10 min average
intervals (experiments conducted in Switzerland) using the R (R Core
Team, 2018) package bLSmodelR (https://github.com/ChHaeni/bLSmodelR, last access: 15 December 2022; version 4.3) as described by Häni et al. (2018). The simulation was performed with six million backward trajectories (N) and the source area defined as 24 individual circles of 5 cm radius as described by
Häni et al. (2018) with a high-performance computer cluster (PRIME –
Programming Rig for Modern Engineering, Aarhus University).
The emissions rate (Q) is proportional to the difference between measured
concentration downwind (Cdownwind) from the
source and the measured background concentration (Cupwind) and the dispersion factor (D):
Q=Cdownwind-CupwindD.
The dispersion factor (D) is calculated as
D=1N∑TDinside2wTD,
where N is the number of backward trajectories from the downwind analyzer
location. The summation refers to the trajectories touching inside the
source area (TDinside), taking the vertical velocity (wTD) at touchdown into account. The calculation of D includes determination of wind profiles and turbulence statistics that are based on the Monin–Obukhov similarity theory (MOST).
Surface deposition velocity
Ammonia is a relatively short-lived gas in the atmosphere and can either be
chemically converted or subjected to dry or wet deposition. Dry deposition
of NH3 is a complex mechanism that is restricted by both atmospheric
dispersion towards and uptake at surfaces (thus, it is correlated to several
different conditions indicated by, e.g., wind speed, solar radiation,
vegetation reactivity). The dry NH3 deposition rate is often expressed
with a deposition velocity (υd). In this study, we model dry deposition as a canopy resistance (“big leaf”) approach
where υd takes place unidirectionally, and it is calculated with the canopy resistances (Hicks et al., 1987):
υd=1Ra+Rb+Rc,
where Ra is the aerodynamic resistance, Rb is the quasi-laminar boundary resistance, and Rc is the bulk canopy resistance. Because Ra (the resistance between the concentration measurement height and the notional height z0) is already included in the bLS model calculations, Eq. (3) can be simplified to represent a surface deposition velocity:
υd*=1Rb+Rc.
For each model trajectory calculation, this surface deposition velocity is
acting on each touchdown outside the emitting source to provide individual
dry deposition fluxes Fd from prevailing touchdown concentration
CTD as
Fd=-υd*CTD.
This reduces the modeled trajectory concentration at each touchdown outside
the source by
ΔCd=CTD⋅exp-2⋅υd*wTD-1.
We refer to Häni et al. (2018) for the derivation of the above equation
and a thorough explanation of the implementation of the dry deposition
algorithm in the bLS model.
In this study, Rb was calculated according to Garland (1977), as a function of the roughness length (z0), friction velocity (u*), kinematic viscosity of air ν, and molecular diffusivity of NH3 in air (δNH3):
Rb=1.45z0u*ν0.24νδNH30.8u*.
Regarding Rc, it is related to the chemical characteristics of the studied gas and the characteristics of the leaf (e.g., type, size). There are different models to calculate Rc. Due to the complexity and the uncertainty of the determination of the resistance, Rc was calculated following the same procedure as by Häni et al. (2018) with the
bLSmodelR. It was assumed that QbLS/Q<1 was solely due to dry
deposition. A similar approach is used here, where 12 values of Rc from 0 to 500 s m-1 were tested in the bLS model that includes ammonia deposition to estimate the Rc, giving QbLS/Q=1 in all intervals. This was done with linear interpolation between the two points closest to Q/Q=1. Using this estimated Rc and the calculated Rb value for each interval, υd* was estimated for all intervals with all instruments. The υd* values are compared to previously reported values for NH3.
Another approach for calculating the Rc is with an empirical equation, which will be used for calculating values for υd*. These calculated values will be compared to the values obtained with the bLS model. It is assumed that Rc is
unidirectional and equal to the sum of the stomatal resistance Rs and the cuticular resistance Rw; see Eq. (6).
1Rc=1Rs+1Rw
The stomatal resistance Rs is calculated with equation Eq. (7)
(Wesely, 2007):
Rs=Rs(min)[1+200SR+0.1]2400Ts40-Ts,
where Rsmin is the minimum bulk canopy Rs for water vapor that is assumed to be equal to 250 s m-1 (Lynn and Carlson, 1990), SR is the solar radiation, and Ts is the soil temperature.
The cuticular resistance is calculated with Eq. (8) (Massad et al., 2010):
Rw=Rwminea(100-RH)e0.15TLAI0.5,
where Rwmin is the minimum cuticular resistance, a is an empirical factor, RH is the relative humidity, T is the air temperature, and LAI is the leaf area index. The parameters Rwmin (10 s m-1), a (0.110), and LAI (2 m2 m-2) were obtained from Massad et al. (2010), Table 1.
Results and discussionRecovered fractions of ammonia and methane
The accuracy of the bLS model is evaluated by the recovered NH3 and
CH4 fractions, QbLS/Q, and the standard deviation σQbLS/Q for all the releases without taking NH3 deposition into account. In all experiments except I-DK and II-DK (Table 1), NH3 and CH4 were released simultaneously. The use of these two gases gives us the additional opportunity to assess potential loss of NH3 downwind from the source by deposition or gas-to-particle conversion, processes that will not occur for CH4 due to its inertness. As the average of all releases and measurement systems, the CH4 recovery rate was 0.95±0.08 (n=19) (Fig. 4). This recovery is similar to 0.93±0.14 (n=8) observed by Gao et al. (2008) with a different controlled release configuration and ground-level sources. There was more variation in the recovery of NH3, with an average of 0.66±0.15 (n=10) for all the releases with the line-integrated system. However, the improved sampling lines (Lines 2 and 3) placed at 15 m downwind from the source had an average recovery of 0.82±0.05 (n=3) for NH3 (Fig. 2). Under comparable conditions, the NH3 recovery rate obtained with the miniDOAS (MD) was 0.91±0.07 (n=3). Häni et al. (2018) observed almost the same recovered fraction, 0.91±0.12, at 15 m from the edge of the source with the MD. The recovery rates of all experiments in this study are shown in Fig. 2, whereas climate conditions such as wind direction, friction velocity u∗, air
temperature, relative humidity (RH), soil temperature, and solar radiation
(SR) from each experiment are presented in Table S2 in the Supplement. In addition, the average of the recovery fraction rates of both gases and their relation in each release are shown in Table 2. I-DK and II-DK were conducted during cold conditions (∼5∘C) with RH ranging from 65 % to 71 %, whereas IV-DK and VIII-DK were conducted in warm conditions (14–18 ∘C) with RH between 48 % and 63 %. The other releases were conducted under moderate temperature conditions (10–13 ∘C) with RH between 39 % and 89 %.
Additional information on the atmospheric conditions, weather conditions,
and recovery fraction rates for each average time interval for each release
is shown in Table S1.
The recovered fractions QbLS/Q of ammonia and methane from each release and analyzer. The downwind distance from the source to the
analyzer is indicated in the legend. Line 1 had a length of 16 m, and it was
heated to 40 ∘C. Line 2 had the same temperature as Line 1, but it was 12 m long. Line 3 had the same length as Line 2 but was heated
to 80 ∘C. Error bars represent the standard deviation.
Average of the recovery fractions QbLS/Q of ammonia and
methane for each release and analyzer with standard deviation. Lines 1 to
3 are described in the Fig. 2 caption.
Three different CRDS sampling line systems have been used from III-DK to
VIII-DK. The difference between the lines was the length and the heating
temperature. Line 1 had a length of 16 m, and it was heated to 40 ∘C. Line 2 had the same temperature as Line 1, but it was
12 m long. Line 3 had the same length as Line 2 but was heated to 80 ∘C, and the critical orifices have a higher inflow than
Lines 1 and 2 (see Sect. 2.4 Setup). We expect that decreasing the
length and increasing the heating temperature of the line will improve
QbLS/Q for NH3 (no expected effect for CH4) by avoiding
adsorption and reducing the response time in the sampling line.
Line 1 was used with the source at ground level and elevated (Table 1), whereas the other two lines only with the source elevated. When the source was at ground level, Line 1 had a recovery rate ranging from 0.42±0.17 to 0.60±0.10 and from 0.75±0.13 to 1.01±0.05 for NH3 and CH4, respectively. The lowest and the highest NH3 recovery rate of Line 1 are directly related to the furthest (60 m) and the shortest (15 m) downwind distance measurement from the source. In addition, the standard deviation σQbLS/Q at the furthest position is higher than at the closest position, which is in accordance with the results from Häni et al. (2018). High uncertainty of the QbLS/Q is related to a smaller difference in concentration between downwind and background concentrations and due to smaller D values (Häni et al., 2018). This is also the reason for the low CH4 recovery rate of Line 1 in III-DK at 60 m (0.75±0.13); downwind concentration is only 4 %–10 % higher than upwind concentrations, since this is one of the lowest
CH4 release rates (Table 1). This is in line with Coates et al. (2021), who observed that the bLS model underestimated 49 % of CO2 released at 50 m fetch distance, partially because the measured downwind concentration was close to the background level. Therefore, in this study, the accuracy of QbLS is mainly influenced by the uncertainty of the concentration measurement, hence the downwind distance is limited by the properties of the gas analyzers and the size of the emission strength of the source. This means the system can be limited in use if the emission source has a large height and low emission strength where, as a rule of thumb, measurements should be conducted at a distance from the source at least 10 times the height of the source (Harper et al., 2011).
In VII-CH, Line 1 was used with the source elevated and had a recovery rate
of 0.69±0.12 for NH3 and 0.95±0.10 for CH4. Line 2
had a numerically higher recovery rate than Line 1, ranging from 0.76±0.08 to 0.81±0.16 and from 0.89±0.20 to 0.99±0.12 for NH3 and CH4, respectively, in IV-DK and VIII-DK. The
length of the line appears to affect the NH3 recovery rate; this might
be due to the increased surface area that NH3 can adsorb to, and there
is a lower flow in each of the critical orifices that decreases the response
time of the system (Shah et al., 2006; Vaittinen et al., 2014). Looking at the measured NH3 rates over time (Fig. 3), higher emissions are
reached with Line 3 for the first hour, indicating a faster time response
compared to Line 2. However, after an hour there was not a clear difference
between the lines. The results indicate that increasing the sampling line
temperature to 80 ∘C had a positive effect on the recovery, which
reached 87 % at a distance of 15 m. From the data obtained by the
open-path analyzer (MD), we can conclude that deposition can cause a
reduction in recovery in the order of 2 %–16 % (Fig. 2). Thus, the recovery obtained with the improved line (Line 3) approaches the recovery obtained with the open-path
analyzer. It should be noted that a direct comparison between Line 3 and the
open-path analyzer (MD) has not been made, and further improvement can still
be suggested for the CRDS sampling line system. Specifically, increasing the
flow through the sampling line and critical orifices will reduce NH3
adsorption in the tubing material. However, similar flow rates through the
sampling orifices in the sampling line must still be ensured.
Ammonia (NH3) fluxes measured with Line 2 and Line 3 in 10 min interval averages in VIII-DK.
The point CRDS system had a recovery rate of 0.70±0.22 and 1.01±0.05 for NH3 and CH4, respectively. The benefit of the
point CRDS system is mainly that increasing the flow in the tubing is less
limited, since there are no critical orifices for which equal flow must be
maintained. However, comparing point and line CRDS systems to the modeled
concentration distribution (Fig. 4), the line-integrated measurement system covers a larger part of the emission plume from the source in a higher wind direction range. In addition, a line-integrated measurement system can reduce uncertainty in the IDM (Flesch et al., 2004), since it is less sensitive to error in the measured wind direction. This is in accordance with Ro et al. (2011), who observed an almost double recovery value of a line-integrated
measurement system for CH4 compared to a point measurement system using
a photoacoustic gas monitor.
Contours of the modeled concentration distribution (CbLS)
for CRE I-DK, CRE III-DK, and CRE IV-DK.
Open-path measurement systems
The recovery rates for the GFs (CH4) ranged from 0.87±0.10 to
1.08±0.15. In V-CH to VII-CH, the corresponding standard deviation
σQbLS/Q of GF 15 m varies from 0.07 to 0.18, while Line 1 (placed parallel to GF 15 m) ranges from 0.06 to 0.10. These standard
deviations σQbLS/Q are comparable with those measured by Gao et al. (2009) (1.03±0.16).
In V-CH and VI-CH (source at ground), the MDs (NH3) had recovery rates
ranging from 0.57±0.12 to 0.93±0.03. In VII-CH, MDs exhibit
higher recoveries ranging from 0.64±0.09 to 0.98±0.10, since
the source was elevated. Generally, it is recommended to do a release
experiment above ground level to reduce the probability of deposition close
to the release area (McBain and Desjardins, 2005). As expected, the recovery rate decreased with downwind distance of the sampling position due to NH3 deposition, which will be evaluated in Sect. 3.4. Comparing MD at 15 m and Line 1 (placed in parallel) in V-CH to VII-CH (Fig. 2), the
recovery rates are higher for MD. The highest difference between MD and Line 1 was in V-CH, where there was the highest RH (87 %). However, there are
no clear patterns explaining the difference between emissions from the
different measurement systems based on atmospheric conditions (Fig. S2). Although the improved recovery with Line 2 (0.81±0.16) and Line 3 (0.87±0.11) in IV-DK and VIII-DK could be influenced by the warmer conditions and solar radiation (Table S2), it is plausible that the line improvements caused
the increase. An increased flow through the orifices and higher temperature
of the sampling line will lead to less NH3 adsorption, thereby getting a
better recovery from the release.
These results show the advantage of an open-path instrument compared to a
closed-path instrument to measure NH3 emissions
(Fig. 2), since open-path avoids prolonged
response caused by the adsorption of NH3 to sampling materials
(Shah et al., 2006; Vaittinen et al., 2014). However, quality assurance is more
challenging for open-path instruments, because the use of a closed gas cell
for calibration with a certified gas standard is tedious and means that the
instrument setup deviates significantly from field measurement
(DeBruyn et al., 2020). Intercomparison to a
reference method (as an alternative to gas standard calibration) may also
introduce uncertainties (Häni et al., 2021). In
addition, the closed-path system presented in this study (line CRDS) is more
flexible with respect to moving the sampling line around the source
depending on the predominant wind direction. This is particularly important
in areas with frequently changing wind direction.
Surface deposition velocity
The corresponding surface deposition velocities (υd*) required to have a recovery rate QbLS/Q=QCH4 are presented in Fig. 5. This approach assumes a recovery equals to the measured QCH4 for each of the measurement systems when taking deposition into account, which is not completely correct for closed-path sampling. In the following, therefore, we refer to deposition velocity required to achieve QbLS/Q=QCH4 as the apparent deposition velocity (υad*). This is included to provide
data on deposition velocities for ammonia for which data are currently very
limited. The recovery rates observed in Fig. 2 show that the MD performed best, whereas lower QbLS/Q were seen in the sampling lines, thus the lowest υad* is expected from MD.
Additional information of Rc and υad* for each time interval in each experiment is shown in Table S1. The apparent surface deposition velocities ranged from 0.2 to 1.8 cm s-1 for open-path data and from 0.2 to 4.2 cm s-1 for the line sampling.
Häni et al. (2018) reported υad* in the range from 0.3 to 1.1 cm s-1. In all the releases where downwind concentrations were measured at different positions, υad* appears to increase with distance. For example, in VI-CH, υad* is
0.3±0.1, 0.7±0.4, and 1.1±0.5 cm s-1 at 15, 30, and 60 m, respectively. The reason for this increase with distance is presently unclear and should be investigated further.
In V-CH, VI-CH, and VII-CH, υad* from Line 1 are 4.9, 4.3, and 4.3 times higher than MD at 15 m. As expected υad* was higher for Line 1, as the Q/Q was lower for Line 1 compared to MD in these experiments. Comparing the apparent deposition velocities from these experiments show comparable values for Lines 2 and 3 but higher values for Line 1, i.e., Line 1 (VII-CH) had υad* of 1.7±0.4, whereas Lines 2 and 3 (VIII-DK) had υad* of 1.4±0.4 and 1.2±0.6 cm s-1, respectively, when measuring 15 m from the elevated source. During VII-CH and VIII-DK the temperature differed with 1 ∘C, and the relative humidity was approximately the same, but wind speed and solar radiation differed (Table S2).
Corresponding apparent surface deposition velocities
(υad*) required to
have a recovery rate QbLS/Q closest to one in all the releases. The error bars represent the standard deviations. All values are shown in Table S1.
Many factors affect the deposition velocity, but it is possible to calculate
υad* from empirical models as explained previously (see Sect. 2.6). Figure 6 shows υad* for MDs in V-CH,
VII-CH, and VIII-CH compared to υad* calculated with the empirical models (Eqs. 3–8). Using the empirical models, υad* varies from 0.13 to 1.02 cm s-1, increasing with the relative humidity (87 %, 76 %, and 52 % RH in V-CH,
VI-CH, and VII-CH, respectively). The difference between the two ways of
estimating υad* is not surprising, since (i) bLS-derived deposition may be influenced by methodological uncertainties and therefore deviate from true deposition, and (ii) calculated resistances are associated with uncertainties due to estimations of physical parameters. In addition, an artificial source most likely will have higher υd* than what is expected from a typical agricultural source (Häni et al., 2018). This is expected, because an artificial source with discrete outlets located near the ground will have a larger deposition close to the release, because the ground has a high capacity for NH3 sorption. On the contrary, a constantly and homogeneously emitting source (e.g., storage tank or manure pile) is not expected to have any significant deposition within the source. This is seen with the higher υad* values found in these experiments compared to the calculated values with the empirical models. The height of the source might also have an influence on υad*. This is indicated by
the lowest υad* in VII-CH, where the source was elevated compared to V-CH and VI-CH, where the gas was released on the surface. Placing the source above ground level will reduce the obstacles (crop on the field) for gas dispersion, reducing surface
deposition. However, the bLS model does not consider the height of the
source. For example, evaluating emissions from the application of liquid
animal manure (ground-level source) or a dairy housing (elevated source)
will have different υad*.
Corresponding apparent surface deposition velocities
(υad*) required to have a recovery rate QbLS/Q closest to one for miniDOAS (MD) in release V-CH, VI-CH, and VII-CH and υad* calculated with an empirical model.
Uncertainties and sensitivity analysis
The line-integrated concentration measurement systems with a CRDS analyzer
(Lines 1 to 3; excluding one value at 60 m) had an average recovery rate of CH4 of 96±4 % (1 standard deviation, n=8) (see also Table 2). Based on this, it is concluded that the method is not biased with respect to CH4. The overall precision of CH4 concentration measurements observed for the three versions of the line ranged from 4 % to 22 % at a 15 m distance (see Table 2). The equivalent precision for NH3 concentration measurements was between 13 % and 23 % (Table 2) for the improved sampling lines (Lines 2 and 3) at a 15 m distance (CRE IV and CRE VIII). The 9 % difference between NH3 recovery rates for MD (CRE V-CH–VII-CH) and for Lines 2 to 3 at a 15 m distance is assessed to represent the sampling line adsorption bias related to the line-integrated system under the best conditions (improved sampling lines and moderate to warm temperatures). The most relevant factors affecting the uncertainties are the determination of the wind direction offset, canopy height, and downwind concentration analyzer height. Furthermore, the accuracy of the concentration measurements and data filtering criteria related to the meteorological conditions (e.g., Flesch et al., 2005) are other important factors that influence the uncertainty of the
measurement. It would require a comprehensive study to evaluate the
contributions of these parameters, which is not in the scope of this study.
Therefore, only sensitivity analyses of the wind direction offset and the
sensor height are included below. The individual uncertainty contributions
would be lumped into the overall precision and bias mentioned above.
A sensitivity analysis of two input parameters for the bLS model was based
on the resulting Q/Q ratio when changing the inlet height of the analyzer
and the wind direction offset compared to the valid measured values in
release VIII-DK. This was done for 11 fluxes in average intervals of 15 min,
where all emissions were estimated again with the bLS model. For the
assessment of the influence of the input for the measurement height all
other variables were kept constant. Likewise, for the influence of the wind
direction, all other variables were kept constant, while the wind direction
offset was changed. The results are presented in Fig. S4, where it can be seen that Q/Q was most sensitive to the changes in wind direction offset, stressing the importance of the true offset in wind direction. Therefore, the wind direction must be thoroughly evaluated for the accuracy of emission estimation, since more or less trajectories have touchdowns inside the source area for the dispersion factor (Eq. 2). In addition, the uncertainty of the Q/Q ratio increases as wind direction offset increases. The emission estimation accuracy from point systems is more sensitive to error in the measured wind direction (Flesch et al., 2004).
The accuracy and precision of the emission estimation also depend on the
detection limits of the concentration sensor analyzer, especially when the
downwind concentration is close to the background level, as it was shown
previously (see Sect. 3.2). Therefore, it is recommended to conduct
concentration and turbulence measurements not far from the source but a
minimum of 10 times the source height according to Harper et al. (2011) at a known height to reduce the uncertainty of the calculated emission rates.
Conclusion
For CH4, the average recovery rate for all releases was 0.95±0.08, which demonstrates that line-averaged measurement with a closed-path analyzer is comparable to an open-path system for inert gases. Under comparable conditions, average recovery rates of 0.82±0.05 (n=3) and 0.91±0.07 (n=3) for NH3 and 0.94±0.02 (n=3) and 0.93±0.05 (n=3) for CH4 were obtained with the closed-path and open-path line-integrated system, respectively. The implementation of the new method presented in this study
will enable measurement of fluxes of multiple gases from different types of
sources and evaluate the effects of mitigation strategies on emissions. In
addition, this method allows for continuous online measurements that resolve
temporal variation in NH3 emissions and the peak emissions of CH4.
It is stressed that the wind direction must be thoroughly evaluated for the
accuracy of emission estimation with the bLS model.
A significant fraction of the emitted NH3 is deposited near the source.
Consequently, including the deposition algorithm in the bLS model will have
less bias in the emission evaluation at ground-level sources (e.g.,
application of liquid animal manure) compared to elevated sources (e.g.,
slurry tank). The present study shows that the estimated deposition
velocities are in the same order of magnitude in all the releases with some
variation across the different approaches (instrument, distance, method).
Code and data availability
Code and data are available at 10.5281/zenodo.7695569 (YoLemes, 2023).
The supplement related to this article is available online at: https://doi.org/10.5194/amt-16-1295-2023-supplement.
Author contributions
YML: conceptualization, experiment design, conducting experiment, data analysis, validation, dispersion modeling, deposition modeling, writing (original draft), writing (review and editing), visualization.
CH: experiment design, conducting experiment, data analysis, validation, dispersion modeling, deposition modeling, writing (review and editing).
JNK: conceptualization, experiment design, validation, dispersion modeling, deposition modeling, writing (review and editing).
AF: conceptualization, experiment design, validation, funding, writing (review and editing).
All authors have read and agreed to the published version of the manuscript.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
Thanks is expressed to Simon Bowald for his great ideas and his help with designing and building up Line 3 and to Marcel Bühler for his valuable contribution in creating some figures for this paper. We also thank technicians Martin Häberli-Wyss, Peter Storegård Nielsen, Jens Kristian Kristensen, and Heidi Grønbæk for their invaluable help during the experimental part of the study. We thank the reviewers for valuable comments that helped to improve the article.
Financial support
This study was funded by the Ministry of the Environment and Food of Denmark as service agreement no. 2019-760-001136.
Review statement
This paper was edited by Daniela Famulari and reviewed by Joseph Pitt and one anonymous referee.
ReferencesAneja, V. P., Schlesinger, W. H., and Erisman, J. W.: Effects of Agriculture
upon the Air Quality and Climate: Research, Policy, and Regulations,
Environ. Sci. Technol., 43, 4234–4240, 10.1021/es8024403, 2009.Bai, M., Loh, Z., Griffith, D. W. T., Turner, D., Eckard, R., Edis, R., Denmead, O. T., Bryant, G. W., Paton-Walsh, C., Tonini, M., McGinn, S. M., and Chen, D.: Performance of open-path lasers and Fourier transform infrared spectroscopic systems in agriculture emissions research, Atmos. Meas. Tech., 15, 3593–3610, 10.5194/amt-15-3593-2022, 2022.Baldé, H., VanderZaag, A. C., Burtt, S., Evans, L., Wagner-Riddle, C.,
Desjardins, R. L., and MacDonald, J. D.: Measured versus modeled methane
emissions from separated liquid dairy manure show large model
underestimates, Agr. Ecosyst. Environ., 230, 261–270,
10.1016/j.agee.2016.06.016, 2016a.Baldé, H., VanderZaag, A. C., Burtt, S. D., Wagner-Riddle, C., Crolla,
A., Desjardins, R. L., and MacDonald, D. J.: Methane emissions from
digestate at an agricultural biogas plant, Bioresource Technol., 216,
914–922, 10.1016/j.biortech.2016.06.031, 2016b.Baldé, H., VanderZaag, A. C., Burtt, S. D., Wagner-Riddle, C., Evans,
L., Gordon, R., Desjardins, R. L., and MacDonald, J. D.: Ammonia emissions
from liquid manure storages are affected by anaerobic digestion and
solid-liquid separation, Agr. Forest Meteorol., 258, 80–88,
10.1016/j.agrformet.2018.01.036, 2018.Bühler, M., Häni, C., Ammann, C., Mohn, J., Neftel, A., Schrade, S.,
Zähner, M., Zeyer, K., Brönnimann, S., and Kupper, T.: Assessment of
the inverse dispersion method for the determination of methane emissions
from a dairy housing, Agr. Forest Meteorol., 307, 108501,
10.1016/j.agrformet.2021.108501, 2021.Carozzi, M., Loubet, B., Acutis, M., Rana, G., and Ferrara, R. M.: Inverse
dispersion modelling highlights the efficiency of slurry injection to reduce
ammonia losses by agriculture in the Po Valley (Italy), Agr. Forest Meteorol., 171–172, 306–318,
10.1016/j.agrformet.2012.12.012, 2013.Coates, T. W., Alam, M., Flesch, T. K., and Hernandez-Ramirez, G.: Field testing two flux footprint models, Atmos. Meas. Tech., 14, 7147–7152, 10.5194/amt-14-7147-2021, 2021.DeBruyn, Z. J., Wagner-Riddle, C., and VanderZaag, A.: Assessment of
Open-path Spectrometer Accuracy at Low Path-integrated Methane
Concentrations, Atmosphere, 11, 184, 10.3390/atmos11020184,
2020.Delre, A., Mønster, J., Samuelsson, J., Fredenslund, A. M., and Scheutz,
C.: Emission quantification using the tracer gas dispersion method: The
influence of instrument, tracer gas species and source simulation, Sci. Total Environ., 634, 59–66, 10.1016/j.scitotenv.2018.03.289, 2018.Desjardins, R. L., Denmead, O. T., Harper, L., McBain, M., Massé, D.,
and Kaharabata, S.: Evaluation of a micrometeorological mass balance method
employing an open-path laser for measuring methane emissions, Atmos. Environ., 38, 6855–6866, 10.1016/j.atmosenv.2004.09.008,
2004.EEA: EMEP/EEA air pollutant emission inventory guidebook 2019: technical
guidance to prepare national emission inventories, Publications Office, LU, 10.2800/293657, 2019.
FAO: Food and Agriculture Organization of the United Nations. The future of
food and agriculture: trends and challenges, Food and Agriculture
Organization of the United Nations, Rome, 163 pp., ISBN 978-92-5-109551-5, 2017.Flesch, T., Wilson, J., Harper, L., and Crenna, B.: Estimating gas emissions
from a farm with an inverse-dispersion technique, Atmos. Environ., 39, 4863–4874, 10.1016/j.atmosenv.2005.04.032, 2005.Flesch, T. K., Wilson, J. D., and Yee, E.: Backward-Time Lagrangian
Stochastic Dispersion Model and Their Application to Estimate Gaseous
Emissions, J. Appl. Meteorol., 34, 1320–1332, 10.1175/1520-0450(1995)034<1320:BTLSDM>2.0.CO;2, 1995.Flesch, T. K., Wilson, J. D., Harper, L. A., Crenna, B. P., and Sharpe, R. R.: Deducing Ground-to-Air Emissions from Observed Trace Gas Concentrations:
A Field Trial, J. Appl. Meteorol., 43, 487–502, 10.1175/JAM2214.1, 2004.Flesch, T. K., Wilson, J. D., Harper, L. A., Todd, R. W., and Cole, N. A.:
Determining ammonia emissions from a cattle feedlot with an inverse
dispersion technique, Agr. Forest Meteorol., 144, 139–155,
10.1016/j.agrformet.2007.02.006, 2007.Flesch, T. K., Desjardins, R. L., and Worth, D.: Fugitive methane emissions
from an agricultural biodigester, Biomass Bioenerg., 35, 3927–3935,
10.1016/j.biombioe.2011.06.009, 2011.Flesch, T. K., McGinn, S. M., Chen, D., Wilson, J. D., and Desjardins, R. L.: Data filtering for inverse dispersion emission calculations, Agr. Forest Meteorol., 198–199, 1–6, 10.1016/j.agrformet.2014.07.010, 2014.Fredenslund, A. M., Rees-White, T. C., Beaven, R. P., Delre, A., Finlayson,
A., Helmore, J., Allen, G., and Scheutz, C.: Validation and error assessment
of the mobile tracer gas dispersion method for measurement of fugitive
emissions from area sources, Waste Manage., 83, 68–78,
10.1016/j.wasman.2018.10.036, 2019.Gao, Z., Desjardins, R. L., van Haarlem, R. P., and Flesch, T. K.:
Estimating Gas Emissions from Multiple Sources Using a Backward Lagrangian
Stochastic Model, J. Air Waste Manage., 58, 1415–1421, 10.3155/1047-3289.58.11.1415, 2008.Gao, Z., Desjardins, R. L., and Flesch, T. K.: Comparison of a simplified
micrometeorological mass difference technique and an inverse dispersion
technique for estimating methane emissions from small area sources,
Agr. Forest Meteorol., 149, 891–898, 10.1016/j.agrformet.2008.11.005, 2009.Garland, J. A.: The dry deposition of sulphur dioxide to land and water
surfaces, P. Roy. Soc. Lond. A-Mat., 354, 245–268,
10.1098/rspa.1977.0066, 1977.Grant, R. H., Boehm, M. T., and Bogan, B. W.: Methane and carbon dioxide
emissions from manure storage facilities at two free-stall dairies,
Agr. Forest Meteorol., 213, 102–113,
10.1016/j.agrformet.2015.06.008, 2015.
Hafner, S. D.: The ALFAM2 database on ammonia emission from field-applied
manure: Description and illustrative analysis, Agr. Forest Meteorol., 258, 66–79, 2018.Häni, C., Flechard, C., Neftel, A., Sintermann, J., and Kupper, T.: Accounting for Field-Scale Dry Deposition in Backward Lagrangian Stochastic Dispersion Modelling of NH3 Emissions, Atmosphere, 9, 146, 10.3390/atmos9040146, 2018.Häni, C., Bühler, M., Neftel, A., Ammann, C., and Kupper, T.: Performance of open-path GasFinder3 devices for CH4 concentration measurements close to ambient levels, Atmos. Meas. Tech., 14, 1733–1741, 10.5194/amt-14-1733-2021, 2021.Harper, L. A., Flesch, T. K., Powell, J. M., Coblentz, W. K., Jokela, W. E.,
and Martin, N. P.: Ammonia emissions from dairy production in Wisconsin,
J. Dairy Sci., 92, 2326–2337, 10.3168/jds.2008-1753, 2009.Harper, L. A., Flesch, T. K., Weaver, K. H., and Wilson, J. D.: The Effect
of Biofuel Production on Swine Farm Methane and Ammonia Emissions,
J. Environ. Qual., 39, 1984–1992, 10.2134/jeq2010.0172,
2010.Harper, L. A., Denmead, O. T., and Flesch, T. K.: Micrometeorological
techniques for measurement of enteric greenhouse gas emissions, Anim. Feed
Sci. Tech., 166–167, 227–239, 10.1016/j.anifeedsci.2011.04.013, 2011.Hicks, B. B., Baldocchi, D. D., Meyers, T. P., Hosker, R. P., and Matt, D.
R.: A preliminary multiple resistance routine for deriving dry deposition
velocities from measured quantities, Water Air Soil Poll., 36, 311–330,
10.1007/BF00229675, 1987.Hu, E., Babcock, E. L., Bialkowski, S. E., Jones, S. B., and Tuller, M.:
Methods and Techniques for Measuring Gas Emissions from Agricultural and
Animal Feeding Operations, Crit. Rev. Anal. Chem., 44, 200–219, 10.1080/10408347.2013.843055, 2014.Hu, N., Flesch, T. K., Wilson, J. D., Baron, V. S., and Basarab, J. A.:
Refining an inverse dispersion method to quantify gas sources on rolling
terrain, Agr. Forest Meteorol., 225, 1–7,
10.1016/j.agrformet.2016.05.007, 2016.Kamp, J. N., Chowdhury, A., Adamsen, A. P. S., and Feilberg, A.: Negligible influence of livestock contaminants and sampling system on ammonia measurements with cavity ring-down spectroscopy, Atmos. Meas. Tech., 12, 2837–2850, 10.5194/amt-12-2837-2019, 2019.Kamp, J. N., Häni, C., Nyord, T., Feilberg, A., and Sørensen, L. L.:
Calculation of NH3 Emissions, Evaluation of Backward Lagrangian Stochastic
Dispersion Model and Aerodynamic Gradient Method, Atmosphere, 12, 102, 10.3390/atmos12010102, 2021.Kupper, T., Eugster, R., Sintermann, J., and Häni, C.: Ammonia emissions
from an uncovered dairy slurry storage tank over two years: Interactions
with tank operations and meteorological conditions, Biosystems Engineering,
204, 36–49, 10.1016/j.biosystemseng.2021.01.001, 2021.Lemes, Y. M., Garcia, P., Nyord, T., Feilberg, A., and Kamp, J. N.:
Full-Scale Investigation of Methane and Ammonia Mitigation by Early
Single-Dose Slurry Storage Acidification, ACS Agric. Sci. Technol., 2,
1196–1205, 10.1021/acsagscitech.2c00172, 2022.Lynn, B. H. and Carlson, T. N.: A stomatal resistance model illustrating
plant vs. external control of transpiration, Agr. Forest Meteorol., 52, 5–43, 10.1016/0168-1923(90)90099-R, 1990.Massad, R.-S., Nemitz, E., and Sutton, M. A.: Review and parameterisation of bi-directional ammonia exchange between vegetation and the atmosphere, Atmos. Chem. Phys., 10, 10359–10386, 10.5194/acp-10-10359-2010, 2010.McBain, M. C. and Desjardins, R. L.: The evaluation of a backward Lagrangian
stochastic (bLS) model to estimate greenhouse gas emissions from
agricultural sources using a synthetic tracer source, Agr. Forest Meteorol., 135, 61–72, 10.1016/j.agrformet.2005.10.003, 2005.McGinn, S. M., Flesch, T. K., Crenna, B. P., Beauchemin, K. A., and Coates,
T.: Quantifying Ammonia Emissions from a Cattle Feedlot using a Dispersion
Model, J. Environ. Qual., 36, 1585–1590, 10.2134/jeq2007.0167, 2007.McGinn, S. M., Coates, T., Flesch, T. K., and Crenna, B.: Ammonia emission
from dairy cow manure stored in a lagoon over summer, Can. J. Soil. Sci.,
88, 611–615, 10.4141/CJSS08002, 2008.McGinn, S. M., Turner, D., Tomkins, N., Charmley, E., Bishop-Hurley, G., and
Chen, D.: Methane Emissions from Grazing Cattle Using Point-Source
Dispersion, J. Environ. Qual., 40, 22–27,
10.2134/jeq2010.0239, 2011.McGinn, S. M., Flesch, T. K., Beauchemin, K. A., Shreck, A., and Kindermann,
M.: Micrometeorological Methods for Measuring Methane Emission Reduction at
Beef Cattle Feedlots: Evaluation of 3-Nitrooxypropanol Feed Additive, J.
Environ. Gual., 48, 1454–1461, 10.2134/jeq2018.11.0412, 2019.
NEC Directive 2016/2284/EU: Directive (EU) 2016/2284 of the European Parliament and of the Council of 14 December 2016 on the reduction of national emissions of certain atmospheric pollutants, amending Directive 2003/35/EC and repealing Directive 2001/81/EC, 2016.OECD and FAO: OECD-FAO Agricultural Outlook 2019-2028, OECD,
10.1787/agr_outlook-2019-en, 2019.
Pedersen, J. M., Feilberg, A., Kamp, J. N., Hafner, S., and Nyord, T.:
Ammonia emission measurement with an online wind tunnel system for evaluation of manure application techniques, Atmos. Environ., 230, 117562, 10.1016/j.atmosenv.2020.117562, 2020.
Platt, U. and Stutz, J.: Differential optical absorption spectroscopy:
principles and applications, Springer, Berlin, 597 pp., ISBN: 978-3-540-21193-8, 2008.
R Core Team: R: A language and environment for statistical computing; R
Foundation for Statistical, Computing, Vienna, Austria, 2018.Ro, K. S., Johnson, M. H., Hunt, P. G., and Flesch, T. K.: Measuring Trace
Gas Emission from Multi-Distributed Sources Using Vertical Radial Plume
Mapping (VRPM) and Backward Lagrangian Stochastic (bLS) Techniques,
Atmosphere, 2, 553–566, 10.3390/atmos2030553, 2011.Ro, K. S., Stone, K. C., Johnson, M. H., Hunt, P. G., Flesch, T. K., and
Todd, R. W.: Optimal Sensor Locations for the Backward Lagrangian Stochastic
Technique in Measuring Lagoon Gas Emission, J. Environ. Qual., 43, 1111–1118, 10.2134/jeq2013.05.0163, 2014.Sanz, A., Misselbrook, T., Sanz, M. J., and Vallejo, A.: Use of an inverse
dispersion technique for estimating ammonia emission from surface-applied
slurry, Atmos. Environ., 44, 999–1002, 10.1016/j.atmosenv.2009.08.044, 2010.Shah, S. B., Grabow, G. L., and Westerman, P. W.: Ammonia Adsorption in Five
Types of Flexible Tubing Materials, Appl. Eng. Agric., 22, 919–923, 10.13031/2013.22253, 2006.Sintermann, J., Ammann, C., Kuhn, U., Spirig, C., Hirschberger, R., Gärtner, A., and Neftel, A.: Determination of field scale ammonia emissions for common slurry spreading practice with two independent methods, Atmos. Meas. Tech., 4, 1821–1840, 10.5194/amt-4-1821-2011, 2011.Sintermann, J., Dietrich, K., Häni, C., Bell, M., Jocher, M., and Neftel, A.: A miniDOAS instrument optimised for ammonia field measurements, Atmos. Meas. Tech., 9, 2721–2734, 10.5194/amt-9-2721-2016, 2016.Sommer, S. G., McGinn, S. M., Hao, X., and Larney, F. J.: Techniques for
measuring gas emissions from a composting stockpile of cattle manure,
Atmos. Environ., 38, 4643–4652, 10.1016/j.atmosenv.2004.05.014, 2004.Todd, R. W., Cole, N. A., Rhoades, M. B., Parker, D. B., and Casey, K. D.:
Daily, Monthly, Seasonal, and Annual Ammonia Emissions from Southern High
Plains Cattle Feedyards, J. Environ. Qual., 40, 1090–1095,
10.2134/jeq2010.0307, 2011.Vaittinen, O., Metsälä, M., Persijn, S., Vainio, M., and Halonen, L.: Adsorption of ammonia on treated stainless steel and polymer surfaces, Appl. Phys. B-Laser. O., 115, 185–196, 10.1007/s00340-013-5590-3, 2014.VanderZaag, A. C., Flesch, T. K., Desjardins, R. L., Baldé, H., and
Wright, T.: Measuring methane emissions from two dairy farms: Seasonal and
manure-management effects, Agr. Forest Meteorol., 194,
259–267, 10.1016/j.agrformet.2014.02.003, 2014.van Haarlem, R. P., Desjardins, R. L., Gao, Z., Flesch, T. K., and Li, X.:
Methane and ammonia emissions from a beef feedlot in western Canada for a
twelve-day period in the fall, Can. J. Anim. Sci., 88, 641–649,
10.4141/CJAS08034, 2008.Vechi, N. T., Mellqvist, J., and Scheutz, C.: Quantification of methane
emissions from cattle farms, using the tracer gas dispersion method,
Agr. Ecosyst. Environ., 330, 107885,
10.1016/j.agee.2022.107885, 2022.Voglmeier, K., Jocher, M., Häni, C., and Ammann, C.: Ammonia emission measurements of an intensively grazed pasture, Biogeosciences, 15, 4593–4608, 10.5194/bg-15-4593-2018, 2018.Wesely, M.: Parameterization of surface resistances to gaseous dry
deposition in regional-scale numerical models?, Atmos. Environ., 41,
52–63, 10.1016/j.atmosenv.2007.10.058, 2007.
Wilson, J. D., Flesch, T. K., and Harper, L. A.: Micro-meteorological
methods for estimating surface exchange with a disturbed windflow,
Agr. Forest Meteorol., 107, 207–225,
10.1016/S0168-1923(00)00238-0, 2001.Yang, W., Que, H., Wang, S., Zhu, A., Zhang, Y., He, Y., Xin, X., and Zhang,
X.: Comparison of backward Lagrangian stochastic model with
micrometeorological mass balance method for measuring ammonia emissions from
rice field, Atmos. Environ., 211, 268–273,
10.1016/j.atmosenv.2019.05.028, 2019.YoLemes: AU-BCE-EE/Lemes-2023-ControlReleaseExperiments-Paper: Lemes-2023-ControlReleaseExperiments-Paper, Version v2, Zenodo [code], 10.5281/zenodo.7695569, 2023.