Long-term information of phytoplankton bloom is critical for assessing the processes driving blooms in lakes. A three-decade survey of the phytoplankton blooms was completed for Erhai Lake from 1987 to 2016 with Landsat imagery. A modified three-band model using Landsat broad bands is developed by comparing reflectance data from Landsat imagery to two field datasets. The model is applied to the archived imagery (1987–2016) to predict chlorophyll-a (Chl-a). Predicted ln(Chl-a) and observed ln(Chl-a) measurements are significantly correlated (R² = 0.70; RMSE = 0.13 ug/L). Bloom maps are generated by identifying Landsat pixels that have Chl-a concentrations larger than 20 ug/L as bloom area. Bloom extent and magnitude are estimated. Our study reveals that algal blooms first occurred in 1996 with a bloom area of 150 km². Bloom occurred frequently from 2002 to 2016, with extreme blooms in 2003, 2013 and 2016. Algal blooms were mostly distributed in the northern and southern part of the lake. The proposed method uses one model for all Landsat images for Erhai Lake and can predict past blooms and extend the record to early years when field data is not available. The bloom extent and magnitude produced in this study can be used as the basis for the understanding of the processes that control the bloom outbreak.

Albedo is a key factor in modulating the absorption of solar radiation on ice surfaces. Satellite measurements have shown a general reduction in albedo across the Greenland ice sheet over the past few decades, particularly along the western margin of the ice sheet, a region known as the Dark Zone (albedo < 0.45). Here we chose a combination of Landsat 4–8 and Sentinel 2 imagery to enable us to derive the longest record of albedo variations in the Dark Zone, running from 1984 to 2020. We developed a simple, pragmatic and efficient sensor transformation to provide a long time series of consistent, harmonized satellite imagery. Narrow to broadband conversion algorithms were developed from regression models of harmonized satellite data and in situ albedo from the Program for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather stations. The albedo derived from the harmonized Landsat and Sentinel 2 data shows that the maximum extent of the Dark Zone expanded rapidly between 2005 and 2007, increasing to ~280% of the average annual maximum extent of 2900 km² to ~8000 km² since. The Dark Zone is continuing to darken slowly, with the average annual minimum albedo decreasing at a rate of $\sim \!-0.0006 \pm 0.0004 \, {\rm a}^{-1}$ (p = 0.16, 2001–2020).

Dark (low albedo) surface ice on the Greenland Ice Sheet enhances melting and subsequent runoff, a major mass loss contributor during the ablation season. The accumulation of both biological (e.g., glacier ice algae) and abiotic (e.g., mineral dust) light‐absorbing particulates are important darkening factors, that are potentially influenced by the duration of snow‐free, bare ice (a phenological factor), and other geo‐topographical factors such as elevation, slope, aspect and the distance from the ice margin. Here, we present the first medium‐resolution (30 m) analysis of the phenological and geo‐topographical controls on the distribution of dark ice in SE and SW Greenland from statistical analysis of data derived from a harmonized satellite albedo product and ArcticDEM. The duration of bare ice primarily controls the distribution of dark surface ice, allowing for algae growth on inland ice surfaces in particular, whereas geo‐topographical factors are only secondary controls.

Albedo plays a key role in regulating the absorption of solar radiation within ice surfaces and hence strongly regulates the production of meltwater. A combination of Landsat and Sentinel 2 data provides the longest continuous medium resolution (10–30 m) earth surface observatory records. An albedo product (harmonized satellite albedo, hereafter HSA) has already been developed and validated for the Greenland Ice Sheet (GrIS), using harmonized Landsat 4–8 and Sentinel 2 datasets. In this paper, the HSA was validated for various Arctic and alpine glaciers and ice caps using in situ measurements. We determine the optimal spatial window size in point-to-pixel analysis, the best practices in evaluating remote sensing algorithms with groundtruth data, and cross sensor comparison of the Landsat 9 (L9) and Landsat 8 (L8) data. The impact of the spatial window size on measured ice surface homogeneity and albedo validation was analysed at both local and regional scales. Homogeneity statistics calculated from the grey-level co-occurrence matrix (GLCM) suggest that the ice surface becomes more homogeneous as the image resolution becomes coarser. The optimal spatial window size was found to be 90 m, based on maximizing the statistical and graphical measures while minimizing the root mean square error and bias. HSAs generally agree closely with in situ albedo measurements (e.g. Pearson’s R ranges from 0.68 to 0.92) across various Arctic and alpine glaciers and ice caps. Cross sensor differences between L9 and L8 are minor, and we suggest that no harmonization is necessary to add L9 to our HSA product.