1
|
Rajan M, Karunanidhi D, Gurugnanam B, Subramani T. Assessment of groundwater suitability for drinking and irrigation purposes with probable health threats in a semiarid river basin of South India. WATER ENVIRONMENT RESEARCH : A RESEARCH PUBLICATION OF THE WATER ENVIRONMENT FEDERATION 2025; 97:e70011. [PMID: 39904506 DOI: 10.1002/wer.70011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2024] [Revised: 12/18/2024] [Accepted: 12/30/2024] [Indexed: 02/06/2025]
Abstract
In the semiarid river basin of south India, the present study focuses on the appropriateness of water for drinking and irrigation as well as the risks to human health posed by pollutants. A total of 68 groundwater samples were evaluated for irrigation and consumption purposes. With a high electrical conductivity peaking at 3430 μS/cm and an alkaline composition, the groundwater has a high salinity and poor water quality. Durov's figure displays a trend along the dissolution or mixing line and identifies the geochemical facies of groundwater samples. According to water quality indexes, the majority of samples are categorized as unfit for human consumption (26.47%), extremely bad (36.76%), and poor (26.47%). According to elemental concentrations, the data are grouped into three clusters using hierarchical cluster analysis. According to the geographical distribution, nitrate levels are safe over about 320.25 km2 and dangerous over about 121.10 km2, whereas fluoride levels are safe over about 293.92 km2 and dangerous over about 147.43 km2. About 50.65 km2, 14.70% of the samples, fell into the no restriction category for irrigation, indicating acceptable standards. Low sodium levels in soils are indicated by parameters like SAR, %Na, PI, RSC, MR, and KR; SAR values fall into the C2S1, C3S1, and C4S1 categories. According to Doneen's diagram, 70.5% of samples had a PI >75, indicating suitability; the Wilcox diagram classified 22.05% of samples as excellent and 69.11% as good to permissible for irrigation. According to human health risk assessment, 75% of babies, 63% of children, 75% of teens, and 54% of adults have THI values >1 for fluoride. About 45% of newborns, 42% of kids, 45% of teenagers, and 29% of adults are at risk for nitrate. Infants, kids, and teenagers are at the danger. In order to safeguard human health against fluoride and nitrate, the study emphasizes the necessity of efficiently managing groundwater resources, lowering agricultural pollution, and assuring clean drinking water. PRACTITIONER POINTS: In the area, 79.25 km2 has good drinking water quality based on DWQI. Based on IWQI, 70.33 km2 area is recognized as suitable for agricultural practices. Geogenic and anthropogenic activities contribute to fluoride and nitrate pollution in water. Based on THI, infants and children are more prone to fluoride and nitrate contamination.
Collapse
Affiliation(s)
- Meera Rajan
- Department of Civil Engineering, Hindusthan College of Engineering and Technology (Autonomous), Coimbatore, India
| | - D Karunanidhi
- Department of Civil Engineering, Hindusthan College of Engineering and Technology (Autonomous), Coimbatore, India
| | - B Gurugnanam
- Centre for Applied Geology, The Gandhigram Rural Institute Deemed to be University, Gandhigram, Dindigul, India
| | - T Subramani
- Department of Geology & Department of Mining Engineering, College of Engineering Guindy (CEG), Anna University, Chennai, India
| |
Collapse
|
2
|
Bărbulescu A, Barbeș L. Assessing the Efficiency of a Drinking Water Treatment Plant Using Statistical Methods and Quality Indices. TOXICS 2023; 11:988. [PMID: 38133389 PMCID: PMC10747972 DOI: 10.3390/toxics11120988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 11/25/2023] [Accepted: 11/29/2023] [Indexed: 12/23/2023]
Abstract
This study presents the efficiency of a drinking water treatment plant from Constanța, Romania. Individual and aggregated indices are proposed and built using nine water parameters for this aim. The analysis of individual indices permits the detection of the period of malfunctioning of the water treatment plant with respect to various parameters at various sampling points. In contrast, the cumulated indices indicate the overall performance of the treatment plant during the study period, considering all water parameters. It was shown that the outliers significantly impact the values of some indices. Comparisons between the simple average and weighted average indices (built taking into account the importance of each parameter) better reflect the impact on the water quality of some chemical elements that might harm people's health when improperly removed.
Collapse
Affiliation(s)
- Alina Bărbulescu
- Department of Civil Engineering, Transilvania University of Brașov, 5 Turnului Str., 500152 Braşov, Romania;
| | - Lucica Barbeș
- Department of Chemistry and Chemical Engineering, Ovidius University of Constanța, 124 Mamaia Bd., 900112 Constanţa, Romania
- Doctoral School of Biotechnical Systems Engineering, Politehnica University of Bucharest, 313, Splaiul Independenţei, 060042 Bucharest, Romania
| |
Collapse
|
3
|
Chen Y, Hao C, Yang L, Yao L, Gao T, Li J. Toward understanding the interaction of shale gas-water-carbon nexus in Sichuan-Chongqing region based on county-level water security evaluation. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2023; 30:99326-99344. [PMID: 37610545 DOI: 10.1007/s11356-023-29265-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Accepted: 08/07/2023] [Indexed: 08/24/2023]
Abstract
This study develops a comprehensive framework for understanding the interaction of shale gas-water-carbon nexus in Sichuan-Chongqing region. Within this framework, a county-level water security index (WSI) evaluation system is structured. Spatial autocorrelation model and spatial matching degree model are integrated to illustrate the spatial agglomeration characteristics of water security and the water-carbon relationship, respectively. The impacts of shale gas development on water security and carbon emissions are evaluated based on identification of shale well productivity. Results show that about 25.17% of counties with WSI < 0.4 (unsafe), especially in the eastern region. The central cities (such as Chengdu and Neijiang) should take active steps to reach a safety threshold (WSI ≥ 0.6). Population growth can accelerate water security deterioration through uncertainty analysis. Moreover, the spatial matching degree between WSI and carbon emissions in most cities is extremely poor (< 0.5), implying that these cities should optimize their energy structure and promote green transformation. Water used for shale gas extraction can hardly be ignored from a county-scale perspective, especially in Tongliang, Tongnan, and Jianyang. The future shale gas development would pose a threat to the regional climate.
Collapse
Affiliation(s)
- Yizhong Chen
- School of Economics and Management, Hebei University of Technology, Tianjin, 300401, China.
| | - Can Hao
- School of Economics and Management, Hebei University of Technology, Tianjin, 300401, China
| | - Lingzhi Yang
- School of Economics and Management, Hebei University of Technology, Tianjin, 300401, China
| | - Lan Yao
- School of Economics and Management, Hebei University of Technology, Tianjin, 300401, China
| | - Tianyuan Gao
- School of Economics and Management, Hebei University of Technology, Tianjin, 300401, China
| | - Jing Li
- Hebei Key Laboratory of Environmental Change and Ecological Construction, College of Resource and Environment Science, Hebei Normal University, Shijiazhuang, 050024, China
| |
Collapse
|
4
|
Çankaya Ş, Varol M, Bekleyen A. Hydrochemistry, water quality and health risk assessment of streams in Bismil plain, an important agricultural area in southeast Türkiye. ENVIRONMENTAL POLLUTION (BARKING, ESSEX : 1987) 2023; 331:121874. [PMID: 37230170 DOI: 10.1016/j.envpol.2023.121874] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/22/2023] [Revised: 05/03/2023] [Accepted: 05/22/2023] [Indexed: 05/27/2023]
Abstract
In this study, the water quality of Ambar, Kuruçay, Pamuk and Salat streams, which are the important tributaries of the Tigris River in the Bismil Plain (Diyarbakır, Türkiye) was assessed using 19 physicochemical parameters. Except for a few exceptions, all parameters in the water samples taken from the streams were below the drinking water limit values. Kuruçay Stream had significantly higher TOC, Na+, NO3-, NO2-, Cl- and SO42- levels and lower DO levels than other streams (p < 0.05) due to sewage water discharges, animal manure storage areas near the stream and irrigation return flows. In all streams, Ca-HCO3 was dominant water type. Gibbs diagram indicated that rock weathering is the major factor controlling the hydrochemistry of the streams. According to the results of water quality index (WQI), all sampling stations of Ambar, Pamuk and Salat streams and K1 station of Kuruçay Stream had "good" quality water for drinking purposes, while K2 station of Kuruçay Stream had "poor" quality water. Irrigation indices (permeability index, sodium percentage, magnesium hazard, residual sodium carbonate, Kelley's ratio, sodium adsorption ratio and potential salinity) revealed that all water samples taken from the streams were suitable for irrigation. The water samples from Ambar, Pamuk and Salat streams were in the C2S1 (medium salinity and low alkalinity) category, while the samples from Kuruçay Stream were in the C2S1 and C3S1 (high salinity and low alkalinity) categories. Both hazard quotient and hazard index values of NO3-N, NO2-N and F- for children and adults were found below 1, indicating that adverse health effects are not expected from exposure to these contaminants via water ingestion and dermal contact. The findings of this study showed that the water quality status of Kuruçay Stream is worse as it receives large amounts of irrigation return flows compared to other streams.
Collapse
Affiliation(s)
- Şihap Çankaya
- Dicle University, Department of Biology, Faculty of Science, Diyarbakır, Türkiye
| | - Memet Varol
- Malatya Turgut Özal University, Doğanşehir V.K. Vocational School, Aquaculture Department, Malatya, Türkiye.
| | - Aysel Bekleyen
- Dicle University, Department of Biology, Faculty of Science, Diyarbakır, Türkiye
| |
Collapse
|
5
|
Radmanesh M, Ahmadi SH, Sepaskhah AR. Measurement and simulation of irrigation performance in continuous and surge furrow irrigation using WinSRFR and SIRMOD models. Sci Rep 2023; 13:5768. [PMID: 37031260 PMCID: PMC10082809 DOI: 10.1038/s41598-023-32842-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 04/03/2023] [Indexed: 04/10/2023] Open
Abstract
The SIRMOD and WinSRFR models were used to model and assess the irrigation performance under continuous and surge irrigation strategies with two furrow lengths of 70 m and 90 m and stream sizes of 0.4 l/s and 0.6 l/s for each length. According to the normalized root mean squared error (NRMSE) and the relative error (RE), WinSRFR had, on average, excellent accuracy in the continuous and surge irrigation for simulating advance-recession times (NRMSE: 6.15 and 4.24% for advance time, and 2.20 and 5.20% for recession time), infiltrated water depth (NRMSE: 3.37 and 6.38%), and runoff volume (RE: 6.93 and 2.57%), respectively. SIRMOD had also, on average, excellent simulation in the continuous and surge irrigation for advance-recession times (NRMSE: 3.34 and 2.45% for advance time, and 2.28 and 6.41% for recession time), infiltrated water depth (NRMSE: 2.98 and 5.27%), and runoff volume (RE: 5.31 and 17.49%), respectively. The average of irrigation application efficiency (AE), distribution uniformity (DU), deep percolation (DP), and tail-water ratio (TWR) were 61.50, 90.25, 11.75, and 26.75% in continuous irrigation, and 72.03, 94.09, 8.39, and 19.57% in surge irrigation, respectively, which shows that surge irrigation increased AE (irrigation management performance) and DU (irrigation method performance) and reduced DP and TWR compared to continuous irrigation. Moreover, longer furrow lengths increased AE and DP under both irrigation methods, while it decreased TWR and DU. However, increasing the stream size decreased AE and DP and increased TWR under both continuous and surge irrigations. The higher stream size improved DU in continuous irrigation but reduced DU in surge irrigation. The results confirmed that both SIRMOD and WinSRFR are reliable analytical tools to evaluate furrow irrigation strategies for improving irrigation management. In conclusion, this study showed that surface irrigation models could be employed in practice by irrigation engineers and practitioners to design and define the optimized furrow length and stream size in arid and semi-arid areas where efficient and high performance irrigation strategies are required to save water and reduce water loss.
Collapse
Affiliation(s)
- Mojgan Radmanesh
- Water Engineering Department, School of Agriculture, Shiraz University, Shiraz, Iran
| | - Seyed Hamid Ahmadi
- Water Engineering Department, School of Agriculture, Shiraz University, Shiraz, Iran.
- Drought Research Center, Shiraz University, Shiraz, Iran.
| | - Ali Reza Sepaskhah
- Water Engineering Department, School of Agriculture, Shiraz University, Shiraz, Iran
- Drought Research Center, Shiraz University, Shiraz, Iran
| |
Collapse
|
6
|
Tokatlı C, Varol M, Ustaoğlu F. Ecological and health risk assessment and quantitative source apportionment of dissolved metals in ponds used for drinking and irrigation purposes. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2023; 30:52818-52829. [PMID: 36849683 DOI: 10.1007/s11356-023-26078-2] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Accepted: 02/18/2023] [Indexed: 06/18/2023]
Abstract
In this study, dissolved metal levels of 10 different ponds used as irrigation and drinking water sources in the north of Saros Bay (Türkiye) were evaluated using multivariate statistical methods, contamination and ecological risk indices, and absolute principal component score-multiple linear regression (APCS-MLR). The mean levels of metals in the ponds ranged from 0.045 µg/L (Cd) to 127 µg/L (Mn). Pond 7 used for drinking water source had the lowest total metal level. Only Mn levels in two ponds (P1 and P2) slightly exceeded the critical value set by EU Drinking Water Directive. However, the levels of all metals in all ponds were lower than the critical values set for irrigation water and aquatic life. According to the heavy metal pollution index (HPI), five ponds showed low metal pollution in terms of drinking water quality, four ponds showed moderate metal pollution, and one pond (P1) showed moderate to heavy pollution. According to the Nemerow pollution index (NPI) values (0.26-1.82), ponds P1 and P2 showed slight metal pollution, while other ponds showed insignificant metal pollution. Contamination degree (CD) values of ponds varied between 0.95 and 3.33, indicating that all ponds showed low pollution. In terms of irrigation water quality, all ponds showed low or insignificant metal pollution according to the HPI, NPI, and CD values. According to the ecological risk index (ERI) values, metals in all ponds posed low ecological risks for both drinking and irrigation purposes. Factor analysis identified two potential sources: mixed sources and natural sources. The APCS-MLR model results revealed that mixed sources and natural sources contributed 78.99% and 21.01% to dissolved metals in the ponds, respectively. Health risk assessment results indicated that both individual and combined metals in the ponds would not cause non-carcinogenic risks to both adults and children. Similarly, it was found that Cr and As would not cause carcinogenic risks to the residents of the region.
Collapse
Affiliation(s)
- Cem Tokatlı
- İpsala Laboratory Technology Department, Trakya University, Edirne, Turkey
| | - Memet Varol
- Department of Aquaculture, Doğanşehir V.K. Vocational School, Malatya Turgut Özal University, Malatya, Turkey.
| | - Fikret Ustaoğlu
- Biology Department, Faculty of Arts and Science, Giresun University, Giresun, Turkey
| |
Collapse
|
7
|
Yan Z, Li P, Li Z, Xu Y, Zhao C, Cui Z. Effects of land use and slope on water quality at multi-spatial scales: a case study of the Weihe River Basin. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2023; 30:57599-57616. [PMID: 36971941 DOI: 10.1007/s11356-023-25956-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Accepted: 02/11/2023] [Indexed: 05/10/2023]
Abstract
Exploring the impact of land use and slope on basin water quality can effectively contribute to the protection of the latter at the landscape level. This research concentrates on the Weihe River Basin (WRB). Water samples were collected from 40 sites within the WRB in April and October 2021. A quantitative analysis of the relationship between integrated landscape pattern (land use type, landscape configuration, slope) and basin water quality at the sub-basin, riparian zone, and river scales was conducted based on multiple linear regression analysis (MLR) and redundancy analysis (RDA). The correlation between water quality variables and land use was higher in the dry season than in the wet season. The riparian scale was the best spatial scale model to explain the relationship between land use and water quality. Agricultural and urban lands had a strong correlation with water quality, which was most affected by land use area and morphological indicators. In addition, the greater the area and aggregation of forest land and grassland, the better the water quality, while urban land presented larger areas with poorer water quality. The influence of steeper slopes on water quality was more remarkable than that of plains at the sub-basin scale, while the impact of flatter areas was greater at the riparian zone scale. The results indicated the importance of multiple time-space scales to reveal the complex relationship between land use and water quality. We suggest that watershed water quality management should focus on multi-scale landscape planning measures.
Collapse
Affiliation(s)
- Zixuan Yan
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China
- State Key Laboratory of National Forestry Administration On Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an University of Technology, Xi'an, 710048, Shaanxi, China
| | - Peng Li
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China.
- State Key Laboratory of National Forestry Administration On Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an University of Technology, Xi'an, 710048, Shaanxi, China.
| | - Zhanbin Li
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China
- State Key Laboratory of National Forestry Administration On Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an University of Technology, Xi'an, 710048, Shaanxi, China
| | - Yaotao Xu
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China
- State Key Laboratory of National Forestry Administration On Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an University of Technology, Xi'an, 710048, Shaanxi, China
| | - Chenxu Zhao
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China
| | - Zhiwei Cui
- State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi'an University of Technology, No.5, South Jinhua Road, Xi'an, 710048, Shaanxi, China
- State Key Laboratory of National Forestry Administration On Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an University of Technology, Xi'an, 710048, Shaanxi, China
| |
Collapse
|
8
|
Varol M, Tokatlı C. Evaluation of the water quality of a highly polluted stream with water quality indices and health risk assessment methods. CHEMOSPHERE 2023; 311:137096. [PMID: 36334749 DOI: 10.1016/j.chemosphere.2022.137096] [Citation(s) in RCA: 48] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Revised: 10/21/2022] [Accepted: 10/30/2022] [Indexed: 06/16/2023]
Abstract
The water quality of Çorlu Stream, located in the Thrace region of Türkiye, and exposed to intense industrial pressure, was evaluated by monitoring 10 toxic metals and 13 other water quality variables in the dry and wet seasons of 2021. Seven different water quality indices were applied to determine the pollution level at the sampling stations in the stream. In addition, human health risks from exposure to toxic metals in stream water via ingestion and dermal contact were evaluated. The results showed that the water quality at stations S2 and S3 of Çorlu Stream receiving domestic and industrial discharges are seriously polluted by NH4-N, PO4-P, COD, BOD5 and suspended solids according to surface water quality standards. In addition, these stations were highly polluted and had poor water quality according to the results of the water quality indices. The average Cr level at station S3 exceeded the permissible levels set for the protection of aquatic life due to effluent discharges from the leather factories. Considering the results of the health risk assessment methods, non-carcinogenic risks from ingestion of combined metals in stream water can be expected at station S3 for both children and adults and at station S2 for children. Also, it was estimated that Cr and As at station S3 may cause carcinogenic health risks for residents.
Collapse
Affiliation(s)
- Memet Varol
- Department of Aquaculture, Doğanşehir Vahap Küçük Vocational School, Malatya Turgut Özal University, Turkey.
| | - Cem Tokatlı
- Laboratory Technology Department, Trakya University, İpsala, Edirne, Turkey
| |
Collapse
|
9
|
Varol M, Karakaya G, Alpaslan K. Water quality assessment of the Karasu River (Turkey) using various indices, multivariate statistics and APCS-MLR model. CHEMOSPHERE 2022; 308:136415. [PMID: 36099988 DOI: 10.1016/j.chemosphere.2022.136415] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 08/22/2022] [Accepted: 09/07/2022] [Indexed: 06/15/2023]
Abstract
Determining the water quality status of a river and accurately identifying potential pollution sources threatening the river are pillars in effective control of pollution and sustainable water management. In this study, water quality indices, multivariate statistics and absolute principal component score-multiple linear regression (APCS-MLR) were applied to evaluate the water quality of the Karasu River, the main tributary of the Euphrates River (Turkey). For this, 19 water quality variables were monitored monthly at eight stations along the river during one year. Based on the mean dissolved oxygen (DO), electrical conductivity (EC), nitrate-nitrogen (NO3-N), orthophosphate-phosphorus (PO4-P), total phosphorus (TP), ammonium-nitrogen (NH4-N), chemical oxygen demand (COD) and total nitrogen (TN) levels, most stations of the river had "very good" water status according to surface water quality criteria. Spatial cluster analysis (CA) divided eight stations into three regions as clean region, moderate clean region and very clean region. The mean values of Nutrient Pollution Index indicated that the river was "no polluted". Similarly, Water Quality Index and Organic Pollution Index values indicated that the river water quality was between "good" and "excellent". A minimum water quality index (WQImin) consisted of ten crucial parameters was not significantly different with the WQI based on all the 17 parameters. Discriminant analysis (DA) results showed that water temperature (WT), EC, chlorophyll-a (Chl-a), potassium (K), calcium (Ca), NO3-N and COD are the variables responsible for temporal changes, while WT, total dissolves solids (TDS), Chl-a, K, magnesium (Mg), Ca, NH4-N and COD are the variables responsible for spatial changes in the river water quality. Principal component analysis/factor analysis (PCA/FA) identified four potential sources, including anthropogenic, natural, seasonal and phytoplankton. Source apportionment in the APCS-MLR model revealed that seasonal and anthropogenic sources contributed 35.2% and 25.5% to river water quality parameters, respectively, followed by phytoplankton (21.4%) and natural sources (17.9%).
Collapse
Affiliation(s)
- Memet Varol
- Department of Aquaculture, Doğanşehir V.K. Vocational School, Malatya Turgut Özal University, Turkey.
| | | | | |
Collapse
|