The quality of the surface water along the Kunhar River in Mansehra district, KPK Pakistan was assessed for ten water samples stations. A variety of parameters indicating water quality including physicochemical parameters, nutrients, heavy metals, and antibiotic residues were measured for both the rainy and dry seasons, the two main tropical seasons in Mansehra, using the standard methods. Kunhar River, one of the local drinking water sources, was studied to assess the heavy metal content, health risk assessment, and its suitability for human consumption. Health risk assessment for all the stations indicated that there is no particularly dangerous single heavy metal, but their cumulative effect is indicated by the hazard index. Concentrations of metals in water have a trend of decreasing in water, indicating that the accumulation of metals can affect the water chemistry of aquatic systems due to any physical or chemical process in the system. Greater consideration should be given to the variety of metals in relationship to multivariate analyses, suggesting that the industrial and residential activities were more important contributors to the pollution of the Kunhar River than the agricultural activities in Mansehra district. Such metals pose risks to aquatic life.

  • In the Kunhar River, seasonal fluctuations in the levels of nutrients, toxic substances, as well as certain antibiotics were identified.

  • The study revealed that domestic activities caused significantly more environmental damage to the Kunhar River than agricultural practices.

  • High-pollution sites of the Kunhar River have been affected by human settlement.

Various aspects of development depend significantly on rivers because they are an important source of freshwater (Link et al. 2016; Gleeson & Richter 2018). Grazing by livestock is the main anthropogenic disturbance of native grasslands, and the main factor affecting the accumulation of potentially toxic metals (PTMs) by grazing animals is the presence of the metals (Li et al. 2021). Human activities can alter dissolved organic matter (DOM) in lakes through both direct (i.e., exporting DOM of anthropogenic sources) and indirect effects (i.e., enhancing the autochthonous production of DOM via nutrient loading). Distinguishing between the direct and indirect effects of anthropogenic activities remains highly challenging due to the interdependence of associated environmental variables (Liu et al. 2022). Drinking water is typically acquired from surface waters such as surface water bodies (lakes, rivers, streams, etc.) and groundwater (springs, underground streams, etc.) (Asfaw et al. 2018; Coffer et al. 2021; Ali & Muhammad 2022; Muhammad & Usman 2022). Also, floods bring about immediate material damage caused by the direct effects of high water levels and accelerated flows (Soomro et al. 2021b). Agricultural overflow and effluents from settlements and factories make their way to rivers in a watershed, where pollutants are either diluted or transported downstream. Plants and factories, the hallmarks of industrialization, are essential to a country's economic growth (Yao et al. 2015). To protect the ecological environment, biodiversity, and natural resources, one must protect the integrity and authenticity of typical and unique ecosystems from damage, and leave a valuable natural heritage for future generations (Pan & Chen 2021; Zhao et al. 2022). A variety of physical and geochemical proxies for weathering intensity have been explored to understand loess provenances, environmental, and climatic variations in the regions of deposition (He et al. 2021). Physical and biogeochemical heterogeneity dramatically impacts fluid flow and reactive solute transport behaviors in geological formations across scales (Zhang et al. 2021). Rivers susceptible to residential wastewater discharge and other human influences and drivers of different factors were identified by utilizing a geographical detector, especially the influences of human activities in the region (Amin et al. 2021; Quan et al. 2022). The residential environment has been suggested as an important source of environmental chemicals to human exposure. House dust has been broadly employed to investigate the occurrences of environmental chemicals in an indoor environment (Yang et al. 2022). River ecosystems across the globe have already suffered from the impacts of climate change, and these impacts are only expected to worsen (Soomro et al. 2022). However, the raw materials for the production of wastes they release are harmful to the environment and cause water and soil contamination at multiple levels (López-Delgado et al. 2020; Silva et al. 2022; Verma & Sharma 2022). As more countries pursue plans for geologic disposal of accumulated high-level nuclear waste, quantifying the fate and transport of radionuclides in complex fractured rocks has become an important issue ((Zhang et al. 2022). Excessive reactive N in terrestrial and marine ecosystems from large applications of fertilizer has impacted the balance of the global N cycle and contributed to numerous eco-environmental problems (Lin et al. 2021). Precipitation is a major determinant of vegetation growth and the impact of precipitation variability is more pronounced in ecosystems where sensitive vegetation is apparent in the Kunhar River basin (Soomro et al. 2021a). Deficiency of high efficiency and cost-effective treatment technology prevents industrial effluents from obtaining safely treated and organizational control and monitoring inadequacy. Metals and alkaline earth metals with concentrations of more than 5 g cm3 were generally referred to as ‘heavy metals,’ a term that includes any metallic component with a high density (Liu et al. 2019; Tyopine et al. 2020). Human activities have accelerated global habitat loss and fragmentation, transforming continuous habitats into smaller, isolated habitat fragments surrounded by disturbed landcover (Yan et al. 2022). It's clear that these elements are significant for critical metabolic activities throughout all species of life, but at concentrations greater than any of those found naturally, they can be toxic (Brudvig et al. 2007; Mirrione et al. 2014; Muhammad & Ahmad 2020). Among the most key problems facing modern humanity is the pollution of earth and water by toxic substances connected to characteristics of confidentiality, stability, and permanence (Ernst 2012; Briffa et al. 2020).

Even while toxic substances could be detected in minerals, their occurrence is generally the direct result of human activities (Bloise et al. 2016; Rebello et al. 2021). The quality of drinking water may deteriorate when heavy metals affect the supply of water at various depths, especially surface water, reservoirs, and subsurface (Khatri & Tyagi 2015; Bhateria & Jain 2016; Gyamfi et al. 2019). Measuring the amounts of heavy metals in freshwater is one way to examine water pollution in the Kunhar River (Muhammad et al. 2021). Health complications and the incapability of people to get along with each other are two visible signs of the widespread poverty that have been recorded (Braithwaite & Mont 2009). The locals consume the water for all kinds of domestic purposes, including drinking, cleaning, bathing, and washing. The purpose of this research was to investigate whether or not it was safe to drink water from the Kunhar River, as well as to quantify the number of heavy metals in it.

The Kunhar River watershed travels across the southern Himalayan plateau from 34.2–35.1°N to 73.3–74.1°E; the area of the catchment is 2,632 km2. It originates in the Kaghan Valley and afterward travels via the Lulusar Lake to reach us. It passes through Batakundi, Naran, Kaghan, Mahandri, Jared, Paras, Kawai, Balakot, Kassi Tarana, and Bisian as shown in Figure 1. It joins the Jhelum River near Rara before continuing to its final destination. Since our sample site is at a higher elevation, it is necessary to highlight the present climate trend, which is extremely dynamic in respect to seasonal fluctuations, since this is one of the most distinguishing characteristics of any study area. Catchments at high elevations throughout Pakistan are all under increasing stress caused by both contaminants and climatic condition changes, and studies in such regions around Pakistan are challenging, and are mentioned in Figure 2. A total of 56 water samples were collected along the river from March to September (summer, wet season) and November 2021 to February 2022 (winter, dry season).
Figure 1

Location of sampling sites (Batakundi, Naran, Kaghan, Mahandri, Jared, Paras, Kawai, Balakot, Kassi, Tarana, and Bisian).

Figure 1

Location of sampling sites (Batakundi, Naran, Kaghan, Mahandri, Jared, Paras, Kawai, Balakot, Kassi, Tarana, and Bisian).

Close modal
Figure 2

Variation of seasonal discharge (m3/s), temperature (°C), precipitation (mm/day), and evaporation for a period in the Kunhar River (2010–2021). Source: Pakistan Meteorological Department.

Figure 2

Variation of seasonal discharge (m3/s), temperature (°C), precipitation (mm/day), and evaporation for a period in the Kunhar River (2010–2021). Source: Pakistan Meteorological Department.

Close modal

Analytical method

Using an HTC-2 hygrometer, the temperature and relative humidity were measured instantly. The pH and EC of Kunhar River water samples were measured using a pH meter (HI 2211, Hanna, USA) (Kandhro et al. 2016) with a portable electrical conductivity tester (TET) (HI-98302, Hanna, USA) (Koukina et al. 2021). The total dissolved solids (TDS), turbidity, total suspended solids (TSS), total hardness (TH), biological oxygen demand (BOD), phosphate, fluoride, nitrate, sulfate, and chloride ratio in the river samples were determined and tested according to a standard procedure (Gorde & Jadhav 2013). Trace metals in water were analyzed using a 0.45-μm syringe filter and 1% HNO3 acidification (Finch et al. 2018). A suspension of (100 g/l) sediments was used to perform pH and EC measurements. Losses on ignition (LOI) techniques are used by paleolimnologists, geochemists, and sedimentologists to figure out the configuration of sediment and rock samples (Terasmaa et al. 2013). It is a cheap and low-tech method. In the current research, oven drying and the LOI technique were used to evaluate the sediment's moisture levels and estimated its carbon content. To estimate how much cation exchange capacity (CEC) sediments have quantitatively, barium is used to get rid of soil-exchangeable cations. To get rid of extra barium, four deionized water rinses are used. A known amount of calcium is then exchanged for barium, and the amount of calcium in the solution without barium is found. The CEC is based on the difference between how much calcium was added and how much calcium was found in the solution that came out of it (Johnson et al. 2014). The method can detect about 2.0 mEq/100 g of water. The total metal content of the sediment samples was determined by the digestion of the sediment samples. Additionally, a value control of sandy loamy mud was employed to ensure the technique's effectiveness. Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-OES) is a well-known analytical technique that has a wide linear dynamic range, the limit of detection, high precision, and accuracy (Biata et al. 2018). It can be used to quickly and easily determine the major, minor, and trace elements in river water and also widely used to determine heavy metal contents in water. Heavy metals like Zn, As, Cu, Cd, Hg, Cr, Fe, Pb, and Se were determined in the remainders from acid breakdown and acidic water samples using the ICP-OES method (Alansi et al. 2021). The methodology of the entire experiment is shown in Figure 3.
Figure 3

Flow diagram of the research methodology.

Figure 3

Flow diagram of the research methodology.

Close modal

Health risk assessment

The metals determined to be at harmful levels were examined regarding their potential effect on human life (Fe, Mn, Zn, Cd, and Cr). The analysis followed a healthcare risk-based approach of non-carcinogenic living beings, defined by Qiu et al. (2021). The toxic substances in the Kunhar River water were determined using the following formula (Muhammad et al. 2011).

(1)
where CDI is the daily intake of heavy metals (mg/l) to which consumers might also be highly susceptible, CW (mg/l) is probably the presence of heavy metals in the riverside water, IR is the average daily intake, EF is how often subjected to said absorption, ED is the frequency of absorption, BW is the weight of the body, and AT is the mean time.

Hazard quotient

We were able to identify the hazard quotient (HQ) for non-carcinogenic harm by USEPA (1999).
(2)
where chronic daily intake is the daily intake of heavy metals (mg/l)) where people's intake might also be highly susceptible, and is the daily starting point, or the minimum effective dose at which an individual can be treated for an extended period without adverse consequences. Reference dose (RfD) levels for ingested effects of heavy metals are shown in Figure 5: If HQ > 1, it was an indication of possibly hazardous quasi-effects which is a risk, while HQ < 1 represents a satisfactory amount.

Hazard index

Due to the presence of numerous toxic substances, their interactions are taken into account. It is typically assumed that the risk of the harmful metals caused by toxic metals co-existing within the same vicinity is cumulative. The HQs, when combined collectively, constitute the hazard index (HI), a measurement of aggregate toxic risk (Paustenbach 2015).
(3)
where HI is the ratio of harmful hazards, and n is the number of metals taken into account. It is believed that the non-carcinogenic adverse effect due to this environmental mechanism or toxin is minimal.

Analysis of antibiotic residues

Intake of antibiotics from the environment by drinking water and food may disturb the microbiome, especially the gut microbiota in the human body. A solid-phase extraction accompanied by an excessive performing separation system that used a mixture of liquid chromatography and mass spectrometry in tandem was performed to determine the chemicals in the surface water samples (HPLC-MS/MS) (Hong et al. 2021). Ethylenediaminetetraacetic acid (EDTA; 100 mg/200 ml) was combined with such a water sample (Madadrang et al. 2012). Two different solid-phase extraction (SPE) cylinders (3 ml, 60 mg) have been used, and the device was set to operate at a fluid rate of 5 ml/min. Before incorporating water samples into the instrument, 3 ml of acetonitrile–water was added into the concurrent SPE tubes, followed by 3 ml of Milli-Q at a continual fluid percentage of 3 ml/min. SPE containers are, therefore, placed through with a suction drying for 2 min. Consequently, the preferred antibiotic concentration levels contained in the SPE cylinder were retrieved using a combination of acetonitrile and water (2.5 ml) and methanol and formic acid, accompanied by 1 ml of Milli-Q water. Untreated sewage specimens were reduced from 10 to 1 ml by rinsing the molecules in a stream of nitrogen at 35 °C. The analysis was performed on a 1290 HPLC system utilizing HPLC-MS/MS to use a ZORBAX 1.8 m RRHD (502.1 mm) section (Matuszewski et al. 2003). In an intensity mode, 0.1% formic acid (A) was disintegrated in acetonitrile (B) at a flow rate of 0.4 ml/min, with the following criteria: 95% (A) for 1 min; from 95 to 50% (A) for 2 min; maintains at 50% (A) for 2.6 min; from 50 to 95% (A) for 2 min. A positive development stage evaluation was carried out on a Triple Quadruple 6400 ESI mass spectrometer by Ali et al. (2015). The MRM peak value transition was utilized for the monitoring, and the eventual change in results was controlled for quality control.

Data analyses and statistical methods

The mean observations categorized in the same human activities have been presented in graphs to assess the influence of various human activities in water samples using Origin 2018 software. Analysis of variance (ANOVA) was used to determine the statistical relationship between the influence of multiple human activities on the quality of water (ANOVA). To detect a statistical difference between the four different categories of human-caused activities, T-tests were performed to compare the amount that each pair of variables changed throughout the seasons in regards to the river water's quality using Microsoft Excel 2016 (Charfi et al. 2013). The criterion for evaluating whether or not a variation was substantial was a p-value of lower than 0.017. Each sample spot to measure physicochemical attributes, nutrition levels, and heavy metal content in factor analysis and principal component analysis (PCA) was executed using Origin 2018. The method detection limit for all variables below the measurement thresholds within the input data for the multivariate analysis was set at 50%.

The concentration of heavy metals

One-way ANOVA showed that there was a significant difference; Station 7 was the source of the variation. Electrical conductivity measurements fluctuated around 0.09 and 0.26 mS/cm. From February to June, some levels reported at Station 5 were too excessive, while in January, certain values recorded at Station 7 were too high. An ANOVA showed statistically significant variability, with Stations 4 and 7 being the major contributors. Between 7.33 and 7.72, the pH concentrations were recorded. The outcomes of a one-way ANOVA testing exhibit that it was a statistically considerable variation throughout all locations. Turbidity readings varied between 1.13 and 2.74 mg/l. Each value was within the allowable range. The summary of the heavy metal content is presented in Table 1.

Table 1

The concentration of heavy metals in the Mansehra district at the Kunhar River

ParametersUnitsS1S2S3S4S5S6S7S8S9S10Acceptable limitsLimit of detection
Pb μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 30 
Na μg/l 1,208.6 1,121.8 951.3 1,027.5 1,013.5 1,204.8 1,204.3 1,104.3 1,007.5 981.2 20,000 240 
Cu μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 2,000 30 
Cd μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 10 10 
Zn μg/l 10.8 BDL 9.3 10.6 10.8 9.8 BDL BDL 10.7 BDL 5,000 10 
Co μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL – 10 
Hg μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 10 
Se μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 20 10 
As μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 40 
Fe μg/l 19.8 14.8 15.4 17.6 BDL 18.2 BDL BDL BDL BDL 300 10 
Cr μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 20 
TSS mg/l 04 10 6.5 10.2 9.7 150–200 – 
TDS mg/l 64 78 50 127 141 156 159 148 138 142 500 – 
Turbidity NTU 1.13 1.34 1.26 1.15 0.86 2.49 2.48 2.68 2.56 2.74 ≤5 0.31 
Fluoride mg/l 0.08 0.06 0.05 0.08 0.10 0.13 0.16 0.12 0.11 0.12 0.04 
Hardness mg/l 58 71 98 109 122 148 144 149 142 134 100–300 
pH – 7.33 7.31 7.29 7.46 7.68 7.54 7.59 7.74 7.64 7.72 6.5–8.5 – 
Sulfate mg/l 14 15 17 24 37 31 31 33 26 28 200 0.24 
Phosphate mg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL – 0.05 
Temp. °C 9.7 9.8 10 9.5 12.1 11.3 10.8 11.4 9.8 12.3  – 
BOD mg/l 37 42 84 67 52 86 34 98 94 73 – 
EC mS/cm 0.09 0.12 0.13 0.8 0.10 0.07 0.23 0.24 0.14 0.26 1–2 – 
Chloride mg/l 22 27 20 43 45 38 54 67 65 68 200–250 
Nitrate mg/l 1.5 11 12 11 10 10.2 4.8 10 0.03 
ParametersUnitsS1S2S3S4S5S6S7S8S9S10Acceptable limitsLimit of detection
Pb μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 30 
Na μg/l 1,208.6 1,121.8 951.3 1,027.5 1,013.5 1,204.8 1,204.3 1,104.3 1,007.5 981.2 20,000 240 
Cu μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 2,000 30 
Cd μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 10 10 
Zn μg/l 10.8 BDL 9.3 10.6 10.8 9.8 BDL BDL 10.7 BDL 5,000 10 
Co μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL – 10 
Hg μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 10 
Se μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 20 10 
As μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 40 
Fe μg/l 19.8 14.8 15.4 17.6 BDL 18.2 BDL BDL BDL BDL 300 10 
Cr μg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL 50 20 
TSS mg/l 04 10 6.5 10.2 9.7 150–200 – 
TDS mg/l 64 78 50 127 141 156 159 148 138 142 500 – 
Turbidity NTU 1.13 1.34 1.26 1.15 0.86 2.49 2.48 2.68 2.56 2.74 ≤5 0.31 
Fluoride mg/l 0.08 0.06 0.05 0.08 0.10 0.13 0.16 0.12 0.11 0.12 0.04 
Hardness mg/l 58 71 98 109 122 148 144 149 142 134 100–300 
pH – 7.33 7.31 7.29 7.46 7.68 7.54 7.59 7.74 7.64 7.72 6.5–8.5 – 
Sulfate mg/l 14 15 17 24 37 31 31 33 26 28 200 0.24 
Phosphate mg/l BDL BDL BDL BDL BDL BDL BDL BDL BDL BDL – 0.05 
Temp. °C 9.7 9.8 10 9.5 12.1 11.3 10.8 11.4 9.8 12.3  – 
BOD mg/l 37 42 84 67 52 86 34 98 94 73 – 
EC mS/cm 0.09 0.12 0.13 0.8 0.10 0.07 0.23 0.24 0.14 0.26 1–2 – 
Chloride mg/l 22 27 20 43 45 38 54 67 65 68 200–250 
Nitrate mg/l 1.5 11 12 11 10 10.2 4.8 10 0.03 

Note: Acceptable limits WHO 2011 (World Health Organization 2011).

Each ANOVA revealed a statistically substantial variation between locations, indicating Station 2 as the source of variability. The TSS values were as follows: 0.4–9.7 mg/l. In all the stations, the results were between tolerable limits. According to a one-way ANOVA, Station 4 was the main source of the unique variance. The TDS levels were relatively constant, ranging between 64 and 142 mg/l. Depending on the outcomes of a one-way ANOVA, it is obvious that Stations 2 and 3 are the primary contributors to the variance observed. Chloride amounts were determined to be between 22 and 68 mg/l. Station 2 and 3 data from March 2021 to September 2021 revealed certain readings that were too high. An ANOVA using only the results at Stations 2 and 3 indicated a statistically considerable change. Nitrate concentrations were estimated to be anywhere from 0.02 to 0.09 mg/l. All assessed amounts were found to be within the allowed limits. A considerable difference was found, with Station 1 being the main contributor, according to a one-way ANOVA analysis.

Health risk assessment

Risk assessment of daily chronic intake

Iron (Fe) results were as follows: 0.068 mg/kg/day for adults and 0.157 mg/kg/day for children at Station 4; 0.0482 mg/kg/day for adults and 0.0113 mg/kg/day for children at Station 7; and 0.0430 mg/kg/day for adults and 0.100 mg/kg/day for children at Station 3. Only in Station 2 was the maximum permissible concentration of manganese detected. The recommended daily intake (RDI) for adults was 0.008 mg/kg, while the RDA for children was 0.019 mg/kg. The chronic daily intake (CDI) of the assessed heavy metals is presented in Figure 4.
Figure 4

Chronic intake of the daily heavy metals (mg/kg/day).

Figure 4

Chronic intake of the daily heavy metals (mg/kg/day).

Close modal
Figure 5

HQ values and total HI of the heavy metals.

Figure 5

HQ values and total HI of the heavy metals.

Close modal

Zinc (Zn) CDI results for children and adults in Station 5 were 0.031 and 0.072 mg/kg/day, respectively. Such values were 0.024 and 0.056 mg/kg/day in Station 3 for children, and 0.018, and 0.042 mg/kg/day in Station 3 for adults. In Station 7, the Cr levels for children and adults were 0.002 and 0.006 mg/kg/day, respectively. In Stations 1 and 4, the Cr values for children and adults were 0.001 and 0.003 mg/kg/day, respectively. Cadmium (Cd) values for adults and children in Station 8 were 0.002 and 0.004 mg/kg/day, respectively. In Stations 9 and 6, the values for adults and children were 0.001 and 0.003 mg/kg/day, respectively. Nickel was also only found in Station 2, where it was too high. The daily chronic intake for adults was 0.001 mg/kg/day and for children, it was 0.002 mg/kg/day.

Risk assessment of HQ

For children and adults, the Fe is 9.71 and 22.4 in Station 1; 6.85 and 16.41 in Station 2; and 6.14 and 14.28 in Station 3, respectively. For both adults and children, the Fe was more than 1 in all the stations. In Station 1, the Mn for adults is 0.57 and for children it is 1.36. The children's HQ was greater than 1. In Station 1, the Zn for adults is 0.10 and for children, it is 0.24. In Station 2, it is 0.08 and in Station 3, it is 0.06 and 0.14. In Stations 1, 2, and 3, the Zn for both adults and children was less than 1. In Station 1, the Cr for adults and children is 0.67 and 20, respectively; in Station 2, it is 3.33 and 10; and in Station 3, it is 3.33 and 3.33, respectively. The Cr for children in Station 1 was higher than for adults and children in Stations 2 and 3. The HQ and CD for adults and children is stations 1, 2, and 6 in Station 2, and stations 2 and 2 in Station 3. In all of the stations, the HQ and CD for both adults and children was higher than 1. Ni for adults and children is only 0.05 and 0.1 at Station 1. The HQ was less than 1 for both adults and children. The HQ of the assessed heavy metals is presented in Figure 5.

Concentration of antibiotic residues

Ten locations along the Kunhar River were tested to assess the concentrations of six different antibiotics in the surface water of the river using the LC–MSMS method as shown in Figure 6. For a wide range of highly adequate antibiotics (amoxicillin, amikacin, cefixime, and ciprofloxacin) excluding the sulfamethoxazole and trimethoprim, all study sites might have concentrations below the strategy reference range (Figure 7). Sulfamethoxazole and trimethoprim exhibited very low abundances throughout all sites assessed, apart from industrial areas (0.3 mg/l sulfamethoxazole, 0.01 mg/l trimethoprim) and (0.11 mg/l sulfamethoxazole, 0.013 mg/l for trimethoprim). Six antibiotics examined, sulfamethoxazole, and trimethoprim were most often found in the water (Nguyen Dang Giang et al. 2015). The concentration levels of these antibiotics were 0.021 and 0.01 mg/l , respectively (Xue et al. 2022). Surface water antibiotic residues could provide an evolutionary advantage for the growth and spreading of antibiotic-resistant organisms.
Figure 6

Concentration of heavy metals affected by different anthropogenic activities.

Figure 6

Concentration of heavy metals affected by different anthropogenic activities.

Close modal
Figure 7

The antibiotic concentration of residues in the surface water along the Kunhar River.

Figure 7

The antibiotic concentration of residues in the surface water along the Kunhar River.

Close modal

Heavy metal content

The metal concentration levels at all monitoring sites have been above the adherence to safety thresholds, which may be responsible for the increase in anthropogenic and natural impacts due to human activities, particularly at Stations 1 and 2. When compared to results from other studies, the results obtained here were lower. Values in the range from 0.75 to 15.01 mg/l were found in a few northern parts of Pakistan's water sources, whereas agrarian warmer weather in the northern part measured levels around 1.12 and 5.10 mg/l. The element iron (Fe) is essential for the survival of almost all known living things. Iron is a crucial element of protein molecules, enzymatic, and hemoglobin. As a result of its frequent distribution in the natural environment, iron tends to be found at a higher concentration in Pakistan's freshwater ecosystem than other metals. Fe inadequacy, which can cause anemia and fatigue, is most common in children (those under 5 years of age), pregnant women, and people with compromised immune systems (Clark 2008). An excess of iron can also have harmful consequences for the nervous system, whereas the symptomatic iron toxic effects may induce mitochondrial dysfunction due to metal oxidative stress to the genetic material (Jomova & Valko 2011). Several of the evaluated concentration levels of heavy metals were observed to be incredibly high, particularly at Stations 1 and 2, and this could be due to the combined effects of all the human activities occurring in and around Station 1 and the downstream impacts of the effluent from Station 2 upstream outflow. Overall, the levels detected in this research dropped within the earlier assessed range of 0.06–0.51 mg/l (Iqbal & Shah 2015) in the northern region of Pakistan during the spring and early summer. Due to its contribution as a co-factor in various enzyme functions, manganese is an essential trace mineral for all living things (Hüssy et al. 2021). There is a lot of biochemical importance to heavy metals, and it is not harmful at all. Manganese-rich water is the norm instead of the exception in most regions. Manganism is a chronic disorder that occurs if too much manganese is consumed. Zn limits in drinkable water have been linked to children's recognition and behavior problems (Giri & Singh 2015). Manganese deficiency anemia can generate from persistent exposure to excessive amounts of heavy metals because it impedes the body's capacity to absorb the iron from nutrition. As with copper, manganese ingestion decreases the effectiveness of copper metalloenzymes. The intake of too much manganese can lead to hypertension in individuals over the age of 40 years and symptoms emulating Parkinson's disease (Malavolta & Mocchegiani 2018). Particularly at Stations 1 and 2, all recorded concentration levels of Zn were significantly higher than the upper limit expected. It is also feasible that geomorphological factors or human activities are reprimanded. Zn is also essential for healthy animal and human growth and survival. Its compounds or natural complexes can be found in almost all edible and consumable water bodies (Mandal et al. 1998). The growth of bone fragments and the development and function of the reproductive system are some of the infrequent consequences of zinc toxicity. Among the clinical symptoms of zinc toxic effects are digestive problems, dark urine, organ failure, kidney problems, and metabolic disorders (Cummings & Kovacic 2009). Since all of the measured concentration levels of metals were within the original study limit for acceptance, researchers have found no evidence of an issue. Anthropogenic activities could account for the comparatively high value observed at Stations 1 and 2. Concentration levels as low as 0.02 and 0.15 mg/l were evaluated in the spring and summer time in Pakistan, while (Ekere et al. 2014) few drinking sources of water of Northern Pakistan had already evaluated concentration levels of 0.21–2.65 mg/l. The evaluated concentration levels of metals were considerably higher than the permissible level. All these geomorphological and living thing aspects could be in practice, particularly at Stations 1 and 2. Once chromium thresholds in drinkable water are kept within the normal limit, they aid in the maintenance of healthy glucose uptake, trying to make chromium a vital nutrient for human beings. Disease and inflammation of the skin can generate when chromium concentration levels are too high. Stations 1 and 2 specifically had cadmium amounts that have been substantially greater than the recommended peak value, which could be due to human interference. It is not highly essential, but it can still be toxic because it contains impurities (Murphy et al. 2012). Cadmium, even in smaller concentrations, can harm the arteries of the kidneys. Chemically, Cd can substitute for Zn, resulting in higher blood pressure and kidney problems. Also, it induces human immune disease, which is incredibly painful and has consequences from enzyme intervention. High blood pressure, organ failure, and even the advancement of tumors are all possible scenarios of intense Cd toxicity in living beings (Engwa et al. 2019). Ni limits at Stations 1 and 2 have been observed every once in a while to be too high. Numerous distinct varieties of animals, microbes, and plants depend on Ni, and as a consequence, deficit or toxic effects could take place if there was not a sufficient amount of it (Cempel & Nikel 2006). Despite becoming widespread and essential for the mechanism of several living things, nickel's toxic effects could be precipitated by either human activity release or geological formations changing circumstances in the same regions. Those who are exposed to the toxic effects may encounter problems with their breathing and immune function. The majority of individuals are exposed to Ni largely through their nutrition and drinkable water (Puckett 1995).

Chronic daily intake

Iron CDI levels both for adults and children were greater than that of the RfD (0.007 mg/kg/day) when consumed by humans. It is indeed feasible that the pH level of the river appears to contribute to the increased iron CDI value systems by trying to make the river's naturally high mineral content more drinkable. According to Lenntech (n.d.), a reduction in the pH enhances the water solubility from certain metals (Onyele & Anyanwu 2018); in a similar way, both show increased iron concentrations in the CDI. As an outcome, river water could pose a risk to individuals who drink something due to the presence of iron. The CDI amount of heavy metals for an adult was below the oral standard value (RfD) of 0.014 mg/kg/day, whereas that of children surpassed the RfD. Manganese was not regarded as a possible future risk to adults exposed to the water of Kunhar but not to the children. Zinc concentration levels in Kunhar River water were lower than the oral standard value (RfD) (0.3 mg/kg/day) both for children and adults (Muhammad et al. 2021). Živković et al. (2019) recorded lower CDI values for Zn. CDI values of chromium for adults and children were above the oral reference dose (RfD) (0.0003 mg/kg/day). All these natural and anthropogenic factors could contribute to increased cadmium CDI measurements. Consequently, river water containing cadmium may pose a risk for those who drink it. Nickel concentration levels primarily due to ingestion in Station 1 are significantly below the oral standard value (RfD) of 0.02 mg/kg/day for both adults and children.

HQ and HI

All the stations had a HQ for metals greater than 1. Each of the stations had much greater iron levels (1), and only Station 1 and the children's concentration of manganese were too greater. All stations had zinc levels below the safe threshold for both adults and children. All these children and adults at all stations had excessively high cadmium amounts, and chromium levels were too high except for adults at Station 1. It was discovered that children were especially vulnerable because of the high HQ values from certain age categories (Mohammadi et al. 2019). Some of the metals perceived in this research also had an HQ value greater than 1. Some of the metals examined in this research had extremely high CDI values, which also likely contributed to the large HQ observed values. Water customers in all impacted stations may be exposed to such metals over a lengthy period. Because of the severity of the measured long-term health threat, the non-carcinogenic negative impact could be ignored as superficial.

Anthropogenic activities and water quality parameters in the Mansehra district

All surface water samples gathered from the Kunhar River during the wet and dry seasons were confined to the PCA and correlational analysis (e.g., Pearson and Spearman analysis) to evaluate whether or not there was a correlation between the physical and chemical parameters, essential minerals, toxic substances, and anthropogenic impacts. The PCA identified eight essential aspects, that together took into account 86.67% of the overall variance in the set of data. In total, 75.1% of the variability in parameters of water quality might be clarified by the five fundamental elements that were obtained. Positive maximum load values higher than 0.3 are connected with pH, temperature, BOD5, Cu, Ni, and Zn on the initial singular value (PC1), which helps to explain 28.16% of the overall variability. Organic material and toxic substances (for instance, Cu, Ni, and Zn) in the Kunhar River may well be predominantly polluted by domestic and commercial wastes in the Mansehra district, as shown by the hopeful similarity between the PC1 and many samples taken at the residential and commercial locations. Samples from agrarian and less affected areas were negatively associated, revealing that agricultural practices in the Mansehra district did not result in significant organic material and heavy metal contaminations in the Kunhar River. This may be a result of the fact that the selected sample water was situated in the center of the Mansehra district, where it revealed both industrial effluents and domestic sewage from residential areas. Positive charging values ranging from 0.39 to 0.48 were coupled with conductivity, salinity, and DO on the second principal component, which also took into account 16.95% of the amount of variance. A perfect coefficient between residences and less impacted site specimens recommended that ionic compounds were prevalent in the surface water collected there. Approximately 10 and 7% of the amount of variance of all parameters of water quality were clarified by the constituents. Both residential and commercial wastewater may have an impact on the elements due to their similarity with organic matter and heavy metals.

There might have been an essential assessment of health risks because a few of the toxic substances evaluated were just too intense. There was no single heavy metal found to be a cause of the health risk to the stations, but the combined impact of all of them was fully evident in the HI. So, every station's HI was substantially higher than the threshold of 1. It has serious implications for the health of both children and adults who drink polluted water. To a great extent, anthropogenic activities accelerated naturally produced surface and groundwater sources of contamination with heavy metals. According to the limited working seasons for field involvements in the Himalayas, maintenance of high-altitude areas surrounding the Kunhar River is indeed a massive task for the globe. According to the results, the Kunhar River's water quality is excellent, except for a slightly higher BOD value that can be minimized by preventing the dumping of human-oriented trash around the river or stream. According to the metal concentration in the sediment, sediment particles may act as a drain for metals because it adsorbs particles and carries the largest environmental contamination load. Current quantities of metals in water present no risk to people or aquatic life. Metals in river sediments have been shown to have a minor concentration factor and contamination factor, which might arise with the period if human behaviors along the river are not controlled, which impacts air quality, sanitation and hygiene, and food and nutrition. Predictors of climatic variability effects on water quality and quantity may benefit from our research. As a result of the increased desire for the tourism sector, the land is now under increased stress, increasing the threat of soil degradation, contamination, the destruction of the environmental ecosystems, and the risk to threatened animals. There was also a lot of plastic trash accumulated on the riverbank, which is an increasing environmental danger. These impacts tend to slowly destroy the environmental and social resources that maintain the tourist industry. Climate change has been the major cause of river extinction in the past, so monitoring and management of natural resources are critical to ensure that rivers will not become completely extinct.

It is the author's wish to acknowledge and thank the Pakistan Meteorological Department (PMD) for providing valuable data for this research.

S-e-h.S. conceptualized the study; prepared the methodology; wrote and prepared the original draft; wrote, reviewed, and edited the article. X.S. conceptualized the study, acquired funds, was involved in project administration, collected resources, and supervised the study. J.G. conceptualized the study, acquired funds, wrote, reviewed, and edited the article. C.H. wrote, reviewed, and edited the article. S.J. and H.M.Z. wrote, reviewed, and edited the article, did data curation, and analysed the data.

This work was supported by the Projects of Natural Natural Science Foundation of China (51922065 and 52279069) and the National Natural Science Foundation of China (52179018).

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Alansi
R. Q.
,
Mohammed
A. M.
,
Ali
M. M.
,
Mokbel Ghalib
W. A.
&
Ponnappa
S. C.
2021
Determination of heavy metals in groundwater around Al-Buraihi sewage station in Taiz City, Yemen
.
Journal of Health Pollution
11
(
30
),
210604
.
Ali
W.
&
Muhammad
S.
2022
Spatial distribution of contaminants and water quality assessment using an indexical approach, Astore River basin, Western Himalayas, Northern Pakistan
.
Geocarto International
37, 14005–14026.
Asfaw
A.
,
Maher
K.
&
Shucksmith
J.
2018
Modelling of metaldehyde concentrations in surface waters: a travel time based approach
.
Journal of Hydrology
562
,
397
410
.
Bhateria
R.
&
Jain
D.
2016
Water quality assessment of lake water: a review
.
Sustainable Water Resources Management
2
(
2
),
161
173
.
Bloise
A.
,
Barca
D.
,
Gualtieri
A. F.
,
Pollastri
S.
&
Belluso
E.
2016
Trace elements in hazardous mineral fibres
.
Environmental Pollution
216
,
314
323
.
Cempel
M.
&
Nikel
G.
2006
Nickel: a review of its sources and environmental toxicology
.
Polish Journal of Environmental Studies
15
(
3
).
Clark
S. F.
2008
Iron deficiency anemia
.
Nutrition in Clinical Practice
23
(
2
),
128
141
.
Coffer, M. M., Schaeffer, B. A., Foreman, K., Porteous, A., Loftin, K. A., Stumpf, R. P., P. J. Werdell., Urquhart, E., Ryan J. Albert, R. J. & Darling, J. A.
2021
Assessing cyanobacterial frequency and abundance at surface waters near drinking water intakes across the United States
.
Water Research
201
,
117377
.
Cummings
J. E.
&
Kovacic
J. P.
2009
The ubiquitous role of zinc in health and disease
.
Journal of Veterinary Emergency and Critical Care
19
(
3
),
215
240
.
Ekere
N. R.
,
Hedioha
J. N.
,
Eze
I. S.
&
Agbazue
V. E.
2014
Health risk assessment in relation to heavy metals in water sources in rural regions of South East Nigeria
.
International Journal of Physical Sciences
9
(
6
),
109
116
.
Engwa
G. A.
,
Ferdinand
P. U.
,
Nwalo
F. N.
&
Unachukwu
M. N.
2019
Mechanism and health effects of heavy metal toxicity in humans
.
Poisoning in the Modern World-New Tricks for an Old Dog
10
,
70
90
.
Finch
C.
,
Roldan
R.
,
Walsh
L.
,
Kelly
J.
&
Amor
S.
2018
Analytical Methods for Chemical Analysis of Geological Materials
.
Government of Newfoundland and Labrador, Department of Natural Resources, Geological Survey, Open File NFLD/3316
, St. Johns, NL.
Gorde
S.
&
Jadhav
M.
2013
Assessment of water quality parameters: a review
.
Journal of Engineering Research and Applications
3
(
6
),
2029
2035
.
Gyamfi
E.
,
Appiah-Adjei
E. K.
&
Adjei
K. A.
2019
Potential heavy metal pollution of soil and water resources from artisanal mining in Kokoteasua, Ghana
.
Groundwater for Sustainable Development
8
,
450
456
.
He, M.-Y., Dong, J.-B., Jin, Z., Liu, C.-Y., Xiao, J., Zhang, F., Sun, h., Zhao, Z. Q., Gou, L. F., Liu, W. G., Luo, C. G., Song, Y. G., Ma, L. & Deng, L.
2021
Pedogenic processes in loess-paleosol sediments: clues from Li isotopes of leachate in Luochuan loess
.
Geochimica et Cosmochimica Acta
299
,
151
162
.
Hong
B.
,
Yu
S.
,
Zhou
M.
,
Li
J.
,
Ding
J.
&
Niu
Y.
2021
Development of a pH-paralleling approach of quantifying six-category pharmaceuticals in surface water using SPE-HPLC-MS/MS
.
Watershed Ecology and the Environment
3
,
1
16
.
Hüssy, K., Limburg, K. E., de Pontual, H., Thomas, O. R., Cook, P. K., Heimbrand, Y., Blass, M. & Sturrock, A. M.
2021
Trace element patterns in otoliths: the role of biomineralization
.
Reviews in Fisheries Science & Aquaculture
29
(
4
),
445
477
.
Iqbal
J.
&
Shah
M. H.
2015
Study of selected metals distribution, source apportionment, and risk assessment in suburban soil, Pakistan
.
Journal of Chemistry
2015
(2–3), 1–8.
Johnson
C. E.
,
Driscoll
C. T.
,
Blum
J. D.
,
Fahey
T. J.
&
Battles
J. J.
2014
Soil chemical dynamics after calcium silicate addition to a northern hardwood forest
.
Soil Science Society of America Journal
78
(
4
),
1458
1468
.
Jomova
K.
&
Valko
M.
2011
Advances in metal-induced oxidative stress and human disease
.
Toxicology
283
(
2–3
),
65
87
.
Kandhro
A. J.
,
Samoan
N. A.
,
Laghari
J. H.
,
Chandio
A. M.
&
Yousfani
A. H.
2016
Assessment of arsenic and essential metal ions in the quality of groundwater sources of Taluka Daur, District Shaheed Benazeer Abad, Sindh, Pakistan
.
Journal of Peoples University of Medical & Health Sciences Nawabshah.(JPUMHS)
6
(
1
),
01
08
.
Li
J.
,
Zhao
Y.
,
Zhang
A.
,
Song
B.
&
Hill
R. L.
2021
Effect of grazing exclusion on nitrous oxide emissions during freeze-thaw cycles in a typical steppe of Inner Mongolia
.
Agriculture, Ecosystems & Environment
307
,
107217
.
Lin, X., Lu, K., Hardison, A. K., Liu, Z., Xu, X., Gao, D., Gong, J. & Gardner, W. S.
2021
Membrane inlet mass spectrometry method (REOX/MIMS) to measure 15N-nitrate in isotope-enrichment experiments
.
Ecological Indicators
126
,
107639
..
López-Delgado
A.
,
Robla
J. I.
,
Padilla
I.
,
López-Andrés
S.
&
Romero
M.
2020
Zero-waste process for the transformation of a hazardous aluminum waste into a raw material to obtain zeolites
.
Journal of Cleaner Production
255
,
120178
.
Madadrang, C. J., Kim, H. Y., Gao, G., Wang, N., Zhu, J., Feng, H., Gorring, M., Kasner, M. L. & Hou, S.
2012
Adsorption behavior of EDTA-graphene oxide for Pb (II) removal
.
ACS Applied Materials & Interfaces
4
(
3
),
1186
1193
.
Malavolta
M.
&
Mocchegiani
E.
2018
Trace Elements and Minerals in Health and Longevity
, Vol.
8
.
Springer
, Cham.
Mandal
B.
,
Chowdhury
T. R.
,
Samanta
G.
,
Mukherjee
D.
,
Chanda
C.
,
Saha
K.
&
Chakraborti
D.
1998
Impact of safe water for drinking and cooking on five arsenic-affected families for 2 years in West Bengal, India
.
Science of the Total Environment
218
(
2–3
),
185
201
.
Matuszewski
B. K.
,
Constanzer
M.
&
Chavez-Eng
C.
2003
Strategies for the assessment of matrix effect in quantitative bioanalytical methods based on HPLC− MS/MS
.
Analytical Chemistry
75
(
13
),
3019
3030
.
Mirrione
M. M.
,
Schulz
D.
,
Lapidus
K. A.
,
Zhang
S.
,
Goodman
W.
&
Henn
F. A.
2014
Increased metabolic activity in the septum and habenula during stress is linked to subsequent expression of learned helplessness behavior
.
Frontiers in Human Neuroscience
8
,
29
.
Mohammadi, A. A., Zarei, A., Majidi, S., Ghaderpoury, A., Hashempour, Y., Saghi, M. H., Alinejad, A., Yousefi, M., Hosseingholizadeh, N. & Ghaderpoori, M.
2019
Carcinogenic and non-carcinogenic health risk assessment of heavy metals in drinking water of Khorramabad, Iran
.
MethodsX
6
,
1642
1651
.
Muhammad
S.
,
Ullah
S.
,
Ali
W.
,
Jadoon
I. A.
&
Arif
M.
2021
Spatial distribution of heavy metal and risk indices of water and sediments in the Kunhar River and its tributaries
.
Geocarto International
, 37 (20), 5985–6003.
Murphy
E. A.
,
Post
G. B.
,
Buckley
B. T.
,
Lippincott
R. L.
&
Robson
M. G.
2012
Future challenges to protecting public health from drinking-water contaminants
.
Annual Review of Public Health
33
,
209
.
Nguyen Dang Giang
C.
,
Sebesvari
Z.
,
Renaud
F.
,
Rosendahl
I.
,
Hoang Minh
Q.
&
Amelung
W.
2015
Occurrence and dissipation of the antibiotics sulfamethoxazole, sulfadiazine, trimethoprim, and enrofloxacin in the Mekong Delta, Vietnam
.
Plos one
10
(
7
),
e0131855
.
Onyele
O. G.
&
Anyanwu
E. D.
2018
Human health risk assessment of some heavy metals in a rural spring, southeastern Nigeria
.
Afr J Environ Nat Sci Res
1
(
1
),
15
23
.
Paustenbach
D. J.
2015
Human and Ecological Risk Assessment: Theory and Practice (Wiley Classics Library)
.
John Wiley & Sons
, Hoboken, NJ.
Puckett
L. J.
1995
Identifying the major sources of nutrient water pollution
.
Environmental Science & Technology
29
(
9
),
408A
414A
.
Qiu
H.
,
Gui
H.
,
Fang
P.
&
Li
G.
2021
Groundwater pollution and human health risk based on Monte Carlo simulation in a typical mining area in Northern Anhui Province, China
.
International Journal of Coal Science & Technology
8
(
5
),
1118
1129
.
Rebello
S.
,
Anoopkumar
A.
,
Aneesh
E. M.
,
Sindhu
R.
,
Binod
P.
,
Kim
S. H.
&
Pandey
A.
2021
Hazardous minerals mining: challenges and solutions
.
Journal of Hazardous Materials
402
,
123474
.
Silva, L. F., Oliveira, M. L., Crissien, T. J., Santosh, M., Bolivar, J., Shao, L., Dotto, G. L., Gasparotto, J. & Schindler, M.
2022
A review on the environmental impact of phosphogypsum and potential health impacts through the release of nanoparticles
.
Chemosphere
286
,
131513
.
Soomro
S.-e.-h.
,
Hu
C.
,
Jian
S.
,
Wu
Q.
,
Boota
M. W.
&
Soomro
M. H. A. A.
2021a
Precipitation changes and their relationships with vegetation responses during 1982–2015 in Kunhar River basin, Pakistan
.
Water Supply
21
(
7
),
3657
3671
.
Soomro
S.-e.-h.
,
Hu
C.
,
Boota
M. W.
,
Soomro
M. H. A. A.
,
Jian
S.
,
Zafar
Z.
&
Li
X.
2021b
Mapping flood extend and its impact on land use/land cover and settlements variations: a case study of Layyah District, Punjab, Pakistan
.
Acta Geophysica
69
,
2291
2304
.
Soomro, S.-e.-h., Hu, C., Boota, M. W., Ahmed, Z., Chengshuai, L., Zhenyue, H., Xiang, L. & Soomro, M. H. A. A.
2022
River flood susceptibility and basin maturity analyzed using a coupled approach of Geo-morphometric parameters and SWAT model
.
Water Resources Management
36
(
7
),
2131
2160
.
Tyopine
A. A.
,
Sikakwe
G. U.
,
Obalum
S. E.
&
Okoye
C. O.
2020
Relative distribution of rare-earth metals alongside alkaline earth and alkali metals in rhizosphere of agricultural soils in humid tropical environment
.
Environmental Monitoring and Assessment
192
(
8
),
1
13
.
USEPA
1999
Guidance for Performing Aggregate Exposure and Risk Assessments
.
Office of Pesticide Programs
,
Washington, DC
.
Verma
S. K.
&
Sharma
P. C.
2022
Current trends in solid tannery waste management
.
Critical Reviews in Biotechnology
1
18
.
World Health Organization
2011
Safe Drinking-Water From Desalination
.
Xue, J., Lei, D., Zhao, X., Hu, Y., Yao, S., Lin, K., Wang, Z. & Cui, C.
2022
Antibiotic residue and toxicity assessment of wastewater during the pharmaceutical production processes
.
Chemosphere
291
,
132837
.
Yao
J.
,
Weng
Y.
,
Dickey
A.
&
Wang
K. Y.
2015
Plants as factories for human pharmaceuticals: applications and challenges
.
International Journal of Molecular Sciences
16
(
12
),
28549
28565
.
Zhang, X., Ma, F., Yin, S., Wallace, C. D., Soltanian, M. R., Dai, Z., Ritzi, R. W., Ma, Z., Zhan, C. & Lü, X.
2021
Application of upscaling methods for fluid flow and mass transport in multi-scale heterogeneous media: a critical review
.
Applied Energy
303
,
117603
.
Zhang
X.
,
Wang
Z.
,
Reimus
P.
,
Ma
F.
,
Soltanian
M. R.
,
Xing
B.
, Zang, J., Wang, Y. &
Dai
Z.
2022
Plutonium reactive transport in fractured granite: multi-species experiments and simulations
.
Water Research
224
,
119068
.
Zhao
L.
,
Du
M.
,
Du
W.
,
Guo
J.
,
Liao
Z.
,
Kang
X.
&
Liu
Q.
2022
Evaluation of the carbon sink capacity of the proposed Kunlun Mountain National Park
.
International Journal of Environmental Research and Public Health
19
(
16
),
9887
.
Živković
N.
,
Takić
L.
,
Djordjević
L.
,
Djordjević
A.
,
Mladenović-Ranisavljević
I.
,
Golubović
T.
&
Božilov
A.
2019
Concentrations of heavy metal cations and a health risk assessment of sediments and river surface water: a case study from a Serbian mine
.
Polish Journal of Environmental Studies
28
(
3
),
2009
2020
.

Author notes

These authors contributed equally to this paper

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).