Impact of forest fires, biogenic emissions and high temperatures on the elevated Eastern Mediterranean ozone levels during the hot summer of 2007

International audienceThe hot summer of 2007 in southeast Europe has been studied using two regional atmospheric chemistry models; WRF-Chem and EMEP MSC-W. The region was struck by three heat waves and a number of forest fire episodes, greatly affecting air pollution levels. We have focused on ozone and its precursors using state-of-the-art inventories for anthropogenic, biogenic and forest fire emissions. The models have been evaluated against measurement data, and processes leading to ozone formation have been quantified. Heat wave episodes are projected to occur more frequently in a future climate, and therefore this study also makes a contribution to climate change impact research. The plume from the Greek forest fires in August 2007 is clearly seen in satellite observations of CO and NO2 columns, showing extreme levels of CO in and downwindof the fires. Model simulations reflect the location and influence of the fires relatively well, but the modelled magnitude of CO in the plume core is too low. Most likely, this is caused by underestimation of CO in the emission inventories, suggesting that the CO/NOx ratios of fire emissions should be re-assessed. Moreover, higher maximum values are seen in WRF-Chem than in EMEP MSC-W, presumably due to differences in plume rise altitudes as the first model emits a larger fraction of the fire emissions in the lowermost model layer. The model results are also in fairly good agreement with surface ozone measurements. Biogenic VOC emissions reacting with anthropogenic NOx emissions are calculated to contribute significantly to the levels of ozone in the region, but the magnitude and geographical distribution depend strongly on the model and biogenic emission module used. During the July and August heat waves, ozone levels increased substantially due to a combination of forest fire emissions and the effect of high temperatures. We found that the largest temperature impact on ozone was through the temperature dependence of the biogenic emissions, closely followed by the effect of reduced dry deposiion caused by closing of the plants' stomata at very high temperatures. The impact of high temperatures on the ozone chemistry was much lower. The results suggest that forest fire emissions, and the temperature effect on biogenic emissions and dry deposition, will potentially lead to substantial ozone increases in a warmer climate

Inter-comparison of MAX-DOAS measurements of tropospheric HONO slant column densities and vertical profiles during the CINDI-2 campaign

We present the inter-comparison of delta slant column densities (SCDs) and vertical profiles of nitrous acid (HONO) derived from measurements of different multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments and using different inversion algorithms during the Second Cabauw Inter-comparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2) in September 2016 at Cabauw, the Netherlands (51.97∘ N, 4.93∘ E). The HONO vertical profiles, vertical column densities (VCDs), and near-surface volume mixing ratios are compared between different MAX-DOAS instruments and profile inversion algorithms for the first time. Systematic and random discrepancies of the HONO results are derived from the comparisons of all data sets against their median values. Systematic discrepancies of HONO delta SCDs are observed in the range of ±0.3×1015 molec. cm−2, which is half of the typical random discrepancy of 0.6×1015 molec. cm−2. For a typical high HONO delta SCD of 2×1015 molec. cm−2, the relative systematic and random discrepancies are about 15 % and 30 %, respectively. The inter-comparison of HONO profiles shows that both systematic and random discrepancies of HONO VCDs and near-surface volume mixing ratios (VMRs) are mostly in the range of ∼±0.5×1014 molec. cm−2 and ∼±0.1 ppb (typically ∼20 %). Further we find that the discrepancies of the retrieved HONO profiles are dominated by discrepancies of the HONO delta SCDs. The profile retrievals only contribute to the discrepancies of the HONO profiles by ∼5 %. However, some data sets with substantially larger discrepancies than the typical values indicate that inappropriate implementations of profile inversion algorithms and configurations of radiative transfer models in the profile retrievals can also be an important uncertainty source. In addition, estimations of measurement uncertainties of HONO dSCDs, which can significantly impact profile retrievals using the optimal estimation method, need to consider not only DOAS fit errors, but also atmospheric variability, especially for an instrument with a DOAS fit error lower than ∼3×1014 molec. cm−2. The MAX-DOAS results during the CINDI-2 campaign indicate that the peak HONO levels (e.g. near-surface VMRs of ∼0.4 ppb) often appeared in the early morning and below 0.2 km. The near-surface VMRs retrieved from the MAX-DOAS observations are compared with those measured using a co-located long-path DOAS instrument. The systematic differences are smaller than 0.15 and 0.07 ppb during early morning and around noon, respectively. Since true HONO values at high altitudes are not known in the absence of real measurements, in order to evaluate the abilities of profile inversion algorithms to respond to different HONO profile shapes, we performed sensitivity studies using synthetic HONO delta SCDs simulated by a radiative transfer model with assumed HONO profiles. The tests indicate that the profile inversion algorithms based on the optimal estimation method with proper configurations can reproduce the different HONO profile shapes well. Therefore we conclude that the features of HONO accumulated near the surface derived from MAX-DOAS measurements are expected to represent the ambient HONO profiles well

Overview: On the transport and transformation of pollutants in the outflow of major population centres – observational data from the EMeRGe European intensive operational period in summer 2017

Megacities and other major population centres (MPCs) worldwide are major sources of air pollution, both locally as well as downwind. The overall assessment and prediction of the impact of MPC pollution on tropospheric chemistry are challenging. The present work provides an overview of the highlights of a major new contribution to the understanding of this issue based on the data and analysis of the EMeRGe (Effect of Megacities on the transport and transformation of pollutants on the Regional to Global scales) international project. EMeRGe focuses on atmospheric chemistry, dynamics, and transport of local and regional pollution originating in MPCs. Airborne measurements, taking advantage of the long range capabilities of the High Altitude and LOng Range Research Aircraft (HALO, https://www.halo-spp.de, last access: 22 March 2022), are a central part of the project. The synergistic use and consistent interpretation of observational data sets of different spatial and temporal resolution (e.g. from ground-based networks, airborne campaigns, and satellite measurements) supported by modelling within EMeRGe provide unique insight to test the current understanding of MPC pollution outflows. In order to obtain an adequate set of measurements at different spatial scales, two field experiments were positioned in time and space to contrast situations when the photochemical transformation of plumes emerging from MPCs is large. These experiments were conducted in summer 2017 over Europe and in the inter-monsoon period over Asia in spring 2018. The intensive observational periods (IOPs) involved HALO airborne measurements of ozone and its precursors, volatile organic compounds, aerosol particles, and related species as well as coordinated ground-based ancillary observations at different sites. Perfluorocarbon (PFC) tracer releases and model forecasts supported the flight planning, the identification of pollution plumes, and the analysis of chemical transformations during transport. This paper describes the experimental deployment and scientific questions of the IOP in Europe. The MPC targets – London (United Kingdom; UK), the Benelux/Ruhr area (Belgium, the Netherlands, Luxembourg and Germany), Paris (France), Rome and the Po Valley (Italy), and Madrid and Barcelona (Spain) – were investigated during seven HALO research flights with an aircraft base in Germany for a total of 53 flight hours. An in-flight comparison of HALO with the collaborating UK-airborne platform Facility for Airborne Atmospheric Measurements (FAAM) took place to assure accuracy and comparability of the instrumentation on board. Overall, EMeRGe unites measurements of near- and far-field emissions and hence deals with complex air masses of local and distant sources. Regional transport of several European MPC outflows was successfully identified and measured. Chemical processing of the MPC emissions was inferred from airborne observations of primary and secondary pollutants and the ratios between species having different chemical lifetimes. Photochemical processing of aerosol and secondary formation or organic acids was evident during the transport of MPC plumes. Urban plumes mix efficiently with natural sources as mineral dust and with biomass burning emissions from vegetation and forest fires. This confirms the importance of wildland fire emissions in Europe and indicates an important but discontinuous contribution to the European emission budget that might be of relevance in the design of efficient mitigation strategies. The present work provides an overview of the most salient results in the European context, with these being addressed in more detail within additional dedicated EMeRGe studies. The deployment and results obtained in Asia will be the subject of separate publications

SciPy 1.0: fundamental algorithms for scientific computing in Python.

SciPy is an open-source scientific computing library for the Python programming language. Since its initial release in 2001, SciPy has become a de facto standard for leveraging scientific algorithms in Python, with over 600 unique code contributors, thousands of dependent packages, over 100,000 dependent repositories and millions of downloads per year. In this work, we provide an overview of the capabilities and development practices of SciPy 1.0 and highlight some recent technical developments

Differences in satellite-derived NOx emission factors between Eurasian and North American boreal forest fires

Current fire emission inventories apply universal emission factors (EFs) for the calculation of NOx emissions over large biomes such as boreal forest. However, recent satellite-based studies over tropical and subtropical regions have indicated spatio-temporal variations in EFs within specific biomes. In this study, satellite measurements of tropospheric NO2 vertical columns (TVC NO2) from the GOME-2 instrument and fire radiative power (FRP) from MODIS are used for the estimation of fire emission rates (FERs) of NOx over Eurasian and North American boreal forests. The retrieval of TVC NO2 is based on a stratospheric correction using simulated stratospheric NO2 instead of applying the reference sector method, which was used in a previous study. The model approach is more suitable for boreal latitudes. TVC NO2 and FRP are spatially aggregated to a 1 degree by 1 degree horizontal resolution and temporally averaged to monthly values. The conversion of the satellite-derived tropospheric NO2 columns into production rates of NOx from fire (Pf) is based on the NO2/NOx ratio as obtained from the MACC reanalysis data set and an assumed lifetime of NOx. A global land cover map is used to define boreal forests across these two regions in order to evaluate the FERs of NOx for this biome. The FERs of NOx, which are derived from the gradients of the linear relationship between Pf and FRP, are more than 30% lower for North American than for Eurasian boreal forest fires. We speculate that these discrepancies are mainly related to the variable nitrogen content in plant tissues, which is higher in deciduous forests dominating large parts in Eurasia. In order to compare the obtained values with EFs found in the literature, the FERs are converted into EFs. The satellite-based EFs of NOx are estimated at 0.83 and 0.61 g/kg for Eurasian and North American boreal forests, respectively, which is in good agreement with the value found in a recent emission factor compilation. However, recent fire emission inventories are based on EFs of NOx that are 3-5 times larger, which indicates that there are still large uncertainties in estimates of NOx from biomass burning, especially on the regional scale

Changes in atmospheric aerosol loading retrieved from space-based measurements during the past decade

The role and potential management of short-lived atmospheric pollutants such as aerosols are currently a topic of scientific and public debates. Our limited knowledge of atmospheric aerosol and its influence on the Earth's radiation balance has a significant impact on the accuracy and error of current predictions of future climate change. In the last few years, there have been several accounts of the changes in atmospheric aerosol derived from satellite observations, but no study considering the uncertainty caused by different/limited temporal sampling of polar-orbiting satellites and cloud disturbance in the trend estimates of cloud-free aerosol optical thickness (AOT). This study presents an approach to minimize the uncertainties by use of weighted least-squares regression and multiple satellite-derived AOTs from the space-born instruments, MODIS (onboard Terra from 2000 to 2009 and Aqua form 2003 to 2008), MISR (Terra from 2000 to 2010), and SeaWiFS (OrbView-2 from 1998 to 2007) and thereby provides more convincing trend estimates for atmospheric aerosols during the past decade. The AOT decreases over western Europe (i.e., by up to about −40% from 2003 to 2008). In contrast, a statistically significant increase (about +34% in the same period) over eastern China is observed and can be attributed to the increase in both industrial output and Asian desert dust

Inter-comparison of MAX-DOAS measurements of tropospheric HONO slant column densities and vertical profiles during the CINDI-2 campaign

We present the inter-comparison of delta slant column densities (SCDs) and vertical profiles of nitrous acid (HONO) derived from measurements of different multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments and using different inversion algorithms during the Second Cabauw Inter-comparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2) in September 2016 at Cabauw, the Netherlands (51.97∘ N, 4.93∘ E). The HONO vertical profiles, vertical column densities (VCDs), and near-surface volume mixing ratios are compared between different MAX-DOAS instruments and profile inversion algorithms for the first time. Systematic and random discrepancies of the HONO results are derived from the comparisons of all data sets against their median values. Systematic discrepancies of HONO delta SCDs are observed in the range of ±0.3×1015 molec. cm−2, which is half of the typical random discrepancy of 0.6×1015 molec. cm−2. For a typical high HONO delta SCD of 2×1015 molec. cm−2, the relative systematic and random discrepancies are about 15 % and 30 %, respectively. The inter-comparison of HONO profiles shows that both systematic and random discrepancies of HONO VCDs and near-surface volume mixing ratios (VMRs) are mostly in the range of ∼±0.5×1014 molec. cm−2 and ∼±0.1 ppb (typically ∼20 %). Further we find that the discrepancies of the retrieved HONO profiles are dominated by discrepancies of the HONO delta SCDs. The profile retrievals only contribute to the discrepancies of the HONO profiles by ∼5 %. However, some data sets with substantially larger discrepancies than the typical values indicate that inappropriate implementations of profile inversion algorithms and configurations of radiative transfer models in the profile retrievals can also be an important uncertainty source. In addition, estimations of measurement uncertainties of HONO dSCDs, which can significantly impact profile retrievals using the optimal estimation method, need to consider not only DOAS fit errors, but also atmospheric variability, especially for an instrument with a DOAS fit error lower than ∼3×1014 molec. cm−2. The MAX-DOAS results during the CINDI-2 campaign indicate that the peak HONO levels (e.g. near-surface VMRs of ∼0.4 ppb) often appeared in the early morning and below 0.2 km. The near-surface VMRs retrieved from the MAX-DOAS observations are compared with those measured using a co-located long-path DOAS instrument. The systematic differences are smaller than 0.15 and 0.07 ppb during early morning and around noon, respectively. Since true HONO values at high altitudes are not known in the absence of real measurements, in order to evaluate the abilities of profile inversion algorithms to respond to different HONO profile shapes, we performed sensitivity studies using synthetic HONO delta SCDs simulated by a radiative transfer model with assumed HONO profiles. The tests indicate that the profile inversion algorithms based on the optimal estimation method with proper configurations can reproduce the different HONO profile shapes well. Therefore we conclude that the features of HONO accumulated near the surface derived from MAX-DOAS measurements are expected to represent the ambient HONO profiles well