Constraining the budget of NOx and VOCs at a remote Tropical island using multi-platform observations and WRF-Chem model simulations

Abstract. Volatile organic compounds (VOCs) act as precursors to ozone and secondary organic aerosols, which have significant health and environmental impacts, and they can reduce the atmospheric oxidative capacity. However, their budget remains poorly quantified, especially over remote areas such as the Tropical oceans. Here, we present high-resolution simulations of atmospheric composition over Réunion Island, in the Indian Ocean, using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). The coexistence and spatial heterogeneity of anthropogenic and biogenic emission sources in this region present a valuable but challenging test of the model performance. The WRF-Chem model is evaluated against several observational datasets, including Proton Transfer Reaction Mass Spectrometry (PTR-MS) measurements of VOCs and oxygenated VOCs (OVOCs) at the Maïdo Observatory (2160 m a.s.l.) in January and July 2019. While the primary goal of our study is a better understanding of (O)VOC budget at remote Tropical latitudes, important model refinements are made to improve the model performance, including the implementation of high-resolution anthropogenic and biogenic isoprene emissions, updates to the chemical mechanism, and adjustments to the boundary conditions. These refinements are supported by comparisons with PTR-MS data as well as with meteorological measurements at Maïdo, in situ NOx and O₃ measurements from the air quality Atmo-Réunion network, Fourier Transform Infrared Spectroscopy (FTIR) measurements of O₃, CO, ethane and several OVOCs, also at Maïdo, and satellite retrievals from the TROPOspheric Monitoring Instrument (TROPOMI).

TROPOMI NO₂ data suggests that anthropogenic emissions, particularly from power plants near Le Port, dominate NOx levels over the island. Both TROPOMI and in situ surface NO₂ comparisons are used to adjust the power plant emissions at Le Port. Surface ozone concentrations are overestimated by ~6 ppbv on average, likely due to the neglect of halogen chemistry in the model. Whereas modelled NO₂ over oceans is too low in summer when the lightning source is turned off, the inclusion of this source results in model overestimations corroborated by comparisons with upper tropospheric NO₂ mixing ratios derived from TROPOMI using the cloud-slicing technique (Marais et al., 2021). The model generally succeeds in reproducing the PTR-MS isoprene and its oxidation products (Iox), except for a moderate underestimation (~30 %) of noontime isoprene concentration and for modelled concentration peaks near dawn and dusk, that are not seen in the observations. The ratio of Iox to isoprene (0.8 at noon in January) is fairly well reproduced by the model. The methanol and monoterpenes observations both suggest overestimations of their biogenic emissions, by factors of about 2 and 5, respectively. Acetaldehyde anthropogenic emissions are likely strongly overestimated, due to the lumping of higher aldehydes into this compound. Without this lumping, the modelled acetaldehyde would be underestimated by almost one order of magnitude, suggesting the existence of a large missing source, likely photochemical. The comparisons suggest the existence of a biogenic source of MEK equivalent to about 3 % of isoprene emissions, likely associated with the dry deposition and conversion to MEK of key isoprene oxidation products. A strong model underestimation of the PTR-MS signal at mass 61 is also found, by a factor of 3–5 during daytime, consistent with previously reported missing sources of acetic and peracetic acid.