Eff ects of dissolved oxygen and nutrients from the Kuroshio on hypoxia off the Changjiang River estuary*

Di TIAN , , Feng ZHOU ,,3,**, Wenyan ZHANG , Han ZHANG ,5, Xiao MA , Xinyu GUO ,

1 State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou 310012, China

2 Observation and Research Station of Yangtze River Delta Marine Ecosystems, Ministry of Natural Resources, Zhoushan 316022, China

3 Ocean College, Zhejiang University, Zhoushan 316021, China

4 Institute of Coastal Research, Helmholtz-Zentrum Geesthacht, Geesthacht 21502, Germany

5 Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China

6 Center for Marine Environmental Study, Ehime University, Matsuyama 790-8577, Japan

Abstract The intrusion of the Kuroshio into the East China Sea (ECS) aff ects the development of hypoxia off the Changjiang (Yangtze) River estuary; however, quantitative analysis of its impacts is lacking. In this study, the Regional Ocean Modeling Systems (ROMS) model coupled with the Carbon, Silicate and Nitrogen Ecosystem (CoSiNE) model was used to investigate the relative importance of dissolved oxygen(DO) and diff erent nutrients (silicate, nitrate, and phosphate) in the Kuroshio on hypoxia in the ECS. Results show that changes in DO concentrations in the Kuroshio modify the distribution and intensity of hypoxia through direct onshore transport by hydrodynamic processes. An increase in Kuroshio DO concentration by 25% or 50% would result in a decrease of the maximum hypoxia extent (MHE) in the ECS by 76% or 86%,respectively, while a 25% decrease in Kuroshio DO would increase the MHE by up to 219%. The contribution of DO in the Taiwan Strait is almost negligible. In contrast to Kuroshio DO, nutrients aff ect hypoxia in the ECS through onshore transport by hydrodynamic and biochemical processes. Changes in phosphate and nitrate concentrations by 25% in the Kuroshio would change the MHE by up to 30% and 18%, respectively,accompanied by apparent changes in surface chlorophyll- a concentrations. The eff ect of silicate on hypoxia is negligible because a 25% change in silicate concentrations in the Kuroshio would result in less than 1%change in the MHE. Our results reveal a hierarchical rank of importance for environmental variables in the Kuroshio (i.e., DO ＞ phosphate ＞ nitrate ＞ silicate) in modifying the development of hypoxia in the ECS.

Keyword: hypoxia; Kuroshio; Changjiang River estuary; East China Sea; dissolved oxygen; nutrient

1 INTRODUCTION

Hypoxia, conventionally def ined as dissolved oxygen (DO) concentrations ≤2.0 mg/L (62.5 mmol/m3;1.4 mL/L), has expanded rapidly throughout coastal areas of the world since the mid-20thcentury (Gilbert et al., 2010; Rabalais et al., 2010). More than 500 coastal hypoxic zones have been reported and are typically considered to result from human activities through over-enrichment of nutrients (Diaz and Rosenberg, 2008; Breitburg et al., 2018a). Def icits in DO threaten the marine ecosystem, which lead to reductions in ocean productivity and biodiversity as well as changes to biogeochemical cycles (Vaquer-Sunyer and Duarte, 2008; Levin et al., 2009; Zhang et al., 2010). The biological eff ects of hypoxia are strengthened in conjunction with ocean acidif ication(Melzner et al., 2013). Currently, hypoxia is one of the most challenging global marine environmental issues (Diaz, 2001; Keeling et al., 2010; Breitburg et al., 2018b).

In the East China Sea (ECS), which has one of the largest continental shelves in the world, a seasonal(June-October) bottom water DO def icit with complex spatial and temporal variations, has been frequently observed off shore of the Changjiang (Yangtze) River estuary (CJE) in recent decades (Tian et al., 1993; Chen et al., 2007; Wei et al., 2007; Zhu et al., 2011, 2017; Ni et al., 2016). Hypoxia in the ECS was f irst detected in the late 1950s (Gu, 1980) but did not receive much attention until summer hypoxia of 1999 was welldocumented (Li et al., 2002). In that summer, the hypoxic zone (＜2 mg/L) was estimated to encompass～13 700 km2at an average thickness of 20 m and a minimum DO concentration of 1 mg/L (Li et al., 2002).Over the past decades, the ECS hypoxic area expanded from 10 000 to 20 000 km2(Wei et al., 2007; Wang et al., 2012), comparable to one of the largest coastal hypoxia areas in the world (northern Gulf of Mexico;Rabalais et al., 2002; Bianchi et al., 2010).

The formation and evolution of coastal hypoxia is controlled by multiple factors, which are closely related to both physical and biogeochemical processes (Zhang and Wirtz, 2017; Zhang et al., 2019b). The mechanisms controlling seasonal hypoxia in the off shore area of the CJE have been widely investigated via f ield observations (Wang, 2009; Li et al., 2011; Wang et al.,2012; Luo et al., 2018) and numerical modeling (Zhou et al., 2017; Zhang et al., 2019a, 2020; Zhou et al.,2020). It is generally suggested that this hypoxia formation and evolution is highly modulated by wind(Ni et al., 2016; Guo et al., 2019; Liblik et al., 2020),runoff (notably from the Changjiang River; Zhang et al., 2018), Kuroshio intrusion (Qian et al., 2017; Zhang et al., 2018; Zuo et al., 2019), nutrient loading (Zhu et al., 2017), and tides (Zhu et al., 2017). Bottom topography is also regarded as a factor inf luencing the residence time of hypoxic water (Wang, 2009).

The Kuroshio, a western boundary current of the Pacif ic Ocean, originates east of the Philippines,intrudes into the South China Sea (SCS) through the Luzon Strait and into the ECS mainly through the shelf-break northeast of Taiwan, China (Wu et al.,2017). The annual mean volume transport of the Kuroshio east of Taiwan, China (to 123°E) was 27.5 Sv (1 Sv=106m3/s), and the annual mean volume transport across the 200 m isobath was around 2.6 Sv in 2006 (Zhou et al., 2015). Kuroshio intrusion has a profound inf luence on marine circulation and ecosystems in the ECS (Zhang et al., 2007; Sugimoto et al., 2009; Yang et al., 2011; Tian et al., 2019; Meng at al., 2020). In the ECS, more nutrients are sourced from the ocean than land (Chen and Wang, 1999).Wong et al. (2000) showed that the nutrient supply of the ECS shelf is mostly from subsurface Kuroshio water and, to a lesser extent, from the Changjiang River. The degradation of marine organic matter, as well as plankton respiration, represents the dominant consumption of oxygen in the bottom waters of the CJE hypoxic zone (Wang et al., 2016).

Recently, some studies have reported the eff ects of Kuroshio intrusion on the development of hypoxia in the ECS (Qian et al., 2017; Zhou et al., 2017; Zuo et al., 2019; Große et al., 2020). The north-south movement of the ECS hypoxic zone is related to the northward extension of the Taiwan Warm Current(Wang et al., 2012). Qian et al. (2017) demonstrated that intruded bottom water masses, mixed by subsurface Kuroshio and SCS waters, exert a synergistic control on DO concentrations in the area off shore of the CJE and signif icantly inf luence hypoxia formation. Yang et al. (2018) suggested that relatively cold, oxygen-poor, and nutrient-rich deep water from the Kuroshio may inf luence ECS hypoxia.Meanwhile, subsurface Kuroshio intrusion also can inf luence the extent of bottom water hypoxia by regulating bottom water residence time off shore of the CJE (Zhang et al., 2019a). However, an opposing view was proposed by Zuo et al. (2019), who argued that the Kuroshio carries high concentrations of DO onto the continental shelf and alleviates the ECS oxygen def icit. This same phenomenon (Kuroshio intrusion) has been found to diminish the degree of coastal hypoxia off shore the Zhujiang (Pearl) River estuary in the northern SCS (Lui et al., 2020). Hence,the eff ects of Kuroshio intrusion on DO concentrations in coastal regions are still unclear. Furthermore, the relative contribution of diff erent nutrients and DO in Kuroshio water to ECS hypoxia is less understood.

Fig.1 Model domain and bathymetry

We addressed this lack of knowledge by investigating the inf luence of Kuroshio DO and nutrients (silicate, nitrate, and phosphate) on ECS hypoxia using a coupled physical-biogeochemical model. The eff ect of Taiwan Strait water on hypoxia is also discussed. In general, the conventional def inition of hypoxic water (DO concentrations of 2.0 mg/L) is below the empirical hypoxia threshold for more sensitive taxa. This has implications for the conservation of coastal biodiversity because hypoxia thresholds diff er greatly among marine benthic organisms (Vaquer-Sunyer and Duarte, 2008). To ensure consistency with previous research (e.g., Chen et al., 2007; Zhou et al., 2017, 2020), hypoxia is def ined in this study as DO concentrations less than 3.0 mg/L (2.1 mL/L).

2 METHOD

2.1 Coupled physical-biogeochemical model

We used a coupled physical-biogeochemical coupled model in this study. The physical model is based on Regional Ocean Modelling Systems(ROMS), which is a free surface, terrain-following and three-dimensional primitive equations ocean model (Shchepetkin and McWilliams, 2005;Haidvogel et al., 2008). The model domain covers the Bohai and Yellow seas and the ECS (23.5°N-41.0°N,117.5°E-132.0°E; Fig.1) at a horizontal resolution of 1/24° (～4 km) and 30 vertical layers in a s-coordinate system.

The model spin-up and forcing conditions (surface heat f lux, momentum, and fresh water; lateral boundary conditions; river runoff ) for a realistic run of 2006 were the same as Zhou et al. (2020). The hydrodynamic model was spun up for three years,initialized from a January climatological meanderived from multi-year observations (Chen, 1992)and driven by monthly forcing data from the Comprehensive Oceanographic Atmospheric Data Sets (da Silva et al., 1994). The realistic run of 2006 was initialized from the output of the last iteration of the spin-up. In this run, the surface f luxes of heat and precipitation/evaporation were from ERA-Interim data sets (European Centre for Medium-Range Weather Forecasts; Dee et al., 2011). The wind vector data were from the Blended Sea Winds (https://www.ncdc.noaa.gov/data-access/marineocean-data/blendedglobal/blended-sea-winds; Zhang et al., 2006; Peng et al., 2013). Lateral boundary conditions for temperature, salinity, current velocity, and sea surface elevation were from the Global Hybrid Coordinate Ocean Model realistic simulation for 2006 (Kelly et al., 2007). The monthly mean for Changjiang River runoff was also included in our realistic run. Further details on the setup of the ROMS model and its validation can be found in Zhou et al. (2015), which demonstrated a good model performance in capturing the key physics of the ECS.

Table 1 Sensitivity experiments, including the reference run (Ref-run), and percentages of dissolved oxygen(DO) and nutrients in each

The biogeochemical model is the Carbon, Silicate and Nitrogen Ecosystem (CoSiNE; Chai et al., 2002;Xiu and Chai, 2014) model, which consists of 11 stated variables, including four nutrients (silicate,nitrate, phosphate, and ammonium), two phytoplankton groups (small phytoplankton and diatoms), two grazer groups (microzooplankton and mesozooplankton), two classes of detritus (small/suspended and large/sinking), and DO, coupled online to physical processes. Among many biological parameters, the photosynthetically active radiation saturation for diatoms is 100.0 W/m2, the maximum specif ic growth rate for small phytoplankton and diatoms is 3.0/d and 3.5/d, respectively. The initial and boundary conditions of the CoSiNE model are derived from the output of the CoSiNE-Pacif ic model(Xiu and Chai, 2011). More details on the setup and parameters of the CoSiNE model and its validation for the ECS are in Zhou et al. (2017, 2020). Their results indicated that the development of hypoxia and observed vertical prof iles of DO are reasonably reproduced by this coupled model.

The realistic run for 2006 had four iterations,driven by the same realistic forcing but with initial conditions from the last iteration. This procedure can reduce the bias between model results and observations. Three-day-mean outputs of the last model year were saved for analysis. The realistic run of 2006 was considered the reference run (Ref-run).

The initial DO concentration is important for the numerical experiments. Moreover, the prognostic DO concentration of the model is also constrained by the boundary DO values where the monthly mean output from the large-domain model was prescribed. In the Ref-run, the spatial distribution of the initial DO concentration in the bottom layer and at surface, and the initially vertical distribution of DO along 23.5°N(the south boundary of the model domain) and at specif ic locations (representing the Kuroshio) are shown in Fig.2. The DO concentration in the Kuroshio varies notably with depth, ranging from ～7-～3 mg/L.The vertical structures (Fig.2d) are consistent with Zuo et al. (2019).

2.2 Sensitivity experiment

Fig.2 Initial dissolved oxygen (DO) concentration in the reference run

Nutrients and DO from the open ocean are transported into the ECS continental shelf via either the Taiwan Strait or Kuroshio (Fig.1). We performed several sensitivity experiments using diff erent DO concentrations in the open ocean to demonstrate their inf luence on ECS hypoxia (Table 1). In the f irst experiment (K-O125), DO concentrations in the Kuroshio and Taiwan Strait, represented by the oceanic DO concentration at the southern boundary,were increased by 25% relative to the Ref-run. In the second experiment (K-O75), oceanic DO concentrations at the southern boundary were decreased by 25%. The third experiment (K-O150)represented an extreme case in which oceanic DO concentrations at the southern boundary were increased by 50%. To distinguish whether the water is from the Kuroshio or from Taiwan Strait, additional two experiments were performed in the scenarios that DO concentrations at the southern boundary east of Taiwan, China, increased by 25% (KE-O125) or decreased by 25% (KE-O75), respectively. The diff erences between K-O125 and K-O75, and KEO125 and KE-O75, show the contribution of Taiwan Strait water. Our decision to use 25% changes(increases or decreases) in oceanic DO relative to the Ref-run is consistent with Große et al. (2020) and in accordance with changes in DO projected by Coupled Model Intercomparison Project Phase 5 (CMIP5)Earth system models over the western North Pacif ic Ocean by the end of the 21stcentury under the highemission Representative Concentration Pathway 8.5 scenario (Bopp et al., 2017).

To explore the response of ECS hypoxia to nutrient loading (nitrate, phosphate, and silicate were treated separately here) from oceanic sources (Kuroshio and Taiwan Strait), two experiments were designed for each nutrient (Table 1) with corresponding concentrations at the southern boundary increased or decreased by 25%. The sensitivity experiment setting is based on an analysis of observation data that showed nitrate concentrations in the Kuroshio north of Taiwan, China, increased by ～25% from the 1980s to 2000s (Guo et al., 2012).

Although the Kuroshio intrusion into the ECS mainly occurs at the shelf break northeast of Taiwan,China, it also can enter the Taiwan Strait after intruding the Luzon Strait and become a part of Taiwan Strait water. Xu et al. (2020) suggested that Taiwan Strait water has little inf luence on the production of chlorophyll off shore of the CJE. The inf luence of oceanic DO from the Taiwan Strait water on hypoxia off shore of the CJE is negligible relative to the Kuroshio (as discussed in the Section 3.1).Therefore, the inf luence of DO and nutrients from the open ocean are mainly attributed to the Kuroshio.

2.3 Maximum hypoxia extent

Following Zhou et al. (2017), the maximum hypoxia extent (MHE) was adopted to represent the most probable spatial coverage of hypoxia during a specif ied period (e.g., 1 month or 1 year):

Omin(i,j)=min[O(i,j, 1), ∙∙∙,O(i,j,k)],

Fig.3 Anomalies in dissolved oxygen (DO) concentration in the bottom water averaged from July to September in DO sensitivity experiments relative to the reference run (Ref-run)

Table 2 The variation ratio (%) of the maximum hypoxia extent (MHE) in a specif ied area (28°N-33°N,121°E-125°E) for each sensitivity experiment with respect to the reference run

whereOis the DO concentration (mg/L) in the bottom layer,Ominis the minimum DO concentration (mg/L)within a specif ied time period,kis the time index;mask(i,j) is an index (1=hypoxia; 0=otherwise; a non-dimensional number) at each model grid; ifOmin(i,j) ≤3.0 mg/L, then mask(i,j)=1, otherwise,mask(i,j)=0;Ais the grid area (km2) andiandjare grid indices in the meridional and zonal directions,respectively. Consistent with the MHE (km2), the 3-day-mean hypoxia extent (km2), which shows the synoptic variability of hypoxia, is also represented by the areas with DO concentrations less than 3.0 mg/L.

3 RESULT

3.1 Inf luence of Kuroshio to DO concentration

The DO concentration in the ECS is widely inf luenced by intruded DO from the Kuroshio (Fig.3).With a 50% increase in open ocean DO content, the DO concentration in the core hypoxic area(29°N-32°N, 122°E-124°E) increased by ～23%during the hypoxic season (July-September); with a 25% increase, the resulting DO concentration increased ～11%. Correspondingly, the area of hypoxic water represented by the MHE decreased along with higher DO concentrations in the open ocean (K-O150 and K-O125) compared with the Ref-run (Fig.4a-b).The distinctive contraction of hypoxia was found in the northern, eastern, and southern areas of the CJE.Compared with the Ref-run, the MHE in K-O125 decreased 76% annually and 75% in August when hypoxia was most signif icant (Table 2). The diff erence in MHE between K-O150 and Ref-run was even larger at ～86% in August and annually (Table 2). In K-O75, the experiment with a 25% reduction in open ocean DO concentrations, the MHE increased 219%annually and 128% in August compared with the Refrun. Hypoxic water also spread signif icantly southward and eastward (Fig.4c).

Fig.4 Spatial distribution of annual minimum dissolved oxygen (DO) concentrations in the bottom water from sensitivity experiments

Fig.5 Comparisons of the reference run (Ref-run) and sensitivity experiments

In addition to the spatial extent of hypoxia, its duration and the minimum DO concentration are also sensitive to changes in oceanic DO concentrations, mainly supplied by the Kuroshio.Variations in the time series of the 3-day-mean hypoxic area were similar among the sensitivity experiments, but there were signif icant diff erences in their magnitude and duration, although the onset of hypoxia was similar (Fig.5a). Generally, lower DO concentrations in the Kuroshio prolong hypoxia,while higher DO concentrations shorten its duration.For example, in K-O75, the hypoxia persisted until October, which was 1 month longer than the Ref-run(Fig.5a). Moreover, the diff erence between K-O125 and the Ref-run was much larger than between K-O150 and K-O125 (but smaller than between K-O75 and the Ref-run), which suggests that hypoxia exhibits a nonlinear response to changes in DO concentrations.

The hypoxic area represented by the MHE was largest in August in all runs (Fig.5c), with the lowest minimum DO concentration appearing at the same time (Fig.6) in the area off shore of the CJE (Fig.7). In August, the MHE was most vulnerable to changes in Kuroshio DO concentrations. We also determined that the hypoxic area is more aff ected by changing Kuroshio DO concentrations than hypoxia intensity,considering that the bottom water minimum DO concentrations in August were similar in the diff erent sensitivity experiments.

Fig.6 Time series of minimum dissolved oxygen (DO) concentrations in 2006

The DO concentrations off shore the CJE are dominated by oceanic DO input from the Kuroshio.In KE-O125 and KE-O75, the DO concentrations in the core hypoxic area (29°N-32°N, 122°E-124°E)varied by ～9% relative to the reference run, slightly smaller than the changes (～11%) in K-O125 and K-O75 (Fig.3). The spatial distribution of hypoxic water in KE-O125/KE-O75 were similar to K-O125/K-O75 (Fig.8), which suggests that the hypoxic area is mainly aff ected by Kuroshio intrusion into the bottom water of the shelf-break northeast of Taiwan,China, and any eff ect from the Taiwan Strait is almost negligible.

3.2 Inf luence of nutrients from the Kuroshio on hypoxia

Variations in major nutrient (nitrate, phosphate,and silicate) concentrations of the Kuroshio were not as important as DO variations to the development of hypoxia off shore of the CJE and in ECS bottom waters (Figs.9 & 10). A 25% change in each nutrient only led to a 2% variation in bottom water DO concentrations over the core hypoxic area in July-September (Fig.3). Generally, the MHE expands slightly with increased nutrient concentrations and shrinks westward when nutrient concentrations decrease. The largest change in the monthly MHE in all sensitivity experiments occurred again in August(Table 3). The hypoxia extent was more sensitive to changes in phosphate concentrations than other nutrients because phosphate is usually a limiting factor for phytoplankton growth off shore of the CJE.With a 25% increase in phosphate concentrations, the MHE grew by 19% in August; there was also a 13%growth in the MHE when nitrate was increased by 25%. With a 25% decrease in phosphate and nitrate concentrations, the MHE contracted by 28% and 18%, respectively, in August. These quantitative results suggest that the MHE is more sensitive todecreased nutrient concentrations than to increases in nutrients. The impact of silicate concentrations on the MHE was negligible (Table 3).

Table 3 Change ratio (%) in the maximum hypoxia extent (MHE) over a specif ic area (28°N-33°N, 121°E-125°E) for each experiment with respect to the reference run (MHE=100%)

Fig.10 Spatial distribution of annual minimum dissolved oxygen (DO) concentrations in the bottom water from sensitivity experiments

4 DISCUSSION

4.1 Eff ects on surface chlorophyll a

Our sensitivity experiments indicate that intrusions of the Kuroshio aff ect hypoxia in the area off shore of the CJE in diff erent ways. The DO content in the Kuroshio aff ects the ECS bottom water DO concentrations directly through advection and diff usion, which are hydrodynamic processes associated with Kuroshio intrusion; these processes eventually modify the intensity and distribution of hypoxia. In contrast, the eff ect of this DO intrusion on biochemical processes is limited (Fig.11). A 50%change in the DO in the Kuroshio led to only a slight change (＜2%) in ECS phytoplankton concentrations.Hence, these results support the conclusion that intruded Kuroshio water masses exert a synergistic control on the initial DO concentration that led to the development of hypoxia in the ECS, and future DO decreases in the open ocean may exacerbate hypoxia(Qian et al., 2017).

Nutrient concentrations in the Kuroshio aff ect the ECS hypoxic zone through hydrodynamic and biogeochemical processes. Among nitrate, phosphate,and silicate, phosphate is the chief limiting factor on phytoplankton growth. Our results show that increased concentrations of nitrate and phosphate stimulated phytoplankton growth (represented by chlorophyll-aconcentrations) in Zhejiang coastal water, whereas nutrient decreases constrained growth (Fig.12a-d). In the area off shore of the CJE and its adjacent sea, surface chlorophyll-aconcentrations in K-P75 and K-P125 obviously varied relative to the Ref-run, with changes mainly 6% to 12% more than in K-N75 and K-N125,especially in the region of water depths ＜30 m. Photic zone export and deeper degradation of high amounts of biomass eventually overtake oxygen consumption in the bottom water. Overall, the inf luence of silicate is slight (Fig.12e & f), which is consistent with previous results that silicate limitation on phytoplankton growth in the CJE rarely occurs (Zhou et al., 2008).

4.2 Relative importance of riverine and oceanic nutrients

Fig.11 Changes in Surface chlorophyll- a (Chl- a) concentrations averaged from July to September in DO sensitivity experiments relative to the reference run (Ref-run)

Fig.12 Changes in surface chlorophyll- a (Chl- a) concentrations averaged from July to September in nutrient sensitivity experiments relative to the reference run (Ref-run)

Hypoxic water is associated with both riverine and oceanic nutrients. Zhou et al. (2017) suggested that annual mean MHE would increase by 29% if all nutrients carried by the Kuroshio were elevated by 25%, which is consistent with our results in the nutrient sensitivity experiments (Table 3). Our study further quantif ied the relative importance of three major nutrients in the Kuroshio as phosphate ＞ nitrate＞ silicate. The relative importance of riverine nutrients in aff ecting hypoxia is expected to be smaller compared with the Kuroshio. Specif ically, the MHE should expand by 14% when all nutrient loadings from the Changjiang River are increased by 30%(Zhou et al., 2017). That is partly because the nutrient f lux into the ECS from the open ocean exceeds the Changjiang River (Chen and Wang, 1999; Wong et al., 2000; Zhang et al., 2007). Chen and Wang (1999)suggested that rivers contribute only 7% of the total phosphorus input and about 38% of the total nitrogen input in the ECS. Xu et al. (2020) also pointed out that summer algal blooms in the northern ECS are mainly a result of riverine nutrient input from the Changjiang River (contribution of 59%), while in the southern ECS blooms are supported by oceanic nutrients(contribution of 85%). Moreover, riverine nutrients are quite rich in nitrate and silicate but low in phosphate. Compared with oceanic sources, riverine nitrogen from the Changjiang River has a larger eff ect on oxygen consumption in the ECS (Große et al.,2020). In our study, the nitrate concentration at the mouth of the Changjiang River was about 95 mmol/m3in the simulated year, lower than silicate(～102 mmol/m3) but almost two orders of magnitude higher than phosphate (～1.58 mmol/m3). Nitrate transport from the Changjiang River was ～2.1 kmol/s(66×109mol/a) and phosphate transport was～0.034 kmol/s (1.05×109mol/a). The riverine N/P ratio was ～62, which is smaller than a previous estimate of 111 (Chen and Wang, 1999), but we only considered NO3in our study; however, our N/P ratio was still higher than the Redf ield ratio (16) for phytoplankton growth. According to Chen and Wang(1999), seawater transported to the ECS has an N/P ratio of 21 (N: 301×109mol/a; P: 14.5×109mol/a),which is much closer to the Redf ield ratio. That is also the reason surface chlorophyllain the ECS is more sensitive to variations in Kuroshio phosphate concentrations with respect to nitrate and silicate.

4.3 Future improvement

In the sensitivity experiments of our study, the amplitude of relative change in DO and nutrients is set to be same for all types of Kuroshio waters.Although the intruded Kuroshio Subsurface Water is proved to be more important for the ECS than other types of Kuroshio water (Kuroshio Surface Water,Kuroshio Intermediate Water, and Kuroshio Deep Water). However, further sensitivity experiments with diff erent amplitudes of relative change in DO and nutrients for diff erent types of Kuroshio water would be worth to quantitatively investigate in future work if suffi cient f ield data are collected to support the scenarios.

5 CONCLUSION

The eff ects of DO and diff erent nutrients (silicate,nitrate, and phosphate) in the Kuroshio on the development of hypoxia in the ECS were studied with a validated ROMS-CoSiNE model. Our results suggest that DO concentrations in the Kuroshio modify the distribution and intensity of hypoxia in the ECS mainly through hydrodynamic processes associated with intrusion of Kuroshio water. Kuroshio DO concentrations have a limited eff ect on phytoplankton growth in the ECS. The extent of hypoxia, represented by the MHE, showed a negative correlation with DO concentrations (R=0.96,P＜0.01).Relative to the Ref-run, the MHE in the ECS decreased by 76% when Kuroshio DO concentration increased by 25%, whereas the MHE increased by 219% when Kuroshio DO concentration decreased by 25%. A further increase in Kuroshio DO concentrations by 50% would lead to a decrease of the MHE by 86%compared with the Ref-run. The eff ect of DO from Taiwan Strait water was almost negligible compared with the oceanic source east of the Taiwan Strait.

The eff ect of Kuroshio nutrients on hypoxia in the area off shore of the CJE was opposite to that of Kuroshio DO. Increased nutrients typically stimulate phytoplankton growth, resulting in an expansion of the MHE, whereas decreased nutrients suppress phytoplankton and contract the MHE. The MHE off shore of the CJE is most sensitive to changes in Kuroshio phosphate concentrations (～22%-30%change in the MHE with a 25% change in phosphate),while it is least sensitive to changes in silicate(～0.3%-0.9% change in the MHE with a 25% change in silicate). The MHE changed by ～14%-18% with a 25% change in nitrate. Our modeling results suggest that the DO concentration in the area off shore of the CJE is inf luenced by both intruded DO and nutrients from the Kuroshio, and variations in DO in the Kuroshio seem to play a more important role than nutrients in the development of hypoxia. In reality,decreases in DO and increased nutrients in the Kuroshio always occur simultaneously, so their combined eff ects on ECS hypoxia are nonlinearly dependent on their sum.

6 DATA AVAILABILITY STATEMENT

The datasets generated for this study are available on request to the corresponding author. Model outputs used for the f igures and tables will be available at https://www.soed.org.cn/index.php/scientif ic/data.

7 ACKNOWLEDGMENT

We thank anonymous reviewers for their constructive comments that helped to improve this paper.