Search published articles


Showing 152 results for Analysis

M. Shamaeizadeh, S. Soltani,
Volume 18, Issue 70 (3-2015)
Abstract

Hydrologic drought which usually affects wide regions can be studied through Low flow index. In this study, to predict hydrologic drought in North Karoon watershed, 14 stations with suitable and long enough duration data were recorded in the 1387-88 water year. Then 13 physiographic and climatic characteristics of the chosen stations were used to perform homogeneity test for cluster analysis. 7 day low flow series were calculated in each station and according to chi-square and Kolomogragh smirnov tests and parameter, 2 parameter gamma distribution was selected as the best regional distribution for this region. Therefore, a seven day low flow index was estimated using FREQ for 5,10,20,50,100 return periods. Regional analysis was performed using a multiple regression method. Moreover, flow duration curves were delineated to obtain Q95 index. Then, zoning maps for Q95، Q7,2 ،Q7,10, Q7,100 were prepared. The results of regional analysis indicated that the averages of height and slope were the two most effective parameters in low flow in this watershed. The investigation of zoning maps showed that southeastern part of this watershed experiences severe droughts compared with other parts.


S. Dowlatabadi, S. M. A. Zomorodian,
Volume 19, Issue 71 (6-2015)
Abstract

One of the most essential and appropriate groundwater model components is accurate information of the recharge values among input data often introduced to the model as the percentage of rainfall of aquifers. The recharge values are influenced by many temporal and spatial factors. Firoozabad plain is one of the suitable plains for agriculture in the Fars province in which utilization of groundwater resources has been banned since 23 September 2002, due to the declining water level and negative balance. The main purpose of this study was to estimate the recharge values of groundwater aquifer by using SWAT in the MODFLOW model. Firstly, surface water was simulated via SWAT model, and sensitivity analysis, calibration, validation and uncertainty analysis of results were performed by SWAT-CUP software. After extraction of aquifer recharge values from the calibrated model, the groundwater of basin was simulated via MODFLOW model in both steady and unsteady conditions. Following the model calibration, the hydrodynamic coefficients of plain were determined and sensitivity of model was checked in terms of hydraulic conductivity and discharge rate of pumping wells. As for the confidence, the model was revalidated, which proved in simulating the behavior of the aquifer very well.


S. Akhavan, A. Jodi Hameze Abad,
Volume 19, Issue 72 (8-2015)
Abstract

Urmia Lake, located in north-west of Iran, has been exposed to various threats such as drought, construction of dams, land use changes and increased global temperature. Due to the importance of Urmia Lake, it is feasible to conduct different kinds of studies to identify the problems of its watershed. The main objective of this study was to evaluate SWAT program’s ability to simulate runoff in Urmia Lake watershed with an area of 52000 km2. The model was run for the 1980-1997 period. Calibration and validation periods were from 1980 to 1991 and from 1992 to 1997, respectively. The results of calibration for 10% and 85% of hydrometric stations were very good and suitable, respectively. Also, validation results for 25% and 45% of hydrometric stations were very good and suitable, respectively. These results show the high ability of SWAT model to simulate discharges in Urmia Lake watershed. Moreover, some factors influencing inflow to the lake in recent years were evaluated. The outcomes revealed that recent changes (dam cconstructions, climate change and land use change) in the watershed have caused inflow volume to the lake to decrease by 80%. So, if natural management conditions had prevailed in the watershed, the Lake’s conditions would have been much better.
H. Beigi Harchegani, S. S. Heshmati,
Volume 19, Issue 72 (8-2015)
Abstract

The aim of this paper is to adapt a water quality index for individual samples and to compare the results with that of the original GIS-based approach. Thirteen water quality parameters observed in 97 wells from the Shahrekord aquifer were used. In GIS-based method, quality parameters maps are difference-normalized, ranked and GWQI map is drawn. In derived method, observations from individual wells were separately and similarly treated to obtain WQI for each well. Both GWQI maps displayed similar trends and were highly correlated (R=0.91). While the minimum and mean GWQI for both methods were identical (respectively 81 and 84) the derived method estimated the maximum GWQI slightly lower (7%) and showed up to 6% difference in water quality class coverage. Overall, the derived method GWQI is more correlated with observations and performs better than the GIS-based method, and therefore, can be used for determining the overall quality of individual water samples and without the requirement of samples being spatially distributed.


M. Goodarzi, J. Abedi Koupai, M. Heidarpour, H. R. Safavi,
Volume 19, Issue 73 (11-2015)
Abstract

Due to the time and space changes of hydrological events in the arid and semi-arid regions, recharge measurement in these areas is very difficult. Hence, groundwater recharge is a complicated phenomenon for which there is not a fixed method to determine. The aim of this research was to develop a method for estimation of groundwater recharge based on a hybrid method. In this study, a hybrid method for calculating recharge was presented by combining empirical methods with a mathematical model, MODFLOW, and AHP analysis. The results showed that the most important parameters affecting groundwater recharge are soil properties, unsaturated thickness, land cover, land slope, irrigation and precipitation, from which the soil properties and precipitation are most important. The results showed that the overall impact of small changes in precipitation and temperature significantly affect the groundwater recharge, and heavy soils are much more sensitive to these changes than light soils. By changing 10% precipitation, the recharge rate is changed between 16% and 77% and by changing 1ºC temperature, the recharge rate is changed between 6% and 42%. Also, results showed that precipitation and evapotranspiration changes in four months including December, January, February and March had significant effects on annual recharge rate. Using the results of this research, the vulnerable areas of the plain, appropriate places and time for artificial recharge could be identified. Overall, the results of this study can be useful in various aspects of groundwater management.


M. Khoshravesh, J. Abedi-Koupai, E. Nikzad-Tehrani,
Volume 19, Issue 74 (1-2016)
Abstract

During the past few decades, the southern part of the Caspian Sea has more frequently experienced extreme climatic events such as drought and flood. Trend analysis of hydro-climatic variables was conducted using non-parametric Mann-Kendall test and regression test for Neka basin in the north of Iran.       Trends of precipitation and stream flow characteristics including maximum flow, mean flow and low flow indices were analyzed at the annual, seasonal and monthly time scales from 1358 to 1391 (34 years). Results showed a general decrease in annual and winter precipitation and decrease in daily maximum precipitation, with an increased trend in daily maximum precipitation of spring season. A decreasing trend was observed in 7-day low flow in summer for all sub-basins. Annual and monthly mean flows specifically in winter in all sub-basins decreased, but annual maximum flow increased from upstream to downstream. Land use changes showed that deforestation and urbanization increased during 34 years in the mid and downstream sub-basins. The analysis showed that low flow indices and mean flows are strictly sensitive to climate change. Overall, from hydrological perspective, these results indicate that the study region is getting dryer and facing more severe drought events. The results of this study can predict future droughts to make better decisions for irrigation planning and management of water resources.


L. Kashi Zenouzi, Sh. Banej Shafiee, A. A. Jafari,
Volume 20, Issue 76 (8-2016)
Abstract

In this study the effect of temperature, evaporation or evapotranspiration, precipitation, hillside direction and altitudinal classes, texture and acidity of soil on organic carbon content in the depths of 15 and 45 cm were evaluated. Paired t-test results showed that there is a significant difference between measured parameters in two soil depths. After preparing required data and processing them, outlier's data were removed. Then, base maps for each of the information layers were prepared by Arc GIS9.3 software and all relatd information fit together by overlapping them. Pearson correlation between environmental factors and soil organic carbon values were calculated and it was found that in the depth of 15 cm, the correlation between soil organic carbon values and two environmental factors including temperature and altitude were significant at the level 0.01. As well the results of statistical analysis by using principal component analysis (PCA) method showed that the factors temperature, evaporation (1%), and silt and clay (5%) have had a significant effect on the amount of soil organic carbon. The first, second, and third axes with eigenvalues of 98/4, 78/3 and 92/1, respectively, explained the values 0.33, 0.25, and 0.13 % of correlation between organic carbon and environmental data.


Sh. Ahmadi Doabi, M. Afyuni, H. Khademi, M. Karami,
Volume 20, Issue 76 (8-2016)
Abstract

Heavy metals in dust can directly enter to the human body through ingestion and inhalation. They can pollute the water and soil resources via atmospheric precipitation and accumulate in the plant tissues and then enter human body through water and food. This research aimed to study the heavy metals concentration in dust in Kermanshah province and to identify their sources. 49 samples of dust were collected in the cities of Kermanshah, Songhor, Gilangharb, Ghasre-Shirin, Sahneh, Sarpolzahab, Kangavar, Paveh and Javanrood during the spring 2013. The concentration of Zn, Cu, Ni, Cr, Mn and Fe were determined using an atomic absorption spectrometer following the sample extraction with a mixture of HCL and HNO3 (3:1 ratio). The average concentrations of Zn, Cu, Ni, Cr, Mn and Fe were 182.3, 48.6, 115.3, 73.9, 428.1 and 23161 mg kg-1, respectively. Correlation, cluster and principal component analyses were used to identify probable natural and anthropogenic sources of contaminants, and the enrichment factor was used to identify probable effects of human activity on the concentration of heavy metals. The results indicated that metal concentrations, except for Fe and Mn, were higher in comparison with the world soils. Zn, Cu, Ni and Cr are mainly of anthropogenic origin, while Fe and Mn are mainly of natural origin. Zn and Cu are mainly of traffic sources and partly of industrial sources, and Ni and Cr are mainly derived from industrial sources, combustion processes, combined with traffic sources. The analysis of EF revealed moderate enrichment for Mn and Cr, and significant enrichment for Zn, Cu and Ni. Based on the results of this study, more attention should be paid to identifying and controlling the sources of contaminants such as heavy metals in dust in order to prevent their associated pollution.


H. Saghi,
Volume 21, Issue 1 (6-2017)
Abstract

The water supply networks have always been of significance to researchers as a hydraulic system of transferring and distributing water. The pressure gradient is the main reason of water transfer in networks, and in case of non-standard pressure increase, the undesirable phenomenon of leakage occurs in the network. Leakage in urban water distribution networks causes water waste and enormous financial losses. Therefore, there exists the need to manage and minimize the amount of leakage. In this study, a water distribution network is modeled using the potentialities of hydraulic analysis model, the EPANET 2.0, and, by presenting a new model, the leakage location is recognized. In order to do this, we, firstly, entered all the network parameters into the software. Then, the network was analyzed supposing the non-occurrence of leakage and the amounts of nodal pressures were measured. Moreover, the nodal pressures were estimated by creating a hypothetical leakage in one of the network nodes and analyzing the network. Finally, the position of leakage was determined by defining the leakage index and comparing it in various nodes. The results show that the suggested method is efficiently capable of predicting the leakage position in the network so it can safely replace other methods, especially destructive methods used in recognizing the leakage position in the network.

A. H. Boali, H. Bashari, R. Jafari, M. Soleimani,
Volume 21, Issue 2 (8-2017)
Abstract

Appropriate criteria and methods are required to assess desertification potential in various ecosystems. This paper aimed to assess desertification levels in Segzi plain located in east part of Isfahan, with a focus on soil quality criteria used in MEDALUS model. Bayesian Belief Networks (BBNs) were also used to convert MEDALUS model into a predictive, cause and effects model. Soil samples were collected from 17 soil profiles in all land units and some of their characteristics such as texture, soluble sodium and chlorine, organic material, Sodium Absorption Ratio (SAR), Electrical Conductivity (EC) and CaSo4 of all soil samples were determined in soil laboratory. The effects of measured soil quality indicators on desertification intensity levels were assessed using sensitivity and scenario analysis in BBNs. Results showed that the used integrated method can appropriately accommodate uncertainty in the desertification assessments approaches created as a result of the influence of different soil characteristics on desertification. According to the results of MEDALUS model, 28.28 % and 71.72 % of the study area were classified as poor and moderate areas in terms of soil quality respectively. Sensitivity analysis by both models showed that soil organic matter, SAR and EC were identified as the most important edaphic variables responsible for desertification in the study area. Evaluating the effects of various management practices on these variables can assist managers to achieve sound management strategies for controlling desertification.
 


F. Jalilian, B. Behmanesh, M. Mohammad Esmaeili, P. Gholami,
Volume 21, Issue 2 (8-2017)
Abstract

In this study, different indices of vegetation cover variations and different physicochemical properties of soil in three treatments of flood spreading, enclosure and grazing (control) were investigated and compared in in the region of Peshert in Mazandaran province. In order to measure different soil characteristics, 18 soil samples (six withdrawals at any treatment) from a depth of zero to 30 cm were taken from the desired treatments. In order to investigate different vegetation indices, a total of 90 plots (nine transects of 100 m) were run using systematic random sampling in the studied treatments and the necessary measurements were done (30 plots at any treatment). Then, in each of these plots, canopy coverage percentage was determined separately for each species and to evaluate and assess the diversity and richness in all three treatments, Shannon-Wiener and Simpson diversity indices and Menhink and Margalef richness indices were used. Finally, the data obtained from both sections of soil and vegetation in three studied treatments were compared and analyzed using one-way ANOVA and Duncan test. The results showed that floodwater spreading and enclosure significantly increased the percentage of sand and total Nitrogen, and significantly reduced the percentage of silt and potassium compared to control treatment. Also, percentage of clay and organic matter, soil pH levels, conductivity and soil phosphorus showed no significant differences in the treatments under study. The results of variance analysis of various indices of diversity, richness and species evenness showed that all indicators had significant responses in three treatments and the highest diversity and species richness were observed in flood spreading and enclosure treatments. Due to changes in soil properties and vegetation in flood spreading and enclosure treatments compared to the control treatment, it can be stated that operations of floodwater spreading and enclosure in the studied region has had positive effect on modification of soil texture, increasing the permeability of the soil and ultimately improvement of the vegetation.


M. Khoshravesh, M. Valizadeh,
Volume 21, Issue 2 (8-2017)
Abstract

Construction of a dam in the upstream of Tajan basin, as one of the human activities in order to provide drinking water, hydropower, agriculture, flood control or other purposes, can be effective on other sectors such as water resources. The purpose of this study is investigating the impact of Rajaei dam on groundwater resources in Sari-Neka plain during 26 years period (1985-2011) using geostatistical method, parametric and non-parametric tests. The results of statistical tests showed that groundwater table variation during 26 years has experienced insignificant upward trend. Quality factors such as EC, TDS, TH, Cl, SO4, Mg, Na and SAR had significant downward trend at 5% confidence level for Pearson and Kendall test. The Quality changes of groundwater resources in the region had good condition due to construction of dam. The results of spatial analysis for the study area using Kriging interpolation method in the three years for 1985, 1999, 2011 showed that in 1999 (the year of exploitation of the dam) most of the quality factors had reached the highest level from 1985 until the dam operation (1999) and then decreased to 2011. The results of this study showed that, agricultural activities after dam construction, operation of irrigation system and drainage network design of Rajaei dams, have raised the water table in northern part of the region and consequently had a significant impact on the quality of water resources.
 


A. H. Boali, R. Jafari, H. Bashari,
Volume 21, Issue 3 (11-2017)
Abstract

This paper aimed to assess the severity of desertification in Segzi plain located in the eastern part of Isfahan city, focusing on groundwater quality criteria used in MEDALUS model. Bayesian Belief networks (BBNs) were also used to convert MEDALUS model into a predictive, cause and effects model. Different techniques such as Kriging and IDW were applied to water quality data of 12 groundwater wells to map continuous variations of the CL, SAR, EC, TDS, pH and decline in water table indices in GIS environment. The effects of measured water quality indicators on desertification severity levels were assessed using sensitivity and scenario analysis in BBNs model. According to the results of the MEDALUS, the desertification of the study area was classified as severe class due to its low quality of groundwater. Sensitivity analysis by the both models showed that decline in waater table, water chloride content and electrical conductivity were the most important parameters responsible for desertification in the region from ground water condition standpoint. The determination coefficient between the outputs of the MEDALUS and BBNs models (R2>0.63) indicated that the results of both models were significantly correlated (α=5 %). These results indicate that the application of BBNs model in desertification assessment can appropriately accommodate the uncertainty of desertification methods and can help managers to make better decision for upcoming land management projects.
 


K. Nosrati, M. Majdi,
Volume 21, Issue 4 (2-2018)
Abstract

The soil pollution especially in urban soils is projected to increase drastically and its effects on chemical cycles are yet to be known. Approaches to measure air and water quality are well established, but urban soil quality assessment has received little attention. Soil quality assessment can help as a way to better understand the pollution increase outcomes in urban environments and to establish approaches and integrated soil quality assessment protocols in urban planning and landscape management. Considering lack of information in urban soil quality of Iran, the objective of this study was to assess soil quality under urban land use effect using minimum data set in western part of Tehran. In view of this, 56 soil samples were collected in three land use types of agricultural, parks and urban landscapes, and vacant urban lots and 12 physicochemical properties were measured. The results of analysis of variance (one-way ANOVA) showed that under influence of the land use types, organic carbon, total nitrogen, lime, bulk density and sodium have significant differences. The factor analysis was used to select minimum data set and the results showed that two factors with eigenvalues more than one, explaining more than 68% of total variance, have the most loading factors on organic carbon and sodium. Finally, soil quality indicator (SQI) was determined and compared in different land use types. The results showed that SQI has significant difference in urban land use types and the least soil quality is related to vacant urban lots.
 


S. Esmailian, A. Talebi, M. Esmailian,
Volume 22, Issue 1 (6-2018)
Abstract

This research was aimed to simulate and prioritize the effective factors on water erosion using USLE-M in the system dynamic model. In this integrated model, by using the system dynamic simulation software (Vensim), all variables and factors involved in erosion and soil loss were considered according to the USLE-M model. After model implementation, the estimated values and observations were compared and then sensitivity analysis was done to determine the sensitive parameters. Then, calibration was performed on the sensitive parameters. This study found that that the results of the model were acceptable for soil erosion simulation due to considering all the effective factors in soil erosion. The results of the sensitivity analysis also indicated high model sensitivity to the slope and vegetation cover in high and low slopes, respectively. By investigating the changes in various parameters such as vegetation cover and slope on erosion, the optimal vegetation cover with 67 and 40% slope, was estimated to be 20 and 60%, respectively.

M. Nouri, M. Homaee, M. Bannayan,
Volume 22, Issue 1 (6-2018)
Abstract

In this study, the trends of changes of the standardized precipitation index in a 12-month timescale (SPI-12) and seasonal and annual precipitation were investigated in 21 humid and semi-arid stations of Iran during the 1976-2014 time period. After removing the serial correlation of some series, the trend of precipitation and SPI-12 was detected using the Mann-Kendall nonparametric trend test. The results revealed that the trends of annual precipitation had been declining in all stations over the past 39 years.  The seasonal precipitation trend in winter, spring, autumn and summer was downward in approximately 90, 95, 47 and 37% of the studied stations, respectively. In addition, the descending trend of wintertime precipitation was significant in Sanandaj, Khoy, Urmia, Hamedan, Mashhad, Torbat-e-heydarieh, Nozheh and Qazvin. Also, the temporal trend of SPI-12 was decreasing in all surveyed stations except Shahrekord. Furthermore, SPI-12 showed a significant downward trend only in Sanandaj and Fasa. Moreover, the most severe meteorological drought occurred in the period 1999-2000, in Ramsar, Urmia and Hamedan, and in the period 2008-2009, in Tabriz, Sanandaj, Shiraz, Fasa, Qazvin, Mashhad, Torbat-e-heydarieh, Shahrekord, Gorgan and Kermanshah stations. Overall, the results of this study indicated that the trend of precipitation in most studied sites, particularly in semi-arid parts of the northeast and southwest of Iran, has changed due to the severe and long metrological drought that has occurred in the recent decade (2005-2015).
 


S. P. Mousavi, M. A. Asghar Mokhtari, Y. Khosravi, A. Rafiee, R. Hoseinzade,
Volume 22, Issue 2 (9-2018)
Abstract

In this study, the distribution of heavy metals pollution including arsenic, antimony, nickel, copper, cadmium, cobalt, bismuth, lead and zinc in the stream sediments of Zarshuran- Aghdarreh area was investigated by using statistical techniques and the geometric integration of each sample basin. For this purpose, the degree of pollution in 154 stream sediment samples was analyzed and the distribution maps for enrichment factors were prepared by using a combination technique, pixel estimation, and statistical and geostatistical methods. The results of calculating the enrichment factors indicated that the higher enrichment was related to arsenic, antimony, bismuth, cadmium and lead. Furthermore, the concentration of zinc, copper, lead, arsenic, antimony, cadmium and bismuth in the stream sediment samples was higher than the global average. Application of the principal component analysis on the data led to the recognition of 9 main components for the dataset; the first 5 were components with eigen values greater than 1 and a cumulative percentage more than 85%. Arsenic, antimony, cadmium, lead and zinc in the first component, cobalt in the second component, bismuth in the third component, copper in the fourth component and nickel in the fifth component had the highest values.

M. Jahan, B. Amiri,
Volume 22, Issue 3 (11-2018)
Abstract

Factor analysis is one of the multivariate statistical techniques that considers the interrelationships between apparently irrelevant variables and helps researchers to find the hidden reasons for the occurrence of an event. In order to evaluate the effects of different irrigation levels and humic acid foliar application and identify the factors affecting water use efficiencies of sesame (Sesamum indicum L.), maize (Zea mays L.) and common bean (Phaseolus vulgaris L.), a split plots experiment based on RCBD design with three replications was conducted during the 2014-15 growing season, at the Research Farm of Ferdowsi University of Mashhad, Iran. Irrigation levels (50 and 100% of water requirement) and foliar application and non-application of humic acid were assigned to main and sub plots, respectively. The results showed that in sesame, the highest seed yield and biological yield were obtained from 100% of water requirement and humic acid spraying treatment. In maize, humic acid spraying under condition of supplying 50% of water requirement increased seed weight per plant, plant height, and leaf area index and soil pH In bean, the highest seed weight per plant, plant height, leaf area index, crop growth rate and soil phosphorous content were observed in the treatment of 100% of water requirement and humic acid spraying. Factor analysis results also showed that in sesame, the variables of seed yield, biological yield, seed weight per plant, plant height, leaf area index, crop growth rate, soil phosphorous and water use efficiency were assigned to the first factor and the variables of soil nitrogen, soil pH and EC were assigned to the second one. In maize, seed yield was assigned in the same group with the variables of biological yield, leaf area index, crop growth rate, soil phosphorous, EC and pH and water use efficiency; in bean, this was with the variables of seed yield and water use efficiency. In general, the research results revealed that identifying the effective variables in each factor and those logical nominations according to Eco physiological knowledge can lead to the direct management of effective variables with regard to associated factor, thereby leading to water efficiency improvement.

M. A. Geranmehr, M. R. Chamani, K. Asghari,
Volume 22, Issue 3 (11-2018)
Abstract

A water distribution network (WDN) may not be able to satisfy all required demands when it’s in the pressure deficit mode or under over-loaded demand conditions. Analysis of the network in this mode requires pressure dependent analysis (PDA). Unlike demand driven analysis (DDA), PDA needs an extra equation for every node to relate the nodal demand and the nodal pressure; so it should be solved with the other network’s equations simultaneously. In this paper, based on the Particle Swarm Optimization (PSO) algorithm, a decision support system has been developed by using MATLAB and EPANET for PDA simulation in WDNs. A four-loop network selected from the literature was analyzed using different scenarios and different pressure dependent functions presented by the previous investigations. The results showed that the proposed model (PSO-PDA) was as accurate as the previous ones and provided better convergence. The results of the nodes’ pressure and discharge also indicated minor differences obtained by different PDA functions. However, the differences between the results of PDA and DDA were considerable.

A. Shahbaee Kotenaee, M. Foroumadi, O. Ahmadi,
Volume 22, Issue 3 (11-2018)
Abstract

One of the major issues in the contemporary world is climate change. The behavior and characteristics of parameters affecting climate change can cause them to be seen and hidden. As one of the effective ways to detect overt and covert behaviors for periodic climatic data series, Spectral analysis can be used. It is the analysis of each of the wavelengths series, making this behavior clear. Accordingly, the present study was an attempt to use the method of spectral analysis, data cycles in the minimum temperature, maximum temperature and precipitation in Ramsar station (located in the western regions of Mazandaran province) an nd Babolsar (located in the central parts of this province) in a period from1961 to 2014. For this purpose, temperature and precipitation data were obtained from these stations; MATLAB software environment and the environment for the software were logged for each of the variable in the stations. The results revealed that the minimum temperature at both stations had significant cycles, with the return period being 2 to 5 years; Remote Link could be fit into the cycle parameters such as NAO, AO and ENSO. Analysis of the period gram showed cycles 8 and 5/13-year-old and 5-year-old period in Ramsar and Babolsar. During the rainy cycles, the difference between the two stations and the difference in the geographical position affected systems, and rain accounted for the difference in speed dual-zone climate indicator for Remote Link.


Page 6 from 8     

© 2024 CC BY-NC 4.0 | JWSS - Isfahan University of Technology

Designed & Developed by : Yektaweb