Next Article in Journal
Spatiotemporal Assessment of Water Pollution for Beira Lake, Sri Lanka
Previous Article in Journal
Modelling Trends in Urban Flood Resilience towards Improving the Adaptability of Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Assessment of Ecological Quality Combining Biological and Environmental Data in the Yellow River Estuary

1
College of Oceanography and Ecological Science, Shanghai Ocean University, Shanghai 201306, China
2
State Environmental Protection Key Laboratory of Coastal Ecosystem, National Marine Environmental Monitoring Center, Dalian 116023, China
3
Observation and Research Station of Bohai Strait Eco-Corridor, First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(11), 1615; https://doi.org/10.3390/w16111615
Submission received: 22 April 2024 / Revised: 1 June 2024 / Accepted: 3 June 2024 / Published: 5 June 2024

Abstract

:
The integrated assessment of ecological quality in estuarine ecosystems holds significant importance for environmental management. Previous monitoring programs predominantly focused on environmental data, lacking a comprehensive quality assessment approach. To address this gap, this study aimed to integrate environmental factors with macrofaunal community information to evaluate the ecological quality status of the Yellow River Estuary. A total of 13 stations were routinely monitored in August for four consecutive years to collect environmental and biological data. Candidate indicators were screened based on variation coefficients, distribution ranges, and redundancy analysis, identifying 16 indicators belonging to three categories (i.e., seawater, sediment, and biology). The model fit and the interrelationship of the components were determined using structural equation modelling (SEM). The main results were as follows. (1) A total of 144 macrofaunal taxa, belonging to eight animal phyla and 98 families, were identified, with a dominance of Annelida (37.8%) and Mollusca (33.3%). The environmental variables most strongly correlated with the macrofaunal community were TOC, DO, Cd, and Md. (2) NO2 and heavy metals represented the two most direct factors of environmental pollution, while the factor load of biodiversity indices (H’, J, and D) was large in the biology category. (3) The evaluation results indicated that 78.85% of the total samples were between the average and upper levels of ecological quality, but only 7.69% of samples were at the “high” level. The framework system for the evaluation of ecological quality constructed in this study provides a theoretical and practical basis for the evaluation of the effectiveness of conservation management of the Yellow River Estuary.

1. Introduction

Estuaries and coastal ecosystems are specific habitats situated at the interface of land and sea, characterized by highly dynamic natural disturbances [1,2]. They serve a critical function in delivering a variety of ecosystem services, including regulatory (e.g., as carbon sinks and feeding areas), provisioning (e.g., as fishery and aquaculture resources), and cultural (e.g., as areas for recreational activities and sources of local ecological knowledge) services [3,4]. However, in recent years, estuarine ecosystems have been increasingly threatened by a variety of anthropogenic activities, leading to habitat deterioration and biodiversity degradation [5]. Consequently, routine assessments of ecological quality are essential for effective management aimed at enhancing deteriorating estuarine ecosystems worldwide [6,7,8,9,10]. These assessments provide important information for the evaluation of the effectiveness of nature-based solutions (NBSs) [11]. As both natural and anthropogenic stressors can affect ecosystems simultaneously, it is essential to distinguish their interrelation with observed effects and properly account for the inherent complexity of ecological data [12,13,14].
Traditional evaluations of habitat quality in estuarine waters have predominantly depended on chemical parameters, such as nitrogen, phosphorus, heavy metals, and organic pollutants [15,16,17]. This approach entails the comparison of the content of these substances with environmental quality standards to determine the quality level [18]. However, conducting a comprehensive evaluation solely based on chemical data presents significant challenges, because such data fail to capture the biological impacts resulting from the mixing of multiple pollutants. Fortunately, the European Water Framework Directive (WFD) provides a roadmap for the establishment of integrative standards, combining hydro-morphological elements, physicochemical elements, and biological quality elements [19,20]. It requires that member states establish specific ecological quality standards and targets, thus ensuring the attainment of a “good status” for water bodies [21]. Accordingly, numerous tools have been developed to assess the coastal ecological quality status (EcoQs). Among these, macrobenthos-based indices have been widely accepted by the scientific community [22,23,24,25].
Macrofauna can be categorized into different ecological groups based on their sensitivity and tolerance to disturbance and pollution. Several benthic biotic indices (BBIs), based on the relative abundance of these ecological groups, have been developed, including AZTI Marine Biotic Index (AMBI) [26], Multivariate AZTI Marine Biological Index (M-AMBI) [27], and BENTIX [28]. Although they were initially designed to assess the ecological quality of coastal areas in Europe, they have been widely applied, either individually or in combination, to coastal areas in other countries. In China, AMBI and M-AMBI have been utilized in the estuaries of the Yellow River, Yangtze River, and Xiaoqing River [29,30,31]; similarly, AMBI and BENTIX have found application in the intertidal zone of subtropical islands [32].
The monitoring programs described in most ecological quality assessment guidelines for coastal waters primarily concentrate on environmental parameters. When biological data are included, a comprehensive quality assessment method based on multiple parameters has yet to be developed. To fill this gap, this study aimed to explore the importance of both abiotic and biotic information, for assessing the ecological quality of the Yellow River Estuary in an integrative way. We aim to (1) characterize the temporal patterns of macrofauna for four consecutive years; (2) identify the set of environmental variables that better explains biological community patterns; and (3) develop an integrative method to assess ecological quality status using structural equation modeling (SEM). We hypothesized that the biological data contribute equally to the ecological quality evaluation compared with the seawater and sediment variables using SEM.

2. Materials and Methods

2.1. Study Area

The Yellow River, which flows westward into the Bohai Sea, is China’s second longest river and holds the position of the world’s second largest river by sediment load [33]. Over the past several thousand years, it has accounted for 6% of the worldwide sediment flux into the ocean [34]. Like many large rivers in other countries, the Yellow River has a series of dams and reservoirs, which were constructed to support local industrial, urbanized, and economic developments [35]. Consequently, the Yellow River has experienced a rapid decline in water discharge and sediment load. Since 2002, the water-sediment regulation scheme (WSRS) has been employed to tackle the issue of inadequate water supply for sediment transport in the Yellow River [36]. This scheme releases a substantial volume of water over a short period, representing over 60% of the yearly sediment load and 50% of the annual nutrient flux [37].

2.2. Sampling and Laboratory Analyses

Thirteen stations were routinely sampled for macrofauna investigation in August from 2019 to 2022 in the Yellow River Estuary (Figure 1). Environmental and biological samples were concurrently collected. Bottom seawater samples were collected at each station by using Niskin bottles. Seawater environmental factors including salinity (Sal), dissolved oxygen (DO, mg/L), and pH were measured in situ using a YSI 6600 system (YSI ProPlus Incorporated, USA). Seawater samples were taken to the laboratory to analyze the nitrite nitrogen (NO2-N), nitrate nitrogen (NO3-N), ammonia nitrogen (NH4-N), and phosphate (PO4-P) [15]. In compliance with the China National Standard [38], two replicate samples were obtained using a Van Veen grab sampler (0.1 m2). Upon sieving with a 0.5 mm mesh size, the biological samples were immediately preserved in a 5% buffered formaldehyde solution. Within the laboratory, all macrofaunal organisms were sorted using a stereomicroscope. Individuals were identified to the most specific taxonomic level achievable. The species names were validated against the database of the World Register of Marine Species (http://www.marinespecies.org/) accessed on 5 August 2023. Taxon density was recorded for each station by counting all individuals collected. Meanwhile, 100 g of sediment were taken to measure the total organic carbon (TOC), sulfide (Sul), and heavy metal contents (Cu, Pb, Zn, Cd, Hg, and Cr). Sediment grain size was analyzed using a laser particle size analyzer (Mastersizer 2000; Malvern, UK) [18].

2.3. Data Analyses

Univariate analysis was performed using the PRIMER v7 package [39]. The main macrofaunal community parameters were the species richness (S: number of taxa), abundance (A: number of individuals/m2), dominance, Shannon–Wiener (H’, log2), Simpson (D), and Pielou (J) indices. The AMBI and M-AMBI indices were determined using AZTI’s AMBI software (ver. 5.0) (http://ambi.azti.es/) and the species list from December 2022, following the guidelines. Additionally, the BENTIX index was calculated based on the relative percentages of the ecological groups to which the species belonged. Regarding the temporal variation, differences among different years were evaluated using Student’s t-test.
For multivariate analysis, macrofaunal abundances underwent square-root transformation to diminish the influence of dominant taxa. Significant differences in macrofaunal community composition in different years were evaluated via principal coordinate analysis (PCoA) based on the Bray–Curtis distance using the “Vegan” packages. The influence of environmental factors on macrofaunal community structures was analyzed through redundancy analysis (RDA) with the “Vegetation” packages. Concurrently, distance-based linear modeling (DistLM) was employed to evaluate the relative contributions of abiotic parameters to the observed variability in macrofaunal community structures [34]. The model construction employed the best combination of predictors according to the “Best” procedure, utilizing AIC (Akaike information criteria).

2.4. Development of the Integrated Ecological Quality Index

Initially, three important components of the coastal ecological integrity framework system, i.e., water, sediment, and biology, were selected from the perspective of the basic characteristics of estuarine ecosystems. Then, the ecological quality indicators were screened and determined. Next, the standardized method was used to unify the indicator dimensions, and the factors loading in the structural equation model were used to calculate the indicator weights. Finally, an integrated ecological quality index was calculated via the weighted summation method, and the results were graded according to the equidistant method.
(1)
Screening and identification of evaluation indicators
To ensure the uniqueness of the information provided by the selected indicators, this study employed variability analyses, normality tests, collinearity diagnostics, and redundancy tests. These analyses validated the representativeness and non-redundancy of each selected indicator. The analysis of candidate indicators’ variability and distribution range revealed that most water quality indicators exhibited high variability (CV > 50%), signifying significant fluctuations. Conversely, sediment indicators, with the exception of Sul and Hg, displayed low variability, indicative of minor fluctuations. Given the non-normal distribution of most indicators (p < 0.05), spearman correlation analysis and collinearity diagnostics were conducted for redundancy testing in the candidate indicators. The results indicated that the M-AMBI index demonstrated high correlations with several indicators (p > 0.7), and significant collinearity was observed among the H’, J, D, and M-AMBI indices (VIF > 10). Consequently, the M-AMBI indicator was excluded from the candidate list. Consequently, the integrated ecological quality index in this study consisted of 16 indicators.
(2)
Standardization of evaluation indicators
To comprehensively compare how each indicator represents ecological quality, we categorized the indicators into two types: positive indicators and negative indicators. Positive indicators imply that the higher the original value, the better the estuarine ecological quality, while negative indicators indicate that the higher the original value, the worse the estuarine ecological quality. Each indicator is standardized according to the following formulae.
Positive indicator:
z = x i x m i n + x m a x x m i n 5 × z 1 , x m i n + x m a x x m i n 5 × z
Negative indicators:
z = x i x m a x x m a x x m i n 5 × z , x m a x x m a x x m i n 5 × z 1
xi represents any given data point for each indicator, xmax is the maximum value for each indicator, xmin is the minimum value for each indicator, and z is 1, 2, 3, 4, 5.
(3)
Weight calculation for evaluation indicators
According to the factor load in the structural equation model, the weight wa of each indicator was calculated as follows:
w i = q i q t
w j = q j q t                        
w a = w i × w j                  
wi and wj correspond to the first-order and second-order weights, respectively; qi is the loading of the observed variable on the first-order latent variable, qj is the loading on the second-order latent variable, and qt is the sum of the absolute values of all loadings on that latent variable; wa represents the total weight of each indicator (Table 1).
(4)
Determination of ecological quality status
We used the weighted sum method to synthesize multiple indicators to obtain a comprehensive value that could reflect the ecological quality status. The equidistant scoring method was used to divide the scores obtained from small to large into five grades, namely: high, good, moderate, poor, and bad (Table 2).

3. Results

3.1. Environmental Characteristics

In this study, the percentage of silt and clay was high in the majority of stations in the Yellow River Estuary. There was no obvious temporal variation in Sal, NH4-N, NO3-N, PO4-P, Sul, Cu, or Zn, while DO, NO2-N, Md, TOC, Cd, Cr, Hg, and Pb varied significantly (Figure 2, Table S1). The sediment properties (e.g., TOC, slit, and sand) exhibited higher spatial variations than the water variables (e.g., DO and Sal). In general, the farther offshore, the higher the silt and clay content, while the TOC exhibited similar patterns.

3.2. Macrofaunal Community

A total of 144 macrofaunal species were collected from the Yellow River Estuary during the four years (Table S2). Annelida had the highest number of species, with 54 species, accounting for 37.8% of the total species number, followed by Mollusca, with 48 species, accounting for 33.3% of the total species number. Others included twenty-nine Arthropoda, five Echinodermata, four Chordata, two Nemertea, one Platyhelminthes, and one Brachiopoda. The average species number was 78, and the average density ranged from 120 to 172 ind./m2 (Figure 3a). The community was dominated by Mollusca in 2020, whereas it was dominated by Annelida in the other three years (Figure 3b). The composition of dominant species showed great temporal variation (Figure 3b, Table 3). A Venn plot showed that there were 39 species in common between 2019 and 2020, 45 species in common between 2020 and 2021, and 65 species in common between 2021 and 2022 (Figure 3c). To further clarify the variability of community species composition across years, PCoA analysis was used to measure the similarity of the species composition over the four years, and the results showed that communities from 2020 distinctly clustered compared to other years, and they were separated by the PC2 axis (Figure 3d).

3.3. Influence of Environmental Data on Macrofaunal Assemblages

The relationship between the macrofaunal communities and the environmental factors is demonstrated in the RDA ordination diagram (Figure 4). Monte Carlo tests for the all ordination axes were highly significant (p < 0.05). The first two axes of the RDA explained 37.23% of the total variation. The first axis was positively correlated with Md, TOC, Hg, and Pb but negatively correlated with DO and NO2-N. It was obvious that the macrofaunal communities structure of the Yellow River Estuary in different years was affected by different environmental factors. According to DistLM, the environmental variables which significantly correlated with the whole macrofaunal community (identified at the species level) were mainly TOC (5.2% of variation explained), DO (4.8%), and Cd (4.6%) (Table 4), while the Md, NO2-N, NO3-N, NH4-N, Hg, and Sal were also significant. The best model (AIC = 418.5, R2 = 13.5) was obtained for the combination of DO, NH4-N, and Cd (Table 3). When the three main phyla (Annelida, Arthropoda, and Mollusca) were analyzed separately, DistLM showed that variables explained a higher proportion of variability in Arthropoda, followed by Annelida and, finally, Mollusca (Table S3); therefore, the Annelida phylum can be concluded to best reflect the environmental conditions in the Yellow River Estuary.

3.4. Integrated Ecological Quality Evaluation

In this study, we utilized the SEM to construct an estuarine ecological quality assessment model, validating its fit and obtaining the intrinsic connections among various secondary indicators. The results of the data verification indicated that the primary indicators of water, sediment, and biology all had Cronbach’s alphas above 0.6 (Table S4), confirming data reliability and fulfilling the modeling requirements for constructing the integrated ecological quality model. The model fit results (Table 5) indicated optimal model fitting according to absolute fit indices like chi-square/degrees of freedom (CMIN/DF = 1.126) and root-mean-square error of approximation (RMSEA = 0.05). Incremental fit indices like the comparative fit index (CFI = 0.968) also reached optimal standards, suggesting that the constructed ecological quality evaluation model for the estuary closely resembles the real model, with a high overall fit. These validation results indicated that the model outputs could effectively explain the interrelations between ecological quality and secondary indicators in the Yellow River Estuary. The path analysis results for these secondary indicators and primary indicators indicated significant relationships, demonstrating that the selected secondary indicators could be effectively applied to ecological quality assessment (Figure 5). Using the weighted sum method for assessing ecological quality based on SEM, our results showed that the ecological quality values for the study area ranged from 1.70 to 4.59. Although 78.85% of the total samples were between the average and upper levels of ecological quality, only 7.69% of samples were at the “high” level (Table 2, Table S5).

4. Discussion

In this research, we conducted a routine four-year investigation of macrobenthic communities in the Yellow River Estuary. In total, 144 species were collected (62 species in 2019, 75 species in 2020, 81 species in 2021, and 97 species in 2022). These species belong to 125 genera, distributed across 98 families, 39 orders, 14 classes, and 8 phyla. These results are consistent with those of Li et al. [40], who identified 88 species of macrofauna at 26 stations in this area. Our findings also indicated an increasing trend year on year in the species number of macrofauna, which likely reflects the successful impact of ecological protection and restoration initiatives within the Yellow River Estuary. In most years, the macrofaunal communities were dominated by polychaetes in this study, consistent with findings from other studies in the Yellow River Estuary and adjacent sea areas such as Laizhou Bay [30,41]. Compared with historical data, the macrofaunal communities have shown a trend toward miniaturization in the Bohai Sea. Large individual echinoderms, crustaceans, and mollusks have gradually been replaced by smaller individual polychaetes and mollusks [42,43]. This trend may be associated with marine eutrophication, overfishing, and continuous anthropogenic disturbances [44]. However, in 2020, mollusks were dominant in relative abundance (71%), primarily due to the influence of Arcuatula senhousia collected at S13, which exhibited a high abundance of up to 6520 ind./m2. Overall, the sediment in this area is characterized by a higher clay component and a smaller median particle size, conducive to the habitation of molluscan species. Additionally, the explosive reproduction and gregarious habituation of Arcuatula senhousia made the large quantity collected at S13 an expected phenomenon. Similar findings were reported by Joel et al. [45], in whose study a significant number of Arcuatula senhousia were also collected at a certain station.
The patterns of macrofaunal communities in this study were primarily explained by TOC, DO, Cd, and Md. These variables are also supported by the findings of previous research [46,47]. Sediment granularity is considered one of the main factors regulating soft-bottom macrofaunal communities. Blanchet et al. (2008) [48] also reported granularity, TOC, and heavy metals as significant variables in explaining the distribution of mollusks, crustaceans, and polychaetes in estuarine environments. The importance of sediment grain size distribution to the structure of macrobenthic animal communities underscores the need to consider the impact of human activities on sediment transport systems, as well as terrestrial organic matter and nutrient loads, in the management of estuarine ecosystems.
Estuarine ecological quality evaluation is a complex system project, and although many scholars have made significant progress in this area, considerable controversy remains over the objectivity and accuracy of their evaluation results [8]. Our study aimed to optimize the monitoring and evaluation of estuarine ecological quality by establishing an evaluation framework, selecting evaluation indicators, and adopting objective quantification methods. The construction of an index system represented the initial stage in the assessment of ecological quality, with the selection of indices that accurately mirrored the critical elements of the estuarine ecosystem being imperative. Commonly, macrobenthic communities are recognized as robust indicators of ecosystem health due to their significant associations with sediment processes, which are concurrently connected to water column dynamics [7,10,49,50]. Hence, in this study, five indicators (H’, J, D, AMBI, and Bentix) were combined to describe the characteristics of the benthic macroinvertebrate community rather than a universal index of other indicators. Given the scope of national routine monitoring programs implemented across China, the index system devised for ecological quality assessment is appropriate for application across the general coastal waters of China. Future validations of the ecological quality evaluation model should incorporate both spatial and temporal dimensions to enhance its robustness and applicability.
Structural equation modeling, which is used in this research, has been the fastest growing branch of statistics over the past 30 years. Its major advantage lies in the analysis of latent variables that cannot be directly measured [51], allowing for the use of multi-scale indicators, including both biological and abiotic variables, to better represent the hierarchy and heterogeneity of complex ecosystems [52,53,54]. Indicator weights form a critical basis for evaluating the ecological quality of estuaries. SEM derives path coefficients from the intrinsic logical relationships in original data to establish weights, providing a new perspective for determining weights for various assessment indicators of estuarine ecological quality [55,56]. In this study, model fit indices such as χ2/df and RMSEA meet excellent standards, indicating that the estuarine ecological quality index model is well adapted to the 16 selected indicators. Among environmental factors, the sediment weights are higher than that of the water. Specifically, the heavy metal factor load of the sediment was high, and this represents the main pollution of estuaries. Coastal and estuarine ecosystems in China are now facing increasing metal pollution pressure [57,58], which affects the growth and reproduction of benthic organisms. The weights of the four water quality indicators (NO2-N > NO3-N > NH4-N > PO4-P) suggest that nitrogenous input might be more serious than phosphorous input in the study area. Urban industrial and agricultural nitrogen and phosphorus sewage inputs, along with sediment nitrogen and phosphorus releases, jointly affect the estuarine water quality [59,60]. Biological indices have been proven to be effective indicators of ecological quality. In this study, the high weights of the biodiversity indices D, H’, and J indicate that species richness and evenness can greatly influence the biological quality of the study area. This may be due to the uneven distribution characteristics of the macrofaunal community in the Yellow River Estuary. Furthermore, high sedimentation rates in the Yellow River Estuary lead to deteriorating habitat conditions for benthic organisms, reduced food availability, thus affecting biodiversity and ecosystem stability. Overall, the results obtained in this study indicate that effective measures should be taken to control nutrient and heavy metal pollution, to restore biodiversity in the Yellow River Estuary.

5. Conclusions

In conclusion, this research examines the macrofaunal communities and ecological quality in the Yellow River Estuary. Generally, Polychaeta was identified as the dominant taxonomic group within the macrofaunal community, and TOC, DO, Cd, and Md were significant environmental factors that influence the community. The structural equation modeling reveals that environmental factors are highly weighted, particularly emphasizing the substantial disruption caused by nutrients and heavy metals on estuarine ecological quality. The factor load of biodiversity indices (H’, J, and D) was high, indicating the important role of macrofaunal biodiversity in the studied area. The evaluation results indicate that 78.85% of the total samples were between the average and upper levels of ecological quality, but only 7.69% of samples were at the “high” level. Future efforts should focus on the continuous monitoring of environmental factors and macrobenthic community dynamics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w16111615/s1, Table S1. Environmental factors; Table S2. Macrobenthos community abundance at 13 sampling stations in the Yellow River Estuary from 2019 to 2022; Table S3. Results of DistLM (distance-based linear modelling) analysis for macrofauna. The variables explaining macrofaunal community are included. Analyses were conducted for Annelida, Arthropoda, and Mollusca separately. The table shows significance levels for each predictor variable and the proportion of variation explained for any faunal assemblage. Only variables with significant results (p < 0.05) in marginal tests are shown. Overall best models are also included. AIC: Akaike information criterion; SS (trace): sum of squares of sequential test, Prop.: proportion of variability explained by each factor (in the marginal tests without coaction of factors), Sal: salinity, DO: dissolved oxygen, TOC: total organic carbon; Table S4. The statistics for data reliability and validity test; Table S5. The integrated ecological quality index for each sample.

Author Contributions

Conceptualization, X.G. and W.L.; methodology, X.G.; software, X.G.; validation, H.L. and Y.Z.; formal analysis, X.G.; investigation, H.S.; resources, H.S.; data curation, X.G.; writing—original draft preparation, X.G.; writing—review and editing, X.G.; visualization, Y.L.; supervision, H.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Open Fund from the Observation and Research Station of Bohai Strait Eco-Corridor grant number No. BH202201.

Data Availability Statement

The original contributions presented in the study are included in the Supplementary Material, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kim, H.G.; Song, S.J.; Bae, H.; Noh, J.; Lee, C.; Kwon, B.O.; Lee, J.H.; Ryu, J.; Khim, J.S. Natural and anthropogenic impacts on long-term meiobenthic communities in two contrasting nearshore habitats. Environ. Int. 2020, 134, 105200. [Google Scholar] [CrossRef]
  2. Feebarani, J.; Joydas, T.V.; Damodaran, R.; Borja, Á. Benthic quality assessment in a naturally- and human-stressed tropical estuary. Ecol. Indic. 2016, 67, 380–390. [Google Scholar] [CrossRef]
  3. Barbier, E.B.; Hacker, S.D.; Kennedy, C.; Koch, E.W.; Stier, A.C.; Silliman, B.R. The value of estuarine and coastal ecosystem services. Ecol. Monogr. 2011, 81, 169–193. [Google Scholar] [CrossRef]
  4. Fulford, R.S.; Russell, M.; Hagy, J.D.; Breitburg, D. Managing estuaries for ecosystem function. Glob. Ecol. Conserv. 2020, 21, e00892. [Google Scholar] [CrossRef]
  5. Halpern, B.S.; Walbridge, S.; Selkoe, K.A.; Kappel, C.V.; Micheli, F.; d’Agrosa, C.; Bruno, J.F.; Casey, K.S.; Ebert, C.; Fox, H.E. A global map of human impact on marine ecosystems. Science 2008, 319, 948–952. [Google Scholar] [CrossRef]
  6. Pelletier, M.; Gillett, D.; Hamilton, A.; Grayson, T.; Hansen, V.; Leppo, E.; Weisberg, S.; Borja, A. Adaptation and application of multivariate AMBI (M-AMBI) in US coastal waters. Ecol. Indic. 2018, 89, 818–827. [Google Scholar] [CrossRef]
  7. Mulik, J.; Sukumaran, S.; Dias, H.Q. Can the ecological status of three differentially impacted monsoonal tropical estuaries in NW India, be adequately assessed by a common estuarine benthic index? Ecol. Indic. 2020, 119, 106807. [Google Scholar] [CrossRef]
  8. Borja, A.; Dauer, D.M. Assessing the environmental quality status in estuarine and coastal systems: Comparing methodologies and indices. Ecol. Indic. 2008, 8, 331–337. [Google Scholar] [CrossRef]
  9. Borja, A.; Dauer, D.; Diaz, R.; Llansó, R.; Muxika, I.; Rodríguez, J.; Schaffner, L. Assessing estuarine benthic quality conditions in Chesapeake Bay: A comparison of three indices. Ecol. Indic. 2008, 8, 395–403. [Google Scholar] [CrossRef]
  10. Yang, C.; Song, H.; Sun, Y.; Xie, P.; Liu, Y.; Li, H. Seasonal implications for taxonomic sufficiency to simplify M-AMBI methodology in the coastal area adjacent to a eutrophic estuary. Acta Oceanol. Sin. 2023, 42, 108–116. [Google Scholar] [CrossRef]
  11. Fang, X.; Li, J.; Ma, Q.; Zhou, R.; Du, S. A quantitative review of nature-based solutions for urban sustainability (2016–2022): From science to implementation. Sci. Total Environ. 2024, 927, 172219. [Google Scholar] [CrossRef]
  12. Couce, E.; Engelhard, G.H.; Schratzberger, M. Capturing threshold responses of marine benthos along gradients of natural and anthropogenic change. J. Appl. Ecol. 2020, 57, 1137–1148. [Google Scholar] [CrossRef]
  13. Mulik, J.; Sukumaran, S.; Srinivas, T.; Vijapure, T. Comparative efficacy of benthic biotic indices in assessing the Ecological Quality Status (EcoQS) of the stressed Ulhas estuary, India. Mar. Pollut. Bull. 2017, 120, 192–202. [Google Scholar] [CrossRef]
  14. Mulik, J.; Sukumaran, S.; Dias, H.Q. Is the benthic index AMBI impervious to seasonality and data transformations while evaluating the ecological status of an anthropized monsoonal estuary? Ocean Coast. Manag. 2020, 186, 105080. [Google Scholar] [CrossRef]
  15. Lin, H.; Li, H.; Yang, X.; Xu, Z.; Tong, Y.; Yu, X. Comprehensive investigation and assessment of nutrient and heavy metal contamination in the surface water of coastal Bohai Sea in China. J. Ocean Univ. China 2020, 19, 843–852. [Google Scholar] [CrossRef]
  16. Li, H.; Lin, L.; Ye, S.; Li, H.; Fan, J. Assessment of nutrient and heavy metal contamination in the seawater and sediment of Yalujiang Estuary. Mar. Pollut. Bull. 2017, 117, 499–506. [Google Scholar] [CrossRef]
  17. Li, H.; Ye, S.; Ye, J.; Fan, J.; Gao, M.; Guo, H. Baseline survey of sediments and marine organisms in Liaohe Estuary: Heavy metals, polychlorinated biphenyls and organochlorine pesticides. Mar. Pollut. Bull. 2017, 114, 555–563. [Google Scholar] [CrossRef]
  18. Li, H.; Gao, X.; Gu, Y.; Wang, R.; Xie, P.; Liang, M.; Ming, H.; Su, J. Comprehensive large-scale investigation and assessment of trace metal in the coastal sediments of Bohai Sea. Mar. Pollut. Bull. 2018, 129, 126–134. [Google Scholar] [CrossRef]
  19. Borja, A.; Chust, G.; Muxika, I. Forever young: The successful story of a marine biotic index. Adv. Mar. Biol. 2019, 82, 93–127. [Google Scholar]
  20. Hering, D.; Borja, A.; Carstensen, J.; Carvalho, L.; Elliott, M.; Feld, C.K.; Heiskanen, A.-S.; Johnson, R.K.; Moe, J.; Pont, D.; et al. The European Water Framework Directive at the age of 10: A critical review of the achievements with recommendations for the future. Sci. Total Environ. 2010, 408, 4007–4019. [Google Scholar] [CrossRef]
  21. Reyjol, Y.; Argillier, C.; Bonne, W.; Borja, A.; Buijse, A.D.; Cardoso, A.C.; Daufresne, M.; Kernan, M.; Ferreira, M.T.; Poikane, S.; et al. Assessing the ecological status in the context of the European Water Framework Directive: Where do we go now? Sci. Total Environ. 2014, 497–498, 332–344. [Google Scholar] [CrossRef]
  22. Dong, J.; Wang, X.; Bidegain, G.; Sun, X.; Bian, X.; Zhang, X. Assessment of the benthic ecological quality status (EcoQs) of Laizhou Bay (China) with an integrated AMBI, M−AMBI, BENTIX, BO2A and feeding evenness index. Ecol. Indic. 2023, 153, 110456. [Google Scholar] [CrossRef]
  23. Dong, J.; Sun, X.; Zhang, Y.; Zhan, Q.; Zhang, X. Assessing benthic habitat ecological quality using four benthic indices in the coastal waters of Sanshandao, Laizhou Bay, China. Ecol. Indic. 2021, 129, 107980. [Google Scholar] [CrossRef]
  24. Liang, J.; Ma, C.-W.; Kim, S.-K.; Park, S.-H. Assessing the Benthic Ecological Quality in the Intertidal Zone of Cheonsu Bay, Korea, Using Multiple Biotic Indices. Water 2024, 16, 272. [Google Scholar] [CrossRef]
  25. Dreujou, E.; McKindsey, C.W.; Grant, C.; Tréau de Coeli, L.; St-Louis, R.; Archambault, P. Biodiversity and Habitat Assessment of Coastal Benthic Communities in a Sub-Arctic Industrial Harbor Area. Water 2020, 12, 2424. [Google Scholar] [CrossRef]
  26. Borja, A.; Franco, J.; Pérez, V. A Marine Biotic Index to Establish the Ecological Quality of Soft-Bottom Benthos Within European Estuarine and Coastal Environments. Mar. Pollut. Bull. 2000, 40, 1100–1114. [Google Scholar] [CrossRef]
  27. Muxika, I.; Borja, Á.; Bald, J. Using historical data, expert judgement and multivariate analysis in assessing reference conditions and benthic ecological status, according to the European Water Framework Directive. Mar. Pollut. Bull. 2007, 55, 16–29. [Google Scholar] [CrossRef]
  28. Simboura, N.; Zenetos, A. Benthic indicators to use in Ecological Quality classification of Mediterranean soft bottom marine ecosystems, including a new Biotic Index. Mediterr. Mar. Sci. 2002, 3, 77–111. [Google Scholar] [CrossRef]
  29. Yan, J.; Sui, J.; Xu, Y.; Li, X.; Wang, H.; Zhang, B. Assessment of the benthic ecological status in adjacent areas of the Yangtze River Estuary, China, using AMBI, M-AMBI and BOPA biotic indices. Mar. Pollut. Bull. 2020, 153, 111020. [Google Scholar] [CrossRef]
  30. Li, B.; Li, X.; Bouma, T.J.; Soissons, L.M.; Cozzoli, F.; Wang, Q.; Zhou, Z.; Chen, L. Analysis of macrobenthic assemblages and ecological health of Yellow River Delta, China, using AMBI & M-AMBI assessment method. Mar. Pollut. Bull. 2017, 119, 23–32. [Google Scholar]
  31. Luo, X.; Lin, S.; Yang, J.; Shen, J.; Fan, Y.; Zhang, L. Benthic habitat quality assessment based on biological indices in Xiaoqing River estuary and its adjacent sea of Laizhou Bay, China. J. Ocean Univ. China 2017, 16, 537–546. [Google Scholar] [CrossRef]
  32. Sun, Y.; Chen, B.; Wu, H.; Huang, H.; Ma, Z.; Tang, K. Assessing benthic ecological status in subtropical islands, China using AMBI and Bentix indices. Estuar. Coast. Shelf Sci. 2018, 207, 345–350. [Google Scholar] [CrossRef]
  33. Bi, N.; Sun, Z.; Wang, H.; Wu, X.; Fan, Y.; Xu, C.; Yang, Z. Response of channel scouring and deposition to the regulation of large reservoirs: A case study of the lower reaches of the Yellow River (Huanghe). J. Hydrol. 2019, 568, 972–984. [Google Scholar] [CrossRef]
  34. Liu, S.; Li, L.; Zhang, G.; Liu, Z.; Yu, Z.; Ren, J. Impacts of human activities on nutrient transports in the Huanghe (Yellow River) estuary. J. Hydrol. 2012, 430-431, 103–110. [Google Scholar] [CrossRef]
  35. Wang, H.; Wu, X.; Bi, N.; Li, S.; Yuan, P.; Wang, A.; Syvitski, J.P.M.; Saito, Y.; Yang, Z.; Liu, S.; et al. Impacts of the dam-orientated water-sediment regulation scheme on the lower reaches and delta of the Yellow River (Huanghe): A review. Glob. Planet. Chang. 2017, 157, 93–113. [Google Scholar] [CrossRef]
  36. Zhang, J.; Li, F.; Lv, Q.; Wang, Y.; Yu, J.; Gao, Y.; Ren, Z.; Zhang, X.; Lv, Z. Impact of the water–sediment regulation scheme on the phytoplankton community in the Yellow River estuary. J. Clean. Prod. 2021, 294, 126291. [Google Scholar] [CrossRef]
  37. Li, G.; Sheng, L. Model of water-sediment regulation in Yellow River and its effect. Sci. China Technol. Sci. 2011, 54, 924–930. [Google Scholar] [CrossRef]
  38. GB/T 12763-2007; Specifications for Oceanographic SurVey-Part 3: Marine Meteorological Observations. China Standards Press: Beijing, China, 2017.
  39. Clarke, K.; Gorley, R. Getting Started with PRIMER v7; PRIMER-e: Plymouth, Plymouth Marine Laboratory: Boston, MA, USA, 2015; p. 20. [Google Scholar]
  40. Li, S.; Zhang, X.; Su, H.; Chen, S.; Chen, W.; Li, F.; Wang, T.; Sun, Z.; Weng, X. Taxonomic diversity of the macrobenthos and its relationship with environmental factors in the Huanghe River (Yellow River) estuary and adjacent seas, China. J. Sea Res. 2023, 194, 102402. [Google Scholar] [CrossRef]
  41. Wang, L.; Luo, X.; Yang, J.; Zhang, J.; Fan, Y. Assessing benthic habitat quality using biotic indices in the Laizhou Bay, China. Acta Oceanol. Sin. 2020, 39, 53–62. [Google Scholar] [CrossRef]
  42. Zhou, H.; Zhang, Z.; Liu, X.; Hua, E. Decadal change in sublittoral macrofaunal biodiversity in the Bohai Sea, China. Mar. Pollut. Bull. 2012, 64, 2364–2373. [Google Scholar] [CrossRef]
  43. Zhou, H.; Zhang, Z.; Liu, X.; Tu, L.; Yu, Z. Changes in the shelf macrobenthic community over large temporal and spatial scales in the Bohai Sea, China. J. Mar. Syst. 2007, 67, 312–321. [Google Scholar] [CrossRef]
  44. Xu, Y.; Sui, J.; Ma, L.; Li, X.; Wang, H.; Zhang, B. Temporal variation of macrobenthic community zonation over nearly 60 years and the effects of latitude and depth in the southern Yellow Sea and East China Sea. Sci. Total Environ. 2020, 739, 139760. [Google Scholar] [CrossRef]
  45. Harding, J.M.S.; Segal, M.R.; Reynolds, J.D. Location Is everything: Evaluating the effects of terrestrial and marine resource subsidies on an estuarine bivalve. PLoS ONE 2015, 10, e0125167. [Google Scholar] [CrossRef]
  46. Dong, J.; Zhao, L.; Sun, X.; Hu, C.; Wang, Y.; Li, W.-T.; Zhang, P.; Zhang, X. Response of macrobenthic communities to heavy metal pollution in Laoshan Bay, China: A trait-based method. Mar. Pollut. Bull. 2021, 167, 112292. [Google Scholar] [CrossRef]
  47. Maghsoudlou, A.; Momtazi, F.; Hashtroudi, M.S. Ecological Quality Status (EcoQs) of Chabahar sub-tropical bay based on multimetric macrobenthos-indexes approach: Response of bio-indexes to sediment structural/pollutant variables. Reg. Stud. Mar. Sci. 2020, 40, 101524. [Google Scholar] [CrossRef]
  48. Blanchet, H.; Lavesque, N.; Ruellet, T.; Dauvin, J.C.; Sauriau, P.G.; Desroy, N.; Desclaux, C.; Leconte, M.; Bachelet, G.; Janson, A.L.; et al. Use of biotic indices in semi-enclosed coastal ecosystems and transitional waters habitats—Implications for the implementation of the European Water Framework Directive. Ecol. Indic. 2008, 8, 360–372. [Google Scholar] [CrossRef]
  49. Mulik, J.; Sukumaran, S.; Tatiparthi, S. Factors structuring spatio-temporal dynamics of macrobenthic communities of three differently modified tropical estuaries. Mar. Pollut. Bull. 2019, 150, 110767. [Google Scholar] [CrossRef]
  50. Bae, H.; Lee, J.-H.; Song, S.J.; Ryu, J.; Noh, J.; Kwon, B.-O.; Choi, K.; Khim, J.S. Spatiotemporal variations in macrofaunal assemblages linked to site-specific environmental factors in two contrasting nearshore habitats. Environ. Pollut. 2018, 241, 596–606. [Google Scholar] [CrossRef] [PubMed]
  51. Tuominen, L.S.; Helle, S.; Helanterä, H.; Karell, P.; Rapeli, L.; Richmond, D.; Vuorisalo, T.; Brommer, J.E. Structural equation modeling reveals decoupling of ecological and self-perceived outcomes in a garden box social-ecological system. Sci. Rep. 2022, 12, 6425. [Google Scholar] [CrossRef] [PubMed]
  52. Capmourteres, V.; Anand, M. Assessing ecological integrity: A multi-scale structural and functional approach using Structural Equation Modeling. Ecol. Indic. 2016, 71, 258–269. [Google Scholar] [CrossRef]
  53. Fan, Y.; Chen, J.; Shirkey, G.; John, R.; Wu, S.R.; Park, H.; Shao, C. Applications of structural equation modeling (SEM) in ecological studies: An updated review. Ecol. Process. 2016, 5, 19. [Google Scholar] [CrossRef]
  54. Fernandes, A.C.P.; Sanches Fernandes, L.F.; Moura, J.P.; Cortes, R.M.V.; Pacheco, F.A.L. A structural equation model to predict macroinvertebrate-based ecological status in catchments influenced by anthropogenic pressures. Sci. Total Environ. 2019, 681, 242–257. [Google Scholar] [CrossRef]
  55. Garrido, M.; Hansen, S.K.; Yaari, R.; Hawlena, H. A model selection approach to structural equation modelling: A critical evaluation and a road map for ecologists. Methods Ecol. Evolut. 2022, 13, 42–53. [Google Scholar] [CrossRef]
  56. Grace, J.B. Instrumental variable methods in structural equation models. Methods Ecol. Evolut. 2021, 12, 1148–1157. [Google Scholar] [CrossRef]
  57. Wang, S.; Xu, X.; Sun, Y.; Liu, J.; Li, H. Heavy metal pollution in coastal areas of South China: A review. Mar. Pollut. Bull. 2013, 76, 7–15. [Google Scholar] [CrossRef]
  58. Gao, X.; Chen, C.-T.A. Heavy metal pollution status in surface sediments of the coastal Bohai Bay. Water Res. 2012, 46, 1901–1911. [Google Scholar] [CrossRef]
  59. Wang, B.; Xin, M.; Wei, Q.; Xie, L. A historical overview of coastal eutrophication in the China Seas. Mar. Pollut. Bull. 2018, 136, 394–400. [Google Scholar] [CrossRef] [PubMed]
  60. Rabalais, N.N.; Turner, R.E.; Díaz, R.J.; Justić, D. Global change and eutrophication of coastal waters. ICES J. Mar. Sci. 2009, 66, 1528–1537. [Google Scholar] [CrossRef]
Figure 1. Map showing the locations of macrofauna sampling stations in the Yellow River Estuary.
Figure 1. Map showing the locations of macrofauna sampling stations in the Yellow River Estuary.
Water 16 01615 g001
Figure 2. Boxplots illustrating the temporal variation in environmental factors in the Yellow River Estuary. Different letters above the bars indicate significant differences (t-test, p < 0.05), while the same letters indicate no significant difference.
Figure 2. Boxplots illustrating the temporal variation in environmental factors in the Yellow River Estuary. Different letters above the bars indicate significant differences (t-test, p < 0.05), while the same letters indicate no significant difference.
Water 16 01615 g002
Figure 3. The characteristics of macrofaunal communities from 2019 to 2022 in the Yellow River Estuary. Temporal variation of density, number of species (a), and relative abundance (b). Venn plot of shared and unique species (c). Principal coordinate analysis (PCoA) scatter plot based on Bray–Curtis distance matrix (d).
Figure 3. The characteristics of macrofaunal communities from 2019 to 2022 in the Yellow River Estuary. Temporal variation of density, number of species (a), and relative abundance (b). Venn plot of shared and unique species (c). Principal coordinate analysis (PCoA) scatter plot based on Bray–Curtis distance matrix (d).
Water 16 01615 g003
Figure 4. Graphic representation of the dbRDA, showing the ordination of the macrofaunal community and variables with respect to the first two axes.
Figure 4. Graphic representation of the dbRDA, showing the ordination of the macrofaunal community and variables with respect to the first two axes.
Water 16 01615 g004
Figure 5. The results of the model fitting for the structural equation of ecological quality index in the Yellow River Estuary (chi-square: 109.2, df: 97). The ovals, boxes, and circles represent latent, observed, and error variables, respectively. The numbers on the long arrows represent path coefficients (factor loadings), which use the standardized coefficient (Std. Coeff.) *** p < 0.001, ** p < 0.01, * p < 0.05.
Figure 5. The results of the model fitting for the structural equation of ecological quality index in the Yellow River Estuary (chi-square: 109.2, df: 97). The ovals, boxes, and circles represent latent, observed, and error variables, respectively. The numbers on the long arrows represent path coefficients (factor loadings), which use the standardized coefficient (Std. Coeff.) *** p < 0.001, ** p < 0.01, * p < 0.05.
Water 16 01615 g005
Table 1. The weight distribution of each indicator at all levels.
Table 1. The weight distribution of each indicator at all levels.
GoalComponentFirst-Order Weight (wi)IndicatorSecond-Order
Weight (wj)
Total Weight (wa)
Ecological qualityWater0.314NO20.3530.111
NO30.2780.087
NH40.2040.064
PO40.1650.052
Sediment0.386TOC0.1060.041
Sul0.0990.038
Cu0.2170.084
Zn0.1860.072
Cr0.1130.044
Cd0.1520.059
Pb0.1280.049
Biology0.300H’0.2410.072
J0.2620.079
D0.2980.089
AMBI0.0830.025
Bentix0.1160.035
Table 2. Classification of the grades of the ecological quality of the Yellow River Estuary.
Table 2. Classification of the grades of the ecological quality of the Yellow River Estuary.
GradeHighGoodModeratePoorBad
Ecological quality indexx ≥ 4.03.4 ≤ x < 4.02.8 ≤ x < 3.42.2 ≤ x < 2.8x < 2.2
Number of samples4122565
Proportion (%)7.6923.0848.0811.549.61
Table 3. The characteristics of dominant species in the macrofauna in the Yellow River Estuary.
Table 3. The characteristics of dominant species in the macrofauna in the Yellow River Estuary.
Dominant SpeciesEcological GroupDominance
2019202020212022
Aricidea fragilis WebsterEGI0.039+0.020+
Ringicula doliaris A. GouldEGI+0.035++
Goniada maculata ÖrstedEGII0.057+
Glycinde bonhourei gurjanovaeEGII0.0960.167
Arcuatula senhousiaEGIII0.116++
Sternaspis scutata RanzaniEGIII+0.025
Heteromastus filiformis ClaparèdeEGIV0.0760.0250.0770.089
Sigambra bassi HartmanEGIV0.031+0.033+
Note: the species is defined as a dominant species (Y ≥ 0.02). + and − indicate the dominant species present or absent during the different years.
Table 4. Results of DistLM (distance-based linear modelling) analysis. The variables explaining macrofaunal community (species abundances) are included. Table shows significance levels for each predictor variable and the proportion of variation explained. Only variables with significant results (p < 0.05) in marginal tests are shown. Overall best models are also included. AIC: Akaike information criterion; SS (trace): sum of squares of sequential test, Prop.: proportion of variability explained by each factor (in the marginal tests without coaction of factors), Sal: salinity, DO: dissolved oxygen, TOC: total organic carbon.
Table 4. Results of DistLM (distance-based linear modelling) analysis. The variables explaining macrofaunal community (species abundances) are included. Table shows significance levels for each predictor variable and the proportion of variation explained. Only variables with significant results (p < 0.05) in marginal tests are shown. Overall best models are also included. AIC: Akaike information criterion; SS (trace): sum of squares of sequential test, Prop.: proportion of variability explained by each factor (in the marginal tests without coaction of factors), Sal: salinity, DO: dissolved oxygen, TOC: total organic carbon.
VariablesSS (Trace)Pseudo-FpProportion (%)
TOC83072.720.0015.2
DO7668.22.500.0014.8
Cd7491.52.440.0014.6
Md6432.72.080.0074.0
NO2-N5833.81.880.0143.6
NO3-N5881.81.900.0123.6
NH4-N5583.31.800.0143.5
Hg5202.21.670.0353.2
Sal5031.61.610.0453.1
AICR2Variables
418.513.5DO, NH4-N, Cd
418.816.2DO, NO2-N, NH4-N, Cd
418.916.1Sal, DO, NH4-N, Cd
Table 5. Model fitting criteria and results.
Table 5. Model fitting criteria and results.
Classify Fit MeasuresAbsolute Fit MeasuresIncremental Fit Measures
CMIN/DFRMSEACFITLIIFI
Reference standard<2<0.08>0.9>0.9>0.9
Initial fit index1.6650.1140.8000.7620.810
Modified fit index1.1260.0500.9630.9550.966
Notes: CMIN/DF is the ratio of likelihood ratio to degrees of freedom; RMSEA is the root-mean-square error of approximation; CFI is the comparative fit index; TLI is the Tucker–Lewis index; IFI is the incremental fit index.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, X.; Li, W.; Zhang, Y.; Song, H.; Li, Y.; Li, H. Integrated Assessment of Ecological Quality Combining Biological and Environmental Data in the Yellow River Estuary. Water 2024, 16, 1615. https://doi.org/10.3390/w16111615

AMA Style

Gao X, Li W, Zhang Y, Song H, Li Y, Li H. Integrated Assessment of Ecological Quality Combining Biological and Environmental Data in the Yellow River Estuary. Water. 2024; 16(11):1615. https://doi.org/10.3390/w16111615

Chicago/Turabian Style

Gao, Xin, Wen Li, Yunlei Zhang, Hongjun Song, Ying Li, and Hongjun Li. 2024. "Integrated Assessment of Ecological Quality Combining Biological and Environmental Data in the Yellow River Estuary" Water 16, no. 11: 1615. https://doi.org/10.3390/w16111615

APA Style

Gao, X., Li, W., Zhang, Y., Song, H., Li, Y., & Li, H. (2024). Integrated Assessment of Ecological Quality Combining Biological and Environmental Data in the Yellow River Estuary. Water, 16(11), 1615. https://doi.org/10.3390/w16111615

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