Blue oak landscape classification outputs based on multi-temporal satellite observations and species distribution models
Data files
Jun 23, 2026 version files 22.22 MB
-
00_custom_functions.R
3.12 KB
-
00_functions.R
10.86 KB
-
01_landsat_stack_subset.R
19.36 KB
-
01_sentinel_unzip.R
2.76 KB
-
02_landsat_unmix_dyecreek.R
35.90 KB
-
02_landsat_unmix_randall.R
38.61 KB
-
02_sentinel_process.R
16.32 KB
-
03_sentinel_unmix_dyecreek.R
17.96 KB
-
03_sentinel_unmix_randall.R
34.26 KB
-
04_sdm_integration_env_analysis.R
55.25 KB
-
05_dyecreek_results.R
24.13 KB
-
06_veg_trends.R
10.90 KB
-
DyeCreek_Landsat_ecological_classes.tif
116.80 KB
-
DyeCreek_Landsat_TMR3.tif
653.23 KB
-
DyeCreek_Sentinel_ecological_classes.tif
777.12 KB
-
DyeCreek_Sentinel2_TMR3.tif
5.68 MB
-
landsat_functions.R
12.99 KB
-
Randall_Landsat_ecological_classes.tif
229.46 KB
-
Randall_Landsat_TMR3.tif
1.33 MB
-
Randall_Sentinel_ecological_classes.tif
1.58 MB
-
Randall_Sentinel2_TMR3.tif
11.56 MB
-
README.md
8.18 KB
Abstract
This dataset contains spatial outputs from an approach that integrates vegetation resilience metrics with climate-based projections of future habitat suitability for blue oak (Quercus douglasii) in California. The approach combines temporal mixture residual (TMR) analysis of multi-decadal satellite time series with species distribution models (SDMs) forecasting mid-21st century habitat suitability under a high-emissions climate scenario (RCP 8.5, 2040–2069), averaged for two global climate models (CNRM-CM5 and HadGEM2-ES).
Data were generated for two conservation properties managed by The Nature Conservancy: Randall Preserve (Kern County) and Dye Creek Preserve (Tehama County).
The outputs are provided as GeoTIFF raster files (.tif), grouped into two categories:
- Land classification rasters (4 files):
- Landsat-based classification for Randall Preserve
- Landsat-based classification for Dye Creek Preserve
- Sentinel-2–based classification for Randall Preserve
- Sentinel-2–based classification for Dye Creek Preserve
Each classification raster contains five discrete classes, representing combinations of historical vegetation resilience (derived from TMR residuals) and projected climatic suitability for blue oak habitat. These classes can be used to identify zones of high future potential, low stability, or mismatch between historical resilience and future climate suitability.
- Temporal mixture residual (TMR) rasters (4 files):
-
Landsat-based TMR residuals for Randall Preserve
-
Landsat-based TMR residuals for Dye Creek Preserve
-
Sentinel-2–based TMR residuals for Randall Preserve
-
Sentinel-2–based TMR residuals for Dye Creek Preserve
The TMR rasters provide continuous values that quantify deviations from expected multi-decadal vegetation dynamics. Higher residuals indicate greater instability in vegetation trajectories, while lower residuals indicate more stable, resilient vegetation dynamics.
-
Together, these datasets enable users to explore how historical vegetation resilience interacts with projected climate suitability at multiple spatial resolutions (30 m for Landsat, 10 m for Sentinel-2). They provide spatially explicit guidance for conservation planning and management in blue oak woodlands, highlighting areas where conservation investment may be most effective under future climate change.
Dataset DOI: 10.5061/dryad.v6wwpzh8s
Description of the data and file structure
To classify oak woodlands at the Randall and Dye Creek Preserves, we combined satellite-derived vegetation resilience metrics with climate-driven forecasts of blue oak (Quercus douglasii) habitat suitability (Rose et al., 2023). Future suitability was projected for 2040–2069 under a high-emissions scenario (RCP 8.5), averaged across two climate models (CNRM-CM5 and HadGEM2-ES). Vegetation resilience was quantified using temporal mixture residual (TMR) analysis of multi-decadal Landsat (30 m) and Sentinel-2 (10 m) imagery.
We applied k-means clustering to integrate future habitat suitability with the TMR metrics, producing five landscape classes per site and satellite platform. The optimal number of clusters was determined by evaluating within-cluster sum of squares (WSS), and cluster labels were reordered to maximize interpretability between data sources within each site.
Files and variables
File: DyeCreek_Sentinel_ecological_classes.tif
Description: Raster map of five ecological landscape classes for the Dye Creek Preserve derived from 10 m Sentinel-2 imagery. Classes integrate historical vegetation resilience with projected mid-21st century blue oak habitat suitability.
- Class 1: High suitability; low cover, high interannual variability with recent increases
- Class 2: Low–moderate suitability; low cover, high interannual variability with recent increases
- Class 3: Moderate suitability; stable cover
- Class 4: Mixed suitability; high historical cover with losses during drought and Park Fire
- Class 5: Low suitability; mixed historical conditions
File: DyeCreek_Landsat_ecological_classes.tif
Description: Raster map of five ecological landscape classes for the Dye Creek Preserve derived from 30 m Landsat imagery. Classes integrate historical vegetation resilience with projected mid-21st century blue oak habitat suitability.
- Class 1: High suitability; interannual variability, increasing recently
- Class 2: Low–moderate suitability; interannual variability, increasing recently
- Class 3: Moderate suitability; stable vegetation cover
- Class 4: Mixed suitability; drought-affected with some recovery
- Class 5: Low suitability; mixed historical conditions
File: Randall_Landsat_ecological_classes.tif
Description: Raster map of five ecological landscape classes for the Randall Preserve derived from 30 m Landsat imagery. Classes integrate historical vegetation resilience with projected mid-21st century blue oak habitat suitability.
- Class 1: High suitability; increasing vegetation cover
- Class 2: High suitability; drought-affected vegetation cover but recovering
- Class 3: Moderate suitability; stable or moderately increasing vegetation cover
- Class 4: Moderate suitability; drought-affected with recovery
- Class 5: Low suitability; mixed historical conditions
File: Randall_Sentinel_ecological_classes.tif
Description: Raster map of five ecological landscape classes for the Randall Preserve derived from 10 m Sentinel-2 imagery. Classes integrate historical vegetation resilience with projected mid-21st century blue oak habitat suitability.
- Class 1: High suitability; high vegetation cover and recent increases
- Class 2: High suitability; instability during drought
- Class 3: Moderate suitability; high, increasing cover
- Class 4: Moderate suitability; instability during drought
- Class 5: Low suitability; mixed historical conditions
File: Randall_Landsat_TMR3.tif
Description: Raster of the third principal component of the temporal mixture model residuals (TMR) from Landsat imagery for the Randall Preserve, representing multi-decadal vegetation change and resilience patterns (30 m, 1982-2025).
File: DyeCreek_Landsat_TMR3.tif
Description: Raster of the third principal component of the temporal mixture model residuals (TMR) from Landsat imagery for the Dye Creek Preserve, representing multi-decadal vegetation change and resilience patterns (30 m, 1984-2025).
File: DyeCreek_Sentinel2_TMR3.tif
Description: Raster of the third principal component of the temporal mixture model residuals (TMR) from Sentinel-2 imagery for the Dye Creek Preserve, representing vegetation change and resilience patterns (10 m, 2015-2025).
File: Randall_Sentinel2_TMR3.tif
Description: Raster of the third principal component of the temporal mixture model residuals (TMR) from Sentinel-2 imagery for the Randall Preserve, representing vegetation change and resilience patterns (10 m, 2015-2025).
Usage notes
All raster files (.tif) are standard GeoTIFFs and can be opened with any GIS or remote sensing software, including:
- QGIS (free, open-source): https://qgis.org
- R with the
terraorrasterpackages (free, open-source): https://cran.r-project.org - Python with
rasterioorGDAL(free, open-source) - ArcGIS (proprietary, licensed)
All .R scripts require R (https://www.r-project.org, free, open-source) and the packages listed in the Code/software section below. Scripts are numbered in the order they should be run, as outputs from earlier scripts are used as inputs to later scripts.
Code/software
All code used for the analysis can be found in the following files:
- 00_functions.R: General-purpose helper functions used across multiple downstream scripts (e.g., file I/O, raster handling utilities).
- 00_custom_functions.R: Project-specific custom functions supporting the unmixing and classification workflow.
- 01_landsat_stack_subset.R: Builds multi-temporal Landsat image stacks and subsets them to each preserve's spatial extent.
- 01_sentinel_unzip.R: Unzips and organizes raw Sentinel-2 imagery archives prior to processing.
- 02_landsat_unmix_randall.R: Performs temporal mixture residual (spectral unmixing) analysis on the Landsat time series for Randall Preserve.
- 02_landsat_unmix_dyecreek.R: Performs temporal mixture residual (spectral unmixing) analysis on the Landsat time series for Dye Creek Preserve.
- 02_sentinel_process.R: Pre-processes raw Sentinel-2 imagery (e.g., cloud masking, mosaicking, band stacking) prior to unmixing.
- 03_sentinel_unmix_randall.R: Performs temporal mixture residual analysis on the Sentinel-2 time series for Randall Preserve.
- 03_sentinel_unmix_dyecreek.R: Performs temporal mixture residual analysis on the Sentinel-2 time series for Dye Creek Preserve.
- 04_sdm_integration_env_analysis.R: Integrates species distribution model (SDM) habitat suitability projections with TMR outputs and conducts associated environmental analyses for both preserves.
- 05_dyecreek_results.R: Generates the k-means landscape classification, summary statistics, and figures for Dye Creek Preserve results.
- 06_veg_trends.R: Plots and analyzes long-term vegetation trends from the TMR time series across both preserves.
- landsat_functions.R: Helper functions specific to Landsat data processing and unmixing steps.
The majority of the analyses and visualizations were conducted using open source R packages (e.g., terra, dplyr, and ggplot2).
Selection of the temporal end members used to construct the temporal mixture models was done using ENVI Classic software.
Access information
Other publicly accessible locations of the data:
- Species distribution models used in the analysis are published here: https://datadryad.org/dataset/doi:10.5061/dryad.gf1vhhmw6
Data was derived from the following sources:
- Landsat data can be downloaded free of charge using USGS Earth Explorer https://earthexplorer.usgs.gov
- Sentinel-2 data can be downloaded free of charge using the Copernicus browser
