Glacier Discharge Research Articles

Ice flow modelling to constrain the surface mass balance and ice discharge of San Rafael Glacier, Northern Patagonia Icefield

ABSTRACTWe simulate the ice dynamics of the San Rafael Glacier (SRG) in the Northern Patagonia Icefield (46.7°S, 73.5°W), using glacier geometry obtained by airborne gravity measurements. The full-Stokes ice flow model (Elmer/Ice) is initialized using an inverse method to infer the basal friction coefficient from a satellite-derived surface velocity mosaic. The high surface velocities (7.6 km a−1) near the glacier front are explained by low basal shear stresses (<25 kPa). The modelling results suggest that 98% of the surface velocities are due to basal sliding in the fast-flowing glacier tongue (>1 km a−1). We force the model using different surface mass-balance scenarios taken or adapted from previous studies and geodetic elevation changes between 2000 and 2012. Our results suggest that previous estimates of average surface mass balance over the entire glacier (Ḃ) were likely too high, mainly due to an overestimation in the accumulation area. We propose that most of SRG imbalance is due to the large ice discharge (−0.83 ± 0.08 Gt a−1) and a slightly positiveḂ(0.08 ± 0.06 Gt a−1). The committed mass-loss estimate over the next century is −0.34 ± 0.03 Gt a−1. This study demonstrates that surface mass-balance estimates and glacier wastage projections can be improved using a physically based ice flow model.

Palaeoceanographic evolution of the SW Svalbard shelf over the last 14 000 years

The palaeoceanographic evolution of the SW Svalbard shelf west of Hornsund over the last 14 000 years was reconstructed using benthic foraminiferal assemblages, stable oxygen and carbon isotopes, and grain‐size and ice‐rafted debris data. The results reveal the complexity of the feedbacks influencing the shelf environment: the inflow of Atlantic and Arctic waters (AW and ArW, respectively), and the influence of sea ice and tidewater glaciers. The inflow of subsurface AW onto the shelf gradually increased with the first major intrusion at the end of the Bølling‐Allerød. During the Younger Dryas, the shelf was affected by fresh water originating from sea ice and glacier discharge. Glaciomarine conditions prevailed until the earliest Holocene with the intense deliveries of icebergs and meltwater from retreating glaciers and the occasional penetration of AW onto the shelf. Other major intrusions of AW occurred before and after the Preboreal oscillation (early Holocene), which resulted in more dynamic and open‐water conditions. Between 10.5 and 9.7 cal. ka BP, the shelf environment transformed from glaciomarine to open marine conditions. Between c. 9.7 and 6.1 cal. ka BP the AW advection reached its maximum, resulting in a highly dynamic and productive environment. At c. 6.1 cal. ka BP, the inflow of AW onto the Svalbard shelf decreased due to the intensification of the Greenland Gyre and the subduction of AW under the sea‐ice‐bearing ArW. Bioproductivity decreased over the next c. 5500 years. During the Little Ice Age, bioproductivity increased due to favourable conditions in the marginal sea‐ice zone despite the effects of cooling. The renewed advection of AW after AD 1850 started the climate warming trend observed presently. Our findings show that δ18O can be used to reconstruct the dominances of different water‐masses and, with some caution, as a proxy for the presence of sea ice in frontal areas over the northwestern Eurasian shelves.

Effects of increase glacier discharge on phytoplankton bloom dynamics and pelagic geochemistry in a high Arctic fjord

Arctic fjords experience extremely pronounced seasonal variability and spatial heterogeneity associated with changes in ice cover, glacial retreat and the intrusion of continental shelf’s adjacent water masses. Global warming intensifies natural environmental variability on these important systems, yet the regional and global effects of these processes are still poorly understood. In the present study, we examine seasonal and spatial variability in Kongsfjorden, on the western coast of Spitsbergen, Svalbard. We report hydrological, biological, and biogeochemical data collected during spring, summer, and fall 2012. Our results show a strong phytoplankton bloom with the highest chlorophyll a (Chla) levels ever reported in this area, peaking 15.5µg/L during late May and completely dominated by large diatoms at the inner fjord, that may sustain both pelagic and benthic production under weakly stratified conditions at the glacier front. A progressively stronger stratification of the water column during summer and fall was shaped by the intrusion of warm Atlantic water (T>3°C and Sal>34.65) into the fjord at around 100m depth, and by turbid freshwater plumes (T<1°C and Sal<34.65) at the surface due to glacier meltwater input. Biopolymeric carbon fractions and isotopic signatures of the particulate organic material (POM) revealed very fresh and labile material produced during the spring bloom (13C enriched, with values up to −22.7‰ at the highest Chl a peak, and high in carbohydrates and proteins content – up to 167 and 148µg/L, respectively-), and a clear and strong continental signature of the POM present during late summer and fall (13C depleted, with values averaging −26.5‰, and high in lipid content – up to 92µg/L-) when freshwater melting is accentuated. Our data evidence the importance of combining both physical (i.e. water mass dominance) and geochemical (i.e. characteristics of material released by glacier runoff) data in order to understand the timing, intensity and characteristics of the phytoplankton bloom in Kongsfjorden, a continuously changing system due to sustained warming. In a scenario where glacial retreat is predicted to increase the impacts of meltwater discharge and associated delivery of organic and inorganic material to the surrounding waters, special attention is required in order to predict the consequences for Arctic fjords and adjacent shelf ecosystems.

Emerging impact of Greenland meltwater on deepwater formation in the North Atlantic Ocean

The Greenland ice sheet has experienced increasing mass loss since the 1990s1, 2. The enhanced freshwater flux due to both surface melt and outlet glacier discharge is assuming an increasingly important role in the changing freshwater budget of the subarctic Atlantic3. The sustained and increasing freshwater fluxes from Greenland to the surface ocean could lead to a suppression of deep winter convection in the Labrador Sea, with potential ramifications for the strength of the Atlantic meridional overturning circulation4, 5, 6. Here we assess the impact of the increases in the freshwater fluxes, reconstructed with full spatial resolution3, using a global ocean circulation model with a grid spacing fine enough to capture the small-scale, eddying transport processes in the subpolar North Atlantic. Our simulations suggest that the invasion of meltwater from the West Greenland shelf has initiated a gradual freshening trend at the surface of the Labrador Sea. Although the freshening is still smaller than the variability associated with the episodic ‘great salinity anomalies’, the accumulation of meltwater may become large enough to progressively dampen the deep winter convection in the coming years. We conclude that the freshwater anomaly has not yet had a significant impact on the Atlantic meridional overturning circulation.

Seasonal copula models for the analysis of glacier discharge at King George Island, Antarctica

Modelling glacier discharge is an important issue in hydrology and climate research. Glaciers represent a fundamental water resource when melting of ice and snow contributes to runoff. Glaciers are also studied as natural global warming sensors. GLACKMA association has implemented one of their Pilot Experimental Catchment areas at the King George Island in the Antarctica which records values of the liquid discharge from Collins glacier. In this paper, we propose the use of time-varying copula models for analyzing the relationship between air temperature and glacier discharge, which is clearly non constant and non linear through time. A seasonal copula model is defined where both the marginal and copula parameters vary periodically along time following a seasonal dynamic. Full Bayesian inference is performed such that the marginal and copula parameters are estimated in a one single step, in contrast with the usual two-step approach. Bayesian prediction and model selection is also carried out for the proposed model such that Bayesian credible intervals can be obtained for the conditional glacier discharge given a value of the temperature at any given time point. The proposed methodology is illustrated using the GLACKMA real data where there is, in addition, a hydrological year of missing discharge data which were not possible to measure accurately due to problems in the sounding.

Melt season hydrological characteristics of the Parlung No. 4 Glacier, in Gangrigabu Mountains, south‐east Tibetan Plateau

AbstractWith global warming, hazards relating to glacial melt, such as glacial lake outburst floods, are becoming progressively more serious. However, glacial melt processes and their hydrological consequences are very poorly understood. This study collected glacier discharge data from the terminus of the Parlung No. 4 Glacier throughout the melt season (May–October) during 2008, 2010, 2011 and 2012 to study its specific hydrological characteristics. Time series and multivariate regression analyses were employed to investigate the relationships between discharge and meteorological factors involved, as well as their correlation to discharge estimations. The 0‐ to 3‐day time series analysis showed that discharge rates were highly autocorrelated and that discharge was significantly positively correlated to air temperature, vapour pressure and daily incoming shortwave radiation as well as weakly positively correlated to precipitation. A multiple‐regression exponential model using the independent variables of the daily mean temperature and the vapour pressure exclusively was applied to simulate daily discharge in the basin with a high degree of accuracy. On average, July yielded the maximum monthly mean discharge, followed by August. Discharge in July and August accounted for 53% of the total discharge during the main melt season. The daily cycle of discharge changed as the melt season progressed, reflecting hydrological processes and characteristics of snow melt and glacier ice/snow melt, as well as their transitional periods. Subsequently, regular variations in the characteristics of the diurnal cycle of discharge, storage and delay were observed as the melt season progressed. In addition, the reasons behind the inter‐annual variation in the characteristics of discharge and glacier discharge from the Tibetan Plateau and its surrounding areas are compared and discussed. Copyright © 2015 John Wiley &amp; Sons, Ltd.

Greenland Ice Sheet Mass Balance Reconstruction. Part III: Marine Ice Loss and Total Mass Balance (1840–2010)

Abstract Greenland ice sheet mass loss to the marine environment occurs by some combination of iceberg calving and underwater melting (referred to here as marine ice loss, LM). This study quantifies the relation between LM and meltwater runoff (R) at the ice sheet scale. A theoretical basis is presented explaining how variability in R can be expected to govern much of the LM variability over annual to decadal time scales. It is found that R enhances LM through three processes: 1) increased glacier discharge by ice warming–softening and basal lubrication–sliding; 2) increased calving susceptibility through undercutting glacier front geometry and reducing ice integrity; and 3) increased underwater melting from forcing marine convection. Applying a semiempirical LM f(R) parameterization to a surface mass balance reconstruction enables total ice sheet mass budget closure over the 1840–2010 period. The estimated cumulative 171-yr net ice sheet sea level contribution is 25 ± 10 mm, the rise punctuated by periods of ice sheet net mass gain (sea level drawdown) (1893–1900, 1938–47, and 1972–98). The sea level contribution accelerated at 27.6 mm yr−1 century−1 over the entire reconstruction, reaching a peak sea level rise contribution of 6.1 mm decade−1 during 2002–10.

Modelling the effects of land use and climate changes on hydrology in the Ursern Valley, Switzerland

AbstractWhile many studies have been conducted in mountainous catchments to examine the impact of climate change on hydrology, the interactions between climate changes and land use components have largely unknown impacts on hydrology in alpine regions. They need to be given special attention in order to devise possible strategies concerning general development in these regions. Thus, the main aim was to examine the impact of land use (i.e. bushland expansion) and climate changes (i.e. increase of temperature) on hydrology by model simulations. For this purpose, the physically based WaSiM‐ETH model was applied to the catchment of Ursern Valley in the central Alps (191 km2) over the period of 1983−2005. Modelling results showed that the reduction of the mean monthly discharge during the summer period is due primarily to the retreat of snow discharge in time and secondarily to the reduction in the glacier surface area together with its retreat in time, rather than the increase in the evapotranspiration due to the expansion of the “green alder” on the expense of grassland. The significant decrease in summer discharge during July, August and September shows a change in the regime from b‐glacio‐nival to nivo‐glacial. These changes are confirmed by the modeling results that attest to a temporal shift in snowmelt and glacier discharge towards earlier in the year: March, April and May for snowmelt and May and June for glacier discharge. It is expected that the yearly total discharge due to the land use changes will be reduced by 0.6% in the near future, whereas, it will be reduced by about 5% if climate change is also taken into account. Copyright © 2013 John Wiley &amp; Sons, Ltd.

Organic carbon export from the Greenland ice sheet

Glacial meltwater exports a unique type of organic carbon to marine systems, distinct from non-glacially derived riverine export, potentially capable of stimulating downstream marine primary productivity. Here, we describe for the first time the bulk-level dissolved organic carbon (DOC) and particulate organic carbon (POC) isotopic composition of glacial meltwater draining the Greenland ice sheet (GrIS). These data, in conjunction with an earlier study that investigated the molecular-level composition of GrIS dissolved organic matter, collectively describe the concentration, radiocarbon content, and lability of organic carbon in subglacial discharge from a land-terminating outlet glacier during a melt season. By scaling up our measurements across the ice sheet, we estimate that the annual DOC flux from the GrIS (0.08Tg/y) is equivalent to that from a small Arctic river (discharge (Q)<50km3/y), and that the annual POC flux from the GrIS (0.9Tg/y) may be comparable to that of a large Arctic river (Q>200km3/y). The DOC flux is derived primarily from beneath the glacier (subglacial) (>75%) in the early season, and from surface ice-melt (up to 100%) transmitting through the subglacial environment at the peak of the meltseason. The POC flux is primarily derived from the subglacial environment throughout the meltseason. The early season (low flow) glacier discharge contains higher DOC concentrations (0.5–4.1mgL−1), and exports more enriched carbon (Δ14CDOC ∼ −250‰) compared to the peak season (high flow) discharge, when the concentrations are lower (0.1–0.6mgL−1) and the Δ14C is more depleted (Δ14CDOC ∼ −400‰). Conversely, the POC export (1.4–13.2mgL−1, Δ14CPOC ∼ −250‰) shows no temporal variation in either concentration or radiocarbon content throughout the meltseason. Dissolved ion loads in concomitant samples reflected the seasonal evolution of the subglacial drainage system, confirming the influence of subglacial hydrology on the composition of the bulk carbon pools. Based on this work, we conclude that (1) different mechanisms control the DOC and POC flux from glacial systems; (2) chemically-distinct DOC pools are accessed by seasonally-evolving hydrological flow-paths; and (3) the GrIS can deliver labile carbon, which may also be 14C-depleted, to downstream proglacial and marine environments.

Glacier dynamics over the last quarter of a century at Helheim, Kangerdlugssuaq and 14 other major Greenland outlet glaciers

Abstract. The Greenland ice sheet is experiencing increasing rates of mass loss, the majority of which results from changes in discharge from tidewater glaciers. Both atmospheric and ocean drivers have been implicated in these dynamic changes, but understanding the nature of the response has been hampered by the lack of measurements of glacier flow rates predating the recent period of warming. Here, using Landsat-5 data from 1985 onwards, we extend back in time the record of surface velocities and ice-front position for 16 of Greenland's fastest-flowing tidewater glaciers, and compare these to more recent data from Landsat-7 and satellite-borne synthetic-aperture radar. Climate re-analysis data and sea surface temperatures from 1982 show that since 1995 most of Greenland and its surrounding oceans have experienced significant overall warming, and a switch to a warming trend. During the period from 1985 to 1995 when Greenland and the surrounding oceans were not warming, major tidewater outlet glaciers around Greenland, including Kangerdlugssuaq and Helheim, were dynamically stable. Since the mid-1990s, glacier discharge has consistently been both greater and more variable. Together, these observations support the hypothesis that recent dynamic change is a rapid response to climate forcing. Both air and ocean temperatures in this region are predicted to continue to warm, and will therefore likely drive further change in outlet glacier discharge.

Greenland ice sheet albedo feedback: thermodynamics and atmospheric drivers

Abstract. Greenland ice sheet mass loss has accelerated in the past decade responding to combined glacier discharge and surface melt water runoff increases. During summer, absorbed solar energy, modulated at the surface primarily by albedo, is the dominant factor governing surface melt variability in the ablation area. Using satellite-derived surface albedo with calibrated regional climate modeled surface air temperature and surface downward solar irradiance, we determine the spatial dependence and quantitative impact of the ice sheet albedo feedback over 12 summer periods beginning in 2000. We find that, while albedo feedback defined by the change in net solar shortwave flux and temperature over time is positive over 97% of the ice sheet, when defined using paired annual anomalies, a second-order negative feedback is evident over 63% of the accumulation area. This negative feedback damps the accumulation area response to warming due to a positive correlation between snowfall and surface air temperature anomalies. Positive anomaly-gauged feedback concentrated in the ablation area accounts for more than half of the overall increase in melting when satellite-derived melt duration is used to define the timing when net shortwave flux is sunk into melting. Abnormally strong anticyclonic circulation, associated with a persistent summer North Atlantic Oscillation extreme since 2007, enabled three amplifying mechanisms to maximize the albedo feedback: (1) increased warm (south) air advection along the western ice sheet increased surface sensible heating that in turn enhanced snow grain metamorphic rates, further reducing albedo; (2) increased surface downward shortwave flux, leading to more surface heating and further albedo reduction; and (3) reduced snowfall rates sustained low albedo, maximizing surface solar heating, progressively lowering albedo over multiple years. The summer net infrared and solar radiation for the high elevation accumulation area approached positive values during this period. Thus, it is reasonable to expect 100% melt area over the ice sheet within another similar decade of warming.

A theory of glacial quarrying for landscape evolution models

Numerical models are used commonly to study the topographic evolution of glacially eroded landscapes. These models are grounded on a power-law rule that relates bedrock erosion rate to either the sliding speed or discharge of glaciers. This rule, however, is poorly linked to the principal process of glacial erosion, quarrying, in which bedrock blocks are dislodged from the bed by sliding ice. A new model of quarrying allows this erosion rule to be evaluated. As in past quarrying models, ice-bed separation during sliding controls deviatoric stresses in the rock that cause crack growth. Unlike past models, bedrock strength heterogeneity resulting from preglacial fractures is included using a Weibull statistical distribution of rock strength. This strength distribution is predicated on the observation that larger rock bodies have lower strengths because they have a greater probability of containing a large fracture. Results can, indeed, be closely fitted with a power-law erosion rule, but its nonlinearity, the range of sliding speed over which it applies, and erosion rates depend sensitively on bedrock strength heterogeneity and effective pressure. This theory anchors large-scale models of glacial erosion to the primary small-scale process that these models hope to simulate and reinforces recent emphasis on the role of bedrock fractures in accelerating geomorphic processes. Moreover, by linking basal water pressure to erosion rate, the theory can improve efforts with numerical models to study feedbacks between subglacial hydrology and landscape evolution.

Late‐summer thermal regime of a small proglacial lake

AbstractThis study was motivated by an interest in understanding the potential effects of climate change and glacier retreat on late summer water temperatures in alpine areas. Fieldwork was carried out between July and September 2007 at Place Lake, located below Place Glacier in the southern Coast Mountains of British Columbia. Place Lake has an area of 72 000 m2, a single inlet and outlet channel, and an approximate residence time of 4 days. Warming between the inlet and outlet of the lake ranged up to 3 °C and averaged 1.8 °C, which exceeds the amount of warming that occurred over the 1 km reach of Place Creek between the lake outlet and tree line. Over a 23‐day period, net radiation totalled about 210 MJ·m–2, with sensible heat flux adding another 56 MJ m‐2. The latent heat flux consumed about 8% of the surface heat input. The dominant heat sink was the net horizontal advection associated with lake inflow and outflow. Early in the study period, temperatures between the surface and 6‐m depth were dominantly at or above 4 °C and were generally neutral to thermally stable, whereas temperatures decreased with depth below 6 m and exhibited irregular sub‐diurnal variations. The maximum outflow temperature of almost 7 °C occurred in this period. We hypothesize that turbidity currents associated with cold, sediment‐laden glacier discharge formed an underflow and influenced temperatures in the deeper portion of the lake but did not mix with the upper layers. Later in the study period, the lake was dominantly well mixed with some near‐surface stability associated with nocturnal cooling. Further research is required to examine the combined effects of sediment concentrations and thermal processes on mixing in small proglacial lakes to make projections of the consequences of glacier retreat on alpine lake and stream temperatures. Copyright © 2011 John Wiley &amp; Sons, Ltd.

Geological constraints on glacio-isostatic adjustment models of relative sea-level change during deglaciation of Prince Gustav Channel, Antarctic Peninsula

The recent disintegration of Antarctic Peninsula ice shelves, and the associated accelerated discharge and retreat of continental glaciers, has highlighted the necessity of quantifying the current rate of Antarctic ice mass loss and the regional contributions to future sea-level rise. Observations of present day ice mass change need to be corrected for ongoing glacial isostatic adjustment, a process which must be constrained by geological data. However, there are relatively little geological data on the geometry, volume and melt history of the Antarctic Peninsula Ice Sheet (APIS) after Termination 1, and during the Holocene so the glacial isostatic correction remains poorly constrained. To address this we provide field constraints on the timing and rate of APIS deglaciation, and changes in relative sea-level (RSL) for the north-eastern Antarctic Peninsula based on geomorphological evidence of former marine limits, and radiocarbon-dated marine-freshwater transitions from a series of isolation basins at different altitudes on Beak Island. Relative sea-level fell from a maximum of c. 15 m above present at c. 8000 cal yr BP, at a rate of 3.91 mm yr −1 declining to c. 2.11 mm yr −1 between c. 6900–2900 cal yr BP, 1.63 mm yr −1 between c. 2900–1800 cal yr BP, and finally to 0.29 mm yr −1 during the last c. 1800 years. The new Beak Island RSL curve improves the spatial coverage of RSL data in the Antarctic. It is in broad agreement with some glacio-isostatic adjustment models applied to this location, and with work undertaken elsewhere on the Antarctic Peninsula. These geological and RSL constraints from Beak Island imply significant thinning of the north-eastern APIS by the early Holocene. Further, they provide key data for the glacial isostatic correction required by satellite-derived gravity measurements of contemporary ice mass loss, which can be used to better assess the future contribution of the APIS to rising sea-levels.

Analysis of seasonal variations in mass balance and meltwater discharge of the tropical Zongo Glacier by application of a distributed energy balance model

[1] A distributed energy balance model was applied to Zongo Glacier, Bolivia (16°S, 6000–4900 m above sea level, 2.4 km2), to investigate atmospheric forcing that controls seasonal variations in the mass balance and in meltwater discharge of glaciers in the outer tropics. Surface energy fluxes and melt rates were simulated for each 20 × 20 m2 grid cell at an hourly resolution, for the hydrological year 1999–2000, using meteorological measurements in the ablation area. Model outputs were compared to measurements of meltwater discharge, snow cover extent, and albedo at two weather stations set up on the glacier. Changes in melt rate in three distinct seasons were related to snowfall and cloud radiative properties. During the dry season (May to August), the low melt rate was mainly caused by low long-wave emission of the cloudless thin atmosphere found at these high altitudes. From September to December, meltwater discharge increased to its annual maximum caused by an increase in solar radiation, which was close to its summer peak, as well as a decrease in glacier albedo. From January on, melt was reduced by snowfalls in the core wet season via the albedo effect but was maintained thanks to high long-wave emission from convective clouds. The frequent changes in snow cover throughout the long ablation season lead to large vertical mass balance gradients. Annual mass balance depends on the timing and length of the wet season, which interrupts the period of highest melt rates caused by solar radiation.

A Southern Greenland Ice Sheet Glacier Discharge Reconstruction: 1958-2007

The Greenland ice sheet has been observed to discharge more ice volume between 1996 and 2005 in apparent response to surface climate warming at glaciers south of 70° N latitude. We find that melt rates explain 60-90% of the ice discharge increase, with highest apparent sensitivity evident for the southwestern region. An empirical melt sensitivity model is then used to reconstruct total ice discharge variations due to apparent melt sensitivity from the southern Greenland ice sheet over 50 years (1958-2007). Recent increases in reconstructed southern glacier discharge exceed one standard deviation of the 1958-2007 mean in two years: 2003 and 2007. The 2007 estimated ice discharge value nearly exceeds the reconstruction uncertainty determined using Monte Carlo methods. In combination with climatological surface mass balance estimates, the reconstructed ice discharge suggests that Greenland's sea level contribution during 1961 to 2003 resolves over half of the unexplained IPCC [2007] global sea level rise.

Reconstructions of annual discharge and equilibrium line altitude of glaciers at Qilian Shan, northwest China, from 1978 to 2002

AbstractDischarge was calculated from a mountainous area, including discharge from glaciers, in the Qilian Shan (Qilian Mountains) of northwest China. The studied Yingluoxia basin is 9983 km2 in area, with glaciers making up 0·3% of the basin. The calculation method was based on the heat balance, requiring only daily temperature and precipitation. Calculated annual discharge from the basin corresponded well with the observed data. Calculated annual discharge from glaciers was 3·6% of the total discharge from the basin.The temporal trend of the calculated equilibrium line altitude (ELA) at the July 1st Glacier (western side of the Yingluoxia basin) was similar to that of the observed ELA. The calculated annual mass balance of glaciers within the Yingluoxia basin has a larger negative value than the other glaciers in China, as the ratio of accumulation area to the total glacier area in the Yingluoxia basin is much lower than in neighbouring basins to the west. Copyright © 2010 John Wiley &amp; Sons, Ltd.

Estimated impact of black carbon deposition during pre-monsoon season from Nepal Climate Observatory – Pyramid data and snow albedo changes over Himalayan glaciers

Abstract. The possible minimal range of reduction in snow surface albedo due to dry deposition of black carbon (BC) in the pre-monsoon period (March–May) was estimated as a lower bound together with the estimation of its accuracy, based on atmospheric observations at the Nepal Climate Observatory – Pyramid (NCO-P) sited at 5079 m a.s.l. in the Himalayan region. A total BC deposition rate was estimated as 2.89 μg m−2 day−1 providing a total deposition of 266 μg m−2 for March–May at the site, based on a calculation with a minimal deposition velocity of 1.0×10−4 m s−1 with atmospheric data of equivalent BC concentration. Main BC size at NCO-P site was determined as 103.1–669.8 nm by correlation analyses between equivalent BC concentration and particulate size distributions in the atmosphere. The BC deposition from the size distribution data was also estimated. It was found that 8.7% of the estimated dry deposition corresponds to the estimated BC deposition from equivalent BC concentration data. If all the BC is deposited uniformly on the top 2-cm pure snow, the corresponding BC concentration is 26.0–68.2 μg kg−1, assuming snow density variations of 195–512 kg m−3 of Yala Glacier close to NCO-P site. Such a concentration of BC in snow could result in 2.0–5.2% albedo reductions. By assuming these albedo reductions continue throughout the year, and then applying simple numerical experiments with a glacier mass balance model, we estimated reductions would lead to runoff increases of 70–204 mm of water. This runoff is the equivalent of 11.6–33.9% of the annual discharge of a typical Tibetan glacier. Our estimates of BC concentration in snow surface for pre-monsoon season is comparable to those at similar altitudes in the Himalayan region, where glaciers and perpetual snow regions begin, in the vicinity of NCO-P. Our estimates from only BC are likely to represent a lower bound for snow albedo reductions, because we used a fixed slower deposition velocity. In addition, we excluded the effects of atmospheric wind and turbulence, snow aging, dust deposition, and snow albedo feedbacks. This preliminary study represents the first investigation of BC deposition and related albedo on snow, using atmospheric aerosol data observed at the southern slope in the Himalayas.

Subglacial modulation of the hydrograph from glacierized basins

AbstractThe extent of basin glacierization has important implications for the hydrograph in part, because snow, firn and ice impart different delays in water transport through the system. Here, the significance of subglacial drainage morphology in modulating the hydrograph is examined with a one‐dimensional physically based model. The conceptual model of subglacial drainage comprises both ‘fast’ and ‘slow’ elements, respectively associated with summer and winter drainage regimes. The additional possibility of a permeable glacier substrate is taken into account by allowing water transport in a subglacial aquifer. Forced by prescribed rates of melt‐water delivery to the glacier bed, the model predicts glacier discharge by drainage system provenance. The effects of (1) ‘hard’ versus ‘soft’ glacier beds, (2) subsurface permeability and groundwater flow, and (3) glacier geometry are then investigated. Hydrograph character, in the form of peak timing and amplitude, symmetry with respect to the forcing, and the amplitude of diurnal fluctuations, is affected by the partitioning of water through the various flow elements. Hard beds and impermeable substrates maximize the discharge routed through the fast‐drainage system in the simulations, generally resulting in higher seasonal discharge maxima and stronger diurnal variations in discharge. High hydraulic transmissivities, either at the glacier bed or in underlying strata, hinder the development of the fast‐drainage system in the simulations, producing hydrographs of lower amplitude. Glacier geometry has a modest effect, with adverse bed slopes, very thick or very thin ice and short glacier lengths favouring prolonged drainage through the slow system. These results suggest that the morphology and evolution of the subglacial drainage system may play a significant role in determining the character of the hydrograph from glacierized basins. Copyright © 2008 John Wiley &amp; Sons, Ltd.

Pleistocene glacimarine sedimentation on the continental slope off Vancouver Island, British Columbia

Cordilleran Ice Sheet activity has had a significant influence on glacimarine sediment transport and deposition on the western Canadian continental slope resulting in distinct sediment facies. Here we present lithogenic sediment grain size distributions (grains > 250 µm g − 1 , % > 63 µm, % 10 to 60 µm, and % 2 to 4 µm) and mass accumulation rates (MAR) from a marine core (MD02-2496) collected ~ 50 km offshore of Vancouver Island, British Columbia. The lithogenic fraction of sediment deposited on the slope increased as grain size decreased and sediment sorting increased at ~ 30 ka when the Cordilleran Ice Sheet grew into southern British Columbia during the Fraser Glaciation. As the ice sheet advanced onto the continental shelf near core MD02-2496 between 19 and 17 ka, the proximal ice margin dramatically increased lithogenic sediment MARs. Glacimarine sedimentation at this time was associated with cyclic deposition of generally coarse, poorly sorted sediment. An ~ 80-year cyclicity suggests the influence of 1) outburst floods, 2) lunar or 3) solar cycles on tidewater outlet glacier discharge. Rapid ice sheet retreat at ~ 17 ka is indicated by ice rafted debris (IRD) deposition, followed by the deposition of fine, well-sorted silt suggestive of iceberg discharge followed by rock flour in a meltwater plume. Lithogenic sediment MARs decreased dramatically at the Bølling with a brief return to coarse, poorly sorted sediments during the Younger Dryas. Glacimarine sedimentation decreased by ~ 10 ka as the Cordilleran Ice Sheet extent diminished to the areal extent of modern ice fields.