Global variability and controls on the accumulation of fallout radionuclides in cryoconite.

The accumulation of fallout radionuclides (FRNs) from nuclear weapons testing and nuclear accidents has been evaluated for over half a century in natural environments; however, until recently their distribution and abundance within glaciers have been poorly understood. Following a series of individual studies of FRNs, specifically 137Cs, 241Am and 210Pb, deposited on the surface of glaciers, we now understand that cryoconite, a material commonly found in the supraglacial environment, is a highly efficient accumulator of FRNs, both artificial and natural. However, the variability of FRN activity concentrations in cryoconite across the global cryosphere has never been assessed. This study thus aims to both synthesize current knowledge on FRNs in cryoconite and assess the controls on variability of activity concentrations. We present a global database of new and previously published data based on gamma spectrometry of cryoconite and proglacial sediments, and assess the extent to which a suite of environmental and physical factors can explain spatial variability in FRN activity concentrations in cryoconite. We show that FRNs are not only found in cryoconite on glaciers within close proximity to specific sources of radioactivity, but across the global cryosphere, and at activity concentrations up to three orders of magnitude higher than those found in soils and sediments in the surrounding environment. We also show that the organic content of cryoconite exerts a strong control on accumulation of FRNs, and that activity concentrations in cryoconite are some of the highest ever described in environmental matrices outside of nuclear exclusion zones, occasionally in excess of 10,000 Bq kg-1. These findings highlight a need for significant improvements in the understanding of the fate of legacy contaminants within glaciated catchments. Future interdisciplinary research is required on the mechanisms governing their accumulation, storage, and mobility, and their potential to create time-dependent impacts on downstream water quality and ecosystem sustainability.

Satellite-derived multivariate world-wide lake physical variable timeseries for climate studies.

A consistent dataset of lake surface water temperature, ice cover, water-leaving reflectance, water level and extent is presented. The collection constitutes the Lakes Essential Climate Variable (ECV) for inland waters. The data span combined satellite observations from 1992 to 2020 inclusive and quantifies over 2000 relatively large lakes, which represent a small fraction of the number of lakes worldwide but a significant fraction of global freshwater surface. Visible and near-infrared optical imagery, thermal imagery and microwave radar data from satellites have been exploited. All observations are provided in a common grid at 1/120° latitude-longitude resolution, jointly in daily files. The data/algorithms have been validated against in situ measurements where possible. Consistency analysis between the variables has guided the development of the joint dataset. It is the most complete collection of consistent satellite observations of the Lakes ECV currently available. Lakes are of significant interest to scientific disciplines such as hydrology, limnology, climatology, biogeochemistry and geodesy. They are a vital resource for freshwater supply, and key sentinels for global environmental change.

Complementary water quality observations from high and medium resolution Sentinel sensors by aligning chlorophyll-a and turbidity algorithms.

High resolution imaging spectrometers are prerequisite to address significant data gaps in inland optical water quality monitoring. In this work, we provide a data-driven alignment of chlorophyll-a and turbidity derived from the Sentinel-2 MultiSpectral Imager (MSI) with corresponding Sentinel-3 Ocean and Land Colour Instrument (OLCI) products. For chlorophyll-a retrieval, empirical 'ocean colour' blue-green band ratios and a near infra-red (NIR) band ratio algorithm, as well as a semi-analytical three-band NIR-red ratio algorithm, were included in the analysis. Six million co-registrations with MSI and OLCI spanning 24 lakes across five continents were analysed. Following atmospheric correction with POLYMER, the reflectance distributions of the red and NIR bands showed close similarity between the two sensors, whereas the distribution for blue and green bands was positively skewed in the MSI results compared to OLCI. Whilst it is not possible from this analysis to determine the accuracy of reflectance retrieved with either MSI or OLCI results, optimizing water quality algorithms for MSI against those previously derived for the Envisat Medium Resolution Imaging Spectrometer (MERIS) and its follow-on OLCI, supports the wider use of MSI for aquatic applications. Chlorophyll-a algorithms were thus tuned for MSI against concurrent OLCI observations, resulting in significant improvements against the original algorithm coefficients. The mean absolute difference (MAD) for the blue-green band ratio algorithm decreased from 1.95 mg m-3 to 1.11 mg m-3, whilst the correlation coefficient increased from 0.61 to 0.80. For the NIR-red band ratio algorithms improvements were modest, with the MAD decreasing from 4.68 to 4.64 mg m-3 for the empirical red band ratio algorithm, and 3.73 to 3.67 for the semi-analytical 3-band algorithm. Three implementations of the turbidity algorithm showed improvement after tuning with the resulting distributions having reduced bias. The MAD reduced from 0.85 to 0.72, 1.22 to 1.10 and 1.93 to 1.55 FNU for the 665, 708 and 778 nm implementations respectively. However, several sources of uncertainty remain: adjacent land showed high divergence between the sensors, suggesting that high product uncertainty near land continues to be an issue for small water bodies, while it cannot be stated at this point whether MSI or OLCI results are differentially affected. The effect of spectrally wider bands of the MSI on algorithm sensitivity to chlorophyll-a and turbidity cannot be fully established without further availability of in situ optical measurements.

Dynamics of glacier calving at the ungrounded margin of Helheim Glacier, southeast Greenland.

During summer 2013 we installed a network of 19 GPS nodes at the ungrounded margin of Helheim Glacier in southeast Greenland together with three cameras to study iceberg calving mechanisms. The network collected data at rates up to every 7 s and was designed to be robust to loss of nodes as the glacier calved. Data collection covered 55 days, and many nodes survived in locations right at the glacier front to the time of iceberg calving. The observations included a number of significant calving events, and as a consequence the glacier retreated ~1.5 km. The data provide real-time, high-frequency observations in unprecedented proximity to the calving front. The glacier calved by a process of buoyancy-force-induced crevassing in which the ice downglacier of flexion zones rotates upward because it is out of buoyant equilibrium. Calving then occurs back to the flexion zone. This calving process provides a compelling and complete explanation for the data. Tracking of oblique camera images allows identification and characterisation of the flexion zones and their propagation downglacier. Interpretation of the GPS data and camera data in combination allows us to place constraints on the height of the basal cavity that forms beneath the rotating ice downglacier of the flexion zone before calving. The flexion zones are probably formed by the exploitation of basal crevasses, and theoretical considerations suggest that their propagation is strongly enhanced when the glacier base is deeper than buoyant equilibrium. Thus, this calving mechanism is likely to dominate whenever such geometry occurs and is of increasing importance in Greenland.