Quantifying components of the hydrological cycle in glacierised catchments is important for the assessment of the temporal distribution, quantity and quality of water resources available to downstream regions, especially under a changing climate. However, this assessment requires long time series of observations, which are typically unavailable for remote catchments, such as those in mountainous areas. In this study, we leverage a unique ∼40 year time series of hydrological data recorded in the subarctic glacierised Tarfala catchment (Sweden) to explore temporal trends in the components of the catchment water balance (precipitation, runoff, change in storage, and evaporation), and to assess if water balance residuals are associated with specific hydro-climatic conditions. No significant temporal trends were found in precipitation and storage changes of the glacierised area, but significant increases were found in evaporation and summer discharge (in part attributed to glacier volume losses). The annual water balance could not be perfectly closed, and water losses were on average 112 mm y⁻¹ larger than the water inputs over the study period. Among the water balance components, discharge contributed most to the total water balance uncertainty, and storage surplus due to antecedent meteorological conditions could explain why water losses in specific years exceeded the uncertainty bounds. It is therefore essential to consider legacy effects from previous years when applying water balance calculations in mountainous and/or glacierised catchments.

Geophysical surveys provide an efficient and non-invasive means of studying subsurface conditions in numerous sedimentary settings. In this study, we explore the application of three geophysical methods to a proglacial environment, namely ground penetrating radar (GPR), seismic refraction and multi-channel analysis of surface waves (MASW). We apply these geophysical methods to three glacial landforms with contrasting morphologies and sedimentary characteristics, and we use the various responses to assess the applicability and limitations of each method for these proglacial targets. Our analysis shows that GPR and seismic (refraction and MASW) techniques can provide spatially extensive information on the internal architecture and composition of moraines, but careful survey designs are required to optimise data quality in these geologically complex environments. Based on our findings, we define a number of recommendations and a potential workflow to guide future geophysical investigations in analogous settings. We recommend the initial use of GPR in future studies of proglacial environments to inform (a) seismic survey design and (b) the selection of seismic interpretation techniques. We show the benefits of using multiple GPR antenna frequencies (e.g., 25 and 100 MHz) to provide decimetre scale imaging in the near surface (e.g., < 15 m) while also enabling signal penetration to targets at up to similar to 40 m depth (e.g., bedrock). This strategy helps to circumvent changes in radar signal penetration resulting from variations in substrate conductivity or abundant scatterers. Our study also demonstrates the importance of combining multiple geophysical methods together with ground-truthing through sedimentological observations to reduce ambiguity in interpretations. Implementing our recommendations will improve geophysical survey practice in the field of glacial geology and allow geophysical methods to play an increasing role in the interpretation of glacial landforms and sediments.

Estimations of peatland carbon stocks often use generalised values for peat thickness and carbon content. Ground penetrating radar (GPR), a rapid technique for field data collection, has been increasingly demonstrated as an appropriate method of mapping peat thickness. Light Detection and Ranging (LiDAR) data as a method for understanding peatland surface elevation are also becoming more widely available. Reliable mapping and quantification of site-specific carbon stocks (e.g. upland raised bogs) is therefore, becoming increasingly feasible, providing a valuable contribution to regional, national and potentially global carbon stock assessments. This is particularly important because raised bogs, such as those found in South Wales are considerable carbon stores. They are, however, susceptible to climate warming owing to their southerly location within the UK. Accurate estimates of peatland carbon stocks has broader importance because world-wide peatland carbon stores are significant and threatened by climate change, posing a substantial challenge not only due to climate feedbacks if this stored carbon is released into the atmosphere, but also the impact on the other ecosystem services that they provide.

Here, we assess the value of an integrated GPR, LiDAR and Geographic Information System (GIS) approach to improve estimation of regional carbon stocks. We apply the approach to three ombrotrophic raised bogs in South Wales, UK, selected for their conservation value and their topographically-confined raised bog form.

GPR and LiDAR are found to be well suited, respectively, to mapping peat thickness at bog scale and surface elevation, thus allowing surface and basal topographies to be evaluated using GIS. In turn, this allows peat volumes to be estimated. For the first time, we record values between 55,200 m³ and 163,000 m³ for the sites considered here.

The greater confidence in these peat volume estimates results from the ability to calibrate the GPR velocity using a depth-to-target calibration with peat cores extracted at locations encompassing the deepest bog area. Peat thickness is mapped at the bog scale with near centimetre precision, improving the robustness of subsequent volume calculations and our understanding of the contribution of these small but numerous sites to regional carbon stocks. Our evaluation shows that GPR corresponds well with conventional manual probing but is minimally invasive and therefore less disturbing of sensitive peatland sites, while also offering improved coverage and spatial resolution with less time and cost.

In combination with measured bulk density and organic carbon contents, these peat volumes allow carbon stocks to be estimated with greater confidence compared to conventional approaches, having values between 2181 ± 122 tonnes carbon and 6305 ± 351 tonnes carbon at our three sites.

Because climate change can affect the carbon balance and hydrology in permafrost peatlands, a better understanding of their sensitivity to changes in temperature and precipitation is needed. In Tavvavuoma, northernmost Sweden, meteorological parameters and ground thermal properties have been monitored in a peat plateau from 2006 to 2013. During this time period, the air temperature record shows no warming trend, and the late-season thaw depth has been relatively stable at around 55-60cm. Meanwhile, the mean annual ground temperature at 1m depth has increased by 0.06 degrees C/yr and at 2-5m depth the permafrost is currently warmer than -0.3 degrees C. Statistical analyses suggest that interannual changes in thaw depth and ground temperatures are affected by different meteorological factors. Summer air temperatures and annual thawing degree-days control thaw depth (p0.05), whereas winter precipitation/snow depth affects ground temperatures (p0.1). The permafrost in this peat plateau is likely relict and not in equilibrium with current climatic conditions. Since the early 20(th) century, there has been a regional increase in air temperature and snow depth. If the ongoing permafrost warming in Tavvavuoma is a result of these long-term trends, short-term variability in meteorological parameters can still have an impact on the rate of permafrost degradation, but unless pronounced climate cooling occurs, thawing of the peat plateau is inevitable.

In this study, snow particle size variability was investigated along a transect in Dronning Maud Land from the coast to the polar plateau. The aim of the study was to better understand the spatial and temporal variations in surface snow properties. Samples were collected twice daily during a traverse in 2007-08 to capture regional variability. Local variability was assessed by sampling in 10 x 10m grids (5m spacing) at selected locations. The particle size and shape distributions for each site were analysed through digital image analysis. Snow particle size variability is complex at different scales, and shows an internal variability of 0.18-3.31 mm depending on the sample type (surface, grid or pit). Relationships were verified between particle size and both elevation and distance to the coast (moisture source). Regional seasonal changes were also identified, particularly on the lower elevations of the polar plateau. This dataset may be used to quantitatively analyse the optical properties of surface snow for remote sensing. The details of the spatial and temporal variations observed in our data provide a basis for further studies of the complex and coupled processes affecting snow particle size and the interpretation of remote sensing of snow covered areas.

Over 11 000 L of kerosene was deposited on the surface of Rabots glaciar on the Kebnekaise Massif, northern Sweden, following the crash of a Royal Norwegian Air Force aircraft in March 2012. An environmental monitoring programme was subsequently commissioned, including a series of dye tracing experiments during the 2013 melt season, conducted to investigate the transport of pollutants through the glacier hydrological system. This experimental set-up provided a basis from which we could gain new insight into the internal hydrological system of Rabots glaciar. Results of dye tracing experiments reveal a degree of homogeneity in the topology of the drainage system throughout July and August, with an increase in efficiency as the season progresses, as reflected by decreasing temporary storage and dispersivity. Early onset of melting likely led to formation of an efficient, discrete drainage system early in the melt season, subject to decreasing sinuosity and braiding as the season progressed. Four distinct meltwater flow regimes are identified to summarize the temporal and spatial evolution of the system. Analysis of turbidity-discharge hysteresis further supports the formation of discrete, efficient drainage, with clockwise diurnal hysteresis suggesting easy mobilization of readily available sediments in channels. Dye injection immediately downstream of the pollution source zone reveals prolonged storage of dye followed by fast, efficient release. Twinned with a low dye recovery, and supported by sporadic detection of hydrocarbons in the proglacial river, we suggest that meltwater, and thus pollutants in solution, may be released periodically through an efficient, and likely pressurized, hydrological system within the upper reaches of the glacier.

Ground-penetrating radar has been widely used to map the thermal structure of polythermal glaciers. Hitherto, the cold temperate transition surface (CTS) in radargrams has been identified by a labour-intensive and subjective manual picking method. We introduce a new automatic approach for picking the CTS that uses the difference in signal power exhibited by the cold and temperate ice layers. We compare our automatically computed CTS depths with manual picks. Our results show very good agreement between the two methods in most areas (r(2) > 0.7). RMSEs computed at each trace in two-way travel-time from three test sites range from 14 to 19 ns (2.4-3.2 m). The proposed automated method mostly fails in areas showing a rather gradual transition in signal power at the CTS. In some areas, high power originating from non-water sources is misinterpreted by the automatic picking method as 'temperate ice'.

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of similar to 12 000 km(2) and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Orkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at similar to 1600m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of similar to 500m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology. The combined data set is freely available at doi:10.1594/pangaea.830314.

In this study, summer rainfall contributions to streamflow were quantified in the sub-arctic, 30% glacierized Tarfala (21.7km(2)) catchment in northern Sweden for two non-consecutive summer sampling seasons (2004 and 2011). We used two-component hydrograph separation along with isotope ratios (O-18 and D) of rainwater and daily streamwater samplings to estimate relative fraction and uncertainties (because of laboratory instrumentation, temporal variability and spatial gradients) of source water contributions. We hypothesized that the glacier influence on how rainfall becomes runoff is temporally variable and largely dependent on a combination of the timing of decreasing snow cover on glaciers and the relative moisture storage condition within the catchment. The results indicate that the majority of storm runoff was dominated by pre-event water. However, the average event water contribution during storm events differed slightly between both years with 11% reached in 2004 and 22% in 2011. Event water contributions to runoff generally increased over 2011 the sampling season in both the main stream of Tarfala catchment and in the two pro-glacial streams that drain Storglaciaren (the largest glacier in Tarfala catchment covering 2.9km(2)). We credit both the inter-annual and intra-annual differences in event water contributions to large rainfall events late in the summer melt season, low glacier snow cover and elevated soil moisture due to large antecedent precipitation. Together amplification of these two mechanisms under a warming climate might influence the timing and magnitude of floods, the sediment budget and nutrient cycling in glacierized catchments.