جستجوی مقالات مرتبط با کلیدواژه « شبیه سازی » در نشریات گروه « آب و خاک »

Water scarcity and the need for optimal water utilization in arid and semi-arid regions, including Iran, have encouraged water authorities and farmers to adapt modern irrigation systems likedrip irrigation, to make optimal use of water resources. The most important advantage of drip irrigation over other irrigation methods is its ability to control the amount of water applied to each plant. New irrigation methods focus on plant irrigation and not on land irrigation. In arid and semi-arid regions, a drip irrigation system is used to use water optimally and prevent wastage and evaporation. Factors such as soil texture, type of cultivated plant, amount of available water, distance of drippers and laterals, the wetted surface, and the dimensions of the moisture bulb under the soil surface are involved in the design of the drip irrigation system. Due to the variety of soil textures in the earth, the movement of water under the soil surface is different in all kinds of textures, therefore, knowing exactly how water moves in the soil and how the moisture bulb is distributed under the soil surface is of particular importance. The purpose of this study is to investigate the movement of moisture bulbs, check their dimensions under the soil surface in different soil textures and flow rates, and evaluate the capability of the Moment analysis method to simulate this process under various conditions.

To simulate the moisture bulb in different soil textures, detailed information on the physical properties of the soil, including the percentage of particles that make up the soil texture, bulk and real density, porosity, and saturated hydraulic conductivity, is required. In this research, the simulation of the moisture front in different soil texture was conducted using Rosetta software, which defines 12 types of soil textures. In these tests, the source of soil power was considered as surface and point. The total feeding volume of each type of soil texture is 24 L, and this volume was used with different flow rates of two, four, six, and eight L s-1. To numerically simulate the progress of the moisture front, Hydrus software was used. Then the analytical simulation of the moisture front was done using the equations of the Moment analysis method. In this study, an ellipse was drawn to represent the moisture bulb simulated by Hydrus software at different times for the applied flow rates. Coefficient k was used to draw the ellipse, and its appropriate value was determined by minimizing the difference between the model and Hydrus results.

To calculate the moments, the first step is to obtain the values of M00 According to the applied flow rates of two, four, six, and eight L s-1 and the amount of volume intended to feed all types of soil texture, i.e., 24 L, the duration of irrigation is 12, 6, 4, and 3 hr, respectively. The comparison of moisture distribution over all periods and soil textures showed acceptable results, and the distributed subsurface moisture values were similar. In the study of clay texture, with time from the start of irrigation, the difference in the total amount of distributed moisture increased, and the reason for this result is the decrease in the permeability of the clay due to the filling of fine pores. The results indicated that σx2 values changed with the increase in irrigation duration. The highest variance was found in sandy clay with a flow rate of 8 L s-1 (1503.3 cm2), while the lowest variance was observed in clay texture with a flow rate of 4 L s-1 (368.6 cm2). By increasing the amount of applied discharge, σz2 increases and the slope of this increase is different in each soil texture, according to the characteristics of that texture. Also, the effect of irrigation duration on the value of σz2 is evident. In other words, the longer the duration of irrigation, the more the amount of variance changes.

In this research, the accuracy of the Moment analysis method in predicting moisture distribution from drip irrigation was evaluated using results from Hydrus and Moment analysis. The Hydrus results demonstrated that the moisture bulb expanded over time in both the horizontal and vertical directions. The results also indicated higher flow rates increased the horizontal expansion of the moisture bulb, while the duration of irrigation affected both horizontal and depth expansions. Using the moment analysis method, the center of mass of water distribution in the soil and the changes in the moisture front along the x and z axes were determined. By examining and comparing the dimensions of the moisture front resulting from Hydrus and ovals, it was observed that there is a suitable compatibility between the two methods. Therefore, the Moment analysis method can be relied upon to estimate the dimensions of the moisture bulb in drip irrigation. It also provides an efficient and accurate approach to reducing the time and cost of field experiments.

Crop modeling is one of the widely used methods for simulating the crops yield under different farm management. However, paying attention to the sensitivity of these models to the input parameters can improve the calibration operation and reduce the simulation error. For this reason, in the current research, we analyzed the sensitivity of the AquaCrop model under planting date treatments (D1: March 5 and D2: March 15) and irrigation water salinity (S1: Karon water with an average salinity of 2 dS.m-1, S2: combination of Karon water and drains with an average salinity of 4 dS.m-1 and S3: drain with an average salinity of 6 dS.m-1) and cotton cultivars (C1: Khursheed, C2: Golestan, C3: Sajedi and C4: Khordad). For sensitivity analysis, the normalized water productivity (WP*), maximum crop transpiration coefficient (KCTrx), initial crop coefficient (CCo), crop growth coefficient (CGC), crop decline coefficient (CDC) and harvest index (HI) were used. Sensitivity analysis was done using Beven method. The results showed that the highest sensitivity was in parameters of HI (0.69) and WP* (0.6) and the lowest sensitivity was in parameters CDC (0.02) and CCo (0.3). The increase in planting date from D1 to D2 increases the sensitivity of the AquaCrop model to the HI (19%) and WP* (15%), and the increase in salinity from S1 to S3 increases the sensitivity of AquaCrop to the HI (42%) and WP* (50%). Therefore, in the farm management similar to the current research, paying attention to the two aforementioned parameters can cause the accuracy of the results in the calibration (and validation) stage.

Due to the lack of water resources (quantitatively and qualitatively) and application of different irrigation management in terms of interval and depth of irrigation water (irrigation scheduling) in agriculture and orchards, it is necessary to apply proper irrigation management to prevent the accumulation of salt in the root zone. Pistachio trees are resistant to drought and can produce an average yield even with low amounts of water. The irrigation interval of pistachio trees depends on factors such as soil texture, evapotranspiration, water and soil salinity, irrigation method, tree age and available water amount. Various irrigation methods are used for pistachio trees, the most important of which are surface irrigation, surface drip irrigation, subsurface drip irrigation, and bubbler irrigation. Given that there is no comprehensive information on irrigation management and quality of irrigation water for pistachio orchards in the drip irrigation method of main pistachio growing areas of Semnan (Damghan), identifying different pistachio irrigation managements provides this possibility by using different irrigation scenarios and valid and practical models, suitable irrigation management can be performed in order to increase yield and simulate soil salinity. Three orchards were selected from Damghan city of Semnan province (with a focus on pistachio orchards and surface drip irrigation systems) in order to investigate performance, irrigation management and soil salinity changes. In this regard, information including longitude and latitude, area of orchards, number and arrangement of drippers, irrigation scheduling and soil sampling performance for measuring salinity was done.

Using the most accurate methods and models to simulate the impact of climate change on meteorological variables in different regions of the world is of utmost importance. In this study, the accuracy of 10 AOGCM models related to the sixth assessment report of the IPCC (CMIP6) was investigated for simulating temperature and precipitation in the Sefidrood Basin. For this purpose, observational data of temperature and precipitation from 16 weather stations located in the basin during the time period from 1980 to 2014 were compared with the output of AOGCMs. The Kling-Gupta Efficiency (KGE) index was utilized for this comparison. The comparison was conducted on both annual and monthly time scales, and the more accurate models were identified for each time period. The results indicated that the accuracy of AOGCM models in estimating temperature in the study area was higher than their accuracy in estimating precipitation. Additionally, different models exhibited varying capabilities in simulating these variables across different months. Based on the results obtained, the MIROC6 and MRI-EMS2-0 models performed better than other models in estimating the temperature of different months. Furthermore, the HadGEM3-GC31-LL model showed a better performance than other models in estimating historical precipitation for most months of the year. Based on the results obtained, it is necessary to select and use the best AOGCM models for each month before conducting climate change simulation studies in the study area.

In many fields of water sciences, the moisture movement in the soil profile is very important. Due to the non-linearity of the governing equation in moisture movement in the soil profile, the comprehensive analytical solution of this equation is impossible, and there is a need to solve it numerically. The abundance of data and simulation steps reduces the flexibility in using existing software to solve the Richards equation and using this data as initial conditions. Therefore, having a numerical solution plan for the Richards equation will be very useful. Among the existing methods, the fully implicit finite difference method can continue the simulation unceasingly due to its unconditional stability in providing a non-iteration solution. However, using the finite difference method to solve the Richards equation leads to significant mass balance error. Some usual methods to reduce the mass balance error could not be applied to the non-reverse solution. By creating the Richards equation numerical solution in MATLAB software, this research introduces a method to reduce mass balance error, which is a simple method to apply to all kinds of numerical solutions.

In this research, a Non-Iteration implicit solution for the Richards equation was developed. After creating the Richards equation numerical solution, the difference between the entering flux to the soil column and the added water to the same soil column, is added to each node to reduce the mass balance error. This amount is added to the nodes that their matrix pressure head difference (PD) compared to the initial state, is more than a threshold value. This amount applies as a coefficient proportional to the humidity of each node. To investigate the effect of various factors on the performance of the created numerical solution, a series of 11 experiments were designed with different characteristics in terms of infiltration rate, infiltration duration, duration of moisture redistribution in the soil profile, and initial volumetric water content. The optimal value for PD in each experiment was determined in such a way as to minimize the mass balance error. For the integrity of the results, the average PD values equal to 10 cm were used to simulate all cases. These series of experiments were simulated by two-dimensional HYDRUS software, and the results were used to verification of the created numerical solution.

The results showed that the most mass balance error is forced in the infiltration phase. Using the proposed method to preserve the mass balance caused a significant reduction in the mass balance error. The absolute value of the mass balance error varied in a narrow range and was always less than 3.5 percent. It is expected that time and space steps have a considerable effect on the mass balance error. The absolute mass balance error decreased linearly from 3.36 to 3.33, with the infiltration duration increasing from 10 to 30 minutes at a constant infiltration rate of 1 mm 〖min〗^(-1). The absolute mass balance error increased linearly from 3.36 to 3.37% by the initial volumetric water content increasing from 13 to 17% and decreased from 3.40 to 3.35 percent by the increase in the duration of the moisture redistribution in the soil profile from 300 to 900 minutes. Also, the Mean Absolute Relative error (MAE) increased from 1.17 to 9.53% by the infiltration time increasing from 10 to 30 minutes; the MAE increased from 1.17 to 8.89% by the infiltration rate increasing from 1 to 3 mm 〖min〗^(-1); and the MAE increased from 1.17 to 6.89% by the volumetric water content increasing from 13 to 17%.

The results showed that this method leads off the simulation to a significant reduction in the mass balance error; so that, in all experiment cases, the mass balance error ever fluctuated in a small range, less than 3.5%. The results comparison of the created simulation and the HYDRUS software showed that this method provided acceptable results, in cases such as summer rains where the infiltration duration is much less than the water redistribution time in the soil profile. Also, the highest value of the MAE in comparing the created simulations with the obtained results from the HYDRUS software is less than 10%. It proves a good agreement between the created simulation and the HYDRUS software.

The amount of sediment production, the manner and time of sedimentation, the size and composition of sediment grains, and transport among the waterways network are important features of the sedimentation regime of Watersheds; Because changes in each of these factors cause changes in Watershed performance. Therefore, sediment production is a reflection of the importance and amount of erosion processes and sediment sources in the upstream parts of the Watershed and how sediment is transported and deposited from the moment of movement of erosion materials from the point of separation to the exit of it. Sediment connectivity indicators indicate the spatial changes of connectivity patterns in different parts of the Watershed and provide a suitable estimate of the contribution of sediment sources and sediment transport routes. For this purpose, investigating the spatial pattern of sediment flux at the Watershed scale is of particular importance in developing comprehensive management plans and measures to control erosion and sedimentation. Although various indicators and models have been developed in this field, their performance has not been evaluated based on observational data and statistical methods. This research aims to analyze the sediment flux pattern of the Neyriz Watershed located in the east of Fars Province, based on sediment connectivity and sediment transfer capacity indicators, and compare their performance based on field sediment evidence.

At first, a digital model of the ground elevation of the Neyriz Watershed with a spatial resolution of 12.5 m was prepared and its drainage network was extracted. Then, the sediment connectivity index was calculated by considering the upstream and downstream features of the Watershed and by considering the roughness factor as the sediment movement resistance factor, and the sediment connectivity map of the Watershed was made. Then, the sediment transport capacity index map was also prepared based on the concept of topography and using the digital model layer of land height. For this purpose, the sediment transport capacity index was calculated for each pixel of Neyriz Watershed in SAGA-GIS software, and a sediment transport map was prepared. Based on field visits to different parts of the Watershed, 30 positions that had evidence of sediment transfer were recorded as observation points of the sediment transfer event. Also, 30 other positions that did not have signs and evidence of sediment transfer were added to the validation database as observation points of no sediment transfer event. The corresponding values of each index were also extracted in the geographic information system and based on the available information, using evaluation methods based on the error matrix including true skill statistic (TSS), efficiency (E), and F score (F-score), the validation of the mentioned indicators was done quantitatively.

The value of the sediment connectivity index of the Neyriz Watershed in Fars Province varied from -7.24 to 2.43 and its median value was -4.36. The spatial pattern of the sediment connectivity index in this Watershed is such that the middle and western parts have a low amount and the northern, southern, and eastern parts have more amounts. In this research, the drainage network of the Watershed was introduced as the target of receiving the sediment; Parts of the slopes of the Watershed level, which had the conditions of sediment production and transfer were connected to the drainage network in terms of the transfer path, have shown a higher value of connectivity index. This index provides valuable information for land management by considering the upstream characteristics of each point as well as the characteristics of the transfer path to the sediment-receiving target. The value of the sediment transfer capacity index varied from 7.2 to 23.16 and the average value was 10.19. The value of this index is high in the marginal parts of the Watershed where there are sloping lands, and the middle parts of the Watershed have a small value. Based on the findings, the index of sediment connectivity with the true skill statistic (TSS) of 0.833, the efficiency value (E) of 0.916, and the F score of 0.915 is a better performance than the sediment transport capacity index (TSS= 0.633, E=0.816, F-score=0.825). In addition, based on the values of the false positive component in the error matrix, the sediment transport capacity index predicts high sediment flux potential in many situations; While in the field observations, it was not true.

Based on statistical evaluation criteria, the sediment connectivity index has been able to better describe the state of sediment flux and is more consistent with field realities. So that the sediment connectivity index, in addition to considering the characteristics of the upstream area of each point of the Watershed, is possible to consider the path of the sediment particle to the target location of the receiver (such as the nearest branch of the network drainage) has made it possible. Based on the findings obtained in this research, although the sediment transfer capacity index in some parts of the Neyriz Watershed of Fars Province was in line with the sediment connectivity index; However, due to the energy-based nature of the flow, this index only considers the local conditions of the points and ignores the features of the upstream area as well as the process of transporting sediment particles to the downstream. Therefore, it is suggested to use the sediment connectivity index in the erosion and sedimentation studies of Watersheds. Because when information about the amount of sediment production is not available; the Sediment connectivity index can provide useful information about the sediment transfer process in the Watershed.

Soil contamination due to heavy metals is a global environmental issue. One vital aspect for understanding the impact of a contaminant in porous media is to describe their transport behavior using appropriate models. The governing equations for solute transport in soil consist of the convection–dispersion equation (CDE) and the mobile–immobile model (MIM). Mathematical models are usually used to evaluate solute transport in porous media. The first model used to express the transport of solutes and pollutants in porous media is CDE it provides acceptable and satisfactory results in homogeneous soils in laboratory tests. Hydrus-1D is a modeling environment for simulating water, heat, and solute movement in one-dimensional variably saturated media. Sensitivity analyses and model identification are standard approaches in modeling applications to investigate the relative importance of model components that control the system’s behavior. The sensitivity analysis is applied to identify the parameters that influence the model performance most. The sensitivity analysis is defined as the rate of variation in the model outputs due to changes in the input parameters. This study is a fundamental practice for analyzing the behavior of a model under different conditions of an application. The sensitivity analysis could be a practical and powerful tool for investigating the role and importance of model components, such as parameters and forcing data on the model responses.

The loamy soil samples were collected in both disturbed and undisturbed forms from a farm in the Qaramalek area with appropriate humidity located in western Tabriz, Iran, at 38º 5' 59.89'' north and 45º 12' 38.57'' east. To determine and present breakthrough curves, concentration values are required throughout the laboratory columns at different times. To simulate the CDE model, Hydrus software was used. Solute transport parameters such as diffusion coefficient (D), distribution coefficient (Kd), and dispersion coefficient (β) were estimated using soil hydraulic parameters and data related to the metal concentration of cadmium, nickel, and zinc by an inverse modeling method. A sensitivity analysis was carried out for the identification of the most influential factors on the model output. This method examines the impact of input data on a given model and its actual conditions. In line with this purpose, in each run, one input data was changed to a value equal to Positive and negative five to 15%, and the other input data was kept constant. To identify the effect of the input parameters of a given model on its output, the sensitivity analysis for the Hydrus model was utilized. The parameters of hydrodynamic dispersion coefficient (D), distribution coefficient (Kd), and spreading parameter (β) were changed between five to 15 %. Sensitivity analysis was carried out on cadmium, nickel, and zinc metals with densities equal to 50, 100, and 150 mg.l-1 in two disturbed and undisturbed soils.

Examining the breakthrough curves of cadmium in disturbed and undisturbed soils shows that the fitted curves using the Hydrus model and the measured curve almost coincide with each other, which is more obvious in disturbed soils. It should be noted that the model fits better in the disturbed soil than in the undisturbed soil. This may be due to the disruption of the structure the increase in the contact surface of the particles in the disturbed soil and the presence of heterogeneity in the undisturbed soil column. The simulation results show the transport of heavy metals (Zn, Ni, Cd) and Hydrus output have the highest and the lowest sensitivity to dispersion coefficient β and diffusion coefficient (D), respectively. In general, the impact of input parameters can be reported as follows: spreading parameter (β) > distribution coefficient (Kd) > dispersion coefficient (D). Therefore, it can be observed that D has a negligible effect on the model results; and consequently, measurement errors can be ignored.

Sensitivity analysis is used to analyze model behavior under different conditions. This analysis is used to investigate the relative importance of model components that control the system’s behavior. In this research, the transfer of hydraulic parameters of heavy metals Cd, Ni, and Zn in disturbed and undisturbed loam soil columns with initial concentrations of 50, 100, and 150 mg.l-1 was performed under the simulation of the Hydrus-1D model. The comparison of the simulated BTCs of the Hydrus-1D model and the measured data indicates a high agreement between the simulation curves and the measured data. Solute transport parameters such as hydrodynamic dispersion coefficient (D), distribution coefficient (Kd), and spreading parameter (β) were estimated using soil hydraulic parameters and data related to Cd, Ni, and Zn metal concentration by inverse modeling method. Based on the results of sensitivity analysis, the spreading parameter (β) and hydrodynamic dispersion coefficient (D) had the highest and lowest sensitivity, respectively. In other words, due to the significant effect of β changes on the output values of the model, this parameter should be measured more accurately and on the other hand, the measurement errors of parameter D can be ignored. The degree of sensitivity of the parameters was independent of the initial concentration of the elements.

Investigating the effects of climate change on improving agricultural productivity is of great importance. It is predicted that the rising temperature trend will continue in the 21st century, leading to changes in the climatic conditions of  regions. Changes in precipitation distribution ، temperature and water resources are considered as harmful consequences of climate change، which could have detrimental effects on agricultural production. . In this study، the impact of climate change and cultivation on different dates was examined in terms of the yield of the 'Parsi' spring wheat variety in the Qazvin Plain. This study covered the period from 2021to 2100, comparing two information sources, LARS-WG and DKRZ, for producing annual climate change data and using the Aquacrop model to simulate the plant’s response to the mentioned changes. During the periods 2021-2040، 2041-2060، 2061-2080، and 2081-2100، the best planting date (February 4, February 20, March 5, March 20, and April 4) for increasing wheat yield was evaluated.  According to the model results in both scenarios 4.5 and 8.5، in all four periods (2021-2040، 2041-2060، 2061-2080 and 2081-2100), the spring wheat yield will increase compared to the baseline.. The highest yield in all of these periods and models is predicted for the period 2081-2100 under the LARS-WG climate model of in scenario 8.5,  assuming planting on February 4, with an estimated yield of 12.43 tons per hectare and a standard deviation of 0.31 tons per hectare. the lowest yield is expected for the period 2021-2040 under LARS-WG climate conditions in  scenario4.5, with planting on April 4, with an estimated yield of 7.87 tons per hectare and a standard deviation of 0.36 tons per hectare. In all four periods (2021-2040، 2041-2060، 2061-2080 and 2081-2100), February 4 is recommended as the most suitable planting date to increase wheat yield in the Qazvin plain.

To Estimate Parameters and Evaluation of Cropsyst Model to Predict the Growth and Yield of Two Wheat Cultivars in Northeast Khuzestan, an field experiment was conducted in the form of randomized complete block design in a research farm in Ramhormoz with three replications. Factors were included irrigation regime at two levels (full irrigation and stop irrigation after flowering), cultivation date (six dates including: October 30, November 15, November 30, December 15, December 30, and January 15) and two varieties of wheat. (Chamran 2 and Sirvan). Based on the results, the critical day length, the ceiling day length, the number of basic vernalization days, the minimum growth rate in the absence of vernalization, and the saturation point of vernalization for Chamran 2 variety are equal to 14, 8, zero, and 0.9, respectively. and 24, and for the number of Sirvan, it was 13.8, 7.8, zero, 0.93 and 25, respectively. The results obtained from the simulation of flowering and physiological ripening stages based on RMSE, d, CRM and EF statistics showed that the Cropsyst model simulated the mentioned stages with good accuracy for different figures. The simulation results of dry matter of aerial organs for different figures showed that this model was able to simulate 95 to 98 percent of this parameter in different treatments of planting date and irrigation. The error value of the model based on the RMSE statistic for the dry matter of aerial organs in Chamran 2 and Sirvan cultivars was 627 and 351 kg per hectare, respectively.

Evaporation, the process by which water molecules escape a surface after absorbing sufficient energy to overcome vapor pressure, is a major contributor to water scarcity, especially in arid and semi-arid regions where heat readily facilitates this escape. Accurately estimating evaporation losses is crucial for effective water resource management, crop water demand prediction, and irrigation scheduling. Machine learning (ML) has emerged as a powerful tool for tackling the complex and stochastic nature of environmental problems. ML models excel at identifying relationships between predictor variables and outcomes (predictands), often surpassing traditional methods. However, their performance can vary depending on input factors and climatic conditions. Recently, hybrid techniques that combine multiple models have gained traction in climate and hydrology studies. These techniques leverage the strengths of different approaches within a single algorithm, potentially capturing more complex patterns in data series. This research will explore the potential of various individual ML models and propose a novel hybrid approach for estimating pan evaporation in Sistan and Baluchistan Province.

This study investigates pan evaporation simulation and prediction in Sistan and Baluchistan Province, Iran. Synoptic station data (1980-2019) served as model inputs, while pan evaporation measurements from these stations provided the observed values. In this research, in the approach of individual performance of data mining models, eight data mining models were used to simulate and predict evaporation from the pan. In addition to the individual performance approach, the combined VEDL approach was used to provide a hybrid model (a combination of the mentioned eight individual models of deep learning). In this hybrid approach to regression issues, the estimators of all models are averaged to obtain an estimate for a set called vote regressors (VRs). There are two approaches to awarding votes: average voting (AV) and weighted voting (WV). In the case of AV, the weights are equivalent and equal1. A disadvantage of AV is that all of the models in the ensemble are accepted as equally effective; however, this situation is very unlikely, especially if different machine learning algorithms are used. WV specifies a weight coefficient for each ensemble member. The weight can be a floating-point number between Zero and one, in which case the sum is equal to one, or an integer starting at one denoting the number of votes given to the corresponding ensemble member. the weight of each model was selected based on the accuracy of the model's performance using the evaluation criteria obtained from the training implementation section of individual models. the model’s performance was assessed using statistical measures, including R2, RMSE, MAE, and Taylor diagram.

The results showed that all the models had very good results in both the training and testing stages. All models exhibited excellent performance during training and testing. The Artificial Neural Network (ANN) achieved the highest accuracy in both phases at the Zahedan station (R² = 0.89, RMSE = 45.95 in training; R² = 0.96, RMSE = 44.18 in validation). It emerged as the best model for monthly pan evaporation prediction at this station. Other models also performed well, with the Support Vector Machine (SVM) and Random Forest (RF) models achieving R² values of 0.89 and 0.88 in training, respectively. Notably, the BART model ranked second in validation (R² = 0.96). The Tree Model (TM) had the lowest accuracy (R² = 0.84 and 0.93 in training and validation, respectively). Across all stations, ANN, SVM, and RF consistently delivered the best results in both training and testing. In the test phase, the SVM model outperformed others in Khash, Iranshahr, and Chabahar stations (R² = 0.94, 0.96, and 0.94, respectively). At the Saravan station, the RF model achieved the highest R² (0.94) during testing. To develop a hybrid data mining model, the Voting Ensemble for Deep Learning (VEDL) technique was employed with weighted voting in the training stage. The combined model significantly improved upon the best individual model. RMSE decreased from 45.95 to 33.1, R² increased from 0.89 to 0.94, and MAE improved from 32.92 to 23.9. Evaluation using the Taylor diagram further confirmed the superior performance of the VEDL model compared to the individual ANN model.

The results showed that among all the models, ANN, SVM, and RF models had the best performance in the two stages of training and verification. In the validation stage, the SVM model with R2 values equal to 0.94, 0.96, and 0.94 performed best in the Khash, Iranshahr, and Chabahar stations. At the Saravan station, in the Sensji validity stage, the RF model with an R2 value of 0.94 had the best performance among the models. The excellent performance of the models in the two stages of training and validation is another finding of the research, These results are consistent with the results of researchers who have expressed the appropriate efficiency of machine learning models in estimating evaporation/evaporation and transpiration in different climatic regions of Iran. The results of the combined model showed that the combined model improved the results compared to the best individual model so that the RMSE values increased from 45.95 to 33.1, the R2 values increased from 0.89 to 0.94, and the MAE value improved from 32.92 to 23.9. The use of the VEDL approach to estimate evaporation from the pan was a new approach that has not been used in past studies. Therefore, according to the results of this research, the proposed deep sensing model is proposed to estimate the evaporation of arid and semi-arid areas for water resources management and agricultural planning.

Climate change is one of the most important issues in the world in the 21st century which affects various sectors of agriculture, forestry, water and financial markets, and has serious economic consequences (Reidsma et al., 2009). In recent years, the management of agricultural water consumption has always been considered as one of the important issues in water resources management. Koochaki and colleagues (Koochaki and Kamali, 2006) by evaluating the climatic indicators of Iran's agriculture showed that during the next 20 years, the average monthly temperature will increase in almost all regions of the country, and the increase in evaporation and transpiration is one of the most important consequences of this warming. Simulated climate parameters can be obtained through different general GCM atmospheric models. Due to the low spatial resolution of these models, its output should be downscaled using dynamic or statistical methods.

The LARS-WG model predicts meteorological variables for a period of time in the future by using a series of basic and fine-scale meteorological data, output of one of the GCM models. Research has shown that the LARS-WG model has the necessary accuracy for this task. Calculating the amount of evapotranspiration and yield of very complex plants are time-consuming and dependent on spending a lot of money and limited to the tests performed, the shortness of the test time and also the limitation in the number of scenarios that are checked by the test. Therefore, plant models are considered and evaluated by researchers. The AquaCrop model has demonstrated commendable accuracy in various regions of Iran and globally for forecasting plant growth, water consumption efficiency, and evapotranspiration requirements. These predictions hold significant potential for optimizing irrigation strategies across different agricultural settings. AquaCrop is one of the applied agricultural models that was obtained from the modification and revision of FAO publication No. 33 by prominent experts from all over the world. In this study, the values of green water footprint of winter wheat plant (Pishgam) were estimated in climatic conditions obtained from LARS-WG model and DKRZ database under scenarios 4.5 and 8.5 and at different planting dates (15 October, 1 November, 15 November, 30 November and 15 December), in the next 4 periods (2021-2040, 2041-2060, 2061-2080 and 2081-2100) and by Aquacrop model.

The results showed that if planting date is on October 15, in the climatic conditions obtained from the LARS-WG model and under scenarios 4.5 and 8.5, in all future periods, the footprint of green water will increase compared to its value in the base period, and if planting is the rest of the dates, in each of the next 4 periods, the average green water footprint will decrease compared to its value in the base period. The results obtained for the DKRZ database show that the green water footprint attained for the dates of cultivation and periods investigated in scenarios 4.5 and 8.5 does not have a particular trend. On the planting dates of October 15 and November 1 for the periods of 2061-2080 and 2081-2100, the green water footprint will decrease and on the other three dates (15 November, 30 November, and 1 November) for these periods, there will be an increasing trend. On 15 December, for the DKRZ database, in both scenarios defined for all periods, an increase in green water footprint compared to the base period is reported. However, in the period of 2081-2100 in scenario 8.5, a decrease compared to the base period will be observed. The highest amount of green water footprint in all these periods and models for the period 2041-2060 under the climatic conditions of the DKRZ database in scenario 4.5, if the planting date is 15 October, it is estimated that the amount of water consumed is equal to 4272 cubic meters per ton with a standard deviation of 5018 cubic meters per ton is predicted. The lowest footprint of green water for the period 2081-2100 under the climatic conditions obtained from the LARS-WG model in scenario 8.5, if the planting date is on 15 December, is reported to be 232 tons per hectare with a standard deviation of 52.3 tons per hectare.

Due to population growth and Iran's location in arid and semi-arid regions of the world, the need for water and food has increased and as a result, the pressure on water and soil resources will be more than before. On the other hand, the risk of drying up Lake Urmia, which causes environmental problems in the region, requires macro-water planning for the region and the use of optimal cultivation pattern to deal with water scarcity. Therefore, optimal use of water preserves water resources and increases the quality of products. Today more than ever, increasing the production of strategic crops such as wheat requires the proper use of water resources. The main source of food for the Iranian people is wheat and related products, and any action that increases the yield of wheat due to limited soil resources, especially water resources, is important and necessary at the same time. In recent years, significant advances have been made in modeling product growth and development using mechanical models. Plant growth models are increasingly used in the analysis of agricultural systems and simulate the plant's response to growth factors using mathematical equations. The AquaCrop model is one of the dynamic and user-friendly models developed by the FAO. The AquaCrop model receives information about farm, plant, soil, irrigation and climate, and ultimately predicts important parameters such as crop. Wheat yield simulation allows efficient management and better planning under various environmental inputs such as soil and water. To achieve higher accuracy and less model error, field parameters must be properly calibrated by the model to achieve proper performance. Also, calibration of the model, if not done correctly, causes a high error prediction by the model, which leads to incorrect management, water loss, plant drought and other cases. Therefore, using a model that has accurate and close prediction to the AquaCrop model and requires fewer input parameters is essential, which saves time, reduces costs and eliminates calibration errors. However, this model requires relatively large input parameters and is a time-consuming model in the presence of multiple scenarios.3

In recent years, smart models have been able to show high accuracy and become reliable models. Therefore, in the present study, to solve this problem and develop a model with less input data, using the ANN, SVR and SVR-FFA intelligent models and creating 440 scenarios in 2 farms, the performance of the AquaCrop model was compared.99WestW2 farm is located in Miandoab city and has a yield of 6.588 (ton ha-1) and WestW10 farm is located in Mahabad city and has a yield of 5.05 (ton ha-1).

The results of the model are performed using 5 evaluation criteria of Correlation coefficient, Root mean square error, Nash-Sutcliffe coefficient, Wilmot’s index of agreement and, Mean absolute percentage error. The results of this study showed that for both 99WestW2 and WestW10 farms, the SVR-FFA3 model could have the lowest error rate, so that for the yield index, the RMSE value for the mentioned farms was 0.033 and 0.069 (ton ha-1), respectively. The use of three models SVR, SVR-FFA, and ANN and their comparison with the AquaCrop model to predict wheat yield has been done for the first time in this study. The SVR model was able to show the highest accuracy after the SVR-FFA model. For 99WestW2 farm, it can reduce the error rate to 0.043 (ton ha-1) and for WestW10 farm to 0.077 (ton ha-1) and show good performance. The ANN model, after the SVR model, was able to show acceptable accuracy. The ANN model for 99WestW2 farm was able to reduce the error rate to 0.123 (ton ha-1) and for WestW10 farm to 0.094 (ton ha-1). Finally, the ANN model had a relatively higher error than the SVR-FFA and SVR models, respectively, and showed a relatively lower performance than the two models.

Finally, the intelligent SVR-FFA, SVR and ANN models, despite having the least number of inputs, were able to predict yield values in the shortest time and with the highest accuracy. However, the results showed that the lower the model inputs, the weaker the model prediction. For further studies, it is suggested that the ANN model be combined using the firefly algorithm (MLP-FFA) to increase the accuracy of the ANN model and make more accurate predictions of wheat yield.

Most of the regions of the globe have faced changes in the precipitation regime and CO2 concentration with the increase in temperature (especially the minimum temperature). Considering that two-thirds of Iran is considered to be arid and semi-arid regions, which will cause climate changes and successive droughts, and this process will affect agricultural ecosystems and their performance. The aim of this research is to simulate the effects of climate change on growth, biomass yield and corn grain under climate change conditions in Shushtar and Safi Abad cities in the northern part of Khuzestan province in 2020-2050. For this purpose, the production data of precipitation, minimum temperature, maximum temperature and sunshine hours of the LARS-WG statistical model using HADCM3 atmospheric general circulation model and under emission scenarios A2 and B1 were used in the periods of 2020-2050. AquaCrop model was used to simulate plant performance and growth. Before use, the AquaCrop model was calibrated and verified by field data collected in the region. Then, the values of seed yield and biomass in future periods were simulated for the treatment of 50% of the water requirement of the plant. According to the results, under A2 emission scenarios, the evaluation results of the LARS-WG model showed that it had a suitable forecast for the climate parameters and the simulation of the future growing season. And the increase in minimum and maximum temperature will be the average increase in precipitation. The length of the growth period for each station will decrease with the increase of GDD according to the climatic parameters of the region under the influence of climate change. Biomass and grain yield will also increase by one to two tons in different horizons and scenarios assuming the current planting date and full irrigation remain constant.

An artificial neural network (ANN) is a powerful data-driven tool capable of learning the linear and nonlinear relationships governing different systems. However, determining the best-performing algorithm in terms of convergence speed and accuracy for a particular problem is still a fundamental challenge for users of artificial neural networks.

We investigated the most effective algorithm among widely used processes to simulate and estimate nonlinear water quality parameters. For this purpose, we constructed 42 models combining artificial neural network topology (single or multilayer) and training processes. The quality parameters’ data acquired at 107 wells throughout the aquifer of Qorveh-Dehgolan plain were used for training (data from 1996 to 2013) and to test (data from 2014 to 2016) each model.

The results showed that artificial neural networks with a hidden layer that benefits from the optimal number of neurons could simulate the aquifer behavior with high accuracy and in less time. Also, increasing the number of hidden layers while increasing the response accuracy increases the number of optimal network neurons and the duration of the problem analysis. Finally, artificial neural networks based on the Broyden-Fletcher-Goldfarb (BFG) method had the highest efficiency in simulating aquifer behavior, although the performance of the Levenberg Marquart (LM) method is very close to BFG. BFG is more efficient than LM due to its lower Mean Square Error and Standard Deviation (3.46 and 3.09, respectively).

The aim of the study is to evaluate the changes of mean temperature and precipitation parameters using IPCC Sixth Report models (CMIP6) under SSP scenarios (SSP1-2.6, SSP2-4.5 and SPP5-8.5) during the period of 2022-2100 in Kashan Plain. The mean temperature and precipitation data was obtained from 7 stations (Kashan, Kavir-e-Hosseinabad, Kamu, Ardestan, Alavi, Noushabad and Sensen) in Kashan plain considering the base period of 30 years (1984-2014). Also, 7 models were selected from the models of the sixth report (CMIP6). The post-processing of the output of the models was carried out using the linear ratio method. Nash-Sutcliffe indices (NSE), root mean square error (RMSE) and correlation coefficient (r) were used to determine the accuracy of the models. The annual trend of changes was investigated using Mann-Kendall test. Finally, the mean of IPSL-CM6A-LR and BCC-CSM2-MR models outputs was used to simulate mean temperature and precipitation changes in the future period. According to the results, in all of the studied stations, precipitation in the coming period will have a decreasing trend compared to the base period. The mean temperature will also increase in the future period compared to the base period. In Kavir-e-Hossein Abad, Ardestan, Noush Abad and Sen Sen stations, the intensity of temperature increase will be higher than Kashan, Kamu and Alavi. According to the predicted conditions, it is necessary to pay attention to comprehensive policies in the field of adapting to climate change in Kashan Plain.

Due to the upcoming climate threats, challenge of water shortage and its impact on the food security of the growing population in Iran, the optimal use of soil and water resources is very important. With the development of remote sensing technologies, free access to a variety of field data has become widely available, which can be used to reduce the uncertainty of simulation models. The aim of this study was to simulate actual evapotranspiration (ETa) and rice crop coefficients (Kc) during its growth stages using the SWAP model updated with satellite data and evaluate the accuracy of the results with/without updating. This research was conducted at the National Rice Research Institute of Iran in Rasht in the year of 2017. Based on the obtained results, total ETa measured by lysimeter and simulated by SWAP model with and without updating were 395.4, 373.2 and 363.6 mm, respectively. The average crop coefficients during the growth stages of vegetative, reproductive and ripening were estimated as 1.13, 1.49, 1.21, respectively. The crop coefficients for the proposed stages estimated by SWAP model without using satellite data were 1.02, 1.39, 1.04, respectively. After updating with satellite data, the crop coefficients were modified as 1.05, 1.43 and 1.07, respectively. Finally, the statistical analysis indicated that the SWAP model has a reasonable performance in estimation of ETa (RMSE=0.89; EF=0.98; R2=0.74) and rice crop coefficients (RMSE=0.53; EF=0.96; R2=0.63). The results indicate that the SWAP model combined with satellite data improved the accuracy of ETa estimation (RMSE=0.75; EF=0.99; R2=0.86) and rice crop coefficient (RMSE=0.40; EF=0.99; R2=0.74) at field scale.

The use of hydrological models in watersheds has always been of interest to water resources researchers. Hydrological simulation models are valuable tools for investigating challenging issues related to watershed management, such as the effect of climate change on water resources and the effect of urbanization on floods and droughts. Spatial distribution hydrological model WetSpa is used to simulate river flow at basin scale. The model uses the observed topography, land use, soil map, and daily meteorological time series (rainfall, evaporation and temperature) to predict hydrographs and distributional-spatial hydrological parameters of the basin. In this article, the object-oriented, modular and process-oriented model of WetSpa, which is prepared based on the flexible modeling approach, is applied to simulate the daily hydrograph in Khorramabad basin.

The inputs of the model include digital elevation maps, soil type, land use, and time series of precipitation, temperature, and potential evaporation and transpiration, which are from the statistics of 6 meteorological stations in a ten-year period (water year 84-85 until 93-94) is used. After preparing the inputs of the model, at first the maps of the distribution parameters are automatically generated in the map format by the GIS pre-processing component of the model. After that, the model is calibrated using a 5-year statistical period (water year 84-85 to 89-88) of precipitation, temperature, and potential evaporation and transpiration data. The model uses Thiessen polygons to apply precipitation, temperature, and evaporation data. For this purpose, the daily discharges of Jam Anjir hydrometric station located at the outlet of the studied watershed are used. Model calibration is done manually by determining the values of 11 global (general) parameters of the model, so that the best match between simulation and observational hydrograph is obtained. And finally, the validation of the model is carried out based on a 5-year statistical period (water year 89-88 to 94-93) and the values of the global parameters obtained in the calibration stage.

The maps of distributed parameters are produced, which after preparing the inputs of Mashdand's production model showed that the average potential runoff coefficient of the area is 63% and the concentration time of the area is 17 hours. In the following, according to the 11 global parameters, which symbol and range of changes are specified in table (3), the model global (general) parameters values  are obtained in the calibration stage. Comparing the simulated hydrograph by the model and the observed hydrograph in the calibration stage shows that the best match between the observed and simulated data is established with a correlation coefficient of 0.39. Validation of the model is also based on a 5-year statistical period (water year 89-88 to 94-93) and the values of global parameters. The output files of the model illustrate that 26.15% of the precipitation becomes runoff during the calibration period. During the validation period, the share of total runoff from precipitation is 26.42%. Moreover, the simulation results of the model demonstrate the ratio of evaporation to precipitation in the calibration and validation periods is 57.18 and 69.20%, respectively. Additionally, the results of the evaluation of the model based on the Kling-Gupta index (KGE) present the value of 0.68 for the calibration period and 0.74 for the validation period.

In this article, the effectiveness of WetSpa model is investigated in order to simulate the daily flow of Khorram Abad River at Cham Anjir hydrometric station. According to the results obtained from this research, it can be said that the Wetspa spatial distribution model has the ability to simulate the hydrological behavior of the basin with acceptable accuracy. The graphical comparison of the calculated and observed hydrographs for the calibration and evaluation period also shows a relatively good match between the two hydrographs. Examining the results of calculating of the water balance components by the model demonstrates that the outflow in the calibration and validation period accounted for 26.15 and 26.42% of the total precipitation respectively, seems logical considering the major land use of mountains and pastures in the irrigation basin.

Managing and monitoring the quantity of groundwater resources, along with assessing the intensity and direction of their flow, constitute fundamental principles for sustainable utilization of such resources. This study aims to utilize GMS software to simulate and analyze groundwater resources in Nobandegan plain, Fars, Iran. The objectives include deriving long-term variations in these resources, evaluating the water balance, and ascertaining the direction and intensity of groundwater flow.

Data related to pumping wells, observation wells, topographic features of the plain, and rainfall, among other pertinent information, were gathered using reports from Fars Regional Water Company. Subsequently, a comprehensive hydraulic model was developed to simulate, calibrate, and validate the groundwater flow within the plain. Following this, the overall trajectory and flow intensity within the plain aquifer were simulated.

The hydrograph pertaining to groundwater within the plain reveals a ~50 meter water table drop over the course of the investigation period.  The predominant groundwater flow within this plain is directed from the northern mountainous regions toward the southwestern expanses. The maximum hydraulic head is recorded at approximately 1283 meters, while the minimum stands at 1175 meters above mean sea level, indicating a 108-m drop in the aquifer water table elevation. The simulation results of the groundwater flow within the plain reveal an increased velocity in the northern regions and, to some extent, in the southern reaches of the plain. The findings of this study may be employed for groundwater quality studies on this plain.

The knowledge of crops response to salinity stress application methods in growth stages, can lead to better management of stresses. This research was done on the S.C 704 maize, in the mini-lysimeter space and in Qazvin region. The experiment was performed factorial and in a completely randomized design. The salinity treatments of soil saturated extract (the main factor) were applied at four levels of 1.7(S1), 3(S2), 5(S3) and 7(S4) dS.m-1. The crop growth stage treatments (sub-factor) were defined as one-stage of 6 leaves (C1), flowering (C2), seeds milking (C3) and two-stages of C1C2, C1C3 and C2C3. The research target was to simulate the maize yield by AquaCrop model, in conditions of pulsed salinity stress application in crop growth stages. The stress application data in one and two growth stages, were used for calibration and evaluation the AquaCrop model, respectively. From S1 to S4 treatment, the crop dry matter was decreased from 157.2 g. plant-1 to 115.9, 53.2, 77.7, 86.1, 97 and 46.5 g. plant-1 in the C1, C2, C3, C1C2, C1C3 and C2C3 treatments, respectively. The sudden application of salinity stress in a sensitive growth stage (C2 and C3 treatments), was caused the more damage relative to C1C2, C1C3 treatments. Because once application of salinity stress in the 6-leaf stage (C1) has caused the crop adaptation to future stresses (increasing the tolerance threshold). In AquaCrop model evaluation, the statistical parameters of CRM, EF, R2, RMSE, NRMSE and ME were equal to -0.084, 0.833, 0.91, 12.05, 11.34% and 18.32, respectively. Those showed good accuracy of AquaCrop model in simulation the maize yield. As a result, by management the salinity stress in crop growth stages, will be reduced the negative effects on the maize yield. With simulation the crop yield by AquaCrop model, different stress application states can be evaluated.

Nowadays, numerical modeling is known as a powerful method and tool in investigating practical phenomena. So that this method is used in many fields of engineering. Although the numerical model of an engineering problem is usually prepared based on several simplifying assumptions, it is possible that this numerical model does not reflect all aspects of the real problem and is not able to show the real behavior in reality. Direct updating methods use the analytical solution of the problem for this purpose, although these methods specify the necessary corrections without the need to repeat calculations and in one step, but in most cases such corrections do not have a physical meaning in reality but iteration-based updating methods require sensitivity analysis of the effective parameters in the problem in order to find the impact of each of them. One of the solutions to deal with the time-consuming problem of multiple re-analyses in numerical models prepared by commercial software during model updating based on sensitivity analysis is to replace the numerical model with an approximate representative model known as meta-models. The response surface method is one of the common methods for building such meta-models. The response surface technique is actually a test design method to select the design parameters in the experiment with the aim of optimizing some system response functions. From this point of view, the research regarding the response of the network of open channels to the arrival of the flood wave and predicting the characteristics of the flow in the branches of the network can be the solution to minimize the damages caused by it.

In this research, in order to update the numerical model of unsteady flow in the network of open channels by using the combined method of response surface and meta-exploration methods, a branch of the Garmabadr River (Ziquon) has been selected, and the desired data has been prepared and sorted. have been Hec-Ras software has been used as a numerical method. Neural network is also chosen for simulating the numerical method so that the desired parameters can be optimized in the process of using the genetic algorithm method. This area is located in the northern part of Iran and in the northeast of Tehran province. The geographical coordinates include longitude 51 degrees 32 minutes to 51 degrees 38 minutes and latitude 35 degrees 51 minutes to 35 degrees 58 minutes. MATLAB software has been used to build a meta-model or an alternative model. MATLAB is a software that can be called the language of mathematics and modern engineering sciences.

As it is clear from the results of the Hec-Ras program, the discharge of stations number 13 and 22 is equal to 308.03 and 27.89 cubic meters per second, respectively. The estimation error is equal to the difference between the estimated flow rate and the observed flow rate, which is 18.03 and 2.11 cubic meters per second, respectively. The percentage of the estimation error is equal to the ratio of the estimation error to the observed value is equal to 6.2% and 7%, respectively, which is an acceptable value for this research, and this result shows the appropriate performance of the genetic algorithm for optimizing the flow simulated by It shows the neural network, on the other hand, the low percentage of error also determines the accuracy of the intended program. As a result, the error of the genetic program is at an optimal level. As it is known, Manning's coefficients are very important and influential parameters in the behavior of the flow, the change in which causes significant changes in the simulation results. The effect of this uncertainty on the performance of numerical methods was clearly identified in this research. As a result, examining inputs that have inherent uncertainty and determining their appropriate values is a suitable solution for improving numerical methods.

 Genetic algorithm is a powerful method to determine the values of effective parameters in simulated flow performance, and its combination with neural network is a powerful tool for engineering purposes, including the analysis of the behavior of open channels, whose fast convergence is a strong point to choose in these problems. in similar cases.
 Manning's coefficients are very important and influential parameters in the simulated flow behavior, with a small change in it, the simulated model undergoes significant dispersion, which should be considered in the analysis process so that the results are not significantly different from reality.
 Optimizing an effective parameter in simulated flow behavior is an erosive and time-consuming process, and the use of neural network becomes much faster and more reliable due to the high convergence speed of this process.