1. Introduction
One of the deadliest and largest global pandemics in history, the Coronavirus disease 2019 (COVID-19), firstly appeared in November 2019 in China, Wuhan and causes respiratory illness [
1,
2]. The first case of COVID-19 in Europe was confirmed on February 21, 2020, in France (Spiteri et al., 2020). Until the end of April 2023 around 764.500.000 cases have been confirmed globally, with 6.915.000 deaths [
3]. Europe declared 276.000.000 confirmed cases. To combat the virus and to reduce the infections and mortality, governments put in place numerous measures such as travel restrictions, school- and workplace closures, even complete lockdowns [
4,
5]. In consequence, unprecedented in history, the virus reduced various human and economic activities for several months [
6]. This exceptional situation changed the environment in many ways. For example, changes of the anthropogenic heat release related to road traffic emissions and energy consumption for heating and cooling buildings, modified the air- and land surface temperatures (LST) of cities [
7,
8,
9]. Numerous studies investigated the LST changes during the first lockdown [
9,
10,
11,
12]. They observed a general decline in LST compared to the previous year’s average LST. For Andalusia a LST decline of -4.6 K (-19.3%) during the March to June 2020 compared to the same period in 2019 based on Sentinel 3 data was recorded [
13]. Further, Liu et al. [
9] found that during the lockdown the surface urban heat island intensity (SUHII) in China decreased by 0.25 K during the day and 0.23 K at night, and the canopy-layer UHII by 0.42 K at day and 0.39 K at night, respectively. Also, regarding the air pollutants many studies have already shown that the lockdown restrictions affected anthropogenic-related air pollution [
4,
6,
13,
14,
15,
16]. This is especially important because indoor and outdoor air pollution is one of the greatest health risks for people nowadays, claiming about seven million lives annually [
17]. During the lockdown spatial differences in the intensity of changes were recorded. They are mainly explained due to different strict measures imposed by each government, the prevailing sources of the emissions and the weather [
18]. Strongest air pollution drops were seen in Asia and then in Europe. Less strong drops were registered in North America and smallest changes in Africa due to less strict measures [
19].
In detail, looking at nitrogen dioxide (NO₂), the European Space Agency (ESA) [
20] noted a 40-50% reduction over Asia and Europe, derived from the Sentinel-5P satellite between the end of January and the beginning of February 2020 compared to the same period in 2019. The study of Tobías et al. [
21] based on ground measurements showed similar signals. Here, Barcelona had reductions of -45% to -51%. Both, ground measurements and satellite-based studies concluded that the main contributors for the NO₂ reduction are the decline of road transport and industrial emissions. However, after easing the restrictions, concentrations were approximately as high as before the lockdown [
18].
Observations of fine inhalable particles with diameters 2.5 µm and smaller (PM
2.5) were inconsistent. For example, in Chinese cities drops were generally greater in more industrialized cities [
22].In comparison to 2017-2019 reductions of -42% were noted for Wuhan [
16]. Rural areas, where agriculture is the main activity, or places where PM is more prevalent due to natural sources, PM remained at higher levels. Also, the strictness of the lockdown affected the PM
2.5 reduction. Compared to 2019, reductions of -7.1 µg/m
3 without strict measures and -21.1 µg/m
3 with strict measures were reported [
22]. In South European cities there were only slight PM reductions (-8%) compared to 2017-2019 [
16]. The drops were recorded especially at traffic stations and hence attributed to transport and fuel combustion reductions. However, increased domestic heating and garden activities like biomass burning compensated those declines.
A widespread increase was seen regarding ozone (O
3). For example, in Barcelona and in Andalusia higher O
3 concentrations were noted, with +33% to +57%, and +5.9% respectively, obtained from meteorological ground stations, compare to pre-covid levels [
13,
21]. Another study recorded an O
3 increase of +17% compared to 2017–2019 for Europe [
16]. The increase is explained due to reductions in NO
x emissions resulting in lower O
3 titration leading to higher concentrations of O
3. Further, it must be considered that O
3 formation is weather dependent i.e., photochemical sensitive. The sunny weather in this period lead to a higher O
3 formation.
Thus, the emergence of COVID-19 offers a unique opportunity to comprehend and quantify the human impact on the environment. However, most studies focus only on individual cities and single variables which may not be sufficiently representative, e.g. [
10,
13,
22,
23,
24]. Considerably fewer studies analyzed patterns within a continent or at global level, e.g. [
16,
18,
25]. In fact, there is a lot of annual variation and thus differences between cities. Hence, cities may show significant changes in different directions. Thus, the main objectives of this article are first to perform a multiparameter analysis and to comprehensively document the spatial and temporal LST and air pollutant variations, namely NO
2, O
3, and PM
2.5, during the first lockdown period for 43 cities over Europe, compared to the reference period in 2015–2019; and secondly to determine the influence of altered anthropogenic activity on those variables.
3. Results
3.1. Multiparameter overview
Figure 2 shows the results of the individual variables of the KS-test combined with the relative changes in 2020 compared to the reference period. For the KS-test in most cases the null hypothesis can be rejected i.e., the LST, NO
2 and O
3 differences between the lockdown and reference periods are statistically significant. The statistically significant results are shown as a circle. If there is no statistically significant change in the variable, this is indicated by a rhombus as a symbol. This is the case in Eastern Europe for the O
3 values and the LST night-time values. For PM the results are variable with no clear signal regarding the spatial distribution.
In general, a predominant decline for ground-level NO2 was recorded. Although most of the cities showed the same signals, the magnitude of change differed city-wise. The three cities with the largest percentual NO2 reductions are Luxembourg (-54.0% relative change, -12.3 µg/m3 absolute change), Riga (-50.7%, -5.1 µg/m3), and Belgrade (-50.4%, -6.8 µg/ m3). Hamburg (+10.5%, +1.1 µg/ m3), Tirana (+18.4%, +0.7 µg/m3) and Naples (+11.9%, +1.7 µg/m3) show the largest positive anomalies. For O3 a widespread increase is evident. The largest positive anomalies are found in Luxembourg (+41.1%, +20.5 µg/m3), Cologne (+35%, +15.7 µg/m3) and Paris (+27.1%, +13.1 µg/m3). In contrast to the other cities, cities in the Iberian Peninsula, Italy, and Southern France show lower O3 concentrations in 2020. The greatest changes correspond to Naples (-12.1%, -8.9 µg/m3), Valencia (-10.3%, -7.5 µg/m3) and Madrid (7.2%, -5.0 µg/m3). PM2.5 anomalies are inconsistent over space. They decreased predominantly in Northern Europe. Strongest reductions can be observed in Tallin (-38.4%, -2.9 μg/m3) Vilnius (-37.4%, -4.6 μg/m3), Oslo (-28.8%, -2.8 μg/m3). The highest increases were observed for PM2.5 in Dublin (+43.9%, +1.5 μg/m3), Turin (+31.4%, +5.5 μg/m3), Milan (+19.7%, +4.1 μg/m3).
3.2. Nitrogen Dioxide
Figure 3 shows the change in the distribution of the daily NO
2 levels averaged for all cities. Here it is apparent that the mean of the NO
2 concentration during the lockdown 2020 (dark red line) and the binned observations for the individual concentrations are left-shifted, which means that they have decreased. Furthermore, the density is reduced as well, which underlines the flatter continuous density curve for 2020. The density is calculated by dividing the frequency by the class width. Thus, it represents the frequency per unit for the data in each class. At this point it must be emphasized that the individual cities differ a lot. For some cities, like Barcelona the change is high, and in others, like Valetta the change is low (
Figure 4).
The anomalies over the curfew period are shown in more detail in
Figure 5. Here the differences for each day and city are shown compared to the values for each city in the reference period. In general, the changes in the Eastern European and Scandinavian cities are rather small. It is noticeable that there is almost no NO
2 concentration difference in Sarajevo. In contrast, a much lower NO
2 concentrations in 2020 (i.e., a strong negative anomaly) can be seen in Athens and Luxembourg. Some cities show an inconsistent behavior like Brussels, London, Milan, or Paris. On the other hand, Hamburg, Tirana, and Naples predominantly show a positive change in NO
2 concentrations compared to the reference period. Considering all cities, strong concentration declines in 2020 from day 81 of the year (March 22
nd) approximately until day 91 (April 1
st) are noticeable. Thereon, several cities show concentration increases compared to the reference period, interrupted again from negative anomalies during day 104-106 of the year (April 14
th – 16
th).
Looking at the average weekly pattern of all cities (
Figure 6), the absolute NO
2 concentrations are not only lower but the pattern itself changed. In general, measured NO
2 concentrations increase, peaking on Fridays and decrease towards the weekend. A comparatively stronger increase in concentrations from Monday to Friday can be seen in 2020 (~3.5 µg/m
3 in 2020 vs. ~2 µg/m
3 in the reference period). On Sundays and Mondays, the concentrations are the same in 2020. In contrast, in the reference period, higher concentrations on Mondays compared to Sundays are observed.
3.3. Ozone
Examining the O
3 anomalies over the whole period, lower concentrations seem to prevail at the beginning of the lockdown. Around day 80 to 90 of the year (March 22
nd -March 31
st), however, almost all cities show an increase in O
3 concentration (
Figure 7). Very striking is the strong increase in Budapest, Berlin, Cologne, Luxembourg, and Paris. This is followed by a phase which tends to have less O
3 in 2020 that lasts until day 95 (April 5
th). Until day 108 (April 18
th), most cities show higher O
3 concentrations than in the reference period. The cities Lyon, Marseille, Naples, Rome, Madrid, and Valencia, show conspicuous strong negative anomalies, especially towards the end of the study period from day 109 (April 19
th).
In addition to the fact that surface O
3 concentrations are mostly higher in 2020, concentrations tend to be higher on weekends (
Figure 8). The pattern itself in 2020 compared to the reference period is similar. However, there is a smaller concentration increase in 2020 (~1.9 µg/m
3) towards the weekend than for 2015-2019 (~4.2 µg/m
3). In addition, O
3 concentrations decrease in the reference period from Monday to Friday whereas for 2020 they do only until Wednesday. Further, the pattern is inverse in comparison to NO
2.
3.4. Particulate Matter
PM
2.5 shows only a small shift of the mean concentration (
Figure 9). Looking at the histograms, the bins are distributed in a greater range in 2020. Thus, both, observations with higher and also with lower daily values are recorded in 2020, even though it has a comparatively lower density curve. The anomalies are not only inconsistent over space but also over time (
Figure 10). However, it is noticeable that between day 77 and 80 (March 18
th - March 21
st) in some cities there was a simultaneous increase of PM levels in 2020, followed by a drop until day 86 (March 27
th) approximately. From day 87, with a few exceptions (e.g., Valetta, Naples), there was a positive anomaly that lasted about three days.
3.5. Land Surface Temperature
For both, LST and SUHII there is no clear signal of change (
Figure 11). In total there was only a very slight LST and SUHII reduction. Weather uncorrected data show both, higher and lower LST in 2020 compared to the reference period and diverge strongly temporally and spatially (
Figure 12).
3.6. Mobility Data
Figure 13 shows the percentual change of anthropogenic mobility during the study period (red box: March 15
th to April 30
th, 2020) and the time after the first lockdown in Madrid and Stockholm. The cities of Madrid and Stockholm were selected as examples of a strict lockdown and a lockdown with very few ordinances, respectively. Spain and Italy in general had the most stringent measures, while the government of Sweden imposed almost no restrictions [
34]. Scandinavian and Baltic cities had less strict measures, which was also the case in Germany, where schools have been closed, but the industrial sector remained largely open [
35].
After the restrictions have been imposed in mid-March, the mobility was clearly reduced. All categories were visited less, except of “residential” because more time was spent at home. In general, the greatest changes can be observed in the categories retail & recreation and transit stations; and the least in parks & grocery, and pharmacy. After easing the restrictions, mobility gradually increased again, with parks being visited more than usual. It must be stated that the baseline is in February and that fewer parks are visited in winter than in summer anyway. However, the mobility behavior has still not returned to the levels seen in summer 2020. Furthermore, the data show a weekly pattern. People were more outside and less at home on weekends, but in total less than usual. Regarding the workplaces and transit stations more people worked from home and did not use the public transport. In contrast on weekends most of the employees do not work and changes are rather marginal.
Figure 13.
Percentual mobility change from the baseline, Madrid and Stockholm during 15.03-01.08.2020.
Figure 13.
Percentual mobility change from the baseline, Madrid and Stockholm during 15.03-01.08.2020.
Comparing the two cities, in Madrid, mobility was much more restricted, public life almost came to a halt and thus the mobility data show a greater decline. Transit stations, park visits, and retail & recreation had a change of almost -90%. Furthermore, significantly more time was spent at home. Buying groceries was the only allowed opportunity to leave the house [
34]. In contrast, in Stockholm the surplus of time spent at home was significantly lower (~+20%) than for Madrid (~+35%). Also, changes at transit stations and workplaces were only about 50%. From the end of March and ongoing, parks were visited more than during the baseline period, reaching a maximum of +75%. At the beginning of June, the visit to parks increased considerably compared to the baseline in Stockholm (the large daily variations are because park visits are influenced by the weather conditions). Thus, the two exemplary selected cities clearly illustrate how different strictness levels are reflected in changes in mobility behavior.
In
Figure 14 we examine how the daily mobility is related to the variables under study, where each point corresponds to the value of one city on one day. As an example, the change of the number of visitors at transit stations is selected. The other categories follow the same pattern except of "residential" which shows an inversed pattern (not shown). The data distribution of the change of visitors at transit stations and the absolute NO
2 values is clearly not linearly correlated (
Figure 14 a). Furthermore, it is useful to look not only at the absolute NO
2 values, but also at the change in NO
2 concentrations compared to the change at the transit stations. The representation of two rates of change allows inferences about how the NO
2 values change related to the change of the number of visitors at the transit stations. In
Figure 14 (b) it can be seen that there were less people at the transit stations and that the NO
2 concentration has predominantly decreased during COVID-19. More importantly it becomes evident that the change of visitors at the transit stations has no visible influence on the NO
2 change. Overall, we could not establish a clear relationship between the mobility and the variables under study.
Author Contributions
Conceptualization, P.G., B.B. and P.S.; methodology, P.G., B.B. and P.S.; software, P.G.; validation, P.G.; formal analysis, P.G.; investigation, P.G.; resources, P.G. and P.S.; data curation, P.G. and P.S.; writing—original draft preparation, P.G.; writing—review and editing, P.S. and B.B.; visualization, P.G.; project administration, P.G.; All authors have read and agreed to the published version of the manuscript.