The South China Sea Operational Forecasting System (SCSOFS) is an operational system developed and operated by the National Marine Environment Forecasting Center (Zhu et al., 2020, Zhu et al., 2022). The model setups, including the atmospheric forcings, river forcings, and boundary conditions, can be found in Zhu et al. (2022). The system utilizes a three-dimensional hydrodynamic-biogeochemical coupled model with a horizontal resolution of 1/30° and a vertical resolution of 50 layers of the σ layer. The coverage area of the system spans from 4.5°S to 28.3°N and from 99°E to 144°E. The ecosystem model of the SCSOFS is based on the carbon, silicate, and nitrogen ecosystem (CoSiNE) model, which is coupled with the Regional Ocean Modeling System (ROMS) (Shchepetkin and McWilliams, 2005). The CoSiNE model incorporates 15 state variables, including four nutrients [e.g., ammonium (NH4), nitrate (NO3), phosphate (PO4), and silicate (SiOH4)]; two phytoplankton groups [e.g., small phytoplankton (S1) and diatoms (S2)]; the corresponding Chla components [e.g., Chl1 and Chl2); two grazers [e.g., microzooplankton (Z1) and mesozooplankton (Z2)]; two classes of detritus (e.g., nitrogen (DN) and silicon (DSi)); and three carbonate system state variables: dissolved oxygen (DOX), dissolved inorganic carbon (DIC), and total alkalinity (TALK). These components are influenced by physical processes such as advection and diffusion. The CoSiNE model was initially developed to simulate marine geochemical processes in the equatorial Pacific Ocean (Chai et al., 2002) and has been applied at various spatial scales, including the North Pacific Ocean (Xiu and Chai, 2012), SCS (Pan et al., 2017) and PRE (Fang et al., 2022). Additionally, the model has been applied to study the responses of marine ecosystems to typhoon events (Pan et al., 2017; Fang et al., 2022). The parameter settings of the biogeochemical model for this study can be found in Table 1 . The study area lies within the range of 19.2°N to 24.2°N, 110.2°E to 119.2°E and includes the Qiongdong upwelling and Guangdong upwelling areas ( Figure 1 ). The two upwelling regions are predominantly influenced by the summer southwest monsoon, which typically occurs from late April to October. These upwellings also have important biogeochemical impacts on the region (Wu and Li, 2003; Gan et al., 2010).

Information regarding the two typhoons that affected the PRE in 2021 was obtained from the Shanghai Typhoon Research Institute, China Meteorological Administration (https://tcdata.typhoon.org.cn). These data provide a comprehensive overview of typhoons, including detailed information on important parameters such as the maximum wind speeds and minimum central pressures. The data are available at temporal resolutions of three or six hours, allowing for a more detailed analysis.

Typhoon Cempaka originated as a tropical depression in the northern SCS ( Figure 1 , purple dashed line). The typhoon underwent intensification and reached tropical storm status at 00:00 on July 19 (UTC). Subsequently, the typhoon further strengthened and became a strong tropical storm by 9:00 on the same day. By 15:00, it had reached the typhoon level and made landfall off the coast of Guangdong Province at approximately 13:50 on July 20. It weakened to a tropical depression in western Guangdong at 08:00 on the 21st and then moved through western Guangdong into southern Guangxi. On the morning of the 22nd, the typhoon changed direction to a southwesterly trajectory and entered the Gulf of Tonkin while continuously weakening. The second typhoon, Lupit, originated as a tropical depression on the eastern surface of the Leizhou Peninsula at 9:00 on August 2. By 00:00 on August 4, it had intensified into a tropical storm. Gradually moving northeastward, it eventually made landfall near the east coast of Guangdong Province at approximately 3:00 on August 5, a clear distance away from the PRE ( Figure 1 , blue dashed line). More information on both typhoons is shown in Table 2 .

The atmospheric forcing data used in this study were obtained from the Climate Forecast System Reanalysis (CFSR) from the National Centre for Atmospheric Research (https://psl.noaa.gov/data). The temporal resolution of the data is 6 hours, and the spatial resolution is 0.2-0.3°.

The study utilized satellite remote sensing data for Chla that were obtained from MODIS-Aqua, which were provided by the National Aeronautics and Space Administration (NASA) Ocean Color Processing Center (OCDPS, http://oceancolor.gsfc.nasa.gov), with a horizontal resolution of 4 km. Due to the interference from clouds during typhoons, the availability of local remote sensing data is limited (Zheng and Tang, 2007a). To overcome this limitation, this study employs an 8-day average synthetic product derived from the satellite data, which provides a temporal resolution of 8 days before and after a typhoon event. By analyzing these data, this study aimed to gain insights into the impact of typhoons on the Chla distributions in marine ecosystems.

During the analyses of the ecological environment before and after the typhoons, the period was divided into four phases to better understand the changes. The analysis period started one week before the arrival of Typhoon Cempaka and extended to three weeks after its passage. The four distinct phases are as follows:

The wind stress curl-induced EPV was calculated as follows (Price, 1981) (Equations 1–3):

where f=2ωsinθ is the Coriolis parameter, ω is the angular velocity of the earth’s rotation, ρ ₀ is the sea water density, τ → is the wind stress, ρ _a is the air density, U → is the wind speed at 10 m, and C _D is the drag coefficient.

Calculation of Chla budget by using the following model (Equation 4):

where ∂ Chla ∂ t   represents the temporal rate of Chla change (RATE), − ( u · ∂ Chla ∂ x + v · ∂  Chla  ∂ y ) represents the horizontal advection term (HADV), − w · ∂  Chla  ∂ z represents the vertical advection term, ∂ ∂ z ( K D · ∂ Chla ∂ z ) represents the vertical diffusion term (VDIFF), and Bio represents the term caused by biological processes. The Bio term of Chla is directly influenced by the abundance of the two phytoplankton groups, which are controlled by the growth, grazed, and mortality terms. The diatom also includes the sedimentation term. The horizontal diffusion term is neglected in this context due to its much smaller magnitude compared to the other terms.

To assess the ability of the ROMS-CoSiNE coupled model to simulate the changes in the seawater environment during the typhoons, field survey data from July 27 to August 9 (Guo et al., 2023) were utilized for salinity verification ( Figure 2 ). The model outputs exhibited slightly greater values than did the measured data, and the correlation coefficient was determined to be 0.85, with a root mean square error of 2.88. Overall, the simulation effectively captured the salinity trends during this period.

To assess the impact of Typhoon Cempaka on the Chla distributions, we calculated the differences between the composite data from July 20 to 27 and the composite data from July 12 to 19 ( Figure 3A ). The results showed that a distinct Chla bloom area was located west of the PRE, where the Chla concentrations increased by more than 2 mg·m^-3 ( Figure 3A ). Conversely, the region near the shoreline exhibited a decrease in Chla concentrations. In addition, a reduction in the Chla concentrations was observed in the vicinity of the upwelling area in Qiongdong. Similarly, the impact of Typhoon Lupit was determined by subtracting the data from August 5 to 12 from the data between July 28 and August 4 ( Figure 3B ). Unfortunately, the impact of Typhoon Lupit could not be observed using the satellite remote sensing data because of the cloud cover ( Figure 3B ). The model results were subsequently used to estimate the effects of the typhoons on the distributions of the surface Chla concentrations in the same way ( Figures 3C, D ). The Chla distributions before and after the storms are shown in Supplementary Figure S1 . The variation in simulated chlorophyll concentrations are lower than the satellite observations ( Figures 3A, C ), but there is a significant correlation (R=0.43, P<0.01). The RMSE for the Chla differences was 1.22 mg·m^-3, indicating that the simulation of the distribution of the chlorophyll concentrations in the surface layer during this period was satisfactory overall. Moreover, considering the poor observation results available from satellite remote sensing in nearshore turbid waters, the simulation discrepancy is acceptable (Ye et al., 2020; Feng et al., 2022). Based on the model results, the Chla concentrations in the estuary decreased under the influence of Typhoon Lupit, and the Chla concentrations in the algal bloom area that were caused by Typhoon Cempaka also decreased. The Chla concentrations increased on the eastern side of the estuary, while a decrease was observed on the western side ( Figure 3D ), potentially attributed to the mixing effect induced by Typhoon Lupit. For further details on the simulation performance and additional verification of the SCSOFS data, please refer to Zhu et al. (2022) & Guo et al. (2023).

By decomposing the analysis period into four distinct phases, we can observe and study the changes in the ecological environment before, during, and after the typhoons, as well as the impacts of subsequent events such as algal bloom dissipation and the occurrence of Typhoon Lupit. Before and after the typhoons, southwesterly winds prevailed in the northern SCS during summer (Zhu et al., 2019). First, the southwesterly winds were enhanced ( Figures 4A–F ). During the pre-Cempaka period, the region experienced a period with a break in the monsoon with wind speeds mostly less than 6 m·s^-1. This phase was characterized by relatively calm conditions of the southeasterly winds ( Figure 4A ). During the transit-Cempaka period, the typhoon brought a cyclonic wind field and high wind speeds. The wind velocities near the typhoon center reached more than 15 m·s^-1 ( Figure 4B ). During the post-1-Cempaka period, the wind field was gradually dominated by southwesterly winds, with obviously higher wind speeds than those during the pre-typhoon conditions ( Figure 4C ). During the post-2-Cempaka period, including the passage of Typhoon Lupit, the wind field predominantly shifted in a southerly direction ( Figures 4D–F ).

​Heavy rainfall accompanied the typhoons (Jacob and Koblinsky, 2007), resulting in increased precipitation ( Figure 4G ) and decreased sea surface salinity (Wang and Zhang, 2021). Notably, the rainfall amounts averaged over the study region increased by more than 7 times during the typhoon period. Under the combination of enhanced southwesterly winds and increased rainfall, dynamic and biochemical processes may undergo extraordinary changes (Lin and Oey, 2016; Zhang et al., 2021). Therefore, we further investigated the responses of the marine environment to these ocean-atmospheric interactions.

To better understand the physical-biogeochemical changes that are caused by typhoons, we analyzed the changes in temperature, salinity, Chla concentration, and NO3 concentration from the perspective of three dimensions during the four stages ( Figure 5 ). A depth of typhoon-induced mixing greater than the climatological mixing layer depth (MLD) may affect the variable distributions (Jiang et al., 2023). During the transit-Cempaka period, significant negative anomalous sea temperatures were evident over most regions ( Figure 5B ). However, under the influence of the northeasterly wind located west of the typhoon pathway, the Qiongdong upwelling process was suppressed, with a positive anomalous value. Typhoons typically induce deepening of the upper mixed layer in the ocean, resulting in surface cooling and subsurface warming (Price, 1981). The heat dissipation through air-sea heat flux from the sea surface is less important than mixing (Zhang et al., 2021). During post1-Cempaka, an obvious cooling effect was observed in most areas, ranging from 0.1 to 2.8°C ( Figure 5C ). The upwelling regions exhibited different temperature variations, with a maximum temperature increase of 1.9°C, possibly due to the disruption of the original upwelling structure and mixing by the typhoon (Zhao et al., 2009). During the post-2-Cempaka period, the maximum cooling temperature recorded during this period reached 4.5°C, which can be attributed to the persistent cooling effect of Typhoon Cempaka, as well as the additional cooling caused by Typhoon Lupit, alongside the subsequent recovery of the upwelling system ( Figure 5D ).

Salinity serves as an effective tracer for identifying the extent of river freshwater (Zhou et al., 2015). In the case of typhoons, salinity analysis was used to study the impact of the typhoons on the estuary and surrounding areas. During the per-Cempaka period, the sea surface salinities in the estuary were low due to the influence of freshwater from the river. Under the influence of easterly winds, which induce westward flow, the freshwater extent on the western side surpassed that on the eastern side ( Figure 5E ). During the transit-Cempaka and post1-Cempaka periods, obvious decreases in surface salinity were observed on the western side of the estuary, with a maximum decrease of 8.3 ( Figures 5F, G ). During the post-2-Cempaka period, influenced by both southwesterly winds and Typhoon Lupit, the low-salinity waters became more concentrated on the eastern side of the PRE and indicated the movement and redistribution of the freshwater plume ( Figure 5H ).

The Chla and NO3 distributions exhibited similar patterns ( Figures 5I–P ). An increase in NO3 concentrations is fundamental and essential for the proliferation of phytoplankton in nutrient-poor SCSs (Pan et al., 2017). Before Typhoon Cempaka, the surface distributions of Chla and NO3 in the northern SCS demonstrated decreasing trends from the estuary toward the offshore regions ( Figures 5I, M ), which were mainly caused by the offshore decrease in the large amount of nutrients carried by rivers into the sea (Harrison et al., 2008). However, after Typhoon Cempaka, there was an obvious increase in Chla concentrations (>2.0 mg·m^-3) near the coastal areas and a notable decrease in the estuary and the eastern plume region ( Figure 5K ). The changes in Chla and NO3 concentrations exhibited opposite patterns compared to those in the salinity, suggesting an underlying correlation with the transport of water from rivers. Typhoon Cempaka caused high-value areas on the western side of the estuary in the surface layer ( Figures 5K, O ). With the influence of the southwesterly winds, these high Chla and NO3 concentrations gradually shifted toward the eastern side of the estuary ( Figures 5L, P ). In addition, the increased salinity in the estuary corresponded to the decreases in Chla and NO3 concentrations. In the Qiongdong upwelling and Yuedong upwelling areas, consistent warming was observed at different depths, accompanied by decreases in the Chla and NO3 concentrations. According to the modeling results, typhoon transit disrupted the structure of the original upwelling areas, which caused the water to mix with the surrounding water at different temperatures and nutrient levels.

The changes in marine variables along the cross section were further examined to improve the observations of the vertical variations. The modeling results of the element distributions along section DE ( Figure 3A , section DE) revealed several changes that occurred during and after the typhoon, from estuary point E to typhoon center point A ( Figure 6 ). Such a cross section allows effective observations of the freshwater effects.

During the transit-Cempaka period, a cold center appeared in the upper layer of the DE section ( Figure 6B ). This led to a decrease in temperature of 1-2°C in the upper layer. During the post-1-Cempaka period, the BD section warmed with the sinking of isotherms, while the CE section cooled with increasing isotherms ( Figure 6C ). The maximum cooling of approximately 2°C occurred between 10 and 20 m in the central part of section DE, and the temperatures in the CE section decreased during the subsequent two weeks ( Figure 6D ). During the post1-Cempaka period, there was an obvious decrease in salinity in the upper layer of the BC section ( Figure 6G ). The distributions of Chla, NO3, and salinity anomalies in the upper layers exhibited similar but opposite patterns ( Figures 6G, K, O ). During the post-2-Cempaka period, the salinity in the upper layer increased, while the Chla and NO3 concentrations decreased ( Figures 6H, L, P ). The section also revealed higher Chla and NO3 values near the bottom, but these values were not sufficient to account for the post-typhoon increase. These changes in the marine environment during the typhoon cycle indicate complex interactions between various factors, such as winds, river freshwater, upwelling, nutrient distributions, and phytoplankton growth. The factors contributing to these changes will be further discussed.

​The typhoons caused obvious changes in both the physical and ecological elements on the western side of the PRE. These changes included cooling, decreased salinity, and increased Chla and NO3 levels in the surface layer, followed by a subsequent recovery phase. Additionally, the distribution of environmental variables along section DE ( Figure 6 ) revealed that Point A (112.6°E, 21.2°N) underwent a very representative process. Specifically, the changes at Point A, which was positioned along the typhoon track on July 19 and representing the typhoon center, provide a comprehensive analysis of these changes ( Figure 7 ). The vertical temperature profile at point A demonstrates the impact of the typhoons, which exhibit a noticeable increase in the amount of colder water from greater depths and a maximum temperature decrease of 2.5°C at a depth of approximately 20 m ( Figure 7B ). This phenomenon is likely attributed to a combination of water mixing and air-sea heat flux (Zhang et al., 2021). Previous studies have shown that the range of influence of the PRP is mainly within a depth range of 0-15 m (Zhi et al., 2022). The monthly average runoff is utilized in the model because of the lack of daily measured data, thereby excluding the influence of high runoff resulting from rainfall events. Notably, the changes observed above 15 m indicate that obvious alterations in salinities, Chla concentrations, and NO3 concentrations occurred after Typhoon Cempaka ( Figures 7C, E, G ). The anomaly patterns indicate a decrease in salinity in the upper layer, which corresponds to increases in NO3 and Chla, suggesting that the controlling factors for the changes in the upper layer originate away from point A rather than from local upwelling ( Figures 7F, H ). Considering the change in the lower layers is not enough to provide the increment of the upper layer, it is likely that the vertical changes did not dominate the overall trend. The plume carried large amounts of nutrients and a high quantity of phytoplankton biomass, and these changes may be attributed to the influence of the typhoon on the redistribution of plume (Gan et al., 2010; Liu et al., 2021; Wang and Zhang, 2021). Moreover, a noteworthy observation is the temporal synchrony between the increase of Chla and NO3 concentrations in the upper layer, which deviates from previous research findings (Pan et al., 2017). In previous studies, it was believed that the vertical transport and mixing of nutrients led to the peak of nitrate, followed by the peak of chlorophyll a due to phytoplankton proliferation. The simultaneous peak of nitrate and chlorophyll a suggests that this is likely a result of direct water mass exchange by advection and mixing rather than phytoplankton proliferation.

To further determine the vertical transport contribution, the Ekman pumping velocity (EPV) is calculated to characterize the vertical transport, which can characterize the strength of ocean upwelling or subsidence flow caused by wind stress (Wang and Zhang, 2021). The depression at the center of the typhoon was caused by strong wind stresses due to the Ekman pumping effect in the upper ocean. As shown in Supplementary Figure S2 , the EPV had upward velocities near the typhoon center and downward velocities at the periphery.

The initial screening of the impact factors was conducted by using statistical correlation analysis methods to explore and understand the relationships among the ocean state variables under the influence of typhoons ( Figure 8 ). The variations in SSTs during the typhoon period had a negative relationship with the EPV, which was mainly attributed to the upwelling of cold water by Ekman pumping, a process driven by the strong winds associated with typhoons (Yang and Hong, 2021). Moreover, most of the ecological elements exhibited stronger correlations with salinity than with EPV ( Figure 8 ). The decreased salinity is mainly attributed to the extension of the PRP. These findings indicate that vertical variability may not play a dominant role in determining phytoplankton anomalies, while horizontal advection has a more important impact. Moreover, the alterations in ecological variables exhibited some degree of dissimilarity. The absence of a significant correlation between PO4 and salinity may be attributed to phosphorus limitation in the PRP and the impact of biological absorption (Gan et al., 2014). This suggests that the mixing caused by typhoons may have a positive impact on nutrient availability in certain areas. Another noteworthy finding is that there was a significant positive correlation between salinities and diatoms, as well as a rapid increase in Chla and small phytoplankton concentrations in the upper layer of point A after the typhoon, while the diatom concentrations decreased. Previous studies have shown that under sufficient light conditions and without nutrient limitations, diatoms tend to grow faster and dominate the small phytoplankton after typhoons (Chai et al., 2007; Ma et al., 2013; Pan et al., 2017). However, different observations have been made in coastal areas, with some studies reporting increases in dinoflagellate abundance and others reporting increases in diatom abundance and even cases of biological degradation (Tsuchiya et al., 2014; Anglès et al., 2015; Zhou et al., 2021; Thompson et al., 2023). The low diatom concentrations in the PRP set in the model likely influenced the change in algal species composition, further highlighting the importance of horizontal advection (Anglès et al., 2015).

In summary, the correlation analysis revealed a significant relationship between plume expansion, represented by salinity, and ecological variables. To further quantify the impact of the horizontal advection term, we examined the various contributions of the chlorophyll concentrations in the water column.

Shown as in Figure 9 , the results indicate high spatial and temporal consistency between the horizontal advection term and the total variability. The horizontal advection term is larger than the vertical term, which thus allows us to quantitatively determine the dominant role of horizontal advection, and even whether the horizontal advection exceeds the total variability. In general, the changes in the upper layer after typhoon Cempaka occurred were more obvious than those in the lower layer, and the effect of Typhoon Lupit was not obvious in the upper layer at this point. The RATEs of NO3 and Chla also exhibited considerable spatial and temporal consistency, which was also attributed to the dominant role of horizontal advection. Previous studies have noted that the post-typhoon, peak Chla levels typically occur after the peak NO3 levels, reflecting the proliferation of phytoplankton due to nutrient enrichment (Pan et al., 2017). However, in this bloom event, NO3 variability revealed synchronous peaks in both the Chla and NO3 concentrations, maintaining a nearly identical temporal pattern. Additionally, our examination of the NO3 uptake at this point indicated minimal changes, with a smaller degree of variability compared to the RATE of Chla. These findings suggest that the redistribution of water masses and nutrients due to advection is a mechanism that cannot be neglected during typhoon events and that complex interactions occur between physical and ecological processes.

After considering the dominant role of the advection process in driving the changes, the specific change process will now be examined. To better understand the factors that contribute to the increased Chla concentrations, the Chla concentration budget of a specific area within the bloom region ( Figure 3A , black box) was examined. Based on the previous findings, our primary focus was on the variations within the upper 0-15 m.

The results of the flux calculations reveal an obvious impact of typhoon transit on the changes in Chla concentrations. Consistent with the results for point A, the horizontal advection term dominates, the Chla RATE is highly consistent with the change in horizontal advection, and the contribution from vertical changes is relatively minor. Both typhoon events caused the changes in Chla concentrations in the region to be dominated by horizontal advection, which represented the impact of the typhoons on the PRP. However, according to the different movement directions and intensities of the two typhoons, the increase caused by Lupit’s cyclonic wind field was limited and was smaller than the increase caused by Cempaka; additionally, the increase was quickly dominated by the restored southwest monsoon. The northeast current, influenced by the summer southwest monsoon, causes the flow of the PRP toward the northeast of the estuary (Gan et al., 2010), and the passage of typhoons disrupts the flow dynamics, resulting in an amplified impact of the PRP on the western side of the estuary ( Supplementary Figure S3 ). These findings highlight the dominance of the advection process driven by typhoons in the occurrence of algal blooms, which dissipate as the flow field changes (Zhao et al., 2009; Fang et al., 2022; Feng et al., 2022). However, unlike the findings from previous studies, the temporal fluctuations in the Chla and NO3 concentrations appeared nearly synchronous, emphasizing direct Chla transport rather than a process of nutrient-induced phytoplankton growth. We observed that biological processes contributed negatively to the variability in Chla concentrations ( Figure 10A ), which persisted until the occurrence of Typhoon Lupit. Subsequently, we aimed to elucidate this occurrence based on the behavior of the ecological elements that were simulated by the model.

Before Typhoon Cempaka, small phytoplankton accounted for approximately 63% of the total phytoplankton biomass and reached a maximum proportion of 83% in the box region. This finding is consistent with the findings of Li et al. (2009) and Pan et al. (2017), who reported that under the influence of typhoons, small phytoplankton occupied 90% of the biomass on the coast and established initial dominance. Considering the negative correlation between S1 and S2 ( Figure 8 ), it is likely that there was potential competition for nutrients between the two species. The dynamics of the phytoplankton community are influenced by a combination of temperature, light, and nutrient concentrations (Xiao et al., 2018). Typhoon transit also induces a reduction in radiation flux ( Figure 11 ), and decreasing light intensity also impacts biological processes. Phytoplankton enrichment leads to an increase in light extinction, which is also not conducive to phytoplankton growth in the lower layers (Chai et al., 2002). The variations in NO3 uptake during the occurrence of the typhoons reflect this process ( Figure 10A , pink dashed line). The combined impact of radiation flux attenuation and light extinction that is caused by increased phytoplankton concentrations results in a reduction in the uptake rate of NO3 (Lao et al., 2023a). Until the passage of Typhoon Lupit, there was a notable increase in nitrate uptake, aligning with the positive contribution of biological effects on Chla. Specifically, this increase primarily resulted from the growth of small phytoplankton ( Figure 11A ). The negative biological term for diatoms persists throughout the period. Despite the fact that the model considers the sinking of large phytoplankton, this factor does not primarily account for this trend. As previously mentioned, prior to the typhoon, the region experienced a period with a break in the monsoon, which was not conducive to generation of upwelling and had an impact on the nutrient levels in the upper layer. This could account for the adverse contribution of biological terms to the Chla before the typhoon occurred. After Typhoon Cempaka, following the recovery of radiation flux, there was an increase in NO3 uptake along with PO4 absorption ( Figure 11B ), but the Bio values remained negative ( Figure 10A ). When phytoplankton are enriched, zooplankton are also enriched, and the increased feeding pressure of zooplankton could inhibit phytoplankton growth, although the excretions of zooplankton also serve as an important source of nutrients (Chen et al., 2013a; Zhang et al., 2018; Feng et al., 2022). In the CoSiNE model, grazing is the only source term for zooplankton and an important sink term for phytoplankton. Biological term of zooplankton is positive during Cempaka ( Figure 10 , yellow line), representing the grazing of zooplankton on phytoplankton. In general, the negative biological contribution can be attributed to phytoplankton’s inefficient use of increased nutrients, leading to their loss through grazing.

In brief, the advection induced by typhoons triggered complex biogeochemical processes. Following the typhoon, there was an obvious increase in nutrient levels, but the biogeochemical process did not result in a positive contribution of biological processes to Chla. This can be primarily attributed to the combined effects of light limitation and zooplankton feeding pressure. Therefore, this bloom was a simple physical process, and the increase in Chla was mainly due to direct phytoplankton transport. This finding contrasts with previous studies, which primarily attributed blooms to biological alterations (Zheng and Tang, 2007a; Zhao et al., 2009; Liu et al., 2013; Ning et al., 2019). It is believed that future improvements in observational techniques can help to verify these inferences and improve the performance of simulations. Moreover, due to the redistribution process caused by typhoons, the remineralization process of detritus is bound to be affected as well. The role of the changes in nearshore water turbidity during typhoon events is still unknown, and the resuspension process is also an important process that affects phytoplankton growth (Wang and Zhang, 2021; Zhang et al., 2021). Under the strong stirring effect of typhoons, a large amount of oxygen is injected into the interior of the ocean, which can cause the degradation of organic matter in the water column, thereby consuming a large amount of oxygen and releasing a large amount of nutrients, leading to the deterioration of the water, which is not conducive to the growth of phytoplankton (Wang et al., 2017; Lao et al., 2023a). Coupling these processes is believed to be helpful in accurately simulating the biogeochemical impacts of typhoons (Moriarty et al., 2017; Zang et al., 2020). Accurate simulation and prediction of biogeochemical effects caused by typhoons are of great significance to human production activities in coastal areas, such as mariculture (Lin and Lin, 2022; Zhang et al., 2023; Lao et al., 2023b).