skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Landsat derived patterns of lake drainage in northern Alaska between 1975-2019
Forty-five years (i.e. 1975-2019) of Landsat observations were used to map the spatiotemporal patterns of lake drainage in northern Alaska. All Landsat data was pre-processed by the United States Geological Survey and downloaded by google earth engine in a radiometrically, atmospherically, and geometrically terrain-corrected state. We used Landsat surface reflectance products acquired from the Multispectral Scanner (MSS), Terrestrial Mapper (TM), Enhanced Terrestrial Mapper Plus (ETM+), and Operational Land Imager (OLI) sensors to compute eight image mosaics for the ice-free period (June 15 to September 1) at five-year time-periods. This data product represents the change in lake area between time-periods or epochs. All data used to spatially identify patterns of lake change are presented in a map (LakeDrainChg.tif), where the associated morphometric controls on drainage are summarized for each of the ~33,000 lake boundaries (Lake_Drainage.shp). All data can also be viewed within Google Earth Engine (https://code.earthengine.google.com/?accept_repo=users/mjlara71/LakeDrainage; accessed on 13 September 2021) or clone Git repository (git clone https://earthengine.googlesource.com/users/mjlara71/LakeDrainage ; accessed on 13 September 2021).  more » « less
Award ID(s):
1928048
PAR ID:
10311727
Author(s) / Creator(s):
; ;
Publisher / Repository:
NSF Arctic Data Center
Date Published:
Subject(s) / Keyword(s):
Lake drainage gradual lake drainage catastrophic lake drainage talik permafrost gully tundra Alaska
Format(s):
Medium: X Other: text/xml
Sponsoring Org:
National Science Foundation
More Like this
  1. Nearly 25% of all lakes on earth are located at high latitudes. These lakes are formed by a combination of thermokarst, glacial, and geological processes. Evidence suggests that the origin of periglacial lake formation may be an important factor controlling the likelihood of lakes to drain. However, geospatial data regarding the spatial distribution of these dominant Arctic and subarctic lakes are limited or do not exist. Here, we use lake-specific morphological properties using the Arctic Digital Elevation Model (DEM) and Landsat imagery to develop a Thermokarst lake Settlement Index (TSI), which was used in combination with available geospatial datasets of glacier history and yedoma permafrost extent to classify Arctic and subarctic lakes into Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes, respectively. This lake origin dataset was used to evaluate the influence of lake origin on drainage between 1985 and 2019 in northern Alaska. The lake origin map and lake drainage datasets were synthesized using five-year seamless Landsat ETM+ and OLI image composites. Nearly 35,000 lakes and their properties were characterized from Landsat mosaics using an object-based image analysis. Results indicate that the pattern of lake drainage varied by lake origin, and the proportion of lakes that completely drained (i.e., >60% area loss) between 1985 and 2019 in Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes were 12.1, 9.5, 8.7, and 0.0%, respectively. The lakes most vulnerable to draining were small thermokarst (non-yedoma) lakes (12.7%) and large yedoma lakes (12.5%), while the most resilient were large and medium-sized glacial lakes (4.9 and 4.1%) and Maar lakes (0.0%). This analysis provides a simple remote sensing approach to estimate the spatial distribution of dominant lake origins across variable physiography and surficial geology, useful for discriminating between vulnerable versus resilient Arctic and subarctic lakes that are likely to change in warmer and wetter climates. 
    more » « less
  2. We developed a fire perimeter dataset for eastern Siberian taiga and tundra zones from 2001-2020 based on Landsat imagery. Our study area spanned 7.8 million square kilometers across eight ecozones of eastern Siberian taiga and tundra from approximately 61-72.5°North (N) and 100°East (E)-176°West (W). We used the cloud computing power of Google Earth Engine to access the Landsat archive. We generated composite images for the annual growing season (May - September), which allowed us to mitigate missing data from snow-cover, cloud-cover, and the Landsat 7 scan line error. We used annual composites to calculate the difference Normalized Burn Ratio (dNBR) for each year. Finally, we converted the annual dNBR images to binary burned or unburned imagery that was used to vectorize fire perimeters. We map 22,110 fires burning 150.5 million hectares over 20 years. 
    more » « less
  3. The Landsat satellites provide decades of near‐global surface reflectance measurements that are increasingly used to assess interannual changes in terrestrial ecosystem function. These assessments often rely on spectral indices related to vegetation greenness and productivity (e.g. Normalized Difference Vegetation Index, NDVI). Nevertheless, multiple factors impede multi‐decadal assessments of spectral indices using Landsat satellite data, including ease of data access and cleaning, as well as lingering issues with cross‐sensor calibration and challenges with irregular timing of cloud‐free acquisitions. To help address these problems, we developed the ‘LandsatTS' package for R. This software package facilitates sample‐based time series analysis of surface reflectance and spectral indices derived from Landsat sensors. The package includes functions that enable the extraction of the full Landsat 5, 7, and 8 records from Collection 2 for point sample locations or small study regions using Google Earth Engine accessed directly from R. Moreover, the package includes functions for 1) rigorous data cleaning, 2) cross‐sensor calibration, 3) phenological modeling, and 4) time series analysis. For an example application, we show how ‘LandsatTS' can be used to assess changes in annual maximum vegetation greenness from 2000 to 2022 across the Noatak National Preserve in northern Alaska, USA. Overall, this software provides a suite of functions to enable broader use of Landsat satellite data for assessing and monitoring terrestrial ecosystem function during recent decades across local to global geographic extents. 
    more » « less
  4. Abstract Lakes represent as much as ∼25% of the total land surface area in lowland permafrost regions. Though decreasing lake area has become a widespread phenomenon in permafrost regions, our ability to forecast future patterns of lake drainage spanning gradients of space and time remain limited. Here, we modeled the drivers of gradual (steady declining lake area) and catastrophic (temporally abrupt decrease in lake area) lake drainage using 45 years of Landsat observations (i.e. 1975–2019) across 32 690 lakes spanning climate and environmental gradients across northern Alaska. We mapped lake area using supervised support vector machine classifiers and object based image analyses using five-year Landsat image composites spanning 388 968 km2. Drivers of lake drainage were determined with boosted regression tree models, using both static (e.g. lake morphology, proximity to drainage gradient) and dynamic predictor variables (e.g. temperature, precipitation, wildfire). Over the past 45 years, gradual drainage decreased lake area between 10% and 16%, but rates varied over time as the 1990s recorded the highest rates of gradual lake area losses associated with warm periods. Interestingly, the number of catastrophically drained lakes progressively decreased at a rate of ∼37% decade−1from 1975–1979 (102–273 lakes draining year−1) to 2010–2014 (3–8 lakes draining year−1). However this 40 year negative trend was reversed during the most recent time-period (2015–2019), with observations of catastrophic drainage among the highest on record (i.e. 100–250 lakes draining year−1), the majority of which occurred in northwestern Alaska. Gradual drainage processes were driven by lake morphology, summer air and lake temperature, snow cover, active layer depth, and the thermokarst lake settlement index (R2adj= 0.42, CV = 0.35,p< 0.0001), whereas, catastrophic drainage was driven by the thawing season length, total precipitation, permafrost thickness, and lake temperature (R2adj= 0.75, CV = 0.67,p< 0.0001). Models forecast a continued decline in lake area across northern Alaska by 15%–21% by 2050. However these estimates are conservative, as the anticipated amplitude of future climate change were well-beyond historical variability and thus insufficient to forecast abrupt ‘catastrophic’ drainage processes. Results highlight the urgency to understand the potential ecological responses and feedbacks linked with ongoing Arctic landscape reorganization. 
    more » « less
  5. Landsat 5 has produced imagery for decades that can now be viewed and manipulated in Google Earth Engine, but a general, automated way of producing a coherent time series from these images—particularly over cloudy areas in the distant past—is elusive. Here, we create a land use and land cover (LULC) time series for part of tropical Mato Grosso, Brazil, using the Bayesian Updating of Land Cover: Unsupervised (BULC-U) technique. The algorithm built backward in time from the GlobCover 2009 data set, a multi-category global LULC data set at 300 m resolution for the year 2009, combining it with Landsat time series imagery to create a land cover time series for the period 1986–2000. Despite the substantial LULC differences between the 1990s and 2009 in this area, much of the landscape remained the same: we asked whether we could harness those similarities and differences to recreate an accurate version of the earlier LULC. The GlobCover basis and the Landsat-5 images shared neither a common spatial resolution nor time frame, But BULC-U successfully combined the labels from the coarser classification with the spatial detail of Landsat. The result was an accurate fine-scale time series that quantified the expansion of deforestation in the study area, which more than doubled in size during this time. Earth Engine directly enabled the fusion of these different data sets held in its catalog: its flexible treatment of spatial resolution, rapid prototyping, and overall processing speed permitted the development and testing of this study. Many would-be users of remote sensing data are currently limited by the need to have highly specialized knowledge to create classifications of older data. The approach shown here presents fewer obstacles to participation and allows a wide audience to create their own time series of past decades. By leveraging both the varied data catalog and the processing speed of Earth Engine, this research can contribute to the rapid advances underway in multi-temporal image classification techniques. Given Earth Engine’s power and deep catalog, this research further opens up remote sensing to a rapidly growing community of researchers and managers who need to understand the long-term dynamics of terrestrial systems. 
    more » « less