Cleaner Skies during the COVID-19 Lockdown

: During spring 2020, the COVID-19 pandemic caused massive reductions in emissions from industry and ground and airborne transportation. To explore the resulting atmospheric composition changes, we conducted the BLUESKY campaign with two research aircraft and measured trace gases, aerosols, and cloud properties from the boundary layer to the lower stratosphere. From 16 May to 9 June 2020, we performed 20 flights in the early COVID-19 lockdown phase over Europe and the Atlantic Ocean. We found up to 50% reductions in boundary layer nitrogen dioxide concentrations in urban areas from GOME-2B satellite data, along with carbon monoxide reductions in the pollution hot spots. We measured 20%–70% reductions in total reactive nitrogen, carbon monoxide, and fine mode aerosol concentration in profiles over German cities compared to a 10-yr dataset from passenger aircraft. The total aerosol mass was significantly reduced below 5 km altitude, and the organic aerosol fraction also aloft, indicative of decreased organic precursor gas emissions. The reduced aerosol optical thickness caused a perceptible shift in sky color toward the blue part of the spectrum (hence BLUESKY) and increased shortwave radiation at the surface. We find that the 80% decline in air traffic led to substantial reductions in nitrogen oxides at cruise altitudes, in contrail cover, and in resulting radiative forcing. The light extinction and depolarization by cirrus were also reduced in regions with substantially decreased air traffic. General circulation–chemistry model simulations indicate good agreement with the measurements when applying a reduced emission scenario. The comprehensive BLUESKY dataset documents the major impact of anthropogenic emissions on the atmospheric composition.

studies emphasize the need for comprehensive atmospheric composition measurements from the boundary layer to the stratosphere in different parts of the world in order to determine and better understand atmospheric composition changes caused by human activities and distinguish the anthropogenic impact from natural factors.
Travel restrictions resulted in more than 80% reductions in air traffic worldwide during the early lockdown phase (Guevara et al. 2021), and aviation experienced significantly stronger reductions compared to other transport sectors. While air traffic generally recovered within a year in Asia and America, European air traffic lagged behind and showed a significant decrease throughout summer 2020 with a slight recovery toward the end of the year (see Fig. 2). One year after the initial lockdown, European air traffic was still reduced by about 30% compared to pre-COVID-19 levels.
Aircraft emit carbon dioxide (CO 2 ), nitrogen oxides (NO x ), water vapor (H 2 O), and aerosols. The aircraft emissions can modify cirrus clouds and lead to the formation of contrails in cold and humid areas at cruise altitudes (Lee et al. 2010). The recent, comprehensive assessment of air traffic effects on the atmosphere (Lee et al. 2021) shows that aviation up to 2018 contributed about 3.5% to the total anthropogenic effective radiative forcing, with about one-third coming from its CO 2 emissions and two-thirds resulting from the non-CO 2 effects. In fact, the major contributor to effective radiative forcing from aviation is caused by contrail cirrus (57%) (Burkhardt et al. 2018; Lee et al. 2021). The effects of aviation generated aerosol on natural clouds remains uncertain (Lee et al. 2021). Properties of contrail cirrus have been measured from aircraft (Heymsfield et al. 2010;Voigt et al. 2010Voigt et al. , 2017Voigt et al. , 2021Schumann et al. 2017;Bräuer et al. 2021a,b) and satellites (Minnis et al. 2013;Vázquez-Navarro et al. 2015). These observations were used to evaluate the Contrail Cirrus Prediction model (COCIP) (Schumann 2012;Schumann et al. 2017), with which the radiative forcing from contrail cirrus can be calculated. Reduced contrail cirrus optical thickness and radiative forcing caused by diminished air traffic over Europe has been calculated for a 9-month period in 2020 (Schumann et al. 2021a). Significant reductions in seasonal and regional effective radiative forcing from contrail cirrus caused by reduced air traffic emissions during the lockdown 2020 were also found by Gettelman et al. (2021), while the annual mean effective radiative forcing was less affected. Also, the difficulty in separating the impact of meteorology from the reduced aircraft emissions impact on cloud properties induces considerable uncertainties in the calculations (Schumann et al. 2021b). Quaas et al. (2021) used satellite-based cloud retrievals in regions more and less affected by air traffic to derive changes in the aviation impact on cirrus clouds in 2020. Recent studies (Urbanek et al. 2018;Li and Groß 2021) also suggest changes in optical cirrus properties (extinction, depolarization) caused by aged air traffic emissions. Still, many research questions remain with respect to the derivation of changes in cirrus properties and radiative forcing for reduced air traffic in 2020.
Combined, these challenges motivated the BLUESKY mission (www.dlr.de/content/en/articles/ news/2020/02/20200522_bluesky-examines-the-atmosphere-during-the-coronavirus-lockdown.html), with the objective to advance our understanding of the anthropogenic impact on atmospheric composition by acquiring a unique dataset on trace gas, aerosol, and cloud properties measured in the early lockdown phase over central Europe and the northern Atlantic flight corridor. Two aircraft were equipped for the BLUESKY campaign: the High-Altitude and Long-Range Research Aircraft (HALO), a Gulfstream 550 with a range of about 8,000 km and 14.5 km cruise altitude (Tadic et al. 2021;Voigt et al. 2017), and the DLR Falcon with 3,000 km range and up to 12 km cruise altitude (Voigt et al. 2011). For BLUESKY, the aircraft were equipped with instruments that measure long-and short-lived trace gases, aerosol, and cloud properties. The aircraft measurements were combined with satellite data of tropospheric NO 2 column densities from the Global Ozone Monitoring Experiment 2 (GOME-2; Munro et al. 2016) and Sentinel-5P/Tropospheric Monitoring Instrument (TROPOMI; Veefkind et al. 2012), and with cloud data retrieved from the SEVIRI Imager on the Meteosat Second Generation (MSG) satellite, as well as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) (Winker et al. 2010). The measurements were also used to evaluate the global chemistry climate model EMAC (Jöckel et al. 2010) and the global-regional chemistry-climate model MECO(n) (Kerkweg and Jöckel 2012a,b;Mertens et al. 2016).
Due to travel restrictions, BLUESKY was conducted out of the home base in Oberpfaffenhofen (48°5′N, 11°17′E) in southern Germany. The observations were performed in May and June 2020 during the early COVID-19 lockdown phase in Europe. This paper describes the aircraft instrumentation and the campaign strategy supported by satellite and modeling activities, as well as first results and highlights from the BLUESKY mission.

Aircraft preparation and instrumentation
During the first days of the lockdown, we developed the idea to deploy the two aircraft HALO and Falcon to measure the atmospheric composition changes during the COVID-19 lockdown. This required adaptation of the Falcon instrumentation to the new science goals by adding instruments to measure O 3 and NO y , which are expected to be influenced by anthropogenic emissions. For HALO, the payload that was planned for the mission CAFE-Brazil, originally scheduled for May 2021, was used for the BLUESKY mission. The preparation of both aircraft was performed under special health, distance, and safety regulations, as there was no COVID-19 vaccine available at that time. Instrument certification mainly was done from home, and remote access to instruments was promoted; thereby, instrument automatization was advanced for future campaigns. The flights were performed with a minimum crew and with the same personnel to minimize contacts; see Fig. 3.
Also, the scientific payloads of the two aircraft were complementary. The HALO payload included a comprehensive set of trace gas and aerosol instruments and the Falcon was equipped with a smaller trace gas and aerosol payload and in addition included cloud probes. Basic meteorological sensor systems to measure meteorological parameters and humidity were operated on both aircraft. The aircraft instrumentations for BLUESKY are described in more detail below.

HALO instrumentation
The heavily instrumented HALO aircraft focused on the detection of short-and long-lived trace gases as well as aerosol properties and composition ( Fig. 4 and Table 1). To this end, several mass spectrometers were operated on HALO: The chemical ionization-atmospheric pressure interface time-of-flight (CI-APITOF) mass spectrometer was deployed to measure nucleating vapors such as sulfuric acid, methanesulfonic acid, highly oxygenated organic molecules, and naturally occurring ions.
The thermal dissociation iodide-chemical-ionization mass spectrometer (I-CIMS) (Slusher et al. 2004;Dörich et al. 2021) monitors peroxyacetylnitric anhydride (PAN), peroxypropylnitric anhydride (PPN), and other species. A proton transfer time of flight mass spectrometer system monitors several volatile organic compounds (Derstroff et al. 2017) including species like isoprene and monoterpenes, acetonitrile, acetone, and methanol. The gas chromatography mass spectrometer (GC-MS) system collects and separates volatile organic components (VOCs) by gas chromatography and detects them with a quadrupole mass spectrometer with a time resolution of 3.5 min (Bourtsoukidis et al. 2017). The instrument provides information on organohalogen compounds such as methyl chloride, selected alkanes, alkenes, aromatics, and alkyl nitrates as well as isopropylnitrate (IPN).
Measurements of NO and NO 2 were made with the Nitrogen Oxides Analyzer for HALO (NOAH) instrument, a modified instrument from ECO-Physics (Hosaynali Beygi et al. 2011). The Airborne Tropospheric Tracer In situ Laser Absorption (ATTILA) spectrometer measured CO and methane (CH 4 ) using dual quantum cascade midinfrared laser absorption spectroscopy. The hydrogen oxide radicals OH and HO 2 were measured with the Hydroxyl Radical measurement Unit (HORUS) instrument (Marno et al. 2020) based on the laser-induced fluorescence technique. The Tracer In situ quantum cascade laser absorption spectrometer/ Hydrogen and Organic Peroxide (TRIHOP) monitor consists of an infrared-laser absorption spectrometer (Schiller et al. 2008) for in situ measurements of formaldehyde (HCHO) and a dual-enzyme monitor (Hottmann et al. 2020) for the detection of hydrogen peroxides H 2 O 2 and ROOH.
In addition, the aerosol composition was measured using a compact time-of-flight aerosol mass spectrometer (C-ToF-AMS; Drewnick et al. 2005;Schmale et al. 2010;Schulz et al. 2018). The C-ToF-AMS analyzes aerosol particles in a diameter range of approximately 50-800 nm and provides quantitative mass concentrations of organic matter, sulfate, nitrate, and ammonium. For aircraft operation, it is equipped with a constant pressure inlet that ensures a steady mass flow and operation pressure of the aerodynamic lens. The Fast Aerosol Size Distribution (FASD) system contains a condensation particle counter (CPC) battery, an Ultra-High Sensitivity Aerosol Spectrometer (UHSAS) instrument, and an optical particle size spectrometer. The CCN-Rack includes a cloud condensation nuclei (CCN) counter (Holanda et al. 2020), a single-particle soot photometer (SP2; Krüger et al. 2022;Holanda et al. 2020), and an impaction aerosol sampler with the aim to analyze the physical properties and the chemical composition of aerosol particles. Near isokinetic aerosol sampling was achieved with the HALO Aerosol Submicron Inlet (HASI) mounted on the fuselage outside of the aircraft boundary layer (Andreae et al. 2018).
Finally, upward and downward spectral actinic flux densities in the range 280-650 nm were measured by combinations of CCD spectroradiometers and optical receivers on the aircraft top and bottom fuselage (Bohn and Lohse 2017), and O 3 was measured with the Fast Airborne Ozone (FAIRO) instrument by a combination of two techniques: a UV photometer and a chemiluminescence detector (Zahn et al. 2012).

Falcon instrumentation
The Falcon instrumentation ( Fig. 5 and Table 2) included aircraft tracers, stratospheric and tropospheric tracers and, in addition, a dedicated aerosol and cloud payload. The sum of gas-phase reactive nitrogen species (NO y ) was sampled through a rear-facing inlet tube and detected as NO by chemiluminescence after reduction of the reactive odd-nitrogen compounds in a heated gold converter (Ziereis et al. 2022). O 3 was measured using an UV photometer (Schulte and Schlager 1996;Ziereis et al. 2000). CO, CH 4 , and CO 2 were measured by cavity ring down spectroscopy (Klausner et al. 2020). The instruments were calibrated using standard mixtures which can be traced back to reference standards of the National Institute of Standards and of Global Atmosphere Watch. The accuracies of the measurements are 15% for NO y and CO, 5% for O 3 , 0.1% for CH 4 , and 0.02% for CO 2 . The atmospheric chemical ionization mass spectrometer (AIMS) uses SF 5 − reagent ions for the detection of upper tropospheric and stratospheric concentrations of gaseous SO 2 , hydrogen chloride (HCl), nitric acid (HNO 3 ), and chlorine nitrate (ClONO 2 ) Jurkat et al. 2016Jurkat et al. , 2017Marsing et al. 2019).
Water vapor distribution was measured with the accurate frost point hygrometer CR2 Kaufmann et al. 2014;Heller et al. 2017) and the liquid or ice water content with the tunable diode laser water vapor analyzer (WARAN) instrument (Kaufmann et al. 2018). Aerosol number concentrations were detected with a set of CPCs with different cutoff diameters of 5, 18, and 50 nm facilitated by different temperature settings and diffusion screen separators (Fiebig et al. 2005;Feldpausch et al. 2006) and a Particle Soot Absorption Photometer (PSAP) measuring aerosol optical properties at three wavelengths of 467, 530, and 660 nm (Virkkula et al. 2005;Virkkula 2010).
Two additional aerosol instruments and two cloud probes were installed in wing stations on the Falcon, including the cloud, aerosol, and precipitation probe (CAPS) (Baumgardner et al. 2004;Voigt et al. 2017) and the cloud and aerosol spectrometer with polarization (CAS-DPOL) (Baumgardner et al. 2004;Taylor et al. 2019;Kleine et al. 2018), as well as two optical particle spectrometers, the UHSAS and Passive Cavity Aerosol Spectrometer Probe (PCASP) (Voigt et al. 2021).

BLUESKY flight strategy and scope
After approximately 6 weeks of aircraft and payload preparation, the Falcon took off for the first instrument test flight on 16 May 2020. The first HALO instrument test flight followed 5 days later. Both aircraft flew five missions together, then the Falcon had to return for maintenance while HALO continued with three more mission flights until 9 June 2020. Altogether 20 flights were performed during BLUESKY, 12 with the Falcon and 8 with HALO, albeit with different flight lengths of up to 2,860 km with the Falcon and up to 6,734 km with HALO, see also Table 3. Extensive measurements were performed over Germany, western Europe, and the Atlantic from 13°E to 14°W and from 38° to 55°N as shown in Fig. 6. During survey flights, many profiles of atmospheric composition were made near large cities.
The strict travel restrictions led to strongly reduced air traffic and the DLR research aircraft were for the first time allowed to perform low-approach maneuvers down to 3 m altitude above ground at large international passenger airports like Frankfurt, Berlin-Tegel, and Milan, as well as Amsterdam, Barcelona, Madrid, Marseille, and Rome.
Many profiles were measured near large cities and the composition of the atmosphere was probed from the boundary layer to the lower stratosphere. We also measured spatially highly resolved data of a multitude of trace species in the upper troposphere and lower stratosphere, in order to investigate changes in atmospheric composition caused mainly by aviation. Observations of the reactive nitrogen and ozone distributions in the northern Atlantic flight corridor and over Europe were performed. According to Eurocontrol, during normal operation in 2019, up to 25,000 flights were performed per day in Europe (Fig. 2). The number of flights was greatly reduced, by more than 80%, in the early lockdown phase and gradually recovered to a reduction of about 30% in summer 2021. This also allowed investigation of aviation induced changes in properties of cirrus cloud as observed from satellite. A third specific focus was the effect of shutdown on strongly populated industrial areas like the Ruhr area in Germany and the Milan region in Italy. In both of them, NO 2 levels were strongly enhanced, as urban emissions are captured in a river valley and therefore the exchange with the free-tropospheric air masses is reduced. Satellite observations of tropospheric NO 2 from GOME-2 and TROPOMI are presented in the next section.
Tropospheric NO 2 changes in early lockdown phase during BLUESKY 2020 Tropospheric NO 2 is a representative short-lived tracer for anthropogenic emissions from transport, industrial processes, and energy production (Müller et al. 2022). While the transport sector and in particular diesel engines are often the prevailing source of nitrogen dioxide in the boundary layer, there are also natural emissions from soil (Lu et al. 2021) and lightning (Pérez-Invernón et al. 2022) that can contribute to the boundary layer NO 2 concentrations. Here we use satellite observations from the GOME-2 (Munro et al. 2016) and the TROPOMI (Veefkind et al. 2012) instruments to motivate the BLUESKY campaign, to give a general overview on the changes of the short-lived tracer NO 2 and to support the planning of individual flight tracks. TROPOMI aboard Sentinel-5 Precursor (2017 to the present) is a nadir-scanning spectrometer with an equator crossing time at 1330 local solar time and daily global coverage. Retrievals of tropospheric NO 2 vertical column densities have a spatial resolution of 3.5 km × 5.5 km (van Geffen et al. 2019). Each orbit (level 2 products of version 1.2) was sampled onto a grid of 0.01° × 0.01° longitude-latitude resolution to enable a consistent averaging and data

CAS-DPOL Cloud and Aerosol Spectrometer
Particle size distribution, polarization, 0.6-50 μm UHSAS Ultra-High-Sensitivity Aerosol Spectrometer Dry particle size distribution 60-1,000 nm Sauer, DLR Voigt et al. (2017) comparison (Müller et al. 2022). Although observations from TROPOMI currently exhibit the highest spatial resolution to study regional effects on the atmospheric composition the data are available from April 2018 onward only and not yet suitable to examine longer-term variability.
To better analyze the deviation of tropospheric NO 2 levels in 2020 from previous years, GOME-2 data were examined from 2015 to 2020. GOME-2 aboard MetOp-B (2012 to the present) is a nadir-scanning spectrometer, which measures at around 0930 local solar time equator crossing time and has a ground pixel size of 80 km × 40 km (Munro et al. 2016). Level 2 offline products of version 4.8 were sampled onto a grid of 0.25° × 0.25° geographical resolution.
Long-term analyses of global satellite data show that NO 2 pollution can vary strongly from year to year due to weather conditions (Zhou et al. 2012;Georgoulias et al. 2019). To reduce this source of variability and uncertainty, observations from GOME-2B for 2020 are compared to the mean from 2015 to 2019 as reference (Fig. 7). For each year the BLUESKY period was considered accordingly. Although the multiannual mean provides a better statistical baseline to quantify possible lockdown effects in 2020, the 3-week sampling period is still short regarding meteorological variability Reductions in tropospheric NO 2 are evident throughout Europe, the United States, and Asia (Fig. 7). The global mean values for the BLUESKY period suggest a 12% reduction in global NO 2 in 2020 with respect to the 2015-19 average. Large cities and conurbations show the highest reductions in NO 2 of up to 55%. In Germany, the measures to limit the spread of the coronavirus started with strong social restrictions and a first lockdown on 22 March 2020. On 15 April 2020, the German government decided to gradually reopen public life. In terms of NO 2 levels, a strong reduction was detected from 22 March to 15 April 2020, and from 15 April the NO 2 levels increased gradually (Erbertseder and Loyola 2020).  Tropospheric NO 2 densities monitored by Sentinel-5P/TROPOMI overpasses in 2019 and 2020 were used to support the flight planning. Figure 8 shows mean tropospheric NO 2 densities for 2020 and 2019 averaged over the respective BLUESKY period from 16 May to 9 June 2020. For the selected area, the mean NO 2 values are 26.8 ± 8.9 μmol m −2 and 31.7 ± 11.2 μmol m −2 for 2020 and 2019, respectively. The reduction in 2020 is evident throughout Germany and corresponds to an overall decrease of tropospheric NO 2 in 2020 by about 15%. The decrease is pronounced in the Ruhr area and in urban areas such as Cologne, Frankfurt, Munich, and Berlin with reductions between 21% and 38%, attributable to strongly reduced road traffic and industrial production, but also to lower NO 2 levels in 2020 before the lockdown and meteorological conditions that are not further disentangled here. A detailed analysis of the lockdown effects on tropospheric NO 2 and their statistical significance considering source attributions of different sectors (Petit et al. 2021;Feng et al. 2021;Putaud et al. 2021), adjustments to multiannual trend effects (Bekbulat et al. 2020), or corrections for meteorological influences (Zhou et al. 2012;Goldberg et al. 2020) is out of scope of this paper. However, the presented findings are in the range of more detailed studies that consider wind corrections and adjustments to trend effects, but focus on different periods of the lockdown (e.g., Liu et al. 2020;Goldberg et al. 2020).

Highlights from the BLUESKY mission
We briefly present first results from the BLUESKY mission and show tropospheric profiles of trace species and aerosol affected by anthropogenic emissions, investigate the blue color of the sky and address effects from reduced air traffic on cirrus.
Reduced reactive nitrogen and carbon monoxide in Frankfurt city profiles and in the industrial boundary layer. Profiles taken over Frankfurt can be compared to climatological datasets measured by instrumented in-service passenger aircraft starting or landing in Frankfurt.   (Clark et al. 2021) are also shown.
The measured NO y and CO mixing ratios in the lower and middle troposphere are significantly lower than the median MOZAIC values by 40%-70% and 20%-40%, respectively, and the CO mixing ratios are also lower than the average IAGOS CO profile. The observed O 3 concentrations agree with the median MOZAIC and IAGOS profiles. BLUESKY NO y and CO concentrations over Frankfurt, however, are close to the 5th percentiles, and only 5% of the MOZAIC data are lower than the BLUESKY profiles. These differences can in part be caused by the meteorological conditions prevailing in May and June 2020 and by COVID-19 emissions reductions.
To better distinguish emission reductions from other factors, we add simulated vertical profiles on the right side of Fig. 9, which show the relative difference between a business as usual simulation and a COVID-19 simulation with a reduced emission scenario [see Mertens et al. (2020) for a full description of the simulations] of NO y , CO, and O 3 mixing ratio around Frankfurt as simulated by MECO(n). Besides the total difference, the differences attributed to the emission sources land transport, anthropogenic nontraffic, and aviation are shown. In fact, the COVID-19 emission reductions assumed in the model simulations lead to a reduction of O 3 of up to 10%, for CO of up to 13%, and for NO y of up to 34%. The NO y reduction near ground are mainly due to reduction of land transport emissions, while reductions above 10 km are mainly attributed to changes of aircraft emissions. This indicates that the observed lower NO y and CO mixing ratios over Frankfurt can to a large extent be attributed to the reduction of emissions from the transport sector (land transport and aviation) and to anthropogenic nontraffic COVID-19 emission reductions.
It is interesting to note, that the mixing ratios of the IAGOS average CO profile for the years 2016-19 are about 7% lower than the MOZAIC average profiles for the years 1994-2005, which corresponds to the observed rate of decrease in CO in the Northern Hemispheric troposphere (Clark et al. 2021). There may be a variety of reasons for the simulated 10% O 3 reduction derived by the MECO(n) model for the O 3 profile on 28 May 2020. Uncertainties in the dynamics of the model, temporal or spatial shifts of patterns in the model compared to the observations, uncertainties in the emissions may contribute to this model result. A deeper analysis is required to better assess these deviations using more extensive model simulations as shown in the "Highlights from the BLUESKY mission" section.
Multiple Falcon and HALO flights were also performed in the industrial Po valley in Italy, a pollution hot spot in Europe, where pollution from large cities, e.g., Milan, are trapped in the river valley and by the Alps. As an example, Fig. 10 shows a Falcon measurement transect through the boundary layer in the western part of the Po valley between 46.6°N, 9°E and 45°N, 10.4°E on 1 June 2020, covering the major industrial area of Milan. Measured CO volume mixing ratio are compared to HYSPLIT simulations interpolated to the Falcon flight path. The HYSPLIT simulations are based on the EDGAR emission inventory for CO of the Po valley region with business-as-usual emissions for the year 2017 adjusted to CO observations during the intense HALO measurements in the Italian Po valley during the Effect of Megacities on the transport and transformation of pollutants on the Regional to Global scales (EMERGE) campaign (Andrés Hernández et al. 2022). The observed CO mixing ratios in the major polluted areas during BLUESKY are lower by about 30% compared to the HYSPLIT simulations with no emission reductions for traffic and industry. The HYSPLIT trajectories enable a detailed allocation of emission inventory sources to the observed pollutant concentrations in the boundary layer. In future studies this will also comprise other pollutants such as nitrogen oxides with origins from several sectors. The example in Fig. 10 shows how concentrations of measured NO y and HNO 3 vary by about one order of magnitude when descending from 2.5 to 1.2 km altitude. Still the HNO 3 relative fraction remains stable, as can be expected from rapid airmass exchange. This instrument performance allows the investigation of the propagation of (reduced) trace gas and aerosol emissions into the free troposphere.
Profiles of other anthropogenic and natural tracers. During BLUESKY, the I-CIMS and GC-MS instruments on board HALO measured two C3-organic nitrates: peroxypropionylnitrate, CH 3 CH 2 C(O)OONO 2 , (PPN), and isopropylnitrate, CH 3 CH(ONO 2 )CH 3 , (IPN). Both nitrates are formed during the OH-initiated oxidation of propane and are thus closely tied to anthropogenic emissions. PPN is formed subsequent to H-atom abstraction by OH from the CH 3 groups in propane and further oxidation. IPN is formed when the H abstraction takes place at the central C atom from the reaction of the isopropylperoxy radical with NO, albeit at low yield. The vertical profiles of PPN and IPN (Fig. 11, campaign average) indicate maximum values in the boundary layer, which reflects the fact that both are formed in the chemical degradation of propane emitted during anthropogenic activity at ground level. PPN has a broad, second maximum at around 8 km while IPN decreases quasi monotonically through the free and upper troposphere. The divergence in their vertical profiles can be understood in terms of their chemically different loss processes and different altitude dependent atmospheric lifetimes. PPN is thermally unstable at low altitudes, but reacts only slowly with OH so that its lifetime (mainly controlled by photolysis) is several months at the low temperatures of the upper troposphere. In contrast, IPN is stable with respect to both thermal decomposition and photolysis and is lost mainly by reaction with OH.
Also, organic compounds like acetonitrile and methanol and sulfur species like dimethylsulfide (DMS) were measured by the GC-MS and the PTR-Mass Spectrometer system on board HALO. Figure 11 in addition shows tropospheric concentrations of acetonitrile, methanol, and DMS for the HALO flight on 4 June 2020 over Germany, France, Italy, and the Mediterranean. Stratospheric data were separated and are treated separately in light of extensive ozone tests that were performed in the laboratory following BLUESKY. During HALO flight 6 (violet color in Fig. 6) the aircraft entered the marine boundary layer over the Mediterranean Sea and therefore DMS is enhanced at low altitudes caused by marine emissions of DMS, while acetonitrile is depleted as it is taken up at the sea surface. Methanol, a volatile organic compound of primarily biogenic origin, shows considerable structure caused by cloud outflow.
On the same flight, gaseous sulfuric acid was measured with the newly developed chemical ionization-atmospheric pressure interface time-of-flight mass spectrometer (CI-API-TOF) with a constant pressure inlet by reaction with nitrate ions. The inlet system was designed in a way to enable a high flow rate and to minimize sampling losses of nonvolatile compounds like sulfuric acid and highly oxygenated organic molecules. Low concentrations of gaseous H 2 SO 4 were measured in the boundary layer and the free troposphere over central Europe.
Here we show count rates respective for H 2 SO 4 as the calibration of the pressure dependence of the mass spectrometer signal is ongoing. The data will be used to investigate the budget and the partitioning of the sulfur species, as also SO 2 and sulfate aerosol were measured on Falcon and HALO, respectively. SO 2 was measured at cruise altitudes in the upper troposphere and lower stratosphere (UTLS) by the AIMS mass spectrometer on board the Falcon in order to investigate the contribution of anthropogenic and natural sources to the UTLS SO 2 budget. In 10-14 km altitude in the midlatitudes, SO 2 has a short life time of 13 ± 2 days (Höpfner et al. 2015). Figure 12 shows the median and quartiles of SO 2 profiles in the UTLS during the entire BLUESKY campaign for flight altitudes between 7.5 and 12.5 km. The maximum height is restricted by the Falcon flight altitude, while the lower limit is caused by a water vapor interference in the SO 2 measurement. The median SO 2 mixing ratio decreases with altitude from around 0.22 to 0.12 μg m −3 . SO 2 has been measured at similar altitudes above Europe in previous experiments:  (Jurkat et al. 2010). In the stratosphere, the SO 2 concentrations are much smaller than sulfate in aerosol, which is still enhanced from volcanic eruptions (see also Fig. 13).
Reductions of aerosol number concentration and aerosol mass in the troposphere. To study the possible influence of reduced ground and aircraft emissions during lockdown on aerosol mass concentrations, we compared the BLUESKY data to previous datasets obtained over Europe, namely, the EMERGE-EU campaign in July 2017 (Andrés Hernández et al. 2022) and the Chemistry of the Atmosphere-Field Experiment in Africa (CAFE-Africa) test and ferry flights conducted over Europe between July and August 2018. Figure 13 shows medians and quartiles of the parameters measured with the C-TOF-AMS on board HALO. For all four measured species, the concentrations measured in 2020 were lower, while in 2018 they were higher. For organics and sulfate, the EMERGE concentrations in 2017 were significantly higher than in 2020, while nitrate and ammonium show comparable values. Although the variability in the datasets is large, the consistent reductions indicated from the 2020 data for all four aerosol components suggests that they were at least partly due to diminished emissions during the lockdown. By combining our black carbon measurements with EMAC simulations we found a 40% reduction in black carbon mass related to the lockdown effects (Krüger et al. 2022). For organic aerosol mass, this finding holds also for the free troposphere above 5 km altitude, whereas for sulfate, nitrate, and ammonium the free-tropospheric data are comparable for all three datasets.
A recent model study using the EMAC model with reduced emission scenarios is able to reproduce the observations during the BLUESKY campaign (Reifenberg et al. 2022). Anthropogenic aerosol precursor gases include SO 2 , NO x , and VOCs. Although the emissions reductions had slightly leveled off in May and June compared to April 2020 (Guevara et al. 2021), still the NO x reduction was highest with about 10%-20%, followed by SO 2 (9%-12%) and VOCs (2%-6%). Ammonia originates mainly from agriculture; therefore, no significant reduction was observed. However, the amount of NH 4 in the aerosol is determined by available acids (HNO 3 and H 2 SO 4 ) to react with NH 3 to ammonium nitrate and sulfate. Direct PM 2.5 emissions were reduced only by 6%-9%. The finding that the observed differences in aerosol mass concentrations are much larger that the reduction in direct emissions suggests that secondary aerosol formation is the most important aerosol source over Europe. Below 2 km, sulfate and organics show a higher reduction compared to the previous campaigns than nitrate, indicating that the anthropogenic emissions of SO 2 and VOCs play a larger role for secondary aerosol formation than anthropogenic NO x plays for secondary nitrate formation.
As discussed above, the strong increase of sulfate aerosol above 10 km and the higher sulfate mass concentrations observed in 2020 compared to 2018 may in part be due to the Raikoke volcanic eruption in June 2019 (Muser et al. 2020;de Leeuw et al. 2021), but also other, minor volcanic eruptions are likely to play a role. Figure 14 shows profiles of fine mode aerosol number concentrations for particle sizes above 18 nm from the Falcon measurements during BLUESKY and the long-term observations from the In-service Aircraft for a Global Observing System-Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container (IAGOS-CARIBIC) project. The IAGOS-CARIBIC data shown here only include flights during the months May and June between 2005 and 2015. Profile data are restricted to central Europe between 8.8°W and 28°E and at latitudes between 40° and 66°N.  40°-66°N, 28°E-8.8°W).
The median BLUESKY profile shows a clear reduction of the aerosol number concentration of 30%-70% in the free troposphere above 4 km altitude. At the tropopause above 11 km, the difference becomes less prominent, although the BLUESKY median concentrations are still mostly below the CARIBIC values, even though the IAGOS-CARIBIC data were measured 5-15 years before BLUESKY with respective lower anthropogenic emissions. It is not possible to quantify an effect in the boundary layer although it has to be kept in mind that the BLUESKY flights represent a larger variety of low-level flights in different areas of Europe and the CARIBIC dataset is composed of primarily the ascents and descents into and out of Frankfurt international airport.
Blue sky during BLUESKY. The solar radiation at Earth's surface reached a maximum during the COVID-19 lockdown period in spring 2020 over western Europe. Van Heerwaarden et al. (2021) found that surface irradiance was highest since 1928 for the area of the Netherlands. Reasons were an overall very low cloud fraction, several exceptionally dry days, and low aerosol optical thickness (AOT). Regardless of whether AOT concentrations are ultimately due to an anthropogenic influence or related to meteorological variability, the subjective impression was that in 2020 the sky appeared in a deeper blue during the lockdown period. To substantiate this perception, we carried out radiative transfer calculations in the cloud-free atmosphere.
Simulations of the atmospheric radiative transfer are based on the 1D radiative transfer model UVSPEC from the program package libRadtran (Mayer and Kylling 2005;Emde et al. 2016). In this study, UVSPEC is used to compute downward-directed spectral radiances and irradiances at the surface in the visible spectral range (380-780 nm) for a cloud-free midlatitude summer standard atmosphere (Anderson et al. 1986). Regarding the aerosol, a model after Fenn et al. (1985) was used, which assumes a rural type for the boundary layer. To examine changes in sky color, calculated spectra are convolved with the spectral sensitivity of the human eye in the red, green, and blue spectral range. A subsequent spectral integration gives RGB color components and enables the representation of the sky color. Input to the model are different vertical distributions of the aerosol optical thickness AOT(z). Profiles close to reality have been generated by scaling a normalized profile of the volume extinction coefficient (Fenn et al. 1985) with a total AOT provided by the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the Terra and Aqua satellites.
The total aerosol optical thicknesses derived from MODIS data products MOD08_D3 v6.1 (Terra) and MYD08_D3 v6.1 (Aqua) were statistically analyzed over Germany and Europe. Here we consider the BLUESKY period from 23 May to 9 June 2020 as a reference to the years 2015-19. For 2020 the median total aerosol optical thickness AOT MODIS is 0.156 at 0.55 μm, whereas it is 0.247 for the period 2015-19. For comparison, the individual years 2015-19 give an AOT MODIS range between 0.187 and 0.309. Also, the statistics for Europe shows a median of 0.156 for the BLUESKY period in 2020 and a range of 0.175 to 0.265 for 2015-19. Vertical profiles of the aerosol optical thickness based on the AOT MODIS values of 0.156 and 0.247 are presented in Fig. 15a. Figure 15b shows the spectra of the diffuse component of downward directed spectral irradiances for a solar zenith angle of 45°. Calculating the ratio r of downward irradiances at blue and red wavelengths (I 460nm /I 750nm ), for example, results in r = 2.4 and r = 2.1 for AOT MODIS of 0.156 and 0.247, respectively, indicating a shift of the simulated spectrum toward shorter wavelengths. Figures 15c and 15d reveal that the spectra are associated with a contrast in sky color. Transferring angular resolved radiances and hemispherically integrated irradiances into RGB colors clearly illustrates that in case of a reduced AOT MODIS (Figs. 15c and 15d, left halves of circles), the sky appears in a deeper blue color. This result also applies to other solar zenith angles. Figure 15e displays the year-to-year variability of the AOT MODIS in the period from 2015 to 2020 highlighted with calculated sky colors. It is obvious that the blue in 2020 forms a contrast to all previous years, albeit differently pronounced depending on the AOT MODIS differences. Although it is not clear whether changes in the total aerosol optical thickness are mainly due to an anthropogenic or a meteorological influence, the radiative transfer simulations give a clear indication that the lower aerosol optical thickness during the lockdown resulted in a deeper blue color of the sky during BLUESKY compared to the same period in the previous five years. Global model simulations in the "Highlights from the BLUESKY mission" section further investigate the impact of aerosol changes on the radiation budget at the surface.
Effect of pollution on low clouds. To contrast aerosol-cloud interactions in highly polluted regimes against less polluted regimes, intensive measurements of aerosol and clouds were made in the boundary layer over the densely populated Ruhr metropolitan region versus recurrent sampling of the boundary layer over the remote sparsely populated Hohenpeissenberg regions in southern Germany. Further aerosol and cloud measurements were taken above the cities of Berlin, Hamburg, and Munich, during individual arrival and departure routes to and from the respective airports. Cloud droplet number concentrations (CDNC) in low-level clouds (<2,500 m altitude) were measured with the Cloud and Aerosol Spectrometer (CAS) probe aboard the DLR Falcon and more than 8,000 individual clouds were intercepted. We implemented a cloud threshold considering only droplets larger than 3 μm and with a liquid water content of more than 0.01g m −3 , in order to avoid diffuse cloud edges. Only measurements contributed to the remote Hohenpeissenberg statistics from campaign days without a large-scale westerly flow, to avoid an influence of the Munich urban emission plume. Figure 16 shows the cloud droplet number concentrations measured in the remote and the industrial regions. A clear trend is observed between the CDNC in the rural Hohenpeissenberg area with median cloud drop number concentrations of a few tens per cubic centimeters at the lower side of the spectrum and high median CDNC of 385 cm −3 in the densely populated Ruhr metropolitan area. The cities of Hamburg (195 cm −3 ), Berlin (274 cm −3 ), and Munich (311 cm −3 ) range in between.
Reduced contrail cirrus cover and radiative impact. The COCIP (Schumann 2012) simulates the life cycle of contrails (Schumann et al. 2017;Schumann and Heymsfield 2017) that form for given aircraft flight track data and numerical weather predictions and computes their local radiative forcing (Schumann 2012). COCIP has been used for mission planning, for simulation of contrail properties comparable with in situ and satellite observations, and for aviation climate impact mitigation studies (Schumann and Graf 2013;Voigt et al. 2017;Teoh et al. 2020;Schumann et al. 2021a). Figure 17 presents example results for 16 April 2020. Travel restrictions in response to COVID-19 caused an 89% decrease in air traffic flight distance over Europe this day. Despite the low air traffic density, MSG/SEVIRI data indicate contrail occurrence on that day, as shown in the optical thickness of ice clouds, derived using the algorithm described in Strandgren et al. (2017), and from brightness temperature differences.
Contrail and cirrus optical thickness (COT) derived from COCIP on that day are shown in the middle panel in Fig. 17, as well as COCIP COT calculated with the same meteorology but with air traffic data from 16 April 2019. Higher contrail cirrus optical thickness is calculated for the factor-of-5-higher air traffic on 16 April 2019. COCIP was also used to calculate the respective positive and negative radiative forcing (RF) from contrail cirrus for both scenarios. The resulting daily net forcing is positive on average in both years. The computed longwave and shortwave RF components locally exceed 1 W m −2 in magnitude despite low air traffic in 2020. On 16 April 2020, the contrail RF in the presented area is also strongly reduced to about 20% of its value for the same meteorology but for air traffic data in 2019, implying a significant reduction in RF from contrails due to reduced air traffic in 2020 on that day. Annually averaged changes might be different due to the warming and cooling components of the contrail cirrus forcing (Gettelman et al. 2021). A more detailed discussion of the contrail effects including an uncertainty analysis is given in Schumann et al. (2021a) and the extension to longer time periods is discussed in Schumann et al. (2021b).
Changes in cirrus properties from CALIPSO data. Civil aviation may contribute to climate change by inducing contrail cirrus formation and by changing natural cirrus cloud properties (Tesche et al. 2016). Urbanek et al. (2018) showed that cirrus clouds formed in regions highly affected by air traffic have higher mean particle linear depolarization ratios than those formed in pristine regions. To study the effect of reduced aviation on cirrus cloud properties we compared CALIPSO satellite measurements of cirrus during the BLUESKY campaign with cirrus data measured at the same time period in the previous years . The measurements were performed with CALIOP instrument on board the CALIPSO satellite (Winker et al. 2007(Winker et al. , 2010Stephens et al. 2018). CALIOP uses three channels at 1,064 and 532 nm for measuring the total backscatter and the orthogonal component of depolarization. For this study we used the level 2 5 km cloud profile products with a vertical resolution of up to 30 m. To compare the CALIOP cirrus observations with the measurements during the BLUESKY campaign, we analyzed the datasets covering the Falcon cruise area from 50° to 55°N and from 15°W to 5°E as well as latitudes from 45° to 55°N and longitudes from 5° to 15°E. CALIPSO passes over the study area three to four times each day resulting in 90 overpasses in the period 16 May-9 June each year. Since all BLUESKY measurements were performed during the day, we only used daytime satellite lidar observations for the intercomparison.
To investigate the changes of cirrus cloud properties as possible consequence of the significantly reduced aviation during the COVID-19 pandemic, we used the particle linear depolarization ratio (PLDR), the extinction coefficients (σ ci ), and the calculated effective optical thickness of cirrus clouds ( where r t and r b are the cloud-top and cloud-bottom heights of cirrus, respectively). As cirrus cloud properties strongly depend on temperatures (e.g., Urbanek et al. 2018), and as a temperature T of −50°C is one of the threshold conditions for contrail formation (Schumann 1996), we divided the data into two subsets with temperatures from −38° to −50°C and temperatures below −50°C (Li and Groß 2021). Figure 18 shows the probability density function (PDF) of the PLDR (left) for cirrus clouds at T < −50°C (top panel) and −50° < T < −38°C (bottom panel), as well as of the extinction coefficient (top right) and the effective optical depth (bottom right) for T < −50°C of all cirrus cloud cases measured during the day with CALIOP during the BLUESKY period in the different years. The resulting distributions of cirrus PLDR show a non-Gaussian shape for both temperature ranges with larger values at the lower temperature range (T < −50°C). The PDFs of the PLDR of cirrus clouds in the higher temperature range show almost no deviation from each other, and also the cirrus cloud measurements in 2020 are well within the variability of the previous years. The median of the PLDR distribution is at about 0.35-0.38 for all the years. For the cirrus PLDR at the lower temperature range, the situation looks quite different. The years 2010-19 show only limited variability with a median value between 0.42 and 0.45. However, the PLDR for the cirrus clouds measured during BLUESKY in 2020 show a significantly lower median value of 0.39. The σ ci and τ ci distributions show the same features as PLDR; σ ci and τ ci values for cirrus clouds in 2020 are on average smaller at lower temperatures than in the previous years. The median value of the extinction coefficient for the years 2010-19 varied between 0.20 and 0.26 km −1 while it was about 0.14 in 2020. A similar but not as significant picture is found for the cirrus cloud effective optical thickness. Most median values in the years 2010-19 were within a range of 0.10 and 0.16, whereas the median value in 2020 was about 0.09. Low values are, however, also found for cirrus clouds in 2013 and 2015.
Global chemistry-climate modeling. The input and the comparison of measurement data with global models are essential to enhance our understanding of atmospheric processes and their relation to climate. The EMAC atmospheric chemistry-climate model was used both for forecasting and the postcampaign data analysis. In both cases EMAC was nudged by Newtonian relaxation toward data from the European Centre for Medium-Range Weather Forecasts (ECMWF). The EMAC atmospheric model core is the fifth-generation European Centre Hamburg general circulation model (ECHAM5). For higher resolution in the field campaign region, the forecasts of atmospheric composition were performed with the MECO(n) model, which couples the global EMAC model online with the enhanced resolution regional model COSMO/MESSy, allowing for seamless zooming into regions of interest in the global model context (Kerkweg and Jöckel 2012a,b;Mertens et al. 2016). A continuous analysis simulation is performed from which 5-day forecasts are branched of every 12 h, with nudging toward the ECMWF operational analysis/forecast data. MECO(n) was configured with one refinement ranging from eastern North America to eastern Europe with 0.44° resolution in latitude and longitude. Comprehensive tropospheric gas phase chemistry was calculated as described by Mertens et al. (2016), in addition the source apportionment method by Grewe et al. (2017) and Mertens et al. (2020) was applied. Here, it uses a business as usual air traffic scenario to help identify the regions where large aviation signals onto the NO y mixing ratios are predicted.
As an example, Fig. 19 shows the NO y mixing ratios at 250 hPa averaged for 0700-1600 UTC 2 June 2020 and the relative contribution of aviation nitrogen oxide emissions to the NO y mixing ratios. The model forecast shows regions with relative contributions from air traffic of more than 50% predicted west of Ireland for a business as usual air traffic scenario. The predictions were used to guide both aircraft into regions with aircraft emissions for aviation in 2019, being reduced during the same period in 2020. Accordingly, HALO and Falcon performed measurement flights in the regions west of Ireland (Fig. 19). The NO x and NO y measurements on the HALO and the Falcon during their flights into this area show individual spikes in NO x and NO y indicating the presence of a few aircraft plumes. In addition, NO x and NO y plumes from aircraft flying on 2 June 2020 as well as for air traffic in 2019 were calculated with the COCIP model and folded onto the Falcon and HALO flight paths. Due to reduced air traffic in 2020 only few individual plumes are derived, while calculations for emissions from the denser air traffic in 2019 show many more spikes.
Postcampaign data analysis was also performed with the EMAC model, adopting reduced emissions resulting from the lockdown in Europe as estimated by Guevara et al. (2021). EMAC was nudged toward the ERA-Interim data (Berrisford et al. 2011) to reproduce the observed synoptic weather conditions. Figure 20 shows a comparison of observed temperature, specific humidity, CO, NO, O 3 , and photolysis frequencies j(NO 2 ) from HALO measurements and EMAC results from simulations with the reduced emission (COVID-19) scenario (Reifenberg et al. 2022). The good agreement between observed and simulated data for temperature, which is a nudged variable and therefore shows the same variance as the ERA-Interim dataset, indicates the quality of the reproduced meteorology. As expected, somewhat larger deviations from the observations are found for unconstrained variables, such as O 3 , CO, and NO, although the overall comparison shows a very good agreement of the low-emission scenario simulations with the HALO measurements. In particular, the model is able to reproduce the observed NO, which was strongly reduced in the entire tropospheric column, because of the strong reductions in ground and air traffic (e.g., Schumann et al. 2021a). Furthermore, also parameters which are strongly influenced by parameterizations are satisfactorily represented in the model, such as the photolysis frequency of j(NO 2 ), whose variability is associated with clouds. In addition, O 3 is in general agreement between model and observations, albeit with modeled O 3 concentrations tending to slightly higher values. The EMAC model has also been used to investigate the impact of reduced emissions on direct and indirect aerosol radiative forcing during COVID-19 lockdown in Europe. Reifenberg et al. (2022) find large differences of selected tracers and aerosol between the reduced emission (COVID-19) scenario and the business as usual scenario at aircraft cruise altitudes in the upper troposphere mainly due to reduced air traffic. In addition, noticeable differences are found in the boundary layer affected mainly by ground transportation and industry. The reduction in aerosol surfaces leads to an increase in incoming solar radiation at the surface during the BLUESKY period (Reifenberg et al. 2022), in addition to the blue sky color.

Conclusions and outlook
From 16 May to 9 June 2020 the HALO and DLR Falcon performed 20 flights over Europe during the early COVID-19 lockdown to investigate the impact of reduced anthropogenic emission on the atmospheric composition during the BLUESKY mission. A comprehensive and unique dataset of trace gases, aerosols, and clouds was measured from the boundary layer to the lower stratosphere and profiles of atmospheric constituents were derived and compared to measurements in pre-COVID-19 times and to satellite data and models. This paper presents first highlights from BLUESKY: 1) Significant (10%-50%) tropospheric NO 2 reductions were observed over industrialized continental areas as well as in major city outflows during the early lockdown phase, indicated by TROPOMI and GOME satellite data. 2) Tropospheric NO y and CO profiles over Frankfurt showed significant (20%-70%) reductions on 28 May 2020 compared to the MOZAIC climatology from 2004 to 2015 and to IAGOS measurements from 2016 to 2019. The measured PPN and IPN profiles are in line with reduced anthropogenic influence on atmospheric chemistry. Falcon measurements of CO and NO y in the industrial boundary layer of the Po valley in Italy show reductions up to 30%. EMAC model results confirm the impact of reduced emissions on the NO y and CO profiles over Frankfurt. A comparison of EMAC results to BLUESKY O 3 data suggests that tropospheric O 3 is in general agreement with the observations, but tends to be slightly elevated compared to the HALO observations during the lockdown phase in Europe.
3) A suite of sulfur species was measured, such as DMS, SO 2 , and H 2 SO 4 as well as sulfate aerosol, which makes it possible to investigate the sulfur budget in the UTLS. The stratospheric sulfate was still perturbed by the aged emissions from the volcanic Raikoke eruption in June 2019 and by smaller eruptions thereafter. 4) The aerosol fine mode number concentrations and mass were substantially reduced in continental profiles with respect to observations from previous European summer campaigns and from the MOZAIC dataset. In the lower troposphere blow 5 km, we found strong aerosol mass reductions. For black carbon aerosol, the lockdown-related reduction was as high as 40% (Krüger et al. 2022). Reduced organic aerosol particulates aloft imply that reduced VOC emissions transported from the surface might have led to a reduction of secondary organic aerosol production in the free troposphere. 5) The perceived deeper blue sky during the BLUESKY period in central Europe can be explained by reduced light scattering from aerosol profiles with a 40% lower AOT than usual. 6) Low-level clouds measured in remote parts of Bavaria contained lower cloud droplet number concentrations compared to clouds over large cities and the industrial Ruhr region. 7) The 80% reduced air traffic led to significant reductions in contrail cover, contrail optical thickness and to a reduced, but still positive radiative forcing from contrail cirrus for a case study relating to 16 April 2020. Also, observed reductions in NO y at cruise altitudes can be explained by simulated reductions in aircraft NO x emissions. 8) The extinction and the depolarization ratio of cirrus clouds below −50°C measured by CALIOP on daytime between 16 May and 9 June 2020 exhibits reductions at aircraft cruise altitudes above Europe implying a potential effect from reduced aircraft aerosol emissions on cirrus properties. 9) The global model EMAC has been used in its nudged version for flight planning. First postcampaign simulation results indicate good agreement for a scenario with emissions that were reduced during the lockdown period. Reifenberg et al. (2022) derive a higher fraction of incoming solar radiation at the surface for the reduced COVID-19 emission scenario compared to the business as usual simulations, mainly due to reduced aerosol surface areas.
Altogether, BLUESKY provides a comprehensive dataset on trace gases, aerosols, and clouds over Europe, which documents the anthropogenic impact on atmospheric composition and climate.