The LandTrendr (LT) algorithm has been used widely for analysis of change in Landsat spectral time series data, but requires significant pre-processing, data management, and computational resources, and is only accessible to the community in a proprietary programming language (IDL). Here, we introduce LT for the Google Earth Engine (GEE) platform. The GEE platform simplifies pre-processing steps, allowing focus on the translation of the core temporal segmentation algorithm. Temporal segmentation involved a series of repeated random access calls to each pixel’s time series, resulting in a set of breakpoints (“vertices”) that bound straight-line segments. The translation of the algorithm into GEE included both transliteration and code analysis, resulting in improvement and logic error fixes. At six study areas representing diverse land cover types across the U.S., we conducted a direct comparison of the new LT-GEE code against the heritage code (LT-IDL). The algorithms agreed in most cases, and where disagreements occurred, they were largely attributable to logic error fixes in the code translation process. The practical impact of these changes is minimal, as shown by an example of forest disturbance mapping. We conclude that the LT-GEE algorithm represents a faithful translation of the LT code into a platform easily accessible by the broader user community.

There is a pressing need to map changes in forest structure from the earliest time period possible given forest management policies and accelerated disturbances from climate change. The availability of Landsat data from over four decades helps researchers study an ecologically meaningful length of time. Forest structure is most often mapped utilizing lidar data, however these data are prohibitively expensive and cover a narrow temporal window relative to the Landsat archive. Here we describe a technique to use the entire length of the Landsat archive from Multispectral Scanner to Operational Land Imager (M2O) to produce three novel outcomes: (1) we used the M2O dataset and standard change vector analysis methods to classify annual forest structure in northwestern Montana from 1972 to 2015, (2) we improved the accuracy of each yearly forest structure classification by applying temporal continuity rules to the whole time series, with final accuracies ranging from 97% to 68% respectively for two and six-category classifications, and (3) we demonstrated the importance of pre-1984 Landsat data for long-term change studies. As the Landsat program continues to acquire Earth imagery into the foreseeable future, time series analyses that aid in classifying forest structure accurately will be key to the success of any land management changes in the future.