Satellite-Derived Summer Albedo Variations on the Greenland Ice Sheet from 1979 to 2024 Linked with Climatic Indices

CLARA-A3 currently provides the longest temporal coverage among available albedo products, with improvements in both retrieval algorithms and product coverage compared to earlier versions. This study first evaluates the performance of the CLARA-A3-SAL product over Greenland Ice Sheet (GrIS) and subsequently applies it to investigate spatiotemporal trends in summer albedo from 1979 to 2024. Validation against 32 in situ observation sites indicates negligible bias in the interior regions, with RMSE values ranging from 0.01 to 0.07. Although larger errors exist in the coastal ablation zone due to unresolved sub-grid surface heterogeneity, the product successfully captures observed spatiotemporal variability and long-term trends, demonstrating that CLARA-A3-SAL provides a generally reliable representation of surface albedo. Since 1979, the summer surface albedo averaged over the entire ice sheet has decreased at a rate of −0.24% decade−1. Albedo in the dry snow area has remained relatively stable and showed no significant correlation with most climate variables, except for the North Atlantic Oscillation (NAO) and the Greenland Blocking Index (GBI). Conversely, the marginal zone has undergone substantial darkening (−0.66% decade−1), which is strongly correlated with temperature, snowfall and melt, with meltwater showing the highest correlation (r = −0.90, p < 0.01). This suggests that meltwater-driven grain growth and exposure of bare ice are the primary drivers of albedo reduction over the non-dry snow zone. Large-scale atmospheric circulation also plays a key role: the GBI exhibits the strongest association with albedo (r = −0.63, p < 0.05), underscoring the importance of persistent blocking in amplifying surface warming and darkening. Furthermore, decadal-scale variability associated with the Atlantic Multidecadal Oscillation (AMO) and the Pacific Decadal Oscillation (PDO) modulates both the magnitude and spatial pattern of albedo changes across GrIS, with AMO+ generally linked to reduced albedo and PDO+ tending to enhance it.