JWSS - Journal of Water and Soil Science

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.

The river flow prediction is a key aspect of hydrology that plays a significant role in water resources management, flood risk reduction, and agricultural planning. This study simulates the monthly flow of the Razavar River, located in western Iran, using an extreme learning machine (ELM) model enhanced by the Whale (WOA) Optimization Algorithm and Grasshopper Optimization Algorithm (GOA) metaheuristic optimization algorithms. The data used include river flow, precipitation, evaporation, and temperature, which were collected for 10 years with a monthly time step and normalized in the numerical range of zero to one. 80% of the data is used for training, and the remaining 20% for model evaluation. The performance of the models is measured with the statistical indices RMSE, NSE, and R². First, the basic ELM model is developed using the trial-and-error method to adjust the weights between the hidden and output layers. Then, the WOA and GOA algorithms are used to optimize the weights. The results show that the basic ELM model performs worse than the optimized models (Train: RMSE=0.1427, NSE=0.7795, R²=0.7911, Test: RMSE=0.1406, NSE=0.7811, R2=0.7916). While the WOA-ELM and GOA-ELM models provide similar results, the WOA-ELM model shows better performance in complex conditions (Train: RMSE=0.1215, NSE=0.7869, R2=0.7932, Test: RMSE=0.1165, NSE=0.7872, R2=0.7933). The results of this research show that meta-heuristic optimization algorithms play an important role in improving the performance of river flow prediction models due to their ability to search comprehensively and avoid getting stuck in local optima. The findings of this study emphasize the importance of applying these techniques in water resources management and sustainable planning and will pave the way for future research in this area.

One of the most significant hydraulic issues in determining the opening of river bridges is the lack of flow choking due to a reduction in the width of the flood passage. In this paper, determining the required opening for flow passage at a bridge location has been investigated using the concept of specific energy, one-dimensional, and three-dimensional flow modeling. First, the maximum encroachment of the embankments on the sides of the bridge in the river has been determined in such a way that it does not change the flow situation upstream of the bridge, using the concept of specific energy. The dimensions obtained for the bridge opening have been simulated numerically in two one-dimensional and three-dimensional models, and the flow condition at the bridge site and upstream has been evaluated and compared. The results showed that the one-dimensional numerical model predicts, on average, 67 percent higher amount of afflux than the three-dimensional model, while the maximum shear stress obtained from the one-dimensional model is, on average, 33 percent lower than that of the three-dimensional model. The effect of the bridge skewness on the amount of afflux and other hydraulic parameters of the flow, including bed shear stress and maximum velocity, has also been investigated using a three-dimensional model. The afflux was obtained at a 19.2 percent of normal depth at a skew of 40 degrees.

The Helmand River is the main source of water for the Sistan Plain and underpins the provision of key ecosystem services, including water supply, agricultural production, forage provision, soil protection, and wildlife habitat. The operation of the Kamal Khan Dam and the diversion of river flow toward Goud-e Zereh have reduced downstream water availability and accelerated land degradation in Sistan. This study aimed to assess land-use and land-cover changes associated with Helmand River water availability before and after the operation of the dam. Land-use/land-cover maps for 2019, 2020, 2024, 2025, and 2026 were produced using Landsat imagery, supervised classification, and the Random Forest algorithm, and their accuracy was evaluated using field-based reference samples. The results showed that, after dam operation, the Hamoun-e Hirmand water body completely dried up, whereas the water-covered area around Kamal Khan Dam and in the Goud-e Zereh region increased from 3,039 ha in 2019 to 62,794 ha in 2026. In the Sistan Plain, agricultural land reached its maximum extent of 69,025 ha in 2020 but declined to 5,776 ha in 2025. In addition, the area of saline lands in the Sistan Plain and Hamoun increased from a minimum of 110,414 ha in 2020 to a maximum of 274,708 ha in 2026. These changes were accompanied by a reduction in rangelands, an expansion of barren lands, increased soil salinity, and greater soil susceptibility to wind erosion. The findings indicate that the diversion of water toward Goud-e Zereh, together with associated declines in agricultural land, rangelands, and soil moisture and the expansion of saline and barren lands, has weakened the ecosystem services provided by the Helmand River.