Detecting early warning signals of tree mortality in boreal North America using multiscale satellite data.

Increasing tree mortality from global change drivers such as drought and biotic infestations is a widespread phenomenon, including in the boreal zone where climate changes and feedbacks to the Earth system are relatively large. Despite the importance for science and management communities, our ability to forecast tree mortality at landscape to continental scales is limited. However, two independent information streams have the potential to inform and improve mortality forecasts: repeat forest inventories and satellite remote sensing. Time series of tree-level growth patterns indicate that productivity declines and related temporal dynamics often precede mortality years to decades before death. Plot-level productivity, in turn, has been related to satellite-based indices such as the Normalized difference vegetation index (NDVI). Here we link these two data sources to show that early warning signals of mortality are evident in several NDVI-based metrics up to 24 years before death. We focus on two repeat forest inventories and three NDVI products across western boreal North America where productivity and mortality dynamics are influenced by periodic drought. These data sources capture a range of forest conditions and spatial resolution to highlight the sensitivity and limitations of our approach. Overall, results indicate potential to use satellite NDVI for early warning signals of mortality. Relationships are broadly consistent across inventories, species, and spatial resolutions, although the utility of coarse-scale imagery in the heterogeneous aspen parkland was limited. Longer-term NDVI data and annually remeasured sites with high mortality levels generate the strongest signals, although we still found robust relationships at sites remeasured at a typical 5 year frequency. The approach and relationships developed here can be used as a basis for improving forest mortality models and monitoring systems.

Automatic sub-pixel co-registration of Landsat-8 OLI and Sentinel-2A MSI images using phase correlation and machine learning based mapping.

This study investigates misregistration issues between Landsat-8/OLI and Sentinel-2A/MSI at 30 m resolution, and between multi-temporal Sentinel-2A images at 10 m resolution using a phase correlation approach and multiple transformation functions. Co-registration of 45 Landsat-8 to Sentinel-2A pairs and 37 Sentinel-2A to Sentinel-2A pairs were analyzed. Phase correlation proved to be a robust approach that allowed us to identify hundreds and thousands of control points on images acquired more than 100 days apart. Overall, misregistration of up to 1.6 pixels at 30 m resolution between Landsat-8 and Sentinel-2A images, and 1.2 pixels and 2.8 pixels at 10 m resolution between multi-temporal Sentinel-2A images from the same and different orbits, respectively, were observed. The non-linear Random Forest regression used for constructing the mapping function showed best results in terms of root mean square error (RMSE), yielding an average RMSE error of 0.07±0.02 pixels at 30 m resolution, and 0.09±0.05 and 0.15±0.06 pixels at 10 m resolution for the same and adjacent Sentinel-2A orbits, respectively, for multiple tiles and multiple conditions. A simpler 1st order polynomial function (affine transformation) yielded RMSE of 0.08±0.02 pixels at 30 m resolution and 0.12±0.06 (same Sentinel-2A orbits) and 0.20±0.09 (adjacent orbits) pixels at 10 m resolution.

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product.

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Large carbon release legacy from bark beetle outbreaks across Western United States.

Warmer conditions over the past two decades have contributed to rapid expansion of bark beetle outbreaks killing millions of trees over a large fraction of western United States (US) forests. These outbreaks reduce plant productivity by killing trees and transfer carbon from live to dead pools where carbon is slowly emitted to the atmosphere via heterotrophic respiration which subsequently feeds back to climate change. Recent studies have begun to examine the local impacts of bark beetle outbreaks in individual stands, but the full regional carbon consequences remain undocumented for the western US. In this study, we quantify the regional carbon impacts of the bark beetle outbreaks taking place in western US forests. The work relies on a combination of postdisturbance forest regrowth trajectories derived from forest inventory data and a process-based carbon cycle model tracking decomposition, as well as aerial detection survey (ADS) data documenting the regional extent and severity of recent outbreaks. We find that biomass killed by bark beetle attacks across beetle-affected areas in western US forests from 2000 to 2009 ranges from 5 to 15 Tg C yr(-1) and caused a reduction of net ecosystem productivity (NEP) of about 6.1-9.3 Tg C y(-1) by 2009. Uncertainties result largely from a lack of detailed surveys of the extent and severity of outbreaks, calling out a need for improved characterization across western US forests. The carbon flux legacy of 2000-2009 outbreaks will continue decades into the future (e.g., 2040-2060) as committed emissions from heterotrophic respiration of beetle-killed biomass are balanced by forest regrowth and accumulation.

Land-cover change and avian diversity in the conterminous United States.

Changes in land use and land cover have affected and will continue to affect biological diversity worldwide. Yet, understanding the spatially extensive effects of land-cover change has been challenging because data that are consistent over space and time are lacking. We used the U.S. National Land Cover Dataset Land Cover Change Retrofit Product and North American Breeding Bird Survey data to examine land-cover change and its associations with diversity of birds with principally terrestrial life cycles (landbirds) in the conterminous United States. We used mixed-effects models and model selection to rank associations by ecoregion. Land cover in 3.22% of the area considered in our analyses changed from 1992 to 2001, and changes in species richness and abundance of birds were strongly associated with land-cover changes. Changes in species richness and abundance were primarily associated with changes in nondominant types of land cover, yet in many ecoregions different types of land cover were associated with species richness than were associated with abundance. Conversion of natural land cover to anthropogenic land cover was more strongly associated with changes in bird species richness and abundance than persistence of natural land cover in nearly all ecoregions and different covariates were most strongly associated with species richness than with abundance in 11 of 17 ecoregions. Loss of grassland and shrubland affected bird species richness and abundance in forested ecoregions. Loss of wetland was associated with bird abundance in forested ecoregions. Our findings highlight the value of understanding changes in nondominant land cover types and their association with bird diversity in the United States.

Conservation of forest birds: evidence of a shifting baseline in community structure.

BACKGROUND: Quantifying changes in forest bird diversity is an essential task for developing effective conservation actions. When subtle changes in diversity accumulate over time, annual comparisons may offer an incomplete perspective of changes in diversity. In this case, progressive change, the comparison of changes in diversity from a baseline condition, may offer greater insight because changes in diversity are assessed over longer periods of times. Our objectives were to determine how forest bird diversity has changed over time and whether those changes were associated with forest disturbance. METHODOLOGY/PRINCIPAL FINDINGS: We used North American Breeding Bird Survey data, a time series of Landsat images classified with respect to land cover change, and mixed-effects models to associate changes in forest bird community structure with forest disturbance, latitude, and longitude in the conterminous United States for the years 1985 to 2006. We document a significant divergence from the baseline structure for all birds of similar migratory habit and nest location, and all forest birds as a group from 1985 to 2006. Unexpectedly, decreases in progressive similarity resulted from small changes in richness (<1 species per route for the 22-year study period) and modest losses in abundance (-28.7 - -10.2 individuals per route) that varied by migratory habit and nest location. Forest disturbance increased progressive similarity for Neotropical migrants, permanent residents, ground nesting, and cavity nesting species. We also documented highest progressive similarity in the eastern United States. CONCLUSIONS/SIGNIFICANCE: Contemporary forest bird community structure is changing rapidly over a relatively short period of time (e.g., approximately 22 years). Forest disturbance and forest regeneration are primary factors associated with contemporary forest bird community structure, longitude and latitude are secondary factors, and forest loss is a tertiary factor. Importantly, these findings suggest some regions of the United States may already fall below the habitat amount threshold where fragmentation effects become important predictors of forest bird community structure.