National-scale mapping of tree species composition in Switzerland from high-resolution aerial imagery and multimodal data
Organizations: Photogrammetry and Remote Sensing Group, ETH Zurich, Zurich, 8093, Switzerland · Forest Resources Management Group, ETH Zurich, Zurich, 8093, Switzerland · Ecosystem Ecology Group, Swiss Federal Research Institute WSL, Zürcherstrasse 111, Birmensdorf, 8902, Switzerland · ETH Swiss Geo Lab, D4 Business Village, Platz 1a, Root, 6039, Switzerland ·
Abstract
Tree-species composition shapes forest productivity and how forests respond to climate disturbances such as drought, windthrow and insect outbreaks. Yet most national forest maps assign each pixel a single dominant species and cannot show the share of accompanying species. We mapped the basal-area shares of ten main tree species and a residual Other class across Switzerland at 30 m resolution. A deep-learning model combining 0.1 m aerial RGB imagery, airborne laser scanning (ALS), annual AlphaEarth satellite embeddings, and climate variables was trained on basal-area observations from 3,992 Swiss National Forest Inventory (NFI) sites. In spatial cross-validation, the model reached a macro- of 0.485 across the 11 classes. On an independent forest-reserve dataset, the map reached a macro- of 0.552 for individual 30 m pixels and 0.733 for plots, and identified the dominant class in 73.2% of pixels and 81.1% of plots. The map shows maple and ash as widespread minor components that dominant-class maps omit. It also reproduces the elevational sequence from broadleaved lowland stands through spruce to larch and Arolla pine near the treeline. By resolving the species mixture within each 30 m pixel, the map complements field inventory statistics with spatially continuous information that can support forest monitoring, the management of mixed stands, and climate-adaptation planning across Switzerland.
Figures & tables
| Study/product | Prediction target | Predictor modalities | Scale |
| Categorical | |||
| Francini et al. (2024) | 7 species groups | S-2 | Netherlands |
| Mu et al. (2025a) | 8 species/species groups | S-2 | Germany |
| Waser et al. (2017) | 2 leaf types; broadleaf probability | Aerial multispectral + terrain | Switzerland |
| Blickensdörfer et al. (2024) | 11 species/species groups | S-1 + S-2 + environment | Germany |
| Koch et al. (2024) | 15 species | S-2 + terrain | Switzerland |
| Modality | Dataset | Spatial resolution | Acquisition time |
| Aerial RGB | SWISSIMAGE | 0.1 m | Leaf-on acquisition nearest the NFI survey year |
| ALS | swissSURFACE3D | 15–20 points m -2 | Single available acquisition |
| Satellite embeddings | AlphaEarth Foundations | 10 m | NFI survey year |
| Climate | CHELSA BIOCLIM+ | 1 km | 1981–2010 long-term mean |
| Model | Backbone | Leaf level | Species level | |||
| MAE | MAE | |||||
| ResNet-50 | ResNet-50, ImageNet | 0.109 0.004 | 0.757 0.020 | 0.058 0.002 | 0.380 0.043 | 0.515 0.027 |
| Random Forest | – | 0.167 0.004 | 0.669 0.017 | 0.083 0.001 | 0.313 0.012 | 0.342 0.031 |
| AlphaEarth MLP | – | 0.149 0.008 | 0.667 0.027 | 0.084 0.002 | 0.157 0.014 | 0.194 0.012 |
| Ours | ViT-L/16, SAT-493M | 0.103 0.006 | 0.781 0.020 | 0.052 0.002 | 0.446 0.062 | 0.570 0.043 |
| ViT-L/16, LVD-1689M | 0.097 0.004 | 0.794 0.013 | 0.051 0.002 | 0.444 0.063 | 0.566 0.051 | |
| Modality | Leaf level | Species level | ||||||
| RGB | ALS | AEF | CHELSA | MAE | MAE | |||
| ✓ | ✓ | ✓ | ✓ | 0.093 0.003 | 0.803 0.017 | 0.048 0.001 | 0.485 0.038 | 0.585 0.014 |
| ✓ | ✓ | ✓ | 0.094 0.002 | 0.808 0.013 | 0.049 0.002 | 0.465 0.056 | 0.581 0.031 | |
| ✓ | ✓ | ✓ | 0.093 0.002 | 0.814 0.013 | 0.048 0.001 | 0.483 0.052 | 0.593 0.026 | |
| ✓ | ✓ | ✓ | 0.099 0.004 | 0.791 0.017 | 0.050 0.002 | 0.464 0.070 | 0.580 0.029 | |
| ✓ | ✓ | 0.099 0.003 | 0.793 0.018 | 0.051 0.001 | 0.461 0.044 | 0.581 0.027 | ||
| Metric | Spatial 5-fold CV | Independent, plot-level | Independent, pixel-level |
| NFI ( ) | Forest reserve ( ) | Forest reserve ( ) | |
| 0.483 | 0.733 | 0.552 | |
| MAE | 0.048 | 0.038 | 0.042 |
| RMSE | 0.125 | 0.088 | 0.109 |
| Dominant-class accuracy | 0.734 | 0.811 | 0.732 |
| Dominant-class | 0.585 | 0.704 | 0.596 |
Appendix figures & tables18 assets
Supplementary material from the paper’s appendix.
Appendix
| Step | Screening criterion | Description | Tree points | Sites |
| 1 | Raw NFI tree observations | All available tree records | 63,969 | 5,998 |
| 2 | Recent inventory observations | Recorded in 2017 or later | 42,831 | 4,052 |
| 3 | Canopy-layer screening | Exclude the lowest stratum | 39,863 | 4,043 |
| 4 | Quality control and completeness | Exclude disturbed sites and sites with degraded imagery; require complete multimodal records | 39,380 | 3,992 |
| Pixel side length | 10 m | 14 m | 18 m | 20 m | 24 m | 30 m (used here) |
| Sites fully contained (%) | 10.9 | 24.0 | 49.2 | 65.5 | 90.6 | 99.0 |
| Tree points contained (%) | 34.6 | 61.2 | 82.7 | 90.2 | 98.1 | 99.9 |
| No. | English | Latin | Included species |
| 1 | Spruce | Picea spp. | Picea abies ; Picea spec. |
| 2 | Fir | Abies alba | Abies alba |
| 3 | Pine | Pinus spp. | Pinus sylvestris ; Pinus nigra ; Pinus strobus ; Pinus mugo arborea |
| 4 | Larch | Larix spp. | Larix decidua ; Larix kaempferi |
| 5 | Arolla pine | Pinus cembra | Pinus cembra |
| 6 | Beech | Fagus sylvatica | Fagus sylvatica |
| Taxon | Leaf type | Taxon | Leaf type |
| Abies alba | Coniferous | Pinus sylvestris | Coniferous |
| Acer campestre | Broadleaved | Platanus spec. | Broadleaved |
| Acer opalus | Broadleaved | Populus alba | Broadleaved |
| Acer platanoides | Broadleaved | Populus nigra s.l. | Broadleaved |
| Acer pseudoplatanus | Broadleaved | Populus spec. | Broadleaved |
| Aesculus hippocastanum | Broadleaved | Populus tremula | Broadleaved |
| Variable | Interpretation | Unit |
| bio01 | Annual mean temperature | ∘ C |
| bio02 | Mean diurnal temperature range | ∘ C |
| bio03 | Isothermality (bio02/bio07 100) | % |
| bio04 | Temperature seasonality (standard deviation 100) | ∘ C 100 |
| bio05 | Maximum temperature of the warmest month | ∘ C |
| bio06 | Minimum temperature of the coldest month | ∘ C |
| Modality | Model input | Encoder | Trainable |
| Aerial imagery | 280 280 3 crop, resized to 256 256 (28 m, 0.1 m GSD) | DINOv3 ViT-H+/16, blocks {7, 15, 23, 31}, 8 8 patch pooling | No |
| ALS | Points within 14 m radius, 0.1 m voxels | Sparse 3D U-Net (channel widths 32–160) + 64 learned queries | Yes |
| Satellite embedding | 3 3 64 AlphaEarth patch (10 m pixels, 30 m footprint) | Shared 2-layer MLP + positional embedding | Yes |
| Climate | 19 CHELSA BIOCLIM+ variables | 2-layer MLP | Yes |
| Setting | Value |
| Optimizer | AdamW |
| Initial learning rate | |
| Weight decay | 0.05 |
| Adam betas | (0.9, 0.999) |
| Scheduler | Cosine decay |
| Minimum learning rate |
| Variant | Leaf level | Species level | |||
| MAE | MAE | ||||
| AlphaEarth embeddings † | 0.093 0.003 | 0.803 0.017 | 0.048 0.001 | 0.485 0.038 | 0.585 0.014 |
| Sentinel-2, July composite | 0.098 0.004 | 0.791 0.015 | 0.050 0.001 | 0.464 0.044 | 0.588 0.029 |
| Sentinel-2, growing season | 0.099 0.003 | 0.786 0.018 | 0.050 0.001 | 0.458 0.052 | 0.579 0.041 |
| Sentinel-2, all year | 0.098 0.003 | 0.790 0.011 | 0.050 0.001 | 0.480 0.042 | 0.593 0.020 |
| Variant | Leaf level | Species level | |||
| MAE | MAE | ||||
| CHELSA BIOCLIM+ † | 0.093 0.003 | 0.803 0.017 | 0.048 0.001 | 0.485 0.038 | 0.585 0.014 |
| Climplicit embedding | 0.092 0.003 | 0.813 0.015 | 0.049 0.001 | 0.487 0.042 | 0.593 0.027 |
| Location MLP | 0.093 0.003 | 0.815 0.014 | 0.049 0.002 | 0.473 0.044 | 0.585 0.023 |
| Class set | Spatial 5-fold CV | Independent, plot | Independent, pixel |
| NFI ( ) | Forest reserve ( ) | Forest reserve ( ) | |
| All 11 classes (including chestnut) | 0.485 | a | n.e. b |
| Ten evaluable classes (excluding chestnut) | 0.483 | 0.733 | 0.552 |
| Class | MAE | RMSE | ||||||
| Plot | Pixel | Plot | Pixel | Plot | Pixel | Plot | Pixel | |
| Spruce | 0.085 | 0.115 | 0.172 | 0.237 | 0.668 | 0.677 | 89 | 264 |
| Fir | 0.034 | 0.038 | 0.071 | 0.092 | 0.829 | 0.643 | 48 | 100 |
| Pine | 0.044 | 0.063 | 0.135 | 0.194 | 0.718 | 0.421 | 50 | 102 |
| Larch | 0.020 | 0.029 | 0.058 | 0.085 | 0.869 | 0.728 | 23 | 61 |
| Arolla pine | 0.031 | 0.039 | 0.118 | 0.138 | 0.700 | 0.610 | 18 | 56 |
| Class | MAE | |||
| CV | Pixel-level | CV | Pixel-level | |
| Spruce | 0.111 | 0.115 | 0.746 | 0.677 |
| Fir | 0.061 | 0.038 | 0.668 | 0.643 |
| Pine | 0.024 | 0.063 | 0.599 | 0.421 |
| Larch | 0.039 | 0.029 | 0.742 | 0.728 |
| Arolla pine | 0.013 | 0.039 | 0.339 | 0.610 |
Tree Species Classification
Tree species classification identifies tree taxa from images, spectral measurements, 3D structure, and satellite or aerial observations.
See all 20 papers