GOSAT CH 4 Vertical Proﬁles over the Indian Subcontinent: Effect of a Priori and Averaging Kernels for Climate Applications

: We examined methane (CH 4 ) variability over different regions of India and the surrounding oceans derived from thermal infrared (TIR) band observations (TIR CH 4 ) by the Thermal and Near-infrared Sensor for carbon Observation—Fourier Transform Spectrometer (TANSO-FTS) onboard the Greenhouse gases Observation SATellite (GOSAT) for the period 2009–2014. This study attempts to understand the sensitivity of the vertical proﬁle retrievals at different layers of the troposphere and lower stratosphere, on the basis of the averaging kernel (AK) functions and a priori assumptions, as applied to the simulated concentrations by the MIROC4.0-based Atmospheric Chemistry-Transport Model (MIROC4-ACTM). We stress that this is of particular importance when the satellite-derived products are analyzed using different ACTMs other than those used as retrieved a priori. A comparison of modeled and retrieved CH 4 vertical proﬁles shows that the GOSAT/TANSO-FTS TIR instrument has sufﬁcient sensitivity to provide critical information about the transport of CH 4 from the top of the boundary layer to the upper troposphere. The mean mismatch between TIR CH 4 and model is within 50 ppb, except for the altitude range above 150 hPa, where the sensitivity of TIR CH 4 observations becomes very low. Convolved model proﬁles with TIR CH 4 AK reduces the mismatch to less than the retrieval uncertainty. Distinct seasonal variations of CH 4 have been observed near the atmospheric boundary layer (800 hPa), free troposphere (500 hPa), and upper troposphere (300 hPa) over the northern and southern regions of India, corresponding to the southwest monsoon (July–September) and post-monsoon (October–December) seasons. Analysis of the transport and emission contributions to CH 4 suggests that the CH 4 seasonal cycle over the Indian subcontinent is governed by both the heterogeneous distributions of surface emissions and the inﬂuence of the global monsoon divergent wind circulations. The major contrast between monsoon, and pre- and post-monsoon proﬁles of CH 4 over Indian regions are noticed near the boundary layer heights, which is mainly caused by seasonal change in local emission strength with a peak during summer due to increased emissions from the paddy ﬁelds and wetlands. A strong difference between seasons in the middle and upper troposphere is caused by convective transport of the emission signals from the surface and redistribution in the monsoon anticyclone of upper troposphere. TIR CH 4 observations provide additional information on CH 4 in the region compared to what is known from in situ data and total-column (XCH 4 ) measurements. Based on two emission sensitivity simulations compared to TIR CH 4 observations, we suggest that the emissions of CH 4 from the India region were 51.2 ± 4.6 Tg year − 1 during the period 2009–2014. Our results suggest that improvements in the a priori proﬁle shape in the upper troposphere and lower stratosphere (UT/LS) region would help better interpretation of CH 4 cycling in the earth’s environment.


Introduction
The South Asia region, consisting of India, Pakistan, Bangladesh, Nepal, Bhutan, and Sri Lanka, play an import part in the global methane (CH 4 ) budget, because the regional total emissions contributed 8% (about 500 Tg) of the CH 4 global total emissions during the 2000s [1]. The recent economic growth of India has led to a significant increase in industrial emissions [2][3][4], especially in the Indo-Gangetic Plain (IGP), encompassing northern regions of the Indian subcontinent. The IGP is one of the most densely populated regions with strong emissions of various tracers due to human activity [5]. However, there are little or no long-term measurements of CH 4 and the other greenhouse gases to evaluate the inventory emissions [6][7][8][9].
The Asian Summer Monsoon Anticyclone (ASMA) is a dominant circulation pattern in the Upper Troposphere and Lower Stratosphere (UT/LS) in the northern hemisphere in summer, which extends from Southeast Asia to the Middle East [10,11]. The Asian monsoons may be classified into a few sub-systems, such as the South Asian Monsoon, which affects the Indian subcontinent and surrounding regions. The monsoon is associated with persistent strong convection over India and the Bay of Bengal, elevated surface heating over the Tibetan Plateau, and orographic uplifting at the southern/south-western slopes of the Himalayas, which contributes to overall ascension of boundary layer air to the upper troposphere (up to~200 hPa) [12]. The deep convection and associated circulation patterns of the monsoon provides an important pathway for polluted boundary layer air to reach UT/LS [13,14]. The atmospheric compounds can then be advected horizontally to other regions, or further uplifted in to the stratosphere [1,8,15,16]. Due to the influence of deep convection and long-range transport, the chemical tracers such as CH 4 , CO, and ozone often show extreme values [17,18]. Global observations from satellite instruments can complement and extend the information available from the in situ observations at sparse surface networks and aircraft campaigns to improve our knowledge of the processes controlling the emission and distribution of CH 4 , and to monitor its variability on different scales. The shortwave infrared (SWIR) and thermal infrared (TIR) bands are available for measurements of CH 4 from space. Observations in the SWIR, such as those from the Thermal Additionally, Near-infrared Sensor for carbon Observation-Fourier Transform Spectrometer (TANSO-FTS) onboard the Greenhouse gases Observation SATellite (GOSAT) [19][20][21], provide information on column-averaged CH 4 (XCH 4 ), in cloud-free conditions during the daytime satellite overpass. On the other hand, TIR CH 4 observations provide much greater geographical and temporal coverage, and more importantly, the measurements of vertical profiles would allow a better understanding of the CH 4 cycle over a region. The sensitivity of TIR CH 4 observations is stronger in the mid to upper troposphere and relatively low near the surface, because its spectral signatures depend on thermal contrast between the atmosphere and surface [22,23].
The heterogeneity in surface emissions variations of XCH 4 governed by complex atmospheric transport mechanisms during the southwest monsoon season in July-September and northeast monsoon season in October-December is observed by GOSAT. The work by Chandra et al. [8] has highlighted the difficulty in interpreting the emissions from the surface by columnar CH 4 measurements from SWIR spectra, without using an atmospheric chemistry-transport model. At the same time, Ricaud et al. [24] investigated the impact of the Asian Monsoon Anticyclone to the distribution of the mid-to-upper CH 4 in the Mediterranean Basin using a wide variety of datasets including GOSAT/TANSO-FTS TIR observations.
Critical data for interpretation of satellite retrievals of trace gas are the a priori information on the shape of the vertical profile and the averaging kernel (AK) functions. Both should be implemented to the simulated profile to resample the simulated profile in a manner that mimics the satellite vertical sensitivity [25]. A priori profiles are derived from a chemical transport model or a representation which has different physical processes, emissions, or spatial resolutions. As a result, comparisons between satellite-retrieved and model-simulated column abundances are greatly influenced by the a priori profile shape [26].
This study analyzes the vertical distributions of CH 4 over the Asian monsoon region. We used CH 4 mixing ratios observed from GOSAT/TANSO-FTS TIR (hereafter referred to as "GOSAT-TIR") Version 1 (Level 2 V1) and simulated by MIROC4.0-based Atmospheric Chemistry-Transport Model (MIROC4-ACTM) for the period from 23 April 2009, until 24 May 2014. We discuss challenges in using GOSAT-TIR CH 4 (hereafter referred to as "TIR CH 4 ") data to determine relative contributions of surface emissions and transport in the formation of CH 4 seasonal cycles over different parts of India and the surrounding oceans. The paper is structured as follows. In Section 2, we briefly describe the spaceborne instrument GOSAT-TIR and vertical profile retrievals of CH 4 , the MIROC4-ACTM simulation setup, the study domain, and data processing. The meteorology and climatology of CH 4 inferred from the different datasets over the study domain, variability of CH 4 vertical profiles, and the impact of the Asian Monsoon Anticyclone to the distribution of the tropospheric CH 4 are discussed in Section 3. Major conclusions are given in Section 4.

Study Domain
This work follows the setup described by Chandra et al. [8] and uses 10 regions ( Figure 1a  model-simulated column abundances are greatly influenced by the a priori profile shape [26]. This study analyzes the vertical distributions of CH4 over the Asian monsoon region. We used CH4 mixing ratios observed from GOSAT/TANSO-FTS TIR (hereafter referred to as "GOSAT-TIR") Version 1 (Level 2 V1) and simulated by MIROC4.0-based Atmospheric Chemistry-Transport Model (MIROC4-ACTM) for the period from 23 April 2009, until 24 May 2014. We discuss challenges in using GOSAT-TIR CH4 (hereafter referred to as "TIR CH4") data to determine relative contributions of surface emissions and transport in the formation of CH4 seasonal cycles over different parts of India and the surrounding oceans. The paper is structured as follows. In Section 2, we briefly describe the spaceborne instrument GOSAT-TIR and vertical profile retrievals of CH4, the MIROC4-ACTM simulation setup, the study domain, and data processing. The meteorology and climatology of CH4 inferred from the different datasets over the study domain, variability of CH4 vertical profiles, and the impact of the Asian Monsoon Anticyclone to the distribution of the tropospheric CH4 are discussed in Section 3. Major conclusions are given in Section 4.

Study Domain
This work follows the setup described by Chandra et al. [8] and uses 10 regions  The map of the regional divisions for the analysis (a), depicted by unicolor boxes with the name abbreviations defined in Section 2.1. The black lines in-land show major rivers in the plot domain; the human settlements along the rivers are major sources of CH4. The TIR CH4 a priori with 1-σ SD uncertainty (red line with error bars) and TIR CH4 profile with the retrieval error (blue line with shaded area) are shown in panels (b,d,f), and the 22 AK functions corresponding to the 22 retrieval layers are in panels (c,e,g). Data in panels (b-g) are shown for the regions Northeast India, Southern India, and the Bay of Bengal, averaged over the period July-September 2011. Figure 1. The map of the regional divisions for the analysis (a), depicted by unicolor boxes with the name abbreviations defined in Section 2.1. The black lines in-land show major rivers in the plot domain; the human settlements along the rivers are major sources of CH 4 . The TIR CH 4 a priori with 1-σ SD uncertainty (red line with error bars) and TIR CH 4 profile with the retrieval error (blue line with shaded area) are shown in panels (b,d,f), and the 22 AK functions corresponding to the 22 retrieval layers are in panels (c,e,g). Data in panels (b-g) are shown for the regions Northeast India, Southern India, and the Bay of Bengal, averaged over the period July-September 2011.

GOSAT-TIR Retrievals
GOSAT was the first satellite dedicated to global observations of greenhouse gases CO 2 and CH 4 from space. Since the launch on 23 January 2009, observations are made on a 666 km sun-synchronous orbit with a 3-day revisit cycle, and 12-day operation cycle, with overpasses twice a day at local noon 13:00 ± 15 min and near-midnight [27]. During daytime, GOSAT/TANSO-FTS can simultaneously observe CH 4 column-averaged dry-air mole fractions and CH 4 profiles in the same field of view (a nadir footprint diameter of 10.5 km) using SWIR and TIR spectral bands, respectively [28,29]. Therefore, daytime observations are screened more thoroughly by TANCO-Cloud and Aerosol Imager (CAI). At night, observations and filtering are carried out only by the TIR algorithms. Therefore, the number of scenes of daytime and night-time observations can differ significantly, especially during the monsoon period. However, the average difference between the retrieved CH 4 is usually within ±20 ppb when data within one day are compared. We leave a more detailed analysis of the daily variation of CH 4 outside the scope of this work.
The retrieval algorithm for the GOSAT-TIR V1 CH 4 adopted a nonlinear maximum a posteriori (MAP) method with linear mapping [22]. A priori covariance matrix for CH 4 in the V1 CH 4 retrieval is set to be a diagonal matrix with vertically fixed diagonal elements with a standard deviation of 4%. The first retrieval version of the TIR CH 4 product (Level 2 V00.01) and its validation analysis showed that the TIR XCH 4 profiles were within 0.5% agreement with the aircraft XCH 4 values over the tropical ocean [30]. Work by Holl et al. [31] compared CH 4 data from the Canadian Space Agency's Atmospheric Chemistry Experiment FTS (ACE-FTS), ground-based FTS, and the current released version of GOSAT-TIR (Level 2 V1) in the Canadian high Arctic, although TIR CH 4 measurement information content is too low for true profile retrieval because of the low thermal contrast and the low signal-to-noise ratio over the Arctic. Zou et al. [32] performed global comparisons with the Atmospheric Infrared Sounder (AIRS) retrievals for the levels of 300-600 hPa, where both AIRS and GOSAT-TIR V1 CH 4 have peak sensitivities. Mean mismatch in CH 4 (GOSAT-AIRS) were 10.3-31.8 and −16.2 ± 25.7 ppbv for the levels of 300 and 600 hPa, respectively. Comparison of the XCH 4 shows that TIR CH 4 agrees with AIRS to within 1% in the mid-latitude regions of the Southern Hemisphere and in the tropics, and tõ 1-2% in the mid to high latitudes [32]. Olsen et al. [33] performed global comparisons of GOSAT-TIR V1 CH 4 vertical profiles to ACE-FTS (version 3.5) and the European Space Agency's Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), as well as 16 ground stations with the Network for the Detection of Atmospheric Composition Change (NDACC). In the overlapping altitude ranges of the three satellite data products there is a small, but consistent, positive bias of around 20 ppbv, or 1% in TIR CH 4 data. In the upper troposphere, no bias between TIR CH 4 and NDACC was found.
In a more recent comparison, the average bias in CH 4 profile retrieved from GOSAT-TIR spectra with a spectral correction scheme was less than 2% over the full altitude range, when compared with data from the Monitoring Atmospheric Composition and Climate (MACC) scaled to the total column measurements of the Total Carbon Column Observing Network (TCCON) [23].
Global comparisons are also conducted based on HIPPO, CARIBIC, JMA, and CON-TRAIL/ASE aircraft observations (Saitoh, N. et al.: Intensive validation analysis of GOSAT/ TANSO-FTS thermal infrared CH 4 data (version 1) based on aircraft observations, in preparation). In low latitudes corresponding to the India location, compared datasets include CARIBIC profiles over MAA (Chennai, India), BOG (El Dorado, Colombia), and CCS (Venezuela) airports and CONTRAIL/ASE over GUAM (U.S.A.) airport. The validations show that TIR V1 CH 4 and aircraft CH 4 profiles agreed with each other within 10-15 ppb and there was no evident seasonal dependence in the CH 4 differences.

TIR CH 4 Profile Properties
The observation by GOSAT-TIR band enables analysis of the vertical structure of atmospheric CH 4 . This band has relatively high spectral resolution of~0.2 cm −1 and Remote Sens. 2021, 13, 1677 5 of 25 provides CH 4 vertical profiles in 22 layers from the surface to 0.1 hPa. The degrees of freedom (DOF) of signal for TIR CH 4 observations (V1 algorithm version) is around 1 over low-latitude parts of India. DOF means the number of purely (mathematically) independent pieces of information. However, the real atmospheric layers correlate with each other, so even if its DOF is close to 1, TIR has the ability to derive new knowledge about CH 4 concentrations. Figure 1 shows a typical TIR CH 4 profile with retrieval errors and a priori profile and its corresponding AK profiles for the southern and northeast regions of India, and over the Bay of Bengal during the monsoon season. The a priori profiles (Figure 1b,d,f) used in the CH 4 retrieval were provided by the National Institute for Environmental Studies (NIES) transport model using 47 sigma levels. The stratospheric part of the NIES model was adjusted to the observed age of air for CO 2 and long-term satellite observations from HALOE for CH 4 [34].
Analysis of the AK profiles of the TIR CH 4 products (Figure 1c,e,g) show that the GOSAT-TIR measurement signal is beyond the retrieval error at pressure levels of 500-150 hPa with the sensitivity maximum at the levels of 200-400 hPa decreasing towards the surface. The spectra sensitivity does not change significantly between the different parts of the study domain but becomes low in the lower troposphere (below 800 hPa), emphasizing the major limitation of TIR CH 4 measurements.

MIROC4-ACTM Simulations
The measurements were compared to results of simulations by the MIROC4-ACTM simulation [35,36]. The MIROC4-ACTM runs at a horizontal resolution of T42 spectral truncations (~2.8 ×~2.8) with 67 sigma-pressure vertical levels. The MIROC4-ACTM simulated horizontal winds (U and V) and temperature (T) are nudged to the Japan Meteorological Agency reanalysis fields (JRA-55) at all the vertical levels [37]. The model uses an optimal OH field based on a scaled version of the seasonally varying OH field [38,39].
Two simulations were performed using combinations of inverted fluxes at monthly intervals, based on two a priori emission scenarios [40]:
The anthropogenic emissions for an individual sector were taken from EDGAR (the Emission Database for Global Atmospheric Research, version 4.3.2 inventory [4]) on an annual basis for the period 1970-2012, and the global total emissions were extrapolated to 2017 using the 2011-2012 rate of change. GFED (Global Fire Database version 4s [43]) provides interannually varying monthly emissions for biomass burning. The other emissions (including ocean, termites, mud volcano, etc.) were taken from the TransCom-CH 4 inter-comparison experiment [44]. While all other emissions were common, the VISIT (Vegetation Integrated Simulator of Trace gases [45,46]) provided two distinct sets of wetland and rice emissions on a monthly basis using two different schemes developed by [41,42], which are referred to as "ACTM Cao " and "ACTM WH ", respectively ( Figure 2). The main difference between these approaches in most cases can be explained by the combined effect of changes in soil temperature and the position of the water table. In general, the emissions WH scheme are about 5-10% larger than those by the Cao scheme, except for the WIGP, EIGP, and NEI regions of India and Bangladesh where the maximum difference reaches 20-40% ( Figure 2). Moreover, there are small number of hotspots in Southeast Asia (e.g., Mekong River Delta region), in the case of Cao scheme compared to WH scheme.

Data Processing
The MIROC4-ACTM simulations were sampled and collocated in space and time to the TIR CH 4 observation points. The criteria for the collocation were the nearest model grid cell in space, and the nearest hour in time. For vertical profile comparison, the MIROC4-ACTM data were interpolated on the retrieval pressure levels of the TIR CH 4 product, i.e., from 67 model levels to 22 retrieval levels.

Data Processing
The MIROC4-ACTM simulations were sampled and collocated in space and time to the TIR CH4 observation points. The criteria for the collocation were the nearest model grid cell in space, and the nearest hour in time. For vertical profile comparison, the MI-ROC4-ACTM data were interpolated on the retrieval pressure levels of the TIR CH4 product, i.e., from 67 model levels to 22 retrieval levels.

AK Functions and the Retrieval Sensitivity
AK functions are defined to provide a simple characterization of the relationship between the retrieval and the true state. The retrieval sensitivity can be obtained from the sum of the columns of the AK matrix, which is also referred to as "the area of the averaging kernel" [25]. Along with "raw" model simulation results (ACTM , ), we analyzed (ACTM , ) profiles convolved with retrieval a priori and the TIR CH4 AK matrix using the following vector equation [25,30]: Here, A is an AK matrix, Xa priori represents a vector of an a priori vertical profile, and , and , are vectors of "raw" and convolved model simulated profiles, respectively.

The Prophet Analysis and Forecasting Model
TIR CH4 observations show large temporal and spatial heterogeneity, which complicate the derivation of seasonal cycle variation. The Prophet is a novel time-series analysis and forecasting model [47], which performs the smoothing of time-series data based on a

AK Functions and the Retrieval Sensitivity
AK functions are defined to provide a simple characterization of the relationship between the retrieval and the true state. The retrieval sensitivity can be obtained from the sum of the columns of the AK matrix, which is also referred to as "the area of the averaging kernel" [25]. Along with "raw" model simulation results (ACTM Cao,WH ), we analyzed (ACTM AK Cao,WH ) profiles convolved with retrieval a priori and the TIR CH 4 AK matrix using the following vector equation [25,30]: Here, A is an AK matrix, X a priori represents a vector of an a priori vertical profile, and X ACTM Cao,WH and X ACTM AK Cao,WH are vectors of "raw" and convolved model simulated profiles, respectively.

The Prophet Analysis and Forecasting Model
TIR CH 4 observations show large temporal and spatial heterogeneity, which complicate the derivation of seasonal cycle variation. The Prophet is a novel time-series analysis and forecasting model [47], which performs the smoothing of time-series data based on a generalized additive model with three main properties: trend, seasonality, and holidays. Compared to traditional exponential smoothing, Prophet can handle temporal patterns with multiple periods and has no requirements regarding the regularity of measurement spacing. The model has a robust performance in the presence of missing data and trend shifts, and typically handles outliers well while working with time-series that have several seasons of historical data with strong seasonal patterns. Prophet managed to use all data points for the study period, thereby increasing accuracy and reducing the sensitivity to random outliers [48].

Prevailing Atmospheric Conditions over the Indian Subcontinent
Key components of the climatology in the Indian Ocean and the surrounding areas are the annual migration of the Intertropical Convergence Zone (ITCZ) and seasonal development of monsoon winds [10,11,49,50]. In boreal spring, the ITCZ migrates northward across the Indian Ocean and reaches its northernmost position at approximately 35 • N during summer. A strong pressure gradient between the low-pressure zone over the Tibetan Plateau and a high-pressure zone over the Southern Indian Ocean generates a strong near-surface monsoonal airflow from July to September (Figure 3c1-c3). In autumn, the ITCZ retreats southward and reaches its southernmost position in January. The reversed pressure gradient during the winter months generates the moderate and dry northeast monsoon (Figure 3d1-d3). Devasthale and Fueglistaler [51] showed that a significant fraction of high opaque clouds reaches and penetrates the tropical tropopause layer during active ASMA conditions ( Figure 3c). The overall frequency of convective clouds (reaching at least 200 hPa) is higher in July and August. Most of the deep convection occurs over the Bay of Bengal and central northeast India [52], while very deep convection over the Tibetan plateau is comparatively weak, and may play only a secondary role in troposphere-to-stratosphere transport.

CH 4 over India Observed by GOSAT-TIR and Simulated by MIROC4-ACTM
We analyzed CH 4 distributions from GOSAT-TIR and MIROC4-ACTM regrided on the horizontal grid with a resolution of 3.0 • × 3.0 • and interpolated at the levels of constant pressure of 800, 500, and 300 hPa, which represent the top of the boundary layer, the free troposphere, and the upper troposphere heights of the Earth's atmosphere, respectively. CH 4 concentrations were averaged over time for four periods: pre-monsoon April-June (AMJ), monsoon July-September (JAS), post-monsoon October-December (OND), and winter season January-March (JFM) with focus on the first three periods, because no strong difference was found between JFM and AMJ. Due to a lack of TIR CH 4 data in cloudy scenes and the influence of the complex orography of the studied area, the number of points used for averaging in each grid cell varies with height over land (Figure 4a1-a3). This is especially noticeable for the northern regions of India, because a significant part of Tibet and the Himalayas are above the level of 800 hPa (Figure 4a1). Northern India also has large sources of CH 4 with different types of emission. These two factors cause large standard deviations (SD) in CH 4 (Figure 4b1-b3). For South India and the marine regions, the SD values are much lower compared to those over the land. In the middle and upper troposphere, the perturbations from the heterogeneity of the emissions are smoothed out, the density of observation points increases; therefore, the averaging errors decrease. At a height of 300 hPa, the average SD for TIR CH 4 is approximately 25 ppb. Remote Sens. 2021, 13, x FOR PEER REVIEW 8 of 26   Tibet and the Himalayas are above the level of 800 hPa (Figure 4a1). Northern India also has large sources of CH4 with different types of emission. These two factors cause large standard deviations (SD) in CH4 (Figure 4b1-b3). For South India and the marine regions, the SD values are much lower compared to those over the land. In the middle and upper troposphere, the perturbations from the heterogeneity of the emissions are smoothed out, the density of observation points increases; therefore, the averaging errors decrease. At a height of 300 hPa, the average SD for TIR CH4 is approximately 25 ppb. The density of observation points decreases with the onset of the monsoon season ( Figure S1); however, it remains sufficient to detect significant changes in CH4 concentrations even considering the relatively large SD values there. A significant increase in concentration values is noticeable, primarily in the middle and upper troposphere, which is due to the repeatedly confirmed effect of convective transport from surface sources upward. After reaching a level near the tropopause, the increased concentrations are distrib- The density of observation points decreases with the onset of the monsoon season ( Figure S1); however, it remains sufficient to detect significant changes in CH 4 concentrations even considering the relatively large SD values there. A significant increase in concentration values is noticeable, primarily in the middle and upper troposphere, which is due to the repeatedly confirmed effect of convective transport from surface sources upward. After reaching a level near the tropopause, the increased concentrations are distributed by three jets: the lateral (the cross-equatorial circulation) and transverse (flows between the arid regions of north Africa and the Near East and south Asia) monsoons, and the Walker Circulation is extended across the Pacific Ocean East [10]. The CH 4 concentration in the eastern jets is higher, because it is formed over more northern areas with larger emission. The influence of the third component (the cross-equatorial circulation) is more noticeable in the post-monsoon period ( Figure S2).

CH4 over India Observed by GOSAT-TIR and Simulated by MIROC4-ACTM
TIR CH 4 retrievals are constrained to the a priori CH 4 data (panels b1-b3 of Figures 5  and 6) especially in lower pressure levels due to the relatively low signal-to-noise ratio of the TIR spectra at CH 4 absorption bands [30,32]. Nevertheless, the TIR CH 4 product shows clear differences in CH 4 from the a priori values, even in the lower part of the atmosphere where sensitivity is weak (panels c1-c3 of Figures 5 and 6 for AMJ and JAS, and Figures  S3 and S4 for JFM and OND). These differences are mainly located in the regions with high flux uncertainty. Apparently, the a priori CH 4 simulated by NIES TM with the Emission Database for Global Atmospheric Research (EDGAR) version 32FT2000 [53] for anthropogenic CH 4 , and GISS emissions [44,54] for natural CH 4 shows large negative bias in this area. This implies that an additional signal of CH 4 concentration could be captured by the GOSAT-TIR measurements.
where sensitivity is weak (panels c1-c3 of Figures 5 and 6 for AMJ and JAS, and Figures  S3 and S4 for JFM and OND). These differences are mainly located in the regions with high flux uncertainty. Apparently, the a priori CH4 simulated by NIES TM with the Emission Database for Global Atmospheric Research (EDGAR) version 32FT2000 [53] for anthropogenic CH4, and GISS emissions [44,54] for natural CH4 shows large negative bias in this area. This implies that an additional signal of CH4 concentration could be captured by the GOSAT-TIR measurements. Figure 5. Latitude-longitude distributions of CH 4 at the levels of 800, 500, and 300 hPa (the left, middle, and right panels, respectively) observed by GOSAT-TIR for the season AMJ 2011. The top (a1-a3), middle (b1-b3) and bottom (c1-c3) rows show TIR CH 4 , TIR CH 4 a priori, and difference between the TIR CH 4 and a priori distributions, respectively.
As explained in Section 2.3, MIROC4-ACTM simulations were performed with two flux combinations reflecting different approaches for estimation of the wetland CH 4 emissions. In the pre-and post-monsoon seasons (AMJ and OND), the excess concentration due to additional emission is locked in the boundary layer (as seen from MOPITT CO) [5]; therefore, we can only detect a slight increase in concentration at the levels selected for the analysis. CH 4 simulated using both emission schemes was consistent with the GOSAT-TIR retrieval with an averaged mismatch within ±2%, the heterogeneity of which is apparently caused by transport regimes. By analogy to the CH 4 distribution from GOSAT-TIR, the increased scatter was found in modeled CH 4 over IGP, wherein the enhanced values extended up to the level of 300 hPa (see Supplementary Figures S5 and S6).
Remote Sens. 2021, 13, x FOR PEER REVIEW 11 of 26 Figure 5. Latitude-longitude distributions of CH4 at the levels of 800, 500, and 300 hPa (the left, middle, and right panels, respectively) observed by GOSAT-TIR for the season AMJ 2011. The top (a1-a3), middle (b1-b3) and bottom (c1-c3) rows show TIR CH4, TIR CH4 a priori, and difference between the TIR CH4 and a priori distributions, respectively. As explained in Section 2.3, MIROC4-ACTM simulations were performed with two flux combinations reflecting different approaches for estimation of the wetland CH4 emissions. In the pre-and post-monsoon seasons (AMJ and OND), the excess concentration due to additional emission is locked in the boundary layer (as seen from MOPITT CO) [5]; therefore, we can only detect a slight increase in concentration at the levels selected for the analysis. CH4 simulated using both emission schemes was consistent with the GOSAT-TIR retrieval with an averaged mismatch within ±2%, the heterogeneity of which is apparently caused by transport regimes. By analogy to the CH4 distribution from GOSAT-TIR, the increased scatter was found in modeled CH4 over IGP, wherein the enhanced values extended up to the level of 300 hPa (see Supplementary Figures S5 and S6). During the monsoon, the difference between emission scenarios becomes significant, as additional CH 4 -rich airmass is carried to the middle and upper atmosphere (Figure 7). The large mismatch in comparison with GOSAT-TIR (panels b1-b3 of Figure 7) emphasizes the overestimation of CH 4 emission in the WH scheme. Apparently, production of CH 4 due to wetland and rice along the Ganges river in the WH scheme of the VISIT model is excessively high. More detailed flux strength estimations are in Section 3.6.
During the monsoon, the difference between emission scenarios becomes significant, as additional CH4-rich airmass is carried to the middle and upper atmosphere (Figure 7). The large mismatch in comparison with GOSAT-TIR (panels b1-b3 of Figure 7) emphasizes the overestimation of CH4 emission in the WH scheme. Apparently, production of CH4 due to wetland and rice along the Ganges river in the WH scheme of the VISIT model is excessively high. More detailed flux strength estimations are in Section 3.6.  Figure 8 depicts seasonal mean CH4 vertical profiles observed by GOSAT-TIR and simulated by the model for pre-monsoon (April-June) and monsoon (July-September) seasons of 2011. The variation of GOSAT-TIR sensitivity is taken into account by the implementation of the a priori profiles and AK functions (Figure 1c,e,g) to the modelled datasets (ACTM and ACTM ). The variabilities shown in Figure 8 are larger than GO-SAT-TIR retrieval errors, derived here as the diagonal elements of the posteriori error covariance matrices based on the MAP method, which include random error components of the retrieval. Therefore, GOSAT-TIR and the model show good agreements (mismatch is inside 1-σ SD uncertainty) within both errors (natural variabilities and retrieval random errors). The vertical CH4 profiles have a characteristic curved shape with double peak. The first peak near the surface is associated with emissions from local sources; the second one at the level of 150-200 hPa is caused by the vertical updraft [55,56]. The vertical gradient between the near-surface and upper troposphere levels increases in the direction from south-west to north-east, reflecting the increase in CH4 surface flux intensities (Figure 2).   Figure 8 are larger than GOSAT-TIR retrieval errors, derived here as the diagonal elements of the posteriori error covariance matrices based on the MAP method, which include random error components of the retrieval. Therefore, GOSAT-TIR and the model show good agreements (mismatch is inside 1-σ SD uncertainty) within both errors (natural variabilities and retrieval random errors). The vertical CH 4 profiles have a characteristic curved shape with double peak. The first peak near the surface is associated with emissions from local sources; the second one at the level of 150-200 hPa is caused by the vertical updraft [55,56]. The vertical gradient between the near-surface and upper troposphere levels increases in the direction from south-west to north-east, reflecting the increase in CH 4 surface flux intensities (Figure 2). Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 26 Monsoons cause a powerful perturbation of concentration along the entire vertical profile up to the level of the tropopause. Two southern regions (the Arabian Sea and Southern India; Figure 8a,h) are located near the entry point of Somali Jet-atmospheric masses with low CH4 concentrations coming from the Indian Ocean [7,49]. These regions do not have significant sources of CH4; therefore, concentrations in the vertical profiles increase with height due to transport from other regions. The third southern region (Bay of Bengal; Figure 8c) has similar properties, but at the same time, it is under the influence of transport from neighboring regions (i.e., East India, EIGP), as evidenced by a large spread near the surface.

CH 4 Vertical Profiles
The use of AK takes into account the relatively low vertical resolution of satellite measurements and the change in the sensitivity of the retrieval with altitude. Convolved model profiles with TIR CH4 AK (Equation (1)) along the a priori profile smooth the model profiles to fit the GOSAT-TIR vertical resolution and reduce the mismatch at the levels Monsoons cause a powerful perturbation of concentration along the entire vertical profile up to the level of the tropopause. Two southern regions (the Arabian Sea and Southern India; Figure 8a,h) are located near the entry point of Somali Jet-atmospheric masses with low CH 4 concentrations coming from the Indian Ocean [7,49]. These regions do not have significant sources of CH 4 ; therefore, concentrations in the vertical profiles increase with height due to transport from other regions. The third southern region (Bay of Bengal; Figure 8c) has similar properties, but at the same time, it is under the influence of transport from neighboring regions (i.e., East India, EIGP), as evidenced by a large spread near the surface.
The use of AK takes into account the relatively low vertical resolution of satellite measurements and the change in the sensitivity of the retrieval with altitude. Convolved model profiles with TIR CH 4 AK (Equation (1)) along the a priori profile smooth the model profiles to fit the GOSAT-TIR vertical resolution and reduce the mismatch at the levels where the sensitivity of satellite sensors is weak, as shown in Figure 9.The use of TIR CH 4 AK functions have a significant smoothing effect, approaching the MIROC4-ACTM model profiles to a priori so much that the difference between the calculations for the Cao and WH emission scenarios becomes barely distinguishable. This is especially visible above the level of 150 hPa, where the sensitivity of GOSAT-TIR there drops sharply, and the satellite retrievals and the AK convolved model profiles strongly follow the a priori profiles. where the sensitivity of satellite sensors is weak, as shown in Figure 9.The use of TIR CH4 AK functions have a significant smoothing effect, approaching the MIROC4-ACTM model profiles to a priori so much that the difference between the calculations for the Cao and WH emission scenarios becomes barely distinguishable. This is especially visible above the level of 150 hPa, where the sensitivity of GOSAT-TIR there drops sharply, and the satellite retrievals and the AK convolved model profiles strongly follow the a priori profiles.

CH4 Time-Altitude Variation
The monsoon anticyclone shows substantial intra-seasonal oscillations, which are connected to variable forcing from transient deep convection over the Indian subcontinent and the Bay of Bengal. This variability is typically associated with active/break cycles of the monsoon with timescales of ~10-20 days [57,58]. Significant correlations exist between outgoing longwave radiation (OLR; Figure 3a2,d2) and circulation within the monsoon region, such that the entire balanced anticyclone varies in concert with convective heating: enhanced convection leads to warmer tropospheric temperatures, stronger anticyclonic circulation, and colder lower stratospheric (and tropopause) temperatures [14]. This causes a significant heterogeneity of the flux transported upward and CH4 concentration in the upper layers during ASMA (Figure 10).
The IGP region experiences intense agricultural activity and use of traditional biofuels. In the winter months, the IGP is often enveloped by thick fog and haze. The prevailing winds at low altitudes (surface to ~850 hPa) are northerly to northwesterly with low wind

CH 4 Time-Altitude Variation
The monsoon anticyclone shows substantial intra-seasonal oscillations, which are connected to variable forcing from transient deep convection over the Indian subcontinent and the Bay of Bengal. This variability is typically associated with active/break cycles of the monsoon with timescales of~10-20 days [57,58]. Significant correlations exist between outgoing longwave radiation (OLR; Figure 3a2,d2) and circulation within the monsoon region, such that the entire balanced anticyclone varies in concert with convective heating: enhanced convection leads to warmer tropospheric temperatures, stronger anticyclonic circulation, and colder lower stratospheric (and tropopause) temperatures [14]. This causes a significant heterogeneity of the flux transported upward and CH 4 concentration in the upper layers during ASMA (Figure 10).
The IGP region experiences intense agricultural activity and use of traditional biofuels. In the winter months, the IGP is often enveloped by thick fog and haze. The prevailing winds at low altitudes (surface to~850 hPa) are northerly to northwesterly with low wind speeds (<5 m/s), and the eastern parts of the IGP are impacted by a localized area of strong subsidence in winter [59]. These conditions tend to trap the pollution at low altitudes [5]. speeds (<5 m/s), and the eastern parts of the IGP are impacted by a localized area of strong subsidence in winter [59]. These conditions tend to trap the pollution at low altitudes [5]. Figure 10. Time-altitude cross-section of CH4 from GOSAT-TIR retrieval, GOSAT-TIR a priori, and ACTM (the left, middle, and right panels, respectively) for the Arabian Sea, Arid, and EIGP regions (for all regions please see Figure S7). Note that the profiles are shown for the tropospheric altitudes because the GOSAT-TIR retrieval system is not sensitive to the stratospheric altitudes (see Figure 8 and the associated text).
Apparently, the a priori profiles from the NIES model show insufficient vertical transport due to inefficient convective parameterization required to simulate tracer transport under monsoon conditions ( Figure 10). This problem forced the transition to a more sophisticated reanalysis JRA-25/JCDAS (Japanese 25-year Reanalysis/Climate Data Figure 10. Time-altitude cross-section of CH 4 from GOSAT-TIR retrieval, GOSAT-TIR a priori, and ACTM AK Cao (the left, middle, and right panels, respectively) for the Arabian Sea, Arid, and EIGP regions (for all regions please see Figure S7). Note that the profiles are shown for the tropospheric altitudes because the GOSAT-TIR retrieval system is not sensitive to the stratospheric altitudes (see Figure 8 and the associated text).
Apparently, the a priori profiles from the NIES model show insufficient vertical transport due to inefficient convective parameterization required to simulate tracer transport under monsoon conditions ( Figure 10). This problem forced the transition to a more sophisticated reanalysis JRA-25/JCDAS (Japanese 25-year Reanalysis/Climate Data Assimilation System developed by the Japan Meteorological Agency (JMA)) and the adaptation of the new parameterization as described by [55]. The GOSAT-TIR retrieval tries to compensate for such a concentration deficit in the upper troposphere.

Seasonal Variation of CH 4
The Prophet analysis model was implemented to derive mean seasonal cycles of CH 4 from GOSAT-TIR and MIROC4-ACTM for the levels of 800 and 300 hPa for 2009-2014 ignoring year-to-year variations ( Figure 11). Seasonal changes are controlled primarily by meteorological parameters; therefore, the effect of the summer monsoon is most noticeable because it causes enhanced redistribution of CH 4 from lower layers to the middle and upper troposphere. As a result, in summer, the minimum CH 4 is found in the lower troposphere (800 hPa), while the maximum is in the upper part ( Figure 11). The amplitude of seasonal changes is determined by the net amount of the sources; therefore, it increases from south to north from marine regions to the most densely populated areas. In the same direction, the difference between the Cao and WH fluxes increases and reached maxima in the three most northern regions (WIGP, EIGP and NEI). This inequality determines the difference between seasonal variability for host regions and spreads to the nearest neighbors. It is especially noticeable for AI and WI, where intrinsic fluxes are much smaller.
With the onset of autumn, the deep convective transport is suppressed; therefore, under the influence of the Hadley cell circulation, the slow outflow of air masses starts in the opposite south-west direction. This situation is marked by the peak of concentration at 800 hPa, which slowly moves from the northern regions (over EIGP in October) to the southern (over the Arabian Sea in late November). Despite a significant (about 50%) difference in Cao and WH emissions, the concentrations smoothed by Equation (1) differ significantly less than without implementation of the AK functions ( Figure S8).
The numerical values of the seasonal cycle derived by the Prophet model ( Figure 11) and averaged for individual seasons and levels of 800 and 300 hPa are presented in Table 1. Based on these results, we found that the amplitude of seasonal changes at the level of 800 hPa was maximum for the regions in the southern and central regions of the domain (SI, WI, CI, EI), where it reached values of 100-140 ppb. At the 300 hPa level, the maximum amplitude was much lower (up to 50 ppb) in regions with high convective activity and adjacent to them (WIGP, BB, AI, WI). The seasonal cycles at the levels of 500 and 300 hPa are very similar; therefore, the free troposphere is not given in the table. Table 1. Mean and 1-σ SD uncertainty of multi-year (2009-2014) seasonal variation of CH 4 at the levels of 800 and 300 hPa for JFM, AMJ, JAS, and OND derived from GOSAT-TIR by the Prophet model as shown in Figure 11. West IGP -24.9 ± 6.0 -11.9 ± 2.9 -20.2 ± 14.3 -10.0 ± 11.3 42.7 ± 8.5 29.5 ± 9.8 2.7 ± 26.0 -7.7 ± 5.3

Regional CH 4 Emission Estimation
India occupies a large region of South Asia, where fewer observations limit the chance to reduce the uncertainty in the greenhouse gases (including CH 4 ) flux. Used in this work, the Cao and WH flux combinations for the South Asia region for the period 2009-2014 account for 65.7 ± 5.8 and 82.4 ± 7.3 Tg year −1 , respectively. Figure 12 shows the time variation of modelled and observed CH 4 averaged over the area of South Asia for three levels (300, 500, and 800 hPa). A significant difference in emissions appears during the monsoon, with almost equal effects in the middle and upper troposphere. Inter-annual change of trend is weak but increasing towards the top of the boundary layer (Figure 12c), where GOSAT-TIR variability can be greatly influenced due to the large spread (SD) of individual samplings (Figure 4). Another important factor is the GOSAT-TIR retrieval a priori profile derived from NIES TM with the coarse vertical grid and simplified scheme for modeling of the boundary layer height, which shows strong diurnal and seasonal variations [60]. Good consistency in the phase is found, because the ACTM AK Cao , and ACTM AK WH are overlapped during non-monsoon periods, and concentration mismatch is strongly associated with the flux difference.

Comparison with Ground Observations
GOSAT-TIR measurements are less sensitive to CH 4 concentrations near the surface, which makes it difficult to discuss based on direct comparison with ground-based observations. However, such comparisons are useful for evaluating the MIROC4-ACTM model and a priori distributions. The number of ground observation stations in India is limited; therefore, we also included observations from the surrounded regions for an overall assessment of the comparison results.  -axis), respectively. The same results, but without implementation of the AK functions is shown in Figure S9.

Comparison with Ground Observations
GOSAT-TIR measurements are less sensitive to CH4 concentrations near the surface, which makes it difficult to discuss based on direct comparison with ground-based observations. However, such comparisons are useful for evaluating the MIROC4-ACTM model  AK Cao , and ACTM AK WH , respectively. Solid and dashed lines are for monthly and yearly averaged concentrations (left y-axis), and dotted line shows the difference between the model simulations (right y-axis), respectively. The same results, but without implementation of the AK functions is shown in Figure S9.
Selected sites represent a diverse variety of CH4, which has different interpretations. The remote marine site (SEY) and high-mountain sites (WLG) are influenced by largescale processes. Low bias (~5 ppb) and SD (~20 ppb) confirm the ability to reproduce the global CH4 budget. The BKT, CRI, PBL, and PON sites represent coastal regions with variable local meteorology controlling CH4. Judging by the correlation coefficient (r), the performance of the model is quite high (r = 0.9), although it decreases in mountainous and densely populated areas (r = 0.6).

Discussion
Analysis of CH4 from GOSAT-TIR and MIROC4-ACTM has several important aspects underlining GOSAT-TIR retrieval itself, but also its implementation on studies over India. Despite essential progress, the development of satellite methods for studying at-

Discussion
Analysis of CH 4 from GOSAT-TIR and MIROC4-ACTM has several important aspects underlining GOSAT-TIR retrieval itself, but also its implementation on studies over India. Despite essential progress, the development of satellite methods for studying atmospheric CH 4 is obstructed by a few limitations. The launch rate of new orbital instruments is significantly ahead of the development of a ground-based and aircraft measurement network for their validation. Due to the scarcity of suitable aircraft observations over India, validation of GOSAT-TIR profiles cannot cover all seasons and land regions studied in this work.
In comparison with AIRS satellite observations [61], the finer vertical resolution of GOSAT-TIR allows capturing more detailed features in CH 4 vertical profiles and therefore more prominent patterns related to different regions and seasons. The significant difference in the CH 4 profile gradient, and its seasonal variability (winter and summer) between a priori and the MIROC4-ACTM model, was revealed in the UT/LS zone (150-200 hPa). During monsoon season, daily variation in a priori profiles is also found in the middle troposphere. While [8] mainly used the model simulations to understand the vertical transport (after validating the model using GOSAT-SWIR measurements), we showed the seasonal evolution of transport and emissions on CH 4 at different layers of the troposphere using both the model and GOSAT-TIR measurements. Unlike work conducted by [7], we found a strong difference between seasons in the middle and upper troposphere caused by variability in atmospheric circulation and vertical convection.
A priori profiles play an essential role in processing a satellite signal, especially for the CH 4 , which has a significant change in a lifetime with altitude. The choice of such a profile (usually provided by model calculations) is critical because of small DOF in observed spectra; the retrieval algorithms cannot overcome large errors in input a priori data. The TransCom-CH 4 experiment [44] showed a significant scatter between the participated models, including the NIES model later selected for calculating the GOSAT-TIR a priori profile (described in Section 2.2). In this study, the difference in the CH 4 profile gradient, its seasonal variability (winter and summer) between a priori and the MIROC4-ACTM model was revealed in the UT/LS zone (levels of 150-20 hPa). This emphasizes the difference in vertical resolution and reanalyzes for the calculation of vertical profiles (described in Sections 2.3 and 2.4). However, the reason for the misfits (GOSAT-TIR a priori vs. ACTM) is difficult to reconcile without additional studies with the use of custom a priori profiles in retrieval.
In the case of long-term projects, the updating of a priori data in accordance with the current progress becomes important. From the moment of the GOSAT launch, the calculations of a priori profiles are carried out according to the same scheme. This is important for the long-term consistency of the TIR CH 4 product but does not take into account the recent improvements (e.g., new reanalysis data, higher vertical resolution, and convective parameterizations) implemented in MIROC4-ACTM [36], as well as further understanding of the CH 4 budget [62]. The release of a new version of retrieval algorithms designed to improve and update satellite products is not regular; therefore, perhaps, for further progress, the retrieval algorithms may use some custom a priori information for checking initial results. Therefore, this will require greater transparency of technical information from satellite projects and significant optimization of retrieval calculations, because such tasks require large computational resources. However, this a priori profile issue remains hidden (but no less relevant) in the field of the main efforts of the scientific community working with XCH 4 derived from the SWIR band, apart from applying the a priori correction.
Another issue of concern, as we found, was the implementation of AK functions. Without AK processing, mean MIROC4-ACTM and GOSAT-TIR mismatches were within 50 ppb, except for the level of 150 hPa and upward, where the GOSAT-TIR sensitivity becomes very low. Convolved model profiles with retrieval a priori and AK functions reduce the mismatch to below uncertainty, making the influence of the a priori profiles too large.
The lack of reliable data for the estimation of regional CH 4 budgets can be mitigated by TIR CH 4 . Relying on our comparison (Figure 12), we suggest the Cao flux combination with the annual mean South Asian emission of 65.7 ± 5.8 Tg year −1 for the period 2009-2014 as more plausible. This confirms the assessment made by [63], indicating that the EDGAR inventory (version 4.2FT2010) with a value of 73.3-83.2 Tg year −1 overestimated the South Asia regional emission by 10-15 Tg year −1 . A significant part of the extra fluxes was concentrated in a few relatively small regions in Northern India ( Figure 2). However, our best estimate emission of 51.2 ± 4.6 Tg year −1 over India is much greater than the 19.6-24.3 Tg year −1 estimated by [6], who combined in situ data of different time coverage and SWIR CH 4 retrievals in the trajectory-based modelling framework, but is in good agreement with the mean emissions of 48.9 Tg year −1 calculated by [40]. Simulation with two scenarios showed that during the monsoon, significant CH 4 amounts due to extra fluxes can quickly propagate to the UT/LS zone and not be detected by ground-based measurements.
Overall, the MIROC4-ACTM simulations of CH 4 in the Indian regions are consistent with the GOSAT-TIR samplings, in terms of seasonality and global variability. Although we found noticeable error in the model data in phase and amplitude at the end of the summer-fall period, the model performance in CH 4 transport in the troposphere and the lower stratosphere was improved due to the use of MIROC4.0 as the meteorological model. Inconsistencies seen in the model-observation comparisons provide opportunities for further flux optimization with inverse modeling methods. More insights are expected after the extension of the released data period of the TIR CH 4 product.

Conclusions
CH 4 from GOSAT-TIR observations and MIROC4-ACTM simulations, sampled at the location and time of the satellite overpass, were analyzed over India and surrounding oceanic regions for the period 2009-2014. The area of our analysis was subdivided into eight land and two ocean regions. We studied the seasonal variation of CH 4 at the atmospheric levels of 800 and 300 hPa over Indian regions and found a large time-spatial inhomogeneity. The major factors controlling it included change in local emission strength, variability in atmospheric circulation, and vertical convection caused by the Asian summer monsoon.
GOSAT-TIR observations are less affected by cloud conditions and could be performed during night; therefore, suitable data coverage and density (10-200 observation samples per 3.0 • × 3.0 • grid cell per month) are obtained to study CH 4 over India. By capturing signals from the top of the atmospheric boundary layer (except for high mountain regions) up to UT/LS, detailed horizontal and vertical features of CH 4 are revealed. Differences from the a priori values suggest that the CH 4 concentration signal was captured even in the lower part of the troposphere, where sensitivity of the GOSAT-TIR sensor is relatively weak compared to the middle and upper troposphere. The mean TIR CH 4 model mismatch was within 50 ppb, except for the altitude range above 150 hPa, where the sensitivity of GOSAT-TIR observations becomes very low.
The MIROC4-ACTM simulations were conducted using two surface fluxes optimized by the inverse analysis describing different approaches of modelling CH 4 emissions from wetlands. Among these scenarios, the Cao scheme seems to be more balanced than WH for individual regions and the whole of South Asia during the monsoon season. In the other periods, no strong difference was found. The flux difference is effective to show the strength of surface concentration propagation to the middle and upper troposphere. Using comparison with GOSAT-TIR, we identified the annual mean emission for the period 2009-2014 for South Asia and India as of 65.7 ± 5.8 and 51.2 ± 4.6 Tg year −1 , respectively.
Despite substantial improvements and testing in various validation projects [44,63], CTMs are still considered as an insufficiently reliable method of studies. However, for many regions/cases, it is only one available source of new knowledge, that would identify general trends, indicate optimal methods of development, and motivate the making of accurate field/aircraft measurements. Our results suggest that an additional comparison with aircraft observations is necessary to analyze the GOSAT-TIR and MIROC4-ACTM mismatch found above the level of 150 hPa. Furthermore, the selection of the a priori model for satellite data retrieval could play a significant role and should be addressed in the developments of future retrieval systems.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/rs13091677/s1, Figure S1. Latitude-longitude distributions of TIR CH 4 properties at the levels of 800, 500, and 300 hPa (the left, middle, and right panels respectively) for the season JAS 2011. The first upper panels (a1-a3) show GOSAT-TIR CH 4 observation points numbers, the bottom panels (d1-d3) show GOSAT-TIR CH 4 1-σ SD, respectively. Figure S2. Same as Figure S1, but for JFM 2011. Figure S3. Latitude-longitude distributions of CH 4 at the levels of 800, 500, and 300 hPa (the left, middle, and right panels, respectively) observed by GOSAT-TIR for the season JFM 2011. The top (a1-a3), middle (b1-b3) and bottom (c1-c3) rows show TIR CH 4 , TIR CH 4 a priori, and difference between the TIR CH 4 and a priori distributions, respectively. Figure S4. Same as Figure S3, but for OND 2011. Figure S5. Latitude-longitude difference in CH 4 distributions simulated by MIROC4-ACTM and observed by GOSAT-TIR at the levels of 800, 500, and 300 hPa (the left, middle, and right panels, respectively) for JFM 2011. The top row (a1-a3) and bottom row (b1-b3) show the difference in CH 4 between GOSAT-TIR, and ACTM Cao and ACTM WH , respectively. Figure S6. Same as Figure S5, but for OND 2011. Figure S7. Time-altitude cross-section of CH 4 form GOSAT-TIR retrieval, GOSAT-TIR a priori and ACTM AK Cao (the left, middle, and right panels respectively) for considered regions. Figure S8. Multi-year (2009-2014) seasonal variation of CH 4 (right y-axis) derived by implementation of the Prophet model for levels of 800 (red lines) and 300 hPa (blue lines) over considered regions from GOSAT-TIR (solid line), ACTM Cao (dashed line), and ACTM WH (dotted line), respectively. At the background bar plots represent Cao (dark grey) and WH (light grey) CH 4 fluxes (left y-axis), respectively. Please note the different scale of y-axes (left) for fluxes. Figure S9. Time series of CH 4 averaged over the area of South Asia for levels of (a) 300, (b) 500 and (c) 800 hPa. Symbols state the GOSAT-TIR observations, red and blue lines are for ACTM Cao , and ACTM WH , respectively. Solid