From imagery to ecology: Leveraging time series of all available Landsat observations to map and monitor ecosystem state and dynamics

https://zslpublications.onlinelibrary.wiley.com/doi/full/10.1002/rse2.24Published versio

United States Forest Disturbance Trends Observed Using Landsat Time Series

Disturbance events strongly affect the composition, structure, and function of forest ecosystems; however, existing U.S. land management inventories were not designed to monitor disturbance. To begin addressing this gap, the North American Forest Dynamics (NAFD) project has examined a geographic sample of 50 Landsat satellite image time series to assess trends in forest disturbance across the conterminous United States for 1985-2005. The geographic sample design used a probability-based scheme to encompass major forest types and maximize geographic dispersion. For each sample location disturbance was identified in the Landsat series using the Vegetation Change Tracker (VCT) algorithm. The NAFD analysis indicates that, on average, 2.77 Mha/yr of forests were disturbed annually, representing 1.09%/yr of US forestland. These satellite-based national disturbance rates estimates tend to be lower than those derived from land management inventories, reflecting both methodological and definitional differences. In particular the VCT approach used with a biennial time step has limited sensitivity to low-intensity disturbances. Unlike prior satellite studies, our biennial forest disturbance rates vary by nearly a factor of two between high and low years. High western US disturbance rates were associated with active fire years and insect activity, while variability in the east is more strongly related to harvest rates in managed forests. We note that generating a geographic sample based on representing forest type and variability may be problematic since the spatial pattern of disturbance does not necessarily correlate with forest type. We also find that the prevalence of diffuse, non-stand clearing disturbance in US forests makes the application of a biennial geographic sample problematic. Future satellite-based studies of disturbance at regional and national scales should focus on wall-to-wall analyses with annual time step for improved accuracy

LANDSCAPE SCALE SPECTRAL-TEMPORAL MODELLING OF BAMBOO-DOMINATED FOREST SUCCESSION WITHIN THE ATLANTIC FOREST OF SOUTHERN BRAZIL

Tropical and subtropical ecosystems have become vulnerable to biological invasion (i.e., bamboo) due to human induced forest fragmentation. Bamboo ecological processes have been found to impede forest development, resulting in a state of arrested succession, which has been found to significantly reduce biodiversity, thus contributing to biotic homogenization. In this study we use a semi-empirical approach to develop a community-level spatially explicit ecological process model (hybrid model) using a time-series of Landsat imagery to describe single-landscape scale ecological processes of a pervasive bamboo species (Merostachys skvortzovii) found throughout the Araucaria forest, a critically threatened subtype of Atlantic forest of southern Brazil. The model is subsequently used to map bamboo spatial distribution at a multiple-landscape scale to examine patch pattern throughout a portion of the Araucaria forest. It was determined that the M. skvortzovii lifecycle is a synchronized process occurring at single and multiple-landscapes scale and is comprised of four broad lifecycle phases: pioneer predominance, mature bamboo, dieback and pioneer regeneration. Bamboo patch pattern was found to be associated with human settlement and geographic features, with clusters of patches sharing the same shape and size observed at multiple scales

Multi-temporal Forest Cover Change and Forest Density Trend Detection in a Mediterranean Environment

The loss of forests along with the various types of shrubs in the Mediterranean region is seen as an important driver of climate change and has been repeatedly related with the observed land degradation and desertification in the region. Nevertheless, the extent of woody perennial vegetation cover (WPVC) and its density remain largely unclear. Here, we apply a series of algorithms and methods operationally used in Australia for large-scale WPVC mapping and monitoring and demonstrate their applicability in the Mediterranean region using a Spanish area as the trial site. Five Landsat TM and ETM+ images from various dates spanning 14 years are used to map changes in the extent of WPVC and to identify areas with a declining, stabilising or recovering trend. Results show that the applied methodology, which incorporates (i) preprocessing of the Landsat imagery, (ii) a canonical variate analysis to spectrally discriminate between woody and non-woody land cover types, (iii) a conditional probability network and (iv) spectral indices for mapping woody cover and density trend, is highly successful and well suited for use in Mediterranean environments. A rigorous accuracy assessment is undertaken producing overall accuracies above 97% for both woody and non-woody cover types and all dates. Results also show that in the area of study, the majority of WPVC disturbances were due to forest fires, which represent the region's most frequent natural and anthropogenic disturbance. This raises significant concerns about the future of the area's WPVC. Regeneration compensated to some degree for the high disturbance rates. Copyright © 2015 John Wiley & Sons, Ltd

Impacts of disturbance history on forest carbon stocks and fluxes: Merging satellite disturbance mapping with forest inventory data in a carbon cycle model framework

Forest carbon stocks and fluxes are highly dynamic following stand-clearing disturbances from severe fire and harvest and this presents a significant challenge for continental carbon budget assessments. In this work we use forest inventory data to parameterize a carbon cycle model to represent post-disturbance carbon trajectories of carbon pools and fluxes for specific forest types growing in high and low site productivity class settings. We then apply these trajectories to landscapes and regions based on forest age distributions derived from either the FIA data or from Landsat time series stacks (1985–2006) for 54 representative scenes throughout most of the conterminous United States.Weestimate the net carbon uptake in forests caused by post-disturbance growth and decomposition (“regrowth sink”) for forested regions across the country. At the landscape scale, the prevailing condition of positive net ecosystem productivity (NEP) is in stark contrast to local patcheswith large sources, particularly in the west where fires and clear cuts create contiguous disturbed patches. At the continental scale, regional differences in disturbance rates reflect management patterns of high disturbance rates in the Southeastern and South Central states, and lower disturbance rates in the Northeast andNorthern Lakes States. Despite low contemporary disturbance rates in the Northeast and Northern Lakes States (0.61 and 0.74% y−1), the regrowth sink there remains of moderate to large strength (88 and 57 g C m−2 y−1) owing to the continued legacy from historical clearing. Large regrowth sinks are also found in the Southeast, South Central, and Pacific Southwest regions (85, 86, and 95 g C m−2 y−1) where disturbance rates also tend to be higher (1.59, 1.38, and 0.93% y−1). Overall, the Landsat-derived disturbance rates are elevated relative to FIA-derived rates (1.19 versus 0.93% y−1) particularly for western regions. The differences only modestly adjust regional- and continental-scale carbon budgets, reducing NEP from forest regrowth by about 8%

Using an OBCD approach and Landsat TM data to detect harvesting on nonindustrial private property in Upper Michigan

Forest dynamics influence climate, biodiversity, and livelihoods at multiple scales, yet current resource policy addressing these dynamics is ineffective without reliable land use land cover change data. The collective impact of harvest decisions by many small forest owners can be substantial at the landscape scale, yet monitoring harvests and regrowth in these forests is challenging. Remote sensing is an obvious route to detect and monitor small-scale land use dynamics over large areas. Using an annual series of Landsat-5 Thematic Mapper (TM) images and a GIS shapefile of property boundaries, we identified units where harvests occurred from 2005 to 2011 using an Object-Based Change Detection (OBCD) approach. Percent of basal area harvested was verified using stand-level harvest data. Our method detected all harvests above 20% basal area removal in all forest types (northern hardwoods, mixed deciduous/coniferous, coniferous), on properties as small as 10 acres (0.4 ha; approximately four Landsat pixels). Our results had a resolution of about 10% basal area (that is, a selective harvest removal of 30% could be distinguished from one of 40%). Our method can be automated and used to measure annual harvest rates and intensities for large areas of the United States, providing critical information on land use transition

Improving the precision of dynamic forest parameter estimates using Landsat

The use of satellite-derived classification maps to improve post-stratified forest parameter estimates is well established.When reducing the variance of post-stratification estimates for forest change parameters such as forest growth, it is logical to use a change-related strata map. At the stand level, a time series of Landsat images is ideally suited for producing such a map. In this study, we generate strata maps based on trajectories of Landsat Thematic Mapper-based normalized difference vegetation index values, with a focus on post-disturbance recovery and recent measurements. These trajectories, from1985 to 2010, are converted to harmonic regression coefficient estimates and classified according to a hierarchical clustering algorithm from a training sample. The resulting strata maps are then used in conjunction with measured plots to estimate forest status and change parameters in an Alabama, USA study area. These estimates and the variance of the estimates are then used to calculate the estimated relative efficiencies of the post-stratified estimates. Estimated relative efficiencies around or above 1.2 were observed for total growth, total mortality, and total removals, with different strata maps being more effective for each. Possible avenues for improvement of the approach include the following: (1) enlarging the study area and (2) using the Landsat images closest to the time of measurement for each plot. Multitemporal satellite-derived strata maps show promise for improving the precision of change parameter estimates

Improving the precision of dynamic forest parameter estimates using Landsat

The use of satellite-derived classification maps to improve post-stratified forest parameter estimates is well established.When reducing the variance of post-stratification estimates for forest change parameters such as forest growth, it is logical to use a change-related strata map. At the stand level, a time series of Landsat images is ideally suited for producing such a map. In this study, we generate strata maps based on trajectories of Landsat Thematic Mapper-based normalized difference vegetation index values, with a focus on post-disturbance recovery and recent measurements. These trajectories, from1985 to 2010, are converted to harmonic regression coefficient estimates and classified according to a hierarchical clustering algorithm from a training sample. The resulting strata maps are then used in conjunction with measured plots to estimate forest status and change parameters in an Alabama, USA study area. These estimates and the variance of the estimates are then used to calculate the estimated relative efficiencies of the post-stratified estimates. Estimated relative efficiencies around or above 1.2 were observed for total growth, total mortality, and total removals, with different strata maps being more effective for each. Possible avenues for improvement of the approach include the following: (1) enlarging the study area and (2) using the Landsat images closest to the time of measurement for each plot. Multitemporal satellite-derived strata maps show promise for improving the precision of change parameter estimates

Utilization of Landsat Imagery to Assess the Impacts of Oil and Gas Extraction on the Tazovsky Peninsula, Siberia

Climatic warming of the Arctic is leading to landscape change through cascading biophysical feedbacks; development, such as oil and gas exploration and extraction, can accelerate or worsen these impacts. Due to restricted access to oil and natural gas fields, in situ environmental impact studies are only allowed in some regions. Satellite imagery analysis provides a mean for assessing impacts in areas with limited access. The Yamburg oil and gas field in western Siberia serves as a case study to assess the effects of infrastructure on an Arctic landscape. This project quantifies the land-cover disturbance that occurred during the development and expansion of the Yamburg field. Google’s recently developed, cloud-based image processing platform, Google Earth Engine, was used in conjunction with traditional Geographic Information System (GIS) analysis to detect, map, and quantify the impacts of infrastructure on the Tazovsky Peninsula between 1983 and 2016, utilizing imagery from the Landsat 4, 5, and 8 satellites. Landscape fragmentation metrics, the Normalized Difference Vegetation Index (NDVI), and change analysis quantified the impacts of extraction infrastructure on the surrounding landscape. As distance from the infrastructure and time since field establishment increased, the associated impacts decreased