Vortex-wide chlorine activation by a mesoscale PSC event in the Arctic winter of 2009/10

In the Arctic polar vortex of the 2009/10 winter temperatures were low enough to al-low widespread formation of Polar Stratospheric Clouds (PSC). These clouds occurred during the initial chlorine activation phase which provided the opportunity to investigate the impact of PSCs on chlorine activation. Satellite observations of gas-phase species 5 and PSCs are used in combination with trajectory modeling to assess this initial activation. The initial activation occurred in association with the formation of PSCs over the east coast of Greenland at the beginning of January 2010. Although this area of PSCs covered only a small portion of the vortex, it was responsible for almost the entire initial activation of chlorine vortex wide. Observations show HCl mixing ratios decreased 10 rapidly in and downstream of this region. Trajectory calculations and simpliﬁed heterogeneous chemistry modeling conﬁrmed that the initial chlorine activation continued until ClONO 2 was completely depleted and the activated air masses were advected throughout the polar vortex. For the calculation of heterogeneous reaction rates, surface area density is estimated from backscatter observations. Modeled heterogeneous 15 reaction rates along trajectories intersecting with the PSC indicate that the initial phase of chlorine activation occurred in just a few hours. These calculations also indicate that chlorine activation on the binary background aerosol is signiﬁcantly slower than on the PSCs and the observed chlorine activation can only be explained by an increase in surface area density due to PSCs. Furthermore, there is a strong correlation between the the depletion hand, showed that nearly complete chlorine activation can be achieved in a mountain-wave PSC at su ﬃ ciently low temperatures and large surface area densities. This study investigates the inﬂuence of mesoscale PSCs on the chemical compo- 20 sition of the entire vortex. Mesoscale PSCs are larger in spatial scale than mountain-wave induced PSCs but still only cover a small fraction of the polar vortex. We use data from the space-borne lidar CALIOP (Cloud-Aerosol LIdar with Orthogonal Polarization) instrument for studying PSCs and data from the Microwave Limb Sounder (MLS) in combination with chemistry transport model calculations to examine the im- 25 pact of a mesoscale PSC on the chemical composition of the Arctic vortex in January 2010. CALIOP backscatter observations are used to derive particle surface area density (SAD) which is then used to calculate chlorine activation. Through these anal-33733


Introduction
Heterogeneous chemistry on stratospheric aerosol and polar stratospheric clouds (PSCs) plays a crucial role in the formation of the Antarctic ozone hole (Solomon, 1999).While the stratospheric aerosol layer is present globally at all times (Junge et al., 1961), PSCs only form over the polar regions of the winter hemisphere (Lowe and MacKenzie, 2008).PSCs are ubiquitous in the Antarctic winter but wide-spread PSC occurrence over the Arctic is only observed during very cold winters (Pitts et al., 2009(Pitts et al., , 2011)).PSCs have a substantial influence on the chemical composition of the polar stratosphere: on the one hand heterogeneous chemistry on PSC particles impacts the partitioning of inorganic chlorine by chlorine activation (Solomon et al., 1986), i.e., by converting inorganic chlorine reservoir species (mainly HCl and ClONO 2 ) to photo-active species (ClO x = ClO + 2 • Cl 2 O 2 ).On the other hand, PSC particles can grow large enough to effectively denitrify (Fahey et al., 2001) and dehydrate (Kelly et al., 1989) the lower stratosphere.Although PSCs have a pivotal role in determining stratospheric chemical composition and have been subject to extensive research since the 1980s, some details are still subject to uncertainty.
Two of the key questions are the following: (1) to which extent is heterogeneous chemistry on PSC particles responsible for chlorine activation and (2) what are the timescales for this processing?Several studies (e.g.Drdla and Müller, 2012;Wohltmann et al., 2013;Kirner et al., 2015) have suggested that the influence of PSCs on chlorine activation in subordinate to that of cold binary aerosol on a vortex-wide scale.When temperatures approach the frost point, hetero-Published by Copernicus Publications on behalf of the European Geosciences Union.
geneous reaction rates are large enough to activate chlorine on timescales of minutes to hours regardless of surface type (Wegner et al., 2012).On the other hand, Carslaw et al. (1998) and Kühl et al. (2004) showed that nearly complete chlorine activation can be achieved in a mountain-wave PSC at sufficiently low temperatures and large surface area densities.
This study investigates the influence of mesoscale PSCs on the chemical composition of the entire vortex.Mesoscale PSCs are larger in spatial scale than mountain-wave induced PSCs but still only cover a small fraction of the polar vortex.We use data from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP, Winker et al., 2009) instrument for studying PSCs and data from the Microwave Limb Sounder (MLS, Waters et al., 2006) in combination with model calculations of the Chemical Lagrangian Model of the Stratosphere (CLaMS, McKenna et al., 2002) to examine the impact of a mesoscale PSC on the chemical composition of the Arctic vortex in January 2010.CALIOP backscatter observations are used to derive particle surface area density (SAD) which is then used to calculate chlorine activation.Through these analyses, we will examine the impact of PSC SAD enhancements on chlorine activation compared with the cold binary background aerosol.

Instrumentation and model description
The Cloud-Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO) satellite is member of the A-train satellite constellation (Stephens et al., 2002) orbiting at an inclination of 98.2 • that provides coverage up to 82 • latitude in each hemisphere.Its primary instrument is CALIOP which measures backscatter at 1064 and 532 nm, with the 532 nm channel separated into orthogonal polarization components parallel and perpendicular to the polarization plane of the outgoing laser beam.The CALIOP classification scheme distinguishes between supercooled ternary solutions (STS), mixtures of STS and nitric acid trihydrate (NAT) and ice.Mixtures of STS and NAT are further divided into three groups: MIX 1, MIX 2 and MIX 2 enhanced, in order of increasing NAT number density.CALIOP also distinguishes between synoptic and wave ice.PSCs are identified in the CALIOP measurements as enhancements in backscatter at 532 nm with an altitude-dependent threshold of β 532 between 2 and 4 × 10 −5 km −1 sr −1 in the algorithm described in Pitts et al. (2013).We show sensitivity studies that use different threshold backscatter ratios between 3 × 10 −5 and 1.5 × 10 −4 km −1 sr −1 , where the latter backscatter threshold is high enough to indicate the presence of ice PSCs.Details about the CALIOP PSC classification algorithm can be found in Pitts et al. (2009Pitts et al. ( , 2011Pitts et al. ( , 2013)).
The microwave limb sounder (MLS) is an instrument on the Aura satellite, which is also part of the A-train constellation and has provided nearly continuous measurements since 2004.MLS provides about 3500 profiles of gas-phase species and temperature from Earth's surface to 90 km altitude between 82 • N and 82 • S per day.In this study, we use MLS HCl data from retrieval version 3.3 (Livesey et al., 2011) which are comparable to version 2.2.Retrieval version 2.2 for HCl has been validated in Santee et al. (2008).The vertical resolution and precision of the HCl observations are 3 km and 0.2-0.3ppbv, respectively.
CALIOP and MLS observations are linearly interpolated to a 20 • by 2 • (longitude × latitude) grid on fixed potential temperature surfaces each day.Potential temperature is calculated from the original MLS pressure levels and temperature data from the modern-era retrospective analysis for research and applications (MERRA, Rienecker et al., 2011).
Air mass trajectories are calculated with CLaMS which uses a fourth-order Runge-Kutta scheme.The wind fields for these trajectories are taken from ERA-INTERIM (Dee et al., 2011).Threshold temperatures T NAT (NAT existence temperature) and T ACl (chlorine activation temperature) which can serve as indicators for chlorine activation (Drdla and Müller, 2012) are calculated with temperatures and trace gas concentration from a CLaMS simulation described in Grooß et al. (2014).The vortex edge is defined according to Nash et al. (1996).
A simple algorithm was defined to determine the locations of air masses that had previously passed through the PSC event.Trajectories are used to track air masses downwind of a PSC that is defined by the CALIOP detection threshold.These air masses are called "Processed Air" since we assume that chlorine activation predominantly occurred in such air masses.On a given day a grid box is labeled as filled with "Processed Air" if a processed air mass has passed through this grid box on that day.A grid box is then labeled as filled with "Unprocessed Air" if on a given day no processed air masses have passed through this grid box even if this grid box was marked with "Processed Air" on previous days.This way we can describe for each day the fraction of the vortex that is filled with air masses which have previously encountered a PSC.
The determination of heterogeneous reaction rates requires an estimation of the particle surface area density (SAD) of PSCs and the background aerosol as realistic as possible.Following the approach of Gobbi (1995), we performed Mie calculations for unimodal lognormal particle size distributions representative of the polar stratosphere to derive a relationship between CALIOP measurements of particulate backscatter at 532 nm and the liquid particle SAD. Figure 1a shows calculated liquid particle SAD and backscatter (β 532,liquid ) as a function of temperature relative to the frost point (T − T frost ) for the following conditions: pressure = 30 hPa, HNO 3 = 15 ppbv, H 2 O = 5 ppmv, particle number density N = 10 cm −3 (Curtius et al., 2005), and lognormal geometric standard deviation σ = 1.6 (Lambert et al., 2012) which describes the width of the distribution.For these calculations, liquid particle volume was prescribed as a func-tion of temperature according to Carslaw et al. (1995); mode radius (for the Mie calculations) and SAD were then calculated from particle volume using standard relationships for a unimodal lognormal.
Calculated values of liquid SAD and β 532,liquid from Fig. 1a are plotted against each other in Fig. 1b.The dashed curve shows the least-squares quadratic fit in log-log space between the two parameters, which we will use in this study.With fixed number density N = 10 cm −3 and width of lognormal distribution σ = 1.6, calculated points fall along the same least-squares curve for different pressure levels (50 and 70 hPa), HNO 3 mixing ratios (2, 5, and 10 ppbv), and H 2 O mixing ratio (2 ppmv).Values calculated for different number densities (N = 5 and N = 15 cm −3 ) and geometric standard deviations (σ = 1.3 and σ = 1.8) fall along different curves, but all points lie within the limits depicted by the solid gray curves.Given the small envelope defined by those parameters the relationship between SAD and β 532,liquid is robust and covers a wide range of possible conditions in the polar vortex.Measured values of β 532,liquid larger than 1.5 × 10 −4 km −1 sr −1 likely signify ice PSCs.The SAD estimated using the least-squares relationship for liquid particles can be considered a lower limit of SAD for ice PSCs.The relationship between backscatter and surface area density for STS particles can then be described by and NAT particles.Figure 2 also shows the limitations of our current observations with gaps of hundreds of kilometers between the orbits.During this PSC event, the vortex was shifted away from the pole over northeast Canada and northern Russia.The vortex remained stable during this time frame with only limited dynamic disturbances.
Figure 3 shows backscatter over a 6-day period from 31 December 2009 to 5 January 2010 and the prevailing winds.From the meteorological data we estimate the wind speed in the vicinity of the PSC to be 20 m s −1 which yields a residence time of about 1 day for air parcels traveling through the PSC.Chlorine activation is limited by the availability of the reservoir species and will cease once either HCl or ClONO 2 are depleted (Portmann et al., 1996).Heterogeneous chemistry on stratospheric aerosols is primarily the Therefore, the observations of gas-phase HCl by MLS can serve as an indicator of heterogeneous processing with low values of HCl indicating air masses where inorganic chlorine has been activated.Since HCl is more abundant than ClONO 2 at the beginning of winter, ClONO 2 is the limiting component in chlorine activation (Portmann et al., 1996).
Figure 3a shows that low HCl mixing ratios are co-located with the PSCs over northern Greenland on 31 December and that HCl mixing ratios upstream of the PSC are more than twice as high compared to air exposed to the PSC.CALIOP observed backscatter values above 1.5 × 10 −4 km −1 sr −1 which corresponds to a tenfold increase of SAD relative to background conditions, or to around 10 µm 2 cm −3 .During all subsequent days, HCl is always lower in the vicinity of PSCs compared to areas upstream with HCl mixing ratios decreasing to less than 0.5 ppbv during this period, about a quarter of their values in a chemically unperturbed vortex.
The low HCl mixing ratios persist downstream of the PSCs and do not decrease further, in agreement with current knowledge that chlorine activation is highly dependent on SAD, temperature and the ratio of HCl / ClONO 2 .Model calculations (not shown) indicate that ClONO 2 is totally depleted downstream of the PSC.Additional activation can only occur when ClONO 2 is regenerated (Portmann et al., 1996).ClONO 2 controls the total amount of chlorine that can be activated while PSCs determine the area in which activation occurs.Low HCl concentrations observed in regions where our calculations indicate that only unprocessed air is present, could originate from either physical mixing with processed air or averaging along the MLS field of view.We also need to stress that such trajectory calculations rely on the quality of the wind fields used and are subject to errors which are difficult to quantify.
The air which was exposed to PSCs on 31 December 2009 (Fig. 3a, enclosed by red and blue contour lines) has completed a full circumnavigation of the vortex after 5 days (5 January 2010, Fig. 3f) and re-entered the region where it first encountered PSCs.Air which was not exposed to PSCs tends to have higher HCl mixing ratios, indicating that little chlorine activation occurred outside of PSCs, as shown clearly, for example, on 2 January in Fig. 3c.Air between Norway and Novaya Zemlya followed its own cyclonic circulation separate from the polar vortex.Air in this region did not encounter PSCs nor sunlight, leading to constant HCl mixing ratios.This situation persisted for several days and the trajectories passing through PSCs actually flowed around this region.After 3 January (Fig. 3d-e) air with low HCl mixing ratios flows into the PSC covered area from the unobserved area north of 82 • so we do not have a complete history of HCl concentrations in relation to PSC exposure for these air masses.
As Fig. 3 indicates, it takes about 5 days for air masses to circumnavigate the polar vortex.This means that a PSC with a sufficiently large meridional extent can activate chlorine throughout the entire vortex even when the PSC itself only covers a small fraction of the vortex.Such a PSC can serve as a "processing reactor" for chlorine activation with heterogeneous chemistry basically limited to the time air masses spend inside the PSC.Vortex-averaged HCl mixing ratios which are commonly used to estimate chlorine activation would not accurately represent the true nature of activation within the vortex, under these conditions.
While chlorine activation is generally observed within the boundary of elevated backscatter values as observed by CALIOP, the lowest HCl mixing ratios are correlated with the highest backscatter values.Figure 4 shows the ratio of mean HCl mixing ratios in processed air which has been ex- posed to backscatter values above certain thresholds to the vortex average HCl mixing ratio.The backscatter thresholds in Fig. 4a-c indicate the formation of STS. Figure 4d suggests the presence of ice since backscatter values of this magnitude correspond to higher surface area densities than physically possible with only the condensation of HNO 3 .
The ratio of mean HCl mixing ratios to the vortex average HCl mixing ratio is also shown for air which has been exposed to temperatures below either T NAT or T ACl .After 5 days the HCl mixing ratios within processed air are indistinguishable from the vortex mean (indicated by the white contour on the right side of each panel) because all air within the vortex has been exposed to PSCs.With smaller threshold backscatter values chosen for the area describing the "processing reactor" (Fig. 4d to a), HCl mixing ratios in air which has been exposed to PSCs gradually approach the vortex average.This shows that chlorine activation occurs locally in the area covered by PSCs but not on a vortex-wide scale.
Temperature thresholds for chlorine activation do not capture this localized effect.The air exposed to an area enclosed by a temperature threshold is always very similar to the vortex average since those areas already encompass the main part of the vortex.Therefore, neither T NAT nor T ACl provide realistic information on the location where chlorine activation occurs for the period considered in this study.Only thresholds based on backscatter provide this information.
In addition, Fig. 4 indicates that at the beginning of the activation phase under study, the average HCl mixing ratio in the area above a chosen backscatter threshold decreases as the backscatter threshold is increased.Since backscatter can be used as a proxy for surface area density, the surface area provided by PSCs appears to have a direct affect on gasphase HCl concentrations.While HCl mixing ratios above a backscatter threshold indicative of STS (Fig. 4a) have a minimum of 70 % of vortex average HCl mixing ratios, HCl mixing ratios above a threshold suggestive of ice PSCs (Fig. 4d) have a minimum of 50 % of the vortex average HCl mixing ratios.Therefore, chlorine activation processes visibly faster with increasing surface area density.The magnitude of chlorine activation appears to be correlated with available surface area density provided by PSCs.
The mesoscale PSC event leaves a visible mark on HCl mixing ratios.When separating the polar vortex into regions of processed and unprocessed air, we find lower HCl mixing ratios in processed air (Fig. 5) between 31 December and 9 January.Observations by MLS show that for all days between 20 and 25 km, HCl is lower in processed air compared to unprocessed air.However, the large uncertainties indicate that the distinction between processed and unprocessed air is difficult and that the mean values are calculated from heterogeneous air masses.Still, we see a clear trend that chlorine activation has progressed further in air masses which have been exposed to PSCs.
The analysis for Fig. 5 was repeated with T NAT and T ACl (not shown) as indicators for processed air.Here, processed air describes those grid boxes with temperatures less than the respective threshold temperatures and grid boxes where trajectories are present which have encountered temperatures less than the temperature thresholds in their past.Two main differences emerge.First, the daily mean mixing ratios for processed and unprocessed air are more similar, making a distinction between them difficult.Second, processed air covers the entire vortex more quickly because both temperature thresholds cover a larger area of the vortex than our PSC threshold.This leads to fewer available data points.Using either of the temperature thresholds to describe processed air  The ratio is calculated for six areas, four encompassed by backscatter thresholds (a-d) and two by temperature thresholds (T NAT , e and T ACl , f).A ratio of less than one indicates air inside the contour has a lower HCl mixing ratio than the vortex average.
does not show any significant difference between processed and unprocessed air.

Modeling heterogeneous chemistry
For a quantitative understanding of the shown HCl decrease we use the trajectories combined with a very simplified heterogeneous chemistry scheme.Using the relationship in Eq. ( 1), we can use CALIOP backscatter to calculate PSC SAD along trajectories passing through the maximum PSC backscatter.The trajectories are initialized with H 2 O, HNO 3 and HCl from MLS observations and it is assumed that at the start of the trajectories no chlorine has been activated yet.Hence, ClONO 2 is initialized as the difference between total inorganic chlorine (Cl y ) and HCl with Cl y derived from the Cl y -N 2 O tracer-tracer relationship (Grooß et al., 2002).
The simplified chemistry consists of three heterogeneous reactions which are modeled along the trajectories: No additional reactions are included in the calculations since they are not relevant to assess chlorine activation on timescales of up to 1 day.The uptake coefficients for all three reactions are calculated for STS even when backscatter values indicate the presence of ice.However, no difference in HCl mixing ratios is evident when the uptake coefficients for ice are used once backscatter values suggest the presence of ice PSCs.
Figure 6a-c shows backscatter, temperature, time below T NAT -3 K and time below T ICE for three different trajectories.These trajectories are initialized on 30 December at different locations so each of them would encounter their first maximum in backscatter on another day.Trajectory 1 (Fig. 6a) encounters its maximum on 31 December, Trajectory 2 (Fig. 6b) on 1 January and Trajectory 3 (Fig. 6c) on 2 January.Figure 6d-f then shows the calculated SAD, modeled HCl and ClONO 2 for two different cases and observations by MLS along these trajectories.
The maximum backscatter suggests the presence of ice clouds when the trajectories encounter their first maximum in backscatter; however, temperatures only decrease below the frost point on 6 January for trajectory 1 which encountered its backscatter maximum on 31 December (Fig. 6a).This suggests that temperatures below the frost point might not be resolved by the used meteorological data.The evolution of modeled HCl and observations (Fig. 6df) shows very good agreement for the first activation phase.Nakajima et al. (2016) also reported very good agreement of model and observations during this phase.HCl along the trajectories deviates from the observations after a couple of days.However, chlorine activation in the first half of the trajectories can still be analyzed.A decrease in HCl coincides with an increase in surface area density and ClONO 2 is totally depleted after the initial activation.This confirms that ClONO 2 is indeed the limiting factor for chlorine activation because observed HCl also does not decrease further once the model shows complete removal of ClONO 2 .The complete removal of ClONO 2 also occurs before the first maximum in surface area density (Fig. 6d-f), but activation starts when SAD increases above background values.Neither observations nor modeled chemistry indicate significant chlorine activation in the absence of PSCs.In fact, when heterogeneous chemistry is only modeled on a binary sulfate aerosol (without uptake of HNO 3 as temperatures decrease) chlorine activation does not only progress slower (Fig. 6d-f, dashed lines) but also does not completely deplete ClONO 2 .
While previous studies have shown that for an entire winter season heterogeneous chemistry on cold binary aerosol is sufficient to achieve complete chlorine activation (Drdla and Müller, 2012;Wohltmann et al., 2013), in the case of a mesoscale event as described in this study chlorine activation on PSCs is faster than on the background aerosol.The observed activation could not have occurred without the additional surface area provided by PSCs.However, Fig. 6f also shows that after 5 days, heterogeneous chemistry on only the binary aerosol has depleted all the initially available ClONO 2 .In agreement with Drdla and Müller (2012) the degree of chlorine activation on only the background aerosol would become indistinguishable from chlorine activation on PSCs after about 10 days, or the time it takes for an air parcel to fully circumnavigate the vortex twice.

Conclusions
We have analyzed CALIOP and MLS observations in combination with modeled trajectories to quantify the initial chlorine activation phase for the winter 2009/10 and constrain the spatial and temporal scales on which chlorine activation occurred; thus, this answers the question regarding the extent to which heterogeneous chemistry on PSC particles is responsible for chlorine activation and the timescales for this T. Wegner et al.: Chlorine activation on mesoscale PSCs processing.Our analysis has shown that mesoscale PSCs can have a substantial effect on chlorine chemistry throughout the polar vortex, even though the PSCs themselves only cover a small fraction of the polar vortex.A substantial decrease in HCl is observed in air masses exposed to PSCs for about 24 h.MLS observations indicate that air masses with low HCl mixing ratios occur downstream of this PSC event in agreement with the trajectory calculations.The modeled trajectories provide a solid approximation of the path the air has taken after it encountered the PSC and allow to distinguish between processed and unprocessed air masses.The average daily HCl mixing ratio shows substantially smaller HCl values in processed than in unprocessed air.We also show that chlorine activation in the polar vortex is not always a uniform process but can occur in mesoscale PSC "processing reactors".Heterogeneous chemistry occurs in these "processing reactors" and the air is subsequently advected and mixed throughout the vortex.Chlorine activation does not occur homogeneously throughout the vortex and HCl mixing ratios can vary significantly, especially during this initial activation phase.Therefore, a vortex average point of view does not provide an accurate representation during this phase.The trajectory calculations show that the availability of ClONO 2 limits the extent of chlorine activation.For the first time, CALIOP backscatter observations were utilized to estimate surface area density and model heterogeneous chemistry along trajectories.Results from these simulations are in good agreement with observed mixing ratios of HCl.These calculations also show that the SAD enhancements from PSCs lead to faster chlorine activation than would occur on the background aerosol.While the background aerosol could eventually activate the same amount of chlorine as the PSCs, over the time and spatial scales considered in this study, the observed rate of chlorine activation can only be explained by the additional surface area provided by PSCs.While this study focuses on the Arctic, similar conditions like the situation over Greenland in January 2010 can also occur over the Antarctic.Optically thick PSCs over Antarctica predominantly occur over the Antarctic peninsula; therefore, this area can also serve as a processing reactor for chlorine activation.
Figure 1.(a) Temperature dependence of STS surface area density (top x scale) and the corresponding theoretical backscatter signal at 532 nm for typical stratospheric conditions.(b) Relationship between STS surface area density and particulate backscatter color coded according to their theoretical temperature.Gray curves mark the deviation from the fitted function for varying ambient conditions.

Figure 3 .
Figure 3.The meteorological situation over a 6-day period from 31 December 2009 to 5 January 2010 for the Arctic winter 2009/2010 on the 510 K isentrope (approximately 22 km).Arrows show the prevailing wind direction and speed.The vortex edge is calculated according to Nash et al. (1996) and shown in solid black.Solid magenta indicates regions of the vortex with air which has passed through a PSC.Blue and red contours show the particulate backscatter from CALIOP where blue represents the detection limit for CALIOP (3 × 10 −5 km −1 sr −1 ) and red values are indicative of PSCs with a surface area density of more than 10 times the background values (1.5 × 10 −4 km −1 sr −1 ).Gray contours show HCl mixing ratios from MLS.

Figure 4 .
Figure 4. Ratio of average HCl mixing ratio of air that passed through a specified area to the vortex average HCl mixing ratio.The ratio is calculated for six areas, four encompassed by backscatter thresholds (a-d) and two by temperature thresholds (T NAT , e and T ACl , f).A ratio of less than one indicates air inside the contour has a lower HCl mixing ratio than the vortex average.

Figure 5 .
Figure 5. Correlation of daily average MLS HCl in processed (Backscatter threshold: 3 × 10 −5 km −1 sr −1 ) and unprocessed air from 31 December to 9 January for six different altitudes.Error bars indicate the 1σ standard deviation of calculating the daily average HCl mixing ratios.

Figure 6 .
Figure 6.Temporal evolution of three trajectories on 510 K potential temperature which were initialized at the same time at different locations.The trajectories encounter their first maximum in backscatter on 31 December (a, d), 1 January (b, e) and 2 January (c, f), respectively.The period from 30 December 2009 to 8 January 2010 is shown for all three trajectories.Left column (a-c) shows backscatter (blue), temperature (red), time below T NAT -3 K (red shading) and time below T ICE (blue shading).Right column (d-f) shows surface area density (blue), HCl (red) and ClONO 2 (green).Red stars are the daily mean HCl mixing ratios observed by MLS interpolated on the position of the trajectory and circles are CLaMS mixing ratios of HCl (red) and ClONO 2 (green), respectively.Solid red and green lines represent HCl and ClONO 2 where SAD for calculating heterogeneous reaction rates is calculated from the observed backscatter.Dashed red and green lines represent HCl and ClONO 2 where SAD is calculated for a binary aerosol.