Romania Forest Lab · Whole massif study · 28 October 2026

Forest area and canopy give different answers in Retezat.

ForTy maps less forest in both massifs during 2017–2025. Copernicus maps a small crown-area decrease in Făgăraș and an increase in Retezat during 2018–2024. These are different periods and measurements.

Mapped estimates · EPSG:3035 · fixed geographical boundaries · no causal attribution

The same annual beta product measures both directions. Net change equals gains minus losses; the two massifs share one vertical scale.
Fig. 1 The same annual beta product measures both directions. Net change equals gains minus losses; the two massifs share one vertical scale.
Massif2017 forest (ha)Lost (ha)Gained (ha)Net (ha)Net / baseline
Făgăraș116,075.721,310.62292.18-1,018.44-0.877%
Retezat31,839.06131.5833.96-97.62-0.307%

These totals measure changes in modelled forest classes. They are not validated inventory estimates, canopy-density measurements, or proof of permanent deforestation.

Endpoint change and temporary disturbance answer different questions.

Lines show mapped area changes from each massif’s 2017 baseline. Solid lines include planted and plantation forest; dashed lines include only natural forest.
Fig. 2 Lines show mapped area changes from each massif’s 2017 baseline. Solid lines include planted and plantation forest; dashed lines include only natural forest.
MassifNatural loss (ha)Natural gain (ha)Natural → planted (ha)Forest recovered by 2025 (ha)
Făgăraș1,196.40275.690.01185.88
Retezat102.2025.830.0011.25

A pixel forested in both endpoints can lose its forest label and recover between them. Recovery totals follow cells forested in 2017; later forest gains are represented in annual turnover. Recovery here means a return to a modelled forest class. It does not establish mature canopy, ecological restoration, or recovery of biodiversity. Natural-to-planted conversion is included in natural-forest loss, but is not total-forest loss.

The measurements follow the massif polygons.

Făgăraș: forest classes 1/2/3/4/7, clipped to published relief polygon 552. Loss and gain refer to the 2017 and 2025 endpoints; recovery records intervening model changes.
Fig. 3 Făgăraș: forest classes 1/2/3/4/7, clipped to published relief polygon 552. Loss and gain refer to the 2017 and 2025 endpoints; recovery records intervening model changes.
Retezat: forest classes 1/2/3/4/7, clipped to published relief polygon 576. Loss and gain refer to the 2017 and 2025 endpoints; recovery records intervening model changes.
Fig. 4 Retezat: forest classes 1/2/3/4/7, clipped to published relief polygon 576. Loss and gain refer to the 2017 and 2025 endpoints; recovery records intervening model changes.

The figures use analytical rasters, not browser map tiles. Download the 10 m GeoTIFFs to inspect small changes that may be difficult to see at report scale.

Changing the screening rule changes which hectares qualify.

MassifRuleScore ≥Lost (ha)Gained (ha)Net (ha)Excluded (ha)
FăgărașRaw endpoints0.01,310.62292.18-1,018.440.00
FăgărașRaw endpoints0.51,275.85286.55-989.30854.40
FăgărașRaw endpoints0.7534.08183.16-350.9226,226.75
FăgărașRaw endpoints0.9No eligible cells——155,814.93
FăgărașStable endpoints0.01,023.81288.97-734.84356.89
FăgărașStable endpoints0.5992.07277.76-714.311,239.11
FăgărașStable endpoints0.7354.77129.50-225.2726,567.59
FăgărașStable endpoints0.9No eligible cells——155,814.93
RetezatRaw endpoints0.0131.5833.96-97.620.00
RetezatRaw endpoints0.5129.4231.54-97.88214.06
RetezatRaw endpoints0.769.3614.46-54.907,772.76
RetezatRaw endpoints0.9No eligible cells——41,297.73
RetezatStable endpoints0.096.2732.29-63.9840.79
RetezatStable endpoints0.592.5328.84-63.69261.86
RetezatStable endpoints0.737.418.02-29.397,824.73
RetezatStable endpoints0.9No eligible cells——41,297.73

Raw endpoints use all common valid cells. Stable endpoints require the same binary forest status in 2017/2018 and in 2024/2025. Score filters apply to the maximum native-class score at both endpoints, or all four years for stable endpoints. Excluded cells contribute to neither loss nor gain. Each row has its own eligible support; the variations are not statistical confidence intervals.

Canopy disturbance stays separate from forest-type balance.

MassifHansen loss (ha)Both products (ha)ForTy loss only (ha)Hansen loss only (ha)
Făgăraș2,545.691,102.83207.791,442.86
Retezat180.1285.3446.2494.78

Hansen years 2018–2025 are compared on an aligned 30 m equal-area grid using fractional 10 m support. The baseline filter is ForTy forest in 2017; no year-2000 canopy-density threshold was applied. These counts therefore differ from default Global Forest Watch statistics. Loss can reflect harvest, mortality, fire, or other disturbance. Repeated losses at one pixel are not fully represented.

Hansen overlap measures cross-product agreement, not independent validation. Hansen detects canopy disturbance while ForTy predicts land-use types; differences can be legitimate. Their loss and gain totals cannot be mixed into a net figure.

Copernicus compares canopy during 2018–2024.

Changes in crown-covered area mostly cancel in Făgăraș; Retezat shows a net increase. These are density-weighted hectares, not hectares of newly forested land.
Fig. 5 Changes in crown-covered area mostly cancel in Făgăraș; Retezat shows a net increase. These are density-weighted hectares, not hectares of newly forested land.
MassifPeriodCanopy lost (ha)Canopy gained (ha)Canopy net (ha)ForTy forest net (ha)
fagaras2018–2021501.6039.88-461.72-452.45
fagaras2021–2024239.4477.12-162.32-259.41
retezat2018–202123.522.52-21.00-28.72
retezat2021–20249.00181.84+172.84-33.51
MassifCrown decrease (ha)Crown increase (ha)Net crowns (ha)
fagaras2,974.292,899.19-75.10
retezat483.49778.46+294.97

Presence change measures tree-covered extent on a 20 m grid. The ForTy column uses matching years, but different definitions: forest land use can persist through temporary canopy removal. Density measures vertical crown projection on 10 m cells: sum(density/100 × 0.01 ha). Positive and negative density differences are spatial changes between 2018 and 2024, not cumulative annual events. The presence-change product is marked not validated by Copernicus.

Every total has fixed support and a reproducible definition.

Boundaries follow the 2008 Romanian relief dataset after Posea and Badea’s geographical regionalization. They are generalized massif polygons, not cadastral or park boundaries. The WGS84 download omits its projection file; the publisher’s GCS_WGS84 designation is used explicitly.

Forest comprises native beta classes 1, 2, 3, 4, and 7. Natural forest comprises 1, 2, and 7. Tree crops remain separate. Each 10 m cell represents 0.01 ha. The common mask requires valid observations in all nine years, with at least 99% coverage. Pixel centres decide polygon membership; raster and vector areas differ slightly at edges.

All years use ForTy v2 beta. No v1 baseline is spliced into the series. Sorted first-valid tile ownership prevents duplicate area. Sources are cached with hashes and generations where supplied. Scores remain uncalibrated model confidence; they are not canopy fractions. Classification changes and temporarily unstocked forest can resemble loss or recovery.

Rerun: uv run --project research forest-lab balance. After all source files are cached: FOREST_LAB_OFFLINE=1 uv run --project research --offline forest-lab balance. Focused checks: uv run --project research --offline pytest -q research/tests/test_balance.py.

Change rasters: 0 outside/excluded; 1 nonforest retained; 2 forest retained; 3 endpoint loss; 4 endpoint gain; 5 forest lost and recovered by 2025. Each annual raster band preserves native classes 1–7.

Sources and interpretation limits

No field verification, causal attribution, legal determination, or local accuracy estimate is included. Layout uses the lab’s bundled Visual Explainer shell under its existing MIT license.