Do alternative inventories converge on the spatiotemporal representation of spring ammonia emissions in France?

. Agriculture is the main source of ammonia (NH 3 ) in France, an important gaseous precursor of atmospheric particulate matter (PM). National and global emission inventories are known to have difﬁculty representing the large spatial and temporal variability inherent to atmospheric NH 3 . In this study, we compare NH 3 emissions in France during spring 2011 from one reference inventory, the TNO inventory, and two alternative inventories that account in different manners for both the spatial and temporal variabilities of the emissions: (i) the NH 3 SAT satellite-derived inventory based on IASI NH 3 columns and (ii) the CADASTRE-CIT inventory that combines NH 3 emissions due to nitrogen fertilization calculated with the mechanistic model VOLT’AIR on the database of the CADASTRE_NH 3 framework and other source emissions from the CITEPA. The total spring budgets, from March to May 2011, at the national level are higher when calculated with both alternative inventories than with the reference, the difference being more marked with CADASTRE-CIT. NH 3 SAT and CADASTRE-CIT inventories both yield to large NH 3 spring emissions due to fertilization on soils with high pH in the northeastern part of France (65 and 135 ktNH 3 , respectively, vs. 48 ktNH 3 for TNO-GEN), while soil properties are not accounted for by the TNO-GEN methodology. For the other parts of France, the differences are smaller. The timing of fertilization and associated ammonia emissions is closely related to the nitrogen requirements and hence the phenological stage of the crops, and therefore to the crop year’s speciﬁc weather conditions. Maximum emissions are observed in March for 2011 for some regions for both alternative inventories, while April is the period with maximum emissions for the reference inventory regardless of the region or the year. Comparing the inventories at ﬁner temporal resolutions, typically at daily scale, large differences are found. The convergence of alternative, independent and complementary methods on the spatiotemporal

Abstract. Agriculture is the main source of ammonia (NH 3 ) in France, an important gaseous precursor of atmospheric particulate matter (PM). National and global emission inventories are known to have difficulty representing the large spatial and temporal variability inherent to atmospheric NH 3 . In this study, we compare NH 3 emissions in France during spring 2011 from one reference inventory, the TNO inventory, and two alternative inventories that account in different manners for both the spatial and temporal variabilities of the emissions: (i) the NH 3 SAT satellite-derived inventory based on IASI NH 3 columns and (ii) the CADASTRE-CIT inventory that combines NH 3 emissions due to nitrogen fertilization calculated with the mechanistic model VOLT'AIR on the database of the CADASTRE_NH 3 framework and other source emissions from the CITEPA. The total spring budgets, from March to May 2011, at the national level are higher when calculated with both alternative inventories than with the reference, the difference being more marked with CADASTRE-CIT. NH 3 SAT and CADASTRE-CIT inventories both yield to large NH 3 spring emissions due to fertilization on soils with high pH in the northeastern part of France (65 and 135 kt NH 3 , respectively, vs. 48 kt NH 3 for TNO-GEN), while soil properties are not accounted for by the TNO-GEN methodology. For the other parts of France, the differences are smaller. The timing of fertilization and associated ammonia emissions is closely related to the nitrogen requirements and hence the phenological stage of the crops, and therefore to the crop year's specific weather conditions. Maximum emissions are observed in March for 2011 for some regions for both alternative inventories, while April is the period with maximum emissions for the reference inventory regardless of the region or the year. Comparing the inventories at finer temporal resolutions, typically at daily scale, large differences are found. The convergence of alternative, independent and complementary methods on the spatiotemporal representation of the spring NH 3 emissions, particularly over areas where the contribution of mineral fertilizer spreading to the spring budget is strong, encourages further developments in both prospective complementary directions, as this will help improve national NH 3 emission inventories.
Published by Copernicus Publications on behalf of the European Geosciences Union.

Introduction
France is a major crop producer and a major exporter of agricultural and food products. In 2014, it produced 2 %, 4 %, 5 %, 8 %, 8 % and 14 % of the global production of maize, sunflower, wheat, barley, rapeseed and sugar beet, respectively (Food and Agriculture Organization of the United Nations FAO;Schauberger et al., 2018). Through this food cultivation and also due to animal husbandry, agriculture is the main source of ammonia (NH 3 ) in the country. As an important gaseous precursor of particulate pollution, harmful to human life (Lelieveld et al., 2015;WHO, 2016), ammonia plays an important role in the regulation of inorganic aerosol concentrations Schaap, 2004, Bauer et al., 2016) and contributes to N deposition and potential exceedance of critical loads of ecosystems (Erisman et al., 2007;EEA European Environment Agency, 2014). In order to limit air pollution, also responsible for acidification and eutrophication, the new European National Emission Ceilings Directive 2016/2284, replacing the Directive 2001/81/EC, has now set ambitious national reduction commitments for ammonia. Ammonia emissions indeed have to be reduced by 19 % in 2030, compared with the 2005 levels (OJEU, 2016).
At the European scale, total NH 3 emissions are provided by the European Monitoring and Evaluation Program (EMEP) (Vestreng, 2005) or by the TNO-MACCIII (Kuenen et al., 2014) inventories that rely on national annual declarations and estimates of emission factors. Emissions are accounted for without separating fertilization and livestock. These reference inventories are widely used by the scientific community to study the impact of pollutant emissions on the chemical composition of the troposphere and on air quality. Nevertheless, uncertainties on the quantification of the NH 3 emissions are usually estimated to be between 100 % and 300 % of the annual budgets in the reference inventories (EMEP/EEA, 2016; Kuenen et al., 2014). In addition, the temporal and spatial variability may not be well represented in the reference inventories, as the temporal profiles used do not account for meteorology, soil properties and other local conditions. Moreover, fertilizer spreading is of particular interest, as these are applied during small periods, especially during a few weeks at the end of winter and early spring. However, the exact timing of fertilizer spreading is difficult to predict, as it depends on agricultural practices and meteorological conditions, which is not taken into account in the temporal disaggregation of the reference emission inventories. Both the inaccurate temporal resolutions and the misrepresentation of the spreading emissions largely explain the difficulty encountered by models to represent seasonal or daily patterns of NH 3 concentrations (Menut et al., 2012), and consequently particulate matter levels (Fortems-Cheiney et al., 2016).
To reduce these uncertainties, a better quantification of agricultural ammonia emissions and their temporal and spatial evolution is necessary. In particular, one of the challenges is to capture the right timing of fertilizer spreading at the weekly or even at the daily scale in order to reflect the effect of environmental and agronomic conditions on ammonia emissions. To this end, mechanistic models taking into account meteorological conditions, soil properties and agricultural practices have been developed (e.g., for Denmark (Skjøth et al., 2004), for the UK (Hellsten et al., 2008), and for mineral fertilization in springtime in France, Hamaoui-Laguel et al., 2014). Limitations for such approaches come from the fact that detailed agricultural data needed as input for such models are not available for most of the European countries. Moreover, agricultural practices of a specific country cannot be extrapolated to another country (Skjøth et al., 2011).
As an alternative to direct emission modeling, attempts have been made to constrain ammonia emissions through inverse approaches, based on ammonium wet deposition data (Paulot et al., 2014) or on satellite-observed atmospheric ammonia distributions, e.g., from the Tropospheric Emission Spectrometer (TES) (Zhu et al., 2013;Zhang et al., 2018), from the Infrared Atmospheric Sounding Interferometer (IASI) (Fortems-Cheiney et al., 2016;Van Damme et al., 2018;Adams et al., 2019) or from the Cross-track Infrared Sounder (CrIS) (Adams et al., 2019;Dammers et al., 2019). In principle, such emission estimates can be available shortly after observation. The advantage of satellite-derived estimates is also that these can be derived globally, at a high temporal scale (e.g., daily scale under clear sky). The downside of these, however, is that they do not provide information on the underlying sources of the emissions (fertilizers vs. husbandry), or for example the date of fertilization, the type of fertilizers used, the fertilization rates, etc., that could be important for the regulation of NH 3 emissions.
In this context, we compare ammonia emissions in France from inventories using the different approaches mentioned above: (i) the reference, in the sense that this is the method currently used, hereafter called TNO-GEN, is the European inventory based on the annual budgets provided by the TNO-MACCIII inventory (Kuenen et al., 2014) and on seasonal profiles from GENEMIS (Ebel et al., 1997); (ii) a first alternative inventory, hereafter called NH 3 SAT, is based on a top-down approach starting from the IASI-derived NH 3 columns; (iii) the other alternative inventory, hereafter called CADASTRE-CIT, is based on a bottom-up approach quantifying NH 3 emissions due to nitrogen fertilization combining spatiotemporal data and calculations performed within the CADASTRE_NH 3 framework with the mechanistic model VOLT'AIR Génermont et al., 2018) completed with livestock and other source emissions from the French Interprofessional Technical Centre for Studies on Air Pollution (CITEPA). This study aims at assessing the potential contribution of better spatial and temporal representation of fertilization-related ammonia emissions to the quality of ammonia emission inventories. The improvement is assessed in terms of total budget, spatial distribu-tion and timing of the emissions. The study period, spring 2011 (from 1 March to 31 May 2011), was chosen following three criteria. Firstly, because at the time of the study, the last French agricultural data were available from AGRESTE (AGRESTE, 2014) for the agricultural year 2010-2011, allowing the application of the CADASTRE_NH 3 framework for the quantification of the spatiotemporal distribution of NH 3 emissions due to nitrogen fertilization for this crop year Génermont et al., 2018). Secondly, ammonia emissions are enhanced during spring in accordance with crop requirements (Skjøth et al., 2004;Ramanantenasoa et al., 2018;Génermont et al., 2018). Finally, unlike autumn and winter months, the NH 3 spring levels are detectable with a better confidence in the IASI satellite observations , allowing the extension of the preliminary work of Fortems-Cheiney et al. (2016) to deduce NH 3 emissions from the IASI satellite instrument.
The three inventories and methods to build them used for this study are presented in Sect. 2, and the results of the comparison are given and discussed in Sect. 3.

Inventories
The three inventories TNO-GEN, NH 3 SAT and CADASTRE-CIT compared in this study are described in Table 1 and in the following sections. It is worth noting that only the CADASTRE-CIT inventory provides information on the respective contribution of fertilization and livestock emissions. The spatial and temporal resolutions of the inventories are also shown in Table 1.
The inter-comparison is made at the 0.5 • (longitude) × 0.25 • (latitude) resolution. The outputs of the TNO-GEN and the CADASTRE-CIT inventories could have consequently been aggregated or disaggregated.

TNO-GEN
In this study, TNO-GEN combines the annual budgets provided by the TNO-MACCIII inventory and the seasonal profiles to deduce the monthly variability of NH 3 emissions. This inventory is based on official annual emission data submitted by countries to EMEP/CEIP (European Monitoring and Evaluation Programme/Centre on Emission Inventories and Projections) for air pollutants. It is the update of the TNO-MACCII inventory (Kuenen et al., 2014). It is an inventory at 0.125 • × 0.0625 • resolution providing annual emissions of NH 3 from the agricultural sector, without separating the contributions from fertilization and livestock. Hereafter, we use the TNO-MACIII emissions of the year 2011. The seasonal profile of these emissions is prescribed according to the typical national factors provided by GENEMIS. This seasonal temporal profile used for the temporalization of emissions -the same one applied to the entire country -leads to a maximum in NH 3 emissions systematically in April over France (Ebel et al., 1997). The emissions remain constant between days in each month and between hours in each day.

NH 3 SAT
As a first alternative, a mass-balance approach, which is a common method for the quantification of the surface fluxes of short-lived species (Palmer et al., 2003;Jaeglé et al., 2004;Boersma et al., 2008;Lin et al., 2010) was set up. We used it to deduce NH 3 emissions from differences between NH 3 total columns observed by the IASI instrument and simulated by the CHIMERE regional chemical transport model (CTM) using the TNO-GEN inventory as inputs data.
2.2.1 The regional CTM CHIMERE CHIMERE simulates concentrations of gaseous and particulate chemical species (Menut et al., 2013;Mailler et al., 2017). For this study, we used the CHIMERE version 2013a. The horizontal resolution is given as follows: The vertical grid contains 17 layers from the surface to 200 hPa. This model is driven by the European Centre for Medium-Range Weather Forecasts global meteorological fields (Owens and Hewson, 2018). Climatological values from the LMDz-INCA global model (Szopa et al., 2008) are used to prescribe concentrations at the lateral and top boundaries and the initial atmospheric composition in the domain. For inorganic species, aerosol thermodynamic equilibrium is achieved using the ISORROPIA model (Nenes et al., 1998).
As in most of the models in the world, the parameterization of NH 3 dry deposition is unidirectional in CHIMERE. The parameterization of a bidirectional exchange with surfaces in Wichink Kruit et al. (2012) increased their yearly mean modeled LOTOS-EUROS European ammonia concentrations almost everywhere, and particularly over agricultural source areas. However, Zhu et al. (2015), with the Goddard Earth Observing System-Chemistry (GEOS-Chem) global CTM, estimated a decrease in NH 3 European concentrations in April, when the inclusion of a compensation point for vegetation is included. Further work needs to be done to better investigate the sensitivity of NH 3 concentrations to the deposition. Nevertheless, without such parameterization for bi-directional exchange, Azouz et al. (2019) assessed that regional models such as CHIMERE usually operating with large grid cell sizes simulate quite well the average NH 3 dry deposition flux over a large domain of simulation.
The evaluation of CHIMERE NH 3 and NH + 4 concentrations should be done against NH 3 (as done in Fortems-Cheiney et al., 2016) and NH + 4 measurements. Nevertheless, to our knowledge, there is no available NH 3 measurement over France for the focused period here. There is interpretable NH + 4 surface measurements at only one site, making the interpretation of the results difficult. NH 3 and NH +  (Petetin et al., 2016). From June to September 2010, 83 % of modeled total NH x was gaseous, while in the model, it was only 50 %; coherent with this, NH 3 was underestimated especially during warmer days. Thus, it is concluded for one site and season that particulate NH + 4 has a low to medium impact on NH 3 .

The IASI observations
We use data from the IASI-A instrument, flying on a low Sun-synchronous polar orbit aboard the Metop satellite since October 2006, with Equator crossing times of 09:30 (descending mode) and 21:30 (ascending mode) local sidereal time (LST) . The spatial resolution of its observations is about 12 km × 12 km at nadir. The algorithm used to retrieve NH 3 columns from the radiance spectra is described in Van Damme et al. (2017). Several improvements have been introduced since the description of Van Damme et al. (2014) and the version v1 used in our previous study (Fortems-Cheiney et al., 2016). In this study, we use the reanalyzed dataset ANNI-NH3-v2.2R, relying on ERA-Interim meteorological input data from the European Centre for Medium-Range Weather Forecasts (ECMWF) rather than the operationally provided Eumetsat IASI Level 2 (L2) data used for the standard nearreal-time version (Van Damme et al., 2017). We only consider land measurements from the morning overpass, as IASI is more sensitive at this time to the boundary layer, owing to more favorable thermal conditions (Clarisse et al., 2010, Van Damme et al., 2014. The IASI total columns are averaged into "superobservations" (average of all IASI data within the 0.5 • × 0.25 • resolution and for the given CHIMERE physical time step of about 5-10 min). As suggested by Van Damme et al. (2017), we no longer use weighted averages for this purpose. We performed a sensitivity test by selecting the IASI pixels for which the retrieval error does not exceed 100 %: the results of the temporal and spatial variability of the NH 3 French emissions presented in Sect. 3 did not significantly change, showing the robustness of the IASI NH 3 product (not shown).
The resulting monthly means of IASI NH 3 columns from March to May 2011 are shown in Fig. 1a-c. The spatiotemporal variability -with the highest values over northeastern France in March, and over northwestern France in Aprilis confirmed by the IASI 10-year and by the CrIS 5-year monthly means shown in Viatte et al. (2020). Note that the potential of IASI to provide information at high temporal resolution, up to daily scale, can be hampered by the cloud coverage as only observations with a cloud coverage lower than 10 % are delivered (Van Damme et al., 2017). To evaluate the impact of this limitation, the number of IASI superobservations used to calculate these monthly means, which represents the number of days over a month covered by IASI, is shown in Fig. 1d-f. On average, more than half of the month is sampled by IASI during spring, except in May in the northwestern part of France. The regions showing large IASI NH 3 values are consequently well sampled.

Deducing NH 3 SAT emissions
Relative differences between simulated columns by the CHIMERE regional CTM (described in Sect. 2.2.1, using the TNO-GEN emissions for the year 2011, described in Sect. 2.1) and observed IASI total columns (described in Sect. 2.2.2) are applied as a corrective factor to the reference emissions at daily and at grid-cell resolutions over France, from February to May 2011. As IASI "super-observations" provide one piece of information per day, the diurnal time profile of reference emissions cannot be improved: we apply the same daily correction factor to all hourly NH 3 emissions. When IASI is not available (i.e., observations with a cloud coverage higher than 10 %), the correction is not applied and the emissions remain equal to the TNO-GEN ones. To compare the emissions with the CADASTRE-CIT inventory, and their respective simulations with CHIMERE, the correction is only applied over France here.
With the mass-balance approach, the transport to neighboring cells is assumed to be negligible following Palmer et al. (2003). This approach has been debated by Turner et al. (2012), who found that non-local sources contribute substantially to columns of short-lived species. Li et al. (2019) evaluated the ability of both a mass-balance approach and a variational assimilation to recover known NH 3 emissions at different spatial resolutions. At a 2 • × 2.5 • resolution, they found that both methods yielded similar values. At a 0.25 • × 0.3125 • , the mass-balance approach led to values about 20 % higher compared to the variational ones. With our 0.5 • × 0.25 • resolution, the use of a mass-balance approach would lead to additional errors of less than about 20 % for the quantification of NH 3 emissions. This uncertainty is acceptable and much lower than the uncertainty existing in the annual and national budgets provided by emission inventories (EMEP/EEA, 2016). In this context, we choose to perform such a mass-balance approach to deduce NH 3 emissions from IASI ANNI-NH 3 -v2.2R super-observations. Additional uncertainty comes from the IASI observations. The IASI minimum detection limit is about 2-3 ppbv (∼ 4-6.10 15 molec. cm −2 ) (Clarisse et al., 2010). The signal-tonoise ratio therefore presents better performance for regions with high local concentrations (e.g., northern part of France, Fig. 1) than over low local concentration areas (e.g., southern parts of France, Fig. 1). There is no evaluation available yet for the IASI ANNI-NH 3 -v2.2R product used here.

CADASTRE-CIT
As a second alternative, a bottom-up approach was set up based on the finest national inventories available for anthropogenic sources of ammonia. The CADASTRE_NH 3 framework provides such an inventory for organic and mineral fertilization practices. This is, however, not the case for the other anthropogenic sources. For livestock emissions, with the exception of the stage of effluent spreading in the field, the less detailed inventory of the CITEPA is used. To meet the objectives of better specialization and temporalization, specific procedures are applied. These inventories are completed by the TNO-GEN inventory for the emissions of the other sectors. CADASTRE_NH 3 was implemented in order to represent in a realistic way spatiotemporal variability of French NH 3 emissions due to mineral and organic N fertilization and is fully described in Ramanantenasoa et al. (2018) and in Génermont et al. (2018). It has been constructed through the combined use of two types of resources: the process-based VOLT'AIR model and geo-referenced and temporally explicit databases for soil properties, meteorological conditions and N fertilization. VOLT'AIR is a 1D process-based model predicting NH 3 emissions from N fertilizers on bare soils, from physical, chemical and biological processes (Le Cadre, 2004;Garcia et al., 2012). It incorporates current knowledge on NH 3 volatilization after application of the main types of organic manure and mineral N fertilizers in the field. It takes into account the major factors known to influence NH 3 volatilization in the field, i.e., soil properties, weather conditions, cultural practices, and properties of mineral fertilizers and organic products. It runs at an hourly time step at the field scale for a period of several weeks, thus covering the entire volatilization duration of fertilization events.
Local features are attributed to each simulation unit, the small agricultural regions ( ; and nitrogen fertilization management practices, derived from data of the national AGRESTE survey of cultural practices for arable crops and grassland (Department of Statistics and Forecasting of the French Ministry of Agriculture) (AGRESTE, 2014). All input data required by VOLT'AIR are geographically overlaid and intersected with a geographical information system to generate input combinations in each SAR. Each input combination is used as the input data for a virtual 300 m × 300 m field for a simulation using VOLT'AIR. Exact times and dates of fertilizations are required to run VOLT'AIR, but for the sake of robustness, the statistical analysis of the survey data has been performed on the basis of 2-week intervals for the date of fertilization. Fertilizations are thus randomly distributed within these 2week intervals in proportion to their respective representation following Ramanantenasoa et al. (2018). Each simulation of NH 3 emissions is run at an hourly time step for a period of 2 months, starting 1 month before the fertilization in order to calculate soil water content at the time of application, and ending one month after fertilization, in order to cover the whole volatilization event.
About 160 000 runs with the VOLT'AIR model have been performed over the crop year 2010-2011 to produce ammo-nia emissions per hour, per hectare, per crop type, per SAR. Emissions can be aggregated at different spatial and temporal scales. At the spatial scale, they are weighted with the contribution of (i) each N fertilization management applied to each crop in each SAR and (ii) the area of the crop cultivated in the SAR. A procedure allows ammonia emissions to be produced at the required grid scale for the inventory comparison: it is based on cultivated areas for each crop as the key of desegregation-reaggregation from the SAR to the 0.015625 • × 0.03125 • grid. At the temporal scale, emissions are aggregated over daily, weekly or monthly bases for the sake of comparison with TNO-GEN and NH 3 SAT inventories. Volatilization taking place over several days, from a few days to several weeks, one fertilization in one field contributes to ammonia emissions over several days or weeks. Weather condition effects on overall ammonia emissions are thus the result of both their effects on fertilization timing and their effects on volatilization intensity and dynamics over 30 d from fertilizer application.
As they are not available in the agricultural practice survey, N fertilizations of vegetables, fruits and vines are not accounted for in CADASTRE_NH 3 : their contribution is minor for France overall, only accounting for 5 % of the total agricultural area (AGRESTE, 2010), but is important in particular regions. As the agricultural practice survey does not provide information over Corsica, this inventory is completed by the TNO-GEN inventory over this region.

Livestock emissions
As for the TNO-GEN inventory, French NH 3 emissions from livestock for the CADASTRE-CIT inventory are generated by using annual national emissions provided by CITEPA for 2011 (Fig. 3a). Nevertheless, here these emissions have been spatially distributed differently than for the TNO-GEN inventory. This has been done by using the FAO Gridded Livestock of the World database with a resolution of 30 arcsec. The temporalization of the emissions has been performed as a function of temperature and wind speed with the parameterizations of Skøjth et al. (2011) for the different subsectors.

Results and discussion
First, we analyze the different contributions of livestock and of fertilization to the spring budget in the CADASTRE-CIT inventory. Then, the comparison of the two alternative inventories NH 3 SAT and CADASTRE-CIT versus the reference inventory TNO-GEN and their inter-comparison are made at different temporal and spatial resolutions. We evaluate the inventories at the national scale and at the scale of the different French administrative regions (the administrative division in France on level 2 of the unified NUTS territory classification, NUTS2, shown in Fig. 2). We also analyze their spatial variability at the 0.5 • (longitude) × 0.25 • (latitude) resolution in order to draw a first picture of the consistency of the inventories in terms of the spring NH 3 total budget and to identify regions of interest. Finally, we focus on the temporal variability of the identified regions and discuss the agricultural practices that can influence the variability but also down to which temporal resolution the comparison of the inventories is relevant.

Respective contributions of different sources
The different contributions of livestock and of fertilization to the annual and to the spring budget in the CADASTRE-CIT inventory are shown in Fig. 3. This figure shows that the contribution of fertilization to the annual French budget can be strong, with emissions occurring mainly during spring.
The different contributions of livestock and of fertilization to the spring budget in the CADASTRE-CIT inventory highlight four different domains of interest. We can see that the contribution of fertilization on the high emissions of the northeastern part of France (box A in Fig. 2) is strong. For example, the contribution of fertilization is about 99 % in the region Champagne-Ardennes, and about 85 % in the region Picardie (Table 2). These emissions are mainly due to the use of mineral fertilizer on barley, sugar beet and potato Génermont et al., 2018). In particular, the use of urea or nitrogen solution and the high soil pH (Hamaoui-Laguel et al., 2014;Ramanantenasoa et al., 2018;Génermont et al., 2018) -parameters not taken into account by the TNO-GEN inventory -seem to be the factors responsible for the high emissions in this domain.
The second domain of interest is the northwestern part of France (box B in Fig. 2). Over this domain with high emissions, the NH 3 emissions are due in roughly equal parts to livestock (including animal housing, manure storage and grazing) and to fertilization, with a high use of organic manure . Livestock farming indeed produces large amounts of livestock manure available for application on grassland and on arable crops.
The third domain of interest is the southeastern part of France (box C in Fig. 2), showing the smallest spring NH 3 emissions. Finally, the contribution of fertilization on the emissions of the southwestern part of France (box D in Fig. 2) is strong.

French spring NH 3 total budget and its main spatial features
The spring NH 3 total budget is shown in  higher than the TNO-GEN reference one (234 kt NH 3 ).The CADASTRE-CIT inventory estimates higher NH 3 spring emissions, by about 30 %, than NH 3 SAT. The relative agreement on national budget between TNO-GEN and NH 3 SAT must be nuanced as total budget values from NH 3 SAT and TNO-GEN are close, but large differences in the spatial distribution of the French NH 3 emissions between TNO-GEN and both NH 3 SAT and CADASTRE-CIT can be observed (Fig. 4).
Indeed, over France and for the spring budget, the spatial correlation compared to the CADASTRE-CIT inventory, which should better represent the agricultural practices and their spatial distribution, is improved when using the NH 3 SAT inventory instead of using TNO-GEN (Pearson correlation coefficient r = 0.78 with NH 3 SAT against r = 0.72 with TNO-GEN).
The northwestern part of France presents the largest NH 3 emissions according to the TNO-GEN inventory (Fig. 4a). The TNO-GEN, NH 3 SAT and CADASTRE-CIT inventories lead to similar spring budget (68, 73 and 71 kt NH 3 , respectively), over this domain.
Over the southeastern part of France, CADASTRE-CIT is about 23 % lower than NH 3 SAT (28 and 37 kt NH 3 , respectively, Table 2). One hypothesis to explain the lower NH 3 emissions in CADASTRE-CIT is that market gardening is important in this area and not taken into account in the CADASTRE-CIT inventory Génermont et al., 2018). Nevertheless, market gardening is not included, to our knowledge, in the TNO-GEN inventory. TNO-GEN and NH 3 SAT inventories being in quite good agreement in terms of budget (35 and 37 kt NH 3 , respectively, Table 2. French national and regional budgets of NH 3 spring emissions, from March to May 2011, in kilotons of NH 3 (kt NH 3 ). The relative differences compared to the TNO-GEN are presented between brackets, in percentage. Regions for which the inventories NH 3 SAT and CADASTRE-CIT propose the same sign of relative differences are marked in bold. The contributions of the fertilization emissions to the NH 3 regional spring budget in the CADASTRE-CIT inventory are shown in percentage.  Table 2), further work is required to understand these discrepancies. Finally, over the southwestern part of France (box D in Fig. 2), IASI observations only trigger slight corrections to the TNO-GEN inventory over this area (29 and 30 kt NH 3 , respectively), and CADASTRE-CIT is 36 % higher than TNO-GEN and NH 3 SAT (Table 2).

Temporal variability of the NH 3 emissions at the subseasonal scale
Monthly regional budgets have been calculated for the three inventories. Figure 5 presents the monthly variability of the NH 3 emissions from February to May 2011 for the four domains of interest presented above. February is only displayed here as a baseline to show the sharp peak of NH 3 emissions in March over some domains. The contribution of the emissions due to livestock in the CADASTRE-CIT spring budget is also given. The TNO-GEN inventory shows rather similar NH 3 emissions from March to May for all regions, with a slight maximum in April (Fig. 5a), imposed by the GENEMIS monthly profiles used for the temporalization of emissions (Ebel et al., 1997). On the contrary, over the northeastern part of France, both NH 3 SAT and CADASTRE-CIT inventories show a maximum in March, and a decrease until May by about a factor of 1.5 to 2 (Fig. 5a). This maximum in March is also noticed by Tournadre et al. (2020), providing 9 years of total column observations from ground-based infrared remote sensing over the Paris megacity. We calculated the monthly contribution of livestock to the NH 3 emissions based on CADASTRE-CIT, which allows one to separate this contribution from the fertilization one. As in Fig. 3, Fig. 5a confirms that NH 3 emissions are mainly due to fertilization in the northeastern part of France and shows that the seasonal variation is mainly driven by this contribution. CADASTRE-CIT shows larger values than NH 3 SAT, which might be partly due to a possible low bias in the IASI observations (Dammers et al., 2017;Tournadre et al., 2020), combined with the fact that the TNO-GEN inventory (negatively biased compared to CADASTRE-CIT) is used as a prior for the mass balance approach with no correction applied to this prior when IASI observations are not available.
To go further in the comparison, we analyzed the daily variability of the NH 3 SAT and CADASTRE-CIT inventories, with TNO-GEN representing no daily variability (Fig. 6). To interpret the results, some limitations of the inventories have to be considered. For NH 3 SAT, the corrections applied to the TNO-GEN emissions are only applied for clear-sky conditions when IASI observations are available. In the CADASTRE-CIT, fertilization days are randomly selected within two-week intervals of application extracted from the farm survey analysis ; thus the actual day of fertilization is unknown. However, the NH 3 volatilization is continuous over several days after spreading, reducing the uncertainty introduced by this random selection. Moreover, the random selection is made at the field scale (see Sect. 2.3), and then spatially averaging at the CHIMERE resolution should also smooth the random selection effect. Both effects of weather conditions on fertilization timing, on the one hand, and on volatilization intensity  Fig. 2. The contribution of the emissions due to livestock in the CADASTRE-CIT monthly budgets is also given. Units are kilotons of NH 3 per month (kt NH 3 month −1 ). The monthly regional coverage of the IASI super-observations is given in percentage (in grey). and dynamics at the time of application and after application, on the other hand, are realistically produced. Hence, the NH 3 volatilization is continuous over several days after spreading, introducing additional smoothing.
The CADASTRE-CIT inventory presents a high day-today variability from March to May 2011 ( Fig. 6a) with several strong maxima of emissions -characteristic of emissions due to fertilizer application and a significant effect of the varying meteorological conditions. Day-to-day variability is also large for NH 3 SAT. However, NH 3 SAT and CADASTRE-CIT maxima are not very well correlated. Over the first 15 d in March 2011, the high emissions occurring from 9 to 16 March in CADASTRE-CIT are not reproduced by the NH 3 SAT inventory, potentially because of a lack of IASI coverage for this period (only about 40 % of the domain - Fig. 6a). NH 3 SAT shows an emission maximum 1 week earlier from 1 to 7 March. This time gap in emission  Fig. 2. Note that the scale is different between (a, b) and (c, d). Units are kilotons of NH 3 per day (kt NH 3 d −1 ). The contribution of the emissions due to livestock in the CADASTRE-CIT daily variability is also given. The daily regional coverage of the IASI super-observations is given in percentage (in grey). maxima could be explained by the random selection of fertilization days in CADASTRE-CIT. Over the last 15 d in March 2011 and over the first 15 d in April 2011, the maxima of emissions estimated by NH 3 SAT and CADASTRE-CIT are more correlated (e.g., from 7 to 11 April in CADASTRE-CIT versus from 9 to 11 March in NH 3 SAT). Finally, CADASTRE-CIT still shows high emissions in the last 15 d of April and in May 2011, particularly related to the use of fertilizer over corn. Despite a good coverage with IASI observations, no specific high emissions are derived from IASI for the same periods.
Over the northwestern part of France, the CADASTRE-CIT inventory is in agreement with TNO-GEN, with a slight maximum in April (Fig. 5b). The NH 3 SAT inventory shows a maximum in March, related to the high emissions seen during the first week of March 2011 (Fig. 6b). However, the monthly differences are much smaller than for the north-eastern part of France. Also for this domain, the day-to-day variabilities provided by CADASTRE-CIT and NH 3 SAT are mostly uncorrelated: a strong sharp maximum of emissions during the first week of March 2011 is seen in NH 3 SATbut not reproduced in CADASTRE-CIT (Fig. 6b). The highest daily emissions in CADASTRE-CIT occur during the two last weeks of April 2011 and are not reproduced by NH 3 SAT.
Over the southeastern and southwestern parts of France, month-to-month variations and emission amounts are much smaller than for the two previous domains (Figs. 5c, 6c and 5d, 6d, respectively). In the southeastern part of France, NH 3 SAT emissions are slightly larger than TNO-GEN and CADASTRE-CIT at the beginning of the period (February and March, Fig. 5c). The fertilization contribution to CADASTRE-CIT emissions decreases slightly at the end of the period. In the southwestern part of France, NH 3 SAT and TNO-GEN are very similar (Fig. 5d), and CADASTRE-CIT is slightly larger at the end of spring. This increase is mainly related to fertilization emissions, the livestock contribution being stable. Again, daily time series between both NH 3 SAT and CADASTRE-CIT inventories are uncorrelated ( Fig. 6c and d).

Conclusions
In this study, we performed an inter-comparison of two alternative inventories with the TNO-GEN reference inventory that quantify the French NH 3 emissions during spring 2011. One of the main conclusions of this study is that over regions with large mineral fertilizer use, like over the northeastern part of France, induced NH 3 emissions are probably considerably underestimated by the TNO-GEN reference inventory, as both the NH 3 SAT (constrained by IASI observations) and CADASTRE-CIT (process level oriented) show much larger values. For instance, over northeastern France, NH 3 SAT and CADASTRE-CIT show, respectively, a factor 1.4 and 2.8 larger spring 2011 emissions than TNO-GEN. Over the whole of France, NH 3 emissions are still more than 50 % larger in CADASTRE-CIT than in TNO-GEN. Average French NH 3 SAT emissions are about 10 % larger than TNO-GEN ones. Over the southern part of France, with more diverse agriculture as compared to the crop-intensive one in northeastern France, differences between the inventories are on the whole lower, and signs between CADASTRE-CIT/TNO-GEN and NH 3 SAT/TNO-GEN corrections are often opposite for different regions.
Month-to-month variations are again much more pronounced over northeastern France and show a maximum in March for both CADASTRE-CIT and NH 3 SAT. Day-today variations are large in CADASTRE-CIT and NH 3 SAT, roughly a factor of 5 between minimal and maximal values. This shows the interest in evaluating NH 3 emissions at a daily scale because this input is required for chemistry transport modeling of particulate matter formation and thus air quality. However, time series delivered by CADASTRE-CIT and NH 3 SAT are uncorrelated for all regions considered. This result can be partly explained by the lack of data in IASI NH 3 column observations under partially cloudy conditions, and by the fact that available information on agricultural practices is resolved at a 2-week scale.
Yet, current results of our study have important implications for air quality modeling over Europe. The important changes in the spatial distribution of NH 3 emissions as a function of soil properties are of general concern not only for France, but for all of Europe. Soils are alkaline or neutral (pH > 6) not only over northeastern France, but also over large parts of Italy, eastern Spain or eastern Germany (Reuter et al., 2008). Over these regions, our study suggests potentially larger NH 3 emissions than with a constant emission factor treatment, with impacts on fine-particle formation. These features should be included in "operational" emission inventories used for air quality modeling.
Thus, as a general conclusion, use of observational constrained or process-oriented emission inventories is clearly of added value for estimating better monthly averages over French areas with intensive mineral fertilizer use, but capacity for delivering day-to-day variations is not yet proven. This warrants further studies, including both those that refine hypotheses on days chosen by farmers for fertilizer spread out as a function of meteorological conditions and those that acquire and use continuous surface NH 3 measurements for validating satellite-or process-derived NH 3 emission variability.
Author contributions. All authors have contributed to writing the paper (main authors: AFC, GD, SG and MB). AFC has performed the mass-balance approach to deduce NH 3 emissions from NH 3 total columns observed by the IASI satellite instrument. KD, JMG and SG have performed the bottom-up approach providing the CADASTRE_NH 3 inventory for organic and mineral fertilization practices. FC has compiled this CADASTRE_NH 3 inventory with livestock emissions from the CITEPA. LC, PFC, MVD and CC are the principal investigators of the IASI NH 3 product. All authors discussed the results and contributed to the final paper.
Belgian F.R.S-FNRS. The authors are indebted to all those who provided input data for the CADASTRE_NH3 2010-2011 crop year: the Service de la Statistique et de la Prospective (SSP), Department of Statistics and Forecasting of the French Ministry of Agriculture, for the French cultural practices survey, supported by a public grant overseen by the French National Research Agency (ANR) as part of the "Investissements d'avenir" program (reference: ANR-10-EQPX-17 -Centre d'Accès Sécurisé aux Données -CASD); Météo-France for the data on weather conditions; the Agence de Service et de Paiement (ASP) and the Observatoire du Développement Rural (ODR) service unit for data on the Land Parcel Identification System (LPIS) for France. The authors are also grateful to Maharavo Marie Julie Ramanantenasoa for her technical assistance, funded by the Cadastre_NH3 project under agreement number ADEME 1081C0031. IASI is a joint mission of EUMETSAT and the Centre National d'Etudes Spatiales (CNES, France). The authors acknowledge the AERIS data infrastructure for providing access to the IASI data in this study. This research has been supported by the Belgian State Federal Office for Scientific, Technical and Cultural Affairs (Prodex arrangement IASI.FLOW). This work was granted access to the HPC resources of TGCC under the allocation A0070107232 made by GENCI.
Financial support. For this study, Audrey Fortems-Cheiney and Karine Dufossé were funded by the Amp'Air and the PolQA Primequal projects, respectively, under agreement numbers ADEME 1660C0013 and 1662C0023.
Review statement. This paper was edited by Drew Gentner and reviewed by three anonymous referees.