Monitoring and classifying forest disturbance using Landsat time series has improved greatly over the past decade, with many new algorithms taking advantage of the high-quality, cost free data in the ...archive. Much of the innovation has been focused on use of sophisticated workflows that consist of a logical sequence of processes and rules, multiple statistical functions, and parameter sets that must be calibrated to accurately classify disturbance. For many algorithms, calibration has been local to areas of interest and the algorithm's classification performance has been good under those circumstances. When applied elsewhere, however, algorithm performance has suffered. An alternative strategy for calibration may be to use the locally tested parameter values in conjunction with a statistical approach (e.g., Random Forests; RF) to align algorithm classification with a reference disturbance dataset, a process we call secondary classification. We tested that strategy here using RF with LandTrendr, an algorithm that runs on one spectral band or index. Disturbance detection using secondary classification was spectral band- or index-dependent, with each spectral dimension providing some unique detections and different error rates. Using secondary classification, we tested whether an integrated multispectral LandTrendr ensemble, with various combinations of the six basic Landsat reflectance bands and seven common spectral indices, improves algorithm performance. Results indicated a substantial reduction in errors relative to secondary classification based on single bands/indices, revealing the importance of a multispectral approach to forest disturbance detection. To explain the importance of specific bands and spectral indices in the multispectral ensemble, we developed a disturbance signal-to-noise metric that clearly highlighted the value of shortwave-infrared reflectance, especially when paired with near-infrared reflectance.
•A Landsat disturbance signal-to-noise ratio (DSNR) metric was developed.•DSNR shown to be strongly correlated to forest disturbance detection error•Secondary classification used to improve forest disturbance detection accuracy.•SWIR bands and SWIR-based indices most important•NIR band important when used in combination with SWIR
The 2001 Forest Service Roadless Rule prohibits roadbuilding in forests across an area equivalent to the combined states of New York and Maine (236 000 km2). There have been recent assertions that ...roads are needed to prevent fire and to keep forests healthy. Despite twenty years of ongoing forest health monitoring and the unique scope and ecological significance of this network of roadless areas, there has to date been no integrated assessment of the relationship between roads and forest health. Here, this question was addressed by synthesizing different sources of nationally consistent, longitudinal monitoring data. Agency management records show that a lack of roads has not stopped fire prevention measures; fuel management activities in roadless areas have actually been more numerous on a per-square kilometer basis than elsewhere in the National Forest System, although activities in areas with roads cover larger areas. Historical fire maps indicate that forests with and without roads have burned at similar rates since the Rule took effect. The apparent neutrality of roads with respect to fire occurrence may be due to higher rates of human caused ignition near roads offsetting advantages related to more agile positioning of fire-fighting assets. Beyond the fire dimension of forest health, analysis of over 15 000 inventory plots showed that while tree root disease is only weakly correlated with proximity to roads, roads are strongly associated with the spread of invasive plant species in national forests. Non-native plants are twice as common within 152 meters (500 feet) of a road as farther away. Speculation that eliminating road prohibitions would improve forest health is not supported by nearly twenty years of monitoring data.
While Landsat has proved to be effective for monitoring many elements of forest condition and change, the platform has well-documented limitations in measuring forest structure, the vertical ...distribution of the canopy. This is important because structure determines several key ecosystem functions, including: carbon storage; habitat suitability; and timber volume. Canopy structure is directly measured by LiDAR, and it should be possible to train Landsat structure models at a highly local scale with the dense, global sample of full waveform LiDAR observations collected by NASA’s Global Ecosystem Dynamics Investigation (GEDI). Local models are expected to perform better because: (a) such models may take advantage of localized correlations between structure and canopy surface reflectance; and (b) to the extent that models revert to the mean of the calibration data due to a lack of discrimination, local models will revert to a more representative mean. We tested Landsat-based relative height predictions using a new GEDI asset on Google Earth Engine, described here. Mean prediction error declined by 23% and important prediction biases at the extremes of the range of canopy height dropped as model calibration became more local, minimizing forest structure signal saturation commonly associated with Landsat and other passive optical sensors. Our results suggest that Landsat-based maps of structural variables such as height and biomass may substantially benefit from the kind of local calibration that GEDI’s dense sample of LiDAR data supports.
The United States (U.S.) federal government provides imagery obtained by federally funded Earth Observation satellites typically at no cost. For many years Landsat was an exception to this trend, ...until 2008 when the United States Geological Survey (USGS) made Landsat data accessible via the internet for free. Substantial increases in downloads of Landsat imagery ensued and led to a rapid expansion of science and operational applications, serving government, private sector, and civil society. The Landsat program hence provides an example to space agencies worldwide on the value of open access for Earth Observation data and has spurred the adaption of similar policies globally, including the European Copernicus Program. Here, we describe important aspects of the Landsat free and open data policy and highlight the importance and continued relevance of this policy.
•Free and open data policy is key to the ongoing success of the Landsat program.•Described important aspects of the Landsat free and open data policy•Highlighted the importance and continued relevance of Landsat open data policy
The U.S. Geological Survey Land Change Monitoring, Assessment and Projection (USGS LCMAP) initiative is working toward a comprehensive capability to characterize land cover and land cover change ...using dense Landsat time series data. A suite of products including annual land cover maps and annual land cover change maps will be produced using the Landsat 4-8 data record. LCMAP products will initially be created for the conterminous United States (CONUS) and then extended to include Alaska and Hawaii. A critical component of LCMAP is the collection of reference data using the TimeSync tool, a web-based interface for manually interpreting and recording land cover from Landsat data supplemented with fine resolution imagery and other ancillary data. These reference data will be used for area estimation and validation of the LCMAP annual land cover products. Nearly 12,000 LCMAP reference sample pixels have been interpreted and a simple random subsample of these pixels has been interpreted independently by a second analyst (hereafter referred to as “duplicate interpretations”). The annual land cover reference class labels for the 1984–2016 monitoring period obtained from these duplicate interpretations are used to address the following questions: 1) How consistent are the reference class labels among interpreters overall and per class? 2) Does consistency vary by geographic region? 3) Does consistency vary as interpreters gain experience over time? 4) Does interpreter consistency change with improving availability and quality of imagery from 1984 to 2016? Overall agreement between interpreters was 88%. Class-specific agreement ranged from 46% for Disturbed to 94% for Water, with more prevalent classes (Tree Cover, Grass/Shrub and Cropland) generally having greater agreement than rare classes (Developed, Barren and Wetland). Agreement between interpreters remained approximately the same over the 12-month period during which these interpretations were completed. Increasing availability of Landsat and Google Earth fine resolution data over the 1984 to 2016 monitoring period coincided with increased interpreter consistency for the post-2000 data record. The reference data interpretation and quality assurance protocols implemented for LCMAP demonstrate the technical and practical feasibility of using the Landsat archive and intensive human interpretation to produce national, annual reference land cover data over a 30-year period. Protocols to estimate and enhance interpreter consistency are critical elements to document and ensure the quality of these reference data.
•A subset of pixels with duplicate interpretations quantifies consistency.•Interpreter agreement was 88% overall ranging from 46% (Disturbed) to 94% (Water).•Regional variation in class-specific agreement was observed.•Agreement stayed about the same as interpretations were finished over time.•Agreement was greater for data after 2000 coinciding with increased data density.
Accurate estimation of aboveground forest biomass stocks is required to assess the impacts of land use changes such as deforestation and subsequent regrowth on concentrations of atmospheric CO2. The ...Global Ecosystem Dynamics Investigation (GEDI) is a lidar mission launched by NASA to the International Space Station in 2018. GEDI was specifically designed to retrieve vegetation structure within a novel, theoretical sampling design that explicitly quantifies biomass and its uncertainty across a variety of spatial scales. In this paper we provide the estimates of pan-tropical and temperate biomass derived from two years of GEDI observations. We present estimates of mean biomass densities at 1 km resolution, as well as estimates aggregated to the national level for every country GEDI observes, and at the sub-national level for the United States. For all estimates we provide the standard error of the mean biomass. These data serve as a baseline for current biomass stocks and their future changes, and the mission’s integrated use of formal statistical inference points the way towards the possibility of a new generation of powerful monitoring tools from space.
Recent developments in remote sensing (RS) technology have made several sources of auxiliary data available to support forest inventories. Thus, a pertinent question is how different sources of RS ...data should be combined with field data to make inventories cost-efficient. Hierarchical model-based estimation has been proposed as a promising way of combining: (i) wall-to-wall optical data that are only weakly correlated with forest structure; (ii) a discontinuous sample of active RS data that are more strongly correlated with structure; and (iii) a sparse sample of field data. Model predictions based on the strongly correlated RS data source are used for estimating a model linking the target quantity with weakly correlated wall-to-wall RS data. Basing the inference on the latter model, uncertainties due to both modeling steps must be accounted for to obtain reliable variance estimates of estimated population parameters, such as totals or means. Here, we generalize previously existing estimators for hierarchical model-based estimation to cases with non-homogeneous error variance and cases with correlated errors, for example due to clustered sample data. This is an important generalization to take into account data from practical surveys. We apply the new estimation framework to case studies that mimic the data that will be available from the Global Ecosystem Dynamics Investigation (GEDI) mission and compare the proposed estimation framework with alternative methods. Aboveground biomass was the variable of interest, Landsat data were available wall-to-wall, and sample RS data were obtained from an airborne LiDAR campaign that produced simulated GEDI waveforms. The results show that generalized hierarchical model-based estimation has potential to yield more precise estimates than approaches utilizing only one source of RS data, such as conventional model-based and hybrid inferential approaches.
Disturbance is a critical ecological process in forested systems, and disturbance maps are important for understanding forest dynamics. Landsat data are a key remote sensing dataset for monitoring ...forest disturbance and there recently has been major growth in the development of disturbance mapping algorithms. Many of these algorithms take advantage of the high temporal data volume to mine subtle signals in Landsat time series, but as those signals become subtler, they are more likely to be mixed with noise in Landsat data. This study examines the similarity among seven different algorithms in their ability to map the full range of magnitudes of forest disturbance over six different Landsat scenes distributed across the conterminous US. The maps agreed very well in terms of the amount of undisturbed forest over time; however, for the ~30% of forest mapped as disturbed in a given year by at least one algorithm, there was little agreement about which pixels were affected. Algorithms that targeted higher-magnitude disturbances exhibited higher omission errors but lower commission errors than those targeting a broader range of disturbance magnitudes. These results suggest that a user of any given forest disturbance map should understand the map’s strengths and weaknesses (in terms of omission and commission error rates), with respect to the disturbance targets of interest.
When characterizing the processes that shape ecosystems, ecologists increasingly use the unique perspective offered by repeat observations of remotely sensed imagery. However, the concept of change ...embodied in much of the traditional remote-sensing literature was primarily limited to capturing large or extreme changes occurring in natural systems, omitting many more subtle processes of interest to ecologists. Recent technical advances have led to a fundamental shift toward an ecological view of change. Although this conceptual shift began with coarser-scale global imagery, it has now reached users of Landsat imagery, since these datasets have temporal and spatial characteristics appropriate to many ecological questions. We argue that this ecologically relevant perspective of change allows the novel characterization of important dynamic processes, including disturbances, long-term trends, cyclical functions, and feedbacks, and that these improvements are already facilitating our understanding of critical driving forces, such as climate change, ecological interactions, and economic pressures.
There are several new and imminent space-based sensors intended to support mapping of forest structure and biomass. These instruments, along with advancing cloud-based mapping platforms, will soon ...contribute to a proliferation of biomass maps. One means of differentiating the quality of different maps and estimation strategies will be comparison of results against independent field-based estimates at various scales. The Forest Inventory and Analysis Program of the US Forest Service (FIA) maintains a designed sample of uniformly measured field plots across the conterminous United States. This paper reports production of a map of statistical estimates of mean biomass, created at approximately the finest scale (64,000-hectare hexagons) allowed by FIA’s sample density. This map may be useful for assessing the accuracy of future remotely sensed biomass estimates. Equally important, fine-scale mapping of FIA estimates highlights several ways in which field- and remote sensing-based methods must be aligned to ensure comparability. For example, the biomass in standing dead trees, which may or may not be included in biomass estimates, represents a source of potential discrepancy that FIA shows to be particularly important in the Western US. Likewise, alternative allometric equations (which link measurable tree dimensions such as diameter to difficult-to-measure variables like biomass) strongly impact biomass estimates in ways that can vary over short distances. Potential mismatch in the conditions counted as forests also varies greatly over space. Field-to-map comparisons will ideally minimize these sources of uncertainty by adopting common allometry, carbon pools, and forest definitions. Our national hexagon-level benchmark estimates, provided in Supplementary Files, therefore addresses multiple pools and allometric approaches independently, while providing explicit forest area and uncertainty information. This range of information is intended to allow scientists to minimize potential discrepancies in support of unambiguous validation.