Dynamics of O2 and pCO2 in a Southeast Asia seagrass meadow: Metabolic rates and carbon sink capacity

Dissolved oxygen (DO) and partial pressure of CO2 (pCO2) were measured at half-hourly intervals from June 29 to September 9, 2019, in a seagrass meadow in the Southeast Asia archipelagos region. The open water mass balance of the O2 approach was used to calculate metabolic rates (i.e., gross primary production (GPP), community respiration (CR), and net community production (NCP). The calculations show that GPP and CR rates in the seagrass meadow of Dongsha Island were approximately 2.5 times higher than the global means (GPP, 507 ± 173 vs. 225 ± 11 mmol O2 m-2 d-1; CR, 497 ± 171 vs. 188 ± 10 mmol O2 m-2 d-1), while NCP was similar to the global mean (8 ± 61 vs. 27 ± 6 mmol O2 m-2 d-1), suggesting that seagrass meadows with high GPP may not necessarily hold high potential for carbon sequestration. The current data set also reveal that NCP tended to increase with GPP only at lower GPP levels, while NCP did not increase with GPP anymore at higher GPP levels. Moreover, the autotrophic/heterotrophic status did not correspond well to the sink/source behavior of CO2, suggesting that organic carbon metabolism could not be the only dominant factor in determining the sink/source status in a typical seagrass meadow underlain by carbonate sediments, which was further supported by the observed decrease in the trend of pCO2 with a relatively stable NCP level over the study period. These results demonstrate that the metabolism and the relationship between NCP and pCO2 in the seagrass meadows of Dongsha Island may deviate greatly from the global mean condition. To obtain a better assessment of the global potential of seagrass meadows as a nature-based solution for carbon sequestration, more regional-specific studies are still needed in the key regions, such as Indonesia and the Pacific archipelagos, that support extensive seagrass meadows but have not been charted.


Introduction
Seagrass meadows are highly productive habitats in shallow coastal waters and rank among the most productive ecosystems (Duarte and Chiscano, 1999). Recently, they have received substantial attention as a 'blue carbon' sink due to their high rates of primary production and carbon burial in sediments (Duarte et al., 2010;Fourqurean et al., 2012). As such, seagrass conservation and restoration are considered potential climate change mitigation and adaptation strategies (Nellemann et al., 2009;McLeod et al., 2011;Duarte et al., 2013). Yet, the high production and carbon input also render the seagrass meadow as sites of elevated abundance and activity of heterotrophs, which support high community respiration rates (Middelburg et al., 2005). As a result, the status of seagrass meadows as a source or a sink of carbon may largely depend on the metabolic balance of the primary producers (including seagrass itself, epiphytic algae, and microphytobenthos) and the associated heterotrophs (including the pelagic and benthic compartments) in the community.
The community metabolic status in seagrass meadows can be described by net community production (NCP), referring to the difference between gross primary production (GPP) and community respiration (CR). Autotrophic meadows (NCP > 0) fix organic carbon in excess to local demand (GPP > CR) and may, therefore, either store or export organic carbon to adjacent communities, while heterotrophic meadows (NCP< 0) require additional organic carbon and that the community may be sustained by allochthonous organic carbon (Duarte et al., 2010). Determining the metabolic status of seagrass meadows is therefore crucial for a better understanding of their potential in "blue carbon" sink and thus in climate change mitigation (Champenois and Borges, 2021).
Seagrass meadows support a large range of metabolic rates from autotrophy (Barroń et al., 2006) to heterotrophy (Van Dam et al., 2019), which vary significantly among different community types and species along with seasonal variations (Duarte et al., 2010). In fact, earlier efforts at assessing seagrass metabolism existed a large imbalance in the geographic distributions of the data available. In particular, there is a considerable gap in reports of seagrass community metabolism in the Indo-Pacific, Africa, and South America (Duarte et al., 2013). To gain a better sink potential of seagrass meadows on a global scale, observations are required from the key regions where supporting extensive seagrass meadows have not been reported.
Despite being a hot spot for the diversity of seagrass species and habitat types, the seagrass meadows in Southeast Asia archipelagos have been poorly studied (Fortes et al., 2018), and to our knowledge, metabolic rates in this region have not yet been reported. Lee et al. (2017) have reported the community-scale metabolism of seagrass meadows in East Asia (Korea) for the first time as autotrophic state, but their study site was located in temperate area. In the present study, we report a data set of O 2 , partial pressure of CO 2 (pCO 2 ), and pH measurements at half-hourly intervals, which have been recorded continuously from June 29 to September 9, 2019, in a seagrass meadow on the northern shore of Dongsha Island, located in the largest tropical marginal sea in Southeast Asia (i.e., the South China Sea). Community metabolic rates (GPP, CR, and NCP) were estimated on a daily basis from this data set using the open water mass balance of the O 2 approach (Odum, 1956;Champenois and Borges, 2021). We also investigated the relationship between metabolism and concurrent concentrations of pCO 2 in the water column. Overall, the objectives of this study are (i) to document the metabolism rates and the relationships between them for the first time in the seagrass meadows in the Southeast Asia archipelagos region, (ii) to compare the present result with previously global synthesis, and (iii) to provide new insight into the relationship between metabolism and pCO 2 dynamics for tropical seagrass meadows.

The study site
Dongsha Island (also known as Pratas Island) is a coral island located on the western side of Dongsha Atoll, a ring-shaped coral reef ecosystem in the northern South China Sea ( Figure 1). The island is quite small around 2.80 km long and 0.87 km wide with an area of approximately 1.74 km 2 , and it is made up of coral atolls and reef flats. Vines and bushes cover some of the island, and the rest is surrounded by white coral sand (Dai, 2006). The East Asia Monsoon prevails over Dongsha Island with distinct seasonality: the southwesterly monsoon starts in June and lasts until August in the summer, while the northeasterly monsoon starts in October and predominates over the winter and early spring. The transition periods are April-May for the southwesterly monsoon and September for the northeasterly monsoon (Chai et al., 2009). As a military control zone, there is no permanent inhabitants on the island, and it is visited only by military personnel, researchers and fishermen. Furthermore, there is no freshwater flow into the sea from the island, and most of freshwater on the island is produced by seawater desalination and recycled. Therefore, the inflow from the artifact at Dongsha Island nearly could be negligible.
Despite the small area of Dongsha Island, massive seagrass meadows extend from the intertidal to the subtidal zones with a total coverage area of 11.85 km 2 around the entire island . A total of seven seagrass species from six genera and two families have been identified at Dongsha Island, including Thalassia hemprichii, Halophila ovalis, Cymodocea rotundata, C. serrulata, Halodule uninervis, Syringodium isoetifolium, and Thalassodentron ciliatum (Lin et al., 2005), and the coverage of seagrass meadows was generally >75% around the island (Lee et al., 2021). Among them, the most abundant species around Dongsha Island are T. hemprichii, H. uninervis, C. rotundata, and C. serrulata , and which are also popular species in Indo-Pacific (Short et al., 2007). Besides, T. hemprichii and H. uninervis are widely appeared in Hengchun Peninsula and surrounding islands of Taiwan, but C. rotundata, and C. serrulata are only recorded in Dongsha and Penghu Islands (Lin et al., 2005).This study was conducted in the seagrass meadow on the northern shore (NS, Figure 1) of Dongsha Island, which is a multispecies seagrass meadow with three dominant species of T. hemprichii, C. rotundata, and C. serrulata. The average total seagrass biomass, total seagrass production, seagrass cover, and shoot density were previously reported to be 949.4 ± 62.1 (g DW m −2 ), 8.47 ± 2.92 (g DW m −2 ), 81.91 ± 2.13 (%), and 2677 ± 485 (shoots m -2 ), respectively, based on four seasonal surveys in April, August, October 2010, and February 2011. The large area of the seagrass meadow, the high shoot density and coverage of seagrasses indicate that the seagrass on Dongsha Island is in a good heath condition (Lin et al., 2005). The epiphytic algae cover on seagrass leaves was previously reported to be approximately 4%, 5%, 9%, and 5% in April, June, August, and October 2010, respectively (Chang, 2012). The seasonal variations in salinity and chlorophyll a concentration in water column around the Dongsha Island were relatively minor, which were 33.8 ± 0.4, 33.8 ± 0.5, 33.8 ± 0.5, 34.6 ± 0.3 for salinity (Chou et al., 2021), and 0.81 ± 0.37, 0.70 ± 0.20, 0.60 ± 0.25, 0.41 ± 0.13 mg m -3 for Chlorophyll a (Liang, 2012) in spring, summer, autumn, and winter, respectively. Furthermore, only 4 rainy days were recorded during the entire study period with tiny precipitation less than 2 mm d -1 , suggesting that salinity could remain within a relatively stable level in the study site during the study period.

Dissolved oxygen, pCO 2 , pH, and light intensity measurements
High-resolution dissolved O 2 (DO), pCO 2 , and pH were measured simultaneously at intervals of 30 min, 1 h, and 30 min using HOBO U26 Dissolved Oxygen Data Logger (Onset, MA, USA) and Submersible Autonomous Moored Instrument sensors (SAMI-CO 2 and SAMI-pH; Sunburst Sensors, LLC, MT, USA), respectively. To avoid biofouling, SAMI-CO 2 and SAMI-pH are designed to isolate the sensor in an enclosed cell that is periodically flushed with sample that has been exposed to a biocide such as tributyl tin or copper, which has been proven to be useful in keeping biofouling to a minimum (Lai et al., 2018). Furthermore, the sensors using in this study were covered by copper mesh to protect from fish bites and other sources of damages (Lai et al., 2018). Before deploying, the DO logger was first calibrated to 100% saturation by placing it in watersaturated air and then to 0% saturation by placing it in sodium sulfite solution. The accuracy and resolution of DO measurement were 6 and 0.6 mM, respectively. SAMI-CO 2 and SAMI-pH sensors were returned to Sunburst Sensors LLC for re-calibration and refurbishment just before their deployments. The accuracy and resolution of pCO 2 and pH measurements were ±3 matm and<1 matm and ±0.003 and<0.001 pH units, respectively. The long-term drift of SAMI-CO 2 and SAMI-pH sensors was<1 matm and<0.001, respectively, over six months. Light intensity was recorded by a HOBO UA-002-64 Data Logger (Onset, MA, USA) at an interval of 30 min, and the measured lux was converted to photosynthetically active radiation (PAR, mmol/m 2 /s) by dividing lux by 52.63 (Langhans and Tibbitts, 1997). These sensors were deployed within the seagrass canopy 0.1 m above the bottom Map showing (A) the location of Dongsha Atoll, (B) the distribution of dominant seagrass (studies by Lin et al., 2005) and the study site on the northern shore of Dongsha Island (NS, asterisk), and (C) a photo of the sensors deployed within the seagrass canopy.
( Figure 1C), for a period of 73 days, from 6/29 to 9/9 2019. The water column depth is approximately 1 m at the study site. Furthermore, the wind, wave and tide data during the study period were collected by Central Weather Bureau of Taiwan (Buoy Observation Data Annual  Report 2019 and Tide Tables 2019). 2.3 Gross primary production, community respiration, and net community production calculations The metabolic rates (i.e., gross primary production, GPP; community respiration, CR; net community production, NCP) were computed using the open water mass balance of O 2 (Odum, 1956). First, the hourly respiration was calculated as the temporal excursion in DO during the nighttime hours (18:00-06:00) with corrections for air-water gas exchange using the following equation: where DDO Dt is the rate of change in DO between two consecutive measurements (0.5 h); h is the water depth (m); and air-water O 2 flux is the air-water exchange flux of O 2 , calculated with a bulk transfer approach using the formula of Wanninkhof (1992): where k 600 was the gas transfer velocity estimated using the parametrization of Ho et al. (2006), in which wind data was obtained from the meteorology station on Dongsha Island; S c was Schmidt numbers calculated from in situ salinity and temperature data (Wanninkhof, 1992); O 2(water) and O 2(air) were measured and saturated DO concentrations, respectively. Time-series variations in k 600 were shown in Figure S1, and an error analysis on the calculations of CR, GPP and NCP resulted from the uncertainty in k 600 were given in Figure S2. Likewise, the change in the rate of DO between two consecutive measurements corrected for the exchange of O 2 with the atmosphere during daytime (06:00-18:00) corresponds to the hourly rate of NCP. Assuming a constant hourly rate of CR through the diel cycle, the hourly rate of GPP was calculated as NCP-CR (Champenois and Borges, 2012). Daily integrated GPP was computed by multiplying the average hourly rate by 12 (i.e., the daytime hours); daily integrated CR was calculated by multiplying the average hourly rate by 24 (i.e., the daytime + nighttime hours); and daily integrated NCP was computed as the sum of daily integrated GPP and daily integrated CR.

Site characteristics: Wind, wave, and tide
The daily wind, wave and tide variations are showed in Figure 2. Daily mean wind speed ranged from 1.6 to 11 m s -1 with an average speed of 6.0 m s -1 during the study periods. Southwesterly wind direction was dominant before Julian day 232 and wind direction was variated after that day. These results reflect that southwesterly monsoon was prevailed during the summer (July and August) but ceased in the early autumn (September), which is consistent with the general pattern of the East Asia Monsoon. Daily mean wave was within a range of 0.56 -3.37 m with an average level of 1.9 m, and most of daily wave were southwesterly direction. Tidal range was from 0.04 m of neap tide to 1.42 m of spring tide with an average level of 0.73 m.
3.2 In-situ monitoring of temperature, light intensity, dissolved oxygen, pCO 2 , and pH The diel patterns of temperature, light intensity, DO, pCO 2 , and pH are shown in Figure 3. Over the study period, water temperature varied between 27.5 and 35.0°C and averaged 30.4 ± 1.4°C with a regular diurnal trend (i.e., increasing during daytime and decreasing during nighttime). The lowest temperature period occurred on Julian days 213-214 and 232-233, which corresponded well with the lowest PAR intensity period, suggesting light intensity may exert a strong control over water temperature variation. The variation ranges and the means ( ± SD) of DO, pCO 2 , and pH were 16-529 mM and 195 ± 102 mM, 39-1565 matm and 438 ± 252 matm, and 7.55-8.38 and 7.97 ± 0.16 pH units, respectively. Generally, DO, pCO 2 , and pH revealed distinct diurnal patterns, with DO and pH increasing but pCO 2 decreasing during the day and DO and pH decreasing but pCO 2 increasing during the night. As shown in Figure S3, both DO saturation and pCO 2 closely followed the diurnal cycle of photosynthesis and respiration, suggesting that organic carbon metabolism could play an important role in regulating the CO 2 dynamics in the seagrass meadows of Dongsha Island (Chou et al., 2018). Furthermore, similar to the temperature variation, the smallest amplitude of diel variation in DO corresponded to the lowest PAR intensity period, implying a possible important role of light intensity in regulating metabolic activities (i.e., photosynthesis and respiration) in seagrass meadows.
Provided that the average atmospheric pCO 2 with 100% humidity was 395 matm during the study period, seawater pCO 2 was lower and higher than the atmospheric pCO 2 in up to 52% and 48% of the time, respectively, indicating contrasting times when seagrass meadows could either act as a CO 2 sink or source for the atmosphere. Likewise, assuming that the pH value of the source water was 8.0, pH was lower and higher than the source water pH in up to 55% and 45% of the time, respectively, during the entire study period, demonstrating divergent times when seagrass meadows would either enhance or mitigate ocean acidification.
3.3 Daily variability of gross primary production, community respiration, and net community production The daily variability of GPP, CR, and NCP is presented in Figure 4. Over the study period, GPP was characterized by strong day-to-day variability, ranging between 82 and 832 mmol O 2 m −2 d −1 with an average of 507 ± 173 mmol O 2 m −2 d −1 . The lower GPP corresponded well with the lower PAR intensity, suggesting that GPP variations could be closely related to the daily variability of irradiance. CR ranged between 99 and 900 mmol O 2 m −2 d −1 with an average of 497 ± 171 mmol O 2 m −2 d −1 , and its variations closely tracked those of GPP, implying that high GPP could also support high CR. NCP also showed a strong day-to-day variability but did not track those of either GPP or CR. During the study period, NCP ranged between −138 and 197 mmol O 2 m −2 d −1 and averaged 8 ± 61 mmol O 2 m −2 d −1 , indicating the range of ecosystem metabolism from a distinctly autotrophic to a clearly heterotrophic status and might be near metabolic balance, in general.

Comparison of seagrass community metabolism between the present study and global synthesis
Based on a data set containing 403 individual estimates derived from a total of 155 different sites, Duarte et al. (2010) estimated the global mean of GPP, CR, and NCP (mean ± SE) in seagrass meadows to be 225 ± 11, 188 ± 10, and 27 ± 6 mmol O 2 m −2 d −1 , respectively, while these metabolic rates were 507 ± 173, 497 ± 171, and 8 ± 61 mmol O 2 m −2 d −1 , respectively, in the present study (Table 1). This Frontiers in Marine Science frontiersin.org comparison clearly demonstrates that the averaged GPP and CR in the seagrass meadows on the northern shore of Dongsha Island were nearly 2.5 times higher than the global means, in contrast, the average NCP was similar to the global mean. This finding may not be surprising because the previous syntheses have shown that GPP, CR, and NCP reveal significant differences across different biogeographic locations and that tropical seagrass meadows tended to support higher metabolic rates, with somewhat lower NCP than temperate ones (Duarte et al., 2010). As Dongsha Island is located in the largest tropical marginal sea (the South China Sea) in the northwest Pacific, higher GPP and CR but lower NCP could be expected. Additionally, it is worth noting that the average total biomass of seagrass meadows on Dongsha Island (947-1047 g DW m -2 , Huang et al., 2015) is within the high end of the global dataset (716 ± 140 g DW m -2 , Fourqurean et al., 2012; Table 1). Furthermore, previous study has found that Thalassia hemprichii has the highest specific GPP to seagrass biomass ratio for meadows dominated by different species, including Cymodocea nodosa, Cymodocea rotundata, Enhalus acoroides, Halodule uninervis, Halodule wrightii, Halophila ovalis, Posidonia oceanica, Ruppia maritime, Syringodium filiforme, Syringodium isoetifolium, Thalassia hemprichii, Thalassia testudinum and Zostera marina (Duarte et al., 2010), and Thalassia hemprichii is one of the dominant species in the seagrass meadows on the northern shore of Dongsha Island . Therefore, the high biomass and high species-specific GPP to biomass ratio may partially explain the observed high GPP in the present study.
With the increasing recognition of their importance for carbon sequestration, seagrass conservation and restoration have been widely proposed as blue carbon strategies to mitigate and adapt to climate change (Nellemann et al., 2009;McLeod et al., 2011;Macreadie et al., 2019). Nevertheless, the present study suggests that seagrass meadows with high GPP in the tropical northwest Pacific do not necessarily hold a high potential for carbon sequestration due to the lower NCP. Accordingly, to obtain a better assessment of the potential of seagrass meadows as a nature-based solution for carbon sequestration, more regional-specific metabolic studies are imperatively needed before implementing any seagrass conservation and restoration plans, particularly in the northwestern Pacific archipelagos, where are hot spots of seagrass meadows but with very few record of community metabolism. 4.2 Relationships between community respiration, net community production, and gross primary production A strong quadratic relationship was observed between CR and GPP (CR = 129 + 0.42xGPP + 0.0006xGPP 2 , r=0.94; Figure 5A), indicating that high production rates may support high community respiration rates. In seagrass meadows, CR comprises both autotrophic respiration (AR) and heterotrophic respiration (HR). Generally, AR is linked to GPP at the time scale of minutes to hours, whereas HR is linked to GPP at the time scale of hours to days for bacteria and days to weeks for metazoans (Del Giorgio and Williams, 2005;Champenois and Borges, 2021). We observed a strong coupling between daily CR and GPP ( Figure 5A), thus suggesting that HR from metazoans is less likely the dominant contributor to the observed strong coupling of GPP and CR in seagrass meadows of Dongsha Island. In contrast, seagrasses could release a substantial amount of DOC (Barroń and Duarte, 2009), which can support bacterial HR, so enhanced HR in response to DOC release associated with high GPP might contribute to the close coupling of GPP and CR. It is then likely that the tight coupling between CR and GPP is mainly due to a close coupling between AR and GPP and/or HR for bacteria and GPP. Furthermore, the Time-series variations in (A) temperature, (B) PAR, (C) DO and DO saturation level (DO%), (D) pCO 2 (red solid line) and temperature normalized pCO 2 (TpCO 2 , blue dashed line), and (E) pH (red solid line) and temperature normalized pH (TpH, blue dashed line). TpCO 2 and TpH were normalized to a constant temperature of 30.4°C, the average temperature during the study period.
quadratic increasing trend between CR and GPP suggest that the rate of increase in CR with GPP would be lower at lower GPP than that at higher GPP. The different CR increasing rates at various GPP levels implies that the individual contribution of AR and HR to CR could respond differently to the variation in GPP. In fact, the separation of the individual contribution of AR and HR to CR remains a challenge to date in the metabolism studies of seagrass meadows (Champenois and Borges, 2021). The current results further suggest that how the Daily (A) community respiration (CR), (B) gross primary production (GPP), and (C) net community production (NCP) from June 29 to September 9, 2019 (Julian day 180-252) over a seagrass meadow on the northern shore of Dongsha Island. The superimposed red line represents the daily photosynthetically active radiation (PAR).
TABLE 1 Comparison of gross primary production (GPP), community respiration (CR), net community production (NCP), and seagrass biomass in dry weight units of Dongsha Island (this study) with those for global, temperate, and tropical seagrass meadows as studied by Duarte et al. (2010) for GPP,CR and NCP, and Fourqurean et al. (2012)  AR and HR may respond differently to GPP variations is also an open question that needs to be carefully examined for a better understanding of the relationship between CR and GPP. Unlike the strong correlation between CR and GPP, variability was observed in the relationship between NCP and GPP ( Figure 5B). Although previous studies have suggested that NCP tends to increase with GPP (Duarte and Agustı, 1998;Champenois and Borges, 2021), the present dataset revealed that the increasing trend of NCP with GPP only held true at lower GPP (GPP<~500 mmol O 2 m -2 d -1 ), and NCP did not increase with GPP anymore at higher GPP (GPP>~500 mmol O 2 m -2 d -1 ; Figure 5B). Interestingly, the lower rate of increase in CR at low GPP but a higher rate of increase in CR at high GPP can reasonably explain the observed divergent relationships between NCP and GPP at different GPP levels.
An average threshold GPP of 186 mmol O 2 m -2 d -1 on a global scale was proposed by Duarte et al. (2010), above which the seagrass meadows tended to be autotrophic (NCP>0) and shifted to heterotrophy (NCP<0) at lower GPP values. The average GPP in the current study was higher than 186 mmol O 2 m -2 d -1 , and NCP was positive and thus qualitatively in agreement with the global assessment. In spite of this agreement, the present dataset revealed a quite different GPP threshold from the global average; the present results show that seagrass meadows in the study area tended to be heterotrophic when GPP< 347 and GPP > 620 mmol O 2 m -2 d -1 , but they tended to be autotropic when GPP between 347 and 620 mmol O 2 m -2 d -1 ( Figure 5B). These results demonstrate again that the metabolic relationship in the seagrass meadows of Dongsha Island may deviate greatly from the global mean condition reported in previous synthesis studies, and thus more regional-specific studies are still needed in the key regions where support extensive seagrass meadows but have not yet been charted, such as Indonesia and Pacific archipelagos. The relationship between (A) community respiration (CR) and gross community production (GPP) and (B) net community production (NCP) and GPP in seagrass meadows of Dongsha Island.

Relationship between net community production and pCO 2
Although NCP and daily average pCO 2 showed an expected negative correlation (Champenois and Borges, 2021; Figure 6A), the autotrophic/heterotrophic status did not correspond well to the sink/ source behavior of CO 2 in the seagrass meadow of Dongsha Island on a daily scale. As shown in Figure 6A, when NCP<0 (heterotrophic), pCO 2 could fluctuate from a strong sink status (pCO 2 < 200 matm) to a strong source status (pCO 2 > 900 matm), providing air pCO 2 = 395 matm. Similarly, when NCP>0 (autotrophic), pCO 2 could also demonstrate either a sink or a source status. These results suggest that the organic carbon metabolism could not be the only dominant factor in determining the sink/source status of seagrass meadows (Macreadie et al., 2017;Saderne et al., 2019). This proposition can be further supported by the observed divergent temporal trends of NCP and pCO 2 in the present dataset. As shown in Figure 6B, NCP did not show either a decrease or an increasing trend over the study period, while pCO 2 revealed a significant decreasing trend with time.
Several recent studies have highlighted that in addition to organic metabolism, inorganic carbon processing (i.e., calcium carbonate precipitation/dissolution; Macreadie et al., 2017;Saderne et al., 2019) and the external sources of land-based pollutants (Banerjee et al., 2018) may also play a pivotal role in regulating the capacity of CO 2 sequestration of seagrass meadows. Since Dongsha Island is a military control zone, there is no permanent inhabitants neither no freshwater flow into the sea from the island, the impact of land-based pollution sources could be negligible. In terms of carbonate precipitation/ dissolution, Van Dam et al. (2021) recently reported that calcification-induced CO 2 emissions may exceed organic carbon sequestration in the seagrass-dominated estuary in central Florida Bay, and they argued that net carbon sequestration of seagrasses might be overestimated if calcification-induced CO 2 emissions are not accounted for. On the contrary, our previous studies around Dongsha Island revealed that the enhancement of total alkalinity (TA) coupled with carbonate dissolution may enhance the CO 2 sequestration capacity in seagrass meadows (Chou et al., 2018;Chou et al., 2021). We suggest that the lower pCO 2 in late August and early September (Julian day 230-250) compared with that in late June and early July (Julian day 180-200) could be associated with the favorable condition for TA accumulation in water column deriving from sedimentary carbonate dissolution and/or anaerobic metabolism. During the monsoon transition period, the weak wind results in less energetic hydrodynamics and longer residence time for the waters around Dongsha Island, making it favorable for the occurrence of sedimentary metabolic carbonate dissolution and/or anaerobic TA generation and the subsequent TA accumulation in the overlying waters (Chou et al., 2018;Chou et al., 2021). The elevated TA may lead to a lower pCO 2 level at a relatively stable NCP level. Unfortunately, TA data were unavailable during the monsoon transition period for northeasterly in September, but Chou et al. (2021) indeed found significantly higher TA during the monsoon transition period for southwesterly in April than in summer. This result implies that the different hydrodynamic regimes may lead to divergent sedimentary organic and inorganic metabolic processes, which may exert an important influence on regulating the carbon sequestration capacity of seagrass meadows. Therefore, the focus on organic carbon metabolism alone may not be enough, and hydrodynamics and calcium carbonate cycling remain key uncertainties that will need to be addressed to more accurately estimate the net carbon sequestration potential of seagrass meadows, especially for seagrass meadows with underlying carbonate sediments.
In summary, the present study is the first to report the metabolic rates of seagrass meadows in the Southeast Asia region. The results show that although daily NCP varied largely from autotrophy to heterotrophy, the average NCP over the entire study period was nearly in balance. Furthermore, the calculated GPP and CR in the present study were significantly higher, but NCP was similar to the global means, suggesting that seagrass meadows with higher GPP may not necessarily hold higher blue carbon potential. We also found that organic carbon metabolism was not the only dominant factor in determining the sink/source status in seagrass meadows of Dongsha Island, and hydrodynamics and calcium carbonate cycling may exert an important influence on regulating the carbon sequestration capacity of carbonate seagrass meadows. Since the present work was done mainly during the southwesterly monsoon period, more studies with better temporal coverage, in particular during the northeasterly monsoon and monsoon transition period are still needed to elucidate the linkage between monsoon-driven change in hydrodynamic and biogeochemical processes controlling CO 2 variation in seagrass meadow on Dongsha Island. Overall, the current results in the Southeast Asia region deviate greatly from the global mean condition, and thus regional-specific studies are imperatively needed before implementing any seagrass conservation and restoration as a nature-based solution for carbon sequestration.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions
W-CC conceived the study design, performed the field work, and wrote the manuscript. L-FF, C-CH, Y-YS, W-JH, H-KL, and T-YC contributed to the data collection and interpretation. All authors contributed to the article and approved the submitted version.