Constraining China's land carbon sink from emerging satellite CO2 observations: Progress and challenges

Land carbon sink is a vital component for the achievement of China's ambitious carbon neutrality goal, but its magnitude is poorly known. Atmospheric observations and inverse models are valuable tools to constrain the China's land carbon sink. Space-based CO₂ measurements from satellites form an emerging data stream for application of such atmospheric inversions. Here, we reviewed the satellite missions that is dedicated to the monitoring of CO₂ , and the recent progresses on the inversion of China's land carbon sink using satellite CO₂ measurements. We summarized the limitations and challenges in current space platforms, retrieval algorithms, and the inverse modeling. It is shown that there are large uncertainties of contemporary satellite-based estimates of China's land carbon sink. We discussed future opportunities of continuous improvements in three aspects to better constrain China's land carbon sink with space-based CO₂ measurements.

Toward High Precision XCO2 Retrievals From TanSat Observations: Retrieval Improvement and Validation Against TCCON Measurements

TanSat is the 1st Chinese carbon dioxide (CO2) measurement satellite, launched in 2016. In this study, the University of Leicester Full Physics (UoL-FP) algorithm is implemented for TanSat nadir mode XCO2 retrievals. We develop a spectrum correction method to reduce the retrieval errors by the online fitting of an 8th order Fourier series. The spectrum-correction model and its a priori parameters are developed by analyzing the solar calibration measurement. This correction provides a significant improvement to the O2 A band retrieval. Accordingly, we extend the previous TanSat single CO2 weak band retrieval to a combined O2 A and CO2 weak band retrieval. A Genetic Algorithm (GA) has been applied to determine the threshold values of post-screening filters. In total, 18.3% of the retrieved data is identified as high quality compared to the original measurements. The same quality control parameters have been used in a footprint independent multiple linear regression bias correction due to the strong correlation with the XCO2 retrieval error. Twenty sites of the Total Column Carbon Observing Network (TCCON) have been selected to validate our new approach for the TanSat XCO2 retrieval. We show that our new approach produces a significant improvement on the XCO2 retrieval accuracy and precision when compared to TCCON with an average bias and RMSE of -0.08 ppm and 1.47 ppm, respectively. The methods used in this study can help to improve the XCO2 retrieval from TanSat and subsequently the Level-2 data production, and hence will be applied in the TanSat operational XCO2 processing.

Effect of Structural Uncertainty in Passive Microwave Soil Moisture Retrieval Algorithm

Passive microwave sensors use a radiative transfer model (RTM) to retrieve soil moisture (SM) using brightness temperatures (T_B) at low microwave frequencies. Vegetation optical depth (VOD) is a key input to the RTM. Retrieval algorithms can analytically invert the RTM using dual-polarized T_B measurements to retrieve the VOD and SM concurrently. Algorithms in this regard typically use the τ-ω types of models, which consist of two third-order polynomial equations and, thus, can have multiple solutions. Through this work, we find that uncertainty occurs due to the structural indeterminacy that is inherent in all τ-ω types of models in passive microwave SM retrieval algorithms. In the process, a new analytical solution for concurrent VOD and SM retrieval is presented, along with two widely used existing analytical solutions. All three solutions are applied to a fixed framework of RTM to retrieve VOD and SM on a global scale, using X-band Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) T_B data. Results indicate that, with structural uncertainty, there ensues a noticeable impact on the VOD and SM retrievals. In an era where the sensitivity of retrieval algorithms is still being researched, we believe the structural indeterminacy of RTM identified here would contribute to uncertainty in the soil moisture retrievals.

Assessment of Retrieved N₂O, NO₂, and HF Profiles from the Atmospheric Infrared Ultraspectral Sounder Based on Simulated Spectra

The Atmospheric Infrared Ultraspectral Sounder (AIUS), the first high-resolution (0.02 cm^−1) solar occultation sounder, aboard GF5, was launched in May 2018 from China. However, relevant studies about vertical profiles of atmospheric constituents based on its operational data were not conducted until half a year later. Due to an urgent need for Hin-orbit tests, the real spectra (called reference spectra hereafter) were substituted with simulated spectra calculated from the reference forward model (RFM) plus different random noises at different altitudes. In the generation process of the reference spectra for N₂O, NO₂, and HF species, ACE-FTS (Atmospheric Chemistry Experiment⁻Fourier Transform Spectrometer instrument on the SCISAT satellite) level 2 products replace corresponding profiles included in the atmospheric background profiles. The optimal estimation method is employed to extract N₂O, NO₂, and HF profiles in this study. Comparing the retrieved results with ACE-FTS level 2 products, the relative deviations for these three species are calculated. For N₂O, the average relative deviation is less than 6% at altitudes below 25 km, while larger deviations are observed in the range of 25⁻45 km, with the maximum being at ~25%. Additionally, the difference for NO₂ is less than 5% in the 20⁻45 km range, with a larger discrepancy found below 20 km and above 45 km; the maximum deviation reaches ±40%. For HF, the relative deviation is less than 6% for all tangent heights, implying satisfactory retrieval. The vertical resolution, averaging kernel, and number of degrees of freedom are used to assess the retrieval algorithm, which indicate that the retrieved information content is much more attributable to the reference spectra contribution than to the a priori profile. Finally, a large number of retrieval tests are performed for N₂O, NO₂, and HF in selected areas covering the Arctic region, northern middle latitude, tropics, southern middle latitude, and Antarctic region, and reliable results are obtained. Thus, to a great extent, the algorithm adopted in the AIUS system can process retrievals reliably and precisely.

Reducing the Discrepancy Between ASTER and MODIS Land Surface Temperature Products

Human-induced global warming has significantly increased the importance of satellite monitoring of land surface temperature (LST) on a global scale. The MODerate-resolution Imaging Spectroradiometer (MODIS) provides a 1-km resolution LST product with almost daily coverage of the Earth, invaluable to both local and global change studies. The Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) provides a LST product with a high spatial resolution of 90-m and a 16-day recurrent cycle, simultaneously acquired at the same height and nadir view as MODIS. ASTER and MODIS are complementary in resolution, offering a unique opportunity for scale-related studies. ASTER and MODIS LST have been widely used but the errors in LST were mostly disregarded. Correction of ASTER-to-MODIS LST discrepancies is essential for studies reliant upon the joint use of these sensors. In this study, we compared three correction approaches: the Wan et al.'s approach, the refined Wan et al.'s approach, and the generalized split window (GSW) algorithm based approach. The Wan et al.'s approach corrects the MODIS 1-km LST using MODIS 5-km LST. The refined approach modifies the Wan et al.'s approach through incorporating ASTER emissivity and MODIS 5-km data. The GSW algorithm approach does not use MODIS 5-km but only ASTER emissivity data. We examined the case over a semi-arid terrain area for the part of the Loess Plateau of China. All the approaches reduced the ASTER-to-MODIS LST discrepancy effectively. With terrain correction, the original ASTER-to-MODIS LST difference reduced from 2.7±1.28 K to -0.1±1.87 K for the Wan et al.'s approach, 0.2±1.57 K for the refined approach, and 0.1±1.33 K for the GSW algorithm based approach. Among all the approaches, the GSW algorithm based approach performed best in terms of mean, standard deviation, root mean square root, and correlation coefficient.