# Weekly derived top-down volatile-organic-compound fluxes over Europe from TROPOMI HCHO data from 2018 to 2021

Glenn-Michael Oomen, Jean-François Müller, Trissevgeni Stavrakou, Isabelle de Smedt, Thomas Blumenstock, Rigel Kivi, Maria Makarova, Mathias Palm, Amelie Röhling, Yao Té, et al.

► To cite this version:

Glenn-Michael Oomen, Jean-François Müller, Trissevgeni Stavrakou, Isabelle de Smedt, Thomas Blumenstock, et al. Weekly derived top-down volatile-organic-compound fluxes over Europe from TROPOMI HCHO data from 2018 to 2021. Atmospheric Chemistry and Physics, 2024, 24, pp.449-474. 10.5194/acp-24-449-2024. insu-04822497

HAL Id: insu-04822497

https://insu.hal.science/insu-04822497v1

Submitted on 6 Dec 2024

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

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

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

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

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

Weekly derived top-down volatile-organic-compound fluxes over Europe from TROPOMI HCHO data from 2018 to 2021

Glenn-Michael Oomen¹, Jean-François Müller¹, Trissevgeni Stavrakou¹, Isabelle De Smedt¹,
Thomas Blumenstock², Rigel Kivi³, Maria Makarova⁴, Mathias Palm⁵, Amelie Röhling², Yao Té⁶,
Corinne Vigouroux¹, Martina M. Friedrich¹, Udo Frieß⁷, François Hendrick¹, Alexis Merlaud¹,
Ankie Piters⁸, Andreas Richter⁵, Michel Van Roozendael¹, and Thomas Wagner⁹

¹Royal Belgian Institute for Space Aeronomy (BIRA-IASB), Brussels, Belgium

²Karlsruhe Institute of Technology (KIT), IMK-ASF, Karlsruhe, Germany

³Finnish Meteorological Institute (FMI), Sodankylä, Finland

⁴Atmospheric Physics Department, Saint Petersburg State University, St. Petersburg, Russia

⁵Institute of Environmental Physics, University of Bremen, Bremen, Germany

⁶LERMA-IPSL, Sorbonne Université, CNRS, Observatoire de Paris, PSL Université, 75005 Paris, France

⁷Institute of Environmental Physics, Heidelberg University, Heidelberg, Germany

⁸Royal Netherlands Meteorological Institute (KNMI), De Bilt, the Netherlands

⁹Max Planck Institute for Chemistry (MPI-C), Mainz, Germany

**Correspondence:** Glenn-Michael Oomen (glenn-michael.oomen@aeronomie.be)

Received: 31 August 2023 – Discussion started: 19 September 2023

Revised: 10 November 2023 – Accepted: 13 November 2023 – Published: 11 January 2024

**Abstract.** Volatile organic compounds (VOCs) are key precursors of particulate matter and tropospheric ozone. Although the terrestrial biosphere is by far the largest source of VOCs into the atmosphere, the emissions of biogenic VOCs remain poorly constrained at the regional scale. In this work, we derive top-down biogenic emissions over Europe using weekly averaged TROPOMI formaldehyde (HCHO) data from 2018 to 2021. The systematic bias of the TROPOMI HCHO columns is characterized and corrected for based on comparisons with FTIR data at seven European stations. The top-down fluxes of biogenic, pyrogenic, and anthropogenic VOC sources are optimized using an inversion framework based on the MAGRITTEv1.1 chemistry transport model and its ad joint. The inversion leads to strongly increased isoprene emissions with respect to the MEGAN–MOHYCAN inventory over the model domain (from 8.1 to 18.5 Tg yr⁻¹), which is driven by the high observed TROPOMI HCHO columns in southern Europe. The impact of the inversion on biomass burning VOCs (+13 %) and anthro pogenic VOCs (−17 %) is moderate. An evaluation of the optimized HCHO distribution against ground-based remote sensing (FTIR and MAX-DOAS) and in situ data provides generally improved agreement at stations below about 50° N but indicates overestimated emissions in northern Scandinavia. Sensitivity inversions show that the top-down emissions are robust with respect to changes in the inversion settings and in the model chem ical mechanism, leading to differences of up to 10 % in the total emissions. However, the top-down emissions are very sensitive to the bias correction of the observed columns, as the biogenic emissions are 3 times lower when the correction is not applied. Furthermore, the use of different a priori biogenic emissions has a significant impact on the inversion results due to large differences among bottom-up inventories. The sensitivity run using CAMS-GLOB-BIOv3.1 as a priori emissions in the inversion results in 30 % lower emissions with respect to the optimization using MEGAN–MOHYCAN. In regions with large temperature and cloud cover variations, there is strong week-to-week variability in the observed HCHO columns. The top-down emissions, which are optimized at weekly increments, have a much improved capability of representing these large fluctuations than an inversion using monthly increments.
# 1 Introduction

The emissions of non-methane volatile organic compounds (NMVOCs) play an important role for both air quality and climate. They are precursors of secondary organic aerosols (SOAs; Claeys et al., 2004; Hallquist et al., 2009) and also impact tropospheric ozone levels through their interaction with $NO_x$ (Houweling et al., 1998; Churkina et al., 2017). Globally, the largest source of NMVOCs is of biogenic origin, making up about 85% of the total emissions as compared to 12% and 3% for anthropogenic and pyrogenic emissions, respectively (Stavrakou et al., 2018; Granier et al., 2019). Among the biogenic volatile organic compounds (BVOCs), isoprene ($C_5H_8$) makes up about half of the global emissions, making it the most important NMVOC (Guenther et al., 2012; Sindelarova et al., 2014). Even though isoprene emission models agree on several important features such as temperature and light density dependence, seasonality, and interannual variability, there are still large inconsistencies among different inventories regarding the annual global emission of isoprene, ranging from ~300 to ~600 Tg yr$^{-1}$ (Sindelarova et al., 2022). The main drivers of these discrepancies come from differences in the climate inputs, the emission factors, and the vegetation fields entering the model (Arneth et al., 2011).

Isoprene has a strong influence on tropospheric chemistry due to its large emissions and its fast photochemical sink, mainly through its reaction with OH radicals. The latter leads to the formation of various organic compounds through a long cascade of chemical reactions (Wennberg et al., 2018; Müller et al., 2019), some of which remain poorly characterized, especially in remote environments characterized by low levels of nitrogen oxides ($NO_x$, defined as $NO + NO_2$). Isoprene has an average lifetime of ~1 h in the troposphere (Bates and Jacob, 2019). One of its highest-yield oxidation products is formaldehyde (HCHO), a compound that can be readily measured using both ground- and space-based detectors. The main source of HCHO is the oxidation of volatile organic compounds (VOCs), whereas direct emission provides only a small contribution. Methane oxidation accounts for approximately 60% of the global HCHO budget (Stavrakou et al., 2009), but due to its long lifetime, it functions as a background source. Indeed, short-lived VOCs, such as isoprene and monoterpenes, have a larger impact on the spatial distribution of HCHO. Additionally, HCHO is also relatively quickly destroyed through photolysis and reaction with OH and has a lifetime of approximately 5 h (Stavrakou et al., 2015). Transport processes are relatively unimportant for its measured distribution because of its short lifetime and of the rapid chemical processing leading to HCHO production (except at night, during winter, or in $NO_x$-poor environments; see Marais et al., 2012). Consequently, HCHO is an

excellent tracer, and its total columns can be directly linked to the emissions of reactive NMVOCs.

Since the advent of satellite HCHO data products, many studies have been conducted on the derivation of top-down biogenic emissions in an inversion framework. This includes studies performed using satellite HCHO data from GOME (e.g., Palmer et al., 2006), SCIAMACHY (e.g., Stavrakou et al., 2009), OMI (e.g., Millet et al., 2008; Stavrakou et al., 2018), and GOME-2 (e.g., Stavrakou et al., 2015). HCHO is a trace gas with a low optical density in the atmosphere, and its absorption features in the ultraviolet (UV) are heavily blended with BrO and $O_3$ lines. Consequently, HCHO vertical columns from satellite measurements are relatively noisy, and substantial spatial and temporal averaging is needed in order to derive reliable top-down emissions.

In this study, we use the HCHO product of the Tropospheric Ozone Monitoring Instrument (TROPOMI; De Smedt et al., 2018) mounted on the Sentinel-5P (S5P) satellite as observational constraints in the top-down inversion. Due to its unprecedented spatial resolution and high signal-to-noise ratio, we use weekly averages of the data at 0.5° × 0.5° spatial resolution, whereas previous studies used averages at monthly or longer timescales. Because the emissions of biogenic VOCs are strongly temperature-dependent, HCHO volume mixing ratios can show strong temporal variability in the course of 1 month. We can expect that limiting the extent of temporal averaging will significantly improve the constraints in the inversion.

We focus our study on the European domain as depicted in Fig. 1. This domain contributes only a small amount to the global biogenic VOC budget, since the highest isoprene emissions take place in warmer, more densely vegetated regions such as tropical forests. Nonetheless, the quantification of the BVOC emissions in Europe is important, since the high levels of anthropogenic $NO_x$ and sulfate particles make BVOCs more prone to the production of tropospheric ozone (Seinfeld and Pandis, 2016) and secondary organic aerosols (Schwantes et al., 2019; Bryant et al., 2020). Furthermore, the high amount of validation data available in the domain, as is shown in Fig. 1, makes this region an excellent target for our inversion study. Previous top-down VOC emission studies over Europe using HCHO data include Dufour et al. (2009) and Curci et al. (2010), which used the CHIMERE model in combination with SCIAMACHY and OMI observations, respectively. Their inversion studies resulted in slightly reduced isoprene emissions with respect to the a priori estimates obtained from the MEGAN bottom-up inventory (Guenther et al., 1995, 2006). More precisely, Curci et al. (2010) found a decrease in isoprene emissions from 3.16 to 3.04 Tg in the domain (35–75° N, 15° W–25° E) for the months May to September in 2005. In a global-scale inversion study by Bauwens et al. (2016) using OMI HCHO data from 2005 to 2013 and the IMAGESv2 model, a moder-
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Locations of ground-based HCHO column stations and in situ HCHO measurement sites in Europe presented in Table 1.</figcaption></figure>

ate increase in isoprene emissions was inferred in both eastern and western Europe, with a total increase over the European domain from 6.8 to 8.4 Tg yr⁻¹.

The inversion of VOC emissions is sensitive to constraints from satellite observations. Consequently, it is important to carefully validate the satellite measurements using ground-based and aircraft measurements. Several validation studies have shown that OMI and TROPOMI HCHO columns are too low over source regions with respect to validation datasets (Zhu et al., 2020; Vigouroux et al., 2020; De Smedt et al., 2021). This bias was not accounted for in previous HCHO emission inversion studies. The origin of the bias is unknown, although some known difficulties in the HCHO retrieval include aerosol effects, uncertainties in the treatment of clouds, spectral interferences with BrO and ozone, and a coarse albedo climatology (Vigouroux et al., 2020). The bias can be corrected for based on a linear regression of the ground-based and satellite vertical columns. In this work, we derive a bias correction based on ground-based FTIR measurements at stations located in the European domain.

We describe the different types of data (ground-based and satellite) used in this work in Sect. 2. Section 3 describes the setup of the inversion, as well as the chemistry transport model MAGRITTEv1.1 around which the inversion framework is built. A detailed characterization of the biases of the TROPOMI HCHO data product is presented in Sect. 4 based on FTIR and MAX-DOAS data over Europe. Based on this analysis, a procedure for correcting the biases is proposed and used in the emission inversions. The inversion results are shown in Sect. 5, including a comparison of the top-down model with ground-based data. In Sect. 6, we provide an analysis of the emission uncertainties by means of sensitivity inversions. Finally, we present the conclusions in Sect. 7.

## 2 HCHO data

In this section, we provide an overview of the measurement datasets used in this work. TROPOMI HCHO column

data are used to constrain the VOC emissions in the model (Sect. 2.1). We perform a careful validation of the satellite data using ground-based FTIR vertical columns from a range of stations across Europe (Sect. 2.2). We further validate the top-down results by including European MAX-DOAS stations in our ground-based remote sensing dataset (Sect. 2.3). Finally, in situ HCHO concentration data are also used to evaluate the model before and after inversion (Sect. 2.4). The full list of ground-based stations used in this work is shown in Table 1 and Fig. 1.

### 2.1 TROPOMI HCHO column densities

The S5P satellite has a nadir-viewing orientation in a low-Earth, Sun-synchronous polar orbit. It has a daily equatorial overpass at 13:30 LT. The TROPOMI instrument, onboard the S5P platform, is an imaging spectrometer covering ultraviolet, visible, near-infrared, and short-wavelength infrared ranges. In the UV-VIS range, the initial spatial resolution of the instrument was 3.5 × 7 km², which was further improved to 3.5 × 5.5 km² after 6 August 2019.

The HCHO retrieval relies on the differential optical absorption spectroscopy (DOAS) method, which allows for the isolation of narrow trace gas absorption features. The TROPOMI retrieval algorithm is based on the QA4ECV algorithm developed for its predecessor OMI (De Smedt et al., 2018). The first step is the derivation of the slant column density in the wavelength range 328.5–359 nm. The DOAS reference spectrum is updated daily with an average of Earth reflected radiances measured in a remote sector of the equatorial Pacific Ocean. The slant columns are converted into tropospheric vertical columns using altitude-dependent air mass factors, which are obtained from a look-up table calculated with the VLIDORT v2.6 radiative transfer model (Spurr et al., 2008), for a large number of atmospheric conditions depending on the observation angles, surface reflectivity, surface pressure, and cloud properties. A priori vertical profiles are provided by the TM5-MP daily analysis on 34 vertical levels at a spatial resolution of 1° (Williams et al., 2017). The surface albedo is taken from the monthly OMI albedo climatology at 0.5° (minimum Lambertian equivalent reflectivity, LER; Kleipool et al., 2008). The final step of the algorithm consists of a background correction of the columns, using again the Pacific Ocean as a reference sector. The background HCHO columns, mostly due to methane oxidation, are replaced by the TM5-MP model columns in the same sector (De Smedt et al., 2021).

We use the operational offline level 2 data product (https://doi.org/10.5270/S5P-tjlxdf2), which provides HCHO tropospheric column data together with uncertainty estimates, as well as a priori profiles and averaging kernels (De Smedt et al., 2018, 2021). The latter provide information on the vertical sensitivity of the satellite measurement, which is necessary to make meaningful comparisons with our model. We only select data with a quality value higher
<table><caption>**Table 1.** List of stations providing ground-based remote sensing measurements (FTIR and MAX-DOAS) and in situ measurements.</caption><thead><tr><th colspan="8">FTIR stations</th></tr><tr><th>Station</th><th>Code</th><th>Latitude (° N)</th><th>Longitude (° E)</th><th>Altitude (m)</th><th>Years</th><th>No. of days</th></tr></thead><tbody><tr><td>Bremen</td><td>BRE</td><td>53.10</td><td>8.85</td><td>30</td><td>2018–2020</td><td>23</td></tr><tr><td>Garmisch</td><td>GAR</td><td>47.48</td><td>11.06</td><td>743</td><td>2018–2021</td><td>141</td></tr><tr><td>Karlsruhe</td><td>KAR</td><td>49.10</td><td>8.42</td><td>110</td><td>2018–2020</td><td>136</td></tr><tr><td>Kiruna</td><td>KIR</td><td>67.84</td><td>20.40</td><td>420</td><td>2018–2021</td><td>116</td></tr><tr><td>Paris</td><td>PAI</td><td>48.97</td><td>2.37</td><td>60</td><td>2018–2021</td><td>94</td></tr><tr><td>Sodankylä</td><td>SOD</td><td>67.37</td><td>26.63</td><td>190</td><td>2018–2021</td><td>203</td></tr><tr><td>St. Petersburg</td><td>SPE</td><td>59.88</td><td>29.83</td><td>20</td><td>2018–2021</td><td>94</td></tr></tbody><thead><tr><th colspan="8">MAX-DOAS stations</th></tr><tr><th>Station</th><th>Code</th><th>Latitude (° N)</th><th>Longitude (° E)</th><th>Altitude (m)</th><th>Years</th><th>No. of days</th></tr></thead><tbody><tr><td>Athens</td><td>ATH</td><td>38.05</td><td>23.86</td><td>532</td><td>2018–2021</td><td>388</td></tr><tr><td>Bremen</td><td>BRE</td><td>53.10</td><td>8.85</td><td>46</td><td>2018–2021</td><td>301</td></tr><tr><td>Cabauw</td><td>CAB</td><td>51.97</td><td>4.93</td><td>3</td><td>2021</td><td>67</td></tr><tr><td>De Bilt</td><td>DEB</td><td>52.10</td><td>5.18</td><td>22</td><td>2020–2021</td><td>87</td></tr><tr><td>Heidelberg</td><td>HEI</td><td>49.42</td><td>8.67</td><td>145</td><td>2018–2021</td><td>314</td></tr><tr><td>Mainz</td><td>MAZ</td><td>49.99</td><td>8.23</td><td>150</td><td>2018–2021</td><td>220</td></tr><tr><td>Uccle</td><td>UCC</td><td>50.80</td><td>4.36</td><td>95</td><td>2018–2021</td><td>139</td></tr></tbody><thead><tr><th colspan="8">In situ stations</th></tr><tr><th>Station</th><th>Code</th><th>Latitude (° N)</th><th>Longitude (° E)</th><th>Altitude (m)</th><th>Years</th><th>No. of days</th></tr></thead><tbody><tr><td>Lahemaa</td><td>LAH</td><td>59.50</td><td>25.9</td><td>32</td><td>2018–2019</td><td>201</td></tr><tr><td>Peyrusse Vieille</td><td>PEY</td><td>43.62</td><td>0.18</td><td>200</td><td>2018–2020</td><td>280</td></tr><tr><td>San Pablo de los Montes</td><td>PAB</td><td>39.55</td><td>-4.35</td><td>917</td><td>2018–2021</td><td>209</td></tr><tr><td>La Tardière</td><td>TAR</td><td>46.65</td><td>-0.75</td><td>133</td><td>2018</td><td>55</td></tr><tr><td>Birkenes</td><td>BIR</td><td>58.38</td><td>8.25</td><td>190</td><td>1994–2005</td><td>613</td></tr><tr><td>Brotjacklriegel</td><td>BRO</td><td>48.82</td><td>13.22</td><td>1016</td><td>1999–2004</td><td>477</td></tr><tr><td>Campisábalos</td><td>CAM</td><td>41.27</td><td>-3.14</td><td>1360</td><td>2004–2009</td><td>551</td></tr><tr><td>Donon</td><td>DON</td><td>48.50</td><td>7.13</td><td>775</td><td>1994–1999</td><td>481</td></tr><tr><td>Ispra</td><td>ISP</td><td>45.80</td><td>8.63</td><td>209</td><td>1994–1997</td><td>288</td></tr><tr><td>Košetice</td><td>KOS</td><td>49.57</td><td>15.08</td><td>535</td><td>1994–2005</td><td>1097</td></tr><tr><td>Montelibretti</td><td>MON</td><td>42.10</td><td>12.63</td><td>48</td><td>1994–1996</td><td>122</td></tr><tr><td>Schmücke</td><td>SCH</td><td>50.65</td><td>10.77</td><td>937</td><td>1999–2004</td><td>516</td></tr><tr><td>Waldhof</td><td>WAL</td><td>52.80</td><td>10.76</td><td>74</td><td>1999–2004</td><td>535</td></tr><tr><td>Zingst</td><td>ZIN</td><td>54.44</td><td>12.72</td><td>1</td><td>1999–2004</td><td>533</td></tr></tbody></table>

than 0.5, as recommended by the S5P HCHO product user manual (https://sentinel.esa.int/documents/247904/2474726/Sentinel-5P-Level-2-Product-User-Manual-Formaldehyde, last access: 21 August 2023). We use the clear-sky product, meaning that the cloud information is only used to filter the satellite observations, and the cloud correction to the air mass factors is ignored (De Smedt et al., 2021). To ensure that cloud contamination has a minimal effect on the total columns, we impose that the cloud fraction for each pixel is lower than 0.2, which is stricter than the standard cloud filter (fraction of 0.4 when using a quality value of 0.5).

Figure 2 shows the TROPOMI HCHO vertical columns over Europe for 4 weeks in July 2019, gridded spatially onto the model resolution (0.5° × 0.5°). High temperatures and high levels of solar radiation generally favor biogenic

emissions and HCHO production from long-lived VOC precursors. On average, this leads to higher HCHO columns in southern Europe due to higher temperatures and radiation levels with decreasing latitude. Local minima generally correspond to mountain ranges, oceans, and arid or semi-arid areas. Meteorology plays an important role in the spatio-temporal variability in the observed columns (Fig. 2). Most notably, the high HCHO columns in central Europe in the last week of July 2019 are due to a record-breaking heat wave (Sousa et al., 2020). Weekly averaged TROPOMI HCHO columns are used as observational constraints in our inversion framework. At the resolution of 0.5°, and using weekly averaged columns, the precision of the TROPOMI HCHO columns is estimated to range between 1 and 3 × 10¹⁵ molec.cm⁻², from background to continental condi-
tions (De Smedt et al., 2021). The bias of the HCHO vertical columns is estimated to amount to -30% of the column for high HCHO levels and to 30% of the column for low HCHO levels (Vigouroux et al., 2020). A detailed derivation of the bias over the target domain is presented in Sect. 4.

## 2.2 FTIR column observations

The FTIR stations used in this study (Table 1) are part of the international networks TCCON (Total Carbon Column Observing Network; http://www.tccon.caltech.edu/, last access: 21 August 2023) and/or NDACC (Network for the Detection of Atmospheric Composition Change, https://www2.acom.ucar.edu/irwg, last access: 21 August 2023). They record remote sensing solar absorption measurements in the infrared domain, under clear-sky conditions, using high-resolution spectrometers from the same manufacturer (Bruker 120/5HR). The retrievals of atmospheric total columns and low-resolution vertical profiles are performed using either one of the two algorithms available in the community (PROFITT9, Hase et al., 2006, and SFIT4, Pougatchev et al., 1995), both based on similar line-by-line forward models and optimal-estimation methods (Rodgers, 2000) and therefore providing very consistent outputs.

To build a coherent dataset useful for satellite validation or model evaluation, the retrieval settings of HCHO have been harmonized within the FTIR community (Vigouroux et al., 2018). This dataset has been used for the validation of TROPOMI (Vigouroux et al., 2020, updated in quarterly reports at https://mpc-vdaf.tropomi.eu, last access: 21 August 2023), GOME-2 (https://acsaf.org/valreps.php, last access: 21 August 2023), or more recently OMPS (Kwon et al., 2023). Currently 28 stations provide HCHO column data, among which 7 are in the European domain (Table 1, Fig. 1).

The details of the harmonized retrieval settings are documented in Vigouroux et al. (2018). In summary, the fitted spectral signatures lie in the 3.6 µm region and belong to the $v_1$ and $v_5$ bands. The spectroscopic database is a compilation from Geoffrey Toon, the atm16 line list (http://mark4sun.jpl.nasa.gov/toon/linelist/linelist.html, last access: 21 August 2023), which corresponds to HITRAN 2012 (Rothman et al., 2013) for HCHO and is optimized for each of the interfering species. The Tikhonov regularization (Tikhonov, 1963) is used for constraining/optimizing the retrievals of low-resolution vertical HCHO profiles. The degrees of freedom of signal (DOFSs; trace of the averaging kernel matrix) amount to only 1 to 1.4 for the stations used in this study, implying that only total column information can be retrieved. The averaging kernel shape shows that the sensitivity is highest in the free troposphere (although still present at the surface).

The FTIR dataset is well characterized in terms of uncertainty budget, which is calculated following Rodgers (2000). The systematic uncertainty on the HCHO columns is about 13% (Vigouroux et al., 2018) and is dominated by the uncertainty on the spectroscopic parameters. The random uncer-

tainties, dominated by the measurement noise, can be as low as $1.0 \times 10^{14}$ molec.cm$^{-2}$ (8%) for the pristine Eureka site and up to $5.3 \times 10^{14}$ molec.cm$^{-2}$ (7%) for the polluted Paris site. In addition, the FTIR data have significant random and systematic smoothing uncertainties (Vigouroux et al., 2018), but this component vanishes when the model profiles are smoothed using the FTIR averaging kernel matrix (Rodgers and Connor, 2003), as done in this study.

## 2.3 MAX-DOAS column observations

The MAX-DOAS (Multi-AXIs Differential Optical Absorption Spectroscopy) technique (Hönninger et al., 2004; Platt and Stutz, 2008) involves recording UV-visible spectra of the scattered skylight at different elevation angles between the horizon and the zenith. The different elevation angles yield different optical paths and thus different absorption strengths for trace gases present in the troposphere. By combining these multi-axis measurements with appropriate inversion schemes, it is possible to infer the vertical distribution and the integrated vertical column of several species absorbing in the UV-visible range, such as $NO_2$, HCHO, or ozone. MAX-DOAS measurements also allow the absorption of the oxygen collision complex ($O_4$) to be quantified, which yields information on the tropospheric profile of aerosol extinction. The latter is crucial for the assessment of the optical path and the quantitative interpretation of the measured trace gas absorptions.

Table 1 presents the MAX-DOAS stations used in this work. They rely on various types of MAX-DOAS instruments differing in size, geometry, and signal-to-noise ratio and are located in Germany, the Netherlands, Belgium, and Greece. Depending on the spectral ranges of the spectrometers, the absorption of HCHO is measured either in the 324 to 359 nm or 336 to 359 nm spectral window, while $O_4$ is fitted in the 338 to 370 nm window. The consistency of the HCHO and $O_4$ absorption measurements has been demonstrated by regular intercomparison studies (e.g., Kreher et al., 2020). For stations that have multiple viewing directions, we treat the different directions equally.

Regarding the MAX-DOAS data analysis, this work uses the facilities developed within the FRM4DOAS project (Fiducial Reference Measurements for Ground-Based DOAS Air-Quality Observations; https://frm4doas.aeronomie.be, last access: 23 August 2023). FRM4DOAS aims at harmonizing the data processing from MAX-DOAS instruments operated within the International Network for the Detection of Atmospheric Composition Change (NDACC) by incorporating existing retrieval algorithms into a fully traceable, automated, and quality-controlled processing chain.

The FRM4DOAS evaluation starts with the production of differential slant column densities by applying the QDOAS analysis tool (Danckaert and Fayt, 2017). The QDOAS settings for FRM4DOAS are described in Hendrick et al. (2018). The FRM4DOAS system implements two MAX-
<figure><img src="image_5.png" /><figcaption>**Figure 2.** TROPOMI vertical HCHO columns over Europe for the different weeks of July 2019, gridded at 0.5° × 0.5°. White regions correspond to grid cells without valid data over the week (usually due to clouds or snow/ice). Week 1 corresponds to days 1 to 8, week 2 corresponds to days 9 to 16, week 3 corresponds to days 17 to 24, and week 4 corresponds to days 25 to 31.</figcaption></figure>

DOAS retrieval algorithms: MAPA (Mainz Profile Algorithm; Beirle et al., 2019), which is based on a parameterization of the retrieval profile shape and a Monte Carlo approach for the inversion, and MMF (Mexican MAX-DOAS Fit; Friedrich et al., 2019), an optimal-estimation-based algorithm using the radiative transfer code VLIDORT v2.7 (Spurr et al., 2008) as a forward model. In this work, analogously to the official NO₂ product, only HCHO MMF data consistent with MAPA HCHO are retained. Both inversion algorithms have been extensively tested and validated using synthetic (Frieß et al., 2019) and real data (Tirpitz et al., 2021; Karagkiozidis et al., 2022). Vigouroux et al. (2009) quantified the uncertainties in individual MAX-DOAS measurements of HCHO vertical column densities (VCDs). The systematic uncertainties mainly originate from the HCHO cross-section and amount to 9% of the HCHO VCD. The smoothing error (20%) dominates the random error. After taking into account the retrieval noise (9%) and the errors in the forward model parameters (12%), this leads to a total random error of around 25%.

## 2.4 In situ

Additional data to evaluate our model are provided by ground-based in situ measurement stations. We use data from the EMEP network (Tørseth et al., 2012), available at the

EBAS database (https://ebas nilu.no). The list of stations with available HCHO concentration measurements is provided in Table 1. All in situ stations use high-performance liquid chromatography (HPLC) to derive HCHO concentrations at daily or 3 d intervals.

For the majority of in situ stations in Table 1, measurements were taken before the time span covered in this study (i.e., before 2018). For those stations, we compute a climatological average of the measurements by averaging the data per month across the different years. However, caution is required when comparing a climatological average to a recent model, since HCHO concentrations are sensitive to meteorology and land-use changes.

# 3 Inversion methodology

## 3.1 Model description

The inversion method is built around the Model of Atmospheric composition at Global and Regional scales using Inversion Techniques for Trace gas Emissions (MAGRIT-TEv1.1; Müller et al., 2019). This model is based on its predecessor IMAGES (Müller and Stavrakou, 2005; Stavrakou et al., 2012; Bauwens et al., 2016) and simulates the atmospheric composition at either the global or regional scale. Its chemical mechanism includes 141 long-lived and 41 short-
lived compounds. In particular, it contains a state-of-the-art chemical mechanism for the atmospheric degradation of isoprene, based on the Leuven Isoprene Mechanism (LIM1; Peeters et al., 2014), and includes updates from recent experimental studies (Wennberg et al., 2018; Berndt et al., 2019). This chemical mechanism is thoroughly described in Müller et al. (2019). The only important update with respect to the mechanism of Müller et al. (2019) concerns the isomerization rate of δ-hydroxyperoxy radicals from isoprene, which was based on a determination by the Caltech group (Wennberg et al., 2018). In this study, faster isomerization rates are used, based on the LIM1 determination of the 1,6 H-shift rate, combined with a 5-fold enhancement of the O₂ addition and elimination rates (also from LIM1). This enhanced rate was found necessary by Novelli et al. (2020) to reproduce their chamber measurements of HOₓ radicals and isoprene oxidation products in dedicated oxidation experiments in the simulation chamber SAPHIR.

Meteorological fields are obtained from the ECMWF ERA5 reanalysis data (Hersbach et al., 2020). The MAGRITTE model is first run in its global configuration at a spatial resolution of 2° × 2.5° (latitude, longitude). This simulation provides the boundary conditions for the MAGRITTEv1.1 model runs at 0.5° × 0.5° over the European domain delimited by 34–72° N, 12° W–45° E. The model has 40 equally spaced vertical levels, extending from the surface up to a pressure of about 44 hPa.

Anthropogenic VOC emissions are obtained from the CAMS-GLOB-ANT inventory (Granier et al., 2019, version 4.3); biomass burning emissions are obtained from the Quick Fire Emissions Dataset (QFED) version 2.4 (Darnenov and da Silva, 2015) combined with emission factors from Andreae (2019); and biogenic emissions of isoprene, monoterpenes, and methanol are obtained from the Model of Emissions of Gases and Aerosols from Nature (MEGAN) coupled to the multi-layer canopy environment model (MOHYCAN; Müller et al., 2008; Guenther et al., 2012; Stavrakou et al., 2018). The anthropogenic VOC emissions display a slow declining trend, from 16.8 Tg VOC in 2018 to 16.6 Tg VOC in 2021. The annual biomass burning fluxes over the European domain amount to 1.6 and 1.7 Tg VOC in 2018 and 2021, whereas the fluxes are much higher in 2019 (2.5 Tg VOC) and 2020 (2.6 Tg VOC). The annual a priori isoprene fluxes amount to 8.2, 7.9, 7.8, and 8.4 Tg in 2018, 2019, 2020, and 2021, respectively.

The biogenic emission fluxes are dependent on a range of environmental and phenological factors. In the MEGAN model, the biogenic emission flux F is calculated as

$$F = \epsilon \cdot C_{CE} \cdot \gamma_{A} \cdot \gamma_{CO_{2}} \cdot \gamma_{SM} \cdot \sum_{k=1}^{8} [(\text{LAI})_{k} \cdot (\gamma_{PT})], \quad (1)$$

where $\epsilon$ represents the standard emission rates obtained from the gridded distribution of Guenther et al. (2012), and $C_{CE} = 0.52$ is an adjustment factor ensuring that F equals $\epsilon$ un-

der standard conditions (Müller et al., 2008). LAI is the leaf area index, and the different $\gamma$ factors describe the emission activity and are related to photon density and leaf temperature (PT), leaf age (A), CO₂ concentration, and soil moisture (SM) (see Guenther et al., 2006). In this work, we neglect the effects of drought ($\gamma_{SM} = 1$), and we use the parameterization of Possell and Hewitt (2011) for $\gamma_{CO₂}$. The summation in Eq. (1) is performed over eight layers of the canopy, for which the MOHYCAN model calculates the direct and diffuse radiative fluxes and leaf temperature in each layer. For the above-canopy meteorological parameters, the ECMWF (European Centre for Medium-range Weather Forecasts) ERA5 reanalysis is used (Hersbach et al., 2020).

The MAGRITTE model accounts for diurnal variations in the chemical compounds through correction factors on the chemical reaction rates, photolysis rates, convective fluxes, and boundary layer mixing diffusivities computed via a full simulation of the diurnal cycle using a time step of 20 min (Stavrakou et al., 2009; Bauwens et al., 2016). Diurnal profile shapes calculated with the short time step are also used to estimate HCHO vertical columns at the overpass time of the TROPOMI instrument. The simulations are carried out with a time step of 24 h using daily averaged chemical and dynamical parameters adjusted for the diurnal cycle using the pre-calculated correction factors. This setup provides an optimal balance between model accuracy and computational cost.

## 3.2 Inversion setup

The inversion aims at optimizing the emissions in the model in order to improve the comparison between modeled HCHO columns and TROPOMI observations. The emissions are optimized per emission category (biogenic, pyrogenic, and anthropogenic), per grid cell, and per week. We only use TROPOMI data between May and September, when HCHO columns are highest, biogenic and pyrogenic emissions become important, and the photochemical lifetime of HCHO precursors is short, resulting in a minimal impact of dilution and transport of VOCs away from the emission areas. Even though the assimilation window only includes the summer months, the emissions can be modified throughout the full year due to temporal correlations imposed on the emission parameters (see below). We always consider 4 weeks per month, which are defined as sets of 8, 7, 8, and 7 d or 8, 8, 8, and 7 d for months with 30 and 31 d, respectively. The inversions are performed independently for each year. We can write the top-down emissions $G(x, t)$ as

$$G(x, t) = \sum_{j=1}^{3} \phi_{j}(x, t) \exp(f_{j}(x, t)), \quad (2)$$

where j denotes the emission category, $\phi_j$ denotes the a priori emission fluxes, $f_j$ denotes the emission parameters to be optimized, and x and t are the spatial coordinates (longitude and latitude) and time (week). For the biogenic com-
ponent, the emission parameters are only applied to isoprene and monoterpenes, as they make up the majority of the total biogenic emissions. Furthermore, isoprene and monoterpenes are short-lived and therefore have a strong impact on the HCHO distribution. The fluxes from a pixel and source category are only optimized when the a priori weekly value exceeds $10^{10}$ molec. cm$^{-2}$ s$^{-1}$ at least once throughout the year. The total number of parameters to be optimized by the inversion is of the order of $2 \times 10^5$ for each year.

The modeled HCHO columns are compared to the observed, bias-corrected HCHO columns of TROPOMI after application of the TROPOMI averaging kernels to the modeled profiles. We define the cost function $J$ by

$$J(f) = \frac{1}{2} \left[ (H(f) - y)^T \mathbf{E}^{-1} (H(f) - y) + f^T \mathbf{B}^{-1} f \right], \quad (3)$$

where $f = f_j(x, t)$ is a dimensionless vector consisting of the emission parameters of Eq. (2) to be optimized; $H(f)$ denotes the HCHO columns simulated by the model and smoothed by the TROPOMI averaging kernels; $y$ denotes the bias-corrected TROPOMI HCHO columns; and $\mathbf{E}$ and $\mathbf{B}$ are the error covariance matrices of the HCHO columns and emission parameters, respectively (Müller and Stavrakou, 2005). The first term in Eq. (3) quantifies the difference between the modeled and observed HCHO columns, whereas the second term acts as regularization, preventing the optimized emissions from excessively deviating from the a priori emissions. The error covariance matrix $\mathbf{E}$ is assumed to be diagonal (i.e., errors are uncorrelated), and its elements consist of the retrieval errors in the HCHO columns, to which a model error taken equal to $2 \times 10^{15}$ molec. cm$^{-2}$ is quadratically added. The matrix $\mathbf{B}$, on the other hand, has a diagonal component consisting of the square of the relative errors in the flux parameters $f$, which are set equal to 0.9 (corresponding to an uncertainty factor of 2.5) for biogenic and biomass burning emissions and to 0.7 (factor of 2 uncertainty) for the anthropogenic emission fluxes. The non-diagonal elements of $\mathbf{B}$ are dependent on spatio-temporal correlations of the errors in the fluxes. For biogenic and biomass burning emissions, the spatial correlations are assumed to decrease exponentially with distance $d_{ij}$ between two grid cells $i$ and $j$, such that

$$B_{ij} = \left(\frac{\sigma_i}{\phi_i}\right)^2 \exp\left(-\frac{d_{ij}}{l}\right), \quad (4)$$

where $l$ is the decorrelation length assumed to be 200 km, $\phi_i$ is the total flux emitted from grid cell $i$, and $\sigma_i/\phi_i$ is its relative error. Furthermore, we assume that biogenic emissions emitted by different plant functional types (PFTs) are uncorrelated, such that

$$B_{ij} = \sum_{n} e_{i}^{n} e_{j}^{n} \left(\frac{\sigma_{i}}{\phi_{i}}\right) \left(\frac{\sigma_{j}}{\phi_{j}}\right) \exp\left(-\frac{d_{ij}}{l}\right), \quad (5)$$

with $e_i^n$ being the fraction of flux emitted in cell $i$ by PFT $n$. Similarly, the anthropogenic emissions from different countries

or from different EDGARv4.3.2 categories (road transport, emission processes during production and application, fuel exploitation, energy for buildings, combustion for manufacturing, and the rest) are assumed uncorrelated. This choice is justified by the fact that anthropogenic emission policies are usually taken at the country level and that generally uniform methodologies and emission factors are applied within each country in bottom-up inventories. The spatial correlation coefficient is constant (value of 0.2) for a given category within the same country. No additional dependence on distance is assumed. The temporal correlation of the errors is assumed to decrease linearly from 0.7 for consecutive weeks to zero after 3 months for all emission categories. Finally, after averaging the data weekly, we filter out TROPOMI HCHO columns below $3 \times 10^{15}$ molec. cm$^{-2}$, since these low columns are most uncertain, as discussed in Sect. 4). This filtering of low columns results in too-high top-down emissions in northern Europe (see Sect. 5).

The inversion derives updated emission parameters $f$ that minimize the cost function in Eq. (3). To achieve this, we make use of the adjoint of the MAGRITTEv1.1 model, which calculates the gradient of the cost function with respect to the emission parameters. The cost function is iteratively minimized using a large-scale quasi-Newton method (Gilbert and Lemaréchal, 1989), and a new set of emission parameters $f$ is calculated. This procedure is repeated until the norm of the gradient of the cost function is decreased by a factor of 100 with respect to its initial value. This is achieved after about 40 to 50 iterations.

## 4 Validation of TROPOMI HCHO over Europe

Previous studies have shown that satellite HCHO data are biased low for high HCHO columns and high for low HCHO columns compared to aircraft (Zhu et al., 2020) and ground-based data (Vigouroux et al., 2020; De Smedt et al., 2021). A careful validation of the satellite product is therefore required in order to prevent the propagation of these biases into the derived top-down emissions. The nature of these biases is unknown at this point. To rule out a possible geographical dependence of the bias, we focus on stations in the European domain to derive the bias correction. We use the FTIR station data described in Sect. 2.2 and perform the validation using daily spatio-temporal collocation criteria, following the approach of Vigouroux et al. (2020). More specifically, we average TROPOMI scenes within 20 km of the station location. Additionally, we select ground-based observations acquired within 3 h of the satellite overpass (i.e., from 10:30 to 16:30 LT) in order to minimize the effect of diurnal variations. These criteria provide an optimal balance between data averaging and collocation (see Vigouroux et al., 2020; De Smedt et al., 2021). Days with no valid TROPOMI measurements or no valid ground-based measurement are excluded. For the validation, we only select ground-based and
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Average vertical sensitivity profiles for TROPOMI (blue), FTIR (orange), and MAX-DOAS (red) instruments. Both the averaging kernels (AVKs) and pressure levels are averaged using all measurements at the station locations.</figcaption></figure>

TROPOMI data from May to September, since we use data during this period as constraints for the inversion. Furthermore, we do not consider MAX-DOAS data for TROPOMI validation due to the difference in vertical sensitivities between TROPOMI and MAX-DOAS, shown in Fig. 3. The vertical sensitivities of the FTIR and TROPOMI observations have more overlap, specifically in the higher atmospheric layers, whereas the MAX-DOAS measurement sensitivity is maximum in the boundary layer. FTIR HCHO columns are therefore better suited for validation purposes, whereas MAX-DOAS data are particularly useful to provide complementary information to the satellite columns.

To make a meaningful comparison between space-based and ground-based data, it is necessary to apply a “smoothing” technique to the ground-based data (Rodgers and Connor, 2003). In this regard, we follow the methodology described in Sect. 4.2 of Vigouroux et al. (2020). First, the vertical profile of HCHO at each station is corrected for the different a priori profile that was used in the retrieval of the satellite column. Using the averaging kernel of the ground-based measurement, the a priori profile is substituted by the TM5-MP profile of the TROPOMI retrieval, thereby modifying the retrieved profile at altitudes where the measurement is not very sensitive. The modified station profile is then smoothed by applying the averaging kernel of the collocated clear-sky TROPOMI measurements. Finally, a correction to the total column is applied to account for topographical elevation variations in the TROPOMI footprint around the station. The smoothed columns can then be directly compared to the satellite columns. In the comparisons, we used weekly

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Scatterplot of space-based (TROPOMI) vertical HCHO columns versus ground-based (FTIR) vertical HCHO columns. The data have been collocated and weekly averaged. The linear fit is shown in blue, and its formal uncertainties are depicted by the blue-shaded area.</figcaption></figure>

averages, since weekly averaged TROPOMI HCHO columns are used as top-down constraints in the inversion. A comparison between TROPOMI HCHO columns and the smoothed FTIR columns is illustrated in the scatterplot of Fig. 4.

The linear regression is performed using the Passing-Bablok estimator (Passing and Bablok, 1983, 1984; Bablok et al., 1988), which is commonly used in the comparison between measurements of different methods by virtue of its symmetrical properties and its robustness in the presence of outliers. The comparison of the collocated data shows that TROPOMI HCHO columns are biased with respect to FTIR measurements. The FTIR data suggest that the high TROPOMI HCHO columns are too low, and the low TROPOMI columns are too high with respect to the FTIR measurements. This agrees well with previous FTIR-based bias corrections over the global domain (Vigouroux et al., 2020). We note that for low TROPOMI HCHO columns, i.e., below $3 \times 10^{15} \text{ molec. cm}^{-2}$, the regression line falls below a large majority of FTIR data, indicating that it is not suitable in that range.

The statistical uncertainty in the linear fit calculated and shown inset in Fig. 4 is clearly underestimated, given the substantial column differences between the different stations (see also monthly means for individual stations in Fig. S1 in the Supplement). We note that stations in northern Europe (Kiruna and Sodankylä) show no significant bias with respect to TROPOMI, while stations at mid-latitudes (Germany and France) show a larger bias. The linear scaling of
the bias with respect to the vertical column suggests a large bias in southern Europe. However, due to the lack of ground-based data available in southern Europe, we cannot confirm or refute this claim at this point. An accurate bias correction is particularly important for southern Europe since the highest biogenic emissions are released from the warmer southern European ecosystems.

An additional source of uncertainty originates in data averaging. Averaging the HCHO data inherently introduces a bias in the fit towards weeks with fewer data. We also find that averaging over longer periods (e.g., monthly) of the satellite data generally results in larger slopes and smaller intercepts in the fit compared to the weekly based validation. Furthermore, different assumptions in the retrieval of ground-based and satellite HCHO columns might also skew the dataset, since negative HCHO columns are absent in the FTIR dataset, whereas they are present in the TROPOMI dataset.

In our inversion scheme, we apply the bias correction based on FTIR data (Fig. 4), given by

$$[\text{TROPOMI}]_{\text{BC}} = 1.48 \times [\text{TROPOMI}] - 1.40 \times 10^{15} \text{molec.cm}^{-2}, \quad (6)$$

where [TROPOMI] and [TROPOMI]BC are the retrieved and bias-corrected (weekly averaged) HCHO total columns, respectively. This FTIR-based relation of Eq. (6) is in good agreement with previously published bias corrections of the TROPOMI HCHO product using FTIR (Vigouroux et al., 2020) and MAX-DOAS data (De Smedt et al., 2021). As illustrated in Fig. 5a and b, the FTIR-based correction brings small changes to the HCHO levels in northern Europe but induces substantial increases (30% to 40% over emission sources) in southern Europe.

# 5 Inversion results

## 5.1 Modeled HCHO columns

The a priori and a posteriori model columns are shown in Fig. 5c and d, respectively. The a priori column distribution is calculated based on a priori emission inventories as described in Sect. 3. The model computes the HCHO column at the time of the TROPOMI overpass. The modeled vertical profile is then convolved with the averaging kernel of the satellite observation, and a smoothed total column is derived. We note that this implies that model columns are calculated only whenever valid satellite measurements are available, and the weekly model averages are sampled as the observations. The modeled columns can therefore be directly compared to the TROPOMI HCHO columns of Fig. 5b. The a priori model provides relatively good agreement with the HCHO observations in terms of spatial distribution with, on average, higher columns in southern Europe (~10¹⁶ molec.cm⁻²) compared to the north (~5 × 10¹⁵ molec.cm⁻²). Some anthropogenic

sources appear to be overestimated in the model, such as the Ruhr area in Germany. Moreover, comparison of the a priori model with the bias-corrected data shows that the a priori columns in southern Europe are generally too low (see Fig. 5e).

The optimized HCHO columns inferred through inversion of the bias-corrected TROPOMI HCHO columns are shown in Fig. 5d. Because the bias-corrected columns are much higher than in the a priori model in southern Europe, the inversion leads to large column increases of about 2 to 3 × 10¹⁵ molec.cm⁻² in the Mediterranean region. This can also be seen in Fig. 5e and f, which illustrate the absolute difference between the bias-corrected observations (Fig. 5b) and either the a priori model or the a posteriori model. The a priori difference map shows the model underestimation in southern Europe clearly. This is largely corrected for by the optimization, which shows good agreement between model and observations, especially in mainland Europe. In northern Europe and over the seas, differences of the order of 1 × 10¹⁵ molec.cm⁻² are seen. Even higher discrepancies are found over the northern seas (> 65° N), which is caused by the large uncertainties in the satellite product and by the filtering out of low weekly columns (< 3 × 10¹⁵ molec.cm⁻²). This filter leads to a positive bias of the temporally averaged columns in regions where the columns are frequently lower than this threshold, as is the case over northern seas and Scandinavia. Furthermore, whereas the precision of weekly averaged TROPOMI columns amounts to only a few percent of the column and plays a negligible role in southern Europe, the high cloudiness and low columns found at high latitudes lead to larger random errors, sometimes of the same order as the columns, thereby aggravating the potential bias due to the threshold. The inversion results in these regions should therefore be considered with caution. We also note that oceanic data are excluded from the inversion because the columns over oceans are mostly due to the oxidation of methane (Stavrakou et al., 2009). Consequently, the comparison between model and observations does not change significantly over the northern seas. We note however a substantial reduction in the bias over the Mediterranean and Black seas, due to emission changes over land and subsequent transport of HCHO precursors (or of intermediates in their oxidation mechanisms) to those areas.

Figure 6 shows time series of observed and modeled HCHO columns in 2019 for four regions in the European domain: Germany, central Europe, Iberia, and Anatolia. The optimized model based on our standard inversion described in Sect. 3 uses weekly TROPOMI observations (orange curve). Figure 6 includes the results of an inversion using monthly TROPOMI HCHO columns in order to evaluate the effect of temporal averaging on the top-down emissions (see also scenario 10 in Sect. 6 and Table 3). The emission increments of the monthly inversion run are applied to the weekly a priori HCHO columns in order to allow for comparison with the standard weekly inversion results.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** HCHO columns over the European domain averaged from May to September for the years 2018 to 2021. Panels (a) and (b) show TROPOMI HCHO columns without bias correction and with the FTIR-based bias correction from Eq. (6), respectively. The latter is used as constraints in the inversion framework. Panels (c–d) contain modeled HCHO columns computed by the MAGRITTEv1.1 model. Panel (e) depicts the a priori HCHO columns, and (f) shows the columns based on the optimized emissions resulting from the top-down inversion. Finally, (g) is the absolute difference between the bias-corrected TROPOMI observations (b) and the a priori model (c), and (f) shows the absolute difference between (b) and the optimized model (d).</figcaption></figure>

Over Germany, the TROPOMI HCHO columns have a larger variability than the a priori model. The weekly optimization improves the model by increasing the HCHO columns during peaks and decreasing the model during troughs in the observed time series. The monthly inversion does not affect the weekly variability in the model; therefore the Pearson correlation coefficient of the monthly optimization (R = 0.84) is lower than the standard weekly optimization (R = 0.92). The second region of Fig. 6 focuses on central Europe, more precisely the Pannonian Basin. This

region was chosen for its biogeographical and climatological uniformity, as it is surrounded on all sides by mountain ranges. Similar to Germany, we find a large improvement in the correlation coefficient for the weekly top-down inversion (from 0.92 to 0.98), where the high and low HCHO columns in the time series are much better reproduced. For the regions Iberia and Anatolia, both in southern Europe, the results of the weekly and monthly optimizations are similar. The TROPOMI HCHO columns are much higher than the a priori model, which is partly due to the effect of the bias
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a) Regions used for averaging in the following panels. (b–e) Weekly time series of HCHO columns of TROPOMI observations (black), MAGRITTE a priori model (blue), top-down weekly inversion (orange), and top-down monthly inversion (purple). The HCHO columns cover May to September for the year 2019. In each panel, the Pearson correlation coefficient is indicated for the a priori model (blue) and weekly (orange) and monthly (purple) inversion results.</figcaption></figure>

correction. The top-down HCHO columns follow the high observed columns as expected. Due to the relatively constant meteorology in southern Europe (i.e., high temperatures and few clouds), the variability in HCHO columns is much lower. Consequently, the difference between the weekly and monthly inversion run is relatively small.

## 5.2 Top-down emissions

Figure 7a shows the a priori isoprene emissions of the MEGAN–MOHYCAN inventory, which were used as input for the top-down inversion. High-emission regions generally have a warm climate, along with vegetation consisting of strong isoprene emitters. The a priori emissions are heavily dependent on assumptions regarding the variable and uncer-

tain emission factors of the different vegetation types, as well as on the land-use map that enters the emission model. Figure 7b shows top-down isoprene emissions optimized based on bias-corrected TROPOMI HCHO columns. The emission enhancement is presented in Fig. 7c, calculated as the ratio of the top-down emissions over the a priori emissions. The increment map shows increases in the emissions across the majority of the domain, except for some regions between 50 and 60° N. Especially in southern Europe, the data suggest higher isoprene emissions in order to explain the high observed TROPOMI HCHO columns. Similarly, at high latitudes there are increased top-down isoprene emissions, but this is likely due to the positive bias in this region introduced by the filtering of low TROPOMI HCHO columns.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Results of the inversion of isoprene emissions over Europe. Panel (a) shows isoprene emissions from the MEGAN–MOHYCAN inventory used as input for the inversion. Panel (b) shows the optimized emissions resulting from our top-down inversion. Panel (c) shows the enhancement map derived as the ratio of the optimized and the a priori emissions (b divided by a). The data are averaged over the summer months (May to September) from 2018 to 2021.</figcaption></figure>

Table 2 presents the annual emissions in the European domain. The inversion strongly increases isoprene emissions, by more than a factor of 2 with respect to the a priori emissions (from 8.1 Tg yr⁻¹ in MEGAN–MOHYCAN to 18.5 Tg yr⁻¹ in the top-down emissions). In a previous inversion study by Bauwens et al. (2016) using OMI data, only a mild increase in isoprene emissions was derived over Europe (from 6.8 to 8.4 Tg yr⁻¹, for a slightly different domain). This is mainly due to the absence of a bias correction

<table><caption>**Table 2.** Average annual emissions in Tg yr⁻¹ of the three categories of VOC emissions in the a priori and optimized models for the years 2018 to 2021.</caption><thead><tr><td>Emission category</td><td>A priori</td><td>Optimized</td></tr></thead><tbody><tr><td>Isoprene</td><td>8.1</td><td>18.5</td></tr><tr><td>Biomass burning VOCs</td><td>2.3</td><td>2.6</td></tr><tr><td>Anthropogenic VOCs</td><td>16.9</td><td>14.0</td></tr></tbody></table>

in their inversion study. The results of Bauwens et al. (2016) suggested top-down isoprene emissions across the domain of 8.4 Tg yr⁻¹ over the years 2005 to 2013, i.e., 23% higher than their a priori biogenic emissions. The majority of the biogenic emissions take place in southern Europe, where the high observed HCHO columns lead to a strong increase in the top-down emissions in our inversion. The time series of isoprene emissions in Europe (Fig. 8a) shows that the biogenic emission increases mainly occur during the summer peak from June to August with up to a factor of 3 higher emissions. As expected, the emissions are close to zero during winter. We note that the emissions are optimized using TROPOMI observations only during summer months, but due to temporal correlations imposed in the inversion, the emissions during winter months are also affected.

Besides isoprene, monoterpenes are also impacted by the inversion. A large increase in their emissions is derived, from 3.9 to 6.9 Tg yr⁻¹. The enhancement of the total monoterpenene emissions is smaller than for isoprene because a large fraction of those emissions originate from boreal forests, where the emission increments are much lower (see also Fig. S7). Methanol, another prominent BVOC with a total source of ca. 6.3 Tg yr⁻¹ over the model domain in MEGAN–MOHYCAN, is not optimized by the inversion. Due to its long lifetime of several days, it mainly contributes to the background HCHO column and is therefore less suitable for top-down inversion studies using HCHO column data.

The biomass burning VOC emissions of the top-down model are on average higher than the QFED a priori inventory. The most prominent source of biomass burning VOCs in Europe is agricultural waste burning. This is prominent in the time series in Fig. 8b, which shows peaks around April and September mainly due to agricultural fires in Ukraine and Russia (see Fig. S3). Figure 8b also shows substantial year-to-year variability, in which the top-down inversion does not always lead to changes in the emissions with respect to the QFED a priori inventory. However, the data on average lead to higher top-down fire emissions which are mainly caused by increased agricultural fire VOC emissions in eastern Europe and the Danube river valley, as can be seen in Fig. S3c.

In terms of anthropogenic emissions, the optimized fluxes are lower than the a priori fluxes from CAMS-GLOB-ANT by about 20%. This decrease is driven by generally lower
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Interannual time series of monthly biogenic emissions (a), biomass burning VOC emissions (b), and anthropogenic VOC emissions (c) across the European domain in teragrams. The a priori values are denoted by the solid lines, and the optimized values are shown by the dotted lines.</figcaption></figure>

observed HCHO columns over cities and industrial regions. Furthermore, the TROPOMI HCHO columns are lower than that of the a priori model during colder weeks (see Fig. 6), which suggests a larger relative contribution of biogenic emissions to the total column. Figure S4 shows that the optimized emissions are lower than CAMS-GLOB-ANT for most countries at mid-latitudes (between roughly 45 and 55°N), with the largest differences occurring in Germany. The effect of the imposed country-dependent spatial correlations of the emissions parameters is clearly visible, as the emission increments of anthropogenic VOCs are similar within country borders. The time series in Fig. 8c depicts a weak seasonal cycle with higher emissions in winter compared to summer. This seasonal cycle is preserved in the top-down fluxes, but with lower values during the whole year.

## 5.3 Model evaluation

We evaluate the model against ground-based remote sensing and in situ measurements. For the FTIR and MAX-DOAS comparisons, we select model data at the time of the measurement at the 0.5° × 0.5° pixel in which the station is located. The model profile is smoothed using the averaging kernel of the measurement. Comparison of the optimized model to FTIR station data shows an improvement with respect to the a priori model at Bremen, Garmisch, Karlsruhe, and

Paris, as shown by the reduction in average absolute deviations (AADs) at those stations (see Fig. 9). The opposite effect occurs at the northern stations (Kiruna, Sodankylä, and St. Petersburg), where the inversion worsens the agreement between model and data. More specifically, model overestimations are found at high-latitude Scandinavian sites, invalidating the emission increases in northern Scandinavia. At St. Petersburg, the TROPOMI HCHO columns are generally very low, leading to a decrease in the modeled HCHO columns after the inversion, yet FTIR data suggest higher HCHO columns. The reduction in emissions that follows from the inversion in this region (around 60°N; see Fig. 7) is therefore not corroborated by the FTIR data at St. Petersburg.

The comparison of the model columns with the MAX-DOAS measurements of stations listed in Table 1 is shown in Fig. 10. For all stations (except Athens), the optimized HCHO columns are not improved with respect to the a priori column, and the average absolute difference is increased at several locations. The MAX-DOAS data show a worsening of the optimized model over Germany and the Benelux, whereas the FTIR data suggest an improvement of the model over Germany and France. The MAX-DOAS stations (except Athens) are located in regions where the a priori HCHO model columns are on average larger than the observations, which leads to a decrease in VOC emissions and HCHO columns after optimization. On the other hand, the a priori
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Time series of monthly averages of observed FTIR and a priori and a posteriori modeled HCHO columns from May to September. The modeled columns are calculated using either a priori or optimized emissions. The model columns are smoothed using the averaging kernels of the collocated measurements. The average absolute deviation (AAD) between the observations and model, calculated as the mean of $|\text{HCHO}_{\text{obs}} - \text{HCHO}_{\text{model}}|$ in units of $10^{15} \text{ molec. cm}^{-2}$, is given at the top right of each panel.</figcaption></figure>

ori model underestimates the MAX-DOAS data (Fig. 10). This discrepancy suggests a potential issue with the vertical profile of HCHO in the MAGRITTE model, i.e., too large background concentrations at higher altitudes, to which TROPOMI retrievals are much more sensitive than MAX-DOAS (Fig. 3). This overestimation causes the decrease in anthropogenic VOC emissions over Germany and the Benelux (see Fig. S4), leading to the discrepancy with MAX-DOAS observations. On the other hand, the higher emissions of the optimized model in southern Europe are corroborated by the Athens MAX-DOAS data. This station is located at an

altitude of 500 m, which is higher than the city itself. Since MAX-DOAS observations only consider upward-looking elevation angles, an important part of the HCHO column is not measured. We take this into account by only including part of the model column at 500 m altitude and higher for the comparison. The increased biogenic VOC emissions around Athens lead to increased HCHO mixing ratios of the model above 500 m altitude, which agrees with the MAX-DOAS data at Athens.

As an additional evaluation, we compare the modeled near-surface HCHO concentrations to in situ HCHO data
<figure><img src="image_13.png" /><figcaption>Figure 10. Same as Fig. 9, but with MAX-DOAS HCHO columns.</figcaption></figure>

in Fig. 11 and Table S1. We sample the model at the pixel and vertical layer corresponding to the coordinates and altitude (above sea level) of the station and compute the ratio of monthly averaged daytime (usually 8–12 h) modeled concentration to the corresponding monthly averaged collocated in situ measurements. The ratio is computed using a climatological average if all observations are dated before May 2018 (see Table 1). We find that for the two Spanish stations (San Pablo de los Montes and Campisábalos), the a priori model surface concentration is already higher than the measurements. This discrepancy is increased even further after the optimization. As for the other southern European stations, we find improvements in the south of France (La Tardière

and Peyrusse Vieille) and minor improvements for the Italian stations (Ispra and Montelibretti), where the a priori model is about 25% too low, but the optimized model is 10% to 20% higher than the observations. For the stations located at higher latitudes, we generally see a small effect of the inversion on the optimized model, with the a priori model performing better for the Waldhof and Brotjacklriegel stations. We note that some caution is needed for comparisons that are not temporally collocated (all except LAH, PAB, PEY, and TAR; see Table 1), since meteorology has an important effect on HCHO concentrations, and for the stations located in mountainous terrain (PAB, CAM, BRO, DON, and SCH),
due to local orography-induced transport patterns impacting
the concentrations.

To summarize, the optimized model provides better agree-
ment with the FTIR stations below about 55° N and with the
only southern MAX-DOAS station, Athens. The other MAX-
DOAS stations located in Germany and the Benelux display
higher HCHO columns than the a priori model, yet the (bias-
corrected) TROPOMI HCHO columns do not suggest an in-
crease in the emissions in this area. The three FTIR stations
in northern Europe compare better with the a priori than with
the optimized model, thereby contradicting the emission in-
creases in northern Scandinavia and the emission decreases
in the St. Petersburg area. In southern Europe, we see large
increases for biogenic emissions of the optimized model.
This is corroborated by HCHO column measurements at the
Athens station, although in situ measurements provide mixed
results in Spain, France, and Italy.

6 Sensitivity inversions

The results of the top-down inversion are subject to uncer-
tainties related to the model, the inversion setup, the a pri-
ori inventories, and the observations. These uncertainties are
difficult to quantify. In this section, we present a range of 11
sensitivity inversions aimed at evaluating the effect of dif-
ferent sources of uncertainty on the derived top-down emis-
sions. We provide an overview of the different sensitivity
runs (S1 to S11) in Table 3. Each sensitivity run was per-
formed for the year 2019, and their annual emissions are
shown in Fig. 12.

The first scenario (S1) described in Table 3 is an inver-
sion constrained by TROPOMI HCHO columns without bias
correction. The data used to constrain the inversion in this
scenario are depicted in Fig. 5a, as opposed to Fig. 5b for the
standard inversion scenario (S0). Because the bias correction
leads to much higher HCHO columns over strongly emit-
ting regions (~40% increase), the S1 run has much lower
emissions across Europe for all emission categories. Since
the effect of the bias correction is strongest in southern Eu-
rope, the biogenic emissions are most affected, with the S1
scenario having 3 times lower isoprene emissions compared
to the standard run (6.7 Tg yr⁻¹ versus 18.3 Tg yr⁻¹ for S0).
The S1 results highlight the strong impact of TROPOMI ob-
servations on the top-down emissions, and therefore the im-
portance of a comprehensive and homogeneous validation
dataset.

Next, we evaluate the effect of different a priori in-
ventories on the inversion results. In scenario S2, we
substitute MEGAN-MOHYCAN biogenic emissions with
the CAMS-GLOB-BIOv3.1 inventory (Sindelarova et al.,
2022). Figure 13a shows the isoprene emissions of
CAMS-GLOB-BIOv3.1 gridded to 0.5° × 0.5° and aver-
aged from May to September in 2019. On average, the
CAMS-GLOB-BIOv3.1 isoprene emissions are much lower

than MEGAN-MOHYCAN, with annual total emissions
of 4.1 Tg for CAMS-GLOB-BIOv3.1 versus 7.9 Tg for
MEGAN-MOHYCAN. Furthermore, the distribution of
emissions is very different, especially in Turkey, which is
a major source region of biogenic emissions in MEGAN-
MOHYCAN. Those differences in a priori emissions have
a major effect on the optimized top-down emissions, which
are shown in Fig. 13b. The optimized total isoprene emis-
sions are more than 30% lower in the S2 run, while the ef-
fect on the optimized fire and anthropogenic VOC emissions
is limited. Although the top-down biogenic emissions of the
S2 run are lower than in the standard S0 run, we still see a
large increase with respect to the a priori emissions of about
300%. Similar to the standard case, the emissions are mainly
increased in southern Europe and most heavily in Turkey by
over a factor of 10 (see Fig. 13c).

In scenario S3, the a priori biomass burning emissions are substituted with the GFED4s inventory (van der Werf et al., 2017). The GFED4s VOC emissions are based on a burned area approach and are on average lower than the QFED2.4 emissions used in our standard inversion, based on a fire radiative power approach (Darmenov and da Silva, 2015). However, there is good agreement between the two inventories in terms of temporal and spatial occurrence of fire events, as they both make use of MODIS satellite data (Pan et al., 2020). The lower GFED4s emissions of scenario S3 result in lower top-down fire emissions, totaling 2.1 Tg VOC for the S3 run as opposed to 3.5 Tg VOC for the standard S0 run. Since biomass burning is not a dominant emission source of VOCs in Europe, the effect of differences in a priori emissions has only a limited impact on the top-down biogenic and anthropogenic VOC emissions.

The S4 scenario tests an inversion using biomass burning emission factors from Akagi et al. (2011), instead of emission factors from Andreae (2019) used in the standard inversion. Using the Akagi et al. (2011) emission factors, the total emission of pyrogenic VOCs is decreased by about 20% in the a priori and optimized emissions. This can be mainly attributed to differences in emissions for agricultural waste burning. Specifically, the emission factors of large alkanes (butane or higher) are substantially lower in the S4 run. As for the S3 scenario, the lower a priori emissions lead to lower a posteriori pyrogenic emissions. The impact of this sensitivity run on the optimized biogenic and anthropogenic emissions is small.

Model assumptions represent an important source of un-
certainty in the inversions. In the next scenarios, we evaluate
the effects of some of these assumptions on the top-down
results. The S5 and S6 runs evaluate the effect of the as-
sumed errors in the emission parameters of the a priori inven-
tories by decreasing and increasing the errors, respectively.
As described in Sect. 3, because bottom-up inventories have
no quantified uncertainties, we assume relative errors in the
flux parameters corresponding to a factor of 2.5 for biogenic
and biomass burning emissions and a factor of 2 for anthro-
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Comparison between in situ measurements and near-surface concentrations from the a priori simulation (a) and from the optimized model (b). The data are monthly averaged, and only data from May to September are included in the ratio calculation. The color of the station symbol represents the median ratio of the model concentration and the station concentration. The ratios calculated with the a priori and optimized model are given in Table S1. Stations with data after 2018 use collocated model concentrations for the comparison. For other stations with only older data, climatological averages are used (see Table 1).</figcaption></figure>

<table><caption>**Table 3.** Sensitivity inversions conducted in this study and annual VOC emissions in Tg yr⁻¹ over the model domain. AVOC stands for anthropogenic VOC.</caption><thead><tr><th rowspan="2">Scenario</th><th rowspan="2">Description</th><th colspan="2">Isoprene (Tg)</th><th colspan="2">Fire VOC (Tg)</th><th colspan="2">AVOC (Tg)</th></tr><tr><th>A priori</th><th>Optim.</th><th>A priori</th><th>Optimized</th><th>A priori</th><th>Optimized</th></tr></thead><tbody><tr><td>S0</td><td>Standard inversion run described in Sect. 3</td><td>7.9</td><td>18.3</td><td>2.7</td><td>3.5</td><td>16.9</td><td>13.4</td></tr><tr><td>S1</td><td>As S0, but using TROPOMI HCHO columns without bias correction</td><td>7.9</td><td>6.7</td><td>2.7</td><td>2.0</td><td>16.9</td><td>11.7</td></tr><tr><td>S2</td><td>As S0, but with CAMS-GLOB-BIOv3.1 (Sindelarova et al., 2022) as a priori biogenic emission inventory</td><td>4.1</td><td>12.6</td><td>2.7</td><td>3.6</td><td>16.9</td><td>13.9</td></tr><tr><td>S3</td><td>As S0, but with GFED4s (van der Werf et al., 2017) as a priori biomass burning emission inventory</td><td>7.9</td><td>18.6</td><td>1.3</td><td>2.2</td><td>16.9</td><td>13.1</td></tr><tr><td>S4</td><td>As S0, but with biomass burning emission factors from Akagi et al. (2011)</td><td>7.9</td><td>17.8</td><td>2.0</td><td>2.4</td><td>16.9</td><td>13.1</td></tr><tr><td>S5</td><td>As S0, but with lower errors in the emission parameters (1.5 times lower)</td><td>7.9</td><td>17.0</td><td>2.7</td><td>3.3</td><td>16.9</td><td>14.5</td></tr><tr><td>S6</td><td>As S0, but with higher errors in the emission parameters (1.5 times higher)</td><td>7.9</td><td>18.6</td><td>2.7</td><td>3.4</td><td>16.9</td><td>11.5</td></tr><tr><td>S7</td><td>As S0, but with halved decorrelation lengths for a priori errors in biogenic and fire emissions (100 km instead of 200 km)</td><td>7.9</td><td>16.8</td><td>2.7</td><td>3.4</td><td>16.9</td><td>13.7</td></tr><tr><td>S8</td><td>As S0, but with doubled decorrelation lengths for a priori errors in biogenic and fire emissions (400 km instead of 200 km)</td><td>7.9</td><td>18.6</td><td>2.7</td><td>2.7</td><td>16.9</td><td>13.3</td></tr><tr><td>S9</td><td>As S0, but with lower isomerization rates of isoprene peroxy radicals, consistent with Wennberg et al. (2018)</td><td>7.9</td><td>17.2</td><td>2.7</td><td>3.3</td><td>16.9</td><td>13.1</td></tr><tr><td>S10</td><td>As S0, but using monthly instead of weekly averaged columns and optimizing monthly instead of weekly emission parameters</td><td>7.9</td><td>17.1</td><td>2.7</td><td>3.8</td><td>16.9</td><td>14.9</td></tr><tr><td>S11</td><td>As S0, but using TROPOMI HCHO data with less strict cloud filtering (&lt; 40%)</td><td>7.9</td><td>18.3</td><td>2.7</td><td>3.2</td><td>16.9</td><td>13.3</td></tr></tbody></table>
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Bar chart of annual emissions of the a priori inventory, the standard inversion run (S0), and the sensitivity inversions (S1 to S11) for biogenic emissions (a), biomass burning emissions (b), and anthropogenic VOC emissions (c). The emissions are calculated for the year 2019.</figcaption></figure>

pogenic emissions. By decreasing (increasing) these errors, the inversion will be less (more) inclined to depart from the a priori values in order to match the observations. Indeed, the S6 run, in which the errors are increased by a factor of 1.5, differs more strongly from the a priori inventories compared to the standard S0 run. Similarly, the S5 run with 1.5 times lower errors leads to smaller increases in biogenic emissions and smaller decreases in anthropogenic emissions. The effect of changing the error parameters is strongest for the anthropogenic emissions (~20%), though overall the effect of changing the errors is limited (less than 10% for biogenic emissions).

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Same as Fig. 7, but for scenario S2 using the CAMS-GLOB-BIOv3.1 inventory as a priori biogenic emissions (a). Panel (b) shows the optimized emissions of scenario S2, and (c) shows the ratio of the optimized and the a priori emissions (b divided by a). The data are averaged over the summer months (May to September) in 2019.</figcaption></figure>

Additionally, the correlations between errors in biogenic and pyrogenic emissions are assumed to decay exponentially with distance between the grid cells, for which the decorrelation length is taken to be 200 km (see Eqs. 4 and 5). The effect of the decorrelation length is equivalent to smoothing the emissions by introducing a spatial correlation of the emission parameters. In sensitivity inversion scenarios S7 and S8, we test the effect of a shorter (100 km) and longer (400 km) decorrelation length on the inversion results. Increasing the decorrelation length leads to overall increased biogenic emissions, since the increased emission parameters in southern Europe are forced over longer distances. Conversely, the S7
run with half the decorrelation length leads to about 10% lower biogenic emissions. As expected, the anthropogenic emissions are less affected by these changes.

Besides the inversion setup, the chemical mechanism is an additional source of uncertainty. In particular, the oxidation of isoprene by OH radicals is a complex process, and the initial oxidation steps are still an active research topic (Müller et al., 2019; Møller et al., 2019; Novelli et al., 2020; Schwantes et al., 2020; Medeiros et al., 2022). For example, the isomerization rates of specific isoprene peroxy radicals remain uncertain. In the standard run, we apply 1,6 H-shift rates based on the LIM1 mechanism (Peeters et al., 2014), in agreement with a recent laboratory determination (Novelli et al., 2020). These rates are higher than those given in Müller et al. (2019), which were derived by Wennberg et al. (2018) based on experimental work at Caltech (Teng et al., 2017). The lower isomerization rate leads to a higher HCHO yield from the oxidation of isoprene, which therefore leads to higher HCHO columns. The reason lies in the lower HCHO yields from the further oxidation of isomerization products like the hydroperoxy aldehydes (Peeters et al., 2014; Müller et al., 2019). The optimization of this scenario (S9) results in total isoprene emissions of 17.2 Tg yr⁻¹, i.e., 6% lower than in the standard run (18.3 Tg yr⁻¹). This reduction is largely compensated by higher pyrogenic and anthropogenic emissions, such that the total VOC emission of S9 is only slightly (2%) lower than in S0.

In all previous scenarios (S0 to S9), weekly averaged TROPOMI HCHO columns were used to constrain the emissions, and we update the emissions at weekly increments. In the S10 run, we evaluate the impact of temporal averaging by using monthly averaged HCHO columns and monthly emission parameters in the inversion. In Sect. 5.1, we showed that the monthly inversion is less capable of matching the weekly variability in the HCHO time series (see Fig. 6). Furthermore, the emissions of the S10 run are closer to the a priori model. For example, the annual anthropogenic VOC emissions of the S10 run are less decreased (–12%) with respect to the a priori inventory compared to the weekly S0 run (–21%), whereas the biogenic emission increase shows the opposite effect. This can be expected since the monthly inversion has about 4 times fewer observational constraints than the weekly inversion, which results in smaller departures from the a priori. The number of emission parameters is also lower in the monthly inversion, but this effect is mitigated due to temporal correlations between emission parameters. Additionally, the monthly inversion is less capable of disentangling the biogenic and anthropogenic emission sources, since it averages out meteorological variability at weekly timescales, which is a major driver of biogenic emission variability. The Benelux region provides an interesting example, since the weekly inversion displays increased biogenic emissions, whereas the monthly inversion leads to decreased biogenic emissions. The strong variability in the HCHO columns in the Benelux is displayed in Fig. 2. In order

to explain the large fluctuations in the observations, biogenic emissions are increased by the inversion during warm weeks, and anthropogenic emissions are decreased in order to account for the low columns during cold weeks.

Over northern Europe, the results of the monthly inversion are very different from the standard weekly inversion. In particular, the S10 optimization does not display the increase in the emissions above 60° N (see Figs. S5 and S6). As previously discussed, the increased biogenic emissions in Scandinavia in the standard run are due to the filtering out of low HCHO columns (<3 × 10¹⁵ molec.cm⁻²), implying that only high columns (i.e., positive outliers) remain in the dataset. However, the monthly averaging of the S10 run reduces the noise in the TROPOMI HCHO columns, leading to less filtering of low values and therefore to lower average HCHO columns at high latitudes. Consequently, the monthly inversion does not show an increase in emissions in Scandinavia and provides more reliable results in this region. In conclusion, the weekly inversion provides generally improved results compared to the monthly inversion, except at northern latitudes (above about 60° N).

Another important source of uncertainty in the emission inversion is the effect of clouds. To address this, the standard inversion includes a strict cloud filter which removes scenes with a cloud fraction larger than 0.2. However, this might imply that the emission inversion is biased towards cloud-free conditions. To evaluate this effect, we include a sensitivity run constrained by TROPOMI HCHO columns without the strict cloud filter (run S11). In the unfiltered TROPOMI dataset, the condition for QA (above 0.5) is still applied, which removes scenes with cloud fractions larger than 0.4. Because the observational dataset is modified for this run, we have recalculated the bias correction based on FTIR data, with the same methodology as that used in Sect. 4. The bias correction for this dataset is adjusted to [TROPOMI*]$_{BC}$ = 1.40 × [TROPOMI*] – 0.63 × 10¹⁵ molec.cm⁻², where [TROPOMI*] denotes the TROPOMI HCHO columns without a strict cloud filter.

The S11 sensitivity run leads to very small differences in total emissions compared to the standard run. Even though the cloud filter has a large impact on the TROPOMI HCHO columns, recalibrating the bias correction results in minor changes in the HCHO columns in southern Europe. Consequently, the emission patterns in southern Europe, which dominate the biogenic emissions over the domain, do not change much with respect to the standard run. In general, the bias-corrected [TROPOMI*] HCHO columns are lower over land masses and higher over seas with respect to the filtered dataset, with the largest differences occurring at midlatitude regions. This leads to on average slightly lower biogenic and anthropogenic VOC emissions in France, Belgium, the Netherlands, and Germany for the S11 run. However, as can be seen in Fig. 12, the total emissions do not change substantially, and overall the effect of the cloud filter on the inversion results is small.
To summarize, by far the largest uncertainty in our top-down emissions arises from the application of the bias correction. This emphasizes the importance of continuous and high-quality validation of the HCHO satellite product. This is especially important in southern Europe, where the bias correction has the biggest impact and where validation data are scarce. The biogenic a priori inventory also plays an important role in the inversion, mainly due to the large differences among the available bottom-up inventories. The range of different inversion settings and model parameters that we have tested results in typical differences of roughly 10% to 20% in the final total emission fluxes, as can be seen in Fig. 12.

## 7 Conclusions

We have performed a comprehensive top-down inversion of VOC emissions over Europe from 2018 to 2021 constrained by TROPOMI HCHO columns. We have used the MAGRIT-TEv1.1 chemistry transport model and its adjoint to optimize biogenic, pyrogenic, and anthropogenic VOC emissions at a weekly temporal resolution. The inversion suggests a strong underestimation of MEGAN emissions, especially in southern Europe, with the total annual isoprene flux increasing from 8.1 to 18.5 Tg yr⁻¹. This large increase is driven by the high (bias-corrected) HCHO columns observed by TROPOMI during the summer months in southern Europe. The biomass burning VOC emissions show a slight increase (from 2.3 to 2.6 Tg yr⁻¹) in the inversion, with a strong interannual variability in the emission increments. The anthropogenic VOC emissions are decreased by the inversion (from 16.9 to 14.0 Tg yr⁻¹), with the main reduction occurring over Germany.

Comparison of the TROPOMI HCHO columns and ground-based FTIR HCHO columns shows that TROPOMI HCHO columns are generally too low, especially for high columns. We derived a linear relationship characterizing the bias against FTIR data at European stations (excluding mountain sites), which agrees with the bias correction reported in previous validation studies using global FTIR and MAX-DOAS networks (Vigouroux et al., 2020; De Smedt et al., 2021). Additionally, we used FTIR, MAX-DOAS, and in situ data to validate the top-down model results and found that the optimized model performs well for FTIR stations below about 55° N. In northern Europe, where the signal-to-noise ratio of TROPOMI data is usually low, our top-down emissions are too high due to the filtering of low HCHO columns in the inversion. At MAX-DOAS stations, the optimized model shows no improvement, or even a worsening of the comparison with ground-based measurements, except at Athens. The MAX-DOAS data suggest that the emissions are too low in the Benelux and Germany, specifically during colder and cloudy days. Finally, the comparison with in situ data shows improved results at several locations (south of France, Italy) but too-high concentrations of the model

at other places, including the Spanish stations. These mixed results might reflect the limited representativeness of point measurements and/or difficulties with the model, e.g., regarding boundary layer mixing.

To evaluate the robustness of the inversion results, we have performed a range of sensitivity inversions. As expected, the top-down emissions are very sensitive to the observed HCHO columns and therefore also to the bias correction of the TROPOMI data. This highlights the importance of comprehensive and extended validation studies, using network remote sensing data and/or airborne validation campaigns. We note a scarcity of validation data in southern Europe, although at the time of writing, there are plans for FTIR measurement sites in Italy and Cyprus. Additionally, we find that the choice of a priori biogenic emissions plays an important role due to large differences between available bottom-up emission inventories. Other model-related effects on the inversion induce differences of up to 10% in the final annual total emissions.

Finally, the inversion using weekly averages was found to provide better results than the monthly inversion run, except in the northernmost part of the continent. Especially regions with variable meteorology are better treated by the weekly inversion due to the strong fluctuations in the HCHO columns. As the HCHO data quality continues to improve with the advent of more advanced satellite missions, the observational constraints on the inversion become more tight, allowing for higher spatio-temporal resolution of the inversion studies. In this European inversion study, we found that weekly averages and a spatial resolution of 0.5° are convenient in terms of signal-to-noise ratio and data availability. However, it may be possible to perform higher-resolution inversions in other regions in the world characterized by strong emissions (hence a strong HCHO signal) and moderate or low cloudiness.

**Data availability.** The TROPOMI HCHO dataset is a Copernicus operational product and is available at https://doi.org/10.5270/S5Ptjlxfd2 (Copernicus Sentinel-5P, 2018). The FTIR and MAX-DOAS datasets can be requested by contacting the PIs of each station. In situ data from the EMEP network are available at the EBAS data portal (https://ebas-nilu.no, EBAS, 2024). The top-down biogenic and biomass burning emission datasets generated in this study are available at the SEEDS project data portal (https://www.seedsproject.eu/data, SEEDS, 2024) and at the BIRA-IASB emission portal (https://doi.org/10.18758/71021088, Oomen et al., 2024). The MEGAN-MOHYCAN bottom-up inventory is also available at https://emissions.aeronomie.be (MEGAN-MOHYCAN, 2024). The FTIR data were made available as part of the PRODEX TROVA-E2 Belspo project (2021–2023).

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/acp-24-449-2024-supplement.
**Author contributions.** GMO carried out the analysis and wrote the manuscript. JFM and TS designed the inversion scheme and performed the optimizations. IDS described the TROPOMI HCHO data in Sect. 2.1. TB, RK, MM, MP, ARö, YT, and CV provided the FTIR measurements. CV provided the description of the FTIR datasets in Sect. 2.2. MMF, UF, FH, AM, AP, ARi, MVR, and TW provided the FRM₄DOAS measurements. AM provided the description of the FRM₄DOAS data (Sect. 2.3). All authors read and commented on the manuscript.

**Competing interests.** At least one of the (co-)authors is a member of the editorial board of *Atmospheric Chemistry and Physics*. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

**Acknowledgements.** We thank the two anonymous referees for their constructive remarks that helped improve this paper.

**Financial support.** This research was performed as part of the SEEDS project funded by the European Commission under the H2020 program (grant no. 101004318, 2021–2023). The Paris TCCON site has received funding from Sorbonne Université, the French research center CNRS, the French space agency CNES, and Région Île-de-France. This study also received funding from the FRM₄DOAS(-2.0) project under ESA contract numbers 4000118181/16/I-EF and 4000135355/21/I-DT-Ir.

**Review statement.** This paper was edited by Qiang Zhang and reviewed by two anonymous referees.

## References

Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072, https://doi.org/10.5194/acp-11-4039-2011, 2011.

Andreae, M. O.: Emission of trace gases and aerosols from biomass burning – an updated assessment, Atmos. Chem. Phys., 19, 8523–8546, https://doi.org/10.5194/acp-19-8523-2019, 2019.

Arneth, A., Schurgers, G., Lathiere, J., Duhl, T., Beerling, D. J., Hewitt, C. N., Martin, M., and Guenther, A.: Global terrestrial isoprene emission models: sensitivity to variability in climate and vegetation, Atmos. Chem. Phys., 11, 8037–8052, https://doi.org/10.5194/acp-11-8037-2011, 2011.

Bablok, W., Passing, H., Bender, R., and Schneider, B.: A General Regression Procedure for Method Transformation. Application of Linear Regression Procedures for Method Comparison Studies in Clinical Chemistry, Part III, J. Clin. Chem. Clin. Bio., 26, 783–790, https://doi.org/10.1515/cclm.1988.26.11.783, 1988.

Bates, K. H. and Jacob, D. J.: A new model mechanism for atmospheric oxidation of isoprene: global effects on oxidants, nitrogen oxides, organic products, and secondary organic aerosol, Atmos. Chem. Phys., 19, 9613–9640, https://doi.org/10.5194/acp-19-9613-2019, 2019.

Bauwens, M., Stavrakou, T., Müller, J.-F., De Smedt, I., Van Roozendael, M., van der Werf, G. R., Wiedinmyer, C., Kaiser, J. W., Sindelarova, K., and Guenther, A.: Nine years of global hydrocarbon emissions based on source inversion of OMI formaldehyde observations, Atmos. Chem. Phys., 16, 10133–10158, https://doi.org/10.5194/acp-16-10133-2016, 2016.

Beirle, S., Dörner, S., Donner, S., Remmers, J., Wang, Y., and Wagner, T.: The Mainz profile algorithm (MAPA), Atmos. Meas. Tech., 12, 1785–1806, https://doi.org/10.5194/amt-12-1785-2019, 2019.

Berndt, T., Hyttinen, N., Herrmann, H., and Hansel, A.: First oxidation products from the reaction of hydroxyl radicals with isoprene for pristine environmental conditions, Commun. Chem., 2, 21, https://doi.org/10.1038/s42004-019-0120-9, 2019.

Bryant, D. J., Dixon, W. J., Hopkins, J. R., Dunmore, R. E., Pereira, K. L., Shaw, M., Squires, F. A., Bannan, T. J., Mehra, A., Worrall, S. D., Bacak, A., Coe, H., Percival, C. J., Whalley, L. K., Heard, D. E., Slater, E. J., Ouyang, B., Cui, T., Surratt, J. D., Liu, D., Shi, Z., Harrison, R., Sun, Y., Xu, W., Lewis, A. C., Lee, J. D., Rickard, A. R., and Hamilton, J. F.: Strong anthropogenic control of secondary organic aerosol formation from isoprene in Beijing, Atmos. Chem. Phys., 20, 7531–7552, https://doi.org/10.5194/acp-20-7531-2020, 2020.

Churkina, G., Kuik, F., Bonn, B., Lauer, A., Grote, R., Tomiak, K., and Butler, T. M.: Effect of VOC Emissions from Vegetation on Air Quality in Berlin during a Heatwave, Environ. Sci. Technol., 51, 6120–6130, https://doi.org/10.1021/acs.est.6b06514, 2017.

Claeys, M., Graham, B., Vas, G., Wang, W., Vermeylen, R., Pashynska, V., Cafmeyer, J., Guyon, P., Andreae, M. O., Artaxo, P., and Maenhaut, W.: Formation of Secondary Organic Aerosols Through Photooxidation of Isoprene, Science, 303, 1173–1176, https://doi.org/10.1126/science.1092805, 2004.

Copernicus Sentinel-5P (processed by ESA): TROPOMI Level 2 Formaldehyde Total Column products, Version 01, European Space Agency [data set], https://doi.org/10.5270/S5P-tjlxfd2, 2018.

Curci, G., Palmer, P. I., Kurosu, T. P., Chance, K., and Visconti, G.: Estimating European volatile organic compound emissions using satellite observations of formaldehyde from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 10, 11501–11517, https://doi.org/10.5194/acp-10-11501-2010, 2010.

Danckaert, T. and Fayt, C.: QDOAS Software user manual, Royal Belgian Institute for Space Aeronomy (BIRA-IASB), Brussels, Belgium, https://uv-vis.aeronomie.be/software/QDOAS/QDOAS_manual.pdf (last access: 5 January 2024), 2017.

Darmenov, A. and da Silva, A.: The quick fire emissions dataset (QFED) – Documentation of versions 2.1, 2.2 and 2.4. NASA//TM-2015-104606, Vol. 38, NASA Global Modeling and
Assimilation Office, 183 pp., https://gmao.gsfc.nasa.gov/pubs/docs/Darmenov796.pdf (last access: 5 January 2024), 2015 (data available at: https://portal.nccs.nasa.gov/datashare/iesa/aerosol/emissions/QFED/v2.4r6/, last access: 5 January 2024).

De Smedt, I., Theys, N., Yu, H., Danckaert, T., Lerot, C., Compernolle, S., Van Roozendael, M., Richter, A., Hilboll, A., Peters, E., Pedergnana, M., Loyola, D., Beirle, S., Wagner, T., Eskes, H., van Geffen, J., Boersma, K. F., and Veefkind, P.: Algorithm theoretical baseline for formaldehyde retrievals from S5P TROPOMI and from the QA4ECV project, Atmos. Meas. Tech., 11, 2395–2426, https://doi.org/10.5194/amt-11-2395-2018, 2018.

De Smedt, I., Pinardi, G., Vigouroux, C., Compernolle, S., Bais, A., Benavent, N., Boersma, F., Chan, K.-L., Donner, S., Eichmann, K.-U., Hedelt, P., Hendrick, F., Irie, H., Kumar, V., Lambert, J.-C., Langerock, B., Lerot, C., Liu, C., Loyola, D., Piters, A., Richter, A., Rivera Cárdenas, C., Romahn, F., Ryan, R. G., Sinha, V., Theys, N., Vlietinck, J., Wagner, T., Wang, T., Yu, H., and Van Roozendael, M.: Comparative assessment of TROPOMI and OMI formaldehyde observations and validation against MAX-DOAS network column measurements, Atmos. Chem. Phys., 21, 12561–12593, https://doi.org/10.5194/acp-21-12561-2021, 2021.

Dufour, G., Wittrock, F., Camredon, M., Beekmann, M., Richter, A., Aumont, B., and Burrows, J. P.: SCIAMACHY formaldehyde observations: constraint for isoprene emission estimates over Europe?, Atmos. Chem. Phys., 9, 1647–1664, https://doi.org/10.5194/acp-9-1647-2009, 2009.

EBAS: Data from European Observatories, EBAS, https://ebas.nilu.no/, last access: 5 January 2024.

Friedrich, M. M., Rivera, C., Stremme, W., Ojeda, Z., Arellano, J., Bezanilla, A., García-Reynoso, J. A., and Grutter, M.: NO₂ vertical profiles and column densities from MAX-DOAS measurements in Mexico City, Atmos. Meas. Tech., 12, 2545–2565, https://doi.org/10.5194/amt-12-2545-2019, 2019.

Frieß, U., Beirle, S., Alvarado Bonilla, L., Bösch, T., Friedrich, M. M., Hendrick, F., Piters, A., Richter, A., van Roozendael, M., Rozanov, V. V., Spinei, E., Tirpitz, J.-L., Vlemmix, T., Wagner, T., and Wang, Y.: Intercomparison of MAX-DOAS vertical profile retrieval algorithms: studies using synthetic data, Atmos. Meas. Tech., 12, 2155–2181, https://doi.org/10.5194/amt-12-2155-2019, 2019.

Gilbert, J. C. and Lemaréchal, C.: Some numerical experiments with variable-storage quasi-Newton algorithms, Math. Program., 45, 407–435, 1989.

Granier, C., Darras, S., Denier Van Der Gon, H., Doubalova, J., Elguindi, N., Galle, B., Gauss, M., Guevara, M., Jalkanen, J. P., Kuenen, J. J. P., Liousse, C., Quack, B., Simpson, D., and Sindelarova, K.: The Copernicus Atmosphere Monitoring Service global and regional emissions (April 2019 version), Copernicus Atmosphere Monitoring Service (CAMS) report, https://doi.org/10.24380/d0bn-kx16, 2019.

Guenther, A., Hewitt, C. N., Erickson, D., Fall, R., Geron, C., Graedel, T., Harley, P., Klinger, L., Lerdau, M., McKay, W. A., Pierce, T., Scholes, B., Steinbrecher, R., Tallamraju, R., Taylor, J., and Zimmerman, P.: A global model of natural volatile organic compound emissions, J. Geophys. Res.-Atmos., 100, 8873–8892, https://doi.org/10.1029/94JD02950, 1995.

Guenther, A., Karl, T., Harley, P., Wiedinmyer, C., Palmer, P. I., and Geron, C.: Estimates of global terrestrial isoprene

emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature), Atmos. Chem. Phys., 6, 3181–3210, https://doi.org/10.5194/acp-6-3181-2006, 2006.

Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492, https://doi.org/10.5194/gmd-5-1471-2012, 2012.

Hallquist, M., Wenger, J. C., Baltensperger, U., Rudich, Y., Simpson, D., Claeys, M., Dommen, J., Donahue, N. M., George, C., Goldstein, A. H., Hamilton, J. F., Herrmann, H., Hoffmann, T., Iinuma, Y., Jang, M., Jenkin, M. E., Jimenez, J. L., Kiendler-Scharr, A., Maenhaut, W., McFiggans, G., Mentel, Th. F., Monod, A., Prévôt, A. S. H., Seinfeld, J. H., Surratt, J. D., Szmigielski, R., and Wildt, J.: The formation, properties and impact of secondary organic aerosol: current and emerging issues, Atmos. Chem. Phys., 9, 5155–5236, https://doi.org/10.5194/acp-9-5155-2009, 2009.

Hase, F., Demoulin, P., Sauval, A. J., Toon, G. C., Bernath, P. F., Goldman, A., Hannigan, J. W., and Rinsland, C. P.: An empirical line-by-line model for the infrared solar transmittance spectrum from 700 to 5000 cm⁻¹, J. Quant. Spectrosc. Ra., 102, 450–463, https://doi.org/10.1016/j.jqsrt.2006.02.026, 2006.

Hendrick, F., Fayt, C., Friedrich, M., Frieß, U., Richter, A., Beirle, S., Wagner, T., and Wang, Y.: Fiducial Reference Measurements for Ground-Based DOAS Air-Quality Observations D6: MAX-DOAS Algorithm ATBD, Tech. rep., Brussels, Belgium, https://frm4doas.aeronomie.be/ProjectDir/Deliverables/FRM4DOAS_D6_MAXDOAS_Algorithm_ATBD_v02_20180130.pdf (last access: 5 January 2024), 2018.

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Bal samo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803, 2020.

Hönninger, G., von Friedeburg, C., and Platt, U.: Multi axis differential optical absorption spectroscopy (MAX-DOAS), Atmos. Chem. Phys., 4, 231–254, https://doi.org/10.5194/acp-4-231-2004, 2004.

Houweling, S., Dentener, F., and Lelieveld, J.: The impact of nonmethane hydrocarbon compounds on tropospheric photochemistry, J. Geophys. Res.-Atmos., 103, 10673–10696, https://doi.org/10.1029/97JD03582, 1998.

Karagkiozidis, D., Friedrich, M. M., Beirle, S., Bais, A., Hendrick, F., Voudouri, K. A., Fountoulakis, I., Karanikolas, A., Tzoumaka, P., Van Roozendael, M., Balis, D., and Wagner, T.: Retrieval of tropospheric aerosol, NO₂, and HCHO vertical profiles from MAX-DOAS observations over Thessaloniki, Greece: intercomparison and validation of two inversion algorithms, Atmos. Meas. Tech., 15, 1269–1301, https://doi.org/10.5194/amt-15-1269-2022, 2022.
Kleipool, Q. L., Dobber, M. R., de Haan, J. F., and Levelt, P. F.: Earth surface reflectance climatology from 3 years of OMI data, J. Geophys. Res.-Atmos., 113, D18308, https://doi.org/10.1029/2008JD010290, 2008.

Kreher, K., Van Roozendael, M., Hendrick, F., Apituley, A., Dimitropoulou, E., Frieß, U., Richter, A., Wagner, T., Lampel, J., Abuhassan, N., Ang, L., Anguas, M., Bais, A., Benavent, N., Bösch, T., Bognar, K., Borovski, A., Bruchkouski, I., Cede, A., Chan, K. L., Donner, S., Drosoglou, T., Fayt, C., Finkenzeller, H., Garcia-Nieto, D., Gielen, C., Gómez-Martín, L., Hao, N., Henzing, B., Herman, J. R., Hermans, C., Hoque, S., Irie, H., Jin, J., Johnston, P., Khayyam Butt, J., Khokhar, F., Koenig, T. K., Kuhn, J., Kumar, V., Liu, C., Ma, J., Merlaud, A., Mishra, A. K., Müller, M., Navarro-Comas, M., Ostendorf, M., Pazmino, A., Peters, E., Pinardi, G., Pinharanda, M., Piters, A., Platt, U., Postylyakov, O., Prados-Roman, C., Puentedura, O., Querel, R., Saiz-Lopez, A., Schönhardt, A., Schreier, S. F., Seyler, A., Sinha, V., Spinei, E., Strong, K., Tack, F., Tian, X., Tiefengraber, M., Tirpitz, J.-L., van Gent, J., Volkamer, R., Vrekoussis, M., Wang, S., Wang, Z., Wenig, M., Wittrock, F., Xie, P. H., Xu, J., Yela, M., Zhang, C., and Zhao, X.: Intercomparison of NO₂, O₄, O₃ and HCHO slant column measurements by MAX-DOAS and zenith-sky UV-visible spectrometers during CINDI-2, Atmos. Meas. Tech., 13, 2169–2208, https://doi.org/10.5194/amt-13-2169-2020, 2020.

Kwon, H.-A., González Abad, G., Nowlan, C. R., Chong, H., Souri, A. H., Vigouroux, C., Röhling, A., Kivi, R., Makarova, M., Notholt, J., Palm, M., Winkler, H., Té, Y., Sussmann, R., Rettinger, M., Mahieu, E., Strong, K., Lutsch, E., Yamanouchi, S., Nagahama, T., Hanningan, J. W., Zhou, M., Murata, I., Grutter, M., Stremme, W., De Mazière, M., Jones, N., Smale, D., and Morino, I.: Validation of OMPS Suomi NPP and OMPS NOAA-20 formaldehyde total columns with NDACC FTIR observations, Earth Space Sci., 10, e2022EA002778, https://doi.org/10.1029/2022EA002778, 2023.

Marais, E. A., Jacob, D. J., Kurosu, T. P., Chance, K., Murphy, J. G., Reeves, C., Mills, G., Casadio, S., Millet, D. B., Barkley, M. P., Paulot, F., and Mao, J.: Isoprene emissions in Africa inferred from OMI observations of formaldehyde columns, Atmos. Chem. Phys., 12, 6219–6235, https://doi.org/10.5194/acp-12-6219-2012, 2012.

Medeiros, D. J., Blitz, M. A., Seakins, P. W., and Whalley, L. K.: Direct measurements of isoprene autoxidation: Pinpointing atmospheric oxidation in tropical forests, J. Am. Chem. Soc., 4, 809–818, https://doi.org/10.1021/jacsau.1c00525, 2022.

MEGAN-MOHYCAN: monthly bottom-up isoprene emissions, Royal Belgian Institute for Space Aeronomy, https://emissions.aeronomie.be, last access: 5 January 2024.

Millet, D. B., Jacob, D. J., Boersma, K. F., Fu, T.-M., Kurosu, T. P., Chance, K., Heald, C. L., and Guenther, A.: Spatial distribution of isoprene emissions from North America derived from formaldehyde column measurements by the OMI satellite sensor, J. Geophys. Res.-Atmos., 113, D02307, https://doi.org/10.1029/2007JD008950, 2008.

Møller, K. H., Bates, K. H., and Kjaergaard, H. G.: The importance of peroxy radical hydrogen-shift reactions in atmospheric isoprene oxidation, J. Phys. Chem. A, 123, 920–930, https://doi.org/10.1021/acs.jpca.8b10432, 2019.

Müller, J.-F. and Stavrakou, T.: Inversion of CO and NOₓ emissions using the adjoint of the IMAGES model, Atmos. Chem. Phys., 5, 1157–1186, https://doi.org/10.5194/acp-5-1157-2005, 2005.

Müller, J.-F., Stavrakou, T., and Peeters, J.: Chemistry and deposition in the Model of Atmospheric composition at Global and Regional scales using Inversion Techniques for Trace gas Emissions (MAGRITTE v1.1) – Part 1: Chemical mechanism, Geosci. Model Dev., 12, 2307–2356, https://doi.org/10.5194/gmd-12-2307-2019, 2019.

Müller, J.-F., Stavrakou, T., Wallens, S., De Smedt, I., Van Roozendael, M., Potosnak, M. J., Rinne, J., Munger, B., Goldstein, A., and Guenther, A. B.: Global isoprene emissions estimated using MEGAN, ECMWF analyses and a detailed canopy environment model, Atmos. Chem. Phys., 8, 1329–1341, https://doi.org/10.5194/acp-8-1329-2008, 2008.

Novelli, A., Vereecken, L., Bohn, B., Dorn, H.-P., Gkatzelis, G. I., Hofzumahaus, A., Holland, F., Reimer, D., Rohrer, F., Rosanka, S., Taraborrelli, D., Tillmann, R., Wegener, R., Yu, Z., Kiendler-Scharr, A., Wahner, A., and Fuchs, H.: Importance of isomerization reactions for OH radical regeneration from the photo-oxidation of isoprene investigated in the atmospheric simulation chamber SAPHIR, Atmos. Chem. Phys., 20, 3333–3355, https://doi.org/10.5194/acp-20-3333-2020, 2020.

Oomen, G.-M., Stavrakou, T., and Müller, J.-F.: Top-down isoprene emissions over Europe based on TROPOMI HCHO data, Royal Belgian Institute for Space Aeronomy [data set], https://doi.org/10.18758/71021088, 2024.

Palmer, P. I., Abbot, D. S., Fu, T.-M., Jacob, D. J., Chance, K., Kurosu, T. P., Guenther, A., Wiedinmyer, C., Stanton, J. C., Pilling, M. J., Pressley, S. N., Lamb, B., and Sumner, A. L.: Quantifying the seasonal and interannual variability of North American isoprene emissions using satellite observations of the formaldehyde column, J. Geophys. Res.-Atmos., 111, D12315, https://doi.org/10.1029/2005JD006689, 2006.

Pan, X., Ichoku, C., Chin, M., Bian, H., Darmenov, A., Colarco, P., Ellison, L., Kucsera, T., da Silva, A., Wang, J., Oda, T., and Cui, G.: Six global biomass burning emission datasets: intercomparison and application in one global aerosol model, Atmos. Chem. Phys., 20, 969–994, https://doi.org/10.5194/acp-20-969-2020, 2020.

Passing, H. and Bablok, W.: A New Biometrical Procedure for Testing the Equality of Measurements from Two Different Analytical Methods. Application of linear regression procedures for method comparison studies in Clinical Chemistry, Part I, Clin. Chem. Lab. Med., 21, 709–720, https://doi.org/10.1515/cclm.1983.21.11.709, 1983.

Passing, H. and Bablok, W.: Comparison of Several Regression Procedures for Method Comparison Studies and Determination of Sample Sizes Application of linear regression procedures for method comparison studies in Clinical Chemistry, Part II, J. Clin. Chem. Clin. Bio., 22, 431–445, https://doi.org/10.1515/cclm.1984.22.6.431, 1984.

Peeters, J., Müller, J.-F., Stavrakou, T., and Nguyen, V. S.: Hydroxyl Radical Recycling in Isoprene Oxidation Driven by Hydrogen Bonding and Hydrogen Tunneling: The Upgraded LIM1 Mechanism, J. Phys. Chem. A, 118, 8625–8643, https://doi.org/10.1021/jp5033146, 2014.
Platt, U. and Stutz, J.: Differential Optical Absorption Spectroscopy, Springer, https://doi.org/10.1007/978-3-540-75776-4, 2008.

Possell, M. and Hewitt, C. N.: Isoprene emissions from plants are mediated by atmospheric CO₂ concentrations, Glob. Change Biol., 17, 1595–1610, https://doi.org/10.1111/j.1365-2486.2010.02306.x, 2011.

Pougatchev, N. S., Connor, B. J., and Rinsland, C. P.: Infrared measurements of the ozone vertical distribution above Kitt Peak, J. Geophys. Res., 100, 16689–16697, https://doi.org/10.1029/95JD01296, 1995.

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding, WORLD SCIENTIFIC, https://doi.org/10.1142/3171, 2000.

Rodgers, C. D. and Connor, B. J.: Intercomparison of remote sounding instruments, J. Geophys. Res.-Atmos., 108, 4116, https://doi.org/10.1029/2002JD002299, 2003.

Rothman, L. S., Gordon, I. E., Babikov, Y., Barbe, A., Chris Benner, D., Bernath, P. F., Birk, M., Bizzocchi, L., Boudon, V., Brown, L. R., Campargue, A., Chance, K., Cohen, E. A., Coudert, L. H., Devi, V. M., Drouin, B. J., Fayt, A., Flaud, J. M., Gamache, R. R., Harrison, J. J., Hartmann, J. M., Hill, C., Hodges, J. T., Jacquemart, D., Jolly, A., Lamouroux, J., Le Roy, R. J., Li, G., Long, D. A., Lyulin, O. M., Mackie, C. J., Massie, S. T., Mikhailenko, S., Müller, H. S. P., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V., Perrin, A., Polovtseva, E. R., Richard, C., Smith, M. A. H., Starikova, E., Sung, K., Tashkun, S., Tennyson, J., Toon, G. C., Tyuterev, V. G., and Wagner, G.: The HITRAN2012 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 130, 4–50, https://doi.org/10.1016/j.jqsrt.2013.07.002, 2013.

Schwantes, R. H., Charan, S. M., Bates, K. H., Huang, Y., Nguyen, T. B., Mai, H., Kong, W., Flagan, R. C., and Seinfeld, J. H.: Low-volatility compounds contribute significantly to isoprene secondary organic aerosol (SOA) under high-NOₓ conditions, Atmos. Chem. Phys., 19, 7255–7278, https://doi.org/10.5194/acp-19-7255-2019, 2019.

Schwantes, R. H., Emmons, L. K., Orlando, J. J., Barth, M. C., Tyn-dall, G. S., Hall, S. R., Ullmann, K., St. Clair, J. M., Blake, D. R., Wisthaler, A., and Bui, T. P. V.: Comprehensive isoprene and terpene gas-phase chemistry improves simulated surface ozone in the southeastern US, Atmos. Chem. Phys., 20, 3739–3776, https://doi.org/10.5194/acp-20-3739-2020, 2020.

SEEDS: daily top-down isoprene and biomass burning emission fluxes from 2018 to 2022, SEEDS, https://www.seedsproject.eu/data, last access: 5 January 2024.

Seinfeld, J. and Pandis, S.: Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, ISBN 9781118947401, https://www.wiley.com/en-ie/Atmospheric+Chemistry+and+Physics:+From+Air+Pollution+to+Climate+Change,+3rd+Edition-p-9781118947401 (last access: 5 January 2024), 2016.

Sindelarova, K., Granier, C., Bouarar, I., Guenther, A., Tilmes, S., Stavrakou, T., Müller, J.-F., Kuhn, U., Stefani, P., and Knorr, W.: Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years, Atmos. Chem. Phys., 14, 9317–9341, https://doi.org/10.5194/acp-14-9317-2014, 2014.

Sindelarova, K., Markova, J., Simpson, D., Huszar, P., Karlicky, J., Darras, S., and Granier, C.: High-resolution biogenic global emission inventory for the time period 2000–2019 for

air quality modelling, Earth Syst. Sci. Data, 14, 251–270, https://doi.org/10.5194/essd-14-251-2022, 2022.

Sousa, P. M., Barriopedro, D., García-Herrera, R., Ordóñez, C., Soares, P. M. M., and Trigo, R. M.: Distinct influences of large-scale circulation and regional feedbacks in two exceptional 2019 European heatwaves, Commun. Earth Env., 1, 48, https://doi.org/10.1038/s43247-020-00048-9, 2020.

Spurr, R., de Haan, J., van Oss, R., and Vasilkov, A.: Discrete-ordinate radiative transfer in a stratified medium with first-order rotational Raman scattering, J. Quant. Spectrosc. Ra., 109, 404–425, https://doi.org/10.1016/j.jqsrt.2007.08.011, 2008.

Stavrakou, T., Müller, J.-F., De Smedt, I., Van Roozendael, M., van der Werf, G. R., Giglio, L., and Guenther, A.: Evaluating the performance of pyrogenic and biogenic emission inventories against one decade of space-based formaldehyde columns, Atmos. Chem. Phys., 9, 1037–1060, https://doi.org/10.5194/acp-9-1037-2009, 2009.

Stavrakou, T., Müller, J. F., Peeters, J., Razavi, A., Clarisse, L., Clerbaux, C., Coheur, P. F., Hurtmans, D., de Mazière, M., Vigouroux, C., Deutscher, N. M., Griffith, D. W. T., Jones, N., and Paton-Walsh, C.: Satellite evidence for a large source of formic acid from boreal and tropical forests, Nat. Geosci., 5, 26–30, https://doi.org/10.1038/ngeo1354, 2012.

Stavrakou, T., Müller, J.-F., Bauwens, M., De Smedt, I., Van Roozendael, M., De Mazière, M., Vigouroux, C., Hendrick, F., George, M., Clerbaux, C., Coheur, P.-F., and Guenther, A.: How consistent are top-down hydrocarbon emissions based on formaldehyde observations from GOME-2 and OMI?, Atmos. Chem. Phys., 15, 11861–11884, https://doi.org/10.5194/acp-15-11861-2015, 2015.

Stavrakou, T., Müller, J. F., Bauwens, M., De Smedt, I., Van Roozendael, M., and Guenther, A.: Impact of Short-Term Climate Variability on Volatile Organic Compounds Emissions Assessed Using OMI Satellite Formaldehyde Observations, Geophys. Res. Lett., 45, 8681–8689, https://doi.org/10.1029/2018GL078676, 2018.

Teng, A. P., Crounse, J. D., and Wennberg, P.: Isoprene peroxy radical dynamics, J. Am. Chem. Soc., 139, 5367–5677, https://doi.org/10.1021/jacs.6b12838, 2017.

Tikhonov, A.: Solution of Incorrectly Formulated Problems and the Regularization Method, Dokl. Acad. Nauk SSSR, 151, 501–504, 1963.

Tirpitz, J.-L., Frieß, U., Hendrick, F., Alberti, C., Allaart, M., Apituley, A., Bais, A., Beirle, S., Berkhout, S., Bognar, K., Bösch, T., Bruchkouski, I., Cede, A., Chan, K. L., den Hoed, M., Donner, S., Drosoglou, T., Fayt, C., Friedrich, M. M., Frumau, A., Gast, L., Gielen, C., Gomez-Martín, L., Hao, N., Hensen, A., Henzing, B., Hermans, C., Jin, J., Kreher, K., Kuhn, J., Lampel, J., Li, A., Liu, C., Liu, H., Ma, J., Merlaud, A., Peters, E., Pinardi, G., Piters, A., Platt, U., Puentedura, O., Richter, A., Schmitt, S., Spinei, E., Stein Zweers, D., Strong, K., Swart, D., Tack, F., Tiefengraber, M., van der Hoff, R., van Roozendael, M., Vlemmix, T., Vonk, J., Wagner, T., Wang, Y., Wang, Z., Wenig, M., Wiegner, M., Witrock, F., Xie, P., Xing, C., Xu, J., Yela, M., Zhang, C., and Zhao, X.: Intercomparison of MAX-DOAS vertical profile retrieval algorithms: studies on field data from the CINDI-2 campaign, Atmos. Meas. Tech., 14, 1–35, https://doi.org/10.5194/amt-14-1-2021, 2021.
Tørseth, K., Aas, W., Breivik, K., Fjæraa, A. M., Fiebig, M., Hjellbrekke, A. G., Lund Myhre, C., Solberg, S., and Yttri, K. E.: Introduction to the European Monitoring and Evaluation Programme (EMEP) and observed atmospheric composition change during 1972–2009, Atmos. Chem. Phys., 12, 5447–5481, https://doi.org/10.5194/acp-12-5447-2012, 2012.

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720, https://doi.org/10.5194/essd-9-697-2017, 2017.

Vigouroux, C., Hendrick, F., Stavrakou, T., Dils, B., De Smedt, I., Hermans, C., Merlaud, A., Scolas, F., Senten, C., Vanhaelewyn, G., Fally, S., Carleer, M., Metzger, J.-M., Müller, J.-F., Van Roozendael, M., and De Mazière, M.: Ground-based FTIR and MAX-DOAS observations of formaldehyde at Réunion Island and comparisons with satellite and model data, Atmos. Chem. Phys., 9, 9523–9544, https://doi.org/10.5194/acp-9-9523-2009, 2009.

Vigouroux, C., Bauer Aquino, C. A., Bauwens, M., Becker, C., Blumenstock, T., De Mazière, M., García, O., Grutter, M., Guarin, C., Hannigan, J., Hase, F., Jones, N., Kivi, R., Koshelev, D., Langerock, B., Lutsch, E., Makarova, M., Metzger, J.-M., Müller, J.-F., Notholt, J., Ortega, I., Palm, M., Paton-Walsh, C., Poberovskii, A., Rettinger, M., Robinson, J., Smale, D., Stavrakou, T., Stremme, W., Strong, K., Sussmann, R., Té, Y., and Toon, G.: NDACC harmonized formaldehyde time series from 21 FTIR stations covering a wide range of column abundances, Atmos. Meas. Tech., 11, 5049–5073, https://doi.org/10.5194/amt-11-5049-2018, 2018.

Vigouroux, C., Langerock, B., Bauer Aquino, C. A., Blumenstock, T., Cheng, Z., De Mazière, M., De Smedt, I., Grutter, M., Hannigan, J. W., Jones, N., Kivi, R., Loyola, D., Lutsch, E., Mahieu, E., Makarova, M., Metzger, J.-M., Morino, I., Murata, I., Nagahama, T., Notholt, J., Ortega, I., Palm, M., Pinardi, G., Röhling, A., Smale, D., Stremme, W., Strong, K., Sussmann, R., Té, Y., van Roozendael, M., Wang, P., and Winkler, H.: TROPOMI-Sentinel-5 Precursor formaldehyde validation using an extensive network of ground-based Fourier-transform infrared stations, Atmos. Meas. Tech., 13, 3751–3767, https://doi.org/10.5194/amt-13-3751-2020, 2020.

Wennberg, P. O., Bates, K. H., Crounse, J. D., Dodson, L. G., McVay, R. C., Mertens, L. A., Nguyen, T. B., Praske, E., Schwantes, R. H., Smarte, M. D., St Clair, J. M., Teng, A. P., Zhang, X., and Seinfeld, J. H.: Gas-Phase Reactions of Isoprene and Its Major Oxidation Products, Chem. Rev., 118, 3337–3390, https://doi.org/10.1021/acs.chemrev.7b00439, 2018.

Williams, J. E., Boersma, K. F., Le Sager, P., and Verstraeten, W. W.: The high-resolution version of TM5-MP for optimized satellite retrievals: description and validation, Geosci. Model Dev., 10, 721–750, https://doi.org/10.5194/gmd-10-721-2017, 2017.

Zhu, L., González Abad, G., Nowlan, C. R., Chan Miller, C., Chance, K., Apel, E. C., DiGangi, J. P., Fried, A., Hanisco, T. F., Hornbrook, R. S., Hu, L., Kaiser, J., Keutsch, F. N., Permar, W., St. Clair, J. M., and Wolfe, G. M.: Validation of satellite formaldehyde (HCHO) retrievals using observations from 12 aircraft campaigns, Atmos. Chem. Phys., 20, 12329–12345, https://doi.org/10.5194/acp-20-12329-2020, 2020.