Plankton Metabolism in Coastal Waters of the Guangdong-Hong Kong-Macao Greater Bay: Regional Variance and Driving Factors

Metabolisms of field plankton community, including gross primary production (GPP), community respiration (CR), and net community production (NCP), usually indicate the status of the health, resource production, and carbon budget of marine ecosystems. In this study, we explored the regional variance and driving forces of plankton metabolism in coastal waters of the Guangdong-Hong Kong-Macao Greater Bay Area (GGBA), a fast-developed area with complex hydrological and environmental states. The results showed that the maximum GPP and CR occurred in the estuarine plume of the GGBA in summer, while in winter the more active metabolisms of plankton community occurred in the Daya Bay, with the GPP and CR being respectively mediated by the nutrient level and temperature. Moreover, four regional zones were divided on the base of the environments and biological factors in surface water of the GGBA i.e., the river-runoff zone, river-plume zone, nearshore and far-offshore zones. The metabolic states in these zones varied significantly due to the regional and seasonal variations of, for example, the nutrient level, temperature, and turbidity driven by multiple factors including land-derived runoffs, anthropogenic activities, the Yuedong Coastal Current, and offshore seawater-intrusions. On the whole, the GGBA areas exhibited the weak heterotrophic processes in both summer (NCP = -24.9 ± 26.7 mg C m-3 d-1) and winter (NCP = -51.2 ± 8.51 mg C m -3 d-1). In addition, we found that the higher CR occurred to the bottom layers of the river plume and nearshore zones wherein hypoxia happened, indicating a possible attribution of plankton community respiration to the hypoxia in the GGBA.


INTRODUCTION
Rising atmospheric CO 2 causes global climate change, which brings unprecedented challenges to human's sustainable development. As an important carbon reservoir on Earth, the ocean has been estimated to absorb 31 ± 4% of anthropogenic CO 2 from 1994 to 2007 (Gruber et al., 2019), wherein the plankton metabolisms mediate the exchange of CO 2 between air and sea (Odum, 1956;Falkowski, 1998;Duarte et al., 2013). The dynamics of gross primary production (GPP) and community respiration (CR) can directly reflect the carbon metabolism balance in marine ecosystems, as indicated by the net community production (NCP) (Duarte and Regaudie-de-Gioux, 2009;Cai, 2011). When the GPP is over the CR (i.e., NCP>0), marine ecosystem has a potential to absorb CO 2 and export organic matters and is thus autotrophic. When NCP <0, the marine ecosystem emits CO 2 and becomes a carbon source and is thus heterotrophic, and when NCP ≈ 0, the CO 2 absorption and emission tends to be in equilibrium (Duarte and Regaudie-de-Gioux, 2009;Caffrey et al., 2014). The carbon metabolism of the plankton community that often regulates the distribution and variation of carbon flux in marine ecosystem usually influences the evaluations of the carbon cycle and nutrient dynamics in the oceans.
The metabolism status of the marine plankton community is susceptible to multiple factors and shows greatly spatial and temporal changes (Karl et al., 2003). Therefore, it is not surprising that the controversial results often occur among different regions and seasons depended on whether the marine ecosystem is autotrophic or heterotrophic (Ducklow and Doney, 2013;Serret et al., 2015). A famous controversy is about the judgment of autotrophy or heterotrophy in oligotrophic oceans Williams et al., 2013). Such a debate of metabolic status also occurs in estuaries or bays. Many studies suggested that the estuarine and nearshore areas are heterotrophic and are the sources of CO 2 due to the landderived inputs of excess organic carbon, while the continental shelves are autotrophic and the sinks of CO 2 (Heip et al., 1995;Caffrey, 2004;Borges et al., 2005;Duarte and Prairie, 2005;Chen et al., 2013;Cloern et al., 2014). However, some recent studies showed that some estuaries, coastal and coral reef areas are autotrophic (Bauer et al., 2013;McKinnon et al., 2017;Agusti et al., 2018;Garcia Corral et al., 2020). Such a controversy demonstrates that the metabolic-balanced characteristics, and potential driving mechanisms need to be studied further in combination with specific biological and environmental factors because the regional estuarine and coastal ecosystems are usually affected by various human activities and environmental forces with a spatial variance.
As one of major urban agglomerations in the world, the Guangdong-Hong Kong-Macao Greater Bay Area (GGBA or Greater Bay Area) is one of the regions with the strongest economic vitality and fastest development in China. The coastal area of the GGBA covers the Pearl River Estuary (PRE) and its adjacent waters (Figure 1) and forms a coastal ecosystem with diverse habitats (estuary, bay, and offshore waters) and high biodiversity (Yang et al., 2019;Yu, 2019;Wang et al., 2021). Due to the inputs of nutrients and organic matters, the estuarine and nearshore waters are important junction of connecting land to open sea with active biochemical geochemical cycles (Gattuso et al., 1998). Previous studies showed that the greatly regional variances often occurred in physical and chemical environments of the PRE and its adjacent waters due to the Pearl River runoffs, human activities, and other physical processes (such as coastal current and coastal upwelling), which further affected the carbon fixation process of plankton assemblages (Cai et al., 2002;Song et al., 2004;Song et al., 2011). Many studies also showed that phytoplankton primary production and microbial metabolism mainly mediated the airsea carbon flux in the PRE (Cai et al., 2004;Guo et al., 2009;Dai et al., 2014;Li et al., 2019), suggesting the importance of plankton metabolisms in this ecosystem. However, there are few studies on the couplings of the respiration and primary production of the plankton community in coastal waters of the GGBA, therefore, the balance of autotrophic and heterotrophic processes and their variations and regulation mechanisms are still unclear as the environments therein are very complicated.
With rapid urbanization and economic growth, eutrophication, algal bloom, and hypoxia have become serious ecological and environmental problems in coastal waters of the GGBA (Yin et al., 2004a;Harrison et al., 2008;Song et al., 2009;Qian et al., 2018). As both algal bloom and hypoxia are closely related to the processes of plankton production and metabolism, studies on the balance between them and the influencing factors will help us to deeply understand the potential formation mechanism of these environmental and ecological events. Therefore, the GPP, CR, and NCP of the plankton community, as well as their associated environmental and biological factors were characterized in coastal waters of the GGBA during both wet (summer) and dry (winter) seasons, to clarify 1) the metabolic balance and carbon budget of plankton community in the GGBA and how the physicalenvironmental-biological factors shape their spatial and temporal variations; and 2) the relationship of plankton metabolic processes to the ecological risks of e.g., algal blooms and hypoxia events, as well as its potential indicative functions.

Study Area and Sample Collection
In this study, we carried out the experiments in 3 transactions located in the PRE (A1-A14), Daya Bay (B1-B10), and Huangmaohai Bay (C1-C10) during wet (summer; July 14-24, 2019); and dry seasons (winter; December 8-18, 2020) ( Figure 1). During the two investigations, the profiles of temperature, salinity, dissolved oxygen (DO), pH, turbidity, nutrients, and chlorophyll a (Chl a) concentration were obtained at all the 34 stations; and the cell size-fractionated Chl a, GPP, CR, and NCP in surface and bottom layers were measured at 20 of the 34 stations (labeled in red color). The seawater in the surface (~0.5 m depth), 5, 10, 15, 20, 30, 50 m, and bottom layers were collected with an acid-cleaned (1N HCl) 5-L polycarbonate bucket for Chl a and nutrient analysis as described below.

Environments and Chlorophyll a Measurements
At each station, the temperature, salinity, and DO profiles were measured with a multi-parameter water quality monitor Sonde (YSI 6600;Yellow Springs Instruments,Yellow Springs,OH,USA).
For nutrient analysis, the collected seawater from each depth was filtrated through 0.7 mm pore-sized glass fiber filters (25 mm, Whatman GF/F). The filtrate was then dispensed into an 80 ml polycarbonate bottle, frozen immediately, and stored at -20°C for subsequent analysis. Concentrations of NH + 4 , NO − 2 , NO − 3 , PO 3− 4 (DIP), and SiO 2− 3 were measured with an automatic nutrient analyzer (AA3; Seal Instruments, Hamburg, Germany) using the traditional colorimetry methods (Hansen and Koroleff, 1999). In this study, the DIN concentration was calculated as the sum of NH + 4 , NO − 2 , and NO − 3 concentrations. For the Chl a measurement, 250 ml of collected seawater was prefiltered through a 200 mm mesh and then sequentially filtrated through a 20 mm pore-sized nylon-net filter (25 mm, Millipore), 3 mm pore-sized polycarbonate filter (25 mm, Millipore), and 0.7 mm pore-sized glass fiber filter (25 mm, Whatman GF/F) to collect micro-(>20 mm), nano-(3-20 mm), and pico cell-sized phytoplankton cells (<3 mm), respectively. The total Chl a was obtained by directly filtrating the 200 mm pore-sized meshprefiltered seawater onto a 0.7 mm pore-sized GF/F filter. All the filters with phytoplankton cells were stored at -20°C for later analysis. The Chl a concentration was measured fluorometrically using a Turner Designs 10 fluorometer (Parsons et al., 1984) after extracting in the dark for 24 h at 4°C in 10 ml acetone (90% v/v).

NCP, CR, and GPP
To measure NCP and CR, the collected seawater was gently siphoned into 6 narrow-mouth 250 ml borosilicate glass Winkler bottles. Approximately 3 of the 6 bottles were wrapped in aluminum foil for dark incubation, and the remaining 3 replicates were incubated with the neutral cloth covered to achieve an in situ light level. The in situ light conditions were measured with the Secci disk method . After recording the initial DO level with a high-precision fluorescent oxygen probe (Microx 4; PreSens, Regensburg, Germany), all the Winkler bottles were incubated for 24 h in a water tanker with running surface seawater to control temperature. After the incubations, the final DO concentration in each bottle was measured again, and the NCP and CR were calculated as

Statistical Analysis
Before doing analysis, the data were tested for normality (Shapiro-Wilk test) and the homogeneity of variance (Levene's test). Then, the independent t test or Mann-Whitney test was used to identify the seasonal differences of environmental and metabolic factors. Principal component analysis (PCA) for environmental and biological factors was performed with CANOCO 5.0 software (Lepšand Smilauer, 2003). Data were standardized to minimize the influence of the variances of variables, and the individual transformed observations were named scores (Wu et al., 2016). The GGBA was divided into several zones according to the scores of the first two principal components of each station as well as their actual geographical locations (Song et al., 2011); then, the regional differences of environments and metabolisms were identified between different zones. Tukey's honestly significant difference (HSD) test was used to verify the significance of regional differences of environmental and biological factors, and the results were visualized using "ggplot2" package in R v4.0.5 (Ginestet, 2011;Zhao et al., 2021). Spearman correlation was used to determine the correlation of environment factors to plankton community metabolisms, and the results were visualized using the "ggcor" package in R v4.0.5 (Huang et al., 2020). An optimal multiple regression model was used to evaluate the contribution of environmental factors to the variation of plankton community metabolisms in R v4.0.5 (Grömping, 2006;Field et al., 2012). To estimate the importance of the variables, percentage increases in the mean squared error (MSE) were used: higher MSE% values indicate more important variables, according to the result of Spearman correlation analysis and the optimal multiple regression model with variance decomposition analysis in R v4.0.5 (Jiao et al., 2020). The map of sampling stations, spatial distribution, and transection plots was drawn using Ocean Data View (Schlitzer, 2020).

Environmental Characteristics
In summer, the seawater stratification generally existed in the GGBA according to vertical distributions of salinity and temperature, and all the investigated areas in the PRE, Daya Bay, and Huangmaohai Bay showed clear stratification, especially in nearshore waters (Figures 2A, B). The isohaline of 20 psu extended to the distances of over 100 and 40 km in transections A and C, respectively ( Figure 2B-1, -5), and the salinity was higher than 29 psu in transection B ( Figure 2B-3).
High DO and Chl a concentration generally appeared in the surface water , as well as in the offshore areas of transections A and B. Hypoxia events occurred in bottom of all three transections, with the lowest DO level of less than 2.2 mg L -1 in transection A ( Figure 2C-1), less than 3.3 mg L -1 in transection B ( Figure 2C-3), and less than 1.7 mg L -1 in transection C, respectively ( Figure 2C-5). Moreover, the surface water with DIN levels of over 25 mmol L -1 extended to the distances of over 100 and 50 km far from the land in transections A and C, respectively ( Figure 2E-1, -5). In winter, the water column was well mixed than that in summer. The isohaline of 20 psu shrank to the distances of less than 20 and 10 km far from the land in transections A and C, respectively ( Figure 2B-2, -6).
Consistently, the DO concentration distributed evenly in vertical scale and no hypoxia happened in the bottom layer, and Chl a concentration was lower than summer (Mann-Whitney test, p < 0.05) especially in nearshore areas ( Figure 2D). In addition, the surface water with higher DIN level (>25 mmol L -1 ) in winter just extended to~50 and~30 km far from the land in transections A and C, respectively ( Figure 2E-2, -6).
In general, nano-phytoplankton (3-20 µm) dominated (~60%) both the surface and bottom layers of nearshore waters of all three transections during summer and winter periods, while the pico-phytoplankton (<3 µm) dominated (~50%) both layers of the offshore waters. Moreover, contribution of nano-Chl a to total Chl a decreased offshoreward, while that of pico-Chl a increased. The contribution of micro-Chl a (>20 µm) to total Chl a (i.e., varying from 3.17% to 38%) showed no clear spatial change from nearshore to offshore waters ( Figure 3).

Plankton Community Metabolisms
Coinciding with Chl a biomass, the GPP varied from 52.3 to 661 mg C m -3 d -1 in surface water during the summer period and decreased from nearshore to offshore waters ( Figure 4A-1). The GPP reached as high as 661 mg C m -3 d -1 in estuarine-offshore waters, and such high GPP waters extended from southeast to northwest of the investigated area. In winter, however, the GPP ranged from 19.2 to 303 mg C m -3 d -1 and was significantly lower than that in summer (Mann-Whitney test, p < 0.05). The highest GPP (i.e., 303 mg C m -3 d -1 ) was observed in nearshore waters of the Daya Bay ( Figure 4A-2). Moreover, the GPP in the bottom layer was very low in both summer and winter periods (i.e., summer, 19.5 ± 4.5 mg C m -3 d -1 ; winter, 8.9 ± 1.5 mg C m -3 d -1 ), and displayed less spatial variation than that in surface water ( Figure 4A-3, -4).
Similarly, surface CR varied from 65.0 to 1,213 mg C m -3 d -1 in summer, and the active CR co-occurred with high GPP and Chl a, e.g., the CR reached up to 1,213 mg C m -3 d -1 in the  interface of plume-and saline-waters, where the GPP and Chl a reached as high as 659 mg C m -3 d -1 and 21.7 mg L -1 , respectively ( Figure 4B-1). In winter, the surface CR varied from 20.0 to 210 mg C m -3 d -1 ( Figure 4B-2). As compared to surface water, the CR in the bottom layer was markedly lower in summer (Mann-Whitney test, p < 0.05); however, such a difference did not occur in winter ( Figure 4B-3, -4) and no significant difference occurred between summer and winter periods either (i.e., 131 ± 79.9 vs. 110 ± 51.4 mg C m -3 d -1 ).
The NCP in surface water ranged from -554 to 377 mg C m -3 d -1 in summer and showed a similar spatial change as the GPP and CR ( Figure 4C-1). The highest NCP presented in inner PRE, while the lowest appeared in the interface of plume-and saline-water. In winter, surface NCP varied from -113 to 141 mg C m -3 d -1 and was lower than that in summer (Mann-Whitney test, p < 0.05) ( Figure 4C-2). The NCP in the bottom layer appeared to be negative in both summer and winter periods ( Figure 4C-3, -4), with the lowest values of -245 mg C m -3 d -1 in summer and of -241 mg C m -3 d -1 in winter presented in the Daya Bay. Moreover, the surface water of the GGBA was generally autotrophic in summer , and was close to the metabolic balance in winter; while the bottom water was heterotrophic in both investigated periods ( Figure 4C).

Principal Component Analysis Results on Regional Variances
Based on the PCA results in surface water as well as the geographical locations, the investigated area of the GGBA was divided into four regional zones (Figures 5A, B). In summer, Zone I that is located in the PRE had lower salinity and higher DIP, DIN, and SiO 2− 3 concentrations ( Figures 6B, D Table S1). Moreover, pico-Chl a allocation in Zone III was higher than Zones I and II but lower than Zone IV ( Figure 6I-1). The proportion of nano-phytoplankton in Zone IV was significantly lower, as compared to other zones ( Figure 6H-1). In addition, Chl a, GPP, and CR in Zone IV were lowest among all zones of the GGBA (Figures 6F, J, K; Table S1).
In winter, the surface water of Zone I is located near the PRE and Huangmaohai Bay and covered a narrower area than summer, and was closer to the upper reaches of the PRE. The salinity in Zone I was significantly lower than that in other zones ( Figure 6B-2), while the DIP, DIN, and SiO 2− 3 concentrations were significantly higher (Figures 6D, E; Table S1) but with moderate Chl a and lower GPP and NCP ( Figures 6F, J, L). Zone II covered the interface area of fresh-and saline-water and was closer to the coast, as compared to summer. The salinity in Zone II was lower than that of Zone III or IV, and the DIN was higher than that of Zone III and IV; however, the DIP was lower, leading to the higher N:P ratio (~82) than that in Zone III (~13) or IV (~31) (Figures 6B, D, E; Table S1). Surface Chl a in Zone II appeared to be moderate, just higher than Zone IV but lower than Zone III ( Figure 6F-2). The GPP and NCP in Zone II were lower either, while the CR was similar to the other zones ( Figures 6J, K, L). As compared to summer, Zone III in winter just covered the  nearshore area of the Daya Bay, wherein the salinity, temperature, and DIP were higher than Zone II or IV ( Figures 6A, B, E). The Chl a, GPP, and NCP in Zone III were significantly higher than other zones ( Figures 6F, J, L). Zone IV that covered a wider range of the offshore area in winter than summer had higher temperature and salinity than other zones but lower DIP, DIN and SiO 2− 3 concentrations. Moreover, Zone IV had the lowest Chl a, GPP, and NCP among all zones in winter ( Figure 6; Table S1).
In contrast to surface water, the bottom water that generally had higher salinity and lower temperature ( Figures 6A, B), as well as lower Chl a, GPP, CR, and NCP ( Figures 6F, J, K, L), was divided into three zones in both summer and winter periods ( Figures 5C,  D). In summer, Zone I that is located in the inner parts of the PRE and Huangmaohai Bay had higher temperature and nutrients, as well as lower salinity ( Figures 6A, B, D, E; Table S1). Zone II covered the inner part of the Daya Bay and the interface area of fresh-and saline-water, wherein it had the lowest DO among all three zones ( Figure 6C-3). Nano-phytoplankton dominated Zone II, wherein GPP and CR were highest and NCP was most negative ( Figures 6H, J, K, L). Zone III covered the offshore waters off Zone II and was dominated by pico-phytoplankton and characterized with lower nutrients and Chl a ( Figures 6I, D, E, F). Moreover, the FIGURE 6 | Regional variance of temperature (A,°C), salinity (B), DO (C, mg L -1 ), DIN (D, mmol L -1 ), DIP (E, mmol L -1 ), Chl a concentration (F, mg L -1 ), Micro-% (G), Nano-% (H), Pico-% (I), GPP (J, mg C m -3 d -1 ), CR (K, mg C m -3 d -1 ) and NCP (L, mg C m -3 d -1 ) in surface (-1, -2) and bottom (-3, -4) layer along transections of the GGBA in summer (-1, -3) and winter (-2, -4) periods. Different letters in each box plot indicate the significant differences among different groups (Turkey' HSD test, p < 0.05) GPP and CR in Zone III were low too, and the NCP was even closer to zero ( Figures 6J, K, L).
In winter, the bottom water of Zone I covered similar area as in summer, and had the highest nutrients concentration and lowest temperature and salinity among all three zones ( Figures 6A, B, D, E). Moreover, the Chl a and NCP in Zone I were higher than Zone II or III, and nano-phytoplankton dominated therein ( Figures 6F, H, L). Zone II is located in the Daya Bay and covered its offshore area, but did not cover the coastal waters of the PRE and Huangmaohai Bay. Zone III covered the offshore waters of the investigated area except that off the Daya Bay. Both the Zones II and III were characterized with the higher temperature and salinity than Zone I, as well as lower nutrients (Figures 6A, B, D, E). The pico-Chl a allocation was higher in these two zones ( Figure 6I-4). Finally, the Chl a biomass in Zone II was higher than that in Zone III, while the proportion of pico-Chl a was lower ( Figure 6; Table S1).

Regulation of Plankton Community Metabolisms
The environmental factors that influence the production and consumption of plankton community include temperature, light, nutrients, and turbidity (Rochelle-Newall et al., 2007;Kemp and Testa, 2011). Previous studies showed that temperature and nutrients are key factors to regulate the primary production and community respiration of the estuaries and their adjacent coastal areas (Caffrey, 2004;Caffrey et al., 2014;Hung et al., 2014). Consistently, in this study temperature and nutrients were identified as the key environmental factors for regulating the metabolisms of the plankton community of the GGBA, although the extent of these regulations on the GPP and CR differed greatly. According to the optimal regression model, field temperature explained more variance of the surface CR as compared to other variables, while the nutrients explained more of the GPP (Figure 7). Furthermore, the surface NCP was positively correlated to GPP (p < 0.01), but not to CR (p > 0.05) ( Table S2), again indicating that the nutrients regulate the GPP and thus vary the NCP of the surface waters of GGBA. However, such differences in regulatory factors between primary production and respiratory metabolism were not obvious in the bottom waters, since the multiple regression analysis result suggests the similar key environmental factors explaining the variation of GPP and CR, both of which are significantly correlated with NCP (p < 0.01) (Figure 7; Table S2).
Coastal ecosystems of the GGBA are usually influenced by the Pearl River runoffs, Yuedong Coastal Current, and intruded oligotrophic oceanic water from the South China Sea (Yin   , 2004a;Harrison et al., 2008;Song et al., 2011;Zu et al., 2014;Zu and Gan, 2015). Together with the influences of monsoon, tide, topography, and anthropogenic inputs, the metabolic processes of plankton community in this area varied greatly in spatial and seasonal scales (Figure 4). In this case, the regional analysis would be helpful to get a clear understanding about the complex relationship between community metabolism and driving forces in each regions. According to the PCA results, the GGBA was divided into different zones ( Figure 5), corresponding to the different extent of the comprehensive effects of driving factors like the Pearl River runoffs, anthropogenic activities, and offshore water currents. These regional differences are more obvious in the upper water layer due to the generally higher and more various impacts of the driving factors therein ( Figure 6; Table S1; Song et al., 2011;Dai et al., 2014;Zu et al., 2014). The upper reach of the PRE (Zone I) was mainly influenced by the runoffs of Pearl River, with the influence range extending more seaward in summer. This region had the highest nutrient level but not highest Chl a and GPP ( Figure 6; Table S1), probably due to the restriction of solar irradiation caused by the high turbidity and low transparency and disturbance therein (Lu and Gan, 2015). However, the metabolic state of this region was generally autotrophic in summer, inconsistent to other eutrophic estuaries where they were usually reported to be heterotrophic (Ducklow and McCallister, 2004;Caffrey et al., 2014), which may be explained by the greater negative effect of high disturbance upon the CR than the GPP (Alizadeh et al., 2018).
The river plume-dominated area (Zone II) often expands wider in summer due to the increased land-derived runoffs (Figure 2; Zu and Gan, 2015;Zhao et al., 2020) and was characterized with high Chl a and primary production ( Figure 6; Table S1), consistent with previous studies (Yin et al., 2000;Yin et al., 2004b;Harrison et al., 2008;Song et al., 2011;Lu and Gan, 2015). Nutrients in this area are plentiful, and the restrictions of temperature, turbidity, and so on are less, probably leading to the higher primary production and plankton respiration ( Figure 4A, B). In particular, the surface CR in this region is close to or even higher than the GPP, which may be attributed to the active metabolic decomposition during the postalgal bloom period (Figures 2, 4; Table S1; Buchan et al., 2014;Wemheuer et al., 2014), thus resulting in the low NCP and even a partially heterotrophic state.
The eastern nearshore area of the GGBA (Zone III) was less affected by the runoffs of Pearl River. So, the nutrient sources of this region are mainly taken from anthropogenic activities and terrigenous inputs, as well as the vertical upwelling in a certain period of rainy season (Song et al., 2015;Wu et al., 2017;Zhang et al., 2019), which may increase the GPP and CR to some extent. Moreover, the metabolic process in this zone was more active than the upper reaches of the PRE as well, possibly due to low turbidity. This coastal water showed an obvious autotrophic state in its surface layer in both summer and winter, being consistent with other coastal waters of e.g., the coastal waters of the Red Sea and the northern Australia (McKinnon et al., 2017;Loṕez-Sandoval et al., 2019). The far-offshore region (Zone IV), where had low GPP and CR, is understandable due to the offshore oligotrophic water with low nutrients and organic matters. In addition, a weak heterotrophy occurred in this region, indicating the primary production could not satisfy the requirements of community metabolism. It is consistent with the results from the similar oligotrophic oceanic waters Regaudie-de-gioux and Duarte, 2013).
The bottom water in the GGBA was less influenced by the Pearl River runoffs, and thus more stable than surface water, with the environmental factors varying less in both temporal and spatial scales (Figures 2, 6; Zu et al., 2014). Due to lacking solar light, the GPP therein was low; so, the regional variation of NCP could mainly be dependent on the CR (Table S2). Moreover, efficient transport of organic and inorganic nutrients via terrestrial inputs or downward depositions from the high-GPP upper water may make the community metabolism more active in nearshore areas (Ke et al., 2019;Dan et al., 2021), while in bottom of the offshore waters, the community metabolism was less active ( Figure 6; Table S1), probably resulting from the lacking nutrient inputs.
The Greater Bay Area covers the coastal areas with spatially environmental diversity. The corresponding spatial variations of key environmental factors, such as nutrient supply, turbidity and temperature lead to the regional complexity of plankton metabolic processes. According to the water-column integrated data, the plankton community metabolisms in winter shifted from heterotrophy in upper reaches of the estuary to autotrophy in nearshore waters and back to heterotrophy again in offshore waters, consistent with other coasts (Gattuso et al., 1998;Ducklow and McCallister, 2004). In summer however, the expansion of river plume led to a large-scale heterotrophy in coastal waters off the PRE accompanying with an algal bloom, which has seldom been reported before. This is possibly because that the heterotrophic processes during the algal bloom are more active than primary production, thus leading to the negative state of metabolic balance .

Plankton Metabolisms as Indicators of Ecosystem Health and Carbon Budget
The metabolic status of the plankton community is an important indicator for the health and carbon budget of an ecosystem. Liu et al. (2012) reported the close relationship between high primary production and algal blooms in nearshore waters of the Daya Bay and suggested the primary production as an index of coastal eutrophication. Similarly, here, an obviously high GPP was associated with an algal bloom [featured as maximum Chl a >20 mg L -1 , more than twice of the criterion value described by He et al. (2013)] in the Pearl River plume during summer period (Figures 2, 4). Previous studies showed that this area is conducive to the primary production process and biomass accumulation due to the large amount of nutrient inputs from river runoffs, as well as the more stable water column in summer, resulting in the frequent records of algal blooms (Yin et al., 2000;Yin et al., 2004b;Harrison et al., 2008;Song et al., 2011;Lu and Gan, 2015;Alizadeh et al., 2018). In contrast, the upper reaches of the estuary are less likely to form algal blooms, since phytoplankton production therein are usually limited by turbidity and river-flow disturbance.
Hypoxia often occurs in the bottom water of the GGBA in summer, such as in the PRE and Daya Bay (Yin et al., 2004a;Song et al., 2015;Li et al., 2018;Lu et al., 2018;Qian et al., 2018). In summer, obvious stratification often forms in estuarine and coastal waters, which may block the transportation of oxygen from upper to bottom layers and thus hydrologically increase the probability of hypoxia (Figure 2; Lu and Gan, 2015;Zhao et al., 2020). On the other hand, high production in the upper layer usually enhances the downward transportation of organic carbon, which may have favored the community respiration in the bottom layer, and thus led to the negative NCP value and heterotrophy therein (Li et al., 2019). Consistently, the higher bottom CR, as compared to its adjacent waters, occurred around the high surface GPP area within the river plume ( Figure 4). In addition, the benthic respiration (mainly microbial metabolism) may also benefit from the downward transportation of excess organic matters that are associated with high primary production in surface water (Aller et al., 2008;Zhu et al., 2011). Therefore, the high oxygen consumption in summer due to the active community metabolism and long-term stratification may result in the low DO in bottom water ( Figure 2). It is also suggested that the formation of hypoxic zone is closely related to the massive growth and reproduction of surface phytoplankton cells (Qian et al., 2018). More interestingly, the bottom hypoxia zone in the estuary was not rightly located under the areas with high surface Chl a and primary production but under its adjacent area (Figure 2). This can be explained by the saltwater wedge effect that has transported the sinking organic matters to the nearby area through the current intrusion in lower or bottom layer (Hu and Li, 2009). However, this is not the case in the Daya Bay because there is no such wedge effect (Figure 2). The hypoxia seldom happened in winter either (Figure 2), possibly due to the strong vertical mixing caused by northeastern monsoon that has transported the high DO surface water to the bottom layer and thus enhanced the DO level therein. The lower temperature in winter may be another cause through enhancing oxygen solubility and lowering respiration metabolic rate.
A previous study suggested that annual sea-air CO 2 exchange near the PRE was close to equilibrium, and the CO 2 -absorbing capacity was enhanced in summer (Zhai et al., 2013). In this study, the primary production and community respiration in the GGBA is close to be balanced in both winter and summer, and NCP is also higher in summer (Mann-Whitney test, p < 0.05), indicating the high similarity and potential correlation between the community metabolic carbon flux and carbon budgets. Moreover, the specific distribution of plankton community carbon budget varied greatly in spatial and temporal scales, making it difficult to exactly evaluate the detailed carbon flux. It is generally believed that the upstream of estuaries is the source of CO 2 , while the waters controlled by a river plume is the sink due to the high primary production (Guo et al., 2009;Dai et al., 2014). In this study however, the metabolic respiration appeared to be more active than primary production in river plume in summer, probably contributing to the carbon export of this area. In the eastern coastal part of the GGBA, the plankton community metabolisms were much less affected by the Pearl River runoffs, making this area appear to be a carbon source in both winter and summer (Figure 4). This heterotrophy also existed in the nearshore waters of the Daya Bay during a certain summer period , suggesting the complication of ecosystem carbon flux even in a small-scaled area like the GGBA due to its temporal and spatial variations of the impacts of anthropogenic activities and other environmental factors.
In conclusion, we found that the field plankton community metabolisms in the GGBA showed a great regional difference in both summer and winter periods, interactively influenced by the physical-environmental-biological factors. The surface water of the GGBA exhibited autotrophic in summer, while such an autotrophic in summer status was weaker in winter; however, the bottom water exhibited heterotrophic in both summer and winter periods. Our results also showed that the high surface GPP and bottom CR can be considered as the important indicators for the algal blooms and hypoxia events that happened during that time in the GGBA. Further studies on the dynamic variations of plankton metabolism and its driving mechanisms are necessary for an in-depth understanding of the function of coastal ecosystem and better assessment and management of potential ecological risks.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
LZ, XS, and GL designed the experiments and analyzed the data. CX, YH, XF, and CZ conducted the shipboard sampling and environmental measurements. LZ and YH performed the in situ plankton metabolism measurements. ZW and ZO conducted nutrient measurements. LZ, XS, GL, and CX wrote, discussed, and edited the paper. All authors contributed to the article and approved the submitted version.