Laboratoire de l'Atmosphère et des Cyclones

Abstract. The Lagrangian particle dispersion model FLEXPART in its original version in the mid-1990s was designed for calculating the long-range and mesoscale dispersion of hazardous substances from point sources, such as those released after an accident in a nuclear power plant. Over the past decades, the model has evolved into a comprehensive tool for multi-scale atmospheric transport modeling and analysis and has attracted a global user community. Its application fields have been extended to a large range of atmospheric gases and aerosols, e.g., greenhouse gases, short-lived climate forcers like black carbon and volcanic ash, and it has also been used to study the atmospheric branch of the water cycle. Given suitable meteorological input data, it can be used for scales from dozens of meters to global. In particular, inverse modeling based on source–receptor relationships from FLEXPART has become widely used. In this paper, we present FLEXPART version 10.4, which works with meteorological input data from the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) and data from the United States National Centers of Environmental Prediction (NCEP) Global Forecast System (GFS). Since the last publication of a detailed FLEXPART description (version 6.2), the model has been improved in different aspects such as performance, physicochemical parameterizations, input/output formats, and available preprocessing and post-processing software. The model code has also been parallelized using the Message Passing Interface (MPI). We demonstrate that the model scales well up to using 256 processors, with a parallel efficiency greater than 75 % for up to 64 processes on multiple nodes in runs with very large numbers of particles. The deviation from 100 % efficiency is almost entirely due to the remaining nonparallelized parts of the code, suggesting large potential for further speedup. A new turbulence scheme for the convective boundary layer has been developed that considers the skewness in the vertical velocity distribution (updrafts and downdrafts) and vertical gradients in air density. FLEXPART is the only model available considering both effects, making it highly accurate for small-scale applications, e.g., to quantify dispersion in the vicinity of a point source. The wet deposition scheme for aerosols has been completely rewritten and a new, more detailed gravitational settling parameterization for aerosols has also been implemented. FLEXPART has had the option of running backward in time from atmospheric concentrations at receptor locations for many years, but this has now been extended to also work for deposition values and may become useful, for instance, for the interpretation of ice core measurements. To our knowledge, to date FLEXPART is the only model with that capability. Furthermore, the temporal variation and temperature dependence of chemical reactions with the OH radical have been included, allowing for more accurate simulations for species with intermediate lifetimes against the reaction with OH, such as ethane. Finally, user settings can now be specified in a more flexible namelist format, and output files can be produced in NetCDF format instead of FLEXPART's customary binary format. In this paper, we describe these new developments. Moreover, we present some tools for the preparation of the meteorological input data and for processing FLEXPART output data, and we briefly report on alternative FLEXPART versions.

Abstract. This paper reports on consolidated ground-based validation results of the atmospheric NO2 data produced operationally since April 2018 by the TROPOspheric Monitoring Instrument (TROPOMI) on board of the ESA/EU Copernicus Sentinel-5 Precursor (S5P) satellite. Tropospheric, stratospheric, and total NO2 column data from S5P are compared to correlative measurements collected from, respectively, 19 Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS), 26 Network for the Detection of Atmospheric Composition Change (NDACC) Zenith-Scattered-Light DOAS (ZSL-DOAS), and 25 Pandonia Global Network (PGN)/Pandora instruments distributed globally. The validation methodology gives special care to minimizing mismatch errors due to imperfect spatio-temporal co-location of the satellite and correlative data, e.g. by using tailored observation operators to account for differences in smoothing and in sampling of atmospheric structures and variability and photochemical modelling to reduce diurnal cycle effects. Compared to the ground-based measurements, S5P data show, on average, (i) a negative bias for the tropospheric column data, of typically −23 % to −37 % in clean to slightly polluted conditions but reaching values as high as −51 % over highly polluted areas; (ii) a slight negative median difference for the stratospheric column data, of about −0.2 Pmolec cm−2, i.e. approx. −2 % in summer to −15 % in winter; and (iii) a bias ranging from zero to −50 % for the total column data, found to depend on the amplitude of the total NO2 column, with small to slightly positive bias values for columns below 6 Pmolec cm−2 and negative values above. The dispersion between S5P and correlative measurements contains mostly random components, which remain within mission requirements for the stratospheric column data (0.5 Pmolec cm−2) but exceed those for the tropospheric column data (0.7 Pmolec cm−2). While a part of the biases and dispersion may be due to representativeness differences such as different area averaging and measurement times, it is known that errors in the S5P tropospheric columns exist due to shortcomings in the (horizontally coarse) a priori profile representation in the TM5-MP chemical transport model used in the S5P retrieval and, to a lesser extent, to the treatment of cloud effects and aerosols. Although considerable differences (up to 2 Pmolec cm−2 and more) are observed at single ground-pixel level, the near-real-time (NRTI) and offline (OFFL) versions of the S5P NO2 operational data processor provide similar NO2 column values and validation results when globally averaged, with the NRTI values being on average 0.79 % larger than the OFFL values.

HyMeX-SOP1 collected unprecedented observations of atmosphere, ocean, land, and rivers

Abstract. This paper presents the Meso-NH model version 5.4. Meso-NH is an atmospheric non hydrostatic research model that is applied to a broad range of resolutions, from synoptic to turbulent scales, and is designed for studies of physics and chemistry. It is a limited-area model employing advanced numerical techniques, including monotonic advection schemes for scalar transport and fourth-order centered or odd-order WENO advection schemes for momentum. The model includes state-of-the-art physics parameterization schemes that are important to represent convective-scale phenomena and turbulent eddies, as well as flows at larger scales. In addition, Meso-NH has been expanded to provide capabilities for a range of Earth system prediction applications such as chemistry and aerosols, electricity and lightning, hydrology, wildland fires, volcanic eruptions, and cyclones with ocean coupling. Here, we present the main innovations to the dynamics and physics of the code since the pioneer paper of Lafore et al. (1998) and provide an overview of recent applications and couplings.

Emissions of gases and particles from the combustion of fossil fuels and biofuels in Africa are expected to increase significantly in the near future due to the rapid growth of African cities and megacities. There is currently no regional emissions inventory that provides estimates of anthropogenic combustion for the African continent. This work provides a quantification of the evolution of African combustion emissions from 2005 to 2030, using a bottom-up method. This inventory predicts very large increases in black carbon, organic carbon, CO, NOx, SO2 and non-methane hydrocarbon emissions if no emission regulations are implemented. This paper discusses the effectiveness of scenarios involving certain fuels, specific to Africa in each activity sector and each region (western, eastern, northern and southern Africa), to reduce the emissions. The estimated trends in African emissions are consistent with emissions provided by global inventories, but they display a larger range of values. African combustion emissions contributed significantly to global emissions in 2005. This contribution will increase more significantly by 2030: organic carbon emissions will for example make up 50% of the global emissions in 2030. Furthermore, we show that the magnitude of African anthropogenic emissions could be similar to African biomass burning emissions around 2030. © 2014 IOP Publishing Ltd.

We compare Tropospheric Emission Spectrometer (TES) version 2 (V002) nadir ozone profiles with ozonesonde profiles from the Intercontinental Chemical Transport Experiment Ozonesonde Network Study, the World Ozone and Ultraviolet Data Center, the Global Monitoring Division of the Earth System Research Laboratory, and the Southern Hemisphere Additional Ozonesonde archives. Approximately 1600 coincidences spanning 72.5°S–80.3°N from October 2004 to October 2006 are found. The TES averaging kernel and constraint are applied to the ozonesonde data to account for the TES measurement sensitivity and vertical resolution. TES sonde differences are examined in six latitude zones after excluding profiles with thick high clouds. Values for the bias and standard deviation are determined using correlations of mean values of TES ozone and sonde ozone in the upper troposphere (UT) and lower troposphere (LT). The UT biases range from 2.9 to 10.6 ppbv, and the LT biases range from 3.7 to 9.2 ppbv, excluding the Arctic and Antarctic LT where TES sensitivity is low. A similar approach is used to assess seasonal differences in the northern midlatitudes where the density and frequency of sonde measurements are greatest. These results are briefly compared to TES V001 ozone validation work which also used ozonesondes but was carried out prior to improvements in the radiometric calibration and ozone retrieval in V002. Overall, the large number of TES and sonde comparisons indicate a positive bias of approximately 3–10 ppbv for the TES V002 nadir ozone data set and have helped to identify areas of potential improvement for future retrieval versions.

We present validation studies of MLS version 2.2 upper tropospheric and stratospheric ozone profiles using ozonesonde and lidar data as well as climatological data. Ozone measurements from over 60 ozonesonde stations worldwide and three lidar stations are compared with coincident MLS data. The MLS ozone stratospheric data between 150 and 3 hPa agree well with ozonesonde measurements, within 8% for the global average. MLS values at 215 hPa are biased high compared to ozonesondes by ∼20% at middle to high latitude, although there is a lot of variability in this altitude region. Comparisons between MLS and ground‐based lidar measurements from Mauna Loa, Hawaii, from the Table Mountain Facility, California, and from the Observatoire de Haute‐Provence, France, give very good agreement, within ∼5%, for the stratospheric values. The comparisons between MLS and the Table Mountain Facility tropospheric ozone lidar show that MLS data are biased high by ∼30% at 215 hPa, consistent with that indicated by the ozonesonde data. We obtain better global average agreement between MLS and ozonesonde partial column values down to 215 hPa, although the average MLS values at low to middle latitudes are higher than the ozonesonde values by up to a few percent. MLS v2.2 ozone data agree better than the MLS v1.5 data with ozonesonde and lidar measurements. MLS tropical data show the wave one longitudinal pattern in the upper troposphere, with similarities to the average distribution from ozonesondes. High upper tropospheric ozone values are also observed by MLS in the tropical Pacific from June to November.

Abstract The Southern Hemisphere ADditional OZonesonde (SHADOZ) network was assembled to validate a new generation of ozone‐monitoring satellites and to better characterize the vertical structure of tropical ozone in the troposphere and stratosphere. Beginning with nine stations in 1998, more than 7,000 ozone and P ‐ T ‐ U profiles are available from 14 SHADOZ sites that have operated continuously for at least a decade. We analyze ozone profiles from the recently reprocessed SHADOZ data set that is based on adjustments for inconsistencies caused by varying ozonesonde instruments and operating techniques. First, sonde‐derived total ozone column amounts are compared to the overpasses from the Earth Probe/Total Ozone Mapping Spectrometer, Ozone Monitoring Instrument, and Ozone Mapping and Profiler Suite satellites that cover 1998–2016. Second, characteristics of the stratospheric and tropospheric columns are examined along with ozone structure in the tropical tropopause layer (TTL). We find that (1) relative to our earlier evaluations of SHADOZ data, in 2003, 2007, and 2012, sonde‐satellite total ozone column offsets at 12 stations are 2% or less, a significant improvement; (2) as in prior studies, the 10 tropical SHADOZ stations, defined as within ±19° latitude, display statistically uniform stratospheric column ozone, 229 ± 3.9 DU (Dobson units), and a tropospheric zonal wave‐one pattern with a 14 DU mean amplitude; (3) the TTL ozone column, which is also zonally uniform, masks complex vertical structure, and this argues against using satellites for lower stratospheric ozone trends; and (4) reprocessing has led to more uniform stratospheric column amounts across sites and reduced bias in stratospheric profiles. As a consequence, the uncertainty in total column ozone now averages 5%.

Abstract Reprocessed ozonesonde data from eight SHADOZ (Southern Hemisphere ADditional OZonesondes) sites have been used to derive the first analysis of uncertainty estimates for both profile and total column ozone (TCO). The ozone uncertainty is a composite of the uncertainties of the individual terms in the ozone partial pressure ( P O3 ) equation, those being the ozone sensor current, background current, internal pump temperature, pump efficiency factors, conversion efficiency, and flow rate. Overall, P O3 uncertainties (Δ P O3 ) are within 15% and peak around the tropopause (15 ± 3 km) where ozone is a minimum and Δ P O3 approaches the measured signal. The uncertainty in the background and sensor currents dominates the overall Δ P O3 in the troposphere including the tropopause region, while the uncertainties in the conversion efficiency and flow rate dominate in the stratosphere. Seasonally, Δ P O3 is generally a maximum in the March–May, with the exception of SHADOZ sites in Asia, for which the highest Δ P O3 occurs in September–February. As a first approach, we calculate sonde TCO uncertainty (ΔTCO) by integrating the profile ΔP O3 and adding the ozone residual uncertainty, derived from the McPeters and Labow (2012, doi:10.1029/2011JD017006) 1 σ ozone mixing ratios. Overall, ΔTCO are within ±15 Dobson units (DU), representing ~5–6% of the TCO. Total Ozone Mapping Spectrometer and Ozone Monitoring Instrument (TOMS and OMI) satellite overpasses are generally within the sonde ΔTCO. However, there is a discontinuity between TOMS v8.6 (1998 to September 2004) and OMI (October 2004–2016) TCO on the order of 10 DU that accounts for the significant 16 DU overall difference observed between sonde and TOMS. By comparison, the sonde‐OMI absolute difference for the eight stations is only ~4 DU.

Abstract. Aerosol particles are essential constituents of the Earth's atmosphere, impacting the earth radiation balance directly by scattering and absorbing solar radiation, and indirectly by acting as cloud condensation nuclei. In contrast to most greenhouse gases, aerosol particles have short atmospheric residence times, resulting in a highly heterogeneous distribution in space and time. There is a clear need to document this variability at regional scale through observations involving, in particular, the in situ near-surface segment of the atmospheric observation system. This paper will provide the widest effort so far to document variability of climate-relevant in situ aerosol properties (namely wavelength dependent particle light scattering and absorption coefficients, particle number concentration and particle number size distribution) from all sites connected to the Global Atmosphere Watch network. High-quality data from almost 90 stations worldwide have been collected and controlled for quality and are reported for a reference year in 2017, providing a very extended and robust view of the variability of these variables worldwide. The range of variability observed worldwide for light scattering and absorption coefficients, single-scattering albedo, and particle number concentration are presented together with preliminary information on their long-term trends and comparison with model simulation for the different stations. The scope of the present paper is also to provide the necessary suite of information, including data provision procedures, quality control and analysis, data policy, and usage of the ground-based aerosol measurement network. It delivers to users of the World Data Centre on Aerosol, the required confidence in data products in the form of a fully characterized value chain, including uncertainty estimation and requirements for contributing to the global climate monitoring system.

Abstract The underwater Hunga Tonga-Hunga Ha-apai volcano erupted in the early hours of 15th January 2022, and injected volcanic gases and aerosols to over 50 km altitude. Here we synthesise satellite, ground-based, in situ and radiosonde observations of the eruption to investigate the strength of the stratospheric aerosol and water vapour perturbations in the initial weeks after the eruption and we quantify the net radiative impact across the two species using offline radiative transfer modelling. We find that the Hunga Tonga-Hunga Ha-apai eruption produced the largest global perturbation of stratospheric aerosols since the Pinatubo eruption in 1991 and the largest perturbation of stratospheric water vapour observed in the satellite era. Immediately after the eruption, water vapour radiative cooling dominated the local stratospheric heating/cooling rates, while at the top-of-the-atmosphere and surface, volcanic aerosol cooling dominated the radiative forcing. However, after two weeks, due to dispersion/dilution, water vapour heating started to dominate the top-of-the-atmosphere radiative forcing, leading to a net warming of the climate system.

Large volcanic eruptions, although rare events, can influence the chemistry and the dynamics of the stratosphere for several years after the eruption. Here we show that the eruption of the submarine volcano Hunga Tonga-Hunga Ha'apai on 15 January 2022 injected at least 50 teragrams of water vapor directly into the stratosphere. This event raised the amount of water vapor in the developing stratospheric plume by several orders of magnitude and possibly increased the amount of global stratospheric water vapor by more than 5%. This extraordinary eruption may have initiated an atmospheric response different from that of previous well-studied large volcanic eruptions.

Abstract. Since the nineties, atmospheric measurement systems have been deployed at Reunion Island, mainly for monitoring the atmospheric composition in the framework of NDSC/NDACC (Network for the Detection of Stratospheric Change/Network for the Detection of Atmospheric Composition Change). The location of Reunion Island presents a great interest because there are very few multi-instrumented stations in the tropics and particularly in the southern hemisphere. In 2012, a new observatory was commissioned in Maïdo at 2200 m above sea level: it hosts various instruments for atmospheric measurements, including lidar systems, spectro-radiometers and in situ gas and aerosol measurements. This new high-altitude Maïdo station provides an opportunity: 1. to improve the performance of the optical instruments above the marine boundary layer, and to open new perspectives on upper troposphere and lower stratosphere studies; 2. to develop in situ measurements of the atmospheric composition for climate change surveys, in a reference site in the tropical/subtropical region of the southern hemisphere; 3. to offer trans-national access to host experiments or measurement campaigns for focused process studies.

Abstract. Formaldehyde (HCHO) columns have been retrieved from ground-based Fourier transform infrared (FTIR) campaign measurements in 2004 and 2007 and from UV-Visible MAX-DOAS measurements in 2004–2005 at the NDACC site of Réunion Island (21° S, 55° E). The FTIR and MAX-DOAS daily mean formaldehyde total columns are intercompared in their common measurement period, from August to October 2004. The ground-based data are also compared to correlative SCIAMACHY data. The comparisons account for the vertical sensitivity differences of the data sets, by including their respective averaging kernels. Complete error budgets are also presented. The FTIR and MAX-DOAS daily mean total columns agree very well: no significant bias is observed and the standard deviation of the comparisons is only 8%. Both FTIR and MAX-DOAS HCHO total columns are in good agreement with SCIAMACHY values in the 2004–2005 period, with standard deviations of 21% and 31%, respectively. The same seasonal cycle is observed by the different instruments, with a minimum in austral winter and a maximum in February–March. The FTIR and MAX-DOAS data are confronted with HCHO columns calculated by a global CTM, the IMAGES model. The model underestimates the HCHO columns by 23–29% in comparison with FTIR, and by 15% in comparison with DOAS. This bias might have multiple causes, including an underestimation of OH concentrations in the model (as indicated by a sensitivity study using prescribed OH fields) and/or an underestimated contribution of large-scale transport of HCHO precursors from Madagascar. The latter hypothesis is comforted by the large observed day-to-day variability of HCHO columns, and by the observation that the peak values of FTIR columns can often be associated with free tropospheric transport patterns from source regions over Madagascar to Réunion Island, according to simulations performed with the Lagrangian particle dispersion model FLEXPART.

Here we present extensive observations of stratospheric and upper tropospheric water vapor using the balloon‐borne Cryogenic Frost point Hygrometer (CFH) in support of the Aura Microwave Limb Sounder (MLS) satellite instrument. Coincident measurements were used for the validation of MLS version 1.5 and for a limited validation of MLS version 2.2 water vapor. The sensitivity of MLS is on average 30% lower than that of CFH, which is fully compensated by a constant offset at stratospheric levels but only partially compensated at tropospheric levels, leading to an upper tropospheric dry bias. The sensitivity of MLS observations may be adjusted using the correlation parameters provided here. For version 1.5 stratospheric observations at pressures of 68 hPa and smaller MLS retrievals and CFH in situ observations agree on average to within 2.3% ± 11.8%. At 100 hPa the agreement is to within 6.4% ± 22% and at upper tropospheric pressures to within 23% ± 37%. In the tropical stratosphere during the boreal winter the agreement is not as good. The “tape recorder” amplitude in MLS observations depends on the vertical profile of water vapor mixing ratio and shows a significant interannual variation. The agreement between stratospheric observations by MLS version 2.2 and CFH is comparable to the agreement using MLS version 1.5. The variability in the difference between observations by MLS version 2.2 and CFH at tropospheric levels is significantly reduced, but a tropospheric dry bias and a reduced sensitivity remain in this version. In the validation data set a dry bias at 177.8 hPa of −24.1% ± 16.0% is statistically significant.

Abstract The eruption of the submarine Hunga volcano in January 2022 was associated with a powerful blast that injected volcanic material to altitudes up to 58 km. From a combination of various types of satellite and ground-based observations supported by transport modeling, we show evidence for an unprecedented increase in the global stratospheric water mass by 13% relative to climatological levels, and a 5-fold increase of stratospheric aerosol load, the highest in the last three decades. Owing to the extreme injection altitude, the volcanic plume circumnavigated the Earth in only 1 week and dispersed nearly pole-to-pole in three months. The unique nature and magnitude of the global stratospheric perturbation by the Hunga eruption ranks it among the most remarkable climatic events in the modern observation era, with a range of potential long-lasting repercussions for stratospheric composition and climate.

Abstract. This study presents the principles of the new coupling interface based on the SURFEX multi-surface model and the OASIS3-MCT coupler. As SURFEX can be plugged into several atmospheric models, it can be used in a wide range of applications, from global and regional coupled climate systems to high-resolution numerical weather prediction systems or very fine-scale models dedicated to process studies. The objective of this development is to build and share a common structure for the atmosphere–surface coupling of all these applications, involving on the one hand atmospheric models and on the other hand ocean, ice, hydrology, and wave models. The numerical and physical principles of SURFEX interface between the different component models are described, and the different coupled systems in which the SURFEX OASIS3-MCT-based coupling interface is already implemented are presented.

Abstract. Spectral measurements of BrO using zenith-sky and off-axis viewing geometries are combined in a linear multiple regression retrieval algorithm to provide stratospheric and tropospheric BrO vertical columns. One year of measurement data are investigated over Reunion-Island (20.9° S, 55.5° E), from August 2004 to June 2005. A comparison between the stratospheric columns retrieved at 45°, 80°, 85°, 87.5° and 92.5° solar zenith angles and photochemical simulations initialized by chemical fields from the 3-D-CTM SLIMCAT and further constrained by observed NO2 profiles shows a good agreement only by considering a contribution from the very short-lived organic bromine substances to the stratospheric inorganic bromine budget, of 6 to 8 pptv. Furthermore, stratospheric BrO profiles retrieved from late twilight zenith-sky observations are consistent with a total inorganic bromine (Bry) loading of approximately 23 pptv. This represents 6 to 7 pptv more than can be supplied by long-lived organic bromine sources, and therefore supports an added contribution from very short-lived organic bromine substances as recently suggested in several other studies. Moreover strong evidences are presented for the existence of a substantial amount of BrO in the tropical free-troposphere, around 6 km altitude, possibly supplied by the decomposition of short-lived biogenic bromine organic compounds. Tropospheric BrO vertical columns of 1.1±0.45×1013 molec/cm² are derived for the entire observation period. Comparisons between ground-based BrO vertical columns and total BrO columns derived from SCIAMACHY (onboard the ENVISAT satellite) nadir observations in a latitudinal band centered around 21° S present a good level of consistency, which further strengthens the conclusions of our study.

Abstract Tropical cyclones are some of the most devastating natural hazards and the “three beasts”—Harvey, Irma, and Maria—during the Atlantic hurricane season 2017 are recent examples. The European Centre for Medium-Range Weather Forecasts (ECMWF) is working on fulfilling its 2016–25 strategy in which early warnings for extreme events will be made possible by a high-resolution Earth system ensemble forecasting system. Several verification reports acknowledge deterministic and probabilistic tropical cyclone tracks from ECMWF as world leading. However, producing reliable intensity forecasts is still a difficult task for the ECMWF global forecasting model, especially regarding maximum wind speed. This article will put the ECMWF strategy into a tropical cyclone perspective and highlight some key research activities, using Harvey, Irma, and Maria as examples. We describe the observation usage around tropical cyclones in data assimilation and give examples of their impact. From a model perspective, we show the impact of running at 5-km resolution and also the impact of applying ocean coupling. Finally, we discuss the future challenges to tackle the errors in intensity forecasts for tropical cyclones.

Abstract. Reunion Island (21° S, 55° E), situated in the Indian Ocean at about 800 km east of Madagascar, is appropriately located to monitor the outflow of biomass burning pollution from Southern Africa and Madagascar, in the case of short-lived compounds, and from other Southern Hemispheric landmasses such as South America, in the case of longer-lived species. Ground-based Fourier transform infrared (FTIR) solar absorption observations are sensitive to a large number of biomass burning products. We present in this work the FTIR retrieval strategies, suitable for very humid sites such as Reunion Island, for hydrogen cyanide (HCN), ethane (C2H6), acetylene (C2H2), methanol (CH3OH), and formic acid (HCOOH). We provide their total columns time-series obtained from the measurements during August–October 2004, May–October 2007, and May 2009–December 2010. We show that biomass burning explains a large part of the observed seasonal and interannual variability of the chemical species. The correlations between the daily mean total columns of each of the species and those of CO, also measured with our FTIR spectrometer at Reunion Island, are very good from August to November (R ≥ 0.86). This allows us to derive, for that period, the following enhancement ratios with respect to CO: 0.0047, 0.0078, 0.0020, 0.012, and 0.0046 for HCN, C2H6, C2H2, CH3OH, and HCOOH, respectively. The HCN ground-based data are compared to the chemical transport model GEOS-Chem, while the data for the other species are compared to the IMAGESv2 model. We show that using the HCN/CO ratio derived from our measurements (0.0047) in GEOS-Chem reduces the underestimation of the modeled HCN columns compared with the FTIR measurements. The comparisons between IMAGESv2 and the long-lived species C2H6 and C2H2 indicate that the biomass burning emissions used in the model (from the GFED3 inventory) are probably underestimated in the late September–October period for all years of measurements, and especially in 2004. The comparisons with the short-lived species, CH3OH and HCOOH, with lifetimes of around 5 days, suggest that the emission underestimation in late September–October 2004, occurs more specifically in the Southeastern Africa-Madagascar region. The very good correlation of CH3OH and HCOOH with CO suggests that, despite the dominance of the biogenic source of these compounds on the global scale, biomass burning is their major source at Reunion Island between August and November.