Air quality and meteorological patterns of an early spring heatwave event in an industrialized area of Attica, Greece.

Heatwaves-excessively hot ambient conditions that are considered a serious threat to human health-are often associated with poor air quality. The aim of this study was to examine the impact of an early heatwave episode in an industrialized plain in the eastern Mediterranean region (Thriasio, Greece) on human thermal discomfort and urban air quality. The heatwave occurred in mid (15-20) May 2020, shortly after some of the restrictions that were improsed to halt the spread of coronavirus disease 2019 (COVID-19) in Greece were lifted (on 4 May). The discomfort index (DI) and the daily air quality index (DAQI) were calculated on an hourly basis throughout spring 2020 (March, April, May) using data from two stations that measure meteorological parameters and air pollutant concentrations in the Thriasio Plain. The analysis showed that the air temperature increased during 7-17 May to levels that were more than 10 °C above the monthly average value (25.8 °C). The maximum measured air temperature was 38 °C (on 17 May). The results showed a high level of thermal discomfort. The DI exceeded the threshold of 24 °C for several hours during 13-20 May. Increased air pollution levels were also identified. The average DAQI was estimated as 0.83 ± 0.1 and 1.14 ± 0.2 at two monitoring stations in the region of interest during the heatwave. Particulate matter (diameter < 10 μm) appeared to contribute significantly to the poor air quality. Significant correlations between the air temperature, DI, and AQSI were also identified. © Springer Nature Switzerland AG 2021.

Introduction

Globally, heatwaves are getting hotter, longer, and more frequent, particularly in regions affected the most by climate change. Moreover, these heatwave trends are accelerating (Perkins-Kirkpatrick and Lewis 2020). For instance, in the Mediterranean Basin, heatwave frequency increased by two days per decade from 1950 to 2017, and by 6.4 days per decade from 1980 to 2017 (Perkins-Kirkpatrick and Lewis 2020).

Extreme temperature events have a major impact on various aspects of human life, especially when these events are unexpected (Fontana et al. 2015; Sailor et al. 2019). Extreme temperatures may cause wildfires (Mavrakis and Salvati 2015) and marine heatwaves (Mavrakis and Tsiros 2018), can increase the energy demand for cooling (Añel et al. 2017) and affect productivity (Xia et al. 2018), are directly linked to loss of life (Xu et al. 2016), and indirectly affect the prevalence of infectious diseases (Naval et al. 2017). Several studies have examined heatwave and drought events through a unified approach (Mavrakis et al. 2015b, 2016). In many cases, these two environmental hazards are related, so it can be unclear whether effects of such events on human well-being are due to a short-duration heatwave or a longer-term drought. Previous studies focusing on the Thriasio Plain suggest that the local climatic regime has become warmer and drier (Mavrakis et al. 2015a, 2015b, 2016). This represents an emblematic example of the urban–rural relationship.

The Thriasio Plain is located close to the city of Athens—the capital of Greece and the main industrial hub of the country. This plain has the highest industrial activity, fuel consumption, and industrial pollution levels in Greece (Mavrakis et al. 2008). The region is characterized by the spontaneous development of various sectors, including residential areas, agriculture, and a range of industries such as oil refineries and the steel, steelworking, chemical, product, nonmetal mineral material, and cement industries. There are nine large quarries on the mountains surrounding the plain. The seashore is occupied by various industries that are performing a variety of activities (both legislated and unlegislated). The 13 port facilities of Elefsis Harbor handle 5500 ships per year and a total cargo load that is 2.5 times larger than that handled by Piraeus Harbor. All of these activities coexist within the same environment, with the road network (three freeways) and two rail networks crossing the urban poles in this region (Mavrakis et al. 2016, 2020; Cecchini et al. 2019).

During 15–20 May 2020, Greece experienced the earliest heatwave ever recorded in the country according to the criteria of Metaxas and Kallos (1980). The air temperature during the heatwave exceeded the highest value that had been recorded during the second ten-day period of May (35.8 °C in 11th May 2003) for over 150 years, according to the historical climatic archive of the National Observatory of Athens (Founda 2020). The Thriasio Plain is particularly vulnerable to temperature extremes. When combined with poor air quality due to the intense industrial activity in the area, these temperature extremes can cause unfavorable environmental conditions for human health (Theoharatos et al. 2010; Mavrakis et al. 2012).

In the present study, we examined the thermal and air quality conditions in the Thriasio Plain during the heatwave experienced in Greece during 15–20 May 2020. Two commonly used indices, the discomfort index (DI) and the daily air quality (stress) index (DAQI), were estimated. This early heatwave occurred shortly after some of the restrictions that were imposed to halt the spread of coronavirus disease 2019 (COVID-19) in Greece were lifted (on 4 May 2020), when the process of gradually restoring normal industrial and transportation activities in the area began. The results of this study can be used by policy makers to formulate appropriate resilience measures that support the health of the population during heatwaves (McElroy et al. 2020).

Materials and methods

Study area

Greece has a Mediterranean climate (Kottek et al. 2006). The warmest months are July and August (the mean maximum temperature ranges between 29 °C and 35 °C; HNMS 2020). Spring (March–May) is short (HNMS 2020). In Athens, the monthly mean maximum temperature in May has increased by about 0.2 °C per decade over the last 130 years (Founda 2020). The highest monthly mean maximum temperature recorded was 29 °C (in 1945), followed by 28.8 °C (in 2003). The lowest monthly mean maximum temperature was 20.2 °C (in 1919; Founda 2020). Based on data recorded at the Thissio Climatic Station of the National Observatory of Athens (NOA) from 1860 onwards, the highest hourly maximum temperature (37.6 °C) was recorded on 29 May 2008. Hourly maximum temperatures of over 35 °C have been recorded on a total of 16 days, all of which were in late May except for a temperature of 35.8 °C recorded on 11 May 2003. Hourly maximum temperatures of over 36 °C have been recorded on four days in late May (25 May 1945, 29 May 1950, 26 May 1994, 29 May 2008).

In mid (15–20) May 2020, Greece experienced a heatwave with extremely high temperatures (39–40 °C); indeed, it set many maximum temperature records (HNMS 2020; Vougioukas et al. 2020). In Athens, the highest air temperature (37 °C, Thissio Station) was recorded on 16 and 17 May 2020, which was the highest in the last 160 years for mid-May (Founda 2020). At the same time, a low-intensity but long-duration dust transport episode from Africa was recorded.

The Thriasio Plain (Fig. 1) extends over an area of 120 km2 in the prefecture of Western Attica, northwest of Greater Athens. There are three municipalities in the Thriasio Plain: Elefsis, Aspropyrgos, and Mandra, with a resident population of 78,038 according to the 2011 census (HSA 2020). The Thriasio Plain has a smooth surface morphology that is slightly inclined towards the sea (3%) and is surrounded by areas with high relief. Its close proximity to the sea and the local climatic conditions produce local regressive atmospheric circulation patterns that greatly inhibit atmospheric self-cleaning through dispersion and transport mechanisms. Air temperature inversion heights are especially low during the cold period of the year; they are frequently lower than the surrounding hills and comparable to the heights of the highest chimneys of the large industrial complexes. Thus, air pollutants can be trapped within a shallow layer, resulting in high (daily and hourly) pollution concentrations (Lykoudis et al. 2008; Mavrakis et al. 2008, 2015b; Toumpos et al. 2017).

Fig. 1

Map of Attica Prefecture

Data sources

Hourly values of the air temperature (Ta, °C), the relative humidity (RH, %), and the concentrations of sulfur dioxide (SO2), nitrogen dioxide (ΝΟ2), ozone (Ο3), and particulate matter less than 10 μm in diameter (PM10) for the time period between 1 March and 31 May 2020 were used in this study. These data were obtained from the meteorological stations in Elefsis Air Base (LGEL), Elefsis City (ELE), and Alonistra (AL). LGEL is operated by the Hellenic National Meteorological Service (HNMS), World Meteorological Organization (WMO) code LGEL (167180). ELE is located in the center of Elefsis City (1000 m from the seashore), away from large constructions. It is operated by the Bureau of Pollution for the Municipality of Elefsis, and provides air quality data (NO2, O3, PM10). AL is located in Aspropyrgos (500 m south of the city center and 1000 m from the seashore). It is operated by the Bureau of the Environment for the Municipality of Aspropyrgos (BEMA), and provides air quality (SO2, NO2, O3, PM10) and meteorological (Ta, RH) data.

Assessment of the thermal conditions and air quality

Thermal conditions were assessed using the discomfort index (DI, °C; Thom 1959), which was estimated based on an equation in Paliatsos and Nastos (1999), modified for Greece:where Ta is the dry bulb temperature (°C) and RH is the relative humidity (%). The DI categorizes thermal environments into levels of thermal discomfort according to the assessment scale presented in Table 1.

Table 1

Assessment scale for the discomfort index (DI; Giles et al. 1990) and the daily air quality index (DAQI; Katsoulis and Kassomenos 2004)

DI (°C)	Thermal conditions	DAQI	Air quality
18–21	No discomfort	< 0.2	Very low
21–24	< 50% of the population feels discomfort	0.2–0.4	Low
24–26	> 50% of the population feels discomfort	0.4–0.6	Moderate
27–29	Most of the population suffers discomfort	0.6–0.8	Distinct
29–32	Everyone feels severe stress	> 0.8	Strong
Over 32	State of medical emergency	Independent of DAQI	Extreme

Air quality was assessed based on the daily air quality index (DAQI). The DAQI is a measure of the air quality and the health risk from air pollution that is derived from the average daily concentration of air pollutants using the following formula (Matzarakis and Mayer 1991; Kassomenos et al. 2004; Katsoulis and Kassomenos 2004):where i is the air pollutant under consideration (SO2, NO2, O3, or PM10) and “MI-24 h value” is the threshold concentration of that air pollutant according to European Community standards (hourly SO2: 350 μg/m3; hourly NO2: 200 μg/m3; mean 8-hourly O3: 180 μg/m3; daily PM10: 50 μg/m3). Values were calculated for hourly data, which were then used to derive daily averages. The assessment scale for the DAQI is presented in Table 1.

Statistical analysis

Pearson’s correlation test was applied to examine the possible relationships between air temperature, relative humidity, air pollutant concentrations, DI, and DAQI. The confidence level was set at 0.05.

Results

The average monthly Ta varied between 13.1 °C and 21.8 °C during the period from 1 March to 31 May 2020 in the Thriassio Plain (Table 2). The estimated monthly Ta values from the AL station were slightly higher than those from the LGEL station. The hourly air temperature ranged between 3 °C and 38 °C.

Table 2

Descriptive statistics for parameters of interest

		March 2020				April 2020				May 2020				March–May 2020				Lockdown 2020				Heatwave May 2020
		Average	Stdev	Max	Min	Average	Stdev	Max	Min	Average	Stdev	Max	Min	Average	Stdev	Max	Min	Average	Stdev	Max	Min	Average	Stdev	Max	Min
LGEL	Ta (°C)	13.1	3.9	24.0	3.0	15.0	4.5	27.0	7.0	21.4	5.3	38.0	12.0	16.5	5.8	38.0	3.0	14.9	4.6	27.0	5.0	27.1	5.1	38.0	18.0
	RH (%)	67	17	100	26	58	17	100	20	50	16	94	14	58	18	100	14	60	18	100	20	43	15	83	14
	DI (°C)	13.1	3.2	20.6	4.2	14.6	3.3	22.2	7.3	19.1	3.3	27.1	12.1	15.6	4.1	27.1	4.2	14.5	3.4	22.2	6.3	22.8	2.6	27.1	17.2
AL	Ta (°C)	13.9	3.6	24.6	5.6	15.7	4.1	25.0	8.1	21.8	4.6	36.1	12.8	17.0	5.3	36.1	5.6	15.6	4.1	26.2	8.1	27.2	3.8	36.1	20.4
	RH (%)	64	19	100	25	55	19	100	18	48	15	100	14	56	19	100	14	59	20	100	18	43	13	88	15
	DI (°C)	13.9	2.8	20.5	7.3	15.2	2.9	21.1	8.6	19.5	2.9	26.3	12.8	16.1	3.7	26.3	7.3	15.1	3.0	21.7	8.4	23.0	2.0	26.3	18.8
	SO2 (μg/m3)	12	9	59	0	11	9	110	0	14	16	123	0	12	12	123	0	10	8	110	0	29	26	123	0
	NO2 (μg/m3)	44	23	145	7	37	24	144	5	49	29	160	8	43	25	160	5	36	23	144	5	72	31	160	22
	O3 (μg/m3)	41	26	130	0	62	28	125	1	60	35	220	1	54	31	220	0	58	29	125	1	60	52	220	1
	PM10 (μg/m3)	41	27	209	4	35	22	156	1	47	37	295	3	41	29	295	1	35	22	156	1	90	41	295	4
	DAQI	0.41	0.10	0.61	0.23	0.40	0.12	0.65	0.19	0.54	0.32	0.94	0.23	0.45	0.21	0.94	0.19	0.39	0.12	0.65	0.19	0.83	0.10	0.94	0.71
ELE	NO2 (μg/m3)					65	52	207	1	86	69	345	2	77	63	345	1	57	33	118	5	133	75	345	6
	O3 (μg/m3)	9	7	35	3	49	25	126	3	57	28	156	5	45	29	156	3	40	23	79	5	52	39	156	5
	PM10 (μg/m3)	27	19	85	2	26	16	107	1	43	35	493	7	31	25	493	1	25	13	55	4	72	44	493	24
	DAQI	0.33	0.17	0.63	0.12	0.53	0.22	1.11	0.20	0.82	0.29	1.39	0.40	0.59	0.30	1.39	0.12	0.48	0.24	1.11	0.12	1.14	0.19	1.39	0.91

Hourly data were calculated and then used to derive daily averages

AL Alonistra, ELE Elefsis City, LGEL Elefsis Air Base, Max maximum, Min minimum, Stdev standard deviation, T air temperature, RH relative humidity, SO sulfur dioxide, NO nitrogen dioxide, O ozone, PM particulate matter less than 10 μm in diameter

Rapid increases in the mean daily air temperature, from 16.7 °C to 28.7 °C at the LGEL station and from 18.3 °C to 28.6 °C at the AL station, were observed between 7 May and 17 May. The maximum daily Ta occurred on 16 and 17 May at both LGEL (38 °C on both days) and AL (35.1 °C and 36.1 °C, respectively); see Fig. 2a. On each day, Ta exceeded the threshold of 36 °C for four consecutive hours (from 12:00 to 16:00 LST) at LGEL. On 17 May, Ta was 38 °C for four consecutive hours at LGEL, while the daily range of Ta was 20 °C (minimum Ta = 18 °C). At AL, Ta reached its maximum value (36.1 °C) for 1 h at 15:00 LST.

Fig. 2

Variation in a the maximum daily air temperature and b the maximum discomfort index (DI) and daily air quality index (DAQI) values

The average daily DI ranged between 13.1 °C and 19.5 °C. The DI values indicated that conditions leading to thermal discomfort (DI ≥ 24 °C) were present between 13 and 20 May 2020 (Fig. 2b), mainly during the hours immediately after midday (12:00 to 16:00 LST). Conditions leading to discomfort (24.1 °C to 25.9 °C) were also noted during the early evening hours (17:00 to 20:00 LST) at AL. A DI value of > 26 °C occurred each day from 16 to 20 May at LGEL. The maximum DI was recorded at 15:00 LST on 17 May at both LGEL (27.1 °C) and AL (26.3 °C) (Fig. 2b).

SO2, NO2, O3, and PM10 reached their maximum concentrations during the heatwave (Fig. 3a and b). The maximum daily concentration of SO2 (123 μg/m3) was recorded on 19 May (Fig. 3a), and the maximum daily concentration of NO2 ranged between 99 μg/m3 (20 May) and 160 μg/m3 (16 May) at AL (Fig. 3a) and between 164 μg/m3 (20 May) and 345 μg/m3 (16 May) at ELE station. The maximum daily concentration of O3 varied between 113 μg/m3 (20 May) and 219 μg/m3 (18 May) at AL, and between 84 μg/m3 (20 May) and 156 μg/m3 (17 May) at ELE. The maximum daily concentration of PM10 ranged from 156 μg/m3 (15 May) to 295 μg/m3 (20 May) at AL (Fig. 3a), and from 93 μg/m3 (16 May) to 493 μg/m3 (17 May) at ELE (Fig. 3b). The mean PM10 concentration was 90 ± 41 μg/m3 and 72 ± 44 μg/m3 at AL and ELE, respectively (Table 2).

Fig. 3

Variation in the daily maximum concentrations of air pollutants at a Alonistra station (AL) and b Elefsis station (ELE)

The daily air quality index (DAQI) was estimated using the concentrations of SO2, NO2, O3, and PM10 recorded at ELE and AL (Table 2). The highest estimated average monthly DAQI occurred in May (0.82 at ELE). The estimated average DAQI was below 0.5 during the lockdown and 0.83 at AL and 1.14 at ELE during the heatwave (15–20 May). Figure 2b shows that poor air quality was present during the latter period, with DAQI values exceeding the threshold of 0.8 (strong health risk) for several days. The maximum daily DAQI value reached 1.39 at ELE and 0.94 at AL (18 May). The mean daily DAQI value exceeded the threshold of 0.8 at both ELE and AL during the heatwave.

Correlations between the meteorological and air pollution parameters were analyzed using data from the stations in the region of interest, based on Pearson’s linear correlation coefficients (Table 3). A significant positive correlation was found between Ta and DI; the estimated correlation coefficients ranged between 0.75 and 1.0. RH and DI were negatively correlated, with correlation coefficients ranging between − 0.49 and − 0.55. A significant relationship was identified between DI and DAQI; the estimated coefficients varied between 0.57 and 0.82. The PM10 concentration showed the strongest correlation with DAQI (correlation coefficients: 0.84–1.00).

Table 3

Pearson correlation coefficients for the associations between air temperature (Ta), relative humidity (RH), air pollutant concentrations (SO2, NO2, O3, PM10), discomfort index (DI), and daily air quality index (DAQI)

Station		LGEL			AL								ELE
		Ta	RH	DI	Ta	RH	DI	SO2	NO2	O3	PM10	DAQI	NO2	O3	PM10	DAQI
LGEL	Ta	1
	RH	− 0.58	1
	DI	1.00	− 0.57	1
AL	Ta	0.99	− 0.59	0.99	1
	RH	− 0.50	0.91	− 0.50	− 0.49	1
	DI	0.99	− 0.57	0.99	1.00	− 0.49	1
	SO2	0.45	− 0.18	0.42	0.45	− 0.23	0.43	1
	NO2	0.45	− 0.08	0.42	0.45	− 0.13	0.43	0.68	1
	O3	0.27	− 0.57	0.28	0.29	− 0.45	0.29	− 0.19	− 0.53	1
	PM10	0.56	− 0.10	0.54	0.56	− 0.14	0.56	0.69	0.72	− 0.17	1
	DAQI	0.59	− 0.15	0.57	0.60	− 0.18	0.59	0.70	0.72	− 0.12	1.00	1
ELE	NO2	0.63	− 0.22	0.60	0.61	− 0.18	0.59	0.66	0.86	− 0.43	0.74	0.75	1
	O3	0.44	− 0.62	0.46	0.45	− 0.58	0.47	0.02	− 0.11	0.70	− 0.03	0.01	− 0.68	1
	PM10	0.67	− 0.29	0.66	0.68	− 0.30	0.67	0.50	0.63	0.05	0.85	0.86	0.55	0.08	1
	DAQI	0.82	− 0.46	0.81	0.82	− 0.45	0.82	0.59	0.69	0.16	0.84	0.86	0.75	0.30	0.92	1

Statistically significant coefficients are shown in bold

AL Alonistra, ELE Elefsis City, LGEL Elefsis Air Base, T air temperature, RH relative humidity, SO sulfur dioxide, NO nitrogen dioxide, O ozone, PM particulate matter less than 10 μm in diameter

Discussion

The present study examined the thermal conditions and air quality during the spring of 2020 in an industrialized plain in Attica, Greece. This period was marked by the imposition of restrictions aimed at suppressing the spread of COVID-19 and the subsequent gradual easing of these restrictions in mid-May, which coincided with an early heatwave. The maximum air temperature during the heatwave was more than 10 °C above the monthly normal climatic average (25.8 °C, LGEL station; HNMS 2020). The results showed a high level of thermal discomfort and poor air quality. Significant correlations between the air temperature, DI, and DAQI were identified. DAQI was strongly correlated with PM10 concentration, indicating that the PM10 concentration is an important influence on air quality.

Heatwaves are associated with the large-scale circulation of air masses in the upper atmosphere, as this circulation favors the occurrence of such phenomena. Depending on the season, the subtropical jet stream plays a crucial role in heatwave generation at northern latitudes, such as the 2018 and 2019 heatwaves in Central and Northern Europe, because it transports hot air masses at these latitudes. The early heatwave in May seems to have occurred due to an unusually early shift of the jet stream to northern latitudes (above the northern border of Greece), leading to the circulation and transportation of very hot air masses from North Africa, as documented in previous studies (Flocas et al. 2009; Theoharatos et al. 2010; Mavrakis and Tsiros 2018). Figure 4 shows the results of back-trajectory analysis of the air masses using the FLEXTRA model. This model uses meteorological data provided by ECMWF (the European Centre for Medium-Range Weather Forecasts) (www.nilu.no/trajectories; Stohl 1998). The results show the flow of air masses from the Sahara Desert and their circulation above the central Mediterranean Sea and Greece during the aforementioned time period. The circulation of hot, dry air masses above the Mediterranean Sea humidified the masses (Theoharatos et al. 2010). At the same time, due to the low speed of the wind above the sea, the sea surface temperature increased, yielding a marine heatwave event, in agreement with Mavrakis and Tsiros (2018). Also, according to Table 3, the PM10 values were strongly correlated with the DAQI values. This is an indication that the high PM10 concentrations are due to a Sahara dust transport episode. The calm conditions in the area favored high PM10 concentrations and prevented the dispersion of particulates.

Fig. 4

Back trajectories of air masses for the period 15–20 May 2020 at 12 UTC. The trajectories were calculated using the FLEXTRA model and meteorological data were provided by the European Centre for Medium-Range Weather Forecasts (ECMWF)

The abrupt increases in the air temperature, DI, DAQI, and PM10 came soon after a lockdown period lasting almost two months (from 23 March to 4 May 2020) that was prompted by the COVID-19 pandemic (Fig. 2). This lockdown was implemented despite being strongly criticized (Ioannidis et al. 2020). During the lockdown, the values of both DI and DAQI were relatively low according to all of the criteria mentioned above. The early spring heatwave led to abrupt increases in DI and DAQI (i.e., poor air quality). The time evolutions of these two indices during this period took the form of a shock wave. Sudden large increases in DI and DAQI in a particular area can be lethal, as noted in previous studies, due to thermal shock in the general population and poor air quality (Pantavou et al. 2011; Theoharatos et al. 2010).

Conclusions

This study showed an unusual and rapid increase in air temperature in the Thriasio Plain (Greece) that was indicative of an early spring heatwave, with the air temperature rising to more than 10 °C above the monthly average value for this region. Conditions of thermal discomfort were found to be positively correlated with poor air quality, mainly due to high levels of PM10 that may have been transported from the Sahara Desert. High ambient temperatures combined with high humidity and low air velocity (calm conditions) led to high thermal discomfort and contributed to poor air quality for a relatively long period of time.

Such episodes are not expected to persist in early spring, so the occurrence of this heatwave event should encourage authorities to plan precautionary measures for future events.