Abstract
Intensive aquaculture usually produces large volumes of nutrient-rich wastewater, which is essential to treat to avoid eutrophication. This study aimed to evaluate the performance of five continuously aerated biofilm reactors treating simulated high-strength aquaculture wastewater under different dissolved oxygen (DO) levels, and the effects of DO increments on simultaneous nitrification-denitrification (SND). Continuous aeration was beneficial to complete nitrification. Total inorganic nitrogen (TIN), principally ammonium, was mainly removed by SND. The SND rate response to different DO levels was fitted well by the power function of y = 54.81 + 371.58/(1 − 0.16*x)^(−1/0.24) (R2 = 0.897, P = 0.000). When the TIN was removed completely, the optimal SND rate was defined and corresponded to a value of 121.8%. Accordingly, the optimal DO concentration was calculated as 2.10 mg/L, close to the actual level of 1.83 mg/L, at which the highest proportional removals of total nitrogen (58.0%) and TIN (57.3%) were obtained. Phosphorus was also removed by denitrifying polyphosphate-accumulating organisms.
HIGHLIGHTS
Dissolved ammonium was mainly removed by simultaneous nitrification-denitrification (SND) process.
The optimal SND rate for the current wastewater was 121.8%.
The optimal DO corresponding to the best SND rate was 2.10 mg/L.
Phosphorus was mainly removed by denitrifying phosphate-accumulating organisms (PAOs).
Graphical Abstract
INTRODUCTION
With rising living standards, demand for fish protein is increasing. Promoted by the demand for aquatic products, aquaculture expansion, usually involving intensive fish production, is inevitable (Calone et al. 2019). Huge amounts of wastewater are still generated during production (Hesni et al. 2020). Compared to natural surface waters, intensive aquaculture wastewaters usually contain high concentrations of organics, nitrogen (mainly ammonium), and phosphorus (von Ahnen et al. 2018). If discharged directly without proper treatment, such wastewater will cause serious eutrophication. Under increasing pressure for environmental protection, a more effective way to deal with this kind of mixed component wastewater is needed urgently.
Nitrogen removal by the conventional activated sludge process is usually achieved in sequential continuous reactors performing nitrification and denitrification separately (Gao et al. 2020). This requires a large area and has the disadvantages of complex operation (Velasco-Garduno et al. 2019), poor impact resistance (Li et al. 2019), and easy appearance of filamentous bulking (Li et al. 2020a). On the contrary, biofilm reactors might be an effective option for intensive aquaculture wastewater treatment (Davidson et al. 2019), offering the advantages of convenience, economy, large biomass, strong shock load resistance, and relatively lower sludge yield (Zhang et al. 2018; Chai et al. 2019). Due to diffusion limitations, anaerobic, anoxic, and aerobic areas may be distributed within the biofilm, forming numerous, similar ‘microreactors’ under aeration conditions (Cao et al. 2017). These might promote simultaneous nitrification-denitrification (SND) in a single reactor, eliminating the need for a separate anoxic compartment (Gao et al. 2020) and improving nitrogen removal efficiency (Chai et al. 2019). At the same time, due to the role of endogenous denitrification and phosphorus removal, such a reactor has a high carbon utilization rate and can achieve simultaneous nitrogen and phosphorus removal (Wang et al. 2015).
Many factors affect the SND rate of biofilm reactors, including the DO concentration, carbon-to-nitrogen (C/N) ratio, packing filler, salinity, etc (Machat et al. 2019; Xia et al. 2019; Zhang et al. 2020a). DO is one of the most critical factors because different DO concentrations affect the number and activity of various functional microbes in the reactor directly. For instance, when DO increases, the number of ammonium-oxidizing bacteria (AOB) and nitrite-oxidizing bacteria (NOB) increases (Cao et al. 2017), and the nitrifying and denitrifying capacities increase and decrease respectively (Zhang et al. 2020a). DO can also affect the proportions of aerobic and anaerobic areas inside the biofilm directly, and then the microbial community structure (Cao et al. 2017; Wang et al. 2018), making it one of the most important factors for enhancing the SND rate.
Recently, there have been inconsistent reports on the optimal DO to SND rate. Wang et al. (2020) and Cao et al. (2017) state that maximum TN/TIN removal was achieved in a moving-bed sequencing batch reactor (SBR) at a DO of 2.5 mg/L when treating synthetic ammonium wastewater. Peng et al. (2020) reported that the SND rate in an anaerobic granular SBR reached the maximum of 53.86% at a DO of 2 mg/L. Zhang et al. (2020b) reported that the SND rate achieved a maximum of 92.3% in an aerobic biofilm system at a DO of 2 mg/L, when treating real industrial/domestic wastewater. Subtil et al. (2019) found that TN removal achieved a maximum of 66% in a single-stage membrane bioreactor at a DO of 0.8 ± 0.1 mg/L. However, there is still little information about the response of SND rate to different DO levels in biofilm reactors for the treatment of high-strength aquaculture wastewater.
In this study a biofilm reactor was used to treat simulated high-strength aquaculture wastewater through SND, to explore the response of SND rate to different DO levels. Using linear fitting between DO and the SND rate, the aim was to determine the optimal level of DO in a biofilm reactor to treat this kind of organic-rich aquaculture wastewater through SND. Ordination analysis between performance and environmental factors was used to obtain the response mechanism of SND to different factors in the biofilm reactor, which would help guide practical application.
MATERIALS AND METHODS
System construction
For this study, five identical SBRs (R1, R2, R3, R4, R5) were constructed of Plexiglas – 20 cm (L and W) × 25 cm (H). They were supported by a 3.5 mm thick PVC plate with a 4.5% perforation rate, above which was 12.5 cm of lightweight, spherical ceramsite (Φ: 3–5 mm; porosity: 0.433). Two aeration sand trays (Φ: 8 cm) were fixed at the bottom of the supporting layer, where a 20 mm drainage hole was drilled. Liquid-level probes were installed near the bottom and top of the reactor, to control liquid volume (Figure 1).
Each SBR had its own, time-controlled aeration system – the two aeration trays in the supporting layer were connected to the same two-chamber micro-air pump (supply rate: 2 × 200 L/h; power: 4 W; maximum lift: 2 m; brand: EHEIM; model: EM3704) through 4 mm internal diameter rubber hoses fitted with check and regulating valves. During feeding, identical influent pumps simultaneously transferred the synthetic high-strength aquaculture wastewater from the same feed tank to the SBRs for subsequent reaction.
Experimental design
In the study, the SBRs were operated in sequential batch mode with the help of the time controllers and liquid-level relays. The four cycles per day each comprised 3 minutes feeding, 6 hours aeration, and 2 minutes drainage. Unlike conventional SBR operation, no sedimentation or idle stage was set. To explore the impact of DO on reactor performance, the DO levels in R1 to R5 were controlled at around 1, 2, 3, 4, and 5 mg/L (He et al. 2009), respectively, by adjusting the valve on the rubber aeration hose.
During start-up, each reactor was inoculated intermittently over 3 days with 3.5 L of activated sludge. The activated sludge came from a previous biofilter operated for simulated high-strength aquaculture wastewater treatment over a year. The reactors' mixed liquid suspended solids (MLSS) concentrations were maintained at approximately 2,000 mg/L. After inoculation, all reactors began treating the aquaculture wastewater for microbial acclimation until maturation, which was characterized by stable carbon, nitrogen, and phosphorus removal in the systems. The hydraulic loading rate (HLR) for the reactors was maintained at 0.35 m3/(m2·d) throughout.
Most pollutants in intensive aquaculture wastewater come from unused feed/excreta (Kokou & Fountoulaki 2018), so the aquaculture wastewater in this study was simulated by dissolving 20 g fish feed (product series: 156, from Tongwei Feed Co., Ltd) in 3 L tap water in a 5 L beaker with magnetic stirring for 5 minutes. After sedimentation for 15 minutes, the supernatant was supplemented with anhydrous sodium acetate (4.90 g), ammonium chloride (2.44), and potassium dihydrogen phosphate (1.09) as carbon, ammonium, and phosphorus sources, respectively, in a 75 L feed tank. The simulated wastewater's final nutrient compositions are shown in Table 1. Trace elements, as described by He et al. (2016), were also added to the simulated wastewater and are listed in Table 2.
Indices . | Sources . | Concentration in this study (mg/L) . | Grade V Standard (mg/L) . |
---|---|---|---|
DCOD | Fish feed + CH3COONa | 84 ± 13.1 | / |
TCOD | Fish feed + CH3COONa | 111 ± 12.4 | ≤40 |
IP | Fish feed + KH2PO4 | 4.1 ± 0.2 | / |
TP | Fish feed + KH2PO4 | 4.4 ± 0.2 | ≤0.4 |
NO3−-N | Fish feed | 1.9 ± 0.0 | / |
TAN | Fish feed + NH4Cl | 8.7 ± 0.2 | ≤2.0 |
NO2−-N | Fish feed | 0.0 ± 0.0 | / |
TIN | Fish feed + NH4Cl | 10.6 ± 0.2 | / |
TN | Fish feed + NH4Cl | 12 ± 0.6 | ≤2.0 |
Org-N | Fish feed | 1.4 ± 0.5 | / |
DCOD/TIN | 7.9 | / |
Indices . | Sources . | Concentration in this study (mg/L) . | Grade V Standard (mg/L) . |
---|---|---|---|
DCOD | Fish feed + CH3COONa | 84 ± 13.1 | / |
TCOD | Fish feed + CH3COONa | 111 ± 12.4 | ≤40 |
IP | Fish feed + KH2PO4 | 4.1 ± 0.2 | / |
TP | Fish feed + KH2PO4 | 4.4 ± 0.2 | ≤0.4 |
NO3−-N | Fish feed | 1.9 ± 0.0 | / |
TAN | Fish feed + NH4Cl | 8.7 ± 0.2 | ≤2.0 |
NO2−-N | Fish feed | 0.0 ± 0.0 | / |
TIN | Fish feed + NH4Cl | 10.6 ± 0.2 | / |
TN | Fish feed + NH4Cl | 12 ± 0.6 | ≤2.0 |
Org-N | Fish feed | 1.4 ± 0.5 | / |
DCOD/TIN | 7.9 | / |
Note: The Grade V Standard is the environmental quality standard for surface water in China (GB3838-2002) (Xu et al. 2021).
Item . | EDTA . | FeCl3 . | CaCl2 . | H3BO3 . | KI . | CuSO4·5H2O . | MnCl2·4H2O . | ZnSO4·7 H2O . | CoCl2·6H2O . | Na2MoO4·2 H2O . |
---|---|---|---|---|---|---|---|---|---|---|
Concentration (mg/L) | 10.0 | 0.9 | 5.0 | 0.15 | 0.18 | 0.03 | 0.06 | 0.12 | 0.15 | 0.06 |
Item . | EDTA . | FeCl3 . | CaCl2 . | H3BO3 . | KI . | CuSO4·5H2O . | MnCl2·4H2O . | ZnSO4·7 H2O . | CoCl2·6H2O . | Na2MoO4·2 H2O . |
---|---|---|---|---|---|---|---|---|---|---|
Concentration (mg/L) | 10.0 | 0.9 | 5.0 | 0.15 | 0.18 | 0.03 | 0.06 | 0.12 | 0.15 | 0.06 |
Data collection
Daily samples were collected of the raw influent (i.e. the simulated wastewater) and the 6-hour effluents from the reactors. The samples were analyzed for dissolved chemical oxygen demand (DCOD), total chemical oxygen demand (TCOD), nitrite nitrogen (NO2−-N), nitrate nitrogen (NO3−-N), total ammoniacal nitrogen (TAN), total nitrogen (TN), inorganic and total phosphorus (IP and TP) using standard methods (APHA 2017). TCOD and TP were determined directly, while DCOD and IP were measured after filtration through a 0.45 μm fiber membrane. Real-time temperature (T), DO, alkalinity (pH), redox potential (ORP), specific conductivity (SC), total dissolved solids (TDS), and salinity (Sal) were measured with a YSI 6600 V2 multiparametric sonde (Yellow Spring Instruments, USA).
Statistical analysis
One-way ANOVA (analysis of variance) was used to test the impact of DO level on reactor performance, with LSD (equal variances assumed) or Tamhane's T2 test (equal variance not assumed) for multiple comparisons. Linear fitting between the DO and SND rates was performed to explore the SND rate's response to different DO levels. Redundancy analysis (RDA) was also employed to evaluate performance and environmental factor relationships. All significance was set at P < 0.05.
RESULTS AND DISCUSSION
Comparison of online physicochemical parameters
Table 3 shows that the physicochemical parameters with obvious distinctions between the influent and effluent include T, DO, and ORP, while those with limited change include pH, EC, TDS, and salinity. Because the experiment was done in a greenhouse, the wastewater's temperature was significantly below that indoors, so, after a period of retention, the effluent temperature was significantly higher than that of the influent. The significant decrease in effluent DO concentration was attributed to its consumption arising from the oxidation of organics and ammonium (Machat et al. 2019). ORP concentration generally reflected those of DO and TAN (Lackner & Horn 2012). The effluent ORP from R1 and R2 was negative – that is, the reactor DO was relatively low and the TAN content relatively high.
Site . | T (°C) . | DO (mg·L−1) . | pH . | ORP (mV) . | EC (μs·cm−1) . | TDS (g·L−1) . | Sal (/wt‰) . |
---|---|---|---|---|---|---|---|
Influent | 11.5 ± 0.8 a | 11.41 ± 0.32 a | 7.58 ± 0.25 a | 206 ± 92 a | 447 ± 11 a | 0.29 ± 0.01 a | 0.22 ± 0.01 a |
Eff._R1 | 16.8 ± 0.8 b | 0.96 ± 0.30 b | 7.36 ± 0.11 b | −46 ± 37 b | 439 ± 9 a | 0.29 ± 0.01 a | 0.21 ± 0.01 a |
Eff._R2 | 17.1 ± 0.8 bc | 1.83 ± 0.39 c | 7.34 ± 0.12 b | −11 ± 27 bc | 396 ± 59 ab | 0.26 ± 0.04 ab | 0.19 ± 0.03 ab |
Eff._R3 | 17.7 ± 1.0 c | 2.93 ± 0.40 d | 7.33 ± 0.12 b | 7 ± 24 c | 368 ± 54 b | 0.24 ± 0.04 b | 0.18 ± 0.03 b |
Eff._R4 | 17.8 ± 0.9 c | 4.22 ± 0.29 e | 7.32 ± 0.12 b | 18 ± 21 c | 393 ± 13 b | 0.26 ± 0.01 b | 0.19 ± 0.01 b |
Eff._R5 | 17.8 ± 0.8 c | 5.19 ± 0.39 f | 7.32 ± 0.12 b | 26 ± 20 c | 408 ± 6 b | 0.26 ± 0.00 b | 0.20 ± 0.00 b |
P value | 0.000 | 0.000 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 |
Site . | T (°C) . | DO (mg·L−1) . | pH . | ORP (mV) . | EC (μs·cm−1) . | TDS (g·L−1) . | Sal (/wt‰) . |
---|---|---|---|---|---|---|---|
Influent | 11.5 ± 0.8 a | 11.41 ± 0.32 a | 7.58 ± 0.25 a | 206 ± 92 a | 447 ± 11 a | 0.29 ± 0.01 a | 0.22 ± 0.01 a |
Eff._R1 | 16.8 ± 0.8 b | 0.96 ± 0.30 b | 7.36 ± 0.11 b | −46 ± 37 b | 439 ± 9 a | 0.29 ± 0.01 a | 0.21 ± 0.01 a |
Eff._R2 | 17.1 ± 0.8 bc | 1.83 ± 0.39 c | 7.34 ± 0.12 b | −11 ± 27 bc | 396 ± 59 ab | 0.26 ± 0.04 ab | 0.19 ± 0.03 ab |
Eff._R3 | 17.7 ± 1.0 c | 2.93 ± 0.40 d | 7.33 ± 0.12 b | 7 ± 24 c | 368 ± 54 b | 0.24 ± 0.04 b | 0.18 ± 0.03 b |
Eff._R4 | 17.8 ± 0.9 c | 4.22 ± 0.29 e | 7.32 ± 0.12 b | 18 ± 21 c | 393 ± 13 b | 0.26 ± 0.01 b | 0.19 ± 0.01 b |
Eff._R5 | 17.8 ± 0.8 c | 5.19 ± 0.39 f | 7.32 ± 0.12 b | 26 ± 20 c | 408 ± 6 b | 0.26 ± 0.00 b | 0.20 ± 0.00 b |
P value | 0.000 | 0.000 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 |
Note: Significant differences in columns are marked with different letters.
There was no significant difference in pH among the effluents, perhaps because of the SND taking place. The alkalinity consumed by nitrification could be supplemented by that produced by denitrification (Wang et al. 2020). Nevertheless, the pH of all effluents decreased significantly compared with the influent, for which there are two possible reasons. Firstly, a small amount of sodium acetate was added to the influent as carbon source, which would increase the pH (Li et al. 2020b), and its exhaustion – indicated by the low effluent COD concentration, Figure 2(f) – in the reactor, made the effluent pH slightly lower than that of the influent. Secondly, nitrification consumed alkalinity and denitrification produced it, and nitrification in the reactor might be more efficient than denitrification, resulting in a slight decrease in pH (Li et al. 2014). This was especially obvious among reactors with high DO levels (Table 3).
Comparison of treatment performance
As shown in Figure 2, the proportional removal of most pollutants was significantly different among the reactors except for the Org-N index (Figure 2(c)). The latter probably because no stage was set for sludge discharge and suspended sludge was occasionally discharged with the effluent, which was influenced heavily by sampling flow speed. The suspended sludge was composed mainly of organics involving Org-N (Zhang et al. 2019).
The effluent nitrite concentrations were all extremely low (no more than 0.05 mg/L) (Figure 2(c)), so the statistical difference in accumulation rate was meaningless. Nevertheless, the lack of obvious nitrite accumulation in the effluent indicated that complete nitrification occurred in the reactors even under the relatively low DO conditions – e.g., in R1 and R2. Several studies show that continuous aeration is helpful to complete nitrification acclimation (Boog et al. 2018; Hou et al. 2018; Chai et al. 2019) even under low DO conditions (Subtil et al. 2019), which accords with the present study. Complete nitrification was also demonstrated by the obvious nitrate accumulation at 3 to 5 mg/L DO levels (Figure 2(b)).
With increasing DO, the proportional removal of TAN increased significantly (Figure 2(a)). When the DO concentration was 1 mg/L, almost no TAN was removed. When it was 2 mg/L, the TAN removal rate was below 60%. When it was increased to 3 to 5 mg/L, however, TAN was removed almost completely. Correspondingly, the NO3-N concentration increased gradually with increasing DO (Figure 2(b)). When DO was 1 to 2 mg/L, the NO3-N removal rate was close to 100%, indicating that both the influent NO3-N and the proportion of it arising from TAN oxidization were denitrified completely, which also implied that SND was occurring in the reactors.
The TN and TIN concentrations increased significantly at first and then decreased gradually as DO increased from 1 to 5 mg/L (Figure 2(d)). Because the TIN removal denoted the extent of denitrification (Zhang et al. 2020a), it was obvious that denitrification efficiency was relatively high when DO was between 1 and 2 mg/L. Since effluent TN was mainly TAN, nitrification was the dominant factor controlling TN removal. Meanwhile, combined with organics removal in R1 and R2 (Figure 2(f)), it could be inferred that increasing the DO concentration to an appropriate level should improve TN removal given the presence of sufficient organics, which could be stored temporarily as an endogenous carbon source (Wang et al. 2015). When DO increased from 3 to 5 mg/L, TAN was removed completely and nitrate accumulation increased, while the proportional removal of TN and TIN decreased slowly, indicating that complete nitrification and aerobic denitrification occurred inside these systems. TN and TIN could still be removed when TCOD and DCOD removal was similar and their effluent concentrations were low (Figure 2(f)). This suggested that the reactors could remove nitrogen by endogenous denitrification when the C/N ratio was low (Wang et al. 2015), hinting that the carbon source was not the only inhibiting factor for denitrification. It is thought that the decrease of the hypoxic zone due to increased oxygen penetration thickness in the biofilm caused by the increase in DO might be another decisive factor inhibiting denitrification (Cao et al. 2017; Wang et al. 2018).
In this study, phosphorus displayed a low treatment performance (<40%, Figure 2(e)). Because of the continuous aeration, conditions outside the reactor biofilms were always aerobic (i.e., in the liquid phase – Table 3), perhaps causing constant anoxic/anaerobic conditions inside the biofilm (Cao et al. 2017; Wang et al. 2018; Jiang et al. 2020), thus preventing alternating anaerobic/aerobic conditions over time. This would be detrimental to phosphorus release by conventional PAOs and restrain subsequent phosphorus uptake (Zhao et al. 2018). When treating high phosphate concentrations, the packing filler of the ceramsite used in the reactors could easily reach phosphate adsorption saturation (Zhang et al. 2019), so the positive phosphorus removal shown in Figure 2(e) arose from microbial assimilation and absorption (Rahimi et al. 2011). Moreover, Spearman correlations showed that the TIN/TN removal was positively correlated with that of IP/TP (R = 0.900, P = 0.019), implying that phosphorus was mainly removed by denitrifying PAOs.
Response of SND rate to DO variation
Taking the optimal SND rate as the expected value, the optimal DO concentration, which would meet the assumption of complete dissolved nitrogen removal, could be calculated. Using the optimal SND rate in Equation (3) gave DO = 2.10 mg/L, which suggested that increasing DO from 1.83 to 2.10 mg/L in R2 could improve TIN and/or TN removal efficiency. This is consistent with the trend shown in Figure 2(d).
RDA analysis between performance and effluent characteristics
For the on-effluent line parameters (Figure 4(a)), the changes in EC, TDS, and salinity were consistent, and all three were related negatively to the proportional removal of TIN, indicating that NOx-N removed by denitrification was an important contributor to them. The positive correlation between effluent pH and NOx-N removal was due to the alkalinity produced by denitrification (Wang et al. 2020), while the negative correlation between effluent pH and ammonia removal was due to the alkalinity consumed by nitrification.
The negative correlation between effluent pH and organics removal was attributed to the fact that the sodium acetate is alkaline (Li et al. 2020b), its degradation leading to a pH decrease. Meanwhile, the negative correlation between pH and TN/TIN removal might be explained by the large proportion of TAN in the influent, as nitrification, which consumes alkalinity, was the main factor affecting TN/TIN removal.
DO, T and ORP correlate negatively with the proportional removal of nitrate, but positively with those of TAN, COD, and TN. These correlations suggest that higher DO concentrations could provide more electron acceptors for heterotrophs and nitrifiers to oxidize organics and TAN, respectively (Machat et al. 2019; Jiang et al. 2020), and higher temperature could intensify their activity (Ma et al. 2020). These conditions would produce more nitrate, facilitating and stronger denitrification processes in the reactors.
For the effluent nutrient indices (Figure 4(b)), proportional TAN removal was positively related to nitrate concentration, because TAN was mainly oxidized to nitrate. TAN removal was also negatively related to effluent TCOD and DCOD because the excess effluent organic content inhibited nitrifier activity (Friedman et al. 2018). The obvious positive correlation between TN/TIN removal and nitrite suggested that the ammonium removal path in the reactors might include short-cut nitrification and denitrification (Jiang et al. 2020). Proportional TN/TIN removal was negatively correlated with effluent TCOD, DCOD, and TAN, indicating low nitrogen removal efficiency when oxidation in the reactor was insufficient. They also implied both insufficient ammonia oxidation capacity (more ammonia remained in the effluent), and insufficient carbon source utilization by denitrification due to the lack of electronic acceptors (not enough nitrate generated). Finally, the proportional removal of TP and IP was significantly negatively correlated with effluent nitrate, TN, and TIN, indicating that denitrifying phosphorus removal was the main path for simultaneous nitrogen and phosphorus removal in the reactors (Wang et al. 2015), which was consistent with the previous discussion (section 3.2).
CONCLUSIONS
Dissolved nitrogen was largely removed by SND from the aquaculture wastewater. The SND rate's response to different DO concentrations was described by a power function, through which the optimal SND rate or DO concentration could be calculated. Although dissolved nitrogen might be removed completely to meet the stringent environmental quality standard for surface water (GB3838-2002), the reactor displayed a low phosphorus removal performance. The phosphorus was removed mainly by denitrifying PAOs because of the lack of alternating anaerobic/aerobic conditions. Phosphorus could be removed to a greater extent by other processes, however, such as chemical coagulation/precipitation, to meet the environmental standard.
ACKNOWLEDGEMENTS
This work was financially supported by the Key Laboratory of Recreational Fisheries, Ministry of Agriculture and Rural Areas (ZJK201905), the Open Research Program from the Agriculture Ministry Key Laboratory of Healthy Freshwater Aquaculture (ZJK201902), and the Fundamental Research Funds for the Central Universities (WUT: 2019III107CG).
DECLARATION OF INTEREST STATEMENT
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
CREDIT AUTHORSHIP CONTRIBUTION STATEMENT
Shiyang Zhang: Investigation, Data curation & Writing. Jing Chen: Data curation & Analysis. Julin Yuan: Conceptualization, Funding acquisition. Guangjun Wang: Conceptualization, Funding acquisition.
DATA AVAILABILITY STATEMENT
All relevant data are available from the corresponding author.