The current rate of transformation of ecosystems and its consequences on the provision of benefits or Ecosystem Services (ES) make necessary the need for information on the quantification of their provision and identifying areas with the most capacity to provide ES. The main goal of this study is to develop a method to locally measure four ES (forage, timber, firewood and carbon storage) provided by plant diversity and to explore which biophysical variables can predict the services at larger scales in order to map ES provision for a whole watershed. The study was undertaken in the Cuitzmala watershed (1082 km2), located in the pacific coast of
The most explicative model for forage (r²=70%) was explained mainly for the canopy closure index, slope, elevation and distance to stream variables. The best model for timber (r²=50%) was explained for the canopy closure, normalized vegetation index for the wet season, normalized infrared index for the dry season, distance to road, to towns and to streams. The model for the firewood (r²=62%) was explained for canopy closure, normalized vegetation index for both seasons, normalized infrared index for the dry season and slope. The best model for the carbon storage (r²=65%), was explained for the canopy closure, normalized vegetation indexes for both seasons, slope and distance to stream. These models were used to extrapolate the ES values to the watershed and finally map the services. The resulting forage map showed a higher provision in the middle part of the watershed, while the timber, firewood and carbon storage maps showed higher values in the lower and upper parts of the watershed. These maps can be used jointly with biodiversity maps to establish conservation priorities in the watershed.