Advanced Search
Article Contents

A Differential Optical Absorption Spectroscopy Method for X CO2 Retrieval from Ground-Based Fourier Transform Spectrometers Measurements of the Direct Solar Beam


doi: 10.1007/s00376-015-4213-9

  • A differential optical absorption spectroscopy (DOAS)-like algorithm is developed to retrieve the column-averaged dry-air mole fraction of carbon dioxide from ground-based hyper-spectral measurements of the direct solar beam. Different to the spectral fitting method, which minimizes the difference between the observed and simulated spectra, the ratios of multiple channel-pairs——one weak and one strong absorption channel——are used to retrieve X CO2 from measurements of the shortwave infrared (SWIR) band. Based on sensitivity tests, a super channel-pair is carefully selected to reduce the effects of solar lines, water vapor, air temperature, pressure, instrument noise, and frequency shift on retrieval errors. The new algorithm reduces computational cost and the retrievals are less sensitive to temperature and H2O uncertainty than the spectral fitting method. Multi-day Total Carbon Column Observing Network (TCCON) measurements under clear-sky conditions at two sites (Tsukuba and Bremen) are used to derive X CO2 for the algorithm evaluation and validation. The DOAS-like results agree very well with those of the TCCON algorithm after correction of an airmass-dependent bias.
  • 加载中
  • Bovensmann H., J. P. Burrows, M. Buchwitz, J. Frerick, S. No\"el, V. V. Rozanov, K. V. Chance, and A. P. H. Goede, 1999: SCIAMACHY: Mission objectives and measurement modes. J. Atmos. Sci., 56, 127- 150.
    Buchwitz M., Coauthors, 2005: Carbon monoxide, methane and carbon dioxide columns retrieved from SCIAMACHY by WFM-DOAS: Year 2003 initial data set. Atmospheric Chemistry and Physics, 5, 3313- 3329.
    Canadell, J. G., Coauthors, 2010: Interactions of the carbon cycle, human activity, and the climate system: A research portfolio. Current Opinion in Environmental Sustainability, 2, 301- 311.
    Kurucz R. L., 1995: The solar spectrum: Atlases and line identifications. Workshop on Laboratory and Astronomical High Resolution Spectra, ASP Comference Series, No. 81, A. J. Sauval, R. Blomme and N. Grevesse, 17- 31.
    Deutscher N.M., Coauthors, 2010: Total column CO2 measurements at Darwin, Australia——site description and calibration against in situ aircraft profiles. Atmospheric Measurement Techniques, 3, 947- 958.
    Easterling D. R., M. F. Wehner, 2009: Is the climate warming or cooling? Geophys. Res. Lett., 36, L08706, doi: 10.1029/ 2009GL037810.
    Kuze A., H. Suto, M. Nakajima, and T. Hamazaki, 2009: Thermal and near infrared sensor for carbon observation Fourier-transform spectrometer on the Greenhouse Gases Observing Satellite for greenhouse gases monitoring. Appl. Opt., 48, 6716- 6733.
    Liu Y., D. X. Yang, and Z. N. Cai, 2013: A retrieval algorithm for TanSat XCO2 observation: Retrieval experiments using GOSAT data. Chinese Science Bulletin, 58, 1520- 1523.
    Min Q.-L., M. Z. Duan, 2005: Simultaneously retrieving cloud optical depth and effective radius for optically thin clouds. J. Geophys. Res., 110,D21201, doi: 10.1029/2005JD 006136.
    Min Q. L., E. Joseph, and M. Z. Duan, 2004: Retrievals of thin cloud optical depth from a multifilter rotating shadowband radiometer. J. Geophys. Res., 109,D02201, doi: 10.1029/2003 JD003964.
    Morino I., Coauthors, 2011: Preliminary validation of column-averaged volume mixing ratios of carbon dioxide and methane retrieved from GOSAT short-wavelength infrared spectra. Atmospheric Measurement Techniques, 4, 1061- 1076.
    Reuter M., Coauthors, 2011: Retrieval of atmospheric CO2 with enhanced accuracy and precision from SCIAMACHY: Validation with FTS measurements and comparison with model results .J. Geophys. Res., 116,D04301, doi: 10.1029/ 2010JD015047.
    Solomon S., Coauthors, 2007: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press,996 pp.
    Stephens B.B., Coauthors, 2007: Weak northern and strong tropical land carbon uptake from vertical profiles of atmospheric CO2. Science, 316, 1732- 1735.
    Wunch D., Coauthors, 2011a: The total carbon column observing network. Philosophical Transactions of the Royal Society A Mathematical Physical and Engineering Sciences, 369, 2087- 2112.
    Wunch D., Coauthors, 2011b: A method for evaluating bias in global measurements of CO2 total columns from space. Atmospheric Chemistry and Physics, 11, 12 317- 12 337.
    Yoshida Y., Y. Ota, N. Eguchi, N. Kikuchi, K. Nobuta, H. Tran, I. Morino, and T. Yokota, 2011: Retrieval algorithm for CO2 and CH4 column abundances from short-wavelength infrared spectral observations by the Greenhouse gases observing satellite. Atmospheric Measurement Techniques, 4, 717- 734.
    Yoshida Y., Coauthors, 2013: Improvement of the retrieval algorithm for GOSAT SWIR XCO2 and XCH4 and their validation using TCCON data. Atmospheric Measurement Techniques, 6, 1533- 1547.
  • [1] Minqiang ZHOU, Qichen NI, Zhaonan CAI, Bavo LANGEROCK, Jingyi JIANG, Ke CHE, Jiaxin WANG, Weidong NAN, Yi LIU, Pucai WANG, 2023: Ground-Based Atmospheric CO2, CH4, and CO Column Measurements at Golmud in the Qinghai-Tibetan Plateau and Comparisons with TROPOMI/S5P Satellite Observations, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 223-234.  doi: 10.1007/s00376-022-2116-0
    [2] WU Xue*, WANG Xin, and LÜ Daren, 2014: Retrieval of Vertical Distribution of Tropospheric Refractivity through Ground-Based GPS Observation, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 37-47.  doi: 10.1007/s00376-013-2215-z
    [3] Denghui JI, Minqiang ZHOU, Pucai WANG, Yang YANG, Ting WANG, Xiaoyu SUN, Christian HERMANS, Bo YAO, Gengchen WANG, 2020: Deriving Temporal and Vertical Distributions of Methane in Xianghe Using Ground-based Fourier Transform Infrared and Gas-analyzer Measurements, ADVANCES IN ATMOSPHERIC SCIENCES, 37, 597-607.  doi: 10.1007/s00376-020-9233-4
    [4] Denghui JI, Minqiang ZHOU, Pucai WANG, Yang YANG, Ting WANG, Xiaoyu SUN, Christian HERMANS, Bo YAO, Gengchen WANG, 2024: Erratum to: Deriving Temporal and Vertical Distributions of Methane in Xianghe Using Ground-based Fourier Transform Infrared and Gas-analyzer Measurements, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 377-377.  doi: 10.1007s00376-023-3009-6
    [5] WANG Yuesi, WANG Yinghong, 2003: Quick Measurement of CH4, CO2 and N2O Emissions from a Short-Plant Ecosystem, ADVANCES IN ATMOSPHERIC SCIENCES, 20, 842-844.  doi: 10.1007/BF02915410
    [6] Bozhen LI, Gen ZHANG, Lingjun XIA, Ping KONG, Mingjin ZHAN, Rui SU, 2020: Spatial and Temporal Distributions of Atmospheric CO2 in East China Based on Data from Three Satellites, ADVANCES IN ATMOSPHERIC SCIENCES, 37, 1323-1337.  doi: 10.1007/s00376-020-0123-6
    [7] Bo LIU, Juan HUO, Daren LYU, Xin WANG, 2021: Assessment of FY-4A and Himawari-8 Cloud Top Height Retrieval through Comparison with Ground-Based Millimeter Radar at Sites in Tibet and Beijing, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 1334-1350.  doi: 10.1007/s00376-021-0337-2
    [8] Lingyun ZHANG, Yanfang SONG, Jialin SHI, Qun SHEN, Deng HU, Qiang GAO, Wei CHEN, Kien-Woh KOW, Chengheng PANG, Nannan SUN, Wei WEI, 2022: Frontiers of CO2 Capture and Utilization (CCU) towards Carbon Neutrality, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1252-1270.  doi: 10.1007/s00376-022-1467-x
    [9] XU Yongfu, LI Yangchun, 2009: Estimates of Anthropogenic CO2 Uptake in a Global Ocean Model, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 265-274.  doi: 10.1007/s00376-009-0265-z
    [10] Xu Yongfu, 1992: The Buffer Capability of the Ocean to Increasing Atmospheric CO2, ADVANCES IN ATMOSPHERIC SCIENCES, 9, 501-510.  doi: 10.1007/BF02677083
    [11] LIU Liping, ZHANG Zhiqiang, YU Danru, YANG Hu, ZHAO Chonghui, ZHONG Lingzhi, 2012: Comparison of Precipitation Observations from a Prototype Space-based Cloud Radar and Ground-based Radars, ADVANCES IN ATMOSPHERIC SCIENCES, 29, 1318-1329.  doi: 10.1007/s00376-012-1233-6
    [12] Xiaotong ZHU, Qingqing LI, Jinhua YU, Dan WU, Kai YAO, 2018: Geometric Characteristics of Tropical Cyclone Eyes before Landfall in South China Based on Ground-Based Radar Observations, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 592-603.  doi: 10.1007/s00376-017-7144-9
    [13] L.S. Hingane, 1989: Effect of Increasing CO2 on the Stratospheric Level of CO and O3, ADVANCES IN ATMOSPHERIC SCIENCES, 6, 390-392.  doi: 10.1007/BF02661544
    [14] Dongxu YANG, Janne HAKKARAINEN, Yi LIU, Iolanda IALONGO, Zhaonan CAI, Johanna TAMMINEN, 2023: Detection of Anthropogenic CO2 Emission Signatures with TanSat CO2 and with Copernicus Sentinel-5 Precursor (S5P) NO2 Measurements: First Results, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 1-5.  doi: 10.1007/s00376-022-2237-5
    [15] B. S. K. REDDY, K. R. KUMAR, G. BALAKRISHNAIAH, K. R. GOPAL, R. R. REDDY, V. SIVAKUMAR, S. Md. ARAFATH, A. P. LINGASWAMY, S. PAVANKUMARI, K. UMADEVI, Y. N. AHAMMED, 2013: Ground-Based In Situ Measurements of Near-Surface Aerosol Mass Concentration over Anantapur: Heterogeneity in Source Impacts, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 235-246.  doi: 10.1007/s00376-012-1234-5
    [16] MIAO Qun, and Bart GEERTS, 2013: Airborne measurements of the impact of ground-based glaciogenic cloud seeding on orographic precipitation, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1025-1038.  doi: 10.1007/s00376-012-2128-2
    [17] Debashis NATH, CHEN Wen, 2013: Investigating the Dominant Source for the Generation of Gravity Waves during Indian Summer Monsoon Using Ground-based Measurements, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 153-166.  doi: 10.1007/s00376-012-1273-y
    [18] HUO Juan, ZHANG Wenxing, ZENG Xiaoxia, Lü Daren, LIU Yi, 2013: Examination of the Quality of GOSAT/CAI Cloud Flag Data over Beijing Using Ground-based Cloud Data, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1526-1534.  doi: 10.1007/s00376-013-2267-0
    [19] Wei Chong, Xue Yongkang, Zhu Xiaoming, Zou Shouxiang, 1984: DETERMINATION OF ATMOSPHERIC PRECIPITABLE WATER AND HUMIDITY PROFILES BY A GROUND-BASED 1.35 cm RADIOMETER, ADVANCES IN ATMOSPHERIC SCIENCES, 1, 119-139.  doi: 10.1007/BF03187623
    [20] Yun HE, Fan YI, Fuchao LIU, Zhenping YIN, Jun ZHOU, 2022: Ice Nucleation of Cirrus Clouds Related to the Transported Dust Layer Observed by Ground-Based Lidars over Wuhan, China, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 2071-2086.  doi: 10.1007/s00376-021-1192-x

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 22 September 2014
Manuscript revised: 08 January 2015
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

A Differential Optical Absorption Spectroscopy Method for X CO2 Retrieval from Ground-Based Fourier Transform Spectrometers Measurements of the Direct Solar Beam

  • 1. Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, College of Atmospheric Sciences, Lanzhou University, Lanzhou 730000
  • 2. Key Laboratory of Middle Atmosphere and Global Environment Observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 10029
  • 3. Wuhan University, Wuhan 430000
  • 4. Atmospheric Science Research Center, State University of New York, Albany, NY 12203, USA

Abstract: A differential optical absorption spectroscopy (DOAS)-like algorithm is developed to retrieve the column-averaged dry-air mole fraction of carbon dioxide from ground-based hyper-spectral measurements of the direct solar beam. Different to the spectral fitting method, which minimizes the difference between the observed and simulated spectra, the ratios of multiple channel-pairs——one weak and one strong absorption channel——are used to retrieve X CO2 from measurements of the shortwave infrared (SWIR) band. Based on sensitivity tests, a super channel-pair is carefully selected to reduce the effects of solar lines, water vapor, air temperature, pressure, instrument noise, and frequency shift on retrieval errors. The new algorithm reduces computational cost and the retrievals are less sensitive to temperature and H2O uncertainty than the spectral fitting method. Multi-day Total Carbon Column Observing Network (TCCON) measurements under clear-sky conditions at two sites (Tsukuba and Bremen) are used to derive X CO2 for the algorithm evaluation and validation. The DOAS-like results agree very well with those of the TCCON algorithm after correction of an airmass-dependent bias.

1. Introduction
  • Carbon dioxide (CO2) is considered to be the main greenhouse gas causing current global warming (Solomon et al., 2007). However, (Easterling and Wehner, 2009) reported that records of surface air temperature show no warming trend or even a slight cooling trend, while greenhouse gas levels are still increasing. The disagreement about climate change is mostly due to the lack of long-term records of CO2 measurements, especially for large area measurements and CO2 sources and sinks (Stephens et al., 2007; Canadell et al., 2010).

    It is advantageous to use satellite remote sensing to monitor atmospheric CO2 globally. However, at present, only the satellite datasets of column-averaged dry-air mole fraction of CO2 (X CO2) from the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) on board the Environmental Satellite (ENVISAT) (Bovensmann et al., 1999) and the Thermal and Near-infrared Sensor for Carbon Observation-Fourier Transform Spectrometer (TANSO-FTS) on board the Greenhouse Gases Observing Satellite (GOSAT) (Kuze et al., 2009), are used to estimate regional CO2 fluxes. Both instruments use the reflected solar radiation in the shortwave infrared (SWIR) spectral region, making them sensitive to the variation of near-surface CO2 concentrations. Unfortunately, the low spectral resolution of SCIAMACHY limits the inversion accuracy, with a single retrieval precision of about 2.5 ppm, as compared to ground-based Fourier transform spectrometer (FTS) measurements (Buchwitz et al., 2005; Reuter et al., 2011). The biases and standard deviations of the column-averaged dry-air mole fraction of carbon dioxide (X CO2) from the SWIR L2 V02.xx GOSAT retrieval algorithm reach -1.48 and 2.09 ppm, respectively (Morino et al., 2011; Yoshida et al., 2011; Yoshida et al., 2013).

    Figure 1.  (a) Temperature and (b) water vapor profiles at 0000 (black), 0600 (blue) and 1200 UTC (red) in one day, and (c) inversion errors with same coefficients a and b.

    Compared with satellite observations of reflected light, ground-based observations of the direct solar beam are less influenced by surface albedo, aerosols etc. Therefore, ground-based observations can achieve higher accuracy and precision in determining the CO2 total column amount. However, at present, the Total Carbon Column Observing Network (TCCON) is the only existing network that retrieves the total column concentration of CO2 from ground-based FTS measurements for satellite validation. TCCON achieves a network-wide uncertainty of X CO2 of better than 0.8 ppm, with after correcting for an airmass-dependent bias and calibrating to aircraft vertical profiles (Wunch et al., 2011a, b).

    A Chinese satellite for CO2 monitoring is planned for launch in 2015 (Liu et al., 2013). To validate the satellite retrievals, a surface observation network has been set up to measure the hyper-spectrum of the direct solar beam in the SWIR bands. To derive the total column amount of CO2 from these spectral measurements, a retrieval algorithm is needed. In this paper, a new DOAS-like algorithm is developed, in which multiple pairs of CO2 absorption ratios (one in the weak CO2 absorption channel and one in the strong CO2 absorption channel) are used to derive the column CO2. More importantly, both channels in the pair are carefully selected to reduce their sensitivity to the surface pressure, air temperature, water vapor, noise and frequency shift. Compared with the spectral fitting method, DOAS-like retrievals are less sensitive to temperature and H2O uncertainty.

    Figure 2.  Channels used in the retrieval algorithm.

    Figure 3.  Error analysis of X CO2 retrieval when the strong absorption channel is located at the (a) far-wing (plus symbol), (b) in the middle (triangle symbol), and (c) center (cross symbol) of the absorption line. Five-hundred retrievals with random noise under the assumption of SNR of 500 are illustrated. Left-hand panels represent the spectrum of a single absorption line, and right-hand panels the errors of X CO2 retrievals.

2. Retrieval algorithm
  • Our retrieval algorithm is based on the fact that the photon path lengths within a narrow spectral range are equal. Therefore, the ratio of the channel pair is proportional to X CO2 if the surface pressure, temperature profile and water vapor are known. Based on the Lambert-Beer law, a ground-based measurement of the direct solar beam for a fixed wavelength can be expressed as \begin{equation} \label{eq1} I_\lambda=I_{{\rm sca},\lambda}+I_{0,\lambda}e^{-\tau m} , \end{equation} where IΛ is the downward radiance measured at the bottom of the atmosphere for wavelength Λ, I sca,Λ is the forward scattering contribution in the incident direction, and I0 is the incoming solar radiance at the top of the atmosphere. m is the air mass factor. τ is the optical depth in the vertical optical path, which can be written as \begin{equation} \label{eq2} \tau=\tau_{{\rm CO}_2}+\tau_{{\rm H}_2{\rm O}}+\tau_{\rm aer}+\tau_{\rm Ray} , \end{equation} where the right-hand terms represent optical depth of CO2 absorption, water vapor absorption, aerosol extinction, and Rayleigh scattering, respectively. The scattering term of I sca,Λ in Eq. (2) is negligible due to a very small field of view (FOV) (2.4 mrad) of the spectrometer, particularly for small aerosol particles and small aerosol optical depths (Min et al., 2004; Min and Duan, 2005; Wunch et al., 2011a). Therefore, the radiance can be simplified as \begin{equation} \label{eq3} I_\lambda=I_0e^{-(\tau_{{\rm CO}_2}+\tau_{{\rm H}_2{\rm O}}+\tau_{\rm aer}+\tau_{\rm Ray})m} . \end{equation} In a very limited spectral range, the variation of τ aer and τ Ray across the spectral range can be ignored. Therefore the ratio of the selected channel pair is insensitive to the loading of aerosol and Rayleigh scattering. Hence, we have \begin{equation} \label{eq4} \dfrac{I_{\lambda_1}}{I_{\lambda_2}}=\dfrac{I_{0,\lambda_1}}{I_{0,\lambda_2}}e^{-[(\tau_{{{\rm CO}_2},\lambda_1}- \tau_{{\rm CO}_2,\lambda_2})-(\tau_{{\rm H}_2{\rm O},\lambda_1}-\tau_{{\rm H}_2{\rm O},\lambda_2})]m} . \end{equation} Letting r=IΛ1/IΛ2 and r0=I0,Λ1/I0,Λ2, Eq. (5) can be rewritten as \begin{equation} \label{eq5} r=r_0e^{-[(\tau_{{\rm CO}_2,\lambda_1}-\tau_{{\rm CO}_2,\lambda_2})-(\tau_{{\rm H}_2{\rm O},\lambda_1}-\tau_{{\rm H}_2{\rm O},\lambda_2})]m} . \end{equation} By taking the logarithm of Eq. (6), we have \begin{equation} \label{eq6} \ln(r)=\ln(r_0)-[(\tau_{{\rm CO}_2,\lambda_1}-\tau_{{\rm CO}_2,\lambda_2})-(\tau_{{\rm H}_2{\rm O},\lambda_1}-\tau_{{\rm H}_2{\rm O},\lambda_2})]m . \end{equation} The optical depth τ CO2 is proportional to the total number of molecules of CO2 per surface area, which is positively correlated to X CO2 if the surface pressure, air temperature and CO2 volume mixing ratio (VMR) profile are assumed to be known. Furthermore, only channel pairs with weak H2O absorption interference are selected. Therefore, the difference associated with water vapor is small and can be treated as a correction coefficient. Then, Eq. (7) is simplified as \begin{equation} \label{eq7} X_{{\rm CO}_2}=a\ln(r)+b . \end{equation} Through the about pair selection procedure, the retrieval, i.e., Eq. (8), is weakly sensitive to the atmospheric state uncertainty (temperature and water vapor). Nonetheless, the coefficients of a and b are weakly dependent on temperature and water vapor in the atmosphere. To further reduce the error associated with the atmospheric state, we can calculate both coefficients with the surface pressure of in-situ measurements collected by automatic meteorological stations and reanalysis/forecasting atmospheric profiles. The profiles can be fixed for multiple measurements within some specific time period because only channels that are independent of temperature and water vapor are used in our retrieval algorithm. To illustrate the feasibility of fixed profiles, several inversions calculated by the different coefficients a and b at 0000, 0600 and 1200 UTC are shown in Fig. 1. All the errors are less than 0.15 ppm. The DOAS-like algorithm of Eq. (8) only has one unknown parameter. Hence, no iteration is needed.

    Figure 4.  The X CO2 errors of each channel for +1 hPa bias of the surface pressure when the channel is used as the only strong absorption channel in the super channel-pair.

    Figure 5.  The X CO2 errors of each channel for a +1 K bias of temperature profile when the channel is used as the only strong absorption channel in the super channel-pair.

    Figure 6.  The strong CO2 absorption channels distributed on different sides of an absorption line (upper panel) and the inversion errors (lower panel) caused by frequency shift when channel 1 (plus symbol), channel 2 (triangle symbol), and the average of channel 1 and channel 2 (cross symbol) are used as the strong absorption channel in the super channel-pair.

    Figure 7.  The prior atmosphere profiles of temperature (left panel), water vapor (middle panel) and CO2 (right panel) used in the sensitivity tests.

  • The Atmospheric Environment Research (AER) Line By Line Radiative Transfer Model (LBLRTM) is used to calculate the gases' absorption coefficients. The line parameters database, aer_v_3.2, is also from AER and a Voigt line shape is assumed. The extra-terrestrial solar spectrum is provided by Kurucz (Kurucz, 1995). The atmosphere is divided into 70 layers in 1 km intervals. The surface pressure and a priori VMR profiles at TCCON stations are from TCCON auxiliary data. In the forward model, LBLRTM is used to calculate the gas absorption optical depth with a step of 0.001 cm-1, and these absorptions along with the Rayleigh scattering are used to calculate the fine structure of direct radiance. Finally, the fine spectra are convolved with a resolution and sample rate of 0.02 cm-1 and 7.5 kHz for measurements in Tsukuba, and 0.014 cm-1 and 7.5 kHz for measurements in Bremen (Wunch et al., 2011a).

    Figure 8.  Left panels: inversion X CO2 of TCCON (black) and the DOAS-like (blue) method at Tsukuba and Bremen. Right panels: the corresponding SZA time series.

  • The DOAS-like method could reduce computational cost, but the super channel-pair must be carefully selected to reduce the impacts of H2O absorption, the solar Fraunhofer lines, and other factors such as instrument noise, temperature, pressure, frequency shift etc. In our channel-pair, the mean of 430 channels with very weak CO2 absorption is regarded as the weak absorption channel in the super channel-pair, which is applied to the following analyses, and the mean of some strong CO2 absorption channels is regarded as the strong absorption channel in the super channel-pair, as shown in Fig. 2. The selection of the strong absorption channel in the super channel-pair is presented in the following paragraphs.

    The effects of random noise are also analyzed in the strong absorption channel selection to avoid large errors. Figure 3 illustrates the errors due to instrument noise in differently positioned strong CO2 absorption channels if only one strong CO2 absorption channel is used in X CO2 retrieval. It is clearly shown that when the strong CO2 absorption channel located at the far wing is used, large errors could be introduced due to the reduced information content of CO2 (Fig. 3a); while in the line center, low signal-to-noise ratio (SNR) results in large uncertainty (Fig. 3c).

    The line strength and absorption coefficients depend on pressure and temperature. For the ultra-high spectral resolution measurements, an inaccurate pressure and temperature profile will introduce extra errors in the retrieval of CO2. Figure 4 shows the X CO2 errors of each channel for a +1 hPa bias of surface pressure, which is calculated by comparing inversions with and without a 1 hPa change, when the channel is regarded as the only strong CO2 absorption channel. In the error calculations, the coefficients a and b are calculated under the surface pressure, while the "measurements" are given under the +1 hPa bias of the surface pressure. Similarly, Fig. 5 shows the X CO2 errors for a +1 K shift of the temperature profile. As shown in Fig. 4, inversion errors caused by the +1 hPa pressure bias of most channels are positive, except for some channels in the weak absorption area. To reduce the impact of pressure, only channels with an inversion error of less than 1 ppm are selected to be the component of the strong absorption channels in the super channel-pair. Different from that of pressure, the error due to the +1 K offset of the temperature profile could be either positive or negative, and it could be minimized by careful channel selection in real retrievals.

    Figure 9.  The variation of X CO2 anomalies between the DOAS-like method and TCCON results with SZA at Tsukuba and Bremen.

    Inaccurate wavelength registration is another source of error in retrievals. As shown in Fig. 6, the errors in X CO2 retrieval due to a frequency shift of 0.003 cm-1 could be up to 25 ppm if only one strong CO2 absorption channel located on one side of the line center is used. But if strong CO2 absorption channels located on both sides of the line center are used, the errors due to the frequency shift tend to be very small, or even zero.

    Based on the above sensitivity studies, and the additional removal of the channels with strong H2O absorption and Fraunhofer lines, the final 588 strong CO2 absorption channels are used in our retrieval. In order to evaluate the dependence of the DOAS-like method on the atmospheric state uncertainty, one-year prior profiles in Tsukuba, as shown in Fig. 7, are used in simulated inversions. The results for both the DOAS-like and spectral fitting methods are listed in Table 1, in which errors of 1 K for the temperature profile, +5% for water vapor, +1 hPa for surface pressure, and 0.001 cm-1 for frequency offset are assumed for solar zenith angles (SZAs) at 20° and 70°. For specific atmospheric parameter analysis, both the spectral fitting and DOAS-like method have one unknown state vector. Relatively, the DOAS-like retrievals are less sensitive to the temperature and H2O uncertainties, especially for large SZAs and high H2O amounts. The effects of surface pressure and frequency shift to being slightly better in the spectral fitting method.

    Figure 10.  The inversion X CO2 of TCCON (black) and the DOAS-like method (blue) in Tsukuba and Bremen after correction.

3. Case studies and comparisons
  • To validate the DOAS-like algorithm, TCCON data in Tsukuba, Japan (36.0513°N, 140.1215°E) and Bremen, Germany (53.10°N, 8.85°E) are used. The spectra at both stations are measured with an FTS (IFS 125HR, Bruker Optics GmbH, Germany). The absorption spectrum is calculated by a Fourier transform of the interferogram, which is formed by beams reflected from a moving mirror and a static mirror. The resolution and sample rate of the FTS are determined by the maximum optical path differences (MOPDs) and speed of the moving mirror. The MOPDs of the FTS in Tsukuba and Germany are 45.01 and 64.29 cm, respectively. The retrieved X CO2 using the DOAS-like method are illustrated in Fig. 8 (left panels), and the results of the official TCCON algorithm are also included for comparison. At first sight, the X CO2 of the DOAS-like method is smaller than that of the official TCCON algorithm. After comparing the difference between the TCCON and DOAS-like retrievals with the SZA (Fig. 8, right panels), we find that the difference is linearly dependent on the SZA (Fig. 9). Moreover, the linear relationship does not vary with time and place. For the TCCON results, a post-retrieval algorithm is used to correct an airmass-dependent bias based on the assumption that any symmetric variability within a day should be an artifact (Deutscher et al., 2010; Wunch et al., 2011a). Through a simple correction process in which the linear dependency on the SZA is removed, the DOAS-like and TCCON results agree well with each other, as shown in Fig. 10. The standard deviation of the difference between the TCCON and DOAS-like methods is less than 0.8 ppm, both in Tsukuba and Bremen (Fig. 11). This suggests that the DOAS-like algorithm provides retrievals with similar precision to TCCON. However, the temporal variability of the atmospheric state in Fig. 12 limits the possibility of a higher inversion accuracy. Certainly, there could be many other factors for the low values of DOAS-like retrievals. For example, the solar lines provided by Kurucz used in our algorithm are not so good (Yoshida et al., 2013), and the FTS only focuses on the center of the solar disk due to its very small FOV. This inaccurate extra-terrestrial solar spectrum may be a factor for our lower value of X CO2.

4. Conclusions and future directions
  • A new algorithm using a channel-pair ratio to derive X CO2 is presented in this paper. The algorithm is similar to that of the DOAS method. For the purpose of channel selection, the effects of solar lines, water vapor, air temperature, pressure, instrument noise and wavelength registration shift on the retrieval error are analyzed through a series of sensitivity tests. One super channel-pair is used in the retrieval algorithm. FTS measurements at the TCCON stations in Tsukuba and Bremen are used to validate the new algorithm by comparing our results with the official TCCON product. Our X CO2 results are lower than those of TCCON with airmass correction. Taking the TCCON data as a reference, our results are further corrected using an SZA-dependent method. After the correction, our corrected results agree well with those of the TCCON products, suggesting that this new algorithm is useful. However, due to insufficient ground measurements, the new retrieval method is validated by observations at only two stations. Clearly, a thorough validation with extensive

    Figure 11.  Comparison between the results of TCCON and the DOAS-like method after correction at Tsukuba and Bremen. The standard deviations are shown in each panel.

    observation is warranted for our DOAS-like algorithm.

    Figure 12.  The time series of surface pressure (left panels), surface temperature (middle panels), and H2O column (right panels) for each day.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return