The Influence of Instrumental Line Shape Degradation on the Partial Columns of O 3 , CO , CH 4 and N 2 O Derived from High-Resolution FTIR Spectrometry

High resolution Fourier transform infrared (FTIR) measurement of direct sunlight does not only provide information of trace gas total columns, but also vertical distribution. Measured O3, CO, CH4, and N2O can be separated into multiple partial columns using the optimal estimation method (OEM). The retrieval of trace gas profiles is sensitive to the instrument line shape (ILS) of the FTIR spectrometer. In this paper, we present an investigation of the influence of ILS degradation on the partial column retrieval of O3, CO, CH4, and N2O. Sensitivities of the partial column, error, and degrees of freedom (DOFs) of each layer to different levels of ILS degradation for O3, CO, CH4, and N2O are estimated. We then evaluate the impact of ILS degradation on the long-term measurements. In addition, we derive the range of ILS degradation corresponding to the acceptable uncertainties of O3, CO, CH4, and N2O results. The results show that the uncertainties induced by the ILS degradation on the absolute value, error, and the DOFs of the partial column are altitude and gas species dependent. The uncertainties of the partial columns of O3 and CO are larger than those on CH4 and N2O. The stratospheric partial columns are more sensitive to the ILS degradation compared to the tropospheric part. Our result improves the understanding of the ILS degradation on the FTIR measurements, which is important for the quantification of the measurement uncertainties and minimizes the bias of the inter-comparison between different measurement platforms. This is especially useful for the validation of satellite observations, the data assimilation of chemical model simulations, and the quantification of the source/sink/trend from the FTIR measurements.

Within the NDACC network, the solar spectra are acquired with similar instruments and are then processed with dedicated algorithms [11], which ensures consistent results between different FTIR (Fourier transform infrared) sites.However, the instrumental line shape (ILS) of the FTIR spectrometers may drift slowly due to mechanical degradation over time or may change abruptly because of operator intervention [12].Moreover, Sun et al. (2017) found that the ILS status is dependent on the optical settings, because the mechanical errors between different field stops may be inconsistent [13].The routine observation may change the entrance field stop size if incident radiation changes, which may introduce inconsistency into ILS.All of these misalignments will result in biases if not properly characterized.
Depending on the degrees of freedom (DOFs) of the retrievals, the total columns of certain NDACC gases, like O 3 , CH 4 , CO, and N 2 O, can be divided into multiple independent partial columns [8,9,14].In certain circumstances, such as validation of chemical model simulations or satellite observations or source/sink/trend estimations, the partial column is more useful than the total column because it is only integrated over the most relevant layers and eliminates the influence from the layers below and/or above [8].Vigouroux et al. (2008Vigouroux et al. ( , 2015) ) divided the O 3 measurements at eight NDACC FTIR stations, namely, Ny-Ålesund (79 • N), Thule (77 • N), Kiruna (68 • N), Harestua (60 • N), Jungfraujoch (47 • N), Izaña (28 • N), Wollongong (34 • S), and Lauder (45 • S), into four independent partial columns, one in the troposphere and three in the stratosphere up to about 45 km.These O 3 partial columns were validated with the coincident Sondes or LIDAR data and were then used for the investigation of O 3 trends and variabilities [8,15].Zhou et al. (2018) studied the atmospheric CO and CH 4 variability with the partial column times series measured with FTIR spectrometers at two sites (St Denis and Maïdo) on Reunion Island (21 • S) in the Indian Ocean.Meanwhile, these partial column times series were compared to the in situ measurements, the GEOS-Chem model simulations, and the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) satellite measurements [16].García et al., (2018) used the NDACC FTIR measurements made between 2007 and 2017 in subtropical Izaña (28 • N), the mid-latitude station of Karlsruhe (49 • N) and the Kiruna polar site (68 • N) for the validation of the IASI (Infrared Atmospheric Sounding Interferometer) CH 4 and N 2 O partial column products [17].Sun et al (2018) investigated the influence of ILS degradation on the total column of all current NDACC mandatory gases and deduced the range of ILS deviations that can be acceptable [18].However, these deductions do not apply to the partial column that is integrated over each independent layer because the sensitivity of the profiles of all NDACC gases to ILS degradation is not constant at different altitudes (Figures S1 and S2).This altitude-dependent sensitivity influences the partial columns differently.Accurate knowledge of this influence is beneficial to minimize the time-to-time and site-to-site biases, which can minimize the propagation error in either the validation of chemical model simulations or satellite observations or source/sink/trend estimations [19].Currently, the influence of ILS degradation on the partial columns of the NDACC gases is still not fully quantified.
This paper investigates the influence of ILS degradation on the partial columns of O 3 , CO, CH 4 , and N 2 O derived from high-resolution FTIR spectrometry and deduces the range of the ILS degradation corresponding to the acceptable uncertainties of O 3 , CO, CH 4 , and N 2 O results.

The Influence of ILS Degradation
The measured spectra were the Fourier transform of the interferogram recorded by the FTIR instrument.The ILS is the Fourier transform of the weighting applied to the interferogram.This weighting consists of two parts: an artificially applied part to change the spectrum and an unavoidable part which is due to the fact that the interferogram is finite in length (box car function), the divergence of the beam is non-zero (due to the non-zero entrance aperture), and several other effects which are due to misalignment ( [20], Chapter 9).We used the latest version of the SFIT4 (version 0.9.4.4) algorithm to perform profile retrieval [1].The basic principle of SFIT4 is the use of an Optimal Estimation Method (OEM) to fit the calculated-to-measured spectra with an iterative Gauss-Newton scheme [21].The SFIT4 includes a nonlinear least-squares spectral fitting subroutine that iteratively generates forward-modelled spectra until the best fit to the measured spectrum is achieved.This iteration process can be expressed as where y C n and y M n are the calculated and measured spectra at the n spectral point, respectively; y n top is the atmospheric top layer spectrum; τ(ν n ) is the atmospheric transmittance spectrum; C is the continuum level; S is the continuum tilt across a window; ν is the spectral frequency drift; ν 0 is the center frequency of the chosen spectral window; ILS(ν n , ν) is the instrument line shape function; z offset is the offset at zero path difference; m is the number of spectral sampling points; and y M n − y C n is called fitting residuals.The τ(ν n ) is related to the absorption cross section α k (ν n , z ) by where k refers to the k-th absorber and ρ k (z ) is its concentration.According to Rodgers (2000), the optimal estimation solution for profile x, the DOF d s , and the total error T err with the iterative Gauss−Newton scheme can be expressed as where F(x i ) is the forward model calculation; S a is the a priori covariance matrix; x a is the a priori state vector; S ε is the measurement noise error covariance matrix; K i is the Jacobian matrix which links the measurement vector y to the state vector x i : ∆y = K i ∆x i ; E s is the smoothing error calculated via Equation ( 7); E m is the measurement error calculated via Equation (8); and E b is the forward model parameter error calculated via Equation ( 9) ( [21], Chapter 3).
For a perfectly aligned spectrometer, the ILS might be close to the theoretical limit, and the measured spectra can be well reproduced by the synthetic spectra calculated with an ideal ILS.However, if an FTIR spectrometer is subject to misalignment, the resulting ILS degradation could affect the calculated spectra and therefore cause the retrieval to reach convergence with the wrong atmospheric status.The x a and S a are independent of the ILS because they are a priori parameters.
The S ε and S b depend on the ILS, since the ILS is part of the forward model calculation and can introduce correlations among the different spectral points, causing non-diagonal matrices of S ε and S b .The Jacobians also depend on the ILS because they represent the sensitivity of the spectra to the atmospheric status, and the analyzed spectra depend on the ILS.
If we assume that the number of iterations is j, and S ε , S b , and K b turn into S ε , S b , and K b , respectively, when using an error ILS, the resulting biases for x, the DOFs, and the total error can be expressed as where E s , E m , and E b are calculated via Equations ( 7)-( 9) using S ε , S b , and K b .In the partial profiles of the DOFs and error calculation, the elements of the corresponding matrices are set to zero for the altitudes outside of the concerned partial layer boundaries.The partial column between altitudes z 1 and z 2 is calculated by integrating profile x by Equation (14), where A m is the air-mass profile derived from the FTIR retrievals:

The Influence of Reference Selection
In this study, the fractional difference (%D) is defined by Equation ( 15): where X represents either the partial column, partial DOFs, or partial error deduced with a degraded ILS.X ref is the same as X but is deduced with a nominal ideal ILS.We used the results with an ideal ILS as the reference.If the ILS of the actual measurements deviates from ideal conditions, it will cause an offset to both X and X ref .
Assuming the offset is ∆, then Equation (15) turns into Equation (16).The NDACC network regularly performs cell measurements to diagnose misalignment of the spectrometer and to realign the instrument when indicated [11].For a perfectly aligned spectrometer, (X ref -∆) is close to X ref , and thus, the influence of ideal condition selection as the reference is of secondary significance (%D ≈ %D ).

Experimental Scheme
We used the program ALIGN60 to simulate ILS degradation in a high resolution FTIR spectrometer.The ALIGN60 is a raytracing model for FTIR spectrometers following the classical Michelson design, assuming one fixed and one movable arm, and using cube corners instead of plane mirrors.It calculates the resulting phase distortions in the recombined beam and from these deduces the variable intensity observed by the detector.The ALIGN60 takes into account the lateral shear error of the movable retro-reflector as a function of the OPD (optical path difference), a decenter of the field stop with respect to the optical axis, an unsharp boundary line, or deformation of the field stop image (possibly caused by a defocused collimator), and vignetting effects with increasing OPD.It can generate trustworthy results with respect to all types of misalignment [22][23][24].The ILS degradation can be regarded as positive modulation efficiency (ME) degradation if its amplitude is larger than unity or negative ME degradation if it is less than unity [22][23][24].Sun et al. (2018) compared the influences of typical types of ILS degradation (denoted as ILS a to j in Table 1 in Sun et al. ( 2018)) on the total columns of all NDACC gases, and found that, for positive ME, the ILS j generally causes the maximal influence, and for negative ME, the ILS i causes the maximal influence [18].The ILS i, as a consequence of cosine and sine bending of the moveable arm, causes a chordal increase in misalignment within the first half of OPD and causes a chordal decrease in misalignment within the second half of OPD.The ILS j, as a consequence of constant shear plus cosine and sine bending of the movable arm, causes a chordal decrease in misalignment within the first half of OPD and causes a chordal increase in misalignment within the second half of OPD.This study selected ILS i and j to represent negative and positive ME, respectively, to investigate the influence of ILS degradation on partial columns of O 3 , CO, CH 4 , and N 2 O.By perturbing the input parameters included into the ALIGN60, we simulated seven levels of ILS i (ranging from −4% to −28% as shown in Figure S3) and ILS j (ranging from 4% to 28% as shown in Figure S4).We then incorporated them into the SFIT forward model in sequence and calculated the fractional difference (D%) in the partial column, partial error, and partial DOFs of each gas relative to the retrieval with an ideal ILS.We first performed a case study by selecting only one typical spectrum for each gas, and then we evaluated the consistency of the resulting deduction with one year of measurements.Table 1 reports the characteristics of the spectra used in the case study.The spectra were randomly selected from routine measurements in clear sky on February 16, 2016.The spectra used in the consistency evaluation were from August 2015 to August 2016, and spanned a large range of gas concentrations, solar zenith angles (SZAs), atmospheric humidity, surface pressures, temperatures, wind speeds, and wind directions (Figure S5).All actual ILS degradations of the FTIR spectrometer within the selected period were less than 2% and can be regarded as ideal.The tropopause height for Hefei derived from the National Centers for Environmental Protection (NCEP) database was 15.1 km with a standard deviation (1σ) of 1.1 km.

Retrieval Strategy
The retrieval settings for O 3 , CO, CH 4 , and N 2 O followed the NDACC standard conventions [25].All spectroscopic line parameters were adopted from HITRAN 2008 [1].A priori profiles of all gases except H 2 O were from a dedicated WACCM (Whole Atmosphere Community Climate Model) run.A priori profiles of pressure, temperature and H 2 O were interpolated from the National Centers for Environmental Protection and National Center for Atmospheric Research (NCEP/NCAR) reanalysis [25].We assumed S ε to be diagonal, and set its diagonal elements to the inverse square of the signal-to-noise ratio (SNR) of the fitted spectra and its non-diagonal elements to zero.For all gases, the diagonal elements of S a were set to the standard deviation of a dedicated WACCM run from 1980 to 2020, and its non-diagonal elements were set to zero.
We defined the layers by summing up the diagonal entries of the averaging kernel matrix from the surface up to higher altitudes.Whenever 1.0 was reached a new layer was started; i.e., for each layer the DOFs were at least 1.0.This has been used in several previous studies in order to get a product that is mainly determined by the measurements and only weakly affected by the a priori information [8,15,19].Table 2 summarizes typical total DOFs and the sensitivity ranges (the sensitivity is larger than 0.5) of the retrieved O 3 , CO, CH 4 , and N 2 O profiles in Hefei, China.The partial DOFs and altitude ranges obtained at each partial layer are also added in Table 2.The retrieved profiles of O 3 , CO, CH 4 , and N 2 O were separated into four, three, three, and four partial layers, respectively.The thus defined layers were independent on the basis of the resolution of the averaging kernels, as can be seen in Figure 1, where the partial column averaging kernels (PAVKs) of the four gases are given.Their PAVKs were resolved at their FWHM (Full Width at Half Maximum), where the averaging kernels peaked at the right altitude; i.e., at the middle of the chosen layer ranges.  of the FTIR spectrometer within the selected period were less than 2% and can be regarded as ideal.
The tropopause height for Hefei derived from the National Centers for Environmental Protection (NCEP) database was 15.1 km with a standard deviation (1σ) of 1.1 km.

Retrieval Strategy
The retrieval settings for O3, CO, CH4, and N2O followed the NDACC standard conventions [25].All spectroscopic line parameters were adopted from HITRAN 2008 [1].A priori profiles of all gases except H2O were from a dedicated WACCM (Whole Atmosphere Community Climate Model) run.A priori profiles of pressure, temperature and H2O were interpolated from the National Centers for Environmental Protection and National Center for Atmospheric Research (NCEP/NCAR) reanalysis [25].We assumed Sε to be diagonal, and set its diagonal elements to the inverse square of the signalto-noise ratio (SNR) of the fitted spectra and its non-diagonal elements to zero.For all gases, the diagonal elements of Sa were set to the standard deviation of a dedicated WACCM run from 1980 to 2020, and its non-diagonal elements were set to zero.
We defined the layers by summing up the diagonal entries of the averaging kernel matrix from the surface up to higher altitudes.Whenever 1.0 was reached a new layer was started; i.e., for each layer the DOFs were at least 1.0.This has been used in several previous studies in order to get a product that is mainly determined by the measurements and only weakly affected by the a priori information [8,15,19].Table 2 summarizes typical total DOFs and the sensitivity ranges (the sensitivity is larger than 0.5) of the retrieved O3, CO, CH4, and N2O profiles in Hefei, China.The partial DOFs and altitude ranges obtained at each partial layer are also added in Table 2.The retrieved profiles of O3, CO, CH4, and N2O were separated into four, three, three, and four partial layers, respectively.The thus defined layers were independent on the basis of the resolution of the averaging kernels, as can be seen in Figure 1, where the partial column averaging kernels (PAVKs) of the four gases are given.Their PAVKs were resolved at their FWHM (Full Width at Half Maximum), where the averaging kernels peaked at the right altitude; i.e., at the middle of the chosen layer ranges.

Sensitivity Study
Figures 2-5 show the sensitivities of the partial column, partial error, and partial DOFs at each partial layer for O 3 , CO, CH 4 , and N 2 O with respect to different levels of ILS degradation, respectively.The partial column, partial error, and partial DOFs were calculated via Equations ( 5), (6), and ( 14), respectively, by setting the elements of the corresponding matrices to zero for the altitudes outside of the concerned partial column boundaries.The fractional differences in the partial column, partial error, and partial DOFs were calculated via Equation (15).

Sensitivity Study
Figures 2, 3, 4 and 5 show the sensitivities of the partial column, partial error, and partial DOFs at each partial layer for O3, CO, CH4, and N2O with respect to different levels of ILS degradation, respectively.The partial column, partial error, and partial DOFs were calculated via Equations ( 14), ( 6) and (5), respectively, by setting the elements of the corresponding matrices to zero for the altitudes outside of the concerned partial column boundaries.The fractional differences in the partial column, partial error, and partial DOFs were calculated via Equation (15).

O 3
PC1 (0-9 km) and most PC2 (9-18 km) are in the troposphere, and PC3 (18-27 km) and PC4 (27-44.4km) are in the stratosphere.Generally, for all partial layers, the partial column is more sensitive to ILS degradation than the partial error and partial DOFs.Positive ME causes an overestimated partial column at PC1, PC2, and PC3 and causes an underestimated partial column at PC4.In contrast, negative ME causes an underestimated partial column at PC1, PC2, and PC3, and causes an overestimated partial column at PC4.For the tropospheric PC1 and PC2, the influence of positive ME on the partial column is generally larger than that of negative ME, while for the stratospheric PC3 and PC4, the influence of negative ME on the partial column is generally larger than that of positive ME.

CO
PC1 (0-3.5 km) and PC2 (3.5-12 km) are in the troposphere, and most PC3 (12-68 km) is in the stratosphere.Generally, for all partial layers, the partial column, partial error and partial DOFs are more sensitive to positive ME degradation than negative ME degradation.Positive ME causes an overestimated partial column at PC1 and PC3 and causes an underestimated partial column at PC2.In contrast, negative ME causes an underestimated partial column at PC1 and PC3 and causes an overestimated partial column at PC2.The influence of ILS degradation on the partial error at each partial layer is opposite to that on the partial column.For all partial layers, positive ME overestimated partial DOFs and negative ME underestimated partial DOFs.

O3
PC1 (0-9 km) and most PC2 (9-18 km) are in the troposphere, and PC3 (18-27 km) and PC4 (27-44.4km) are in the stratosphere.Generally, for all partial layers, the partial column is more sensitive to ILS degradation than the partial error and partial DOFs.Positive ME causes an overestimated partial column at PC1, PC2, and PC3 and causes an underestimated partial column at PC4.In contrast, negative ME causes an underestimated partial column at PC1, PC2, and PC3, and causes an overestimated partial column at PC4.For the tropospheric PC1 and PC2, the influence of positive ME on the partial column is generally larger than that of negative ME, while for the stratospheric PC3 and PC4, the influence of negative ME on the partial column is generally larger than that of positive ME.

CO
PC1 (0-3.5 km) and PC2 (3.5-12 km) are in the troposphere, and most PC3 (12-68 km) is in the stratosphere.Generally, for all partial layers, the partial column, partial error and partial DOFs are more sensitive to positive ME degradation than negative ME degradation.Positive ME causes an overestimated partial column at PC1 and PC3 and causes an underestimated partial column at PC2.In contrast, negative ME causes an underestimated partial column at PC1 and PC3 and causes an overestimated partial column at PC2.The influence of ILS degradation on the partial error at each partial layer is opposite to that on the partial column.For all partial layers, positive ME overestimated partial DOFs and negative ME underestimated partial DOFs.

N 2 O
PC1 (0-5 km) and most PC2 (5-11.5 km) are in the troposphere, PC3 (11.5-20 km) is a mixture of troposphere and stratosphere, and PC4 (20-35 km) is in the stratosphere.The partial columns at PC3 and PC4 are more sensitive to ILS degradation than those at PC1 and PC2.For PC1, PC2, and PC3, the influence of positive ME on the partial column is generally larger than that of negative ME.At PC4, the positive and negative ME show comparable influences on the partial column.Positive ME causes an overestimated partial column at PC1 and causes an underestimated partial column at all other partial layers.The influence of negative ME on the partial columns at PC1 and PC2 is typically less than 1%, regardless of ILS degradation, while negative ME degradation underestimates the partial column at PC3 and overestimates the partial column at PC4.The changes in partial error and partial DOFs do not vary monotonously with ILS degradation.

Consistency Evaluation
The simulated ILS with an ME amplitude of −4% and ILS j with an ME amplitude of 4% were incorporated in the spectra from August 2015 to August 2016 and then the results were compared to those deduced with an ideal ILS.The variability (mean ± standard deviation) of fractional differences in partial column, partial error, and partial DOFs calculated with −4% of ILS i and 4% of ILS j were used to evaluate the consistency of the influences of negative and positive ME degradations.
With respect to positive ME degradation, Figures S6-S9 show the variability of fractional differences in the partial column, partial error, and partial DOFs for O 3 , CO, CH 4 , and N 2 O as a function of SZA, respectively.Figures S10-S13 are the same as Figures S6-S9 but for negative ME degradation.The resulting statistics are summarized in Table 3.Generally, for both positive and negative ME degradation, the variability of the partial column, partial error, and partial DOFs at each partial layer for O 3 and CO is larger than that for CH 4 and N 2 O.For O 3 , the variability of the partial column, partial error, and partial DOFs at each partial layer are comparable.The variability of these quantities at PC1 and PC2 is larger than that at PC3 and PC4.With 4% of ILS degradation, the variability of fractional difference in the partial column at PC1, PC2, PC3, and PC4 was shown to be 3.3 ± 13.2%, 5.2 ± 21.1%, 0.7 ± 9.7%, and −3.2 ± 8.5%, respectively.With −4% of ILS degradation, the variability at PC1, PC2, PC3, and PC4 was shown to be 0.92 ± 13.0%, −3.3 ± 23.4%, −1.7 ± 9.4%, and 2.6 ± 9.1%, respectively.
For CO, the partial column and partial error show larger variability at each partial layer than the partial DOFs.Large variability of the partial column for CO is observed at all partial layers.With 4% of ILS degradation, the variability of the fractional difference in the partial column at PC1, PC2, and PC3 was shown to be −3.9 ± 41.2%, 3.7 ± 23.5%, and 3.4 ± 37.8%, respectively.With −4% of ILS degradation, the variability at PC1, PC2, and PC3 was shown to be −4.2 ± 38.7%, 3.3 ± 25.7%, and 6.5 ± 36.5%, respectively.

Discussion and Recommendations
The influence of ILS degradation on the partial quantities of O 3 , CO, CH 4 , and N 2 O is generally larger than that on the total quantities presented in Sun et al. (2018) [18].The influence of ILS degradation on the partial quantities of all gases is altitude and gas species dependent, and the level of the influence varies across ILS degradation levels.Positive and negative ILS degradation generally have opposite influences on partial quantities.The influence of ILS degradation on partial quantities of O 3 and CO is larger than that on CH 4 and N 2 O.For all gases, the partial columns in the troposphere are more sensitive to positive ILS degradation than negative ILS degradation.Generally, the partial columns in the stratosphere are more sensitive to ILS degradation than those in the troposphere.
The sensitivity curves of O 3 , CO, and CH 4 are smoother than that of N 2 O, because ILS degradation alters the S ε , S b , and Jacobian matrices of these gases in different ways.For O 3 , CO, and CH 4 , ILS degradation does not alter the number of iterations, and the sensitivity varies monotonously across ILS degradation levels, while for N 2 O, the number of iterations is occasionally altered, and the sensitivity is also altered.
The fact that the influence of ILS degradation on the partial quantities is altitude and gas species dependent can be explained by Equation ( 17) which shows how the ILS affects the calculated spectra in the forward model calculation.
where S k , L k , and M k refer to the absorption intensity, the ideal line shape, and the molecular mass for the k-th absorber, respectively.Typically, L k can be expressed by a Voigt function as a consequence of translational effects and collisional effects.P(z ) and T(z ) refer to the pressure and temperature at altitude z , respectively.The absorption intensity, molecular mass, and fitting frequency are gas dependent, and the line shape, pressure, and temperature are altitude dependent.The influence of ILS degradation on partial quantities is thus altitude and gas species dependent.Specifically, the absorption line shape in the stratosphere is narrower than that in the troposphere, and stratospheric quantities are more sensitive to ILS degradation.For all gases, we set the acceptable fractional difference in the partial column as 2% and assumed a proportional variability at each partial layer.The permitted maximum ILS degradation for each gas at each partial layer is recommended in Table 4 as: 1.
For O 3 and CO, a precise ILS should be used because the variability in ILS degradation influence at each partial layer is very large.We recommend the incorporation of the measured ILS in partial column retrieval.One option is to use the output of routine cell measurements with LINEFIT code [11].

Conclusions
We assessed the influence of instrumental line shape (ILS) degradation on the partial columns of O 3 , CO, CH 4 , and N 2 O derived from high-resolution FTIR spectrometry.The sensitivities of the partial column, partial error, and partial DOFs (degrees of freedom) with respect to different levels of ILS degradation for O 3 , CO, CH 4 , and N 2 O were first investigated, and then the consistency of the resulting deductions was evaluated with one year of measurements.Finally, the maximum ILS deviations allowable for limiting the influence within 2% were deduced.
The influence of ILS degradation on partial quantities is altitude and gas species dependent.The influence of ILS degradation on the partial columns of O 3 and CO is larger than that for CH 4 and N 2 O.For all gases, the partial columns in the stratosphere are more sensitive to ILS degradation than those in the troposphere.In order to limit the fractional difference in the partial column within 2%, it is recommended that a precise ILS deduced from cell measurements is used for O 3 and CO; for CH 4 , the influence of positive ME (modulation efficiency) degradation at PC2 (7-16 km) and negative ME degradation at PC1 (0-7 km) can be regarded as negligible.The maximum positive ME degradation at PC1 (0-7 km) and PC3 (16-37 km) should be less than 13% and 1.5%, respectively.The maximum negative ME degradation at PC2 (7-16 km) and PC3 (16-37 km) should be less than 12% and 4.5%, respectively; for N 2 O, the maximum positive ME degradation at PC1 (0-5 km), PC2 (5-11.5 km), PC3 (11.5-20 km) and PC4 (20-35 km) should be less than 6%, 4%, 1.5%, and 1.5%, respectively.The influence of negative ME degradation at PC1 (0-5 km) and PC2 (5-11.5 km) can be regarded as negligible.The maximum negative ME degradation at PC3 (11.5-20 km) and PC4 (20-35 km) should be less than 2% and 1.5%, respectively.
Our results provide improved knowledge of the ILS degradation in ground-based FTIR measurements, which is important for the quantification of measurement uncertainties and minimization of the bias of inter-comparisons between different measurement platforms.This is especially useful for the validation of satellite observations, the data assimilation of chemical model simulations, and the quantification of the source/sink/trend from FTS measurements.

Figure 2 .
Figure 2. Sensitivity of the partial column, partial error, and partial degrees of freedom (DOFs) at each partial layer for O3 with respect to the modulation efficiency (ME) deviation.The solid red line (P_PC), solid blue line (P_Perr), and solid black line (P_Pdofs) represent the sensitivity of the partial column, partial error, and partial DOFs with respect to positive ME deviation, respectively.The dashed red line (N_PC), dashed blue line (N_Perr), and dashed black line (N_Pdofs) represent the sensitivity of the partial column, partial error, and partial DOFs with respect to negative ME deviation, respectively.The ME deviations were modelled by ALIGN60.The results were deduced from the spectra recorded at Hefei on 16 February 2016.

Figure 2 .
Figure 2. Sensitivity of the partial column, partial error, and partial degrees of freedom (DOFs) at each partial layer for O 3 with respect to the modulation efficiency (ME) deviation.The solid red line (P_PC), solid blue line (P_Perr), and solid black line (P_Pdofs) represent the sensitivity of the partial column, partial error, and partial DOFs with respect to positive ME deviation, respectively.The dashed red line (N_PC), dashed blue line (N_Perr), and dashed black line (N_Pdofs) represent the sensitivity of the partial column, partial error, and partial DOFs with respect to negative ME deviation, respectively.The ME deviations were modelled by ALIGN60.The results were deduced from the spectra recorded at Hefei on 16 February 2016.

Figure 3 .
Figure 3.The same as Figure 2 but for CO.

Figure 4 .
Figure 4.The same as Figure 2 but for CH4.

Figure 3 . 15 Figure 3 .
Figure 3.The same as Figure 2 but for CO.

Figure 4 .
Figure 4.The same as Figure 2 but for CH4.Figure 4. The same as Figure 2 but for CH 4 .

Figure 4 .
Figure 4.The same as Figure 2 but for CH4.Figure 4. The same as Figure 2 but for CH 4 .

Figure 5 .
Figure 5.The same as Figure 2 but for N2O.

Figure 5 .
Figure 5.The same as Figure 2 but for N 2 O.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2072-4292/10/12/2041/s1,FigureS1: Sensitivity of profile with respect to positive modulation efficiency (ME) deviation for current NDACC mandatory gases; FigureS2: The same as Figure S1 but for negative ME deviation; Figure S3: Simulated negative ME deviations along with optical path difference (OPD); Figure S4: Simulated positive ME deviations along with OPD; Figure S5: The meteorological data and solar zenith angles (SZAs) record at Hefei, China; Figure S6: Variability of fractional difference in partial column, partial error, and partial DOFs of O 3 at each partial layer as function of SZA from August 2015 to August 2016 where ILS j with a maximum ME deviation of 4% is used; Figure S7: The same as Figure S6 but for CO; Figure S8: The same as Figure S6 but for CH 4 ; Figure S9: The same as Figure S6 but for N 2 O; Figure S10: The same as S6 but for a maximum ME deviation of −4%; Figure S11: The same as S7 but for a maximum ME deviation of −4%; Figure S12: The same as S8 but for a maximum ME deviation of −4%; Figure S13: The same as S9 but for a maximum ME deviation of−4%.

Table 1 .
Description of the spectra selected for the case study.

Table 2 .
Typical degrees of freedom (DOFs) for the signal and sensitivity ranges of the retrieved O 3 , CO, CH 4 , and N 2 O profiles in Hefei, China.

Table 3 .
The variability (mean ± standard deviation) of fractional differences in the partial column, partial error, and partial DOFs for O 3 , CO, CH 4 , and N 2 O as a function of the solar zenith angle (SZA).

Table 4 .
Recommendations for limiting the difference in the partial columns of all gases within 2%.
* The influence is less than 2%.Even an excessively large instrumental line shape (ILS) degradation of 20% can normally be regarded as being negligible.