1. Introduction
On March 9
th, 2020, the Italian Government issued a DPCM (
Decreto del Presidente del Consiglio dei Ministri, Decree of the President of the Council of Ministers) introducing the first full scale lockdown (LD) in the world, outside China, due to the Covid-19 outbreak [
1]. The first LD was extremely strict and effectively prohibited non-essential activities, thus causing most of them to come to a complete stop or be vastly reduced. The LD persisted until May 18
th, when a second DPCM lifted most restrictions [
2]. In the following year and a half, other minor LDs were introduced, though none of them were as strict as the first, thus making that LD a unique circumstance for the evaluation of atmospheric concentrations of select parameters. In fact, following similar LDs throughout the globe, researchers quickly performed studies on the unprecedented environmental conditions brought by LD measures [
3,
4,
5,
6,
7,
8,
9,
10].
The environmental response to reduced anthropic activities was not constant among the main parameters, greenhouse gases (GHGs), aerosols, and other key compounds observed globally across international networks such as the WMO/GAW (World Meteorological Organization – Global Atmosphere Watch). For example, despite generally lower anthropogenic emissions [
11] and improved air quality [
12,
13], methane experienced a global surge in 2020 [
14,
15]. Stevenson et al. (2022) [
16] reported a 3.8 to 5.8 ppb increase in the 2020 annual growth rate of methane linked to the large reduction in anthropogenic NO
x (nitrogen oxides) releases to the atmosphere, which was partially counterbalanced by fewer CO (carbon monoxide) and non-methane volatile organic compounds (NMVOC) emissions; the net increase, accounting for NO
x, CO, and NMVOC counterbalance effects, has been reported as 2.9 (1.7 to 4.0) ppb. NO
x compounds are known to be key regulating factors of methane [
17,
18], accounting for an approximate 1000 ppb total reduction in atmospheric concentration levels [
19]. It is also worth noting that LDs have led to an increase in methane emissions related to the energy sector [
20]. In the first year of the Covid-19 pandemic, a 5.3-5.5 ppb increase in the annual methane growth rate up to the value of 15.0 ppb yr−1 has been observed [
21]. Locally, the first multi-year evaluation of LMT methane cycles and trends shown in D’Amico et al. (2024a) [
22] highlighted a 2020 surge which was in accordance with the global trend observed by NOAA (National Oceanic and Atmospheric Administration) [
21]. Conversely, several studies reported a decline in other compounds such as CO, as a direct consequence of LD restrictions and reduced anthropogenic emissions [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24].
Research studies in Italy primarily focused on the impact of LD restrictions on atmospheric particle concentrations [
25,
26], urban air quality [
27,
28,
29], specific GHGs in a variety of environments [
30,
31,
32], multiparameter evaluations [
33,
34]. Globally, a review by Addas and Maghrabi (2021) [
35] reported that most studies published within one year from the first LDs focused, in order, on NO
2, PM
2.5, PM
10, SO
2, and CO as key indicators of changes in air quality due to LD restrictions. In the context of Europe, as early as 2020 several studies reported increase in air quality throughout the continent [
36,
37,
38]. Following research provided information on reduced emissions on a continental scale, in particular with respect to nitrogen compounds [
39,
40,
41].
With respect to the Lamezia Terme (LMT) observation site in Calabria, Southern Italy, a previous study assessed the effect of LD restrictions on nanoparticle concentrations via a comparison with data gathered at Lecce (Apulia, Southern Italy) [
42]. However, no other research has ventured deeper into a multiparameter evaluation of the first LD period at LMT, accounting for multiple gases and aerosols. This study is therefore a first attempt at analyzing CO, CO
2, CH
4, BC, and NO
x trends as well as key meteorological data observed at LMT during that period and provide new insights with respect to source apportionment in an area characterized by multiple anthropogenic and natural emission outputs [
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43].
Each of the analyzed parameters is an indicator of certain sources. Carbon monoxide (CO) is a common byproduct of combustion processes and its atmospheric trends have been frequently linked to the application of sustainable policies and new technologies to combustion engines [
44]; carbon dioxide (CO
2) is a primary output of fossil fuel burning and has been known for decades to be a primary driver of climate change [
45,
46,
47,
48,
49,
50]; methane (CH
4) is also a byproduct of the fuel combustion, for instance, of regular vehicles [
51] and airplane engines [
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43,
44,
45,
46,
47,
48,
49,
50,
51,
52], though natural outputs such as wetlands [
54] and specific anthropogenic emissions such as livestock [
55] are responsible for a considerable fraction of the annual global output; eBC (equivalent black carbon) is an excellent indicator of combustion processes [
56], known to pose both health [
57] and climate [
58,
59,
60,
61] hazards; nitrogen oxides (NO
x) also serve as indicators of anthropogenic activities, with various studies highlighting their reduction during LDs [
62,
63]. Fossil and biomass burning, as well as the use of fertilizers in agriculture, are among the main anthropogenic sources of atmospheric NO
x [
64,
65].
Overall, the analysis of these parameters, combined with the peculiar configuration of the observation site and local wind circulation, are expected to provide new insights on the first 2020 LD in Southern Italy, which in turn can be used to gather new information on local source apportionment and, by extension, the development of new sustainable policies at a local scale.
This research is organized as follows: section 2 will describe and characterize the observation site of LMT; section 3 will show the results of performed analyses; sections 4 and 5 will be focused on discussion and the conclusions of this study.
4. Discussion
In this study, Italy’s pioneeristic role in introducing a nationwide lockdown (LD) [
1] during the Covid-19 pandemic in early 2020 and Lamezia Terme’s (LMT) coastal location (
Figure 1) and local wind circulation (
Figure 2) have both been exploited to assess CO, CO
2, CH
4, eBC and NO
x parameters between February and July 2020, a time span exceeding the first nationwide LD (March 9
th – July 31
st) in order to account for pre- and post-LD trends and late environmental responses to LD restrictions. The first Italian LD, which was more strict than similar measures meant to contain Covid-19 in late 2020 and in 2021, therefore provided a unique circumstance to analyze key parameters with reduced anthropic activities. Among the observed parameters, CH
4 alone has so far been subject to a detailed multi-year study accounting for most of LMT’s observation history [
22], while the others are as of today not characterized on an observation site level.
Generally speaking, the Mediterranean Basin is considered as a hotspot for air quality assessments [
81] as well as a context where multiple air mass transport mechanisms combine [
82,
83,
84]. The first LD occurred during a seasonal transition between Mediterranean winter and spring, then summer, and resulted into most anthropic activities either being completely interrupted or significantly reduced. The most documented reduction in anthropic activity in the LMT area is that of the nearby international airport (
Figure 1), where only one daily flight to/from Rome-Fiumicino (IATA: FCO; ICAO: LIRF) operated during the LD to ensure the most basic connections to air travel networks [
85]. The LD has also affected other activities, such as the education sector, in a country where schools are normally open from Monday to Saturday and urban rush hours in the Lamezia Terme municipality area are linked to that form of commuting. The A2 highway (
Figure 1) has also experienced a sharp reduction in vehicular traffic, though the toll free policy applied to it does not allow a detailed estimate on traffic reduction.
The seasonal transition during which the first LD has occurred allowed to test several correlations with parameters usually linked to domestic heating and, consequently, daily temperatures. In fact, although a full-scale census and assessment on the phenomenon does not exist, many households in the rural areas nearby rely on biomass burning (in particular, wood) as domestic heating during the winter season.
Figure 3 confirms the seasonal change, which has resulted in a nearly constant increase in daily temperatures since mid-April, with the exception of a heat wave slightly before the restriction lift, which occurred on May 18
th [
2]. This allowed a direct comparison between the above-mentioned daily temperatures and daily-aggregated data concerning key parameters, as seen in
Figure 4 and
Table 2, which allowed to determine two different types of response to rising temperatures. In particular, CO (
Figure 4A), eBC (
Figure 4D), and NO
x (
Figure 4E), as shown by 90% confidence ellipses and R
2 values reported in
Table 2, are characterized by a correlation between daily temperatures and observed concentrations. CO in particular seems the most affected parameter, as winter concentrations are higher due to more emissions from domestic heating as well as lower OH (hydroxyl radical) concentrations, the latter being a notable sink [
86]. Conversely, CO
2 (
Figure 4B) and CH
4 (
Figure 4C) show no correlation, as the range of daily concentrations does not seem to be affected by temperature changes and the R
2 values reported in
Table 2 are low.
Different trends among observed parameters have also been reported in the evaluation of hourly averages throughout the entire observation period (FEB-JUL 2020), as shown in
Figure 5. In the case of CO (
Figure 5A), hourly data have added extra detail to the correlation with daily temperatures, as peaks fall considerably from mid-April onwards, reflecting an increase in daily temperatures (
Figure 3). This is also observed in the case of eBC (
Figure 5D), though the pattern is not as prominent as that observed for CO. CO
2 (
Figure 5B) and CH
4 (
Figure 5C) both lean towards a stabilization of observed hourly trends but show an upward trend in the PLD period. CH
4 in particular yields the notable hourly peaks in the PLD period, a finding that is not in accordance with the general seasonal trend of CH
4 at LMT reported in D’Amico et al. (2024a) [
22] but does reflect the peak experienced by this compound in 2020. Previous studies linked CH
4 concentrations and their relative stability to local sources, such as landfills (
Figure 1), waste management, and livestock [
43], which are all believed to have been constant during the LD. NO
x have experienced a decline (
Figure 5E), though peaks in hourly averages are present during the LD.
Following the daily cycle analysis seen in D’Amico et al. (2024a) [
22], ALD, LD, and PLD trends have been analyzed for the purpose of this research study, with results shown in
Figure 6. At LMT, the daily cycle is heavily influenced by local wind circulation (
Figure 2): daytime flows generally come from the sea and yield low values in all parameters, while nighttime flows come from the northeast and are enriched. Unlike the trends seen in D’Amico et al. (2024a),
Figure 6 show different responses of daily cycle between ALD, LD, and PLD periods: CO (
Figure 6A) mole fractions are very low during the PLD period, but follow similar trends during the ALD and LD periods and nighttime peaks are consistent with previous research which attributed them to stable layer conditions [
43]; CO
2 (
Figure 6B) and CH
4 (
Figure 6C) both show the actual daily cycle expected at LMT due to local wind circulation, with ALD values being generally lower than their LD and PLD counterparts; eBC (
Figure 6D) and NO
x (
Figure 6E) both show a daily cycle affected by a perturbance in the early morning, which is linked to the transition between NE-continental and W-seaside winds. The observed NO
x pattern in early morning hours is compatible with the findings of a previous study [
43], which linked them to rush hour emissions within the shallow coastal PBL. In fact, the ALD peak is significantly greater than LD and PLD concentrations observed at the same hours, highlighting a reduction in anthropic activities (specifically, transportation).
The influence of local wind circulation was further analyzed in
Figure 7.
Table 3 shows that, once the filters are applied, the W to NE ratio in terms of gathered hourly data is 1.67-1.71. D’Amico et al. (2024a) [
22], in analyzing CH
4 data at LMT, already reported a clear differentiation between W and NE mole fractions. The same study also reported a correlation between NE concentrations and wind speed, with lower speeds being correlated with the highest concentrations observed at the site, attributable to nearby sources. The study found evidence of a “hyperbola branch pattern” (HBP) in wind speeds and mole fractions. This finding from previous research is not only confirmed for CH
4 in the context of the 2020 LD period but is also extended to all other observed parameters (CO, CO
2, eBC and NO
x) and allows to better constrain some of the local sources of pollution. CO (
Figure 7A1, 7A2, 7A3) shows a clear seasonal trend, with ALD and PLD mole fractions being well differentiated; the main source of CO is from the NE sector and shows little influence of wind speeds, which may point to remote sources. The highest reported values from the W sector have been observed during the LD at low to moderate wind speeds, which also point to contributions from other locations in the context of the Mediterranean Basin. CO
2 (
Figure 7B1, 7B2, 7B3) shows a clear W/NE differentiation, with LD data yielding generally low values regardless of wind speed, though peaks linked to the LD period are reported. CH
4 (
Figure 7C1, 7C2, 7C3) shows patterns similar to those of CO
2, which is consistent with other evaluations in this research paper, as well as finding from previous research [
43], and may indicate contributions from remote and local sources alike. eBC (
Figure 7D1, 7D2, 7D3) shows a seasonal trend, plus a peculiar behavior during the LD period: contributions from nearby and remote locations are prominent in the NE sector alone. NO
x (
Figure 7E1, 7E2, 7E3) concentrations show reduced influence from wind speeds, which results into high values even from the W sector, which previous research attributed to offshore ship emissions related to traffic to and from the Gioia Tauro port, located ≈55 km S-SW of LMT, which is a key logistical hub in the central Mediterranean [
43].
A final comparison between ALD, LD, and PLD data has been performed on a weekday
latu sensu (Monday to Sunday) basis, as shown in
Figure 8, following the findings from D’Amico et al. (2024a, 2024b) [
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43,
44,
45,
46,
47,
48,
49,
50,
51,
52,
53,
54,
55,
56,
57,
58,
59,
60,
61,
62,
63,
64,
65,
66,
67,
68,
69,
70,
71,
72]. By averaging gathered hourly data per weekday, anthropic influences have been tested under the assumption that no natural phenomenon would lead to a proper weekly cycle, unlike anthropogenic emissions which are subject to such cycles [
87,
88,
89,
90]. CO (
Figure 8A) shows no LD and PLD weekly cycles, while a weekly cycle during the ALD period is present and is consistent with changes in anthropic activities across the week, a pattern that is believed to have been reduced to a bare minimum during the LD. The LD period does not show a proper weekly cycle for CO, which is consistent with constant domestic heating emissions during the lockdown; in addition to that, the PLD period does not show a weekly cycle, which is also consistent with a shift from CO emissions from domestic heating (and comparable sources of emission) to wildfires, which are assumed to be spread equally over the course of a standard week. In the case of CO, fluctuations in daily concentrations differ up to ≈50%. CO
2 (
Figure 8B), CH
4 (
Figure 8C), eBC (
Figure 8D), and NO
x (
Figure 8E) all show well defined weekly changes in the ALD period. In the case of CO
2, it is worth noting that the observed fluctuations are in the order of a few ppm around the average of 427.5 ppm and are therefore not comparable to the major fluctuations seen for other parameters. CH
4 shows a clear pattern in the ALD period, which is consistent with a weekly cycle of anthropogenic emissions; LD and PLD CH
4 weekly behaviors tend towards a normalization, which supports the hypothesis of a constant, or nearly constant, output from sources such as landfills and livestock. eBC’s flattening in the PLD period is consistent with a shift from anthropogenic emissions to other outputs such as wildfires (which can be of anthropic origin, but do not have a clear weekly cycle). NO
x has a prominent ALD cycle, with weekly fluctuations reaching a peak of ≈250%, a threshold not observed in any of the other parameters.
In addition to weekly fluctuations, averaged values per weekday (dotted horizontal lines) in
Figure 8 also highlight the differences between absolute ALD, LD, and PLD concentrations, which reflect clear seasonal changes in the case of CO and eBC and major shifts in emission sources. Both parameters have much lower PLD averages compared to their LD and ALD counterparts.
An assessment of the environmental response to LDs also needs to consider restrictions and regulations enforced in other countries. As described in section 2.1, in fact, LMT’s location in the central Mediterranean area makes the site subject to the influence of several European and African outputs [
74]. The Mediterranean Basin itself is also known to be subject to air masses originating in continental Europe and enriched in various pollutants [
84,
85,
86,
87,
88,
89,
90,
91]. At LMT, the W sector is assumed not to have physical obstacles for hundreds of kilometers, as the Italian island of Sardinia is in fact located ≈600 km W-NW from LMT. Continental Spain and France are located ≈1300 and ≈1000 km in those directions, respectively. In the case of the NE corridor however, it is worth noting that both Greece and Albania are in a ≈350 km radius from the LMT observation site in the E/NE direction, meaning that at least some of the observed emission outputs could be attributable to such “remote” regions as well as other countries in the Balkans. For instance, Greece had a different LD policy in early 2020 compared to Italy (March 23
rd – May 4
th in Greece, 42 days [
92]) which does not overlap with the definitions of ALD, LD, and PLD used in this paper.
Overall, a cross analysis of the ALD, LD, and PLD periods at LMT has allowed to better constrain, for the first time in its operational history, the seasonal cycles of parameters such as CO and eBC in conditions with extremely low anthropogenic emissions, as well as provide new insights on the sources of CO2, CH4, and NOx emissions in the NE sector. These new findings significantly expand the knowledge on source apportionment at LMT and provide baseline data for future research relying on additional atmospheric tracers.
5. Conclusions
For the first time, the first Italian lockdown (LD) of 2020 has been used to assess the behavior of CO (carbon monoxide), CO2 (carbon dioxide), CH4 (methane), eBC (equivalent black carbon), and NOx (nitrogen oxides, NO + NO2) at the WMO/GAW regional coastal site of Lamezia Terme in Calabria, Southern Italy. The study exploited the station’s location in the context of the Mediterranean Basin, where local wind circulation patterns allow to discriminate western-seaside corridors yielding generally low concentration values with northeastern-continental corridors, enriched in most pollutants and agents of climate change. The study considered a period longer than the actual first nationwide LD, which was introduced on March 9th and consequently lifted on May 18th. Data gathered between February and July therefore allowed to assess the behavior of observed parameters before (ALD), during (LD), and after (PLD) the first lockdown.
The research has allowed to better constrain the correlation between key parameters and seasonal changes, as the observation period is characterized by a nearly constant increase in daily temperatures. CO and eBC have proven to be the most affected by seasonal changes, thus confirming that local sources of emissions may in fact be related to domestic heating and similar outputs.
Other evaluations considered hourly averages and, particularly, their patterns in daily cycles. The analysis of these cycles, differentiated by ALD, LD, and PLD, allowed to better constrain local and potentially remote emission sources, as well as the influences of wind circulation on observed data. Consequently, the integration of wind data (direction and speed) has also allowed to find notable differences between the W and NE corridors, with the former being depleted in most parameters across the entire observation period. The NE corridor is characterized by the same hyperbola branch pattern (HBP) observed in a previous study, with low wind speeds being correlated with high concentrations and, vice versa, high speeds yielding low values. The finding has allowed to identify local sources in the NE sector as responsible for most of these peaks.
Finally, following the implementation of new methods seen in previous research, all parameters have been evaluated with respect to their weekly cycles, under the assumption that no natural phenomenon (unlike anthropic activities) would lead to a proper weekly cycle and trend. This analysis has proved the ALD period to be affected by such cycles, while the LD and PLD periods were much less affected, further corroborating the hypothesis proposed in other research by which anthropogenic emissions are responsible for a significant fraction of the peaks observed at LMT station.
Author Contributions
Conceptualization, F.D. and C.R.C.; methodology, F.D., C.R.C., D.G., and T.L.F.; software, F.D.; validation, C.R.C., I.A., D.G., T.L.F. and P.C.; formal analysis, F.D.; investigation, F.D.; data curation, F.D., I.A., D.G., E.A., T.L.F., P.C., L.M., D.P., S.S. and G.D.B.; writing—original draft preparation, F.D.; writing—review and editing, F.D., CR.C., I.A., D.G., E.A., T.L.F., M.D.P., P.C., L.M., D.P., S.S. and G.D.B.; visualization, F.D., C.R.C., D.G. and T.L.F.; supervision, C.R.C. and P.C.; funding acquisition, C.R.C. and M.D.P. All authors have read and agreed to the published version of the manuscript.