Seasonal Variations of Marine Environment and Primary Production in the Taiwan Strait

The first data set of seasonal marine environment and euphotic zone integrated primary production (IP) variations in the Taiwan Strait was reported. The measured annual IP was 123 ± 86 gC m–2 year–1 (338 ± 235 mgC m–2 day–1), and its seasonal variations can be described with a left-skewed normal distribution curve. The average seasonal IP values from the highest to the lowest were summer (664 ± 270 mgC m–2 day–1), autumn (350 ± 118 mgC m–2 day–1), spring (202 ± 110 mgC m–2 day–1), and winter (137 ± 68 mgC m–2 day–1). The lowest IP was during the nutrient-rich winter because it had a short insolation duration, low incident photosynthetic active radiation (PAR) and low light transmission (shallow euphotic zone depth) due to strong vertical mixing. In contrast to the winter, the highest IP was during the nutrient-depleted summer, which had a long insolation duration, high incident PAR and high light transmission (deeper euphotic zone depth). In addition, the heterotrophic nutrients from upwelling in the south might also support the highest IP in summer. As three primary water masses exist in the Taiwan Strait and three of them have different characteristics, different mixing ratios of water masses may cause different chemical and hydrographic conditions, which leads to different levels of Chl a concentrations and primary production. It is worth to mention that offshore wind farm (OWF) construction in the Changyun Rise (CYR) of the Taiwan Strait is on-going. As primary production is the foundation for a marine ecosystem and supports the food web and fish stock, the results of this research can not only be used as the baseline for evaluating the OWF impact on the marine ecosystem but also be used for assessing their influence on fishery resources.


INTRODUCTION
The Taiwan Strait, with an average depth of 60 m, is located in the western Pacific Ocean and is an important channel in transporting water and chemical constituents between the East China Sea and the South China Sea (Chen et al., 2010;Hong et al., 2011;Huang et al., 2015). It is also an important path way for ship navigation and fish migration (Hong et al., 2011). The Taiwan Strait is strongly influenced by the East Asian monsoon. Wind stress explains a majority of the transport reversals (Ko, 2003). As the wind directions have changed, the transportations and circulations of water masses in the Taiwan Strait have changed (Jan et al., 2002(Jan et al., , 2006Chen, 2003Chen, , 2011Chen and Sheu, 2006;Wu et al., 2007;Chen et al., 2010). The wind fields over the strait are southwesterly dominated in summer and northeasterly dominated in winter. Seasonal changes of the wind directions have resulted in different currents and circulations, which brings different water masses into the Taiwan Strait (Jan et al., 2006). Jan et al. (2010) revealed three primary water masses in the Taiwan Strait: the Kuroshio Branch Water (KBW) with high temperature and high salinity, the Mixed China Coastal Water (MCCW) with low temperature and low salinity, and the South China Sea Water (SCSW) with intermediate temperature and salinity. Additionally, tropical storms and typhoons also influence this area during the summer via induced heavy rainfall, strong vertical water mixing or upwelling (Ko et al., 2016). Jan et al. (2002) pointed out that the deep Penghu Channel guides the intruding KBW northward into the Changyun Rise (CYR) and partially blocks along-strait transport on the eastern side of the strait in winter. On the other hand, the MCCW is conveyed by the northeast monsoon-driven flow to the northern Taiwan Strait during the winter (Jan et al., 2002). In addition, Jan et al. (2010) found that the source water in the Taiwan Strait comes from the SCSW in summer. Huang et al. (2018) also reported that during the summer and autumn, the main water mass in the Penghu Channel is the SCSW. Moreover, since precipitation is heavier in summer, water in the top layer is warmer and fresher during the summer than during the autumn . As three primary water masses in the Taiwan Strait have different characteristics, different mixing ratios of water masses may cause different chemical and hydrographic conditions. This dynamic will lead to differing levels in chlorophyll a (Chl a) and primary production, which may be influenced by the phytoplankton community structure (Hong et al., 2011).
Since the wind energy potential of the offshore areas of Taiwan's west coast is excellent (Archer and Jacobson, 2005;Fang, 2014;Chen et al., 2015) and Taiwanese government has sought to increase the use of renewable energy, 14 offshore wind farms (OWFs) are planned for installation in the CYR of the Taiwan Strait, offshore of Changhua county (Figure 1). They will have a wind power capacity of approximately 2400 MW, which will occupy 62.6% of the total offshore wind power capacity in Taiwan. Although wind energy is commonly recognized as a clean and environmentally friendly renewable energy resource, the potential environmental impacts due to installation, operation, maintenance and removal phases of the wind turbines cannot be ignored (Wilhelmsson and Malm, 2008;Smyth et al., 2015). Therefore, marine environmental investigations before OWF installation are needed to ensure and maintain the sustainability of both the social economy and marine environment. To understand the seasonal variations of marine environment and primary production in the Taiwan Strait and possible future impact of OWF on marine ecosystem, four seasonal cruises in 2017 were conducted. Hydrography, nutrients, Chl a and primary production conditions were investigated during the cruises.

Cruise Information and Sampling Instruments
Seasonal water samples were taken at five sampling stations in the Taiwan Strait (Figure 1) during four cruises on board an R/V Ocean Researcher 2: OR2-2212 (January 2017), OR2-2229 (April 2017), OR2-2253 (July 2017), and OR2-2269 (October 2017). Most of the sampling stations were in the CYR and all of them were located on the east of 119.5 • E in order to follow the "median line" agreement. The vertical profiles of seawater temperature, salinity, downwelling scalar irradiance, and light transmission (TM) were measured with a conductivity-temperature-depth instrument (CTD; SBE 9/11, Seabird Scientific, United States). Downwelling scalar irradiance was measured with a photosynthetic active radiation (PAR) sensor (Chelsea Technologies Group Ltd., United Kingdom). TM was measured with a 25 cm path length transmissometer (C-Star, Seabird Scientific, United States). Seawater samples from various depths were collected using a carousel water sampler (SBE32, Seabird Scientific, United States) that was fitted with 20 L Tefloncoated Go-Flo bottles (General Oceanic, United States) mounted on a CTD assembly. The daily cycle of incident PAR above the sea surface was measured on deck using a 4p scalar irradiance sensor (QSR-2100, Biospherical Instruments Inc., United States). The euphotic zone depth was defined as the depth of 1% of the surface PAR remains (Kirk, 1994). On deck, incubations of photosynthetic-irradiance (P B -E curve) experiments were carried out during the daytime at 2-3 depths to obtain the primary productivity at each depth. Primary productivity obtained from the on deck P B -E curve experiment is described in detail elsewhere (Gong et al., 2003). Euphotic zone integrated primary production (IP) was calculated using a trapezoid rule. The satellite sea surface temperature data were from the group for high resolution sea surface temperature (GHRSST) project, NASA and the Ocean Data Bank of the Ministry of Science and Technology, R.O.C. 1

Chemical Analysis
Water samples for nutrient analysis were placed in 100 ml polypropylene bottles and frozen immediately with liquid nitrogen. The methods employed for nutrient analysis are described in detail elsewhere (Gong et al., 1995a(Gong et al., , 1996(Gong et al., , 1997. The precision in the measurements of nitrate (NO 3 − ), phosphate (PO 4 3− ), and silicate (SiO 4 4− ) was ± 0.3, ± 0.01 and ± 0.1 µmol L −1 , respectively. Water samples for Chl a analysis were immediately filtered through GF/F filter paper (Whatman, 47 mm) and stored at -20 • C. A Turner Designs model 10-AU-005 fluorometer (Turner Inc., United States) was utilized to measure the Chl a concentration following extraction by 90% acetone (Strickland and Parsons, 1972;Gong et al., 1993Gong et al., , 1995b.

RESULTS
Marine hydrographic and chemical data as well as Chl a and primary production were measured during four seasonal cruises in the Taiwan Strait. Detailed information of the cruises and the euphotic zone results for averaged hydrographic data, nutrient concentrations, Chl a concentrations and IP at each station are listed in Table 1. Figure 2 presents the daily wind observation data from the Central Weather Bureau of Taiwan at Penghu Weather Station in 2017. The incident PAR daily cycles for the four seasonal cruises are shown in Figure 3. The results of marine environmental conditions and primary production of the four seasons in 2017 are described below.

Marine Environmental Conditions and Primary Production in Winter
The prevailing northeast monsoon was observed in winter (Figure 2). The PAR daily cycle in winter (Figure 3) indicated that the insolation duration (also known as day length) was approximately 10.5 h and the daily total PAR was 33.5 Einstein m −2 day −1 ( Table 1). Seawater temperature and salinity ranged 1 | Seasonal variation of day length, total photosynthetic active radiation (PAR) above the sea surface, euphotic zone depth (Ze), euphotic zone averaged marine environmental data, and primary production for each station.

Season
St. St., station no.; Long, longitude ( • E); Lat, latitude ( • N); BD, water depth (m); DL, day length (hour); E T , total PAR above sea surface (Einstein m −2 day −1 ); Ze, euphotic zone depth (m); T, temperature ( • C); S, salinity; TM, light transmission (%); NO 3 − , nitrate (µmol L −1 ); PO 4 3− , phosphate (µmol L −1 ); SiO 4 4− , silicate (µmol L −1 ); Chl a, chlorophyll a (mg m −3 ); IB, euphotic zone integrated Chl a (mg m −2 ); IP, euphotic zone integrated primary productivity (mgC m −2 day −1 ). from 19.02 to 23.29 • C and 32.41 to 34.40, respectively ( Table 1). The lowest temperature and salinity values were observed in station 3 while the highest temperature and salinity values were observed near the shore (station 4; Figure 4A). The cross-sections of the Taiwan Strait (Figures 5a-d) display as vertically uniform, indicating a strong vertical mixing within the prevailing northeast monsoon. The euphotic zone depth ranged between 5 and 22 m and decreased from offshore to inshore (Figure 5b). Station 5 had the lowest TM, which was under 30% throughout the whole water column. A low temperature and salinity with high nutrients were observed at stations 2 and 3 while a high temperature and salinity with relatively low nutrients were observed at station 4 (Figures 5a,c,d). Chl a concentrations ranged between 0.51 and 0.80 mg m −3 (Figure 5d) and had the same pattern as the nutrient distributions. The euphotic zone integrated Chl a concentration (IB) and IP varied ranges were 3.4-17.6 mg m −2 and 65-234 mgC m −2 day −1 , respectively ( Table 1). The highest IB and IP were at station 2 while the lowest ones were at station 5 ( Figure 6).

Marine Environmental Conditions and Primary Production in Spring
Wind speed began to slow and changed direction in early April (Figure 2). The PAR daily cycle in spring (Figure 3) indicated that the insolation duration was approximately 12.5 h and the daily total PAR was 59.5 Einstein m −2 day −1 ( Table 1). Seawater temperature and salinity ranged from 23.11-24.10 • C to 34.57-34.61, respectively ( Table 1). The cross-sections of the Taiwan Strait (Figures 5e-h) show that water columns were thoroughly mixed vertically. Seawater in all sampling stations was of a high salinity with similar temperatures (Figure 4B). The euphotic zone depth ranged between 5 and 17 m and decreased from offshore to inshore, similar to TM (Figure 5f). There were lower NO 3 − concentrations in the water columns, except for some NO 3 − available in the bottom water of station 2, and relatively high PO 4 3− concentrations were observed at station 5 (Figure 5g). SiO 4 4− concentrations decreased from offshore to inshore stations, and higher concentrations of Chl a were observed at station 2 (Figure 5h). IB and IP varied between 5.0-20.1 mg m −2 and 108-383 mgC m −2 day −1 , respectively ( Table 1). The highest IP was found at station 2 whereas the lowest IP was at station 5 because of its low IB, TM and shallow euphotic zone depth (Figure 6).

Marine Environmental Conditions and Primary Production in Summer
In summer, wind switched to the southwest direction and the wind speed slowed (Figure 2). The PAR daily cycle in summer (Figure 3) indicated that the insolation duration was approximately 13.5 h and the daily total PAR was 70.2 Einstein m −2 day −1 ( Table 1). Seawater temperature and salinity ranged from 27.98-29.91 • C to 33.43-34.14, respectively ( Table 1). Water columns were thoroughly mixed vertically in stations 3, 4, and 5 ( Figure 4C). On the other hand, the water columns were stratified and the seawater temperature decreased as the salinity increased with depth at stations 1 and 2 (Figure 5i). The euphotic zone depth ranged between 10 and 39 m and TM values were high in all the water columns; both of them decreased from offshore to inshore (Figure 5j). Relatively high nutrient concentrations were found in the bottom water of station 1 (Figure 5k). Higher concentrations of Chl a were observed at stations 1, 2, and 3, especially in the subsurface layer (Figure 5l). IB and IP showed ranges of 6.2-39.8 mg m −2 and 331-985 mgC m −2 day −1 , respectively ( Table 1). The highest IP was at station 2 whereas the lowest IP was at station 4 ( Figure 6).

Marine Environmental Conditions and Primary Production in Autumn
The northeast monsoon picked up after October (Figure 2). The PAR daily cycle in autumn (Figure 3) indicated that the insolation duration was approximately 12 h and the daily total PAR was 69.2 Einstein m −2 day −1 ( Table 1). Seawater temperature and salinity ranged from 27.34-29.55 • C and 33.61-34.02, respectively ( Table 1). The water columns were stratified at stations 1 and 2 but were thoroughly mixed vertically in stations 3 and 5 (Figures 4D, 5m). The euphotic zone depth varied between 23 and 52 m, and the TM values were high in all the water columns (Figure 5n). In most of the water columns, nutrient concentrations were low during the autumn, except for the bottom water of station 1 (Figure 5o). Although Chl a concentrations were relatively low in autumn (Figure 5p), the TM values were high and the euphotic zone depths were deep in all the sampling stations. As the result, IB and IP were high with ranges of 16.2-28.5 mg m −2 and 175-472 mgC m −2 day −1 , respectively ( Table 1). The highest IP was at station 1 while the lowest IP was at station 4 ( Figure 6).

DISCUSSION
The objective of this research is to understand the seasonal variations of marine environment and primary production in the Taiwan Strait. As primary production is a limiting factor of fisheries' catch rates (Chassot et al., 2010), the results of this research can not only be used as the baseline for evaluating the OWF impact on the marine ecosystem but also be used for assessing their influence on fishery resources. Three environmental factors, which are light availability, water column stability and nutrients availability, influence primary production the most. Although assemblage composition of phytoplankton as well as zooplankton grazing pressure (Tsai et al., 2018) also influence the primary production, these were not included in this study. According to the research results, IP not only varied between each sampling station in the same season but also differed among the four seasons ( Table 1). The highest averaged IP (664 ± 270 mgC m −2 day −1 ) was observed in summer and was approximately 5 times higher than that in winter (137 ± 68 mgC m −2 day −1 ; Table 2). Consequently, three environmental factors of the Taiwan Strait during four seasons were analyzed to understand the key factors that influenced the IP seasonal variation.

Light Availability
Solar radiation is the main external source of energy for marine ecosystem (Chassot et al., 2010). Primary production has fixed the solar energy via phytoplankton in the euphotic zone. Therefore, insolation duration, light intensity and transparency of the water columns influence the primary production. The seasonal insolation durations (Figure 3) ranging from the highest to the lowest were summer (13.5 h), spring (12.5 h), autumn (12 h), and winter (10.5 h). In addition, the seasonal mean values of the daily total PAR ranging from the highest to the lowest were 70.2, 69.2, 59.5 and 33.5 Einstein m −2 day −1 in summer, autumn, spring and winter (Table 1), respectively. Moreover, the seasonal mean TM ranging from the highest to the lowest were 95.2, 85, 69.5, and 65.5% in autumn, summer, spring, and winter ( Table 2), respectively. The TM in summer being lower than in autumn might be due to more phytoplankton present in the water columns (higher Chl a concentrations) in summer. Hence, light availability was better in summer and autumn than in winter and spring ( Table 3).

Water Column Stability
Generally, wind stress is strong during the winter in the Taiwan Strait (Jan et al., 2002), and it induces violent waves and mixing while reducing TM and the euphotic zone depth. The euphotic zone depth value is a measure of water clarity, which is an  important parameter regarding the health of ecosystem (Lee et al., 2007). The inshore station (station 5) had the shallowest water depth, so the water column was easily affected by terrestrial inputs or resuspend particles from the sediment, especially during the winter (Naik and Chen, 2008). As a result, the relatively shallow euphotic zone depth was observed at station 5 throughout all seasons. Figure 2 shows that strong northeastern winds had started since January and slowed down after April. The wind changed its direction to southwestern between June and August, switched back to a northeastern direction and then increased speed after October. According to the data, which is presented as a dashed line in the cross-section figures (Figures 5b,f,j,n), a lower TM and shallower euphotic zone depth were found in winter and spring due to strong vertical mixing, which indicated that the water columns were unstable. On the other hand, the water columns were thoroughly stratified in summer and autumn (Figures 5i,m). Therefore, water column stability was higher in summer and autumn than in winter and spring ( Table 3).

Water Masses Variety and Nutrients Availability
As Taiwan Strait is located in the East Asian monsoon region, seasonal changes of the wind direction have caused different currents and circulations. This have brought different water masses into the Taiwan Strait (Jan et al., 2002(Jan et al., , 2006Chen, 2003Chen, , 2011Chen and Sheu, 2006;Wu et al., 2007;Chen et al., 2010). In winter, the strong northeastern winds had brought the MCCW with cold, brackish and high nutrients to the western part of the Taiwan Strait and the CYR blocks part of the MCCW and forces a U-shaped flow pattern in the northern Strait (Jan et al., 2002(Jan et al., , 2006Wu et al., 2007). On the other hand, the intruding KBW flowed northward into the CYR through the Penghu Channel and partially blocks along-strait transport on the eastern side of the strait (Jan et al., 2002). The low temperature, salinity and high NO 3 − concentrations water existed in station 3 while high temperature, salinity and low NO 3 − concentrations water was in station 4 ( Figures 4A, 7). In winter, water west of 119.8 • E (station 1, 2, and 3) had high nutrient concentrations with a low temperature and salinity, which might have come from MCCW; the water east of 119.8 • E (stations 4 and 5) had relatively low nutrient concentrations with a high temperature and salinity (Figures 4A, 6A, 7 and Table 2), which might have come from KBW. The MCCW confronted the KBW in the CYR in winter. Water was well mixed vertically (Figures 4A, 5a) and it was the mixture of the MCCW and the KBW horizontally. As the northeastern wind began to slow and changed direction in spring, the MCCW retreated and the KBW expanded northward. Water in the Taiwan Strait had a higher salinity in spring than in other seasons (Figure 4) and it had low NO 3 − concentrations (Figures 6A, 7), meaning that it presumably came from the KBW. In summer and autumn, water in the Taiwan Strait had medium salinity and low nutrient concentrations (Figures 4C,D, 6A, 7), meaning that it presumably came from SCSW. The seasonal variations of nutrients are strongly influenced by the seasonal change of different water mass occupation. Consequently, the seasonal change of water masses in the Taiwan Strait in this research is in accordance with previous research (Jan et al., 2002(Jan et al., , 2006Chen, 2003Chen, , 2011Chen and Sheu, 2006;Wu et al., 2007;Chen et al., 2010;Huang et al., 2018). Hong et al. (2011) revealed that there are three upwelling areas in the southern Taiwan Strait in summer. Upwelling near the Chinese coast occurs during the summer and depends on the prevailing southwesterly winds, while upwelling around the Taiwan Bank and the Penghu Islands occurs throughout the year, which is due to the interaction between the bottom topography and the northward flow (Hong et al., 2011). According to Jan   2010), the acoustic Doppler current profiler (ADCP) data indicated that the northeastward transport increases in spring and summer, peaks in early August and weakens thereafter. In addition, typhoon Roke was formed southeast of Taiwan on July 21 and traveled northwest and passed through the Luzon Strait, reaching Hong Kong on July 23. The typhoon might have induced upwelling, and the upwelling may have been enhanced by the cyclonic vorticity (Ko et al., 2016). Southwesterly winds changed direction after October (Figure 2), which may have slowed down the northward flow and upwelling. According to the satellite sea surface temperature (SST) images (Figure 8), the area near Penghu Island and Taiwan Bank had lower SST in summer than in autumn. The temperature gradient between low SST areas and the sea areas nearby was larger in summer than in autumn, which indicates that the upwelling might be stronger in summer than in autumn. The northward currents carried the nutrientrich seawater from the upwelling area to our sampling area. As a result, water in summer had higher nutrient concentrations than in autumn ( Figure 6A). Consequently, nutrients' availability was high in winter, medium in summer and low in autumn and spring ( Table 3).

Seasonal Distribution of the Euphotic Zone Integrated Chl a Concentration and Primary Production
The IB were high in spring and summer, but they were low in winter and lower in autumn ( Table 2). During the spring and summer, Chl a concentrations were higher in the western part than the eastern part of the Taiwan Strait (Figures 5h,l), and the highest IB was revealed at station 2 (Table 1 and Figure 6B).
Although the average Chl a concentrations in autumn were the lowest, the euphotic zone depths were the deepest among all seasons, so the values of IB were relatively high ( Table 2 and Figure 6B). The seasonal mean IB values from highest to lowest were summer (21.2 mg m −2 ), autumn (19.9 mg m −2 ), spring (10.2 mg m −2 ) and winter (9.7 mg m −2 ; Table 2).
The highest IB and IP values were observed in station 2 in summer and the lowest IP value was found at station 5 in winter ( Figure 6B). The lowest IP values were observed in winter ( Figure 6B). Even if the water in winter was rich in nutrients (Figures 5c,d), a short insolation duration (Figure 3), low incident PAR and shallow euphotic zone depth (Table 1) due to strong vertical mixing resulted in the lowest IP values among all seasons. The highest IP values were observed in summer ( Table 2) at approximately five times higher than that in winter. Even though the nutrient concentrations were relatively low in summer (Figures 5k,l), Chl a concentrations were high (Table 1), especially at station 1, 2, and 3 (Figure 5l).  Hu et al. (2015) also found the water in the coastal upwelling system in the Taiwan Strait had low nutrients but high Chl a concentrations. Because The insolation duration (Figure 3), incident PAR, Chl a concentrations and IB values were high in summer (Table 1), the highest IP values were revealed. In addition, the high values of IP in summer were presumably supported by the heterotrophic nutrients, which might have been supplied by the enhanced upwelling in the south. Overall, seasonal IP values ranged from highest to lowest in the following order: summer > autumn > spring > winter.
If we assume that the maximum IP occurred in mid-July (i.e., day 199) and that the winter value represented the base of the IP in the Taiwan Strait, the annual IP variations can be described with a left-skewed normal distribution curve (Figure 9). The fitted left-skewed normal distribution curve is: where µ = 199, σ = 52 and X is the Julian day. The mean absolute relative error between the observed value and the value computed from the equation was approximately 3%. The average annual IP determined from the equation for the Taiwan Strait is approximately 120 gC m −2 year −1 (or 328 mgC m −2 day −1 ), which is comparable with the average annual IP 123 ± 86 gC m −2 year −1 (338 ± 235 mgC m −2 day − 1) from the measurement.
Nevertheless, the wind speed slowed down (Figure 2) after our spring sampling day, which means that the signal for spring might show up after our sampling day. Thus, the IP values for spring might be underestimated. With more stability of the marine environment, the IP values in spring may have increased. In addition, typhoons may not only enhance the upwelling in the sea but also bring the rainfall on the land, which increases the freshwater input and nutrient supplies from the land to the sea. Conclusively, the IP values for summer may increase when more typhoons approach.

Potential Impacts of Offshore Wind Farm Installations on Primary Production
In recent decades, the number of OWFs has increased rapidly. With increasing OWF constructions, environmental impact assessments and monitoring programs have been deployed worldwide. However, they mainly focus on sea mammals, seabirds, benthos or demersal fish but generally ignore the potential effects that OWFs may have on the marine ecosystem (Floeter et al., 2017). As the marine environment is home for marine creatures, a change in the marine ecosystem will have major impacts on them. Primary production is the foundation for a marine ecosystem and supports the food web and fish stock. Therefore, potential OWF installation impacts on the marine environment, especially primary production, should be addressed.
According to Smyth et al. (2015), OWF installation might cause several environmental pressures associated with foundations, piles, cabling routes, waste from vessels and disturbance from vessels. Carpenter et al. (2016) indicated that as tidal currents move past the OWF foundation structures, they generate a turbulent wake that will contribute to a mixing of the stratified water column. In addition, van der Molen et al. (2014) suggested that primary production may be decreased by higher turbulence levels induced by turbine foundations. Vanhellemont and Ruddick (2014) also found that the underwater light field will be affected by increased suspended particles and the turbid wakes may significantly impact sediment transport and downstream sedimentation. Even though OWF installation may increase the vertical mixing and bring nutrients to the upper layer, it also increases suspended particles and decreases light penetration. In addition, surface sediments of the CYR are composed of sands (Liao et al., 2008), which may be disturbed during the OWF installation and thus increase the suspended particles in the seawater. As vertical mixing in the Taiwan Strait was strong and water stability was already low during the winter and spring (Table 3), the impacts on IP due to OWF installation during these seasons may be limited. On the other hand, the impacts of OWF installation on IP during the summer and autumn may be greater. Figure 9 shows the annual IP variation curve in the Taiwan Strait. If we assume the OWF installation takes 90 days and it begins on the 1st of June (day 152), approximately 46% (55.6 gC m −2 year −1 ) of the total annual IP may be influenced. Since the highest IP was in summer, the potential impact of OWF installation on IP reduction may be larger during the summer than during other seasons. Although it can be argued that this IP variation curve only came from data in one year, it can be improved when the data set accumulated. This annual IP variation curve not only gives an insight into how the seasonal IP changed but also takes a glimpse into how the IP might be affected by OWF installation, which could provide suggestive information for governance, managers, operators and stakeholders to reduce potential impacts of OWF installation on the marine environment.
To reduce potential impacts on the marine environment and gain social benefits from ecosystem services, the investigation or monitoring of the potential effects of OWFs on the marine ecosystem is needed. Marine environmental investigations and/or OWF monitoring can provide the necessary scientific data for managers to reduce the potential impacts on the marine ecosystem as well as inform governmental authorities, operators, operating companies and other stakeholders about the stage changes of the marine environment.

CONCLUSION
During four seasonal research cruises in the Taiwan Strait, the average measured IP was 338 ± 235 mgC m −2 day −1 . As three primary water masses exist in the Taiwan Strait and three of them have different characteristics, different mixing ratios of water masses may cause different chemical and hydrographic conditions, which leads to different levels of Chl a concentrations and primary production. Although water in winter was rich in nutrients, it had a short insolation duration, low incident PAR and shallow euphotic zone depth. As a result, winter had the lowest IP values (137 ± 68 mgC m −2 day −1 ) among the four seasons. Owing to the strong vertical mixing in winter and spring, the euphotic zone depths were shallower in winter and spring than in summer and autumn, and it caused low IP values during winter and spring. The highest values of IP were observed in summer, which was 664 ± 270 mgC m −2 day −1 . This might be supported by the heterotrophic nutrients from upwelling around the Taiwan Bank and the Penghu Islands. The primary impact of an OWF installation on the marine environment is increasing disturbance of the seawater in addition to noise. An increase in disturbance may decrease the euphotic zone depth, which reduces IP. As naturally vertical mixing is strong during the winter and spring, the impact on IP from an OWF installation is minor during these seasons than others. On the other hand, the potential impact of an OWF installation on IP reduction may be larger during the summer than in other seasons. Consequently, long-term marine environmental monitoring can provide information and may reduce the potential environmental impacts as well as maintain the sustainability of both the social economy and healthy marine environment.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.

AUTHOR CONTRIBUTIONS
G-CG was the project leader and developed the idea of the study. C-CC, A-YT, and T-YC led the research cruises and did the sampling. W-LY and WH did the samples and data analysis. K-WL obtained and calculated data related to satellite. H-CT wrote the first draft of the manuscript. All authors contributed equally to the interpretation of the results and in writing the final manuscript.