1. Introduction
Air pollution has been one of Europe’s leading political concerns since the late 1970s. It continues to significantly impact the health of the European population today, particularly in urban areas. In terms of harm to human health, Europe’s most serious pollutants are particulate matter (PM), NO2, and ground-level O3. Italy is among the European countries most affected by air pollution, with 58,600 premature deaths attributed to PM2.5 exposure, 14,600 to NO2 exposure, and 3000 to O3 exposure [1].
Meteorological parameters play an important role in determining air-pollution concentrations. Generally, PM concentrations decrease with an increase in precipitation rate, wind speed, and temperature [2]. The NO2 in our atmosphere is primarily a function of the magnitude of nitrogen oxide (NOx) emissions and weather factors, such as sun angle, wind speed, and temperature. Meteorological variations between years can cause column NO2 differences of ~15% over monthly timescales [3]. O3 generally increases with increasing temperature and decreases with increasing relative humidity [4].
Road traffic is the leading cause of air pollution in cities, and it is responsible for almost a quarter of Europe’s greenhouse gas emissions [5]. According to a recent report published by Legambiente, a recognized environmental organization in Italy, each year, several Italian cities exceed air-pollution limits, especially PM and O3 pollutant levels [6]. In 2020, as in previous years, Italy’s largest cities were forced to ban hundreds of thousands of vehicles from the road after days of persistent smog.
Despite the societal and economic impacts of the COVID-19 pandemic [7,8], the lockdown measures put in place by the Italian government provided an unprecedented opportunity to increase our knowledge on the contribution of transportation and industrial emissions on air quality in our cities. As expected, after Italy adopted the mandatory measures, there was a significant reduction in transport- and industry-related emissions. Recent research clearly discloses that the restrictive measures adopted during March–April 2020 brought about a significant reduction (−64.6%) in Rome’s personal vehicle usage [9].
On the evening of 9 March 2020, Italian Prime Minister Giuseppe Conte imposed a national quarantine in response to the growing COVID-19 pandemic restricting people’s movement, except for reasons related to basic human needs, work, and health [10]. On 22 March 2020, a new decree closed down all non-essential industries throughout the country and restricted inter-city movement further by requiring travelers to provide justification and documentation to the authorities, when moving between and within cities [11]. The restrictions lasted three months up to 3 June 2020.
This study aimed to assess the effects of the reduced emissions during the COVID-19 period on the air quality in three Italian cities, Florence, Pisa, and Lucca. For this study’s purpose, we compared the concentrations of PM10, PM2.5, NO2, and O3 measured before, during, and after the lockdown, with the corresponding values.
2. Methods
2.1. Area of the Study
We identified three cities for this study (Figure 1), all located in Tuscany’s northern region:
-
Florence, the capital city of the Tuscany region and the Province of Florence, is the most populated, with about 372,000 inhabitants that live in the municipality, but over 1,520,000 that live in the greater metropolitan area;
-
Pisa, the capital city of the Province of Pisa, has over 91,000 residents living in the municipality and about 200,000 living in the surrounding area;
-
Lucca, the capital city of the Province of Lucca, which has more than 90,000 residents.
Each of these cities has problems related to air quality, with the measured pollutant values often exceeding EU law’s daily limits. In contrast, the selected cities differ amongst several characteristics, such as the number of inhabitants (
2.2. Air-Quality and Meteorological Data Collection and Processing
We based the data used in the following analyses from ARPAT (Agenzia Regionale per la Protezione Ambientale della Toscana, the Regional Agency for the Environmental Protection of Tuscany) website (
We collected the data from the three cities (Firenze, Lucca, and Pisa). For each city and each pollutant, we specify the number of measuring stations in Table 1. The number of stations in each city varies according to the contaminant measured and the actual city. ARPAT provided data verification and validation by monitoring the instrumental performance and applying quality-control procedures.
For each pollutant, from 1 January 2019, we focused on the two standard periods dividing it further into the following subperiods:
-
[1 January–8 March 2019] vs. [1 January–8 March 2020] → pre-lockdown period;
-
[9 March–3 June 2019] vs. [9 March–3 June 2020] → lockdown period;
-
[4 June–12 August 2019] vs. [4 June–12 August 2020] → post-lockdown period.
We used a paired t-test to determine whether the mean difference between the two sets of observations was significant. The null hypothesis (H0) assumed that the true mean difference between the paired samples was zero, while the alternative hypothesis (H1) assumed that the true mean difference between the paired samples was not equal to zero. We chose the cutoff value of 0.05 to determine statistical significance. This value corresponds to a 5% (or less) chance of obtaining a result like the one that was observed if the null hypothesis was true. Regarding the main assumptions of the paired t-test, it is reasonable to assume that the observations were independent of one another. Furthermore, we tested the assumption of normality and found that the data were approximately normally distributed. Finally, we examined the relative percentage variation of the average concentrations throughout the subperiods.
For this analysis, we downloaded data on wind speed, rainfall, relative humidity, temperature, and solar irradiance for the entire period from the Meteoblue AG–Switzerland website (
To give further details concerning the meteorological data, we used the following data measured at the Pisa San Giusto, Firenze Peretola, and Lucca meteorological stations: average daily air temperature (as Celsius degrees), average daily air humidity (as a percentage value), average daily wind speed (as km/h), and average daily sea-level air pressure (as mbar). For both 2019 and 2020, we did not derive statistically significant differences. In particular, we considered the wind speed as average daily wind speed and not the wind directions during each period for the two years, 2019 and 2020, since the measuring stations are distant from each other, and since we are not interested in the directions in which the pollutants are spread but only in their values as measured by every measuring station. In all the three cases of Pisa San Giusto, Firenze Peretola, and Lucca, we found that the average daily wind speed differences between the same periods of the two years 2019 and 2020 were not statistically significant at a α = 0.05 level.
3. Results and Discussion
3.1. Particulate Matter
Compared to the same period in 2019, PM10 and PM2.5 concentrations during the pre-lockdown period are not significantly different between the air-monitoring stations across the three cities, Florence, Pisa, and Lucca. Table 2 and Table 3 show that the FI-GRAMSCI monitoring station is the only air-monitoring station in which there was a statistically significant reduction of PM10 and PM2.5 during the lockdown period, with a variation of about 30.8% and 50.1%, respectively. On the other hand, there was no significant difference between the remaining air-monitoring stations. This result is in line with a recent report on outdoor air pollution changes in Scotland during the COVID-19 lockdown, which showed that PM concentrations did not decline, despite reducing vehicular traffic [13]. It is very likely that the environmental response, measured as the reduction in PM concentrations, varies according to the dominant source of emission in each country and because of the specific meteorological conditions [14].
Unexpectedly, we also observed an increase in the PM2.5 concentration at the PI-BORGHETTO monitoring station, the Pisa city center traffic station, with a variation of about 33.3%. Changes in meteorological conditions do not explain these variations, given the absence of a significant difference between the two years. We hypothesize that because people spent more time indoors during the lockdown period, the higher levels of PM2.5 are associated with the increased use of heating systems in buildings. Moreover, we hypothesize there was only a small change in the volume of traffic in this area of the city. This hypothesis highlights that other PM sources are significantly present and that holistic source-control measures are needed for improved air quality in urban environments [15].
Figure 2 and Figure 3 graph the PM10 and PM2.5 daily mean concentrations for the three stations in which we note a statistically significant difference. Although previous studies reported a decrease of PM10 and PM2.5 concentrations [16,17,18,19], we noticed that they had compared the two periods of the same year: the pre-lockdown and lockdown period. In our view, this methodology represents a bias since, every year, it is possible to observe a decrease of particulate matter concentrations due to seasonal variations [20]. However, other studies reported a reduction in PM, using a methodology similar to ours, for example, by comparing the same periods of the previous year or with averages of several prior years [21,22,23,24]. These differences should be interpreted in light of the aforementioned variations in meteorological conditions and sources of emissions between different countries and cities. Regarding the FI-GRAMSCI measuring station, we must consider that it is the critical point (hot spot) in Tuscany, with much higher values than any other air-monitoring station. Therefore, this result should be interpreted with this peculiarity in mind and confirmed with further analyses by comparing this station with similar hot-spot air-monitoring stations in other cities than Pisa, Lucca, and Florence itself. Overall, our results show no direct evidence of a relationship between the implementation of the lockdown measures and the reduction of particulate matter in urban centers, except heavy traffic areas. Further research is needed to clarify this relationship, and it would also be interesting to examine the citizens’ risk perception of air pollution in the context of the current Covid-19 emergency [25].
3.2. Nitrogen Dioxide
NO2 primarily gets in the air from burning fuel from cars, trucks, buses, power plants, and off-road equipment. It reacts with other chemicals in the air to form both PM and O3, and it is harmful when inhaled due to the effect it has on the respiratory system. While we noted fluctuating NO2 levels in different areas, the concentration of NO2, during the lockdown was significantly less among all the air-monitoring stations across the cities included in this study. Table 4 indicates statistically significant reductions (p < 0.05) of NO2 concentration during the lockdown period, compared to the same period of 2019. These results clearly indicate that reducing traffic emissions contributes to lower NO2 concentration levels and represents a significant improvement in public health. They should motivate us to adopt new urban mobility policies to reduce pollution in our cities and protect human health even after the COVID-19 crisis [26].
Figure 4 exhibits the daily mean of NO2 concentration during the three periods.
In detail, compared to the levels in 2019, concentrations of NO2 decreased by 38.5%, 32.1%, 39.4%, 40.1%, 41.6%, 35.0%, and 39.4%, respectively, at the FI-GRAMSCI, FI-BASSI, FI-MOSSE, PI-BORGHETTO, PI-PASSI, LU-MICHELETTO, and LU-CARIGNANO measuring stations. These results are consistent with the findings of previously published works. For example, Nakada et al., 2020 [22], observed a decrease of up to 54.3% in NO2 concentrations (μg∙m−3) in the urban roads of São Paulo state, Brazil. Moreover, Collivignarelli et al., 2020, observed a drastic drop in the concentration of NOx and NO2 in all the areas covered by that study on air quality in Milan, Italy. Researchers also noted an apparent decrease in NO2 levels in Rio de Janeiro, Brazil [19], in 22 cities covering different regions of India [21], the megacity Delhi, India [18], the Yangtze River Delta Region [23], and 44 cities in Northern China [24]. Satellite measurements also capture regional NO2 concentration reductions in air quality by the TROPOspheric Monitoring Instrument (TROPOMI) on the European Space Agency’s Sentinel-5 satellite [27].
As shown in Table 5, we must point out that the NO2 concentration during 2020, is statistically lower than that of 2019, even before the lockdown. However, the variations during the lockdown are higher, leading us to deduce that the decrease in vehicular traffic positively affected the air quality. Road traffic is the principal outdoor source of nitrogen dioxide [28], and the relationship between traffic and NO2 has been examined in many studies [29,30,31].
3.3. Ozone
Tropospheric O3 is a harmful air pollutant created by chemical reactions between NOx and volatile organic compounds (VOC) in the presence of sunlight. A recent study showed that the lockdown measures caused an increase of O3 concentrations at urban stations of four Southern European cities and Wuhan (China). This phenomenon’s explanation is mainly based on the unprecedented reduction in NOx emissions leading to a lower O3 titration by NO [32]. On the other hand, our results show that the O3 concentration did not show significant reductions during the period of lockdown (Table 6). Figure 5 shows the increasing ozone trend across the three periods, a typical phenomenon observed during the spring season due to the higher solar radiation [33]. However, we must point out that our data are related to background and suburban stations, and, for this reason, we cannot compare them with the aforementioned study, and we cannot draw any significant conclusions about O3 general trends in urban areas.
4. Conclusions
In this article, we assessed the impact of the Italian nationwide lockdown due to the COVID-19 outbreak on air quality in three medium-sized cities. The strict travel restrictions and limitations on the movement between and within the cities during the Italian lockdown period provided an unprecedented opportunity to assess anthropogenic activities’ effect on urban air quality. More specifically, the adopted restrictive measures resulted mainly in a significant reduction in vehicular traffic and industrial activities. We aimed to evaluate the effects of these restrictive measures on the concentration of four primary air pollutants collected by the Regional Agency for the Environmental Protection of Tuscany, which are PM10, PM2.5, NO2, and O3. Our results did not show significant reductions in PM levels during the lockdown period, except at one monitoring station in an area of heavy traffic, which is considered a hot spot because of the consistently high levels recorded for all of Tuscany. These results are not consistent with the published studies on this topic, and we can assume that the environmental response varies according to the dominant emission source and the specific meteorological conditions. Our research field should further explore the interrelationship between the different sources of emissions and particulate matter concentration. On the other hand, the reduction in NO2 pollution levels, consistent with other studies, was statistically significant at all the air-monitoring stations across the cities used in this study, showing a relevant traffic volume relationship. These results should motivate the politicians to adopt new urban policies to reduce pollution in our cities and protect human health even after the COVID-19 crisis. Finally, for O3 pollutant levels, we did not observe a significant reduction during the lockdown period. However, we should take into account that none of these monitoring stations are placed in the traffic areas, and, therefore, we cannot draw significant conclusions on the trend of O3 during the lockdown period. In conclusion, we believe it is necessary to adopt holistic source-control measures for improved air quality in urban environments.
Author Contributions
Conceptualization, M.M.M.S.-V. and G.D.; methodology, G.D., L.C., M.C., and A.L.M.; software, G.D., L.C., M.C., and A.L.M.; validation, G.D., L.C., and A.L.M.; formal analysis, G.D., L.C., and A.L.M.; G.D., L.C., and A.L.M.; writing—original draft preparation, G.D., L.C., and A.L.M.; writing—review and editing, M.M.M.S.-V. and G.D.; supervision, M.M.M.S.-V. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Acknowledgments
We wish to thank Maria Grazia De Agazio for proofreading, editing, and assisting with translating the English manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figures and Tables
Figure 2. The daily mean concentration of PM10 (μg∙m−3). In each graph, the lockdown period is depicted between the two vertical bars.
Figure 3. The daily mean concentration of PM2.5 (μg∙m−3). In each graph, the lockdown period is depicted between the two vertical bars.
Figure 4. The daily mean concentration of NO2 (μg∙m−3). The area between the vertical bars in each graph indicates the duration of the lockdown.
Figure 5. The daily mean concentration of O3 (μg∙m−3). The area between the vertical bars in each graph indicates the duration of the lockdown duration.
Characteristics of the air-monitoring stations included in the study.
Air-Monitoring Station | City | Station Type | Area Type | Pollutants | |||
---|---|---|---|---|---|---|---|
PM10 | PM2.5 | NO2 | O3 | ||||
FI-GRAMSCI | Florence | Traffic | Urban | ✓ | ✓ | ✓ | |
FI-BASSI | Florence | Background | Urban | ✓ | ✓ | ✓ | |
FI-MOSSE | Florence | Traffic | Urban | ✓ | ✓ | ||
FI-SETTIGNANO | Florence | Background | Suburban | ✓ | |||
PI-BORGHETTO | Pisa | Traffic | Urban | ✓ | ✓ | ✓ | |
PI-PASSI | Pisa | Background | Urban | ✓ | ✓ | ✓ | ✓ |
LU-MICHELETTO | Lucca | Traffic | Urban | ✓ | ✓ | ||
LU-CARIGNANO | Lucca | Background | Rural | ✓ | |||
LU-SAN-CONCORDIO | Lucca | Background | Urban | ✓ | ✓ |
Station type. Traffic-measuring stations: air-monitoring stations located in areas where the pollution level is most influenced by traffic emissions from neighboring roads with medium-high traffic intensity. Background (or general) measuring stations: air-monitoring stations located where the pollution level is not influenced mostly by emissions from specific sources (industries, traffic, residential heating, etc.) but rather from a combination of all sources, upwind from the station concerning the predominant wind directions at the site. Area type. Urban fixed sampling site: a fixed site placed in a growing or predominantly built-up area. Suburban fixed sample site: a fixed site placed in largely built-up areas where there are both built-up and non-urbanized areas. Rural fixed sampling site: a fixed site selected in areas other than those identified for urban and suburban locations. The site is defined as remote rural if it is more than 50 km away from the emission sources. PM = particulate matter.
Table 2Output of the two-sample t-test to evaluate the lockdown effects on PM10 air concentration. The PM10 values are expressed in μg∙m−3.
First Period |
Second Period (Lockdown) |
Third Period |
|||||||
---|---|---|---|---|---|---|---|---|---|
Mean |
Mean (SD) |
p-Value | Mean (SD) |
Mean (SD) |
p-Value | Mean (SD) |
Mean (SD) |
p-Value | |
FI-GRAMSCI | 37 |
36 |
0.635 | 23 |
16 |
<0.001 | 24 |
17 |
<0.001 |
FI-BASSI | 23 |
28 |
0.300 | 14 |
15 |
0.186 | 19 |
15 |
0.002 |
FI-MOSSE | 30 |
32 |
0.427 | 14 |
15 |
0.269 | 21 |
15 |
<0.001 |
PI-BORGHETTO | 39 |
36 |
0.237 | 19 |
18 |
0.255 | 26 |
18 |
<0.001 |
PI-PASSI | 35 |
33 |
0.223 | 16 |
17 |
0.982 | 22 |
15 |
<0.001 |
LU-MICHELETTO | 42 |
42 |
0.996 | 17 |
19 |
0.141 | 20 |
18 |
0.131 |
LU-SAN-CONCORDIO | 40 |
37 |
0.322 | 15 |
18 |
0.037 | 21 |
15 |
<0.001 |
Output of the two-sample t-test to evaluate the lockdown effects on PM2.5 air concentration. The PM2.5 values are expressed in μg∙m−3.
First Period |
Second Period (Lockdown) |
Third Period |
|||||||
---|---|---|---|---|---|---|---|---|---|
Mean (SD) |
Mean |
p-Value | Mean (SD) |
Mean |
p-Value | Mean (SD) |
Mean |
p-Value | |
FI-GRAMSCI | 24 |
21 |
0.207 | 12 |
10 |
0.045 | 14 |
10 |
<0.001 |
FI-BASSI | 19 |
20 |
0.552 | 9 |
10 |
0.056 | 12 |
9 |
<0.001 |
PI-BORGHETTO | 25 |
22 |
0.179 | 8 |
11 |
<0.001 | 11 |
8 |
<0.001 |
PI-PASSI | 29 |
26 |
0.245 | 11 |
12 |
0.437 | 15 |
10 |
<0.001 |
Output of the two-sample t-test to evaluate the effects of lockdown on NO2 air concentrations. The NO2 values are expressed in μg∙m−3.
First Period |
Second Period (Lockdown) |
Third Period |
|||||||
---|---|---|---|---|---|---|---|---|---|
Mean (SD) |
Mean (SD) |
p-Value | Mean (SD) |
Mean (SD) |
p-Value | Mean (SD) |
Mean (SD) |
p-Value | |
FI-GRAMSCI | 117 |
110 |
0.081 | 98 |
60 |
<0.001 | 94 |
76 |
<0.001 |
FI-BASSI | 61 |
60 |
0.856 | 36 |
25 |
<0.001 | 31 |
22 |
<0.001 |
FI-MOSSE | 89 |
79 |
0.005 | 59 |
36 |
<0.001 | 46 |
34 |
<0.001 |
PI-BORGHETTO | 85 |
76 |
0.003 | 58 |
35 |
<0.001 | 48 |
34 |
<0.001 |
PI-PASSI | 64 |
55 |
0.001 | 34 |
20 |
<0.001 | 24 |
16 |
<0.001 |
LU-MICHELETTO | 69 |
58 |
<0.001 | 45 |
29 |
<0.001 | 35 |
27 |
<0.001 |
LU-CARIGNANO | 77 |
62 |
<0.001 | 46 |
28 |
<0.001 | 36 |
26 |
<0.001 |
Variations NO2 (μg∙m−3) during the three periods.
First Period 2020 vs. 2019 |
Second Period 2020 vs. 2019 |
Third Period 2020 vs. 2019 |
|
---|---|---|---|
FI-GRAMSCI | 5.7% | 38.5% | 19.6% |
FI-BASSI | 1.0% | 32.1% | 28.3% |
FI-MOSSE | 11.3% | 39.4% | 26.2% |
PI-BORGHETTO | 9.7% | 40.1% | 29.4% |
PI-PASSI | 13.2% | 41.6% | 33.1% |
LU-MICHELETTO | 17.1% | 35.0% | 23.3% |
LU-CARIGNANO | 19.6% | 39.4% | 25.9% |
Output of the two-sample t-test to evaluate the effects of lockdown on O3 air concentrations. The O3 values are expressed in μg∙m−3.
Pre-Lockdown | Lockdown | Post-Lockdown | |||||||
---|---|---|---|---|---|---|---|---|---|
Mean (SD) |
Mean |
p-Value | Mean (SD) |
Mean (SD) |
p-Value | Mean (SD) |
Mean |
p-Value | |
FI-SETTIGNANO | 70 |
66 |
0.067 | 100 |
100 |
0.939 | 128 |
116 |
0.007 |
PI-PASSI | 61 |
62 |
0.886 | 92 |
95 |
0.198 | 107 |
99 |
0.008 |
LU-CARIGNANO | 80 |
70 |
<0.001 | 100 |
98.42 (17) | 0.423 | 129 |
108 |
<0.001 |
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020 by the authors.
Abstract
Despite the societal and economic impacts of the COVID-19 pandemic, the lockdown measures put in place by the Italian government provided an unprecedented opportunity to increase our knowledge of the effect transportation and industry-related emissions have on the air quality in our cities. This study assessed the effect of reduced emissions during the lockdown period, due to COVID-19, on air quality in three Italian cities, Florence, Pisa, and Lucca. For this study, we compared the concentration of particulate matter PM10, PM2.5, NO2, and O3 measured during the lockdown period, with values obtained in the same period of 2019. Our results show no evidence of a direct relationship between the lockdown measures implemented and PM reduction in urban centers, except in areas with heavy traffic. Consistent with recently published studies, we did, however, observe a significant decrease in NO2 concentrations among all the air-monitoring stations for each city in this study. Finally, O3 levels remained unchanged during the lockdown period. Of note, there were slight variations in the meteorological conditions for the same periods of different years. Our results suggest a need for further studies on the impact of vehicular traffic and industrial activities on PM air pollution, including adopting holistic source-control measures for improved air quality in urban environments.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details



1 Department of Preventive Medicine and Public Health, Food Sciences, Toxicology, and Legal Medicine, School of Pharmacy, University of Valencia, Avenida Vicente Andres Estellés s/n, Burjassot, 46100 Valencia, Spain;
2 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
3 Hygiene and Public Health Unit, Department of Public Health, AUSL Imola, 40026 Bologna, Italy;
4 Department of Preventive Medicine and Public Health, Food Sciences, Toxicology, and Legal Medicine, School of Pharmacy, University of Valencia, Avenida Vicente Andres Estellés s/n, Burjassot, 46100 Valencia, Spain;
5 Department of Preventive Medicine and Public Health, Food Sciences, Toxicology, and Legal Medicine, School of Pharmacy, University of Valencia, Avenida Vicente Andres Estellés s/n, Burjassot, 46100 Valencia, Spain;