Scalogram habitat measures as predictors of bird abundance

Birds select habitat characteristics, such as variability in habitat structure, across multiple spatial scales (grain and extent). Measuring habitat variability at multiple scales can better capture factors that influence avifauna communities than focusing on one scale only. One valuable tool in assessing habitat heterogeneity is the cumulative dynamic habitat index (DHI), which is derived from satellite data and captures temporal variability in vegetation productivity. Our goals were to develop new habitat measures from the cumulative DHI at multiple scales based on scalograms, and to test their performance in models of bird abundance. We counted birds at 188 plots during three breeding seasons (2007–2009) at Fort McCoy military installation, USA, to assess the abundance of forest (ovenbird), shrubland (indigo bunting), and grassland (grasshopper sparrow) bird specialists. We then calculated NDVI based on PlanetScope (3 m), Sentinel-2 (10 m), Landsat-8 (30 m), and MODIS (250 m) data to quantify cumulative DHI. We summarized the averaged NDVI cumulative DHI within multiple extents around each bird survey and developed 11 new habitat measures to test their predictive power in models of bird abundance. We found positive relationships between cumulative DHI at different extents and the abundances of both ovenbirds and indigo buntings, a forest and a shrubland specialist, respectively; and a negative relationship with grasshopper sparrows, a grassland specialist. In multiple linear regression models that incorporated single- and multi-grain predictors, the scalogram habitat measures explained moderate to high levels of variability in bird abundance, with R2 = 0.77, 0.37, and 0.75 for our forest, shrubland, and grassland specialists, respectively. Our results show that scalograms are an effective tool for capturing multiscale habitat configuration, because they capture the variability of habitat conditions in forests, shrublands, and grasslands. The scalogram habitat measures that we developed can be computed using our new R package ‘scalogram’.

File: Ecography-2025-Silveira-Scalogram-habitat-measures-as-predictors-of-bird-abundance.pdf

Substantial increases in burned area in circumboreal forests from 1983 to 2020 captured by the AVHRR record and a new autoregressive burned area detection algorithm

Wildfire maintains boreal forest health by catalyzing nutrient cycling and forest succession. However, increased annual burned area due to climate warming may facilitate forest loss and soil carbon release, which makes it important to monitor circumboreal burned area. Our goal was to characterize regional changes in circumboreal burned area from 1983 to 2020 using Advanced Very High Resolution Radiometer (AVHRR) data, and to identify the ecoregions where increases in burned area represent significant trends. We accomplished this by developing and applying a new burned area mapping algorithm that is based on an autoregressive analysis of the AVHRR and MODIS MOD09CMGv061 time series. Our algorithm worked well, and resulting burned area totals were similar to those of the MODIS MCD64A1v61 burned area product for years where both were available (2001􀀀 2020); however, the advantage of our AVHRR burned area dataset is that it extends back to 1983. Based on the resulting burned area maps, we evaluated circumboreal burned area changes and tested for significant trends. Net changes were substantial: while only 5.37 % of the circumboreal biome burned in the 1980s, 8.22 % did during the 2010s, an increase of 2.85 % (= 8.22–5.37 %) that corresponds to a proportional increase in area burned of 0.53 (= 2.85/5.37). In one ecoregion, Muskwa Slave Lake Forests, burned area more than quadrupled from the 1980s to the 2010s, and in three it more than tripled (Northern Canadian Shield Taiga, Yukon Interior Dry Forests, and Northeast Siberian Taiga). Furthermore, despite interannual variability in burned area typically being high, we found statistically significant increasing trends in burned area in seven of the twenty-three boreal ecoregions, corresponding to 19.6 % of boreal forests (35 % of North American and 11 % of Eurasian boreal forests), while only one ecoregion (Eastern Canadian Shield Taiga) had a decreasing trend. By analyzing the longterm AVHRR record, we were able to capture much larger increases in burned area than from the shorter MODIS record, allowing us to quantify how widespread and substantial these increases have been. By analyzing ecoregions, we found that north-eastern Siberian, north-western Canadian, and Alaskan boreal forests have experienced the most increases in burned area. These increases in burned area may have implications for future forest persistence and carbon storage within Eurasian and North American boreal forests.

File: Stephens-et-al.-2025-AVHRR-Boreal-Burned-Area-Mapping-1983-2020.pdf

Potential fence density in central and Western North America and implications for Bison (Bison bison) restoration

Fences serve multiple purposes, including livestock management, agriculture, property delineation, and conservation. However, fences often act as ecological barriers, limiting wildlife movement and access to resources, particularly for species like bison (Bison bison) in North America. Despite the substantial impacts of fencing, large-scale datasets on fence densities are lacking. Our goal was to create potential fence density maps for the western and central U.S. and Canada using GIS modelling and freely accessible data. Specifically, we aimed to: (1) map potential fence density and identify high density of fence, (2) contrast the potential fence density map with the patterns of high human influences, and (3) identify areas with high bison habitat suitability and low density of potential fences. Using GIS modelling, we generated potential fence density maps by integrating data on land parcels, croplands, roads, and railroads. Subsequently, we identified regions with high and low potential fence density and compared them with patterns of human influence and bison habitat suitability. We found high total potential fence density in central regions of Canada and the U.S., mainly due to agriculture and transportation corridors. Interestingly, areas with high potential density of fence in the western U.S., often had low other human influence, suggesting that human influence maps may underestimate impacts if they miss fences. We also identified large areas with high bison habitat suitability and low fence density, which are promising for bison restoration. Our findings highlight the importance of assessing fences for wildlife conservation and supporting bison restoration in the Great Plains.

File: silveira_FenceDensity_BioCons_2025.pdf

Impact of data density and endmember definitions on long-term trends in ground cover fractions across European grasslands

Long-term monitoring of grasslands is pivotal for ensuring continuity of many environmental services and for supporting food security and environmental modeling. Remote sensing provides an irreplaceable source of information for studying changes in grasslands. Specifically, Spectral Mixture Analysis (SMA) allows for quantification of physically meaningful ground cover fractions of grassland ecosystems (i.e., green vegetation, nonphotosynthetic vegetation, and soil), which is crucial for our understanding of change processes and their drivers. However, although popular due to straightforward implementation and low computational cost, ‘classical’ SMA relies on a single endmember definition for each targeted ground cover component, thus offering limited suitability and generalization capability for heterogeneous landscapes. Furthermore, the impact of irregular data density on SMA-based long-term trends in grassland ground cover has also not yet been critically addressed. We conducted a systematic assessment of i) the impact of data density on long-term trends in ground cover fractions in grasslands; and ii) the effect of endmember definition used in ‘classical’ SMA on pixel- and map-level trends of grassland ground cover fractions. We performed our study for 13 sites across European grasslands and derived the trends based on the Cumulative Endmember Fractions calculated from monthly composites. We compared three different data density scenarios, i.e., 1984–2021 Landsat data record as is, 1984–2021 Landsat data record with the monthly probability of data after 2014 adjusted to the pre-2014 levels, and the combined 1984–2021 Landsat and 2015–2021 Sentinel-2 datasets. For each site we ran SMA using a selection of sitespecific and generalized endmembers, and compared the pixel- and map-level trends. Our results indicated no significant impact of varying data density on the long-term trends from Cumulative Endmember Fractions in European grasslands. Conversely, the use of different endmember definitions led in some regions to significantly different pixel- and map-level long-term trends raising questions about the suitability of the ‘classical’ SMA for complex landscapes and large territories. Therefore, we caution against using the ‘classical’ SMA for remotesensing- based applications across broader scales or in heterogenous landscapes, particularly for trend analyses, as the results may lead to erroneous conclusions.

File: Lewinska-et-al.-2025-Impact-of-data-density-and-endmember-definitions.pdf

Landscape scale effects of primary productivity on forest bird species occurrence and abundance in Argentina

Context Approaches estimating landscape effects
on biodiversity frequently focus on a single extent,
finding one ‘optimal’ extent, or use narrow extents.
However, species perceive the environment in different
ways, select habitat hierarchically, and respond to
multiple selection pressures at extents that best predict
each pressure.
Objective We aimed to assess multi-scale relationships
between primary productivity and species
occurrences and abundances.
Methods We used a multi-scale approach, called
‘scalograms’, to assess landscape level effects of primary
productivity, in the form of Dynamic Habitat
Indices (DHIs) on the occurrences and abundances of
100 Argentinian forest bird species. We used average
DHI values within multiple extents (3 × 3 to 101 ×
101 pixels; 30 m resolution), and 11 ‘scalogram’ metrics
as environmental inputs in occurrence and abundance
models.
Results Average cumulative DHI values in extents
81 × 81 to 101 × 101 pixels (5.9 – 9.2 km2)
and maximum
cumulative DHI across extents were in the top
three predictors of species occurrences (included in
models for 41% and 18% of species, respectively).
Average cumulative DHI values in various extents
contributed ~ 1.6 times more predictive power to
occurrence models than expected. For species abundances,
average DHI values and scalogram measures
were in the top three predictors for < 2% of species and contributed less model predictive power than expected, regardless of DHI type (cumulative, minimum, variation). Conclusions Argentinian forest bird occurrences, but not abundances, respond to high levels of primary productivity at multiple, broad extents rather than a single ‘optimal’ extent. Factors other than primary productivity appear to be more important for predicting abundance.

File: Olah-et-al_2025_Argentina_DHIs_scalograms_bird-occurence_abundance.pdf

Medium-resolution Dynamic Habitat Indices from Landsat and Sentinel-2 satellite imagery

Biodiversity science requires effective tools to predict patterns of species diversity at multiple temporal and
spatial scales. The Dynamic Habitat Indices (DHIs) are remotely sensed indices that summarize aboveground
vegetation productivity in a way that is ecologically relevant for biodiversity assessments. Existing global DHIs,
derived from MODIS at 1-km resolution, predict species richness at broad scales well, but that resolution is coarse
relative to the grain at which many species perceive their habitat. With the much finer spatial resolution of
Sentinel-2 and Landsat data, plus Landsat’s longer data record, it is possible to track potential changes of
vegetation and its impacts on biodiversity at a finer grain over longer periods. Here, our main goals were to
derive the DHIs from 10-m Sentinel-2, 30-m Landsat, and 250-m MODIS data for the conterminous US and
compare all DHIs at two spatial extents, and to evaluate the ability of these DHIs to predict bird species richness
in 25 National Ecological Observatory Network terrestrial sites. In addition, we derived the Landsat DHIs for
1991–2000 and investigated how they changed by 2011–2020. We found that the Sentinel-2, Landsat, and
MODIS DHIs were highly correlated when summarized by ecoregion (Spearman correlation ranging from 0.89 to
0.99), indicating good agreement between them and that we were able to overcome the lower temporal resolution
of Sentinel-2 and Landsat. Sentinel-2 and Landsat DHIs outperformed MODIS in modeling species richness
for all bird guilds, explaining up to 49% of variance of grassland affiliates in linear regression models.
Furthermore medium-resolution DHIs (10–30 m resolution) captured spatial heterogeneity much better than
MODIS DHIs. We observed considerable changes in Landsat DHIs from 1991–2000 to 2011–2020, such as
increased cumulative DHI along the West Coast, in mountain ranges, and in the South, but lower cumulative DHI
in the Midwest. Our newly derived DHIs for the conterminous US have great potential for use in biodiversity
science and conservation.

File: Razenkova-etal-2025-DHIs-from-Landsat-and-Sentinal2.pdf

Conservation value and ecosystem service provision of Nothofagus antarctica forests based on phenocluster categories

Traditional approaches of forest classifications were based on tree species composition, but recently combine phenology
and climate to characterise functional (cyclic and seasonal greenness) rather than structural or compositional
components (phenoclusters). The objective was to compare the conservation value (capacity to support more native
biodiversity) and provision of ecosystem services (ES) in different phenocluster categories of Nothofagus antarctica
forests in Tierra del Fuego (Argentina). We used available models (ES, potential biodiversity) and ground-truth data of
145 stands, comparing phenocluster values using uni- and multivariate analyses. Conservation value and capacity to
supply ES significantly varied among phenocluster categories: (i) cultural, regulating, and provisioning ES and potential
biodiversity at landscape level, (ii) soil carbon and nitrogen, (iii) dominant height, crown cover, basal area, total volume,
and domestic animal stock, and (iv) understory plant richness and cover at stand level. These differences are linked to
the forest capacity to support more native biodiversity and ES. Besides, multivariate analyses supporting the split of
this forest type into four phenocluster subtypes (coast, highland, ecotone with other types, and degraded or secondary
forests). Our findings suggest the needs of specific management and conservation proposals, based on phenoclusters
rather than forest types defined by tree canopy-cover composition.

File: Martinez-Pastur-et-al.-2025_Conservation-value_Nothofagus-antarctica_based-on-phenoclusters.pdf

A half-century of land cover changes in the Caucasus derived from Corona spy satellite and Landsat images

Land cover change substantially affects ecosystems and leaves long-lasting legacies. Unfortunately, land cover analyses typically begin in the mid-1980s, when 30-m Landsat data became available, missing major global changes that occurred in the 1960s and 1970s. We aimed to quantify long-term land cover changes in the Caucasus (240,000 km2) comparing the magnitude of Soviet-era (1965–1987) versus post-Soviet changes (1987–2015). We (a) mapped land cover based on 1965 Corona spy satellite imagery and (b) quantified long-term changes by comparing 1965 Corona with 1987 and 2015 Landsatbased classifications while accounting for the differences in sensors’ spatial and spectral resolutions. Our Corona-derived map accuracy was 74.4 ± 3.7%, and change accuracies were 66.0 ± 4.2% for 1965–1987 and 61.6 ± 2.8% for 1965–2015. Overall, 30% of the land changed during the Soviet era compared to 20% during the post-Soviet era, highlighting the importance of mapping those early changes. Change trajectories differed considerably during the Soviet era and thereafter. For example, forests were lost during the Soviet era (− 6%) but gained area post-1987 (+ 5%). Croplands were often lost (− 18%) due to grassland gains (+ 11%), which were continuous, but at different rates (4% versus 7%), whereas croplands were lost in both eras, especially post-1987 (3% versus 16%). There were stark differences among countries: Azerbaijan underwent post-Soviet cropland gains, while the Russian Caucasus and Georgia experienced forest gains. Our results highlight the feasibility and value of early spy satellite data for long-term land cover change analyses, particularly in regions with substantial land cover changes then.

File: s10113-025-02360-6.pdf

Avian diversity across guilds in North America versus vegetation structure as measured by the Global Ecosystem Dynamics Investiation (GEDI)

Avian diversity, a key indicator of ecosystem health, is closely related to canopy structure. Most avian diversity models are based on either optical remote sensing or airborne lidar data, but the latter is limited to small study areas. The launch of the Global Ecosystem Dynamics Investigation (GEDI) instrument in 2018 has opened new avenues for exploring the influence of vegetation structure on avian diversity. To examine how direct measurements of canopy structural characteristics explain bird diversity across North America, we analyzed 18 GEDI metrics from 2019 to 2022, along with corresponding Breeding Bird Survey (BBS) counts and AVONET morphological data, analyzing effects across broad regions and at varying spatial extents. We grouped 440 bird species into 20 ecological guilds under six guild categories and employed random forest algorithms to model avian diversity across eight spatial extents (1, 2, 3, 4, 5, 10, 20, and 39.2 km). The models predicted six diversity indices, including species richness (sRich), functional richness (fRich), evenness (fEve), dispersion (fDis), divergence (fDiv), and redundancy (fRed) across eight spatial extents. The best-predicted guilds varied for each diversity index. The most accurate models were sRich (pseudo-R2 = 0.71, RMSE = 4.28) and fRed (pseudo-R2 = 0.60, RMSE = 0.13) for forest specialists guilds; fRich (pseudo-R2 = 0.55, RMSE = 0.18) for urban guilds; fEve (pseudo-R2 = 0.28, RMSE = 0.08) for insectivore guilds; and fDiv (pseudo-R2 = 0.38, RMSE = 0.12) and fDis (pseudo-R2 = 0.53, RMSE = 0.87) for short distance migrants guilds. Our results highlight the critical role of canopy structure, including its horizontal and vertical distribution and variation, in predicting avian diversity, as measured by the mean number of detected modes (num_detectedmodes), the standard deviation of foliage height diversity (FHD), num_detectedmodes, canopy cover, and plant area index (PAI) across the spatial extents centered on BBS routes. Therefore, we recommend incorporating the GEDI metrics into avian diversity modeling and mapping across North America, thereby potentially enhancing bird habitat management and conservation efforts.

File: Xu-1-s2.0-S0034425724004723-main.pdf

Croplands abandoned between 1986 and 2018 across the U.S.: spatiotemporal patterns and current land uses

Knowing where and when croplands have been abandoned or otherwise removed from cultivation is fundamental to evaluating future uses of these areas, e.g. as sites for ecological restoration, recultivation, bioenergy production, or other uses. However, large uncertainties remain about the location and time of cropland abandonment and how this process and the availability of associated lands vary spatially and temporally across the United States. Here, we present a nationwide, 30 m resolution map of croplands abandoned throughout the period of 1986–2018 for the conterminous United States (CONUS). We mapped the location and time of abandonment from annual cropland layers we created in Google Earth Engine from 30 m resolution Landsat imagery using an automated classification method and training data from the U.S. Department of Agriculture Cropland Data Layer. Our abandonment map has overall accuracies of 0.91 and 0.65 for the location and time of abandonment, respectively. From 1986 to 2018, 12.3 (±2.87) million hectares (Mha) of croplands were abandoned across CONUS, with areas of greatest change over the Ogallala Aquifer, the southern Mississippi Alluvial Plain, the Atlantic Coast, North Dakota, northern Montana, and eastern Washington state. The average annual nationwide abandoned area across our study period was 0.51 Mha per year. Annual abandonment peaked between 1997 and 1999 at a rate of 0.63 Mha year−1 , followed by a continuous decrease to 0.41 Mha year−1 in 2009–2011. Among the abandoned croplands, 53% (6.5 Mha) changed to grassland and pasture, 18.6% (2.28 Mha) to shrubland and forest, 8.4% (1.03 Mha) to wetlands, and 4.6% (0.56 Mha) to non-vegetated lands. Of the areas that we mapped as abandoned, 19.6% (2.41 Mha) were enrolled in the Conservation Reserve Program as of 2020. Our new map highlights the long-term dynamic nature of agricultural land use and its relation to various competitive pressures and land use policies in the United States.

File: Xie_2024_Environ._Res._Lett._19_044009.pdf