JWSS - Journal of Water and Soil Science

One of the factors leading to the contamination of water resources is human activity, producing waste materials. In this paper, the effects of contamination on the water quality of Taleghan River, was simulated using of Qual2k model and the seasonal changes were evaluated. The qualified data collected during two months, August (as the dry season) and February (as the wet season), were used in the modelling. The results showed that the dissolved oxygen change was in the range of 4.5-6.52 mg/L in August. However, it changed between 4.8-5.3 mg/l in February and this reduction in the wet season was due to the run off deposition and the seepages of farmland near the river. Furthermore, BOD in the wet season and the dry season changed by 6-31 and 10-26 mg/l, respectively. These changes were due to the sewage dilution in flow during the wet season. Evaluation of the pH values in wet and dry seasons also showed that water of the river was more alkaline in the wet season, which is due to the effect of non-point resources or the fertilizer entrance after farm land bleaching by rain. Evaluation of Taleghan river EC variation also showed these changes were not significant during the wet and dry seasons. Water temperature was altered by 3-10⁰C and 19-25⁰C, respectively, in February and August.

In the present study, we used 27 precipitation average monthly data from synoptic, climatologic, rain-guage and evaporative stations located in Zayandeh-Rud river basin for the period of 1970-2014. Before interpolating, the missing data in the time series of each station was reconstructed by the normal ratio method. Also, for the data quality control, the Dickey-Fuller and Shapiro-Wilk tests were used to check the data stationarity and normality. Then, these data were interpolated by six interpolation methods including   Inverse Distance Weighting, Natural Neighbor, Tension Spline, Regularized Spline, Ordinary Kriging and Universal Kriging; then each method was evaluated using the cross-validation technique with MAE, MBE and RMSE indices. The results showed that among the spatial interpolation methods, Natural Neighbor method with MAE of 0.24 had the best performance for interpolating precipitation among all of the methods. Also, among Ordinary Kriging, Universal Kriging, Spline and Inverse Distance Weighting methods, respectively, Exponential Kriging with MAE 0.54, Quadratic Drift Kriging with MAE of 0.5, Tension Spline with the MAE of 0.54 and Inverse Distance Weighting with the power of 4 with MAE of 0.57 had the least error compared to other IDW methods.

In general, rivers are one of the best and most accessible water resources at the disposal of mankind.  So, given the effect of the force of water and changes on the  flow patterns and consequently, on river morphology changes, the analysis of the flow in the river is important and necessary to organize projects, flood control and water supply structures downstream. In this study, by using numerical models CCHE2D hydraulic conditions Dimeh River Bridge between Oregon Bridge Sudjan was investigated. CCHE Model is a mathematical model capable of simulating the flow patterns and sediment transport in rivers and canals laboratory network. The numerical model in 1998, based on the calculations by the National Centre for Water Science and Engineering, University of Mississippi (NCCHE), was developed and has been applied in many research projects related to water engineering. At the outset, the input data required model provides and numerical model was implemented. In the next step, the results of the model were calibrated and validated using field data measurements; eventually, they were extracted and their model results were compared; it was confirmed that CCHE model could still simulate the flow pattern.

In this research, conjunctive and integrated operation of surface and ground water resources of Behbahan plain (Maroon dam's reservoir and existing wells, respectively) was investigated. Simulation of allocation of water demands in this basin was performed by four scenarios, using WEAP software: 1) current conditions (M1), 2) reference scenario for the next 16 years (M2), 3) land development scenario (M3), and 4) optimal scenario (M4). The optimal scenario was performed with multi-purpose linear programming. Based on the results, drinking water demands was satisfied completely in all scenarios. Under the scenario of current conditions, all agricultural demands, except the traditional rights, supplied more than 50% in the low-flow months. In the reference scenario, water supply for agricultural demands in some months was less than 100% and even in June and July, the water supply for North and South Irrigation networks of Behbehan plain was less than 10%. In the land development scenario, agricultural demands of all irrigation networks, except Ramhormoz network, satisfied more than 90% in all months. The optimal scenario performed better than other scenarios for minimum Maroon River flow and volume of storage in the reservoir. Comparison of the four scenarios in satisfying the environmental needs also revealed that the optimal scenario performed better than the other three scenarios in the spring months. However, it provided less than 100% of water needs in the whole year. Comparison of the four scenarios also showed that the first two scenarios had the highest reliability percent in the Jayzan-Fajr, South Behbahan and North Behbahan Irrigation Networks and traditional water rights. Frequency of storage-time-probability from the storage volume in the optimal scenario also showed that maximum storage lifetime of the lasting storage volume was 558 million m3 (which was equal to half of the volume of Maroon dam’s reservoir) with the highest probability (60%).

After Hamidieh Diversion Dam near the city of Hamidieh, Karkheh River is divided into two streams known as Hufel and Nissan. At the lower flow rates, Nissan makes up a greater share than Hufel due to the steeper slope of the former. This study attempted to construct a hydraulic structure to appropriately divide water flow in Hufel. In a laboratory experiment, a flume with a 90-degree bend was used at Islamic Azad University of Ahvaz. Various experiments were conducted at different widths and heights. Furthermore, this model was simulated through CCHE2D, the results of which were compared against those of physical and mathematical models. The results indicated that the weir height increased the deviation flow percentage to the Hufel stream due to rising water level. Moreover, the deviation flow percentage to Hufel was declined as the weir width was increased due to falling water level. At Hufel, the installation of rectangular weir in different dimensions yielded the minimum of 34.3% and the maximum of 61.5% increase in the flow rate. In the normal mode without any weirs installed, however, there would be an increase in the flow rate, as compared to the mode where a weir has been installed. This can be associated with the flow controlled by the weir. On average, the deviation flow rate was increased by 2.8% in the weir mode and 7.7% in the weir-less one. An increase in the Froude number from 0.21 to 0.38 led to a lower average deviation flow rate by 19.3%. Moreover, the results of the simulation through CCHE2D were demonstrated to be largely similar to those of physical model experiments. However, an increase in the Froude number did not lead to a decline in the deviation flow rate (i.e. it remained constant). This trend was inconsistent with the results of the physical model.

Occurrence of heavy floods in rivers causes a lot of damages and losses. In this research, to highlight the river-training reaches in 10.9 km of Babolrud River, first, using topographic map of the area, the Tin layer was created in GIS software. Then, using the HEC-GeoRAS extension, the main route and cross sections of the river were prepared and introduced to the HEC-RAS model. River discharge with return periods of 2 to 200 years was calculated. Flow analysis in the agricultural and urban areas was performed and the areas which need training measures such as flood-retaining walls and levees were specified and designed. The structural design, stability control and sliding was performed based on the Standard No. 518, using RetainWall software, and design and control of levee stability was performed based on the Standard No. 214, using GeoStudio software. The cost of project implementation was estimated based on the Price List of 2017. The output of HEC-RAS software showed that height of flood-retaining wall in 3 urban reaches ranged between 1.73 to 2.8 m and in 5 agricultural reaches ranged between 1.46 to 2.25 m. It was concluded that the overall cost of levee implementation is about 9.01 billion Rials, of gravity concrete flood-retaining wall is about 9.26 billion Rials and of concrete cantilever inverted T shape flood-retaining wall is about 10.05 billion Rials. Thus, using flood-retaining levee is the most economical option.

Hydrodynamic models proposed for simulating flow hydraulic in rivers assume the flow in one direction and simulate the hydraulic parameters based on the one-dimensional Saint-Venant equations. In this research, a two-dimensional HEC-RAS model was used to simulate the flow in the Karun River, between Mollasani and Farsiat stations. Geographic information system (GIS) and river cross sections were used to prepare the altitude map using the satellite image of the study area. Modeling results in river bends showed that the maximum velocity occurred in the outer bend, which coincided with the flow mechanism in the bends. Based on the results, grid type and density have little effect on flow depth modeling. However, the characteristics of the mesh used had a great influence on the velocity distribution, so that the regular high-density mesh had the best accuracy in simulating the flow velocity. Statistical analysis showed that the RMSE for the flow discharge and flow depth were 17.95 m³/s and 0.05 m, respectively. In addition, the Nash–Sutcliffe efficiency index was calculated to be above 0.9 for the discharge and flow depth, which could be considered as a desirable value.

Nowadays, the prediction of river discharge is one of the important issues in hydrology and water resources; the results of daily river discharge pattern could be used in the management of water resources and hydraulic structures and flood prediction. In this research, Gene Expression Programming (GEP), parametric Linear Regression (LR), parametric Nonlinear Regression (NLR) and non-parametric K- Nearest Neighbor (K-NN) were used to predict the average daily discharge of Karun River in Mollasani hydrometric station for the statistical period of 1967-2017. Different combinations of the recorded data were used as the input pattern to predict the mean daily river discharge. The obtained esults  indicated that GEP, with R²= 0.827, RMSE= 59.45 and MAE= 26.64, had a  better performance, as compared to LR, NLR and K-NN methods, at the  validation stage for daily Karun River discharge prediction with 5-day lag, at the Mollasani station. Also, the performance of the models in the maximum discharge prediction showed that all models underestimated the flow discharge in most cases. 

Investigating and understanding river change issues is one of the important factors in sediment hydraulic sciences and river engineering. These studies can be done with the help of physical, mathematical models, or both, but due to financial and time constraints, mathematical models are more general and often used. In this study, the GSTARS model was used to investigate erosion and sedimentation and select the most appropriate function in 12.5 km in length from the Talar river in Mazandaran Province. Simulation using the 55 sections taken in 2006, the daily flow data of the hydrometric station of the Shirgah, located at the beginning of the rich and characteristics of the river sediment, was done. The calibration and validation of the model with cross sections taken in 2012 showed that Yang's sediment transport equation has the highest correlation with reality and can be used to predict river change. The amount of sediment depleted from the case study using the Yang equation is estimated at 8590 tons per year. Also, the study of longitudinal profiles of the river with different sediment transfer functions showed that the study reach at the end range has an erosion trend and is not capable of sand and gravel mining.

Reducing surface water resources and successive droughts and consequently excessive use of groundwater resources, especially for agricultural purposes, have caused irreparable damage to the natural resources of the country. In the meantime, knowing the status of the water balance of the plain can help to effective management of water resources in the region. Samalqan plain is located in a semi-arid climate in North Khorasan Province. Since the surface water resources for water supply are not very reliable, so, the main source of water supply in the region is well. Due to the existence of rivers in the plain, the low thickness of the alluvium, groundwater level fluctuations, and the high uncertainty in the calculation of hydrodynamic coefficients, the need for careful hydrogeological studies and determining the role of each parameter affecting groundwater is necessary. This study was conducted to simulate the Samalqan aquifer and analysis of water balance for the years 2003 to 2013 using the MODFLOW model. To identify the groundwater recharge rate, this component was estimated by the SWAT model. Calibration and validation of the model with an error of 1.1% and 1.2%, respectively, indicated that an appropriate estimation between the simulated and observed heads. Assessment of the groundwater hydrograph in the observation wells showed that the groundwater level in most places has many monthly and seasonal fluctuations. After drawing the potential lines of the plain, the inputs and outputs were identified, and using the reserve volume changes, the water balance was determined. The results showed that the water balance of the plain was negative and the reservoir deficit was estimated at 9.14 million cubic meters. Therefore, this model can be used to predict the future situation of aquifer and the management of water resources in the region.

Bed shear stress is one of the most important hydraulic parameters to determine the amount of bed and suspended load and the bed and bank scouring in rivers. Bed shear stress depends on bedforms (ripples, dunes, and anti-dunes) in alluvial rivers. In this study, the effect of artificial ripple bedforms on bed shear stress has been investigated. Two types of uniform granulation with average sizes (d50) of 0.51 and 2.18 mm were used to roughen the surface of the artificial ripples. The bedform length and height were 20 and 4 cm, respectively. The angles of its upstream and downstream to the horizon were selected equal to 16.4 and 32 degrees, respectively. Different flow rates (Q= 10, 15, 20, 25, and 30 l/s) and different bed slopes (S= 0, 0.0001, 0.0005, 0.001, and 0.0015) were examined. The results showed that by increasing the particle size on the bed surface, total shear stress (t_b ), grain-related bed-shear stress (t^¢_b ), and form-related bed-shear stress ( t^²_b )  increase. The value of t_b , t^¢_b , and t^²_b in bed form roughened by sediment size of 2.18 mm were, on average, 22.38, 30.86, and 22.3% more than the bed form roughened by sediment size of 0.51 mm, respectively.

The use of structures has economic and safety advantages compared to other energy-consuming structures. In this research, to investigate the effect of the length of the sill of the flip bucket spillway on the scour downstream, experiments were conducted in a rectangular laboratory flume made of Plexiglas. The scouring downstream of the flip bucket spillway was investigated using a flip bucket spillway with four relative sill lengths and four threshold angles at four current intensities in the channel in this research. The results of this research showed that by the increase in the length of the sill in the flip bucket spillway, the energy consumption in the spillway increased and the scour depth downstream decreased. Also, increasing the relative length of the sill by 70% at the sill angle of 45 degrees, the scouring depth is reduced by about 88%. Also, a relationship was presented to determine the maximum depth of relative scour, and the correlation coefficient of the results obtained from this equation with the laboratory results is about 0.92.

The rapid and complex movement of sediments in rivers and coastal areas with highly erosive and unsteady flows presents river engineers with numerous problems in the geomorphology of alluvial rivers. Accurately predicting these complex processes in the water-sediment system (a multiphase, dense, granular flow system) is still a major challenge for mesh-based models. Due to the ability of meshless Lagrangian methods to model large deformations and discontinuities, meshless Lagrangian methods can provide a unique way to deal with this complexity. In the current research, the capabilities of the weakly compressibility moving particle semi-implicit (WC-MPS) model in soil-fluid interaction modeling are developed to enable the modeling of sediment transport and erosion effects behind coastal walls. In this method, granular material is considered a non-Newtonian and viscoplastic fluid. The 𝜇(I) rheological model has been used to predict the non-Newtonian behavior of the granular phase. To verify the application of the present model in simulating the interaction of liquid and solid phases, first, the widely used problem of dam break on an erodible bed was modeled. The NRMSE model was calculated to be approximately 6%, which indicates the efficiency and accuracy of the target model in this problem. At the end, the scouring of coastal walls was simulated by the WC-MPS method using 𝜇(I) rheology model. Investigations show that the processes related to erosion and scouring can be well modeled using the current Lagrangian method. The numerical results show excellent agreement with the laboratory measurements. It should be noted that the mean error of the mentioned model is estimated to be 10%.

Drought represents one of the most critical natural disasters, exerting profound impacts on agriculture, society, the economy, and water resources. Various indices are used to monitor drought and its effects. This study aims to monitor drought in the Zayandeh-Rud Basin using the Standardized Precipitation Index (SPI), the Standardized Precipitation-Evapotranspiration Index (SPEI), the Evaporative Demand Drought Index (EDDI), the Palmer Drought Severity Index (PDSI), and the Reconnaissance Drought Index (RDI). All these indices are based on potential evapotranspiration, incorporating parameters such as precipitation, temperature, relative humidity, wind speed, and sunshine duration. These five indices were calculated and evaluated during the statistical period of 1993–2023 for meteorological stations in Isfahan, East Isfahan, Kabootarabad, Daran, Shahreza, Najafabad, and Mobarakeh. After calculating the indices and using spatial zoning maps, the studied stations were compared in terms of these indices. The continuity of dry and wet periods, as well as the intensity of droughts and wet spells, was analyzed. Subsequently, drought intensities during different years in these stations were ranked using the TOPSIS model based on factors such as precipitation, potential evapotranspiration, and station elevation. The results showed that in stations with a dry climate (such as Isfahan, East Isfahan, and Shahreza), drought occurrences (as indicated by higher rankings) have been consecutive over multiple years. Comparing the performance of the indices in the studied stations using spatial zoning maps revealed that the intensity of droughts and wet spells in regions with dry and semi-dry climates was not very significant. However, in areas with humid climates, the fluctuations in drought and wet spell intensities were quite substantial. The findings indicate that the PDSI and EDDI indices are more suitable for evaluating drought in dry climates.