We present observations of a gravity wave (GW) event over the Netherlands and Belgium in the night and early morning of 30 June 2022. This event was captured by instruments that are part of the meteorological observational network of The Netherlands and Belgium, including the operational weather radars and automatic lidar ceilometers, and by the (tower) in situ and Doppler lidar at the Cabauw atmospheric research station. Together they provide 3D information on the GW event, with detailed information on the vertical profiles from Cabauw, and the spatial extent and evolution from the observational synoptic network. Two distinct GW trains were observed, both interpreted as ducted GW that are trapped in the lowermost 500 m of the stable nocturnal boundary layer. The GWs showed large vertical velocity amplitudes up to 3 m/s, resulting in strong modulations of wind, temperature, humidity and pressure. By comparing our observations to the Copernicus European Regional ReAnalysis (CERRA) wind and temperature fields, we interpret that these GWs were generated by convergence lines over Belgium and the North Sea that was linked to an occlusion front.

Atmospheric composition plays an important role in present and near-future climate change. Airborne particles exert direct and indirect radiative impacts and can serve as cloud condensation and ice nuclei, having therefore a strong influence on cloud formation and precipitation. Furthermore, a detailed understanding of present-day atmospheric transport pathways of particles from source to deposition in Antarctica remains essential.Since 2010, the aerosol total number and size distribution, aerosol absorption coefficient and mass concentration of light-absorbing aerosols and the aerosol total scattering coefficient have been monitored at the Belgian research station Princess Elisabeth Antarctica (PEA). The station is situated in Dronning Maud Land, East Antarctica (71.95° S, 23.35° E, 1390 m asl). Besides these instruments, a cloud condensation nuclei counter was operated during three austral summers. Meteorological data come from an automatic weather station. In this work, we investigate the climatology of the particle properties with respect to the air mass origin. To that end, we used the FLEXTRA trajectory model to investigate transport pathways into Antarctica. The model was driven with ECMWF ERA-5 meteorological fields. 10-days 3D backward trajectories, starting from PEA, were calculated for the period 01/01/2010 to 31/12/2020, in 3-hour-intervals. A k-means cluster analysis has been done based on latitude, longitude and altitude, resulting in four clusters of air mass origin.We will present results for the climatology of particle properties and the air mass origin. In addition, the backward trajectories have been combined with measured atmospheric particle properties and parameters like potential vorticity and exposure to sunshine duration, showing the distribution of the measured atmospheric particle properties between and within the air mass origin clusters. Some distinct features could be seen in the air mass origin clustering. Source regions from South America, Southern Africa and Australia, New Zealand were limited and the Southern Ocean was a main source region, as was the Antarctic continent itself. For each season, the dominating cluster represented mainly air masses of Antarctic continental origin with a large influence of upper tropospheric air. We will show further results of our analysis on air mass origin and atmospheric and particle properties, with respect to differentiations between seasons, clusters, continental and maritime origin and source altitude compartments.

The Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) is employed as an intercomparison tool for validating TROPOspheric Monitoring Instrument (TROPOMI) satellite NO2 retrievals against high-resolution Airborne Prism EXperiment (APEX) remote sensing observations performed in June 2019 in the region of Antwerp, a major hotspot of NO2 pollution in Europe. The model is first evaluated using meteorological and chemical observations in this area. Sensitivity simulations varying the model planetary layer boundary (PBL) parameterization were conducted for a 3 d period in June 2019, indicating a generally good performance of most parameterizations against meteorological data (namely ceilometer, surface meteorology, and balloon measurements), except for a moderate overestimation (∼ 1 m s−1) of near-surface wind speed. On average, all but one of the PBL schemes reproduce the surface NO2 measurements at stations of the Belgian Interregional Environmental Agency fairly well, although surface NO2 is generally underestimated during the day (between −4.3 % and −25.1 % on average) and overestimated at night (8.2 %–77.3 %). This discrepancy in the diurnal evolution arises despite (1) implementing a detailed representation of the diurnal cycle of emissions (Crippa et al., 2020) and (2) correcting the modeled concentrations to account for measurement interferences due to NOy reservoir species, which increases NO2 concentrations by about 20 % during the day. The model is further evaluated by comparing a 15 d simulation with surface NO2, NO, CO, and O3 data in the Antwerp region. The modeled daytime NO2 concentrations are more negatively biased during weekdays than during weekends, indicating a misrepresentation of the weekly temporal profile applied to the emissions obtained from Crippa et al. (2020). Using a mass balance approach, we determined a new weekly profile of NOx emissions, leading to a homogenization of the relative bias among the different weekdays. The ratio of weekend to weekday emissions is significantly lower in this updated profile (0.6) than in the profile based on Crippa et al. (2020; 0.84). Comparisons with remote sensing observations generally show a good reproduction of the spatial patterns of NO2 columns by the model. The model underestimated both APEX (by ca. −37 %) and TROPOMI columns (ca. −25 %) on 27 June, whereas no significant bias is found on 29 June. The two datasets are intercompared by using the model as an intermediate platform to account for differences in vertical sensitivity through the application of averaging kernels. The derived bias of TROPOMI v1.3.1 NO2 with respect to APEX is about −10 % for columns between (6–12) × 1015 molec. cm−2. The obtained bias for TROPOMI v1.3.1 increases with the NO2 column, following CAPEX=1.217Cv1.3-0.783 × 1015 molec. cm−2, in line with previous validation campaigns. The bias is slightly lower for the reprocessed TROPOMI v2.3.1, with CAPEX=1.055CPAL-0.437 × 1015 molec. cm−2 (PAL). Finally, a mass balance approach was used to perform a crude inversion of NOx emissions based on 15 d averaged TROPOMI columns. The emission correction is conducted only in regions with high columns and high sensitivity to emission changes in order to minimize the errors due to wind transport. The results suggest that emissions increase over Brussels–Antwerp (+20 %), the Ruhr Valley (13 %), and especially Paris (+39 %), and emissions decrease above a cluster of power plants in western Germany.

Starting in 1969 and comprising three launches a week, the Uccle (Brussels, Belgium) ozonesonde dataset is one of longest and densest in the world. Moreover, as the only major change was the switch from Brewer-Mast (BM) to electrochemical concentration cell (ECC) ozonesonde types in 1997 (when the emissions of ozone-depleting substances peaked), the Uccle time series is very homogenous. In this paper, we briefly describe the efforts that were undertaken during the first 3 decades of the 50 years of ozonesonde observations to guarantee the homogeneity between ascent and descent profiles, under changing environmental conditions (e.g. SO2), and between the different ozonesonde types. This paper focuses on the 50-year-long Uccle ozonesonde dataset and aims to demonstrate its past, present, and future relevance to ozone research in two application areas: (i) the assessment of the temporal evolution of ozone from the surface to the (middle) stratosphere, and (ii) as the backbone for validation and stability analysis of both stratospheric and tropospheric satellite ozone retrievals. Using the Long-term Ozone Trends and Uncertainties in the Stratosphere (LOTUS) multiple linear regression model (SPARC/IO3C/GAW, 2019), we found that the stratospheric ozone concentrations at Uccle have declined at a significant rate of around 2 % per decade since 1969, which is also rather consistent over the different stratospheric levels. This overall decrease can mainly be assigned to the 1969–1996 period with a rather consistent rate of decrease of around −4 % per decade. Since 2000, a recovery of between +1 % per decade and +3 % per decade of the stratospheric ozone levels above Uccle has been observed, although it is not significant and is not seen for the upper stratospheric levels measured by ozonesondes. Throughout the entire free troposphere, a very consistent increase in the ozone concentrations of 2 % per decade to 3 % per decade has been measured since both 1969 and 1995, with the trend since 1995 being in almost perfect agreement with the trends derived from the In-service Aircraft for a Global Observing System (IAGOS) ascent/descent profiles at Frankfurt. As the number of tropopause folding events in the Uccle time series has increased significantly over time, increased stratosphere-to-troposphere transport of recovering stratospheric ozone might partly explain these increasing tropospheric ozone concentrations, despite the levelling-off of (tropospheric) ozone precursor emissions and notwithstanding the continued increase in mean surface ozone concentrations. Furthermore, we illustrate the crucial role of ozonesonde measurements for the validation of satellite ozone profile retrievals. With the operational validation of the Global Ozone Monitoring Experiment-2 (GOME-2), we show how the Uccle dataset can be used to evaluate the performance of a degradation correction for the MetOp-A/GOME-2 UV (ultraviolet) sensors. In another example, we illustrate that the Microwave Limb Sounder (MLS) overpass ozone profiles in the stratosphere agree within ±5 % with the Uccle ozone profiles between 10 and 70 hPa. Another instrument on the same Aura satellite platform, the Tropospheric Emission Spectrometer (TES), is generally positively biased with respect to the Uccle ozonesondes in the troposphere by up to ∼ 10 ppbv, corresponding to relative differences of up to ∼ 15 %. Using the Uccle ozonesonde time series as a reference, we also demonstrate that the temporal stability of those last two satellite retrievals is excellent.

The purpose of the EUNADICS-AV (European Natural Airborne Disaster Information and Coordination System for Aviation) prototype early warning system (EWS) is to develop the combined use of harmonised data products from satellite, ground-based and in situ instruments to produce alerts of airborne hazards (volcanic, dust, smoke and radionuclide clouds), satisfying the requirement of aviation air traffic management (ATM) stakeholders (https://cordis.europa.eu/project/id/723986, last access: 5 November 2021). The alert products developed by the EUNADICS-AV EWS, i.e. near-real-time (NRT) observations, email notifications and netCDF (Network Common Data Form) alert data products (called NCAP files), have shown significant interest in using selective detection of natural airborne hazards from polar-orbiting satellites. The combination of several sensors inside a single global system demonstrates the advantage of using a triggered approach to obtain selective detection from observations, which cannot initially discriminate the different aerosol types. Satellite products from hyperspectral ultraviolet–visible (UV–vis) and infrared (IR) sensors (e.g. TROPOMI – TROPOspheric Monitoring Instrument – and IASI – Infrared Atmospheric Sounding Interferometer) and a broadband geostationary imager (Spinning Enhanced Visible and InfraRed Imager; SEVIRI) and retrievals from ground-based networks (e.g. EARLINET – European Aerosol Research Lidar Network, E-PROFILE and the regional network from volcano observatories) are combined by our system to create tailored alert products (e.g. selective ash detection, SO2 column and plume height, dust cloud, and smoke from wildfires). A total of 23 different alert products are implemented, using 1 geostationary and 13 polar-orbiting satellite platforms, 3 external existing service, and 2 EU and 2 regional ground-based networks. This allows for the identification and the tracking of extreme events. The EUNADICS-AV EWS has also shown the need to implement a future relay of radiological data (gamma dose rate and radionuclides concentrations in ground-level air) in the case of a nuclear accident. This highlights the interest of operating early warnings with the use of a homogenised dataset. For the four types of airborne hazard, the EUNADICS-AV EWS has demonstrated its capability to provide NRT alert data products to trigger data assimilation and dispersion modelling providing forecasts and inverse modelling for source term estimate. Not all of our alert data products (NCAP files) are publicly disseminated. Access to our alert products is currently restricted to key users (i.e. Volcanic Ash Advisory Centres, national meteorological services, the World Meteorological Organization, governments, volcano observatories and research collaborators), as these are considered pre-decisional products. On the other hand, thanks to the EUNADICS-AV–SACS (Support to Aviation Control Service) web interface (https://sacs.aeronomie.be, last access: 5 November 2021), the main part of the satellite observations used by the EUNADICS-AV EWS is shown in NRT, with public email notification of volcanic emission and delivery of tailored images and NCAP files. All of the ATM stakeholders (e.g. pilots, airlines and passengers) can access these alert products through this free channel.

The particle composition is important for air quality studies in the urban atmosphere. In particular, particles smaller than 1 μm are increasingly in the focus with respect to human health because of they are inhalable deeply into the human lung. A relevant part of such (ultra-)fine particles are light-absorbing particles. Important sources for such particles in the cities and residential areas are emissions from traffic (mainly soot), but also emissions from wood-burning stoves in private households. The relative contributions of these sources to the atmospheric particle load are important to know in order to be able to reduce hazardous emissions by specific measures. The Royal Meteorological Institute of Belgium (RMI) gathers ambient aerosol data in Brussels with a 7-wavelengths aethalometer (mass concentration and absorption coefficient of light-absorbing aerosol). The measurement site is located in the sub-urban, rather residential southern part of Brussels and measurements are representative for the urban background. At the same site, the boundary layer height and atmospheric stability parameter are derived from a co-located ceilometer and eddy-covariance system. Further, RMI uses the multi-scale chemical transport model CHIMERE, coupled to the high resolution regional numerical weather prediction limited area model ALARO. In this context, CHIMERE will be used to evaluate emission scenarios. The wavelength dependency of the measured aerosol parameters revealed distinct variations of aerosol composition, both on a daily, weekly, and seasonal scale. The Absorption Angstrom Exponent (AAE; spectral dependency of the absorption coefficient) revealed values around 1.3 during winter and around 1.0 during summer months. As fresh soot has a flat spectrum, (i.e., AAE around 1.0), higher AAE values during winter meant that other light-absorbing compounds, which absorb stronger in the UV, increased in importance (e.g., wood burning aerosol). The multi-wavelengths data is exploited to derive the relative contribution of fresh soot (thus traffic emissions) and other sources to the amount of light-absorbing aerosol. The influence of boundary layer height and atmospheric stability on the aerosol data will be presented.

The Royal Meteorological Institute of Belgium (RMI) gathered ambient aerosol data in Brussels with a 7-wavelengths aethalometer and a 3 wavelengths integrating nephelometer. The amount of light-absorbing particles showed a clear daily and weekly cycle, with a sharp peak in the morning rush hour time and a broader peak in the evening. During weekends, the rush hour peak diminished. The spectral dependency of the absorption coefficient revealed peak contributions of traffic emissions to the amount of light-absorbing particles of up to 90%. Other sources (like wood burning from households) showed peak contributions of up to 35%. These percentages showed in addition a clear daily, weekly and seasonal cycle, with higher contributions of these other sources during night time, weekends and summer.

Ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements of aerosols and tropospheric nitrogen dioxide (NO2) were carried out in Uccle (50.8∘ N, 4.35∘ E), Brussels, during 1 year from March 2018 until March 2019. The instrument was operated in both the UV and visible wavelength ranges in a dual-scan configuration consisting of two sub-modes: (1) an elevation scan in a fixed viewing azimuthal direction (the so-called main azimuthal direction) pointing to the northeast and (2) an azimuthal scan in a fixed low elevation angle (2∘). By applying a vertical profile inversion algorithm in the main azimuthal direction and a parameterization technique in the other azimuthal directions, near-surface NO2 volume mixing ratios (VMRs) and vertical column densities (VCDs) were retrieved in 10 different azimuthal directions. The dual-scan MAX-DOAS dataset allows for partly resolving the horizontal distribution of NO2 around the measurement site and studying its seasonal variations. Furthermore, we show that measuring the tropospheric NO2 VCDs in different azimuthal directions improves the spatial colocation with measurements from the Sentinel-5 Precursor (S5P), leading to a reduction of the spread in validation results. By using NO2 vertical profile information derived from the MAX-DOAS measurements, we also resolve a systematic underestimation in S5P NO2 data due to the use of inadequate a priori NO2 profile shape data in the satellite retrieval.

Abstract. Ground-based Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) measurements of aerosols and tropospheric nitrogen dioxide (NO2) were carried out in Uccle (50.8° N, 4.35° E) Brussels, during one year from March 2018 until March 2019. The instrument was operated in both UV and visible (Vis) wavelength ranges in a dual-scan configuration consisting of two sub-modes: (1) an elevation scan in a fixed viewing azimuthal direction (the so-called main azimuthal direction) pointing to the Northeast and (2) an azimuthal scan in a fixed low elevation angle (2°). By applying a vertical profile inversion algorithm in the main azimuthal direction and a parameterization technique in the other azimuthal directions, near-surface NO2 concentrations (VMRs) and vertical column densities (VCDs) were retrieved in ten different azimuthal directions. The dual-scan MAX-DOAS dataset allows partly resolving the horizontal distribution of NO2 around the measurement site and studying its seasonal variations. Furthermore, we show that measuring the tropospheric NO2 VCDs in different azimuthal directions improves the spatial colocation with measurements from the Sentinel-5 Precursor (S5P), leading to a reduction of the spread in validation results. By using NO2 vertical profile information derived from the MAX-DOAS measurements, we also resolve a systematic underestimation in S5P NO2 data due to the use of inadequate a-priori NO2 profile shape data in the satellite retrieval.

Ground-based Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) measurements of aerosols and tropospheric nitrogen dioxide (NO2) were carried out in Uccle (50.8° N, 4.35° E) Brussels, during one year from March 2018 until March 2019. The instrument was operated in both UV and visible (Vis) wavelength ranges in a dual-scan configuration consisting of two sub-modes: (1) an elevation scan in a fixed viewing azimuthal direction (the so-called main azimuthal direction) pointing to the Northeast and (2) an azimuthal scan in a fixed low elevation angle (2°). By applying a vertical profile inversion algorithm in the main azimuthal direction and a parameterization technique in the other azimuthal directions, near-surface NO2 concentrations (VMRs) and vertical column densities (VCDs) were retrieved in ten different azimuthal directions. The dual-scan MAX-DOAS dataset allows partly resolving the horizontal distribution of NO2 around the measurement site and studying its seasonal variations. Furthermore, we show that measuring the tropospheric NO2 VCDs in different azimuthal directions improves the spatial colocation with measurements from the Sentinel-5 Precursor (S5P), leading to a reduction of the spread in validation results. By using NO2 vertical profile information derived from the MAX-DOAS measurements, we also resolve a systematic underestimation in S5P NO2 data due to the use of inadequate a-priori NO2 profile shape data in the satellite retrieval.

Since 2010, aerosol number, size and optical have been monitored at the Belgian Antarctic research station Princess Elisabeth. Atmospheric particles have been collected during the last three austral summers for chemical analysis by high-volume sampling. In addition, passive samplers have been installed at seven locations along a 250 km N-S transect from the coast to the Antarctic plateau, and the collected samples have been analysed on organics, VOCs and inorganics, respectively. Results will be presented for the multi-year seasonality of aerosol phyiscal properties. Further, first results of the chemical analyses for organics, inorganics and VOCs will be presented. E.g., the mineralogy of some samples point to distinct sources in South-America or southern Africa, respectivley. Simulations with the FLEXPART model will show how particles from these lower latitudes could reach Antarctica. FLEXTRA and FLEXPART results will also be shown for the air mass origin of specific events.

For three austral summer seasons (2013–2016, each from December to February) aerosol particles arriving at the Belgian Antarctic research station Princess Elisabeth (PE) in Dronning Maud Land in East Antarctica were characterized. This included number concentrations of total aerosol particles (NCN) and cloud condensation nuclei (NCCN), the particle number size distribution (PNSD), the aerosol particle hygroscopicity, and the influence of the air mass origin on NCN and NCCN. In general NCN was found to range from 40 to 6700 cm−3, with a median of 333 cm−3, while NCCN was found to cover a range between less than 10 and 1300 cm−3 for supersaturations (SSs) between 0.1 % and 0.7 %. It is shown that the aerosol is dominated by the Aitken mode, being characterized by a significant amount of small, and therefore likely secondarily formed, aerosol particles, with 94 % and 36 % of the aerosol particles smaller than 90 and ≈35 nm, respectively. Measurements of the basic meteorological parameters as well as the history of the air masses arriving at the measurement station indicate that the station is influenced by both marine air masses originating from the Southern Ocean and coastal areas around Antarctica (marine events – MEs) and continental air masses (continental events – CEs). CEs, which were defined as instances when the air masses spent at least 90 % of the time over the Antarctic continent during the last 10 days prior to arrival at the measurements station, occurred during 61 % of the time during which measurements were done. CEs came along with rather constant NCN and NCCN values, which we denote as Antarctic continental background concentrations. MEs, however, cause large fluctuations in NCN and NCCN, with low concentrations likely caused by scavenging due to precipitation and high concentrations likely originating from new particle formation (NPF) based on marine precursors. The application of HYSPLIT back trajectories in form of the potential source contribution function (PSCF) analysis indicate that the region of the Southern Ocean is a potential source of Aitken mode particles. On the basis of PNSDs, together with NCCN measured at an SS of 0.1 %, median values for the critical diameter for cloud droplet activation and the aerosol particle hygroscopicity parameter κ were determined to be 110 nm and 1, respectively. For particles larger than ≈110 nm the Southern Ocean together with parts of the Antarctic ice shelf regions were found to be potential source regions. While the former may contribute sea spray particles directly, the contribution of the latter may be due to the emission of sea salt aerosol particles, released from snow particles from surface snow layers, e.g., during periods of high wind speed, leading to drifting or blowing snow. The region of the Antarctic inland plateau, however, was not found to feature a significant source region for aerosol particles in general or for cloud condensation nuclei measured at the PE station in the austral summer.

Local surface mass balance (SMB) measurements are crucial for understanding changes in the total mass of the Antarctic Ice Sheet, including its contribution to sea level rise. Despite continuous attempts to decipher mechanisms controlling the local and regional SMB, a clear understanding of the separate components is still lacking, while snowfall measurements are almost absent. In this study, the different terms of the SMB are quantified at the Princess Elisabeth (PE) station in Dronning Maud Land, East Antarctica. Furthermore, the relationship between snowfall and accumulation at the surface is investigated. To achieve this, a unique collocated set of ground-based and in situ remote sensing instrumentation (Micro Rain Radar, ceilometer, automatic weather station, among others) was set up and operated for a time period of 37 months. Snowfall originates mainly from moist and warm air advected from lower latitudes associated with cyclone activity. However, snowfall events are not always associated with accumulation. During 38 % of the observed snowfall cases, the freshly fallen snow is ablated by the wind during the course of the event. Generally, snow storms of longer duration and larger spatial extent have a higher chance of resulting in accumulation on a local scale, while shorter events usually result in ablation (on average 17 and 12 h respectively). A large part of the accumulation at the station takes place when preceding snowfall events were occurring in synoptic upstream areas. This fresh snow is easily picked up and transported in shallow drifting snow layers over tens of kilometres, even when wind speeds are relatively low (< 7 ms−1). Ablation events are mainly related to katabatic winds originating from the Antarctic plateau and the mountain ranges in the south. These dry winds are able to remove snow and lead to a decrease in the local SMB. This work highlights that the local SMB is strongly influenced by synoptic upstream conditions.

With the remnants of hurricane Ophelia reaching the British Isles, a rare pollution event crossed Europe in October 2017, with transport of dust from the Sahara and of thick smoke from disastrous wildfires in northern Portugal/Spain. On October 15-16, both dust and smoke were rapidly advected by Ophelia-energized strong southerly winds into the Bay of Biscay, north-western France and across the UK. Two unusual events were intersecting: a tropical storm and high intensity late-season fires. We now benefit from a range of observations above Europe allowing us to track the transport of such a pollution plume. We used an integrated observational approach for analysing this event in near-real-time, combining both observations and model analysis. Using near real time radiance observations from the IASI thermal infrared remote sensors onboard the two Metop satellites, focused maps of dust, carbon monoxide and ammonia (two proxies for fire plumes) could be derived. The European lidar/ceilometer network established the overpass time for both dust and the smoke plumes at different locations, along with altitude information. The IAGOS network on passenger aircraft also reported large peaks of CO during departures and landings, with a measured CO concentration matching well with the IASI observations. Finally, the smoke and dust particle-size distributions are compared with aerosol forecasts and analyses produced by the Copernicus Atmosphere Monitoring Service (CAMS) using the ECMWF IFS, constrained by MODIS and PMAp satellite observations of aerosol optical depth, to inform on the different aerosol sources.