Abstract. Surface meltwater is widespread around the margin of the Antarctic Ice Sheet and has the potential to influence ice-shelf stability, ice-dynamic processes and ice-albedo feedbacks. Whilst the general spatial distribution of surface meltwater across the Antarctic continent is now relatively well known, our understanding of the seasonal and multi-year evolution of surface meltwater is limited. Attempts to generate robust time series of melt cover have largely been constrained by computational expense or limited ice surface visibility associated with mapping from optical satellite imagery. Here, we implement an existing meltwater detection method alongside a novel method for calculating visibility metrics within Google Earth Engine. This enables us to quantify uncertainty induced by cloud cover and variable image data coverage, allowing us to automatically generate time series of surface melt area over large spatial and temporal scales. We demonstrate our method on the Amery Ice Shelf region of East Antarctica, analysing 4,164 Landsat 7 and 8 optical images between 2005 and 2020. Results show high interannual variability in surface meltwater cover, with mapped cumulative lake area totals ranging from 384 km2 to 3898 km2 per melt season. However, by incorporating image visibility assessments into our results, we estimate that cumulative total lake areas are on average 42 % higher than minimum mapped values, highlighting the importance of accounting for variations in image visibility when mapping lake areas. In a typical melt season, total lake area remains low throughout November and early December, before increasing, on average, by an order of magnitude during the second half of December. Peak lake area most commonly occurs during January, before decreasing during February as lakes freeze over. We show that modelled melt predictions from a regional climate model provides a good indication of lake cover in the Amery region, and that annual lake coverage is strongly associated with phases of the Southern Annular Mode (SAM); surface melt is typically highest in years with a negative austral summer SAM index. Furthermore, we suggest that melt-albedo feedbacks modulate the spatial distribution of meltwater in the region, with the exposure of blue ice from persistent katabatic wind scouring influencing the susceptibility of melt ponding. Results demonstrate how our method could be scaled up to generate a multi-year time series record of surface water extent from optical imagery at a continent-wide scale.