Distribution and Discharge of Dissolved Methane in the Middle Okinawa Trough, East China Sea

Widespread seepage of methane from seafloor sediments on continental margins are released into seawater, a portion of which may escape to the atmosphere. To assess the water column distribution characteristics of methane and its input to the atmosphere, we investigated methane emissions from the shelf and west slope of the back-arc Okinawa Trough (OT), East China Sea. Our results showed a heterogeneity distribution of methane within the water column. The highest value, which was more than 10 times of the background concentration, occurred near a cold seep in the north of the study area which was discovered by a remotely operated underwater vehicle (ROV). Other sources of methane to the water column of the OT, besides cold seepage input, probably also include in situ aerobic methane production, advective transport from the continental shelf, and/or hydrothermal venting. Furthermore, the sea-to-air flux of methane throughout the study area was up to 116 μmol m–2d–1, noticeably higher than that in many other continental shelf waters and seep sites globally, indicating that this region is an active CH4 emission area. Our findings demonstrate that methane discharged from both cold seeps and hydrothermal vents have a significant influence on the methane cycle in the OT, providing a new insight for the methane budget of back-arc basins.


INTRODUCTION
Methane (CH 4 ), a significant greenhouse gas, has a global warming potential about 20-40 times higher than that of CO 2 on average over a century timescale (St. Louis et al., 2000;IPCC, 2013). Although CH 4 emissions from the oceans account for only a small portion (2%) of the global CH 4 budget (Reeburgh, 2007), it still plays an important role in the atmospheric chemistry, airsea interaction, and distribution of chemosynthetic communities. Previous studies estimated that marine CH 4 flux to the atmosphere ranges from 0.4 to 1.8 Tg yr −1 in the open ocean (Rhee et al., 2009), while it reaches 13 Tg yr −1 for continental shelves (Bange et al., 1994). This observation suggested that the continental shelf would be the dominant place for the oceanic CH 4 emissions.
The major sources of dissolved CH 4 in the ocean are terrestrial inputs, sediment emission (including geological sources), and biological metabolism (Fleischer et al., 2001;Jayakumar et al., 2001;Skarke et al., 2014;Mau et al., 2017). Most of this CH 4 gas could be dissolved in seawater (McGinnis et al., 2006), generating patches of high CH 4 concentration (Leifer et al., 2000), but only a small amount of them can reach the atmosphere (Reeburgh, 2007). In this case, the amount of CH 4 estimated based on sea-air gas flux is markedly lower than that discharged from sediment (Milkov and Sassen, 2003). Nevertheless, due to high variability in hydrochemical properties and release rate of CH 4 from sediment, the oceanic CH 4 budget is poorly documented up to date. Therefore, it is highly essential to carry out more measurements of CH 4 production, dissolution, migration, and emissions from the seabed in order to better understand the oceanic CH 4 budget and potential influence on climate change.
The first studies on the characteristics of dissolved CH 4 in the East China Sea (ECS) have been conducted since the early 1990s (Tsurushima et al., 1996). Previous research demonstrated that terrestrial input and hydrographic circulation systems [Kuroshio, the Taiwan Warm Current Water (TWCW), as well as the East China Sea Coastal Current (ESCC)] were the main controlling factors on CH 4 concentration and distribution patterns (Lee et al., 2000;Rehder and Suess, 2001;Liu et al., 2003;Zhang et al., 2004Zhang et al., , 2008. However, given that a variety of mud volcanoes, pockmarks, and cold seeps have been discovered in recent decades in the OT (Yin et al., 2003;Li et al., 2015;Sun et al., 2015;Xu et al., 2018), it is reasonable to infer that they may contribute a substantial amount of CH 4 to seawater. In addition, as a universal extreme system within the OT, modern hydrothermal activity is pervasive and vigorous from north to south, which can also discharge a large amount of CH 4 along with hydrothermal fluids (Sakai et al., 1990;Ishibashi et al., 1995Ishibashi et al., , 2014Inagaki et al., 2006;Miyazaki et al., 2017). Previous studies found that the distance between the site of the newly discovered cold seeps Xu et al., 2018) and the hydrothermal vents found in the past (Inagaki et al., 2006;Ishibashi et al., 2014) are less than several tens of kilometers apart. The proximity of the cold seeps to the hydrothermal fluid complicates the material circulation in the OT . Recently conducted research has confirmed that both hydrothermal and cold seep fluids could contribute a considerable amount of carbon to the seawater Zhang et al., 2019). In such situations, due to the complicated supply of methane-rich fluids derived from cold seeps and hydrothermal vents, the aforementioned region has been an important target to explore the characteristics and sources of CH 4 in the OT.
In this study, we measured dissolved CH 4 concentrations and related physical oceanographic parameters from the fulldepth water, and then estimated CH 4 fluxes at the seaair interface in the OT (Figure 1). Our study aimed to determine the fate of the discharged methane and evaluated the regional contribution to atmospheric methane, it would shed light on better understanding the methane budget and the biogeochemical carbon cycle in the OT.

GEOLOGICAL SETTING
The OT is an active semi-depth back-arc basin formed behind the Ryukyu arc-trench system in the West Pacific (Sibuet and Olu, 1998;Shinjo et al., 1999). It is approximately 1200 km in length with the average width of approximately 104 km and characterized by active rifting structures and magmatism along the depression. In this trough, a large section of the seawater is deeper than 1000 m, and the maximum depth is about 2300 m. The thickness of sediments in the northern section of the OT is up to 8 km, due to a huge supply of terrigenous matter from the ESC continental shelf (Sibuet et al., 1987). A topographical study showed that the OT has a typical U-shape, and submarine canyons, fans, and turbidite deposits have extensively developed on the western slope. Modern hydrothermal activities are pervasive and vigorous in the middle and southern sections of the OT (Glasby and Notsu, 2003). In addition, seafloor cold seeps have been discovered in the northern and central OT recently Sun et al., 2015;Xu et al., 2018).
The hydrological and circulation systems in the OT are complex, and are mainly controlled by the force from the Kuroshio Current and coastal current (Yu et al., 2009). Amongst, the Kuroshio Current, carrying oligotrophic waters northward with low concentrations of nutrients and high temperature and salinity, is the most important factor controlling the environment and hydrological characteristics. It flows above the OT as a water mass up to 100 km wide and 800-1,000 m deep and at a velocity ranging from 45 to 150 cm/s (Qin, 1987). It is the only source of heat and ocean material on the shelf of ECS.

Sampling and Regional Survey
Forty sites of full-depth water column samples were selected in this study from the western slope of the middle OT and continental shelf during the integrated environmental and geological expedition of R/V Zhang Jian from June to September 2016 (Figure 1). Seawater samples from different sites and depths were collected with cleaned Go-Flo bottles (20-L) mounted on a Seabird 911 plus CTD/rosette. Surface waters were collected at a depth of ∼3 m and near-bottom waters were at ∼50 m above the seafloor. Once on board, water was collected into 40 mL glass serum vials using a silicone tube. The vials were thrice overfilled without introducing any air bubbles in the sample. Then, those were sealed with a butyl rubber stopper and aluminum cap. After that, methane analyses of these samples were conducted in the laboratory immediately. Temperature, salinity, and oxygen data were obtained from the CTD profiles.
A high-quality in situ giant clam beds map was obtained by a remotely operated underwater vehicle (ROV) Beaver (Shanghai Jiao Tong University) in the cruise of R/V Zhang Jian in 2017 (Figure 1). The dive sites were characterized by several acoustically detected flares in the water column, which were first discovered by the Kongsberg EM302 multi-beam echosounder on-board.

Dissolved CH 4 Concentration
Headspace gas chromatography was applied for on-board measurements of dissolved gases, as well as high resolution mapping of CH 4 concentrations. Briefly, a 40 mL glass vial with half a sample of seawater was immersed into a 25 • C constant temperature water bath to balance the CH 4 distribution between headspace and seawater. After that, 0.5 mL of headspace gas was injected into a gas chromatograph (GC910 produced by Ke Chuang Company, Shanghai) with a flame-ionization detector. High-purity helium (99.999%) was used as the carrier gas at a flow rate of 30 mL/min. The precision of the analytical method was ± 2.5%.

Calculations of Saturation and Sea-Air Flux
Saturation values R, expressed in%, were calculated as the ratio of the concentration of dissolved gas to the expected equilibrium water concentration. The saturation and sea-air CH 4 fluxes (F, µmol·m −2 ·d −1 ) were calculated as follows: Where C obs is the observed dissolved CH 4 concentration (mol L −1 ); C eq is the CH 4 concentration in the seawater equilibrated with air, which can be calculated using in situ temperature, salinity, and the solubility equation of Wiesenburg and Guinasso (1979). However, atmospheric CH 4 concentration was not determined in the present study. Therefore, a global mean atmospheric CH 4 mixing ratio of 1831.16 ppb of 2013, taken from the NOAA/ESRL Global Monitoring Division in situ program 1 , was used for calculations. k (cm h −1 ) is the gas transfer velocity, which is a function of wind speed and the Schmidt number (Sc). In the current research, the formula presented by Wanninkhof (1992Wanninkhof ( , 2014) was used to calculate the gas transfer coefficient k. The average wind speed during the cruises is 8.3 m s −1 .

Horizontal Distribution of CH 4
We collected the methane data from the shelf and west slope of the OT, which ranged from 2.7 to 24.7 nM, with an average of 6.5 nM. The horizontal distributions of CH 4 in water are shown in Figure 2, presenting a great variation both in the surface and bottom, from 2.0 to 19.9 nM and 2.0 to 24.7 nM, respectively. Careful study found that the CH 4 concentrations in the bottom at about 60% of the sites are approximately 10-270% higher than that in the shallower and surface water (Figure 2 and Supplementary Figure S2). This finding suggests there may be an important benthic CH 4 originating from the underlying sediment, which can be identified through sediment core incubations (Chronopoulou et al., 2017). However, the high CH 4 concentration anomalies are observed on the lower slope of the north study area, with the maximum concentration of 24.7 nM, which is approximately more than 10 times the background concentrations of 2-3 nM in the OT (Tsurushima et al., 1996) and the typical open ocean (Bates et al., 1996;Kock et al., 2008;Tseng et al., 2017). This irregular distribution of CH 4 is possibly pointing to seafloor methane seep from sediment.

Distribution of CH 4 Across the Shelf and Slope
We divided our data into six transects. Transects P1-P5 were conducted across the shelf and slope while P6 followed the contour line of the middle slope (Figure 1). Along the P1 transect (Figure 3), CH 4 concentrations changed in a wide range, from 2.0 to 20.3 nM. The most conspicuous feature observed in the majority of our sites for the P1 transect was that CH 4 concentrations noticeably elevated from the surface to the bottom water. This was particularly pronounced at CTD4 and CTD5, where CH 4 concentrations decreased from the surface to 600 m, then increased from 600 m to the seabed (Figure 3). Finally, they reached the maximum value of 20.3 nM for CTD4 and 10.1 nM for CTD5 at the bottom water, which are far more than that in the typical open ocean (Tilbrook and Karl, 1995;Watanabe et al., 1995).
The CH 4 concentrations of the P2 and P3 transects unveiled a noticeable shelf-to-slope trend, which increased and then decreased with distance from the shelf. On the continental shelf with shallow water depth, all sites were characterized by relatively high CH 4 concentration, with the value > 10 nM, especially in surface waters, which may be influenced by shelf mixed water (Ye et al., 2016;Sun et al., 2018). However, the maximum values in both transects appeared on the slope. In the P3 transect, the highest concentration was 12 times higher than the background value, which appeared at a water depth of approximately 900 m at the site of CTD10. Whereas, the maximum value of P2 appeared at a depth of 1,100 m (near bottom water) at the site of CTD21, with a concentration of 18.4 nM.
Moreover, the characteristics of CH 4 concentration in the P4 and P5 transects showed relatively high value in the shelf area and on the lower slope, while it was relatively low on middle slope. This scene can be observed at the CTD26 site with a water depth of 165 m on the shelf area, where CH 4 concentration in the whole water depth profile was relatively high, with the range from 8.6 to 10.8 nM. That might be associated with the input of high concentration CH 4 from land sources or the production and diffusion from sediment (Ye et al., 2016). However, in the P5 transect located at the canyon in the southern part of the OT, CH 4 concentration was lower than that in other sections and was close to the background value, which may be due to the complex geological structure or hydrography (Nakamura et al., 2013). Previous study showed that the internal waves and tides, combined with the Kuroshio current, enhanced the disturbance and mixing of sea water in submarine canyons (Gao, 1996), that resulted in a rapidly dilution of emitted CH 4 by surrounding waters. Thus, it was one of the important reasons for the relatively low CH 4 concentration in this area.
In the present research, we summarized the distribution characteristic of CH 4 in the P6 transect along the contour line on the west middle slope of the OT. Our results showed that the CH 4 concentration ranged from 1.9 to 24.7 nM, and the maximum and minimum values were at CTD10 in the middle OT and at CTD37 in the southern area (Figure 3), respectively. However, at the FIGURE 3 | The six transects distribution of CH 4 in the OT. The plots were generated using ODV Version 5.1.0 (Schlitzer, 2018). The "gray spaces" in P1 represent areas where there is too little information for ODV to interpolate.
CTD10 site, a high concentration of CH 4 was noted throughout the water column, which may be attributed to methane seepage and methane bubble plumes (Figure 1). Similarly, an abnormally high concentration of CH 4 was found at the depth water of the CTD4 site, with a distance of only 16 km from the CTD10 site. Moreover, the CH 4 concentrations in the near-surface waters of the CTD16 and CTD17 sites were relatively high, with the values of 12.5 and 14.5 nM, respectively.
According to the analysis results of the above-mentioned six transects, it was disclosed that higher CH 4 concentration was observed on the shelf and lower slope area within the depth range of 900-1000 m. In order to identify the CH 4 influence factors of the main water masses and currents, we presented the CTD profile data of all sites in Figure 4 and Supplementary Figure S1. Most data points indicated that the shelf mixed water (SMW) and Kuroshio (salinity > 34 , with gradual reduction of temperature with depth) played substantial roles in regulating the distribution of CH 4 concentration. The SMW distributed in the shelf edge zone included Changjiang Diluted Water (CDW), TWCW, and Kuroshio Surface Water (KSW) (Zhang et al., 2007;Ye et al., 2016). Previous studies reported that a small number of large rivers in East Asia (Sieburth, 1987) were transporting huge amounts of freshwater along with high concentrations of dissolved CH 4 to the ECS, e.g., about 70.6 × 10 6 mol yr −1 CH 4 , were imported from the Changjiang (112-190 nM in the surface water in spring) to the ECS (Tsurushima et al., 1996;Zhang et al., 2004;Ye et al., 2016;Sun et al., 2018). It is noteworthy that dissolved CH 4 in the river water is almost 2 orders of magnitude higher than that found in the continental shelf edge , and in summer, this river water mass extends further to the continental shelf, which may justify the relatively high CH 4 concentration in this area. However, the main stream of the Kuroshio flowing northeastward along the 200 m isoline (Qin, 1987), is the most important factor that affects the study area (Figure 4), with a relatively low CH 4 concentration (Rehder and Suess, 2001;Zhang et al., 2004). This is coupled with our finding in the slope area, and pronounced in the southern area. However, the CH 4 concentration is abnormally high at the sites of the lower slope (e.g., CTD4, 10) even though it is affected by the Kuroshio Intermediate Water (KIW) or Kuroshio Deep Water (KDW)  (Tang et al., 1997;Li et al., 2012;Qi et al., 2014)  ( Figure 4). This may be attributed to the observed methane cold seepage (Figure 1).
However, CH 4 concentration at the CTD10 and CTD4 sites was found anomalously high in the full depth water column, which was pronounced at a depth of 800-1,100 m near the seafloor. These irregular phenomena are closely related to the cold seepage with methane bubble plumes (Figure 1). CH 4 produced in sediments is consumed as it is pushed upward by anaerobic oxidation of methane (AOM) or aerobic CH 4 oxidation (Barnes and Goldberg, 1976). However, if the upward CH 4 is not fully exhausted by these processes, the extra CH 4 is emitted to the ocean, and generates patches of high CH 4 concentration (Leifer et al., 2000;Mau et al., 2017). Furthermore, the presence of gas hydrates is inferred from the bottom-simulating reflections (BSRs) in the OT (Fan and Yang, 2004;Xiwu et al., 2008). Geological models calculated that the water depth of the gas hydrate stability zone of the OT was approximately 500 m (Fan and Yang, 2004). Fittingly, the water depth of the CTD10 and CTD4 sites is about 900 m located in the gas hydrate stability zone. In addition, other evidence of pore-water geochemical anomalies Xu et al., 2018) and geochemical characteristics of authigenic pyrite associated with AOM  have been documented with cold seepages activities appearing near our study sites. Although CH 4 concentration is significantly lower than that in the samples collected by a ROV at the cold seepage sites here, previous studies still supported our augment that anomalously high CH 4 concentrations are attributed to the existence of cold seepage. However, the low concentration of CH 4 in our samples compared with that in the samples collected by a ROV was attributed to the greater distance above the seafloor during sampling. After emission from the seafloor, CH 4 concentration may be noticeably reduced due to continuous oxidation, dissolution, diffusion, and dilution of a large amount of seawater (Leifer et al., 2000;McGinnis et al., 2006;Konn et al., 2009;Faure et al., 2010). Not only that, when the CTD seawater was retrieved, CH 4 was partly lost due to a decrease of pressure (Konno et al., 2006). Furthermore, near-surface waters at the CTD16 and CTD17 sites showed high CH 4 concentrations, which were more than six times of the background value of the lower slope area, suggesting that there may be methane seepage. Nevertheless, compared with other areas, e.g., west Spitsbergen (up to 524 nM; Gentz et al., 2014), the Gulf of Mexico (∼600 nM; Solomon et al., 2009), and the Hikurangi Margin of New Zealand (up to 3,500 nM; Faure et al., 2010), it seemed that CH 4 concentration in this study was lower, and the CH 4 seepage was not considerably active in our study area. It is possible that CH 4 seepage is episodic, or the seepage is on a small scale or only very close (decimeters) to the seafloor at active seep sites (in the current study, sampled bottom water depth was about 50 m above seafloor). Therefore, it is difficult to catch the abnormality of methane with CTD sampling. Generally, cold seep in a finer scale is spatially associated with canyons, which incision downward into older strata and landward into the upper slope could physically disrupt gas hydrate or free gas deposits. On the contrary, no abnormal CH 4 concentration in P5 was detected in the canyon of the study area in spite of the existence of pore water anomalies associated with cold seep (Zhang et al., 2004;Xu et al., 2018). It is also possibly attributed to the episodic characteristic of CH 4 seepage or that emitted CH 4 is diluted swiftly, distributed by surrounding waters, or is actively consumed in oxygenated seawater.
Notably, our sampling site was relatively close to the hydrothermal vents of the OT (Figures 1, 2), thus, the hydrothermal activity may be an influence factor on the CH 4 distribution.

Vertical Distribution of CH 4
The vertical concentration profiles of CH 4 were irregular due to the complex geography and hydrography of the ECS (Zhang et al., 2004 ; Figures 3, 5). Nonetheless, CH 4 concentrations at nearly half of our study sites exhibited obvious peaks in surface or subsurface water. The largest maximum generally occurred at 50-100 m and coincided with the pycnocline. It was likely the result of in situ aerobic CH 4 production caused by microbial biogeochemistry (Reeburgh, 2007). At depths > 200 m, CH 4 peaks were present at some sites; for example, at sites CTD5, 13, and 25 (Figure 5), one broad CH 4 peak was present at the depth of ∼400 m, which might be associated with the advection and sinking of shelf water across the shelf break (Patra et al., 1998;Ye et al., 2016). Below the depth of ∼ 400 m, CH 4 concentration in the majority of profiles decreased with increasing depth and with a concurrent decrease in salinity, indicating the influence of the KIW (Figure 4). Finally, CH 4 concentrations increased again when approaching the sea bottom ( Figure 5). For instance, CH 4 concentrations in the water columns of CTD4,5,8,10,11,12,15,16,17,20,24,and 28 showed a steady increase to the seafloor, with the maximum value of higher than 10 nM. This may be attributed to CH 4 emissions from the underlying sediments discussed above. According to previous studies, biological production of CH 4 is expected to be widespread on the shelf due to a copious supply of terrestrial organic matter (Floodgate and Judd, 1992;Lin et al., 1992;Judd and Hovland, 2009). As this organic matter accumulates faster than the supply rate of oxidizing agents, methanogenesis sets in and causes organic carbon concentrations to be partly oxidized to CO 2 and partly reduced to CH 4 which supports the argument that CH 4 in bottom water is maintained by emission from the underlying organic-rich sediments. The water column increase in CH 4 concentration when approaching the seafloor was particularly pronounced at the CTD4 and CTD10 sites, with the maximum value of higher than 20 nM at the depth of several meters above the seafloor. This can be attributed to CH 4 seepage already reported in this area Xu et al., 2018).
Contrary to other sites, the CTD9 site showed positive CH 4 anomalies in the entire water column. The CTD10 site showed relatively high CH 4 concentration in the depths ranging from ∼450 to 1,000 m, while at CTD4, elevated CH 4 concentrations were observed in the depth range from ∼ 850 to 1,000 m. This may indicate that the depth range affected by the seabed CH 4 seepage, or the rising height of the methane plume (Figure 1). This is also consistent with the observations of the CH 4 bubble plume analyzed through multi-beam data (Sun et al., unpublished data), which revealed that the maximum height of the methane plume can reach 578 m above the seafloor near the CTD10 site where the water depth is about 1,000 m. In general, CH 4 can be transported in the form of dissolved or free gas in water (Faure et al., 2010). Bubble dissolution models tell us that a 9 mm diameter pure CH 4 bubble released into the water column from a depth of 1,000 m will be dissolved completely when it reaches a depth of 495 m, while a bubble of 15 mm will fully dissipate at 366 m (Greinert and McGinnis, 2009;McGinnis et al., 2006). Assuming that methane bubbles from cold seeps rise from 1,000 to 400 m depth and the bubble dissolution model is applicable to our study, we infer that the diameter of the CH 4 bubbles released into the water column must be between 9 and 15 mm.

Methane Lost to the Atmosphere
In this study, we calculated CH 4 saturations in the surface seawater (Table 1), which showed obvious spatial variation and was supersaturated throughout the study areas (Supplementary Figure S3), ranged from 108 to 1078%, with the high CH 4 saturation mainly noted at lower slope sites surrounding CTD10 and CTD4. This indicated that surface seawater is a net source of atmosphere CH 4 . Spatially, CH 4 saturation ranged from 465 to 975% in the shelf edge area, with an average of 676%. Whereas, in the upper slope area (water depth < 800 m), the saturation varied in a comparatively large range, from 110 to 606%. Different from the shelf and upper slope areas, saturation remarkably changed on the lower slope (108-1078%). Spatial variability among different sites was associated with the characteristics of their geological and physical environment, e.g., seep intensity, water depth, and currents.
To quantify the CH 4 lost to the atmosphere in this area, we estimated the sea-air flux. Overall, the sea-air flux of CH 4 throughout the study area presented heterogeneity in space (Figure 6), which ranged from 0.74 to 116 µmol m −2 d −1 , with an average of 26.2 µmol m −2 d −1 . Compared with the CH 4 seaair flux reported by Tsurushima et al. (1996) and Sun et al. (2018) in the shelf, our estimation showed a higher value, which may be due to different sampling and detection methods. In this study, we conducted CH 4 analyses immediately after sampling in a laboratory on-board, which made our detection closer to the actual value than measurements taken after months of store time. Furthermore, the range of the sea-air CH 4 flux on the continental shelf is from 34.8 to 104 µmol m −2 d −1 , with an average of 61.7 µmol m −2 d −1 , obviously higher than that on the global continental shelf (22-37 µ mol m −2 d −1) (Bange et al., 1994). Compared with other continental shelves worldwide, the sea-air flux was also higher than that in other continental shelf areas of the Western Pacific, such as the South China Sea (Tseng et al., 2017). However, an irregular wide range from 0.74 to 116 µmol m −2 d −1 on the slope was noted, with the maximum sea-air CH 4 flux value presented around CTD10, which may result from the direct influence of seafloor cold seepage studied previously. Globally, CH 4 spillovers on the lower slope are significantly higher than those in the Gulf of Cádiz (Ferrón et al., 2009), the Beihai (Bange et al., 1994), the Philippine Sea (Tseng et al., 2017),  Table 1). In summary, our study area is an active area of CH 4 emissions.

Possible Effects of Hydrothermal Activities
In addition to the contribution from cold seepage to CH 4 in the water of the OT, hydrothermal fluid was noted as an efficient carrier of CH 4 . Previous studies demonstrated that approximately more than 20 hydrothermal vents have been discovered since the 1980s in the OT (Halbach et al., 1989;Miyazaki et al., 2017). Significant enrichment of hydrothermal fluids of the OT in gas species is one of the geochemical features compared with those in the typical sediment-starved Mid Ocean Ridge (Sakai et al., 1990;Konno et al., 2006;Kawagucci et al., 2011;Ishibashi et al., 2014). Furthermore, anomalously high CH 4 concentrations were reported in the venting fluid of Iheya North (up to 7.6 mM; Kawagucci et al., 2011) in the  middle OT, and documented CH 4 concentration in hydrothermal fluids of the Minami Ensei nearest to our study sites was approximately 2.1 mM (Kawagucci et al., 2013). These fluids, with CH 4 concentrations generally a thousand times more than the background value, are a significant contributor to the seawater CH 4 (Glasby and Notsu, 2003).
In the current research, CH 4 concentration at different sites was anomalously high in the deep-water of the study area. As shown in Figures 1, 2, distribution of high CH 4 concentrations in the depth water partly appeared near known hydrothermal sites (Inagaki et al., 2006). For instance, in depth of 1050 mbsl (not the bottom water), the data from CTD12 showed a CH 4 concentration of 15.2 nM, and this high concentration may originate from the nearest seep site of CTD10 as discussed earlier.
However, it was found to be very close to the Minami Ensei Knoll hydrothermal vent with the high temperature fluid enriched in CH 4 . The distance measured between them is approximately only 19 km. Although Cowen et al. (2002) reported that CH 4 in hydrothermal plumes can be dramatically decayed with dispersed distances by the microbial oxidation process along with dilution of the plume with background waters; CH 4 evidence in the hydrothermal plume occurred in a distal station 15 km off the hydrothermal vents of the Juan de Fuca Ridge with a CH 4 concentration of 600 nM. However, CH 4 concentration in the Minami Ensei Knoll (∼2.1 mM) and the velocity of the current in the OT were significantly higher than that in the Juan de Fuca Ridge (Cowen et al., 2002;Nakamura et al., 2008;Kawagucci et al., 2013), thus, CH 4 in the hydrothermal plume could be further transferred. Therefore, we speculate that the CH 4 concentration in the CTD12 site may be influenced by the Minami Ensei Knoll hydrothermal fluid. In addition, CTD17 showed a high concentration of CH 4 in deep rather than bottom water, with the maximum value of higher than 10 nM, which may be due to the horizontal extra input by non-buoyant plumes (German and Syefried, 2014). Similar to CTD12, a hydrothermal vent was found near the site named Iheya Ridge (Figure 1) with an extremely high CH 4 concentration of 7.6 mM (Ishibashi et al., 1995;Kawagucci et al., 2011Kawagucci et al., , 2013. In this case, this might effect the CH 4 concentration of CTD17.

Significance of Methane Seepage for the CH 4 Budget
Methane seepage from seafloor sediments remains poorly quantified, although it is widespread on continental margins (Judd and Hovland, 2009;Boetius and Wenzhöfer, 2013). Most previous studies demonstrated that the contribution of the oceans to the global atmospheric CH 4 budget is minor (Schmale et al., 2005;Yvon-Lewis et al., 2011), with inputs originating from surface seawater, only occurring in regions of supersaturation. In contrast, new sites of seafloor CH 4 seepage are being discovered continually (Skarke et al., 2014), and recent studies on CH 4 fluxes at sites of seepage have been conducted. However, previous results based on traditional indirect sampling techniques and modeling suggested that bubble plumes emitted from cold seepages at depths greater than 200 m do not reach the surface mixed layer due to bubble dissolution and CH 4 oxidation (Valentine et al., 2001). This case is present in deep-water hydrocarbon seepages in offshore Svalbard (Graves et al., 2015) and the northern Gulf of Mexico (Hu et al., 2012), where seafloor makes an insignificant local contribution to the atmosphere due to CH 4 oxidation. Similar to these areas, the sea-air flux that we estimated from the west slope of the OT was relatively small. Based on the mean flux, the preliminarily estimated emission rate of CH 4 from our study area in the ECS was 2.9 × 10 −3 Tg yr −1 . This indicated that although the study area only accounted for 2.47% of the surface area in the ESC (∼7.7 × 10 5 km 2 ), it accounted for about 4.5% of CH 4 emissions in the ESC (6.4 × 10 −2 Tg yr −1 ) (Ye et al., 2016). Nevertheless, as a unique back-arc basin, the characteristic of CH 4 can be influenced by both the cold seepage and the pervasive hydrothermal fluids as discussed above. Our study area only represented 2.47% of the ESC and was away from the majority of hydrothermal sites; thus, the contribution of emitted CH 4 from sediment to atmosphere in the ESC was remarkably larger than that estimated previously (Ye et al., 2016).
Additionally, these CH 4 emissions from the seafloor extreme systems may have a significant influence on the carbon cycle. Methane seepage allows buried old or deep sourced carbon to be transported to the seafloor (Boetius and Wenzhöfer, 2013), which contributes allochthonous forms of methane carbon to water (Pohlman et al., 2011;Zhang et al., 2019). At the same time, part of the methane carbon is bio-absorbed and then redeposited in carbonate (Sun et al., 2015), which became the methane sink. On the other hand, Zachos et al. (2004) reported that large releases of methane could have contributed to ocean acidification and with that resulting in the dissolution of large amounts of sedimentary carbonate in historic events. In this case, that would speed up the circulation of methane carbon in the ocean.

CONCLUSION
CH 4 concentrations in the seawaters on the west slope of the OT were in the range of 2.7-24.7 nM, with an average of 6.5 nM. The special distributions revealed that CH 4 concentration in the majority of the transections showed a high to low trend from the shelf to the slope. However, CH 4 concentration in the P1 and P2 transections was subjected to high variability, showing high values in the bottom water of the CTD4 and CTD10 sites. The vertical distributions indicated that dissolved CH 4 can be influenced by in situ aerobic CH 4 production caused by microbial biogeochemistry or physical advective supply from the shallower depth water of the continental shelf.
CH 4 in the surface water was, in all cases, supersaturated with respect to the atmosphere. Besides, the sea-air flux of CH 4 throughout the study area showed high spatial variations, which ranged from 0.74 to 116 µmol m −2 d −1 , demonstrating that surface seawater is a net source of CH 4 to the atmosphere. Based on the mean CH 4 flux, the preliminary estimated rate of CH 4 from our study area in the ECS was 2.9 × 10 −3 Tg yr −1 , which was higher than the value speculated previously. Therefore, the previously estimated CH 4 emission rate in the ESC may be underestimated according to our findings. The coexisted contributions of the cold seepage and hydrothermal vents expand our view of the CH 4 budget of the back-arc basin worldwide.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

AUTHOR CONTRIBUTIONS
ZS and NW organized and designed sampling strategy. XiaZ, LW, XilZ, and BZ collected the sample. XY analyzed the CH 4 concentrations. XiaZ wrote the manuscript. CX, WG, and HC contributed to the data interpretation and discussion of the manuscript at different stages. All authors contributed to the article and approved the submitted version.

ACKNOWLEDGMENTS
We thank the captains and crew of Zhang Jian for their assistance in recovering the samples during the ECS marine geological Expedition during 2016 and 127 ROV cruises during 2017. We also thank Dr. Meng Wang and his team members from Shanghai Jiao Tong University for their help in exploration and sampling in cold seeping areas by the ROV Beaver.