1 Introduction

The North Atlantic Oscillation (NAO) is the normalised sea level pressure difference between the Azorean high and the Icelandic low, which is responsible for a considerable share of variability in European climate. Surface air temperature, the precipitation and evaporation balances in the Atlantic Ocean and the surrounding land masses are influenced asymmetrically both in time and space (Stenseth et al. 2002).

The impacts of the NAO include a meridional sea level pressure gradient in the North Atlantic, with positive phases characterised by an intense north–south gradient and negative phases when the gradient decreases (Hurrell 1995). The strength of each NAO phase is given by the NAO index, which can be a positive or negative number, according to the phase. Its fluctuations induce considerable variations in the mean direction and speed of the wind patterns over the Atlantic (Trigo et al. 2002; Hurrell et al. 2003), the heat and moisture transport to the neighbouring continents (Dayan 2011; Christoudias et al. 2012), and the intensity and amount of storms (Trigo 2006). Consequently, given NAO’s hemispheric character (Greatbatch 2000), other associated phenomena have been reported as suffering from its influence. For instance, the influxes of dust from the Sahara desert towards the Atlantic and Mediterranean (Moulin et al. 1997), and photochemical processes that have an effect on the transport of atmospheric contaminants from North America to Europe (Li et al. 2002). Ultimately, all these impacts can result in direct human health-related parameters, such as hospital admissions due to diseases linked to air pollution episodes (Almendra et al. 2017). Such effects are likely to influence the composition of the atmosphere, either via local processes or through a significant transport of air pollutants (Izquierdo et al. 2014; Camarero et al. 2017).

The atmospheric concentration of persistent organic pollutants (POPs), and of polycyclic aromatic hydrocarbons (PAHs) in particular, is closely related to climate change and climatic fluctuations, which are likely to influence contaminant transport pathways and transfer processes. Predicting how climate variability alters PAHs concentrations in the atmosphere still poses an exceptional challenge (Wang et al. 2010), since the existing literature on the climatological framework (long-term modelling) of the influence of NAO or other climate indexes is scarce (Eckhardt et al. 2003; Christoudias et al. 2012; Bacer et al. 2016), and practically non-existent with the aim of assessing local and large-scale approaches separately.

Regarding semi-volatile organic contaminants (SVOCs), the references to the influence of NAO on their levels and distribution in the atmosphere are scarce, and mainly focusing on polychlorinated biphenyls (PCBs) and organochlorine pesticides (OCPs). Ma et al. (2004a, b) studied the influence of NAO fluctuation on the distribution and deposition of some PCBs and OCPs in the Great Lakes region and the whole of North America, respectively, while Hung et al. (2005) did a similar work while assessing the temporal and spatial trends of the same pollutant families in the Arctic region. MacLeod et al. (2005) reported variations in the atmospheric levels of PCBs in the Northern Hemisphere due to the NAO, using a global-scale model (BETR-Global) and later Wang et al. (2010) investigated ice core concentration records of OCPs DDT and HCH to assess the effect of climate parameters (including NAO) on their presence in mountain glaciers. Finally, Octaviani et al. (2015) studied the meridional transport of persistent organic pollutants and their correlations with the NAO and Arctic Oscillation (AO) indicating that the net burden of PCB153 reaching the Arctic will increase under a future climate scenario. To our knowledge, there are no studies with reference to PAHs, a family of SVOCs of extreme importance and concern (ATSDR 1995).

PAHs are generated by the incomplete combustion or burning of organic carbon-containing matter, and can be found in cigarette smoke, traffic exhaust, natural fuels and as by-products of many industrial processes and natural phenomena like forest fires or volcanoes (Ravindra et al. 2008). They are ubiquitous in the environment in gas and particulate form and can have an effect on human health through several forms of exposure (inhalation, diet, etc.). Among them, benzo(a)pyrene (BaP), with a five-aromatic ring molecular structure (thus present in the atmosphere mainly as particulate matter) (Ratola and Jiménez-Guerrero 2015) is possibly the most studied congener of this group of compounds and is often considered as the PAHs “representative” in environmental research (Nisbet and LaGoy 1992; Ratola and Jiménez-Guerrero 2017).

PAHs have been present in our planet since the beginning and will remain as long as a fossil fuel exists, and a combustion is made. This means that they will always play a crucial role in the assessment, control and mitigation of air pollution (Dat and Chang 2017). Europe is particularly affected by PAHs, given the level of industrialization and use of fossil resources. In fact, the 2019 Air Quality Report of the European Environment Agency (EEA 2019) states that an estimated 17% of the EU-28 urban population was exposed to annual mean BaP concentrations above the target value of 1 ng \(\hbox {m}^{-3}\) for the total content of BaP in the PM10 fraction (EU 2004). Although being the lowest value since 2008, Guerreiro et al. (2016) stated from a comprehensive modeling study over Europe that BaP needs to be continuously monitored as a proxy for all PAHs and the corresponding health effects estimated. And obviously, the transboundary character of airborne pollutants is a challenge for national-based policies. For instance, Schreiberová et al. (2020) estimated 46% and 70% of foreign contribution in the presence of BaP for most of the Czech Republic. A pan-European action is advisable under this scenario, since the progress in air quality mitigation is uneven within Europe. Viana et al. (2020) reported the management strategies put forward by 10 European cities (Antwerp, Berlin, Dublin, Madrid, Malmö, Milan, Paris, Plovdiv, Prague, Vienna) from 2013 to 2018 and identified several knowledge gaps that pointed towards a systemic thinking approach to the future of air quality management.

Climate change also plays a very important part in this equation since the increased frequency and intensity of extreme events complicates the establishment of monitoring strategies and mitigation measures. As an example, forest fires have been growing in number and intensity in recent years, but not only in the southern part of the domain (where they were already frequent), also in the north of Europe as in the recent wildfire events in Scandinavian countries in 2018 (Turco et al. 2019). This means that a tight control of BaP (and the rest of PAHs) is needed to face their atmospheric presence adequately. Consequently, indexes like NAO can inform about real-time changes in BaP concentrations in Europe and, with the help of CTMs, project the near future and long-term incidence. These data will definitely help air quality managers and legislators to decrease the environmental fingerprint of PAHs

Consequently, the objective of the current study is to assess the influence of the NAO fingerprint on the mean concentration of BaP in a European domain and focusing on the effect of regional-scale processes. The target domain comprises the whole Mediterranean basin and its areas of influence, since this is considered one of the most vulnerable areas to air pollution and climate change worldwide (Linares et al. 2020). Moreover, in recent years the influence of climate variations in air pollution over the target area has become evident (e.g. Dayan et al. (2017)), especially the influence of NAO phase in central and southern Europe and the Mediterranean basin (Christoudias et al. 2012).

2 Materials and methods

Taking advantage of the formulations previously reported in literature on the impact of the NAO on European aerosol ground levels (e.g. (Jerez et al. 2013b)), this approach follows an analogous methodology, applied to the widespread contaminant BaP. The novelty of this work is a contribution to the understanding of the regional processes controlled by NAO that affect BaP levels in a climatological timescale (30 years), which to our knowledge is still not established. The analysis was focused on winter and summer seasons using composites that highlight the changes in the mean fields of modelled BaP concentrations between NAO’s high and low phases (NAO+ and NAO−, respectively) (Fig. 1). This figures indicates a slight trend of the NAO index observed for the summer but not any evident trend for winter, which is the main season of interest in NAO-related studies.

Fig. 1
figure 1

Winter (top) and summer (bottom) original CPC/NOAA NAO index series for the period 1971–2000. The former are December–January–February averages and the latter are June–July–August averages. The red (blue) lines denotes the 70th (30th) percentile value of each series. The years/values with a NAO index above (below) that percentile are selected as NAO+ (NAO−) events for the composites analysis

In brief, the NAO indexes available at the Climate Prediction Center (CPC) from the National Ocean and Atmospheric Administration (NOAA) for winter and summer months (December, January, February—DJF and June, July, August—JJA, respectively) are employed to identify positive and negative NAO stages throughout the period of study (1971–2000). The CPC/NOAA NAO index was selected in this case since it is readily available (Pey et al. 2013). Following a common approach (e.g. (Trigo et al. 2002; Jerez et al. 2013a, b)), we defined classes of NAO positive and negative phases using the detrended monthly series of the NAO index. In fact, detrending the original NAO series shown in Fig. 1 in the computational code ensures that both classes have a similar occurrence and neither NAO+ nor NAO− phase prevail at the beginning or the end of the period.

It should be noticed that different data sources and/or methods to compute the NAO index may provide different values. In particular, we are aware that the CPC/NOAA NAO index is not based on the ERA-20C data (Poli et al. 2016) used to drive WRF in our modelling system. However, at the timescale addressed here, this influence is negligible (Kumar et al. 2016; Comas-Bru and Hernández 2018). Hence, this choice is based on its easy access and up-to-date computation approach, which in fact makes it widely used nowadays in studies focused on the recent past (Bladé et al. 2012; Pey et al. 2013; Grams et al. 2017; Alizadeh-Choobari and Adibi 2019).

NAO+ and NAO− phases are defined according to the 70th and 30th percentiles of the NAO series, respectively (Jerez et al. 2013a). The influence of the NAO on BaP concentrations is then evaluated through composites estimating the differences in the mean concentrations of BaP simulated for the periods with negative and positive NAO phases. The significance checks performed on the difference between BaP levels during NAO+ and NAO− phases follow those of Jerez et al. (2013a), and can be summarised in two requirements: if they have statistical significance above the 90% confidence level and if they are backed by temporal correlations between the NAO and the corresponding BaP series that are statistically significant also above the 90% confidence level. Statistical significance is evaluated by performing two-tailed t-tests for the null hypothesis of equal means or zero correlation, respectively (Snedecor and Cochran 1989). Moreover, the signal-to-noise (S/N) ratio (defined as the ratio between the absolute value of the NAO+ minus NAO− differences in BaP concentrations and the standard deviation of the whole BaP seasonal series) is considered a relative measure of the importance of the NAO-impact. In fact, areas where the S/N ratio is below the unit were blurred. In the other areas, the NAO-impact exceeds one standard deviation of the series and hence we are likely moving from one tail to the other of the probability density function describing the variability in the BaP levels when the NAO phase changes.

2.1 Models and domains

BaP levels are obtained from a simulation of the chemical composition of the atmosphere from 1971 to 2000, in a domain comprising most of Europe (Mediterranean basin and central Europe). The modelling system is WRF+CHIMERE, coupled with the European Monitoring and Evaluation Program (EMEP) emissions, and run with a horizontal spatial resolution of 25 km. Thirty-three vertical levels are spaced unevenly up to 10 hPa, with higher resolution near the ground in order to characterize the particularities (12 model layers within the lowest \(\sim\) 2000 m and 15 layers up to \(\sim\) 5500 m). This time period and resolution are significant compared to those commonly applied in climate runs for air pollution purposes. Further description can be found in the Supplemental Material and in Table SM1.

The simulation was designed to disregard the signals from the NAO impact on the long-range atmospheric transport (LRAT) by using constant climatological boundary conditions for the BaP concentrations. Therefore, the contribution of LRAT to changes in BaP during different NAO phases can be considered negligible, and the influence of climatological boundary conditions limited and overwhelmed by local processes.

Although the LRAT of pollutants is recognized and has been studied for some decades, its mechanisms are intricate and, for chemicals like PAHs, are far from being fully understood (Lammel et al. 2015). For BaP, the degradation and fate processes are the most uncertain (Shrivastava et al. 2017). Being almost completely associated with particulate matter, the distribution and LRAT of BaP are mainly ruled by degradation processes, particularly with ozone (Keyte et al. 2013). But also here the assessment of reaction rates and half-lives is difficult, as they are dependent on other parameters like temperature (Mu et al. 2018). Nevertheless, laboratory measurements indicate that ozone degrades particulate BaP rapidly (Zhou et al. 2019). Moreover, Lammel et al. (2009) indicated a modelled value of 48 h for BaP atmospheric lifetime, but some correction in the parameterizations in the meantime pointed towards a value of 6 h (Friedman and Selin 2012). For all these reasons, it is extremely difficult to quantify the percentage of the PAH burden attributable to LRAT but, unless in remote areas, this contribution tends to be very small compared with the local emissions (Callén et al. 2008; Lhotka et al. 2019), especially in highly populated and industrialized areas, that generally characterize our domain. Therefore, for the purposes of this study LRAT can be considered a negligible parameter.

Constant anthropogenic emissions are considered throughout the simulation period to isolate the possible impact of the NAO on the ground levels of atmospheric contamination. Further information on emissions can be found in the Suplemental Material (Fig. SM1). Moreover, the influence of climatology on airborne hazardous chemicals is estimated neglecting potential changes on anthropogenic pollutant emissions, land use, vegetation, or any feedbacks from those substances to meteorology parameters, although allowing variations regarding natural emissions (e.g. Meleux et al. (2007); Jiménez-Guerrero et al. (2012)). The NAO effects on LRAT of the target contaminant were not taken into consideration. Such settings enhance our knowledge on the role of the NAO-ruled underlying mechanisms on a regional level.

The robustness of this simulation was evaluated in other works, which rely on the same data and include the assessment of the integrations (Ratola and Jiménez-Guerrero 2015, 2016), and hence, only a short description is shown here. In brief, the validation of the modelling system used was based on the climatologies of BaP available in the target period, obtained from daily means coming from 13 stations belonging to the EMEP air sampling network (set at least 10 km off large emission sources, this data is assumed to represent regional background levels) (Tørseth et al. 2012). Ratola and Jiménez-Guerrero (2015) found that the modelling system provided consistently mean fractional biases (MFB) and mean fractional errors (MFE) in the range of 20–30% and 30–60%, respectively, when evaluated against EMEP stations in Europe. Moreover, the same authors reported in a subsequent study that the MFB of this modelling system varies from − 19% for spring (March–April–May) to a minor overestimation during winter (December–January–February, +9%), with biases showing a slight annual underestimation around − 3 picograms per cubic meter (pg \(\hbox {m}^{-3}\)) (Ratola and Jiménez-Guerrero 2016). These results support the use of the modelling system as a good depiction of BaP air concentrations and as a valid tool to calculate their changes under different NAO phases.

3 Results

3.1 NAO impact on atmospheric conditions

This Section does not intend to provide novel insights, but essentially to evaluate the ability of our WRF simulations (driving the CHIMERE run) to simulate the extensively reported NAO impact on the European atmospheric conditions. As this analysis was widely conducted by Jerez et al. (2013a), we mainly summarize this work to give the reader an appropriate context for the interpretation of the NAO impact on BaP concentrations, since the analysis is specifically focused on those variables driving air pollution levels. Fig. 2 further confirms the expected NAO impact on the wind field, showing asymmetric responses in winter and summer.

Fig. 2
figure 2

NAO impact in winter (left) and summer (right) on atmospheric conditions. (First to bottom rows) 10-m wind direction (WD10) during NAO+ (red) and NAO− (blue) phases (arrow proportional to the wind speed); NAO+ minus NAO− for the mean 10 m-wind module (WM10) (m \(\hbox {s}^{-1}\)), precipitation (PREC) (mm \(\hbox {month}^{-1}\)), column integrated cloud water (CICW) (mm) and 2-m temperature (T2M) (K). Results are only represented only if significant (p<0.1). Dots blur the values not representing a S/N ratio >1, and contours depict differences expressed as percentage (adapted from Jerez et al. (2013b))

In winter, NAO+ (NAO−) promotes a westerly flow in northern (southern) Europe, with the largest significant differences in the wind speed (above 10%) appearing northward in the western sector of the domain. Conversely, in summer, NAO− phases are associated with a stronger westerly flow in northern Europe than in NAO+ phases, with differences in the wind module exceeding 5% and extending more eastward than in winter. Windier conditions associated to the NAO− phase prevail also in the southwest of the domain during summer, mainly around the Iberian Peninsula (differences about 5%).

Consequences in precipitation, cloudiness and temperature derived from the former conditions are also provided in Fig. 2, as these variables play a key role on the BaP levels. NAO+ minus NAO− differences in precipitation are larger in winter than in summer, when they are overall negative and affect mainly southern Europe ranging west to east from 40% (50 mm \(\hbox {month}^{-1}\)) to 20% (20 mm mm \(\hbox {month}^{-1}\)). In summer, these differences reach 20% (up to 30 mm mm \(\hbox {month}^{-1}\)) over large sectors of Europe, being negative in the north and positive (and less important) in the south. Consistently, the column integrated cloud water (a variable representative of cloudiness) shows negative values in southern Europe during winter (resembling the west to east gradient of the precipitation signal, with differences ranging from 30 to 10%) and in northern areas during summer (differences around 10%), and positive values in north-western areas during winter (i.e., NAO+ enhances cloudiness there by about 20%). Regarding the NAO impact on mean 2 m temperature, the most important signal consists of positive NAO+ minus NAO− differences (up to 3 K) largely spread over northern Europe in winter. Smaller negative (positive) differences (up just to 1 K) appear also over some southernmost (northern) areas in winter (summer).

These results support the capacity of the climate simulation to reproduce the expected responses to the NAO phase, providing a meaningful representation of how the atmospheric conditions vary between NAO+ and NAO− events, that will be used to understand the NAO impact on BaP levels in the next Section.

3.2 NAO impact on BaP levels

The difference between Icelandic Low and the Azores High pressures is intensified throughout the NAO+ phases, originating a stronger separation of the flow over Europe. Under these conditions, a higher northward transport over northern latitudes and an enhanced southbound transport over the south are seen. Also, in winter, NAO+ promotes a windy west flow in north Europe whereas NAO− has a similar effect in the south. NAO+ blocks the zonal circulation in the Mediterranean basin, reducing precipitation and cloudiness there (e.g. (Trigo 2006; Jerez et al. 2013a), among others). High NAO index winters (NAO+) lead to drier weather over the Mediterranean, the majority of central and southern Europe and some areas of the east, while a precipitation rate above average falls over northern Europe (Hurrell 1995). In consequence, higher ground-level concentrations of BaP (almost exclusively in particulate form in air (Bidleman 1988)) are more frequently modelled during NAO+ than during NAO− (Figs. 3 and 4), since the cleaning effect of both enhanced winds and wet deposition and the inhibitory effect of diminished solar radiation are reduced. In this sense, the 10-m wind speed in winter is enhanced in the NAO+ phase by over 10% when compared to NAO− phase in the northern part of the domain.

Fig. 3
figure 3

Absolute ground-level concentrations (ng \(\hbox {m}^{-3}\)) of BaP averaged during periods with NAO+ and NAO− phases for winter (left) and summer (right) months; 1971–2000

The regions holding the highest modelled increases in BaP concentrations during NAO+ and NAO− phases (Fig. 4) are predominantly located in the Mediterranean basin during winter (changes around 0.05 to 0.10 ng \(\hbox {m}^{-3}\) between NAO+ and NAO− phases). For instance, the absolute concentrations in Lisbon and in the Po Valley in northern Italy exceed 0.70 ng \(\hbox {m}^{-3}\) in winter (Fig. 3) during NAO+ phases, being 0.20 ng \(\hbox {m}^{-3}\) higher than during NAO− in January (Fig. 4). This increase in concentration in BaP during NAO+ phase is particularly evident over the Iberian Peninsula, northern Italy and Greece, where the NAO+ minus NAO− signal is clearly negative for 10-m wind speed, precipitation and cloud water content. Thus, NAO+ raises over 100% the concentrations of BaP in southern Europe (Fig. 5) (especially in December, where the changes become more evident) as a consequence of the drier weather and increased stagnation associated to the NAO+ phase over the Mediterranean area.

Fig. 4
figure 4

Absolute ground-level differences (ng \(\hbox {m}^{-3}\)) in BaP concentrations between NAO+ and NAO− phases for winter (left) and summer (right) months; 1971–2000. Results are only shown where statistically significant (p < 0.1)

Fig. 5
figure 5

Differences (%) in ground-level BaP concentrations between NAO+ and NAO−phases for winter (left) and summer (right) months; 1971–2000. Results are only shown where statistically significant (p < 0.1)

Conversely, in the north of France and Germany, BaP absolute background concentrations are about 0.15 ng \(\hbox {m}^{-3}\) during NAO+ positive phase (Fig. 3), and reductions up to − 0.10 ng \(\hbox {m}^{-3}\) (Fig. 4) are associated to NAO+ phases vs. NAO− phases (particularly in January); that is, a − 80% change (Fig. 5). On the other hand, during conditions associated with NAO− phase, the air masses are more stagnant on the centre and south of Europe, allowing the accumulation of BaP, as shown in Figs. 4 and 5 (left column). In fact, the concentrations over the Po Valley (northern Italy), Switzerland and Austria may change over 0.20 ng \(\hbox {m}^{-3}\) (changes around 100%) in NAO+ vs. NAO− phases.

In summer, NAO− phases usually display more intense west flows in the north of Europe than NAO+ phases, with the wind module (10-m wind speed NAO+ minus NAO−) showing differences above 0.4 m \(\hbox {s}^{-1}\) (corresponding to a 5% change) (not shown) and extending more eastwards than in winter. The scenario of stronger winds matching the NAO− phase also occurs more frequently in the southwest of the domain in the warmer months, particularly around the Iberian Peninsula (10-m wind speed NAO+ minus NAO− resulting in differences of − 0.3 m s\(^-1\)). During summer, BaP changes are similar for the entire period (increases of 0.10–0.15 ng \(\hbox {m}^{-3}\) during NAO+ phases vs. NAO− phases, +60%, since the NAO+ concentration is about 0.25 ng \(\hbox {m}^{-3}\) and during NAO− is about 0.15 ng \(\hbox {m}^{-3}\)) (Fig. 3), but displaced towards central and northern Europe as the summer advances (Belgium, the Netherlands and Poland in June–July; northeast France and the British Isles, especially in August), since the role of NAO+ in reducing wind, precipitation and cloudiness is more noticeable there.

In line with the winter findings, the strong effect of the NAO in this area and season is indicated by the wide ranges in the BaP levels between both NAO phases, which can rise up to 100% (Fig. 5) especially in Belgium, the Netherlands and the British Isles. In the latter, however, results have to be taken carefully: BaP concentrations are below 0.05 ng \(\hbox {m}^{-3}\) during both NAO+ and NAO− phases (Fig. 3). In north-central Europe, the corresponding increments associated with NAO+ phase can reach or in some cases even surpass 0.20 ng \(\hbox {m}^{-3}\) in absolute value during August (Fig. 4).

Although wind is a key parameter for atmospheric presence of BaP, other factors can also be decisive on the incidence of this kind of mostly-particulate airborne pollutants, like precipitation and temperature. Rainfall promotes the wet deposition of particulate-bound compounds including PAHs (Birgü et al. 2011), and consequently, the decrease of their air concentrations. The NAO+ minus NAO− variances in precipitation are higher in winter and negative overall (Jerez et al. 2013a), affecting mainly the south of the domain and ranging from 40 to 20% from west to east, respectively.

In summer, the differences only reach 20% over most of Europe, being negative in the North and positive (and weaker) in the South. Observing Fig. 5, the matches are not evident with the differences in BaP concentrations. Although stronger in winter as well, they are not always negative (particularly in southern Europe) and in the summer months the negative trends in the North are also not seen.

This suggests that the influence of dry deposition may not be the most decisive for BaP atmospheric levels. Instead, temperature has a strong influence on the emission patterns, transport behaviour, and environmental distribution of SVOCs. For example, high temperatures promote the revolatilisation of chemicals like PAHs from soil and the ocean into the atmosphere and enhance the gas-phase distribution, contributing to the decrease of the respective half-lives (Delgado-Saborit et al. 2013). When applying the GEOS-Chem model to atmospheric PAHs, Friedman and Selin (Friedman and Selin 2012) identified the inclusion of secondary sources like revolatilization from soils (and to a much lesser extent from water bodies) as a priority for future modelling strategies. To address this issue, Zhang et al. (2017) mentioned that multi-year simulations and multi-compartment modules are mandatory, along with a dry-deposition parameterization (Zhang et al. 2015). In any case, the authors refer that this phenomenon is prone to affect essentially the more volatile 3 and 4-aromatic ring PAHs and not so much the modelled concentrations of 5-ringed BaP.

In turn, lower temperatures prompt the population to increase the use of PAH-emitting sources like traffic or building heating and for this reason are commonly linked with a higher atmospheric presence of PAHs (Morville et al. 2011). Nonetheless, this fact is not included in this contribution since road emissions in this study are not sensitive to temperature and the influence of temperature on BaP levels is established through changes in transport patterns and photochemistry instead, and indirectly through emissions temporal profiles.

The impact of the NAO on the mean 2-m air temperature is also reported by Jerez et al. (2013a), and the authors refer that the most significant signal is the positive NAO+ minus NAO− differences in winter (reaching 3 K) over the north of Europe, together with lower negative (positive) variations (only up to 1 K) found in some southern (northern) territories in the colder (warmer) months. In this case, the BaP concentration differences (Fig. 4) seem to be inversely correlated with those of the temperature (especially in winter), which is in line with the aforementioned tendency to have higher (lower) PAHs incidence with lower (higher) temperatures (Garrido et al. 2014).

As for other types of airborne particulate pollutants (Jerez et al. 2013a), the variation of several climatic indicators such as temperature (mainly in the north), the distribution of rain or the wind patterns associated with shifts in the NAO parameter are also closely related to the influence of NAO on BaP air levels. A summary of the variables involved can be found in Fig. 6.

Fig. 6
figure 6

Summary of the processes involved in the NAO impacts on the concentration of BaP during (top) NAO+ phase and (bottom) NAO− phase (base image adapted from Wanner et al. (2001))

4 Conclusions

The effect of climate phenomena like the North Atlantic Oscillation (NAO) is evident not only on meteorological parameters like temperature, wind or precipitation, but also on the distribution and transport of macropollutants (substances that comprise 90–95% of the mass fraction of environmental pollutant release) and micropollutants (substances with a small share, both with an organic—PAHs, pesticides, organohalogen compounds, etc.—or inorganic—metalloids, heavy metals and radionuclides—origin) (Konstantinova et al. 2020). This entails impacts from the local to the global scale which are challenging to assess.

This contribution sheds some light on the process-understanding related to the atmospheric concentration of PAHs, one of the most studied organic micropollutant families. In fact, when the concentrations of BaP (the most commonly used PAH congener) during NAO+ minus NAO− phases are estimated over Europe, strongly marked effects can be perceived temporally and spatially, namely yielding differences between seasons and north and south, respectively. While the Mediterranean basin displays a greater incidence of BaP in colder months during the NAO+ phase, summer brings the same signals up north. As a consequence, the mean BaP concentrations can vary 100% or in some cases even more, being of the same order of magnitude as the variations found for different daily meteorological conditions. Lighter evidences of opposite trends can also be detected, such as the NAO− phase inducing an increase in BaP levels are observed in northern locations in winter and southern ones in the warmer months. Even taking into consideration that our approach focuses on local-governing processes of air pollution, these results agree with those of Christoudias et al. (2012) in the sense that a stronger Azores High in the NAO+ phase increases the north-to-south transport in Europe and therefore higher BaP incidence.

Apart from enhancing the state-of-the-art on key correlations between climate and air quality, the current study sets the foundations towards (1) improving the predictability of BaP pollution episodes (a considerable body of work is under way to improve the accuracy in the prediction of large-scale circulation modes (e.g. Brands et al. (2012)), and (2) understanding NAO future projections on the long-term (still uncertain) regarding air quality. In addition, the present contribution highlights that meteorological conditions under different phases of the NAO can increase the ground levels of BaP independently of the emissions. Therefore, from a regulatory point of view, it should be born in mind that efforts devoted to reduce PAHs emissions in order to improve air quality can be hampered by the climate penalty associated to the enhanced levels caused by the natural climatic variability.

New chemicals are being produced and released into the environment and their hazardous potential and worldwide distribution need to be assessed as fast as possible. Modelling approaches need detailed information on the emission patterns and substantial field-based data prior to their development. In this sense, field sampling becomes an essential tool to determine the occurrence of hazardous compounds in the environment. Active and especially passive air sampling can provide information on the levels of many pollutants at different locations including remote areas, with a relative low-cost and easy maintenance setup (Domínguez-Morueco et al. 2019). Nevertheless, to help the transformation of these valuable datasets onto a complete understanding of spatial, temporal and chemical transport patterns it is essential to combine the field observations with modelling tools (Ratola and Jiménez-Guerrero 2016). If even for BaP, which has been studied for decades, we have important uncertainties, there is a need to reduce the time it takes to have reliable simulations for both legacy and emerging threats and the influence phenomena like NAO have on their presence and fate.