Impacts of wastewater discharge from Kalar city on Diyala-Sirwan river water quality, Iraq: pollution evaluation, health risks of heavy metals contamination

In this work, a comprehensive assessment has been performed to understand impact of wastewater discharge from Kalar city on water quality parameters of Diyala-Sirwan river, Iraq. Levels of physicochemical parameters of Ca, K, Na, Mg, pH, electrical conductivity, dissolved oxygen, total dissolved solid, turbidity, biological oxygen demand, and total hardness, and heavy metals of Fe, Mn, Ni, Cu, Zn, Hg, Al, Ba, Se, and As have been determined by inductively coupled plasma optical emission spectroscopy and other devices. Pollution indices, health risk assessment, and multivariate analysis have been applied to evaluate pollution intensity originated from wastewater discharge on river water quality at different locations, before and after the impact point. Pollution levels of slight to moderate have been identified using different indices for the river water at all sampling locations after the wastewater discharge point. Cluster and correlation analyses showed that the impact of wastewater discharge on the river water quality was occurring along the river within the study area. Non-carcinogenic health risk assessment for heavy metals in the river, indicating a potential risk, might on people's health, especially after discharge point. Furthermore, a significant effect of carcinogenic risk by heavy metals of As and Hg was recognized in the same area. This study thus helps to understand clearly the alteration that happened in the water quality of Diyala-Sirwan river due to the discharge of untreated wastewater from Kalar city.


Introduction
Surface water quality is a key factor in water supply evaluation for authorities in urban and rural areas (Ismaiel et al. 2018). Usually, surface water quality is evaluated depending upon various physical, chemical, and biological parameters (Issa 2017). Since a great number of industrial effluents and sewage are likely to have been discharged into surface water sources, heavy metals are taken into consideration for effective quality evaluation of any surface water body (Kaushik et al. 2009;Mu et al. 2015;Ullah et al. 2019). Levels of heavy metals and other chemical quality parameters in surface waters are continuously changing due to various inputs, involving natural and anthropogenic sources (Ali et al. 2016;Muhammad and Ahmad 2020). Surface runoff and drainage water from rains have main contributions in carrying plenty of natural loads such as muds, soil, and humus into rivers (Begum et al. 2009). Industrial and urban discharges are the principal anthropogenic inputs of heavy metals and other chemicals in surface water bodies (Salah et al. 2012). River contamination with heavy metals is increasingly becoming an important issue, because of toxic, non-biodegradable behavior of heavy metals in aquatic biological systems Varol and Şen 2012). Subsequently, heavy metals are likely to be transferred to humans, animals, and plants through consuming and using this contaminated water, posing serious health and life complications (Protano et al. 2014). Surface water quality assessment is commonly conducted to identify natural and anthropogenic sources causing pollution of surface water within a particular study area. Quality assessment for cases when a pollution source has been identified gives further illustrations about how pollutants were readily dispersed and distributed in the water body.
Many former studies have been performed on river pollution caused by city wastewater discharges. Several of these works found that river pollution in cities is mainly attributed to heavy metal pollution (Abdel-Satar et al. 2017). For river water quality assessment, studies showed that mostly various heavy metals in water samples have high concentrations than the normal levels (Bhuiyan et al. 2014;Diagomanolin et al. 2004;Reza and Singh 2010). Physicochemical parameters are also being used for water quality assessment. Combining heavy metals with other water quality parameters helps to have a general view on surface water quality condition, as physicochemical parameters water may reveal significant pollution situations as well (Özgür et al. 2012). Various ways have been adopted to evaluate heavy metals and other parameters for surface water quality: the degree of contamination index (C d ) to determine combined effects of risky quality parameters in water (Bhuiyan et al. 2010); water quality index (WQI); heavy metals evaluation index (HEI) (Singaraja et al. 2015); pollution index (PLI); heavy metals pollution index (HPI) to identify heavy metal pollution in water (Cengiz et al. 2017;Prasad and Bose 2001;Reza and Singh 2010); and multivariate statistics (Giri and Singh 2014;Nasrabadi 2015;Simeonov et al. 2003;Wu et al. 2018).
Nearly all rivers in Iraq are suffering from low-quality level and pollution (Hassan et al. 2010;Ibrahim 2012;Kamil and Adel Abdulrazzaq 2010), as a result of industrial and rural wastewater discharges. Several attempts have been established to define potential risks on surface water resources from these discharges for different parts of Iraq (Aziz et al. 2018;Aziz and Fakhrey 2016;Issa 2014;Razzak and Sulaymon 2009;Wissam and Isam 2017). Till now, few studies have been conducted to determine heavy metals and other physicochemical concentration levels in Diyala-Sirwan river (Abdullah 2013;Hassan et al. 2017;Issa and Alshatteri 2018). However, the previous studies were generally investigating the water quality of Diyala-Sirwan river without focusing on an identified source of pollution for the river. This work demonstrates the impact of wastewater discharge from Kalar city on increasing heavy metal concentrations and physicochemical parameters in Diyala-Sirwan river, and their consequences on health and life of the people of the region. The evaluation was made by using various heavy metals pollution and risk assessment indices with different statistical and multivariate methods such as analysis of variance (ANOVA), Friedman test, Pearson correlation analysis matrix (CM), and cluster analysis (CA) and box plots to identify pollution levels. This work aids to establish a trustworthy evaluation of wastewater discharges impact on Diyala-Sirwan river quality, by which further development can be achieved in remediation policies that applied to improve the environmental and health programs in the region.

Description of the study area
The study area Kalar city (34° 37′ N, 45° 19′ E) is located in the south of Sulaimaniyah Province, east of Iraq ( Fig. 1), and is about 200 km northeast of Baghdad. Kalar city has become increasingly populated in the last two decades. The population expansion happened due to both natural economic development and to the unusual ongoing political situation in Iraq. Kalar city of about 281 K population (Kurdistan Region Statistics Office 2020), and comprising many commercial and industrials activities with agricultural in suburban areas as well, which lastly are discharged into Diyala-Sirwan river without treatment. The annual rainfall of the city is 273 mm with no precipitation in the summer season (Garmian Region Agriculture Department 2017).

Water sample collection and preparation
In the study area, water samples were collected from sampling locations from four locations along Diyala-Sirwan river in January 2019: L1 before impact point of Kalar city wastewater discharge, by 10 km; L2 directly at the impact point of wastewater discharge; L3 after 5 km from impact point; and L4 at 10 km from the impact point. The selected water sampling locations cover a stretch of about 20 km ( Fig. 1). At each site, three water samples were collected in clean polyethylene containers of a 250 ml volume, in which water samples were allowed to remain in containers only for a while before taking to analysis. All water samples have been prepared for analysis of heavy metals were acidified with 2% nitric acid (pH ˂ 2), and refrigerated and transferred to the instrumental research laboratory, within 1 to 2 h from the time of collection and stored in a dark place. For dilutions and glassware washing, distilled deionized water always has been used (Marcovecchio et al. 2007). Standard solutions used by ICP-OES analysis were prepared by diluting with several dilutions of concentrations of 0.1, 0.5, 2 ppm into 0.5% nitric acid, which was used as a diluent (Aris et al. 2013).

Water samples analysis
In the laboratory, all returned water samples for heavy metals analysis and other chemical elements were acidified by adding concentrated nitric acid HNO 3 and stored at 25 0 C. The analysis was done at the University of Garmian laboratories by using inductively coupled plasma optical emission spectroscopy, ICP-OES, (Spectro across Germany). The standard solutions were prepared by serial dilutions of the 1000 mg/L used to analysis ten heavy metals of Fe, Mn, Ni, Cu, Zn, Hg, Al, Ba, Se, and As and four chemical elements of Ca, K, Na, and Mg. Instrument conditions of ICP-OES were: RF power/W equals 1400; pump speed was about 30 rpm; coolant flow equals 14 L/min; measure time was 28 s; replicate measurements were 3.
Water samples were also analyzed over the same period for physicochemical water quality parameters of pH, electrical conductivity (EC), dissolved oxygen (DO), total dissolved solid (TDS), turbidity, biological oxygen demand (BOD), and total hardness (TH) at the University of Garmian laboratories by following analysis procedures according to APHA (Rice et al. 2017). These parameters are chosen to evaluate Kalar city wastewater discharge, as they are highly relevant to the occurring Diyala-Sirwan river quality variation concerning the suitability for drinking and various purposes by people of the study area. The accuracy of the analysis method was verified, randomly selected samples their analysis was duplicated.

Heavy metal pollution index (HPI)
HPI as proposed by Mohan et al. (Mohan et al. 1996), by which water quality was assessed in terms of heavy metals importance in water samples. HPI calculation is based on an arithmetic mean a method that converts water existing data into a sole number in terms of heavy metals presence impact on water quality as following equations, M i , S i , and I i are monitored, standard, and ideal values of i-parameter for the investigated heavy metals, Q i is the sub-index of i-parameter, W i is the weight of i-parameter, and n is the total number of parameters in a test. For each parameter, W i is inversely proportional to the recommended standard. (1)

Heavy metals evaluation index (HEI)
HEI is applied to identify water contamination by heavy metals, and heavy metal evaluation index is calculated as following (Boateng et al. 2015), where H c and H mac are the observed and maximum permissible level concentrations for each i-parameter, respectively.

Contamination index (C d )
For C d , a sum of individual contamination parameters above the upper allowable limits as the following (Edet and Offiong 2002), C fi , C Ai , and C Ni are concentration factor, analytical value, and the upper allowable concentration of the i-parameter, respectively.

Water Pollution index (PI) and overall index of pollution (OIP)
Pollution index (PI) evaluates water pollution by metals that originated from different origins. PI is a normalized factor, is a concentration ratio of quality parameters, and their background values as the following equation (Goher et al. 2014).
C i is parameter concentration, and S i is standard (background) concentration (mg/kg) of the same parameter in water sample. For PI, five classifications were considered: no effect pollution (PI < 1.0); slightly affected (1.0 ≤ PI < 2.0); moderately affected (2.0 ≤ PI < 3.0); strongly affected (3.0 ≤ PI ≤ 5.0); and extremely affected (PI > 5.0). To estimate various pollution sources affecting surface water quality, an average of all individual water quality parameters was proposed by (Sargaonkar and Deshpande 2003) as the following equation where Pi is the ith parameter pollution index, and n is the number of parameters included in pollution evaluation. According to measured values of water quality parameters, integer values have been assigned to these parameters (from 1 to 16) to classify surface water quality by indicating the pollution level for each parameter (Issa 2014). OIP classifies surface water quality into five main classes: 0-1 excellent; 1-2 acceptable; 2-4 slightly polluted; 4-8 polluted; and 8-16 heavily polluted (Sharda and Sharma 2013).

Health risk assessment of river samples
Risk assessment of water samples for heavy metals was made using previous methods for two types of exposure: ingestion and dermal absorption (Li and Zhang 2010). The dose received by human pathways is determined by the following two equations that simulate a chemical daily intake (µg/kg/day) by ingestion (CDI ingestion ) and dermal exposure (CDI dermal ) according to US Environmental Protection Agency (US EPA 2004;Wu et al. 2009): where C W is the average concentration of heavy metals in water (µg/L), IR is ingestion rate (in this work it was considered to be 2.5 (L/day), ABS g is gastrointestinal adsorption factor, ABS d is dermal adsorption factor, EF is exposure frequency (in this work 365 days/year was adopted), ED is exposure duration (in this work 65 years was assumed), BW is the average body weight (in this work 75 kg was assumed), AT is averaging time (for non-carcinogens and carcinogens is 25550 days), SA is exposed area of skin (in this work 2800 cm 2 was adopted from (De Miguel et al. 2007)), K p is skin adherence factor (cm/h), ET is the exposure time (in this work 0.5 h/day was assumed), CF is a unit conversion factor, for water equals 1 L/1000 cm 3 .
For non-carcinogen risk, the risk on human life and health was assessed by applying hazard quotients (HQ) for heavy metals and hazard index (HI) (Tripathee et al. 2016;Wang et al. 2017).
where RfD is the reference dose for a toxic element (μg/kg d), the RfD ingestion and RfD dermal for non-carcinogen toxic heavy metals investigated in this work were adapted from (De Miguel et al. 2007;Li and Zhang 2010;US EPA 2004;Wu et al. 2009), while K p values were adapted from (US EPA 2004;Wang et al. 2017). SF value for carcinogenic heavy metal for As was taken from (Giri and Singh 2015).

Statistical analysis
River water samples analysis from the different sampling locations was subjected to several statistical analyses: ANOVA; CM; Friedman test; and CA. These statistical analysis helps to determine spatial differences of physicochemical and heavy metals qualities for Diyala-Sirwan river for pre-and post-discharging of wastewater of Kalar city into the river. ANOVA and Friedman tests were used to recognize the differences in the mean values of tested parameters in river water samples before and after the impact point at different distances. CM identifies the strength of relationships between the investigated heavy metals and water physicochemical parameters. CA was performed to categorize river water samples due to variation of physicochemical and heavy metal parameters based on Ward-algorithmic linkage and Euclidean distance methods. Cluster analysis was conducted using XLSTAT, version 2014 for Excel 2013 software.

River water quality and wastewater discharge impacts
The concentrations of water quality parameters taken at different locations from Diyala-Sirwan river in January 2019 within the study area are shown in Table 1. At each site, 21 heavy metals and physicochemical properties were measured or analyzed for river water samples. Four sampling locations (L1, L2, L3, and L4) have been chosen to demonstrate the impact of wastewater discharge from Kalar city on the water quality of the river. However, from Table 1, it can be seen that EC, TDS, TH, temperature, K, Na, and Mg at the four investigated locations are lower than the maximum admissible limit (MAL). Also, concentrations of heavy metals of As, Ba, Cu, Fe, Ni, Se, and Zn at the four locations are lower than MAL. At impact point sampling location, L2, it can be (11) HI = ∑ HQs seen there is a spike in concentrations for all investigated river water quality parameters of physicochemical and heavy metals. Water samples collected from the location L2 reveal a significant deterioration in water quality, all physicochemical parameters, and heavy metals were at elevated levels, close to or higher than MAL, showing the harmful effect of wastewater discharging without treatment into the river water body.
DO was effectively lowered from 9.56 to 4.5 mg/L at locations L1 and L2, respectively. Turbidity was highly elevated due to wastewater discharge as it can be observed that it was 6.88 NTU at L1 and then was increased to 78.2 NTU at L2. At location L2, pH, Hg, and Mn increased from acceptable levels to be higher than MAL, but pH was returned to be normal at the locations L3 and L4, while Hg and Mn concentrations remained to be higher than MAL at the locations L3 and L4. Additionally, Table 1 shows that Ca and Al concentrations were already elevated in the river before expelling the wastewater of Kalar city, indicating that natural weathering of soils and rocks of river banks might be the source of high concentrations of Ca and Al in the river beside the wastewater of Kalar city. In general, it has been noticed that the water of Diyala-Sirwan river is significantly influenced by Kalar city wastewater discharge especially at the impact point and subsequently at locations after various distances. Discharging untreated wastewater to the river imposes a considerable decline in water properties, suggesting a disturbing in the aquatic life and water quality for drinking, agricultural and urban uses with the study area. Therefore, to some extent, the direct use of Diyala-Sirwan river without treatment at the study area for domestic and agricultural purposes is undesirable. As seen in Table 1, wastewater impact concerned appears to be within an acceptable range, when comparing with some physicochemical rates influenced by wastewater discharges for similar rivers cases in Iraq.
A noticeable increase in heavy metals concentration at/ after impact point of wastewater discharge shows that many heavy metals were originated by anthropogenic activities in Kalar city and then transported in wastewater such as As, Cu, Fe, Hg, Mn, and Zn. Sources of such metals are likely to be tap water and roofs, galvanized material and car washes, the amalgam in teeth for Cu, Zn, and Hg, respectively (Sörme and Lagerkvist 2002).

Analysis of variance
ANOVA at 95% confidence level for the sampling locations showed no significant difference in Diyala-Sirwan river water quality caused by Kalar city wastewater discharge and among waster sampling location, as F value, F critical, and p-value were 0.273, 2.71, and 0.845, respectively.
Friedman test resulted in a ranked table for the investigated variables according is presented as Table 2. The calculated Chi-square was 188.91 which is highly greater than the value of 3.84146 found from the state decision rule list for alpha = 0.05, k = 2, showing that a significant difference is occurring between the two conditions of pre-and postimpact point.

Friedman test analysis
Friedman test analysis is a non-parametric statistical tool that is mostly used to identify the variations in group sets of obtained data. In this work, the difference between the two conditions of Diyala-Sirwan river, pre-and post-impact point of wastewater discharge, has been investigated for alpha equals 0.05. Similar ANOVA hypotheses of null (H 0 ) and alternative hypotheses (H 1 ) were used to double-check the results of ANOVA.

Correlation analysis matrix (CM)
CM gives numbers, if they are closer to -1 or + 1, indicating strong correlation, whereas numbers closer to zero are meaning a weak correlation. In this work, Pearson correlation analysis was performed for river water samples of the post-impact point of wastewater discharge at a significant level of 0.05, as illustrated in Table 3. CM of the 22 physicochemical and heavy metals parameters in river water samples after wastewater discharge shows various strong (higher than 0.7,) positive or negative correlations. Table 3 exhibits many strong positive relationships, between EC, DO, TDS, turbidity, BOD, TH, temperature, High positive relationships between physicochemical parameters and heavy metals suggest that similar inputs sources of chemical pollutants come in urban wastewater discharge into the river (Bastami et al. 2012). Also, strong relationships between heavy metals in the river after the impact point suggest similar input sources and transport ways of these metals (Jiang et al. 2014).
Negative strong correlations between pH and most of the investigated physicochemical and heavy metals parameters in river water refer to the fact that low pH in the river water leads to an increase in heavy metals and chemical parameters solubility (Gonzalez et al. 1990). Lowered concentrations of heavy metals in river water gradually from impact point location were due to the formation of insoluble higher oxides and hydroxides of the heavy metals resulting in their precipitation (Suresh et al. 2012).

Cluster analysis (CA)
CA analysis was conducted to identify similarities among clustered results that were established by this method by presenting internal clusters homogeneity and significant external heterogeneity between clusters. CA in this study was applied based on spatial similarity among river water samples of sampling locations in the study area.
From the results presented in Fig. 2, a dendrogram of cluster analysis result established three distinct clusters depending on the similarity of physicochemical parameters and heavy metals in river water samples for sampling locations. Cluster 1 which represents the location L2, the wastewater discharge impact point, is noticeably different from other clusters, at which a high pollutant input is occurring and leading to a significant deterioration in the water quality of the river.
Cluster 3 was for locations L3 and L4, water samples in this cluster were quite similar and, at the same time, were different from cluster 2 (L1 location), meaning that the river did not fairly overcome the high input of pollutants due to wastewater discharge even though at 10 km from the impact point. It can be assumed that due to the high loud of pollutants that discharged at point L2 the river  Table 3 Pearson correlation matrix for physicochemical parameters and heavy metals in Diyala-Sirwan river after the impact point of wastewater discharge self-purification process takes longer time and distance than 10 km to retain the same water quality of location L1, before impact point, therefore, as shown from CA, the influence of pollutants in the river might stay for long distance from their sources (Cukrov et al. 2008).

Pollution evaluation indices of Diyala-Sirwan River
Calculations of pollution evaluation indices HPI, HEI, C d , PI, and OIP in this study were based on river water quality standards taken from the literature (as illustrated in Table 1). For HPI, HEI, and C d indices, mean concentrations of heavy metals in water samples were implemented for the locations before and after the impact point of wastewater discharge. While for PI and OIP indices, physicochemical and heavy metals parameters dataset were used for the same locations mentioned before. In HPI calculations, ideal concentrations of heavy metals were adapted from literature (Cengiz et al. 2017;Herojeet et al. 2015;Issa and Alshatteri 2018). Herojeet et al. (Herojeet et al. 2015) identified three classes for HPI pollution results: low (< 15); medium (15-30); and high (> 30). While, the classification by HEI results for river pollution are; low (< 1.24), medium (1.24-2.48), and high (> 2.48) as stated by (Khoshnam et al. 2017). C d classification for river pollution are three classes: low (C d < 1), medium (C d = 1-3) and high (C d > 3) (Clesceri et al. 1998). Figure 3 displays values of HPI, HEI, C d , and OIP, at location L2 of the impact point of wastewater discharge, there is a spike in all pollution indices: HPI increased from 89.31 to 314.69; HEI increased from 3.11 to 9.29; Cd increased from -6.89 to -0.71; OIP increased from 3.89 to 4.84. According to HPI and HEI, Diyala-Sirwan river was already polluted and then with Kalar city wastewater discharge, the river has  been more polluted. Cd index indicates that the river water quality in terms of heavy metals pollution was not considerably impacted by the wastewater discharge as the river was slightly polluted before and after the impact point. Regarding OIP, Diyala-Sirwan river condition was changed from slightly polluted to be polluted at location L2. For locations L3 and L4, at 5 and 10 km after the impact point, generally the river could not regain the water quality condition to be similar to location L1 (before impact point). The reason for that might be the river was already carrying considerable concentrations of heavy metals before reaching the impact point and a high amount of pollutants was discharged by Kalar city wastewater into the river, which takes more time and longer distance to be naturally recovered. Table 4 shows the PI results for the most important water quality parameters. Before the impact point, PI values range between 0.001 and 0.761 indicating no pollution effect in the river, except PI for Al was 1.038 meaning that the river was slightly affected by this parameter. After the impact point along with a distance of 10 km, PI values range between 0.004 and 2.687, even though the values of all tested parameters were elevated, but the highest value was for Hg for which the river is considered a moderately polluted concerning the concentration of this heavy metal. Slight pollution effects were also noticed for Al and Ca.
High PI value for Al before the impact point most likely comes from the water treatment plant of Kalar city, in which considerable amounts of aluminum sulfate are used and the resulted treatment sludge expelled without treatment into Diyala-Sirwan river. Otherwise, the elevated PI of Hg and Ca after impact point suggests a release of some industrial effluents into the city wastewater of these two parameters (Hammes et al. 2003;Wagner-Döbler et al. 2000).

Health risk assessment of Diyala-Sirwan River due to wastewater discharge
Three HI classes were proposed for non-carcinogenic effects by (Cui et al. 2015): HI is higher than 1, the risk is significant, HI is between 1 to 0.1, the risk is moderately significant, and HI is less than 0.1, the risk is neglected. From (Fig. 4), it can be seen that the total non-carcinogenic effects represented by HI on human health within the study area for people exposed to Diyala-Sirwan river by different pathways have variable risk potentials for the studied heavy metals. In general, all the heavy metals have HI values less than unity, suggesting that no significant effect is occurring on the people in the area. For non-carcinogenic risk, both As and Hg showed that they are moderately significant at sampling locations L2, L3, and L4. Additionally, Hg was not significant at L1, meaning that the risk was originated due to wastewater discharge in the river, while for As, HI results were increased after the impact point, showing that the wastewater discharge in the river contributed in increasing the risk of Hg existence in river water.
Sources of As in Diyala-Sirwan river and wastewater of Kalar city most probably are mixed from natural origins such as soil and rocks weathering, and from anthropogenic activities, which is mainly for the study area are fossil combustion and application of arsenic-containing fertilizers and pesticides for agricultural purposes (Ungureanu et al. 2015). Figure 4 shows that carcinogenic risk of As in water samples at sampling location was 1.134*10 -4 , and after impact point was ranging from 1.323*10 -4 to 1.89*10 -4 , presenting significant carcinogenic risk. These results show that direct exposure to Diyala-Sirwan river at the study area is unsafe, as these results are fallen within a suggested range, between 10 -4 and 10 -6 , of potential carcinogenic risk (USEPA 2010).
The carcinogenic risk of As, which is determined as an incremental probability in lifetime risk that might cause cancer for any individual as a result of drinking water consumption or through dermal exposure (Iqbal and Shah 2013), might cause many serious kinds of cancer in the human body like cancer of liver, lung, kidney, skin, and bladder (Banerjee et al. 2011). Regarding health risk assessment, it is worth mentioning that the approach applied in this study involves parameters of SF, K p , and RfD might not be exactly compatible with Iraqi conditions. Therefore, the current study is only preliminary evaluation of Diyala-Sirwan river condition after wastewater discharge has occurred, so more studies are required to figure out more precisely the risk arising from mixing Kalar city wastewater with the river water and subsequent effects on people and living creatures in the area.

Conclusions
In this work, different evaluation methods have been performed for physicochemical and heavy metals water quality parameters of Diyala-Sirwan river impacted by wastewater discharge from Kalar city, east of Iraq. Results of pollution indices, HIP, HEI, Cd, Pi, and OIP, indicated that the river was considerably impacted by the wastewater discharge, showing that heavy metals concentrations were remarkably elevated within the study area. Multivariate statistics, CA and CM, revealed that the impact of wastewater discharge on river water quality remained significant in water samples collected from all sampling locations after the wastewater impact point. The health risk assessment showed wastewater discharge into the river was most likely responsible for considerable non-carcinogenic of the investigated heavy metals, and carcinogenic risk of both As and Hg on individuals in the study area. The results obtained in this study are important to build a reliable system for environmental monitoring and measuring pollution of Diyala-Sirwan river.
Funding No funding was received for conducting this study.

Conflict of interest
The authors declare that there is no conflict of interest regarding the publication of this paper.

Human and animal rights
This article does not contain any studies involving human participants or animals performed by any of the authors.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.