Zarvan All stories

Forests

Bark beetle prediction

We ranked forest in three areas of southern Sweden by how likely spruce bark beetles were to attack in the 2024 season, using data up to 2023. In forward tests, plots in the highest risk class were attacked about 3.8 times as often as the average plot.

Open the risk map
Satellite view of Östergötland, Jönköping and Kalmar with high-risk forest in red and the study areas outlined
Part of the map. Red marks the forest ranked riskiest within each study area; towns, farmland and water have no rating.
21%of field plots in the very high class attacked per year
3.8×the average attack rate (95% range 2.4–5.2)
4,243field plot-years used as labels, 288 of them attacked

The goal: where to look first

Nobody can inspect every stand of spruce. The map ranks forest so that inspection can start where an attack is most likely. It says where to look first, not how much forest is threatened.

How we built it

  1. Field truth

    SLU's inventory recorded freshly attacked spruce on field plots of 20–25 m radius from 2020 to 2024. Plot positions are confidential and do not appear on the map.

  2. What the forest looks like from above

    Sentinel-2 imagery, land cover, tree height, soil wetness, terrain, climate, nearby clear-cuts, beetle traps and 9.2 million AI-detected dead spruces from Skogsstyrelsen.

  3. More examples to learn from

    We added 4,200 dead-spruce events and 3,216 spruce controls at low weight. This raised the score on unseen Östergötland forest from 0.17 to 0.23. We call this model B5+.

  4. Tests on data the model had not seen

    Train outside Östergötland and test there; train on earlier years and test the next; compare with Skogsstyrelsen's published risk index and with later dead-spruce detections.

  5. The map

    Forest is ranked within each area and split into five classes, at 100 m resolution, using the October forecast for the 2024 season.

The classes separate risk

Attack rates fall steadily from the very high class to the low class. The percentages shift between years, so the multiples are the steadier guide.

Very hightop 5% of forest21% · 3.8×Highnext 5%12% · 2.1×Elevatednext 15%8% · 1.4×Averagenext 25%5% · 0.9×Lowlowest 50%3% · 0.6×
Share of plots attacked per year, and multiple of the average. Measured on 2,568 plot-years with 141 attacks in forward tests 2022–2024, outside Östergötland and without photos of the test year.

How it compares

B5+ clearly beat Skogsstyrelsen's risk index in the forward 2023 test. In Östergötland the two cannot be told apart.

Field testRisk indexB5+
Östergötland: AUPRC0.1820.227
Östergötland: attacks in the riskiest 10%24%29%
Forward 2023: AUPRC0.1100.248
Forward 2023: attacks in the riskiest 10%19%32%

AUPRC measures how well attacked plots rank at the top; higher is better. Östergötland has 418 plot-years with 38 attacks, and forward 2023 has 994 with 63. With no skill, the riskiest 10% would hold about 10% of attacks.

Simple field records and attack history did about as well as anything on Östergötland, and with only 38 attacks the ranges overlap widely. In forward tests, B5+ pulled ahead of the field model in 2022 and 2023 (2.5× and 3.9× better than chance, against 1.5× and 1.8×) but not in 2024, with its 12 attacks (2.9× against 3.6×).

0.00.10.20.30.40.5No skill (0.091)Field model + B5+ combined0.30Field records and history0.30B5+, the model behind the map0.23Sentinel-2 alone0.16Climate alone0.10
Test score (AUPRC) with 95% range, trained outside Östergötland and tested there.

How far to trust it

What comes next

Later forecast dates help only slightly: the riskiest 10% held 29% of attacks for an October forecast, 30% for April and 34% for June. That supports the main goal, a forecast map for 2027 issued before the season and scored when the 2027 field data arrive.

Open the risk map

Data

DataFromWhat it gives
SLU bark-beetle inventory (TFDI/NRS)SLUField labels: fresh attacks per plot
Sentinel-2 L2AESA / Planetary ComputerVegetation and moisture indices, 3×3 pixel window around each plot
Sentinel-1 RTCESA / Planetary ComputerRadar backscatter; tested, then left out of B5
ESA WorldCover 2021ESATree cover share
Döda granarSkogsstyrelsenAI-detected dead spruces from aerial photos, all of Sweden 2017–2024 (9.2 million trees): counts within 250 m, 1 km, 5 km
ERA5-LandCopernicus CDSApril–September temperature, degree-days, precipitation, water balance, summer soil moisture vs 1991–2020
NMD 2018 land coverNaturvårdsverketSpruce, conifer and clear-felled shares within 50 m and 250 m
NMD tree height (laser)NaturvårdsverketTree height classes (5–45 m) within 20 m; share of trees over 20 m
Barkborre svärmningsrapportSkogsstyrelsenWeekly beetle trap catches; only from 2021, so it can help the 2022–2024 labels only
Utförd avverkningSkogsstyrelsenFinal fellings 2015–2024 (620,000): clear-cut area within 500 m and count within 1 km
Copernicus DEM GLO-30ESA / Planetary ComputerElevation, slope, south exposure
MarkfuktighetsindexSLU / Metria via NaturvårdsverketSoil wetness index (10 m) within 20 m

Boundaries from SCB. Satellite background: contains modified Copernicus Sentinel data 2023, via Microsoft Planetary Computer.