Greenhouse gases Observing SATellite 2 (GOSAT-2): mission overview

The Greenhouse gases Observing SATellite 2 (GOSAT-2) was launched in October 2018 as a successor to GOSAT (launched in 2009), the first satellite to specialize in greenhouse gas observations. Compared to the GOSAT sensors, the sensors of GOSAT-2 offer higher performance in most respects. The quality and quantity of data from observations are expected to be improved accordingly. The signal-to-noise ratio (SNR) is better in both the SWIR and TIR bands of TANSO-FTS-2, which is the main sensor of GOSAT-2. This improvement ultimately enhances the accuracy of greenhouse gas concentration analysis. Furthermore, because of the improved SNR in the SWIR band, the northern limit at which data are obtainable in high-latitude regions of the Northern Hemisphere in winter, where observation data have remained unavailable because of weak signal strength, has moved to higher latitudes. As better data are obtained in greater quantities, progress in carbon cycle research for high-latitude regions is anticipated. Moreover, the improvement of SNR in the TIR band is expected to be considerable. Particularly, the resolutions of the vertical concentration distributions of CO2 and CH4 have been improved drastically. The first function introduced for GOSAT-2 that is not in GOSAT is an intelligent pointing mechanism: a cloud area avoidance function using the in-field camera of TANSO-FTS-2. This function can increase the amounts of observation data globally and can improve the accuracy of CO2 emissions estimation and measurements of uptake intensity. The effects are expected to be strong, especially for the tropics because cumulus clouds are the most common cloud type. The intelligent pointing system can avoid the clouds effectively. Another important benefit of TANSO-FTS-2 is that the wavelength range of Band 3 of SWIR has been expanded for measuring carbon monoxide (CO). Because CO originates from combustion, it is used to evaluate some effects of human activities in urban areas and biomass burning in fields. Particularly, black carbon-type aerosols can be measured by the sub-sensor, TANSO-CAI-2, to assess biomass burning along with CO2 and CO by TANSO-FTS-2.


Introduction
The Greenhouse gases Observing SATellite (GOSAT), the first satellite dedicated to observing greenhouse gases, was launched in 2009 (Hamazaki et al. 2005;Kuze et al. 2009). It continues to make observations normally. The satellite's main sensor, Thermal And Near infrared Sensor for carbon Observation-Fourier Transform Spectrometer (TANSO-FTS), has observed dry air column average concentrations of carbon dioxide and methane and their upper air concentrations for over a decade. Based on these data, inversion analyses of these gases on a subcontinental scale have been performed (e.g., Takagi et al. 2011), contributing to the improvement of the estimation accuracy of the surface fluxes of the gases. Furthermore, observations of solar-induced fluorescence (SIF), an index related to the activity of terrestrial vegetation, were first conducted (Joiner et al. 2011;Frankenberg et al. 2011aFrankenberg et al. , 2011b. The TANSO-Cloud and Aerosol Imager (CAI), an image sensor onboard the satellite, has begun to analyze land aerosols that include black carbon (e.g., Hashimoto and Nakajima 2017).
Based on these achievements, GOSAT-2, a successor of GOSAT, was launched in October 2018. The signal-tonoise ratio of its main sensor, TANSO-FTS-2, has been improved considerably in both the Short Wavelength InfraRed (SWIR) and Thermal InfraRed (TIR) bands. A new spectral band has been added for observing carbon monoxide (CO). Because GOSAT has the important shortcoming that the probability of contamination caused by clouds is closely attributable to its large instantaneous field of view (IFOV), we adopted a new mechanism called intelligent pointing to concentrate the observation points in the clear sky region.
In addition, two observation directions (forward and backward) became available with CAI-2 to better avoid sun glint areas as CAI suffered from them to perform the cloud detection.
Furthermore, the addition of new spectral channels has improved the accuracy of observations of highly lightabsorbing aerosols such as black carbon.
Although GOSAT and GOSAT-2 sensors have large IFOV, they have been accumulating long-term global observation data on greenhouse gas concentrations in the atmosphere. On the other hand, for Orbiting Carbon Observatory-2 (OCO-2) ) and OCO-3 (Eldering et al. 2019), although the observation period is short, high spatial resolution data have been acquired. Combining such different types of sensors is expected not only to compare detection sensitivities, but also to create data sets with well-defined spatiotemporal resolution. On the other hand, regarding atmospheric pollutants such as black carbon, the synergy between GOSAT and GOSAT-2 on-board sensors with a small instantaneous field of view and sensors such as The TROPOspheric Monitoring Instrument (TROPOMI) (e.g., Torres et al. 2020) and Ozone Monitoring Instrument (OMI) (e.g., Jethva and Torres 2019), which can simultaneously measure many chemical substances and aerosols, is expected to bring more detailed knowledge of atmospheric chemistry.
Chapter 2 of this article presents detailed specifications of satellites and sensors. The data analysis algorithms of TANSO-FTS-2 and CAI-2 are described in Chapter 3. In Chapter 4, calibration of the radiance data observed by sensors and verification of level 2 data such as gas concentration and aerosol content are described. In Chapter 5, the scientific goals related to the global and regional carbon cycles, SIF measurements, and air pollution monitoring using GOSAT sensors are described based on the project's science plan (Imasu 2019).

Satellite bus system
Like its predecessor, GOSAT, the GOSAT-2 (known as "IBUKI-2") satellite is dedicated to the observation of greenhouse gases. It carries two mission instruments that provide higher performance than those of GOSAT: Thermal And Near infrared Sensor for carbon Observation-Fourier Transform Spectrometer-2 (TANSO-FTS-2) ("FTS-2") and Cloud and Aerosol Imager-2 (TANSO-CAI-2) ("CAI-2"). GOSAT-2 is flying in a sun-synchronous sub-recurrent orbit. Based on the operation of GOSAT and the mission requirements for GOSAT-2, some functions and performances have been improved as explained below.
The recurrence cycle of GOSAT is 3 days. The distance to the next orbit is about 900 km on the equator. However, the swath of TANSO-CAI is 1,000 km, so areas lying somewhere between the two orbits are observed with large off-nadir angles, which prevents improvement of the aerosol measurement accuracy. Therefore, it was decided that GOSAT-2 should pass between the GOSAT orbits so that the areas lying somewhere between two GOSAT orbits would be observed with a small off-nadir angle. Furthermore, the recurrence cycle which depends on the orbit height and inclination was set as 6 days. Two candidates of orbit heights and lower height orbit of 613 km were selected because the footprint size of FTS-2 would thereby be smaller. The GOSAT orbit was set in such manner that GOSAT flies over Tsukuba, Japan, where the Total Carbon Column Observing Network (TCCON) can acquire validation data. However, because many data acquired at Tsukuba are contaminated by clouds, the GOSAT-2 orbit was set to fly over another TCCON site, Lamont, where it is possible to anticipate sunny days.
In cases where only low power is available, as might occur with the breakdown of one solar array paddle, GOSAT was designed to halt observations in the TIR band of FTS and CAI. By contrast, GOSAT-2, with its larger-sized solar array paddle than that of GOSAT, is designed to continue operations of all mission instruments even after the breakdown of one solar array paddle. Table 1 and Fig. 1, respectively, present the principal specifications of GOSAT-2 and an artistic image of GOSAT-2 on orbit.

TANSO-FTS-2 2.2.1 Overview
TANSO-FTS-2 is a Fourier transform spectrometer (FTS), as is TANSO-FTS on GOSAT. The pointing mirror steers the line of sight to the observation area. After the scene flux passes the interferometer, it is divided into five observation bands by the dichroic filters. Additionally, bands other than the TIR are divided into two polarizations: P and S. They form an image on the detectors by the imaging optics. Table 2 presents the principal specifications of the TANSO-FTS-2. Also, Figs. 2, 3, and 4, respectively, present an outline view of the sensor unit of TANSO-FTS-2, a ray diagram, and details of the division of the band.   Spectral ranges (μm) (cm−1) Band 1: 0.75-0.77 (12,250) Band 1: 0.75-0.77 (12,200) Band 2: 1.56-1.69 (5900-6400) Band 2: 1.56-1.72 (5800-6400)

Improvements from TANSO-FTS on GOSAT
The increase in signal-to-noise ratio (SNR) and the addition of carbon monoxide measurements were necessary to satisfy mission requirements such as anthropogenic emission estimation over mega cities. Furthermore, the effective aperture size was extended from TANSO-FTS. The band range of the 2.0 μm band was extended to 2.3 μm.
To minimize the change of design from TANSO-FTS to the greatest degree possible, the basic design of the interferometer was retained. The aperture size was extended to the maximum size, which is available with the basic design of TANSO-FTS to increase the signal level. Based on GOSAT operations and the mission requirements, the spectral range of each band was optimized to decrease the noise level.
The range of each band was changed from that of GOSAT. The range of Band 1 was shifted slightly based on results of GOSAT data utilization. That of Band 2 was reduced to increase SNR by cutting the range which has not been used. That of Band 3 was extended to 2.3 μm to observe the carbon monoxide. TIR band was divided into two bands, Bands 4 and 5, to increase the SNR. Details about SNR for each spectral band are described by Suto et al. (2021Suto et al. ( , 2022 as a function of radiance intensity. For reference, comparison of SNR of detectors in the same wavenumber range for TANSO-FTS and TANSO-FTS-2 is shown in Tables 2, 3. Basically, apodization processing has not been applied to the spectra of all bands in the retrieval analysis. The IFOV is 15.8 mrad, which is the same as that of GOSAT, but the footprint size was decreased from 10.5 to 9.7 km because of the lower GOSAT-2 altitude. The Instrument Spectral Response Function (ISRF) is almost same as that of TANSO-FTS on GOSAT. The ISRF was evaluated using the Instrument Line Shape (ILS) target on TANSO-FTS-2, which is an integrated sphere with laser diodes of 768 nm and 1543 nm, respectively, for Bands 1 and 2, and a gas cell on the ground. The ISRFs of Bands 1 and 2 were evaluated on orbit using the ILS target. The FWHMs of Bands 1 and 2 were proved to be maintained at almost the same value as that obtained during ground testing by evaluation done for the commissioning phase. The ISRF of Band 2 measured on orbit is shown in Fig. 5.
Except for Band 1, for which the sampling resolution is different from those of other bands, all bands have almost the same FWHM. Details about ILS and its spectral dependences are described by Suto et al. (2021Suto et al. ( , 2022.
In addition, the pointing angle in the along track direction increased from ± 20 degrees to ± 40° to extend the sun glint observation area. GOSAT observes the standard mesh points, which are arranged with regularity. Up to 1500 points, a day can be replaced with specific points.
By contrast, GOSAT-2 is capable of observing all points as specific points. One can select the observation points freely within the capabilities of the pointing mirror moving speed.
The observation points are stored in the ground system as an observation table that is assigned to one recurrence cycle of 6 days. The ground system can hold up to five tables. The observation points in the table will be replaced with new specific points such as points required by users including PIs. The table of one-day observation points is uploaded to GOSAT-2: GOSAT-2 can then hold tables for two days. TANSO-FTS on GOSAT can select the scan duration, which is the time to select one interferogram from a menu of choices: 1.1 s, 2 s, and 4 s. However, the SNR of the data taken with the durations of 1.1 s and 2 s is relatively low. Therefore, TANSO-FTS-2 on GOSAT-2 makes observations with only 4 s scan duration. It is difficult to optimize the gain to fit the scene because TANSO-FTS has only two steps of the gain. TANSO-FTS-2 has 16 steps of the gain to take advantage of the dynamic range, accommodating very different radiance.

Intelligent pointing
Although TANSO-FTS with a point detector observes discretely distributed points, any clouds in the scene, however small, render the data invalid. TANSO-FTS-2 has a function of avoiding the clouds in the   Fig. 6; also, Fig. 7 presents actual pictures acquired by the monitor camera in orbit. It is possible to switch the intelligent pointing on and off. Detailed sequence of intelligent pointing. After completion of the preceding observation, the line of the sight (LOS) is steered to the next scheduled observation point. After the steering was stopped completely, after about 250 ms of the start of the steering, the image around the scheduled observation point is taken by the monitor camera (about 150 ms duration of ). Then, the image is transferred to the processing part; the clouds are detected with about 50 ms duration). Based on these results, the cloud-free area is determined. The LOS is steered to that area (about 50 ms duration of ). The next measurement starts at 650 ms after completion of the preceding observation This pointing mechanism including the algorithm was confirmed by trial manufacturing and flight test, combined with a FTS for airplanes. The intelligent mechanism was confirmed to have functioned as expected at about 45% of the approximately 3600 observation points. However, the RGB camera used for this experiment on the airplane had had about one-hundredth of the IFOV of the camera equipped with TANSO-FTS-2. Therefore, cloud detection was easy. For that reason, 45% seems to be larger than the actual FTS-2 case.
According to the simulation using the monitor camera data obtained by TANSO-FTS on GOSAT, the number of valid data is expected to increase to about two times that of GOSAT. The model used for the flight test of the intelligent pointing and photographs taken before and after cloud avoidance are presented, respectively, in Figs. 8 and 9.

TANSO-CAI-2
TANSO-CAI-2 is a push-broom-type imaging radiometer using a linear array sensor and band pass filter, similarly to TANSO-CAI. The LOS of TANSO-CAI-2 is tilted ± 20° in the along track direction relative to TANSO-CAI, for which the LOS points to the nadir direction. Also, TANSO-CAI-2 consists of five lens systems. Each lens system has two FOVs of forward and backward viewing.  Based on requirements to improve the measurement accuracy of aerosols, three new observation bands were added to TANSO-CAI-2 in addition to the four observation bands of TANSO-CAI. Furthermore, three bands are available for observation by forward and backward viewings. The swath of the 1.6 μm band of TANSO-CAI is 750 km, which is 25% narrower than other bands. Therefore, the swath of the 1.6 μm band of TANSO-CAI-2 was matched with those of other bands. Table 4 presents the principal specifications of TANSO-CAI-2. Also, Figs. 10 and 11 depict an outline view of the sensor unit of TANSO-CAI-2 without insulation, multilayer insulation (MLI), and a cross section of the lens system.
The TANSO-FTS-2 SWIR L2 Chlorophyll Fluorescence and Proxy-method Product includes the six individual retrieval results processed under the assumptions of cloud-free and aerosol-free conditions. Six wavenumber ranges, two each from the three TANSO-FTS-2 SWIR bands, correspond to six independent retrievals. From TANSO-FTS-2 Band 1, data of solar-induced chlorophyll fluorescence (SIF) and surface pressure are retrieved, respectively, using Fraunhofer lines in the wavenumber range, where the O 2 absorption and the O 2 A-band are negligibly small. From TANSO-FTS-2 Band 2, XCO 2 and XCH 4 are retrieved, respectively, from the 1.6 µm CO 2 band and 1.67 µm CH 4 band. From TANSO-FTS-2 Band 3, XCO 2 and XH 2 O, and XCO, XCH 4 , and XH 2 O are retrieved, respectively, using the 2.08 µm CO 2 band and 2.3 µm CO band. SIF is weak radiation emitted by chlorophylls during photosynthesis. It appears as the variation in fractional depth (filling in) of the Fraunhofer Fig. 9 Photographs taken before and after cloud avoidance for two cases (left and right). Green circle represents observation points. a Scheduled pointing, b processing the cloud and cloud free area (white circles are scheduled pointing position and red presents the cloud area), and c re-pointing to the cloud free area lines in the satellite-observed spectrum. Using this filling-in signal, SIF can be retrieved from satellite observed spectra (Frankenberg et al. 2011a(Frankenberg et al. , 2011b. Aside from SIF retrieval, five other retrievals use atmospheric absorption bands. Because these retrievals are processed under the assumption of cloud-free and aerosol-free conditions, the retrieved column-averaged dry-air mole fractions include some effects from optical path length modifications. If the column-averaged dry-air mole fraction of some gas species (X gas ) area is obtained with high accuracy using a model-based calculation or other methods, then that gas is useful as an optical path length proxy by comparing the retrieved X gas to prior knowledge. Because the degree of optical path length modification is regarded as similar to the nearby wavelength, the evaluated optical path length proxy is then useful to target gas retrieval. This method is designated as a "proxy method. " The CO 2 retrieved from the 1.6 µm CO 2 band is often used as the optical path length proxy to obtain XCH 4 from the 1.67 µm CH 4 band (e.g., Frankenberg et al. 2011aFrankenberg et al. , 2011bParker et al. 2011). Furthermore, the difference between the retrieved and a priori values of surface pressure obtained under the assumption of cloud-free and aerosol-free conditions is expected to reflect the degree of optical path length modification because of clouds and aerosols. Therefore, the TANSO-FTS-2 SWIR L2 Chlorophyll Fluorescence and Proxy-method Product includes retrieval results for all TANSO-FTS-2 SWIR observations. The TANSO-FTS-2 SWIR L2 Column-averaged Dryair Mole Fraction Product includes simultaneously retrieved XCO 2 , XCH 4 , XH 2 O, and XCO considering clouds and aerosols and using the several wavenumber ranges. This retrieval is designated as the "full-physics method" after the mode of the treatment of the optical path length modification. In the full-physics method, the optical path length modification attributable to clouds and aerosols is described by atmospheric scattering processes using the forward model. Being different from the proxy method, the full-physics method has the important benefit that no limitation exists for the gas species to retrieve. However, its retrieval error depends strongly on how accurately the forward model can represent the optical path length modifications attributable to clouds and aerosol scattering. Therefore, TANSO-FTS-2 SWIR observations that are affected strongly by clouds and/or aerosols are not processed using the full-physics method. Retrieval results for converged cases in the retrieval processing are stored in the TANSO-FTS-2 SWIR L2 Column-averaged Dry-air Mole Fraction Product.
Differences in the wavenumber range of observed spectra, state vector elements, and treatment of clouds and aerosols can be treated as differences in retrieval settings. The Maximum A Posteriori (MAP) retrieval method (Rodgers 2000) is used for all retrieval processing. Figure 12 presents examples of column-averaged dryair mole fraction products of CO 2 , CH 4 , and CO as averaged values in each 2.5° mesh area for March 2019.

Maximum posterior (MAP) analysis method
According to the MAP analysis method, the optimal estimate of state vector x, which represents unknown parameters to be retrieved, is the value of x which minimizes cost function J(x), as expressed in Eq. (1).
Therein, y represents the observed radiance spectrum, F(x, c) denotes a forward model, i.e., the spectrum computed with the state vector x and other necessary variables c as inputs, x a is the a priori state vector, S a represents the variance-covariance matrix of x a , and S ε is the variance-covariance matrix of observation error and forward model.
In TANSO-FTS-2 SWIR L2 processing, solution x is obtained iteratively from Eq. (2). (1) (2) In that equation, subscript i denotes the number of iterations, K = ∂F(x)/∂x is a Jacobian, D is a diagonal matrix introduced to stabilize the iteration process, and λ represents a coefficient of which the value is adjusted so that the value of the cost function decreases in each step of the iteration.
Here, the posterior error covariance matrix S and the averaging kernel matrix AK are given, respectively, by Eqs. (3) and (4).
As described earlier, the elements of state vector x and the wavenumber range of the observed spectrum y change for each retrieval process, but the same a priori values are used for each type of state vector element ( Table 5). The column-averaged dry-air mole fraction of target gas species X x is obtained by the vertical weighted integral of the gas concentration profile x x obtained using MAP analysis (Eqs. 5 and 6). Here, subscript x denotes the elements corresponding to the gas species to be retrieved. To describe the gas concentration profile, we use the layer-averaged (partial column-averaged) dry-air mole fraction of a gas species in a vertical grid that evenly divides the atmosphere by atmospheric pressure.
(3) Therein, subscript j denotes a vertical grid; also, w dry is a partial column amount of dry air.
The posterior error σ x 2 , column averaging kernel a x , and degrees of freedom for signals (DFS) with respect to X x (DFS x ) are given, respectively, by Eqs. (7-9).
In those equations, S xx and AK xx , respectively, represent the sub-matrices of S and AK, with regard to the retrieved gas species.

Pre-launch evaluation results based on the designed value of TANSO-FTS-2
Compared to TANSO-FTS, TANSO-FTS-2 SWIR bands show higher SNR values. Its spectral coverage is expanded and optimized for CO and SIF observation (Table 6). This section presents our investigation of the effects of these differences in TANSO-FTS-2 SWIR design on the posterior error (Eq. 7), a column averaging kernel (Eq. 8), and DFS (Eq. 9).
For the evaluation, the three cases shown in Table 7 are considered. The effectiveness of the improved SNR and spectral coverage is evaluated, respectively, by comparing cases 1 and 2, and cases 2 and 3. To make the comparison conditions uniform, the actual observation geometric conditions of TANSO-FTS are used; also, the state vector elements and a priori values used to derive the TANSO-FTS-SWIR L2 product version 02  are used. However, the variance-covariance matrix S a for the a priori value of the gas profile is a diagonal matrix. Also, the diagonal elements are set as (10%) 2 for CO 2 , CH 4 , and H 2 O, and (100%) 2 for CO. Furthermore, the variance-covariance matrix for observation error S ε is a diagonal matrix. For simplicity, its diagonal elements are assumed to consist only of shot noise, which is proportional to the square root of the total number of photons incident on the detector. Under these conditions, evaluation is conducted using 15,751 observation geometric conditions of TANSO-FTS in January, April, July, and October 2011. Figure 13 presents an example of the averaging kernel matrix (Eq. 4) and the column averaging kernel (Eq. 8).
The averaging kernel matrix and the column averaging kernel, respectively, represent the sensitivity of the retrieved profile and X gas to the actual ("true") profile of target gas. With the improvement of SNR (cases 1 and 2) and the expansion of the wavenumber range (cases 2 and 3), the averaging kernel matrix is changed, but its effect on the column averaging kernel is small. As the SNR improves, the column averaging kernel tends to approach unity, albeit slightly. In any case, the column averaging kernel shows a value that is close to unity over the whole atmospheric layer, indicating that it has sufficient sensitivity to X gas . Figure 14 portrays histograms of DFS and posterior error. With the improvement of SNR and the expansion of the wavenumber range, DFS increases and posterior error decreases. Regarding XCO 2 , almost no difference exists between cases 2 and 3 because the wavenumber range added to case 3 does not include  the CO 2 absorption band. Because DFS exceeds unity in most cases, it can be confirmed that the design result of TANSO-FTS-2 has performance that is sufficient for retrieving X gas . The improvement of SNR contributes greatly to reduction of the posterior error. The posterior errors of XCO 2 , XCH 4 , and XH 2 O are suppressed, respectively, to about 0.85, 0.55, and 0.80 times for GOSAT-2 (case 3) compared to GOSAT (case 1). However, because this evaluation is based on the designed value of TANSO-FTS-2, there is a simplified aspect to the evaluation. It is appropriate to regard it as a result obtained under ideal conditions. It is necessary to make a more realistic evaluation using actual characteristics of TANSO-FTS-2.

Greenhouse gases and temperature
We have adopted a nonlinear MAP method with linear mapping (Rodgers 2000) to retrieve vertical profiles of greenhouse gases and temperature from the thermal infrared (TIR) bands (Bands 4 and 5) of TANSO-FTS-2 on GOSAT-2, as is true also for the operational data processing of the data of TANSO-FTS on GOSAT (Saitoh et al. 2009. Regarding retrieval, we have introduced linear mapping between vertical layer grids for radiative transfer calculation (full grids) and grids for retrieval calculation (retrieval grids). The retrieval grids, which are defined as a subset of the full grids, are fixed for all retrieval processing. The main target gases of the TIR bands of TANSO-FTS-2 are CO 2 , CH 4 , H 2 O, O 3 , and N 2 O. Vertical profiles of CO 2 , O 3 , and temperature are retrieved simultaneously from TIR Band 5 spectra at around 15 μm (690-750 cm −1 and 790-795 cm −1 ), 10 μm (930-990 cm −1 ), and 9 μm (1040-1090 cm −1 ) absorption bands. Vertical profiles of CH 4 , H 2 O, and N 2 O are retrieved simultaneously from TIR Band 4 and 5 spectra at around 7 μm and 8 μm absorption bands (1140-1370 cm −1 ).
A priori and initial values of meteorological parameters (temperature and water vapor) and the ozone field used in TIR retrieval were referred from the Japan Meteorological Agency (JMA) Japanese 55-year Reanalysis (JRA-55). Profiles of a priori and initial concentrations of CO 2 and CH 4 were taken from outputs of the Nonhydrostatic ICosahedral Atmospheric Model (NICAM)based Transport Model (TM) (Niwa et al. 2011), which are the same as those used for SWIR retrieval. Profiles Fig. 13 Example of the averaging kernel matrix (colored lines) and the column averaging kernel (black line). They are calculated under the diffusive reflectivity of 0.3 and the solar zenith angle of 30°. It corresponds to cases 1, 2, and 3 from the top, and CO 2 , CH 4 , H 2 O, and CO from the left. The vertical axis corresponds to the layer number of the atmospheric layer. It is counted from the top-of-the-atmosphere toward the surface. Colors of the averaging kernel matrix also correspond to the layer number of the atmospheric layer, in the order of purple (top-of-the-atmosphere), blue, green, orange, and red (surface) of a priori and initial concentrations of N 2 O are taken from outputs of the Atmospheric Chemistry Transport Model (ACTM) (Ishijima et al. 2010). A priori covariance matrices for these target parameters are set as diagonal matrices with fixed values as elements throughout a vertical profile. Similarly to GOSAT, variation with a standard deviation of approximately 2.5% for CO 2 is assumed for an entire vertical profile . For TIR retrieval, we used the Line-By-Line Radiative Transfer Model (LBLRTM) (Clough et al. 2005) and the Highresolution Transmission Molecular Absorption Database (HITRAN) (Rothman et al. 2005(Rothman et al. , 2009(Rothman et al. , 2013 for lineby-line calculations of optical thickness for each atmospheric layer of the full grids. Figure 15a shows typical averaging kernel functions (Rodgers 2000) of CO 2 and CH 4 retrievals from TANSO-FTS-2 TIR spectra for atmospheric conditions at low latitudes. For calculating the averaging kernels, we have assumed signal-to-noise ratio (SNR) values for Bands 4 and 5 of TANSO-FTS-2, which were evaluated at the TANSO-FTS-2 Critical Design Review (CDR). The CO 2 and CH 4 averaging kernel functions for the TANSO-FTS TIR band calculated under the same atmospheric conditions are depicted in Fig. 15b for reference. As depicted in Fig. 15, the TANSO-FTS-2 TIR spectra include more information related to CO 2 and CH 4 concentrations than TANSO-FTS. The much higher sensitivity to CH 4 concentration of TANSO-FTS-2 is attributable to the high SNR values of TANSO-FTS-2 Band 4 (1045 at 1300 cm −1 ).
The abundances of the species are much less than the main absorbers (H 2 O, CO 2 , CH 4 , and O 3 ). Moreover, their absorption signals are generally weaker than those of the main absorbers. Retrievals of the species are subject to interference from these absorbers and from atmospheric temperature, surface temperature, and surface emissivity because of the overlaps of the absorption bands and water vapor continuum. Therefore, the retrievals must incorporate consideration of these influences. The retrieval algorithms used for earlier studies particularly addressing the species are typically based on the MAP solution (Rodgers 2000), which is spectral fitting using the Bayesian theorem, similar to that used for greenhouse gases. In these algorithms, the interferences are retrieved simultaneously or separately retrieved in advance to input the results into the retrieval of the main target. Because of weak absorption and interference, the information contents of the minor species within the observed spectra are small; also, the retrievals have low ability to resolve vertical profiles. In fact, TANSO-FTS-2 has the highest spectral resolution among the current operating space-borne nadir-viewing TIR sounders, with an SNR much higher than that of TANSO-FTS (Suto et al. , 2022. Therefore, the information contents and the abilities to resolve vertically for the species are expected to be enhanced greatly by TANSO-FTS-2.

Clouds and aerosols
The cloud screening process is important for analyzing gas components such as CO 2 . For analysis of the TANSO-FTS-2 (SWIR) band data, the cloud flag of CAI-2 and the cloud information obtained from the SWIR band data itself are used together. Similarly, for analysis of TANSO-FTS-2 (TIR) band data, the cloud flag of CAI-2 and the cloud information obtained from the TIR band data itself are used together. However, because CAI-2 is not operated during nighttime, the cloud information obtained from the TIR data itself will be used alone.
For cloud screening based on the TIR band data, independent methods of three types are used together. The first is the threshold method for brightness temperature. In this case, when the value of the difference between the observed brightness temperature and the ground surface temperature in the atmospheric window region exceeds a certain threshold value, it is identified that the field of view is covered with clouds. At this time, meteorological analysis data are used as the ground surface temperature, but the surface emissivity estimated from data such as the vegetation index and snow and ice occupancy in the field obtained from MODIS data is used. For this method, the cloud top temperature and height in the field of view are derived as products.
The second method is a slicing method using a CO 2 absorption band. Here, the method described by Someya et al. (2016) is used, which employs a pseudo-spectral channel re-constructed by combining spectral channels having the same peak altitude of the weighting function. In addition to clouds, this method can detect optically thick aerosol layers such as yellow sands (Someya et al. 2019). Using this method, the cloud top height in the field of view is derived as a product.
The third method is the so-called split window method proposed by Inoue (1985). Using this method, the cloud top temperature, the effective cloud particle size, thermodynamic phase of the particles, and their optical thicknesses are derived as products.

TANSO-CAI-2 3.3.1 Cloud discrimination
The presence of clouds in the instantaneous field of view (IFOV) of the TANSO-FTS-2 engenders incorrect estimates of greenhouse gas concentrations. Consequently, an important role of TANSO-CAI-2 is to perform cloud discrimination to identify and reject cloud-contaminated TANSO-FTS-2 data. Improving the cloud discrimination accuracy necessitates refinement of the existing cloud discrimination algorithm for TANSO-CAI on GOSAT: the CLoud and Aerosol Unbiased Decision Intellectual Algorithm used with CAI (CLAUDIA1-CAI) (Ishida and Nakajima 2009). A new cloud discrimination algorithm (CLAUDIA3) using a support vector machine (SVM) (Vapnik and Lerner 1963) was developed (Ishida et al. 2018). Primary verification results show that the accuracy of CLAUDIA3-CAI is greater than that of CLAU-DIA1-CAI (Oishi et al. 2016). CLAUDIA3-CAI is in the improvement stage because difficulties to be resolved were also given in the article. A summary of CLAUDIA3-CAI is presented below. CLAUDIA1 performs cloud discrimination using thresholds that are set based on experience. CLAUDIA3 uses SVM to decide the thresholds objectively using multivariate analysis. SVM is a supervised pattern recognition method using TANSO-CAI-2 Bands 3-5 (forward) or Bands 8-10 (backward).
First, it determines the following items using training samples of typical clear and cloudy pixels: (1) a decision function to discriminate between two classifications (clear and cloudy), (2) the thresholds, and (3) the support vectors, which are training samples specified by the decision function. The support vectors are decided in a high-dimensional feature space of the training samples. Next, it performs cloud discrimination using the decision function, thresholds, and support vectors it determined. The flow of CLAUDIA3-CAI is depicted in Fig. 16. The results are expressed as integrated clear-sky confidence levels (integrated-CCL). Integrated-CCL of 0 expresses that the pixel is cloudy; 1 signifies that the pixel is cloudfree. Ambiguous pixels between cloudy and cloud-free are described using numerical values from 0 to 1. For CLAUDIA3-CAI, the integrated-CCL of 0.5 corresponds to the separating hyperplane of clear support vectors and cloudy support vectors (Oishi et al. 2017.
Users can adjust the boundaries which separate cloud and clear-sky areas with the integrated-CCL thresholds. They can also discriminate the presence of clouds. Figure 17 presents examples of an input CAI Level 1B product image and its cloud discrimination result obtained using CLAUDIA3-CAI.
TANSO-CAI-2 has forward-(+ 20°) and backwardlooking (− 20°) bands. Because cloud discrimination will be performed in each direction, the difference between forward-looking and backward-looking bands of cloud discrimination has been discussed (Oishi et al. 2017. The use of the forward and backward cloud discrimination during processing for calculating greenhouse gas concentrations (Oishi et al. 2015) is examined.

Aerosol products
TANSO-CAI-2 takes two-directional observations in forward and backward directions. Ten observation bands and seven wavelengths exist, as presented in Table 3. The spatial resolution (IFOV) is 920 m at a 1630 nm wavelength; it is 460 m at the other wavelengths. Actually, TANSO-CAI-2 is characterized by having two wavelengths in the ultraviolet (UV) region of 340 nm and 380 nm that are sensitive to aerosol light absorption characteristics. Therefore, the retrieval accuracy of aerosol absorption is expected to be enhanced. Particle size information can also be estimated using data of the visible-near-infrared band.
The TANSO-CAI-2 Level 2 aerosol products are intended to derive the aerosol optical thickness at a 1630 nm wavelength (AOT1630), which is necessary for improving the CO 2 derivation accuracy by TANSO-FTS-2. In addition, the retrieved aerosol products are the aerosol optical thickness at a 550 nm wavelength (AOT550) as an amount of the atmospheric aerosols, the aerosol Ångström exponent (AE) as particle size information, ePM2.5 as an equivalent amount of PM2.5 to elucidate the anthropogenic air pollutants, and the black carbon volume fraction (BCF). The derived resolution of the product is 2 km × 2 km in the Asian region, and 5 km × 5 km otherwise (Fig. 18). Table 8 presents the TANSO-CAI-2 L2 aerosol property products and the required accuracy.

Aerosol retrieval algorithm
The multi-wavelength and multi-pixel method (MWPM) (Hashimoto and Nakajima 2017) is applied as the algorithm to retrieve TANSO-CAI-2 Level 2 aerosol products. The MWPM was developed for accurate derivation of aerosol characteristics, even in areas with inhomogeneous surface reflectance, such as urban areas for which it was difficult to derive aerosol characteristics from satellite images. This method expands the principle of the Kaufman Neutral method (Kaufman 1987) and simultaneously retrieves aerosol optical thickness and light absorption characteristics. It does so by simultaneously analyzing the observation points in the region where reflectance in the spatial direction changes.
The MWPM uses a combination of surfaces with different reflectance on the assumption of an aerosol that is distributed smoothly over the surfaces, as depicted in Fig. 19. Equation 10 expresses the apparent reflectance (R) by a single scattering approximation on a Lambertian ground surface as shown below.
Therein, A g stands for the surface reflectance; m and m 0 , respectively, represent inverses of the cosines of incident and scattered zenith angles. τ, ω, and P(Θ), respectively, denote the AOT, single scattering albedo (SSA), and phase function of the atmospheric air mass at scattering angle Θ. Also, γ + and γ − , respectively, represent the downward and upward fractions of light scattering. The neutral reflectance (R c ) is a case in which R c = R = A g in Eq. 10 corresponding to the middle case of Fig. 19 when a path radiance (first term on the right-hand side, RHS) and surface-reflected radiance (second term of the RHS) are canceled out, irrespective of the value of τ. For low-reflectance and high-reflectance surfaces, aerosols in the atmosphere, respectively, increase and decrease the apparent reflectance at the top of the atmosphere (TOA) with R c as a boundary, depending on the light absorption characteristics of the aerosol particles. Presuming that AOT and other aerosol properties in the atmosphere do not change considerably in the horizontal direction, an optimal set of aerosol parameters (e.g., AOT, SSA, A g ) is found. Therefore, AOT and SSA are obtained by solving a set of radiative transfer equations (RTEs) using the spatial differences of surface reflectance (Hashimoto and Nakajima 2017).
The optimal method and smoothing constraint for atmospheric parameters are combined to solve the inverse problems for aerosol properties. Bayes' theorem is applied as an optimal method of problem solving to maximize the posterior probability density function P(u|R) as shown below.
Therein, P(u) is a probability density function of event u. Also, P(u|R) is a conditional probability density function where event u occurs under event R. Then, the smoothing constraint for atmospheric aerosol properties (i.e., a condition of minimizing the square of the second derivative of the aerosol spatial distribution with smoothing coefficients) is added to the optimal method.
(10) In those equations, γ denotes the smoothing factor, x stands for a certain pixel, and ∆x is the distance between two adjacent pixels to calculate the parameter's smoothness. Also, R is the measurement vector (i.e., observed apparent reflectance) at each measurement band. In addition, u and u a , respectively, stand u ij = lnτ fine lnτ seasalt lnτ coarse lnf BC lnA g,band1 lnA g,band2 · · · lnA g,band10 T ij , (16) u a,ij = lnτ a,fine lnτ a,seasalt lnτ a,coarse lnf BCa lnA ga,band1 lnA ga,band2 · · · lnA ga,band10 T ij , for the state and a priori vectors. Actually, τ fine , τ seasalt , and τ coarse , respectively, denote fine, sea, and coarse mode of aerosol optical thickness at 550 nm wavelength: AOT550fine, AOT550sea, and AOT550coarse. Total AOT550 is the sum of these AOT550s. f BC is the volume fraction of black carbon (BCF). A g,bandN is the surface albedo at band N. Subscript "a" represents the a priori value. Over the ocean, the wind speed at 10 m above sea level is used instead of the surface albedo. Five bands of forward or backward viewing are used, respectively, for analysis to avoid the sun glint region. S e and S a represent the covariance matrices of measurement errors and a priori values. B is the boundary condition for the state vector. f(u) denotes a forward model (i.e., radiative transfer model). The forward model incorporates the radiative transfer model, accelerated RSTAR using a neural network method (Takenaka et al. 2011). Developed at the University of Tokyo, RSTAR (Nakajima andTanaka 1986, 1988;Sekiguchi and Nakajima 2008) can simulate radiation fields in the atmosphere-land-ocean system. Both AOT550 and BCF are derived by solving Eq. (12).

Fig. 18
Image of CAI-2 aerosol product with resolution of 2 km × 2 km for the Asian region corresponding to GOSAT-2 passes 4-35 and frames 13-20, and 5 km × 5 km for other areas For analyses, the region to be analyzed is divided into several small regions designated as subdomains (e.g., 5 × 5 pixels). Aerosol properties at each pixel in the subdomains are solved sequentially with a boundary condition.

Aerosol retrieval procedure
Figure 20 portrays the procedure for deriving TANSO-CAI-2 aerosol characteristics. First, using the observed data for one month before and after the product retrieval date (accounting for a total of 11 days), the surface reflectance is derived by the minimum and second minimum method (Fukuda et al. 2013). Then, using the surface reflectance and meteorological field data, aerosol properties such as total AOT550 and AOT1600 and BCF are retrieved by the MWPM. The AOT at each wavelength including 1600 nm is also derived from the retrieved AOT550s and the assumed aerosol model. Then, the Ångstrom exponent (AE) and ePM2.5 are calculated using the retrieved AOTs and BCF.
Using the wavelength characteristics of retrieved AOT at wavelengths of 441, 546, 672, and 865 nm, AE is calculated as where α is the Ångström exponent, τ i stands for the total amount of AOT at wavelength i, expresses the wavelength, and n represents the number of wavelengths used to calculate AE.
ePM2.5 is derived from the retrieved BCF, fine mode of AOT550, and the assumed aerosol particle size as Therein, m is the aerosol height information parameter at unit meter. Also, ρ sulfate and ρ BC , respectively, denote the sulfate and black carbon particle densities. In addition, r stands for the particle radius; dv fine /dlnr represents the volume size distribution of aerosol particles in fine mode.

Overview
Sensor calibration is performed by pre-flight test (PFT) before the launch, and again on-board, with vicarious calibrations after launching. The TANSO-FTS-2 SWIR band and TANSO-CAI-2 observe sunlight reflected from the Earth's surface. The PFT is performed to evaluate the sensor output to radiance by the integrating sphere light source, which is calibrated using a standard radiometer. Consequently, the observation outputs can be converted into radiance by application of the radiometric conversion factors determined by the PFT results. Observation radiance of the TANSO-FTS-2 TIR band is achieved by conversion from the sensor output calibrated with deep space as a cold target and the onboard blackbody as a hot target on orbit. In the PFT, it is necessary to obtain a nonlinear response of detectors and polarization characteristics of the optics, which is difficult to determine on orbit.

Onboard calibration
The onboard calibration plan for TANSO-FTS-2 is presented in Table 9. Integrating sphere lasers are carried to confirm the instrument line shape (ILS) of Band 1 and Band 2. Thereby, the temporal stability on orbit can be tracked because of the PFT. Double-sided diffusers for solar irradiance calibration are operated to evaluate the degradation of the radiometric response of the SWIR band by observing both sides alternately. Lunar calibration for radiometric calibration for the SWIR band is planned for observation at the phase angle of + 7° waning (17)  near a Full Moon. The appearance of the Moon is half the size of the TANSO-FTS-2 field of view (FOV). Radiometric calibration of the TIR band is performed on orbit by deep space and onboard blackbody two times dayside and four times nightside. The phase transit cell is carried to confirm the performance of the thermometer attached to the onboard blackbody. The FOV camera is operated to take photographs before and after intelligent pointing for the detection of cloud information of each sounding. The onboard calibration plan for the CAI-2 is shown in Table 10. The TANSO-CAI-2 onboard calibration is operated by lunar calibration and nighttime calibration; it is not equipped with onboard calibrator. The lunar calibration for radiometry is planned to observe at phase angle of + 7° waning near Full Moon by along-track (AT) scanning with center pixels. The GOSAT-2 satellite is targeted to the Moon at a pitch angle of + 20° for forward-looking bands, and at − 20° for backward-looking bands. The lunar calibration for pixel-to-pixel response correction is operated by crosstrack (CT) scanning with all pixels at a phase angle of − 7° waxing near the Full Moon. Nighttime calibration for dark-time output is planned to observe the dark night ocean near a New Moon.

Vicarious calibration
Vicarious calibration after the launch is shown in Table 11 for the TANSO-FTS-2, and in Table 12 for the TANSO-CAI-2. A field campaign is conducted in summer at Railroad Valley (RRV) Playa, Nevada, USA,   by measuring the surface reflectance coincident with GOSAT-2 target observation for vicarious radiometric calibration of the TANSO-FTS-2 SWIR band and TANSO-CAI-2. The RRV experiment is a collaboration with NASA OCO-2/-3 team and other greenhouse gas projects. Cross-calibration for the TANSO-CAI-2 VNIR band is evaluated at the dark ocean target using AER-ONET Ocean Color (OC), and at a bright land desert target using VIIRS observed radiance. Those for UV bands 1 and 6 are evaluated using Marine Optical BuoY (MOBY) water leaving radiance (nLw). Inter-satellite comparisons with data from similar sensors, especially those of GOSAT, are planned that the TANSO-FTS-2 targets at some GOSAT TANSO-FTS observation points during 1 or 2 days in the GOSAT-2 6-day revisit, both at close orbits. The TANSO-FTS-2 SWIR band is compared with OCO-2 and TROPOMI, and the TIR band with AIRS, IASI and CrIS using simultaneous observation data. The TANSO-CAI-2 VNIR band is compared with VIIRS, whereas the UV band is compared with OMPS for Band 1, and with OMI for Band 6. The temporal stability of the radiometric response is also evaluated by monitoring desert sites recommended by Committee on Earth Observation Satellites (CEOS). Geolocation is evaluated by comparison with a reference image after launching. The TANSO-FTS-2 geolocation is evaluated using a FOV camera image. The TANSO-CAI-2 geolocation is evaluated as a standard in Band 4 for forward-looking bands, and in Band 9 for backward-looking bands. Other bands are evaluated as band-to-band registration relative to the standard band.

Purpose and outline of validation
The product data observed using satellite remote sensing techniques are necessary for evaluation with other independent observation data for their scientific and political uses. Level 2 and higher level products will be validated and compared, respectively, whereas level 1 products will be calibrated. For validation, the product is evaluated with other independent observation data with better precision and well-known bias. For comparison, the product is evaluated with other independent observation data with equivalent or less precise or model data. During calibration, a physical quantity such as radiance is evaluated with independent observation data having better precision. Well-known bias and characteristics of instruments are assessed using various measurement methods.
A procedure for validation is the following, for example, selection of the validation data, quality check and screening of validation and satellite data, colocation, summary statistics, correlation analysis, and bias correction. For selection of the validation data, similarity, precision, and bias of data should be considered. Similarity signifies a similar physical quantity such as the columnaveraged dry mole fraction such as XCO 2 or a profile of the dry air mole fraction to satellite data. If each group of data differs in quantity for validation and satellite data, then satellite or validation data conversion should be made of into same physical units. During quality checking and validation of satellite data, anomalies caused by instruments, weather conditions, and some phenomena should be discriminated according to the purpose of the validation level, such as preliminary, routine, geophysical, advanced, and special-purpose detection for gas emission using statistical techniques. After quality checking, they should be removed. For colocation, various methods are used, such as geographical (Inoue et al. 2013), T700 (Wunch et al. 2011a), model-based , and geostatistical (Nguyen et al. 2014). They are summarized in a report by Nguyen et al. (2014). Then, validation and satellite data are made summary statistics (mean and median with or without weight) with or without averaging kernel and prior corrections. Next, for correlation analysis, validation and satellite data are compared with or without weights. When weights for both data are considered for correlation analysis, the so-called York method is used. Characteristics such as bias and scatter are eventually estimated. If there are validation criteria for each product as a satellite mission and statistics achieve them, the products are recognized as validated. However, statistical examination using this procedure can be applied for validation. If possible, satellite data are recommended to be corrected against their biases using empirical linear regression. This motivation and various methods are summarized in a report by Inoue et al. (2016). Normally, correlation analyses between the satellite biases and simultaneously retrieved auxiliary parameters are made. These correlations are used to biascorrect the GOSAT data, thereby removing these spurious correlations.

Overview of the GOSAT-2 products validation plan
As validation sources for XCO 2 , XCH 4 , XCO, and XH 2 O retrieved from the SWIR spectra onboard GOSAT-2 TANSO-FTS-2, TCCON is a primary validation source. TCCON is a worldwide network of ground-based FTSs that provide time series of column-averaged abundances of various atmospheric constituents. These constituents including CO 2 , CH 4 , N 2 O, HF, CO, H 2 O, and HDO are retrieved from near-infrared solar absorption spectra similar to GOSAT-2 using a nonlinear least-squares fitting algorithm known as GFIT (Wunch et al. 2010(Wunch et al. , 2011a(Wunch et al. , 2011b. The TCCON data are also useful for the validation of satellite vertical profiles of CO 2 and CH 4 after conversion into the same physical quantities. In addition to the TCCON data, the obtained Network for the Detection of Atmospheric Composition Change (NDACC) FTS data are useful. The Infrared Working Group (IRWG) represents a network of infrared Fourier-transform spectrometers that is part of NDACC. An overview of NDACC is described in a report by De Mazière et al.
: It is a multi-national collection of more than 20 high-resolution spectrometers that regularly record the atmospheric absorption spectrum from sites distributed from pole to pole. The absorption spectra are used to retrieve concentrations of a number of the gaseous atmospheric constituents, including O 3 , HNO 3 , HCl, HF, CO, N 2 O, CH 4 , HCN, C 2 H 6 , and ClONO 2 . Because the NDACC FTS data are column abundances and vertical profiles of constituents, they are useful to validate SWIR and TIR products of GOSAT-2. In addition to FTS data, DIAL (e.g., CO 2 , Sakaizawa et al. 2009;O 3 , Uchino et al. 2017), other portable solar light observation instruments such as portable FTS network called COllaborative Carbon Column Observing Network (COCCON) (Frey et al. 2019), Optical Spectrum Analyzer (OSA, Kobayashi et al. 2010;Kawasaki et al. 2012), fiber-etalon sun-photometer (Kobayashi et al. 2010;Kawasaki et al. 2012), and a laser heterodyne radiometer (Wilson et al. 2014) are used. Profiles of CO 2 , CH 4 , H 2 O, and other gases retrieved from the TIR spectra can be validated using aircraft data such as the Comprehensive Observation Network for TRace gases by AIrLiner project (CONTRAIL) (Machida et al. 2008), In-service Aircraft for a Global Observing System (IAGOS) core (Filges et al. 2015), the NOAA Earth System Research Laboratory/Global Monitoring Division (ESRL/GMD) (Xiong et al. 2008;Sweeney et al. 2015), the US Department of Energy (DOE) (Biraud et al. 2013;Schmid et al. 2014), and CO 2 sonde (Ouchi et al. 2019), and AirCore (Karion et al. 2010). Their profiles can also be useful for validating the SWIR products by conversion into X gas .
Surface measurements of greenhouse gases (e.g., WDCGG 2022) can also provide a useful comparison for SWIR and TIR products.
Solar-induced fluorescence (SIF) provides reliable information of plant physiological functioning at a large spatial scale. Actually, SIF can be validated using groundbased (e.g., Julitta et al. 2016) and aircraft (e.g., Garzonio et al. 2017) SIF data recorded with a hyperspectral spectrometer. Moreover, SIF can be compared with satellite data provided by GOME-2, GOSAT, OCO-2, TanSat, FLEX, and others, in addition to vegetation model data.
In situ data and ground-based remote sensing data from AERONET and SKYNET (Nakajima et al. 2020a(Nakajima et al. , 2020b) are used to validate PM2.5 and black carbon (BC) from TANSO-CAI-2. Satellite data from MODIS, Earth-CARE, and others can also be useful for comparison. As advanced validation to improve the GOSAT data quality further, investigation can be made about the effects of aerosols and thin cirrus clouds on the GOSAT-2 product over Rivulets, Tsukuba, Saga, Burgos, and Lauder using FTS, lidar, and skyradiometer.
Campaign observation is expected to be important where validation data are sparse. Under such circumstances, combinations of portable instruments and aircraft and satellite observation are recommended according to the campaign purposes.

Ongoing validation activities of the GOSAT-2 products
After starting the L2 processing, we have been conducting validation analysis using TCCON data for XCO 2 , XCH 4 , XCO, and XH 2 O. We have made evaluations of the influences of aerosol and thin cirrus on them using data obtained from skyradiometer and lidar. We also have started intercomparisons of XCO 2 , XCH 4 , and XCO with OCO-2, TROPOMI, and model data. The CO 2 profile has been compared with those taken from CONTRAIL and model data. The SIF data have been compared with those from satellite data from OCO-2, etc. Moreover, the AOT data from TANSO-CAI-2 have been validated using the AERONET and SKYNET data; also, BC and PM2.5 have been compared with surface BC and PM2.5 data. Some of the validation analyses have been used for review processes in public release.

Global inversion
At the moment, available estimates of surface fluxes of greenhouse gases (GHGs) still have marked degrees of uncertainty because our understanding of flux mechanisms remains insufficient. Many gaps exist in in situ observation networks. Satellites, including GOSAT-2, that observe atmospheric mole fractions of GHGs are expected to reduce those uncertainties because of their wide spatial coverages. Generally, an atmospheric transport model is used to connect surface fluxes with observed atmospheric mole fractions. Based on those findings, an inverse analysis is performed to estimate surface fluxes. An inverse analysis is useful in that it provides quantitative estimates of spatiotemporal variations of surface fluxes consistently with observed variations of atmospheric mole fractions. In a typical global inverse analysis, surface fluxes are estimated on a much larger scale (a few 100 km) than a footprint of satellite observations (a few to 10 km), because of the computationally demanding simulations of an atmospheric transport model. However, because satellite retrievals represent partially column-averaged mole fractions and because they consequently have large representativeness, satellite observations can be major input data of a global inverse analysis. Such a global inverse analysis is expected to provide valuable implications of large-scale (sometimes interhemispheric scale) spatial patterns and seasonal variations of surface GHG fluxes.
In the GOSAT-2 project, the Nonhydrostatic ICosahedral Atmospheric Model (NICAM)-based Transport Model (NICAM-TM) (Satoh et al. 2014;Niwa et al. 2011) is used to prepare prior estimates of retrieval calculations (i.e., the L2 product). Furthermore, the inverse analysis based on NICAM-TM, which is designated as NICAM-based Inverse Simulation for Monitoring (NIS-MON) CO 2 and CH 4 , provides surface flux estimates of CO 2 and CH 4 using each L2 product of GOSAT-2 (i.e., the L4 products) (Saito et al. 2019). NISMON employs the 4D-Var method, by which forward and backward transport simulations are performed iteratively to optimize fluxes (Niwa et al. 2017a(Niwa et al. , 2017b. The optimizer is based on a quasi-Newtonian method with the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm (Fujii 2005). The L4 products of GOSAT-2 are derived from NISMON analyses with a horizontal resolution of ca. 223 km, which is the mean interval of the "glevel-5" NICAM grid (Fig. 21). It is noteworthy that the NICAM grid shapes are hexagons or pentagons, which are generated by dividing an icosahedron recursively (Satoh et al. 2014). The temporal resolution of optimized fluxes is set as monthly. As elaborated in Niwa et al. (2017b), surface fluxes are optimized with prescribed spatial correlations in the prior errors. In the experiment of reported by Niwa et al. (2017b) with ground-based observations, the error correlation scales are set as 500 km for terrestrial areas and as 1000 km for oceans. These error correlations act as a smoother and improve flux estimates. That would be true even for the case with satellite data. Therefore, although fluxes are optimized at each model Fig. 21 Model grid of the "glevel-5" NICAM. The mean grid interval is ca. 223 km (Satoh et al. 2014) grid, the practical resolution of the fluxes that the satellite data constrain would be several times as large as the model grid. This important point must be recognized before analyzing the L4 products.
Estimated surface fluxes by inversion have theoretical errors named "posterior errors, " which are derived from diagonal elements of the posterior error covariance matrix. Because that matrix size is n × n, where n denotes the number of estimated fluxes (i.e., the number of model grids multiplied by the temporal resolution), it is not a trivial task to obtain posterior errors. However, NISMON uses a state-of-the-art method developed by Niwa and Fujii (2020) to estimate posterior errors. By investigating the degree to which posterior errors are reduced from prior errors, one can evaluate the strengths of observational constraints, an example of which is presented in Fig. 22. Posterior errors in Fig. 22 are derived using the same settings as those reported by Niwa et al. (2017b), but the ground-based observations are replaced by the GOSAT TANSO-FTS SWIR data for 2010. Because they depend strongly on the settings of prior errors and model-observation mismatch errors, which are key tuning parameters of an inversion, the values themselves in Fig. 22 might not be valid for a practical satellite inversion. Nevertheless, spatiotemporal patterns give valuable implications. As shown in Fig. 22, the satellite constraints are marginal in winter high-latitude areas. This is true because satellite retrievals are not available where sunlight does not exist. In addition, prior errors of terrestrial biosphere fluxes are often set as small because of the low biosphere activity in winter, resulting in weak observational constraints. In summer mid-latitudes to high latitudes and tropical areas, the satellite constraints are considerably strong. This fact is especially important for Fig. 22 Error reduction ratios of CO 2 flux estimates for each month. They are derived by the same settings of Niwa et al. (2017b), but with the GOSAT-1 SWIR data for 2010. The error reduction ratio is defined as σpos−σpri σpri × 100 , where σ pos and σ pri , respectively, represent the posterior and prior errors. This represents observational constraints. Here, the observation model mismatch uncertainty is assumed as 2 ppm for individual SWIR data Page 26 of 35 Imasu et al. Progress in Earth and Planetary Science 2023, 10(1):33 Siberia and tropical areas, where ground-based observation networks are extremely sparse. This picture would be true also for GOSAT-2.

Regional targets
Monitoring CO 2 emissions from large-scale CO 2 emission sources such as electric power plants and urban areas with heavy traffic is an important task. To this end, we have conducted specific observations by GOSAT targeting urban areas and power plants, particularly, intensive observations of the Kanto Plain, including Tokyo started in 2010. Similarly, we are conducting target observations with GOSAT-2 over major megacities around the world. An example of the observed data is presented in Fig. 23. These observation data are combined with the ground-based observation data, after which they are used as input data to regional transport models. The CO 2 emission intensity was estimated using inverse analysis.
By improving the SNR and by reducing the radiance biases, the estimation accuracy of the boundary layer (mixed layer) altitude obtained from the vertical distribution analysis of the temperature using the TIR band is improved. Using this information, it is expected that the accuracy of the method used for estimating the CO 2 concentration near the surface layer will be improved through the combined analysis of the columnar CO 2 concentration retrieved from SWIR band data and the upper air concentration obtained from the TIR band (Imasu 2019). This method is particularly effective for urban areas where surface temperatures are high. For such situations, the thermal contrast, which is defined by the difference between the atmosphere and the surface temperature, is high. For that reason, the resolution of the vertical distribution of air temperature in the lower atmosphere is expected to be improved.
Another major improvement of TANSO-FTS-2 is that the wavelength range of SWIR band 3 for carbon monoxide (CO) measurement has been expanded. Because CO is a substance of combustion origin, it is expected that information related to differences between CO 2 of combustion origin and CO 2 of vegetation respiratory origin in urban areas can be obtained (Takekawa et al. 2004).
Hotspots, in a climatological sense, are large carbon stores and large carbon sinks. These are areas that are susceptible to climate change and land use changes. They can be major future sources of greenhouse gas emissions. Pre-identifying and then monitoring these hotspots continuously will support early detection of signs of carbon cycle change attributable to climate change.
Specific areas include land areas where global warming effects are expected to be greater at high latitudes, especially areas where CO 2 absorption increases with extension of the growing season, areas where CO 2 absorption decreases because of dryness during summer, areas where CO 2 and CH 4 emissions increase because of large-scale forest fires, areas where albedo is expected to change because of early snowmelt, areas where frozen soil melts and the accompanying increase in CO 2 and CH 4 emissions from the soil occurs, and areas where CH 4 hydrates exist.

Solar-induced chlorophyll fluorescence (SIF)
Photosynthetic CO 2 assimilation by plant vegetation, known as gross primary production (GPP), is the largest carbon flux on the Earth surface; it is also an important driver of global climate (Global Carbon Project 2017). Its rate is extremely sensitive to environmental variation. For that reason, it varies spatiotemporally to a great degree. To date, many efforts have been undertaken using various techniques, including satellite remote sensing, to monitor the GPP of various forms of vegetation. Earlier satellite remote sensing studies have used reflectance spectral indices (e.g., normalized difference vegetation index, NDVI) that estimate plant status related to structural or pigment properties such as the leaf area index (White et al. 1997) and leaf chlorophyll content (Daughtry et al. 2000). Although these data are useful to predict the potential GPP (Running et al. 2004), they are insufficient for assessment of physiological changes of the photosynthetic system that often occur without changes in leaf area or chlorophyll content. Therefore, such vegetation indices are insufficient to assess spatiotemporal variations in GPP. In recent years, complementing such reflectance-based indices, solar-induced chlorophyll fluorescence (SIF) has been used to estimate GPP and plant stress (e.g., Parazoo et al. 2013Parazoo et al. , 2014Liu et al. 2017). Chlorophyll fluorescence is a radiation in the red and near-infrared regions with two peaks at approximately 690 and 740 nm emitted by chlorophyll-a molecules. A large part of the light energy absorbed by chlorophyll is used for photochemistry. Therefore, carbon assimilation and photorespiration, to the extent that it remains, are lost as heat and fluorescence (Baker 2008;Porcar-Castell et al. 2014). Chlorophyll fluorescence yield changes depending on the status of photochemistry and heat dissipation. In fact, it provides various information related to the photosynthetic machinery (Hikosaka and Tsujimoto 2021). However, detecting SIF in the field has been challenging (Meroni et al. 2009) because the magnitude of SIF is much lower than the reflected radiation from leaves. At a leaf scale, pulse amplitude-modulated (PAM) fluorometers have been able to detect chlorophyll fluorescence. They have been used in many physiological and ecological studies. In remote sensing, including satellite observations, SIF can be retrieved using high-spectral resolution measurements using the Fraunhofer depth method (Meroni et al. 2009).
The first SIF retrieval from satellite observations on a global scale was conducted with TANSO-FTS onboard GOSAT (Frankenberg et al. 2011cJoiner et al. 2011;Guanter et al. 2012). Band 1 of TANSO-FTS covers the overlapping part of SIF and Fraunhofer lines in the near-infrared region with high spectral resolution. Subsequently, SIF has been retrieved from the data measured by satellites aiming at measuring GHGs concentration, such as SCIAMACHY onboard ENVISAT (Joiner et al. 2012), GOME-2 of Op-A and B (Joiner et al. 2013), OCO-2 (Sun et al. 2017), TROPOMI of Sentinel 5P (Guanter et al. 2015;Köehler et al. 2018), and TanSat (Du et al. 2018). These satellite SIF data have been used for mutual validation (e.g., Oshio et al. 2019).
TANSO-FTS-2 onboard GOSAT-2 is also available to retrieve SIF. Those SIF data are provided as one standard product. GOSAT-2 is expected to provide SIF data that have higher validity than those of GOSAT because of the intelligent pointing system. Generally speaking, satellites are unable to observe SIF accurately in cloudcovered areas. The intelligent pointing system is expected to increase the number of clear-sky observations and to enhance our understanding of climate change effects on terrestrial ecosystems, especially in tropical rainforests, which are often cloud-covered. Figure 24 portrays a global map of GOSAT-2 SIF of August 2019. Actually, SIF was greatest in areas with tropical vegetation such as the Amazon, Central Africa, and Southeast Asia, which are known to have the highest GPP on Earth. By contrast, it was lowest in desert areas such as the Sahara and Arabia, and the Arctic and Antarctic regions, where vegetation is scarce. Results demonstrated that SIF is correlated linearly with data-driven GPP estimates at global and/or annual scales (Frankenberg et al. 2011c;Guanter et al. 2012;Sun et al. 2017). The correlations are rationalized to the light use efficiency (LUE) model for GPP by Monteith (1972). Actually, GPP is regarded as the product of proportional to incoming photosynthetically active radiation (PAR), the fraction absorbed PAR (fAPAR) by vegetation, and the LUE for GPP. Similarly, SIF is conceptualized as the product of PAR, fAPAR, LUE for SIF. These equations can be combined as simple correlations between GPP and SIF Guanter et al. 2014;Sun et al. 2017). This approach has been used for various ecosystems such as those of the Amazon rainforest (Parazoo et al. 2013) and the boreal forests of Eurasia and North America (Walther et al. 2016). However, the relation between GPP and SIF differs depending on the season (Thum et al. 2017), vegetation type, and canopy structure Migliavacca et al. 2017). Another approach is using the radiative transfer models, which simulate underlying mechanisms such as the biophysiological processes and scattering and reabsorption effects on SIF signals within a canopy. The SCOPE model, developed by van der Tol et al. (2009), is a 1-D radiative transfer, photosynthesis, and energy balance model. SCOPE has contributed to assessment of GPP from SIF Zhang et al. 2014;Thum et al. 2017).
Using the SCOPE-like model and carbon-cycle models, GOSAT-2's SIF would enable us to improve the accuracy of the L4 product of GOSAT-2, CO 2 flux. Accurate estimation of GPP helps elucidate the consequences between meteorological conditions and CO 2 assimilation functions of vegetation, and to predict future trends. Additionally, used in conjunction with XCO 2 data, it might contribute to estimation of anthropogenic CO 2 emissions.

Air pollution monitoring and related climate studies
One main objective of the GOSAT-2 satellite mission is air pollution monitoring. As reported by the World Health Organization (WHO 2016) and the Organization for Economic Co-operation and Development (OECD 2020), global health effects of air pollution are as great as annual excess mortality of 3-4 million people. Furthermore, the IPCC Special Report on Global Warming of 1.5° and the Sixth Assessment Report (AR6) assessment report have concluded that marked reduction of the Short-Lived Climate Forcer (SLCF) along with CO 2 reduction is necessary to attain the 1.5° reduction target for global warming mitigation (IPCC 2018(IPCC , 2021. Recent studies indicate that a SLCF reduction scenario must be designed carefully by the full analysis of complex negative and positive feedback processes among the constituents (Stohl et al. 2015;Nakajima et al. 2020aNakajima et al. , 2020b. A small but significant surface temperature decrease by the total BC reduction, for example, is explained by counterwarming effects of increases in the upper clouds because of atmospheric stability change (Takemura and Suzuki 2019).
Earth observation by the GOSAT-2 satellite contributes to air pollution monitoring and to collecting useful knowledge for the SLCF climate impacts with their unique simultaneous measurements of CH 4 , CO, and aerosols with the TANSO-FTS2 sensor with a 9.7 km FOV and the TANSO-CAI-2 imager with two UV bands of 340 nm and 380 nm and with a 500 m FOV. The fine FOV of TANSO-CAI-2 is especially useful for pollution monitoring over heterogeneous land surfaces such as city areas and coastal zones when data are analyzed using the MWMP algorithm (Hashimoto and Nakajima 2017;Shi and Nakajima 2018;Shi et al. 2019). Figure 25 shows that the two UV bands can detect smoke aerosols clearly in the case of an Australian forest fire (JAXA 2020; Nakajima et al. 2021). Estimation of the BC fraction is promising, as reported by Shi et al. (2020) for the Australian forest fire case, and for the Indian case (Gogoi et al. 2021(Gogoi et al. , 2022. Recent reports of close correlation of CO 2 concentration with short-lived atmospheric materials, i.e., NO 2 and aerosols, in model simulations and airborne and satellite observations around large-scale emission sources (e.g., Goldberg et al. 2019;Fujinawa et al. 2021) indicate that the combination of GOSAT-series satellites with other satellites for short-lived material observation is a promising method for providing wide-area information that is useful for emission inventory construction and the United Nations Framework Convention on Climate Change (UNFCCC) global stock taking routine. Itahashi et al. (2019) demonstrated an update of REAS NO x inventory through assimilation of OMI-observed NO 2 data with short turn-around time of less than several years. This benefit is expected to be increased in the next-phase global monitoring system by Geostationary Environment Monitoring Spectrometer (GEMS), Tropospheric Emissions: Monitoring Pollution (TEMPO), and Sentinel-4 geostationary satellites (Kim et al. 2020). Therefore, designing a satellite that is able to monitor NOx and GHGs simultaneously is expected to be important for future work.

Conclusions
Compared to the first GOSAT sensors, the GOSAT-2 equipped sensors have demonstrated improved performance in most respects. The quality and quantity of observation data are expected to be improved accordingly. The signal-to-noise ratio (SNR) is improved in both the SWIR and TIR bands of TANSO-FTS-2, which is the main sensor of GOSAT-2. This improvement will ultimately enhance the accuracy of greenhouse gas concentration analysis. In addition, because of the improved SNR in the SWIR band, the northern limit at which data are obtainable in high-latitude regions of the Northern Hemisphere in winter, where observation data were heretofore unavailable because of weak signal strength, has moved to higher latitudes. As more and better data are obtained, progress in carbon cycle research for highlatitude regions is anticipated. Moreover, it is estimated that the improvement of the SNR in the TIR band can be considerable. Particularly, the resolution of the vertical concentration distribution of CO 2 and CH 4 has been improved drastically. Using the CO 2 column concentration obtained from the SWIR band and the upper air concentration of CO 2 obtained from the TIR band in combination, the retrieved CO 2 concentrations near the ground surface, which are important for estimating the amount of CO 2 emission and uptake, are obtainable. This benefit of GOSAT/GOSAT-2 is not provided by satellites equipped with other similar sensors.
The first function introduced into GOSAT-2 which has not been implemented in GOSAT is an intelligent pointing mechanism: a cloud area avoidance function using the in-field camera of TANSO-FTS-2. This function can be expected to increase the amounts of observation data globally and to improve the estimation accuracy of CO 2 emissions and uptake intensity. The effect is expected to be strong, especially in the tropics, because the most common cloud type is cumulus there. The intelligent pointing system can avoid clouds effectively.
Another major change of TANSO-FTS-2 is that the wavelength range of Band 3 of SWIR has been expanded for measuring carbon monoxide (CO). Because CO is a substance originating from combustion, it is used to evaluate the effects of human activities in urban areas and biomass burning in the field. Particularly, black carbontype aerosols can be measured by TANSO-CAI-2 for biomass burning as well as CO 2 and CO by TANSO-FTS-2. The partitioning ratio of these three substances might provide information that characterizes the fire, such as combustion completeness. In fact, GOSAT-2 is the only satellite that can simultaneously measure Short-Lived Climate Forces (SLCF) such as CH 4 , BC, and groundlevel ozone together with CO 2 .
The ultimate goal of most greenhouse gas observation satellite projects, including GOSAT and GOSAT-2, is elucidation of anthropogenic greenhouse gas emissions and their effects on the carbon cycle and climate. However, to extract only the effects of anthropogenic origin from the information obtained using these satellite-based sensors, it is necessary to evaluate and eliminate the effects of natural origin which are included therein. Eliminating these influences is mainly a question of how to achieve accurate evaluation of CO 2 emissions (respiration) and absorption (photosynthesis) by terrestrial ecosystems. One method is cooperation with satellites launched to evaluate the dynamics of terrestrial ecosystems. By combining related to various activities of terrestrial ecosystems such as the vegetation index and heat dissipation (photochemical reflection index) obtained from these satellites with SIF information obtained from GOSAT-2, the physiological and ecological activity of vegetation can be better evaluated. In this way, the combined use of various satellite data makes it possible to acquire more advanced information related to the carbon cycle.
Future research undertaken to elucidate the global carbon cycle including greenhouse gases and research to evaluate the global environment are expected to become more sophisticated and increasingly complicated. Therefore, achieving sufficient scientific results with only a single satellite project is difficult. The only solution is addressing these issues through collaboration with numerous satellite projects and extensive international cooperation.