Observation and modeling of high-7Be concentration events at the surface in northern Europe associated with the instability of the Arctic polar vortex in early 2003

Events of very high concentrations of 7Be cosmogenic radionuclide have been recorded at low-elevation surface stations in the subpolar regions of Europe during the cold season. With an aim to investigate the mechanisms responsible for those peak 7Be events, and in particular to verify if they are associated with the fast descent of stratospheric air masses occurring during sudden stratospheric warming (SSW) events, we analyze 7Be observations at six sampling sites in Fennoscandia during January–March 2003 when very high 7Be concentrations were observed and the Arctic vortex was relatively unstable as a consequence of several SSW events. We use the GEOS-Chem chemistry and transport model driven by the MERRA-2 meteorological reanalysis to simulate tropospheric 7Be over northern Europe. We show that the model reasonably reproduces the temporal evolution of surface 7Be concentrations observed at the six sampling sites. Our analysis of model simulations, surface 7Be observations, atmospheric soundings of ozone and temperature and surface ozone measurements indicates that the 7Be peak observed in late February 2003 (between 20 and 28 February 2003) at the six sampling sites in Fennoscandia was associated with downward transport of stratospheric vortex air that originated during an SSW that occurred a few days earlier (between 18 and 21 February 2003).

The spatial and temporal variability of the 7 Be surface concentrations in Europe and their relationship with meteorological variables was previously analyzed in many studies (e.g., Piñero García et al., 2012;Błażej and Mietelski, 2014). The impact of the 11-year solar modulation on the 7 Be concentrations in the air is well established (e.g., Leppänen et al., 2010). The distinctive spring-summer maximum of 7 Be concentrations is widely described and mainly linked with the increased downward transport from the upper troposphere resulting from the intense convection and higher tropopause height typical of the warm season (Cristofanelli et al., 2006;Gerasopoulos et al., 2001Gerasopoulos et al., , 2003. In addition, cases of high 7 Be surface concentration, some of which occurred over the autumn-winter season, have been analyzed in Europe, e.g., over the Iberian Peninsula  and at high-altitude stations in the Alps and the Apennines (Brattich et al., 2017a, b;Cristofanelli et al., 2006Cristofanelli et al., , 2009. The spring-summer maximum was originally observed with fission products injected into the stratosphere during atmospheric nuclear tests (Dutkiewicz and Husain, 1985;Cristofanelli et al., 2018). A recent study by Salminen-Paatero et al. (2019), who used potential vorticity analysis to gain insights into stratosphere-to-troposphere transport of radionuclides at Rovaniemi (Finnish Lapland) indicated that the transfer of stratospheric air into the upper troposphere was at its maximum in March followed by descent to the ground level during late spring and early summer.
Further, over the last decade many studies have investigated the 7 Be records in northern Europe (Leppänen et al., 2010(Leppänen et al., , 2012Leppänen and Paatero, 2013;Sarvan et al., 2017;Leppänen, 2019), one of the three regions in Europe identified with a distinct 7 Be behavior in the surface air Hernández-Ceballos et al., 2015, 2016. Among these studies, Ajtić et al. (2016) analyzed the 7 Be concentration measured in Helsinki, Finland, over 25 years  and pointed out a relatively high number of 7 Be extremes occurring over autumn and winter: more specifically, 10 % of the highest 7 Be concentrations (above the 90th percentile) were observed in the cold season (October-March). Furthermore, recent studies have also indicated that the polar vortices can have a notable influence on the wintertime 7 Be surface concentrations in both the Northern (Ajtić et al., 2018;Bianchi et al., 2019;Terzi and Kalinowski, 2017) and the Southern Hemisphere (Pacini et al., 2015).
In particular, Ajtić et al. (2018) and Bianchi et al. (2019) employed two different methodologies to identify episodes of extremely high 7 Be surface concentrations in autumn and winter, pointing out a large number of cases over the October-March period of the years investigated. The comparison of the dates identified in both analyses showed an overlap with the events of the so-called sudden stratospheric warming (SSW) of the Arctic vortex, i.e., a sudden rise in the polar temperatures that leads to a highly irregular shape of the vortex and its misalignment from the pole. Ajtić et al. (2018) also noted cases of extremely high 7 Be concentrations occurring right after a very low 7 Be concentration over the Fennoscandian Peninsula during autumn and winter. Overall, this relationship between the SSW of the Arctic vortex and high 7 Be surface concentrations is likely linked to the perturbed stratosphere-troposphere interactions associated with SSW events, which could favor a fast descent of (1) midlatitude air rich in 7 Be, thus increasing this radionuclide's surface abundance, and (2) aged vortex air wherein 7 Be is subjected to radioactive decay and not transported from outside the vortex, thus decreasing the 7 Be surface abundance.
The atmospheric circulation in the Arctic is dominated by the presence of two distinct polar vortices, one in the troposphere and one in the stratosphere (Waugh et al., 2017). The two vortices present well-defined features: first, the vortex in the troposphere is much more extended than the stratospheric one, and second, while the tropospheric vortex is present all-year round, the stratospheric polar vortex exists only from fall to spring (Waugh et al., 2017). The stratospheric polar vortex in winter stems from the large-scale temperature gradients between the midlatitudes and the poles. Therefore, the stratospheric polar vortex begins to form in autumn as a result of the decreasing solar heating in the polar regions; it strengthens during winter and then breaks down in spring when solar radiation returns to the polar region. Larger topographic and land-sea contrasts and the resulting stronger upward-propagating waves in the Northern Hemisphere make the northern stratospheric vortex, or the Arctic vortex, weaker and more distorted than its Southern Hemisphere counterpart, the Antarctic vortex. SSW, a major mode of the temporal variability of the Arctic vortex, is the strongest manifestation of the coupling of the stratospheretroposphere system, with influence on the tropospheric flow lasting for many weeks  and with significant effects on chemical composition in the middle atmosphere (Sofieva et al., 2012;Tao et al., 2015). While major SSW, the so-called vortex split , can even cause the stratospheric vortex to break down during midwinter (Waugh et al., 2017), vortex displacements are instead characterized by a shift of the polar vortex off the pole and its subsequent distortion into a "comma shape" during the extrusion of a vortex filament . Previous works suggested that the occurrence of SSW events is capable of greatly perturbing the polar vortex and hence the stratospheric potential vorticity (PV) distribution (e.g., Matthewman et al., 2009) and the vertical distribution of ozone (e.g., Sonneman et al., 2006;Madhu, 2016). Additionally, the meteorological conditions associated with SSW events in the Arctic have been linked with the occurrences of 7 Be winter extremes, especially in the presence of a very high Scandinavian teleconnection index (Ajtic et al., 2018).
While the initial scientific interest over the stratospheric polar vortex was especially linked to the stratospheric ozone loss over the poles, it is now recognized that the vortices might affect the processes in the troposphere and surface weather (e.g., Mitchell et al., 2013). The present work aims to investigate in detail the atmospheric processes responsible for high 7 Be activities recorded in the cold season over the Fennoscandian Peninsula and its relationship with the Arctic polar vortex through model simulations. For this purpose, we conduct 7 Be simulations for the period of January-March 2003 using the GEOS-Chem global 3-D chemical and transport (CTM) model. The period was selected because of the large number of events with extremely high 7 Be concentrations at the surface in Fennoscandia; some of these events were preceded by very low surface concentrations in the lower troposphere (< 10th percentile). This period thus offers the opportunity to test the hypothesis that SSW events facilitate a fast descent of not only the midlatitude but also vortex air (Ajtić et al., 2018). To achieve this goal, our analysis will therefore focus on the following: investigating the processes responsible for the variability of 7 Be concentrations in surface air in northern Europe; better understanding whether and how SSW and the Arctic vortex wintertime instability influence the surface concentrations of 7 Be in northern Europe; quantifying the rate of air subsidence on the inner and outer side of the vortex during the period of its instability.
To analyze the influence of SSW and of the Arctic polar vortex on 7 Be concentrations, we first assess the performance of the GEOS-Chem model in reproducing the observed 7 Be variability. We then use model simulations together with other supporting measurements from soundings and meteorological datasets to examine the processes responsible for the variability in the 7 Be concentrations over the period of January-March 2003. As opposed to the cosmogenic origin of 7 Be, 210 Pb (halflife 22.3 years) is a nuclide of crustal origin derived from decay of 222 Rn (half-life 3.8 d), which is emitted from soils by decay of 226 Ra. Owing to the contrasting natural origins of the two nuclides, the 7 Be / 210 Pb ratio is often regarded as indicative of vertical transport processes and convective activity in the atmosphere (e.g., Koch et al., 1996;Tositti et al., 2004;Brattich et al., 2017a, b). After being produced by contrasting physical mechanisms, both 7 Be and 210 Pb rapidly attach to ambient submicron-sized particles (e.g., Gaffney et al., 2004) and are removed by wet (mainly) and dry (secondarily) deposition processes of their carrier aerosol. The bias in the simulated 7 Be / 210 Pb ratio due to uncertainties in the model deposition schemes is thus reduced. For this reason, besides 7 Be, the 7 Be / 210 Pb ratio was also analyzed to gain further insights into vertical transport processes during the study period.
The rest of this paper is organized as follows. Section 2 describes the radioactivity ( 7 Be and 210 Pb) and meteorological data used. Section 3 provides a brief description of GEOS-Chem, the HYSPLIT trajectory model, and statistical parameters used to assess the model's performance in reproducing the observations. Section 4 presents an overview of the 7 Be observations made in northern Europe in the boreal winter 2003. Section 4 presents and discusses the main results of this work. Specifically, Sect. 4.1 analyzes the precipitation and transport pattern in the study region, while Sect. 4.2 assesses how well the GEOS-Chem model performs in reproducing the observed variability in the monthly mean surface 7 Be concentrations during the study period. Section 4.3 further evaluates the performance of the model in reproducing the short-time variability of 7 Be in northern Europe, followed by an interpretation of the observed variability using model simulations and additional meteorological and ozone observations in Sect. 4.4. Finally, summary and conclusions are given in Sect. 5.

Data
In this section we briefly describe the 7 Be and 210 Pb radioactivity data as well as the meteorological datasets analyzed in this work.

7 Be data
Since 1988, the Radioactivity Environmental Monitoring data bank (REMdb) (https://data.jrc.ec.europa.eu/collection/ id-0117, last access: 23 November 2021) has brought together and stored in a harmonized way environmental radioactivity data (air, water, milk and mixed diet) measured by the European Union member states (Sangiorgi et al., 2019). Among the set of sample types and measurements recommended in 2000/473/Euratom (European Commission, 2000), measurements of natural radioelements, such as 7 Be in surface air, are required, and hence, it is very closely monitored and widely stored in the REMdb (De Cort et al., 2007).
Within the REMdb, the activity concentration of 7 Be in the surface air in northern Europe (latitude north of 55 • N) is available for six surface sampling sites (Hernández-Ceballos et al., 2015), all located at elevations ranging from 9 to 130 m a.s.l.: Ivalo, Umeå, Helsinki, Kista, Harku and Risø, however, with varying start dates and sampling frequencies (Fig. 1a). The largest dataset is for Helsinki where, since 1999, the sampling has been performed daily or once every 2 d. Datasets for Ivalo, Umeå, Kista and Risø also span more than two decades and have a good temporal coverage (roughly once a week since 1995). In particular, 7 Be activity concentrations were obtained by gamma-spectrometry analysis performed by the European Union Competent Authorities. Aerosol samples were collected on filter papers using air samplers with a flow rate of several hundred cubic meters per day, and then their radioactivity concentrations were analyzed in laboratories.

210 Pb data
Daily aerosol samples were collected in Helsinki on the roof of the Finnish Meteorological Institute's main building (60 • 10 N, 24 • 57 E). Filters (Munktell MGA, diameter Ø = 240 mm) were changed every day at 06:00 UTC. The air volume was about 3500 m 3 d −1 . The filters were assayed for 210 Pb by alpha counting of the in-grown daughter nuclide 210 Po (Mattsson et al., 1996).

Meteorological and ozone data
The 7 Be variability is tightly linked to horizontal and vertical transport of the carrier aerosol, and to precipitation that leads to the radionuclide's removal from the atmosphere. Here we used the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis (Gelaro et al., 2017) to assist in the data analysis and to drive the GEOS-Chem model simulations. MERRA-2 is produced with version 5.12.4 of the Goddard Earth Observing System (GEOS) atmospheric data assimilation system. It assimilates modern observations of the atmosphere, ocean, land and chemistry and includes assimilation of aerosol remote sensing data.
Vertical soundings of air temperature from the Finnish Meteorological Institute's (FMI) Arctic Space Centre (http: //fmiarc.fmi.fi, last access: 23 November 2021) at Sodankylä, northern Finland (67.37 • N, 26.63 • E) were obtained from the University of Wyoming (http://weather. uwyo.edu/upperair/sounding.html, last access: 23 November 2021). To study the effect of downward transport of stratospheric air masses into the troposphere, potential vorticity (PV) values (Holton et al., 1995) were calculated from ERA-Interim wind, temperature, and surface pressure fields (Dee et al., 2011) obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF), Reading, UK.

Methods
In this section we give a brief description of the GEOS-Chem and HYSPLIT models and the statistical parameters used to indicate the model performances.
We used the 7 Be production rates recommended by Lal and Peters (1967) for a maximum solar activity year (1958), which has been shown to produce the best results compared to aircraft 7 Be observations in the stratosphere where 7 Be concentrations mainly result from a balance between production and radioactive decay, and their observations can be used as a constraint on the 7 Be source (Koch et al., 1996;. In this work, the production rates of Lal and Peters (1967) are formulated as a function of latitude and pressure without seasonal variation (Koch et al., 1996). About twothirds of atmospheric 7 Be is generated in the stratosphere. While there are more accurate modern models (e.g., Masarik and Beer, 1999;Webber et al., 2007;Usoskin and Kovaltsov, 2008;Poluianov et al., 2016), the use of the Lal and Peters production rates is sufficient for this study where we mainly focus on atmospheric transport features. Although the Lal and Peters sources are used here for temporal variability studies, they are not applicable for quantitative studies of the 7 Be activities (e.g., Golubenko et al., 2021) 222 Rn emission follows a recent work by Zhang et al. (2021), in which a customized emission map was built upon a few previously published emission scenarios and evaluated against global 222 Rn surface observations and aircraft profiles. 222 Rn emission flux rate is a function of latitude, longitude and month. 7 Be and 210 Pb are assumed to behave like aerosols once formed in the atmosphere and subject to dry and wet deposition . Both wet and dry deposition for 222 Rn are neglected due to its inert nature.
GEOS-Chem simulations in this work are driven by the MERRA-2 meteorological reanalysis. The native resolution of MERRA-2 is 0.667 • longitude by 0.5 • latitude, with 72 vertical layers (top at 0.01 hPa). The meteorological fields are regridded into 2.5 • longitude by 2 • latitude for the GEOS-Chem simulations in this work. GEOS-Chem uses the TP-CORE advection algorithm of Lin and Rood (1996). Convective transport is calculated using archived convective mass fluxes (Wu et al., 2007). Boundary-layer mixing is based on the non-local scheme implemented by Lin and McElroy (2010). The wet deposition scheme follows that of  and includes rainout (in-cloud scavenging) due to stratiform and anvil precipitation, scavenging in convective updrafts (Mari et al., 2000), and washout (below-cloud scav- enging) by precipitation (Wang et al., 2011). Precipitation formation and evaporation fields are archived in MERRA-2 and used directly by the model wet deposition scheme. Dry deposition is based on the resistance-in-series scheme of Wesely (1989).
In addition to the standard model simulations of 7 Be and 210 Pb, we separately transport 7 Be produced in the model layers above the MERRA-2 thermal tropopause (i.e., stratospheric 7 Be tracer) to quantify the stratospheric contribution to 7 Be in the troposphere. This approach was previously used by Liu et al. (2001Liu et al. ( , 2016. Stratospheric fraction of 7 Be is defined as the ratio of the stratospheric 7 Be tracer concentration to the 7 Be concentration from the standard simulation. All model simulations are conducted for the period of January 2002-March 2003 with initial conditions from a previous 5-year simulation. Hourly and monthly mean outputs for January-March 2003 are used for analysis.

HYSPLIT
The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model, developed by the NOAA's Air Resources Laboratory (ARL) (Stein et al., 2015), was used to calculate a set of backward trajectories during the study period. To compute the 96 h 3-D backward trajectories at 00:00, 06:00, 12:00 and 18:00 UTC and with different ending heights, 100, 500, 1000 and 1500 m above ground level, the NCEP (National Centers for Environmental Prediction) FNL Operational Global Analysis (NCEP/NWS/NOAA/U.S. Department of Commerce, 2000) meteorological files were used. Computation used the vertical velocity (w, m/s) field contained in the meteorological input file. While 96 h was considered a sufficiently long period to represent the synoptic air flows, the heights were selected to help us to understand the behavior of the airflows circulating in the atmospheric boundary layer (ABL), just above the ABL and in the free troposphere. We used the cluster methodology implemented in the HYSPLIT model to group the calculated trajectories according to their length and curvature and thus identify the airflow patterns over the whole period of the analysis (Hernández-Ceballos et al., 2013;Brattich et al., 2016). It is worth mentioning that clusters, as well as trajectories, indicate an estimation of the general airflow rather than the exact pathway of an air parcel (e.g., Jorba et al., 2004;Salvador et al., 2008).

Evaluation of the model output
The performance of the model in reproducing observed activity concentrations was evaluated by calculating some basic statistical parameters, such as the mean and standard deviation and other indicators, according to the methodology developed by Hanna (1993) and summarized later by Chang and Hanna (2004). Specifically, the performance of the CTM was evaluated using the following set of indicators, proposed by Carruthers et al. (2000): -The mean bias (MB), a measure of the mean difference between the modeled and observed concentrations: where C m is modeled concentration and C o is observed concentration.
-The normalized mean square error (NMSE), a measure of the mean difference between matched pairs of modeled and observed concentrations: -The fraction of modeled concentrations within a factor of 2 of observations (FA2), i.e., for which 0.5 < C m /C o < 2.  -The Pearson's correlation coefficient (R), a measure of the extent of a linear relationship between the modeled and observed concentrations: Additionally, the statistical significance of the Pearson's correlation coefficient and the Z test were also evaluated to identify the presence of statistically significant differences between observed and simulated concentrations. A perfect model has MB and NMSE values equal to 0 and FA2 value equal to 1, while the R results range from −1 (perfect negative relationship) to +1 (perfect positive relationship), where 0 implies no relationship between the variables. To better understand the quantitative differences between observations and simulations, scatter plots were used.  As indicated by Ajtić et al. (2018) and Bianchi et al. (2019), the winter of 2003 offers a good opportunity to investigate a possible link between SSW events and extreme surface concentrations of 7 Be detected in northern Europe. This period is sufficiently covered by the 7 Be activity concentration measurements at all six monitoring sites.
In particular, very high 7 Be activity concentrations, above the 90th percentile (calculated over the 1995-2011 period) simultaneously at most of the Fennoscandian Peninsula sampling sites, were recorded around 23-24 February 2003 (Ajtić et al., , 2018 (Fig. 1b). During the 2002-2003 boreal winter, the Arctic vortex was relatively unstable, with six SSW events taking place over the whole season (Peters et al., 2010). Two very pronounced episodes, which were both associated with the vortex splitting and fast SSW recovery, occurred in January and February, respectively (Günther et al., 2008). The evolution of the vortex caused vortex filamentation and vigorous mixing of the vortex and midlatitude stratospheric air (Günther et al., 2008;Müller et al., 2003). Several  Huret et al., 2006;Müller et al., 2007). Such disturbances around the pole are expected to affect the troposphere, i.e., weather conditions (Baldwin and Dunkerton, 2001), and air chemical composition (Hsu, 1980;Limpasuvan et al., 2004). Hence, the high 7 Be concentrations that were measured in Fennoscandia around 24 February 2003 could be a result of downward motion of midlatitude stratospheric air. Interestingly, prior to this episode, very low (below the 10th percentile for each site) surface concentrations of 7 Be were measured in Risø, Kista and Ivalo on 3, 10 and 16 February 2003, respectively (Fig. 1). These low values were tentatively linked by Ajtic et al. (2018) with the transport of aged stratospheric vortex air poor in 7 Be, even though they are more likely related to precipitation scavenging that occurred in the days before, as shown by the ECA&D (European Climate Assessment & Dataset, https://www.ecad.eu/, last access: 23 November 2021) records. The reader is referred to Ajtić et al. (2018) for more details.

Analysis of winter precipitation and transport in the Fennoscandian Peninsula: observations vs. model simulations
Before analyzing the temporal pattern of simulated 7 Be concentrations, we analyzed the precipitation and transport pattern in the MERRA-2 meteorological dataset that drives the GEOS-Chem simulations. In particular, the MERRA-2 precipitation was evaluated against the data from Global Precipitation Climatology Project (GPCP) v2.2 (https://rda.ucar. edu/datasets/ds728.2/, last access: 23 November 2021) satellite and surface observations in winter 2003 (Adler et al., 2003). Figure 2 shows the MERRA-2 and GPCP monthly precipitation in winter for the region within 0-90 • N and 90 • W-90 • E. Good agreement is found between the MERRA-2 and the GPCP precipitations averaged over the region. Specifically, the geographical distribution of precipitation in MERRA-2 shows some important features that are consistent with the observed climatology precipitations: the desert climate in northern Africa with very low precipitation through-E.  Be concentration events at the surface in northern Europe out the period, high precipitation over the North Atlantic region during winter and the seasonal pattern of precipitation in Europe which is similar to that in the North Atlantic region.
To assess the capability of the model to correctly capture the trend in precipitation during the observation period at the sampling sites, we examined the normalized differences between the MERRA-2 and the observed precipitation, calculated as a difference between the MERRA-2 and the observed values, normalized over the observed value (Table 1).
Overall, the MERRA-2 precipitation tends to be higher than that of GPCP at all sampling sites (Table 1) except for Harku and Helsinki, and especially in the February-March period. This result is in agreement with the findings of Gelaro et al. (2017), who compared the global precipitation of MERRA-2 and GPCP and reported a general positive bias over northern high latitudes. However, the agreement between MERRA-2 and GPCP precipitation seasonality is reasonable, as indicated by the correlation coefficient values, higher than 0.85 at all sites except for Ivalo (−0.32), and the low NMSE values, in the range of 0-0.42 (Table 1). The low negative correlation at Ivalo is due to the fact that while the GPCP-observed precipitation at this site is similar between January and February with a general tendency towards lower values from January to March 2003, the model simulates a decrease from January to February with a small increase in March. However, the statistical parameters reported in Table  1 indicate a small overall discrepancy between the GPCP and MERRA-2 precipitation at all sites. Figure 3 shows that winter circulation in the Fennoscandian Peninsula is dominated by SW and W winds (Chen, 2000;Linderson, 2001). The analysis of the main circulation in the 3 months in Fig. 3 reveals low wind speeds from S-SW in the study area and period. A region of strong wind speeds, possibly corresponding to the Arctic vortex, is clearly visible at surface level to the west of the study area in all 3 months. In addition, there appears to be a convergence area (opposite wind directions) between 60 and 75 • N. Model-simulated 7 Be / 210 Pb ratios and fraction of stratospheric 7 Be increased over the 3-month period and peaked in March, suggesting increasing stratospheric influence, subsidence or convective mixing in the study region.

Variations of the monthly mean surface 7 Be
concentrations in the Arctic region: model simulations vs. observations Figure 4a shows a scatter plot comparing the simulated and observed monthly mean 7 Be concentrations at the six sampling sites. Table 2 reports the statistical parameters and the normalized differences that indicate the performance of the GEOS-Chem model in reproducing the observed 7 Be monthly means.
In general, the model simulates well the month-to-month trend in 7 Be concentrations measured at the sampling sites, as indicated by the fact that all the values fall within the 95 % confidence levels (Fig. 4a) and the high positive correlation coefficients (>0.7) except for Ivalo and the low MB and NMSE values (Table 2). In fact, the normalized differences are not very high (generally <1), except at Risø. The bias between the model and the observations is partly attributed to the coarse horizontal resolution of the model. Overall, the simulations underestimate the observed values, likely due to uncertainties associated with the deposition schemes and/or precipitation as discussed earlier.
The use of the 7 Be production rate of Lal and Peters (1967) for a solar maximum year (1958) may also partly explain the tendency of simulated 7 Be to be lower than observed. The sunspot number in 2003 (99.3) was rather low (slowly decreasing from 2000, a solar maximum year, and reaching minimum in 2008) compared with the value of 184.8 in 1958. As is known, the galactic cosmic-ray intensity, largely responsible for the production of cosmogenic radionuclides, at the Earth's orbit is inversely related to solar activity (Potgieter, 2013), leading to the well-known phase opposition between sunspot number and 7 Be concentration (e.g., Hernández-Ceballos et al., 2015). Sunspot number data herein used were extracted from the World Data Center for the production, preservation and dissemination of the international sunspot number (Sunspot Index and Long-term Solar Observation, SILSO, Royal Observatory of Belgium, Brussels, https://wwwbis.sidc.be/silso/datafiles, last access: 23 November 2021).

Variations of the 7 Be weekly and daily mean surface concentrations in the Arctic region: observations vs. model simulations
After analyzing the model's performance in reproducing 7 Be monthly mean observations in the previous section, here we compare the simulated and observed weekly (daily in the case of Helsinki) 7 Be activity concentrations at the six sampling sites (Fig. 4b). Table 3 shows the corresponding parameters that indicate the performance of the GEOS-Chem model in reproducing observations. The weekly evolution of simulated versus observed 7 Be concentrations at these sites is shown in Fig. 5. As with the monthly means, the model generally represents adequately the temporal pattern but not the magnitude of weekly mean concentrations, which tend to be lower than those observed (Table 3, Figs. 4b and 5). This bias can arise from the higher precipitation in the model than in the observations and/or from errors in the deposition schemes. However, the correct reproduction of the 7 Be temporal pattern, as indicated by the high correlation values at all sampling sites with the exception of Risø, suggests that the model captures the transport processes leading to the peak in 7 Be concentrations at the end of February 2003 and the preceding very low concentration values. In addition, the low MB and NMSE values calculated at all sites and especially at Ivalo and Umeå suggest that the model reproduces adequately the observed Table 1. Normalized differences between the MERRA-2 and observed precipitation, calculated as differences between the MERRA-2 and the observed values, normalized over the observed value, at each sampling site (left). Positive values indicate that the model tends to overestimate observations while the negative ones indicate underestimation. Statistical parameters (mean ± SD = mean ± standard deviation; MB = mean bias; NMSE = normalized mean square error; R = correlation coefficient; FA2 = factor of 2) indicating the model performance in reproducing GPCP monthly accumulated precipitation at the six sampling sites in northern Europe (right).  (Fig. 5). The significance of the Z test for the differences between mean observed and simulated values indicates the presence of statistically significant differences, reflecting the model's general tendency to underestimate observations. As for the 7 Be / 210 Pb ratio in Helsinki (Table 3), the model tends to underestimate the observed ratio, which could be due to the model underestimating 7 Be and overestimating 210 Pb. Nevertheless, the relatively high correlation between the simulated and observed ratios suggests a reasonable simulation of the temporal pattern of this tracer.

Understanding the 7 Be variations during the 2002-2003 boreal winter
As mentioned earlier (Sect. 4.1), an SSW event occurred at the end of February 2003. We concentrated our analysis on   7 Be / 210 Pb ratio, maps of surface winds and relative humidity, ozone soundings, and vertical cross sections of simulated 7 Be activity concentrations and calculated potential vorticity, and simulated and observed vertical pro- files of air temperature. The results were further supported with the analysis of the clusters of back-trajectories during the two different periods of low and high 7 Be concentrations. Figure 6 presents the temporal (weekly mean) pattern of 7 Be / 210 Pb and of the stratospheric fraction of 7 Be (calculated as the ratio of the stratospheric 7 Be tracer concentration to the total 7 Be concentration in the troposphere) at the six sampling sites, while daily observations of the 7 Be / 210 Pb ratio at Helsinki and Sodankylä (67.367 • N, 26.629 • E; 160 km south of Ivalo) are presented in the Supplement (hereafter SI). At the beginning of February, the ratio was generally quite low at all the sites. In contrast, the week of 19-26 February 2003 was marked by an evident peak in the 7 Be / 210 Pb ratio and a simultaneous increase in the fraction of 7 Be originating in the stratosphere at all sites, which together could be the first indication of a prominent vertical transport from the UT-LS region. The low 7 Be / 210 Pb ratio at Risø results from the simultaneous increase in both 7 Be and 210 Pb concentrations, together with the delayed stratospheric influence at this site as evidenced by the pattern of the stratospheric 7 Be fraction (Fig. 6b).
Hence, we further examined the vertical profiles of temperature with an aim to identify differences in vertical transport near the beginning and end of February. The soundings from the Sodankylä station in the Arctic offer three sets of measurements for each of the investigated periods: on 10 E.  Be concentration events at the surface in northern Europe and 16 February, which fall into the period when very low 7 Be concentrations were recorded in Kista and Ivalo, respectively; 22 and 24 February, the days marked by extremely high 7 Be concentrations over the Fennoscandian Peninsula; and 20 and 21 February in the period of transition to high 7 Be concentration over the Fennoscandian Peninsula. Figure 7 shows air temperature profiles in the MERRA-2 dataset and atmospheric soundings at the Sodankylä station. Besides the very good agreement between the MERRA-2 and observed temperatures, a warming of the stratosphere (20-60 km) and a different vertical temperature structure of the lower stratosphere around 20-24 February as compared to the 10 and 16 February profiles are also evident. The observation of the increase in stratospheric temperatures suggests the link between the SSW and the 7 Be peak observed at the six sampling sites located in northern Europe. In addition, the ozone soundings at the Sodankylä station reveal an ozone mixing ratio peak in the lower troposphere (∼ 1.5-3 km) on 19 February 2003 as compared to those observed during 12, 26 and 28 February 2003 (Fig. 8a), consistent with downward transport from higher altitudes around that day. Despite the chemical ozone loss in the Arctic vortex in the stratosphere in 2003 as observed by ozone soundings (Tilmes et al., 2006), obviously lower-stratospheric ozone was still enhanced rel-ative to tropospheric ozone. In addition, average O 3 values recorded at surface air quality stations located in Denmark, Finland and Sweden, which are available through the saqgetr R package (Grange, 2019), show enhanced O 3 concentrations in late February 2003, consistent with the aforementioned peaks in the 7 Be / 210 Pb ratio as well as stratospheric 7 Be fraction. This further suggests the transport of stratospheric air masses to the surface.
Simultaneously, the analysis of maps of surface transport and relative humidity (Fig. 9) highlights the different winds and relative humidity values in the two periods, with low relative humidity values (∼40 %-50 %) suggesting subsidence from 18 to 21 February, and the transition from a clockwise circulation to the fast and complex wind system typical of the second period corresponding to the high 7 Be peak and the SSW. Together with the increase in stratospheric temperatures observed in Fig. 7, the reversal of zonal winds confirms the link between the 7 Be peak and the SSW event, in agreement with Sofieva et al. (2012). Interestingly, the parcel of the lowest relative humidity values occurs during the 18-21 February period, i.e., a couple of days before the dates of the 7 Be peaks in the measurements and those peaks in the simulated 7 Be / 210 Pb ratio and stratospheric 7 Be fraction (Figs. 1b and 6). Together with enhanced ozone concen- trations observed in the lower troposphere on 19 February, this suggests that the downward transport from the UT-LS was triggered by the SSW occurring a few days before.
To better constrain the stratospheric origin of the air masses arriving at the sampling sites during the two periods, we further analyzed the potential vorticity data from ECMWF during the month of February 2003 at three latitudes (63, 64.5 and 66 • N) along the 21 • E meridian (Fig. 10), while vertical cross sections of MERRA-2 PV are provided in Fig. S3. The data clearly reveal a bubble of high potential vorticity down to the surface at the three latitudes from 18-22 February 2003, particularly at the northernmost latitude where values higher than 1.6 PVU, a value considered as a threshold for stratospheric air in the lower troposphere especially when in conjunction with low relative humidity, high 7 Be / 210 Pb ratios and ozone (Cristofanelli et al., 2006), were observed. Indeed, PV is considered as a quasi-passive tracer, because of its conservation properties under adiabatic conditions (Hoor et al., 2010;Gettelman et al., 2011), and the tropopause level can be identified as regions of strong enhancements in gradients of PV, essentially indicated by distinct values in the troposphere and stratosphere. Recent studies of PV cutoff lows in the Northern Hemisphere extratropics (Pinheiro et al., 2017;Portmann et al., 2021;Muñoz et al., 2020) suggest that the majority of the events are relatively short-lived, persisting for about 2-3 d, a value in good agreement with the vertical motion subject of this study.
The low relative humidity (Fig. 9) and high potential vorticity (Fig. 10) corresponded to high 7 Be descending to lower atmospheric levels, as simulated by the model (Fig. 11). The descending vertical motion from the upper vertical levels during the period is clearly visible in the MERRA-2 vertical pressure velocity (omega) fields sampled at the six sampling sites for the month of February 2003, especially for the northernmost sites (Fig. 12) where omega is largely positive with near-surface values up to 0.3-0.4 Pa s −1 around 18-19 February. Further evidence of this is seen from the maps of vertical pressure velocity (Fig. 13) and of the stratospheric fraction of 7 Be originating in the stratosphere (Fig. 14) at 940 hPa in the boundary layer. Note that positive vertical pressure velocity as seen over both Fennoscandia and Greenland (Fig. 13)  dicates descending motions that can facilitate the transport of stratospherically influenced air, if present, to the groundlevel. This explains why Fennoscandia saw increased stratospheric influence on surface 7 Be concentrations during this period (Fig. 14), but Greenland did not see much. However, on monthly average, Greenland is a region with significant stratospheric influences in February 2003 (middle right panel, Fig. 3).
Analysis and comparison of 4 d back-trajectories at each sampling station allowed a reconstruction of two distinct atmospheric circulation patterns in the two periods. Figure 15 shows the clustering results for three sampling sites (Ivalo, Harku and Risø) during both periods. The stations are ordered as a function of decreasing latitude from high (upper panels) to low (lower panels). Only results at 1000 m (well above the winter PBL height) are shown, and cluster results for other altitudes in the lower troposphere are similar. While the first period (with low 7 Be values) presents the dominance of westerly winds (air masses flowing eastward), as is typical at these latitudes (Zanis et al., 1999), the second period is characterized by a clockwise displacement of airflows with origins at higher altitudes (Fig. 15). This pattern in the second period is better established at lower-latitude stations (Helsinki, Harku, Kista and Risø) compared to higherlatitude ones (Ivalo, Umeå). It likely results from the aged vortex (Günther et al., 2008) and the SSW at the end of February, corresponding to a decrease in the MERRA-2 daily average height of the thermal tropopause on 21-22 February at the Sodankylä station in Finland (Fig. S4) during the SSW (Peethani et al., 2014;Wargan and Coy, 2016). Associated with these processes is the downward transport of stratospheric air previously identified with an independent approach.

Summary and conclusions
We have used a global 3-D model (GEOS-Chem) driven by the MERRA-2 meteorological reanalysis to simulate atmospheric concentrations of 7 Be of cosmogenic origin for the period of January-March 2003. The aim was to verify the mechanisms responsible for the surface 7 Be variabilities in northern Europe, and to test the hypothesis that SSW events may facilitate fast descent of UT-LS vortex air to the surface. The period was selected as it involves two intense SSW events and observations of extremely high 7 Be concentrations at six sampling sites in Fennoscandia. Before using the model's output to investigate the processes responsible for 7 Be variability in northern Europe over the period, we evaluated the MERRA-2 precipitation fields against the GPCP satellite and surface observations. A generally good agreement was found both at regional scale and at the six sampling sites. Analysis of the wind fields in the study period indicated low wind speeds from S-SW in agreement with the major circulation patterns over the Fennoscandian Peninsula in winter, and the presence of a region of strong wind speeds to the west of the study area, likely in connection with the Arctic polar vortex.
The model reproduced efficiently the 7 Be and 7 Be / 210 Pb temporal (i.e., monthly and weekly) patterns at the six sampling sites in the study period, even though it presented a tendency to underestimate the observed surface 7 Be concentrations. The lower modeled values are likely due to its coarse horizontal resolution (2.5 • longitude by 2 • latitude), lack of year-to-year variation in 7 Be production rates and uncertainties associated with precipitation scavenging.
In order to investigate the processes responsible for 7 Be variability at the six sampling sites during the study period, and in particular to test whether the peak 7 Be concentrations measured in Fennoscandia around 24 February 2003 originated from fast descent of stratospheric vortex air facilitated by SSW, we analyzed time-height cross sections of simulated 7 Be and potential vorticity, vertical profiles of air temperature, maps of surface winds and relative humidity, and ozone soundings. The analysis of the temporal vari- ations of simulated 7 Be / 210 Pb ratio and fraction of 7 Be originated in the stratosphere indicates a peak during the week of 19-26 February 2003, suggesting downward transport from the UT-LS region. The latter was corroborated by a layer of ozone mixing ratio enhancements in the lower troposphere recorded by the soundings at the Sodankylä station on 19 February. Furthermore, the vertical profiles of air temperature indicated a warming of the stratosphere and a change in shape in the vicinity of the tropopause region during the period, suggesting the link between the downward transport of the vortex air and SSW.
Our analysis of time-height cross sections of simulated 7 Be concentrations, calculated potential vorticity and MERRA-2 vertical pressure velocity (omega) revealed the vertical downward transport to the surface of a stratospheric air parcel characterized by high potential vorticity, high vertical velocity (in particular at Ivalo on 19 February) and high 7 Be concentrations, further supporting the stratospheric origin of the air masses during the investigated period.
Additionally, low relative humidity and a change in the circulation pattern from slow and clockwise to fast, swirling winds occurred over the study area. The change in the circulation pattern and the downward transport of stratospheric E.  Be concentration events at the surface in northern Europe air was verified by the analysis of the clusters of backtrajectories during the periods of low and high 7 Be concentrations, which showed a change from westerlies to airflows from upper vertical levels.
Altogether, these analyses confirm the link between the SSW and transport of stratospheric air to the surface, resulting in high surface 7 Be concentrations observed in February 2003 in Fennoscandia. Since more frequent SSW events are expected in a warmer climate Tziperman, 2017, 2018;Simpkins, 2017), this link has important implications for the impact of climate change on atmospheric transport, tropospheric composition and air quality in northern highlatitude regions.
Author contributions. JA, MAHC and EB designed the study. HL and BZ conducted the GEOS-Chem model simulations. MAHC led the calculation and analysis of HYSPLIT back-trajectories. EB developed the analysis methodology and led the analysis of observational data and model output, with contributions from all coauthors. JP contributed 210 Pb and meteorological observational datasets. EB wrote the paper with contributions from all coauthors.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. The University of Wyoming, ECMWF, Dr. Rigel Kivi from the Finnish Meteorological Institute Arctic Space Center and Dr. Laura Thölix from the Finnish Meteorological Institute are gratefully acknowledged for providing and helping with temperature profiles, ozone soundings and potential vorticity data used to support the detection of air masses of stratospheric origin in the troposphere. The World Data Center is gratefully acknowledged for the production, preservation and dissemination of the international sunspot number (Sunspot Index and Long-term Solar Observation, SILSO, Royal Observatory of Belgium, Brussels, http://www.sidc.be/silso/ datafiles#total, last access: 23 November 2021). Support for the Twentieth Century Reanalysis Project version 3 dataset is provided by the U.S. Department of Energy, Office of Science Biological and Environmental Research (BER; http://science.energy.gov/ber/, last access: 23 November 2021), by the National Oceanic and Atmospheric Administration Climate Program Office, and by the NOAA Physical Sciences Laboratory. NASA/GSFC/SED/ESD/LA/MAPB and contributors are gratefully acknowledged for developing and computing the GPCP (Global Precipitation Climatology Project) as part of the Global Energy and Water Cycle Exchanges (GEWEX) project. Stuart Grange is gratefully acknowledged for providing access to European air quality data through the saqgetr R package. Review statement. This paper was edited by Amanda Maycock and reviewed by three anonymous referees.