Background and Objectives: The groundwater level of Silakhor plain has been decreasing significantly with the occurrence of successive droughts, industrial growth and increasing water needs. In addition, the cropping pattern of the region in recent years has led to the cultivation of water crops, which sets the need for efficient management in the allocation of limited water resources in the region. In this study, to determine the optimal cropping pattern of major crops in Silakhor plain, with the aim of maximizing farmers' incomes and available water and land constraints, two approaches using Linear Programming and using Multi-Objective Meta Heuristic Algorithms in different exploitation scenarios have been investigated. Materials and Methods: Before using Linear Programming and optimization algorithms, in the first step, 100 different exploitation scenarios with equal intervals were determined for each crop year. Rainfall of the last 10 years in monthly and seasonal conditions was modeled using Artificial Neural Network and Genetic Programming and a better model according to the evaluation criteria of modeling. Then the rainfall of the next three crop years was forecasted and the resulting nutrition was estimated. Due to the need for proper exploitation from aquifers, it is necessary to have less exploitation than recharging in the coming years. In Silakhor plain, 50% of groundwater abstraction is used for horticultural, industrial and drinking products. Therefore, 45% of the feeding volume in each crop year was considered as the minimum exploitation and 140% of the exploitation in 2015 was considered as the maximum exploitation. One approach to solving constrained problems using metaheuristic algorithms is to constrain the problem using the penalty function and define the minimization of the penalty function as a goal. In this regard, in the second step, using Linear Programming, the optimal cropping pattern was followed, by the maximum income of farmers with limited water exploitationin each scenario and available land. Then, by defining the mentioned limitations as different penalty functions, the unresolved issue and maximizing the farmers' income function was considered as the first goal and minimizing the penalty functions as the second goal. The multi-objective optimization algorithm continues to operate until the response obtained from Linear Programming is reached with a maximum error of one percent,It is also an acceptable answer that the amount of the fine is zero. In other words, the answer in question must not exceed the defined limits. In this study, the performance of three types of static, dynamic and classified dynamics penalty functions in three multi-objective algorithms NSGA-II, SPEA-II and PESA-II. Are evaluated. The following equation shows the general form of the objective functions. Cost Function (I): Maximum Net Income Cost Function (II): Minimum Penalty Functions. Results: The results show that along with increasing groundwater exploitation, farmers' incomes also increase,However, in the exploitation of more than 223. 5, 222. 2 and 225. 1 million cubic meters for the cropping years 2020-2021, 2021-2022 and 2022-2023, respectively, the limitation of the total arable land in Silakhor plain prevents the increase of crop cultivation. As a result, the income of farmers in the region will not change. The results of the algorithms also show that the best performance among the algorithms in this issue belongs to the SPEA-II, PESA-II and NSGA-II algorithms with the number of iterations of 12. 1, 14. 5 and 17. 8, respectively. Among the penalty functions, on average in all three algorithms, the best performance belongs to the classified dynamics, dynamic and static penalty functions with the number of iterations of 13. 1, 13. 7 and 17. 5, respectively. Conclusion: Due to the decrease of groundwater level in Silakhor plain, determining different scenarios of groundwater abstraction and optimizing the cropping pattern appropriate to each scenario, in addition to increasing the economic productivity of the region, also facilitates the management of water resources. It is impossible to introduce a single algorithm and penalty function to solve all optimization problems. However, based on the results of this study, to solve the linear constraint problems, the use of the SPEA-II algorithm with a classified dynamics penalty function is recommended.