نشریه علوم آب و خاک (علوم و فنون کشاورزی و منابع طبیعی)
سال بیست و ششم شماره 4 (پیاپی 102، زمستان 1401)

This study determines the amount of irrigation water saved as a result of the subsidy policy to adapt from flood to drip irrigation. We developed a positive mathematical programming model (PMP) to evaluate the effect of economic incentives on farmers’ decisions to choose the type of irrigation technology, cropping pattern, and "water use" and "water consumption" in rural Garkan Shomali district, which is part of the Najafabad aquifer. We collected data through farm surveys, desk research, and expert interviews. The results showed that a reduction in the financial costs of converting flood irrigation into drip irrigation can lead to farmers investing in this technology. In the current water allocation scenario, the subsidy policy increases the water consumption of drip-irrigated crops by 28%, of which 19% is non-consumed water before subsidy payment and the rest is related to the reduction of furrow-irrigated lands. Also, under non-volumetric water delivery conditions, the operating costs reduce and the net income of the farms increases because of the increase in efficiency and the development of the area under cultivation, which increases water consumption while the water use is constant. In the volumetric water delivery scenario, with the increase in subsidies, the net income of the farms will increase without developing the area under cultivation and only because of the increased yield. Therefore, subsidy policy increases irrigation efficiency at both the farm and regional levels and is an effective tool for dealing with drought conditions.

The present study aimed to investigate and compare the laboratory results of energy dissipation and length of vertical Drops equipped with horizontal Screens with the results of standard stilling basins of type one, two, three and four simple vertical Drops. For this purpose, 64 different experiments were performed on vertical Drops equipped with a horizontal Screen at relative distances of 0, 0.25, 0.5, and 0.75 from the edge of Drops, with a porosity of 40 and 50% of the Screen and a height of 20 cm .The results showed that in all experiments and at a constant flow, increasing the distance of the Screen from the edge of Drops does not have much affect the energy dissipation of the current. On average, the downstream energy dissipation for the present study has increased by more than 20% compared to the simple vertical Drop, which can be an excellent alternative to the downstream stilling basin. Among the models of the present study, the most significant reduction in the relative length of the Drops was achieved by the vertical Drops model with a horizontal Screen with a relative distance of 0.75. On average, when using horizontal Screen at four relative distances from the edge of Drops, the relative length of the Drops is reduced by more than 73% compared to the vertical Drops equipped with a standard stilling basin.

The relationship between the concentration of heavy metals and physicochemical factors was studied in the fine-grained sediments of the tidal section of the Mehran delta where mangrove trees have grown. Surface sediments of the tidal zone of the Mehran delta were sampled. The grain size distribution, calcium carbonate, organic matter contents, pH/Eh, and heavy metals concentration were determined in the mud fraction of the sediments. The presence of high calcium carbonate in sediments is an effective factor in the abundance of Mn, Cd, and Pb. Clay and organic matter as adsorptive have contributed to the high concentration of Zn, Cu, while Fe oxy-hydroxides have affected the concentration of Co, Cr, and Ni. According to the index of Enrichment Factor (EF), the average enrichment of the elements in the studied sediments from high to low order is Cr, Ni, Cd, Zn, Mn, Pb, Co, and Cu. According to the Pollution Load Index (PLI), none of the delta sediments including sediments from inside and outside of the mangrove forestare classified as polluted to the heavy metals. According to the Igeo Index, most samples, and only Cd, Cr, and Ni show slight pollution in some samples. Accordingly, there is no current threat of contamination of potentially toxic elements of natural and anthropogenic origins to the mangrove environment and Mehran delta.

Precipitation is one of the most significant climatic parameters; its distribution and values in different areas is the result of complex linear and nonlinear relationships between atmospheric elements-climatic processes and the spatial structure of the earth's surface environment. Classification of data and placing them in small and homogeneous zones can be effective in improving the understanding of these complex relationships and their results. In the present study, zoning and analyzing the distribution of rainfall in Iran concerning environmental factors was performed using the annual precipitation data of 3423 synoptic, climatological, and gauge stations in the country during the period from 1961 to 2015 and the altitude, slope, aspect, and station density data. After standardization and preparation of the data matrix, the optimal number of clusters was determined and the data set was entered into the neural-fuzzy network model (ANFIS-FCM). The results showed that the values of R2  and MAE  indices were 0.76 and 0.23, respectively which indicate the appropriate accuracy of the model. It was also found that in the four output zones of the model, environmental factors have a high impact on the spatial distribution of precipitation. In the first and third zones, the combination of high altitude and slope factors along with geographical proximity to precipitation systems has caused the average annual rainfall in these zones to be 318 and 181 mm, respectively. The mean annual rainfall has decreased to about 100 mm by the weakening of the role of environmental factors in the second and fourth clusters.

Changes in nutrient concentrations of soil can specify optimal management of manure and prevent environmental and water resources pollution. The present study was conducted with the objective of changing macronutrients concentrations of Nitrogen, Phosphorus, and Potassium with amendments application of polyvinyl acetate, bean residual, and a combination of polyvinyl acetate + bean residual for time periods of one, two, and four months. The results showed that the application of soil amendments had various effects on changing Nitrogen, Phosphorus, and Potassium. The maximum amount of Nitrogen related to the treatment of bean residual at the time period of four months before simulation (with a rate of 44.62 percent) and minimum amount of nitrogen related to Polyvinyl acetate treatment at the time period of one month (with a rate of -1.92 percent). The minimum rate of Phosphorus was measured at the treatment of bean residual at the time period of one month before simulation (with a rate of 0.95 percent). The maximum amount of Potassium related to the treatment of Polyvinyl acetate at the time period of four months before simulation (with a rate of 189.35 percent) and the minimum amount of Potassium related to the combination of bean residual + Polyvinyl acetate at the time period of one month after simulation (with a rate of 40.66 percent). Therefore, the application of amendments has various effects on changing soil macronutrients at different time periods.

The pattern of root distribution in layered soils is one of the significant issues in the calculations of soil water and irrigation management and planning. The objective of this study was to determine the pattern of root distribution of soybean in layered soils and its effect on water uptake. The research was conducted in a completely randomized design with 15 treatments consisting of three different textures of soil (light, heavy, and medium) in four replications. The pattern of root distribution was monitored by the sampling of columns at the end of the growth period of the soybean. It was observed that the presence of the layer with medium texture has led to better plant development and growth after comparing the treatments in terms of plant growth. In general, root length density decreased with increasing soil depth, except in cases where there were different layers of soil, and root length density takes place in the following order: root length density in layers with medium texture≥ heavy texture≥ light texture. The rate of root water uptake rate was highest in the sandy layers, intermediate in clay, and lowest in loamy texture. Also, the rate of root water uptake rate increased significantly with increasing depth regardless of treatments. It can be concluded that the pattern of root distribution and plant growth is significantly affected by soil texture and its stratification.

In the present study, the flow rate in flues containing lateral semi-cylinders (SMBF) was simulated and estimated under free and submerged conditions using back vector machine models (SVM), spin multivariate adaptive regression (MARS), and multilayer artificial neural network (MLPNN) model. In free flow mode, the dimensionless parameters extracted from the dimensional analysis include the ratio of upstream flow to throat width and contraction ratio (throat width to channel width), and in the submerged state, in addition to these two parameters, the depth-to-throat width, and bottom-depth parameters upstream depth were used as input and the two-dimensional form of flow rate was used as the output of the models. The results showed that in free flow mode in the validation stage, the MARS model with statistical indices of R2 = 0.985, RMSE = 0.008, MAPE = 0.87%, and the SVM model with statistical indices of  R2 = 0.971, RMSE = 0.0012, MAPE =1.376%, and MLPNN model with statistical indices of R2 = 0.973,  RMSE = 0.011, MAPE = 1.304% have modeled and predicted the flow rate. In the submerged state, the statistical indices of the developed MARS model were R2 = 0.978, RMSE = 0.018, MAPE = 3.6%, and the statistical indices of the SVM model were R2 = 0.988, RMSE = 0.014, 2%. MAPE = 4, and the statistical indicators of the MLPNN model were R2 = 0.966, RMSE = 0.022, and MAPE = 5.7%. In the development of SVM and MLPNN models, radial kernel and hyperbolic tangent functions were used, respectively.

Class A pan evaporation method as one of the most common methods for reference evapotranspiration (ET0) estimation has been widely used in the world due to its simplicity, relatively low cost, and ability to estimate daily ET. In this study, the performance of 8 empirical methods consisting of Allen and Pruitt (1991), Cuenca (1989), Snyder (1992), modified Snyder, Pereira, et al. (1995), Orang (1998), Raghuwanshi and Wallender (1998), and FAO/56 were analyzed to estimate class A pan coefficient and ET0 at Fasa synoptic station located in Fars province. The calculated pan evaporation coefficients from the above equations were compared with measured pan evaporation coefficients which were obtained from the ratio of evapotranspiration calculated by the FAO-Penman-Monteith method to the rate of evaporation from the pan. The results showed that all empirical methods did not predict pan coefficient values well (R2 < 0.3 and NRMSE > 0.25). The comparison results between ET0 from empirical methods and ET0 obtained from FAO-Penman–Monteith indicated that the FAO/56 method had the best performance (R2 = 0.72 and NRMSE = 0.3). To increase the accuracy of empirical pan coefficient equations, these equations were modified with eight years (2007-2015) of meteorological data from the Fasa synoptic station and validated using two years of independent data (2015-2017). The results showed that the accuracy of all empirical models was improved and the Cuenca equation with NRMSE = 0.16 and R2= 0.63 was selected as the best equation for pan coefficient estimation and ET0 (R2 =0.85; NRMSE =0.18) in Fasa region. The sensitivity analysis revealed that the estimated pan coefficient is more sensitive to wind speed, followed by relative humidity, fetch distance, the slope of the saturation vapor pressure curve, sunshine hours, and air pressure. According to statistical results and sensitivity analysis, an equation was expanded for the Fasa region and other areas with the same climate.

The present study was carried out in the Mahyar region of Esfahan Province to determine optimum drip tape spacing for the wheat crops on a silty clay loam soil respecting grain yield as well as yield components, water use efficiency, and variations in the salinity within the soil profile. The experiment was performed for three years from 2017 to 2019 with a randomized complete block design with three replicates and four treatments. The treatments consisted of three tape spacings (A) at 45, (B) at 60, (C) at 75 cm, and the Control (D) was irrigated with the basin method. The same volume of irrigated water was applied to the drip treatments, A, B, and C in every irrigation interval, whereas for treatment D, the local farmers’ practice was followed. Based on the results from compound variance analysis, the treatment effect on both grain yield and biological yield, and on water use efficiency and harvest index was significant at 1% and 5% level of confidence, respectively. The mean water use efficiency in treatments A, B, C, and D was measured as 0.79, 0.79, 0.73, and 0.78 kg m-3; thus, treatments A, B, and D outperformed treatment C. A comparison between the salinity of the soil profile at the beginning and the end of the growing season revealed that the basin irrigation method was more effective on salt leaching than the drip tape system. The results of this study indicated that concerning water use efficiency and crop yield, drip tapes spaced at 45 or 60 cm outperformed those which were 75 cm apart. On the other hand, the work required for irrigation system installation as well as the amount of drip tape residues left on the field at the end of the growing season is larger for tapes spaced at 45 cm compared to those which are 60 cm apart. This will have a significant impact on farmers’ budgets and environmental issues. Therefore, it is recommended to lay the tapes 60 cm apart for the irrigation of wheat crops on silty clay loam soils.

In this research, the trend of spatial changes in extreme indices of temperature related to the health and agriculture sectors such as the number of frost days, number of summer days, number of icing days, number of tropical nights, growing season length, diurnal temperature range, cold spell duration index, and warm spell duration index were investigated for 54 synoptic stations throughout Iran for observational (1976-2005) and future (2025-2054) periods. Daily maximum and minimum temperature data of three regional climate models namely, CCSM4, MPI-ESM-MR, and NORESM1-ME from the CORDEX project under RCP4.5 and RCP8.5 scenarios were downscaled for each station using a developed multiscale bias correction method. Then, trends and changes of extreme temperature indices were investigated using Mann-Kendall and Sen’s trend line slope methods. The results indicated that the warm indices such as the number of summer days and tropical nights indices have had a positive trend at most stations in both observational and future periods. In contrast, cold indices like the number of frost days have had a decreasing trend in most stations. The results of cold and warm spell duration indices showed that most stations have had no trend for both periods. The growing season length has increased in more than 60% of stations (45% having a significant trend) mainly located in the northern, northwestern, and western regions of the country. Based on the results, it can be concluded that without considering thoughtful climate adaptation measures, some parts of the country may face health risks and limited habitability and agriculture in the future.

Reservoir dams have had problems despite all the benefits for humans. one of the most important issues is exposing a large amount of water in contact with the air causing a large amount of water to evaporate. Using chemical methods including heavy alcohols is one of the evaporation suppression methods. In this study, three emulsions of octadecanol, hexadecanol, and a combination of octadecanol, and hexadecanol along with Brij-35 and two physical methods of the canopy and floating balls were used to evaluate the performance of different emulsions. A one-way analysis of variance was applied to compare the mean of evaporation in different chemical and physical methods and a two-way analysis of variance was performed to investigate the main and interaction effects of different meteorological parameters on the value of evaporation. The mean comparison of the evaporation in different methods showed that the two physical methods of the canopy and floating balls had better performance than the chemical methods, and the octadecanol was more efficient than the two other chemical methods. The results of one-way ANOVA showed that among the chemical methods, the octadecanol had no significant difference with floating balls at a 99% probability level (P <0.01). Two-way ANOVA indicated that air temperature and relative humidity had the greatest effect on evaporation. Examination of the effect of different levels of meteorological parameters on the performance of evaporation reduction methods showed that at low temperatures, octadecanol had poor performance than the two physical methods but with increasing temperature, its performance improved. In addition, this monolayer had a suitable performance at low wind speeds compared to physical methods. By increasing wind speed, its performance is severely affected and its efficiency decreases. So, at temperatures above 37° C, an increase in wind speed from 3.5 m/s to above 8.7 m/s has increased evaporation by more than 50%. The effects of monolayers and other evaporation suppression methods on the quality characteristics of the water including dissolved oxygen are significant and should be investigated in future research.

Trend analysis of stream flow provides practical information for better management of water resources on the eve of climate change. Therefore, the present study investigated river flow variations during three decades as well as projections of future discharge in the Gorganrood watershed. The Man-Kendall method has been used to detect the trend and methods of Pettitt, SNHT, and Buishand to identify points of a sudden change in discharge time series in 8 stations of Aq Qala, Galikesh, Gonbad, Haji Ghoshan, Nodeh, Ramyan, Sadgorgan, and Tamar. The Mann-Kendall trend test showed the existence of a significant negative trend (flow reduction) on a daily and annual scale in all stations. Monthly, the strongest negative trend in Aq Qala, Galikesh, Gonbad, Haji Ghoshan, and Ramyan stations was related to July, but in Nodeh and Tamar stations, it was related to August and February, respectively. A decreasing trend was observed in all stations on a seasonal scale, but this trend was not significant in some seasons. The results of the analysis of change points in discharge showed that the change points in the data used in this study are more of a decreasing and in some cases incremental type and some stations, no change points have been identified at all. Therefore, the number of decreasing changes in the studied hydrometric stations is significantly higher than the incremental changes and is more visible from 1993 to 1997 and 2005-2007 in three and four stations, respectively. Also, the most incremental changes among the stations are related to the Aq Qala station in 2017 with a flow rate of 234 cubic meters per second. Investigation of the flow of the basin in the past decades showed significant monotonic and abrupt changes which are mostly toward decreasing the basin’s discharge. The downward trend in discharge values at different time scales for all hydrometric stations of the Gorganrood watershed, which will be more severe in the future due to global climate change, and increasing the region's water needs for various future use due to population growth and the expansion of industries can also be considered as a serious warning for policymakers, planners, and local managers to prevent a possible water crisis in the region in the future with proper planning.

The land use change as well as changes in climatic parameters such as temperature increase affect many natural processes such as soil erosion and sediment production, floods, and degradation of physical and chemical properties of soil. Therefore, it is necessary to pay attention to different aspects of the effect of these changes in studies and macro decisions of the country. In the present study, the SWAT conceptual model was used to test and analyze the existing scenarios in the Marvast basin. After calibrating the model, the two scenarios were tested. The first scenario is in the field of agricultural management and conversion of gardens to agricultural lands and the second scenario is a 0.5-degree increase in temperature by assuming other conditions are constant. The calibration and validation results of the model with the Nash-Sutcliffe test showed 0.66 and 0.68 respectively, which indicate the acceptable performance of the model in the study area. Then, the results of using two scenarios of land use change and heating, especially in recent years showed the effect of 30 percent of the climate scenario on the increase of flooding in the basin. The scenario of changing the use of garden lands to agriculture in two cases of 20% and 50% change of use of 10% and 12% was added to the flooding of the basin. The results indicate that in similar areas of the study area which is located in a dry climate zone, a possible increase in temperature can have a significant effect on flooding in the basin. However, the indirect impact of the human factor in increasing greenhouse gases and flooding in the basin should not be ignored.

One of the most significant environmental crises in arid, semi-arid, sub-humid, and even humid regions is the destructive phenomenon of desertification and in the arid and semi-arid regions is wind erosion. These problems exist in large areas of Iran and it is necessary to use an environmentally friendly and economic method to solve this problem. In this study, calcium bentonite clay was used for the first time in Iran and perhaps in the worlds in the critical region of Sajzi, which covers an area of 65 hectares. Experiments were performed on the crusts after one year of mulching with bentonite clay. The results showed that wind erosion has a negative and significant correlation with the mean weight diameter and geometric weight diameter of aggregate, aggregates with diameters greater than 0.25 mm, shear strength, and penetration resistance. On the other hand, the results of the permeability test using double-ring and by three models (Kostiakov, Horton, and Philip) showed that the lowest mean square error (SSE) and the highest coefficient of determination (R2) belonged to the Kostiakov model in the mulch-applied and control samples. This result indicated the superiority of the Kostiakov model compared to Horton and Philip's models. Wind erosion intensity was also measured in situ using a portable wind tunnel at 20 points in the Sajzi region. The findings showed that mulch application has controlled more than 95% of soil erosion.

In this study, the decade variability of frequency and severity of drought in Iran has been investigated. The one-month scale data from the standardized precipitation-evapotranspiration index (SPEI 01) in the period 1956 - 2015 have been used. Based on the common numerical thresholds, the characteristics of the frequency and severity of drought for each pixel have been calculated and they are the basis for the analysis of the drought situation. Then, the frequency of drought severity classes was calculated and its trend was investigated using the non-parametric Mann-Kendall test. The findings indicated the spatio-temporal variability of drought frequency and intensity patterns in Iran. The frequency of mild droughts has decreased from south to north and from east to west; while the frequency of more severe droughts has increased from north to south and from west to east. The frequency of mild droughts in the southeast, northwest, and northeast has increased by 5 to 40 percent. While the frequency of more severe droughts in most parts of Iran has increased between 10 and 20 percent. Variability in the frequency of more severe droughts is more pronounced in the Central Plateau catchment area as well as in the Persian Gulf-Oman Sea. The trend of drought intensity is decreasing (drought intensification) at the same time as the prevailing rainfall regime in Iran. A significant increase in drought intensity (wet season intensification) is observed only in southeastern Iran at the same time as the monsoon regime. However, extra-arid and arid regions of southeastern Iran are affected by the frequency and severity of drought and have a high degree of vulnerability.

In recent years, indirect methods such as remote sensing and data mining have been used to estimate soil salinity. In this research, the electrical conductivity of 94 soil samples from 0 to 100 cm was measured using the Hypercube technique in the Saveh plain. 23 types of input data were used in the form of topographic and spectral categories. Land area parameters such as the Topographic Wetness Index (TWI), Terrain Classification Index (TCI), Stream Power Index (STP), Digital Elevation Model (DEM), and Length of Slope (LS) were considered as topographic inputs using Arc-GIS and SAGA software. Also, salinity spatial and vegetation indices were extracted from Landsat 8 images and were considered spectral inputs. The GMDH neural network was used to model salinity with a ratio of 70% for training and 30% for validation. The results showed that the soil salinity values were between 0.1 and 18 with mean and standard deviation of 5 and 4.7 dS/m, respectively. Also, the results of modeling indicated that the statistical parameters R2, MBE, and NRMSE in the training step were 0.80, 0.06, and 42.1%, respectively. The same values in the validation step were 0.79, 0.13, and 48.7%, respectively. Therefore, the application of spectral, topographic, and GMDH neural network indices for modeling soil salinity is effective.

Groundwater is always considered one of the important water resources, especially in arid and semi-arid regions of the world, such as Iran. In recent decades, it has decreased drastically due to excessive use. The objective of this study was to determine the best interpolation method and evaluation of the spatiotemporal variations for the groundwater level in the Sahneh-Biston plain of Kermanshah province during three decades from 1991 to 2020. At first, four Gaussian, linear, spherical, and power semi-variograms were obtained for observations. Then, the best semi-variogram and interpolation methods were selected among the evaluated methods for zoning the groundwater level in the region. The lowest value of the sum of RMSE, MBE, and MAE error criteria and the highest coefficient of determination (R2) between observations and estimates in all three decades and the average of the entire period were calculated and considered to evaluate the most appropriate semi-variogram and interpolation methods for spatial distribution. The results showed that the ordinary kriging method with Gaussian semi-variogram is the best method to estimate the groundwater level in the Sahneh-Biston plain. The average difference between the minimum and maximum groundwater levels based on the observation wells of the study area and the zonation method is from 1279 to 1372 meters and 1289 to 1409 meters during the studied period time, respectively. The groundwater level is placed in more depth with the proximity to the central and southern regions. The maximum decrease and increase of groundwater level variations have been 12 and 19 meters during three decades, respectively. Also, the underground water level variations during these three decades showed that both the second and third decades compared to the first decade and the third decade compared to the second decade have increased in more than 50% of the region. This increase can be caused by the optimum management and water use in these years. Therefore, groundwater level monitoring provides effective help for experts and users in planning and optimal management of groundwater for the sustainable development of water resources.

In this study, the accuracy of the Adaptive Network-Based Fuzzy Inference System (ANFIS) in integrating with the Gray Wolf Algorithm (ANFIS-GWO) in predicting groundwater level was evaluated for the first time using unpublished observational data from 1998 to 2018 in the Zarandieh aquifer, central Iran. Three observational wells were randomly selected for analysis. Assessment of evaluation criteria demonstrated that among the proposed scenarios using the hybrid model, the D scenario was selected as the optimal scenario with input data including the previous month's groundwater level, precipitation, temperature, and groundwater extraction. In the D scenario, parameters including MAPE, RMSE, and NASH were 0.29 m, 0.47 m, and 0.99, respectively for the first observational well. Also, C scenario with input data including the previous month's groundwater level, precipitation, and groundwater extraction for the second observational well, for the same parameters mentioned above equal to 0.20 m, 0.26 m, and 0.99. As well for the third observational well, the A scenario with input data including the previous month's groundwater level for the same parameters equal to 0.29 m, 0.41 m, and 0.99 as the optimal scenarios were selected using the ANFIS-GWO model. Based on the results, the Gray wolf algorithm in training the ANFIS model was able to reduce the average forecast error by equal to 0.03 (RMSE) and 0.02 (MAPE) meter and increased the average NASH value equal to 0.01 and increased the accuracy of predictions.

Changes in land quality factors were investigated according to the change in land use of two conventional cropping systems in Khuzestan (Dimcheh region, periodic cultivation system, sugarcane, forest, and deforesting in Zaras region). The results showed that by the change of forest land use, organic carbon from 0.93 to 0.55%, cation exchange capacity (CEC) from 19.6 to 13.3 cmol/kg, C/N from 7.4 to 3.8%, the mean weight diameter of aggregate (MWD) from 1.7 to 1.3%, and microbial respiration from 0.11 to 0.06 mg of CO2 /gr of soil per day decreased and in contrast, the dispersible clay from 4.6 to 19.3% increased. PCA analysis for the parameters showed that five factors justified more than 90% of the variance in the values of FC, PWP, AW, and AF. In the Dimcheh region, the average volumetric moisture content of FC from 31.3% to 27.3%, available water from 12.9% to 9.8%, dispersible clay from 56.1% to 12.3%, and bulk density reduced from 1.6 to 1.4%, organic carbon from 0.45 to 0.78%, C/N from 6.3 to 10.0%, microbial respiration from 0.01 to 0.04 mg of CO2 /gr soil per day and MWD of aggregates increased from 0.77 to 1.3 mm. Five factors including FC, AW, BD, DC, and OM explained more than 90% of the variance.

Runoff is formed by spending some time after rain and significantly depends on rainfall intensity, soil moisture, and slope. One of the fundamental questions about runoff is the time that it starts to create. In this research, the runoff start time in sandy loam soil was evaluated experimentally under different conditions using a precipitation simulator machine. The rainfall intensity parameters of (60, 80, and 100 millimeters per hour) and the slope of (0 and 5 percent) were investigated. The rainfall was created in the three soil treatment types completely dry (Sdry), the dry soil that had been saturated 24 hours before the test (S24hrlag), and the dry soil that had been saturated 48 hours before the test (S48hrlag). Eighteen tests were conducted on this soil. At the end of each test, the soil moisture was measured. The experimental results were compared with the numerical model of Green-Ampt. According to the Kendall and Spearman correlation test results, as the rainfall be intense, the start time of the runoff is lower. Also, the runoff starts at a faster time in the slope of 5 percent for every three types of soil. Also, the results of starting time of the runoff in the soil with a delay of 48 hours in the rain compared to the soil with a delay of 24 hours in the rain are closer in all of the rainfall intensity and slopes compared to the case of dry soil. Therefore, in the experiments related to a delay of 24 hours, the time of the start of runoff decreases. While in tests with a delay of 48 hours, it was not much different from completely dry soil. Also, the Green-Ampt results are close to the experimental results (R2=0.9775), and the maximum difference between the two mentioned methods is 4.8 minutes. Therefore, it can be used with the Green-Ampt method to calculate the start time of runoff in sandy loam soil in different states of rainfall intensity and bed slope.