Plan retention tree groups based on ecosystem service values enriched by hotspot analysis (Getis & Ord 1992). The delineation of the retention tree groups is applied from the highest hotspot values of ecosystem services.

### Input parameters

- Planning area Planning or cutting area as vector polygon format

- ID field Unique field of planning area if there is multiple planning areas

- CHM Canopy height model (m)

- DTW Dept-to-Water -index (m)

- Vegetation zones Vegetation zones of boreal forest as vector polygon format

- Water body Water bodies to consider in Phosphorus retention calculation

- Forest grid Attibute data of forest grid need to have "fertilityclass","MEANHEIGHTPINE","MEANHEIGHTSPRUCE","MEANHEIGHTDECIDUOUS"

- Retention tree count (trees/ha) Count of retention trees. Minimum of the PEFC sertificate is 10 trees / ha

- Advanced parameters Weighting of the ecoystem services. Pret emphasize water protection, DTW biodiversity and water protection, Ds biodiversity, and DWP biodiversity and climate

### Outputs

- Trees Single-trees of planning area with ecosystem services values

- Retention tree groups Proposed retention tree groups based on weight of the indices

Algorithm author: Mikko Kesälä