Abstract
High-resolution solar absorption spectra, observed by ground-based Fourier Transform Infrared spectroscopy (FTIR), are used to retrieve vertical profiles and partial or total column concentrations of many trace gases. In this study, we present the tropospheric CO2 columns retrieved by mid-infrared solar spectra over Hefei, China. To reduce the influence of stratospheric CO2 cross-dependencies on tropospheric CO2, an a posteriori optimization method based on a simple matrix multiplication is used to correct the tropospheric CO2 profiles and columns. The corrected tropospheric CO2 time series show an obvious annual increase and seasonal variation. The tropospheric CO2 annual increase rate is 2.71 ± 0.36 ppm yr-1, with the annual peak value in January, and CO2 decreases to a minimum in August. Further, the corrected tropospheric CO2 from GEOS-Chem simulations are in good agreement with the coincident FTIR data, with a correlation coefficient between GEOS-chem model and FTS of 0.89. The annual increase rate of XCO2 observed from near-infrared solar absorption spectra is in good agreement with the tropospheric CO2 but the annual seasonal amplitude of XCO2 is only about 1/3 of dry-air averaged mole fractions (DMF) of tropospheric CO2. This is mostly attributed to the seasonal variation of CO2 being mainly dominated by sources near the surface.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Carbon dioxide (CO2), as the most important anthropogenic enhanced global greenhouse gas, and has continued to increase since the industrial revolution [1]. This CO2 increase has been deeply influencing global climate, flora and fauna. From 2002 to 2011, the rate of increase in CO2 is approximately 2.0 ± 0.1 ppm yr−1, coinciding with emission reduction efforts worldwide [2]. In recent times, anthropogenic sources have emitted more than 9 GtC yr−1 CO2 into the atmosphere [2,3]. The anthropogenic CO2 emission increase is mainly attributed to anthropogenic activities, like fossil fuel burning, industrial activity and land use-change [4]. Prediction of the evolution of atmospheric greenhouses requires detailed understanding of their sources and sinks [5]. Due to the complex seasonal variability of trace gas sources and transport pathways, the seasonal cycle of trace gases in the troposphere is different from the gases at the surface [6]. Accurate and continued observations of tropospheric CO2 have given valuable insights on climate change, carbon cycle, model simulation and CO2 emission sources [7]. The CO2 concentrations throughout the troposphere can complement very well in situ results.
Recently several methods and techniques have been utilized to measure concentration of atmospheric CO2. The in-situ technique offers high accuracy and precision, usually utilized to monitor tropospheric CO2 [8]. Other examples include, nondispersive infrared gas analyzer (NDIR) and cavity ring-down spectroscopy (CRDS) used to monitor long-time variation of ground-surface CO2, integrated-path differential absorption (IPDA) light detection and ranging (LIDAR) system to measure the column-averaged dry-air mixing ratio of CO2 (XCO2), and a differential absorption lidar (DIAL) system that could provide the CO2 profile variation below the boundary layer [9–11]. However, these measurements are strongly influenced by local sources and small-scale processes [12]. The space-based measurements, such as the Greenhouse gases Observing SATellite (GOSAT), the Greenhouse gases Observing SATellite-2 (GOSAT-2), the Orbiting Carbon Observatory 2(OCO-2) and TanSat are tasked to observe the columns of atmospheric CO2, which provide global atmosphere CO2 concentrations [13–15]. But the satellite data must be validated by ground-based instruments, as satellite observations can be affected by surface properties and aerosols, measurement features that do not impact ground-based instruments. [16–18].
The ground-based high resolution Fourier transform spectrometer (FTIR) is increasingly used to observe total columns and vertical profiles of atmospheric trace gases [20–22]. The global Network for the Detection of Atmospheric Composition Change (NDACC) utilize middle-infrared Fourier transform infrared spectrometers to monitor the profiles of tropospheric and stratospheric trace gases [20,23–26]. Here we use a ground-based high resolution FTIR spectrometer to detect tropospheric CO2.
In this paper, the objective is to use the mid-infrared solar spectra to retrieve vertical profiles and tropospheric CO2 columns. The outline of the paper is as follows. In Sect. 2, we describe the measurement site, observation instrument and simulation model. Sect. 3 presents the retrieval setup for measurement spectra, retrieval results correction and error estimates. In addition, the time series of tropospheric CO2 columns from FTIR observations and GEOS-chem simulations are presented and compared in Sect. 4. Sect. 5 summarizes our findings and conclusions.
2. Methods
2.1 Measurement site
The observation site (31.540 °N, 117.100° E, 35.5 m a.s.l.), adjacent to a reservoir named DongPu, is located in the north-west of Hefei city in eastern China. We installed the ground-based remote sensing system in 2014 (Fig. 1), consisting of a high-resolution FTIR instrument (Bruker IFS 125HR), which has nine scanner compartments (with a maximum spectral resolution of 0.00096 cm−1), and a solar tracker (A547N) that is installed on the roof of the lab building and reflects the solar beam into the FTIR instrument (with the tracking precision of 0.1°), to monitor the atmospheric trace gas profiles and total columns. A meteorological station (ZENO, Coastal Environmental Systems, USA), on the roof of the lab building that is near the solar tracker, has recorded ground pressure, ground temperature, relative humidity, ground wind speed, ground wind direction, ground solar radiation, rain or snow and leaf wetness data since Sep. 2015 [27,28].
2.2 FTIR measurement
The FTIR instrument collects near-infrared and mid-infrared solar absorption spectra alternately on clear days. The FTIR instrument is equipped with a calcium fluoride (CaF2) beamsplitter and an indium gallium arsenide (InGaAs) detector, which can record near-infrared solar spectra (3800–11,000 cm−1) with a spectral resolution of 0.02 cm−1, according to the standard TCCON set up [27]. The mid-infrared spectra (600–4500 cm−1), with a spectral resolution up to 0.005cm-1, are collected by a potassium bromide (KBr) beamsplitter along with indium antimonide (InSb) and mercury cadmium telluride (MCT) detectors. Knowledge of the instrument line shape (ILS) is required to diagnose the alignment of the spectrometer and hence to retrieve total columns of gases from measurements accurately [27,29]. The ILS of the FTIR instrument is monitored by low-pressure HCl and HBr cell measurements [26,27]. The FTIR trace gases uncertainty estimates are computed by the standard deviation of the daily observed results. Aircraft measurements are used to calibrate the FTIR results within TCCON, and has yet to happen at our site, but we will do this work if we have the opportunity later [30,31].
2.3 GEOS-Chem model
The GEOS-Chem is a global 3-D chemical transport model, which has been used in a wide range of research topics on atmospheric composition, including CO2, is driven by the Goddard Earth Observing System (GEOS) assimilated meteorological input data from the NASA Global Modeling and Assimilation Office [32,33]. In this study, we use a spatial resolution of 2° latitude × 2.5° longitude with 47 vertical layers [34]. For atmospheric CO2 concentration, annual fossil fuel emissions in the model are estimated with the Multi-resolution Emission Inventory for China (MEIC) [35,36]. Biomass burning emissions are simulated from the Global Fire Emission Database (GFED2016) [37]. Ocean exchange and annual biofuel flux are acquired according to the paper of Takahashi (2009) [38], and Yevich and Logan (2003) [39], respectively. Balanced ecosystem exchange are calculated by the Carnegie-Ames-Stanford Approach (CASA) biosphere model [40]. It should be noted that the balanced ecosystem exchange contribute no net annual uptake of CO2, but they make the greatest contribution to the seasonal cycle of atmospheric CO2 over most of the globe with the largest impact in the Northern Hemisphere[34]. Ship and plane emissions of CO2 emission are estimated from the International Comprehensive Ocean Atmosphere Data Set (ICOADS) and System for assessing Aviations Global Emissions (SAGE) [40,41]. The concentration of atmospheric CO2 in this study is calculated using these settings. The net ecosystem uptake, which is used to adjust net ecosystem exchange (NEE), are computed from TransCom climatology [42].
3. Tropospheric CO2 retrieval
3.1 CO2 retrieval setup
The tropospheric CO2 concentrations are retrieved using SFIT4 (version 0.9.4.4), which is developed from SFIT2 [43]. The SFIT algorithm is based on the optimal estimation method to obtain the vertical profiles and columns of trace gases [44]. The CO2 micro-windows are selected from mid-infrared spectra which contain strong and well-isolated CO2 absorption lines [21,25,45]. The retrieval parameters and interfering gases for CO2 at the Hefei site is listed in Table 1. Four absorption windows of CO2 were fitted simultaneously. In addition to CO2, the spectroscopic signatures of interfering gases are also considered to reduce fitting residuals. The interfering gases, H2O and HDO, are considered as two different species for fitting. The a priori profile of CO2 and other interfering gases except H2O and HDO are derived from the Whole Atmosphere Community Climate Model (WACCM). The a priori profiles of atmospheric temperature, pressure and H2O are adopted from National Centers for Environmental Protection (NCEP) reanalysis data [46]. The spectroscopic line parameters for CO2 and interfering species are derived from the HITRAN 2012 linelist and recent updated versions [47]. The columns of H2O, HDO, CH4 and O3 are also simulated during the CO2 retrieval. A typical spectrum collected at 04:44 UTC 03 Feb. 2016, with a solar zenith angle of 48.82° is shown in Fig. 2. The fitting residuals of each window fluctuate within ± 0.2% and the averaged root-mean-square error is 0.069%, showing that the spectrum is well fitted.
The averaging kernels are used to depict the sensitivity of height dependence of the retrieved profile to concentration perturbations at the various atmospheric levels [44]. The typical averaging kernel of CO2 is shown in Fig. 3(a). The most sensitive height range of CO2 is mainly located at 5-9 km and 12-25 km. The trace of averaging kernel matrix is characterized as the degrees of freedom for signal (DOFs). The typical DOFs for CO2 over Hefei is about 2.86. This indicates that the retrieved profile can be separated into two independent partial columns, because the criterion for separating partial columns into separate layers is that the partial column DOFs of gases is greater than 1.0 [26,48]. To ensure accurate information from the CO2 retrievals in the troposphere and minimize any stratosphere influence, i.e., the stratosphere–troposphere exchange, the upper altitude limit for the troposphere is set to 12 km, which is lower than the mean value of the tropopause (15 km) [26]. The red lines represent the averaging kernels above 12 km and the blue lines represent the averaging kernels in the troposphere (for altitudes < 12 km) in Fig. 2. The typical DOFs in the troposphere for CO2 is 1.27.
3.2 Posteriori optimization for the tropospheric CO2 retrieval
Rodgers and Connor (2003) provide some mathematical background to a posteriori modification of a retrieval [44]. The profile retrievals are not necessarily perfect estimates of the desired quantity, including noise, but will have an averaging kernel that is not perfect [44]. As discussed in Sect.4 of Rodgers and Connor (2003), what is the best estimate of a function of the state vector. Since the tropospheric CO2 has been considered as a function of the retrieved state vector: in the case of a retrieval that is based on the Bayesian theorem, one can simply apply the function to the retrieved state and get the expectation of the state, given the measurement (see Sect. 4.1 of Rodgers and Connor, 2003). However, in our case the a priori profile and its covariance are not completely known. Instead we use the WACCM simulated profile as the a priori and a constant covariance. We then calculate the tropospheric partial column of CO2 from this approach. However, the tropospheric partial column of CO2 can be better estimated with a posteriori modification. Equation (19) of Rodgers and Connor (2003) shows how this better estimator can be calculated from the retrieved state [44]. Sepulveda et al. (2014) also use this posteriori optimization method to eliminate the upper tropospheric/lower stratospheric (UTLS) influence on the tropospheric partial column of CH4 [20].
As shown in Fig. 3(a), the peak of the tropospheric averaging kernels (blue lines) mainly appear in the troposphere and the peak of stratospheric averaging kernels (red lines) occur mainly in the stratosphere. However, the negative value for the tropospheric averaging kernels in stratosphere means that the stratosphere also contributes to the retrieved tropospheric CO2 [20]. This means that the variation of stratospheric CO2 will significantly affect the retrieved tropospheric CO2. Therefore, we rewrite the retrieved averaging kernel matrix A as:
To reduce the cross-dependencies and separate the tropospheric column, that has a weak dependency on stratosphere, the a posteriori optimization method based on a simple matrix multiplication is used to correct the retrieval results. The matrix C can be written as:
I is the unit matrix. The dimension of the correction matrix C is identical to the averaging kernel matrix. The a posteriori optimized averaging kernel is obtained from C multiplied by the averaging kernel matrix:
The corrected averaging kernels matrix ${A_{corr}}$ is depicted in the Fig. 3(b). The overlap between the tropospheric averaging kernel and the stratosphere averaging kernels decreased between approximately 5 and 20 km, which means the corrected averaging kernel is much less affected by the stratosphere than the original averaging kernel [20]. The column averaging kernels in Fig. 4 are calculated based on the derived partial columns and averaging kernel [49]. The column averaging kernels have the property that they are 1 in the ideal case, but for a real retrieval can be more or less than one depending on the dominance of the apriori profile and other factors. The tropospheric partial column averaging kernel depicts this change in stratospheric cross-sensitivity clearly (Fig. 4). The kernel values improved significantly from negative to positive values in the range from 18-42 km. The shape of the sensitivity of the corrected and uncorrected partial column averaging kernels are similar in tropospheric region. It is apparent that this correction improved the separation between the retrieved tropospheric and stratospheric amounts.
3.3 Error estimations
The retrieved error covariance is calculated following the method of Rodgers and Conor (2003) [44]. The error components for error estimations are listed in Table 2. The square root of diagonal elements of the error covariance matrices, for a typical measurement spectrum (observed in 03:36 UTC, 23 Jun. 2019, with SZA 11.72°), are depicted in Fig. 5(a). The vertical structures of the error profiles reflect the effect of the propagation of different errors in the retrieval. As shown in Fig. 5, the dominant error sources are spectral and measurement error.
The error estimation for the a posteriori optimization is calculated by:
The small cross-dependency of the retrieved lower tropospheric mole fraction on the state of the stratosphere can significantly affect the quality of the lower tropospheric column product [20]. This means that cross-dependencies are one of the leading error sources. This error is referred to as the “sensitivity error”. This error covariance matrix is calculated as:
Table 3 summarizes DOFs values, sensitivity error and the total errors for the tropospheric partial column of CO2 with and without a posteriori optimization. The total error include statistical and systematic errors. For the corrected results we get a sensitivity error that is still a significant fraction of the variability in the raw data. The a posterior optimization is more important for the sensitivity error so that the cross-talk between the partial columns is reduced without affecting the statistical and systematic errors. The DOFs value slightly decreases, but still meets the criterion of separating troposphere and stratosphere (DOFs greater than 1.0). The results show that the total error calculated from the error components list in Table 2 slightly decreases after correction in troposphere. The total error in this study is comparable to the error at the Karlsruhe site, which reported a total error of 4.24% [21]. The a posteriori optimization reduces the stratospheric influence on the tropospheric partial column from the CO2 retrieval. This is reflected in the reduction of the sensitivity error (from 1.33% to 0.97%), but the DOFs also decrease, which is mainly losses in information at and above 12 km that were being folded into the troposphere. Since we separate the troposphere and stratosphere at 12 km, accordingly the tropospheric DOFs reduce (from 1.14 to 1.03). This is the classic trade-off between vertical resolution and profile accuracy: improve one and the other worsens. However in our case, this apparent reduction in DOFs may well be regarded as an improvement as the information loss of 0.11 was mainly a contribution to the tropospheric partial column from the stratosphere.
4. Results and discussion
4.1 Time series of tropospheric CO2
The retrieved state is recalculated by considering the a posteriori optimization matrix C:
The DMF of tropospheric CO2 time series show an obvious annual increase and seasonal variation in Fig. 6. The annual increase of CO2 is due to the contribution of fossil fuel combustion. The CO2 annual increase rate is 2.71 ± 0.36 ppm yr-1. This rate is slightly higher than the values at other observation stations in the Northern Hemisphere. The annual increase rate of XCO2 is 2.3 ppm yr-1 from 2011 to 2015 at the Saga site, which is the TCCON site in Japan near to our site [51]. The annual increase rate of CO2 is about 2 ppm yr-1 from 2005 to 2009 at a rural site near Beijing and 2.05 ppm yr−1 in Mauna Loa with in-situ observation [52], also similar to the value of 1.96 ppm yr-1 from 2003 to 2011 in the Northern Hemisphere observed from SCIAMACHY [53]. This high increase rate is caused by the recent rapid rise of fossil fuel combustion in the Hefei area.
The CO2 seasonality is also captured by the Fourier series fitting (Fig. 6). The annual peak value of CO2 always appears in January, and then to a minimum in August. This seasonal variation amplitude of CO2 is mainly caused by respiration and photosynthesis from plants. To clearly show the seasonal variation, we de-trend the tropospheric CO2 time series. This de-trending is performed by subtracting the long-term trends, and the de-trended daily averaged CO2 of each year are shown in Fig. 7. Table 4 summarizes the seasonal amplitude in each observation year. The seasonal amplitudes are different, i.e., the seasonal amplitude difference is 7.72 ppm between 2016 and 2017. This amplitude discrepancy is mainly due to the different temperature variation in different years because the CO2 seasonal amplitude is strongly influenced by temperature [53–55]. The averaged seasonal amplitude of tropospheric CO2 over observed period is 18.53 ± 3.07 ppm. The total column seasonal amplitude of XCO2 are 8.31 and 13.56 ppm, respectively, over Hefei from 2014 to 2015 and from 2015 to 2016 [27]. The seasonal cycle was mainly driven by biosphere-atmosphere exchange, and it mostly happens near the surface, therefore this seasonal variability is more clearly captured in the troposphere compared with the total column.
4.2 Comparison with GEOS-Chem model
The FTIR DMF of tropospheric CO2 time series was used to compare with GEOS-Chem model simulated results to evaluate the performance of GEOS-Chem model simulation in eastern China. We used the GEOS-Chem model to simulate the CO2 mixing ratio profiles from July 2015 to December 2019 with 1 h time frequency. We sampled the simulation results at the grid cell covering the FTIR site, and interpolated them vertically to the vertical grid of FTIR retrieval [33]. To eliminate the influence of different a priori profiles and averaging kernel matrices, the profile from model simulation is smoothed by the FTIR retrieval averaging kernel [44]:
The daily averaged DMF of tropospheric CO2 obtained by the FTIR instrument and the GEOS-Chem model are presented in Fig. 8. As shown in Fig. 8, the DMF of tropospheric CO2 from GEOS-Chem simulation are in good agreement with the coincident FTS data. Figure 9 shows the GEOS-Chem model results comparisons with coincident FTIR data. The black line is the linear regression curve of the FTIR and GEOS-Chem results and the linear regression method according to York (2004) is used [56]. The GEOS-Chem model data are smoother than the FTIR measurements. Further, the FTIR and GEOS-Chem model result shows high correlation, with the correlation coefficient (r) of 0.89, but also show a bigger bias between FTIR and GEOS-Chem model data, especially in August and September. A positive CO2 annual trend of 2.41 ppm yr-1 is calculated from the GEOS-Chem data, comparable to the FTIR data (2.71 ± 0.36 ppm yr-1). The model also has a good simulation record compared with other FTS and satellite CO2 observation results [57,58]. The seasonal variations of GEOS-Chem model data is similar to the FTIR data, with the maximum and minimum of DMF of tropospheric CO2 occurring in February and August, respectively. The mean bias between FTIR and GEOS-chem data, expressed as FTIR minus GEOS-chem, is -1.66 ± 7.31 ppm. The averaged seasonal amplitude of GEOS-chem is 10.68 ± 0.52 ppm and this is only one half of the FTIR seasonal amplitude (18.53 ± 3.07 ppm), which means GEOS-chem model underestimates the seasonal amplitude of the tropospheric CO2 column. Because the model overestimates the minimum value and underestimates the peak value in each year (Fig. 8), especially the minimum value in August 2016, the average difference between FTIR and GEOS-chem is 13.97 ± 3.63 ppm. The mean biases of the minimum and maximum peaks between simulated and observed values are 8.07 ± 4.78 and 3.51 ± 1.44 ppm, respectively.
The bias in the DMF of tropospheric CO2 between the FTIR observations and GEOS-Chem simulations can be largely attributed to the following reasons: the grid cell and timing of the GEOS-Chem simulated data are not exactly the same as those of FTIR observations; and the accuracy of GEOS-Chem model simulation is affected by several factors, including the emission inventory (with an uncertainty of 31%), the atmosphere-biosphere model, and the meteorological conditions [33, 59–61]. Uncertainties in the model inputs and parameters result in the differences between the GEOS-Chem and the FTIR data.
4.3 Comparison with mole fraction of the total column and surface concentration
The total column abundances of CO2 are retrieved from the near-infrared solar absorption spectra [19]. The retrieval strategies of total column CO2 are harmonized with those of other TCCON sites. We use the GGG2014 software package to analyze the spectra [62]. The core of software is a nonlinear least squares spectral fitting algorithm (GFIT), which scales the a priori profile for the best fit based on the high-resolution transmission molecular absorption (HITRAN) database [63]. The column abundances of gases are then converted into column-averaged dry-air mole fractions (XCO2) by dividing them by the total column of dry air. The division helps to cancel some potential systematic errors, including pointing errors, ILS uncertainty, zero level offsets, and source brightness fluctuation [19]. The method of determining the column of dry air is done by calculating the column of O2 from the spectrum, divided by the dry-air mole fractions of O2 (0.2095).
A wavelength-scanned cavity ring-down spectrometer (CRDS, Picarro model-2014i) has been used to measure surface CO2 concentration since November 2015. To ensure the measurement stability, the CRDS instrument is placed inside the FTIR lab. Standard gases with three different concentrations are used to calibrate the measurement system once a month.
Time series of the daily averaged DMF of tropospheric CO2 and XCO2 are shown in Fig. 10. As show in Fig. 10, these data sets show obvious seasonal variation and an annual increase trend. The annual increase rate of CO2 for the two datasets are in good agreement. The annual increase rate of DMF of tropospheric CO2 and XCO2 are 2.71 ± 0.36 ppm yr-1 and 2.64 ± 0.13 ppm yr-1, respectively. However, as show in Fig. 10, the annual maximum values of DMF of tropospheric CO2 is higher than XCO2 and the annual minima lower than XCO2. The difference between daily averaged DMF of tropospheric CO2 and XCO2 is 17.79 ± 6.62 ppm, while the averaged seasonal amplitude of XCO2 is 7.1 ± 0.98 ppm which is about 1/3 of DMF of tropospheric CO2 seasonal amplitude (18.53 ± 3.07 ppm). The correlation coefficient (R) between the daily averaged DMF of tropospheric CO2 and XCO2 is 0.56 (p < 0.01) (Fig. 11).The seasonal variation of CO2 is mainly dominated by the respiration and photosynthesis of plants, mostly happening at the surface. It would appear that the DMF of tropospheric CO2 has a slightly better ability to capture this seasonality, due to the XCO2 being slightly more diluted and influenced by stratospheric CO2.
The time series of the daily means from the surface concentration, DMF of tropospheric CO2 and XCO2 are shown in Fig. 12. As shown in Fig. 12(a) and (c), the DMF of tropospheric has a smaller bias (8.12 ± 9.28 ppm)with the surface CO2 than the total column derived XCO2 (15.05 ± 9.33 ppm). Further, as show in Fig. 12(b) and (d), the correlation coefficient (R) between the DMF of tropospheric CO2 and surface CO2 is 0.61 (p < 0.01, with the slope of 0.62 ± 0.05), which is also higher than the correlation between XCO2 and surface CO2 (R = 0.54, with the slope of 1.71 ± 0.17). These results still show though that both the DMF of tropospheric CO2 and XCO2 are also influenced by sources in the mid and upper troposphere that are not sampled by the in-situ measurement.
By way of a comparison, a similar test of the bias between the surface CO2 with uncorrected DMF of tropospheric CO2 and the total column XCO2 computed from the MIR (using pressure rather the air column), gives similar biases (6.11 ± 9.70 ppm and 13.51 ± 8.50 ppm, respectively). The correlation coefficient (R) between the uncorrected DMF of tropospheric and surface CO2 is 0.53 while the surface CO2 and XCO2MIR correlations coefficient is R = 0.58. This shows that the a posteriori correction has removed some of the stratospheric influence on the total column products. And the surface concentrations of CO2 measured by the in-situ instrument are higher than XCO2 derived from FTIR data at Anmyeondo station, Korea, which is the nearest TCCON site to our site [64]. Also, the averaged difference between the surface concentrations of CO2 and XCO2 is 13 ppm at Jungfraujoch site [65].
Further work is planned to produce an SFIT4 derived profile from the NIR spectra where this cross-dependency from the stratosphere can be reduced, combined with the know improvements that remove instrumental and other effects.
5. Conclusion
In this study, we present tropospheric CO2 columns retrieved by ground-based high resolution FTIR remote sensing measurements. The tropospheric CO2 concentration were retrieved from mid-infrared solar absorption spectra based on the optimal estimation algorithm, SFIT4.
The negative values for the tropospheric averaging kernels that occur in the stratosphere indicates that stratospheric CO2 can affect the tropospheric CO2 variation. To reduce this dependency on stratosphere variations, a posteriori optimization method based on a simple matrix multiplication was used to correct the retrieved profile. The tropospheric averaging kernel is less negative in the stratosphere after the a posteriori optimization, which means the corrected averaging kernel, and hence the CO2 retrieved profile are much less affected by the stratosphere than the retrieval with the raw averaging kernel. The retrieved error covariances are calculated, and the dominant error sources are spectral and measurement error. The sensitivity of the troposphere partial column is clearly reduced after this correction. The total error slightly decreases from 1.22% to 1.18% after applying the correction in troposphere.
The corrected tropospheric CO2 column during the period from August 2015 to December 2019 are presented. The tropospheric CO2 time series shows an obvious annual increase and seasonal variation. The CO2 annual increase rate is 2.71 ± 0.36 ppm yr-1. The annual peak of CO2 always occurs in January, and then decreases to a minimum in August. The mean annual amplitude of CO2 over the observed period is 18.53 ± 3.07 ppm, and due to temperature variation in different years, the seasonal amplitude for each year has obvious variability.
Further, the observed DMF of tropospheric CO2 is compared with the smoothed simulation data from GEOS-chem model and XCO2 observed from near-infrared solar absorption spectra. The tropospheric CO2 column from GEOS-Chem simulation are in good agreement with the coincident FTIR data, and the correlation coefficient (r) between the two data is 0.89. But the GEOS-chem model underestimates the seasonal amplitude of tropospheric CO2 column, with the averaged seasonality amplitude of 10.68 ± 0.52 ppm that is only one half of the FTIR seasonal amplitudes. This discrepancy is largely attributed to the fact that the grid cell and timing of model data mismatch to FTIR data and the model inputs have high uncertainties. The annual increase rate of XCO2 is in good agreement with the DMF of tropospheric CO2 but the annual seasonal amplitude of XCO2 is only about 1/3 of DMF of tropospheric CO2. This is mostly attributed to the seasonal variation of CO2 being mainly dominated by sources near the surface, therefore this seasonal variability is more clearly captured in the troposphere compared with the total column derived XCO2
Appendix
The O2 total column can be used as an indicator to confirm that the FTIR instrument was stable and prove that the year to year variability in the CO2 drawdown is not related to instrumental problems. The Appendix Fig. 13 presented the time series of O2 total column. As shown in Fig. 13, the O2 total column has an obvious seasonal variation that is stable year to year.
The averaging kernel, partial column averaging kernel and the error covariance matrices observed in 03:36 UTC, 23 Jun. 2019(solar zenith angle: 11.72°) are depicted in Fig. 14 and Fig. 15.
Funding
National Key Technology R&D Program of China (2019YFC0214702, 2018YFC0213201, 2016YFC0200404, 2017YFC0210002, 2018YFC0213104, 2019YFC0214802); National Natural Science Foundation of China (41775025, 41722501, 91544212, 51778596, 41575021, 41977184); Major Projects of High Resolution Earth Observation Systems of National Science and Technology (05-Y30B01-9001-19/20-3); Strategic Priority Research Program of the Chinese Academy of Sciences (XDA23020301); National Key Project for Causes and Control of Heavy Air Pollution (DQGG0102, DQGG0205); Natural Science Foundation of Guangdong Province (2016A030310115).
Acknowledgements
The processing environment of SFIT4 and some plot programs are provided by National Center for Atmospheric Research (NCAR), Boulder, Colorado, USA. The NDACC networks is acknowledged for supplying the SFIT software and advice.
Disclosures
The authors declare no conflicts of interest.
References
1. J. R. Petit, J. Jouzel, and D. Raynaud, “Climate and atmospheric history of the past 420,000 years from the Vostok Ice Core,” Nature 399(6735), 429–436 (1999). [CrossRef]
2. IPCC: Climate Change 2014: Synthesis Report, Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Core Writing Team, Pachauri, R. K., and Meyer, L. A., IPCC, Geneva, Switzerland, 151 pp., (2014).
3. C. Le Quéré, G. P. Peters, and R. J. Andres, “Global carbon budget 2013,” Earth Syst. Sci. Data 6(1), 235–263 (2014). [CrossRef]
4. M. Kiel, D. Wunch, and P. O. Wennberg, “Improved retrieval of gas abundances from near-infrared solar FTIR spectra measured at the Karlsruhe TCCON station,” Atmos. Meas. Tech. 9(2), 669–682 (2016). [CrossRef] .
5. E. Sepúlveda, M. Schneider, and F. Hase, “Longterm validation of tropospheric column-averaged CH4 mole fractions obtained by mid-infrared ground-based FTIR spectrometry,” Atmos. Meas. Tech. 5(6), 1425–1441 (2012). [CrossRef] .
6. B. Pak, R. Langenfelds, and S. Young, “Measurements of biomass burning influences in the troposphere over southeast Australia during the SAFARI 2000 dry season campaign,” J. Geophys. Res. 108(D13), 8480 (2003). [CrossRef] .
7. P. Ciais, C. Sabine, and G. Bala, “Carbon and other biogeochemical cycles, Chapter, in: Climate Change 2013: The Physical Science Basis,” Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK and New York, NY, USA, 465–570, (2013).
8. Y. Sawa, T. Machida, and H. Matsueda, “Seasonal changes of CO2, CH4, N2O, and SF6 in the upper troposphere/lower stratosphere over the Eurasian continent observed by commercial airliner,” Geophys. Res. Lett. 42(6), 2001–2008 (2015). [CrossRef] .
9. M. F. Schibig, M. Steinbacher, and B. Buchmann, “Comparison of continuous in situ CO2 observations at Jungfraujoch using two different measurement techniques,” Atmos. Meas. Tech. 8(1), 57–68 (2015). [CrossRef] .
10. U. N. Singh, T. F. Refaat, and S. Ismail, “Feasibility study of a space-based high pulse energy 2 mu m CO2 IPDA lidar,” Appl. Opt. 56(23), 6531–6547 (2017). [CrossRef]
11. Y. Shibata, C. Nagasawa, and M. Abo, “Development of 1.6 µm DIAL using an OPG/OPA transmitter for measuring atmospheric CO2 concentration profiles,” Appl. Opt. 56(4), 1194–1201 (2017). [CrossRef]
12. S. C. Olsen and J. T. Randerson, “Differences between surface and column atmospheric CO2 and implications for carbon cycle research,” J. Geophys. Res. 109(D2), D02301 (2004). [CrossRef] .
13. C. Frankenberg, R. Pollock, and R. A. M. Lee, “The Orbiting Carbon Observatory (OCO-2): spectrometer performance evaluation using prelaunch direct sun measurements,” Atmos. Meas. Tech. 8(1), 301–313 (2015). [CrossRef] .
14. S. S. Kulawik, C. O’Dell, and V. H. Payne, “Lower-tropospheric CO2 from near-infrared ACOS-GOSAT observations,” Atmos. Chem. Phys. 17(8), 5407–5438 (2017). [CrossRef] .
15. Y. Ran and X. Li, “TanSat: “a new star in global carbon monitoring from China,”,” Sci. Bull. 64(5), 284–285 (2019). [CrossRef] .
16. M. Reuter, H. Bovensmann, and M. Buchwitz, “Retrieval of atmospheric CO2 with enhanced accuracy and precision from SCIAMACHY: Validation with FTS measurements and com-parison with model results,” J. Geophys. Res. 116(D4), D04301 (2011). [CrossRef] .
17. Y. Yoshida, N. Kikuchi, and I. Morino, “Improvement of the retrieval algorithm for GOSAT SWIR XCO2 and XCH4 and their validation using TCCON data,” Atmos. Meas. Tech. 6(6), 1533–1547 (2013). [CrossRef] .
18. D. Wunch, P. O. Wennberg, and G. Osterman, “Comparisons of the Orbiting Carbon Observatory-2 (OCO-2) XCO2 measurements with TCCON,” Atmos. Meas. Tech. 10, 2209–2238 (2017). [CrossRef] .
19. D. Wunch, G. C. Toon, and J. Blavier, “The Total Carbon Column Observing Network,” Phil. Trans. R. Soc. A 369(1943), 2087–2112 (2011). [CrossRef]
20. E. Sepulveda, M. Schneider, and F. Hase, “Tropospheric CH4 signals as observed by NDACC FTIR at globally distributed sites and comparison to GAW surface in situ measurements,” Atmos. Meas. Tech. 7(7), 2337–2360 (2014). [CrossRef] .
21. S. Barthlott, M. Schneider, and F. Hase, “Using XCO2 retrievals for assessing the long-term consistency of NDACC/FTIR data sets,” Atmos. Meas. Tech. 8(3), 1555–1573 (2015). [CrossRef] .
22. Y. Té, P. Jeseck, and B. Franco, “Seasonal variability of surface and column carbon monoxide over the megacity Paris, high-altitude Jungfraujoch and Southern Hemispheric Wollongong stations,” Atmos. Chem. Phys. 16(17), 10911–10925 (2016). [CrossRef] .
23. M. Zhou, B. Langerock, and C Vigouroux, “Atmospheric CO and CH4 time series and seasonal variations on Reunion Island from ground-based in situ and FTIR (NDACC and TCCON) measurements,” Atmos. Chem. Phys. 18(19), 13881–13901 (2018). [CrossRef]
24. M. Schneider, S. Barthlott, and F. Hase, “Groundbased remote sensing of tropospheric water vapour isotopologues within the project MUSICA,” Atmos. Meas. Tech. 5(12), 3007–3027 (2012). [CrossRef] .
25. S. Barthlott, M. Schneider, and F Hase., “Tropospheric water vapour isotopologue data (H162O, H182O and HD16O) as obtained from NDACC/FTIR solar absorption spectra,” Earth Syst. Sci. Data Discuss., doi:10.5194/essd-2016-9, (2016).
26. Y. Sun, C. Liu, and M. Palm, “Ozone seasonal evolution and photochemical production regime in the polluted troposphere in eastern China derived from high-resolution Fourier transform spectrometry (FTS) observations,” Atmos. Chem. Phys. 18(19), 14569–14583 (2018). [CrossRef]
27. W. Wang, Y. Tian, and C. Liu, “Investigating the performance of a greenhouse gas observatory in Hefei, China,” Atmos. Meas. Tech. 10(7), 2627–2643 (2017). [CrossRef] .
28. C. Shan, W. Wang, and C. Liu, “Regional CO emission estimated from ground-based remote sensing at Hefei site, China,” Atmos. Res. 222, 25–35 (2019). [CrossRef]
29. F. Hase, B. J. Drouin, and C. M. Roehl, “Calibration of sealed HCl cells used for TCCON instrumental line shape monitoring,” Atmos. Meas. Tech. 6(12), 3527–3537 (2013). [CrossRef] .
30. D. Wunch, G. C. Toon, and P. O. Wennberg, “Calibration of the Total Carbon Column Observing Network using aircraft profile data,” Atmos. Meas. Tech. 3(5), 1351–1362 (2010). [CrossRef] .
31. T. Tanaka, Y. Miyamoto, and I. Morino, “Aircraft measurements of carbon dioxide and methane for the calibration of ground-based high-resolution Fourier Transform Spectrometers and a comparison to GOSAT data measured over Tsukuba and Moshiri,” Atmos. Meas. Tech. 5(8), 2003–2012 (2012). [CrossRef] .
32. I. Bey, D. J. Jacob, and R. M. Yantosca, “Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation,” J. Geophys. Res. 106(D19), 23073–23095 (2001). [CrossRef] .
33. H. Yin, Y. Sun, and C Liu, “FTIR time series of stratospheric NO2 over Hefei, China, and comparisons with OMI and GEOS-Chem model data,” Opt. Express 27(16), A1225–A1240 (2019). [CrossRef]
34. R. Nassar, D. B. A. Jones, and P. Suntharalingam, “Modeling global atmospheric CO2 with improved emission inventories and CO2 production from the oxidation of other carbon species,” Geosci. Model Dev. 3(2), 689–716 (2010). [CrossRef] .
35. M. Li, Q. Zhang, and J.-I. Kurokawa, “MIX: a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP,” Atmos. Chem. Phys. 17(2), 935–963 (2017). [CrossRef] .
36. F. Liu, Q. Zhang, and D. Tong, “High-resolution inventory of technologies, activities, and emissions of coal-fired power plants in China from 1990 to 2010,” Atmos. Chem. Phys. 15(23), 13299–13317 (2015). [CrossRef] .
37. L. Giglio, J. Randerson, and G. van der Werf, “Analysis of daily, monthly, and annual burned area using the fourth-generation global fire emissions database (GFED4),” J. Geophys. Res. Biogeosci. 118(1), 317–328 (2013). [CrossRef] .
38. T. Takahashi, S. C. Sutherland, and R. Wanninkhof, “Climatological mean and decadal change in surface ocean pCO2, and net sea-air CO2 flux over the global oceans,” Deep Sea Res. Part II 56(8-10), 554–577 (2009). [CrossRef] .
39. R. Yevich and J. A. Logan, “An assessment of biofuel use and burning of agricultural waste in the developing world,” Global Biogeochem. Cycles 17(4), 1095 (2003). [CrossRef] .
40. O. Endresen, E. Sorgard, and H. L. Behrens, “A historical reconstruction of ships’ fuel consumption and emissions,” J. Geophys. Res. 112(D12), D12301 (2007). [CrossRef] .
41. J. T. Wilkerson, M. Z. Jacobson, and A. Malwitz, “Analysis of emission data from global commercial aviation: 2004 and 2006,” Atmos. Chem. Phys. 10(13), 6391–6408 (2010). [CrossRef] .
42. D. F. Baker, R. M. Law, and K. R. Gurney, “TransCom 3 inversion intercomparison: Impact of transport model errors on the interannual variability of regional CO2 fluxes, 1988–2003,” Global Biogeochem. Cycles 20(1), 1988–2003 (2006). [CrossRef] .
43. N. S. Pougatchev, B. J. Connor, and C. P. Rinsland, “Infrared measurements of the ozone vertical distribution above Kitt Peak,” J. Geophys. Res. 100(D8), 16689 (1995). [CrossRef] .
44. C. D. Rodgers and B. J. Connor, “Intercomparison of remote sounding instruments,” J. Geophys. Res. 108(D3), 4116 (2003). [CrossRef] .
45. M. Schneider, S. Barthlott, and F. Hase, “Groundbased remote sensing of tropospheric water vapour isotopologues within the project MUSICA,” Atmos. Meas. Tech. 5(12), 3007–3027 (2012). [CrossRef] .
46. E. Kalnay, M. Kanamitsu, R., and W. Kistler, “The NCEP/NCAR 40-Year Reanalysis Project,” Bull. Amer. Meteor. Soc. 77(3), 437–471 (1996). [CrossRef]
47. L. S. Rothman, I. E. Gordon, and Y. Babikov, “The HITRAN2012 molecular spectroscopic database,” J. Quant. Spectrosc. Radiat. Transfer 130, 4–50 (2013). [CrossRef] .
48. M. Schneider, A. Wiegele, and S. Barthlott, “Accomplishments of the MUSICA project to provide accurate, long-term, global and high-resolution observations of tropospheric {H2O; δD} pairs – a review,” Atmos. Meas. Tech. 9(7), 2845–2875 (2016). [CrossRef] .
49. B. J. Connor, H. Boesch, and G. Toon, “Orbiting Carbon Observatory: Inverse method and prospective error analysis,” J. Geophys. Res. 113(D5), D05305 (2008). [CrossRef] .
50. K. W. Thoning, P. P. Tans, and W. D. Komhyr, “Atmospheric carbon dioxide at Mauna Loa Observatory: 2. Analysis of the NOAA GMCC data,” J. Geophys. Res. 94(D6), 8549–8565 (1989). [CrossRef] .
51. H. Ohyama, S. Kawakami, and T. Tanaka, “Observations of XCO2 and XCH4 with ground-based high-resolution FTS at Saga, Japan, and comparisons with GOSAT products,” Atmos. Meas. Tech. 8(12), 5263–5276 (2015). [CrossRef] .
52. Y. Wang, J. W. Munger, and S. Xu, “CO2 and its correlation with CO at a rural site near Beijing: Implications for combustion efficiency in China,” Atmos. Chem. Phys. 10(18), 8881–8897 (2010). [CrossRef] .
53. O. Schneising, M. Reuter, and M. Buchwitz, “Terrestrial carbon sink observed from space: variation of growth rates and seasonal cycle amplitudes in response to interannual surface temperature variability,” Atmos. Chem. Phys. 14(1), 133–141 (2014). [CrossRef] .
54. Y. Yin, P. Ciais, and F. Chevallier, “Changes in the response of the Northern Hemisphere carbon uptake to temperature over the last three decades,” Geophys. Res. Lett. 45(9), 4371–4380 (2018). [CrossRef] .
55. S. Piao, Z. Liu, and Y. Wang, “On the causes of trends in the seasonal amplitude of atmospheric CO2,” Glob. Change Biol. 24(2), 608–616 (2018). [CrossRef] .
56. N. York, N. M. Evensen, and M. L. Martínez, “Unified equations for the slope, intercept, and standard errors of the best straight line,” Am. J. Phys. 72(3), 367–375 (2004). [CrossRef] .
57. Y. Jing, T. Wang, and P. Zhang, “Global Atmospheric CO2 Concentrations Simulated by GEOS-Chem: Comparison with GOSAT, Carbon Tracker and Ground-Based Measurements,” Atmosphere 9(5), 175 (2018). [CrossRef]
58. Y. Wang, N. M. Deutscher, and M. Palm, “Towards understanding the variability in biospheric CO2 fluxes: using FTIR spectrometry and a chemical transport model to investigate the sources and sinks of carbonyl sulfide and its link to CO2,” Atmos. Chem. Phys. 16(4), 2123–2138 (2016). [CrossRef] .
59. J. Messerschmidt, N. Parazoo, and D. Wunch, “Evaluation of seasonal atmosphere–biosphere exchange estimations with TCCON measurements,” Atmos. Chem. Phys. 13(10), 5103–5115 (2013). [CrossRef] .
60. D. Wunch, P. O. Wennberg, and J. Messerschmidt, “Ground-based FTIR observation of hydrogen chloride (HCl) over Hefei, China, and comparisons with GEOS-Chem model data and other ground-based FTIR stations data,” Opt. Express 28, 8041–8055 (2013). [CrossRef] .
61. H. Yin, Y. Sun, and C. Liu, “Ground-based FTIR observation of hydrogen chloride (HCl) over Hefei, China, and comparisons with GEOS-Chem model data and other ground-based FTIR stations data,” Opt. Express 28, 8041–8055 (2013). [CrossRef] .
62. D. Wunch, G. C. Toon, and V. Sherlook, “The Total Carbon Column Observing Network’s GGG2014 Data Version,” Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA, (2015).
63. G. C. Toon “Telluric line list for GGG2014. TCCON data archive, hosted by the Carbon Dioxide Information Analysis Center,” Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA, (2014).
64. Y. Oh, S. T. Kenea, and T. Goo, “Characteristics of greenhouse gas concentrations derived from ground-based FTS spectra at Anmyeondo, South Korea,” Atmos. Meas. Tech. 11(4), 2361–2374 (2018). [CrossRef]
65. M. F. Schibig, E. Mahieu, and S. Henne, “Intercomparison of in situ NDIR and column FTIR measurements of CO2 at Jungfraujoch,” Atmos. Chem. Phys. 16(15), 9935–9949 (2016). [CrossRef] .