Next Article in Journal
Characterizing Spatiotemporal Patterns of Land Deformation in the Santa Ana Basin, Los Angeles, from InSAR Time Series and Independent Component Analysis
Previous Article in Journal
Strategies for the Storage of Large LiDAR Datasets—A Performance Comparison
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

XCO2 and XCH4 Reconstruction Using GOSAT Satellite Data Based on EOF-Algorithm

1
School of Instrumentation Science and Optoelectronic Engineering, Beihang University, Beijing 100191, China
2
Research Institute of Frontier Science, Beihang University, Beijing 100091, China
3
State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China
4
School of Transportation Science and Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(11), 2622; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14112622
Submission received: 11 February 2022 / Revised: 11 May 2022 / Accepted: 21 May 2022 / Published: 31 May 2022
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
The Greenhouse Gases Observing Satellite (GOSAT) can help to ascertain the global distribution of carbon dioxide (CO2) and methane (CH4), and how the sources and sinks of these gases vary by season, year, and location. However, the data provided by the GOSAT level 2 and 3 products have certain limitations due to their lack of spatial and temporal information; even with the application of the kriging geostatistical method on the level 2 products, the processing algorithms still need further upgrades. In this study, we apply an empirical orthogonal function (EOF)-based method on the GOSAT L3 products (137 images, from January 2010 to May 2021) to estimate the column average of carbon dioxide and methane (XCO2–XCH4) within the entire Earth. The reconstructed results are validated against the Total Carbon Column Observing Network (i.e., TCCON), with 31 in situ stations, and GOSAT L4B column-averaged data, using 107 layers. The results show an excellent agreement with the TCCON data and exhibit an R-squared coefficient of 0.95 regarding the CO2 measurements and 0.86 regarding the CH4 measurements. Therefore, this methodology can be incorporated into the processing steps used to map global greenhouse gases.

1. Introduction

CO2 is the most important greenhouse gas associated with human activities and the second most important gas in global warming after water vapor. Since the industrial revolution, a notable increase in global CO2 concentrations has been observed, rising from 280 ppm (parts per million) in the middle of the 19th century to over 400 ppm as of 2015 [1,2]. Several atmospheric research centers noted this rise; based on preliminary analysis, the global average atmospheric CO2 in 2020 was 412.5 ppm [3], and in 2021, it was 415.7 ppm, based on monthly mean NOAA reports [4,5].
CH4 is a strong greenhouse gas and plays an important role in determining the oxidation capacity of the troposphere, with a global warming potential on a 100-year time-scale 28–34 times that of CO2 [6]. According to current observations, in 2015, the globally averaged content of CH4 was already1845 ppb [7]. Scientific debate continues in the effort to explain the atmospheric CH4 trend observed over the past couple of decades. Records from surface sites reveal a plateau from 2000 to 2007 and a resumed increase after 2007 [8,9]. Amongst the varied surface sources of CH4, the largest are natural wetlands, agriculture, livestock, biomass burning, waste, and fossil fuel production; whereas the primary sink is the OH-radical content in the atmosphere [10,11]. Various hypotheses have been offered that attempt to attribute the behavior in the global growth rate to a particular component or mechanism [12,13,14,15], but currently, there is no consensus within the scientific community. Many of these studies have utilized satellite observations of atmospheric CH4 and have shown the increasing ability of such measurements to characterize global and regional surface CH4 fluxes [16].
Greenhouse gases cause environmental and health issues. They cause climate change by trapping heat, which affects flora and fauna. The climate change caused by greenhouse gas emissions also contributes to wildfires, weather, and droughts. In the last decade, CO2 and CH4 have increased constantly. Localizing the highest concentrations globally and locally is important to alert particular areas in order to determine solutions through which to reduce the emission of these gases. The pandemic situation caused by COVID-19 forced the world to lock down, and this study and others [17] show that greenhouse gases did not reduce significantly; therefore, the search for alternatives to these gases is the next step for humanity. Building this scientific understanding to meet the extraordinary climate mitigation challenge requires frequent, robust, transparent, and traceable data sets and methods that can be scrutinized and replicated [18].
Ground stations, such as TCCON, have been established to accurately measure atmospheric greenhouse gases, such as CO2, CO, N2O, and CH4 [2,19]. A significant amount of effort has been expended and measurements of the greenhouse gases CO2 and CH4 by TCCON have resulted in an accuracy of 0.25% under clear atmospheric conditions [20]. However, due to the limited spatial coverage and uneven distribution of TCCON sites, it is not possible to accurately measure CO2 quantities at subcontinental and regional spatial scales [21]. The reliability of the results from ground stations needs to be enhanced, while reliable TCCON results are only available in North America and Europe due to a lack of information in South America, Africa, and Antarctica. Studies suggest that satellite observations, with their lower accuracy than ground-based measurements but higher spatial coverage, may help to improve CO2 measurements [22,23,24]. Current estimates of CO2 and CH4 from existing satellite instruments have spatially consistent biases that hinder their use for monitoring carbon–climate feedback [25]. This paucity of information has hindered the modeling, prediction, and management of air-quality studies.
GOSAT is a Japan Aerospace Exploration Agency (JAXA) mission designed to monitor the CO2 and CH4 on a subcontinental scale more accurately and to assist environmental administrations in evaluating the carbon balance of the land ecosystem and performing assessments of regional emissions and absorptions. By analyzing the GOSAT data, scientists will accumulate new knowledge on the global distribution and temporal variation of greenhouse gases, as well as the global carbon cycle and its influence on the climate. These new findings will be utilized to predict future climate change and assess its impact.
However, GOSAT has experienced many events in past decade. For example, GOSAT has experienced several incidents resulting in instrument anomalies [26], as follows. In 2014–2015, a solar paddle incident resulted in a temporary instrument shutdown, which caused missing information in December 2014 and January 2015. In 2015 Cryocooler shut down and restarted. In May 2018, there was an incident with and restart of the Command and Data Management System (CDMS). There was a solar-paddle-rotation incident, and a thermal and near-infrared sensor for carbon observation (TANSO) restart, which caused a loss of information, in November 2018. Finally, in May 2020, there was an accelerated zero path difference (ZPD) shift from May 31 until the recovery on June 18. This information was recorded in the report on the GOSAT’s status.
On 23 January 2014, GOSAT/Ibuki had been in orbit for five years. The spacecraft and its payload were operating nominally. With a nominal design life of 5 years and the spacecraft in good health, the mission’s duration was extended (by JAXA, National Institute of Education Sciences (NIES), and Ministry of the Environment (MOE)) for another 4 years [27,28]. In September 2013, the first estimates of the global distribution of CO2 surface fluxes were estimated and published using the total column CO2 measurements retrieved by the SRON-KIT RemoTeC algorithm from GOSAT. The global source-sink estimates of CO2 were based on data from over eighteen months, from 1 June 2009 to 1 December 2010.
It is well known that GOSAT L3 products are generated by interpolating, extrapolating, and smoothing the Fourier transform spectrometer (FTS) level 2 column-averaged CO2 and CH4 mixing ratios. The GOSAT L3 products still have missing data, even with the processing of L2 products with the kriging interpolation method. This method is limited due to the uncertainty generated by large areas of null information; these areas with a high level of uncertainty are therefore ignored. Given the problem of the loss of relevant information, it is necessary to search for alternative solutions. One potential solution is to apply interpolation on a spatial–temporal matrix of the GOSAT data from 2010 until 2021, to retrieve data using information related to space and time. This methodology can be implemented to obtain a reliable amount of interpolated information to study and anticipate future environmental changes and study their effects in areas that have low levels of information concerning GOSAT L3 product records. The space-based remote sensing of the CO2 and CH4 column-averaged dry-air mole fractions XCO2 and XCH4 has the potential to provide global observational constraints on CO2 and CH4 fluxes across the surface–atmosphere boundary and to offer insights into the related biogeochemical cycles. However, the measurement of XCO2 and XCH4 imposes challenging accuracy requirements. The target accuracies for regionally and weekly averaged XCO2 and XCH4 are in the order of a few tenths of a percent based on background concentrations of roughly 385 ppm and 1.8 ppm, respectively, and the requirement for XCH4 is less stringent than for XCO2 [29].
Previous studies applying the EOF algorithm to GOSAT products, such as [30], conclude that the implementation of the EOF-based algorithm for XCO2 retrievals from GOSAT observations should be considered, and the development of very fast XCO2 retrieval algorithms to process the huge amounts of current and future satellite observation data remains of interest [26]. However, separating the relative importance of each factor is often complicated, as some of the influencing factors are interdependent due to scale dependence and other factors. It is important to note that this study should be conducted comprehensively, as it is necessary to analyze more data related to the XCO2 and XCH4 values from the satellite, ground stations, and even climate models, since the values obtained in this study are the products of spatial–temporal interpolation, with a high confidence value, of only the pixels covering the areas where the TCCON ground stations were installed. Several methodologies can be applied to the GOSAT L3 dataset, such as the interpolation of data based on a harmonic trend of each pixel in a time series, but it is important to emphasize that not at all of these methodologies use spatial and temporal information to interpolate data. Therefore, future studies related to this methodology will require the analysis of several databases to improve the reliability of the results.
The objectives of this study are to: (1) evaluate the spatial–temporal patterns of greenhouse gases based on the EOF-GOSAT L3 datasets; (2) analyze the feasibility of the EOF-based algorithm applied to the dataset; and (3) provide a possible approach to integrating an additional processing step into the GOSAT L3 products.

2. Materials and Methods

It is generally known that the Fourier-transform spectrometer–shortwave infrared (FTS-SWIR) Level 3 data products are constructed monthly by interpolating, extrapolating, and smoothing the FTS SWIR Level 2 column-averaged CO2 and CH4 mixing ratios. Kriging, a geostatistical calculating approach, is used in L2 products to obtain L3 products. The data are divided into 2.5-degree cells (GOSAT webpage). The standard error is stored in L3 products and almost all missing information is located in the South and North poles. Due to this problem, a series of stages are developed to obtain missing information. (1) The EOF technique is used to fill in missing data from geophysical fields in image reconstruction due to missing information in GOSAT L3 deliverables. (2) Following image reconstruction, correlation analysis utilizing kernel plots is used to evaluate in situ stations with reconstructed images. Next, the temporal and spatial analysis of CO2 and CH4 is assessed, with the most heavily affected places identified and their growth rates indicated.

2.1. Study Area

Due to the variety of climates, types of vegetation cover, and human activity intensity around the world, it is difficult to derive CO2 concentrations using a single metric. To increase the derived accuracy of the CO2 concentration, we divided the world into six regions, according to the following continents: Africa, Asia, Australia, Europa, North America, and South America. The Antarctic continent was removed from this analysis due to the missing GOSAT information. Moreover, the TCCON project does not have in situ stations in this region and the interpolated results are not sufficient to corroborate the results. For each continent, we focus on two greenhouse gas concentrations (CO2 and CH4).

2.2. Data Collection

Two datasets were obtained. GOSAT satellite image data (which are freely available) from January 2010 to May 2021 were downloaded from GOSAT Data Archive Service (GDAS) website. These global data were collected under cloud-free conditions for estimating column-averaged dry-air mole fractions of major greenhouse gases, such as CO2 and methane. TCCON data were also collected. Both datasets are described in the following sections.

2.2.1. GOSAT

The main role of the GOSAT project is to estimate emissions and absorptions of the greenhouse gases on a subcontinental scale and to help ecological organizations to assess the carbon equilibrium of the land biological system and make appraisals of provincial outflows and retentions. GOSAT orbits at an altitude of approximately 666 km, with a 10.5-km spatial resolution and a three-day temporal resolution [29]. The observation instrument on GOSAT is the thermal and near infrared sensor for carbon observation (TANSO), which is composed of two subunits: the Fourier-transform spectrometer (FTS) and the Cloud and Aerosol Imager (CAI). Thermal and near-infrared sensor for carbon observation–Fourier transform Spectrometer (TANSO-FTS) operates with spectral resolution in three narrow bands in the short-wavelength infrared (SWIR) region (0.76, 1.6, and 2.0 µm) and a wide thermal infrared band (5.5–14.3 µm). Thermal and near-infrared sensor for carbon observation–Cloud and Aerosol Imager (TANSO-CAI) operates with spectral resolution in four narrow bands in the near-ultraviolet-to-near-infrared region (0.38, 0.674, 0.87, and 1.6 µm). TANSO can detect optically thick clouds inside the TANSO-FTS instantaneous field of view (IFOV) and correct the effect of aerosols in the TANSO-FTS spectrum data [24,28,31].
This study used the FTS SWIR Level 3 data products generated from 1 January 2010 to 31 May 2021, freely downloaded from the GOSAT Data Archive Service (GDAS), which was collected under cloud-free conditions for estimating column-averaged dry-air-mole fractions of major greenhouse gases, such as CO2 and methane. GOSAT Level 2 does not have data for January 2015, and GOSAT Level 3 does not have data for June and December 2014 or January and September 2015.

2.2.2. TCCON

The TCCON is a worldwide observation system of ground-based Fourier transform spectrometers (FTSs) that records direct solar spectra in near-infrared. It is associated with the network for the GAW program and the Detection of Atmospheric Composition Change Infrared Working Group (NDACC-IRWG) [2]. It is a ground-based network of FTSs designed to retrieve precise and accurate column abundances of CO2, CH4, N2O, and CO from near-infrared (NIR) solar absorption spectra. It provides accurate measurements of studied greenhouses (CO2 and CH4) at more than a dozen operational sites worldwide. The sites are generally located at low altitudes, where the biosphere/ocean exchange has the most impact on the atmospheric concentrations [32]. The GGG2014 TCCON [33] data used in this study can be obtained from the TCCON Data Archive (http://tccon.ornl.gov/ (accessed on 15 January 2022)) [34], and the data, provided by the Carbon Dioxide Information Analysis Center (CDIAC), range from 2004 in some stations until 2020. The precision of the latest version has been improved by comparing it with that of GGG2014 by correcting laser-sampling errors and adding more aircraft and AirCore profiles as bias correction data [35].
TCCON observations are retrieved to provide a long-term, almost continuous time-series to assist vertical verification as a transfer standard between in situ measurements and satellite observations and to help better understand carbon-cycle processes [36]. The network shares methodologies, technical issues, and membership with the NDACC FTIR Working Group. The network aims to provide temporally dense, globally consistent data coverage network-wide by using a common retrieval technique. The retrieved columns are tied to the ground-based in situ networks by scaling to (in situ) balloon and aircraft profiles [37].
As shown in Figure 1, all of the TCCON stations are represented with blue points, and almost all of the stations are located in the northern hemisphere, especially in developed countries. On the other hand, the developing countries do not have sufficient information, and the Antarctica region has many empty gaps in GOSAT L3 product information. Therefore, the results from this region were excluded, since they are the results of interpolation without corroboration.
The total error for a single measurement of XCO2 at the Lamont site is 0.2 percent or less for SZAs less than 83° [37]. The largest errors are from instrument alignment errors, assuming a linear spectral continuum across a large spectral band, and a priori profile uncertainties.
In Figure 1, all the locations that are part of the in-situ stations of the TCCON network are visible, and Table 1, presents all the coordinates and the names of their respective in situ stations. Furthermore, we estimated the R-squared and b linear slope in order to analyze the variability of results of each in situ station.

2.3. Data Processing

The ideal TCCON–satellite data pair should consist of measurements at the same place during the same time. However, to find a sufficient number of co-located measurements to enable a robust statistical analysis, several spatial–temporal criteria were used in previous validations. Therefore, the statistical analyses were based on the individual data pairs or daily averaged data pairs, and all data pairs were assumed to be of equal weight in this study.
The EOF algorithm was applied to acquire missing greenhouse gas information in GOSAT FTS SWIR Level 3 data. Next, a monthly global mapping of XCO2 and XCH4 was generated based on the reconstructed images, which were validated with TCCON data, as well as density scatter plots to compare the results.
The EOF technique was used to examine the trend of these greenhouse gases over the previous ten years, since CO2 readings are more consistent than CH4 values. The temporal and spatial means of the data were eliminated before the original data set was placed in an n-dimensional spatial–temporal matrix. Next, the first EOF mode was calculated using the singular value decomposition (SVD) technique, which was used to infer a new estimate for the missing data. This procedure was repeated until convergence was obtained for the values given to the missing data with the first EOF mode [27]. Finally, a few charts, such as Hovmöller, power wavelet spectrum, time-series analysis, and decomposition, are shown in Section 3 to represent temporal variation. The spatial analysis, on the other hand, was organized by continent. Despite this, fewer coherence results for the continent analysis were reported in Antarctica due to the limited information provided by GOSAT; they could not be compared to those of any of the TCCON stations because the network of ground in situ stations did not install a climate observatory in this region.

2.4. Data Interpolation

The geostatistical mapping method applied in this study is based on the EOF algorithm used to fill empty gaps using temporal data, this algorithm has been developed at the Geo-Hydrodynamics and Environmental Research (GHER) [38].
Concerning spatial resolution, the L3 product was generated for meshes with dimensions of 2.5° in both latitude and longitude. This grid size seems reasonable when considering the GOSAT lattice observation distances and the present rate at which FTS SWIR L2 data are retrieved [39].
While all interpolation methods have a degree of unreliability, the EOF algorithm exhibits some major advantages compared with many others [40], including the minimization of errors from iterative processing, parameter-free processing, and speed.
A detailed description of the EOF processing is presented in [41], and the libraries and repositories to use are available on GitHub and the Comprehensive R Archive Network (CRAN). EOF algorithm is a process that fills gaps in data fields by repeatedly decomposing them using SVD until the optimal answer when compared to a set of reference values (non-gaps) is discovered. This is accomplished by gradually including more EOFs into the gappy location reconstruction until the error minimization converges. We first show the approach using a synthetic example (Figure 2), in which: (1) a data set is constructed, (2) noise is introduced (based on a noise/signal setting), (3) gaps are made (based on a fraction of gappiness setting), and, lastly, the gaps are interpolated using the DINEOF method. Figure 3 depicts the iterations used to fill in the gaps as the number of EOFs increases. Using the first six EOFs, an optimal solution was found, whereas adding a seventh EOF did not enhance the root mean squared error (RMSE) [42].
Nine mixed signals were used in the synthetic example. The EOFs of the “true” field (black points) and the optimum solution of the gappy data are shown in Figure 4 (red lines).
The comparison of TCCON and GOSAT XCO2 data was based on the same spatially resolved mean monthly data set. The monthly global distributions of TCCON XCO2 (2.5° × 2.5°) were reproduced based on GOSAT Level 3 SWIR products. The reconstructed monthly global distributions of TCCON XCO2 were generated by averaging all measurements in each grid corresponding to the in-situ stations for each month. We mainly used the GOSAT SWIR Level 3 products because the spatial coverage of the reconstructed GOSAT Level 2 global monthly XCO2 distributions (using the same method as TCCON) was smaller than that of the L3 products, as shown in Figure 2. Moreover, the spatial coverage of GOSAT Level 3 XCO2 when applying the EOF algorithm increased in the number of pixels, as can be seen in Figure 5, where each pixel has 2.5 × 2.5 grid degrees as a unit. As can be seen, the number of pixels increased from 2500 to approximately 3500 approximately from 2015, which indicates that the number of pixels with missing information has increased in GOSAT L3 products since 2015, and the months with a high amount of information correspond to July to October in almost every year. It is important to note the total number of pixels in each month layer that is equal to 10,368 pixels; in some cases, the number of reconstructed pixels is approximately 4500, especially for the 2010–2014 layers.

3. Results and Discussion

Three analyses are evaluated in this study. First of all, the GOSAT data validation is presented and discussed in Section 3.1. Next, the data validation is evaluated, comparing the GOSAT interpolated results and the TCCON data; this analysis is discussed in Section 3.2. Finally, the temporal and spatial analysis is assessed in Section 3.3 and Section 3.4.

3.1. Data Validation

We directly compared the monthly mean values of the TCCON in situ station measurements with the L3 products of the GOSAT interpolated data using the EOF algorithm of the XCO2 and XCH4 data in this study. To obtain the GOSAT data, we subtracted the temporal variation in the pixel of each station from the stack raster reconstructed data, as shown in Table 2.
TCCON and GOSAT have different a priori profiles and averaging kernels, which describe the sensitivity of a retrieval algorithm to the true state throughout the atmosphere [33,34]. Recent research has demonstrated that the effect of column averaging kernels on GOSAT XCO2 is relatively negligible when compared with its measurement accuracy, and direct comparison is applicable in the validation of satellite measurements [4,35].
To obtain the results of the validation with TCCON data with GOSAT data, we used the in situ information from 31 stations within the data validation analysis. The time range and geolocation are presented in Table 2, as are the results of the linear fit (linear correlation coefficient k and R-squared coefficient) between the GOSAT L3 XCO2 and XCH4 measurements with the application of the EOF algorithm and the TCCON data in each location. Figure 6 shows the same statistics on the GOSAT measurements as Figure 7, which compares the data of the TCCON GGG2020, GOSAT L3, and GOSAT L4B column averages of XCO2 and XCH4. Figure 6 shows the data correlation results. The RMSE value for the GOSAT L3 and TCCON was 2.74, and the R-squared value was 0.95. The RMSE for the GOSAT L4 and TCCON was 1.88, and the R-squared was 0.96. The RMSE for the GOSAT L3 and GOSAT L4 was 3.02, and the R-squared was 0.94. This allows a comparison of the results of the application of the EOF algorithm on the L3 products with the TCCON station records. On the other hand, Figure 7 shows the same comparison, of the results of the application of the EOF algorithm with the TCCON records. The RMSE value for the GOSAT L3 and TCCON was 1.56 × 10−4, and the R-squared was 0.86. The RMSE for the GOSAT L4 and TCCON was 1.10 × 10−3, and the R-squared was 0.62. The RMSE for the GOSAT L3 and GOSAT L4 was 9.38 × 10−4, and the R-squared was 0.86. Other secondary comparison parameters between the average column results of the GOSAT L4B and TCCON products are visible in Figure 6 and Figure 7.
The TCCON network comprises ground-based high-resolution Fourier-transform spectrometers that perform direct measurements of solar spectra in near-infrared. There are currently 27 operational sites located across North America, Europe, Asia, and Australia. Only one station in South America is registered, and it only contains information from 2015. Some other stations are in the same situation, including “Alesund”, “Izana”, “Manaus”, “Indianapolis”, “Saga”, “Rikubetsu”, “Lauder”, “Orleans”, which only contain a low quantity of data that can be used to compare ground station measurements with satellite data. In this study, the data from 31 stations were compared with the GOSAT data; the total number of training points was 2602, to assess the density scatter plots. We reduced the TCCON information, obtaining the average value of each month, to compare with the monthly information from the GOSAT. The TCCON stations have many time ranges without values because the stations are not currently operating. On the other hand, all the GOSAT pixels of 2.5-degree resolution that contain the geolocation of each station were extracted from a stack of information by using a bilinear method.

3.1.1. Comparison of XCO2 GOSAT and TCCON Data

To validate the CO2 data used in the generation of the final global mapping, we evaluated the daily mean of the TCCON greenhouse gas readings and faced the data with a series of information extracted from pixels that were near 2.5° of the greenhouse GOSAT information.
We plotted 35 charts of TCCON information and overlays with GOSAT information, as shown in Figure 8, where the blue dotted line represents the TCCON readings and the red dotted line represents the GOSAT readings. In certain cases, in situ stations did not record for a period of time, which caused the blank spaces in the plotting of the results.

3.1.2. Comparison of XCH4 GOSAT and TCCON Data

There are signs of a potential hemispheric or latitudinal bias against TCCON in the data, although this is not apparent at all sites. For example, Karlsruhe, Lamont, Tsukuba, and Lauder all have negligible biases but span a large latitude range. It should also be noted that the uncertainty on the TCCON CH4 is approximately 4 (ppb), and for most sites, the GOSAT–TCCON difference is within this range of uncertainty; it is therefore is important to not over-interpret any signals at this scale [43].
As with the XCO2 data obtained by the in situ stations, there were some missing XCH4 readings due to some stations not operating. In the XCH4 readings, only three stations (Zugspitze, Garmisch, and Izana) registered an appreciable difference in line trends, of approximately 0.03 ppm as can be appreciated in Figure 9.

3.1.3. XCO2 and XCH4 Difference between GOSAT L3 and TCCON Data

Along with the mission extension, the reprocessing of the L1B products was completed with the v161161 algorithm for all of the GOSAT data to retrieve the XCO2 and XCH4 data with high accuracy and precision [44]. Accuracy can be defined as the degree of closeness between the information on a map and the values in the real world. On the other hand, precision refers to the level of exactness of the description of data, The expected relative accuracies were 1% for the XCO2 (4 ppm, 3-month average) and 2% for the XCH4 (34 ppb, 3-month average). The inverse algorithms were continuously optimized [45]. A validation study showed that the greenhouse gas observation sensor on the GOSAT/Ibuki, “TANSO-FTS”, observes column-averaged CO2 concentrations (hereafter, “CO2 concentrations”) with a precision (random error) of 0.5% (approximately 2 ppm) in comparison to the data acquired from ground-based observations. This study investigated the detectability of enhanced CO2 concentrations due to anthropogenic emissions by using satellite observations [46]. Therefore, when we analyzed the accuracy of the XCO2 by comparing the TCCON data with the GOSAT L3 by applying the EOF algorithm products, we found that the Bremen, Lamont, Wollongong, and Lamont in situ stations had a difference in values of up to 3 ppm. The difference expressed in Figure 10a, ΔXCO2 = XCO2_TCCON − XCO2_GOSAT_L3, helps to evaluate the precision and accuracy across both data. Moreover, regarding the difference in the XCH4, it is visible that the stations of Bremen and Ålesund had a difference in values up to 0.07 ppm. This analysis and the differences across other stations can be observed in Figure 10b.

3.2. Carbon Dioxide and Methane Growth (Global and Local Analysis)

Figure 11 shows low peaks in June–July of each year. These start to increase and reach a high peak in October–November–December, after which the values remain constant until March–April–May. The growth of XCO2, based on the global average values of the GOSAT L3 products with the application of the EOF algorithm, is determined by the difference of every two consecutive months. The increase in annual averages from 2019 to 2020 (2.5 ppm) was slightly lower than the increase observed between 2018 and 2019, but slightly higher than the average rate of increase over the previous decade (2.40 ppm per year) [47]. Therefore, the variation in the temporal distribution for each month was minimal; however, it presented a constant increase, especially from January to April and from September to December. In 2019, the lowest growth values of CO2 were recorded at the beginning of the year and its end. This information shows, once again, that one of the main greenhouse gases, CO2, continues to increase despite the measures taken worldwide in order to reduce it. This constant increase could cause the effects of global warming and the climatic changes that it entails to become more serious, placing in danger the future of humanity, as well as that of wildlife and flora. Therefore, more measures must be strengthened and implemented to prevent these levels from continuing to rise.
In this case, Figure 12 shows the growth of the XCH4, based on the global mean values of the GOSAT L3 products with the application of the EOF algorithm. The monthly growth was calculated by subtracting the current monthly value from the previous monthly value. The global mean XCH4, calculated from in situ observations, reached a new maximum of 1.889 ± 2 ppm in 2020, an increase of 11 ppm over the value for 2019 [47]. The low CH4 growth peak usually occurred in October–November–December and began to increase to a high peak in April–May–June.
On the other hand, in other studies, it is claimed that the most important factors influencing the temporal variability in CH4 emissions are the presence of mice, the building of burrows, and the connection of reduced-soil regions to aerated regions. Thus, to reduce the CH4 release from landfill cover soils, it is necessary to avoid any direct connection between the refuse body and the aerated soil sections of the atmosphere, which can be caused not only by animal activities, but also by specific soil-structure-forming processes. Moreover, under the given refuse and climatic conditions, a well-aerated soil cover layer of at least 70 cm is recommended [48].
From 2010 to 2021, CO2 and CH4 significantly increased each year, at rates of approximately 2.50 ppm/year and 0.012 ppm/year respectively. These rate estimations are based on linear arrangements in time-series analyses of median values. Significant changes were found in the Northern and Southern Hemispheres. The CO2 and CH4 growth were similar to the value of the global means of 2.4 ppm/year and 8 ppb/year, respectively (for the last decade, 2011–2021), as reported in [47].

3.3. Global Temporal Analysis

Real-time satellite observations can objectively determine changes in atmospheric CO2 concentrations caused by regional emissions. However, the satellite observation mode, cloud cover, atmospheric aerosols, and surface albedo, cam affect the satellite signal. Therefore, monthly data are discontinuous in time and space. The GOSAT satellite acquisition retrieval data does not have global terrestrial coverage, and it only has a small amount of data in some areas. Overall, the data cover only about half of the high latitudes in the southern hemisphere, and Antarctica lacks observational data. Therefore, we gridded the annual GOSAT satellite CO2 data to obtain global coverage.
Figure 13 and Figure 14 show the variation in CO2 and CH4 from January 2010 to May 2020. The plots are arranged by year and meteorological season, based in the Northern Hemisphere, it can be concluded that the highest concentration of CO2 is in the Northern Hemisphere. The highest concentration of CH4 is also in the Northern Hemisphere. The Southern Hemisphere presents lower values in almost all the months.
An important point to highlight is the behavior of greenhouse gases throughout the entire year. Regarding CO2, higher concentrations are common in spring and winter, while the lowest concentrations occur in autumn. Regarding methane, the highest concentrations occur in winter and the lowest occur in summer; the highest concentrations around the world are in India, the south of China, and the middle part of Africa. Additionally, Table 3 presents the CO2 and CH4 growth rates per year, based on the GOSAT and NOAA reports; the NOAA reports in this table are available from Ed Dlugokencky, at NOAA/GML (gml.noaa.gov/ccgg/trends_ch4/ (accessed on 20 January 2022)), and the GOSAT monthly reports are available from www.gosat.nies.go.jp (accessed on 20 January 2022).
Hovmöller plots (latitudinal mean versus time) clarify the zonal propagation of synoptic-scale waves in a latitude band. Thus, the time–latitude cross-sections are the Hovmöller plot of the climate anomalies. The plots various parameters at the Earth’s surface [49]. In Figure 15 and Figure 16 the “y” axis represents time, covering the period from January 2010 to May 2021, and the “x” axis is the latitude. The study developed in [20] shows the global average distribution of atmospheric CO2 with harmonic variation in latitude and highlights that the minimum values of CO2 concentration are in the Southern Hemisphere, as demonstrated by the Hovmöller plots in Figure 15a,b.
Figure 15 shows Hovmöller plots of the EOF GOSAT L3–GOSAT L4B XCO2 total column amount. The GOSAT distribution behaved as expected, showing an increase over time between 2010 and 2019 superimposed over a north–south gradient and a regular seasonal cycle. The simulated results of the L4B products exhibited very similar characteristics and were in very good agreement with the GOSAT data. The difference between the two datasets (lower panel) shows that although there was a small offset between the two, there was very good consistency over time. The GOSAT L3 and L4B data seemed to agree slightly better during the peak of the seasonal cycle, particularly in the tropics; unreliable data were obtained near the 80–90°N latitudinal area, reaching RMSE values of 6 ppm.
On the other hand, Figure 16 shows the Hovmöller plots (latitude mean versus time, respectively) for the GOSAT CH4. The GOSAT distribution behaved as expected, showing an increase over time between 2010 and 2021, increasing at the North and South Poles. At very high northern latitudes, the CH4 concentration was slightly lower at the co-ordinates corresponding to the section 40°S to 90°S; in sections 90°N to 40°S, the higher methane values were concentrated.
The specific area corresponding to Antarctica cannot be analyzed in this section due to missing information in the GOSAT SWIR L3 data, and the EOF temporal interpolation was unable to complete the information from this area because the number of pixels collected from January 2010 to May 2021 was zero. Figure 16a,b shows the temporal analysis of the mean XCH4 values of the EOF–GOSAT SWIR L3 and GOSAT L4B data. On the other hand, Figure 16c shows that the highly reliable data corresponding to the area from the 60°S to 50°N can be used, because the RMSE values are from 0 to 0.02 ppm.
Figure 17 and Figure 18 depict the time-frequency decomposition of the data using Morlet wavelet transform, highlighting the frequency of the events and presenting a harmonic and linear fit of XCO2 and XCH4 to analyze the trend of these greenhouse gases. To construct the harmonic curve fit of a periodic signal, ordinary least squares regressions using sinusoidal and cosine linked curves on the time series data were used. The underlying algorithm is based on [50]. Formula (1) shows the harmonic fit used in Figure 17d and Figure 18d. The linear fit is based on ordinary least squares (OLS), which is a method that is used to quantify the evaluation of the different regression lines.
x t = β 0 + β 1 2 π t + i = 1 n [ β 2 i c o s i 2 π t T + β 3 i s i n i 2 π t T ]
where:
  • t = Time difference in fractal years.
  • T = Length of the time period.
  • n = Total number of harmonics.
  • i = Current harmonic.
  • β x = Independent coefficients derived from OLS regression.
The global mean of the XCO2 during the studied period continuously rose. A seasonal fluctuation in this parameter, visible in Figure 18a, was extracted from the time series XCO2 GOSAT data Figure 18d (black dotted line) to analyze the frequency of these events, which showed a fluctuation between 1.47 to 1.62 ppm. This was evaluated in the Morlet wavelet power spectrum and the global wavelet spectrum Figure 18b,c. In Figure 18b, the Morlet wavelet power spectrum shows a decrease and increment in this gas every six months.
Many meteorological organizations, such as the World Meteorological Organization (WMO), have reported that the global monthly mean CO2 concentration based on data obtained at surface-level monitoring sites has exceeded 400 ppm. However, this was the first time that the whole-atmospheric CO2 mean exceeded 400 ppm, as monitored by the GOSAT, which can observe CO2 concentrations from the Earth’s surface to the top of the atmosphere (about 70 km). This means that CO2 concentrations are increasing not only at the global surface but also in the global atmosphere [51].
Figure 17. Wavelet analysis of the XCO2 based on GOSAT data. (a) Time series (solid black line) and inverse wavelet transform (solid gray line). (b) Normalized wavelet power spectrum of the global mean of XCO2 using the Morlet wavelet (ω0 = 6) as a function of time and Fourier equivalent wave period (in years). The black solid contour lines enclose regions of more than 95% confidence relative to a red-noise random process (α = 0.77). The cross-hatched and shaded area indicates the influence of the cone of influence of the mother wavelet. (c) Global wavelet power spectrum (solid black line) and Fourier power spectrum (solid gray line). The dotted line indicates the 95% confidence level. (d) The global mean of XCO2 in time-series harmonic (red line) and linear (blue line) trends.
Figure 17. Wavelet analysis of the XCO2 based on GOSAT data. (a) Time series (solid black line) and inverse wavelet transform (solid gray line). (b) Normalized wavelet power spectrum of the global mean of XCO2 using the Morlet wavelet (ω0 = 6) as a function of time and Fourier equivalent wave period (in years). The black solid contour lines enclose regions of more than 95% confidence relative to a red-noise random process (α = 0.77). The cross-hatched and shaded area indicates the influence of the cone of influence of the mother wavelet. (c) Global wavelet power spectrum (solid black line) and Fourier power spectrum (solid gray line). The dotted line indicates the 95% confidence level. (d) The global mean of XCO2 in time-series harmonic (red line) and linear (blue line) trends.
Remotesensing 14 02622 g017
Figure 18. Wavelet analysis of the XCH4 based on GOSAT data. (a) Time series (solid black line) and inverse wavelet transform (solid gray line). (b) Normalized wavelet power spectrum of the global mean of XCH4 using the Morlet wavelet (ω0 = 6) as a function of time and Fourier equivalent wave period (in years). The black solid contour lines enclose regions of more than 95% confidence relative to a red-noise random process (α = 0.77). The cross-hatched and shaded area indicates the influence of the cone of influence of the mother wavelet. (c) Global wavelet power spectrum (solid black line) and Fourier power spectrum (solid gray line). The dotted line indicates the 95% confidence level, the wavelet power spectrum of global average seasonal atmospheric methane. (d) Global mean of XCH4 in time-series harmonic (red line) and linear (blue line) trends.
Figure 18. Wavelet analysis of the XCH4 based on GOSAT data. (a) Time series (solid black line) and inverse wavelet transform (solid gray line). (b) Normalized wavelet power spectrum of the global mean of XCH4 using the Morlet wavelet (ω0 = 6) as a function of time and Fourier equivalent wave period (in years). The black solid contour lines enclose regions of more than 95% confidence relative to a red-noise random process (α = 0.77). The cross-hatched and shaded area indicates the influence of the cone of influence of the mother wavelet. (c) Global wavelet power spectrum (solid black line) and Fourier power spectrum (solid gray line). The dotted line indicates the 95% confidence level, the wavelet power spectrum of global average seasonal atmospheric methane. (d) Global mean of XCH4 in time-series harmonic (red line) and linear (blue line) trends.
Remotesensing 14 02622 g018
Figure 18a shows the seasonal variation extracted from the global mean XCH4 GOSAT data measured from January 2010 to May 2021. During the period analyzed, the global mean XCH4 rose continuously, as shown by the harmonic and linear trends in Figure 18d. In Figure 18b, the frequency, increment every six months, and decreased XCH4 values are visible.

3.4. Greenhouse Gas Growth Analysis by Continent

The temporal variation by continent is visibly divided by columns and rows that represent the years (January 2010 to May 2021) and the seasons (spring, summer, autumn, and winter), based on those of the Northern Hemisphere, and the regions analyzed are South America, Africa, Australia, Europe, Asia, and North America. This local analysis highlights the areas with a high concentration of CO2 and methane.
South America presented a range of CO2 of 30 ppm between January 2010 and May 2021. It showed low levels of CO2 in spring, and its maximum values occurred in summer; the lowest value recorded was 385 ppm, registered in 2010. Over time, this gas increased, reaching values of approximately 415 ppm. The most affected country was Brazil, along with its surrounding countries, such as Bolivia, Uruguay, and Paraguay, as well as the north of Argentina. The epicenter of the heat CO2 concentration was at (10°W, 55°S).
Africa had a variation of 40 ppm over the last decade. The highest concentrations of CO2 occurred in winter, and the lowest in autumn. The lowest value recorded was 380 ppm, registered in 2010. Over time, this gas increased to 420 ppm. An area with remarkably high concentrations of CO2 was (20°E, 20°S), near Congo, Gabon, and Cameron. The second most affected area was located near (0°E, 20°N); the nearest countries to this area are Mali, Algeria, Mauritania, and Morocco, as well as the Western Sahara.
Australia presented a range of CO2 of 35 ppm between January 2010 and May 2021, and showed low levels of CO2 in spring. Its maximum values occurred in winter. The lowest value recorded was 380 ppm, registered in 2010. Over the time, this gas increased, reaching values of 415 ppm. The most affected country was Australia, especially its central region, at coordinates of (140°E, 25°S).
Regarding Europe, typically, the highest concentrations of CO2 were the central region. The highest value of CO2 in parts per million recorded in the last 10 years was approximately 420 ppm. The lowest concentrations occurred in autumn, and the maximum concentrations were observed in winter.
Five representative locations in Asia displayed a high CO2 rate: India, Indonesia, Japan, Beijing and the south-western area of China. The coordinates of the representative increments of CO2 are (75°E, 15°N), (45°E, 80°N), (100°E, 25°N), (125°E, 18°N), and (130°E, 45°N). The north-west area of Russia (near Moscow) also presented a high increment of CO2. This region had a variation of 45 ppm over the last decade, and it is one of the most polluted regions on Earth. The highest concentrations of CO2 occurred in winter, and the lowest in autumn; the lowest value recorded was 380 ppm, registered in 2010. Over the time, this gas increased, reaching values of 425 ppm.
North America presented a range of CO2 of 45 ppm between January 2010 and May 2021. The region showed low levels of CO2 in the autumn, and its maximum values occurred in winter. The lowest value recorded was 380 ppm, registered in 2010. Over time, this gas increased, reaching values of approximately 425 ppm. The most affected countries were United States (the central region of the United States was less affected than the east and west) and Canada.
Previous studies were examined to compare the differences between the source data (e.g., land use) of the carbon inventories and the real atmospheric carbon observations [52,53]. The authors of [54] reported that the Orbiting Carbon Observatory-2 (OCO-2) satellite-based XCO2 shows a positive correlation with the Open-Source Data Inventory for Anthropogenic CO2 (ODIAC) for the main pollution regions: the eastern USA, central Europe, and East Asia. There is also literature that features comparisons between ground-based atmospheric observations (eddy flux tower measurements) and the publicly available inventory for urban atmospheric CO2 on the city scale [55,56].
Figure 19 shows the spatial mean tendency of the total columns of CO2 and CH4 in the world, as well as for Asia, Australia, Africa, South America, North America, and Europe. Regarding the XCO2 tendency, Australia shows low variations and an almost linear tendency; Europe, Asia, and North America show variations of approximately 10 ppm per year. However, XCH4 shows low concentrations in South America and, especially, in Australia; higher concentrations of CH4 can be observed in North America, Asia, Europe, and Africa.

4. Conclusions

The main objective of this study was to reconstruct the GOSAT L3 product based on the EOF algorithm and validate it through TCCON measurements. The global spatial–temporal pattern of CO2 and CH4 is illustrated better with the application of this algorithm, as can be seen through the results presented in this paper. The southern part of the Earth, about which there is a low amount of information for spring–summer, is illustrated better, and the exact point in the Southern Hemisphere where the levels of these gases have reduced considerably can be observed. Furthermore, the northern part of the Earth, about which there is little information from the GOSAT L3 product for autumn–winter, is represented better, especially in the ocean areas. However, in the Antarctic area, the information could not be reconstructed, since there were no values recorded in the last 10 years.
The monthly time series shows similitude between the in-situ stations TCCON and GOSAT and CO2 in terms of sufficiency. This assertion is difficult to corroborate in the Southern Hemisphere due to the lack of TCCON stations around the world. Moreover, the TCCON manifested higher CO2 concentrations compared with the CO2 recovered from the GOSAT. The kernel density plots presented in the section on the data validation seem to show higher similitude regarding CO2. The CH4 gas data are more unreliable, especially between 2010 and 2012, due to the lack of information from within those years.
The GOSAT CH4 dataset currently contains over a decade of worldwide CH4 behavior. The ability to assess worldwide and territorial CH4 outflows is imperative to develop our understanding of the worldwide level of methane. We started this work by highlighting the wide assortment of aspects of the past forms of this dataset that have contributed towards this understanding, illustrating the utility of this dataset for analyzing and understanding the worldwide behavior of CH4, which is still rising.
Summarizing the results, it can be concluded that the area most affected by greenhouse gases is the Northern Hemisphere. Moreover, both gases have a harmonic variation over time. In the case of CO2, at the beginning and the end of each year, the gas starts to increase rapidly, and in the middle of the year, its levels reduce slightly (less than the increment). On the other hand, the seasonal variation in CH4 is different because the gas starts to increase in the middle of each year and decreases at the end and the beginning of each year. Both have a similar frequency of events, repeating the same event at approximately the same rate over a given year, which suggests that the seasonal behavior of both gases resembles a sinusoidal wave.
In this study, some shortcomings are evident in the data interpolation, especially for the information obtained from the South America, Africa, and Antarctica regions. The results for these regions were unreliable due to the low quantity or absence of TCCON data in these regions; therefore, GOSAT and TCCON data validation analyses could not be performed. To increase the reliability of the results, in a future study, it will be necessary to compare the interpolated results with simulated results using the atmospheric transport model available in the GOSAT L4 products, or compare them with other models or other satellite data. Furthermore, the empirical nature of the data and their relatively high level of uncertainty indicate the need for future satellite sensors to increase data availability and accuracy.

Author Contributions

F.P.A.L.: investigation, data curation, writing—original draft. G.Z.: conceptualization, methodology, writing—review and editing. G.J.: formal analysis, resources. K.Z.: writing—review and editing. Y.T.: validation, formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant number 2021YFB3901000, 2021YFB3901005 and 2019YFE0126400) and the National Natural Science Foundation of China (grant number 41971320).

Data Availability Statement

The dataset can be downloaded from GitHub and National Earth System Science Data Center, National Science and Technology Infrastructure of China (http://www.geodata.cn, accessed on 12 December 2021).

Acknowledgments

The authors acknowledge the GOSAT Data Archive Service (GDAS) (https://data2.gosat.nies.go.jp/ (accessed on 12 December 2021)) for providing the CO2 and CH4 datasets and the California Institute of Technology (Caltech) for providing the TCCON data (http://tccon.ornl.gov/ (accessed on 12 December 2021)).

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Chevallier, F.; Bréon, F.-M.; Rayner, P.J. Contribution of the Orbiting Carbon Observatory to the Estimation of CO2 Sources and Sinks: Theoretical Study in a Variational Data Assimilation Framework. J. Geophys. Res. 2007, 112, 3–10. [Google Scholar] [CrossRef]
  2. Wunch, D.; Toon, G.C.; Blavier, J.-F.L.; Washenfelder, R.A.; Notholt, J.; Connor, B.J.; Griffith, D.W.T.; Sherlock, V.; Wennberg, P.O. The Total Carbon Column Observing Network. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2011, 369, 2087–2112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Dunn, R.J.H.; Aldred, F.; Gobron, N.; Miller, J.B.; Willett, K.M.; Ades, M.; Adler, R.; Allan, R.P.; Allan, R.; Anderson, J.; et al. Global Climate. Bull. Am. Meteorol. Soc. 2021, 102, S11–S142. [Google Scholar] [CrossRef]
  4. Hall, B.D.; Crotwell, A.M.; Kitzis, D.R.; Mefford, T.; Miller, B.R.; Schibig, M.F.; Tans, P.P. Revision of the World Meteorological Organization Global Atmosphere Watch (WMO/GAW) CO2; Calibration Scale. Atmos. Meas. Tech. 2021, 14, 3015–3032. [Google Scholar] [CrossRef]
  5. Liang, A.; Gong, W.; Han, G.; Xiang, C. Comparison of Satellite-Observed XCO2 from GOSAT, OCO-2, and Ground-Based TCCON. Remote Sens. 2017, 9, 1033. [Google Scholar] [CrossRef] [Green Version]
  6. Etminan, M.; Myhre, G.; Highwood, E.J.; Shine, K.P. Radiative Forcing of Carbon Dioxide, Methane, and Nitrous Oxide: A Significant Revision of the Methane Radiative Forcing. Geophys. Res. Lett. 2016, 43, 12614–12623. [Google Scholar] [CrossRef]
  7. Antonov, K.; Poddubny, V.; Markelov, Y.; Buevich, A.; Medvedev, A. Dynamics of Surface Carbon Dioxide and Methane Concentrations on the Arctic Belyy Island in 2015–2017 Summertime. In Proceedings of the 24th International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, Tomsk, Russia, 2–5 July 2018; Romanovskii, O.A., Matvienko, G.G., Eds.; SPIE: Bellingham, WA, USA, 2018; p. 360. [Google Scholar]
  8. Dlugokencky, E.J.; Bruhwiler, L.; White, J.W.C.; Emmons, L.K.; Novelli, P.C.; Montzka, S.A.; Masarie, K.A.; Lang, P.M.; Crotwell, A.M.; Miller, J.B.; et al. Observational Constraints on Recent Increases in the Atmospheric CH4 Burden. Geophys. Res. Lett. 2009, 36, 2–4. [Google Scholar] [CrossRef] [Green Version]
  9. Rigby, M.; Prinn, R.G.; Fraser, P.J.; Simmonds, P.G.; Langenfelds, R.L.; Huang, J.; Cunnold, D.M.; Steele, L.P.; Krummel, P.B.; Weiss, R.F.; et al. Renewed Growth of Atmospheric Methane. Geophys. Res. Lett. 2008, 35, 2–5. [Google Scholar] [CrossRef] [Green Version]
  10. Kirschke, S.; Bousquet, P.; Ciais, P.; Saunois, M.; Canadell, J.G.; Dlugokencky, E.J.; Bergamaschi, P.; Bergmann, D.; Blake, D.R.; Bruhwiler, L.; et al. Three Decades of Global Methane Sources and Sinks. Nat. Geosci. 2013, 6, 813–823. [Google Scholar] [CrossRef]
  11. Saunois, M.; Stavert, A.R.; Poulter, B.; Bousquet, P.; Canadell, J.G.; Jackson, R.B.; Raymond, P.A.; Dlugokencky, E.J.; Houweling, S.; Patra, P.K.; et al. The Global Methane Budget 2000–2017. Earth Syst. Sci. Data 2020, 12, 1561–1623. [Google Scholar] [CrossRef]
  12. Hausmann, P.; Sussmann, R.; Smale, D. Contribution of Oil and Natural Gas Production to Renewed Increase in Atmospheric Methane (2007–2014): Top–down Estimate from Ethane and Methane Column Observations. Atmos. Chem. Phys. 2016, 16, 3227–3244. [Google Scholar] [CrossRef] [Green Version]
  13. McNorton, J.; Gloor, E.; Wilson, C.; Hayman, G.D.; Gedney, N.; Comyn-Platt, E.; Marthews, T.; Parker, R.J.; Boesch, H.; Chipperfield, M.P. Role of Regional Wetland Emissions in Atmospheric Methane Variability. Geophys. Res. Lett. 2016, 43, 11433–11444. [Google Scholar] [CrossRef] [Green Version]
  14. Nisbet, E.G.; Fisher, R.E.; Lowry, D.; France, J.L.; Allen, G.; Bakkaloglu, S.; Broderick, T.J.; Cain, M.; Coleman, M.; Fernandez, J.; et al. Methane Mitigation: Methods to Reduce Emissions, on the Path to the Paris Agreement. Rev. Geophys. 2020, 58, e2019RG000675. [Google Scholar] [CrossRef]
  15. Schaefer, H.; Fletcher, S.E.M.; Veidt, C.; Lassey, K.R.; Brailsford, G.W.; Bromley, T.M.; Dlugokencky, E.J.; Michel, S.E.; Miller, J.B.; Levin, I.; et al. A 21st-Century Shift from Fossil-Fuel to Biogenic Methane Emissions Indicated by 13CH4. Science 2016, 352, 80–84. [Google Scholar] [CrossRef]
  16. Jacob, D.J.; Turner, A.J.; Maasakkers, J.D.; Sheng, J.; Sun, K.; Liu, X.; Chance, K.; Aben, I.; McKeever, J.; Frankenberg, C. Satellite Observations of Atmospheric Methane and Their Value for Quantifying Methane Emissions. Atmos. Chem. Phys. 2016, 16, 14371–14396. [Google Scholar] [CrossRef] [Green Version]
  17. le Quéré, C.; Jackson, R.B.; Jones, M.W.; Smith, A.J.P.; Abernethy, S.; Andrew, R.M.; De-Gol, A.J.; Willis, D.R.; Shan, Y.; Canadell, J.G.; et al. Temporary Reduction in Daily Global CO2 Emissions during the COVID-19 Forced Confinement. Nat. Clim. Chang. 2020, 10, 647–653. [Google Scholar] [CrossRef]
  18. Friedlingstein, P.; O’Sullivan, M.; Jones, M.W.; Andrew, R.M.; Hauck, J.; Olsen, A.; Peters, G.P.; Peters, W.; Pongratz, J.; Sitch, S.; et al. Global Carbon Budget 2020. Earth Syst. Sci. Data 2020, 12, 3269–3340. [Google Scholar] [CrossRef]
  19. Toon, G.; Blavier, J.-F.; Washenfelder, R.; Wunch, D.; Keppel-Aleks, G.; Wennberg, P.; Connor, B.; Sherlock, V.; Griffith, D.; Deutscher, N.; et al. Total Column Carbon Observing Network (TCCON). In Fourier Transform Spectroscopy; OSA: Washington, DC, USA, 2009. [Google Scholar]
  20. Manning, A.C.; Nisbet, E.G.; Keeling, R.F.; Liss, P.S. Greenhouse Gases in the Earth System: Setting the Agenda to 2030. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2011, 369, 1885–1890. [Google Scholar] [CrossRef]
  21. Koruba, D.; Piotrowski, J.Z.; Latosińska, J. Biomass—Alternative Renewable Energy Source to the Fossil Fuels. E3S Web Conf. 2017, 14, 2015. [Google Scholar] [CrossRef]
  22. Baker, D.F.; Bösch, H.; Doney, S.C.; O’Brien, D.; Schimel, D.S. Carbon Source/Sink Information Provided by Column CO2 Measurements from the Orbiting Carbon Observatory. Atmos. Chem. Phys. 2010, 10, 4145–4165. [Google Scholar] [CrossRef] [Green Version]
  23. Hungershoefer, K.; Breon, F.-M.; Peylin, P.; Chevallier, F.; Rayner, P.; Klonecki, A.; Houweling, S.; Marshall, J. Evaluation of Various Observing Systems for the Global Monitoring of CO2 Surface Fluxes. Atmos. Chem. Phys. 2010, 10, 10503–10520. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, H.; Jiang, F.; Wang, J.; Ju, W.; Chen, J.M. Terrestrial Ecosystem Carbon Flux Estimated Using GOSAT and OCO-2 XCO2 Retrievals. Atmos. Chem. Phys. 2019, 19, 12067–12082. [Google Scholar] [CrossRef] [Green Version]
  25. Bruhwiler, L.; Basu, S.; Butler, J.H.; Chatterjee, A.; Dlugokencky, E.; Kenney, M.A.; McComiskey, A.; Montzka, S.A.; Stanitski, D. Observations of Greenhouse Gases as Climate Indicators. Clim. Chang. 2021, 165, 12. [Google Scholar] [CrossRef] [PubMed]
  26. Kuze, A.; Suto, H.; Shiomi, K.; Kawakami, S.; Tanaka, M.; Ueda, Y.; Deguchi, A.; Yoshida, J.; Yamamoto, Y.; Kataoka, F.; et al. Update on GOSAT TANSO-FTS Performance, Operations, and Data Products after More than 6 Years in Space. Atmos. Meas. Tech. 2016, 9, 2445–2461. [Google Scholar] [CrossRef] [Green Version]
  27. Yokota, T. GOSAT Mission Operation Extended after Nominal Lifetime. NIES GOSAT Project Letter, Issue No 30, March 2014; p. 10. Available online: http://www.gosat.nies.go.jp/eng/newsletter/newsletter30e.pdf (accessed on 17 November 2021).
  28. Suto, H.; Shiomi, K.; Nakajima, M. Space-Based Carbon Monitoring by GOSAT and GOSAT-2. In Proceedings of the IWGGMS-10 (10th International Workshop on Greenhouse Gas Measurements from Space) ESA/ESTEC, Noordwijk, The Netherlands, 5–7 May 2014. [Google Scholar]
  29. Yokota, T.; Yoshida, Y.; Eguchi, N.; Ota, Y.; Tanaka, T.; Watanabe, H.; Maksyutov, S. Global Concentrations of CO2 and CH4 Retrieved from GOSAT: First Preliminary Results. SOLA 2009, 5, 160–163. [Google Scholar] [CrossRef] [Green Version]
  30. Bril, A.; Maksyutov, S. Development of the Empirical Orthogonal Functions-Based Algorithm for the Retrievals of Atmospheric CO2 Total Column Amount from Space-Borne Observations of Reflected Sunlight. J. Appl. Remote Sens. 2018, 12, 046012. [Google Scholar] [CrossRef]
  31. Nakajima, T.; Nakajima, T.Y.; Higurashi, A.; Sano, I.; Takamura, T.; Ishida, H.; Schutgens, N. A Study of Aerosol and Cloud Information Retrievals from CAI Imager on Board GOSAT Satellite. J. Remote Sens. Soc. Jpn. 2008, 28, 178–189. [Google Scholar]
  32. Messerschmidt, J.; Geibel, M.C.; Blumenstock, T.; Chen, H.; Deutscher, N.M.; Engel, A.; Feist, D.G.; Gerbig, C.; Gisi, M.; Hase, F.; et al. Calibration of TCCON Column-Averaged CO2: The First Aircraft Campaign over European TCCON Sites. Atmos. Chem. Phys. 2011, 11, 10765–10777. [Google Scholar] [CrossRef] [Green Version]
  33. Toon, G.C.; Wunch, D. 2015: A stand-alone a priori profile generation tool for GGG2014 release (Version GGG2014. R0), CaltechDATA. Available online: https://data.caltech.edu/records/250 (accessed on 17 November 2021). [CrossRef]
  34. Wunch, D.; Toon, G.C.; Sherlock, V.; Deutscher, N.M.; Liu, C.; Feist, D.G.; Wennberg, P.O. Documentation for the 2014 TCCON Data Release. CaltechDATA; Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory: Oak Ridge, TN, USA, 2015. [CrossRef]
  35. Total Carbon Column Observing Network (TCCON) Team. 2014 TCCON Data Release (Version GGG2014) [Data Set]. CaltechDATA. 2017. Available online: https://data.caltech.edu/records/293 (accessed on 16 August 2019). [CrossRef]
  36. Keppel-Aleks, G.; Wennberg, P.O.; Schneider, T.; Honsowetz, N.Q.; Vay, S.A. Total Column Constraints on Northern Hemisphere Carbon Dioxide Surface Exchange. In Proceedings of the AGU Fall Meeting Abstracts, San Francisco, CA, USA, 15–19 December 2008; Volume 2008. [Google Scholar]
  37. Wunch, D.; Wennberg, P.O.; Toon, G.C.; Connor, B.J.; Fisher, B.; Osterman, G.B.; Frankenberg, C.; Mandrake, L.; O’Dell, C.; Ahonen, P.; et al. A Method for Evaluating Bias in Global Measurements of CO2 Total Columns from Space. Atmos. Chem. Phys. 2011, 11, 12317–12337. [Google Scholar] [CrossRef] [Green Version]
  38. Alvera-Azcárate, A.; Barth, A.; Sirjacobs, D.; Lenartz, F.; Beckers, J.M. Data Interpolating Empirical Orthogonal Functions (DINEOF): A Tool for Geophysical Data Analyses. Mediterr. Mar. Sci. 2011, 12, 5–11. [Google Scholar] [CrossRef] [Green Version]
  39. Watanabe, H.; Hayashi, K.; Saeki, T.; Maksyutov, S.; Nasuno, I.; Shimono, Y.; Hirose, Y.; Takaichi, K.; Kanekon, S.; Ajiro, M.; et al. Global Mapping of Greenhouse Gases Retrieved from GOSAT Level 2 Products by Using a Kriging Method. Int. J. Remote Sens. 2015, 36, 1509–1528. [Google Scholar] [CrossRef] [Green Version]
  40. Taylor, M.H.; Losch, M.; Wenzel, M.; Schröter, J. On the Sensitivity of Field Reconstruction and Prediction Using Empirical Orthogonal Functions Derived from Gappy Data. J. Clim. 2013, 26, 9194–9205. [Google Scholar] [CrossRef] [Green Version]
  41. Alvera-Azcárate, A.; Barth, A.; Rixen, M.; Beckers, J.M. Reconstruction of Incomplete Oceanographic Data Sets Using Empirical Orthogonal Functions: Application to the Adriatic Sea Surface Temperature. Ocean Model. 2005, 9, 325–346. [Google Scholar] [CrossRef] [Green Version]
  42. Beckers, J.M.; Rixen, M. EOF Calculations and Data Filling from Incomplete Oceanographic Datasets*. J. Atmos. Ocean. Technol. 2003, 20, 1839–1856. [Google Scholar] [CrossRef]
  43. Parker, R.J.; Boesch, H.; Wooster, M.J.; Moore, D.P.; Webb, A.J.; Gaveau, D.; Murdiyarso, D. Atmospheric CH4 and CO2 Enhancements and Biomass Burning Emission Ratios Derived from Satellite Observations of the 2015 Indonesian Fire Plumes. Atmos. Chem. Phys. 2016, 16, 10111–10131. [Google Scholar] [CrossRef] [Green Version]
  44. Uno, T.; Oishi, T.; Fujii, Y.; Imaki, M.; Sato, A.; Kataoka, F.; Imai, H.; Hashimoto, M.; Shiomi, K.; Nakajima, M. The Development and On-Orbit Calibration Status of GOSAT-2 TANSO-CAI-2 Instrument. In Proceedings of the International Conference on Space Optics—ICSO 2020, Virtual, 30 March–2 April 2021; Sodnik, Z., Cugny, B., Karafolas, N., Eds.; SPIE: Bellingham, WA, USA, 2021. [Google Scholar]
  45. Shiomi, K.; Suto, H.; Kawakami, S.; Kina, T.; Yoshida, M.; Mitomi, Y.; Sekio, N.; Kataoka, F. Calibration Plan of GOSAT TANSO. In Proceedings of the Fourier Transform Spectroscopy/Hyperspectral Imaging and Sounding of the Environment, Santa Fe, NM, USA, 11–15 February 2007; OSA: Washington, DC, USA, 2007. [Google Scholar]
  46. National Institute for Environmental Studies Ministry of the Environment, Japan Aerospace Exploration Agency. “IBUKI” (GOSAT) Satellite’s Historically Unprecedented Observation of Anthropogenic CO2 Emissions of Mega-City Regions from Space. 2014. Available online: https://www.gosat.nies.go.jp/eng/result/download/GOSAT_L4CO2_20141209_en.pdf (accessed on 23 December 2021).
  47. AK-BHD, M. WMO Greenhouse Gas Bulletin; World Meteorological Organization: Geneva, Switzerland, 2021. [Google Scholar]
  48. Giani, L.; Bredenkamp, J.; Eden, I. Temporal and Spatial Variability of the CH4 Dynamics of Landfill Cover Soils. J. Plant Nutr. Soil Sci. 2002, 165, 205–210. [Google Scholar] [CrossRef]
  49. Hocke, K.; Kämpfer, N. Hovmöller Diagrams of Climate Anomalies in NCEP/NCAR Reanalysis from 1948 to 2009. Clim. Dyn. 2011, 36, 355–364. [Google Scholar] [CrossRef] [Green Version]
  50. Shumway, R.H.; Stoffer, D.S. Time Series Analysis and Its Applications; Springer: New York, NY, USA, 2011; ISBN 978-1-4419-7864-6. [Google Scholar]
  51. The Ministry of the Environment, Japan (MOEJ). Whole-Atmospheric Monthly CO2 Concentration Tops 400 Ppm-Preliminary GOSAT Monitoring Results. 2016. Available online: http://www.gosat.nies.go.jp/newpdf/GOSATpressrelease_20160520_400ppm_en.pdf (accessed on 17 November 2021).
  52. Janardanan, R.; Maksyutov, S.; Oda, T.; Saito, M.; Kaiser, J.W.; Ganshin, A.; Stohl, A.; Matsunaga, T.; Yoshida, Y.; Yokota, T. Comparing GOSAT Observations of Localized CO 2 Enhancements by Large Emitters with Inventory-based Estimates. Geophys. Res. Lett. 2016, 43, 3486–3493. [Google Scholar] [CrossRef] [Green Version]
  53. Hwang, Y.; Um, J.-S. Exploring Causal Relationship between Landforms and Ground Level CO2 in Dalseong Forestry Carbon Project Site of South Korea. Spat. Inf. Res. 2017, 25, 361–370. [Google Scholar] [CrossRef]
  54. Hakkarainen, J.; Ialongo, I.; Tamminen, J. Direct Space-based Observations of Anthropogenic CO2 Emission Areas from OCO-2. Geophys. Res. Lett. 2016, 43, 11400–11406. [Google Scholar] [CrossRef]
  55. Lin, J.C.; Mitchell, L.; Crosman, E.; Mendoza, D.L.; Buchert, M.; Bares, R.; Fasoli, B.; Bowling, D.R.; Pataki, D.; Catharine, D.; et al. CO2 and Carbon Emissions from Cities: Linkages to Air Quality, Socioeconomic Activity, and Stakeholders in the Salt Lake City Urban Area. Bull. Am. Meteorol. Soc. 2018, 99, 2325–2339. [Google Scholar] [CrossRef]
  56. Levin, I.; Hammer, S.; Eichelmann, E.; Vogel, F.R. Verification of Greenhouse Gas Emission Reductions: The Prospect of Atmospheric Monitoring in Polluted Areas. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2011, 369, 1906–1924. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of TCCON stations and quantity of information collected in GGG2014 data.
Figure 1. Location of TCCON stations and quantity of information collected in GGG2014 data.
Remotesensing 14 02622 g001
Figure 2. XCO2 global interpolation procedure corresponding to the GOSAT data of January 2021. (a) FTS SWIR L2 XCO2, (b) FTS SWIR L3 XCO2—interpolating with kriging method, (c) XCO2 temporal interpolation using EOF method.
Figure 2. XCO2 global interpolation procedure corresponding to the GOSAT data of January 2021. (a) FTS SWIR L2 XCO2, (b) FTS SWIR L3 XCO2—interpolating with kriging method, (c) XCO2 temporal interpolation using EOF method.
Remotesensing 14 02622 g002
Figure 3. Record of the iterations performed to fill in the gaps with an increasing number of EOFs (#). (a) XCO2 Temporal matrix, (b) XCH4 temporal matrix.
Figure 3. Record of the iterations performed to fill in the gaps with an increasing number of EOFs (#). (a) XCO2 Temporal matrix, (b) XCH4 temporal matrix.
Remotesensing 14 02622 g003
Figure 4. Comparison of EOF coefficients of the GOSAT L3 applying EOF algorithm and GOSAT L4B column average. Black points represent “true” field and red line represents field gappy data. (a) XCO2 temporal matrix, (b) XCH4 temporal matrix.
Figure 4. Comparison of EOF coefficients of the GOSAT L3 applying EOF algorithm and GOSAT L4B column average. Black points represent “true” field and red line represents field gappy data. (a) XCO2 temporal matrix, (b) XCH4 temporal matrix.
Remotesensing 14 02622 g004aRemotesensing 14 02622 g004b
Figure 5. Quantity of reconstructed pixels applying the EOF algorithm in GOSAT L3 products in a temporal variation from January 2010–May 2021.
Figure 5. Quantity of reconstructed pixels applying the EOF algorithm in GOSAT L3 products in a temporal variation from January 2010–May 2021.
Remotesensing 14 02622 g005
Figure 6. 3-D kernel scatter plot comparing XCO2 GOSAT L3 applying the EOF algorithm, L4B column average, and monthly mean TCCON.
Figure 6. 3-D kernel scatter plot comparing XCO2 GOSAT L3 applying the EOF algorithm, L4B column average, and monthly mean TCCON.
Remotesensing 14 02622 g006
Figure 7. 3-D kernel scatter plot comparing XCH4 GOSAT L3 applying the EOF algorithm, L4B column average, and monthly mean TCCON.
Figure 7. 3-D kernel scatter plot comparing XCH4 GOSAT L3 applying the EOF algorithm, L4B column average, and monthly mean TCCON.
Remotesensing 14 02622 g007
Figure 8. XCO2 TCCON and GOSAT data comparison.
Figure 8. XCO2 TCCON and GOSAT data comparison.
Remotesensing 14 02622 g008
Figure 9. XCH4 TCCON and GOSAT data comparison.
Figure 9. XCH4 TCCON and GOSAT data comparison.
Remotesensing 14 02622 g009
Figure 10. Difference between monthly average TCCON in situ data (January 2010–May 2021) with (a) ΔXCO2, (b) ΔXCH4.
Figure 10. Difference between monthly average TCCON in situ data (January 2010–May 2021) with (a) ΔXCO2, (b) ΔXCH4.
Remotesensing 14 02622 g010
Figure 11. Global mean of XCO2 monthly growth analysis from January 2010 to May 2021.
Figure 11. Global mean of XCO2 monthly growth analysis from January 2010 to May 2021.
Remotesensing 14 02622 g011
Figure 12. Global mean of XCH4 monthly growth analysis from January 2010 to May 2021.
Figure 12. Global mean of XCH4 monthly growth analysis from January 2010 to May 2021.
Remotesensing 14 02622 g012
Figure 13. Spatial and temporal behavior of monthly mean XCO2 L3, with EOF algorithm, from around the world.
Figure 13. Spatial and temporal behavior of monthly mean XCO2 L3, with EOF algorithm, from around the world.
Remotesensing 14 02622 g013
Figure 14. Spatial and temporal behavior of monthly mean XCH4 L3, with EOF algorithm, around the world.
Figure 14. Spatial and temporal behavior of monthly mean XCH4 L3, with EOF algorithm, around the world.
Remotesensing 14 02622 g014
Figure 15. Hovmöller plots. (a) Latitudinal XCO2 mean vs. time (GOSAT L3 with EOF algorithm), (b) latitudinal XCO2 mean vs. time (GOSAT L4B), (c) latitudinal RMSE mean vs. time (GOSAT L4B–GOSAT L3 with EOF algorithm).
Figure 15. Hovmöller plots. (a) Latitudinal XCO2 mean vs. time (GOSAT L3 with EOF algorithm), (b) latitudinal XCO2 mean vs. time (GOSAT L4B), (c) latitudinal RMSE mean vs. time (GOSAT L4B–GOSAT L3 with EOF algorithm).
Remotesensing 14 02622 g015
Figure 16. Hovmöller plots. (a) Latitudinal XCH4 mean vs. time (GOSAT L3 with EOF algorithm), (b) latitudinal XCH4 mean vs. time (GOSAT L4B), (c) latitudinal RMSE mean vs. time (GOSAT L4B–GOSAT L3 with EOF algorithm).
Figure 16. Hovmöller plots. (a) Latitudinal XCH4 mean vs. time (GOSAT L3 with EOF algorithm), (b) latitudinal XCH4 mean vs. time (GOSAT L4B), (c) latitudinal RMSE mean vs. time (GOSAT L4B–GOSAT L3 with EOF algorithm).
Remotesensing 14 02622 g016
Figure 19. (a) Time series of global and continental mean XCO2. (b) Time series of global and continental mean XCH4.
Figure 19. (a) Time series of global and continental mean XCO2. (b) Time series of global and continental mean XCH4.
Remotesensing 14 02622 g019
Table 1. Number of empty pixels in a stack of data from GOSAT L3 products (January 2010–May 2021).
Table 1. Number of empty pixels in a stack of data from GOSAT L3 products (January 2010–May 2021).
ContinentTime RangeNumber of Empty
Pixels
Number of Layers
AfricaJanuary 2010–May 202193,916137
AsiaJanuary 2010–May 2021229,291137
AustraliaJanuary 2010–May 202131,871137
EuropeJanuary 2010–May 202184,418137
North AmericaJanuary 2010–May 2021226,358137
South AmericaJanuary 2010–May 202163,236137
Table 2. TCCON stations coordinate locations, linear b slope coefficient, and R-squared coefficient based on EOF–GOSAT L3 data and mean values per month of TCCON data.
Table 2. TCCON stations coordinate locations, linear b slope coefficient, and R-squared coefficient based on EOF–GOSAT L3 data and mean values per month of TCCON data.
StationDateLat.Long.XCO2 r2XCH4 r2XCO2 bXCH4 b
1Ascension Island2012-05-22 to 2018-10-31−14.33−7.920.920.530.950.62
2Anmeyondo2015-02-02 to 2018-04-18126.3336.540.970.901.061.09
3Bialystok2009-03-01 to 2018-10-0123.0253.230.840.820.941.15
4Bremen2010-01-22 to 2019-08-238.8553.100.920.470.880.55
5Burgos2017-03-03 to 2020-02-28120.6518.530.900.520.870.62
6Caltech2012-09-20 to 2020-10-31−118.1334.140.980.931.081.05
7Darwin2005-08-28 to 2020-02-26130.89−12.430.980.931.031.03
8Edwards2013-07-20 to 2020-10-30−117.8834.960.840.810.971.12
9East trout lake2016-10-07 to 2020-09-07−104.9954.360.940.881.050.92
10Eureka2010-07-24 to 2020-07-07−86.4280.050.900.500.960.57
11Four corners2013-03-16 to 2013-10-04−108.4836.800.920.760.900.96
12Garmisch2007-07-16 to 2019-10-1811.0647.480.970.901.011.03
13Hefei2015-09-18 to 2016-12-31117.1731.900.980.961.050.99
14Indianapolis2012-08-23 to 2012-12-01−86.0039.860.970.770.990.74
15Izana2007-05-18 to 2020-11-25−16.4828.30--0.620.66
16Jet propulsion lab2007-07-31 to 2018-05-14−118.1834.200.970.931.010.95
17Saga2011-07-28 to 2020-09-01130.2933.240.870.470.791.28
18Karlsruhe2010-04-19 to 2020-11-308.4449.100.980.961.001.05
19Lauder2004-06-28 to 2020-08-31169.68−45.050.970.941.091.02
20Nicosia2019-08-31 to 2020-02-2833.3835.140.980.941.021.06
21Lamont2008-07-06 to 2020-10-30−97.4936.600.870.671.200.74
22Orleans2009-08-29 to 2019-09-182.1147.970.930.881.090.94
23Park falls2004-06-02 to 2020-10-31−90.2745.940.980.971.041.07
24Paris2014-09-23 to 2020-02-262.3648.850.960.941.051.04
25Reunion island2011-09-16 to 2020-07-1855.49−20.900.970.961.081.04
26Rikubetsu2013-11-16 to 2019-09-30143.7743.460.720.770.751.04
27Sodankyla2009-05-16 to 2020-10-2026.6367.370.980.960.991.08
28Ålesund2014-04-06 to 2019-09-1511.9278.920.500.240.320.26
29Tsukuba2011-08-04 to 2019-09-30140.1236.050.960.831.020.96
30Wollongong2008-06-26 to 2020-02-28150.88−34.410.910.861.071.04
31Zugspitze2015-04-24 to 2019-10-1710.9847.420.990.951.021.10
Table 3. CO2 and CH4 growth per year, based on GOSAT and NOAA reports.
Table 3. CO2 and CH4 growth per year, based on GOSAT and NOAA reports.
CO2 ppm/YearCH4 ppb/Year
YearNOAAGOSATNOAAGOSAT
20102.42.25.03.3
20111.72.05.112.2
20122.41.94.82.7
20132.42.35.65.2
20142.12.212.910.7
20153.02.110.15.3
20162.83.37.311.2
20172.22.56.68.0
20182.42.38.610.2
20192.52.69.99.4
20202.42.514.89.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lopez, F.P.A.; Zhou, G.; Jing, G.; Zhang, K.; Tan, Y. XCO2 and XCH4 Reconstruction Using GOSAT Satellite Data Based on EOF-Algorithm. Remote Sens. 2022, 14, 2622. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14112622

AMA Style

Lopez FPA, Zhou G, Jing G, Zhang K, Tan Y. XCO2 and XCH4 Reconstruction Using GOSAT Satellite Data Based on EOF-Algorithm. Remote Sensing. 2022; 14(11):2622. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14112622

Chicago/Turabian Style

Lopez, Franz Pablo Antezana, Guanhua Zhou, Guifei Jing, Kai Zhang, and Yumin Tan. 2022. "XCO2 and XCH4 Reconstruction Using GOSAT Satellite Data Based on EOF-Algorithm" Remote Sensing 14, no. 11: 2622. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14112622

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop