Observations of ozone depletion events in a Finnish boreal forest

We investigated the concentrations and vertical profiles of ozone over a 20-year period (1996–2016) at the SMEAR II station in southern Finland. Our results showed that the typical daily median ozone concentrations were in the range of 20–50 ppb with clear diurnal and annual patterns. In general, the profile of ozone concentrations illustrated an increase as a function of heights. The main aim of our study was to address the frequency and strength of ozone depletion events at this boreal forest site. We observed more than a thousand of 10 min periods at 4.2 m, with ozone concentrations below 10 ppb, and a few tens of cases with ozone concentrations below 2 ppb. Among these observations, a number of ozone depletion events that lasted for more than 3 h were identified, and they occurred mainly in autumn and winter months. The low ozone concentrations were likely related to the formation of a low mixing layer under the conditions of low temperatures, low wind speeds, high relative humidities and limited intensity of solar radiation.


Introduction
Ozone (O 3 ) is one of the most important compounds in atmospheric chemistry (e.g.Monks, 2005;Monks et al., 2015).It is a strong greenhouse gas and a key contributor to the oxidation capacity of the atmosphere (Toumi et al., 1994;Simpson et al., 2007).Ozone is also the principal source of the essential atmospheric oxidants, being the major photochemical precursor of the OH radical (Atkinson and Arey, 2003) and a source of the nocturnal NO 3 radical through the NO 2 + O 3 → NO 3 reaction (Brown and Stutz, 2012).Very recently, the importance of ozone in initiating the au-toxidation of volatile organic compounds (VOCs) in the gas phase was recognized, producing rapidly a range of highly-oxidized multifunctional compounds (HOMs, Ehn et al., 2014;Rissanen et al., 2014), which can participate in the formation of secondary organic aerosols and/or contribute significantly on their early growth (Kulmala et al., 1998;Ehn et al., 2014;Trostl et al., 2016).While stratospheric ozone acts to shield the Earth from detrimental ultraviolet radiation, the tropospheric ozone is harmful for human health (e.g.Shang et al., 2013;Turner et al., 2016).It causes negative consequences throughout the biosphere (Gregg et al., 2003;Ashmore, 2005) and reduces crop yields (Emberson et al., 2001;Feng and Kobayashi, 2009), when it is either deposited to the surface of or taken up and then metabolized by plants.
The concentration of tropospheric ozone is governed by its sources and sinks.In the troposphere, O 3 is produced via complex photochemical cycles involving VOCs and nitrogen oxides (NO x consists of NO and NO 2 ) (Sillman et al., 1990;Sillman, 1999).Transport from the stratosphere brings additional ozone to its tropospheric budget.The tropospheric ozone is reduced directly by dry deposition, photolysis and reactions with other trace gases and indirectly via wet deposition of ozone precursors.The dry deposition of ozone can be enhanced under moist conditions (Fuentes et al., 1992).These contributions form an equilibrium state that maintains tropospheric ozone-mixing ratios at concentration levels of about tens of parts per billion (ppb) (Vingarzan, 2004;Monks et al., 2015).Although tropospheric ozone is relatively well mixed in a broad sense, its concentration can vary notably with both time and space depending on the prevailing meteorological conditions.Regional and local surface O 3 con-X.Chen et al.: Observations of ozone depletion events in a Finnish boreal forest centrations are affected, for example, by horizontal and vertical advections of pollutants, entrainment from the free troposphere and other boundary layer dynamic processes (Simpson et al., 2007;Cao et al., 2016;Tang et al., 2017;Wang et al., 2017;Zhou et al., 2017).
The most well-known ozone depletion phenomena are the stratospheric "ozone-hole" episodes taking place above the polar regions caused by the catalytic halogen chemistry (Molina and Rowland, 1974;Solomon, 1999), as well as the frequently reported O 3 depletion events in a polar planetary boundary layer (Simpson et al., 2007).In addition, significant decreases in local surface O 3 concentrations have been observed in polluted urban environments (Gregg et al., 2003).While in the polar areas halogenated compounds are responsible for most of the surface ozone destruction, in urban locations low local ozone concentrations are typically related to complex atmospheric chemistry involving anthropogenic pollutants, mainly nitrogen oxides (Gregg et al., 2003;US National Research Council, 1992).Also, low surface ozone concentrations were reported at different sites in the Amazonian forest (Rummel et al., 2007;Pacifico et al., 2015).Although anormal reductions in surface ozone concentrations have been observed at different geographical locations, there lacks a general characterization of the vertical features of such O 3 depletion phenomena.Thus, it is imperative to study the vertical profile of O 3 to provide insights into observations of surface ozone depletion.
Based on previously reported depletion events in surface ozone, we investigated in detail ozone concentrations and profiles using long-term (> 20 years) data collected between 1 January 1996 and 30 June 2016 at the Hyytiälä SMEAR II station, a boreal forest site in Finland.Our main goals in this study were (1) to characterize the typical seasonal and diurnal behaviour of vertical ozone concentration profiles within and above the boreal forest, (2) to evaluate statistically the low ozone concentration episodes, (3) to assess the atmospheric conditions in association with these lowconcentration episodes, and (4) to identify key atmospheric features during the ozone depletion events.

Experimental setup and data analysis
The SMEAR II (Station for Measuring forest Ecosystem-Atmosphere Relation) infrastructure, located in a boreal forest at Hyytiälä, southern Finland (61 • 51 N, 24 • 17 E, 181 m a.s.l.), offers a unique set of information based on continuous measurements of atmospheric composition and meteorological variables (Hari and Kulmala, 2005).These measurements are located at several heights on a mast, ranging from below-canopy levels (4.2, 8.4, and 16.8 m) to abovecanopy levels (33.6, 50.4, 67.2, 101, and 125 m), capturing the vertical profile of the lower atmosphere at this site.Snow on the ground during winter months is common at our measurement site, and the site can receive transported pollutants from nearby cities (e.g.Tampere and Helsinki) and continental areas (Riuttanen et al., 2013).
The forest stand is dominated by Scots pines (Pinus sylvestris L.) with an average height of 16-18 m (living canopy height ca.10-12 m) and age of 55 years.Other, less abundant species include Norway spruce (Picea abies L. H. Karst.), silver birch (Betula pendula Ehrh.) and trembling aspen (Populus tremula), which together form ca. 25 % of the biomass (Ilvesniemi et al., 2009).The ground vegetation consists of dwarf shrubs (e.g.Vaccinium sp., Calluna vulgaris) and mosses (Pleurozium schreberi, Dicranum sp.).For most of the vertical measurement levels, the data were available since January 1996.But the highest measurement levels (101 and 125 m) were introduced in February 2013.These data altogether form a continuous dataset of atmospheric observations over a 20-year period.

Trace gas measurements
All the trace gas concentrations were recorded in every minute at all the aforementioned measurement levels.O 3 concentrations were determined by a photometric TEI 49C ozone monitor (Thermo Fisher Scientific, Waltham, MA, USA).The instrument had a detection limit of 1 ppb and a relative accuracy of ±3 %.
Carbon monoxide (CO) and carbon dioxide (CO 2 ) concentrations were quantified by their infrared (IR) absorption characteristics.CO was measured with a Horiba APMA 360 monitor (Horiba, Kyoto, Japan) with a detection limit of 20 ppb and relative accuracy of ±3 %.CO 2 was measured with an URAS 4 CO 2 monitor (Mannesmann Hartmann & Braun, Frankfurt am Main, Germany) that had a relative accuracy of ±1 % The CO 2 signal was corrected for H 2 O vapour interference.
Nitrogen oxide (NO x consists of NO and NO 2 ) concentrations were measured with a TEI 42CTL photometric analyser (Thermo Fisher Scientific, Waltham, MA, USA).The measurement had a 0.1 ppb limit of detection during 2006 measurements and this was increased to 0.15 ppb for the measurements conducted in 2009 due to modification of the instrument.The relative accuracy in all of the NO x measurements was ±10 %.

Ancillary measurements
The temperature was monitored from all the heights, up to 67.2 m, with custom-made 4-wired PT-100 sensors.The measurement spots were shielded in order to avoid interferences due to solar radiation and ventilated by fans to ensure a homogeneity in the measured air parcel.The accuracy of the temperature measurement was estimated to ±0.2 • C via a comparison with a mercury-based thermometer.
From April 1998 to mid-June 2012, the relative humidity (RH) values were obtained via measurement of dew point temperature and air temperature at 16.4 m on the sampling tower using a chilled mirror hygrometer (DewTrak Model 200M Meteorological Humidity System, EdgeTech, Marlborough, MA, USA).The measurement of the dew point and ambient temperatures were retrieved with an accuracy of ±0.3 • C given a RH accuracy of ±3 %.Since 16 June 2012, the relative humidity is directly measured with an RH sensor (Rotronic MP102H with Hygroclip HC2-S3, Rotronic AG, Bassersdorf, Switzerland), providing accuracy of ±0.8 %.
The total global radiation was calculated as the sum of direct and diffuse solar radiations.The measurement was placed on a separate tower next to the mast at the height of 18 m.The radiation was measured using pyranometers: a Reemann TP 3 pyranometer (Astrodata, Tõravere, Tartumaa, Estonia) from January 1996 to mid-June 2008, and a Middleton Solar SK08 First Class Pyranometer (Middleton Solar, Yarraville, Australia) since 15 June 2008.These devices have a detection limit of 5 W m −2 .The irradiance of solar radiation on a hemispherical surface (sun disc) was monitored in a wavelength range of 0.30-4.8µm, with an accuracy of ±5 %.
Wind speeds at different levels of the mast were measured either with cup anemometers or ultrasonic anemometers, or both, with an accuracy of ±0.1 m s −1 .Cup anemometers (A101M/L, Vector Instruments, Rhyl, Clwyd, UK) were used for monitoring wind properties at 4.2, 8.4, 16.8, 33.6, 50.4,and 74 m from January 1996 to 4 September 2003.Thereafter, the cup anemometers were replaced at all the measurement levels by ultrasonic anemometers (Ultrasonic anemometer 2-D, Adolf Thies GmbH, Göttingen, Germany).Since 7 July 2011, the wind measurement level at 74 m has been moved down to 67.4 m.Two 3-D ultrasonic anemometers (Metek uSonic-3, Elmshorn, Germany) were installed at 67.4 and 125 m in November 2015.
Friction velocity data were used to characterize the turbulent conditions.The friction velocity (u * ) was determined as 30 min averages of square root of the covariance between horizontal and vertical winds (Monteith and Unsworth, 2013), where u and w are the fluctuations in the horizontal and vertical winds, respectively.Horizontal and vertical winds were measured at 23.3 m with a fast-response acoustic anemometer (Solent HS1199) on the radiation tower.
For assessing the stability of the atmosphere, the potential temperatures were calculated using the following expression (e.g.Monteith and Unsworth, 2013): where T is the ambient temperature (in kelvin), P 0 is 1000 mbar, P is the ambient pressure (in mbar), R is the gas constant, and c p is the isobaric heat capacity.For air, the ratio gas constant to isobaric heat capacity (R/c p ) typically has a value of 0.286.

Data processing and handling
With the exception of wind data, both gas and meteorological data were arithmetically averaged into a 10 min time resolution.The wind data were treated as vectors, where the arithmetic averaging was applied to the two horizontal wind components of the decomposed vectors and the 10 min wind data were restored from these averaged components.
The duration of consecutive low ozone concentrations was determined based on the 10 min data for the concentrations below threshold values of 10, 5, 2, and 1 ppb.Occasionally, the data were available only for every second 10 min interval, mainly between January 1996 and February 1997 for all the heights and at heights of 101 and 125 m from February 2013 onwards.A gap-filling procedure was applied to such data in order to estimate the duration of a low-ozone episode: if ozone concentrations were found below a threshold in every second 10 min interval with the value in between being missing, then the missing value in between two low-concentration records was treated as a low-concentration point.Hereafter, the low-ozone episodes are referred to as ozone concentration drops, and when the ozone concentration remains below a given threshold (10, 5, 2, or 1 ppb) for more than 30 min, it is considered an ozone depletion event.No screening according to snow on the ground was taken into account in our data analysis.

Diurnal and seasonal variations in ozone profiles
A clear seasonality was seen in ozone concentrations at different heights (Fig. 1).The ozone concentrations were the highest in spring (March-April) at all heights and then declined towards the lowest values seen in early winter (November).This observation is similar to the reported feature in surface ozone concentrations at several northern midlatitude sites (Oltmans et al., 2006;Parrish et al., 2013).Atlas et al. (2003) and Vingarzan (2004) pointed out that the spring peak is probably caused by spring recovery of photochemical production (Dibb et al., 2003) and by ozone accumulation through winter (Liu et al., 1987) at the northern mid-to high latitudes, but hardly by transportation from the stratosphere, even though the intrusion of ozone from the stratosphere has been reported to be a major source of springtime ozone appearing in the middle troposphere at high northern latitudes (Dibb et al., 2003).Notably, starting from early summer and continuing until late autumn, ozone concentrations tended to increase with an increasing measurement height, being clearly the lowest close to the surface (Fig. 1).This feature is probably linked to some combination of the biotic ozone uptake by vegetation, ozone deposition on surfaces inside the forest and ozone loss through chemical reactions with forest emissions (Rannik et al., 2012;Monks et al., 2015;Clifton et al., 2017).
The diurnal behaviour of the vertical ozone concentration profile was similar in spring (March-May) and summer (June-August) (Fig. 2): a high ozone concentration was seen throughout the vertical extent of the measurements between about 10:00 and 20:00, likely related to the photochemical production of ozone.The highest ozone concentrations were generally seen above the canopy level (at about 50-60 m) during the daytime, with a clear decline with the decreasing height.This observation can be attributed to the loss of ozone by the forest (Rannik et al., 2012;Monks et al., 2015).During the evening in spring and summer, low ozone concentrations were first (at around 20:00) seen close to the ground and then later on gradually at levels higher up, so that the diurnal minimum ozone concentration was reached between about 05:00 (4.2 m) and 08:00 (125 m).This phenomenon results from the ozone deposition at night, when the photochemical production of ozone is prohibited.Ozone can also be consumed through oxidation reactions with biogenic organic compounds that are accumulated in the nocturnal boundary layer (Hakola et al., 2012).After sunrise, turbulent mixing in the boundary layer is initiated, leading to an increase in the mixed layer height and entrainment of ozone into this layer from above, in addition to which ozone starts to form via photochemical reactions.As a consequence, the ozone concentrations at all levels recovered rapidly, with a tendency of the higher-level concentrations to recover prior to those at lower levels.
The autumn ozone concentration pattern followed that seen in spring and summer, but with a slower ozone concentration recovery in the morning (Fig. 2).Also, the onset of the concentration recovery occurred slightly later in autumn (at 09:00) than in spring and summer (between 06:00 and 08:00).These features are very likely related to the shorter daytime length in autumn compared with spring and summer (Fig. 3).The later sunrise time in autumn facilitated a prolonged deposition-dominant ozone budget and a delayed replenishment.The reduced solar radiation intensity in autumn was unable to sustain a rapid enrichment of ozone from the residual layer.Meanwhile, the photochemical production of ozone was also weakened.As a consequence, the autumn ozone profile showed a gradual concentration recovery from high to low measurement levels in the morning.Moreover, the occurrence of the ozone concentration declination in the afternoon in autumn at around 17:00, well before that in spring and summer, can be possibly ascribed to the earlier sunset time and corresponding increase in the surface layer stability.Finally, the highest ozone concentrations tended to shift upwards in heights from spring to autumn.No clear diurnal pattern in ozone concentrations was seen in winter.

Frequencies of low ozone concentration episodes
As described before, we defined four low-concentration thresholds (10, 5, 2, and 1 ppb) to characterize low ozone concentration episodes.When such an episode lasted for more than 30 min, it is considered as an ozone depletion event, otherwise as an ozone drop.To assist the characterization of long-lasting depletion events longer than 3 h, a gapfilling procedure was utilized as described in Sect.2.3.
The observation frequency, i.e. the number of observations, of ozone concentration drops below 10 ppb exhibited similar diurnal variations at the different measurement levels (Fig. 4).The most pronounced diurnal variation was seen inside the canopy layer (≤ 16.8 m), with low ozone concentrations frequently observed at 04:00-08:00 (Fig. 4a, b and c).Low ozone observations were rare around noon.The diurnal behaviour of ozone concentration drops below 5 ppb resembled those of the below 10 ppb drops.However, excluding the 4.2 m measurement level, ozone concentrations below 5 ppb occurred also often in the afternoon at around 16:00-20:00.For the occurrence of ozone concentration drops below 2 ppb, in addition to the expected noon low, the diur-  nal distribution exhibited a second trough at around 02:00 in the early morning prior to the morning high, again except at 4.2 m.Interestingly at 4.2 m, the moment at which the morning peak observations occurred shifted in time from 04:00 to 06:00 for the sequence of ozone concentration drops below 10, 5, and 2 ppb.This feature was also perceptible at 8.4 and 16.8 m (Fig. 4b and c).No systematic diurnal variations were recognized for the observation frequency of ozone concentration drops below 1 ppb.
Ozone depletion events that lasted for more than 30 min occurred mainly in autumn and winter (Fig. 5).March had the lowest number of ozone depletion events below 10 ppb, after which the number of these events tended to increase over the course of both months (towards December) and measurement heights (from low to high).Except in May and June, ozone depletion events below 10 ppb were observed to reach the measurement level of 67.2 m.Ozone depletion events below 5 ppb were mainly observed below the canopy level (≤ 16.8 m) during June-September, while only during  October-January and in March did they reach 67.2 m.In February, April, and May, no ozone depletion events below 5 ppb were observed at any of the measurement levels.The below-2 ppb ozone depletion events exhibited a similar feature to the below-5 ppb events -however, with the depletion seen only at the lowest measurement level (i.e.4.2 m) in January and summer months.No depletion events below 2 ppb were observed during February-May.The below-1 ppb depletion events were concentrated in the period of October-December.
We found a number of long-lasting (> 3 h) ozone depletion events (Table 1).During the period 1996-2016, there were four such events, in which the ozone depletion below 2 ppb was observed up to the height of 67.2 m.This number increased to 26 in the case of below-5 ppb depletion events and to 161 (or 148 without the gap-filling procedure) in the case of below-10 ppb events.Long-lasting ozone depletion events below 10 ppb occurred all year around, whereas no such events below 5 ppb were observed in February, April, May, or June (Fig. 6).Long-lasting ozone depletion events below 2 ppb were found in October, November, and Decem-Table 1.The number of ozone depletion events seen at different heights that last for more than 3 h.Ozone depletion events were classified as concentration drops below 10, 5, 2, and 1 ppb.The statistics were based on the original 10 min data.Statistics based on gap-filled data are presented in parentheses.The gap-filling procedure resulted in no significant changes in the statistics of depletion events for heights of .Diurnal patterns in median ozone concentrations at 4.2 m in different months, from all the data (black), from below-10 ppb depletion events (red), from below-5 ppb depletion events (orange), from below-2 ppb (purple) and from below-1 ppb (green).A minimum duration threshold of 3 h was used to select days with ozone depletion events, based on the gap-filled data.
ber, whereas below-1 ppb events occurred only in October and November and they were generally associated with high humidity conditions.The variations of median diurnal ozone concentrations for different months derived from all 10 min data ranged from 17 to 45 ppb, with diurnal cycles seen in almost all months, except in January and December (black lines in Fig. 6).The median ozone concentrations obtained from the days with observations of long-lasting ozone depletion events were, in general, lower than those derived from all 10 min data, and displayed similar diurnal patterns to those derived from all 10 min data during May-October.No systematic diurnal behaviour was identified for any type of ozone depletion event between November and April, nor was it found for the below-2 ppb and the below-1 ppb depletion events in October.The difference between the median ozone concentrations determined on all the days and on depletion event days was about 10-15 ppb between October and February.For July-September, this concentration difference was big in the early morning, then diminished towards noon, and enlarged again after 16:00, as a result of boundary layer dynamics.

Features in low-ozone episodes
Ozone concentrations were, on average, lower at higher relative humidities (Fig. 7a), and low-ozone episodes were only www.atmos-chem-phys.net/18/49/2018/Atmos.Chem.Phys., 18, 49-63, 2018 observed in high RH (> 70 %) conditions.According to Altimir et al. (2006), moisture in the air enhances the total ozone deposition, possibly due to the high solubility of ozone in water (Sotelo et al., 1989).At RH levels higher than 70 % almost all surfaces of foliage are covered with thin water film, which efficiently adsorbs ozone and decouples the stomatal conductance from O 3 fluxes, forming a significant non-stomatal sink for ozone (Altimir et al. 2006).This can be responsible for about 40-60 % of the observed decrease in O 3 concentrations in the very humid cases in our data, and we can conclude that moisture content in the air is likely to be one of the key factors controlling these low-ozone episodes.
Low ozone concentrations occurred mostly during the nighttime, associated with little global radiation (Fig. 7b) and low wind speeds (Fig. 7d).Air temperatures have been reported to have a minor effect on the total ozone deposition (Rannik et al., 2012).Although we observed no direct relationship between air temperature and ozone concentration drops (Fig. 7c), the temperature might indirectly affect ozone deposition by influencing e.g. the RH.When ozone concentration dropped lower than 1 or 2 ppb at 67.2 m, there was typically a very small ozone concentration difference between the 67.2 and 4.2 m heights, i.e. no gradient in the vertical ozone concentration profile.These episodes were found in association with friction velocities between 0 and 1 m s −1 , which indicates relatively stable con-ditions with low vertical mixing.Small ozone gradients were mainly observed at low friction velocities (Fig. 8a).This is consistent with the observation, depicted in Fig. 7b, that most of the low ozone concentrations occurred under relatively dark conditions with low solar radiation.Ozone concentration below 5 or 10 ppb observed at 67.2 m occurred mostly when the potential temperature difference between the 67.2 and 4.2 m heights was around 0 • C (Fig. 8b), indicating neutral atmospheric conditions.Under such conditions, the atmosphere does not favour large-scale vertical mixing of ozone, yet small-scale turbulent eddies could be formed, aiding the deposition of ozone and leading to the occurrence of ozone concentration drops.
In addition, there existed a number of cases in which the 4.2 m ozone concentration was low while the concentration difference between 67.2 and 4.2 m was large (Fig. 9a) -i.e. a gradient in the vertical ozone concentration profile.These episodes were also observed at very low friction velocities (Fig. 8c), but they tended to exhibit larger concentration differences at larger differences in the potential temperature between 67.2 and 4.2 m (Fig. 8d).Such observations likely indicate that the gradients in the vertical ozone concentration profile occurred under calm and stably stratified conditions.The narrow range of the low friction velocities (0-0.5 m s −1 ) observed during these cases may suggest a nearly invariant deposition sink for ozone.In such circumstances, the smaller the mixing height, the more effective the loss of surface ozone that can be expected, due to the lack of ozone replenishment from higher up, which possibly results in the occurrence of the large number of ozone depletion events only at lower heights (over 200 below-10 ppb events at the three lowest heights, Table 1).
Based on the above-described features in the vertical ozone concentration profile, two major types of ozone depletion events were identified: (1) the ozone depletion occurred at low measurement levels only with a clear gradient in the vertical ozone profile, termed here as near-surface depletion, and (2) the ozone depletion extended to 67.2 m with a negligible ozone concentration difference between 67.2 and 4.2 m, termed here as full-profile depletion (Fig. 9a, c and d).
A possible third type lies in between the near-surface and full-profile depletion, in which decreases in ozone concentration could be seen at high measurement levels but with a remarkable difference in the ozone concentration between 67.2 and 4.2 m, and a clear gradient can be seen through the vertical concentration profile.This third type is referred to as gradient-profile ozone depletion.However, the boundary between the gradient-profile and the near-surface depletion can often be vague in ambient data, i.e. a depletion event classified as a near-surface event may be re-categorized as the gradient-profile type if a longer period is considered when a decrease trend in high-level ozone concentrations can be seen.During the years 1996-2016, over a hundred fullprofile below-10 ppb depletion events were observed, with 26 below-5 ppb events, four below-2 ppb events, and even one below-1 ppb event (Table 2).We define that a near-surface depletion event should fulfil two criteria: (1) the first quartile of the difference in the ozone concentration between 67.2 m and the height that the depletion reaches has a value larger than 2 ppb, and (2) the third quartile of the ozone concentration measured at 67.2 m during the depletion event is larger than 10 ppb.Then seven below-10 ppb events lasting longer than 3 h with depletion reaching 33.6 m were identified during our study period (Table 2).
The near-surface and the gradient-profile ozone depletion events were observed throughout the year, except in spring, with more pronounced ozone gradients typically seen in later summer and autumn (Fig. 9b).However, the full-profile de- a Full-profile ozone depletion events here are determined using the third quartile of the difference in the ozone concentration between 67.2 and 4.2 m, which has a value smaller than 2 ppb.b Near-surface ozone depletion events here are determined using the first quartile of the difference in the ozone concentration between 67.2 m and the height that the depletion reaches, which has a value larger than 2 ppb.Meanwhile, the third quartile of the ozone concentration measured at 67.2 m during the depletion event is larger than 10 ppb.pletion events were primarily observed in winter months (Fig. 9b).The near-surface depletion events, though often observed at night and early morning hours under relatively dark conditions, were not exclusively restricted in this time window.Occasionally, they also occurred during daytime hours -presumably when a stable condition prevails resulting in a low mixing height.An analysis of atmospheric conditions was carried out for a near-surface depletion event encountered during the night of 14-15 September 2009 (Fig. 9c), and for a full-profile de-pletion event that took place on the morning of 23 October 2006 (Fig. 9d).The near-surface depletion event was initiated before the sunset on 14 September and lasted over the course of the night.However, during the 2-3 h prior to the sunrise, unlike the ozone concentrations at 4.2 and 8.4 m showing a continuous decrease, the ozone concentration at 16.8 m (about the forest canopy height) tended to level out.This observation is likely ascribed to the formation of two decoupling layers above and below the forest canopy.Upon sunrise, the ozone concentrations at 4.2 and 8.4 m started shortly to recover from their lowest levels, while an instantaneous concentration recovery was possibly seen at 16.8 m, accompanied by a decrease in the 33.6 m ozone concentration.These phenomena indicate a rapid development of turbulent mixing vertically above the canopy level.But the break-up of the decoupling and the extension of the mixing height to 67.2 m came with some time lags.Interestingly, after the ozone concentration levels at different heights converged at around 09:00, there seems to be a transient obstruction in the further increase of ozone concentrations below 16.8 m.The ozone concentrations at 16.8 m and higher nearly levelled out until 10:00, while the 8.4 and 4.2 m concentrations showed a drop before a concentration increase was perceptible.Such a feature (in some cases, even a clear reduction in the ozone concentrations measured at all heights) can be seen in more than 30 % of the below-10 ppb near-surface depletion events at 4.2 m.These phenomena might be attributed to the ozone loss via photolysis and other consumptions through chemical reactions.In contrast, the full-profile event on 23 October 2006 occurred before sunrise and lasted through the morning hours, even until almost 14:00 at 4.2 m.The recovery of the ozone concentrations was first seen at 67.2 m before noon and then gradually at lower heights.
Both near-surface and full-profile depletion events displayed distinct features in relation to the air temperature (Fig. 10a-l).While both depletion events were observed at low-temperature conditions, a temperature gradient was only present during the near-surface depletion event, with higher temperatures found at the highest measurement levels.This observation shows that this near-surface depletion event occurred during a temperature inversion.A similar feature was identifiable during another depletion episode on the early morning of 14 September 2010 (Fig. 10a), with very low wind speeds at the low measurement levels and mild turbulent conditions at the high measurement levels (Fig. 10c).In contrast, calm conditions were present during the full-profile ozone depletion event observed on October 2006 (Fig. 10i).
Humid conditions have been previously reported to favour the occurrence of low ozone concentrations (Altimir et al., 2006).During both near-surface and full-profile ozone depletion events, the relative humidity was high (Fig. 10b, h).The results obtained here thus indicate that the ozone depletion at this relatively clean boreal forest site tends to occur under low-mixing calm conditions with a high moisture content in the air (as can be seen in Fig. 10c, d, i, j).Similar low-mixing conditions have been observed previously in association with polar boundary layer ozone depletion events (Simpson et al., 2007;Cao et al., 2016), but no relation to RH, as found in our observations, was reported in the polar region (note that polar ozone depletion is characterized by halogen chemistry, which has not been observed at our boreal forest site).
In addition, it seems that an elevated CO concentration is a potential indicator of a full-profile depletion (Fig. 10k), since low ozone concentrations were found in association with high CO concentrations.The recovery of ozone concentrations from the depletion event was observed first at high measurement levels, and at about the same time the CO concentration began to decrease.However, no clear trace of such a feature was seen in the near-surface depletion case, though the highest CO concentrations seemed to coincide with the lowest ozone concentrations (Fig. 10e).The relation between the ozone depletion and NO x concentration was also inspected (Fig. 10f, l), but no systematic connection that could explain our observations was found.

Conclusions
Here, we have investigated vertical profiles of ozone concentrations inside and above a conifer stand over a period of 20 years  in a boreal forest environment in southern Finland (the SMEAR II station).The variations in ozone concentrations were investigated in relation to other measured trace gases and meteorological variables.Our results show that, in general, ozone concentrations increased with an increasing height above the ground.The daily median ozone concentrations at our measurement site were found to be in the range of about 20-50 ppb, with clear maxima in the springtime.
Temporal significant decreases of ozone concentrations at the ground level have been observed in very polluted urban environments (US National Research Council, 1992;Gregg et al., 2003).However, to our knowledge and based on observations, the occurrence of very low ozone concentrations in relatively clean environments other than the polar regions has not been reported previously.The main aim of our study was to establish the frequency and strength of ozone depletion events at our boreal forest site.By analysing 20 years of ozone measurements, we identified on average more than a thousand of 10 min-periods at 4.2 m, during which ozone concentrations were below 10 ppb, and a few tens of cases of ozone concentration below 2 ppb.These low ozone concentration episodes occurred typically at low temperature, low wind speed, high RH, and limited intensity of solar radiation, which are conditions likely related to the formation of a shallow mixing layer.When such conditions prevail, neither photochemical production nor enrichment from the residual layer can provide sufficient replenishment of ozone into the confined mixing layer, while ozone loss through deposition and chemical consumption continued unperturbed.
Future investigations should relate ozone depletion to ozone chemistry and boundary layer dynamics in order to understand the underlying mechanisms involved in these particular events.In addition, it is worth exploring whether there is a connection between ozone depletion and atmospheric new particle formation.It seems that several unrevealed or at least rarely observed processes are likely to take place during ozone depletion events in clean to moderately polluted environments.Data availability.The data used in this work were collected at Hyytiälä SMEAR II station, which are available in the SmartSMEAR database (https://avaa.tdata.fi/web/smart/smear/).
Competing interests.The authors declare that they have no conflict of interest.Special issue statement.This article is part of the special issue "Pan-Eurasian Experiment (PEEX)".It does not belong to a conference.

Figure 1 .
Figure 1.Daily medians of ozone concentrations at different heights for years 1996-2016.Heights are indicated by different colours.

Figure 2 .
Figure 2. Diurnal profiles of ozone concentrations in different seasons.Spring: March-May, summer: June-August, autumn: September-November, and winter: December-February.Panel (a): ozone concentrations in different scales for different seasons to emphasize the profile information.Note the significantly larger variability of O 3 concentrations during summer months.Panel (b): ozone concentrations in the same scale for different seasons.

Figure 4 .
Figure 4. Hourly distribution of observations of ozone concentration drops below 10, 5, 2, and 1 ppb at different heights: (a) 4.2 m, (b) 8.4 m, (c) 16.8 m, (d) 33.6 m, (e) 50.4 m, and (f) 67.2 m.The statistics were based on the 10 min ozone data and the observation frequency is presented in the logarithmic scale.

Figure 5 .
Figure 5. Vertical extension of ozone depletion events that last longer than 30 min.The numbers of depletion event observations are colour coded in logarithmic scale for ozone concentration drops below 10, 5, 2, and 1 ppb.

Figure 7 .
Figure 7. Ozone concentrations vs. (a) relative humidity (RH), (b) global radiation, (c) temperature (T ), and (d) wind speed.Plots are made with 10 min data collected from the 16.8 m measurement level during the period 1996-2016, with an exception of the global radiation data from 18 m.

Figure 8 .
Figure 8. Ozone concentration difference between 67.2 and 4.2 m measurement levels vs. (a, c) friction velocity determined at 23.3 m and (b, d) potential temperature difference between 67.2 and 4.2 m measurement levels.The colour codes in (a, b) represent concentration drops below 10, 5, 2, and 1 ppb from the ozone measurement at 67.2 m.The green crosses in (c, d) mark the data points fulfilling the conditions that the ozone concentration at 4.2 m is smaller than 2 ppb, while the difference in the ozone concentration between 67.2 and 4.2 m is larger than 2 ppb.

Figure 9 .
Figure 9.The ozone concentration difference between 67.2 and 4.2 m vs. the ozone concentration (≤ 2 ppb) at 4.2 m, (a) colour-coded by ozone concentrations at 67.2 m and (b) colour-coded by month.(c) A near-surface depletion event on 14-15 September, 2009.(d) A full-profile depletion event on 23 October 2006.The vertical black lines indicate the sunrise (SR) and sunset (SS) times.