Climate change, land degradation, and desertification lead to the loss of carbon present in the soil and plants. The carbon dioxide (CO2) concentration in the atmosphere has reached 412 ppm. This is a rise of 47% since the start of the industrial period, when the concentration was close to 280 ppm. Therefore, the sequestration of carbon from the atmosphere to earth is the need of the hour. Many scientists have suggested agroforestry as a potent instrument for climate change mitigation as well as to fetch lucrative benefits.

The Ghodaghodi Lake Complex is a Ramsar site, Nepal’s first bird sanctuary, and has significant ecological and economic values. The lake complex is in the western part of the lowland of the Terai region. Numerous studies indicate a relation between the normalized difference vegetation index (NDVI), land use, and land cover with plant diversity. However, the association between terrestrial plant diversity and NDVI in the Ghodaghodi Lake Complex is unknown but has important implications due to potential land use changes.

This paper addresses a forest harvesting problem with adjacency constraints, including additional environmental constraints to protect wildlife habitats and minimize infrastructure deployment costs. To this end, we propose an integer programming model to include those considerations during the optimization of the harvest regime of a Mexican forest. The model considered was based on the Unit Restriction Model, a benchmark approach that merges the management units before the optimization process.