Introduction
Plastic, a groundbreaking material first developed in the 19th century, has experienced a dramatic production growth since the 1950s, driven by population growth and rising consumption. By 2021, global plastic production had surpassed 390 million tons [1-3]. Plastic waste degrades into microplastics (MPs), defined as plastic particles with diameters <5 mm. Over the past few decades, these MPs have increasingly infiltrated both aquatic and terrestrial environments [1, 4, 5]. MPs can originate from primary sources, such as cosmetic products, and from secondary sources, where larger plastics degrade into smaller pieces through various processes [2, 6]. Notably, most environmental MPs originate from secondary sources, prompting bans on primary MPs, such as microbeads [7].
These small plastic particles, composed of various polymers and additives, pose threats to water quality and ecosystems by carrying harmful chemicals and absorbing/releasing pollutants and heavy metals [8]. Their global environmental concern status arises from their widespread presence and ecological hazards. Despite their ubiquity in the environment, the risks posed by MPs remain incompletely understood, necessitating focused studies of their sources, distribution, and impacts [9-11].
Conducting risk assessments for MPs is essential to identify potential risks to ecosystems and humans, and it requires integrating exposure and hazard assessments in a complex evaluation process due to MPs’ diverse characteristics [12, 13].
MPs pose significant ecological risks, potentially harming ecosystems and human health when they enter food chains. This highlights the critical need to assess their impacts using important indices, including the pollution load index (PLI), risk index (RI), and hazard index (HI) [14, 15]. Human health risks associated with MPs encompass direct and indirect effects, including physical stressors and the transport of environmental contaminants, impacting bodily systems and leading to adverse health outcomes [14, 15]. The ingestion of MPs by organisms and their transfer through the food chain pose serious risks to both human health and ecological systems. This concern has led organizations, such as the World Health Organization (WHO), to stress the urgent need for increased research to better understand the implications of MPs on human health [16].
Given the widespread presence of MPs in various food items consumed by humans, including milk, seafood, honey, table salt, tap water, and bottled water, numerous researchers have increasingly focused on this critical issue. While the intake of certain food products, such as seafood and honey, can be reduced or avoided, exposure to MPs in drinking water is inevitable [17]. The presence of MPs in various food items consumed by humans underscores the urgency of researching their impact on human health, especially through unavoidable exposures, such as drinking water [12, 18].
In aquatic environments, the prevalence of MPs is attributed to inadequate wastewater treatment, transportation activities, and improper waste disposal practices, with industrial and domestic wastewater acting as prominent pathways for MP introduction [19, 20]. Despite global efforts to detect MPs, identifying their diverse sources remains a critical challenge [21]. The growing popularity of exploratory analysis techniques for statistical data evaluation offers opportunities to visualize correlations among multiple variables in a single graph, aiding in the identification of pertinent information for sample classification [22, 23]. Multivariate statistical methods are essential for analyzing large environmental datasets, enabling researchers to accurately identify sources of pollution [24, 25].
While numerous studies have examined MP abundance, distribution, and source identification in aquatic systems using multivariate statistical approaches (e.g. principal component analysis [PCA], cluster analysis [CA], and hierarchical cluster analysis [HCA]), few have placed primary emphasis on comprehensive risk assessment. This systematic review addresses this critical gap by focusing predominantly on ecological and human health risk evaluation, comparatively analyzing indices, such as PLI, polymer HI (PHI), and potential ecological RI (PERI) for ecological risks, alongside estimated daily intake (EDI), target hazard quotient (THQ), and MP pollution load index (MPLI) for health implications, across diverse global aquatic environments. By synthesizing risk data from selected recent studies and integrating insights from source apportionment, this work provides a novel, risk-centered framework to prioritize high-risk areas, inform targeted mitigation strategies, and offer actionable guidance for policymakers and researchers amid escalating MP pollution.
Materials and methods
Search strategy
This systematic review was conducted in accordance with the preferred reporting items for systematic reviews and meta-analyses (PRISMA) 2020 guidelines [26]. A comprehensive literature search was performed in December 2025 across four major scientific databases (PubMed, Scopus, Web of Science, and ScienceDirect) to identify relevant peer-reviewed studies addressing MP source identification and risk assessment in aquatic environments. Grey literature was excluded from the review. The protocol for this systematic review was not registered.
This systematic review was based on the following search strategy and medical subject headings (MeSH), applied to Title and Abstract searches. For MP source identification, the search string was: (“microplastic*” OR “micro-plastic*” OR “micro plastic*”) AND (“source*” OR “origin*” OR “pathway*” OR “emission*” OR “input*” OR “release*” OR “anthropogenic source*” OR “pollution source*” OR “source apportionment” OR “source tracing”) AND (“multivariate” OR “PCA” OR “principal component analysis” OR “CA” OR “cluster analysis” OR “HCA” OR “hierarchical cluster analysis” OR “machine learning” OR “random forest” OR “decision tree”) AND (“aquatic” OR “water*” OR “freshwater” OR “marine” OR “river*” OR “lake*” OR “estuar*” OR “ocean*” OR “sediment*” OR “surface water” OR “wastewater”).
For MP risk assessment, the search string was: (“microplastic*” OR “micro-plastic*” OR “micro plastic*”) AND (“risk assessment*” OR “ecological risk*” OR “health risk*” OR “hazard assessment*” OR “toxicity*” OR “exposure assessment*” OR “polymer hazard*” OR “pollution load index” OR “PLI” OR “PHI” OR “polymer hazard index” OR “PERI” OR “potential ecological risk index” OR “EDI” OR “estimated daily intake” OR “THQ” OR “target hazard quotient” OR “MPLI”) AND (“aquatic” OR “water*” OR “freshwater” OR “marine” OR “river*” OR “lake*” OR “estuar*” OR “ocean*” OR “sediment*” OR “surface water” OR “wastewater” OR “aquaculture”).
All searches were restricted to English-language articles published between January 2015 and December 2025. This time frame was selected to encompass both foundational studies and recent advances in multivariate statistical analyses and risk assessment frameworks related to MPs.
Selection parameters and data extraction
The inclusion criteria comprised original, peer-reviewed research articles with full-text availability that investigated MP sources and/or ecological or human health risks in aquatic environments (e.g. rivers, lakes, reservoirs, estuaries, coastal and marine systems). Eligible source-identification studies were required to apply multivariate statistical approaches (such as PCA, CA, HCA, multidimensional scaling [MDS], multiple correspondence analysis [MCA], or hybrid and machine-learning–assisted techniques), while risk-focused studies needed to report quantitative indices including PLI, PHI, PERI, RI, THQ, or HI.
The exclusion criteria included review articles, letters to the editor, books, book chapters, non-aquatic studies, laboratory-only or modeling-only studies lacking environmental data.
Three authors independently conducted all screening stages to minimize selection bias. Any disagreements regarding study inclusion were resolved through discussion and, where necessary, consensus among the reviewers. For each included study, data were systematically extracted on study location, sampling (aquatic matrix), polymer types, applied statistical methods for source apportionment, and reported ecological and/or human health risk metrics. The reference lists of the included papers were also hand-searched to identify any additional relevant studies. The overall study selection process is illustrated in the PRISMA flowchart (
Figure 1).
Risk of bias assessment
The methodological quality and risk of bias of the included studies were systematically evaluated using the National Heart, Lung, and Blood Institute (NHLBI) quality assessment tool for observational cohort and cross-sectional studies [27]. This tool was selected because most eligible studies employed observational designs using environmental sampling and quantitative analyses relevant to microplastic source identification and ecological or human health risk assessment.
The NHLBI checklist comprises 14 items addressing key domains of internal validity, including clarity of the research question, definition and representativeness of the study population or sampling framework, consistency of eligibility criteria, validity and reliability of exposure and outcome measurements, adequacy of the study timeframe, appropriateness of statistical analyses, and consideration of potential confounding factors. Each item was independently rated as yes, no, cannot determine, not reported, or not applicable in accordance with the NHLBI guidance. An overall methodological quality rating was then assigned to each study based on the collective assessment of internal validity and susceptibility to selection, information, measurement, and confounding bias. Studies were categorized as good (low risk of bias), fair (moderate risk of bias), or poor (high risk of bias). Quality assessment was conducted independently by two reviewers, and any discrepancies were resolved through discussion and consensus. To account for methodological heterogeneity, source identification and risk assessment studies were evaluated separately. Detailed quality assessment outcomes are provided in Supplementary
Tables S1 and
S2.
Table S1 includes 22 ‘fair’ and 2 ‘good’ quality assessments, while
Table S2 contains 14 ‘fair’ and 5 ‘good’ quality assessments. Overall, most of the included studies were rated fair to good, indicating acceptable methodological rigor for synthesis and interpretation.
Results
Statistical methods for identifying MPs sources
Sources and contributions to MP pollution vary significantly across geographic and ecological contexts, underscoring the importance of localized investigations to inform targeted management strategies for aquatic systems [25, 29]. Multivariate statistical methods are instrumental in processing large environmental datasets, converting intricate information into actionable insights by revealing patterns in MP distribution, composition, and origins. Conventional exploratory techniques, such as PCA, CA, and HCA, are widely employed to identify dominant groupings and associations among MP characteristics [30-32]. To move beyond qualitative source interpretation, these techniques are frequently integrated with regression- and receptor-based frameworks, including PCA–multiple linear regression analysis (MLRA), absolute principal component scores–multiple linear regression (APCS-MLR), and positive matrix factorization (PMF), which enable quantitative estimation of relative source contributions [31-33].
In addition, hybrid statistical–machine learning approaches, such as PCA combined with random forest (PCA-RF), along with K-means clustering, decision trees, Naïve Bayes classifiers, MDS, and MCA, facilitate the identification of pollution sources while enhancing water quality assessments [23, 28, 25]. These approaches not only delineate the primary factors influencing aquatic systems but also support integrated risk evaluations, enabling policymakers to prioritize interventions based on empirical evidence.
Table 1 presents a systematic overview of representative studies applying these multivariate statistical techniques, such as PCA, CA, and HCA for MP source identification, highlighting methodological approaches, investigated aquatic environments, and analytical tools employed.

Statistical analyses were predominantly performed using SPSS software, version 26 reflecting its widespread applicability in environmental data interpretation and source apportionment.
Geographical trends in source contributions
Figure 2 shows the regional distribution of studies on MP source identification, highlighting research hotspots and showing that the majority of studies were conducted in Asia, particularly China.
By contrast, European and North American environments demonstrate more localized point-source contributions.
Environmental risk assessment
Ecological risk assessment
To tackle MP pollution, it is crucial to implement targeted control measures in areas at high ecological risk of MP contamination [53]. However, there is currently no systematic and standardized model available to assess the potential ecological risk that is specifically suitable for MP pollution [56]. Consequently, researchers have developed a comprehensive risk assessment model applicable to MP pollution and other pollutants. Three models (PLI, PHI, and PERI) have been used to assess the ecological risk of MPs. This model aimed to improve our understanding and evaluation of the ecological risks associated with MP contamination [57].
PLI of MPs
To evaluate the extent of MP pollution in aquatic settings, we employed an integrated PLI calculation method developed by Qiu et al. (2023) [58]. The assessment model is as follows:

where ci is the MP abundance at each sample site, c0 is the background value of MPs, n is the number of stations. PLI is divided into four levels: I (low, <10), II (medium, 10≤ PLI <20), III (high, 20≤ PLI <30), and IV (extremely high, PLI ≥30) [53, 58-67].
PHI of MPs
The occurrence and effects of MPs in aquatic ecosystems are affected by both their concentration and the specific polymer types they contain, which ultimately influence their toxicity. A risk grading model developed by Lithner et al. (2011) [68] categorizes polymers and calculates their risk scores based on monomer hazard levels. The hazard score values for polypropylene (PP), polyethylene (PE), polyacrylonitrile (PAN), polyethylene terephthalate (PET), polystyrene (PS), polyamide (PA), and polyvinyl chloride (PVC) are 1, 11, 11521, 4, 30, 50, and 5001, respectively [68].
The PHI is then calculated using a specific formula:

where Pi represents the proportion of polymer in each sample, and Sn denotes the risk score associated with the type of MP polymer. The PHI is categorized into four levels: I (low, PHI <10), II (medium, 10≤ PHI <100), III (high, 100≤ PHI <1000), and IV (extremely high, PHI ≥1000) [53, 58-68].
PERI
PERI was utilized to evaluate the extent of MP contamination. The following equations were employed to calculate the PERI:

Where CFi is the enrichment coefficient of MPs in samples; Ci is the abundance of MPs in samples; Cin represents the abundance of MPs in the samples; Tri is the toxicity coefficient of polymer; Pn refers to the abundance of each polymer in the samples; and Sn is the hazard score associated with each polymer. The toxicity coefficient is calculated as the sum of the percentage of specific polymers in the total samples Pn/ Ci multiplied by the toxicity hazard score of the polymers Sn. PERI is divided into five levels: I (low potential risk, PERI <150), II (Medium potential risk, 150≤ PERI <300), III (Considerable potential risk, 300≤ PERI <600), IV (High potential risk, 600≤ PERI <1200), and V (very high potential risk, PERI ≥1200) [53, 58, 59, 65-67, 69].
Health risk assessment
Evaluating exposure and characterizing hazards are crucial components of human health risk assessment (HHRA) [70]. Human health risk assessment (HHRA) is a systematic approach for evaluating the nature and probability of adverse health effects in individuals exposed to environmental hazards. This process considers both current and potential future exposures to pollutants in various environmental media [71].
EDI of MPs
Basaran et al. (2023) [69] utilized the following equation within the deterministic model to estimate an individual’s intake of MPs from consuming branded milk.

EDI is the EDI (particles/mL bw/day), Ma is the amount consumed (mL/day), Mc is the concentration of MPs (particles/mL), and Bw is the body weight (kg). The average body weight is 61.4 kg, and the consumption rate is 130 g/day per individual [69, 70].
THQ
THQ is a ratio comparing an individual’s exposure to a hazardous substance with the level of exposure at which adverse health effects are likely to occur [72]. The following formula was utilized to calculate the THQ of MPs:

Where EF denotes the frequency of exposure (156 days per year, based on three fish meals weekly), ED indicates the duration of exposure (70 years, aligned with a typical lifespan), FIR refers to food consumption (130 g/day), C represents the concentration of HMs in fish muscle (µg/g wet weight), WAB stands for average weight, TA signifies the total exposure time to non-carcinogens (365 days/year multiplied by ED), and RFD specifies the reference dose for each heavy metal. According to this model, if THQ is less than 1, the health effects from oral intake of heavy metals are minimal. Conversely, if THQ exceeds 1, there may be a notable non-carcinogenic risk associated with this exposure pathway [59, 73-75].
MPLI
The MP contamination factor (MCF) and the MPLI in milk were determined using the following equations:

MPi is the number of MPs in the sample, n is the total number of samples, and MPb is the minimum reported average MP concentration (1.68 particles/kg) in processed foods [69].
Table 2 presents the standard grid for this analysis [69, 76].

Ecological risk indices, including PLI, PHI, and PERI, showed a wide range of values across different study areas, reflecting variations in polymer composition and contamination levels.
Table 3 presents representative results from selected studies.
Discussion
Dominant sources identified by multivariate methods
A synthesis of the reviewed literature across diverse aquatic environments indicates that multivariate statistical analyses predominantly attribute MP pollution to anthropogenic sources, with urban wastewater, industrial effluents, agricultural runoff and maritime activities consistently identified as major contributors [36, 39, 49]. PCA-based investigations frequently resolve agricultural and forestry-related inputs as dominant sources in rural river basins. In the Jinjiang River Basin, PE and PP exhibited strong loadings on agricultural components, whereas PVC and PET were associated with household sewage and fishing activities [30]. Similar source structures were reported in Taiwanese rivers, where PCA linked MP variability to vehicle exhaust and petrogenic inputs [35], and in the Huixian Wetland, where effluent discharge emerged as the primary urban driver alongside suburban agricultural runoff [10].
Beyond PCA-based interpretations, quantitative receptor models have provided higher-resolution source apportionment. For instance, PMF-based analysis in southwestern Nigerian rivers resolved multiple contributing factors and identified surface runoff and domestic/industrial wastewater as the dominant MP sources in both surface water and sediments, while packaging-related plastics showed moderate contributions and atmospheric deposition played a minor role, highlighting the predominance of land-based anthropogenic activities such as fishing and inadequate sanitation systems [52]. Similarly, APCS-MLR applied to sediments of the Dan River in China quantified two major anthropogenic sources, with mixed fishery and agricultural activities contributing the majority of MPs, followed by wastewater discharge associated with laundry and personal care products [33].
Cluster-based approaches, including CA and HCA, support these findings. In Laizhou Bay, such analyses identified riverine inputs, plastic recycling activities, and aquaculture as key sources, with PET dominating surface waters [32].
Recent studies extend these interpretations while providing additional resolution in source apportionment [86]. In the Xiangjiang River, multivariate analysis linked specific MP characteristics (e.g. transparent and black colors, 50–100 μm size fractions) to urban wastewater treatment plants and revealed positive associations with microorganisms, such as Proteobacteria, suggesting microbial-mediated transport pathways [45, 53]. In marine environments, the integration of PCA with random forest models achieved high accuracy in identifying commercial polymers (PE, PS, and PP) and indicated contributions from shipping activities and consumer-derived waste [50]. In the Brantas River watershed, PCA and CA confirmed recreational littering as the dominant source at high-abundance sites, with lower contributions from agricultural fields [38]. Similarly, in the Bay of Bengal, PCA and CA grouped PE and PP with anthropogenic factors and indicated medium-to-high ecological risks associated with fishing nets and urban runoff [39]. In Shenzhen’s river systems, the combination of PCA and random forest analysis revealed downstream industrial dominance and demonstrated that correlations based solely on physicochemical parameters tended to underestimate risks without machine-learning integration [51], whereas urban non-point sources are more influential in coastal systems [40, 41]. This variability underscores the need for polymer-specific and hybrid statistical frameworks to disentangle overlapping source pathways.
Geographical trends in source contributions
MP pollution exhibits pronounced geographical variability, influenced by regional differences in population density, industrialization, land-use patterns, and environmental management practices [38, 39]. Studies compiled in
Table 1, indicate that Asian riverine and coastal systems consistently report the highest MP abundance and the most diverse sources. In China, river basins, estuaries, and wetlands, including the Jinjiang River, Laizhou Bay, Dan River, Bohai and Yellow Seas, and Shenzhen rivers, show predominant contributions from land-based anthropogenic activities such as agriculture, aquaculture, urban wastewater, and industrial discharges [30, 32, 33]. Southeast and South Asian regions, including Indonesia and Bangladesh, exhibit elevated MP concentrations in densely populated coastal areas and recreational zones, where urban runoff, fishing, and shipping activities are key drivers [38, 39]. Based on
Figure 2, in Adventfjorden (Norway) and urban wastewater networks in Denmark and Belgium, MPs are primarily associated with effluent discharges, household and industrial wastewater, and combined sewer overflows [34, 48, 49]. Similarly, in Lake Ontario (Canada), MPs are linked to wastewater treatment plant effluents and hydrologically mediated redistribution, highlighting the influence of regulated infrastructure on MP distribution [43].
Evidence from Africa, including southwestern Nigerian rivers, indicates that MP sources are dominated by surface runoff and domestic and industrial wastewater, whereas packaging materials and atmospheric deposition contribute minimally [52].
These patterns indicate that local anthropogenic pressures and regional environmental management largely govern the geographical distribution of MP pollution. Densely populated urban and intensive agricultural areas tend to exhibit diffuse, multi-source contamination. In contrast, well-regulated systems show more localized, point-source contributions, highlighting that broad global generalizations may fail to capture critical site-specific drivers.
Performance and comparison of different statistical techniques
A comparison of statistical techniques employed across the reviewed studies demonstrates that multivariate methods play a central role in resolving MP source patterns, although their effectiveness varies depending on data structure, sample size, and environmental context. PCA remains the most widely applied approach, owing to its capacity to explain a substantial proportion of variance in MP datasets and to identify dominant source-related components, with reported cumulative variance often exceeding 65% [35, 40, 54]. As such, PCA is particularly effective for initial pattern recognition and dimensionality reduction in complex environmental matrices.
However, PCA alone is often insufficient for precise source apportionment, prompting its integration with complementary techniques. Regression-based hybrids such as PCA–MLRA improve quantitative source attribution in marine environments by linking principal components to specific anthropogenic activities [37]. Cluster-based methods, including CA and HCA, are commonly applied to group sampling sites according to similarity in MP composition and abundance, thereby revealing spatial gradients associated with riverine inputs, aquaculture zones, and urban coastal areas [10, 32, 431].
Non-linear ordination methods such as MDS and MCA provide additional insight into complex source structures that may not be adequately captured by linear techniques, particularly in heterogeneous systems influenced by mixed urban, industrial, and agricultural activities [42, 44]. Classification-based approaches, including decision trees and Naïve Bayes classifiers, offer advantages in interpretability and material discrimination, although their performance may be constrained by limited training datasets or highly mixed environmental samples [46-48].
Recent advances have highlighted the superior performance of hybrid statistical–machine learning frameworks. The integration of PCA with random forest models substantially enhances source discrimination and classification accuracy, particularly in marine environments characterized by overlapping spectral and compositional features [50]. Similar hybrid approaches reveal that reliance on traditional correlation analyses alone may underestimate industrial and downstream risks, underscoring the value of machine learning in resolving non-linear relationships between MP properties and environmental drivers [55]. In sediment-focused applications, receptor and hybrid regression models further improve quantitative source resolution. For instance, PMF demonstrated strong performance in resolving multiple mixed anthropogenic sources in surface waters and sediments of southwestern Nigerian rivers, effectively quantifying the dominance of runoff and wastewater-related inputs in complex datasets [52]. Similarly, the combined PCA–APCS–MLR framework applied to Dan River sediments translated statistically derived components into quantitatively meaningful anthropogenic source contributions, illustrating the enhanced interpretability and apportionment capability of regression-based hybrid approaches [33]. In general, receptor models (APCS–MLR) offer more robust quantitative apportionment than PCA alone, while hybrid PCA–RF approaches improve classification accuracy in complex marine datasets.
Despite their demonstrated strengths, all techniques remain subject to methodological limitations, including sensitivity to data pre-processing, normalization procedures, and sampling design. Moreover, statistical associations do not imply causation, necessitating cautious interpretation of inferred contributions from sources.
Comparative analysis of ecological and health risks from MPs
Risk assessment frameworks applied in this review provide critical insights into the multifaceted threats posed by MP pollution in aquatic ecosystems. Although MP-specific standardized models remain limited [56], adapting established indices, such as PLI, PHI, and PERI for ecological evaluation, alongside EDI, THQ, and MPLI for human health assessment, offers a pragmatic approach to quantifying risks.
A comparative analysis across studies, as summarized in
Table 3, reveals pronounced spatial and source-driven heterogeneity in ecological risks. Low PLI values may obscure substantial hazards when highly toxic polymers dominate, as observed in the Pearl River Estuary, where moderate contamination (PLI 1.0–1.43) coincides with extreme PHI (up to 1375) arising from polycarbonate in urban-industrial runoff [61]. Similar patterns are evident in port-dominated systems such as the Persian Gulf, where high PLI (>30) and PHI (>1000) indicate significant ecological vulnerability [77]. In tropical regions, including the Northern Bay of Bengal, land-based sources contribute to elevated PERI values through PVC and mixed polymers, driving ecological risks into categories III–IV [39]. Seasonal variations further modulate risks, exemplified by Thailand’s Ubolratana Reservoir, where dry-season retention amplifies PHI (64–1510) under heightened tourism pressures [60], a pattern corroborated in sedimentary hotspots such as Lake Aba Samuel (PERI >2000) and Yahekou Reservoir (hazardous HI in deeper strata) [79, 80]. Collectively, these findings underscore PHI and PERI as more sensitive indicators of polymer toxicity than abundance-focused PLI, particularly within sediment compartments [82].
Human health implications, although less frequently quantified in MP-specific studies, emerge through EDI and THQ assessments, highlighting ingestion as a primary exposure pathway in urbanized contexts [70, 71]. Notably, children exhibit higher EDI than adults, such as up to 13.03 MP/kg/bw/year via South American bottled water, suggesting heightened susceptibility to oxidative stress and inflammatory responses [82]. Food-based exposure, such as branded milk, demonstrates a moderate MPLI (2.8) with polymer-specific hazards (PU up to 198), resulting in a lifetime exposure of approximately 5289 particles over 70 years [69]. Although THQ values <1 suggest limited acute non-carcinogenic effects, synergistic interactions between MPs and heavy metals in fish from Indonesian streams elevate both total THQ and cancer risks through the Trojan horse mechanism [59]. Comparably low-to-moderate THQ in Thai reservoirs further indicates chronic underestimation of risk due to food chain transfer [60].
The convergence of ecological and human health risks positions MPs as effective vectors for co-contaminants, thereby intensifying threats in polluted systems [59]. Methodological inconsistencies, including variable sampling approaches and lack of uniform background referencing, may bias risk indices [82]. Consequently, the development of MP-tailored frameworks that incorporate temporal variability, synergistic effects, and compartment-specific dynamics is imperative to inform targeted mitigation strategies, including source reduction and advanced treatment interventions [53, 69]. This synthesis advocates proactive, integrated approaches to limit the escalating impacts of MP on aquatic ecosystem integrity and human well-being.
Limitations and challenges
This review is limited to peer reviewed studies published in English, which may have resulted in the exclusion of relevant regional research available in Persian or other local languages. This language constraint may introduce a geographic bias in the synthesis of available evidence. Multivariate statistical methods, such as PCA and CA, are valuable for hypothesis generation in MP source identification but are constrained by their inability to establish causality and often rely on correlational inferences that may be confounded by environmental collinearity and temporal variability [24]. The challenge of overlapping sources, e.g. distinguishing agricultural runoff from industrial discharges, remains unresolved in many applications, as these methods struggle with non-linear interactions without hybrid integrations, such as PCA-RF [50]. These limitations imply a cautious interpretation of the results, particularly in dynamic aquatic systems where unaccounted-for variability could misrepresent pollution pathways, potentially undermining mitigation strategies.
In risk assessment frameworks, the absence of standardized, MP-specific toxicity reference values (RfDs) poses a major hurdle; current adaptations from heavy metal models fail to capture unique MP hazards, such as additive leaching or fragmentation, resulting in questionable THQ and HI estimates [83]. Significant uncertainties in exposure assessment, including bioavailability influenced by weathering and chemical leaching, further complicate evaluations, often underestimating chronic risks from nano-MPs [13]. The ongoing “Trojan horse” debate, which highlights the gap between the risks arising from polymers and those from adsorbed contaminants, such as PAHs, suggests that the latter are amplified in food chains [8]. Such methodological shortcomings can lead to inconsistent risk categorizations across studies, affecting policy decisions; for example, overreliance on borrowed RfDs may downplay health impacts in high-exposure regions. Implications for future research include developing MP-tailored RfDs through toxicological experiments and incorporating advanced modeling to address exposure uncertainties, thereby enhancing the reliability of indices, such as PERI and PHI.
Conclusion
MPs represent a significant and growing threat to aquatic environments, with substantial ecological and potential human health implications. This review demonstrates that multivariate statistical approaches, particularly PCA, CA, and HCA, are effective tools for identifying MP sources and elucidating their strong associations with anthropogenic activities, such as urbanization, industrial discharge, and agricultural runoff.
The findings further highlight that integrating source identification with ecological and human health risk indices, such as PLI, PHI, PERI, and THQ, provides a comprehensive framework for assessing MP distribution and associated risks across diverse environmental settings. Accurate risk evaluation, however, depends on a clear understanding of MP bioaccessibility and bioavailability, which directly influence exposure and toxicity outcomes. Overall, advancing MP research requires developing polymer-specific toxicity data and adopting standardized methodologies for sampling, identification, and analysis to improve comparability and support more robust, reliable risk assessment and management strategies.
Ethical Considerations
Compliance with ethical guidelines
This study was approved by the Research Ethics Committee of Mazandaran University of Medical Sciences, Sari, Iran (Code: IR.MAZUMS.REC.1402.224).
Funding
This study was financially supported by Mazandaran University of Medical Sciences, Sari, Iran.
Authors contributions
Conceptualization: Atefeh Jabari; Methodology: Yalda Hashempour and Afsaneh Fendereski; Data curation and formal analysis: Kosar Kouhi; Writing the original draft: Kosar Kouhi and Atefeh Jabari; Review, editing and supervision: Yalda Hashempour.
Conflict of interest
The authors declared no conflict of interest.
Acknowledgements
The authors thank the Student Research Committee at Mazandaran University of Medical Sciences, Sari, Iran, for their generous financial assistance.