Interannual Variability of Atmospheric CH 4 and Its Driver Over South Korea Captured by Integrated Data in 2019

: Understanding the temporal variability of atmospheric methane (CH 4 ) and its potential drivers can advance the progress toward mitigating changes to the climate. To comprehend interannual variability and spatial characteristics of anomalous CH 4 mole fractions and its drivers, we used integrated data from different platforms such as in situ measurements and satellites (TROPOspheric Monitoring Instrument (TROPOMI) and Greenhouse Gases Observing SATellite (GOSAT)) retrievals. A pronounced change of annual growth rate was detected at Anmyeondo (AMY), Republic of Korea, ranging from − 16.8 to 31.3 ppb yr − 1 as captured in situ through 2015–2020 and 3.9 to 16.4 ppb yr − 1 detected by GOSAT through 2014–2019, respectively. High growth rates were discerned in 2016 (31.3 ppb yr − 1 and 13.4 ppb yr − 1 from in situ and GOSAT, respectively) and 2019 (27.4 ppb yr − 1 and 16.4 ppb yr − 1 from in situ and GOSAT, respectively). The high growth in 2016 was essentially explained by the strong El Niño event in 2015–2016, whereas the large growth rate in 2019 was not related to ENSO. We suggest that the growth rate that appeared in 2019 was related to soil temperature according to the Noah Land Surface Model. The stable isotopic composition of 13 C/ 12 C in CH 4 ( δ 13 -CH 4 ) collected by ﬂask-air sampling at AMY during 2014–2019 supported the soil methane hypothesis. The intercept of the Keeling plot for summer and autumn were found to be − 53.3‰ and − 52.9‰, respectively, which suggested isotopic signature of biogenic emissions. The isotopic values in 2019 exhibited the strongest depletion compared to other periods, which suggests even a stronger biogenic signal. Such changes in the biogenic signal were affected by the variations of soil temperature and soil moisture. We looked more closely at the variability of XCH 4 and the relationship with soil properties. The result indicated a spatial distribution of interannual variability, as well as the captured elevated anomaly over the southwest of the domain in autumn 2019, up to 70 ppb, which was largely explained by the combined effect of soil temperature and soil moisture changes, indicating a pixel-wise correlation of XCH 4 anomaly with those parameters in the range of 0.5–0.8 with a statistical signiﬁcance ( p < 0.05). This implies that the soil-associated drivers are able to exert a large-scale inﬂuence on the regional distribution of CH 4 in Korea.


Introduction
Methane (CH 4 ) is the second most effective greenhouse gas after carbon dioxide (CO 2 ), and therefore a substantial contributor to global warming. CH 4 has multiple sources, including natural sources such as wetlands, oceans, and termites, that account for 35-50% of total emissions, and the rest from anthropogenic sources such as agricultural soils, landfill, rice paddies, ruminants, biomass burning, and energy utilization [1]. According to the Food and Agriculture Organization (FAO), the livestock industry emits 37% of anthropogenic CH 4 , mainly from ruminant enteric fermentation. CH 4 is produced in the soils under anoxic environments by methanogenic microbes during the anaerobic digestion of organic matter or removed from the soils by microbial oxidation, which is carried out in the aerobic zone of methanogenic soils (methanotrophy) and in upland soils that oxidize atmospheric methane [2]. There are environmental factors that affect the process, for example, (i) change of water content in the soil that alters gas diffusion in connection with the oxido-reduction level and CH 4 transfer; (ii) the soil temperature, which influences the microbial activities in the soil, i.e., low soil temperatures decrease CH 4 production by decreasing the activity of methanogens [2][3][4]. Methane is mainly removed through reaction with hydroxyl radical (OH) in the troposphere and stratosphere while its sources are diverse [5]. CH 4 has a relatively short lifetime (~9 years) [6] after emitting into the atmosphere, so we can see the effect in a reasonably short time when we reduce its emission. However, the geographical distribution and sectorial attribution of methane emissions, and the interannual variations of the sources, are uncertain [5,7,8]. This gap hampers the effective formulation of emission-mitigation strategies. Recent global atmospheric CH 4 levels had reached 1867 ppb by the end of the 2010s and are growing faster than at any time in the past two decades [9]. Between 2007 and 2013, global CH 4 has annual growth rates of about 6 ppb yr −1 [10], and has accelerated to 10 ppb yr −1 between 2014 and 2018 [6]. The drivers of this increase are partly addressed and may embrace enhanced emissions and a possible decline in the destruction of methane in the air [6]. Especially, it is known that the growth rate is related to ENSO, e.g., [10,11], while the growth rate increased in 2019 without ENSO over Korea, which is the main focus of this paper for investigation.
Asia is a substantial contributor to climate change, with greenhouse gas emissions not only from CO 2 [12] but also from a CH 4 perspective. East Asia covers highly populated and industrialized areas with large anthropogenic CH 4 sources in the Eurasia continent [13]. Specifically, considering CH 4 distributions over South Korea, 20% of the total emission was derived from the agriculture sector. The other significant sources of CH 4 are enteric fermentation and livestock manure management, which were estimated to have annual growth rates of 1.7% and 2.6%, respectively, from 1990 to 2009 [14]. Complying with the Paris Agreement, the South Korean government recently declared a roadmap to reduce its GHG emissions by 2030, which proposed a reduction in overall greenhouse gas (GHG) emissions from the agriculture sector by 5.2% [15] and by 2050, it will attain a net-zero carbon emission. Therefore, in this context, advancing the understanding of the interannual variations of CH 4 with a focus on indicating anomalous events and their major causes of variations are essential inputs for carbon-mitigation strategies.
To understand the increase of CH 4 in 2019 over South Korea, integrated data comprising of near-surface CH 4 observations from the surface in situ and column-averaged CH 4 retrievals from space-borne instruments (e.g., Greenhouse Gases Observing SATellite (GOSAT) [16], TROPOspheric Monitoring Instrument (TROPOMI) [17] were utilized. The World Meteorological Organization/Global Atmosphere Watch (WMO/GAW) Anmyeondo (AMY) regional station has the longest continuous records of atmospheric CH 4 from 1999 till present in Korea. To explore the possible attributes for the large growth of CH 4 , ancillary data for soil temperature soil moisture content at 10 cm depth parameters from the Land Data Assimilation System (FLDAS) Noah Land Surface Model (LSM) L4 [18] and the El Niño-Southern Oscillation (ENSO) were used to assess the impact of the global climate [11]. We also used 13 C/ 12 C isotopes ratios in CH 4 (δ 13 C-CH 4 ) data from flask-air sampling at AMY to infer the information about the source signature [19].
This paper attempts to address the large growth CH 4 in 2019, characterize their spatial distributions, indicate the source signature, and explore the potential drivers of anomalous atmospheric CH 4 mole fractions with integrated data such as in situ surface and satellite data over South Korea.

AMY Near-Surface In Situ Data
The AMY station is located at the west coast of South Korea at 36.54 • N, 126.33 • E, and 46 m above sea level ( Figure 1). Its climatic conditions are described as winter (December-February) is the coldest season (average temperature around −2 • C) with north-westerly winds while summer (June-August) is the warmest season (average temperature around 25 • C) with southerly winds. For the precipitation, according to each year, the accumulated amount can be different, but its value is highest in summer due to summer monsoon with 362-1267 mm in recent 10 years (2011-2020) while winter observes the lowest level, less than 180 mm. It is a WMO/GAW regional station. Several GHGs have been regularly monitored, of which atmospheric mole fractions of CH 4 from 40 m tower has been continuously measured using Cavity Ring Down Spectroscopy (CRDS 2401, Picarro, CA, USA) based on 5 s intervals from 2016. Before the CRDS, we monitored CH 4 with a Gas Chromatography-Flame Ionization Detector (GC-FID) based on 30 min intervals. This instrument is calibrated against the WMO-X2004A scale with 4 reference tanks every 2 weeks while the GC-FID with 1 reference gas every 6 h. The difference between hourly data from both CRDS and GC-FID and weekly flask-air samples analyzed by National Oceanic and Atmospheric Administration (NOAA) are within the WMO/GAW compatibility goal ±2 ppb. Here, compatibility goal is defined as "property of a set of measurement results for a specified measured, such that the absolute value of the difference of any pair of measured quantity values from two different measurement results is smaller than some chosen multiple of the standard measurement uncertainty of that difference" and also "Meteorological compatibility of measurement results replaces the traditional concept of "staying within the error", as it represents the criterion for deciding whether two measurement results refer to "the same measured or not" [20]. Since the compatibility goal for CH 4 is ±2 ppb, we can consider the differences are almost the same with NOAA. NOAA has a role in the Central Calibration Laboratory, which serves the primary standard for CH 4 in the GAW network; our results were very close to the standard, and our calibration method is also reliable. To highlight the source signatures of CH 4 , flask-air measurements of 13 C/ 12 C in CH 4 (δ 13 C-CH 4 ) [19] (ftp://aftp.cmdl.noaa.gov/data/trace_gases accessed on 4 February 2021) were analyzed. Near-surface data used in this work covered 2014-2020. Further details about the station and instruments can be found in [21]. CH 4 measurements at AMY provide a constraint on the regional methane source strength, which is seen in Tae-ahn Peninsula in close proximity to AMY (about 28 km) [22]. On account of this, the purpose of using the station data in this analysis is to understand whether the behavior of CH 4 growth rate, source signature, and the possible attributes can share the same behavior with the regional average over Korea or not.

TROPOMI Retrievals
The TROPOMI instrument onboard the Sentinel-5 Precursors (S5P) [17] satellite was launched on 13 October 2017 and was put into a near-polar, sun-synchronous orbit with a mean altitude of 824 km above Earth's surface, an Equator crossing time of 13:30 LT, and a revisit time of 1 day. Its daily global coverage was provided at a spatial resolution of 7 × 7 km 2 since its launch in October 2017 and improved to 5.5 × 7 km 2 in August 2019. Data were screened with a retrieval quality flag of 1.0 (QC = 1.0). S5P bias-corrected XCH 4 has been evaluated against TCCON data, and the accuracy and precision were found to be 0.25% and 0.57%, respectively [23]. Data spanning April 2018-December 2020 were utilized to analyze the spatial variations of interannual variability of XCH 4 through finding anomalous distributions with the focus on 2019 over Korea, since it provides much denser spatial sampling as compared to GOSAT.

GOSAT Retrievals
Greenhouse Gases Observing SATellite (GOSAT) was launched into a sun-synchronous orbit on 23 January 2009 by an H-IIA launch vehicle and was positioned in a sun-synchronous orbit at a 666 km altitude. It has a 3-day revisit orbit cycle. The TANSO-FTS onboard GOSAT makes global observations both with the nadir and off-nadir modes, and makes use of the spectral band Short Wavelength InfraRed (SWIR) for deriving CH 4 [24]. Version 9.0 of the University of Leicester (UoL) GOSAT Proxy XCH 4 data [25,26] from 2014 to 2019 were utilized. Parker et al. [26] validated the proxy XCH4 data against ground-based TCCON (Total Carbon Column Observing Network) observations and found an accuracy of 4.8 ppb (~0.27%) and a single precision of 13.4 ppb (~0.74%). Primary use of these data is to analyze the seasonal and interannual variability of XCH 4 of the domain averaged over the region of interest near AMY, though the spatial sparsity of the data coverage is a caveat.

Ancillary Parameters
To understand the attributes for interannual variability and spatial characteristics of atmospheric CH 4 anomalies captured by the instruments, soil temperature and soil moisture content at 10 cm depth with a spatial resolution of 0.1-degree of latitude-longitude and monthly temporal resolution from the Famine Early Warning Systems Network (FEWS NET), the Land Data Assimilation System (FLDAS) and Noah Land Surface Model (LSM) L4 (https://disc.gsfc.nasa.gov/datasets/FLDAS_NOAH01_C_GL_M_001 accessed on 25 January 2021) were used. This model simulation was forced by a combination of the Modern-Era Retrospective analysis for Research 2 (MERRA-2) and Climate Hazard Group Infrared Precipitation with Station (CHIRPS) 6-hourly rainfall data [18] (McNally, Chicago, IL, USA, 2018). FLDAS monthly products of soil moisture were compared to remotely sensed observations and the overall agreement was good [27]. As described in Section 1, soil temperature and soil moisture are the trigger factors for changes in CH 4 emissions coming from agriculture soils through affecting the methanogenic and methanotrophic activities in the soils [2]. South Korea's bottom-up CH 4 surface-emission data from agriculture sector were used to look into the annual variations along with CH 4 growth rate interannual variations.
The El Niño-Southern Oscillation (ENSO) index (https://psl.noaa.gov/data/climat eindices (accessed on 22 January 2021)) was exploited to analyze the influence of ENSO on the interannual variability CH 4 growth rate [11]. ENSO is an indicator of interannual variability in global climate system [28], which affects seasonal to decadal global climate and impacts temperature and precipitation.

Data Analysis Method
To analyze the seasonal variability and long-term trend of the data from the surface and satellites after deseasonalized, we applied Thoning's method [29]. Seasonal variability was extracted from the detrended data, and the long-term trend was derived from the deseasonalized data. Data were selected here to compare the column data. Therefore, the surface data can reflect the activities on a regional/local scale. The annual increase of atmospheric CH 4 in a given year is the increase of mole fraction from January 1 in that year to January 1 of the next year after removing the seasonal cycle. The yearly or monthly anomaly was computed as subtracting the annual means from the individual mean estimates. We also applied Pearson's correlation coefficient to quantify how well the interannual variability of CH 4 co-varied with drivers such as soil temperature and soil moisture. When aggregating TROPOMI overpass soundings within 0.1 • × 0.1 • , we collected all soundings that fall in pixel and made daily averages from which we calculated monthly means. The spatial distributions of TROPOMI data that were made by a cluster average within 0.1 • by 0.1 • latitude-longitude have the same spatial resolution with that of soil temperature and soil moisture data. The gridded data were then averaged on a monthly basis. While selecting satellite soundings at AMY, a spatial coincident criteria of 2.0 • by 2.0 • of latitude-longitude was applied. To infer the source signatures of CH 4 , the flask-air measurements of δ 13 C-CH 4 ratio in CH 4 at AMY were utilized. Furthermore, a Keeling plot [19] was implemented to characterize source signatures responsible for CH 4 enhancement. At the station, seasonality was well captured by both satellites (TROPOMI and GOSAT) and in situ. It was noted that the interannual variability of seasonal maximum was higher than the seasonal minimum but the trend of their periods remained to be consistent.  (Figure 2, right panel). The seasonal cycle of CH 4 in the Northern Hemisphere is more complex [30]. While the destruction of CH 4 due to reaction with OH is expected to be stronger in summer, sources are diverse and also strongly vary with the seasons. As for describing the station locations, for the AMY station, within 35 km to the northeast and southeast, there are the largest thermal power plants fired by coal and heavy oil. In addition to that, agricultural activities such as rice paddies surrounding this station are known to affect the high CH 4 in this region [21]. The ULD station is placed in the eastern part of Ulleungdo Island (an area of 72 km 2 ), which is in the southwest of South Korea, and 155 km away from the mainland. There are farming activities but the station is located at the mountain (220.9 m from the ground), and is less affected by those activities. JGS station is located in the west part of Jeju Island, which is in the southwest of South Korea and 90 km far from the mainland (Figure 1). JGS station from the southwest to the northwest is open to the sea, and the wind speed is the highest among the three stations, so that it is hard to capture the local CH 4 sources. Details about the sampling sites are provided in [20]. Therefore, strong local sources might not affect strong local sources that can influence the seasonal cycle of CH 4 at JGS and ULD, while transport could be playing a major role in modulating the seasonal behavior with a depleted CH 4 signal from the ocean. Unlike JGS and ULD, AMY is located in the proximity of local sources, indicating that the seasonal maximum (2042.7 ppb) and minimum (1982.1 ppb) occurred in the summer and spring, respectively. In ULD and JGS, the minimum values were detected in summer, with 1936.8 ppb and 1933.0 ppb, respectively, whereas the maximum was observed in the winter (ULD) and autumn (JGS), with mean values of 1968.7 ppb and 1977.7 ppb, respectively. Such complex information will be refined further in the subsequent study focused on source investigation of CH 4 and CO 2 over Korea GAW stations. When comparing with the surface data in other global stations located in the Northern Hemisphere, the AMY result is similar to the seasonal cycle observed at Waliguan [31]. The spatial differences on the seasonal cycle of CH 4 depend on the regional sources and sinks, and atmospheric transport. The behavior of the seasonal cycle of XCH 4 integrated over Korea has a similar trend to AMY (Figure 2, right panel) given that the actual magnitudes of amplitude were expected to be varied. This suggests the dominant drivers are the same.

Interannual Variability of CH 4 Growth Rate in Korea
The annual increase was calculated from the trend line (Figure 3a) as the increase from January 1 in 1 year to January 1 in the next year, after removing the seasonal cycle. The successive annual differences of CH 4 at AMY were determined in the range of −16.8 ppb yr −1 -31.3 ppb yr −1 by in situ and 3.9 ppb yr −1 -16.4 ppb yr −1 by GOSAT ( Table 1). The annual averages CH 4 estimated from in situ and GOSAT were 7.55 ± 23.59 ppb yr −1 and 8.92 ± 6.87 ppb yr −1 , respectively ( Table 1). The large growth rates dramatically increased in 2016 (31.3 ppb yr −1 and 13.4, ppb yr −1 as seen by in situ and GOSAT, respectively) and 2019 (27.4 ppb yr −1 and 16.4, ppb yr −1 , in situ and GOSAT, respectively). Though we considered changing instruments from GC-FID to CRDS, this growth rate seems to be related to natural phenomena rather than the measurement system since GOSAT also captured the higher growth rate, and the differences with flasks data are also in the compatibility goal as described in 2.1. Regarding TROPOMI, data in 2018 were started from April so that this large data gap can affect the yearly growth estimate, so we skipped that to discuss the result obtained. Similar to AMY, the spikes of peak growth obtained from the domain averaged were demonstrated in 2016 and 2019. Figure 3 shows deseasonalized long-term trends and the instantaneous growth rate of near-surface and column-averaged CH 4 mole fractions at AMY, as well as integrated over South Korea using in situ and GOSAT. The long-term trend revealed that the increase from 2014 to 2016 corresponding with the cold-to-warm transition from weak La Niña started in the earlier period of 2014 to a strong El Niño event in 2016 (Figure 4). After 2016, the trend has shown to be oscillated, with maximum and minimum growth in 2018 and 2019, respectively. A higher temperature in 2019 was observed in Korea, which could be caused by the presence of atmospheric instability in the region ( Figure S1). A large part of the mid-latitude temperature variability can be described by natural variability, which is essentially chaotic rather than a reflection of a forced signal such as Arctic warming linked to sea ice loss [32]. Geographically, Korea is located in the middle latitudes of the Northern Hemisphere, on the east coast of the Eurasia continent, and also adjacent to the western pacific. Therefore, it unveils complex climate behaviors. In October 2019, high pressure over the north Pacific Ocean was developed, which was larger than the usual in the east-west direction. Such structure affected the temperature and moisture over Korea. After mid-October, anticyclone circulation near the Ural mountain moved to the east. Combined with high atmospheric pressure in the north pacific, it contributed to high temperatures over Korea. Sea-surface temperature was observed higher than normal over the tropical West Indian Ocean. Under this condition, convection was suppressed near the Philippine Sea, whereas active convection was detected near Japan and Korea. In the vicinity, the downdraft of the airmass was noticed and thus contributed to the formation of an anticyclone condition [33].  Table 1. Annual means (ppb) and yearly growth rates (ppb yr −1 ) of in situ CH 4 and GOSAT XCH 4 at AMY derived using curve-fitting method [29] after being deseasonalized. Annual mean is derived from the trend curve, after removing the seasonal cycle, and then yearly growth rate is computed as a difference of successive annual means (i.e., 1 January in that year to 1 January of the next year), and relative increase is given in percentage in open bracket. Growth rate is denoted as GR and relative increase as RI (%).  Quantitatively, the correlation between interannual variabilities of yearly growth rate captured by in situ and soil temperature changes was found to be 0.92 with a statistical significance of a 95% confidence interval. In recent findings, Sun et al. [34] studied the long-term field CH4 flux measurement in the Yangtze River delta, China, and determined that higher temperatures and less rainfall in 2013 caused CH4 emissions to increase from rice paddies by 58-294%. [35] indicated that the rising temperature results in an enhanced CH 4 production per unit of livestock products. In other studies, Lee et al. [36] reported that the forage quality declines with elevating temperature and leads to an increase of enteric CH 4 production, and they indicated the increase is 0.9% • C −1 .

In Situ CH
As we compare the trend and growth rate patterns detected at AMY with the other Korea's in situ stations located in the east (ULD) and south (JGS), it was found to be similar to ULD, but a contrasting pattern with JGS in the period between 2015 and 2017 ( Figure S2), which will be investigated as the reason for disagreement with JGS in the future. The overall temporal characteristics of the growth rate of AMY and the region of interest have similar behaviors of growth rate patterns in the entire period (Figure 3a), and thus imply the significant drivers are the same.
A noteworthy decline of the growth rate in 2020 was registered, which is 1.7 ppb yr −1 . Interestingly, 2020 recognized the COVID-19 pandemic globally, so that many countries have implemented lock-down situations. As a result, this could affect surface emissions for various trace species, and thus can induce the surface temperature change [37]. Despite the reduction of emissions, 2020 was the second warmest year over a period of 2010-2020. As we can see specifically over Korea, 2019 was warmer than 2020 ( Figure S1), which is due to the occurrence of extended summer monsoon Changma in 2020 [38]. Changma is derived from the Korean words for "prolonged rainy period" in summer.

CH 4 Growth Rate Related to Isotopic Ratios in 2019
The isotopic composition of atmospheric methane 13 C/ 12 C (expressed as δ 13 C-CH 4 ) collected by flask-air sampling weekly at AMY was used to understand the source signatures. Figure 5a  We also computed the intercept of the Keeling plot to infer the isotopic signature of the source responsible for enhancements in CH 4 , and the intercept values for summer and autumn were found to be −53.3‰ and −52.9‰, respectively (Figure 5b), which is consistent with biogenic emissions. The aircraft data observed over western Korea in the summer season also revealed biogenic emissions signatures using CH 4 /CO ratio analysis [39]. As reported in [19], high-latitude CH 4 has distinct δ 13 C isotopic signatures, from −60 to −50‰ for CH 4 hydrates; ruminants-digesting C4 plants give off CH 4 at −55 to −50‰, whereas those eating C3 plants give off −65 to −60‰ CH 4 , since C3 and C4 plants have different photosynthetic pathways that result in a different signature. The enhanced depletions found in 2016 and 2019 accounts for enhanced biogenic CH 4 signals that coincided with the high growth rates of CH 4 in those respective years. The observed interannual variability of the CH 4 growth rate associated with the climate variability was not reflected in Korea's CH 4 emission estimate from the agriculture sector from 2014 to 2018, though data for 2019 were not available ( Figure 6).

CH 4 Growth Rate Related to Soil Temperature and Moisture in 2019
Regarding the high growth in 2019 at AMY, we probed the spatial extent of XCH 4 distributions in the same year. The short-term anomalous XCH 4 can contribute to a large growth rate change. In this context, we scrutinized the spatial variations of XCH 4 interannual variability over Korea estimated from the monthly means of TROPOMI data from April 2018 to August 2020. The estimated annual means were predominantly higher over the west, central, and south of most regions (Figure 7a). The illustrated standard deviations highlight the magnitude of interannual variability, and the higher values were clustered over the southwest region, fluctuating within 18-28 ppb (Figure 7b). It was also depicted the number of monthly means data used in Figure 7c. A large increase of XCH 4 was discerned in October 2019, and this was estimated to be a maximum of 70 ppb after subtracting the annual mean ( Figure 8). This change apparently responded to the variations of soil temperature and soil moisture (Figure 9). A higher soil temperature was evident in October 2019. While comparing the actual monthly mean estimate of soil temperature, it was shown above 16 • C in October 2019 over the west and south part of the domain, which is almost 2-3 • C higher than the same month in 2018 (Figure 9a,b). The period of 2019 was characterized by a large positive soil temperature anomaly in the entire domain, which is in contrast to 2018 (Figure 9c,d). Similarly, the spatial moisture variation was also noticeable in those respective periods. Quantitatively, the spatial Pearson correlation of XCH 4 interannual variability with the aforementioned attributes was determined between 0.5 and 0.8 and statistically significant (p < 0.05) (Figure 10). This finding highlighted the influence of the drivers over a large spatial extent in the domain of interest. As discussed in Section 3.3, earlier findings [34][35][36] pointed out that increased temperature causes rising CH 4 emissions from livestock and paddy fields.

Conclusions
In this work, we studied the seasonal and interannual variability of CH 4 , as well as their spatial distributions over South Korea, with the focus on unprecedented growth in 2019. The seasonal cycle of CH 4 at AMY was found to be maximum in summer and minimum in spring, which is similar to the seasonal cycle of XCH 4 averaged over Korea. Thus, the drivers of CH 4 seasonal change over AMY are most likely identical to the regional-scale drivers of CH 4 in Korea. At AMY, the average growth rates computed from in situ and GOSAT were 7.0 ± 19.7 ppb yr −1 and 8.6 ± 5.1 ppb yr −1 , for 2014-2020 and 2014-2019, respectively. The high growth rates were detected in 2016 and 2019, resulting in 27.4-31.4 ppb yr −1 from in situ and 13-16 ppb yr −1 from GOSAT. The high growth rate in 2016 was influenced by a strong El Niño in 2015-2016. On the other hand, the large growth recorded in 2019 was induced by high temperatures due to enhanced regional atmospheric instability. We demonstrate those reasons with two types of data.
(1) The isotopic composition of CH 4 (δ 13 -CH 4 ) collected by flask-air sampling at AMY was analyzed to infer the source signatures, and the result was varied within −47.96‰ to −47.17‰, with an annual average of −47.59 ± 0.19‰. Enhanced depletions appeared in 2016 and 2019 that matched the large growth of CH 4 . These elevated depletions were driven by the summer and autumn seasons. The intercepts of the Keeling plot for summer and autumn were found to be −53.3‰, and −52.9‰, respectively, which is consistent with biogenic emissions that are accountable for the enhancements in CH 4 . The change of soil temperature and soil moisture affected CH 4 emission from biogenic sources.
(2) We also probed the spatial variations of XCH 4 interannual variability with the center of indicating anomalous events over Korea from April 2018 to December 2020 based on TROPOMI observations. We discerned the enhanced spatial distribution of XCH 4 in October 2019. This was characterized by a considerably large anomaly of XCH 4 after removing the annual mean, up to 70 ppb. This suggests that the short-term enhancement can appreciably contribute to the observed yearly growth rate changes. We also found that the likely attributes for such changes were soil temperature and soil moisture content. More importantly, the increase in soil temperature was very pronounced during October 2019. Furthermore, the overall correlation of the interannual variability of XCH 4 and those attributes was determined to be in the range of 0.5-0.8 with the statistical significance of a 95% confidence interval. We can infer that biogenic-related CH 4 with associated drivers such as soil temperature and moisture is relevant in understanding the interannual variability of the growth rate over our domain. Our findings, related to the characteristics of CH 4 and its source characterization derived from the integrated data such as in situ and satellite data, can improve the numerical model uncertainties in the CH 4 budget and source identification (based on isotope information). The quantified relationship of the spatial variability of CH 4 and soil temperature and moisture changes can improve the model uncertainty of CH 4 related to soil variables. When more TROPOMI observations are available, one will also be able to identify temporally stable XCH 4 anomalies and distinguish them from the short-term anomalous CH 4 events in the future. This will open the prospects for providing observational-based contrasts between the local shortterm contributions and long-term (meteorologically driven) exposure of the region to CH 4 transport.