Next Article in Journal
The Infancia y Procesamiento Sensorial (InProS—Childhood and Sensory Processing) Project: Study Protocol for a Cross-Sectional Analysis of Parental and Children’s Sociodemographic and Lifestyle Features and Children’s Sensory Processing
Previous Article in Journal
Infection Spread and High-Resolution Detection of Close Contact Behaviors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation and Analysis of Water Quality of Marine Aquaculture Area

1
Dalian Ocean University, College of Marine Sciences, Dalian 116023, China
2
Key Laboratory of Environment Controlled Aquaculture, Ministry of Education, Dalian Ocean University, Dalian 116023, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(4), 1446; https://doi.org/10.3390/ijerph17041446
Submission received: 5 January 2020 / Revised: 19 February 2020 / Accepted: 20 February 2020 / Published: 24 February 2020

Abstract

:
In the rapid development of marine aquaculture, the water quality of aquatic environments is regarded as a main limiting factor. Therefore, it is necessary to assess the water quality and environmental conditions in marine aquaculture areas and find out the main influencing factors regarding damage to the water quality environment. In the present research, pond aquaculture and cage aquaculture areas were sampled in May, August and November in 2018. Nine water quality indicators were detected, including pH, temperature, salinity, dissolved oxygen, molybdate-reactive phosphorus, chemical oxygen demand, chlorophyll a, inorganic nitrogen and antibiotic resistance genes (ARGs). Principal component analysis (PCA) was used to analyze the water quality conditions, spatial–temporal changes, and the driving factors in pond and cage aquaculture areas. The results showed that three main components were extracted from the pond aquaculture area, which explained 66.82% of the results, the most relevant factors are salinity, dissolved oxygen and ARGs. For the cage aquaculture area, three main components were extracted which can account for 72.99% of the results, the most relevant factors are chlorophyll a, salinity and dissolved oxygen. The comprehensive scores of the principal components indicated that the heaviest polluted months in pond and aquaculture areas were August and November, respectively. The water quality of the pond aquaculture area is mainly limited by the volume of the pond, while aquaculture activities and seasonality are the main factors for cage aquaculture. ARGs in cage culture areas showed more variety and frequency compared with pond culture areas, which indicated that terrestrial input might be one of the sources for ARGs occurrence. The results would be helpful for the relevant authorities to select water quality monitoring parameters in marine aquaculture areas.

1. Introduction

In recent years, with the rapid development of marine aquaculture in China, the increasing scale and breeding density of aquaculture have induced a frequent occurrence of diseases and environmental pollution [1]. A large amount of the nitrogen and phosphorus in feeds cannot be utilized by fish during aquaculture, which makes the nutrient concentration higher in aquaculture areas than in other natural areas [2]. The nutrient input from aquaculture effluent may cause a deterioration of water quality in the aquatic environment [3]. Besides, aquatic product quality is easily influenced by the condition of the water environment [4]. Some antibiotics are considered as “pseudo-persistent” contaminants due to their continuous release and lower degradability in the environment. The persistent organic pollutants in aquatic environment can be enhanced through biological amplification, which may pose a potential hazard to human health. Moreover, the abuse of antibiotics in aquaculture leads to the emergence of drug resistance in aquatic animals, which arouses serious clinical and public health problems on a global level [5,6].
China is the largest producer and user of antibiotics in the world; based on a survey in 2013 [7], veterinary antibiotics consumed 84,240 tons (52%) of the total antibiotics’ usage in China. Antibiotic residues have been detected at a high level in various environmental matrices. In the research of Zhou et al. [8], it was found that 25 antibiotics were detected in the Yangtze River, with a total concentration of 0.22–366 ng/L. A variety of antibiotics were widely detected in Dalian coastal waters, with concentrations ranging from 22.6 to 2402.4 ng/L [9]. The residual antibiotics in the environment could not only induce adverse effects on the nontarget organisms but also are considered as one of the selected stresses for the occurrence of the promotion of bacterial antibiotic resistance. It has been demonstrated that the frequent use of antibiotics in the animal husbandry industry as well as in aquaculture is one of the critical reasons for the high abundance of antibiotic resistance genes (ARGs) [10]. Usually, antibiotics are applied through feed addition in aquaculture; only about 25% antibiotics will be absorbed by aquatic animals, while the excess unabsorbed antibiotics or antibiotics not ingested by fishes will be discharged in the form of urine and feces, then directly dissolved/degraded in water by bacterial biomass, which eventually will come back to aquatic environment. The long-term existence of antibiotics leads to the accumulation of antibiotic-resistant bacteria and ARGs, which may pose a risk of transmission to the general population and to wildlife, causing treatment-resistant illness [11,12].
Zhuanghe, located in the southeast of Liaodong Peninsula, is an important marine aquaculture area in Dalian in China. It has a coast length of 285 km and culture area of 620 km2, with over than 190 aquaculture companies around. The main aquaculture species included sea cucumber (Stichopus japonicus), clam (Meretrix meretrix L), scallop (Patinopecten yessoensis), oyseter (ostrea gigas thunberg), shrimp (Penaeus orientalis), puffer (Takifugu rubripes), etc. Pond culture, cage culture, floating raft culture and indoor industrial culture are the main culturing models. Among them, pond culture is mainly carried out in enclosed ponds constructed with dams by the sea, and the exchange of water relies on the regularity of the tide. A positive aspect is that pond culture is convenient to manage with a relatively small water volume. On the other hand, water exchange is restricted by the water quality from the sea and is vulnerable to being affected by the water environmental conditions, such as organic pollutants, nutrients, temperature, dissolved oxygen and so on [13]. Cage and floating raft cultures are characterized in high-breeding density and high-profit areas and use metal or synthetic fibers to construct the cage or raft; however, these methods are susceptible to waves, temperature and other factors [14]. In recent years, marine water quality has deteriorated due to continuously discharging wastewater from industry and land domestic sewage, as well as untreated effluent from aquaculture. As a consequence, the quality and the yield of aquatic products significantly declined, while the source of contaminants was difficult to determine, which affected the implementation of pollution control. Therefore, it is of great significance to monitor the water quality and make a comprehensive evaluation in aquaculture areas [15].
Currently, the commonly used methods in marine environmental assessment for water quality classification include single factor assessment and multivariate analysis. Single factor assessment is a simple and convenient method for water quality classification, judged according to the most impaired assessment factor and compared to the corresponding water quality standard [16]. The multivariate analysis includes principal component analysis (PCA), factor analysis (FA), the fuzzy comprehensive evaluation method and integrated analysis method, etc. [17]. Compared with single factor assessment, PCA has the advantages of keeping the original information and uncorrelated variables at the same time, with the aim of screening out the independent principal component factors through dimension reduction analysis [18,19]. To effectively address the water quality impairment and pollutant sources, PCA has been extensively used in the evaluation of water quality, eutrophication and pollution degree [20]. At present, there are many studies on the analysis of water quality by PCA [21,22]. In the study by Ouyang, PCA was verified as a useful tool for the identification of important surface water quality monitoring stations and water quality parameters [23]. Besides, PCA was also used in the selection of water quality parameters for the water quality index (WQI), and the number of monitoring parameters was reduced from 28 to 9 after the application of PCA [24]. To ensure the objectivity of results and the effectivity of resources, the local conditions should be considered to reevaluate the PCA results. At the same time, in order to improve the method of water quality assessment, the role of principal component analysis still needs to be compared with other water quality assessment analysis methods that can reduce the workload by reducing the indicators [25].
Although emerging pollutants in the aquatic environment have been frequently detected, they were not included in the scope of water quality monitoring. The types of pollutants varied among different aquaculture models, which affected the distribution and transfer of ARGs in the environment. At present, there is no related research on water quality assessment with ARGs and other water quality parameters. In the present study, principal component analysis was used to assess the temporal and spatial variations of water quality in pond and cage culture areas and to analyze the driving factors among various water quality parameters and the ARGs detection rate in the different aquaculture areas. The results are expected to provide a valuable method for the evaluation of impaired water environments, especially for areas with fixed pollutant input.

2. Materials and Methods

2.1. Sampling Sites and Time

The sampling sites covered both pond culture and cage culture areas in a typical district of marine aquaculture in Zhuanghe. There are nine culture ponds (1–9) featuring the breeding of sea cucumber on the cofferdam shore, and 10 points (10–17) around cage or floating raft with fish and shellfish cultures distributed 15–30 km away from the mainland. The water outside the ponds and cages came from the same source, and water quality differences and effects between them were ignored in the present study. Detailed sampling site information is shown in Figure 1. Considering the growth cycle of cultured species and the fishery activities, samplings was conducted in May, August and November in 2018, which represented the growth, breeding and harvest period of aquaculture organisms, respectively. For sea cucumber culture ponds, water was exchanged depending on the tide once or twice per month. Sea cucumber seedlings were cultured in sampling sites 1–3, while adult sea cucumber was cultured in sampling sites 4–9. For cage culture areas, the bivalves Argopecten irradias, Chlamys farreri, Crassostrea gigas and Scapharca broughtonii were the main breeding species; thereinto fish net cages for Seriola aureovittata breeding were settled near sampling sites 14, 16 and 17. The water depths of the sea cucumber cultural ponds were about 2.0 m, and those of the cage culture area were in the range of 8.6–24.8 m.

2.2. Sampling and Chemical Analysis

The sampling, preservation, transportation and detection of water samples were strictly performed following the National Specification for Marine Monitoring standard methods (GB 17378.4-207) (SEPA, 2002). Water samples were taken at depth of 0.5 and 3.0 m in the pond and cage culture areas, respectively. The parameters of pH (± 0.2), temperature (T, ± 0.02 °C), salinity (SAL, ± 0.1 ppt), dissolved oxygen (DO, ± 0.01 mg/L) were detected by a multi-parameter water quality analyzer (ProPlus, YSI, Yellow Springs, OH, USA) at the sampling sites. Molybdate-reactive phosphorus (MRP) was determined by the methods of phosphor molybdenum blue spectrophotometry. Under acid conditions, molybdic acid reacts with inorganic phosphorus to form yellow ammonium phosphomolybdate, which is then reduced to dark blue molybdenum blue by the original agent. The absorbance was measured at 650 nm by spectrophotometer, and the content of inorganic phosphorus was calculated by a standard curve. Chemical oxygen demand (COD) was determined by the method using alkaline potassium permanganate, which is not affected by salinity [26]. Under the condition of alkaline heating, we oxidized the oxygen demand substance in seawater with a known amount of potassium permanganate. Under sulfuric acid conditions, excess potassium permanganate and manganese dioxide were reduced by potassium iodide. The free iodine was titrated with sodium thiosulfate standard solution. The volume of sodium thiosulfate consumed was determined by a formula and known concentration. Chlorophyll a (Chla) content was determined by Dimethyl Formamide (DMF) extraction methods. The chlorophyll was extracted by DMF, which could be directly determined according to the absorbance after being kept away from light for 12 h [27]. Dissolved inorganic nitrogen (DIN, the sum of ammonium nitrogen (NH4+-N), nitrate nitrogen (NO3--N) and nitrite nitrogen (NO2--N)) were determined by an automatic continuous flow analyzer (AA3, SEAL) [28]. In order to detect the emergence of ARGs in aquaculture wastewater, 1500 mL water samples were filtered by a 0.22 μm filter membrane; then, DNA was extracted by DNA extraction kit (Chengdu Foregene Biotechnology Co., Ltd, Chengdu, China). Seven ARGs (flor, sul1, sul2, tetB, tetM, qnrS and ermB) from five classes of commonly used antibiotics in aquacultures were selected. A Veriti™ 96-Well Thermal Cycler (Applied Biosystems™, ThermoFisher Scientific, Waltham, MA, USA) was employed using 2 × EasyTaq® PCR SuperMix (Beijing Transgen Biotech, Beijing, China) to determine the expression of ARGs. The primer sequences of ARGs are detailed in Table 1. According to the results of questionnaires, no records of antibiotic application in the process of aquaculture were found.

2.3. Data Processing and Principal Component Analysis

The statistical analysis was based on the standardized sample matrix (19 sites × 3 phases × 9 variables) and summarized in Table 2, which presented the mean values ± SD of the results. To avoid misclassification due to wide differences in data dimensionality, data was standardized through z-scale transformation before applying PCA [36]. For the detection rate of AGRs, there is no exact standard range of AGRs that can be used as an index to analyze the other water quality factors. In the present study, when ARGs were detectable, the result was regarded as 1; otherwise, it was regarded as 0. The data were normalized according to the following formula:
Z i j = X i j μ σ
where Z i j represents the normalized value of the i-th station’s j index; X i j represents the measured value of the i-th station’s j indicator; μ represents the mean value of variables; and σ represents the standard deviation of variables.
In this study, all the data statistical analysis was conducted by SPSS 18.0, KMO and Bartlett’s sphericity test were used to verify the applicability of PCA to raw data [37]. Then, the corresponding correlation coefficient matrix, eigenvalues and eigenvectors were calculated. We adopted the principle that the eigenvalue > 1.0 to determine the number of principal components m [21]. The maximum variance was used to normalize the rotation of the original data to rotate the factor [38].
To estimate the degree of contamination, the common factor for each sample needed to be calculated, and the principal component composite score could be used both for model diagnosis and as raw data for further analysis. Firstly, the values of the first m principal components of each sample were calculated as the following formula:
F i g = Z i 1 l g 1 + Z i 2 l i g + Z i g l i g + + Z i g l i g i = 1 , 2 , , n ; g = 1 , 2 , , m .
Then, the weighted summation was performed according to the first m principal components, and the weight coefficient of each principal component was the variance contribution rate of the principal component, i.e., λ g / g = 1 m λ g . For the sample i, the comprehensive score of each month’s standing was calculated using the following formula:
F i = g = 1 m λ g / g = 1 m λ g .
The weighted sum of the previous three principal components was used to get the comprehensive scores of each sample, by which the Box-plot and the sampling sites’ comprehensive score maps were made. Spearman’s rank correlation coefficient analysis was used to analyze the relevance between factors and the principal components, and we obtained the main driving factors [39].

3. Results

3.1. Water Quality Parameter Detection

As shown in Table 2, nine water quality parameters were summarized as mean value ± standard deviation. During the entire detected period, no abnormal raw data for all the parameters were found. The values of pH varied in the range of 7.9–8.4, and no significant differences were observed between seasons and areas. The values of temperature varied in the range of 9.8–24.7 ℃, and the highest water temperature of the pond culture area was shown in May and in August for the cage culture area. The values of SAL in the two aquaculture areas varied in the range of 28.5–32.1 ppt, while the lowest and highest values of SAL were shown in the pond culture area and kept relatively stable in cage culture area. Dissolved oxygen (DO) values varied greatly with seasons, ranging from 7.3 to 10.4 mg/L (DO saturation index 77.18%–99.30%). The contents of Chla varied in the range of 0.43–6.86 μg/L. The amount of COD, MRP and DIN varied in range of 0.56–11.12 mg/L, 0.008–0.328 mg P/L and 0.002–1.796 mg N/L, respectively. Compared with the water quality standard, the dissolved oxygen (<6 mg/L), molybdate-reactive phosphorus (<0.045 mg P/L) and inorganic nitrogen (<0.2 mg N/L) in the present study meet the class I of water quality standard, while contents of COD (> 5 mg/L) exceeding the standard were observed in August in the pond culture area. In contrast, DO, MRP and COD in the cage culture area did not exceed the standard for the entire monitored month.
To further explore the difference of ARGs distribution in different sources, the detection rates of seven ARGs in pond and cage culture areas were analyzed. The detection rates of ARGs varied greatly ranging from 11%–100%, and relatively low detection rates of ARGs were observed in the pond culture area. Six kinds of ARGs were detected in the pond aquaculture area; among them, the detection rates were 11% and 94% for sulfonamides ARGs (sul1 and sul2), 5% for chloramphenicoles ARGs (flor), 6% and 0% for tetracyclines ARGs (tetM and tetB), 11% for quinolones ARGs (qnrS) and 44% for macrolides ARGs (ermB). In the cage culture area, the detection rates of sulfonamides ARGs (sul1 and sul2) and chloramphenicoles ARGs (flor) were 100%, those for tetracyclines ARGs (tetM and tetB) were 81% and 57%, quinolones ARGs (qnrS) and those for macrolides ARGs (ermB) were 87% and 82%, respectively.

3.2. Water Quality Assessment by PCA

Principal component analysis was performed on the data matrix of variables from pond and cage marine aquaculture areas in Zhuanghe using the computed eigenvalues and weights of parameters. After the verification of the data’s validity by Bartlett’s sphericity test (<0.001) and KMO test for pond and cage culture, the nine parameters yielded three principal components (eigenvalues > 1) explaining sample variances of about 66.82% and 72.99% for pond and cage culture areas, respectively (Table 3). Since pond and cage culture areas were combined to calculate the correlation matrix, the correlation coefficients should be interpreted with caution as they are simultaneously affected both by spatial and temporal variations. Various factor loadings in principal components are shown in Table 4. For the pond culture area, PC1 explained 31.79% of total variance, which had positive loadings (> 0.70) on COD (0.72) and Chla (0.74) and significant negative loadings on SAL (–0.86). In the cage culture area, pH (0.71), T (–0.76) and Chla (0.80) occupied higher loading in PC1 (> 0.70).

3.3. Spatial and Temporal Variations of Water Quality

Box and whisker plots of the comprehensive score of principal components were constructed to evaluate water quality in different aquaculture patterns associated with temporal variations (Figure 2). The higher the comprehensive scores of principal components, the more serious the water pollution. In the pond culture area, the highest comprehensive score of principal components was in August (0.88), then in May (–0.46), and lowest in November (–0.53). In the cage culture area, the comprehensive score was highest in November (1.19), second in May (–0.08), and lowest in August (–1.00).
The spatial distributions of comprehensive scores in different months in the pond culture area and cage culture area of Zhuanghe are shown in Figure 3. Obvious spatial distribution trends in different culture stages were observed. In May and August, comprehensive scores in most sampling sites of the pond culture area were significantly higher than that in the cage culture area. In contrast, comprehensive scores in the pond culture area were significantly lower than in the cage culture area in November.
Principal component analysis is an effective tool to figure out the driving factors of water pollution. In order to extract the dominant factors from numerous physicochemical parameters related to water quality, the driving factors of water pollution in different marine aquaculture areas were analyzed. In the present study, a Spearman’s rank correlation coefficient analysis between the principal components and parameters was carried out followed with the research of Lundberg et al. [39]. As shown in Table 5, SAL, DO and ARGs were the main driving factors of PC1, PC2 and PC3 (P < 0.01) in the pond culture area, while Chla, SAL and DO were the main driving factors of PC1, PC2 and PC3 (p < 0.01) in cage culture areas.

4. Discussion

PCA is useful when several related random environmental variables are simultaneously considered. The nine original variables were reduced to three key independent factors that influenced water quality in pond and cage culture areas. In the pond culture area, PC1 had high loadings on SAL, Chla, and COD, which represents the contribution of phytoplankton growth from the natural impacts on the water condition, such as the water exchange rate of pond and rainfall. In early spring, when the ice begins to melt, fluctuations in salinity affect cultured organisms. Even though algae have definite selectivity and tolerability to salinity variability [40], their community structure is still easily influenced by salinity fluctuations [41]. COD served as an important parameter in traditional water quality monitoring to reflect the organic matter in marine waters. To guarantee the rapid growth of sea cucumber, herbicide may be used to suppress the outbreak of macroalgae in sea cucumber culture ponds, and the death of macroalgae may cause the increase of COD [42]. PC2 had high loadings on temperature and DO, which can be regarded as the index of aquaculture organism’s growth conditions in pond culture. Temperature and dissolved oxygen are the critical factors in the process of aquaculture [43], which presented significant influences regarding the growth and development of aquaculture organisms. PC3 had high loadings on DIN and ARGs, which was identified as human activity-induced pollution. Nitrate in ponds is mainly derived from aquaculture activities and wastewater discharges [44], and ARGs occurrences mainly came from microorganisms after antibiotic administration in ponds [45,46,47]. In the cage culture area, PC1 had high loadings on Chla, T, pH which can be regarded as the photosynthesis of phytoplankton. Temperature affects phytoplankton growth by affecting photosynthesis and phosphorus concentration [48], and pH affects the photosynthesis of phytoplankton by affecting fluorescence absorption of chlorophyll [49].In PC2, SAL with positive loading and DIN with negative loading occupied relative high loadings. Salinity variations led to the changes of microbial function and the process of nitrification and denitrification, affecting the contents of DIN in water [50]. Therefore, PC2 can be regarded as the microorganism changes in seawater. In PC3, DO showed a high loading.
Different trends in the comprehensive score box of the pond culture and cage culture were observed, and the season-correlated parameter can be taken as representing the major source of temporal variations in water quality. In the summer of 2018 in Dalian, the abnormally high temperature caused the death of many aquaculture organisms, which may explain the highest comprehensive score of principal components in August in the pond culture area. The volume of water in the pond culture area was limited and easily influenced under high temperature, which would aggravate the pollution and result in a high comprehensive score in August. In the cage culture area, the comprehensive score was highest in November and lowest in August. As regards the better water exchange ability in the cage culture area, high temperature was not the critical factor in the impairment of water quality in August. Moreover, the light intensity and illumination period was suitable for the growth of phytoplankton, which accordingly induced the decrease of nutrient concentration in seawater. In November, the water temperature decreases, which is regarded as one of the factors leading to the decline of primary productivity [51]. Fishery harvest activities and the frequent human activities were also considered as critical reasons for high comprehensive scores and heavy pollution during this period in the cage culture area.
The spatial analysis revealed a significant difference between pond and cage culture areas. April and May are the important cultivation period for sea cucumber seedlings. To achieve the rapid growth rate, soybean and eelgrass (Zostera marina) were fed as the supplemental feeds, the remnants of which may cause a deterioration of water quality and elevate the comprehensive scores of principal components. In cage culture, sampling sites 16 and 17 exhibited high-level contaminants, in which Seriola aureovittata were cultured in fish net cages intensively. The extreme weather played a direct role in the ponds with shallow water depth, which caused the death and hypoxia of aquaculture organisms, as well as the outbreak of epidemic disease. In November, the opposite phenomenon for the scores of the cage culture area was observed, which was significantly higher than that of the pond culture area. This may be due to the pollution degree of the seawater in cage culture areas, which was increased by the influence of human fishery harvest activities and fewer sea cucumber activities in pond culture under low temperature conditions. Considering the descriptions above, the difference in aquaculture water quality was mainly affected by the species and the aquaculture patterns.
The common driving factors in pond and cage culture areas are SAL and DO. Salinity is a key factor in the process of aquaculture and showed a negative correlation in the pond culture area and positive correlation in the cage culture area. From the perspective of original data, the salinity varied within the normal range (28.2–33.0), with the variations probably due to the sudden weather changes [52]. It was found that, in August, Dalian was often affected by heavy rainfall. A cultural pond is a half-open artificial ecological system, which is susceptible to the disturbance of abnormal natural conditions. Heavy rainfall brought a huge amount of freshwater injection and disturbed the sediment in the pond, and the salinity sharply decreased which gave fierce fluctuations in bacteria communities and oxygen concentrations [53], which may also have resulted in the death of sea cucumbers [54]. Thus, this showed the serious pollution and negative correlation in the area with low salinity.
For the cage culture area, rainfall could bring oxygen and other substances and promote the growth of plankton, which is conducive to the absorption of nutrients. Besides, DO is also a crucial indicator which not only directly affects the growth of cultured organisms but also is vital to the decomposition of harmful substances at the bottom of the pond. It takes 4570 g of oxygen to oxidize 1000 g of ammonia to nitrate; the nitrification reaction was blocked under the low DO concentration. The decrease of DO led to the increase of CO2 content and the decrease of the pH value. It has been reported that the oxygen consumption of oxidation decomposition accounts for 40% of the total oxygen consumption in the aquaculture pond, while only 12% was consumed for aquatic organisms, even though this was less acute in seawater than in freshwater [55]. In both cultural areas, DO presented a positive correlation, which indicated that the water quality was better in the area with a high concentration of dissolved oxygen. In aquaculture, DO is influenced by the temperature and density of cultural organisms; it dropped significantly when the temperature and breeding density increased. In addition, Chla showed a positive correlation in the cage culture area. Chla was considered as the key driving factor in the culture area which could represent the density of phytoplankton, which contributed to the utilization of the redundant nutrition. The abundance of Chla indicated the risk of eutrophication, while the weakened sunlight may have decreased the photosynthesis of algae and the production of oxygen. Furthermore, oxygen in the water was to a great extent consumed by the degradation of organic matter from dead phytoplankton cells, which may aggravate the water pollution. Furthermore, ARGs showed a negative correlation in the pond culture area. ARGs were considered as an emerging contaminant in the aquatic environment due to their potential risk. It was reported that 11 ARGs with frequent presence were observed in the sediments of the East China Sea bays, and sulfonamide resistance was the most prevalent ARG type [12]. Consistently, sulfonamide resistance was also detected at a high rate both in pond and cage culture areas. ARGs in the cage culture area showed more variety and frequency compared with the pond area, which indicated the terrestrial input might be one of the sources for ARGs occurrence, such as hospitals, land farms and wastewater treatment plants. ARGs can be widely transferred among bacterial species in the various environmental media. Even no record of antibiotics application in aquaculture, the residues of heavy metals, micro plastics and other organic pollutants in the aquatic environment may promote the transfer of ARGs [56]. Therefore, the risk assessment and removal of ARGs in aquaculture wastewater should be taken into consideration.

5. Conclusions

Water quality assessment with ARGs detection rate and other water quality parameters between different culture models was investigated in 19 sampling sites in the Zhuanghe marine aquaculture region. The PCA results suggest that three principal components explained 66.82% and 72.99% of the results in pond and cage culture areas. According to the time–space analysis of principal component comprehensive scores, August and November were the most polluted months in pond and cage culture areas, and seasonal variation and human fishery harvest activities affect the pollution of these two different types of aquaculture areas seriously. Further analysis showed that Chla, SAL, DO and ARGs are the main driving factors. Because of the limited area of the pond, weather changes exerted critical impacts on its water quality, while the water quality in the cage culture was mainly related to human fishery activities. In the present study, the detection rates of ARGs were varied in pond and cage culture areas. To explore the reason for the occurrence of ARGs in aquaculture, environmental stresses such as antibiotics, heavy metals, etc., should be included in future research.

Author Contributions

Conceptualization, Q.Z.; data curation, X.Z., P.L., J.L. and S.J.; funding acquisition, Y.L.; methodology, Y.Z., L.C. and Z.M.; project administration, Q.Z. and Y.L.; visualization, R.G.; writing—original draft, X.Z.; writing—review and editing, X.Z. and Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key R&D Program of Guangdong Province, grant number 2019B020215001; National Key R&D Program of China, grant number 2017YFB0404000; High-level Innovation Team Overseas Training Project of Liaoning Provincial, grant number 201818; PhD Research Funds of Liaoning Province, grant number 2019-BS-035;” and “The APC was funded by Key R&D Program of Guangdong Province”.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ren, C.; Wang, Z.; Zhang, Y.; Zhang, B.; Chen, L.; Xi, Y.; Xiao, X.; Doughty, R.B.; Liu, M.; Jia, M.; et al. Rapid expansion of coastal aquaculture ponds in China from Landsat observations during 1984–2016. Int. J. Appl. Earth Obs. Geoinf. 2019, 82, 101902. [Google Scholar] [CrossRef]
  2. Funge-Smith, S.J.; Briggs, M.R.P. Nutrient budgets in intensive shrimp ponds: Implications for sustainability. Aquaculture 1998, 164, 117–133. [Google Scholar] [CrossRef]
  3. Falconer, L.; Telfer, T.C.; Ross, L.G. Modelling seasonal nutrient inputs from non-point sources across large catchments of importance to aquaculture. Aquaculture 2018, 495, 682–692. [Google Scholar] [CrossRef] [Green Version]
  4. Saavedra, M.; Pereira, T.G.; Carvalho, L.M.; Pousão-Ferreira, P.; Grade, A.; Teixeira, B.; Quental-Ferreira, H.; Mendes, R.; Bandarra, N.; Gonçalves, A. Wild and farmed meagre, Argyrosomus regius: A nutritional, sensory and histological assessment of quality differences. J. Food Compos. Anal. 2017, 63, 8–14. [Google Scholar] [CrossRef]
  5. Muziasari, W.; Pärnänen, K.; Johnson, T.; Lyra, C.; Karkman, A.; Stedtfeld, R.; Tamminen, M.; Tiedje, J.; Virta, M. Aquaculture changes the profile of antibiotic resistance and mobile genetic element associated genes in Baltic Sea sediments. FEMS Microbiol. Ecol. 2016, 92, 52. [Google Scholar] [CrossRef] [Green Version]
  6. Vikesland, P.; Pruden, A.; Alvarez, P.; Aga, D.; Bürgmann, H.; Li, X.D.; Manaia, C.; Nambi, I.; Wigginton, K.; Zhang, T.; et al. Toward a Comprehensive Strategy to Mitigate Dissemination of Environmental Sources of Antibiotic Resistance. Environ. Sci. Technol. 2017, 51, 13061–13069. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, Q.Q.; Ying, G.G.; Pan, C.G.; Liu, Y.S.; Zhao, J.L. Comprehensive Evaluation of Antibiotics Emission and Fate in the River Basins of China: Source Analysis, Multimedia Modeling, and Linkage to Bacterial Resistance. Environ. Sci. Technol. 2015, 49, 6772–6782. [Google Scholar] [CrossRef]
  8. Zhou, L.-J.; Li, J.; Zhang, Y.; Kong, L.; Jin, M.; Yang, X.; Wu, Q.L. Trends in the occurrence and risk assessment of antibiotics in shallow lakes in the lower-middle reaches of the Yangtze River basin, China. Ecotoxicol. Environ. Saf. 2019, 183, 109511. [Google Scholar] [CrossRef]
  9. Du, J.; Zhao, H.; Wang, Y.; Xie, H.; Zhu, M.; Chen, J. Presence and environmental risk assessment of selected antibiotics in coastal water adjacent to mariculture areas in the Bohai Sea. Ecotoxicol. Environ. Saf. 2019, 177, 117–123. [Google Scholar] [CrossRef]
  10. Yuan, K.; Wang, X.; Chen, X.; Zhao, Z.; Fang, L.; Chen, B.; Jiang, J.; Luan, T.; Chen, B. Occurrence of antibiotic resistance genes in extracellular and intracellular DNA from sediments collected from two types of aquaculture farms. Chemosphere 2019, 234, 520–527. [Google Scholar] [CrossRef]
  11. McKinney, C.; Pruden, A. Ultraviolet Disinfection of Antibiotic Resistant Bacteria (ARBs) and their Antibiotic Resistance Genes (ARGs) in Water and Wastewater. Environ. Sci. Technol. 2012, 46, 13393–13400. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, J.; Su, Z.; Dai, T.; Huang, B.; Mu, Q.; Zhang, Y.; Wen, D. Occurrence and distribution of antibiotic resistance genes in the sediments of the East China Sea bays. J. Environ. Sci. 2019, 81, 156–167. [Google Scholar] [CrossRef] [PubMed]
  13. Fu, J.; Yokoyama, H.; Cui, B.; Zhou, J.; Yan, J.; Ma, X.; Shibata, S. Impacts of sea cucumber farming on biogeochemical characteristics in the Yellow River estuary, Northern China. Phys. Chem. Earth Parts A/B/C 2017, 97, 19–30. [Google Scholar] [CrossRef]
  14. Islam, M.S. Nitrogen and phosphorus budget in coastal and marine cage aquaculture and impacts of effluent loading on ecosystem: Review and analysis towards model development. Mar. Pollut. Bull. 2005, 50, 48–61. [Google Scholar] [CrossRef]
  15. Liao, Y.; Shou, L.; Jiang, Z.; Tang, Y.; Du, P.; Zeng, J.; Chen, Q.; Yan, X.; Chen, J. Effects of fish cage culture and suspended oyster culture on macrobenthic communities in Xiangshan Bay, a semi-enclosed subtropical bay in eastern China. Mar. Pollut. Bull. 2019, 142, 475–483. [Google Scholar] [CrossRef]
  16. Tripathi, M.; Singal, S.K. Allocation of weights using factor analysis for development of a novel water quality index. Ecotoxicol. Environ. Saf. 2019, 183, 109510. [Google Scholar] [CrossRef]
  17. Simeonov, V.; Stratis, J.A.; Samara, C.; Zachariadis, G.; Voutsa, D.; Anthemidis, A.; Sofoniou, M.; Kouimtzis, T. Assessment of the surface water quality in Northern Greece. Water Res. 2003, 37, 4119–4124. [Google Scholar] [CrossRef]
  18. Shine, J.; Ika, R.; Ford, T. Multivariate Statistical Examination of Spatial and Temporal Patterns of Heavy Metal Contamination in New Bedford Harbor Marine Sediments. Environ. Sci. Technol. 1995, 29, 1781–1788. [Google Scholar] [CrossRef]
  19. Loska, K.; Wiechuła, D. Application of principal component analysis for the estimation of source of heavy metal contamination in surface sediments from the Rybnik Reservoir. Chemosphere 2003, 51, 723–733. [Google Scholar] [CrossRef]
  20. Perkins, R.G.; Underwood, G.J.C. Gradients of chlorophyll a and water chemistry along an eutrophic reservoir with determination of the limiting nutrient by in situ nutrient addition. Water Res. 2000, 34, 713–724. [Google Scholar] [CrossRef]
  21. Mandal, U.K.; Warrington, D.N.; Bhardwaj, A.K.; Bar-Tal, A.; Kautsky, L.; Minz, D.; Levy, G.J. Evaluating impact of irrigation water quality on a calcareous clay soil using principal component analysis. Geoderma 2008, 144, 189–197. [Google Scholar] [CrossRef]
  22. Guo, X.; Zhang, X.; Yue, H. Evaluation of hierarchically weighted principal component analysis for water quality management at Jiaozuo mine. Int. Biodeterior. Biodegrad. 2018, 128, 182–185. [Google Scholar] [CrossRef]
  23. Ouyang, Y. Evaluation of river water quality monitoring stations by principal component analysis. Water Res. 2005, 39, 2621–2635. [Google Scholar] [CrossRef]
  24. Tripathi, M.; Singal, S.K. Use of Principal Component Analysis for parameter selection for development of a novel Water Quality Index: A case study of river Ganga India. Ecol. Indic. 2019, 96, 430–436. [Google Scholar] [CrossRef]
  25. Aldieri, L.; Grafstrom, J.; Sundstrom, K.; Vinci, C.P. Wind power and job creation. Sustainability 2019, 12, 45. [Google Scholar] [CrossRef] [Green Version]
  26. Ma, J. Determination of chemical oxygen demand in aqueous samples with non-electrochemical methods. Trends Environ. Anal. Chem. 2017, 14, 37–43. [Google Scholar] [CrossRef]
  27. Cubas, C.; Gloria Lobo, M.; González, M. Optimization of the extraction of chlorophylls in green beans (Phaseolus vulgaris L.) by N,N-dimethylformamide using response surface methodology. J. Food Compos. Anal. 2008, 21, 125–133. [Google Scholar] [CrossRef]
  28. Feng, J.; Li, F.; Zhou, X.; Xu, C.; Fang, F. Nutrient removal ability and economical benefit of a rice-fish co-culture system in aquaculture pond. Ecol. Eng. 2016, 94, 315–319. [Google Scholar] [CrossRef]
  29. Van, T.T.H.; Chin, J.; Chapman, T.; Tran, L.T.; Coloe, P.J. Safety of raw meat and shellfish in Vietnam: An analysis of Escherichia coli isolations for antibiotic resistance and virulence genes. Int. J. Food Microbiol. 2008, 124, 217–223. [Google Scholar] [CrossRef]
  30. Lu, Z.; Na, G.; Gao, H.; Wang, L.; Bao, C.; Yao, Z. Fate of sulfonamide resistance genes in estuary environment and effect of anthropogenic activities. Sci. Total Environ. 2015, 527–528, 429–438. [Google Scholar] [CrossRef]
  31. Looft, T.; Johnson, T.; Allen, H.; Bayles, D.; Alt, D.; Stedtfeld, R.; Sul, W.J.; Stedtfeld, T.; Chai, B.; Cole’s, R.J.; et al. In-feed antibiotic effects on the swine intestinal microbiome. Proc. Natl. Acad. Sci. USA 2012, 109, 1691–1696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Luo, Y.; Mao, D.; Rysz, M.; Zhou, Q.; Zhang, H.; Xu, L.; Alvarez, P. Trends in Antibiotic Resistance Genes Occurrence in the Haihe River, China. Environ. Sci. Technol. 2010, 44, 7220–7225. [Google Scholar] [CrossRef] [PubMed]
  33. Guo, M.-T.; Yuan, Q.-B.; Yang, J. Ultraviolet reduction of erythromycin and tetracycline resistant heterotrophic bacteria and their resistance genes in municipal wastewater. Chemosphere 2013, 93, 2864–2868. [Google Scholar] [CrossRef] [PubMed]
  34. Strahilevitz, J.; Jacoby, G.; Hooper, D.; Robicsek, A. Plasmid-Mediated Quinolone Resistance: A Multifaceted Threat. Clin. Microbiol. Rev. 2009, 22, 664–689. [Google Scholar] [CrossRef] [Green Version]
  35. Chen, J.; Yu, Z.; Michel, F.; Wittum, T.; Morrison, M. Development and Application of Real-Time PCR Assays for Quantification of erm Genes Conferring Resistance to Macrolides-Lincosamides-Streptogramin B in Livestock Manure and Manure Management Systems. Appl. Environ. Microbiol. 2007, 73, 4407–4416. [Google Scholar] [CrossRef] [Green Version]
  36. Finch, B.K.; Beck, A.N. Socio-economic status and z-score standardized height-for-age of U.S.-born children (ages 2–6). Econ. Hum. Biol. 2011, 9, 272–276. [Google Scholar] [CrossRef] [Green Version]
  37. Fávero, L.P.; Belfiore, P. Chapter 12—Principal Component Factor Analysis. In Data Science for Business and Decision Making; Fávero, L.P., Belfiore, P., Eds.; Academic Press: Cambridge, MA, USA, 2019; pp. 383–438. [Google Scholar]
  38. Kaiser, H.F. The varimax criterion for analytic rotation in factor analysis. Psychometrika 1958, 23, 187–200. [Google Scholar] [CrossRef]
  39. Lundberg, C.; Lönnroth, M.; von Numers, M.; Bonsdorff, E. A multivariate assessment of coastal eutrophication. Examples from the Gulf of Finland, northern Baltic Sea. Mar. Pollut. Bull. 2005, 50, 1185–1196. [Google Scholar] [CrossRef]
  40. Singh, R.; Upadhyay, A.K.; Chandra, P.; Singh, D.P. Sodium chloride incites reactive oxygen species in green algae Chlorococcum humicola and Chlorella vulgaris: Implication on lipid synthesis, mineral nutrients and antioxidant system. Bioresour. Technol. 2018, 270, 489–497. [Google Scholar] [CrossRef]
  41. Ryan, K.G.; Tay, M.L.; Martin, A.; McMinn, A.; Davy, S.K. Chlorophyll fluorescence imaging analysis of the responses of Antarctic bottom-ice algae to light and salinity during melting. J. Exp. Mar. Biol. Ecol. 2011, 399, 156–161. [Google Scholar] [CrossRef]
  42. Changzi, G.; Kaijing, Z.; Chunfeng, Q.; College, M.; University, S. Safe Utilization of Prometryne in Sea Cucumber Culture Pond: From Nutrient Diffusion Flux Across Biodeposition-water Interface Aspect. Chin. Agric. Sci. Bull. 2019, 35, 149–153. [Google Scholar]
  43. Qiang, J.; Zhong, C.Y.; Bao, J.W.; Liang, M.; Liang, C.; Li, H.X.; He, J.; Xu, P. The effects of temperature and dissolved oxygen on the growth, survival and oxidative capacity of newly hatched hybrid yellow catfish larvae (Tachysurus fulvidraco♀ × Pseudobagrus vachellii). J. Therm. Biol. 2019, 86, 102436. [Google Scholar] [CrossRef] [PubMed]
  44. Christianson, L.; Lepine, C.; Tsukuda, S.; Saito, K.; Summerfelt, S. Nitrate removal effectiveness of fluidized sulfur-based autotrophic denitrification biofilters for recirculating aquaculture systems. Aquac. Eng. 2015, 68, 10–18. [Google Scholar] [CrossRef] [Green Version]
  45. Páez-Osuna, F.; Guerrero-Galván, S.R.; Ruiz-Fernández, A.C. The environmental impact of shrimp aquaculture and the coastal pollution in Mexico. Mar. Pollut. Bull. 1998, 36, 65–75. [Google Scholar] [CrossRef]
  46. Shen, X.; Jin, G.; Zhao, Y.; Shao, X. Prevalence and distribution analysis of antibiotic resistance genes in a large-scale aquaculture environment. Sci. Total Environ. 2020, 711, 134626. [Google Scholar] [CrossRef] [PubMed]
  47. Luo, Y.; Zhou, Q. Antibiotic resistance genes (ARGs) as emerging pollutants. Huanjing Kexue Xuebao Acta Sci. Circumstantiae 2008, 28, 1499–1505. [Google Scholar] [CrossRef]
  48. Chen, M.; Li, J.; Dai, X.; Sun, Y.; Chen, F. Effect of phosphorus and temperature on chlorophyll a contents and cell sizes of Scenedesmus obliquus and Microcystis aeruginosa. Limnology 2011, 12, 187–192. [Google Scholar] [CrossRef]
  49. Qu, F.; Gong, N.; Wang, S.; Gao, Y.; Sun, C.; Fang, W.; Men, Z. Effect of pH on fluorescence and absorption of aggregates of chlorophyll a and carotenoids. Dye. Pigment. 2020, 173, 107975. [Google Scholar] [CrossRef]
  50. Xia, Z.; Wang, Q.; She, Z.; Gao, M.; Zhao, Y.; Guo, L.; Jin, C. Nitrogen removal pathway and dynamics of microbial community with the increase of salinity in simultaneous nitrification and denitrification process. Sci. Total Environ. 2019, 697, 134047. [Google Scholar] [CrossRef]
  51. Zhou, B.; Wang, Y.; Wang, J.H.; Liu, Y.; Tang, X.X. A multivariate assessment of eutrophication in the coastal waters of the southern Shandong Peninsula. Acta Oceanol. Sin. 2010, 22, 2186–2190. [Google Scholar] [CrossRef]
  52. Bhurtun, P.; Lesven, L.; Ruckebusch, C.; Halkett, C.; Cornard, J.-P.; Billon, G. Understanding the impact of the changes in weather conditions on surface water quality. Sci. Total Environ. 2019, 652, 289–299. [Google Scholar] [CrossRef] [PubMed]
  53. Huang, L.; Xu, Y.; Xu, J.; Ling, J.; Zheng, L.; Zhou, X.; Xie, G. Dissemination of antibiotic resistance genes (ARGs) by rainfall on a cyclic economic breeding livestock farm. Int. Biodeterior. Biodegrad. 2019, 138, 114–121. [Google Scholar] [CrossRef]
  54. Yuan, X.; Yang, H.; Wang, L.; Zhou, Y.; Gabr, H.R. Effects of salinity on energy budget in pond-cultured sea cucumber Apostichopus japonicus (Selenka) (Echinodermata: Holothuroidea). Aquaculture 2010, 306, 348–351. [Google Scholar] [CrossRef]
  55. Li, F.; Sun, Z.; Qi, H.; Zhou, X.; Xu, C.; Wu, D.; Fang, F.; Feng, J.; Zhang, N. Effects of Rice-Fish Co-culture on Oxygen Consumption in Intensive Aquaculture Pond. Rice Sci. 2019, 26, 50–59. [Google Scholar] [CrossRef]
  56. Lu, J.; Zhang, Y.; Wu, J.; Luo, Y. Effects of microplastics on distribution of antibiotic resistance genes in recirculating aquaculture system. Ecotoxicol. Environ. Saf. 2019, 184, 109631. [Google Scholar] [CrossRef]
Figure 1. Sampling sites in the marine aquaculture area of Zhuanghe.
Figure 1. Sampling sites in the marine aquaculture area of Zhuanghe.
Ijerph 17 01446 g001
Figure 2. Temporal distribution of comprehensive scores. (a) Pond aquaculture area; (b) cage aquaculture area. Note: the top, bottom and middle lines of the Box plot represent the upper and lower quartiles and the median, respectively; the vertical part of the box body extending upward and downward represents the range of data distribution.
Figure 2. Temporal distribution of comprehensive scores. (a) Pond aquaculture area; (b) cage aquaculture area. Note: the top, bottom and middle lines of the Box plot represent the upper and lower quartiles and the median, respectively; the vertical part of the box body extending upward and downward represents the range of data distribution.
Ijerph 17 01446 g002
Figure 3. Spatial distributions of the comprehensive scores of principal components in aquaculture area. (a) Principal component comprehensive score in May; (b) principal component comprehensive score in August; (c) principal component comprehensive score in NovemberDriving Factors Analysis
Figure 3. Spatial distributions of the comprehensive scores of principal components in aquaculture area. (a) Principal component comprehensive score in May; (b) principal component comprehensive score in August; (c) principal component comprehensive score in NovemberDriving Factors Analysis
Ijerph 17 01446 g003
Table 1. Primer sequences of antibiotic resistance genes (ARGs).
Table 1. Primer sequences of antibiotic resistance genes (ARGs).
ARGs NamePrimers Sequences (5’-3’)Product Length (bp)Annealing TemperatureReference
florF: TATCTCCCTGTCGTTCCAG
R: AGAACTCGCCGATCAATG
39958 °C[29]
sul1F: CACCGGAAACATCGCTGCA
R: AAGTTCCGCCGCAAGGCT
15960 ℃[30]
sul2F: TCATCTGCCAAACTCGTCGTTA
R: GTCAAAGAACGCCGCAATGT
10556 ℃[31]
tetBF: CGAAGTAGGGGTTGAGACGC
R: AGACCAAGACCCGCTAATGAA
19256 ℃[32]
tetMF: ACAGAAAGCTTATTATATAAC
R: TGGCGTGTCTATGATGTTCAC
17160 ℃[33]
qnrSF: ACGACATTCGTCAACTGCAA
R: TAAATTGGCACCCTGTAGGC
41756 ℃[34]
ermBF: GATACCGTTTACGAAATTGG
R: GAATCGAGACTTGAGTGTGC
36458 ℃[35]
Table 2. Water quality variables in pond and cage culture areas during different sampling seasons.
Table 2. Water quality variables in pond and cage culture areas during different sampling seasons.
Sampling TypeMonthpHT (°C)SALDO (mg/L)MRP
(mg P/L)
COD (mg/L)Chla (μg/L)DIN
(mg N/L)
ARGs (%)
PondMay8.41 ± 0.4526.0 ± 0.731.3 ± 0.87.44 ± 1.440.095 ± 0.0613.0 ± 1.52.03 ± 1.600.100 ± 0.08622% ± 11%
Aug8.13 ± 0.2320.7 ± 0.528.90 ± 0.538.34 ± 1.290.178 ± 0.1507.9 ± 3.14.80 ± 2.060.027 ± 0.02527% ± 8%
Nov8.41 ± 0.1410 ± 0.231.1 ± 0.89.79 ± 0.690.042 ± 0.0492.2 ± 1.63.13 ± 1.610.080 ± 0.57720% ± 7%
CageMay8.10 ± 0.0516.2 ± 0.630.9 ± 0.29.26 ± 0.290.013 ± 0.0052.3 ± 1.62.54 ± 1.880.800 ± 0.99679% ± 12%
Aug7.97 ± 0.0624.3 ± 0.431.0 ± 0.28.24 ± 0.720.032 ± 0.0142.5 ± 1.31.93 ± 1.100.053 ± 0.02777% ± 16%
Nov8.06 ± 0.0313.0 ± 0.131.2 ± 0.08.16 ± 0.220.014 ± 0.0033.3 ± 1.44.26 ± 1.330.042 ± 0.008100% ± 0%
SAL: salinity; DO: dissolved oxygen; MRP: molybdate-reactive phosphorus; COD: chemical oxygen demand; DIN: dissolved inorganic nitrogen.
Table 3. Eigenvalues, variance contribution and accumulated contribution rate of principal components analysis (PCA).
Table 3. Eigenvalues, variance contribution and accumulated contribution rate of principal components analysis (PCA).
Sampling TypePrincipal ComponentEigenvalueVariance ContributionAccumulated Contribution Rate
Pond culture12.8631.7931.79
21.8820.8852.68
31.2714.1466.82
Cage culture13.1635.0835.08
22.3225.8260.90
31.0912.0982.99
Table 4. Various factors loadings in principal components.
Table 4. Various factors loadings in principal components.
FactorPond AquacultureCage Aquaculture
PC1PC2PC3PC1PC2PC3
pH −0.52 −0.37 −0.33 0.71 0.01 0.39
SAL −0.86 −0.10 −0.09 0.24 0.87 −0.28
T −0.04 0.78 0.13 −0.76 −0.30 −0.26
DO 0.03 0.90 0.16 −0.01 0.32 0.83
COD 0.72 0.26 0.27 0.03 0.55 0.36
Chla 0.74 −0.22 −0.10 0.80 −0.19 −0.22
MRP 0.33 0.52 0.10 −0.57 −0.21 −0.63
DIN −0.29 0.37 −0.76 0.47 −0.74 0.21
ARGs −0.02 0.27 0.90 0.57 0.68 −0.06
Table 5. Spearman’s correlation coefficient between the factors and principal components (PCs).
Table 5. Spearman’s correlation coefficient between the factors and principal components (PCs).
FactorPond AquacultureCage Aquaculture
PC1PC2PC3PC1PC2PC3
pH−0.56 **−0.47 *0.330.69 **0.04−0.45 *
T−0.170.73 **0.05−0.75 **−0.61 **0.09
SAL−0.86 **−0.150.150.270.85 **0.25
DO0.020.91 **0.190.020.400.83 **
MRP0.270.57 **0.00−0.49 *−0.180.58 **
COD0.71 **0.40 *−0.320.100.56 **−0.40
Chla0.74 **−0.220.030.77 **0.160.17
DIN−0.260.230.65 **0.19−0.40−0.59 **
ARGs0.130.38 *−0.82 **0.61 **0.83 **0.07
Note: * and ** show significant differences (p < 0.05 and p < 0.01, respectively).

Share and Cite

MDPI and ACS Style

Zhang, X.; Zhang, Y.; Zhang, Q.; Liu, P.; Guo, R.; Jin, S.; Liu, J.; Chen, L.; Ma, Z.; Liu, Y. Evaluation and Analysis of Water Quality of Marine Aquaculture Area. Int. J. Environ. Res. Public Health 2020, 17, 1446. https://doi.org/10.3390/ijerph17041446

AMA Style

Zhang X, Zhang Y, Zhang Q, Liu P, Guo R, Jin S, Liu J, Chen L, Ma Z, Liu Y. Evaluation and Analysis of Water Quality of Marine Aquaculture Area. International Journal of Environmental Research and Public Health. 2020; 17(4):1446. https://doi.org/10.3390/ijerph17041446

Chicago/Turabian Style

Zhang, Xianyu, Yingqi Zhang, Qian Zhang, Peiwu Liu, Rui Guo, Shengyi Jin, Jiawen Liu, Lei Chen, Zhen Ma, and Ying Liu. 2020. "Evaluation and Analysis of Water Quality of Marine Aquaculture Area" International Journal of Environmental Research and Public Health 17, no. 4: 1446. https://doi.org/10.3390/ijerph17041446

APA Style

Zhang, X., Zhang, Y., Zhang, Q., Liu, P., Guo, R., Jin, S., Liu, J., Chen, L., Ma, Z., & Liu, Y. (2020). Evaluation and Analysis of Water Quality of Marine Aquaculture Area. International Journal of Environmental Research and Public Health, 17(4), 1446. https://doi.org/10.3390/ijerph17041446

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop