# Accounting for Field-Scale Dry Deposition in Backward Lagrangian Stochastic Dispersion Modelling of NH3 Emissions

C. Hani, Christophe Flechard, A. Neftel, J. Sintermann, T. Kupper

► To cite this version:

C. Hani, Christophe Flechard, A. Neftel, J. Sintermann, T. Kupper. Accounting for Field-Scale Dry Deposition in Backward Lagrangian Stochastic Dispersion Modelling of NH3 Emissions. Atmosphere, 2018, 9 (4), pp.146. 10.3390/atmos9040146. [hal-02056884](https://hal.science/hal-02056884)

HAL Id: hal-02056884

https://hal.science/hal-02056884

Submitted on 4 Mar 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Accounting for Field-Scale Dry Deposition in Backward Lagrangian Stochastic Dispersion Modelling of NH₃ Emissions

Christoph Häni <fnref n="1" />,* ID, Christophe Flechard <fnref n="2" />, Albrecht Neftel <fnref n="3" />,†, Jörg Sintermann <fnref n="3" />,‡ and Thomas Kupper <fnref n="1" />

Received: 27 February 2018; Accepted: 4 April 2018; Published: 14 April 2018

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** A controlled ammonia (NH₃) release experiment was performed at a grassland site. The aim was to quantify the effect of dry deposition between the source and the receptors (NH₃ measurement locations) on emission rate estimates by means of inverse dispersion modelling. NH₃ was released for three hours at a constant rate of Q = 6.29 mg s⁻<fnref n="1" /> from a grid of 36 orifices spread over an area of 250 m<fnref n="2" />. The increase in line-integrated NH₃ concentration was measured with open-path optical miniDOAS devices at different locations downwind of the artificial source. Using a backward Lagrangian stochastic (bLS) dispersion model (bLSmodelR), the fraction of the modelled release rate to the emitted NH₃ (QbLS/Q) was calculated from the measurements of the individual instruments. QbLS/Q was found to be systematically lower than 1, on average between 0.69 and 0.91, depending on the location of the receptor. We hypothesized that NH₃ dry deposition to grass and soil surfaces was the main factor responsible for the observed depletion of NH₃ between source and receptor. A dry deposition algorithm based on a deposition velocity approach was included in the bLS modelling. Model deposition velocities were evaluated from a 'big-leaf' canopy resistance analogy. Canopy resistances (generally termed R_c) that provided QbLS/Q = 1 ranged from 75 to 290 s m⁻<fnref n="1" />, showing that surface removal of NH₃ by dry deposition can plausibly explain the original underestimation of QbLS/Q. The inclusion of a dry deposition process in dispersion modelling is crucial for emission estimates, which are based on concentration measurements of depositing tracers downwind of homogeneous area sources or heterogeneously-distributed hot spots, such as, e.g., urine patches on pastures in the case of NH₃.

**Keywords:** backward Lagrangian stochastic model; atmospheric surface-layer; micrometeorological techniques; atmospheric ammonia; dry deposition; open-path measurements; differential optical absorption spectroscopy

## 1. Introduction

Estimation of trace gas emission from confined source areas on a local scale (i.e., receptor (measurement) distance to sources less than 500 m) using the combination of inverse dispersion modelling with either concentration or flux measurements is a widespread method, especially in the agricultural sector (e.g., [1–6]).

<footnote n="1">¹ School of Agricultural, Forest and Food Sciences HAFL, Bern University of Applied Sciences, 3052 Zollikofen, Switzerland; thomas.kupper@bfh.ch</footnote>

<footnote n="2">² INRA, Agrocampus Ouest, UMR 1069 SAS, 35042 Rennes, France; christophe.flechard@inra.fr</footnote>

<footnote n="3">³ Climate and Agriculture Group, Agroscope, 8046 Zürich, Switzerland; neftel_a@bluewin.ch (A.N.); sintermann@gmx.ch (J.S.)</footnote>

<footnote>* Correspondence: christoph.haeni@bfh.ch; Tel.: +41-31-910-2269</footnote>

<footnote>† Now at: Neftel Research Expertise, 3033 Wohlen b. Bern, Switzerland.</footnote>

<footnote>‡ Now at: Office of Waste, Water, Energy and Air, Canton of Zurich, 8090 Zurich, Switzerland.</footnote>
In particular, the combination of concentration measurements with backward Lagrangian stochastic (bLS) modelling is a convenient way of emission estimation that has spurred its utilization in the past decade [7]. Flesch et al. [8] present details of a surface-layer bLS model that has been embedded in the framework of the frequently used software 'WindTrax' (http://www.thunderbeachscientific.com/). A variety of trace gases, such as, e.g., methane (CH₄) or ammonia (NH₃), have been investigated using this model (e.g., [9–13]).

While gases like CH₄ satisfy the model assumption of an inert gas, NH₃ shows a high affinity to sorb on any kind of surface (e.g., [14,15]) and measurements of NH₃ downwind of a source area will therefore be affected by the removal of NH₃ through deposition between source and receptor. Thus, the concentration measurement downwind will reflect the net emission that corresponds to the NH₃ exchange in the entire sector between source and receptor. If the aim is the quantification of the emission from the emitting source areas only (i.e., the gross emission flux), the deposition component must be quantified.

For short ranges, i.e., within the first few hundred meters downwind of a source area, the removal of NH₃ is largely dominated by dry deposition [16,17]. Emission estimates without the inclusion of the dry deposition mechanism in the model run have the potential to be biased towards smaller values. The extent of this effect is poorly known and the bias due to dry deposition is usually neglected in results from field measurements.

Artificial sources can be used to validate short-range atmospheric dispersion models. Due to practical reasons, many artificial sources consist of a single or an ensemble of point releases, e.g., critical orifices. If the artificial source is located in or close to the canopy, a significant part of surface-interacting trace gases, such as NH₃ is deposited within and near the source grid. This results in recovered fractions of the emitted tracer that are systematically below 1.

In this paper, we report the extension of a bLS model to include an explicit dry deposition algorithm and we use the results from a controlled release experiment to estimate the magnitude of dry deposition between source and receptor locations. The model environment has been embedded in an R package [18] denominated bLSmodelR. The newest version of bLSmodelR is available at https://www.agrammon.ch/documents-to-download/blsmodelr/.

## 2. Materials and Methods

### 2.1. Backward Lagrangian Stochastic Dispersion Model

The employed bLS model in bLSmodelR is a first-order Lagrangian stochastic dispersion model that is run in backward mode (i.e., backward in time) assuming horizontally homogeneous and vertically inhomogeneous, Gaussian turbulence. The vertical profiles of the wind speed and the turbulence statistics are based on Monin-Obukhov Similarity Theory (MOST).

The relation between the emission rate (Q) of a homogeneous surface area and the corresponding concentration increase (C) at a receptor location, is calculated as:

$$C/Q = \frac{1}{N} \sum_{\text{TD inside}} \left| \frac{2}{w_{\text{TD}}} \right|, \quad (1)$$

where $w_{\text{TD}}$ is the vertical velocity at touchdown and N is the total number of trajectories in the model run that have been released and calculated back-in-time from the receptor location. The sum consists of all touchdowns inside of the emitting area. The concentration increase C represents the difference between the concentration with and without the source present. The latter is often referred to as background concentration.

All the dispersion modelling relevant equations in bLSmodelR were obtained from Flesch et al. [8]. Model aspects that are different from or not specified in Flesch et al. [8] are the following:

* For the numerical discretization, the mixed implicit-explicit Euler scheme described in Flesch et al. [19] is used.
* The trajectories are reflected perfectly at an effective ground level, taken as the level of the roughness length $z_0$ above the displacement height $d$, such that the covariance of the wind components is retained.

* The initialization of the wind components at release from the receptor location is done using an orthogonal projection procedure (http://stats.stackexchange.com/questions/15011/generate-a-random-variable-with-a-defined-correlation-to-an-existing-variable) that has been adapted for three correlated components. This guarantees a Gaussian distribution of initial wind components with a covariance matrix given by:

$$ \operatorname{cov}(\boldsymbol{u}_{\text{ini}}) = \begin{pmatrix} \sigma_u^2 & 0 & -u_*^2 \\ 0 & \sigma_v^2 & 0 \\ -u_*^2 & 0 & \sigma_w^2 \end{pmatrix}, \qquad (2) $$

where $u_*$ is the friction velocity (Section 2.4.3.) and $\sigma_u^2$, $\sigma_v^2$, and $\sigma_w^2$ are the variances of the along-wind, the crosswind, and the vertical wind components.

* The theoretical wind profile in the model is extended by an additional term $\psi_M(\frac{z_0}{L})$ and defined according to MOST:

$$ \bar{u} = \frac{u_*}{k_v} \left\{ \ln\left(\frac{z-d}{z_0}\right) - \psi_M\left(\frac{z-d}{L}\right) + \psi_M\left(\frac{z_0}{L}\right) \right\} \qquad (3) $$

with the stability correction function for momentum given as:

$$ \psi_M(x) = \begin{cases} -4.8x & , \text{ if } L > 0 \\ 2\ln\left(\frac{1+\alpha}{2}\right) + \ln\left(\frac{1+\alpha^2}{2}\right) - 2\operatorname{atan}(\alpha) + \frac{\pi}{2} & , \text{ if } L < 0 \end{cases} \qquad (4) $$

and:

$$ \alpha = (1 - 16x)^{1/4}. \qquad (5) $$

where $x$ is either $(z - d)/L$ or $z_0/L$, $z$ represents the geometric height (i.e., the height above ground), $L$ is the Obukhov-Length, $d$ is the displacement height, $z_0$ is the roughness length, $\bar{u}$ represents the temporal average of the along-wind velocity, and $k_v$ is the von Kármán constant (defined as 0.4).

## 2.2. Dry Deposition Modelling

### 2.2.1. Basic Principle

The deposition modelling is coupled to the dispersion modelling as a concentration modifying process that acts on each touchdown occurring outside of the predefined source areas (Figure 1).

The flux ($F$) between the absorbing surface (at the notional height $z'_0$) and the model effective ground level $z = d + z_0$, the level where touchdowns occur in the bLS model, can be approximated by the product of an exchange velocity $v_{ex}$, active through the quasi-laminar viscous sub-layer, and the difference between the concentration at the absorbing surface $C_{\{d+z_0\}}$ and the concentration at the model surface $C_{\{z'_0\}}$:

$$ F = v_{\text{ex}} (C_{\{z'_0\}} - C_{\{d+z_0\}}). \qquad (6) $$

The super-positioning of concentrations in the bLS model (i.e., the concentration of a single trajectory can be written as the sum of concentration contributions from the background and the different sources: $C_{\text{Traj}} = C_{\text{bgd}} + \sum_i C_{\text{Source}_i}$) allows to reformulate Equation (6) to:

$$ F = v_{\text{ex}} \left( C_{\{z'_0\}}|_{\text{bgd}} - C_{\{d+z_0\}}|_{\text{bgd}} + C_{\{z'_0\}}|_{\text{Sources}} - C_{\{d+z_0\}}|_{\text{Sources}} \right) = F|_{\text{bgd}} + F|_{\text{Source}}. \qquad (7) $$
With respect to the source-receptor relationship, only the flux related to the source emission is needed and the terms that are related to the background concentration are discarded without any loss of information. Further, it is assumed, that the increase in $C_{\{z'_{0}\}|_{Sources}}$ due to already deposited molecules from the source areas can be neglected. Thus, it holds that $C_{\{z'_{0}\}|_{Sources}} = 0$ and Equation (7) can be written as:

$$F = -v_{ex}C_{TD}, \qquad (8)$$

where $C_{TD}$ is used as synonym for $C_{\{d+z_{0}\}|_{Sources}}$, the modelled concentration at touchdown of the trajectory. Since $C_{TD}$ is strictly positive, the right-hand side of Equation (8) is always negative, i.e., a downward directional deposition flux $F_{d}$, and the exchange velocity represents in fact a surface deposition velocity $v_{d}^{*}$:

$$F_{d} = -v_{d}^{*}C_{TD}. \qquad (9)$$

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Illustration of the deposition modelling in the bLSmodelR. The trajectory touchdowns (indicated by circles) at -2.5 s and -2.7 s are inside the source area and not affected by deposition modelling; those at -1.9 s and 0.0 s are outside the source area where deposition occurs in the model.</figcaption></figure>

### 2.2.2. Integration of Dry Deposition in the bLS Model

The integration of Equation (9) into the model, whenever a touchdown occurs outside the source areas, is done by maintaining the conservation of mass for each trajectory. The change in concentration with time in a volume $V$ equals the deposition flux through the area $A$ which is orthogonal to the average flux:

$$V \frac{dC_{TD}}{dt} = AF_{d}. \qquad (10)$$

Dividing Equation (10) by $V$ and treating volume $V$ as notional volume where the deposition takes place with an infinitesimal vertical extension $dz$ (thus, $V = Adz$) and further replacing $F_{d}$ by Equation (9) leads to:

$$\frac{dC_{TD}}{dt} = \frac{F_{d}}{V}A = -\frac{v_{d}^{*}}{dz}C_{TD}, \qquad (11)$$
Thus, at each touchdown, the concentration decrease can be described by integration of Equation (11) over time as:

$$C_{>TD} = C_{<TD} \exp\left(-\frac{v_d^*}{dz} t_{TD}\right), \qquad (12)$$

where $C_{>TD}$ is the concentration after and $C_{<TD}$ the concentration before a touchdown occurs. The 'residence time' $t_{TD}$ (i.e., the amount of time that the trajectory is below $z_0$, in the deposition relevant volume) is given as:

$$t_{TD} = \left| \frac{2dz}{w_{TD}} \right| \qquad (13)$$

and thus:

$$C_{>TD} = C_{<TD} \exp\left(-\left|\frac{2v_d^*}{w_{TD}}\right|\right) \qquad (14)$$

or equivalent, when dividing both sides by the source strength $Q$:

$$\{C/Q\}_{>TD} = \{C/Q\}_{<TD} \exp\left(-\left|\frac{2v_d^*}{w_{TD}}\right|\right). \qquad (15)$$

The deposition pathway inside emitting areas is switched off, i.e., no deposition is modelled. The concentration after the touchdown occurred is only affected by the usual increase due to the source area emission, thus:

$$\{C/Q\}_{>TD} = \{C/Q\}_{<TD} + \left|\frac{2}{w_{TD}}\right|, \qquad (16)$$

if the touchdown occurs inside the source areas.

The modelled $C/Q$ ratio for a single trajectory $i$ is given as the cumulative sum over all touchdowns $j \in \{1, \dots, M\}$, starting at the first (i.e., first in time) touchdown inside the source:

$$\{C/Q\}_{\text{dep}i} = \sum_{j=1}^{M} \begin{cases} \{C/Q\}_{\text{dep},j-1} \left( \exp\left(-\left|\frac{2v_d^*}{w_{TD}}\right|\right) - 1 \right), & \text{if TD outside source} \\ \left|\frac{2}{w_{TD}}\right|, & \text{if TD inside source} \end{cases} \qquad (17)$$

with:

$$\{C/Q\}_{\text{dep},i,j} = \sum_{k<j} \{C/Q\}_{\text{dep},i,k}. \qquad (18)$$

The final modelled $C/Q$ ratio is given as the average value of all trajectories:

$$\{C/Q\}_{\text{dep}} = \frac{1}{N} \sum_{i=1}^{N} \{C/Q\}_{\text{dep}i}. \qquad (19)$$

The deposition algorithm described above is very similar to the deposition implementation described in Wilson et al. [20]. The two approaches differ in their derivation, but provide comparable results.

## 2.3. Structure and Implementation of the bLSmodelR

The bLS model by Flesch et al. [8] together with the deposition algorithm was implemented into an R [18] package called bLSmodelR (Supplement S1, current version available at https://www.agrammon.ch/documents-to-download/blsmodelr/). The deposition algorithm was included in the bLSmodelR package as a post-processing function that allows the modelling of the dry deposition additionally to the dispersion modelling. It provides reduced $C/Q$ ratios due to the user defined dry deposition velocity, where the bLS model needs to be run first in normal mode, i.e., without the dry deposition mechanism included.

The bLSmodelR package makes use of several different available R packages. A complete list of all packages that bLSmodelR depends on can be found on the above-mentioned homepage. Amongst
the most important ones are used mainly due to CPU time, the packages *snowfall* [21] and *snow* [22] for simple parallelization of calculation processes, *Rcpp* [23] for the simple implementation of the core model's C++ code (i.e., the calculation of the individual trajectories) and *data.table* [24] for fast calculation in the model framework.

## 2.4. NH₃ Release Experiment

### 2.4.1. Experimental Site

The experiment took place on 22 September 2014 at a grassland site (46°59′06″ N/07°04′35″ E) in the Seeland region, Switzerland. The release started at 10h40 local time and lasted until 13h42. The sward was treated with a broad-spectrum herbicide (glyphosate) one week before the release experiment. The height of the canopy was estimated as 9.6 cm on average with a rising plate pasture meter (RPM; Filip's Manual Folding Plate Meter, NZ Agriworks Ltd. t/a Jenquip, Feilding, New Zealand).

During the experiment, the wind direction varied between NE to E (45° to 103°), where the only obstacles upwind of the setup was a hedge of lower height at approx. 140 m and on the eastern end, a group of tall trees at approx. 200 m distance. The terrain was horizontally flat and homogeneous within 140 m upwind of the artificial source. There were two farms located within a range of 1 km, one farm at 700 m distance NE (42°) with pasture next to the farm buildings and another farm at a distance of 500 m SSE (163°).

### 2.4.2. Experimental Site

The artificial source included a gas distributor unit at the center where 12 branches consisting of polyamide tubes (6 mm outer diameter × 4 mm inner diameter, Legris, Parker Hannifin, Mayfield Heights, OH, USA) were attached radially forming a hexagonal geometry (Figure 2). Each branch contained three critical flow orifices (100 µm diameter, stainless steel, Lenox Laser Inc., Glen Arm, MD, USA) in series with 3 m distance between neighboring orifices. The source grid was placed at ground level and covered a notional circular area of approximately 250 m². The gas bottle with a nominal NH₃ to N₂ mixing ratio of 4.97% (±2% rel., Carbagas, Bern, Switzerland) was placed 20 m away from the source orifices. The constant release rate of 10 Lₙ min⁻¹ was regulated and measured by a mass flow controller (red-y smart controller, Voegtlin Instruments GmbH, Aesch, Switzerland). The NH₃ total release rate of the source was given as 6.28 mg NH₃ s⁻¹ (174.6 µg NH₃ s⁻¹ per orifice). The response time of the source, i.e., the time between the switching-on of the flow and the reaching of the total NH₃ release rate at the exit (orifices), is in the order of minutes.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Schematic of the artificial source setup. Black/grey circles with crosses: orifices; black lines: polyamide tubes.</figcaption></figure>
### 2.4.3. Dispersion Model Input

The high-frequency fluctuations of 3D wind components along with the sonic temperature were measured at 10 Hz using an ultrasonic anemometer (WindMaster™Pro, Gill Instruments Limited, Lymington, UK). Since the raw data of the WindMaster™Pro series includes a bug produced by the firmware they were corrected according to the indications communicated by Gill instruments (http://gillinstruments.com/data/manuals/KN1509_WindMaster_WBug_info.pdf). The ultrasonic anemometer (denominated 'Sonic') was placed next to the NH₃ measurements at a height of 1.25 m above ground level (a.g.l.) and oriented towards magnetic north using a geological compass. The orientation of the Sonic was corrected by the magnetic declination, i.e., the deviation of the magnetic north to the geographic north. A bubble level was used to align the Sonic in the vertical axis.

The raw data from the Sonic was averaged on a 10-min basis. The wind vector was rotated twice to provide a wind vector with components u, v, and w, representing the along-wind, crosswind and vertical velocity, respectively. The first rotation in the horizontal adjusted the x-axis to the 10-min mean wind direction and the second rotation was done afterwards around the new y-axis until the mean vertical wind equaled zero. All sonic data were linearly detrended and the required model input parameters were calculated from these processed data.

The following set of input parameters must be specified in bLSmodelR: the three MOST parameters {u*, L, z₀}, the wind statistics {σu, σv, σw(zσw)} with corresponding height zσw where the standard deviation of the vertical wind component was estimated, the displacement height d and the wind direction β.

Specifically, the friction velocity was calculated as:

$$u_* = \sqrt{-\overline{u'w'}} \qquad (20)$$

where the overbar represents the temporal average of the product of and the prime denotes the deviation from the average of the rotated wind components u and w, respectively. The displacement height was taken to be 2/3 of the canopy height, resulting in d = 6.4 cm. The value of the roughness length (z₀) was estimated in an iterative procedure solving Equation (3) for z₀ by providing $\bar{u}$, $u_*$ and L from sonic anemometer measurements.

The bLS model was run by supplying direct measurements of the wind statistics σu, σv, σw implying a varying Kolmogorov coefficient (C₀) and related to that, a varying eddy diffusivity equivalent of the model.

### 2.4.4. Concentration Measurement

The NH₃ concentrations in the ambient air were determined with miniDOAS instruments described in detail by Sintermann et al. [25]. They are open-path (thus, sample inlet-free) optical devices that measure the concentration of NH₃, sulfur dioxide (SO₂) and nitric oxide (NO) as line-integrated averages over a measurement path between a light source and a detector (spectrometer) by UV absorption in the wavelength range of 200 to 230 nm. The light path is returned by a retroreflector typically at 10 to 50 m from the light source to keep the light source and the detector at the same location in an actively ventilated, temperature stabilized box (Figure 3).

The four miniDOAS were placed 15 m downwind of the source center at three heights (locations denominated near/bottom at 0.5 m, near/middle at 1.25 m and near/top at 3.0 m) and 70 m downwind of the source center at a height of 1.2 m (location denominated far) (Figure 4). The measurement paths were approx. 70 m (i.e., distance between box and reflector ≈ 35 m) for all miniDOAS instruments.

The positions of all instruments and the individual orifices were recorded by a Global Positioning System (approx. precision 10 cm, GPS Trimble R8 GNSS, Sunnyvale, CA, USA).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Three miniDOAS instruments located at the release site, 15 m downwind of the source center (position *near*). The distances between the boxes containing the light source and the detector (picture foreground) and the reflectors (picture background) were approx. 35 m. The measurement heights were 0.5 m (*near/bottom*), 1.25 m (*near/middle*), and 3.0 m (*near/top*) above ground.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Position of the source (orifices) and the instruments during the artificial source release experiment as measured with the GPS device. Three miniDOAS (*near/top*, *near/middle*, and *near/bottom*) were located close to the artificial source at three different heights, one miniDOAS was located further downwind (*far*). The average wind direction during the release experiment is indicated by the blue arrow.</figcaption></figure>
### 2.4.5. Background Concentration

The background concentration for each individual miniDOAS instrument was linearly interpolated from the concentration measurement before and after the release. The concentration measurements during daytime several days before and after the release experiment with wind from the northeasterly sector (the prevailing wind direction during the release experiment) showed little variation over time and a linear interpolation of the background concentration seems appropriate for the short time of the release. The measured background concentration immediately before and after the release was around 4 µg m⁻³.

### 2.4.6. Recovered Fraction of the Tracer Gas

The ratio of the modelled to the emitted release rate (i.e., the recovered fraction of NH₃ by the model) was calculated as:

$$Q_{bLS}/Q = \frac{C_m - C_b}{Q_{src}\left\{\frac{C}{Q}\right\}_{sim}} \quad (21)$$

where C_b is the contribution of the background concentration to the measured total concentration C_m, both in µg NH₃ m⁻³, Qsrc is the emission rate of the source in µg NH₃ m⁻² s⁻¹ and {C/Q}sim is the modelled C/Q ratio in s m⁻¹ calculated from Equation (1) if dry deposition was not included or Equation (19) if deposition was included in the model run. For the model calculation of {C/Q}sim the source area was defined as 36 circular areas of 5 cm radius, all emitting with equal strength of Qsrc = 22.24 mg NH₃ m⁻² s⁻¹.

### 2.4.7. Surface Deposition Velocity

For each estimate of the recovery rate QbLS/Q (i.e., each interval at each location), the near-surface deposition velocity (v\*d) (Section 2.2) was evaluated by running the model deposition post-processing with different values of v\*d and its optimum was estimated from these results. The values of v\*d that were provided as input to the bLS dry deposition modelling were approximated by a resistances approach [26], where v\*d is modelled as the inverse of the sum of a series of different resistances to deposition:

$$v_{d}^{*} = \frac{1}{R_{b} + R_{c}}. \quad (22)$$

R_b and R_c represent transfer resistances below the aerodynamic surface (d + z₀), with R_b the pseudo-laminar boundary layer resistance and R_c a bulk (big-leaf) canopy resistance, respectively. The aerodynamic resistance (R_a) usually included in resistance schemes is not included in Equation (22), since v\*d is defined between d + z₀ and the absorbing surface (z'₀), and, therefore, R_a is already implicitly accounted for in bLS dispersion modelling.

The pseudo-laminar resistance R_b was calculated following Garland [27] as:

$$R_b = \frac{1.45 \left( \frac{z_0 u_*}{v} \right)^{0.24} \left( \frac{v}{D} \right)^{0.8}}{u_*}, \quad (23)$$

where z₀ is the roughness length, u\* the friction velocity, v the kinematic viscosity of air and D the molecular diffusivity of NH₃ in air. The temperature and pressure dependence of the diffusion coefficient was approximated according to the Chapman-Enskog theory of gas diffusion [28]:

$$D = D_0 \frac{p_0}{p} \left( \frac{T}{T_0} \right)^{3/2}, \quad (24)$$
where *p* is the actual air pressure and *T* is the actual air temperature. The standard diffusion coefficient of ammonia in air (D₀ = 0.20487 cm² s⁻¹ at p₀ = 1 atm and T₀ = 273 K) is calculated as the average of the values summarized in Tang et al. [29].

To estimate the effective magnitude of $v_d^*$, based on the assumption that dry deposition was the main reason for the observed incomplete recovery rates, $R_c$. was varied from 0 to 500 s m⁻¹ in steps of 50 s m⁻¹. In a first step, the optimum of $R_c$ (henceforth referred to as $R_c^1$) was estimated with local polynomial regression (R function *loess*) from the post-processed results such that $Q_{bLS}/Q = 1$. This yields a 'best estimate' of $R_c$ anywhere in the range between 0 and $\infty$ s m⁻¹. In a second step, the corresponding optimum of $v_d^*$ was calculated from Equation (22).

The reason for the indirect analysis of $v_d^*$ via $R_c$ is given in $R_c$ being a surface property that is expected to change only marginally during the release experiment in contrast to the turbulence affected $v_d^*$. Furthermore, $R_c$ can, in theory, be approximated based on certain environmental parameters like air temperature, relative humidity and SO₂ concentration as discussed in Section 4.2.

# 3. Results

## 3.1. *Environmental Data*

### 3.1.1. Meteorological Conditions and Dispersion Parameters

On the day of the experiment, the weather was sunny with a steady easterly wind. During the release, the air temperature increased from 15.5 to 18.9 °C with an average of 17.6 °C, whereas the relative humidity dropped from 61% to 43% (average: 51%) (Figure 5e). The 10-min averaged wind speed at 1.25 meters above ground level from the sonic anemometer measurements ranged from 1.7 to 2.7 m s⁻¹ (average: 2.2 m s⁻¹) with wind directions ranging between 45° and 103° with an average of 74° (Figure 5c).

The friction velocity $u_*$ ranged between 0.16 and 0.26 m s⁻¹. The Obukhov length *L* ranged from -1.3 to -7.1 m (Figure 5d). The roughness length *z*₀ estimated from Equation (3) averaged to 1.0 cm (range: 0.2 to 2.4 cm). The scaled standard deviations for the three components of wind ranged from 3.1 to 4.8 (average: 4.0) for $\sigma_u/u_*$, 3.0 to 6.9 (average: 4.5) for $\sigma_v/u_*$ and 1.15 to 1.90 (average: 1.47) for $\sigma_w/u_*$, resulting in a range of $b_w$ ($\sigma_w/u_*$ at ground level) from 0.99 to 1.35 (average: 1.16) or a corresponding Kolmogorov coefficient $C_0$ of the bLS model between 3.2 and 5.1 (average: 3.9).
<figure><img src="image_8.png" /></figure>

**Figure 5.** Overview on the measured variables during the release experiment. Thin/thick lines show the 1-/10-min averages of the measurements. Black vertical dotted lines indicate the start at 10h40 and the end at 13h42 (local time) of the NH₃ release. Panels (top to bottom): (a) Measured NH₃ concentration increase above background for the miniDOAS instruments at position *near*: light grey: *near/bottom*; grey: *near/middle*; black: *near/top*. (b) Measured NH₃ concentration increase above for the miniDOAS instruments at position *far*. (c) Wind direction β in ° from N (black lines) and wind speed *U* in m s⁻¹ (grey lines) measured with the sonic anemometer at 1.25 m a.g.l. (d) Friction velocity *u*ₗ in m s⁻¹ (black line) and Obukhov Length *L* in m (grey line) as calculated on a 10-min basis. (e) Air temperature *T* in °C (black line) and relative humidity *RH* in % (grey line) measured at 1.25 m a.g.l. (f) Ensemble of the SO₂ concentration measurements from the individual instruments. The SO₂ plume from the Bárðarbunga volcano passed the experimental area during the NH₃ release.
### 3.1.2. SO₂ Plume

Simultaneous to the NH₃ concentration, the concentrations of NO and SO₂ were recorded with the miniDOAS instruments. Co-deposition of SO₂ is likely to influence the magnitude of dry deposition of NH₃ [30]. An SO₂ plume originating from the eruption of the Icelandic volcano Bárðarbunga passed the location from 08h00 (local time) to 17h00 [25], raising the SO₂ concentration from 20 to 70 µg m⁻³ during the NH₃ release with highest concentration just at the end of the release experiment (Figure 5f). There was no detectable difference in the SO₂ concentration between the three miniDOAS systems close to the NH₃ source, but the instrument at position *far* was higher by roughly 10%.

### 3.2. NH₃ Concentration

The average increase in the measured NH₃ concentration during the NH₃ release period amounted to 39.7 µg m⁻³ (*near/bottom*), 22.6 µg m⁻³ (*near/middle*), 6.3 µg m⁻³ (*near/top*), and 2.7 µg m⁻³ (*far*) (Figure 5a,b). The standard deviation σNH₃ of the measured concentration increased from 0.5 µg m⁻³, 0.3 µg m⁻³, 0.2 µg m⁻³, and 0.8 µg m⁻³ before and after the release to 11.5 µg m⁻³, 6.0 µg m⁻³, 3.2 µg m⁻³, and 2.3 µg m⁻³ during the release for the *near/bottom*, *near/middle*, *near/top*, and *far* measurement location, respectively. The standard deviation of the NH₃ concentration at position *far* during the release phase was of the same order as the average increase above background concentration. The instruments' random errors (i.e., the standard error (SE) of the coefficient estimates from the instrument's curve fitting, reflecting the quality of the measurement) averaged to 0.3 µg m⁻³, 0.5 µg m⁻³, 0.3 µg m⁻³, and 0.6 µg m⁻³ during the release and changed only marginally from background concentration measurements. All instruments showed a similar temporal behavior in the measured concentration, with a residual tailing of the concentration above background after the gas flow through the source stopped, i.e., the NH₃ concentration did not instantaneously drop to background levels, unlike the case for CH₄ (data not shown).

### 3.3. Recovered Fractions of NH₃ without Deposition Modelling

The calculated recovered fractions QbLS/Q from the standard model runs are shown in Figure 6. The average QbLS/Q was significantly below 1 for each instrument, estimated as 0.88 (*near/bottom*), 0.91 (*near/middle*), 0.83 (*near/top*) and 0.69 (*far*). The corresponding standard deviations σQ/Q were smaller in the case of larger fractions recovered, i.e., 0.11 (*near/bottom*, number of intervals n = 18), 0.12 (*near/middle*, n = 18), versus 0.19 (*near/top*, n = 18) and 0.19 (*far*, n = 17). One interval (from 12h10 to 12h20) from the measurement at position *far* showed a very high recovery rate (QbLS/Q = 4.4) exceeding the average QbLS/Q by approx. 20 times the standard deviation, indicative of large measurement or modelling errors and, thus, was excluded from the further analysis. The SE of the individual QbLS/Q estimates (i.e., on an interval basis) were calculated from error propagation. They were also smaller in the case of the larger recovered fractions, i.e., 0.10 (*near/bottom*), 0.11 (*near/middle*), versus 0.13 (*near/top*) and 0.32 (*far*). The error propagation assumed independent error contributions from (a) the NH₃ release rate with an SE of the NH₃ release rate assumed to be 5% of the absolute release rate, (b) the concentration measurement (individual SE given by the miniDOAS instruments), (c) the background concentration estimation (SE assumed to be 0.5 µg m⁻³), and (d) the model calculation of the C/Q ratio (SE given by the stochastic error obtained from the model calculation and a systematic model error, which was assumed to be 10% of the model result).
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Calculated recovery rates $Q_{bLS}/Q$ from model runs without accounting for deposition. The order of the instrument position is given (within an interval) from left (darkest grey) to right (lightest grey) as *near/bottom*, *near/middle*, *near/top*, and *far*. Vertical lines indicate the 95% confidence region of the individual estimates. One interval from the measurement at position *far* (12h10 to 12h20, $Q_{bLS}/Q = 4.4$) was excluded from the analysis.</figcaption></figure>

### 3.4. Canopy Resistance and Surface Deposition Velocity

When considering the ensemble of all intervals, the $R_c$ values that were required to raise the uncorrected $Q_{bLS}/Q$ values closest to 1 ($R_c^1$) and the corresponding average surface deposition velocities ($\bar{v}_d^*$) were 280 s m$^{-1}$, ($\bar{v}_d^* = 0.34$ cm s$^{-1}$) for position *near/bottom*, 290 s m$^{-1}$ ($\bar{v}_d^* = 0.33$ cm s$^{-1}$) for position *near/middle*, 120 s m$^{-1}$ ($\bar{v}_d^* = 0.74$ cm s$^{-1}$) for position *near/top* and 75 s m$^{-1}$ ($\bar{v}_d^* = 1.11$ cm s$^{-1}$) for position *far*. Figure 7 shows the numerical impact of the investigated range of $R_c$ (and $\bar{v}_d^*$) values on the calculated $Q_{bLS}/Q$. The maximum deposition possible ($R_c = 0$ s m$^{-1}$) resulted in $Q_{bLS}/Q$ averages between 1.6 and 2.0.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Sensitivity of the $Q_{bLS}/Q$ ratio to the canopy resistance $R_c$ and the average of the corresponding surface deposition velocities $\bar{v}_d^*$ (secondary x-axis on top). Points and vertical lines indicate the average $Q_{bLS}/Q$ over all intervals with the associated 95% confidence intervals. The line shows the local polynomial regression estimate.</figcaption></figure>
Regarding the results from individual intervals, the medians of $v_d^*$ were 0.28 cm s⁻¹ (near/bottom), 0.42 cm s⁻¹ (near/middle), 0.67 cm s⁻¹ (near/top) and 1.06 cm s⁻¹ (far). Values of $R_c^1$ with the corresponding surface deposition velocities for each interval are given in Table A1 of the Appendix A.

# 4. Discussion

## 4.1. Recovered Fraction of NH₃ without Deposition Modelling

The recovered fractions $Q_{bLS}/Q$ for the standard model runs showed similar patterns with time of the release experiment at the four measurement locations (Figure 6) averaging over all intervals between 0.69 and 0.91 times the expected release rate $Q_{src}$. The variation in between different intervals was small ($\sigma_{Q/Q}$ between 0.11 and 0.19) and compared well to previously-published results (e.g., Table A1 in Harper et al., [12]).

The SE in the $Q_{bLS}/Q$ estimate at position *far* is approx. three times higher than at positions *near/bottom*, *near/middle*, and *near/top*. The following factors explain the larger uncertainty at position *far* compared to the other instruments:

1. the higher sensitivity of $Q_{bLS}/Q$ to errors in the calculated average wind direction at 70 m downwind of the source,

2. measurement at the plume edge (Figure 8) and with that:

a. the higher sensitivity of $Q_{bLS}/Q$ due to the smaller increase in the measured concentration above the background concentration (resulting in a higher sensitivity to the interpolation of the background concentration and to the concentration measurement itself); and

b. the higher sensitivity of $Q_{bLS}/Q$ due to the smaller value of $\left\{ \frac{C}{Q} \right\}_{sim}$.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** (Left) Modelled concentration footprint for instrument *far* at time 12h10 to 12h20. The results from this interval were excluded from the analysis for location *far*. The grey areas indicate the spatial distribution of the relative contribution of surface fluxes to the concentration at position *far* (i.e., the concentration footprint). Wind is blowing from east ($\beta = 97^\circ$), i.e., from right to left. (Right) Modelled concentration distribution in the horizontal plane at height 1.2 m a.g.l. between 12h10 and 12h20, assuming an emission rate identical to the source release of 6.29 mg NH₃ s⁻¹.</figcaption></figure>

The latter two also introduce a higher uncertainty in the results of *near/top* compared to *near/middle*, and *near/bottom*, since location *near/top* was most of the time at the edge of the emission plume.
In previous experiments with CH₄ as a passive tracer using the same artificial source grid as in this study, we observed QbLS/Q values with averages of 0.96 and 1.04 under comparable conditions (see Supplement S2). While there are different possible causes for the low values of QbLS/Q in the present study (e.g., the parametrization of the model dispersion that is reflected i.a. by the value of C₀), we focus on the hypothesis that the systematically lower QbLS/Q values in this study are related to the dry deposition removal of NH₃ between source and receptor.

## 4.2. Deposition Modelling

The canopy resistance R_c can be modelled as the inverse of the sum of the inverse of resistances to the stomatal exchange (usually termed R_s) and the uptake to non-stomatal surfaces (e.g., ground surface, plant cuticle, etc.), usually termed R_w [31]. Since the sward of the entire experimental site was treated with glyphosate one week prior to the measurements, an absence of stomatal activity can be assumed (R_s = ∞ s m⁻¹). Therefore, the experimentally-derived R_c values in the context of the present study mostly describe the magnitude of R_w. Average R_c¹ values that explained the observed QbLS/Q are compatible with the R_w parameterization by Massad et al. [32] (several hundred s m⁻¹ at low relative humidity), and the resulting surface deposition velocities v_d* of 0.33 to 1.11 cm s⁻¹ do not disagree with the deposition velocity ranges published in Schrader and Brümmer [15]. It needs to be mentioned that Schrader and Brümmer [15] reviewed deposition velocity values at a certain, unknown reference height above the canopy and their reported values of 0.1 to 1.8 cm s⁻¹ (average: 0.9 cm s⁻¹) for semi-natural land use would (a) increase by a certain amount if corrected down to a height of d + z₀, but also (b) decrease by a certain amount if corrected for the absence of stomatal activity.

The individual estimates of R_c¹ show a comparable inter-quartile range from 127 to 526 s m⁻¹ and from 148 to 1182 s m⁻¹ for positions *near/low* and *near/middle* and, similarly, for positions *near/top* and *far* a comparable inter-quartile range from 35 to 501 s m⁻¹, as well as from 26 to 314 s m⁻¹. These wide ranges reflect the turbulence-induced variability and do not indicate a physical variability of the deposition processes from one interval to the next.

Furthermore, the assumption of a constant deposition velocity over time and space can be questioned. Close to the source, the elevated deposition momentarily enhances the NH₃ loading of plant and soil surfaces, which may lead to a non-zero surface equilibrium concentration. This translates physically into an increased effective surface resistance, which will induce a negative feedback on further deposition over time. Further away from the source, this surface charging-up effect is less pronounced due to a less intense exposure of the absorbing surface to NH₃. As a result, the deposition velocity is expected to rise with the increasing measurement distance to the source, but the effect is not accounted for in the model.

The canopy resistance also depends on the physical and chemical conditions of the absorbing surfaces, e.g., water content, surface temperature or pH, possibly driven by atmospheric acid gas (SO₂, HNO₃, HCl, etc.) and aerosol concentrations [30]. The correlation between the recovery rate and the measured relative humidity, air temperature and SO₂ concentration can be seen (Supplement S3), however, the effects are rather weak and possibly attenuate each other.

After the interruption of the gas flow through the source, the NH₃ concentration did not instantaneously drop to background concentration, in contrast to the immediate response of the CH₄ concentration observed in a different, comparable release experiment (data not shown). As the miniDOAS instruments measure NH₃, NO, and SO₂ open-path, their response to the actual concentration in the air is instantaneous. Therefore, the delay in the NH₃ concentration decrease is not an instrument artefact, but points towards a temporal storage of NH₃ during the release phase [33]. This temporal storage is not considered by the deposition modelling where NH₃ is assumed to absorb irreversibly. Therefore, the modelled values of R_c likely represent lower bounds of R_c.

For any source that consists of small emission hot spots such as the artificial source with 36 point sources (or to be accurate, 36 tiny source areas) at ground level, the fraction of the emission that is deposited very close to the source is high. If we compare this fraction between a setup with the artificial
source and a comparable setup with a homogeneously emitting circular source, the model shows that measurements very close to, or even above the artificial source, will be significantly different for the two model runs (Figure 9). It can then be surmised that in the case of, say, a circular slurry patch in actual field emission experiments, the fraction of dry deposited NH₃ between the source and receptors with the same geometry to the total emitted NH₃ would be much smaller (by a factor of more than three in the example of Figure 9). Conversely, our multiple hotspot artificial source can be compared with a cattle-grazed or sheep-grazed field, in which very significant dry deposition will occur over the unsoiled grass patches in between urine hotspots [34].

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Modelled cumulative, cross-wind integrated fraction of emission (in %) that is deposited up to a distance downwind of the source, in the case of a multiple-orifice artificial source versus the case of a true circular area source (e.g., homogeneous slurry patch) of the same size. The atmospheric conditions were taken from measurement interval 6 (11h30 to 11h40). A model deposition velocity $v_d^*$ of 0.35 cm s⁻¹ was used.</figcaption></figure>

Loubet et al. [35] investigated different release experiments on grassland sites and found cumulative deposition of the released NH₃ within 200 m to be between 4% and 34%. They exhibited similar findings, with largest deposition flux close to the source and a (rather) fast drop in the deposition flux with distance. These findings are consistent with the expected deposition fluxes due to the vertical mixing of the plume with time downwind of a (surface) source.

**5. Conclusions**

Micrometeorological flux measurements always produce net flux estimates. In the case of surface-interacting trace gases like NH₃, a significant part can be deposited downwind of a source. For the interpretation of the flux estimates, it is, hence, important to distinguish between:

* the net emission from the footprint-related area (comprising the field of interest as well as the surrounding area);

* the net emission from confined areas (the field of interest with well-defined boundaries); and

* the net emission from hot spots within such an area (e.g., urine patches in a pasture field as investigated in Bell et al. [34]).

In the latter two cases, when emission measurements are performed downwind of confined source areas such as e.g., pasture areas or farm buildings which include livestock housings or manure stores,
the resulting emission rates excluding dry deposition removal will underestimate the gross emissions, if the target is to evaluate emissions from the individual hot spots such as individual buildings. The estimated emission will represent the net flux from the area defined by the source area (net emission of NH₃) plus the area between source and receptor (deposition of NH₃). Including a dry deposition algorithm in the dispersion modelling is needed to correct for this systematic error. This applies not only for emission estimation in inverse dispersion modelling based on concentration measurements, but also for flux-footprint calculations as used in the eddy-covariance and the flux-gradient method, where the net flux from the entire footprint area is estimated from vertical flux measurements.

Furthermore, the results from this study show that for upscaling the emission from urine patch models (e.g., [36,37]), dry deposition inclusion is mandatory when comparing the modelled net emission from an entire pasture area to flux measurements.

The values of the dry deposition related parameters $R_c$ and $v_d^*$ estimated in this study are specific to the surface (i.e., canopy) conditions prevailing during the experiment (e.g., a grass canopy that has been treated with glyphosate and, thus, has no stomatal activity; low relative humidity; moderate temperatures; etc.) and, in the case of $v_d^*$, additionally on the state of the flow close to the surface (via the extent of the pseudo-laminar layer). Our recommendation for other/future inverse dispersion modelling cases is to use published models or parameterizations of $R_c$ available in the literature (see for example Flechard et al. [38] and references therein, in which different models are being compared). Such readily available models can be used to calculate a range of $R_c$ estimates from the relevant environmental variables and the range of $R_c$ thus obtained by different modelling approaches will provide a range of plausible estimates of the deposited NH₃ between the source and receptor.

**Supplementary Materials:** Supplement S1 (zip): R package *bLSmodelR* ('bLSmodelR_4.0-1.tar.gz', 'bLSmodelR_4.0-1.zip', 'bLSmodelR-manual.pdf', 'Guide2bLSmodelR_v4.r'), processed field data ('Field Data and Geometry.xlsx', 'Field_Data.csv', 'Geometry.csv') and R script to follow the data analysis ('Data_Analysis.r'); Supplement S2 (pdf): CH₄ Release Experiment; Supplement S3 (pdf): Dependence of the Canopy Resistance on Environmental Parameters.

**Acknowledgments:** We thank Markus Jocher (Agroscope INH, Zürich, Switzerland) for technical support with the equipment on site, Bernhard Märki and Alfred Burri (Witzwil farm, Gampelen, Switzerland) for providing the experimental site and their support during the measurement period. The Swiss Federal Office for the Environment FOEN (contract 06.9115.PZ/L225-3604) and the Swiss Federal Office for Agriculture FOAG (contract 09.0084.PJ/K233-1881) are acknowledged for funding this study.

**Author Contributions:** C.H. designed and wrote the paper. C.F., A.N., and T.K. substantially contributed to the paper writing. C.H. wrote the *bLSmodelR* code and executed the dispersion modelling and deposition modelling analysis. T.K. and A.N. conceived and designed the experiment; J.S. and C.H. constructed and tested the artificial source; J.S., C.H., T.K., and A.N. performed the experiment; J.S. and C.H. processed the field data; C.F. contributed to the deposition analysis; A.N. and J.S. conceived, designed, and performed the experiment presented in Supplement S2; and J.S. processed and analyzed the data presented in Supplement S2.

**Conflicts of Interest:** The authors declare no conflict of interest.
Appendix A

<table><caption>**Table A1.** Summary of the bLS input parameters, concentration averages and modelled results for each instrument location and each 10-min interval. $u_*$: friction velocity; $z_0$: roughness length; $L$: Obukhov length; $\sigma_x/u_*$: scaled standard deviation of the along-wind ($u$), crosswind ($v$) and vertical ($w$); at Sonic height 1.25 m) wind velocity component; $\beta$: wind direction; $C_0$: Kolmogorov coefficient; $R_b$: pseudo-laminar boundary layer resistance; $\{C/Q\}_{sim}$: modelled C/Q ratio; $C_m$: measured NH$_3$ concentration; $C_b$: NH$_3$ background concentration; $Q_{bLS}/Q$: recovery rate; $R_c^1$: canopy resistance that provided estimates of $Q_{bLS}/Q$ closest to 1 (see Section 2.4.7.); $v_d^*$: corresponding surface deposition velocity at height $z_0$; fetch: distance to center of the source.</caption><thead><tr><th>Interval</th><th>Time Start (GMT + 1)</th><th>u* (m s<sup>-1</sup>)</th><th>z<sub>0</sub> (m)</th><th>L (m)</th><th>σ<sub>w</sub>/u* (-)</th><th>σ<sub>v</sub>/u* (-)</th><th>σ<sub>w</sub>/u* (-)</th><th>β (°)</th><th>C<sub>0</sub> (-)</th><th>R<sub>b</sub> (s m<sup>-1</sup>)</th><th>{C/Q}<sub>sim</sub> (s m<sup>-1</sup>)</th><th>C<sub>m</sub> - C<sub>b</sub> (µg m<sup>-3</sup>)</th><th>Q<sub>bLS</sub>/Q (-)</th><th>R<sub>c</sub><sup>1</sup> (s m<sup>-1</sup>)</th><th>v<sub>d</sub><sup>*</sup> (cm s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="17">near/bottom, fetch: 15 m, height: 0.5 m above ground level, path length (one way): 36 m</td></tr><tr><td>nb_1</td><td>10h40</td><td>0.17</td><td>0.004</td><td>-2.8</td><td>4.2</td><td>5.1</td><td>1.70</td><td>68</td><td>4.7</td><td>15</td><td>0.00218</td><td>35.7</td><td>0.73</td><td>81</td><td>1.0</td></tr><tr><td>nb_2</td><td>10h50</td><td>0.20</td><td>0.003</td><td>-4.9</td><td>3.9</td><td>3.3</td><td>1.42</td><td>67</td><td>4.0</td><td>13</td><td>0.00242</td><td>40.6</td><td>0.75</td><td>79</td><td>1.1</td></tr><tr><td>nb_3</td><td>11h00</td><td>0.20</td><td>0.007</td><td>-4.5</td><td>3.9</td><td>4.9</td><td>1.50</td><td>76</td><td>4.3</td><td>15</td><td>0.00213</td><td>43.4</td><td>0.92</td><td>399</td><td>0.2</td></tr><tr><td>nb_4</td><td>11h10</td><td>0.22</td><td>0.013</td><td>-4.2</td><td>3.7</td><td>4.4</td><td>1.34</td><td>60</td><td>3.6</td><td>16</td><td>0.00216</td><td>45.7</td><td>0.95</td><td>734</td><td>0.1</td></tr><tr><td>nb_5</td><td>11h20</td><td>0.25</td><td>0.016</td><td>-5.7</td><td>4.0</td><td>4.4</td><td>1.31</td><td>67</td><td>3.7</td><td>16</td><td>0.00196</td><td>43.1</td><td>0.99</td><td>2702</td><td>0.0</td></tr><tr><td>nb_6</td><td>11h30</td><td>0.23</td><td>0.024</td><td>-4.9</td><td>4.6</td><td>4.9</td><td>1.19</td><td>76</td><td>3.2</td><td>18</td><td>0.00192</td><td>38.1</td><td>0.89</td><td>347</td><td>0.3</td></tr><tr><td>nb_7</td><td>11h40</td><td>0.19</td><td>0.005</td><td>-2.7</td><td>4.3</td><td>3.9</td><td>1.60</td><td>73</td><td>4.2</td><td>15</td><td>0.00223</td><td>39.1</td><td>0.79</td><td>114</td><td>0.8</td></tr><tr><td>nb_8</td><td>11h50</td><td>0.18</td><td>0.017</td><td>-5.0</td><td>3.9</td><td>3.4</td><td>1.52</td><td>74</td><td>4.5</td><td>20</td><td>0.00244</td><td>49.0</td><td>0.90</td><td>357</td><td>0.3</td></tr><tr><td>nb_9</td><td>12h00</td><td>0.22</td><td>0.018</td><td>-6.2</td><td>3.1</td><td>4.0</td><td>1.25</td><td>103</td><td>3.5</td><td>18</td><td>0.00154</td><td>26.9</td><td>0.79</td><td>101</td><td>0.8</td></tr><tr><td>nb_10</td><td>12h10</td><td>0.17</td><td>0.002</td><td>-1.9</td><td>4.2</td><td>5.5</td><td>1.90</td><td>97</td><td>5.0</td><td>13</td><td>0.00151</td><td>39.3</td><td>1.17</td><td>∞</td><td>0.0</td></tr><tr><td>nb_11</td><td>12h20</td><td>0.16</td><td>0.004</td><td>-1.3</td><td>4.6</td><td>6.9</td><td>1.86</td><td>82</td><td>4.1</td><td>15</td><td>0.00176</td><td>39.7</td><td>1.02</td><td>∞</td><td>0.0</td></tr><tr><td>nb_12</td><td>12h30</td><td>0.20</td><td>0.010</td><td>-2.7</td><td>4.5</td><td>4.4</td><td>1.51</td><td>60</td><td>3.8</td><td>17</td><td>0.00202</td><td>37.2</td><td>0.83</td><td>166</td><td>0.5</td></tr><tr><td>nb_13</td><td>12h40</td><td>0.17</td><td>0.008</td><td>-1.8</td><td>4.1</td><td>6.8</td><td>1.63</td><td>86</td><td>3.7</td><td>18</td><td>0.00173</td><td>34.0</td><td>0.88</td><td>320</td><td>0.3</td></tr><tr><td>nb_14</td><td>12h50</td><td>0.24</td><td>0.021</td><td>-6.1</td><td>3.4</td><td>3.2</td><td>1.15</td><td>45</td><td>3.2</td><td>17</td><td>0.00250</td><td>46.5</td><td>0.84</td><td>167</td><td>0.5</td></tr><tr><td>nb_15</td><td>13h00</td><td>0.23</td><td>0.006</td><td>-4.1</td><td>3.4</td><td>3.0</td><td>1.32</td><td>58</td><td>3.5</td><td>13</td><td>0.00236</td><td>48.5</td><td>0.92</td><td>425</td><td>0.2</td></tr><tr><td>nb_16</td><td>13h10</td><td>0.16</td><td>0.002</td><td>-1.7</td><td>4.8</td><td>5.1</td><td>1.63</td><td>73</td><td>3.7</td><td>14</td><td>0.00245</td><td>40.3</td><td>0.74</td><td>107</td><td>0.8</td></tr><tr><td>nb_17</td><td>13h20</td><td>0.22</td><td>0.011</td><td>-3.6</td><td>4.1</td><td>4.1</td><td>1.48</td><td>94</td><td>4.0</td><td>16</td><td>0.00157</td><td>30.3</td><td>0.87</td><td>207</td><td>0.4</td></tr><tr><td>nb_18</td><td>13h30</td><td>0.26</td><td>0.015</td><td>-7.1</td><td>3.8</td><td>3.5</td><td>1.16</td><td>65</td><td>3.2</td><td>15</td><td>0.00228</td><td>47.6</td><td>0.94</td><td>560</td><td>0.2</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Interval</th><th>Time Start (GMT + 1)</th><th>u* (m s<sup>-1</sup>)</th><th>z<sub>0</sub> (m)</th><th>L (m)</th><th>σ<sub>u</sub>/u* (-)</th><th>σ<sub>v</sub>/u* (-)</th><th>σ<sub>w</sub>/u* (-)</th><th>β (°)</th><th>C<sub>0</sub> (-)</th><th>R<sub>b</sub> (s m<sup>-1</sup>)</th><th>{C/Q}<sub>sim</sub> (s m<sup>-1</sup>)</th><th>C<sub>m</sub> - C<sub>b</sub> (µg m<sup>-3</sup>)</th><th>Q<sub>bLS</sub>/Q (-)</th><th>R<sub>c</sub><sup>1</sup> (s m<sup>-1</sup>)</th><th>v<sub>d</sub><sup>s</sup> (cm s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="17">near/middle, fetch: 15 m, height: 1.25 m above ground level, path length (one way): 36 m</td></tr><tr><td>nm_1</td><td>10h40</td><td>0.17</td><td>0.004</td><td>-2.8</td><td>4.2</td><td>5.1</td><td>1.70</td><td>68</td><td>4.7</td><td>15</td><td>0.00124</td><td>20.7</td><td>0.75</td><td>63</td><td>1.3</td></tr><tr><td>nm_2</td><td>10h50</td><td>0.20</td><td>0.003</td><td>-4.9</td><td>3.9</td><td>3.3</td><td>1.42</td><td>67</td><td>4.0</td><td>13</td><td>0.00108</td><td>20.3</td><td>0.84</td><td>124</td><td>0.7</td></tr><tr><td>nm_3</td><td>11h00</td><td>0.20</td><td>0.007</td><td>-4.5</td><td>3.9</td><td>4.9</td><td>1.50</td><td>76</td><td>4.3</td><td>15</td><td>0.00110</td><td>23.8</td><td>0.98</td><td>1268</td><td>0.1</td></tr><tr><td>nm_4</td><td>11h10</td><td>0.22</td><td>0.013</td><td>-4.2</td><td>3.7</td><td>4.4</td><td>1.34</td><td>60</td><td>3.6</td><td>16</td><td>0.00132</td><td>26.2</td><td>0.90</td><td>239</td><td>0.4</td></tr><tr><td>nm_5</td><td>11h20</td><td>0.25</td><td>0.016</td><td>-5.7</td><td>4.0</td><td>4.4</td><td>1.31</td><td>67</td><td>3.7</td><td>16</td><td>0.00112</td><td>25.5</td><td>1.02</td><td>∞</td><td>0.0</td></tr><tr><td>nm_6</td><td>11h30</td><td>0.23</td><td>0.024</td><td>-4.9</td><td>4.6</td><td>4.9</td><td>1.19</td><td>76</td><td>3.2</td><td>18</td><td>0.00110</td><td>23.6</td><td>0.96</td><td>925</td><td>0.1</td></tr><tr><td>nm_7</td><td>11h40</td><td>0.19</td><td>0.005</td><td>-2.7</td><td>4.3</td><td>3.9</td><td>1.60</td><td>73</td><td>4.2</td><td>15</td><td>0.00120</td><td>22.7</td><td>0.85</td><td>161</td><td>0.6</td></tr><tr><td>nm_8</td><td>11h50</td><td>0.18</td><td>0.017</td><td>-5.0</td><td>3.9</td><td>3.4</td><td>1.52</td><td>74</td><td>4.5</td><td>20</td><td>0.00155</td><td>27.2</td><td>0.79</td><td>80</td><td>1.0</td></tr><tr><td>nm_9</td><td>12h00</td><td>0.22</td><td>0.018</td><td>-6.2</td><td>3.1</td><td>4.0</td><td>1.25</td><td>103</td><td>3.5</td><td>18</td><td>0.00098</td><td>15.4</td><td>0.70</td><td>38</td><td>1.8</td></tr><tr><td>nm_10</td><td>12h10</td><td>0.17</td><td>0.002</td><td>-1.9</td><td>4.2</td><td>5.5</td><td>1.90</td><td>97</td><td>5.0</td><td>13</td><td>0.00084</td><td>22.7</td><td>1.21</td><td>∞</td><td>0.0</td></tr><tr><td>nm_11</td><td>12h20</td><td>0.16</td><td>0.004</td><td>-1.3</td><td>4.6</td><td>6.9</td><td>1.86</td><td>82</td><td>4.1</td><td>15</td><td>0.00103</td><td>22.9</td><td>1.00</td><td>9126</td><td>0.0</td></tr><tr><td>nm_12</td><td>12h30</td><td>0.20</td><td>0.010</td><td>-2.7</td><td>4.5</td><td>4.4</td><td>1.51</td><td>60</td><td>3.8</td><td>17</td><td>0.00119</td><td>22.2</td><td>0.84</td><td>150</td><td>0.6</td></tr><tr><td>nm_13</td><td>12h40</td><td>0.17</td><td>0.008</td><td>-1.8</td><td>4.1</td><td>6.8</td><td>1.63</td><td>86</td><td>3.7</td><td>18</td><td>0.00103</td><td>21.9</td><td>0.96</td><td>881</td><td>0.1</td></tr><tr><td>nm_14</td><td>12h50</td><td>0.24</td><td>0.021</td><td>-6.1</td><td>3.4</td><td>3.2</td><td>1.15</td><td>45</td><td>3.2</td><td>17</td><td>0.00145</td><td>27.8</td><td>0.87</td><td>172</td><td>0.5</td></tr><tr><td>nm_15</td><td>13h00</td><td>0.23</td><td>0.006</td><td>-4.1</td><td>3.4</td><td>3.0</td><td>1.32</td><td>58</td><td>3.5</td><td>13</td><td>0.00112</td><td>25.9</td><td>1.04</td><td>∞</td><td>0.0</td></tr><tr><td>nm_16</td><td>13h10</td><td>0.16</td><td>0.002</td><td>-1.7</td><td>4.8</td><td>5.1</td><td>1.63</td><td>73</td><td>3.7</td><td>14</td><td>0.00120</td><td>22.5</td><td>0.84</td><td>206</td><td>0.5</td></tr><tr><td>nm_17</td><td>13h20</td><td>0.22</td><td>0.011</td><td>-3.6</td><td>4.1</td><td>4.1</td><td>1.48</td><td>94</td><td>4.0</td><td>16</td><td>0.00096</td><td>19.5</td><td>0.92</td><td>274</td><td>0.3</td></tr><tr><td>nm_18</td><td>13h30</td><td>0.26</td><td>0.015</td><td>-7.1</td><td>3.8</td><td>3.5</td><td>1.16</td><td>65</td><td>3.2</td><td>15</td><td>0.00124</td><td>23.4</td><td>0.85</td><td>147</td><td>0.6</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Interval</th><th>Time Start (GMT + 1)</th><th>u* (m s<sup>-1</sup>)</th><th>z<sub>0</sub> (m)</th><th>L (m)</th><th>σ<sub>u</sub>/u* (-)</th><th>σ<sub>v</sub>/u* (-)</th><th>σ<sub>w</sub>/u* (-)</th><th>β (°)</th><th>C<sub>0</sub> (-)</th><th>R<sub>b</sub> (s m<sup>-1</sup>)</th><th>{C/Q}<sub>sim</sub> (s m<sup>-1</sup>)</th><th>C<sub>m</sub> - C<sub>b</sub> (µg m<sup>-3</sup>)</th><th>Q<sub>bLS</sub>/Q (-)</th><th>R<sub>c</sub><sup>1</sup> (s m<sup>-1</sup>)</th><th>v<sub>d</sub><sup>*</sup> (cm s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="17">near/top, fetch: 15 m, height: 3.0 m above ground level, path length (one way): 37 m</td></tr><tr><td>nt_1</td><td>10h40</td><td>0.17</td><td>0.004</td><td>-2.8</td><td>4.2</td><td>5.1</td><td>1.70</td><td>68</td><td>4.7</td><td>15</td><td>0.00033</td><td>4.9</td><td>0.66</td><td>30</td><td>2.2</td></tr><tr><td>nt_2</td><td>10h50</td><td>0.20</td><td>0.003</td><td>-4.9</td><td>3.9</td><td>3.3</td><td>1.42</td><td>67</td><td>4.0</td><td>13</td><td>0.00021</td><td>3.4</td><td>0.73</td><td>35</td><td>2.1</td></tr><tr><td>nt_3</td><td>11h00</td><td>0.20</td><td>0.007</td><td>-4.5</td><td>3.9</td><td>4.9</td><td>1.50</td><td>76</td><td>4.3</td><td>15</td><td>0.00026</td><td>5.2</td><td>0.88</td><td>149</td><td>0.6</td></tr><tr><td>nt_4</td><td>11h10</td><td>0.22</td><td>0.013</td><td>-4.2</td><td>3.7</td><td>4.4</td><td>1.34</td><td>60</td><td>3.6</td><td>16</td><td>0.00038</td><td>6.4</td><td>0.75</td><td>63</td><td>1.3</td></tr><tr><td>nt_5</td><td>11h20</td><td>0.25</td><td>0.016</td><td>-5.7</td><td>4.0</td><td>4.4</td><td>1.31</td><td>67</td><td>3.7</td><td>16</td><td>0.00031</td><td>7.7</td><td>1.10</td><td>∞</td><td>0.0</td></tr><tr><td>nt_6</td><td>11h30</td><td>0.23</td><td>0.024</td><td>-4.9</td><td>4.6</td><td>4.9</td><td>1.19</td><td>76</td><td>3.2</td><td>18</td><td>0.00034</td><td>6.9</td><td>0.91</td><td>338</td><td>0.3</td></tr><tr><td>nt_7</td><td>11h40</td><td>0.19</td><td>0.005</td><td>-2.7</td><td>4.3</td><td>3.9</td><td>1.60</td><td>73</td><td>4.2</td><td>15</td><td>0.00035</td><td>5.8</td><td>0.74</td><td>51</td><td>1.5</td></tr><tr><td>nt_8</td><td>11h50</td><td>0.18</td><td>0.017</td><td>-5.0</td><td>3.9</td><td>3.4</td><td>1.52</td><td>74</td><td>4.5</td><td>20</td><td>0.00052</td><td>6.2</td><td>0.54</td><td>0</td><td>4.9</td></tr><tr><td>nt_9</td><td>12h00</td><td>0.22</td><td>0.018</td><td>-6.2</td><td>3.1</td><td>4.0</td><td>1.25</td><td>103</td><td>3.5</td><td>18</td><td>0.00025</td><td>2.6</td><td>0.46</td><td>0</td><td>5.7</td></tr><tr><td>nt_10</td><td>12h10</td><td>0.17</td><td>0.002</td><td>-1.9</td><td>4.2</td><td>5.5</td><td>1.90</td><td>97</td><td>5.0</td><td>13</td><td>0.00023</td><td>5.4</td><td>1.04</td><td>∞</td><td>0.0</td></tr><tr><td>nt_11</td><td>12h20</td><td>0.16</td><td>0.004</td><td>-1.3</td><td>4.6</td><td>6.9</td><td>1.86</td><td>82</td><td>4.1</td><td>15</td><td>0.00034</td><td>7.0</td><td>0.93</td><td>555</td><td>0.2</td></tr><tr><td>nt_12</td><td>12h30</td><td>0.20</td><td>0.010</td><td>-2.7</td><td>4.5</td><td>4.4</td><td>1.51</td><td>60</td><td>3.8</td><td>17</td><td>0.00041</td><td>8.4</td><td>0.91</td><td>275</td><td>0.3</td></tr><tr><td>nt_13</td><td>12h40</td><td>0.17</td><td>0.008</td><td>-1.8</td><td>4.1</td><td>6.8</td><td>1.63</td><td>86</td><td>3.7</td><td>18</td><td>0.00033</td><td>7.7</td><td>1.06</td><td>∞</td><td>0.0</td></tr><tr><td>nt_14</td><td>12h50</td><td>0.24</td><td>0.021</td><td>-6.1</td><td>3.4</td><td>3.2</td><td>1.15</td><td>45</td><td>3.2</td><td>17</td><td>0.00043</td><td>8.1</td><td>0.84</td><td>118</td><td>0.7</td></tr><tr><td>nt_15</td><td>13h00</td><td>0.23</td><td>0.006</td><td>-4.1</td><td>3.4</td><td>3.0</td><td>1.32</td><td>58</td><td>3.5</td><td>13</td><td>0.00024</td><td>6.2</td><td>1.15</td><td>∞</td><td>0.0</td></tr><tr><td>nt_16</td><td>13h10</td><td>0.16</td><td>0.002</td><td>-1.7</td><td>4.8</td><td>5.1</td><td>1.63</td><td>73</td><td>3.7</td><td>14</td><td>0.00035</td><td>6.8</td><td>0.88</td><td>283</td><td>0.3</td></tr><tr><td>nt_17</td><td>13h20</td><td>0.22</td><td>0.011</td><td>-3.6</td><td>4.1</td><td>4.1</td><td>1.48</td><td>94</td><td>4.0</td><td>16</td><td>0.00031</td><td>5.0</td><td>0.72</td><td>35</td><td>2.0</td></tr><tr><td>nt_18</td><td>13h30</td><td>0.26</td><td>0.015</td><td>-7.1</td><td>3.8</td><td>3.5</td><td>1.16</td><td>65</td><td>3.2</td><td>15</td><td>0.00029</td><td>4.6</td><td>0.70</td><td>32</td><td>2.1</td></tr></tbody></table>
<table><caption>Table A1. *Cont.*</caption><thead><tr><th>Interval</th><th>Time Start (GMT + 1)</th><th>u* (m s<sup>-1</sup>)</th><th>z<sub>0</sub> (m)</th><th>L (m)</th><th>σ<sub>u</sub>/u* (-)</th><th>σ<sub>v</sub>/u* (-)</th><th>σ<sub>v</sub>/u* (-)</th><th>β (°)</th><th>C<sub>0</sub> (-)</th><th>R<sub>b</sub> (s m<sup>-1</sup>)</th><th>{C/Q}<sub>sim</sub> (s m<sup>-1</sup>)</th><th>C<sub>m</sub> - C<sub>b</sub> (µg m<sup>-3</sup>)</th><th>Q<sub>bLS</sub>/Q (-)</th><th>R<sub>c</sub><sup>1</sup> (s m<sup>-1</sup>)</th><th>v<sub>d</sub><sup>*</sup> (cm s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="17">far, fetch: 70 m, height: 1.2 m above ground level, path length (one way): 33 m</td></tr><tr><td>f_1</td><td>10h40</td><td>0.17</td><td>0.004</td><td>-2.8</td><td>4.2</td><td>5.1</td><td>1.70</td><td>68</td><td>4.7</td><td>15</td><td>0.00020</td><td>2.8</td><td>0.64</td><td>56</td><td>1.4</td></tr><tr><td>f_2</td><td>10h50</td><td>0.20</td><td>0.003</td><td>-4.9</td><td>3.9</td><td>3.3</td><td>1.42</td><td>67</td><td>4.0</td><td>13</td><td>0.00033</td><td>3.7</td><td>0.51</td><td>17</td><td>3.3</td></tr><tr><td>f_3</td><td>11h00</td><td>0.20</td><td>0.007</td><td>-4.5</td><td>3.9</td><td>4.9</td><td>1.50</td><td>76</td><td>4.3</td><td>15</td><td>0.00018</td><td>3.7</td><td>0.93</td><td>548</td><td>0.2</td></tr><tr><td>f_4</td><td>11h10</td><td>0.22</td><td>0.013</td><td>-4.2</td><td>3.7</td><td>4.4</td><td>1.34</td><td>60</td><td>3.6</td><td>16</td><td>0.00021</td><td>2.5</td><td>0.54</td><td>31</td><td>2.1</td></tr><tr><td>f_5</td><td>11h20</td><td>0.25</td><td>0.016</td><td>-5.7</td><td>4.0</td><td>4.4</td><td>1.31</td><td>67</td><td>3.7</td><td>16</td><td>0.00019</td><td>3.3</td><td>0.77</td><td>121</td><td>0.7</td></tr><tr><td>f_6</td><td>11h30</td><td>0.23</td><td>0.024</td><td>-4.9</td><td>4.6</td><td>4.9</td><td>1.19</td><td>76</td><td>3.2</td><td>18</td><td>0.00015</td><td>2.3</td><td>0.72</td><td>109</td><td>0.8</td></tr><tr><td>f_7</td><td>11h40</td><td>0.19</td><td>0.005</td><td>-2.7</td><td>4.3</td><td>3.9</td><td>1.60</td><td>73</td><td>4.2</td><td>15</td><td>0.00017</td><td>2.1</td><td>0.56</td><td>26</td><td>2.5</td></tr><tr><td>f_8</td><td>11h50</td><td>0.18</td><td>0.017</td><td>-5.0</td><td>3.9</td><td>3.4</td><td>1.52</td><td>74</td><td>4.5</td><td>20</td><td>0.00020</td><td>2.5</td><td>0.57</td><td>22</td><td>2.4</td></tr><tr><td>f_9</td><td>12h00</td><td>0.22</td><td>0.018</td><td>-6.2</td><td>3.1</td><td>4.0</td><td>1.25</td><td>103</td><td>3.5</td><td>18</td><td>0.00003</td><td>0.6</td><td>0.94</td><td>314</td><td>0.3</td></tr><tr><td>f_10<sup>1</sup></td><td>12h10</td><td>0.17</td><td>0.002</td><td>-1.9</td><td>4.2</td><td>5.5</td><td>1.90</td><td>97</td><td>5.0</td><td>13</td><td>0.00003</td><td>2.6</td><td>4.37<sup>1</sup></td><td>∞<sup>1</sup></td><td>0.0<sup>1</sup></td></tr><tr><td>f_11</td><td>12h20</td><td>0.16</td><td>0.004</td><td>-1.3</td><td>4.6</td><td>6.9</td><td>1.86</td><td>82</td><td>4.1</td><td>15</td><td>0.00009</td><td>2.1</td><td>1.05</td><td>∞</td><td>0.0</td></tr><tr><td>f_12</td><td>12h30</td><td>0.20</td><td>0.010</td><td>-2.7</td><td>4.5</td><td>4.4</td><td>1.51</td><td>60</td><td>3.8</td><td>17</td><td>0.00019</td><td>2.8</td><td>0.67</td><td>78</td><td>1.1</td></tr><tr><td>f_13</td><td>12h40</td><td>0.17</td><td>0.008</td><td>-1.8</td><td>4.1</td><td>6.8</td><td>1.63</td><td>86</td><td>3.7</td><td>18</td><td>0.00009</td><td>1.8</td><td>0.87</td><td>318</td><td>0.3</td></tr><tr><td>f_14</td><td>12h50</td><td>0.24</td><td>0.021</td><td>-6.1</td><td>3.4</td><td>3.2</td><td>1.15</td><td>45</td><td>3.2</td><td>17</td><td>0.00023</td><td>2.2</td><td>0.44</td><td>6</td><td>4.3</td></tr><tr><td>f_15</td><td>13h00</td><td>0.23</td><td>0.006</td><td>-4.1</td><td>3.4</td><td>3.0</td><td>1.32</td><td>58</td><td>3.5</td><td>13</td><td>0.00039</td><td>3.4</td><td>0.39</td><td>0</td><td>7.2</td></tr><tr><td>f_16</td><td>13h10</td><td>0.16</td><td>0.002</td><td>-1.7</td><td>4.8</td><td>5.1</td><td>1.63</td><td>73</td><td>3.7</td><td>14</td><td>0.00017</td><td>2.7</td><td>0.71</td><td>114</td><td>0.8</td></tr><tr><td>f_17</td><td>13h20</td><td>0.22</td><td>0.011</td><td>-3.6</td><td>4.1</td><td>4.1</td><td>1.48</td><td>94</td><td>4.0</td><td>16</td><td>0.00004</td><td>0.9</td><td>0.93</td><td>380</td><td>0.3</td></tr><tr><td>f_18</td><td>13h30</td><td>0.26</td><td>0.015</td><td>-7.1</td><td>3.8</td><td>3.5</td><td>1.16</td><td>65</td><td>3.2</td><td>15</td><td>0.00029</td><td>3.5</td><td>0.54</td><td>30</td><td>2.2</td></tr></tbody></table>

<footnote n="1">¹ interval f_10 was removed from the analysis.</footnote>
References

1. Coates, T.W.; Flesch, T.K.; McGinn, S.M.; Charmley, E.; Chen, D. Evaluating an eddy covariance technique to estimate point-source emissions and its potential application to grazing cattle. *Agric. For. Meteorol.* **2017**, *234–235*, 164–171. [CrossRef]

2. Felber, R.; Münger, A.; Neftel, A.; Ammann, C. Eddy covariance methane flux measurements over a grazed pasture: Effect of cows as moving point sources. *Biogeosciences* **2015**, *12*, 3419–3468. [CrossRef]

3. Laubach, J.; Kelliher, F.M.; Knight, T.W.; Clark, H.; Molano, G.; Cavanagh, A. Methane emissions from beef cattle—A comparison of paddock- and animal-scale measurements. *Aust. J. Exp. Agric.* **2008**, *48*, 132–137. [CrossRef]

4. Laubach, J.; Taghizadeh-Toosi, A.; Sherlock, R.R.; Kelliher, F.M. Measuring and modelling ammonia emissions from a regular pattern of cattle urine patches. *Agric. For. Meteorol.* **2012**, *156*, 1–17. [CrossRef]

5. Loubet, B.; Génermont, S.; Ferrara, R.M.; Bedos, C.; Decuq, C.; Personne, E.; Fanucci, O.; Durand, B.; Rana, G.; Cellier, P. An inverse model to estimate ammonia emissions from fields. *Eur. J. Soil Sci.* **2010**, *61*, 793–805. [CrossRef]

6. Sintermann, J.; Ammann, C.; Kuhn, U.; Spirig, C.; Hirschberger, R.; Gartner, A.; Neftel, A. Determination of field scale ammonia emissions for common slurry spreading practice with two independent methods. *Atmos. Meas. Tech.* **2011**, *4*, 1821–1840. [CrossRef]

7. Wilson, J.D.; Flesch, T.K.; Crenna, B.P. Estimating Surface-Air Gas Fluxes by Inverse Dispersion Using a Backward Lagrangian Stochastic Trajectory Model. In *Lagrangian Modeling of the Atmosphere*; Lin, J., Brunner, D., Gerbig, C., Stohl, A., Luhar, A., Webley, P., Eds.; American Geophysical Union: Washington, DC, USA, 2012; pp. 149–162.

8. Flesch, T.K.; Wilson, J.D.; Harper, L.A.; Crenna, B.P.; Sharpe, R.R. Deducing ground-to-air emissions from observed trace gas concentrations: A field trial. *J. Appl. Meteorol. 2004*, *43*, 487–502. [CrossRef]

9. Carozzi, M.; Loubet, B.; Acutis, M.; Rana, G.; Ferrara, R.M. Inverse dispersion modelling highlights the efficiency of slurry injection to reduce ammonia losses by agriculture in the Po Valley (Italy). *Agric. For. Meteorol.* **2013**, *171*, 306–318. [CrossRef]

10. Flesch, T.K.; Harper, L.A.; Powell, J.M.; Wilson, J.D. Inverse-dispersion calculation of ammonia emissions from Wisconsin dairy farms. *Trans. ASABE 2009*, *52*, 253–265. [CrossRef]

11. Grant, R.H.; Boehm, M.T.; Bogan, B.W. Methane and carbon dioxide emissions from manure storage facilities at two free-stall dairies. *Agric. For. Meteorol.* **2015**, *213*, 102–113. [CrossRef]

12. Harper, L.A.; Flesch, T.K.; Weaver, K.H.; Wilson, J.D. The effect of biofuel production on swine farm methane and ammonia emissions. *J. Environ. Qual.* **2010**, *39*, 1984–1992. [CrossRef] [PubMed]

13. McGinn, S.M.; Flesch, T.K.; Crenna, B.P.; Beauchemin, K.A.; Coates, T. Quantifying ammonia emissions from a cattle feedlot using a dispersion model. *J. Environ. Qual.* **2007**, *36*, 1585–1590. [CrossRef] [PubMed]

14. Flechard, C.R. Turbulent Exchange of Ammonia above Vegetation (BL). Ph.D. Thesis, University of Nottingham, Nottingham, UK, 1998.

15. Schrader, F.; Brümmer, C. Land use specific ammonia deposition velocities: A review of recent studies (2004–2013). *Water Air Soil Pollut.* **2014**, *225*, 2114. [CrossRef] [PubMed]

16. Asman, W.A.H.; Sutton, M.A.; Schjørring, J.K. Ammonia: Emission, atmospheric transport and deposition. *New Phytol.* **1998**, *139*, 27–48. [CrossRef]

17. Loubet, B.; Asman, W.A.H.; Theobald, M.R.; Hertel, O.; Tang, Y.S.; Robin, P.; Hassouna, M.; Dämmgen, U.; Génermont, S.; Cellier, P.; et al. Ammonia Deposition Near Hot Spots: Processes, Models and Monitoring Methods. In *Atmospheric Ammonia*; Sutton, M.A., Reis, S., Baker, S.M.H., Eds.; Springer: Dordrecht, The Netherlands, 2009; pp. 205–267.

18. R Core Team. *R: A language and environment for statistical computing*; R Foundation for Statistical Computing: Vienna, Austria, 2018.

19. Flesch, T.K.; Wilson, J.D.; Yee, E. Backward-time Lagrangian stochastic dispersion models and their application to estimate gaseous Emissions. *J. Appl. Meteorol. 1995*, *34*, 1320–1332. [CrossRef]

20. Wilson, J.D.; Ferrandino, F.J.; Thurtell, G.W. A relationship between deposition velocity and trajectory reflection probability for use in stochastic Lagrangian dispersion models. *Agric. For. Meteorol.* **1989**, *47*, 139–154. [CrossRef]
21. Knaus, J. snowfall: Easier Cluster Computing (Based on snow). R Package Version 1.84-6.1. Available online: https://CRAN.R-project.org/package=snowfall (accessed on 1 January 2017).

22. Luke, T.; Rossini, A.J.; Li, N.; Sevcikova, H. snow: Simple Network of Workstations. R Package Version 0.4-2. Available online: https://CRAN.R-project.org/package=snow (accessed on 1 January 2017).

23. Eddelbuettel, D.; François, R. Rcpp: Seamless R and C++ Integration. *J. Stat. Softw.* **2011**, *40*, 1–18. [CrossRef]

24. Dowle, M.; Srinivasan, A. data.table: Extension of 'data.frame'. R Package Version 1.10.4. Available online: https://CRAN.R-project.org/package=data.table (accessed on 1 January 2018).

25. Sintermann, J.; Dietrich, K.; Häni, C.; Bell, M.; Jocher, M.; Neftel, A. A miniDOAS instrument optimised for ammonia field measurements. *Atmos. Meas. Tech.* **2016**, *9*, 2721–2734. [CrossRef]

26. Sutton, M.A.; Schjorring, J.K.; Wyers, G.P.; Duyzer, J.H.; Ineson, P.; Powlson, D.S. Plant-atmosphere exchange of ammonia [and discussion]. *Philos. Trans. R. Soc. A* **1995**, *351*, 261–278. [CrossRef]

27. Garland, J.A. The dry deposition of sulphur dioxide to land and water surfaces. *Proc. R. Soc. A* **1977**, *354*, 245–268. [CrossRef]

28. Cussler, E.L. *Diffusion: Mass Transfer in Fluid Systems*, 3rd ed.; 6th printing; Cambridge University Press: Cambridge, UK, 2013.

29. Tang, M.J.; Cox, R.A.; Kalberer, M. Compilation and evaluation of gas phase diffusion coefficients of reactive trace gases in the atmosphere: Volume 1. Inorganic compounds. *Atmos. Chem. Phys.* **2014**, *14*, 9233–9247. [CrossRef]

30. Flechard, C.R.; Massad, R.S.; Loubet, B.; Personne, E.; Simpson, D.; Bash, J.O.; Cooter, E.J.; Nemitz, E.; Sutton, M.A. Advances in understanding, models and parameterizations of biosphere-atmosphere ammonia exchange. *Biogeosciences* **2013**, *10*, 5183–5225. [CrossRef]

31. Flechard, C.R.; Spirig, C.; Neftel, A.; Ammann, C. The annual ammonia budget of fertilised cut grassland—Part 2: Seasonal variations and compensation point modeling. *Biogeosciences* **2010**, *7*, 537–556. [CrossRef]

32. Massad, R.S.; Nemitz, E.; Sutton, M.A. Review and parameterisation of bi-directional ammonia exchange between vegetation and the atmosphere. *Atmos. Chem. Phys.* **2010**, *10*, 10359–10386. [CrossRef]

33. Sutton, M.A.; Burkhardt, J.K.; Guerin, D.; Nemitz, E.; Fowler, D. Development of resistance models to describe measurements of bi-directional ammonia surface-atmosphere exchange. *Atmos. Environ.* **1998**, *32*, 473–480. [CrossRef]

34. Bell, M.; Flechard, C.; Fauvel, Y.; Häni, C.; Sintermann, J.; Jocher, M.; Menzi, H.; Hensen, A.; Neftel, A. Ammonia emissions from a grazed field estimated by miniDOAS measurements and inverse dispersion modelling. *Atmos. Meas. Tech.* **2017**, *10*, 1875–1892. [CrossRef]

35. Loubet, B.; Cellier, P.; Milford, C.; Sutton, M.A. A coupled dispersion and exchange model for short-range dry deposition of atmospheric ammonia. *Q. J. R. Meteorol. Soc.* **2006**, *132*, 1733–1763. [CrossRef]

36. Giltrap, D.; Saggar, S.; Rodriguez, J.; Bishop, P. Modelling NH₃ volatilisation within a urine patch using NZ-DNDC. *Nutr. Cycl. Agroecosyst.* **2017**, *108*, 267–277. [CrossRef]

37. Móring, A.; Vieno, M.; Doherty, R.M.; Milford, C.; Nemitz, E.; Twigg, M.M.; Horváth, L.; Sutton, M.A. Process-based modelling of NH₃ exchange with grazed grasslands. *Biogeosciences* **2017**, *14*, 4161–4193. [CrossRef]

38. Flechard, C.R.; Nemitz, E.; Smith, R.I.; Fowler, D.; Vermeulen, A.T.; Bleeker, A.; Erisman, J.W.; Simpson, D.; Zhang, L.; Tang, Y.S.; et al. Dry deposition of reactive nitrogen to European ecosystems: A comparison of inferential models across the NitroEurope network. *Atmos. Chem. Phys.* **2011**, *11*, 2703–2728. [CrossRef]

<figure><img src="image_13.png" /></figure>