key: cord-0729578-lpu0w3or
authors: Mavrakis, Anastasios; Kapsali, Athanasia; Tsiros, Ioannis X.; Pantavou, Katerina
title: Air quality and meteorological patterns of an early spring heatwave event in an industrialized area of Attica, Greece
date: 2021-01-25
journal: EuroMediterr J Environ Integr
DOI: 10.1007/s41207-020-00237-0
sha: 9a92095fd2d58042fe4a390c6636435a182ca45c
doc_id: 729578
cord_uid: lpu0w3or

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.

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 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 (Mavrakis et al. , 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 (Mavrakis et al. , 2015b (Mavrakis et al. , 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 ). 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 (Mavrakis et al. , 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 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 ).

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 km 2 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 Mavrakis et al. 2008 Mavrakis et al. , 2015b Toumpos et al. 2017 ).

Hourly values of the air temperature (T a , °C), the relative humidity (RH, %), and the concentrations of sulfur dioxide (SO 2 ), nitrogen dioxide (ΝΟ 2 ), ozone (Ο 3 ), and particulate matter less than 10 μm in diameter (PM 10 ) 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 (NO 2 , O 3 , PM 10 ). 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 (SO 2 , NO 2 , O 3 , PM 10 ) and meteorological (T a , RH) data.

Thermal conditions were assessed using the discomfort index (DI, °C; Thom 1959) , which was estimated based on Paliatsos and Nastos (1999) , modified for Greece:

where T a 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 .

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 (SO 2 , NO 2 , O 3 , or PM 10 ) and "MI-24 h value" is the threshold concentration of that air pollutant according to European Community standards (hourly SO 2 : 350 μg/m 3 ; hourly NO 2 : 200 μg/ m 3 ; mean 8-hourly O 3 : 180 μg/m 3 ; daily PM 10 : 50 μg/m 3 ). 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 .

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.

The average monthly T a 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 T a 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.

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 T a 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, T a exceeded the threshold of 36 °C for four consecutive hours (from 12:00 to 16:00 LST) at LGEL. On 17 May, T a was 38 °C for four consecutive hours at LGEL, while the daily

average daily value MI-24 h value , range of T a was 20 °C (minimum T a = 18 °C). At AL, T a reached its maximum value (36.1 °C) for 1 h at 15:00 LST. 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) .

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

The daily air quality index (DAQI) was estimated using the concentrations of SO 2 , NO 2 , O 3 , and PM 10 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 T a 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 PM 10 concentration showed the strongest correlation with DAQI (correlation coefficients: 0.84-1.00). LGEL T a (°C) 13. 

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 PM 10 concentration, indicating that the PM 10 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 Fig. 2 Variation in a the maximum daily air temperature and b the maximum discomfort index (DI) and daily air quality index (DAQI) values 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/traje ctori es; 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 PM 10 values were strongly correlated with the DAQI values. This is an indication that the high PM 10 concentrations are due to a Sahara dust transport episode. The calm conditions in the area favored high PM 10 concentrations and prevented the dispersion of particulates.

The abrupt increases in the air temperature, DI, DAQI, and PM 10 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., Fig. 3 Variation in the daily maximum concentrations of air pollutants at a Alonistra station (AL) and b Elefsis station (ELE) 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) . Table 3 Pearson correlation coefficients for the associations between air temperature (T a ), relative humidity (RH), air pollutant concentrations (SO 2 , NO 2 , O 3 , PM 10 ), discomfort index (DI), and daily air quality index (DAQI) Statistically significant coefficients are shown in bold AL Alonistra, ELE Elefsis City, LGEL Elefsis Air Base, T a air temperature, RH relative humidity, SO 2 sulfur dioxide, NO 2 nitrogen dioxide, O 3 ozone, PM 10 particulate matter less than 10 μm in diameter

LGEL AL ELE 

Impact of cold waves and heat waves on the energy production sector

Urban sprawl and the "olive" landscape: sustainable land management for "crisis

Synoptic and local scale atmospheric circulation associated with air pollution episodes in an urban Mediterranean area

Early heatwaves over Italy and their impacts on durum wheat yields

The temperatures in May in Athens from the 19th century until today and the early heat of 2020. Kosmos Web Mag (Natl Observ Athens) 21

Too hot for comfort: the heatwaves in Greece in 1987 and 1988

Hellenic National Meteorological Service (HNMS) (2020) Climate of Greece

Hellenic Statistical Authority (HSA) (2020) Demographic characteristics 2011

Population-level COVID-19 mortality risk for non-elderly individuals overall and for non-elderly individuals without underlying diseases in pandemic epicenters

Assessment of the air-quality over urban areas by means of biometeorological indices: the case of Athens

Rubel F (2006) World Map of the Köppen-Geiger climate classification updated

Predicting photochemical pollution in an industrial area

The extreme heatwave in Athens July 1987 from the point of view of human biometeorology

Air quality levels in a closed industrialized basin (Thriasion Plain, Greece)

Towards (un)sustainable urban growth? Industrial development, land-use, soil depletion and climate aridity in a Greek agro-forest area

Population dynamics, industrial development and the decline of the traditional agro-forest landscape on the Mediterranean urban fringe: a case study in Greece

Analyzing the behaviour of selected risk indexes during the 2007 Greek forest fires

Mixing ratio as indicator of climate variations at a local scale: trends in an industrial area of the Eastern Mediterranean

Bioregional planning and design, volume II: issues and practices for a bioregional regeneration, part II: regional contexts, practices and projects for a bioregional recovery

Biometeorological and air quality assessment in an industrialized area of eastern Mediterranean-Thriassion Plain

The abrupt increase in the Aegean Sea surface temperature during the

Defining heatwaves using sub-regional meteorological data to maximize benefits of early warning systems to population health

Heatwaves from a synoptic point of view

The impact of global environmental changes on infectious disease emergence with a focus on risks for Brazil

Relation between air pollution episodes and discomfort index

Evaluating thermal comfort conditions and health responses during an extremely hot summer in Athens

Increasing trends in regional heatwaves

The growing threat of heat disasters

Computation, accuracy and applications of trajectories-a review and bibliography

synoptic conditions, bioclimatological assessment, air quality levels and health effects

(eds) Perspectives on atmospheric sciences

Impressive temperature records from the heatwave of

Assessment of the economic impacts of heat waves: a case study of Nanjing

Impact of heatwave on mortality under different heatwave definitions: a systematic review and meta-analysis

The authors would like to thank Dr. Anastasios Christides (Bureau of Environment, Municipality of Elefsis) and Eleni Verouti MSc (Bureau of Environment, Municipality of Aspropyrgos) for providing meteorological and air pollutant data.Funding There were no funding sources for this research.

Conflict of interest The authors declare no conflict of interest.

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 PM 10 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.