Articles | Volume 13, issue 12
https://doi.org/10.5194/amt-13-6837-2020
https://doi.org/10.5194/amt-13-6837-2020
Research article
 | 
16 Dec 2020
Research article |  | 16 Dec 2020

Validation of SMILES HCl profiles over a wide range from the stratosphere to the lower thermosphere

Seidai Nara, Tomohiro O. Sato, Takayoshi Yamada, Tamaki Fujinawa, Kota Kuribayashi, Takeshi Manabe, Lucien Froidevaux, Nathaniel J. Livesey, Kaley A. Walker, Jian Xu, Franz Schreier, Yvan J. Orsolini, Varavut Limpasuvan, Nario Kuno, and Yasuko Kasai
Abstract

Hydrogen chloride (HCl) is the most abundant (more than 95 %) among inorganic chlorine compounds Cly in the upper stratosphere. The HCl molecule is observed to obtain long-term quantitative estimations of the total budget of the stratospheric chlorine compounds. In this study, we provided HCl vertical profiles at altitudes of 16–100 km using the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) from space. The HCl vertical profile from the upper troposphere to the lower thermosphere is reported for the first time from SMILES observations; the data quality is quantified by comparison with other measurements and via theoretical error analysis. We used the SMILES level-2 research product version 3.0.0. The period of the SMILES HCl observation was from 12 October 2009 to 21 April 2010, and the latitude coverage was 40 S–65 N. The average HCl vertical profile showed an increase with altitude up to the stratopause ( 45 km), approximately constant values between the stratopause and the upper mesosphere ( 80 km), and a decrease from the mesopause to the lower thermosphere ( 100 km). This behavior was observed in all latitude regions and reproduced by the Whole Atmosphere Community Climate Model in the specified dynamics configuration (SD-WACCM). We compared the SMILES HCl vertical profiles in the stratosphere and lower mesosphere with HCl profiles from Microwave Limb Sounder (MLS) on the Aura satellite, as well as from the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) on SCISAT and the TErahertz and submillimeter LImb Sounder (TELIS) (balloon borne). The TELIS observations were performed using the superconductive limb emission technique, as used by SMILES. The globally averaged vertical HCl profiles of SMILES agreed well with those of MLS and ACE-FTS within 0.25 and 0.2 ppbv between 20 and 40 km (within 10 % between 30 and 40 km; there is a larger discrepancy below 30 km), respectively. The SMILES HCl concentration was smaller than those of MLS and ACE-FTS as the altitude increased from 40 km, and the difference was approximately 0.4–0.5 ppbv (12 %–15 %) at 50–60 km. The difference between SMILES and TELIS HCl observations was about 0.3 ppbv in the polar winter region between 20 and 34 km, except near 26 km. SMILES HCl error sources that may cause discrepancies with the other observations are investigated by a theoretical error analysis. We calculated errors caused by the uncertainties of spectroscopic parameters, instrument functions, and atmospheric temperature profiles. The Jacobian for the temperature explains the negative bias of the SMILES HCl concentrations at 50–60 km.

1 Introduction

Hydrogen chloride (HCl) is the most abundant species of all total inorganic chlorine compounds Cly in the stratosphere. More than 95 % of Cly exists as HCl at about 1 hPa ( 50 km) (see Froidevaux et al.2008). Thus, HCl has been used for long-term quantitative estimates of the total budget of stratospheric chlorine (WMO2010). Stratospheric HCl vertical profiles have been observed globally by several satellites. The Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite measured stratospheric HCl profiles continuously from late 1991 to November 2005 (Russell III et al.1996). Comparisons of HCl profiles obtained by the Aura Microwave Limb Sounder (MLS) and the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) have been reported by Froidevaux et al. (2008) and Mahieu et al. (2008). These profiles agreed within 5 % (0.15 ppbv) at 0.5 hPa ( 53 km). The Aura/MLS and ACE-FTS HCl concentrations were larger (by 10 %–20 %) than those from HALOE (WMO2010). Balloon-borne remote sensing in the polar winter region was performed using the TErahertz and submillimeter LImb Sounder (TELIS) from 2009 to 2011 (Birk et al.2010). Xu et al. (2018) reported the error analysis and comparison of some trace gas profiles provided by TELIS.

The Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) was launched in 2009 to observe the atmospheric compositions of ozone and species related to the stratospheric ozone destruction cycle, including HCl (Kikuchi et al.2010). The SMILES mission is a joint project of the Japan Aerospace Exploration Agency and the National Institute of Information and Communications Technology (NICT). SMILES is an instrument used to conduct limb observations from the International Space Station (ISS) using supersensitive 4 K superconducting receivers in the submillimeter-wave regions (625 and 649 GHz bands). The non-Sun-synchronous circular orbit of the ISS gave us an opportunity for observations of HCl profiles at different local times. The unprecedented low noise of the SMILES instrument provided a sensitivity 10 times superior to those of previous microwave/submillimeter limb emission instruments used to observe HCl spectra from space (Kikuchi et al.2010). The SMILES observations were used to reveal small abundances of atmospheric species in the stratosphere and mesosphere (e.g., Sato et al.2017; Yamada et al.2020).

A comparison of the HCl profiles inside the Antarctic vortex has been performed for 19–24 November 2009 using the SMILES, Aura/MLS, and ACE-FTS data (Sugita et al.2013). The SMILES HCl values (version 2.1.5 of NICT products) agreed within 10 % (0.3 ppbv) with those of the MLS and ACE-FTS HCl data between the potential-temperature levels of 450 and 575 K, and 425 and 575 K, respectively. In contrast to the previous study, this study used 3.0.0 (v300) SMILES data and provided a global comparison, including the stratosphere and mesosphere. The HCl climatology using SMILES NICT version 2.1.5 was reported by Kreyling et al. (2013).

In this study, we examined HCl vertical profiles from 16 to 100 km, including the upper mesosphere and lower thermosphere (MLT) region. We used the SMILES NICT level-2 product v300, which was released in late 2012 (http://smiles.nict.go.jp/pub/data/index.html, last access: 17 October 2020). HCl vertical profiles from the upper troposphere to the lower thermosphere are reported for the first time. We perform a validation of the SMILES HCl vertical profiles by comparisons with the corresponding global model results of SD-WACCM, satellite observations from MLS and ACE-FTS, and balloon-borne observations from TELIS, and we provide a SMILES HCl error analysis.

The SMILES HCl observations and retrieval procedure are described in Sect. 2. The HCl vertical profile derived by the SMILES product and comparisons of the HCl distribution between 40 S and 60 N with SD-WACCM are shown in Sect. 3. Results of HCl comparisons between SMILES and other instruments are described in Sect. 4. An estimation of the systematic and random errors is presented in Sect. 5. Section 6 describes our conclusions.

2 SMILES HCl observations

The SMILES instrument had been attached to the Japanese Experimental Module (JEM) on the ISS since September 2009. SMILES operational period started on 12 October 2009 and ended on 21 April 2010. The ISS has a non-Sun-synchronous circular orbit with an inclination angle of 51.6 with respect to the Equator. The height of the ISS changed slowly over the observational period ranging from 340 to 360 km. We tilted the line of sight at a 45 angle in the direction of the forward movement of the ISS to observe the northern polar region. The SMILES parameters are summarized in Table 1. The details of the SMILES observations are described in Kikuchi et al. (2010) and in an ozone-validation paper by Kasai et al. (2013), respectively.

Table 1SMILES specification

a Superconductor–insulator–superconductor. b High-electron-mobility transistor.

Download Print Version | Download XLSX

Herein, we briefly describe the SMILES atmospheric observations in terms of the HCl spectra. SMILES had three frequency bands in the submillimeter-wave regions: band A (624.32–625.52 GHz), band B (625.12–626.32 GHz), and band C (649.12–650.32 GHz). Each single scan provided a combination of two out of three of the frequency bands: bands A + B, C + B, or C + A. The combinations changed on a daily basis, as described in Kikuchi et al. (2010). The rotational transitions of two HCl isotopologues were located in bands B and A for H35Cl at 625.9 GHz and H37Cl at 624.9 GHz, respectively. Two acousto-optical spectrometers (AOSs) were equipped in the SMILES instrument. The spectral resolution obtained by AOS1 and AOS2 was 0.8 MHz. The band configuration for each HCl isotopologues had three patterns as follows:

  • H37Cl observation by band A with AOS1,

  • H37Cl observation by band A with AOS2, and

  • H35Cl observation by band B with AOS2.

Band B was always associated with AOS2.

Figure 1 shows the number of observations per day obtained for H35Cl (band B) and H37Cl (band A) for 5 latitude bins during the SMILES observational period. It should be noted that the sampling is not homogeneously distributed for the SMILES observations. The SMILES instrument was sometimes not in operation, for example, when the ISS was boosted up from a low to a high altitude, and the solar panels disturbed the observational line of sight. The SMILES observational latitude range was 65 N–38 S in nominal operations and changed to 38 N–65 S when the ISS performed a yaw maneuver of 180, which happened three times during the 7-month observational period.

We used the SMILES NICT level-2 product v300 for this study. The general method of v300 for the HCl retrieval is similar to that of the level-2 product version 2.1.5 (v215) (Kasai et al.2013). The major updates from the v215 product are as follows: (1) an improvement in the spectrum calibration, particularly for the gain non-linearity of the receiver system (Ochiai et al.2013), (2) an improvement in the accuracy of the tangent height estimation (Ochiai et al.2013), and (3) an update of the temperature retrieval. The temperature profile used in the SMILES v300 retrieval was synthesized, assuming the hydrostatic equilibrium, using the Goddard Earth Observing System version 5 (GEOS-5) reanalysis meteorological datasets and the climatology based on the Aura/MLS measurements (Kuribayashi et al.2017). The GEOS-5 datasets were used in the upper troposphere and stratosphere, and the Aura/MLS datasets were used in the mesosphere and lower thermosphere, respectively. The temperature profile was also retrieved using the ozone transition in the SMILES v300 retrieval process, but it was not applied to the retrieval of atmospheric species including HCl, except for ozone, to avoid systematic error propagation issues (SMILES-NICT2014). Level-1b spectrum data version 008 were used for the SMILES NICT level-2 product v300 (Ochiai et al.2013). The spectroscopic parameters for the HCl retrieval were based on Cazzoli and Puzzarini (2004) and Baron et al. (2011), as summarized in Table 2. We used only the H35Cl data because the intensities of the H37Cl spectra were weaker than those of H35Cl, as shown in Fig. 2.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f01

Figure 1Number of HCl observations for a 5 latitude bin in 1 d for the range of SMILES observational latitudes and period (12 October 2009–21 April 2010). Panels (a) and (b) show the number of observations from band A (H37Cl) and band B (H35Cl), respectively.

Download

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f02

Figure 2(a) An example of the H37Cl spectra observed by SMILES from a single-scan measurement. Tangent heights are approximately 50, 65, and 80 km. Center frequency is 624.98 GHz. (b) An example of the H35Cl spectra for the same condition as for the H37Cl spectra. Center frequency is 625.92 GHz.

Download

Table 2Spectroscopic line parameters measured in the laboratory and used in the processing.

a Laboratory measurements by Cazzoli and Puzzarini (2004). b Baron et al. (2011).

Download Print Version | Download XLSX

The criteria for the selection of the H35Cl profiles used in this study are as follows:

  1. χ2 less than 0.8. χ2 is defined by

    (1) χ 2 = y - F ( x ) T S ϵ - 1 y - F ( x ) + x - x a T S a - 1 x - x a .

    y is the observed spectrum, Sϵ is the covariance matrix of the measurement noise, xa is the a priori state, Sa is the covariance matrix of xa, and F(x) is the forward model depending on the state vector x. We used the US standard atmosphere profiles as the a priori state (xa) (US-Standard1976). They are used separately for polar, equatorial, summer midlatitude, or winter midlatitude regions. The forward model is essentially the atmospheric radiative transfer and the instrument model.

  2. Measurement response (MR) between 0.8 and 1.2 for each vertical grid. MR is defined by

    (2) MR [ i ] = j A [ i , j ] .

    A is the averaging kernel matrix. A value of MR near unity indicates that most of the information in the retrieval results is provided by observations. A low value of MR indicates that the retrieval results are largely influenced by the a priori state and are forced to be identical to a priori values.

An example H35Cl profile and its averaging kernel are presented in Fig. 3 for the single-scan spectrum observation at latitude of 38.3 N, longitude of 12.7 E, and solar zenith angle (SZA) of 160.9 on 15 November 2009. We assumed that the natural isotopic abundance of H35Cl∕HCl was 0.7576 by Berglund and Wieser (2011) for these retrievals. The peak values of the averaging kernels for the SMILES H35Cl profiles exceeded 0.8 in the height range of 16 to 90 km ( 100 to 0.001 hPa) as shown in Fig. 3. The full width at half maximum (FWHM) of the averaging kernels is approximately 3–4 km at altitudes of 18–50 km and becomes greater than 5 km for the altitudes 50–90 km.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f03

Figure 3(a) An example of a HCl profile retrieved from the H35Cl spectrum on 15 November 2009 at a latitude of 38.3 N, longitude of 12.7 E, and SZA of 160.9. The solid red line indicates the retrieved HCl volume mixing ratio (VMR), and the error bar is the root sum square of the smoothing and measurement errors. The dashed black line is an a priori HCl profile used for the retrieval procedure. Panel (b) shows the averaging kernel, the measurement response (indicated by the thick black line), and the FWHM of the averaging kernel profile (indicated by the thick dashed black line).

Download

3 Vertical and latitudinal distribution of SMILES HCl

In this section, we provide the HCl vertical profiles derived from the SMILES NICT level-2 product v300. Figure 4a shows the HCl zonal mean distribution in the latitudinal range from 40 S to 60 N. The time period is between 16 October 2009 and 17 April 2010. The HCl vertical distribution shows an increase with altitude with a maximum below the stratopause ( 45 km), approximately constant values between the stratopause and the upper mesosphere ( 80 km), and a decrease with altitude from the mesopause to the lower thermosphere ( 100 km). In the lower and middle stratosphere, HCl is generated by the reaction of Cl with CH4 and HO2 and transported by circulation (e.g., Brewer–Dobson circulation). The HCl abundance is balanced by production (Cl+HO2HCl+O2) and loss (HCl+OHCl+H2O,HCl+hνH+Cl) in the upper stratosphere and mesosphere. Near the mesopause, the photodissociation becomes the dominant reaction, and the HCl abundance decreases with height (Brasseur and Solomon2005).

This behavior was reproduced by the Whole Atmosphere Community Climate Model version 4 (WACCM4) in the specified dynamics configuration (SD-WACCM) (Lamarque et al.2012; Limpasuvan et al.2016); see Fig. 4b. In the specified dynamics configuration, the simulated meteorological fields in the troposphere and stratosphere are constrained to the Global Modeling and Assimilation Office Modern-Era Retrospective Analysis for Research and Applications (Rienecker et al.2011). The horizontal resolution of SD-WACCM simulation is 1.9 and 2.5 for latitude and longitude, respectively. The vertical grid is from the ground to 150 km with 88 levels, and the time resolution is 3 h. Differences of HCl abundance between SMILES and SD-WACCM are shown in Fig. 4c. For these comparisons, the SD-WACCM profiles are interpolated linearly for each SMILES observation point and time from the nearest model positions and times. Figure 4d shows the number of selected SMILES observations. The SMILES HCl abundance agrees well with the SD-WACCM simulation (within 0.1 ppbv) in the stratosphere and middle mesosphere. At 70–100 km, the difference of the HCl VMR between SMILES and WACCM in the upper mesosphere and lower thermosphere is 0.4–1.0 ppbv, which is larger than the total systematic error in the SMILES HCl VMR estimated by quantitative error analysis (0.1–0.4 ppbv). Thus, the difference is significant. It might be caused by a combination of the SMILES observation and the WACCM model uncertainties. The underestimation of the model could be because of uncertainties in the HCl photodissociation and the reaction with OH radical, which are dominant destruction mechanisms of HCl for the mesosphere/thermosphere (Brasseur and Solomon2005).

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f04

Figure 4Zonal mean distribution of HCl from SMILES, SD-WACCM, and their differences (SMILES – SD-WACCM) are shown in panels (a), (b), and (c), respectively. Each panel displays distributions as a function of latitude (x axis) and altitude (y axis). The latitude range covers between 40 S and 60 N. The number of HCl profiles observed by SMILES is described in panel (d) as a histogram versus latitude. The time period is between 16 October 2009 and 17 April 2010. The data are averaged within latitude bins of 10.

Download

4 Comparison of SMILES HCl profile with those obtained using other instruments

We performed a comparison between SMILES HCl products and other datasets obtained using the MLS on the Aura satellite, the ACE-FTS on the SCISAT satellite, and the balloon-borne TELIS instrument. The characteristics of instruments and datasets are summarized in Table 3. The mean absolute difference, Δabs, for each altitude between SMILES and the other instrument is defined as

(3) Δ abs ( z ) = 1 N ( z ) i = 1 N ( z ) ( x c , i ( z ) - x s , i ( z ) ) ,

where xs,i(z) and xc,i(z) are the HCl volume mixing ratio (VMR) of ith coincidence at an altitude z for SMILES and the other instrument, respectively. N(z) is the number of coincidence at an altitude z. The mean relative difference in percentage, Δrel, for each altitude is given by

(4) Δ rel ( z ) = 1 N ( z ) i = 1 N ( z ) x c , i ( z ) - x s , i ( z ) ( x c , i ( z ) + x s , i ( z ) ) / 2 × 100 .

The coincidence criteria between the SMILES observations and those of the other instruments were set to within 2 for latitude, 8 for longitude, and 5 h for time, following the previous work on MLS observations (Froidevaux et al.2008).

Table 3Characteristics of instruments and datasets used in the comparison.

Download Print Version | Download XLSX

4.1 Comparison with Aura/MLS

The Aura satellite has a Sun-synchronous orbit with an ascending node at 13:45 LT (Waters et al.2006). The measurements cover latitudes from 82 S to 82 N and provide approximately 3500 vertical profiles each day. We used version 4.2 of the Aura/MLS HCl profiles for these comparisons. The Aura/MLS team provided two HCl products: one from spectral band 13 (HCl-640-B13) and one from band 14 (HCl-640-B14). The HCl-640-B13 product is available for the pressure region between 100 and 0.32 hPa. The single profile precision of HCl-640-B13 is less than 0.8 ppbv (25 %) in the stratosphere, and the estimated accuracy is approximately 0.3 ppbv (10 %) (Livesey et al.2018).

The vertical and horizontal resolutions are 3–5 km and 200–400 km, respectively. The vertical resolution of MLS is of the same order as that of SMILES. We used the recommended parameters “status,” “quality,” and “convergence” for screening the Aura/MLS data based on Livesey et al. (2018). The good profiles were selected using (1) quality >1.2, (2) convergence <1.05, and (3) status is even. The SMILES profiles were linearly interpolated onto the MLS pressure grid for the comparisons.

Figure 5 shows a comparison of the HCl vertical profiles derived from SMILES and Aura/MLS in the SMILES latitude range. The total number of coincident profile pairs was 4356. The HCl profiles from the SMILES measurements agreed with those of MLS within 0.25 ppbv between 20 and 40 km (50–3 hPa), while large biases of about 0.4 ppbv were confirmed between 40 and 55 km (3–0.32 hPa). The differences between HCl profiles increase from the stratosphere to the mesosphere. Figure 6 shows the zonal mean absolute differences averaged for 5 latitude bins. The differences did not change with latitude in the altitude region between 35 and 55 km (7–0.32 hPa), while below 35 km altitude, several areas of large differences can be identified, especially in the equatorial region.

Figure 7 displays more detailed features of the differences between the SMILES and MLS profiles. Figure 7a shows the averaged HCl vertical profiles of SMILES (blue) and MLS (red), and the absolute difference in each latitude region (40–20 S, 20 S–20 N, 20–50 N, and 50–65 N). Figure 7b shows the latitudinal distribution averaged in 5 bins for SMILES (blue) and MLS (red) at 30, 40, and 50 km. The error bar shows the 1σ variability for each latitude grid. In the equatorial regions, the SMILES HCl profiles were 0.5 ppbv larger than those of MLS at 30 km (10 hPa). On the other hand, the HCl profiles from the SMILES measurements were 0.4 ppbv lower than those of MLS in all latitude regions at 50 km.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f05

Figure 5Comparison between SMILES and MLS profiles from 24 to 27 January 2010. (a) Mean HCl VMR values (solid lines) and 1σ (shaded areas) for SMILES and MLS. The blue and red lines indicate the SMILES HCl profiles and MLS profiles, respectively. (b) The absolute difference between the SMILES and MLS profiles calculated using Eq (3). (c) The relative difference between the SMILES and MLS profiles calculated using Eq (4).

Download

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f06

Figure 6(a) SMILES zonal mean HCl profiles versus latitude. Profiles that satisfied coincidence criteria between SMILES and MLS for 24–27 January 2010 were used here. (b) Zonal mean HCl profiles from MLS. (c) The absolute difference of the matched pairs versus latitude.

Download

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f07

Figure 7(a) Mean HCl profiles, 1σ uncertainty based on variability (shaded region) and the absolute differences for each latitude region. Upper row: the blue and red lines indicate the mean profiles from SMILES and MLS for 40–20 S, 20 S–20 N, 20–50 N, and 50–65 N (left to right panels). Lower row: the absolute difference between the SMILES and MLS profiles calculated using Eq. (3) for each latitude region. (b) Latitudinal variation of SMILES and Aura/MLS HCl profiles at three altitude levels. The profiles were averaged for each 5 bin. Upper row: mean of SMILES (blue) and MLS (red) profiles for each altitude level. The error bars indicate 1σ uncertainties. Lower row: the means and 1σ of the absolute differences are displayed.

Download

4.2 Comparison with SCISAT/ACE-FTS

The ACE-FTS is an instrument mounted on the Canadian SCISAT satellite. SCISAT moves along an orbit at a 650 km altitude and is inclined at 74 to the Equator (Bernath et al.2005). The ACE comprises two instruments: the Fourier Transform Spectrometer (ACE-FTS) and the Measurement of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation (ACE-MAESTRO). The HCl observation has been performed using ACE-FTS. The ACE-FTS observes solar occultation spectra in the infrared spectral region (750–4400 cm−1) with a high spectral resolution (0.02 cm−1). We used ACE-FTS HCl profiles from the version 4.0 software, which is the latest data version (Boone et al.2020). The vertical resolution of the ACE-FTS HCl retrieval is 3–4 km, and the values from the retrieval grid are interpolated onto the 1 km grid, using a piecewise quadratic method (Bernath et al.2005). The vertical resolution of ACE-FTS is of the same order as that of SMILES. In this study, we used the data within ± 3 times the median absolute deviation (3-MAD) from the median to remove significantly large positive and negative biases based on Boone et al. (2019).

Figure 8 presents the comparison of the HCl vertical profiles between SMILES and ACE-FTS in the latitude range of 40 S–65 N. The total number of coincident data pairs was 935. The SMILES HCl profiles agreed with those of the ACE-FTS by 0.2 ppbv between 20 and 40 km (50–2 hPa), while 0.5 ppbv (15 %) lower than those of ACE-FTS above 40 km. This negative bias in the SMILES HCl concentration was confirmed as in the case of the MLS–SMILES comparison. The results of the comparison with the ACE-FTS for each latitudinal region are presented in Fig. 9a in a manner similar to Fig. 7a. In the tropical region, a large discrepancy was observed at an altitude of 25 km because of relativity poor sampling number. The HCl profiles observed using the SMILES were about 0.5 ppbv (15 %) less than those of ACE-FTS for each latitudinal region above 50 km.

The period of the SMILES H35Cl observation (from 16 October 2009 to 17 April 2010) was covered by the ACE-FTS observation period, while the MLS band-13 observation period had an overlap of only 4 d. We analyzed the difference between SMILES and ACE-FTS for each month and confirmed the seasonal variation of the bias. Figure 9b shows a seasonal variation of the difference between the SMILES and ACE-FTS profiles in the Northern Hemisphere (30–65 N latitude range). The upper panels in Fig. 9b show a mean of the SMILES (blue) and ACE-FTS (green) profiles for each month at three altitude levels. The difference for every month, shown in the lower panels, was represented by the mean for each month at each altitude, and the error bar showed the 1σ standard deviation. No significant seasonal dependence of the difference was observed. The difference of HCl value from these measurements was consequently about 0.5 ppbv at 50 km.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f08

Figure 8Comparison between SMILES and ACE-FTS profiles from 16 October 2009 to 17 April 2010. (a) The mean HCl VMR values (solid lines) and 1σ (shaded areas) for SMILES and ACE-FTS. The blue and green lines are the SMILES and ACE-FTS HCl profiles, respectively. (b) The absolute difference between the SMILES and ACE-FTS profiles calculated using Eq. (3). (c) The relative difference between the SMILES and ACE-FTS profiles calculated using Eq. (4).

Download

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f09

Figure 9(a) The mean of HCl values, 1σ, and the absolute difference for each latitudinal region are shown. Upper row: the blue and green lines indicate the SMILES HCl profiles and ACE-FTS profiles for each latitudinal region (40–20 S, 20–20 N, 20–50 N, and 50–65 N) from the left side. Lower row: the absolute difference between the SMILES and ACE-FTS profiles calculated using Eq. (3) for each latitudinal region. (b) Seasonal variation in the SMILES and ACE-FTS HCl profiles at three altitude levels. Latitudinal range is 30–65 N. Upper row: monthly mean of SMILES (blue) and ACE-FTS (green) HCl profiles for each altitude level. The error bars correspond to 1σ. Lower row: the mean and 1σ uncertainty of the absolute difference.

Download

4.3 Comparison with balloon-borne instrument TELIS

TELIS (Birk et al.2010) is a balloon-borne THz/submillimeter-wave spectrometer with superconductive 4 K technology similar to that of SMILES. TELIS was one of the instruments used in the balloon observation campaigns for three winters (2009–2011) in Kiruna, Sweden. Limb observations were performed from a flight altitude of 30–35 km, toward the stratosphere and the upper troposphere with 1.5–2 km vertical sampling. TELIS was used to observe H37Cl at the transition frequency of 1873.4 GHz in the 1.8 THz channel (level-1b version 3v02), with corrected non-linearity problems in the radiometric calibration using a quadratic term for each frequency segment. The details on the TELIS 1.8 THz channel and its L1 data processing are shown in Suttiwong (2010).

In this study, we compared two H37Cl profiles obtained by TELIS (observation number 20044 and 21537) with the SMILES profiles (band B and identifiers 761 and 762) in terms of geolocation over Kiruna (67.8 N, 20.4 E) and time on 24 January 2010. The total error of the H37Cl profile derived from the TELIS measurement was estimated to be 0.25–0.5 ppbv, with a vertical resolution of 3 km (Xu et al.2018).

Figure 10 shows a comparison of the HCl vertical profiles of SMILES and TELIS. Figure 10a–c show the SMILES and TELIS HCl vertical profiles, their absolute difference, and their relative difference, respectively. We converted the H37Cl amount to HCl using the natural isotopic abundance (H37Cl ∕ HCl = 0.2424). The SMILES and TELIS HCl profiles agreed well, within 0.5 ppbv from 17 to 34 km. A depletion in the HCl profile was observed below 25 km in both the SMILES and TELIS observations. This result is considered to have been caused by chlorine activation in the polar vortex (Webster et al.1993; Wegner et al.2016). Both SMILES and TELIS profiles agree well in general, but the TELIS profiles are larger than the SMILES profiles above 32 km (8 hPa). The H37Cl line is still rather strong at higher altitudes, and the dominant error source of the TELIS data stems from the non-linearity in the calibration process, which shows that even a small uncertainty may result in significant errors in the retrieval (Xu et al.2018). Table 4 summarizes the results of the comparison with other observations.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f10

Figure 10Differences between SMILES and TELIS HCl observations are shown for 24 January 2010. Two of the SMILES observations were taken from closer points, where the latitude, longitude, and SZA were 64.3 N, 30.6 E, and 84.6 SZA and 64.8 N, 38.2 E, 86.2 SZA, respectively. (a) Each profile obtained by SMILES and TELIS instruments; (b) the absolute difference between the SMILES and TELIS profiles calculated using Eq. (3); (c) the relative difference between the SMILES and TELIS profiles calculated using Eq. (4).

Download

Table 4Summary of the HCl comparison study between SMILES and the other instruments.

Download Print Version | Download XLSX

5 Theoretical error analysis

We have evaluated the total error in the HCl vertical profiles observed by SMILES, and we discussed the cause of the bias observed in the comparison study; see Sect. 4.

5.1 Estimation of total error

We employed a perturbation method to estimate the total error of the SMILES HCl profile. The details of the perturbation method of the SMILES error analysis have been described in Kasai et al. (2006) and Sato et al. (2014), Kasai et al. (2013), and Sato et al. (2012) for ozone isotopes, ozone, and ClO, respectively. We assumed an averaged HCl profile within the coincidence with MLS in the southern midlatitude region (20–40 S) as a reference. The total error (Etotal) for each altitude grid was calculated by

(5) E total [ i ] = E noise [ i ] 2 + E smooth [ i ] 2 + E param [ i ] 2 ,

where Enoise is the error due to spectrum noise, Esmooth is the smoothing error, and Eparam is the model parameter error. Enoise and Esmooth were calculated by

(6)Enoise[i]=Snoise[i,i](7)Snoise=DSϵDT

and

(8)Esmooth[i]=Ssmooth[i,i](9)Ssmooth=(U-A)Sa(U-A)T,

where D is the contribution function, and U is the unit matrix.

We calculated the errors due to the uncertainties of the spectroscopic parameters and the instrument functions to calculate Eparam. The error Eparam was calculated by

(10) E param = I y ref , b 0 + Δ b 0 - I y ref , b 0

where I is the inversion function, b0 is the vector of model parameters, and Δb0 is the uncertainty on model parameters. yref is the reference spectrum calculated using the reference profile. The details of the estimation of the total error are described in Sato et al. (2012). The error sources and perturbations for the model parameter used in this study are summarized in Table 5. These parameters and perturbation values are based on Sato et al. (2014). The uncertainties of the spectroscopic parameters of the O3 transition 625.37 GHz were included to estimate the error due to the interference from the line shape of O3 spectrum. The calibration error was not considered in this study because the latest L1b data (version 008) were used, and Sato et al. (2014) reported that the error due to the spectrum calibration in these L1b data was insignificant.

Table 5Error sources and their perturbation values for this study.

The perturbation value was based on Sato et al. (2014).

Download Print Version | Download XLSX

The estimated errors in the SMILES HCl v300 product are presented in Fig. 11. The total model parameter error, labeled as “param”, was calculated by a root sum square (RSS) of all model parameter errors. The three spectroscopic parameters, line intensity, air-pressure-broadening coefficient (γair), and its temperature dependence (nair) were dominant error sources below 30 km. The γair was a major error source between 30 and 60 km, which was about 0.09 ppbv ( 3 %) at 50 km. At altitudes of 60–90 km, the largest error source was the AOS response function, and its peak value reached 0.38 ppbv ( 12 %) at 70 km.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f11

Figure 11Summary of the errors for a single-scan observation. Panels (a) and (b), respectively, show the absolute difference and relative difference, which were estimated using the perturbation method. The green marker indicates an error due to the uncertainty in line intensity (“intensity”). The red and blue squares show errors for the air-broadening coefficient (“γair”) and its temperature dependence (nair) of the H35Cl transition. The red and blue circles show errors for “γair” and nair of the O3 transition (“γair (O3)”, “nair (O3)”). The yellow and cyan symbols show errors for the antenna beam pattern (“antenna”), and the AOS response function (“AOS”). The dashed red and blue lines indicate Enoise (“noise”) and Esmooth (“smooth”), respectively. The solid black lines (“param”) indicate the total model parameter errors, obtained from the RSS of all the estimated parameter errors. The dotted and dashed black lines indicate the difference between SMILES and MLS, and between SMILES and ACE-FTS.

Download

5.2 Discussion: cause of the negative bias of the SMILES HCl vertical profile

In this section, we discuss the cause of the negative bias of the SMILES HCl vertical profile, i.e., approximately 10 % less than those of Aura/MLS and ACE-FTS especially at the altitudes between 40 and 60 km. Such a large bias could not be explained by the total error estimated by the perturbation method described in Sect. 5.1. We further investigated the cause of this bias by looking at the temperature profile used in the retrieval calculation. The temperature profile used for the retrieval procedure of SMILES was lower than those of MLS and ACE-FTS particularly in the upper stratosphere and mesosphere. We estimated the difference of the retrieved HCl vertical profile, Δx, due to the difference of the temperature profile, ΔT, as follows:

(11) Δ x = DK T Δ T , D = x y , K T = y T .

The Jacobian KT indicates the sensitivity of the spectral brightness temperature (y) with reference to changes of the temperature (T). Here, we synthesized the Jacobian with a perturbation of 0.5 K. Figure 12a shows the Jacobian as a function of tangent height. Here, the minimum value of the column of the KT matrix is plotted. A negative Jacobian value means that higher temperatures induce a lower brightness temperature spectrum, thus increasing the HCl abundance in the retrieval to compensate for this underestimation. The temperature profiles used in the retrievals by SMILES, MLS, and ACE-FTS and their differences (MLSSMILES, ACE-FTSSMILES) are shown in Fig. 12b and c, respectively. The vertical profile of temperature of SMILES is approximately 5–10 K lower than those of both MLS and ACE-FTS for altitudes between 50 and 60 km. The a priori temperature profile used in the SMILES retrieval procedure is based on the GEOS-5 profile in the stratosphere and MLS retrieved profile in the mesosphere and above. The altitude limit of the MLS temperature profile is 0.001 hPa, with a vertical resolution of 6–14 km and a precision of 1.2–3.6 K per profile. MLS used GEOS-5 up to 1 hPa as with SMILES. For pressures smaller than 1 hPa, the Committee on Space Research (COSPAR) International Reference Atmosphere (CIRA-86) is used as a priori temperature information (with a loose constraint) in the MLS retrieval procedure (Schwartz et al.2008). The altitude range and vertical resolution of the CIRA-86 profile are ground to 120 and 2 km, respectively (Fleming et al.1990). The MLS temperature profiles in the stratosphere and above were retrieved primarily from bands near the 118 GHz O2 spectral line (Livesey et al.2018). The temperature value retrieved by MLS is 10 K lower than the a priori profile on average in some areas for pressure values smaller than 1 hPa, based on earlier version validation studies (Schwartz et al.2008). The ACE-FTS HCl retrieval procedure uses the temperature profile retrieved from the ACE-FTS measurements between 18 and 125 km. The ACE-FTS temperature profiles were retrieved from CO2 VMR using hydrostatic equilibrium (Boone et al.2020). The vertical resolution of ACE-FTS retrieved temperature is 3–4 km. The temperature values retrieved by ACE-FTS are less than 10 K larger than the MLS derived temperatures (Schwartz et al.2008). These types of differences are also seen in the comparison results performed here.

Figure 12d shows the Δx due to the difference of the temperature profile. This temperature difference caused an increase in SMILES HCl of 0.12 and 0.20 ppbv at 50–60 km for the MLS and ACE-FTS comparisons, respectively. In the comparison study (see Sect. 4), we confirmed the negative bias in the SMILES HCl vertical profile of 0.4 ± 0.38 and 0.5 ± 0.28 ppbv between 40 and 60 km for MLS and ACE-FTS, respectively. We estimated the total error by the perturbation method and investigated the temperature profile used in the retrieval calculation, in order to investigate the cause of this negative bias. The largest error sources were the uncertainties in γair of the H35Cl transition and the temperature profiles used in the retrievals. We assumed a 3 % uncertainty in γair which could lead to a 0.1 ppbv bias in the 40–60 km region. In addition to the error due to the γair coefficients, the effect of temperature differences should be taken into account at altitudes above 50 km. The gradual increase in bias is caused by the difference in altitude at which these two errors become more pronounced. The difference in temperature profiles used in the retrieval between SMILES and ACE-FTS caused a negative bias of about 0.2 ppbv at 50–60 km. In summary, 0.3 in 0.5 ppbv (60 %) negative bias between SMILES and ACE-FTS can be explained by the uncertainty in γair and the temperature profiles used in the SMILES retrievals. The effect of γair error on the negative bias between SMILES and MLS is less than that of SMILES and ACE-FTS, since SMILES and MLS observed the same H35Cl transition lines and the values of the γair are consistent within approximately 1 % (3.39 MHz/Torr for SMILES and 3.42 MHz/Torr for MLS; Drouin2004). A 1 % difference in γair might cause the HCl abundance to increase by about 0.03 ppbv according to our error analysis with a perturbation method. About 40 % (0.15 out of 0.4 ppbv) of the negative bias between SMILES and MLS can be explained. Therefore, our theoretical error analysis shows that the SMILES HCl has a negative bias of at most 0.25 ppbv between 40 and 60 km; the remaining difference between SMILES and MLS or ACE-FTS can be explained by the standard deviation in the comparison result.

https://amt.copernicus.org/articles/13/6837/2020/amt-13-6837-2020-f12

Figure 12(a) The lines represent the Jacobian for the temperature calculated using the mean profile from the SMILES winter midlatitude region (20–40 S). The units of the Jacobian (x axis) are Kelvin. (b) Temperature profile of each instruments. The solid blue and red lines indicate the mean SMILES and MLS temperature profiles within the coincidence criteria used for the HCl comparisons. The dashed blue and green lines indicate the SMILES and ACE-FTS temperature profiles for the same coincidence criteria. (c) The absolute difference of temperature profiles between SMILES and other instruments calculated using Eq. (3). (d) The difference in HCl profiles calculated using Eq. (11).

Download

6 Conclusions

In this study, the HCl vertical profile in a wide range from the upper troposphere to the lower thermosphere was reported for the first time using the SMILES NICT level-2 data product v300. The HCl distribution shows an increase with altitude with a maximum below the stratopause ( 45 km), approximately constant values between the stratopause and the upper mesosphere ( 80 km), and a decrease with altitude from the mesopause to the lower thermosphere ( 100 km). In the lower and middle stratosphere, HCl is generated by the reaction of Cl with CH4 and HO2 and transported by circulation (e.g., Brewer–Dobson circulation). The HCl abundance is balanced by production (Cl+CH4HCl+CH3,Cl+HO2HCl+O2) and loss (HCl+OHCl+H2O, HCl+hνH+Cl) in the upper stratosphere and the mesosphere. Above the mesopause, the photodissociation becomes the dominant reaction and the HCl abundance decreases. This behavior was reproduced by the SD-WACCM model, and the SMILES HCl vertical profile agreed well with the SD-WACCM model within ± 0.1 ppbv for altitudes between 30 and 70 km.

The data quality of the SMILES HCl vertical profile was quantified by comparisons versus other measurements and supported by a theoretical error analysis. We compared the SMILES HCl vertical profiles with well-validated data of two satellite instruments (Aura/MLS and ACE-FTS), as well as a balloon-borne instrument (TELIS) at their temporal–spatial coincidences. The SMILES HCl profiles at 20–40 km showed good agreement, within less than 0.25 ± 0.3 (1σ) and 0.20 ± 0.2 (1σ) ppbv versus MLS and ACE-FTS, respectively. The comparison with TELIS in the polar winter region at 20–34 km showed similar behavior with differences within 0.3 ppbv, which is the same order of magnitude as the systematic error of the TELIS data. A negative bias (< 0.5 ppbv) of the SMILES HCl profiles from 40 to 60 km altitudes was observed in comparisons versus MLS and ACE-FTS HCl profiles.

We estimated the total error for SMILES HCl based on the perturbation method and considering the uncertainties in atmospheric temperature profiles used in the retrievals. The dominant contributions to the systematic errors were from the air-broadening parameter (0.09 ppbv) and the AOS response function (0.38 ppbv) at 30–60 km and 60–100 km altitudes, respectively. The uncertainty in the temperature profile used in the retrieval calculation caused a negative bias of 0.12 to 0.20 ppbv between 50 and 60 km, which was 30 % and 40 % of the HCl abundance difference between SMILES and MLS, and SMILES and ACE-FTS, respectively. The uncertainties of the air-broadening parameter and the temperature profile are capable of contributing a total of 40 %–50 % of the SMILES HCl negative biases at 50–60 km. In summary, our theoretical error analysis showed that the HCl profiles had a negative bias of 0.20–0.25 ppbv at 50–60 km, which is consistent with the observed differences versus MLS and ACE-FTS profiles within 1 standard deviation. The spectroscopic parameters of the HCl transitions and the temperature profile above the stratopause are key parameters for potential improvements in the SMILES retrieval algorithms.

The observation of HCl abundances in the upper atmosphere is important to investigate the long-term total budget of anthropogenic chlorine in the Earth's atmosphere. Further observations and model studies are needed to better understand the sources and sinks, transport processes, and chemical reactions related to HCl.

Data availability

The SMILES data are available at http://smiles.nict.go.jp/pub/data/index.html (last access: 18 November 2020; SMILES-NICT2020). The MLS data are available at https://disc.gsfc.nasa.gov/datasets?page=1&keywords=AURA MLS (last access: 18 November 2020; NASA2020) or see https://mls.jpl.nasa.gov/data/ (last access: 18 November 2020; Aura/MLS2020). The ACE-FTS data are available at http://www.ace.uwaterloo.ca/data.php (last access: 18 November 2020; SCISAT/ACE-FTS2020). The details of the TELIS 1.8 THz channel and its L1 data processing are shown at https://elib.dlr.de/66749/ (last access: 18 November 2020; Suttiwong2010) “Development and characterization of the balloon borne instrument TELIS (TErahertz and submillimeter LImb Sounder): 1.8 THz receiver Suttiwong, Nopporn (2010)” and/or https://elib.dlr.de/97249/ (last access: 18 November 2020; Xu2015) “Inversion for Limb Infrared Atmospheric Sounding Xu, Jian (2015)”. The WACCM data are available at http://www.cesm.ucar.edu/working_groups/Whole-Atmosphere/ (last access: 18 November 2020; WACCM2020).

Author contributions

SN designed the study and performed the analysis. YK designed the study and provided the SMILES data. TOS provided the code for the error analysis and contributed to data analysis and interpretation. TY, TF, and KK contributed to the data analysis and reviewed the manuscript. LF, NJL, KAW, JX, and FS provided MLS, ACE-FTS, and TELIS information and knowledge of HCl in the atmosphere and helped review the manuscript. YJO and VL performed and provided the WACCM simulations, and helped review the manuscript. NK and TM supervised writing of and reviewed the manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “New developments in atmospheric Limb measurements: Instruments, Methods and science applications (AMT/ACP inter-journal SI)”. It is a result of the 10th international limb workshop, Greifswald, Germany, 4–7 June 2019.

Acknowledgements

The JEM/SMILES mission is a joint project of the Japan Aerospace Exploration Agency and the National Institute of Information and Communications Technology. The authors wish to acknowledge the contributions made by their colleagues at JAXA and NICT for managing and supporting the SMILES mission. The authors also thank Kazuya Muranaga (Systems Engineering Consultants Co., Ltd.) and Jan Moller (Molflow Co., Ltd.) for supporting the data processing in the SMILES NICT level-2 product. Work at the Jet Propulsion Laboratory, California Institute of Technology, was performed under a contract with the National Aeronautics and Space Administration. The Atmospheric Chemistry Experiment (ACE), also known as SCISAT, is a Canadian-led mission mainly supported by the Canadian Space Agency.

Financial support

This research has been supported by the Ministry of Internal Affairs and Communications (grant nos. 0155-0285 and 0155-0093). Yasuko Kasai is supported by the Funding Program for Next Generation World-Leading Researchers (NEXT Program) (grant no. GR101).

Review statement

This paper was edited by Miriam Sinnhuber and reviewed by two anonymous referees.

References

Aura/MLS: EOS MLS data, Jet Propulsion Laboratory, available at: https://mls.jpl.nasa.gov/data/, last access: 18 November 2020. a

Baron, P., Urban, J., Sagawa, H., Möller, J., Murtagh, D. P., Mendrok, J., Dupuy, E., Sato, T. O., Ochiai, S., Suzuki, K., Manabe, T., Nishibori, T., Kikuchi, K., Sato, R., Takayanagi, M., Murayama, Y., Shiotani, M., and Kasai, Y.: The Level 2 research product algorithms for the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES), Atmos. Meas. Tech., 4, 2105–2124, https://doi.org/10.5194/amt-4-2105-2011, 2011. a, b

Berglund, M. and Wieser, M. E.: Isotopic compositions of the elements 2009 (IUPAC Technical Report), Pure Appl. Chem., 83, 397–410, https://doi.org/10.1351/PAC-REP-10-06-02, 2011. a

Bernath, P. F., McElroy, C. T., Abrams, M. C., Boone, C., Butler, M., Camy-Peyret, C., Carleer, M., Clerbaux, C., Coheur, P., and Colin, R.: Atmospheric chemistry experiment (ACE): mission overview, Geophys. Res. Lett., 32, L15S01, https://doi.org/10.1029/2005GL022386, 2005. a, b

Birk, M., Wagner, G., Lange, G., Lange, A., Ellison, B., Harman, M., Murk, A., Oelhaf, H., Maucher, G., and Sartorius, C.: TELIS: TErahertz and subMMW LImb Sounder – Project summary after first successful flight, pp. 195–200, 21st international symposium on space terahertz technology, Oxford, 23–25 March, 2010. a, b

Boone, C., Jones, S., and Bernath, P.: Data usage guide and file format description for ACE-FTS level 2 data version 4.0 ASCII format, Tech. rep., Atmospheric Chemistry Experiment Science Operations Center, available at: https://databace.scisat.ca/level2/ace_v4.0/ACE-SOC-0033-ACE-FTS_ascii_data_usage_and_fileformat_for_v4.0.pdf (last access: 18 November 2020), 2019. a

Boone, C., Bernath, P., Cok, D., Jones, S., and Steffen, J.: Version 4 retrievals for the atmospheric chemistry experiment Fourier transform spectrometer (ACE-FTS) and imagers, J. Quant. Spectrosc. Ra., 247, 106939, https://doi.org/10.1016/j.jqsrt.2020.106939, 2020. a, b

Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere: Chemistry and Physics of the Stratosphere and Mesosphere, Springer, P.O. Box 17, 3300 AA Dordrecht, The Netherland, 2005. a, b

Cazzoli, G. and Puzzarini, C.: Hyperfine structure of the J=1<-0 transition of H35Cl and H37Cl: improved ground state parameters, J. Mol. Spectrosc., 226, 161–168, https://doi.org/10.1016/j.jms.2004.03.020, 2004. a, b

Drouin, B. J.: Temperature dependent pressure-induced lineshape of the HCl J=1<-0 rotational transition in nitrogen and oxygen, J. Quant. Spectrosc. Ra., 83, 321–331, https://doi.org/10.1016/S0022-4073(02)00360-6, 2004. a

Fleming, E. L., Chandra, S., Barnett, J., and Corney, M.: Zonal mean temperature, pressure, zonal wind and geopotential height as functions of latitude, Adv. Space Res., 10, 11–59, https://doi.org/10.1016/0273-1177(90)90386-E, 1990. a

Froidevaux, L., Jiang, Y. B., Lambert, A., Livesey, N. J., Read, W. G., Waters, J. W., Fuller, R. A., Marcy, T. P., Popp, P. J., Gao, R. S., Fahey, D. W., Jucks, K. W., Stachnik, R. A., Toon, G. C., Christensen, L. E., Webster, C. R., Bernath, P. F., Boone, C. D., Walker, K. A., Pumphrey, H. C., Harwood, R. S., Manney, G. L., Schwartz, M. J., Daffer, W. H., Drouin, B. J., Cofield, R. E., Cuddy, D. T., Jarnot, R. F., Knosp, B. W., Perun, V. S., Snyder, W. V., Stek, P. C., Thurstans, R. P., and Wagner, P. A.: Validation of Aura Microwave Limb Sounder HCl measurements, J. Geophys. Res.-Atmos., 113, D15S25, https://doi.org/10.1029/2007JD009025, 2008. a, b, c

Kasai, Y. J., Urban, J., Takahashi, C., Hoshino, S., Takahashi, K., Inatani, J., Shiotani, M., and Masuko, H.: Stratospheric ozone isotope enrichment studied by submillimeter wave heterodyne radiometry: the observation capabilities of SMILES, IEEE T. Geosci. Remote, 44, 676–693, https://doi.org/10.1109/TGRS.2005.861754, 2006. a

Kasai, Y., Sagawa, H., Kreyling, D., Dupuy, E., Baron, P., Mendrok, J., Suzuki, K., Sato, T. O., Nishibori, T., Mizobuchi, S., Kikuchi, K., Manabe, T., Ozeki, H., Sugita, T., Fujiwara, M., Irimajiri, Y., Walker, K. A., Bernath, P. F., Boone, C., Stiller, G., von Clarmann, T., Orphal, J., Urban, J., Murtagh, D., Llewellyn, E. J., Degenstein, D., Bourassa, A. E., Lloyd, N. D., Froidevaux, L., Birk, M., Wagner, G., Schreier, F., Xu, J., Vogt, P., Trautmann, T., and Yasui, M.: Validation of stratospheric and mesospheric ozone observed by SMILES from International Space Station, Atmos. Meas. Tech., 6, 2311–2338, https://doi.org/10.5194/amt-6-2311-2013, 2013. a, b, c

Kikuchi, K.-i., Nishibori, T., Ochiai, S., Ozeki, H., Irimajiri, Y., Kasai, Y., Koike, M., Manabe, T., Mizukoshi, K., Murayama, Y., Nagahama, T., Sano, T., Sato, R., Seta, M., Takahashi, C., Takayanagi, M., Masuko, H., Inatani, J., Suzuki, M., and Shiotani, M.: Overview and early results of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES), J. Geophys. Res., 115, D23306, https://doi.org/10.1029/2010JD014379, 2010. a, b, c, d

Kreyling, D., Sagawa, H., Wohltmann, I., Lehmann, R., and Kasai, Y.: SMILES zonal and diurnal variation climatology of stratospheric and mesospheric trace gasses: O3, HCl, HNO3, ClO, BrO, HOCl, HO2, and temperature, J. Geophys. Res.-Atmos., 118, 11,888–11,903, https://doi.org/10.1002/2012JD019420, 2013. a

Kuribayashi, K., Yoshida, N., Jin, H., Orsolini, Y., and Kasai, Y.: Optimal retrieval method to estimate ozone vertical profile in the mesosphere and lower thermosphere (MLT) region from submillimeter-wave limb emission spectra, J. Quant. Spectrosc. Ra., 192, 42–52, https://doi.org/10.1016/j.jqsrt.2017.01.033, 2017. a

Lamarque, J.-F., Emmons, L. K., Hess, P. G., Kinnison, D. E., Tilmes, S., Vitt, F., Heald, C. L., Holland, E. A., Lauritzen, P. H., Neu, J., Orlando, J. J., Rasch, P. J., and Tyndall, G. K.: CAM-chem: description and evaluation of interactive atmospheric chemistry in the Community Earth System Model, Geosci. Model Dev., 5, 369–411, https://doi.org/10.5194/gmd-5-369-2012, 2012. a

Limpasuvan, V., Orsolini, Y. J., Chandran, A., Garcia, R. R., and Smith, A. K.: On the composite response of the MLT to major sudden stratospheric warming events with elevated stratopause, J. Geophys. Res.-Atmos,, 121, 4518–4537, https://doi.org/10.1002/2015JD024401, 2016. a

Livesey, J. N., Read, W. G., Wagner, P. A., Froidevaux, L., Lambert, A., Manney, G. L., Millán Valle, L. F., Pumphrey, H. C., Santee, M. L., Schwartz, M. J., Wang, S., Fuller, R. A., Jarnot, R. F., Knosp, B. W., Martinez, E., and Lay, R. R.: Version 4.2x Level 2 data quality and description document., Tech. Rep., JPL D-33509, Jet Propulsion Laboratory, available at: http://mls.jpl.nasa.gov/data/datadocs.php (last access: 18 November 2020), 2018. a, b, c

Mahieu, E., Duchatelet, P., Demoulin, P., Walker, K. A., Dupuy, E., Froidevaux, L., Randall, C., Catoire, V., Strong, K., Boone, C. D., Bernath, P. F., Blavier, J.-F., Blumenstock, T., Coffey, M., De Mazière, M., Griffith, D., Hannigan, J., Hase, F., Jones, N., Jucks, K. W., Kagawa, A., Kasai, Y., Mebarki, Y., Mikuteit, S., Nassar, R., Notholt, J., Rinsland, C. P., Robert, C., Schrems, O., Senten, C., Smale, D., Taylor, J., Tétard, C., Toon, G. C., Warneke, T., Wood, S. W., Zander, R., and Servais, C.: Validation of ACE-FTS v2.2 measurements of HCl, HF, CCl3F and CCl2F2 using space-, balloon- and ground-based instrument observations, Atmos. Chem. Phys., 8, 6199–6221, https://doi.org/10.5194/acp-8-6199-2008, 2008. a

NASA: National Aeronautics and Space Administration (NASA): GES DISC Data Archive, available at https://disc.gsfc.nasa.gov/datasets?page=1&keywords=AURAMLS, last access: 18 November 2020. a

Ochiai, S., Kikuchi, K., Nishibori, T., Manabe, T., Ozeki, H., Mizobuchi, S., and Irimajiri, Y.: Receiver Performance of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on the International Space Station, IEEE T. Geosci. Remote, 51, 3791–3802, https://doi.org/10.1109/TGRS.2012.2227758, 2013. a, b, c

Rienecker, M. M., Suarez, M. J., Gelaro, R., Todling, R., Bacmeister, J., Liu, E., Bosilovich, M. G., Schubert, S. D., Takacs, L., Kim, G.-K., Bloom, S., Chen, J., Collins, D., Conaty, A., da Silva, A., Gu, W., Joiner, J., Koster, R. D., Lucchesi, R., Molod, A., Owens, T., Pawson, S., Pegion, P., Redder, C. R., Reichle, R., Robertson, F. R., Ruddick, A. G., Sienkiewicz, M., and Woollen, J.: MERRA: NASA's Modern-Era Retrospective Analysis for Research and Applications, J. Clim., 24, 3624–3648, https://doi.org/10.1175/JCLI-D-11-00015.1, 2011. a

Russell III, J. M., Deaver, L. E., Luo, M., Cicerone, R. J., Park, J. H., Gordley, L. L., Toon, G. C., Gunson, M. R., Traub, W. A., Johnson, D. G., Jucks, K. W., Zander, R., and Nolt, I. G.: Validation of hydrogen fluoride measurements made by the Halogen Occultation Experiment from the UARS platform, J. Geophys. Res.-Atmos., 101, 10163–10174, https://doi.org/10.1029/95JD01705, 1996. a

Sato, T. O., Sagawa, H., Kreyling, D., Manabe, T., Ochiai, S., Kikuchi, K., Baron, P., Mendrok, J., Urban, J., Murtagh, D., Yasui, M., and Kasai, Y.: Strato-mesospheric ClO observations by SMILES: error analysis and diurnal variation, Atmos. Meas. Tech., 5, 2809–2825, https://doi.org/10.5194/amt-5-2809-2012, 2012. a, b

Sato, T. O., Sagawa, H., Yoshida, N., and Kasai, Y.: Vertical profile of δ18OOO from the middle stratosphere to lower mesosphere from SMILES spectra, Atmos. Meas. Tech., 7, 941–958, https://doi.org/10.5194/amt-7-941-2014, 2014. a, b, c, d

Sato, T. O., Kuribayashi, K., Yoshida, N., and Kasai, Y.: Diurnal variation of oxygen isotopic enrichment in asymmetric-18 ozone observed by the SMILES from space, Geophys. Res. Lett., 44, 6399–6406, https://doi.org/10.1002/2016GL071924, 2017. a

Schwartz, M. J., Lambert, A., Manney, G. L., Read, W. G., Livesey, N. J., Froidevaux, L., Ao, C. O., Bernath, P. F., Boone, C. D., Cofield, R. E., Daffer, W. H., Drouin, B. J., Fetzer, E. J., Fuller, R. A., Jarnot, R. F., Jiang, J. H., Jiang, Y. B., Knosp, B. W., Krüger, K., Li, J.-L. F., Mlynczak, M. G., Pawson, S., Russell III, J. M., Santee, M. L., Snyder, W. V., Stek, P. C., Thurstans, R. P., Tompkins, A. M., Wagner, P. A., Walker, K. A., Waters, J. W., and Wu, D. L.: Validation of the Aura Microwave Limb Sounder temperature and geopotential height measurements, J. Geophys. Res.-Atmos., 113, D15S11, https://doi.org/10.1029/2007JD008783, 2008. a, b, c

SCISAT/ACE-FTS: Data access for ACE-FTS data products, University of Waterloo, available to: http://www.ace.uwaterloo.ca/data.php, last access: 18 November 2020. a

Smiles-Nict, T.: JEM/SMILES Level2 Research Data (L2r) Products Guide for Version 3.0.0 STRUCTURE OF SMILES L2r PRODUCTS, Tech. rep., National Institute of Information and Communications Technology, available at: https://smiles.nict.go.jp/pub/data/pdf/SMILES-L2r_Guide_v3.0.0.pdf (last access: 18 November 2020), 2014. a

SMILES-NICT: Data access for the SMILES products, National Institute of Information and Communication Technology, available at: https://smiles.nict.go.jp/pub/data/pdf/SMILES-L2r_Guide_v3.0.0.pdf, last access: 18 November 2020. a

Sugita, T., Kasai, Y., Terao, Y., Hayashida, S., Manney, G. L., Daffer, W. H., Sagawa, H., Suzuki, M., Shiotani, M., Walker, K. A., Boone, C. D., and Bernath, P. F.: HCl and ClO profiles inside the Antarctic vortex as observed by SMILES in November 2009: comparisons with MLS and ACE-FTS instruments, Atmos. Meas. Tech., 6, 3099–3113, https://doi.org/10.5194/amt-6-3099-2013, 2013. a

Suttiwong, N.: Development and characterization of the balloon borne instrument TELIS (TErahertz and submillimeter LImb Sounder): 1.8 THz receiver, PhD thesis, University of Bremen, 2010. a, b

US-Standard: U.S Standard Atmosphere 1976, Tech. rep., U.S Government Printing Office, Washington DC, 1976. a

WACCM: WHOLE ATMOSPHERE WORKING GROUP, National Center for Atmospheric Research, available at http://www.cesm.ucar.edu/working_groups/Whole-Atmosphere/, last access: 18 November 2020. a

Waters, J. W., Froidevaux, L., Harwood, R. S., Jarnot, R. F., Pickett, H. M., Read, W. G., Siegel, P. H., Cofield, R. E., Filipiak, M. J., Flower, D. A., Holden, J. R., Lau, G. K., Livesey, N. J., Manney, G. L., Pumphrey, H. C., Santee, M. L., Wu, D. L., Cuddy, D. T., Lay, R. R., Loo, M. S., Perun, V. S., Schwartz, M. J., Stek, P. C., Thurstans, R. P., Boyles, M. A., Chandra, K. M., Chavez, M. C., Chen, G.-S., Chudasama, B. V., Dodge, R., Fuller, R. A., Girard, M. A., Jiang, J. H., Jiang, Y., Knosp, B. W., Labelle, R. C., Lam, J. C., Lee, A. K., Miller, D., Oswald, J. E., Patel, N. C., Pukala, D. M., Quintero, O., Scaff, D. M., Vansnyder, W., Tope, M. C., Wagner, P. A., and Walch, M. J.: The Earth Observing System Microwave Limb Sounder (EOS MLS) on the Aura Satellite, IEEE T. Geosci. Remote, 44, 1075–1092, https://doi.org/10.1109/TGRS.2006.873771, 2006. a

Webster, C. R., May, R. D., Toohey, D. W., Avallone, L. M., Anderson, J. G., Newman, P., Lait, L., Schoeberl, M. R., Elkins, J. W., and Chan, K. R.: Chlorine Chemistry on Polar Stratospheric Cloud Particles in the Arctic Winter, Science, 261, 1130–1134, https://doi.org/10.1126/science.261.5125.1130, 1993.  a

Wegner, T., Pitts, M. C., Poole, L. R., Tritscher, I., Grooß, J.-U., and Nakajima, H.: Vortex-wide chlorine activation by a mesoscale PSC event in the Arctic winter of 2009/10, Atmos. Chem. Phys., 16, 4569–4577, https://doi.org/10.5194/acp-16-4569-2016, 2016. a

WMO: Ozone Depleting Substances (ODSs) and Related Chemicals (Chapter 1), Scientific Assessment of Ozone Depletion 2010, WMO Global Ozone Research and Monitoring Project, Report No. 52, 2010. a, b

Xu, J.: Inversion for Limb Infrared Atmospheric Sounding, Dissertation, Technische Universität München, Munich, Germany, available at: https://elib.dlr.de/97249/ (last access: 18 November 2020), 2015. a

Xu, J., Schreier, F., Wetzel, G., De Lange, A., Birk, M., Trautmann, T., Doicu, A., and Wagner, G.: Performance Assessment of Balloon-Borne Trace Gas Sounding with the Terahertz Channel of TELIS, Remote Sens., 10, 315, https://doi.org/10.3390/rs10020315, 2018. a, b, c

Yamada, T., Sato, T. O., Adachi, T., Winkler, H., Kuribayashi, K., Larsson, R., Yoshida, N., Takahashi, Y., Sato, M., Chen, A. B., Hsu, R. R., Nakano, Y., Fujinawa, T., Nara, S., Uchiyama, Y., and Kasai, Y.: HO Generation Above Sprite-Producing Thunderstorms Derived from Low-Noise SMILES Observation Spectra, Geophys. Res. Lett., 47, e60090, https://doi.org/10.1029/2019GL085529, 2020. a

Download
Short summary
In the atmosphere, more than 80 % of chlorine compounds are anthropogenic. Hydrogen chloride (HCl), the main stratospheric chlorine reservoir, is useful to estimate the total budget of the atmospheric chlorine compounds. We report, for the first time, the HCl vertical distribution from the middle troposphere to the lower thermosphere using a high-sensitivity SMILES measurement; the data quality is quantified by comparisons with other measurements and via theoretical error analysis.