The aim of this research was to identify species suitable for plantation. We first identified species for potentially suitable for plantation based on ecological capabilities regarding soil properties. We determined the area of plantation for different species based on ecological capabilities. Then, we collected relevant data such as growth patterns of different species, labor requirements for plantation and plantation cost. A linear programming model and two integer linear programming models were used for optimization. The appropriate species based on ecological capabilities were ash, elm, maple, oak and bald cypress. A linear programming model was used based on ecological capabilities classification to determine the land area of different species for plantation. Then, two integer linear programming models were employed to select the species for plantation. We set ecological properties unequal for all of the species in the first run of the integer programming model. Two groups were classified: group one included maple and ash; group two included bald cypress,oak and elm. The second integer programming model assumed equal ecological properties for all the species.Results of linear programming showed that maple and bald cypress were appropriate for plantation at the site and their plantation areas should be 151.3 and 355.3 ha, respectively. Results of the first integer linear programming model showed that maple and bald cypress would be economically profitable for plantation. The results of the second integer linear programming model showed that only bald cypress would be appropriate for plantation. The aim of this research was to identify species suitable for plantation. We determined the area of ​​plantation for different species based on ecological capabilities. Then, we collected relevant data such as growth patterns of different species, labor requirements for plantation and plantation cost. A linear programming model and two integer linear programming models were used for optimization. The appropriate species based on ecological capabilities were ash, elm, maple, oak and bald cypress. A linear programming model was used based on ecological capabilities classification to determine the land area of ​​different species for plantation. Then, two integer linear programming models were employed to select the species for plantation. We set ecological properties unequal for all of the species in the first run of the integer programming model. Two groups were cl assified: group one included maple and ash; group two included bald cypress, oak and elm. The second integer programming model assumed equal ecological properties for all the species. Results of linear programming showed that maple and bald cypress were appropriate for plantation at the site and their plantation areas should be 151.3 and 355.3 ha, respectively. Results of the first integer linear programming model showed that maple and bald cypress would be economically profitable for plantation. The results of the second integer linear programming model showed that only bald cypress would be appropriate for plantation.