We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a two multi-band raster stacks. The raster stacks are identified by the vegetation index they represent (i.e., NDVI, TC greenness). Each band within the separate raster stacks represents a month from January 2014 through December 2020 (i.e., band 1 is January 2014 and band 84 is December 2020). Though the full study extends until December 2021, we chose to remove data from 2021 from the phenological analysis due to increased fire activity during 2021. Additionally, we chose not to include the TC transformation metrics of brightness and wetness to focus on metrics to simply measure vegetation greenness.

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a single XML metadata file and four zipped files containing the raster stacks, where each zipped file contains the rasters for each vegetation index, respectively (i.e., NDVI, TC brightness, TC greenness, TC wetness). The raster stacks within each index-specific zipped folder are identified by the season they represent (i.e., spring, late-spring, summer, fall). Each band within the separate raster stacks (n=37) represents a year from 1985 through 2021 (i.e., band 1 is 1985 and band 37 is 2021).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a two multi-band raster stacks. The raster stacks are identified by the vegetation index they represent (i.e., NDVI, TC greenness). Each band within the separate raster stacks represents a month from January 2014 through December 2020 (i.e., band 1 is January 2014 and band 84 is December 2020). Though the full study extends until December 2021, we chose to remove data from 2021 from the phenological analysis due to increased fire activity during 2021. Additionally, we chose not to include the TC transformation metrics of brightness and wetness to focus on metrics to simply measure vegetation greenness.

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a single XML metadata file and four zipped files containing the raster stacks, where each zipped file contains the rasters for each vegetation index, respectively (i.e., NDVI, TC brightness, TC greenness, TC wetness). The raster stacks within each index-specific zipped folder are identified by the season they represent (i.e., spring, late-spring, summer, fall). Each band within the separate raster stacks (n=37) represents a year from 1985 through 2021 (i.e., band 1 is 1985 and band 37 is 2021).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). For each index, we calculate the Sen's Slope across a series of climate periods (see associated Child Item - 1) Database of Trends in Vegetation Properties and Climate Adaptation Variables -- Standardized Precipitation Evapotranspiration Index Timeseries for the Upper Gila River Watershed: 1985 to 2021). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a single XML metadata file and a zipped file containing the raster stacks. The raster stacks are identified by both the season they represent (i.e., spring, late-spring, summer, fall) and the vegetation index (i.e., NDVI, TC brightness, TC greenness, TC wetness). Each band represents Sen's slope value across each climate period, where band 1 represents the Sen's slope across the 1st climate period (i.e., 1985 through 1993), band 2 represents the Sen's slope across the 2nd climate period (i.e., 1993 through 2014), and band 3 represents the Sen's slope across the 3rd climate period (i.e. 2014 through 2021).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). For each index, we calculate the Sen's Slope across a series of climate periods (see associated Child Item - 1) Database of Trends in Vegetation Properties and Climate Adaptation Variables -- Standardized Precipitation Evapotranspiration Index Timeseries for the Upper Gila River Watershed: 1985 to 2021). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a single XML metadata file and a zipped file containing the raster stacks. The raster stacks are identified by both the season they represent (i.e., spring, late-spring, summer, fall) and the vegetation index (i.e., NDVI, TC brightness, TC greenness, TC wetness). Each band represents Sen's slope value across each climate period, where band 1 represents the Sen's slope across the 1st climate period (i.e., 1985 through 1993), band 2 represents the Sen's slope across the 2nd climate period (i.e., 1993 through 2014), and band 3 represents the Sen's slope across the 3rd climate period (i.e. 2014 through 2021).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developing using satellite imagery, including the Normalized Difference Vegetation Index (NDVI) and the Tasseled Cap (TC) Transformation. The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. The TC transformation is an approached used to transform satellite imagery into a collection of spectral metrics that can quantify various aspects of the vegetation and soil surfaces (Kauth and Thomas, 1976). Specifically, the TC transformation develops 6 separate metrics, though we only assess the three primary metrics: (i) brightness (transformation 1), (ii) greenness (transformation 2), and (iii) wetness (transformation 3). The TC transformation metrics are calculated using a series of coefficients multiplied across reflectance values for the suite of Landsat bands, then summed across each metric. No specific range is identified for the TC transformation metrics, though for both brightness and greenness, positive values represent brighter and greener conditions, respectively, while negative values represent wetter conditions for wetness. Because bandwidths differ slightly between Landsat 4, 5, 7 and Landsat 8, we use two sets of coefficients and complete the calculation separately before combining the collections into a single series of images (DeVries et al., 2016; Zhai et al., 2022). For each index, we calculate the Sen's Slope across a series of climate periods (see associated Child Item - 1) Database of Trends in Vegetation Properties and Climate Adaptation Variables -- Standardized Precipitation Evapotranspiration Index Timeseries for the Upper Gila River Watershed: 1985 to 2021). All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a single XML metadata file and a zipped file containing the raster stacks. The raster stacks are identified by both the season they represent (i.e., spring, late-spring, summer, fall) and the vegetation index (i.e., NDVI, TC brightness, TC greenness, TC wetness). Each band represents Sen's slope value across each climate period, where band 1 represents the Sen's slope across the 1st climate period (i.e., 1985 through 1993), band 2 represents the Sen's slope across the 2nd climate period (i.e., 1993 through 2014), and band 3 represents the Sen's slope across the 3rd climate period (i.e. 2014 through 2021).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developed using satellite imagery, including the Normalized Difference Vegetation Index (NDVI). The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a multi-band raster stack of monthly NDVI images from January 2014 through July 2022 covering the area of the Bylas Fire. We included data from 2022, contrasting the full study which only includes data through 2021, to include additional data regarding our post-fire vegetation response analysis. Each band within the raster stack represents a month from 2014 through 2022 (i.e., band 1 is January 2014 and band 103 is July 2022).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. For this study, the vegetation conditions are characterized using a series of remote sensing vegetation indices developed using satellite imagery, including the Normalized Difference Vegetation Index (NDVI). The NDVI is a commonly used vegetation index that quantifies relative greenness of the vegetation based on the plant’s photosynthetic activity, measured as a ratio between the Near Infrared (NIR) and Red bands (Tucker, 1979). The NDVI equation follows: NDVI = (NIR band - Red band) / (NIR band + Red band). NDVI has a range of -1 to 1, though green vegetation theoretically ranges from 0 to 1. Dense green vegetation is represented with values closer to 1 while barren soil, rock, and less-dense surface vegetation has values closer to 0. Values below 0 often represent water due to its unique reflective characteristics. All raster products were developed using the Google Earth Engine (GEE) cloud computing software program for the Upper Gila River watershed. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. This Child Item consists of a multi-band raster stack of monthly NDVI images from January 2014 through July 2022 covering the area of the Bylas Fire. We included data from 2022, contrasting the full study which only includes data through 2021, to include additional data regarding our post-fire vegetation response analysis. Each band within the raster stack represents a month from 2014 through 2022 (i.e., band 1 is January 2014 and band 103 is July 2022).

We apply a research approach that can inform riparian restoration planning by developing products that show recent trends in vegetation conditions identifying areas potentially more at risk for degradation and the associated relationship between riparian vegetation dynamics and climate conditions. The full suite of data products and a link to the associated publication addressing this analysis can be found on the Parent data release. To characterize the climate conditions across the study period, we use the Standardized Precipitation Evapotranspiration Index (SPEI). The SPEI is a water balance index which includes both precipitation and evapotranspiration in its calculation. Conditions from the prior n months, generally ranging from 1 to 60, are compared to the same respective period over the prior years to identify the index value (Vicente-Serrano et al., 2010). Values generally range from -3 to 3, where values less than 0 suggest drought conditions while values greater than 0 suggest wetter than normal conditions. For this study, we are using the 12-month, or 1-year, SPEI to compare annual conditions within the larger Upper Gila River watershed. The SPEI data was extracted into a CSV spreadsheet using data from the Gridded Surface Meteorological (GRIDMET) dataset, which provides a spatially explicit SPEI product in Google Earth Engine (GEE) at a 5-day interval and a spatial resolution of 4-km (Abatzoglou, 2013). In GEE, we quantify overall mean values of SPEI across each 5-day period for the watershed from January 1980 to December 2021. Using R software, we reduced the 5-day values to represent monthly mean values and constrained the analysis to water year 1980 (i.e., October 1980) through water year 2021 (i.e., October 2021). Using the monthly timeseries, we completed the breakpoint analysis in R to identify breaks within the SPEI time series. The algorithm identifies a seasonal pattern within the timeseries. When the seasonal pattern deviates, a breakpoint is then detected. These breaks can be used to pinpoint unique climate periods in the time series. This is a Child Item for the Parent data release, Mapping Riparian Vegetation Response to Climate Change on the San Carlos Apache Reservation and Upper Gila River Watershed to Inform Restoration Priorities: 1935 to Present - Database of Trends in Vegetation Properties and Climate Adaptation Variables. The spreadsheet attached to this Child Item consists of 5 columns, including the (i) month from January 1985 through October 2021, (ii) the 1-year SPEI monthly time series, (iii) the dates identified as breaks within the breakpoint algorithm, (iv) the breakpoint trend identified within the breakpoint algorithm, and (v) the dates that were used as the climate period breaks in this study. The climate periods identified in this spreadsheet using the SPEI data were used as the climate periods in our riparian study.