Evaluation of a new inference method for estimating ammonia volatilisation from multiple agronomic plots
- 1INRA, UMR ECOSYS, INRA, AgroParisTech, Université Paris-Saclay, 78850, Thiverval-Grignon, France
- 2ARVALIS-Institut du Végétal, Station expérimentale de La Jaillière, La Chapelle Saint Sauveur, 44370 Loireauxence, France
- anow at: Agroscope Research Station, Climate and Agriculture, Zurich, Switzerland
Correspondence: Benjamin Loubet (firstname.lastname@example.org)
Tropospheric ammonia (NH3) is a threat to the environment and human health and is mainly emitted by agriculture. Ammonia volatilisation following application of nitrogen in the field accounts for more than 40 % of the total NH3 emissions in France. This represents a major loss of nitrogen use efficiency which needs to be reduced by appropriate agricultural practices. In this study we evaluate a novel method to infer NH3 volatilisation from small agronomic plots consisting of multiple treatments with repetition. The method is based on the combination of a set of NH3 diffusion sensors exposed for durations of 3 h to 1 week and a short-range atmospheric dispersion model, used to retrieve the emissions from each plot. The method is evaluated by mimicking NH3 emissions from an ensemble of nine plots with a resistance analogue–compensation point–surface exchange scheme over a yearly meteorological database separated into 28-day periods. A multifactorial simulation scheme is used to test the effects of sensor numbers and heights, plot dimensions, source strengths, and background concentrations on the quality of the inference method. We further demonstrate by theoretical considerations in the case of an isolated plot that inferring emissions with diffusion sensors integrating over daily periods will always lead to underestimations due to correlations between emissions and atmospheric transfer. We evaluated these underestimations as −8 % ± 6 % of the emissions for a typical western European climate. For multiple plots, we find that this method would lead to median underestimations of −16 % with an interquartile [−8–22 %] for two treatments differing by a factor of up to 20 and a control treatment with no emissions. We further evaluate the methodology for varying background concentrations and NH3 emissions patterns and demonstrate the low sensitivity of the method to these factors. The method was also tested in a real case and proved to provide sound evaluations of NH3 losses from surface applied and incorporated slurry. We hence showed that this novel method should be robust and suitable for estimating NH3 emissions from agronomic plots. We believe that the method could be further improved by using Bayesian inference and inferring surface concentrations rather than surface fluxes. Validating against controlled source is also a remaining challenge.
Tropospheric ammonia (NH3) is mainly emitted by agriculture and has great environmental impacts (atmospheric pollution, eutrophication, reduction of biodiversity), which are increasingly taken into account in European and international regulations (Council, 1996, 2016; UNECE, 2012). Ammonia losses also have great agronomic and economic impacts for farmers, as they reduce nitrogen use efficiency. The varying prices of mineral fertilisers and concerns about environmental and health threats demand improvements in the efficiency of nitrogen utilisation, and especially in recycling nitrogen through organic fertilisation (Sutton et al., 2011). Indeed, NH3 volatilisation during storage of manure and slurry and following their field application is the main source of NH3 in Europe (55 % of the emissions) while farm building emissions represent 45 %. In France, crop farming represents 35 % of the emissions and animal farming represents 65 % (CITEPA, 2017; ECETOC, 1994; EUROSTAT, 2012; Faburé et al., 2011). Reducing NH3 losses from this agricultural sector is therefore a major objective for applied research.
While NH3 emissions from farm buildings and storage can be handled by engineering solutions, losses during organic fertilisation are much more dependent on the combination of application methods (splash plate, band spreading, pressurised injection, open and close slot injection, trailing hose, and trailing shoe), soil type and occupation, and environmental conditions (soil humidity, air temperature, wind speed, solar radiation) (Sommer et al., 2003). For instance, Sintermann et al. (2012) report NH3 losses following cattle and pig slurry application in the field ranging from a few percent to 50 % over large fields and up to 100 % over medium fields. Evaluating ammonia losses from field fertilisation over a range of practices and soil and climatic conditions is therefore key in evaluating the best application methods.
However, characterising these emissions at the field scale requires complex experimental design and most of the time also requires the use of large fields (Ferrara et al., 2016, 2012; Flechard and Fowler, 1998; Loubet et al., 2012; Milford et al., 2009; Sintermann et al., 2011b; Spirig et al., 2010; Sun et al., 2015; Whitehead et al., 2008). Especially useful for measuring ammonia losses are methods that can deal with small- and medium-scale fields (20–50 m on the side) that are commonly used in agronomic trials. Indirect estimation methods (soil nitrogen balance or 15N balance) are not well adapted to evaluate gaseous ammonia losses, mainly because of the soil heterogeneity and also because the method relies on evaluating small variations of large numbers (McGinn and Janzen, 1998). Among existing methods for measuring NH3 emissions, the integrated horizontal flux method (Wilson and Shum, 1992) is well adapted, but is a subject of debate in its practical application since it seems to be systematically biased towards higher estimates (Häni et al., 2016; Sintermann et al., 2012). Alternatively, enclosure methods proved to be not representative for a sticky compound such as ammonia (Pacholski et al., 2006), but more concerning is the fact that ammonia fluxes result from an air-surface equilibrium which is disturbed by the confined environment offered by the chamber. Inverse dispersion modelling approaches either based on backward Lagrangian stochastic models (Flesch et al., 1995) or Eulerian models (Kormann and Meixner, 2001; Loubet et al., 2001) or based on the Philip equation (Philip, 1959) have been demonstrated to be adapted for estimating NH3 volatilisation from strong sources (Loubet et al., 2010; Sommer et al., 2005).
These approaches are well adapted to small or medium fields (≤ 50 × 50 m2) but typically require hourly NH3 concentration measurements. Long-term concentration measurements of NH3 are now well handled by the use of short-path passive samplers developed by Sutton et al. (2001), or active denuders, which have both been used for concentration monitoring for years (Tang et al., 2001, 2009). These active denuders can be adapted for measuring fluxes based on conditional sampling like the conditional time-averaged gradient method (COTAG) (Famulari et al., 2010), which is a useful method but only adapted for large fields (≥ 0.5 ha). The passive samplers have also been shown to be adapted for inverse modelling estimations of NH3 sources for large fields (Carozzi et al., 2013b; Ferrara et al., 2014).
In another field of research, solutions to the multiple source inference problem, which consists of inferring multiple sources based on measured concentrations at multiple points in space and time, have been developed especially since 2008 (Crenna et al., 2008; Gao et al., 2008; Gericke et al., 2011; Mukherjee et al., 2015; Vandré and Kaupenjohann, 1998). They have chiefly been used over regional scales (Flesch et al., 2009; Lushi and Stockie, 2010; Yee and Flesch, 2010), and have been shown to be very dependent on the source-sensor geometry (Crenna et al., 2008; Flesch et al., 2009; Wang et al., 2013). Mukherjee et al. (2015) highlighted the dependency of the inferred source on background concentration and plot disposition by means of an inverse footprint approach. Yee et al. (2008) have shown how to retrieve the number, location and intensity of multiple sources with dispersion models coupled with Bayesian inference methods. Yee and Flesch (2010) have evaluated the inversion and inference methods for determining four point sources using several laser transects. Flesch et al. (2009) have shown that source–receptor geometry is critical in determining whether a multiple-source inversion problem can provide realistic solutions or not. Flesch et al. (2009) have moreover shown that if the geometry is well chosen the accuracy of the method for a 15 min integration time can reach 10 to 20 %. These studies have also shown that the multiple source inference problems can be solved if not ill-conditioned (ill-conditioning depends on the location of sources and concentration sensors and is characterised by a conditioning number κ).
In this study, we pose the following research questions: can inverse dispersion modelling approaches be used for inferring NH3 emissions from multiple small plots (agronomic trials) using passive samplers, and to which degree of accuracy? The answer is given through the investigation of the optimal design in terms of field dimensions, plot location and size, passive sampler locations, and their duration of exposure. Throughout this study, agronomic trials are considered to be multiple small adjacent fields with repetitions of treatments. A typical trial would consist of three repetitions of three treatments. Hence the double challenge that we face in this study is to consider both (i) the multiple-source inference issue (adjacent small fields) and the (ii) time-integration issue (using passive samplers).
To answer these questions, we use a four-step approach: (1) the ammonia emissions are first modelled on each source using prescribed NH3 emissions potential dynamics coupled with a simple soil–vegetation–atmosphere exchange scheme to mimic realistic seasonal, daily and hourly variations in NH3 emissions. (2) These prescribed emissions are then used to estimate the concentration at each target location using short-range atmospheric dispersion modelling over half-hourly periods. (3) The obtained concentrations are then averaged over several integration periods to simulate the behaviour of passive samplers. Finally, (4) the sources are evaluated by inference with dispersion modelling based on the averaged concentrations.
Two dispersion models and several inference methodologies are evaluated. The effect of the size of the source, the locations of targets, the dynamics and magnitude of each source, the meteorological conditions, and the background concentration variability are evaluated and discussed. The feasibility of the method is finally evaluated over a real case with two repetitions of three treatments (slurry spreading, injection and a reference without fertilisation).
At first we present the theoretical background of source inference by optimisation for single and multiple sources with time-averaging concentration sensors. Then the method used to generate a realistic ammonia source is introduced before the description of the dispersion models used for both generating the concentration fields and inferring back the sources. The geometry of the sources, sensor locations and the meteorological data used for this analysis are then shown, and finally the real test case used for evaluating the method is detailed.
2.1 The theory of the source inference method
At first we will recall some important theoretical features of the inverse dispersion modelling approach, which is actually an inference method.
2.1.1 Case of a single area source and a single concentration sampler
We first consider the case of a single area source with a single concentration sampler (target). The source varies with time. The method is based upon the general superimposition principle (Thomson et al., 2007), which relates the concentration at a given location C(x,t) to the source strength S(t) and the background concentration Cbgd(t) using a transfer function D(x,t), which has the dimensions of a transfer resistance (s m−1).
Here x denotes the location of the sensor and t the time. The concentration and source units are in µg N-NH3 m−3 and µg N-, respectively. The superimposition principle implies that the studied tracer must be conservative, which is a reasonable hypothesis for NH3 whose reaction time with acids in the atmosphere is below the transport time for spatial scales below 1000 m (Nemitz et al., 2009). Moreover, in Eq. (1), we assume a spatially homogeneous area source with strength S(t). The spatial homogeneity of the source is less trivial for NH3 than other gas released in agriculture as the source itself depends on the concentration at the surface. However, Loubet et al. (2010) have shown that the heterogeneity of the source can be neglected as long as the dimension of the source is larger than 20 m. Hence, this study is limited to source areas with fetch larger than 20 m and a spread of the concentration samplers over a domain smaller than 1000 m. Moreover, it is interesting to note that for infinitely spread fields, the transfer resistance is linearly linked to the transfer matrix (see Supplement Sect. S1)
2.1.2 Effect of time-averaging sensors on source inference for a single source
Since we consider time-averaging concentration samplers, we develop the time-averaged equation of Eq. (1) over an integration time period τ:
where the overbars denote a time average over the period τ. Similar to turbulent flux calculations, the first part of the right-hand side of Eq. (2) is decomposed using the Reynolds decomposition of a random variable (Kaimal and Finnigan, 1994), giving
where is the time covariance between D(x,t) and S(t). If the averaged background concentration is a known quantity, Eq. (3) can be easily manipulated to give an estimation of the averaged source strength , the quantity we want to infer:
In the right-hand side of Eq. (4), (I) can be calculated from measured and and , which is itself calculated with dispersion models. Conversely, (II) is a priori unknown and depends on the correlation between the source strength and the transfer function . Hence, if (II) is neglected, the inferred source is biased. The relative bias of the method is then
Hence we show in Eq. (5) that time averaging leads to a relative bias which can be quantified by the time covariance between the transfer function and the source strength. However, this quantity is by nature unknown since the dynamics of S(t) is unknown. Determining requires knowledge of the source dynamics, which can be obtained from measurements with a micrometeorological method. It can alternatively be approached by modelling using state-of-the-art ammonia exchange processes as we do here.
In addition to the bias, which is term (II) in Eq. (4), evaluating term (I) is encompassed with errors related to the uncertainties in , and . In particular, cases in which is small may lead to large errors in inferring the source term S. This is linked to the conditioning of the inverse problem and is discussed in Supplement Sect. S2.
2.1.3 Case of multiple sources and multiple concentration samplers with time averaging
If we generalise the approach to multiple sources and multiple receptors, then the transfer function becomes a matrix , which is the contribution of source Sj to concentration at a target located at xi. For reading purposes we simplify the matrix notation to Di,j. Equation (3) then becomes
which in condensed notation gives
If the number of targets is equal to the number of sources, the problem can be solved by inversion of a linear system. If the number of targets is larger than the number of sources, the problem is a multiple linear regression type with unknowns and . The third term on the right-hand side of Eq. (6b) is a bias which is a priori unknown and which we will evaluate in this study.
2.1.4 Source inference methods
The inferred sources, , were derived from Eqs. (3) or (6) assuming the covariance term (last term on right-hand side) was null. The method used to infer the source was either a simple division (Eq. 3) or an optimisation of the linear system using the linear model function lm in R (package stats, R version 3.2.3), with either M = 1 (single source) or M = 9 (multiple sources):
The bias δSi was then evaluated as the difference between the inferred sources and the modelled sources averaged over each period:
As shown in Eqs. (3) and (6) the overall mean bias δSi contains (i) a bias term due to the inference method which is dependent mainly on the conditioning of the matrix Di,j (see Supplement Sect. S2) and (ii) a bias term which is intrinsically linked to the covariance between Di,j and Sj (Eqs. 3 and 6). Thus, with Eq. (8) we evaluate the sum of the two biases without distinction. In order to infer the sources, the elements of the dispersion matrix Di,j need to be determined. The next part details how these were estimated with a dispersion model.
2.2 The dispersion model used for determining the transfer matrix Di,j
The elements of the transfer matrix , which is by definition the concentration at location xi and time t generated by a source Sj of strength Sj(t)=1, were calculated using a dispersion model. The FIDES-3D model (Loubet et al., 2010), based on the analytical solution of the advection–diffusion equation of Philip (1959) was used for that purpose. This model was first compared with a backward Lagrangian stochastic dispersion model (bLS, the WindTrax software, Thunder Beach Scientific, Nanaimo, Canada; Flesch et al., 1995), and successively tuned to mimic the bLS. The two models and how the FIDES model was tuned are briefly described hereafter and detailed in Supplement Sects. S3 and S4.
The FIDES model is based on the Philip (1959) solution of the advection–diffusion equation, which assumes power law profiles for the wind speed U(z) and the vertical diffusivity Kz(z) at height z. This approach also assumes no chemical reactions in the atmosphere and spatial horizontal homogeneity of roughness length (z0), wind speed (U), and vertical and lateral diffusivity (Kz and Ky). The dispersion model is detailed in Huang (1979) and Loubet et al. (2010). The details of the model and the way the transfer function was estimated are detailed in Supplement Sect. S2.
The Schmidt number, which is the ratio of momentum to scalar vertical diffusivity , is key in dispersion modelling, as it determines the vertical diffusion rate of scalars. Wilson (2015) demonstrated that bLS and dispersion models like FIDES give different values of Sc by constitution. In order to assure consistency of the Philip (1959) approach with bLS models, considered as references in dispersion modelling, we chose to tune the Philip (1959) model to get the same Sc number as in WindTrax as described by Flesch et al. (1995). The details are given in Supplement Sect. S4. The comparison showed that the tuned FIDES model gives very similar concentrations to WindTrax at measurement heights lower than 2 m above the source, although slightly overestimated under stable and neutral conditions and slightly underestimated under unstable conditions. The correlation between the two models is however very high (.96), meaning that using the tuned FIDES model to characterise source inference performance will lead to results comparable to WindTrax. Moreover, since in this study the same model is used for predicting and for inferring the fluxes, the results are self-consistent.
2.3 Ammonia sources from simple SVAT modelling and prescribed emissions potentials
In order to evaluate the bias introduced by time averaging the concentrations when inferring single or multiple sources (third term in Eqs. 3 and 6), we generated NH3 emissions patterns mimicking the behaviour of real sources as closely as possible. With that goal, we used the SurfAtm-NH3 model developed by Personne et al. (2009) for two purposes: (i) evaluating the turbulence parameters (the friction velocity u∗ and the Monin Obukhov length L) from the meteorological datasets to parameterise the dispersion models, and (ii) providing the surface temperature T(z0) and the surface resistances in order to calculate ammonia emissions patterns.
The SurfAtm-NH3 model is a one-dimensional, bidirectional surface–vegetation–atmosphere transfer (SVAT) model, which simulates the latent (LE) and sensible (H) heat fluxes, as well as the NH3 fluxes between the biogenic surfaces and the atmosphere. It is a resistance analogue model separately treating the vegetation layer and the soil layer, and coupling a slightly modified (Choudhury and Monteith, 1988) model of energy balance and the two-layer bidirectional NH3 exchange model of Nemitz et al. (2000) with a water balance model. Unless otherwise stated, the surface was considered a bare soil with z0=5 mm, displacement height (d) = 0 m, and leaf area index (LAI) = 0.
The ammonia emissions patterns were modelled using the resistance approach and assuming atmospheric concentration was zero, which is a reasonable assumption following nitrogen application and leads to patterns mimicking reality, which is what we are seeking here:
Here Ra(zref) is the aerodynamic resistance at the reference height zref=3.17 m, and is the soil boundary layer resistance for ammonia as described in Personne et al. (2009). The ground surface compensation point concentration (Cpground) was expressed as a function of Γ, the ratio of to H+ concentrations in the soil liquid phase at the surface, as in Loubet et al. (2012):
where Kh and Kd are the Henry and the dissociation constant for NH3, respectively, and T(z0) is the soil surface temperature. Since we wanted to evaluate the correlation between the transfer function Di,j and the source strength Sj, which is the bias in the inference problem (Eq. 6), the NH3 volatilisation was modelled to reproduce the variety of existing kinetics of NH3 emissions from fields. With that goal, three Γ patterns were simulated:
a constant Γ=Γ0, which would mimic background NH3 emissions from soils;
an exponentially decreasing , which best represents NH3 emissions following slurry application;
a Gaussian , which would represent the typical NH3 emissions following urea application.
Here Γ0 is the maximum Γ during the period, t is the time in days, and τ0 is the duration of the emissions in days. The factor 4.6 was chosen so that when t=τ0, Γ goes down to 1 % of Γ0. The duration of the emissions was chosen to be 4 weeks, τ0=28 days. The timescale of the exponential decrease we used here was around 6 days, which is twice as large as the one reported by Massad et al. (2010) for slurry application (2.9 days). While these Γ patterns gave the weekly trend of NH3 emissions, the daily patterns were produced by the thermodynamical and turbulence drivers of NH3 emissions, which were explicitly taken into account through the compensation point (Eq. 10). To facilitate understanding, in most of the paper only the constant Γ was considered, and the effect of modifying the source strength was evaluated in a sensitivity study.
2.4 Spatial set-up of the sources and concentration sensors
The sources (plots) were considered to be squares with width xplot and aligned south–north. Two configurations were considered: (1) a single-source configuration and (2) a multiple-source configuration, which mimics typical agronomic trials with nine sources (plots) placed next to each other, with three treatments of three repetitions each. Each treatment was assigned a value of Γ0 different from the others, while the three repetitions of the same treatment were assigned the same value of Γ. The concentration sensor (receptors) locations, xi, were set in the middle of each plot, at several heights zi. (Fig. 1).
A number of plot sizes (xplot=25, 50, 100 and 200 m on the side), and receptor heights (zi=0.25, 0.5, 1 and 2 m), were tested successively. Several source strengths and dynamics were also tested: Γ was first considered constant with time (pattern 1) in all the plots, and the Γ0 values of each of the three treatments were either chosen to be significantly different in strength (104, 105, 106), or of the same order of magnitude (1000, 2000, 4000). Then the three Γ patterns (constant, exponential and Gaussian) were randomly assigned to the treatments for each simulation period. The ammonia background concentration, Cbgd, was considered constant and equal to 1 ppb except when studying the sensitivity of the inference method to the background concentration, where it was set as unknown. Throughout this study, an “optimum” block configuration was considered (shown in Fig. 1c), which avoided trivial configurations like aligned blocks and maximised the mean distance between blocks as in a Latin-square design.
2.5 Simulation details
2.5.1 Meteorological data and fertiliser application periods
A range of meteorological conditions were simulated based on the half-hourly meteorological data of the FR-Gri ICOS site in 2008. In total 13 periods of 28 days were considered, which spanned the whole year except the last 2 days of the year. Each period consisted of 1344 half-hourly data.
2.5.2 Concentration sensor integration periods
In order to evaluate the influence of the concentration averaging period on the source inference, several integration periods τ were tested: 0.5 (no integration), 3, 6, 12, 24, 48 and 168 h (7 days). In practice the concentrations were computed at each sensor location using Eq. (6) over 0.5 h: at that timescale, which corresponds to the spectral gap, the covariance term is assumed to be negligible (Van der Hoven, 1957). Then the averaged concentrations were computed for all integration periods.
2.5.3 Sensitivity to inferential method scenarios
Several scenarios were considered and summarised in Table 1.
The background concentration was either supposed known and fixed to the prescribed values (C1–C4) or was inferred (C5–C7).
The three repetitions of each treatment were either supposed to have the same source strength (C2, C4, C5, C6) or they were inferred independently (C1, C3, C7). In C2, C4, C5 and C6, Si=Sm for all i and m values belonging to the same treatment. In practice a new dispersion matrix was calculated by averaging together all columns belonging to the same treatment (matrix dimension N × 3). Three strength values of S were inferred to be tested.
Either one concentration sensor at each source location (zi) was considered (C1, C2, C5) or two sensors positioned at two heights were considered (C3, C4, C6, C7). All the measurement heights and their combinations were considered.
2.6 Statistical indicators
For each run the mean bias (BIAS) and the normalised mean bias (NBIAS) were calculated as , , where Nτ is the number of the time-averaged samples over each 28-day period and cumSi and are the inferred and observed cumulated fluxes over the same period. The medians and interquartiles of these statistical indicators were then calculated over the 13 periods of 28 days for 2008.
2.7 Real experimental test case
In order to evaluate the feasibility of the method we applied it to a real test case (Fig. 2). The trial was located at La Chapelle Saint-Sauveur in France ( N, W) and performed from 5 to 26 April 2011. Soil texture was loamy with a pH in water of 6.2 and a bulk density of 1.4 t m−3 in the first 15 cm. The experimental unit was composed by six squared subplots 20 m wide with two repetitions of three treatments: (1) surface application of cattle slurry, (2) surface application and incorporation of the same slurry, and (3) no application. Slurry pH was 7.5 with a dry matter (DM) content of 6.05 % and C : N ratio of 10.4 and it contained 38.4 g N kg−1 (DM) as total nitrogen and 13.2 g N-NH4 kg−1 (DM) as ammoniacal nitrogen. Slurry was applied on 5 April 2011 at a rate of 49 m3 ha−1, which led to 114 kg N ha−1 and 39 kg N-NH4 ha−1. The application was identical between the two repetitions with a small standard deviation (< 0.2 kg N ha−1). The incorporation was performed in two subplots 1 h after the end of the slurry spreading with a disc harrower at a depth of 0.10 m. The soil humidity between 0 and 5 cm depth was homogeneous over the blocks and decreased from 20 ± 1 to 17 ± 1 % w∕w between the start and the end of the experiment. Meteorological data were measured at less than 50 m from the central plots (Fig. 2). Air temperature, relative humidity, global solar radiation, wind velocity and direction were recorded every 30 min at 2 m height. The turbulence parameters (u∗ and L), input of the dispersion models, were evaluated with a simple energy balance model of Holtslag and Van Ulden (1983) assuming a Bowen ratio of 0.5 and a deep soil temperature equal to the averaged ambient temperature. Ammonia concentration was measured with diffusive samplers (ALPHA), (Sutton et al., 2001; Tang et al., 2001, 2009), which were placed at the centre of each subplot at two heights (0.32 and 0.87 m from the ground) as well as next to the assay at three locations (5 m away from the plots) at 3 m height. The ALPHA samplers were set in place just after slurry application and incorporation (between 14:20 and 14:50 LT) and left exposed subsequently for 3, 22, 23, 23, 71 h (3 days) and 359 h (15 days), hence spanning 21 days. The diffusive samplers were prepared prior to the experiment, stored at 4∘ C in a refrigerator and analysed by colorimetry. Since no background concentrations were measured at a reasonable distance from the field, the background concentration was assumed as the minimum over the whole period of the concentrations measured on the 3 m height masts.
3.1 Meteorological data range and simulated ammonia sources
The meteorological conditions over the 13 periods represented a good sample of temperate climate conditions. The friction velocity u∗ varied between 0.024 and 1.181 m s−1, and the stability parameter z∕L at 1 m height varied between −49 and 21 (Fig. 3). It is noticeable that u∗ showed greater variability during the winter than during the summer, while it was the opposite for z∕L. The surface temperature also showed a structure varying between periods, with a larger temperature range during the summer (from 5.7 to 50.4 ∘C) than during the winter (from −5.2 to 22.9 ∘C). This surface temperature variability is an essential feature to representing real-case ammonia sources (Sutton et al., 2009), which shows a variability reflecting both the surface temperature and the resistance variations (Eqs. 9 and 10).
3.2 Example ammonia concentration dynamics modelled with the tuned FIDES model
The modelled ammonia concentrations reproduced typical patterns measured above field following nitrogen application well, with maximum concentrations during the day and minimum concentrations at night (Fig. 4). These patterns are a consequence of daily variations of the sources driven by surface temperature combined with variations in the aerodynamic transfer function Di,j, which behaves similarly to a transfer resistance (see Supplement Sect. S1). The integration periods are also shown in Fig. 4, which illustrates the progressive loss of information of the pattern structure with integration periods. Particularly, it can be seen that the day-to-night variation is captured up to an integration period of 6 h. Moreover, it should be noted that averaging also means overestimating lower concentrations and underestimating higher concentrations.
3.3 Evaluation of the inference method for a single source and a single sensor
At first we evaluate the bias of the inference method for the simpler case of a single source and a single sensor placed in the centre of the source field at several heights, assuming we know the background concentration (strategy C1; Fig. 1a). This case has the advantage of having a condition number equal to 1 (Supplement Sect. S2 and Eq. S1) and a bias δS which is well defined and equal to (Eq. 8). This section hence focuses on evaluating the influence of sensor height, time integration and source dimension on the bias without dealing with the complexity of the interactions between multiple fields.
3.3.1 Example of inferred source dynamics
Figure 5 reports an example source inference, which shows the progressive smoothing of the source with integration period. We first see that the source strength corresponding to Γ=105 leads to ammonia emissions ranging from 0 to ∼ 1 in the winter, which corresponds to 0.71 kg N ha−1 day−1. Over the entire year, the maximum emissions occur during the hottest days and reach up to 7.1 kg N ha−1 day−1. Regarding the inference method, it can be seen in that example that, up to 24 h, the variability in emissions over the period is captured quite well.
3.3.2 Effect of target height, source dimension and integration period on the bias δS for a single source
In this simpler case shown in Fig. 6, the fractional bias of the inferred emissions is mostly negative for the combination in which the ratio of sensor height to plot dimension is small and integration times are larger than 6 h. According to Eq. (5), this means that the covariance term is negative for these conditions, meaning that any increase in source strength S at a time t is correlated with a decrease in the transfer function D(t)and vice versa. This is expected as S(t) increases with the surface temperature (Eq. 10) and is proportional to (Eq. 9), while D(t) is proportional to the aerodynamic resistance Ra(zref), as shown in Supplement Sect. S1. Hence, over daily periods, S and D are negatively correlated: S increases during the day and decreases at night (due to temperature and wind speed daily patterns), while D decreases during the day and increases at night (mainly due to wind speed patterns). This is expected to be a general feature for NH3 surface fluxes as the daily variability reproduced by the model used in this study is representative of most situations from mineral and organic fertilisation to urine patches or seabird colonies (Ferrara et al., 2014; Flechard et al., 2013; Milford et al., 2001; Móring et al., 2016; Personne et al., 2015; Riddick et al., 2014; Sutton et al., 2013).
The median bias δSi tends to increase in magnitude with the sensor height for large fields (xplot=100 and 200 m), while it decreases for smaller fields (xplot=25 and 50) when sensor height gets close to the field boundary layer height. Furthermore, δSi becomes positive and very large when sensors are above the field boundary layer height (Fig. 6). For large fields, the increase in the magnitude of the bias with lower sensor height is expected as D decreases with height in absolute value. For small fields, the decrease in the bias corresponds to a loss of information as D gets close to zero when the sensor gets closer to the field boundary layer height. For heights above this limit, we observe a change in sign of the bias, which can be explained by the fact that the sensor concentration footprint is not in the source during stable conditions (at night), while it is in the source under unstable conditions during the day. The inference method will hence not work if at least one sensor is not below the plot boundary layer height.
We also note that for integration periods equal to or below 3 h, the fractional bias is slightly positive, which can be explained by the positive correlation between S and D at small timescales. This is because of the influence of u∗ on T(z0): for a given solar radiation and air temperature over small timescales (< 3 h), an increase in u∗ leads to a decrease in T(z0), which leads to an exponential increase in the surface compensation point according to Eq. (10). However, at the same time, Ra(z)−1 decreases, but linearly with u∗. The resulting ammonia emissions calculated with Eq. (9) nevertheless increases because the exponential effect of temperature overcomes the linear effect of the exchange velocity (data not shown). This effect is more visible for large fields than small fields because over small fields an additional effect is that when u∗ decreases, the footprint increases and the source “seen” by the targets hence decreases because it incorporates a fraction of zero emissions sources.
Overall, the median fractional bias for weekly integrated emissions over a 25 m field and sensor heights below 0.5 m was overall −8 % with an interquartile (−14 to −2 %). We can conclude that the bias of the NH3 emissions is reproducible within ±6 %. We can also conclude that it would be better to place the concentration sensor at a low height to minimise the bias of the method.
3.3.3 Effect of surface boundary layer turbulence on the inference method for a single source
The inference method depends on the turbulence at the site and especially on the main drivers of the dispersion, which are the friction velocity and the stability regime. Indeed, Fig. 7 shows that the relative root-mean-square residual of the inferred source (RRMSR) decreases with increasing u∗ at long integration periods and is larger in slightly stable than near-neutral or slightly unstable conditions. Figure 7 also shows that under stable conditions or low u∗ the RRMSR increases by more than an order of magnitude (up to 50 %) when integration periods increase from 6 to 12 h, which catches most of the source variance. We also see that under near-neutral or high u∗ conditions, the third quartile of the RRMSR remains below 10 % for all integration periods. Finally, we also see that the larger third quartiles at short integration periods are obtained with intermediate u∗ values or slightly unstable conditions. A similar response of the bias to u∗ and 1∕L was reported by Fig. 6 in Flesch et al. (2004) and Fig. 3 in Gao et al. (2009) in controlled source experiments. While Gao et al. (2009) attributed the bias of the inference method to parameterisation of the stability dependence of the turbulent parameters (z∕L), in this study this cannot happen since we use the same parameterisation for prescribing the concentration and inferring it. In our case, the interpretation is to be linked with Eq. (5): the smaller u∗ or the most stable conditions also correspond to the larger time derivatives of source strength (driven by surface temperature and surface exchange resistances) as well as the larger time derivatives of transfer function D. We hence expect that under such conditions, the covariance between the transfer function and the source strength will be larger than under near-neutral conditions. In a more heuristic view, under low turbulence, large time derivatives of concentrations are expected above a source due to low mixing (small changes in mixing lead to large variations in concentrations).
We conclude that the inference method with a long integration period will lead to very moderate biases for locations with near-neutral conditions and high wind speed, but may lead to much larger bias under stable conditions and low wind speed as soon as the integration period reaches 12 h.
3.4 Multiple-source case
In contrast to the single-source case, with multiple sources (see Fig. 1b) the inference method leads to biases at small integration times as can be seen in the example reported in Fig. 8. In that specific case, the emissions of treatments 2 (Γ=105) and 3 (Γ=106) are 10 times and 100 times larger than those of treatment 1 (Γ=104), respectively. This leads to concentrations over plots of treatment 1 (and to a lesser extent over those of treatment 2) being highly correlated to emissions from plots of treatment 3 (and hence less with subplots of treatment 1). As a result, inferring emissions of plots of treatment 1 becomes harder as soon as averaging periods become larger or equal to 3 h. This can be viewed as a progressive loss of information of the treatment 1 contribution to concentrations due to the overweighing contribution of treatment 3 plots. However, we also see that treatments 2 and 3 seem quite correctly inferred for integration times smaller than 48 h.
In the following we will first evaluate the influence of the length of integration periods, sensor heights and plot dimensions on the fractional biases made when inferring the source. Each factor will be evaluated independently of the others in order to understand the processes behind it. For these evaluations background concentration was kept constant at 1 µg NH3 m−3. Strategy C1 was used except when testing sensor heights, for which strategy C3, which uses two targets, was also used. These two strategies assume that the background concentration is known, which avoids any compensating effects between source and background concentration inferences. Then the sensitivity of the methodology to the (i) emissions ratios between two of the three treatments and (ii) the variability in the background concentration were evaluated. Finally, seven inversion strategies were compared to determine which was the most robust (Table 1).
3.4.1 Effect of integration periods on the bias
We first consider strategy C1, which is the simplest configuration, in which plots are independent, background concentration is known and one target is used above each plot. Figure 9 shows that for the given treatment range (∼ 1–100 µg NH3 m−2 s−1), the fractional mean bias is lower than 0.2 in magnitude for the treatment emitting the most (treatment 3, Γ=106), lower than 0.4 for the intermediate treatment (treatment 2, Γ=105) and up to 8 for the treatment emitting the least (treatment 1, Γ=105); here we considered the 0.25–0.75 quantiles. The bias of the highest treatment (treatment 3) actually behaves similarly to a single-source case (Fig. 6), with a median bias around 10 % for 48 h integration periods. This is expected because treatment 1 and treatment 2 have a much smaller emissions strength and hence little influence on the concentration above the treatment 3 plots, which therefore behave in a similar manner to a single source. As a consequence, this bias in treatment 3 is mainly due to the anti-correlation between D and S, which increases with integration periods. The fractional mean bias is very large for treatment 1 even for small integration periods. The bias can either be positive or negative, showing that this method does not allow for a correct estimation of the smallest sources.
3.4.2 Effect of target heights on the bias
Figure 10 shows that the bias remains low as long as sensor heights are low enough to catch a sufficient part of the field footprint. When only a single height is used (strategy C1) this means that the sensor should be placed at 0.5 m or below for the field size we have tested here (25 m). The result is similar for a pair of sensors (strategy C3). For the lowest treatment though, the bias (and its variability) remain high whatever the height. It is interesting to notice that the heights which were found to provide an optimal inference of NH3 sources (below 0.5 m) are smaller than ZINST (the height at which the vertical flux can be approximated by the horizontal flux) reported by Wilson et al. (1982) (which were 0.9 m for 40 m diameter circular sources, and which we estimate as 0.65 m based on a power law extrapolation as in Laubach et al., 2012). It is also important to note that this height should vary with both the roughness length z0 and displacement height as was shown by Wilson et al. (1982) for ZINST.
3.4.3 Effect of plot size on the bias
Increasing the plot size from 25 to 200 m in width reduces the bias of the two highest source treatments for which the median bias reaches values around 10 %, while the interquartiles remain stable (Fig. 11). Conversely, in treatment 1 (Γ=104), the bias increases. It is expected that the bias in a multiple-source configuration never becomes smaller than the bias in a single-source problem, which is a limit linked to the time integration (covariance between the source and the concentration; see Eqs. 3 and 6). It is also expected that the biases remain higher than the single-source case until the source size increases sufficiently so that the concentration generated by a block on the neighbour fields becomes negligible compared to the concentration generated by the source below. This is what we observe in treatment 2 (Γ=105) and treatment 3 (Γ=06), with treatment 2 showing a median bias of −13 % (larger than in the single-source case) for the 200 m wide field, while the bias of the largest source tends to be −10 % [−17 %, −1 %], which is the range observed for a single source.
3.4.4 Sensitivity of the method to ratios of emissions potentials among treatments
A central question is the capability of the inference method to resolve small or large differences in emissions from the nearby blocks. Indeed, we can speculate that small differences will be hard to resolve while large differences will lead to large bias. In order to determine the resolution power of the method, we compared the performance of the inference method with a set of three treatments: the first treatment had Γ=0 to mimic a reference field receiving no nitrogen; the second treatment had a constant Γ=1000 corresponding to a small emissions (0.7 kg N ha−1 day−1), and the third treatment Γ was successively set to increasing values from 1500 to 105 (70 kg N ha−1 day−1). In this section we consider the background to be known (sensitivity to the background concentration will be evaluated in the next section).
Figure 12 shows the median and interquartile biases of the cumulated emissions for the longest integration period of 168 h over the ratio of the high-to-low source treatments. The bias of the largest source always remained around 14 %, which is larger than the single-source case. The bias of the lowest source increased with increasing inter-treatment source ratios from 13 to 40 %. In fact we find that the fractional bias increased approximately as a power function of the ratio of the two predicted sources (dotted lines, 0.11 x0.256).
3.4.5 Quality of background concentration estimations
As pointed out by Flesch et al. (2004), the knowledge of the background concentration is essential in a source inference problem. Retrieving the background necessitates having at least Nsources+1 sensors. Hence only strategies with two heights per plot or which assume identical emissions in treatment repetitions can be evaluated in their capacity of retrieving the background (strategy C2 to C7). In order to evaluate the sensitivity of the method when the background concentration varies with time, we set a realistic background concentration as a linear combination of u∗ and air temperature (Ta) with a mean of 6 µg NH3 m−3 and a standard deviation of 0.1 µg NH3 m−3. This test was performed with a range of treatments in order to elucidate the correlations between varying background and varying treatments. We see in Fig. 13 that the concentration, which follows a realistic pattern, is well retrieved even over the longest integration period of 168 h. However, we see that for the treatments with the largest source contrast (Γ=1000 and 105), the background concentration can be overestimated even for small integration periods (6 h). The median residual of the background concentration was smaller in magnitude than 0.05 µg NH3 m−3, except for the case with very large differences among treatments (0, 1000, 10 000), for which the residual reached 0.1 and 0.5 µg NH3 m−3 for the 6 h and 24 h or 168 h integration periods. Furthermore, the background concentrations were overestimated for the largest source ratios and underestimated for the lowest source ratios and longer integration periods (24 and 168 h).
3.4.6 Identifying the most robust strategy
Finally, to identify which strategy is the most suitable for retrieving the emissions from the multi-plot configuration, we compared all strategies for a simulation with a variable background (set as in the previous section) and two source ratios of 2 and 20 between treatments 2 and 3 (treatment 1 being a zero-source reference). We found, as expected, that strategies with known backgrounds have low biases compared to strategies that calculate the background, except for strategy C7, which provided biases similar to strategy C3, which is the strategy equivalent to C7 but with a known background (Fig. 14). We also see that incorporating some knowledge of the sources by assuming plots from the same treatment have the same emissions gave slightly better estimates when the background is known (strategies C2 and C4 compared to C3). This is however not true when the background is unknown, in which case the magnitude of the bias increases up to a median of 0.7 (strategies C5 and C6 compared to C7). It is due to compensation between background concentration and source strength as we have seen in Fig. 14 that the background concentration was overestimated in such cases. We also see, as expected, that the strategies with two sensors placed at different heights above each plot lead to better evaluations of the emissions. Overall, the strategy based on two sensors above each plot, which also assumes that sources are independent, seems to be the most robust (strategy C7). This strategy does not assume the background is known, nor does it assume the plots have similar emissions, which is more adapted to reality. Indeed, even though the same amount of nitrogen is applied in each repetition plot, the emissions may vary due to soil heterogeneity and advection. We finally obtain a median bias for strategy C7 which is −16 % with an interquartile [−8–22 %]. It is important to stress though that the minimums and maximums are further away, which indicates that under some rarer circumstances, the method may overestimate the sources by 12 % or underestimate them by 40 %. These cases correspond to integration periods with very low wind speeds and stable conditions.
3.5 Application of the methodology to a real test case with multiple treatments
The evaluation of the methodology on a real test case is shown in Figs. 15–17. The concentration measured above the surface-applied slurry (up to 200 µg N-NH3 m−3) is much higher than above the two other treatments (below 50 µg N-NH3 m−3) (Fig. 15).
The inference method gives very consistent results both in terms of comparison between repetitions (B1 and B2) of a given treatment and in terms of comparison between treatments (strategy C7 shown in Fig. 16). Surface slurry application showed the largest emissions: 9 ± 0.3 kg N ha−1 in B1 and 10 ± 0.2 kg N ha−1 in B2 (median and confidence interval). This corresponds to an emissions factor around 24 % of the N-NH4 applied and 8 % of the total N applied, which is in line with agronomic references (Sintermann et al., 2011a; Sommer et al., 2006). In contrast, the incorporated slurry showed much smaller emissions: 0.3 ± 0.2 kg N ha−1 in B1 and 0.6 ± 0.2 kg N ha−1 in B2. It is noticeable that no application showed slight deposition, especially in B2: −0.26 ± 0.2 kg N ha−1 in B1 and −1.7 ± 0.2 kg N ha−1 in B2.
Comparing the inference strategies is instructive (Fig. 17). We see that in methods which assume a known background (strategies C3 and C4), the inferred emissions are slightly higher than when background is assumed unknown. We should state that we set the background concentration to the minimum concentration measured on the 3 m height masts because these were located too close to the plots to be considered real background masts. This explains why strategies C3 and C4 lead to higher estimates compared to strategies C6 and C7, as the background may have been underestimated. We also find that all methods consistently infer a deposition flux to the blocks with no application, which is consistent with our knowledge of ammonia exchange between the atmosphere and the ground (Flechard et al., 2013). Indeed, the concentration in the atmosphere, which is enriched by the nearby sources is expected to be higher than near the ground due to a low soil pH (6.1), a low nitrogen content in the soil surface (6–9.5 g N kg−1 DM) and a 20 % humid soil surface, hence leading to a flux from the air to the ground.
From our theoretical study we know that strategy C7 should give a bias around −16 % ± ∼ 7 %. Therefore, we could expect that the real flux is the one measured with C7 times 1.15 (±0.08); hence it would be 10.9 ± 1.3 kg N ha−1. This corresponds to 28 ± 3 % of the N-NH4 applied and ∼ 9 ± 1 % of the total N applied. For the incorporated slurry, the emissions are around 20 times smaller than the emissions from the surface-applied slurry. Under these conditions, the bias on the emissions would be around −20 %, which means that the corrected emissions would range from 0.5 to 2.5 % of the N-NH4 applied and 0.2 and 0.8 % of the total N applied. We should bear in mind that the theoretical correction is based on the median of the simulations performed with the 2008 dataset in Grignon, which had similar meteorological conditions to this trial. It would be much more relevant though for future developments to evaluate the bias based on the same method as developed here but with emissions and meteorological conditions taken from the real case.
3.6 Comparison with previous work
Several studies have reported methodologies for evaluating multiple sources using dispersion models. These were mostly based on backward Lagrangian modelling (Crenna et al., 2008; Flesch et al., 2009; Gao et al., 2008). There were several inference methods reported: the methods based on the inversion of the dispersion matrix Di,j or singular value decomposition of least-square optimisation (Flesch et al., 2009), which optimise the conditioning of the dispersion matrix, and one based on Bayesian inference (Yee and Flesch, 2010). Yee and Flesch (2010) showed that the Bayesian approach would avoid unrealistic source estimates that could appear when the matrix conditioning was poor. Unrealistic source estimates were for instance reported by Flesch et al. (2009), with negative emissions sources.
Ro et al. (2011) evaluated the bLS technique to infer two controlled methane surface sources with laser measurements. They found 0.6 recovery ratios (ratio of inferred to known source) if the fields were not in the footprint of the sensor but with adapted filters; they found a high degree of recovery of 1.1 ± 0.2 and 0.8 ± 0.1 for the two sources. They found that in contradiction to Crenna et al. (2008) and Flesch et al. (2009), even with large conditioning numbers they had high recovery rates.
Misselbrook (2005) compared different methodologies and showed that under high concentrations diffusion samplers may lead to overestimation of up to 70 % of the concentration. They suggest potential issues related to the deformation of the Teflon membrane, which would modify the distance between coated filters and the membrane itself, which could cause sampler saturation. There is hence some concern about the quality of diffusion samplers to measure concentrations at heights close to large sources, which would necessitate field validations.
3.6.1 Sensor positioning and conditioning number
Crenna et al. (2008) have clearly shown that the optimal sensor positioning should be so that each sensor preferentially sees a single source, and reversely, each source should preferentially influence a single sensor. In this study the source–sensor geometry was especially designed in a way that minimises the condition number by placing the sensors in the middle of each plot. For the smallest source (xplot=25 m), the conditioning number ranged from 1.97 to 3.01 (median 2.42) for sensors located at 0.25 m, and increased to 2.6–6.9 (median 3.2) for sensors at 0.5 m, 4.7–150 (median 21) for sensors at 1.0 m, and 40–165 000 (median 640) for sensors at 2 m. This shows that including at least one sensor per block at heights lower than the field width divided by 20 would ensure that the conditioning number remains lower than in most trials reported by Crenna et al. (2008).
By comparing different strategies we have found that the strategies using two sensors over each source systematically led to improved performances (C3 versus C1 and C6 versus C5, Fig. 14). This is also in line with the results of Crenna et al. (2008), who showed that using more sensors separated spatially improves the performance of the inference method. Hence we can conclude that the inference method we used is based on a well-conditioned system which leads to robust results of the least-square optimisation. This is further illustrated by the real-case example (Figs. 15–17), which shows a good reproducibility among block repetitions. Indeed, good reproducibility among repetitions is a check for evaluating the quality of the inference method in real test cases. The use of the Bayesian inference method would however also be valuable in the set-up we propose here.
3.6.2 Effect of time-integrating sensors on the source inference quality
The use of time-averaging sensors for estimating ammonia sources was already reported by Sanz et al. (2010), Theobald et al. (2013), Carozzi et al. (2013a, b), Ferrara et al. (2014) and Riddick et al. (2016a, 2014). All these studies have shown the feasibility of these measurements; however only a few of them allow the estimation of the impact of averaging: Riddick et al. (2014) measured emissions from a bird colony on Ascension Island with WindTrax using both several ALPHA samplers in a transect across the colony and a continuous analyser for ammonia (AiRRmonia, Mechatronics, NL) downwind. They also averaged the continuous sampler concentrations to evaluate the effect of averaging on the emissions estimates. They found as we do here that averaging over monthly periods would lead to systematic underestimations from −9 to −66 %. They also found that estimations from diffusive samplers would lead to average underestimations of −12 %. This is very close to what we find here for a single source over 1 week (Fig. 6). In a similar comparison Riddick et al. (2016b) found that time integration led to slight overestimations with the integration approach, which is within the range of statistics of the bias we have found for the larger area sources (third quartile in Fig. 6).
3.6.3 Dependency on meteorological conditions
We should bear in mind that the use of time-averaging sensors in the inference method is also highly dependent on the surface layer turbulent structure as shown by Fig. 7. We find, as expected, that stable conditions or low wind speed conditions are those that lead to the highest potential bias (as shown by the third quartile under stable conditions at the bottom of Fig. 7). This is a well-known limitation of inverse dispersion modelling which was reported by Flesch et al. (2009, 2004) and which suggested that inverse dispersion would be inaccurate for .15 m s−1 and . However, both our study and the studies of Riddick et al. (2014, 2016b) show that this is not as much of an issue for ammonia emissions. Indeed, this is due to the fact that ammonia emissions follow a daily cycle with low emissions at night and high emissions during the day. This is firstly because the ground surface compensation point concentration (Cpground) has an exponential dependency on surface temperature as assumed in Eq. (10) based on known thermodynamical equilibrium constants (Flechard et al., 2013). This is secondly due to the fact that ammonia emission is a diffusion-based process which is limited by the surface resistances, as modelled in Eq. (9), which leads to small fluxes when Ra(zref) and become large, which happens during low wind speeds (they are both roughly inversely proportional to wind speed) and stable conditions, which also happens at night (Flechard et al., 2013). In real situations, the combination of small turbulence and high surface concentration leads to a further decrease in the flux, which is dependent on the difference between Cpground and the concentration in the atmosphere above (a feature which was not accounted for in this study as this would imply a higher degree of complexity in the modelling approach). This means that the results we found in this study would not apply for species with an emissions pattern with different temporal dynamics (either constant or anti-correlated with surface temperature or wind speed).
In this study we have demonstrated that it is possible to infer, with reasonable biases, ammonia emissions from multiple small fields located near each other using a combination of a dispersion model and a set of passive diffusion sensors which integrate over a few hours to weekly periods. We found that the Philip (1959) analytical model in FIDES gave similar concentrations as the backward Lagrangian stochastic model WindTrax at 2 m above a small source, under neutral and stable stratification as long as the stability correction functions used in both models are similar and the Schmidt number is identical (here set to 0.64). Under unstable conditions FIDES gave 20 % smaller concentrations at 2 m compared to WindTrax.
We demonstrated by theoretical considerations that passive sensors always lead to the underestimation of ammonia emissions for an isolated source because of the negative time correlation between the ammonia emissions and the transfer function. Using a yearly meteorological dataset typical of the oceanic climate of western Europe we found that the bias over weekly integration times is typically −8 ± 6 %, which is in line with previous reports. Larger biases are expected for meteorological conditions with stable conditions and low wind speeds as soon as the integration period is larger than 12 h.
We showed that the quality of the inference method for multiple sources was dependent on the number of sensors considered above each plot. The most essential technique to minimise the bias of the method was to place a sensor in the middle of each source within the boundary layer. The quality of the sensor positioning was evaluated using “condition numbers” which ranged from 2 to 3 for a sensor placed at 25 cm above the ground to much higher values (40–1.6 × 105) for a sensor at 2 m above 25 m width sources. Although the lowest sensors have the best condition number, we would rather recommend using heights of 50 cm above the canopy in order to reduce uncertainty in positioning the sensors close to the ground as well as avoid non-diffusive transfer conditions. Similarly, although the highest sensors had low condition numbers, they were shown to improve the robustness of the sources' inference, especially for evaluating the background concentrations. Using replicates of each treatment was found to be essential for evaluating the quality of the inference and derive robust statistical indicators for each treatment.
When considering a system, characteristic of agronomic trials, composed of a low and a high potential source and a reference with no nitrogen application, we found that the fractional bias remained smaller than around 25 % for ratios between the largest and smallest sources lower than a factor of 5 and increased as a power function of the ratio. Furthermore, the dynamics of the emissions were found not to strongly affect the fractional bias. As expected, we also found that the fractional bias decreased with increasing source dimensions, especially for the lowest source strength in a multiple-source trial.
Finally, a test on a practical trial proved the applicability of the method in real situations with contrasted emissions. We indeed calculated ammonia emissions of around 27 ± 3 % of the total ammoniacal nitrogen applied for surface-applied slurry while we found less than 1 % of emissions for the treatments with incorporated slurry.
This method could also be improved by incorporating knowledge of the surface source dynamics into the inference procedure. Further work is required, however, for validating the method, for instance using prescribed emissions, and to evaluate the method for growing crops using real measurements with diffusion samplers close to the ground.
The model is available as an R package upon request to the authors. The datasets used in this paper can be obtained from the authors upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/bg-15-3439-2018-supplement.
The authors declare that they have no conflict of interest.
This study was supported by EU FP7 NitroEurope-IP (grant number 017841) and
ECLAIRE (grant number 282910), French national projects CASDAR VOLAT'NH3
(grant number 0933), ADEME EVAPRO (grant number 1560C0036) and ADEME EVAMIN
(grant number 1660C0012). We thank Erwan Personne for the use of the
SurfAtm-NH3 model, and the technical team of the ARVALIS research
station of “La Jaillière” for its involvement in the conduct of the
“real test case” experiment.
Edited by: Paul Stoy
Reviewed by: two anonymous referees
Carozzi, M., Ferrara, R. M., Rana, G., and Acutis, M.: Evaluation of mitigation strategies to reduce ammonia losses from slurry fertilisation on arable lands, Sci. Total Environ., 449, 126–133, 2013a.
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, 306–318, 2013b.
Choudhury, B. J. and Monteith, J. L.: A four-layer model for the heat budget of homogeneous land surfaces, Q. J. Roy. Meteor. Soc., 114, 373–398, 1988.
CITEPA: Inventaire des émissions de polluants atmosphériques en France métropolitaine, format CEE-NU, CITEPA 494/Convention MATE 26/2001, Centre Interprofessionnel Technique d'Etudes de la Pollution Atmosphérique, 2017.
Council, E.: Directive 96/61/EC of 24th September 1996 concerning integrated pollution prevention and control, European Council, Brussels, Belgium, 1996.
Council, E.: 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, European Council, Brussels, Belgium, 2016.
Crenna, B. R., Flesch, T. K., and Wilson, J. D.: Influence of source-sensor geometry on multi-source emission rate estimates, Atmos. Environ., 42, 7373–7383, 2008.
ECETOC: Ammonia emissions to air in Western Europe, European Centre for Ecotoxicology and Toxicology of Chemicals, Avenue E Van Nieuwenhuyse 4, Brussels, 1994.
EUROSTAT: Agroenvironmental indicator – ammonia emission, Eurostat, Luxemburg, 2012.
Faburé, J., Rogier, S., Loubet, B., Génermont, S., Saint-Jean, S., Bedos, C., and Cellier, P.: Synthèse bibliographique sur la contribution de l'agriculture à l'émission de particules vers l'atmosphère: identification de facteurs d'émission, ADEME/INRA, 2011.
Famulari, D., Fowler, D., Nemitz, E., Hargreaves, K. J., Storeton-West, R. L., Rutherford, G., Tang, Y. S., Sutton, M. A., and Weston, K. J.: Development of a low-cost system for measuring conditional time-averaged gradients of SO2 and NH3, Environ. Monit. Assess., 161, 11–27, 2010.
Ferrara, R. M., Loubet, B., Di Tommasi, P., Bertolini, T., Magliulo, V., Cellier, P., Eugster, W., and Rana, G.: Eddy covariance measurement of ammonia fluxes: Comparison of high frequency correction methodologies, Agr. Forest Meteorol., 158, 30–42, 2012.
Ferrara, R. M., Loubet, B., Decuq, C., Palumbo, A. D., Di Tommasi, P., Magliulo, V., Masson, S., Personne, E., Cellier, P., and Rana, G.: Ammonia volatilisation following urea fertilisation in an sorghum crop in Italy irrigated, Agr. Forest Meteorol., 195, 179–191, 2014.
Ferrara, R. M., Carozzi, M., Di Tommasi, P., Nelson, D. D., Fratini, G., Bertolini, T., Magliulo, V., Acutis, M., and Rana, G.: Dynamics of ammonia volatilisation measured by eddy covariance during slurry spreading in north Italy, Agr. Ecosyst. Environ., 219, 1–13, 2016.
Flechard, C. R. and Fowler, D.: Atmospheric ammonia at a moorland site. II: Long-term surface-atmosphere micrometeorological flux measurements, Q. J. Roy. Meteor. Soc., 124, 759–791, 1998.
Flechard, C. R., Massad, R.-S., Loubet, B., Personne, E., Simpson, D., Bash, J. O., Cooter, E. J., Nemitz, E., and Sutton, M. A.: Advances in understanding, models and parameterizations of biosphere-atmosphere ammonia exchange, Biogeosciences, 10, 5183–5225, https://doi.org/10.5194/bg-10-5183-2013, 2013.
Flesch, T. K., Wilson, J. D., and Yee, E.: Backward-Time Lagrangian Stochastic Dispersion Models and Their Application to Estimate Gaseous Emissions, J. Appl. Meteorol., 34, 1320–1332, 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, 2004.
Flesch, T. K., Harper, L. A., Desjardins, R. L., Gao, Z. L., and Crenna, B. P.: Multi-Source Emission Determination Using an Inverse-Dispersion Technique, Bound.-Lay. Meteorol., 132, 11–30, 2009.
Gao, Z. L., Desjardins, R. L., van Haarlem, R. P., and Flesch, T. K.: Estimating Gas Emissions from Multiple Sources Using a Backward Lagrangian Stochastic Model, JAPCA J. Air Waste Ma., 58, 1415–1421, 2008.
Gao, Z. L., Mauder, M., Desjardins, R. L., Flesch, T. K., and van Haarlem, R. P.: Assessment of the backward Lagrangian Stochastic dispersion technique for continuous measurements of CH4 emissions, Agr. Forest Meteorol., 149, 1516–1523, 2009.
Gericke, D., Pacholski, A., and Kage, H.: Measurement of ammonia emissions in multi-plot field experiments, Biosyst. Eng., 108, 164–173, 2011.
Häni, C., Sintermann, J., Kupper, T., Jocher, M., and Neftel, A.: Ammonia emission after slurry application to grassland in Switzerland, Atmos. Environ., 125, 92–99, 2016.
Holtslag, A. A. M. and Vanulden, A. P.: A Simple Scheme for Daytime Estimates of the Surface Fluxes from Routine Weather Data, J. Clim. Appl. Meteorol., 22, 517–529, 1983.
Huang, C. H.: A theory of dispersion in turbulent shear flow, Atmos. Environ., 13, 453–463, 1979.
Kaimal, J. C. and Finnigan, J. J.: Atmospheric Boundary Layer Flows, Their structure and measurement, Oxford University Press, New York, 289 pp., 1994.
Kormann, R. and Meixner, F. X.: An analytical footprint model for non-neutral stratification, Bound.-Lay. Meteorol., 99, 207–224, 2001.
Laubach, J., Taghizadeh-Toosi, A., Sherlock, R. R., and Kelliher, F. M.: Measuring and modelling ammonia emissions from a regular pattern of cattle urine patches, Agr. Forest Meteorol., 156, 1–17, 2012.
Loubet, B., Milford, C., Sutton, M. A., and Cellier, P.: Investigation of the interaction between sources and sinks of atmospheric ammonia in an upland landscape using a simplified dispersion-exchange model, J. Geophys. Res.-Atmos., 106, 24183–24195, 2001.
Loubet, B., Genermont, S., Ferrara, R., Bedos, G., Decuq, G., Personne, E., Fanucci, O., Durand, B., Rana, G., and Cellier, P.: An inverse model to estimate ammonia emissions from fields, Eur. J. Soil Sci., 61, 793–805, 2010.
Loubet, B., Decuq, C., Personne, E., Massad, R. S., Flechard, C., Fanucci, O., Mascher, N., Gueudet, J.-C., Masson, S., Durand, B., Genermont, S., Fauvel, Y., and Cellier, P.: Investigating the stomatal, cuticular and soil ammonia fluxes over a growing tritical crop under high acidic loads, Biogeosciences, 9, 1537–1552, https://doi.org/10.5194/bg-9-1537-2012, 2012.
Lushi, E. and Stockie, J. M.: An inverse Gaussian plume approach for estimating atmospheric pollutant emissions from multiple point sources, Atmos. Environ., 44, 1097–1107, 2010.
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, https://doi.org/10.5194/acp-10-10359-2010, 2010.
McGinn, S. M. and Janzen, H. H.: Ammonia sources in agriculture and their measurement, Can. J. Soil Sci., 78, 139–148, 1998.
Milford, C., Hargreaves, K. J., Sutton, M. A., Loubet, B., and Cellier, P.: Fluxes of NH3 and CO2 over upland moorland in the vicinity of agricultural land, J. Geophys. Res.-Atmos., 106, 24169–24181, 2001.
Milford, C., Theobald, M. R., Nemitz, E., Hargreaves, K. J., Horvath, L., Raso, J., Dämmgen, U., Neftel, A., Jones, S. K., Hensen, A., Loubet, B., Cellier, P., and Sutton, M. A.: Ammonia fluxes in relation to cutting and fertilization of an intensively managed grassland derived from an inter-comparison of gradient measurements, Biogeosciences, 6, 819–834, https://doi.org/10.5194/bg-6-819-2009, 2009.
Misselbrook, T. H., Nicholson, F. A., Chambers, B. J., and Johnson, R. A.: Measuring ammonia emissions from land applied manure: an intercomparison of commonly used samplers and techniques, Environ. Pollut., 135, 389–397, 2005.
Móring, A., Vieno, M., Doherty, R. M., Laubach, J., Taghizadeh-Toosi, A., and Sutton, M. A.: A process-based model for ammonia emission from urine patches, GAG (Generation of Ammonia from Grazing): description and sensitivity analysis, Biogeosciences, 13, 1837–1861, https://doi.org/10.5194/bg-13-1837-2016, 2016.
Mukherjee, S., McMillan, A. M. S., Sturman, A. P., Harvey, M. J., and Laubach, J.: Footprint methods to separate N2O emission rates from adjacent paddock areas, Int. J. Biometeorol., 59, 325–338, 2015.
Nemitz, E., Sutton, M. A., Schjoerring, J. K., Husted, S., and Wyers, G. P.: Resistance modelling of ammonia exchange over oilseed rape, Agr. Forest Meteorol., 105, 405–425, 2000.
Nemitz, E., Dorsey, J. R., Flynn, M. J., Gallagher, M. W., Hensen, A., Erisman, J.-W., Owen, S. M., Dämmgen, U., and Sutton, M. A.: Aerosol fluxes and particle growth above managed grassland, Biogeosciences, 6, 1627–1645, https://doi.org/10.5194/bg-6-1627-2009, 2009.
Pacholski, A., Cai, G. X., Nieder, R., Richter, J., Fan, X. H., Zhu, Z. L., and Roelcke, M.: Calibration of a simple method for determining ammonia volatilization in the field – comparative measurements in Henan Province, China, Nutr. Cycl. Agroecosys., 74, 259–273, 2006.
Personne, E., Loubet, B., Herrmann, B., Mattsson, M., Schjoerring, J. K., Nemitz, E., Sutton, M. A., and Cellier, P.: SURFATM-NH3: a model combining the surface energy balance and bi-directional exchanges of ammonia applied at the field scale, Biogeosciences, 6, 1371–1388, https://doi.org/10.5194/bg-6-1371-2009, 2009.
Personne, E., Tardy, F., Genermont, S., Decuq, C., Gueudet, J. C., Mascher, N., Durand, B., Masson, S., Lauransot, M., Flechard, C., Burkhardt, J., and Loubet, B.: Investigating sources and sinks for ammonia exchanges between the atmosphere and a wheat canopy following slurry application with trailing hose, Agr. Forest Meteorol., 207, 11–23, 2015.
Philip, J. R.: The Theory of Local Advection 1, J. Meteorol., 16, 535–547, 1959.
Riddick, S. N., Blackall, T. D., Dragosits, U., Daunt, F., Braban, C. F., Tang, Y. S., MacFarlane, W., Taylor, S., Wanless, S., and Sutton, M. A.: Measurement of ammonia emissions from tropical seabird colonies, Atmos. Environ., 89, 35–42, 2014.
Riddick, S. N, Ward, D., Hess, P., Mahowald, N., Massad, R., and Holland, E.: Estimate of changes in agricultural terrestrial nitrogen pathways and ammonia emissions from 1850 to present in the Community Earth System Model, Biogeosciences, 13, 3397–3426, https://doi.org/10.5194/bg-13-3397-2016, 2016a.
Riddick, S. N., Blackall, T. D., Dragosits, U., Daunt, F., Newell, M., Braban, C. F., Tang, Y. S., Schmale, J., Hill, P. W., Wanless, S., Trathan, P., and Sutton, M. A.: Measurement of ammonia emissions from temperate and sub-polar seabird colonies, Atmos. Environ., 134, 40–50, 2016b.
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, 2011.
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, 2010.
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, https://doi.org/10.5194/amt-4-1821-2011, 2011a.
Sintermann, J., Spirig, C., Jordan, A., Kuhn, U., Ammann, C., and Neftel, A.: Eddy covariance flux measurements of ammonia by high temperature chemical ionisation mass spectrometry, Atmos. Meas. Tech., 4, 599–616, https://doi.org/10.5194/amt-4-599-2011, 2011b.
Sintermann, J., Neftel, A., Ammann, C., Häni, C., Hensen, A., Loubet, B., and Flechard, C. R.: Are ammonia emissions from field-applied slurry substantially over-estimated in European emission inventories?, Biogeosciences, 9, 1611–1632, https://doi.org/10.5194/bg-9-1611-2012, 2012.
Sommer, S. G., Genermont, S., Cellier, P., Hutchings, N. J., Olesen, J. E., and Morvan, T.: Processes controlling ammonia emission from livestock slurry in the field, Eur. J. Agron., 19, 465–486, 2003.
Sommer, S. G., McGinn, S. M., and Flesch, T. K.: Simple use of the backwards Lagrangian stochastic dispersion technique for measuring ammonia emission from small field-plots, Eur. J. Agron., 23, 1–7, 2005.
Sommer, S. G., Jensen, L. S., Clausen, S. B., and SØGaard, H. T.: Ammonia volatilization from surface-applied livestock slurry as affected by slurry composition and slurry infiltration depth, J. Agr. Sci., 144, 229–235, 2006.
Spirig, C., Flechard, C. R., Ammann, C., and Neftel, A.: The annual ammonia budget of fertilised cut grassland – Part 1: Micrometeorological flux measurements and emissions after slurry application, Biogeosciences, 7, 521–536, https://doi.org/10.5194/bg-7-521-2010, 2010.
Sun, K., Tao, L., Miller, D. J., Zondlo, M. A., Shonkwiler, K. B., Nash, C., and Ham, J. M.: Open-path eddy covariance measurements of ammonia fluxes from a beef cattle feedlot, Agr. Forest Meteorol., 213, 193–202, 2015.
Sutton, M. A., Miners, B., Tang, Y. S., Milford, C., Wyers, G. P., Duyzer, J. H., and Fowler, D.: Comparison of low cost measurement techniques for long-term monitoring of atmospheric ammonia, J. Environ. Monitor., 3, 446–453, 2001.
Sutton, M. A., Nemitz, E., Milford, C., Campbell, C., Erisman, J. W., Hensen, A., Cellier, P., David, M., Loubet, B., Personne, E., Schjoerring, J. K., Mattsson, M., Dorsey, J. R., Gallagher, M. W., Horvath, L., Weidinger, T., Meszaros, R., Dämmgen, U., Neftel, A., Herrmann, B., Lehman, B. E., Flechard, C., and Burkhardt, J.: Dynamics of ammonia exchange with cut grassland: synthesis of results and conclusions of the GRAMINAE Integrated Experiment, Biogeosciences, 6, 2907–2934, https://doi.org/10.5194/bg-6-2907-2009, 2009.
Sutton, M. A., Oenema, O., Erisman, J. W., Leip, A., van Grinsven, H., and Winiwarter, W.: Too much of a good thing, Nature, 472, 159–161, 2011.
Sutton, M. A., Reis, S., Riddick, S. N., Dragosits, U., Nemitz, E., Theobald, M. R., Tang, Y. S., Braban, C. F., Vieno, M., Dore, A. J., Mitchell, R. F., Wanless, S., Daunt, F., Fowler, D., Blackall, T. D., Milford, C., Flechard, C. R., Loubet, B., Massad, R., Cellier, P., Personne, E., Coheur, P. F., Clarisse, L., Van Damme, M., Ngadi, Y., Clerbaux, C., Skjoth, C. A., Geels, C., Hertel, O., Wichink Kruit, R. J., Pinder, R. W., Bash, J. O., Walker, J. T., Simpson, D., Horvath, L., Misselbrook, T. H., Bleeker, A., Dentener, F., and de Vries, W.: Towards a climate-dependent paradigm of ammonia emission and deposition, Philos. T. Roy. Soc. B, 368, 1–13, https://doi.org/10.1098/rstb.2013.0166, 2013.
Tang, Y. S., Cape, J. N., and Sutton, M. A.: Development and types of passive samplers for monitoring atmospheric NO2 and NH3 concentrations, Sci. World J., 1, 513–529, 2001.
Tang, Y. S., Simmons, I., van Dijk, N., Di Marco, C., Nemitz, E., Dammgen, U., Gilke, K., Djuricic, V., Vidic, S., Gliha, Z., Borovecki, D., Mitosinkova, M., Hanssen, J. E., Uggerud, T. H., Sanz, M. J., Sanz, P., Chorda, J. V., Flechard, C. R., Fauvel, Y., Ferm, M., Perrino, C., and Sutton, M. A.: European scale application of atmospheric reactive nitrogen measurements in a low-cost approach to infer dry deposition fluxes, Agr. Ecosyst. Environ., 133, 183–195, 2009.
Theobald, M. R., Crittenden, P. D., Tang, Y. S., and Sutton, M. A.: The application of inverse-dispersion and gradient methods to estimate ammonia emissions from a penguin colony, Atmos. Environ., 81, 320–329, 2013.
Thomson, L. C., Hirst, B., Gibson, G., Gillespie, S., Jonathan, P., Skeldon, K. D., and Padgett, M. J.: An improved algorithm for locating a gas source using inverse methods, Atmos. Environ., 41, 1128–1134, 2007.
UNECE (Ed.): Protocol to Abate Acidification, Eutrophication and Ground-level Ozone to the Convention on Long-range Transboundary Air Pollution, available at: http://www.unece.org/env/lrtap/multi h1.html (last access: March 2018), UNECE, Brussels, Belgium, 1999, as amended on 4 May 2012.
Van der Hoven, I.: Power Spectrum of Horizontal Wind Speed in the Frequency Range from 0.0007 to 900 Cycles Per Hour, J. Meteorol., 14, 160–164, 1957.
Vandré, R. and Kaupenjohann, M.: In Situ Measurement of Ammonia Emissions from Organic Fertilizers in Plot Experiments, Soil Sci. Soc. Am. J., 62, 467–473, 1998.
Wang, W., Liu, W. Q., Zhang, T. S., and Ren, M. Y.: Evaluation of backward Lagrangian stochastic (bLS) model to estimate gas emissions from complex sources based on numerical simulations, Atmos. Environ., 67, 211–218, 2013.
Whitehead, J. D., Twigg, M., Famulari, D., Nemitz, E., Sutton, M. A., Gallagher, M. W., and Fowler, D.: Evaluation of laser absorption spectroscopic techniques for eddy covariance flux measurements of ammonia, Environ. Sci. Technol., 42, 2041–2046, 2008.
Wilson, J. D.: Computing the Flux Footprint, Bound.-Lay. Meteorol., 156, 1–14, 2015.
Wilson, J. D. and Shum, W. K. N.: A Reexamination of the Integrated Horizontal Flux Method for Estimating Volatilization from Circular Plots, Agr. Forest Meteorol., 57, 281–295, 1992.
Wilson, J. D., Thurtell, G. W., Kidd, G. E., and Beauchamp, E. G.: Estimation of the rate of gaseous mass transfer from a surface source plot to the atmosphere, Atmos. Environ., 16, 1861–1867, 1982.
Yee, E.: Theory for reconstruction of an unknown number of contaminant sources using probabilistic inference, Bound.-Lay. Meteorol., 127, 359–394, 2008.
Yee, E. and Flesch, T. K.: Inference of emission rates from multiple sources using Bayesian probability theory, J. Environ. Monitor., 12, 622–634, 2010.