Background: Land leveling is one of the most important steps in soil preparation and cultivation. Although land leveling with machines require considerable amount of energy, it delivers a suitable surface slope with minimal deterioration of the soil and damage to plants and other organisms in the soil. Notwithstanding, researchers during recent years have tried to reduce fossil fuel consumption and its deleterious side effects during this operation. The aim of this work was to determine the best linear model using Artificial Neural Network (ANN), Imperialist Competitive Algorithm-ANN, regression, and Adaptive Neural Fuzzy Inference System (ANFIS) to predict the environmental indicators for land leveling and to determine a model to estimate the dependence degree of parameters on each other. Methods: New techniques such as ANN, ICA, GWO-ANN, PSO-ANN, sensitivity analysis, regression, and ANFIS that using them for optimizing energy consumption will lead to a noticeable improvement in the environment. In this research, effects of various soil properties such as embankment volume, soil compressibility factor, specific gravity, moisture content, slope, sand percent, and soil swelling index in energy consumption were investigated. The study was consisted of 350 samples which were collected from 175 regions in two depths. The grid size was set 20 m × 20 m from a 70-ha farmland in Karaj province of Iran. Results: The models that reveals the relationship between the land parameters and the energy indicators were extracted. As it was expected three parameters; density, soil compressibility factor and, embankment volume index had significant effect on fuel consumption. In comparison with ANN, all ICA-ANN models had higher accuracy in prediction according to their higher R 2 value and lower RMSE value. Statistical factors of RMSE and R 2 illustrate the superiority of ICA-ANN over other methods by values about 0.02 and 0.99, respectively. Results also revealed the superiority of integrated techniques over other methods for prediction of complicated problems such as land leveling energy estimation. Conclusion: Results were extracted and statistical analysis was performed, and RMSE as well as coefficient of determination, R 2 , of the models were determined as a criterion to compare selected models. According to the results,