Proglacial River Research Articles

Permafrost as the main factor controlling the fluvial sedimentation style on glaciomarginal fans

This study describes the sedimentary successions of two glaciomarginal fans (=end moraines) deposited during MIS 6 and MIS 2 in the eastern and north-western part of Poland. Sedimentological analyses of the successions indicate they were both deposited close to the ice sheet margin and are dominated by extensive gravelly and sandy sheet beds derived from the upper flow regime with single ice-wedge structures. Qualitative and quantitative analysis of the pancake-like gravelly and sandy sediments was performed using Markov chain analysis. The presence of these, admittedly common, sedimentary successions raises two key points: why sheetfloods are so common on glaciomarginal fans instead of in-channel flows, and whether sedimentary rhythms deposited by sheetfloods represent the full depositional records, or only partial ones due to accompanying erosion. While both points remain unclear, the answers to these questions may be connected to the presence of permafrost, which occupied the foreland of the ice sheet where the glaciomarginal fans were deposited: it can control the type of proglacial river flows on glaciomarginal fans, the depth of erosion and influence the record of sedimentary rhythms.

Sink or source? Methane and carbon dioxide emissions from cryoconite holes, subglacial sediments, and proglacial river runoff during intensive glacier melting on the Tibetan Plateau

High Mountain Asia glaciers are currently ignored in the estimation of global greenhouse gas budgets (e.g., methane (CH4) and carbon dioxide (CO2)). Similar to the Asian Water Tower and Third Pole, the Tibetan Plateau (TP) hosts the largest volume of glaciers outside the polar regions. These glaciers contain large reservoirs of organic carbon that can influence glacial ecosystems under rapid melting. However, no data exist on the current footprint of CH4 and CO2 from glaciers in the TP. Here, we report in situ observations of CH4 and CO2 fluxes for glacial cryoconite holes, subglacial sediments, and proglacial river runoff across the TP. Our results indicate that cryoconite holes and subglacial sediments can accelerate the export of greenhouse gasses during the melting season due to intensive glacier melting. However, to some extent, proglacial river runoff can be a significant sink of atmospheric CO2; this fact was not identified in previous studies. Our findings suggest that variations (source or sink) of greenhouse gasses from TP glacial basins should be considered in regional CH4 and CO2 budgets under climate warming.

Where will widening occur in an outwash braidplain? A new approach to detecting controls on fluvial lateral erosion in a glacierized catchment (north‐western Spitsbergen, Svalbard)

AbstractThe dynamics of fluvial system evolution depend on fluvial processes and their driving forces associated with climatic variations, which affect changes in the morphology of river channels and floodplains. Neither channel slope and morphology, nor the properties of fluvial sediment have previously been considered as determinants of active braidplain widening on outwash plains (formed from valley/alpine glaciers and confined by pre‐existing topography) in the High Arctic region and in the forefields of retreating glaciers. Factors determining widening of braidplain activity of the Waldemar River outwash (north‐western Spitsbergen, Svalbard) were analysed on the basis of geomorphological, sedimentological, glaciological and meteorological research, and indicate significant multiple correlations between meltwater discharge, precipitation, braidplain width, morphology of the braided channel and the textural features of braidplain deposits. The capability of multivariate adaptive regression splines to detect these relationships was described, and threshold values were identified. The results indicate that the rate of active braidplain widening is proportional to the meltwater outflow in proglacial rivers, which can decrease despite a growing rate of glacial ablation. A proposed model enables us to predict zones of braidplain prone to widening activity in the High Arctic (humid) outwash fans and plains as well as fans developed in arid intermountain basins. The necessary conditions to activate braidplain widening processes were (1) spatial changes in the outwash feeding system due to glacier terminus retreat and (2) crossing the thresholds in passive factors (S and d50) controlling lateral erosion intensity. As a result, the braidplain reached a new dynamic equilibrium, in which high‐magnitude–low‐frequency extreme meltwater discharges were of particular importance in terms of braidplain dynamics and are the dominant controls on the pattern of distributary channels.

The Effects of Glacial Cover on Riverine Silicon Isotope Compositions in Chilean Patagonia

Proglacial rivers have been shown to have distinctive silicon (Si) isotope compositions, providing new insights into the mechanisms controlling Si cycling in the subglacial environment and suggesting terrestrial Si isotope exports may have varied between glacial and interglacial periods. However, Si isotope data are currently limited to a small number of glacial systems in the northern hemisphere, and it is unclear how compositions might vary across a spectrum of glacial influence. Using Chilean Patagonia as a unique natural laboratory, we present Si isotope compositions of 0.45 μm filtered (fSi), 0.02 μm filtered (DSi), and reactive amorphous (ASi) fractions from 40 river catchments with variable glacial cover and explore the key controls on Si cycling. The 0.45 μm filtered glacier-fed river samples displayed isotopically light compositions and a positive linear correlation with upstream glacial cover. This relationship was controlled by the inclusion of an isotopically light colloidal-nanoparticulate (0.02–0.45 μm) silicate phase that was only present in glacier-fed rivers and dominated Si budgets in these catchments. This phase was predominately composed of feldspars and its lability in seawater is uncertain, representing a significant unknown in resolving glacial Si isotope exports from this region. When the colloidal-nanoparticulates were removed from solution by ultra-filtration, the resultant DSi isotope compositions of glacier-fed catchments were not isotopically distinct from some non-glacial rivers and exhibited no clear relationship with glacial cover. The colloidal-nanoparticulate concentration of other weathering-sensitive elements (Li, Mg, Ba, Sr) also showed a linear relationship with glacial cover, suggesting that their isotopic compositions could be affected in a similar manner. These findings highlight the benefit of size-fractionated sampling and the need for more research to understand the lability of colloidal-nanoparticulate species, especially in glacier-fed rivers. Finally, we explore the controls on river ASi isotope compositions and show how including these reactive particulate phases is critical to quantifying terrestrial Si isotope budgets, both in Patagonia and other global regions.

Lava–water interaction and hydrothermal activity within the 2014–2015 Holuhraun Lava Flow Field, Iceland

Lava that erupted during the 2014–2015 Holuhraun eruption in Iceland flowed into a proglacial river system, resulting in aqueous cooling of the lava and an ephemeral hydrothermal system. We carried out a monitoring study of this system from 2015 to 2018 to document the cooling of the lava over this time, using thermocouple measurements and data-logging sensors. The heat loss rate from advection through this hydrothermal system in August 2015 was ~5.5 × 108 W; since eruption, aqueous cooling likely accounted for ~1% of the total heat loss from the lava. This estimate excludes steam losses from fumaroles as well as any groundwater that was not released to the surface, and thus is a lower bound. Near the terminus of the flow, advection of heat by flowing water may have locally accounted for tens of percent of the total cooling of that part of the flow. Our data quantify the importance of water cooling for this lava flow and can be compared with models to better understand lava–water interactions more generally. We also provide detailed methods for simple, low-cost monitoring of similar instances in the future.

Ice sheet and palaeoclimate controls on drainage network evolution: an example from Dogger Bank, North Sea

Abstract. Submerged landscapes on continental shelves archive drainage networks formed during periods of sea-level lowstand. The evolution of these postglacial drainage networks also reveals how past climate changes affected the landscape. Ice-marginal and paraglacial drainage networks on low-relief topography are susceptible to reorganisation of water supply, forced by ice-marginal rearrangement, precipitation and temperature variations, and marine inundation. A rare geological archive of climate-driven landscape evolution during the transition from ice-marginal (ca. 23 ka) to a fully submerged marine environment (ca. 8 ka) is preserved at Dogger Bank, in the southern North Sea. In this study, our analysis of high-resolution seismic reflection and cone penetration test data reveal a channel network over a 1330 km2 area that incised glacial and proglacial lake-fill sediments. The channel network sits below coastal and shallow marine sediments and is therefore interpreted to represent a terrestrial drainage network. When mapped out, the channel form morphology reveals two distinct sets. The first set comprises two low-sinuosity, wide (>400 m) channels that contain macroforms of braid and side bars. These channels are interpreted to have originated as proglacial rivers, which drained the ice-sheet margin to the north. The second set of channels (75–200 m wide, with one larger, ∼400 m wide) has higher sinuosity and forms a subdendritic network of tributaries to the proglacial channels. The timing of channel formation lacks chronostratigraphic control. However, the proglacial rivers must have formed as the ice sheet was still on Dogger Bank, before 23 ka, to supply meltwater to the rivers. Ice-sheet retreat from Dogger Bank led to reorganisation of meltwater drainage and abandonment of the proglacial rivers. Palaeoclimate simulations show a cold and dry period at Dogger Bank between 23 and 17 ka. After 17 ka, precipitation increased, and drainage of precipitation formed the second set of channels. The second set of rivers remained active until marine transgression of Dogger Bank at ca. 8.5–8 ka. Overall, this study provides a detailed insight into the evolution of river networks across Dogger Bank and highlights the interplay between external (climate) and internal (local) forcings in drainage network evolution.

Rapid mineral surface weathering beneath the Greenland Ice Sheet shown by radium and uranium isotopes

Glaciated watersheds are regions of intense physical and chemical weathering. In order to gain new insight on subglacial weathering processes, we measured uranium and radium isotopes from a proglacial river draining the Greenland Ice Sheet (GrIS). Time series samples were collected from the spring to mid-summer, a time period during which subglacial drainage pathways are thought to transition from slow-inefficient (distributed) to fast-efficient (channelized) systems. The 228Ra/226Ra activity ratio of the dissolved load varied from 5.2 ± 0.9–16.9 ± 3.6, which was significantly higher than the 228Ra/226Ra ratio of a suspended sediment load sample of 2.1 ± 0.07 and crustal values of ~1. The high 228Ra/226Ra in the dissolved load relative to the source rock material is indicative of mineral surface weathering induced by rapid and continuous flushing of the subglacial drainage network during the course of the melt season and those prior. The 234U/238U ratio (δ234U) varied between 33 and 106‰ with a discharge-weighted mean of 67‰; the seasonal evolution of δ234U did not correlate to geochemical indicators of subglacial meltwater storage time. An experiment designed to measure changes in δ234U with increasing meltwater storage times found that δ234U in the dissolved phase decreased rapidly with increasing storage time. Similarly, samples collected along a transect moving downstream from the ice sheet terminus had decreasing δ234U values from 63 to 15‰ further indicating that with increased weathering, the δ234U of meltwater decreases. Coupled with the relatively low δ234U and high 228Ra/226Ra, U appears to be impacted by rapid chemical weathering of subglacial and suspended sediments. The Leverett River discharge weighted U concentration was 0.13 nM; if this system is considered representative of the broader GrIS, then the total dissolved U flux from the GrIS would be on the order of 6.4 × 104 mol/y. Using a similar set of assumptions, the dissolved 228Ra and 226Ra flux from the GrIS was ~1.1 × 1014 dpm/y and ~ 5.5 × 1013 dpm/y, respectively. These estimates suggest that the 226Ra flux to the ocean from the GrIS is globally significant and that the 228Ra flux in particular is larger than most river inputs.

Direct Observation of Winter Meltwater Drainage From the Greenland Ice Sheet

AbstractMeltwater runoff from the Greenland Ice Sheet (GrIS) significantly contributes to sea level rise and is the dominant driver of enhanced mass loss. While most melt occurs during summer, little is known about its seasonal and/or interannual retention within the GrIS. Here, we document evidence of runoff during winter, ~4 months after summer melt. Ground‐penetrating radar and borehole surveys in the proglacial Isortoq River reveal slowly flowing water beneath >0.5 m of river ice. Geochemical analysis of this water indicates previous contact with the ice sheet bed. Comparable surveys in proglacial rivers draining four neighboring catchments found no winter drainage, despite a brief surface melt event ~10 days prior. We attribute the observed runoff to residual meltwater storage and release enabled by a >600 m deep trough beneath Isunguata Sermia, but not neighboring glaciers. We conclude that the GrIS bed can stay wet and drain small amounts of meltwater year‐round.

The hazardous 2017\u20132019 surge and river damming by Shispare Glacier, Karakoram

In 2017–2019 a surge of Shispare Glacier, a former tributary of the once larger Hasanabad Glacier (Hunza region), dammed the proglacial river of Muchuhar Glacier, which formed an ice-dammed lake and generated a small Glacial Lake Outburst Flood (GLOF). Surge movement produced the highest recorded Karakoram glacier surface flow rate using feature tracking (~18 ± 0.5 m d−1) and resulted in a glacier frontal advance of 1495 ± 47 m. The surge speed was less than reports of earlier Hasanabad advances during 1892/93 (9.3 km) and 1903 (9.7 km). Surges also occurred in 1973 and 2000–2001. Recent surges and lake evolution are examined using feature tracking in satellite images (1990–2019), DEM differencing (1973–2019), and thermal satellite data (2000–2019). The recent active phase of Shispare surge began in April 2018, showed two surface flow maxima in June 2018 and May 2019, and terminated following a GLOF on 22–23 June 2019. The surge likely had hydrological controls influenced in winter by compromised subglacial flow and low meltwater production. It terminated during summer probably because increased meltwater restored efficient channelized flow. We also identify considerable heterogeneity of movement, including spring/summer accelerations.

Intense Chemical Weathering at Glacial Meltwater-Dominated Hailuogou Basin in the Southeastern Tibetan Plateau

Climate warming has caused rapid shrinkage of glaciers in the Tibetan Plateau (TP), but the impact of glacier retreat on the chemical denudation rate remains largely unknown at the temperate glacial basins. The chemical weathering processes were examined at a temperate glacial basin (HLG) in the southeastern TP based on comprehensive data from the supraglacial meltwater, proglacial river water, precipitation and groundwater over two glacier melt seasons in 2008 and 2013. The concentrations of major ions and suspended sediments in river water exhibit a pronounced seasonality and display a close relationship with river discharge, suggesting a strong hydrological control on the chemical and physical weathering processes. Runoff chemistry is dominated by carbonate weathering and sulfide oxidation. HCO3−, Ca2+, and/or SO42− are the dominant ions in meltwater, river water, precipitation and groundwater. For river water, HCO3− and Ca2+ primarily come from calcite weathering, and SO42− is mainly derived from pyrite oxidation. Both solute and sediment fluxes are positively related to river discharge (r = 0.69, p < 0.01 for sediments). The solute flux and yields are 18,095–19,435 t·year−1 and 225–241 t·km−2·year−1, and the sediment load and yields are 126,390 t·year−1 and 1570 t·km−2·year−1, respectively. The solute yields, cationic denudation rate (CDR; 2850–3108 Σ*meq+ m−2·year−1) and chemical weathering intensity (CWI; 616–711 Σ*meq+ m−3·year−1) at HLG are higher than those at most basins irrespective of the lithology, suggesting more intense weathering in the TP in comparison to other glacial basins worldwide.

Contemporary changes of the channel pattern and braided gravel-bed floodplain under rapid small valley glacier recession (Scott River catchment, Spitsbergen)

After the Little Ice Age, high-latitude proglacial river catchments have been observed to transition to paraglacial conditions through the activation of fluvial processes and changes in the river valley and channel forms. This project provides an analysis of contemporary changes in the morphology of the valley bottom and channel development of the proglacial Scott River catchment (NW part of the Wedel Jarlsberg Land, SW Spitsbergen) based on observations conducted during the 2009–2013 melt seasons. Almost half of the catchment (area of 10 km2) is occupied by a valley glacier, currently in the stage of rapid recession. The processes driving the development of the glacier-free valley floor are determined by (1) long-term factors such as geological structure, glaciotectonics, or the morphogenesis of the terminoglacial zone; (2) medium-term factors, including climatic changes and changes in the extent of the glacial terminus and hydrological regime; and (3) short-term factors such as weather-driven flooding events. Geological factors determined the development of two valley narrowings, and a clear division of the valley floor into three morphological zones. The upper gorge in the terminal moraine rampart distinguishes the wide (up to 700 m) upper section of the valley bottom dominated by an outwash plain, where the river system is fed by small tributaries. The lower valley gorge dissecting the elevated marine terrace separates the middle section of the braided river floodplain from the mouth section (alluvial fan). The geomorphological field mapping of the valley floor and riverbed was conducted using geodesic and hydrometric measurements with high-resolution survey techniques based on the global positioning system (GPS). The detail measurements of the Scott River valley and channel were carried out at representative cross-sections and longitudinal profiles in various parts of the catchment. Because of the dominance of the glacial alimentation regime, changes in the discharge rate and bedload transport regime are influenced by the glacial ablation rate, which determines the variability of channel morphology and the contemporary development of the valley bottom. The upper section of the valley is dominated by a multiple-channel braided system fed by sub- and supraglacial waters, while in the middle section the development pattern changes from a meandering single-thread channel to a multiple-channel wandering system to a multiple-channel braided system. The current valley bottom is almost entirely occupied by a braided river floodplain and multiple-channel river bed. Changes in the channel system are related to periods of snowmelt and the related flooding. After the retreat of the snow cover, the channels are highly stable, as their spatial pattern is usually transformed during several days of ablation or ablation-precipitation flood flows. In the river's mouth section, the channel returns to a single-thread meandering pattern before transitioning to a multiple-channel distribution system and finally to a single-thread straight channel. The river below the ebb in the storm bank develops a subaqual prodelta.

Concentrations and loads of DOC, phenols and aldehydes in a proglacial arctic river in relation to hydro-meteorological conditions. A case study from the southern margin of the Bellsund Fjord – SW Spitsbergen

Climate warming accelerates the melting and thawing of cryosphere components. Therefore, it favours the release of contaminants stored in High-Arctic glaciers for many years. The rate of land-based glacier retreat is of particular importance for the hydrological regime of glacial rivers, but also for the chemical composition of their waters. In this study, we examined 84 surface water samples collected during a period of 42 days (melt season 2012) from the Scott River of glacial-hydrological regime (NW Wedel-Jarlsberg Land). The studies of Scott River hydro-meteorological conditions provided indication on the relation between mean air temperatures (T) and mean daily discharge (Q) (r = 0.67). Moreover, they confirmed the role of short-term meteorological events (increasing of air temperature and heavy precipitation) on the acceleration of glacier melt. In this context, the study also demonstrated incidents of the release of a more concentrated pulse of chemical compounds delivered by the Scottbreen valley glacier. The analysis covered fluctuations of physicochemical parameters, and presence of selected organic pollutants (∑phenols, formaldehyde (HCHO)). Based on the correlation matrix analysis, it was found that fluctuations of physicochemical indices (pH, specific electrolytic conductivity (SEC) and dissolved organic carbon (DOC)) were related to changes in T and Q. Principal components analysis (PCA)permitted us to distinguish two principal components (hydro-meteorological and biogeochemical conditions) responsible for 71.1% of the total hydrochemical variation. In terms of the first principal component, it was found that there was a relationship with the intensified ablation of glacier ice and with the second one, the influence of geological and biological conditions of the environment on the chemistry of the river. The combination of these factors directly shaped the hydrochemistry of the proglacial river during the melt season of 2012. The results of organic compound loads, which are presented in this paper provide unique information about the transport of contaminants such as: (DOC) (9.41–579 mgC s−1), ∑phenols (19.7–534 mg s−1), HCHO (34.0–82.3 mg s−1)) in the Scott River. It is worth noting that the differentiation of concentration and mean daily loads of DOC, between the gorge (96.0 mgC s−1) and mouth (99.6 mgC s−1) sections of the river, were related with the activity of a herring gull (Larus argentatus) colony at the river mouth.

Incision and aggradation in proglacial rivers: Post‐Little Ice Age long‐profile adjustments of Southern Iceland outwash plains

AbstractThe retreat of glaciers in response to climate warming leads to substantial changes in meltwater and sediment yield. Glacial shrinkage also induces the emergence and growth of proglacial margin landforms, which strongly affect water and sedimentary transfers from the glacier to the outwash plains. On a decadal timescale, field observations show that outwash plains of retreating glaciers typically exhibit proximal incision, which decreases in magnitude downstream and stops at an inflection point where aggradation begins. Nevertheless, there is a lack of knowledge about the rates and magnitude of this fluvial adjustment and the effects of the proglacial margin configuration on the temperance or the aggravation of this fluvial adjustment to glacier retreat. This paper investigates the proglacial rivers of 14 retreating glaciers in south‐east Iceland over a post‐Little Ice Age timescale, combining fluvial deposits mapping, lichenometric dating, and long‐profile measurements of proglacial fluvial terraces. Our results demonstrate that (a) proximal incision, associated with distal aggradation and downstream migration of the inflection point, is the dominant pattern of proglacial river response to post‐Little Ice Age glacier retreat in Iceland; (b) estimated mean rates of downstream migration of the inflection point range between 5 and 46 m·a−1; and (c) the downstream migration rate of the inflection point is positively correlated with the proportion of proglacial lakes within the glacier foreland. These findings suggest that proglacial margins dominated by proglacial lakes intensify the rates of proximal incision and inflection point migration.

Diurnal and seasonal variation of glacier meltwater hydrochemistry in Qiyi glacierized catchment in Qilian Mountains, Northwest China: implication for chemical weathering

In glacierized catchments, glacier runoff typically shows a strong diurnal cycle in the ablation season (June-September). To elucidate the effect of these processes on the chemical weathering, fresh snowfall and water samples were collected and studied from the supraglacial river, proglacial river, and gauging site in Qiyi glacierized catchment Qilian Mountains, Northwestern China, in the summer of 2011. The pH and electronic conductivity (EC) were determined in the field, and the concentrations of major ions (Na+, K+, Mg2+, Ca2+, Cl-, SO42-, NO3-) were measured. The results indicated that EC linearly increased with increasing distance from the glacial snout, and the concentrations of major ions increased with increasing water-rock interaction time. Along the flow path of the glacier runoff, Na+ and Cl- are more concentrated than other ions in the supraglacial river while Mg2+ and SO42- are more concentrated than other ions at the gauging site. The discharge, pH, EC, and the concentrations of major ions exhibited significant diurnal variation along the flow path. On the other hand, the amplitude of variation diminished from upstream to downstream along the flow path. The chemical weathering rate (Na++K++Mg2++Ca2+) was determined to be 10.9 t/yr/km2. Moreover, further research indicated that the sampling method influenced the assessment of chemical weathering rates. When the sample was collected randomly in one diurnal cycle of hydrography, the estimated ionic flux could deviate -47%~73% based on estimated hourly data. In contrast, if three samples were collected at peak, base flow and the discharge decreasing rate starts to slow down in one diurnal cycle of hydrography, respectively, the deviation would be less than 15%. The smaller the diurnal variation of discharge, the smaller deviation calculated.

Effect of Meteorological Patterns on the Intensity of Streambank Erosion in a Proglacial Gravel-Bed River (Spitsbergen)

Lower parts of proglacial rivers are commonly assumed to be characterised by a multiannual aggradation trend, and streambank erosion is considered to occur rarely and locally. In the years 2009–2013, detailed measurements of channel processes were performed in the Scott River (SW Spitsbergen). More than 60% of its surface area (10 km2) occupies non-glaciated valleys. Since the end of the Little Ice Age, the Scott Glacier has been subject to intensive retreat, resulting in the expansion of the terminoglacial and paraglacial zones. In this area, the Scott River develops an alluvial valley with a proglacial river, which has led to a comparatively low rate of fluvial transport, dominance of suspension over bedload, and the occurrence of various channel patterns. Measurements, performed in the lower course of the valley in two fixed cross-sections of the Scott River channel, document variable annual tendencies with a prevalence of scour over deposition processes in the channel bottom. The balance of scour and fill also differs in particular measurement cross-sections and during the summer season. The maximum erosion indices (1.7 m2) were related to single periods of floods with snow-glacier melt and rainfall origin. The contribution of streambank erosion was usually lower than that of deep erosion both in the annual cycle and during extreme events. The channel-widening index also suggests variable annual (from −1 m to +1 m) and inter-annual tendencies. During a three-day flood from August 2013, in a measurement profile at the mouth of the river, the NNW bank was laterally shifted by as much as 3 m. Annual and inter-seasonal indices of total channel erosion, however, show that changes in the channel-bottom morphology are equalised relatively fast, and in terms of balance the changes usually do not exceed 0.5% of a cross section’s area.

Diurnal, semidiurnal, and fortnightly tidal components in orthotidal proglacial rivers.

The orthotidal rivers are a new concept referring to inland rivers influenced by gravitational tides through the groundwater tides. "Orthotidal signals" is intended to describe tidal signals found in inland streamwaters (with no oceanic input); these tidal signals were locally generated and then exported into streamwaters. Here, we show that orthotidal signals can be found in proglacial rivers due to the gravitational tides affecting the glaciers and their surrounding areas. The gravitational tides act on glacier through earth and atmospheric tides, while the subglacial water is affected in a manner similar to the groundwater tides. We used the wavelet analysis in order to find tidally affected streamwaters. T_TIDE analyses were performed for discovering the tidal constituents. Tidal components with 0.95 confidence level are as follows: O1, PI1, P1, S1, K1, PSI1, M2, T2, S2, K2, and MSf. The amplitude of the diurnal tidal constituents is strongly influenced by the daily thermal cycle. The average amplitude of the semidiurnal tidal constituents is less altered and ranges from 0.0007 to 0.0969m. The lunisolar synodic fortnightly oscillation, found in the time series of the studied river gauges, is a useful signal for detecting orthotidal rivers when using noisier data. The knowledge of the orthotidal oscillations is useful for modeling fine resolution changes in rivers.

Hydrologic exchange and chemical weathering in a proglacial watershed near Kangerlussuaq, west Greenland

The exchange of proglacial river water with active layer pore water could alter water chemical compositions in glacial outwash plains and oceanic solute fluxes. To evaluate effects of this exchange, we sampled Watson River and adjacent pore water during the 2013 melt season at two sandurs in western Greenland; one in Sandflugtdalen and the other near the confluence with Søndre Strømfjord. We measured temperature, specific conductivity, and head gradients between the river and bank over a week-long period at Sandflugtdalen, as well as sediment hydraulic conductivity and chemical compositions of waters from both sites. Specific conductivity of pore water is four to ten times greater than river water as solutes are concentrated from weathering reactions, cryoconcentration, and evaporation. Pore water compositions are predominantly altered by carbonate dissolution and sulfide mineral oxidation. High concentrations of HCO3 and SO4 result from solute recycling and dissolution of secondary Ca-Mg carbonate/sulfate salts initially formed by near-surface evaporation in the summer and at depth by freeze-in of the active layer and cryoconcentration in the winter. High hydraulic conductivity (10−5 to 10−4 m/s) and diurnal fluctuations of river stage during our study caused exchange of river and pore water immediately adjacent to the river channel, with a net loss of river water to the bank. Pore water >6 m from the river continuously flowed away from the river. Approximately 1–8% of the river discharge through the Sandflugtdalen was lost to the river bank during our 6.75 day study based on calculations using Darcy’s Law. Although not sampled, some of this water should discharge to the river during low river stage early and late in the melt season. Elevated pore water solute concentrations in sandurs and water exchange at diurnal and seasonal frequency should impact fluxes of solutes to the ocean, although understanding the magnitude of this effect will require long-term evaluation throughout the melt season.

Meltwater flow through a rapidly deglaciating glacier and foreland catchment system: Virkisjökull, SE Iceland

AbstractVirkisjökull is a rapidly retreating glacier in south-east Iceland. A proglacial lake has formed in the last ten years underlain by buried ice. In this study we estimate water velocities through the glacier, proglacial foreland and proglacial river using tracer tests and continuous meltwater flow measurements. Tracer testing from a glacial moulin to the glacier outlet in September 2013 demonstrated a rapid velocity of 0.58 m s−1. This was comparable to the velocity within the proglacial river, also estimated from tracer testing. A subsequent tracer test from the same glacial moulin under low flow conditions in May 2014 demonstrated a slower velocity of 0.07 m s−1. The glacier outlet river sinks back into the buried ice, and a tracer test from this sink point through the proglacial foreland to the meltwater river beyond the lake indicated a velocity of 0.03 m s−1, suggesting that an ice conduit system within the buried ice is transferring water rapidly beneath the lake. Ground penetrating radar profiles confirm the presence of this buried conduit system. This study provides an example of rapid deglaciation being associated with extensive conduit systems that enable rapid meltwater transfer from glaciers through the proglacial area to meltwater rivers.

Utility of 222Rn as a passive tracer of subglacial distributed system drainage

Water flow beneath the Greenland Ice Sheet (GrIS) has been shown to include slow-inefficient (distributed) and fast-efficient (channelized) drainage systems, in response to meltwater delivery to the bed via both moulins and surface lake drainage. This partitioning between channelized and distributed drainage systems is difficult to quantify yet it plays an important role in bulk meltwater chemistry and glacial velocity, and thus subglacial erosion. Radon-222, which is continuously produced via the decay of 226Ra, accumulates in meltwater that has interacted with rock and sediment. Hence, elevated concentrations of 222Rn should be indicative of meltwater that has flowed through a distributed drainage system network. In the spring and summer of 2011 and 2012, we made hourly 222Rn measurements in the proglacial river of a large outlet glacier of the GrIS (Leverett Glacier, SW Greenland). Radon-222 activities were highest in the early melt season (10–15 dpm L−1), decreasing by a factor of 2–5 (3–5 dpm L−1) following the onset of widespread surface melt. Using a 222Rn mass balance model, we estimate that, on average, greater than 90% of the river 222Rn was sourced from distributed system meltwater. The distributed system 222Rn flux varied on diurnal, weekly, and seasonal time scales with highest fluxes generally occurring on the falling limb of the hydrograph and during expansion of the channelized drainage system. Using laboratory based estimates of distributed system 222Rn, the distributed system water flux generally ranged between 1–5% of the total proglacial river discharge for both seasons. This study provides a promising new method for hydrograph separation in glacial watersheds and for estimating the timing and magnitude of distributed system fluxes expelled at ice sheet margins.

Autogenic incision and terrace formation resulting from abrupt late-glacial base-level fall, lower Chippewa River, Wisconsin, USA

A paucity of research exists regarding the millennial-scale response of inland alluvial streams to abrupt base-level fall. Studies of modern systems indicate that, over short time scales, the response is a diffusion-like process of upstream-propagating incision. In contrast, evidence from the lower Chippewa River (LCR), located in the upper Midwest of the USA, suggests that autogenic controls operating over time scales of several millennia can overwhelm diffusion, resulting in incision that is prolonged and episodic. During the Last Glacial Maximum, the LCR drained the Chippewa Lobe of the Laurentide Ice Sheet to the glacial upper Mississippi River (UMR). As a meltwater stream, it aggraded and filled its valley with glacial outwash, as did its largest tributaries, which were also meltwater streams. Its nonglacial tributaries aggraded, too, filling their valleys with locally derived sediment. During deglaciation, the UMR incised at least twice, abruptly lowering the LCR's base level — ~15m at 16ka or earlier and an additional 40m at ca. 13.4ka. Each of these base-level falls initiated incision of the LCR, led by upstream migrating knickpoints. The propagation of incision has, however, been a lengthy process. The optically stimulated luminescence (OSL) ages of terrace alluvium indicate that, by 13.5ka, incision had advanced up the LCR only 15km, and by 9ka, only 55km. The process has also been episodic, resulting in the formation of fill-cut terraces (inferred from GPR surveys and exposures of terrace alluvium) that are younger and more numerous in the upstream direction. Autogenic increases in sediment load and autogenic bed armoring, the result of periodic tributary-stream rejuvenation and preferential winnowing of fines by the incising river, may have periodically caused knickpoint migration and incision to slow and possibly stop, allowing lateral erosion and floodplain formation to dominate. A decline in sediment flux from stabilizing incised tributary stream systems would have led to renewed knickpoint migration and incision when floods of sufficient magnitude to breach the channel armor occurred. Minimal floodplain development along the upper section of the present-day LCR, along with the channel morphology of an unstable wandering gravel-bed river immediately downstream from it, suggest that the river is still responding to the base-level falls that happened many millennia ago. The autogenic controls on the LCR's response to UMR incision are a direct consequence of the thick fills of noncohesive sediment that accumulated in its valley and the valleys of its tributary streams during the Late Wisconsinan, making the LCR a prime example of a former proglacial river that remains a paraglacial fluvial system.