Tracking down global NH3 point sources with wind-adjusted superresolution

. As a precursor of atmospheric aerosols, ammonia (NH 3 ) is one of the primary gaseous air pollutants. Given its short atmospheric lifetime, ambient NH 3 concentrations are dominated by local sources. In a recent study, Van Damme et al. (2018) have highlighted the importance of NH 3 point sources, especially those associated with feedlots and industrial ammonia production. Their emissions were shown to be largely underestimated in bottom-up emission inventories. The discovery was made possible thanks to the use of oversampling techniques applied to 9 years of global daily IASI NH 3 satellite measurements. Oversampling allows one to increase the spatial resolution of averaged satellite data beyond what the satellites natively offer. Here we apply for the ﬁrst time superresolution techniques, which are commonplace in many ﬁelds that rely on imaging, to measurements of an atmospheric sounder, whose images consist of just single pixels. We demonstrate the principle on synthetic data and on IASI measurements of a surface parameter. Superresolution is a priori less suitable to be applied on measurements of variable atmospheric constituents, in particular those affected by transport. However, by ﬁrst applying the wind-rotation technique, in the superresolution becomes in NH 3 at a very high spatial resolution. We that than using this wind-adjusted superresolution technique, we introduce a new type of NH 3 map that allows tracking down point sources easily than the regular oversampled average. On a subset of known emitters, the source could be located within a median distance of 1.5 km. We subsequently a new global

Abstract. As a precursor of atmospheric aerosols, ammonia (NH 3 ) is one of the primary gaseous air pollutants. Given its short atmospheric lifetime, ambient NH 3 concentrations are dominated by local sources. In a recent study, Van Damme et al. (2018) have highlighted the importance of NH 3 point sources, especially those associated with feedlots and industrial ammonia production. Their emissions were shown to be largely underestimated in bottom-up emission inventories. The discovery was made possible thanks to the use of oversampling techniques applied to 9 years of global daily IASI NH 3 satellite measurements. Oversampling allows one to increase the spatial resolution of averaged satellite data beyond what the satellites natively offer. Here we apply for the first time superresolution techniques, which are commonplace in many fields that rely on imaging, to measurements of an atmospheric sounder, whose images consist of just single pixels. We demonstrate the principle on synthetic data and on IASI measurements of a surface parameter. Superresolution is a priori less suitable to be applied on measurements of variable atmospheric constituents, in particular those affected by transport. However, by first applying the wind-rotation technique, which was introduced in the study of other primary pollutants, superresolution becomes highly effective in mapping NH 3 at a very high spatial resolution. We show that plume transport can be revealed in greater detail than what was previously thought to be possible. Next, using this windadjusted superresolution technique, we introduce a new type of NH 3 map that allows tracking down point sources more easily than the regular oversampled average. On a subset of known emitters, the source could be located within a median distance of 1.5 km. We subsequently present a new global point-source catalog consisting of more than 500 localized and categorized point sources. Compared to our previous catalog, the number of identified sources more than doubled. In addition, we refined the classification of industries into five categories -fertilizer, coking, soda ash, geothermal and explosives industries -and introduced a new urban category for isolated NH 3 hotspots over cities. The latter mainly consists of African megacities, as clear isolation of such urban hotspots is almost never possible elsewhere due to the presence of a diffuse background with higher concentrations. The techniques presented in this paper can most likely be exploited in the study of point sources of other short-lived atmospheric pollutants such as SO 2 and NO 2 .

Introduction
As one of the primary forms of reactive nitrogen, NH 3 is essential in many of the Earth's biogeochemical processes. It is naturally present along with the nitrogen oxides in the global nitrogen cycle (Canfield et al., 2010;Fowler et al., 2013). However, the discovery of ammonia synthesis through the Haber-Bosch process in the early 1900s has made this vital compound available in almost unlimited quantities, supporting the explosive population growth in the last century (Erisman et al., 2008). As a result, the nitrogen cycle is currently perturbed beyond the safe operating space for humanity, which has led to a host of environmental and societal problems (Steffen et al., 2015). The most obvious direct impact of excess NH 3 is that on air quality, as atmospheric NH 3 is one of the main precursors of secondary particulate mat-Published by Copernicus Publications on behalf of the European Geosciences Union.
ter, which has important adverse health impacts (Lelieveld et al., 2015;Bauer et al., 2016). Emissions of the two other important precursors (SO 2 and NO x ), thanks to effective legislation, drastically decreased in the past 20 years in Europe and North America and have started to level off in eastern Asia (Aas et al., 2019;Georgoulias et al., 2019). In contrast, no such decreases are observed or expected in the near future for NH 3 (e.g., Warner et al., 2017;Sutton et al., 2013). Unlike the other precursors, NH 3 emissions are not well regulated, and in fact, the focus on decreasing NO x and SO 2 has already led to increased NH 3 emissions (Chang et al., 2016) and concentrations (Lachatre et al., 2018;Liu et al., 2018).
The lack of a global regulative framework stems in part from the historical relative difficulty in measuring NH 3 concentrations. Satellite-based measurements of NH 3 , which were discovered about a decade ago, offer an attractive complementary means of monitoring NH 3 . Satellite datasets have now reached sufficient maturity to be directly exploitable, even when the individual measurements come with large and variable uncertainties. Using satellite observations we have recently shown the importance of ammonia point sources on regional scales . In total, over 240 of the world's strongest point sources were identified, categorized and quantified. Somewhat expectedly, many of these point sources (or clusters thereof) were found to be associated with "concentrated animal feeding operations" (CAFOs; Zhu et al., 2015;Yuan et al., 2017). However, much more surprisingly was the number of identified industrial emitters and in particular those associated with ammonia and urea-based fertilizer production. An evaluation of the EDGAR inventory additionally showed that emission inventories vastly underestimate the majority of all point-source emissions, even when a conservative average NH 3 lifetime of 12 h is assumed in the calculation of the satellite-derived fluxes. Industrial processes could therefore be extremely important, especially on a regional scale. Altogether, these findings were made possible due to the availability of the large multiyear NH 3 dataset (Whitburn et al., 2016;Van Damme et al., 2017) derived from measurements of the IASI spaceborne instrument (Clerbaux et al., 2009) and the oversampling technique that was applied to sufficiently resolve localized emitters.
Oversampling techniques applied on measurements of satellite sounders allow obtaining average distributions of atmospheric constituents at a higher spatial resolution than the original measurements (Sun et al., 2018). They exploit the fact that the footprint on the ground of satellite measurements varies in location, size and shape each time the satellite samples an area. When pixels partially overlap, some information becomes available on their (subpixel) intersection. High-resolution mapping can, however, only be obtained by combining typically many hundreds of measurements. A crucial condition on which oversampling relies is that the pixel center and ground instantaneous field of view (GIFOV) of satellite measurements is known with a high accuracy (typically < 1 km, as opposed to the coarse spatial resolution of the extent of the satellite pixel, which is typically > 10 km). Practical implementation of oversampling is relatively straightforward once the footprint is known: a fine subgrid is constructed in which the value of each cell of the grid is obtained as the average value of all overlapping GIFOVs. Optionally, the averaging can be weighted to take into account measurement error, total pixel surface area and spatial response function. We refer to Sun et al. (2018) and Van Damme et al. (2018) for comprehensive reference material on averaging and oversampling, detailed algorithmic descriptions, and practical considerations for their implementation.
Oversampling has gradually become commonplace in the field of atmospheric remote sensing, especially in the study of short-lived pollutants such as NO 2 (Wenig et al., 2008;Russell et al., 2010), SO 2 (Fioletov et al., 2011(Fioletov et al., , 2013, HCHO (Zhu et al., 2014) and NH 3 (Van Damme et al., 2014. The increased spatial resolution enables in the first instance a much better identification of emission (point) sources, quantification of their emissions (Streets et al., 2013), and study of transport and plume chemistry (de Foy et al., 2009). Oversampling applied to the study of point sources becomes even more useful when wind fields are taken into account. Beirle et al. (2011) showed that binned averaging per wind direction allows simultaneous estimates of both emission strengths and atmospheric residence times. Valin et al. (2013) and Pommier et al. (2013) introduced the wind-rotation technique whereby each observation is rotated around the presumed point source according to the horizontal wind direction, effectively yielding a distribution in which the winds blow in the same direction. As we will also demonstrate (see Sect. 3), this reduces the overall spread of the transported pollutants and reduces contributions of nearby sources. Combining plume rotation with oversampling has proven to be a very successful technique for the study of NO 2 and SO 2 point sources, leading to massively improved inventories and emission estimates and better constraints on the atmospheric lifetime of these pollutants (Fioletov et al., 2015(Fioletov et al., , 2017Wang et al., 2015;de Foy et al., 2015;Lu et al., 2015;Liu et al., 2016;McLinden et al., 2016).
However, as pointed out in Sun et al. (2018), while oversampling offers an increased resolution, it still yields a smoothed representation of the true distributions. There exists a large field of research, collectively referred to as superresolution (Milanfar, 2010), that attempts to construct highresolution images from several possibly moving or distorted low-resolution representations of the same reality. Oversampling is in essence the simplest way of performing superresolution but in a way that does not fully exploit the spatial information of the measurements. Superresolution has been applied before in the field of remote sensing of land or land cover (e.g., Boucher et al., 2008;Xu et al., 2017), but even though it is theoretically possible, it has not been applied to atmospheric-sounding measurements. In this case, the "images", as taken by sounders, are of the lowest resolution; i.e., Atmos. Meas. Tech., 12, 5457-5473, 2019 www.atmos-meas-tech.net/12/5457/2019/ they correspond to single, uniformly colored pixels. Perhaps the main reason why superresolution has not been attempted before on atmospheric sounders is that these rely on the fact that the low-resolution samples should be derived from an underlying distribution that does not change in time (de Foy et al., 2009). When this is not the case, the smoothing introduced by oversampling is actually desirable. With the arrival of the wind-rotation technique, most of the variability observed for point-source emitters can be corrected for, and therefore superresolution becomes viable for short-lived species as NH 3 . In Sect. 2 we introduce superresolution and demonstrate its effectiveness on measurements of the IASI sounder for a parameter related to (constant) surface emissivity. Next we illustrate the application of what we coin "wind-adjusted supersampling" on an industrial point source of NH 3 . In Sect. 3 we use ideas from McLinden et al. (2016) to provide a new type of NH 3 map, one that is supersampled and wind-corrected at the same time. This map enables the identification of many new point sources in addition to the ones reported in Van Damme et al. (2018). We performed a detailed global analysis of this new map, which led to the identification and categorization of more than 500 point sources and which we present in Sect. 4.

Superresolved oversampling
The general superresolution problem does not have a unique solution, as the available low-resolution measurements typically do not hold all the required information content (i.e., the problem is underdetermined). It is usually also overdetermined because of measurement noise and, for our case especially, because of temporal variability. As a consequence, there is no unique best algorithm, and a myriad of alternatives exist. For this study, we chose the iterative back-projection (IBP) algorithm (Irani and Peleg, 1993), as it takes a particular intuitive and simple form for single-pixel satellite observations and allows addressing the ill-determined nature of the problem. It proceeds as follows. Suppose we have a set of single-pixel measurements M 0 of a spatially variable quantity. For the first iteration, the solution of the algorithm corresponds to the regular oversampling, which we will write as SS 1 = OS 1 = OS(M 0 ) (SS i stands for the solution of the supersampling obtained in iteration i, and OS stands for the oversampling operator). From this oversampled average, we then calculate simulated observations for each of the original individual observations, corresponding to what the instrument would see if the ground truth were SS 1 . The entire set of these simulated measurements will be denoted by M 1 = M(SS 1 ) (with M the operator that simulates the measurements). If the oversampled average OS 1 = SS 1 corresponds to the ground truth, then M 1 would clearly coincide with M 0 . However, as oversampling typically smooths out the observations, this is generally not the case. An im-proved estimate of the average (SS 2 ) can be obtained by adding OS(M 0 − M 1 ) to the oversampled average, therefore correcting (partially) the observed differences. This process then is repeated to obtain increasingly better estimates. The entire algorithm thus reads . . .
The solution converges to an average that is maximally consistent with the observations, i.e., M 0 ≈ M k for sufficiently large k (as shown in Elad and Feuer, 1997, IBP converges to the maximum likelihood estimate whereby M 0 − M k is minimized). Figure 1 illustrates the algorithm on synthetic data with an idealized ground truth made up of nine point sources (Fig. 1a), with a Gaussian spread between 0.5 and 40 pixels. The measurement footprint was assumed to be variable between 7 and 13 pixels. The SS 1 (Fig. 1b), SS 3 (Fig. 1c) and SS 50 (Fig. 1d) averages illustrate the convergence and strengths of the algorithm well, which reproduces most of the point sources near perfectly and even partly resolves the smallest feature (compare also with Sun et al., 2018, Fig. 8). Some small ringing effects are noticeable though after 50 iterations (best visible on a screen), which are the result of the undetermined nature of the problem (Dai et al., 2007). An example applied to real data is shown in Fig. 2, which shows part of the Sahara and Mediterranean Sea. The quantity on which the oversampling is applied is the brightness temperature difference (BTD) between the IASI channels at 1157 and 1168 cm −1 . This BTD, located in the atmospheric window, is sensitive to the sharp change in surface emissivity due to the presence of quartz (see Takashima and Masuda, 1987, who also illustrate that the relevant feature is not seen in airborne dust). Being related to the surface, it can be assumed to be reasonably constant for each overpass of IASI (note that it is not entirely the case: sand dunes do undergo changes over time, and surface emissivities can depend on the viewing angle and can be affected by changes in moist content). Comparison with visible imagery (Fig. 2a) shows, as expected, that the largest BTD values (> 4 K) are associated with the most sandy areas. The other desert areas exhibit widely varying values, and oceans are slightly negative. The oversampled average (Fig. 2b) captures most large features down to about 5 km in size. Recalling that the footprint of IASI is a 12 km diameter circle at nadir, and elongates to an ellipse of up to 20 km×39 km at off-nadir angles, this example illustrates why oversampling is such a powerful tech- Figure 1. Illustration of the supersampling technique on synthetic data. Panel (a) depicts an imaginary ground truth, made up of twodimensional Gaussian distributions, each with a different spread (0.5 to 40 pixels). The rectangles in (a) indicate the assumed footprint size of the measurements, varying between 7 and 13 pixels. Panel (b) shows the results of the common oversampling approach applied to 100 000 measurements scattered randomly over the area. Panels (c, d) provide the results of the supersampling technique after respectively 3 and 50 iterations. nique. However, the additional resolution brought by the supersampling is clear, even after three iterations. The smallest features that can be distinguished are about 3-4 km (after three iterations; Fig. 2c) and 2-3 km (after 30 iterations; Fig. 2d) in diameter. That said, with increasing iterations, artifacts start to appear due to enhancements of noise and the specific sampling of IASI (in particular, stripes parallel to the orbit track become apparent). Such overfitting to the data and a sensitivity to outliers is often seen in maximum likelihood optimizations (Milanfar, 2010). It can therefore be advanta-geous to stop the algorithm after a few iterations (which can also be required for computational reasons).

Wind-adjusted supersampling
In this section we illustrate the previously introduced supersampling on a wind-rotated NH 3 average centered around a point source. The ammonia plant at Horlivka (Gorlovka), Ukraine, was chosen as a test case. This plant made the news in 2013 because of the major NH 3 leak that occurred on Atmos. Meas. Tech., 12, 5457-5473, 2019 www.atmos-meas-tech.net/12/5457/2019/ For each subpanel, the number in the bottom left corner is the total average NH 3 column amount over the entire area.
6 August, killing five people and injuring many more. We refer to the Wikipedia article for a detailed description of the event and a list of related newspaper articles (Wikipedia, 2019). The accident itself was not detected by IASI, but an abrupt drop in the average concentrations after the incident is seen in the satellite observations. In fact, after 2013 NH 3 enhancements are no longer detected at or near the plant. Figure 3 illustrates the processes of oversampling, supersampling and wind rotation on IASI data from 2007 to 2013. Each subpanel depicts the 120 km×60 km area centered around the plant, from top to bottom.
a. Gridded average. In the regular gridded average, each grid cell is assigned the arithmetic average of all observations whose center falls into the grid cell. This method only gives a faithful representation for larger grid-cell sizes, whereas smaller grid sizes provide a higher resolution at the cost of larger noise. Here a grid size of 0.15 • × 0.15 • was chosen. NH 3 enhancements are seen in a wide area around the plant, with a maximum west of the plant of 1 × 10 16 molec. cm −2 .
b. Oversampled average. Oversampling the daily maps before averaging increases the resolution and reveals the point-source nature of the emission, with a maximum close to the plant (around 1.6 × 10 16 molec. cm −2 ).
c. Wind-rotated oversampling. The wind-rotation technique (Valin et al., 2013;Fioletov et al., 2015) consists of rotating the map of daily observations around a presumed point source and along the direction of the wind direction at that point. The rotation was applied here to align the winds in the x direction. Daily horizontal wind fields were taken from the ERA5 reanalysis (ERA5, 2019) and interpolated at an altitude equal to half of the boundary layer height. The NH 3 average shown in this panel was obtained via oversampling applied to all the daily wind-rotated maps. It is important to note that such a distribution can no longer be interpreted as a geographical map, since each pixel is an average of measurements taken at different places. The only map element that is preserved is the distance to the point source. However, looking at the resulting distribution, the advantages brought by wind rotation are obvious. Whereas in the normal oversampled average the NH 3 enhancements are scattered across, aligning the winds significantly enhances both the source and transport (with a maximum of 2 × 10 16 molec. cm −2 ).
d. Wind-adjusted supersampling (i). The figure in this panel was obtained from wind-rotated daily maps, as in the previous panel, but this time the average was calculated with three iterations of the IBP supersampling algorithm. As explained above, supersampling offers most benefits when the underlying distribution can be assumed to be reasonably constant, which is in part achieved by aligning the winds. The resolution is further increased, and as the plume is much less smoothed out, maximum observed columns are also much higher (3.3 × 10 16 molec. cm −2 ). Note that in general for NH 3 , three iterations of the IBP algorithm seem to offer a good compromise of increasing the resolution of the average without introducing artifacts related to overfitting.
e. Wind-adjusted supersampling (ii). In Fig. 3c and d, the point-source location was taken from Van Damme et al. (2018), where the locations were determined based on the location of the maxima in the oversampled averages. In this last panel, the rotation was applied around the center of the presumed source (the chemical plant). The performance of the wind rotation is further enhanced, yielding a distribution fully consistent with that of a single emitting point source whose emissions undergo transport in a fixed direction. The part of the plume located furthest from the source is a bit off axis, which is probably caused by inhomogeneities in the wind fields across the entire scene. This panel also illustrates the sensitivity of the rotation method to small shifts in the location of the center, a fact that we will exploit in the next section.
One useful property of the different procedures is that they all approximately conserve the quantity that is being averaged; i.e., the averaged quantity in each grid is approximately the same as the average quantity in the grid representing the ground truth given that there is a sufficient number of measurements across the entire grid. When the number of measurements is low, this can break down dramatically, as can be seen with the extreme example of a single high-value measurement over an isolated point source. When a gridded average is made from this single measurement onto a coarse grid (e.g., 5 • × 5 • ), the entire grid cell containing the measurement will be associated with this high value, thus overestimating reality. A strict conservation is therefore not possible in general, as not enough information is contained in the original measurements to reconstruct the ground-truth perfectly even on average. That being said, supersampling conserves quantity with respect to the original measurements, when the number of iterations is large enough. This is a consequence of the fact that the back-projected measurements converge to the actual measurements and therefore also their averages. Finally note that wind rotation does not alter the grid average, as rotation simply redistributes the measurements to different locations on the grid. The average total NH 3 columns are indicated in each subpanel of Fig. 3. The average of the ungridded measurements within the considered box equals 5.23 × 10 15 molec. cm −2 . As can be seen, the largest change in the average column is caused by the rotation procedure, but this is simply an artifact caused by limiting the average to a square box around a point source (instead of a circle). This example illustrates that in practice, with differences smaller than 1 %, the different gridding procedures can be assumed to conserve quantity.

An NH 3 point-source map
Having demonstrated the effectiveness of both the windrotation and supersampling approaches in revealing point sources, we are now in a position to introduce a new type of NH 3 map, specifically designed to track down point sources. It is based on a similar map presented in McLinden et al. (2016) for SO 2 , but some important differences were introduced here to make it work for NH 3 . The main idea of McLinden et al. (2016) is to treat each location on Earth as a potential point source and to assign it a value proportional to the downwind (the source) minus the upwind (the background) column. In particular, for a given location, a wind-rotated average is constructed first, similar to Fig. 3c. Representative average columns are then obtained downwind and upwind from the potential source (e.g., in boxes of 10 km×10 km). Finally, the difference of the upwind and downwind average is calculated, and this value is then used to represent the point-source column at that specific location. While the method works nicely for SO 2 , this method proved to be only moderately successful when we applied it to the IASI NH 3 data. In particular, for those places where area sources dominate or where point sources are clustered over too large an area, local variation in the columns produces a noisy map, with many fictitious point sources.
We found that instead of the differences, the downwind average alone produced a more representative point-source map. In addition, applying the method not on the oversampled average but on the supersampled one allows increasing the resolution. There are two key advantages offered by a point-source map constructed in this way as opposed to a regular oversampled average: brighter point sources and smoother (lower) values over the background. The fact that point sources appear brighter is a direct consequence of the plume concentration achieved with wind-rotated supersampling, as shown in the previous section. Smoothing of the background is accomplished by the process of averaging the area downwind. However, by applying the method not on an oversampled average but on a supersampled one, this smoothing is partially offset for point sources. The resulting point-source map has a similar horizontal resolution to the oversampled map, but with increased averaged columns at the point sources and a smoother background distribution.
Examples over two selected regions in North America are shown in Fig. 4b and d. In these examples, the downwind averages were calculated in boxes extending from −5 to 5 km in the y direction and 0 to 20 km in the x direction. Figure 4a and c correspond to the regular oversampled averages. In Fig. 4a, which shows the oversampled average of the southern part of the Saskatchewan province of Canada, no point sources are apparent in the patchy NH 3 distribution. The cor-Atmos. Meas. Tech., 12, 5457-5473, 2019 www.atmos-meas-tech.net/12/5457/2019/ responding point-source map, shown in Fig. 4, is smoother over areas dominated by the diffuse sources, where column variations are close to the measurement uncertainty. In addition, two bright spots are evident, which upon investigation coincide with the location of an ammonia plant (Belle Plaine) and a very large feedlot (> 2 km in length) near the town of Lanigan. Looking back at the oversampled average, even with the advantage of hindsight, these sources can hardly be singled out. Figure 4c and d show the southwestern part of Kansas, US. It is an area well known for its cattle (Harrington and Lu, 2002). In Van Damme et al. (2018) several point sources associated with feedlots were isolated in Kansas and the rest of the High Plains region, but most of the area was found to be too diffuse to allow identification of individual point sources. The NH 3 point-source map facilitates greatly the attribution of these. This is due to the reduction in noise and the fact that the main point sources contrast much more with the background. An added benefit of this is that the location of the maxima in the map is in general closer to the actual emission source than is the case in the oversampled map, making it easier to track down the suspected source with visible imagery and therefore also to assign and identify the point source.
Displaced maxima that are seen in regular averages (windadjusted or not) can also be the result of transport, as noted by Van Damme et al. (2018), who found that especially for coastal sites, the shift can be as much as 20 km. The suspected reason is vertical uplift during transport, which makes NH 3 easier to detect and to measure (as can be seen in Fig. 3c) downwind of the source. The way in which the pointsource map is set up corrects for the effects of transport, as the columns are partially reallocated back to their source by assigning the average downwind column to the point source. We have quantified the ability to locate sources on a careful selection of 36 industrial emitters. These were all chosen to be relatively isolated, with no nearby other industries or other sources, so that the actual emitting source is known with confidence. In addition, only small-to medium-sized plants were chosen (< 1 km across) so that the precise location of the emission is known within a distance of about 500 m. For the regular oversampled map, the sources were found within a median distance of 3.9 km and a mean of 5.4 ± 3.7 km. The furthest distance was 15.2 km. For the point-source map, all but five sites were located within 3 km (with a median of 1.5 km, a mean of 2.1 ± 1.7 km and a maximum of 7.3 km), which confirms its improved performance in geo-allocation of the sources.
A final advantage of the point-source map is its performance in areas mildly affected by fires (e.g., in southeastern Asia, Mexico and parts of South America). Certain hotspots due to fires, with a plume center of around 25 km, can look just like actual point sources. In the point-source map, these often appear less bright and are blurred out over a wider area, with lower columns compared to the oversampled average. On the other hand, as before, actual point sources appear brighter and can emerge from the patchy NH 3 distribution that is characteristic for areas affected by fires. For that reason, comparing the oversampled and the point-source map was found to be very useful for singling out point sources, especially in those areas with larger background values. Example point sources are the ammonia plant in Campana (Argentina) and Bastos (Brazil), an important center of egg production. These were previously difficult to detect but are now easily identified.

Updated point-source catalog
Using the methodology presented in the previous section, NH 3 point-source maps of the world (land only) were constructed at a resolution of 0.01 • × 0.01 • (corresponding to a horizontal resolution of the order of 1-2 km). A few such maps were constructed by varying the size of the averaging box and the applied wind speeds (either in the middle of the boundary layer or at 100 m). While oversampling and back projection are computationally not that demanding, we recall that the construction is based on the treating each grid cell of the 0.01 • × 0.01 • map as a potential point source and therefore relies on the construction of wind-adjusted supersampled maps like Fig. 3e for each grid cell. Therefore, producing a world map at that resolution entails the generation of over 100 000 000 maps similar to Fig. 3e, each at a resolution of 1 km and each using more than 250 000 IASI measurements. A single point-source map therefore takes more than a month of computation. We decided to use all the available 2007-2017 NH 3 data both from IASI/Metop A (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) and IASI/Metop B (2013, which helps in reducing the noise even though it creates averages which are biased towards the last 5 years. The maps were then analyzed to provide an update of the point-source catalog presented in Van Damme et al. (2018). We refer to it for a detailed description of the methodology for the identification and categorization of the point sources, as we used the same method here. In brief, first, the global map is analyzed region per region in search of NH 3 hotspots that are no larger than 50 km across and that exhibit localized and concentrated enhancements compatible with a point source or dense cluster of point sources. Areas dominated by fires are excluded from the analysis. Analysis of areas with many sources or large ambient background concentrations, such as the Indo-Gangetic Plain, is severely hampered and reveals only the very large point sources. Isolated point sources in remote areas, on the other hand, can easily be picked up. The presence of a point source in the catalog should therefore not be seen as a quantitative indicator of its emission strength. Note that in this study we did not attempt to quantify the emission strengths of each. The categorization of the suspected point sources is performed using Google Earth imagery and thirdparty information (mainly inventories of fertilizer plants and online resources). The original categories were as follows: "agriculture", "fertilizer industry", "other industry", "natural" and "non-determined". Here we expanded the number of categories considerably and in particular introduced an urban category and subdivided "other industry" into five new categories, as detailed below.
The new point-source catalog is listed in Table A1 and illustrated on a world map in Fig. 5. Agricultural point sources were found to be invariably associated with CAFOs. Their number more than doubled, from 83 to 216, largely due to the increased attribution in areas of densely located point sources. For many of the previously tagged "source regions", it was possible to resolve large individual emitters. This was the case in the US (particularly in the geographical region that corresponds to the High Plains Aquifer), Mexico and along the coast of Peru. Also notable are several newly exposed large feedlots in Canada and in eastern Australia. For Atmos. Meas. Tech., 12, 5457-5473, 2019 www.atmos-meas-tech.net/12/5457/2019/ the first time, agricultural point sources were also identified in China and Russia. Industrial point sources, as before, are mainly associated with ammonia or urea-based fertilizer production (216, coming from 132) in Europe, northern Africa and Asia. Industrial hotspots were categorized as fertilizer industry as soon as evidence was found of fertilizer production, even when there are clearly other industries present that may contribute. Separate categories were introduced for the previously identified soda ash, geothermal, nickel-mining and coking industries, as additional examples were found for each. One ammonia plant in the US, associated with the manufacturing of explosives, was also assigned a separate category. Emissions over unidentified industries were labeled as a non-determined industry.
An important new category is the urban one. Previously, localized emissions near Mexico City, Bamako (Mali) and Niamey (Niger) were noted. While these hotspots represent diffuse sources, they have been included in the catalog, as the extent of the emissions in the relevant cities is sufficiently local and sufficiently in excess of background values. Thanks to the improved source representation, clear enhancements were found in Kabul (Afghanistan) and 12 African urban agglomerations: Ouagadougou (Burkina Faso), Bamako (Mali), Kano (Nigeria), Niamey (Niger), Maiduguri (Nigeria), Khartoum-Omdurman (Republic of the Sudan), Luanda (Angola), Kinshasa (Congo), Nairobi (Kenya), Addis Ababa (Ethiopia), Bamako (Mali) and Kampala (Uganda). Especially in Asia, atmospheric NH 3 is found in excess over most megacities and with much larger columns than found in these African megacities. However, because of the much larger background columns, and much denser clusters of cities, these could not be singled out as was the case in Africa. Apart from industry, known urban sources of NH 3 include emissions from vehicles, human waste (waste treatment and sewers), biological waste (garbage containers) and domestic fires (including waste incineration) (Adon et al., 2016;Sun et al., 2017;Reche et al., 2015). At least the hotspot at Bamako is consistent with in situ measurements (Adon et al., 2016), which report very high concentrations year-round, between 28 and 73 ppb on a monthly averaged basis. Local conditions surely are key in explaining why some cities in Africa exhibit much larger concentrations than others. Johannesburg (South Africa) for instance blends in completely in the background, with ambient values barely larger than in the rest of South Africa. While this is outside the scope of the paper, there is no doubt that the IASI NH 3 data could be further exploited www.atmos-meas-tech.net/12/5457/2019/ Atmos. Meas. Tech., 12, 5457-5473, 2019 to better understand the driving factors of urban emission on a global scale. Other than at Lake Natron , no other natural NH 3 hotspots have been identified. For a number of presumed point sources, no likely source could be attributed; however given their location (central US, Middle East and eastern Asia), these are most likely anthropogenic.

Conclusions
Oversampling is a technique now commonplace in the field of atmospheric sounding for achieving hyperresolved spatial averages beyond what the satellites natively offer. There is a class of algorithms referred to as superresolution that goes beyond oversampling, but these have until now only been applied to measurements of satellite imagers for surface parameters. Here, we have shown that it is a viable method that can also be applied to the single-pixel images taken by atmospheric sounders for short-lived gases. We demonstrated this with measurements of a quartz emissivity feature over the Sahara, for which a spatial resolution down to 2-3 km could be achieved.
Superresolution is a priori less suitable for measurements of atmospheric gases because of variations in their distribution related to variations in transport. However, by aligning the winds around point-source emitters, much of this variability can be removed. In Sect. 3, we have shown the advantage of applying IBP superresolution on such wind-corrected data. The resulting averaged plumes originating from point sources not only reveal more detail, but maximum concentrations and gradients are also larger and presumably more realistic. Studies of atmospheric lifetime (e.g., Fioletov et al., 2015), which rely on the precise shape of the dispersion, could potentially benefit from this increase in accuracy.
Wind-adjusted superresolution images around point sources form the basis of the NH 3 point-source map, which is an NH 3 average that simultaneously corrects for wind transport, accentuates point sources and smooths area sources. It was inspired by the SO 2 "difference" map presented in McLinden et al. (2016), but as we do not look at differences, the NH 3 map still looks like an NH 3 total column distribution. However, other than for the identification of point sources, such a map is not easily exploitable, as it is a distorted representation of the reality that favors point sources. In-depth analysis allowed us to perform a major update of the global catalog of point sources presented in Van Damme et al. (2018), with more than 500 point sources identified and categorized. As a whole, this study further highlights the importance of point sources on local scales. The world map shows distinct patterns, with agricultural point sources completely dominant in America, in contrast to Europe and Asia, where industrial point sources are prevalent. In Africa, NH 3 hotspots are mainly found near urban agglomerations.
While the point-source catalog was established with a great deal of care, given its size, mistakes will inevitably be present both in the localization of the point sources (due to, for example, noise in the data or NH 3 in transport) and in the categorization. Improvements can probably best be achieved with feedback from the international community, with complementary knowledge on regional sources. For this reason, and to keep track of emerging emission sources, we have set up a website with an interactive global map, visualizing the distribution, type and time evolution of the different point sources (http://www.ulb.ac.be/cpm/NH3-IASI.html, last access: 11 October 2019). With the help of the community, we hope it can become a useful resource for information on global NH 3 point sources.
Data availability. The IASI NH 3 product is available from the Aeris data infrastructure (http://iasi.aeris-data.fr, last access: 11 October 2019). It is also planned to be operationally distributed by EU-METCast under the auspices of the Eumetsat Atmospheric Monitoring Satellite Application Facility (AC-SAF; http://ac-saf.eumetsat. int, last access: 11 October 2019).
Atmos. Meas. Tech., 12, 5457-5473, 2019 www.atmos-meas-tech.net/12/5457/2019/ Appendix A: Point-source catalog Table A1. Updated point-source catalog. The categories are abbreviated as A for agriculture, CI for the coking industry, EI for the explosives industry, FI for the fertilizer industry, GI for the geothermal industry, NDI for a non-determined industry, NI for the nickel industry, SI for the soda ash industry, N for natural, ND for non-determined and U for urban.