Positive Mass Balance Research Articles

Reconstructing the annual mass balance of the Echaurren Norte glacier (Central Andes, 33.5° S) using local and regional hydroclimatic data

Abstract. Despite the great number and variety of glaciers in southern South America, in situ glacier mass-balance records are extremely scarce and glacier–climate relationships are still poorly understood in this region. Here we use the longest (> 35 years) and most complete in situ mass-balance record, available for the Echaurren Norte glacier (ECH) in the Andes at ∼ 33.5° S, to develop a minimal glacier surface mass-balance model that relies on nearby monthly precipitation and air temperature data as forcing. This basic model is able to explain 78 % of the variance in the annual glacier mass-balance record over the 1978–2013 calibration period. An attribution assessment identified precipitation variability as the dominant forcing modulating annual mass balances at ECH, with temperature variations likely playing a secondary role. A regionally averaged series of mean annual streamflow records from both sides of the Andes between ∼ 30 and 37° S is then used to estimate, through simple linear regression, this glacier's annual mass-balance variations since 1909. The reconstruction model captures 68 % of the observed glacier mass-balance variability and shows three periods of sustained positive mass balances embedded in an overall negative trend over the past 105 years. The three periods of sustained positive mass balances (centered in the 1920s–1930s, in the 1980s and in the first decade of the 21st century) coincide with several documented glacier advances in this region. Similar trends observed in other shorter glacier mass-balance series suggest that the Echaurren Norte glacier reconstruction is representative of larger-scale conditions and could be useful for more detailed glaciological, hydrological and climatological assessments in this portion of the Andes.

Comparison of modelled- and remote sensing- derived daily snow line altitudes at Ulugh Muztagh, northern Tibetan Plateau

The ice cap Ulugh Muztagh in the central Kunlun Shan at the northern fringe of the Tibetan Plateau is a very isolated region with arid cold conditions. No observational, meteorological or glaciological ground truth data is available. Using the Moderate-resolution Imaging Spectroradiometer (MODIS) Level 1 radiance Swath Data (MOD02QKM) with a spatial resolution of 250 m, transient snow lines during the months of July to September in 2001 to 2014 are derived. Results are used to calibrate the physical based Coupled Snowpack and Ice surface energy and Mass balance model (COSIMA). The model runs on a representative detail region of Ulugh Muztagh (UM) on a digital elevation model with the same spatial resolution as the MODIS bands. In the absence of field observations, the model is driven solely by dynamically downscaled global analysis data from the High Asia Refined analysis (HAR). We compare remote sensing derived and modelled mean regional transient snow line altitudes in the course of consecutive summer seasons in 2008 to 2010. The resulting snow line altitude (SLA) and annual equilibrium line altitude (ELA) proxy of both methods coincide very well in their interannual variability in accordance with interannual variability of climatic conditions. Since SLAs of both methods do not consistently agree on a daily basis a usage of remote sensing derived SLAs for model calibration in the absence of field observation data is only limitedly feasible for daily analysis. ELA approximation using the highest SLA at the end of ablation period may not be applied to UM because the negative winter mass balance (MB) is not reflected in the summer SLA. The study reveals moderate negative MB for UM throughout the modelling period. The mean regional MB of UM accounts for -523±410 mm w.e. a-1 in the modelling period. Hence UM seems not to belong to the area of the ‘Karakorum anomaly’ comprising a region of positive mass balances in recent years which has its centre presumably in the Western Kunlun Shan.

Climatic controls and climate proxy potential of Lewis Glacier, Mt. Kenya

Abstract. The Lewis Glacier on Mt. Kenya is one of the best studied tropical glaciers and has experienced considerable retreat since a maximum extent in the late 19th century (L19). From distributed mass and energy balance modelling, this study evaluates the current sensitivity of the surface mass and energy balance to climatic drivers, explores climate conditions under which the L19 maximum extent might have been sustained, and discusses the potential for using the glacier retreat to quantify climate change. Multi-year meteorological measurements at 4828 m provide data for input, optimization, and evaluation of a spatially distributed glacier mass balance model to quantify the exchanges of energy and mass at the glacier–atmosphere interface. Currently the glacier loses mass due to the imbalance between insufficient accumulation and enhanced melt, because radiative energy gains cannot be compensated by turbulent energy sinks. Exchanging model input data with synthetic climate scenarios, which were sampled from the meteorological measurements and account for coupled climatic variable perturbations, reveals that the current mass balance is most sensitive to changes in atmospheric moisture (via its impact on solid precipitation, cloudiness, and surface albedo). Positive mass balances result from scenarios with an increase of annual (seasonal) accumulation of 30 % (100 %), compared to values observed today, without significant changes in air temperature required. Scenarios with lower air temperatures are drier and associated with lower accumulation and increased net radiation due to reduced cloudiness and albedo. If the scenarios currently producing positive mass balances are applied to the L19 extent, negative mass balances are the result, meaning that the conditions required to sustain the glacier in its L19 extent are not reflected in today's meteorological observations using model parameters optimized for the present-day glacier. Alternatively, a balanced mass budget for the L19 extent can be achieved by changing both climate and optimized gradients (used to extrapolate the meteorological measurements over the glacier) in a manner that implies a distinctly different coupling between the glacier's local surface-air layer and its surrounding boundary layer. This result underlines the difficulty of deriving palaeoclimates for larger glacier extents on the basis of modern measurements of small glaciers.

Characterization and Mass Balance of PCDD/Fs During the Co-combustion of Sewage Sludge in a Grate-type Municipal Solid Waste Incineration

The emission and distribution characteristics of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) were investigated during the co-combustion of sewage sludge (SS) with municipal solid waste (MSW) in a grate-type furnace, South China. On the output side, the stack gas, slag and fly ash with the dioxin levels of 10.75 pg I-TEQ/Nm3, 4.98 pg I-TEQ/g and 294.9 pg I-TEQ/g, respectively, showed similarities in congener profiles. However, the input side presented a significantly different congener profiles compared to that in the output side, which may attribute to the de novo synthesis. The predominant PCDD/Fs contribution (>90%) in fly ash suggested it as the major sink of dioxins in the output side. A positive mass balance of PCDD/Fs with the value of 12.09μg I-TEQ/t-MF was demonstrated.

Geometric evolution of the Horcones Inferior Glacier (Mount Aconcagua, Central Andes) during the 2002–2006 surge

AbstractThe Central Andes of Chile and Argentina (31–35°S) contain a large number and variety of ice masses, but only two surging glaciers have been studied in this region. We analyzed the 2002–2006 surge of the Horcones Inferior Glacier, Mount Aconcagua, Argentina, based on medium spatial resolution (15–30 m) satellite images and digital elevation models. During the buildup phase the glacier was stagnant, with velocities lower than 0.1 m/d. In the active‐phase velocities reached 14 m/d and the glacier front advanced 3.1 km. At the peak of the active phase (2003–2004), the area‐averaged elevation change was −42 m in the reservoir zone (2.53 km2) and +30 m in the receiving zone (3.31 km2). The estimated ice flux through a cross section located at 4175 meter above sea level was 108 m3 during a period of 391 days, a flux that suggests a mean glacier thickness at this location of ~90 m. The depletion phase showed a recovery of the reservoir zone elevation, the down wasting of the receiving zone (−17 m, 2007–2014), and a return to quiescent velocities. The short active phase, the abrupt change in the velocities, and the high level of the proglacial stream indicate a hydrological switch (Alaska type) trigger. The 2002–2006 and 1984–1990 surges of Horcones Inferior were synchronous with the surges of nearby Grande del Nevado Glacier. These events occurred after periods of positive mass balance, so we hypothesize a climate driver.

A study on the precipitation in Izmir, Turkey: chemical composition and source apportionment by receptor models

The chemistry of rainwater was studied in three sites of a metropolitan city of Turkey to determine the spatial variation, sources affecting composition and the extent of their influence. The acidity of rainwater, Ca, Co, Cr, Zn, Ba, V and Ni showed significant spatial variations. The most of measured species were moderately to extremely enriched compared to the local soil, which is likely arisen from the anthropogenic activities. SO42−, K and Ca were found to be predominantly non-sea origin. Two receptor models, Positive Matrix Factorization (PMF) and Chemical Mass Balance (CMB), agreed well on the source contribution estimations of Al, K, Zn, V, Pb and Cr. The PMF and CMB resolved very similar source contribution profiles to the elemental concentrations of Ca, Zn, V, Pb, K and Na, which the majority of those resulted from mineral industries, bell casting+steel industries, fossil fuel burning, biomass burning+sea salt and sea salt, respectively. The PMF resolved the source contributions of long-transported emissions, particularly for Ba, Sr and Mn, which could not be estimated by the CMB.

Reactivation of Kamb Ice Stream tributaries triggers century‐scale reorganization of Siple Coast ice flow in West Antarctica

Ongoing, centennial-scale flow variability within the Ross ice streams of West Antarctica suggests that the present-day positive mass balance in this region may reverse in the future. Here we use a three-dimensional ice sheet model to simulate ice flow in this region over 250 years. The flow responds to changing basal properties, as a subglacial till layer interacts with water transported in an active subglacial hydrological system. We show that a persistent weak bed beneath the tributaries of the dormant Kamb Ice Stream is a source of internal ice flow instability, which reorganizes all ice streams in this region, leading to a reduced (positive) mass balance within decades and a net loss of ice within two centuries. This hitherto unaccounted for flow variability could raise sea level by 5 mm this century. Furthermore, better constraints on future sea level change from this region will require improved estimates of geothermal heat flux and subglacial water transport.

Rising and falling river flows: contrasting signals of climate change and glacier mass balance from the eastern and western Karakoram

Field observations and geodetic measurements suggest that in the Karakoram Mountains, glaciers are either stable or have expanded since 1990, in sharp contrast to glacier retreats that are prevalently observed in the Himalayas and adjoining high-altitude terrains of central Asia. Decreased discharge in the rivers originating from this region is cited as a supporting evidence for this somewhat anomalous phenomenon. Here, we show that river discharge during the melting season of the glaciers in the eastern and western Karakoram, respectively, exhibits rising and falling trends. We have implemented a statistical procedure involving non-parametric tests combined with a benchmark smoothing technique that has proven to be a powerful method for separating the stochastic component from the trend component in a time series. Precipitation patterns determined from ERA-40 and GPCP data indicate that summer-monsoonal precipitation has increased over the Karakoram Mountains in recent decades. Increasing flows in June and July in the eastern Karakoram are due to an increase in summer-monsoonal precipitation. The rising trend of August discharge is due to an increase in the loss of glacier storage at an approximate average rate of 0.186–0.217 mm d-1 year-1 during the period 1973–2010. Moreover, this rate is higher than the rate of increase in monsoonal snowfall during the months of August and September. Therefore, most plausibly, glacier mass balance in the eastern Karakoram is negative. In the western Karakoram, river flows show declining trends for all summer months for the period 1966–2010, corresponding to a rate of increase of glacier storage by approximately 0.552–0.644 mm d-1 year-1, which is also higher than the rate of increase in summer-monsoonal precipitation. The gain of the cryospheric mass in the western Karakoram is in the form of increased thickness of the glaciers and perennial snowpacks instead of areal expansion. This investigation shows two contrasting patterns of trends of river flows that signify both negative and positive mass balance of the Karakoram glaciers. Trends of river flows are spatially and temporally integrated responses of a watershed to changing climate and thereby are important signals of the conditions of the cryospheric component of a watershed where it is highly significant. However, they cannot unequivocally provide indications of the state and fate of the glaciers in the complex hydrometeorological setting of the Karakoram. Extreme caution and care must be exercised in interpreting trends of river discharge in conjunction with climatic data.

The impact of Saharan dust and black carbon on albedo and long-term mass balance of an Alpine glacier

Abstract. Light-absorbing impurities in snow and ice control glacier melt as shortwave radiation represents the main component of the surface energy balance. Here, we investigate the long-term effect of snow impurities, i.e., mineral dust and black carbon (BC), on albedo and glacier mass balance. The analysis was performed over the period 1914–2014 for two sites on Claridenfirn, Swiss Alps, where an outstanding 100-year record of seasonal mass balance measurements is available. Information on atmospheric deposition of mineral dust and BC over the last century was retrieved from two firn/ice cores of high-alpine sites. A combined mass balance and snow/firn layer model was employed to assess the effects of melt and accumulation processes on the impurity concentration at the surface and thus on albedo and glacier mass balance. Compared to pure snow conditions, the presence of Saharan dust and BC lowered the mean annual albedo by 0.04–0.06 depending on the location on the glacier. Consequently, annual melt was increased by 15–19 %, and the mean annual mass balance was reduced by about 280–490 mm w.e. BC clearly dominated absorption which is about 3 times higher than that of mineral dust. The upper site has experienced mainly positive mass balances and impurity layers were continuously buried whereas at the lower site, surface albedo was more strongly influenced by re-exposure of dust and BC-enriched layers due to frequent years with negative mass balances.

Changes in the firn structure of the western Greenland Ice Sheet caused by recent warming

Abstract. Atmospheric warming over the Greenland Ice Sheet during the last 2 decades has increased the amount of surface meltwater production, resulting in the migration of melt and percolation regimes to higher altitudes and an increase in the amount of ice content from refrozen meltwater found in the firn above the superimposed ice zone. Here we present field and airborne radar observations of buried ice layers within the near-surface (0–20 m) firn in western Greenland, obtained from campaigns between 1998 and 2014. We find a sharp increase in firn-ice content in the form of thick widespread layers in the percolation zone, which decreases the capacity of the firn to store meltwater. The estimated total annual ice content retained in the near-surface firn in areas with positive surface mass balance west of the ice divide in Greenland reached a maximum of 74 ± 25 Gt in 2012, compared to the 1958–1999 average of 13 ± 2 Gt, while the percolation zone area more than doubled between 2003 and 2012. Increased melt and column densification resulted in surface lowering averaging −0.80 ± 0.39 m yr−1 between 1800 and 2800 m in the accumulation zone of western Greenland. Since 2007, modeled annual melt and refreezing rates in the percolation zone at elevations below 2100 m surpass the annual snowfall from the previous year, implying that mass gain in the region is retained after melt in the form of refrozen meltwater. If current melt trends over high elevation regions continue, subsequent changes in firn structure will have implications for the hydrology of the ice sheet and related abrupt seasonal densification could become increasingly significant for altimetry-derived ice sheet mass balance estimates.

Modelling of future mass balance changes of Norwegian glaciers by application of a dynamical–statistical model

The long-term behaviour of Norwegian glaciers is reflected by the long mass-balance records provided by the Norwegian Water Resources and Energy Directorate. These show positive annual mass balances in the 1980s and 1990s at maritime glaciers followed by rapid mass loss since 2000. This study assesses the influence of various atmospheric variables on mass changes of selected Norwegian glaciers by correlation- and cross-validated stepwise multiple regression analyses. The atmospheric variables are constructed from reanalyses by the National Centers for Environmental Prediction and the European Centre for Medium-Range Weather Forecasts. Transfer functions determined by the multiple regression are applied to predictors derived from a multi-model ensemble of climate projections to estimate future mass-balance changes until 2100. The statistical relationship to the North Atlantic Oscillation (NAO), the strongest predictor, is highest for maritime glaciers and less for more continental ones. The mass surplus in the 1980s and 1990s can be attributed to a strong NAO phase and lower air temperatures during the ablation season. The mass loss since 2000 can be explained by an increase of summer air temperatures and a slight weakening of the NAO. From 2000 to 2100 the statistical model predicts predicts changes for glaciers in more continental settings of c. −20 m w.e. (water equivalent) or 0.2 m w.e./a. The corresponding range for their more maritime counterparts is −0.5 to +0.2 m w.e./a. Results from Bayesian classification of observed atmospheric states associated with high melt or high accumulation in the past into different simulated climates in the future suggest that climatic conditions towards the end of the twenty-first century favour less winterly accumulation and more ablation in summer. The posterior probabilities for high accumulation at the end of the twenty-first century are typically 1.5–3 times lower than in the twentieth century while the posterior probabilities for high melt are often 1.5–3 times higher at the end of the twenty-first century than in the twentieth and early twenty-first century.

Reconstructing mass balance of Garabashi Glacier (1800–2005) using dendrochronological data

Выполнена первая реконструкции баланса массы ледника Гарабаши по дендрохронологическим данным с 1800 по 2005 г. Лучше всего модель воспроизводит изменчивость годового баланса ледника при использовании хронологий по ширине годичных колец и по максимальной плотности поздней древесины. Изучен вклад каждой древесно-кольцевой хронологии в модель и её климатический отклик. Установлена высокая согласованность реконструкции баланса массы ледника по дендрохронологическим данным с реконструкциями, полученными ранее по метеорологическим рядам.

Monitoring of lake level changes on the Tibetan Plateau and Tian Shan by retracking Cryosat SARIn waveforms

In this study we present, for the first time, the application of Cryosat SARIn mode data to monitor lakes in mountainous areas and to find their water balance. By applying a novel retracker, tailored for lake level observations, we find at least four useable passes for 125 lakes on the Tibetan Plateau and Tian Shan areas over the period February 2012 to January 2014. From these 125 lakes, 30 are passed at least ten times, for which we compute trends and periodic variations, and 16 lakes more then twenty times, for which we additionally apply a water slope correction. This slope correction accounts for geoid inaccuracies or wind effects. We compared the results over two lakes, Langa Co and Bosten, with Jason-2 measurements. Over Langa Co we find an RMS difference of 0.55 m, while for Lake Bosten this is only 0.26 m. For Lake Bosten, the estimated trends, annual and semi-annual variations from the Cryosat and Jason-2 datasets are compared as well. The estimated amplitudes are comparable, while derived phases differ by a few days. Using the trends of all lakes passed at least ten times, a water volume balance of natural lakes is estimated. A loss of 1.51 ± 0.11 km3 y−1 is observed in the lakes in the Tian Shan area. In Tibet, a positive mass balance is estimated of 1.76 ± 0.24 km3 y−1.

Evaluation of Internal Structure, Volume and Mass of Glacial Bodies by Integrated LiDAR and Ground Penetrating Radar Surveys: The Case Study of Canin Eastern Glacieret (Julian Alps, Italy)

We propose an integrated methodology to image the internal structure, evaluate the volume and estimate the densities of different units within ice bodies, useful for more precise mass estimation of very small glaciers. The procedure encompasses light detection and ranging (LiDAR) and ground penetrating radar (GPR) common offset data. The case study is the Canin Eastern Glacieret (CEG), a very small and maritime glacier in the Eastern Alps, and one of the lowermost glaciers of the European Alps. We calculate both volumetric and mass variations of the analysed ice body by integrating GPR measurements with LiDAR surveys acquired in different years (2006 and 2011). Between 2006 and 2011, the area of the glacieret increased from 8,510 to 17,530 m2 with a gain of 9,016 m2. The observed volume increase has been estimated in 96,350 m3 (+97 %), which corresponds to a positive mass balance of 3.89 m w.e.. This quite unusual finding in the present global warming behaviour is mainly due to the above-average winter accumulation (cw) in the considered period. Moreover, the winter season 2008–2009 represented an exceptional event with a cw equal to 13.38 m, the highest of the available record. Thanks to density estimation, we infer the total mass of the CEG at the time of the geophysical surveys, comparing such results with the ones obtained with available empirical equations, observing an important mass gain in the 5 years considered.

Surface velocity and mass balance of Livingston Island ice cap, Antarctica

Abstract. The mass budget of the ice caps surrounding the Antarctica Peninsula and, in particular, the partitioning of its main components are poorly known. Here we approximate frontal ablation (i.e. the sum of mass losses by calving and submarine melt) and surface mass balance of the ice cap of Livingston Island, the second largest island in the South Shetland Islands archipelago, and analyse variations in surface velocity for the period 2007–2011. Velocities are obtained from feature tracking using 25 PALSAR-1 images, and used in conjunction with estimates of glacier ice thicknesses inferred from principles of glacier dynamics and ground-penetrating radar observations to estimate frontal ablation rates by a flux-gate approach. Glacier-wide surface mass-balance rates are approximated from in situ observations on two glaciers of the ice cap. Within the limitations of the large uncertainties mostly due to unknown ice thicknesses at the flux gates, we find that frontal ablation (−509 ± 263 Mt yr−1, equivalent to −0.73 ± 0.38 m w.e. yr−1 over the ice cap area of 697 km2) and surface ablation (−0.73 ± 0.10 m w.e. yr−1) contribute similar shares to total ablation (−1.46 ± 0.39 m w.e. yr−1). Total mass change (δM = −0.67 &plusmn 0.40 m w.e. yr−1) is negative despite a slightly positive surface mass balance (0.06 ± 0.14 m w.e. yr−1). We find large interannual and, for some basins, pronounced seasonal variations in surface velocities at the flux gates, with higher velocities in summer than in winter. Associated variations in frontal ablation (of ~237 Mt yr−1; −0.34 m w.e. yr−1) highlight the importance of taking into account the seasonality in ice velocities when computing frontal ablation with a flux-gate approach.

Are perfluoroalkyl acids in waste water treatment plant effluents the result of primary emissions from the technosphere or of environmental recirculation?

Wastewater treatment plants (WWTP) have been suggested to be one of the major pathways of perfluoroalkyl acids (PFAAs) from the technosphere to the aquatic environment. The origin of PFAAs in WWTP influents is either from current primary emissions or a result of recirculation of PFAAs that have been residing and transported in the environment for several years or decades. Environmental recirculation can then occur when PFAAs from the environment enter the wastewater stream in, e.g., tap water. In this study 13 PFAAs and perfluorooctane sulfonamide were analyzed in tap water as well as WWTP influent, effluent and sludge from three Swedish cities: Bromma (in the metropolitan area of Stockholm), Bollebygd and Umeå. A mass balance of the WWTPs was assembled for each PFAA. Positive mass balances were observed for PFHxA and PFOA in all WWTPs, indicating the presence of precursor compounds in the technosphere. With regard to environmental recirculation, tap water was an important source of PFAAs to the Bromma WWTP influent, contributing >40% for each quantified sulfonic acid and up to 30% for the carboxylic acids. The PFAAs in tap water from Bollebygd and Umeå did not contribute significantly to the PFAA load in the WWTP influents. Our results show that in order to estimate current primary emissions from the technosphere, it may be necessary to correct the PFAA emission rates in WWTP effluents for PFAAs present in tap water, especially in the case of elevated levels in tap water.

Monitoring of Glacier Changes and Response Time in Chorabari Glacier, Central Himalaya, Garhwal, India

Chorabari Glacier (6.6 sq. km) in the Mandakini River basin, a tributary of the River Alaknanda, Central Himalaya, Garhwal (India) has been monitored in terms of its length and frontal area (snout) changes for the period between 1962 and 2012. Global Positioning System, Survey of India toposheet (1 : 50,000) and ground-based measurements were used to obtain the changes in morphology and size of the glacier. The result shows that the frontal area of the glacier has shrunk by 1% and 344 ± 24 m length loss, with an average rate of 6.8 ± 0.5 m a -1 from 1962 to 2012. The observed terminus records of Chorabari Glacier indicate that the positive mass balance can cause terminus advance in about a 17-year timescale. The lag time of glacier signal transferred from accumulation area to the snout by glacier flow is about 562 years. These observations as well as other studies carried out in the region show a significant reduction in glacier area. The increased retreat rate of the glacier snout is probably a direct consequence of global warming.

Life in the Lizard Slow Lane: Gila Monsters Have Low Rates of Energy Use and Water Flux

Field metabolic rates (FMR) and water influx rates (WIR) were measured in Gila Monsters (Heloderma suspectum) in the Mojave Desert of southern Nevada. Gila Monsters had rates of energy use that were less than half of those expected for lizards of their size, including species that live in arid habitats. Free-living Gila Monsters also had comparatively low water flux rates. The annual energy requirements of adult Gila Monsters averaged approximately 3,766 kJ, with two-thirds of the total (2,533 kJ) allocated to activity while above ground, and one-third to resting metabolism (1,233 kJ; mostly while below ground). Unusually low rates of resting metabolism caused relative costs of activity to be high: compared to other lizards, Gila Monsters have among the highest proportional use of energy for activity, relative to total energy use. The average annual water expenditure of adults was approximately 775 mL, and lizards achieved positive mass balance by drinking rain water during the late summer monsoon season....

Transient bright “halos” on the South Polar Residual Cap of Mars: Implications for mass-balance

Spacecraft imaging of Mars’ south polar region during mid-southern summer of Mars year 28 (2007) observed bright halo-like features surrounding many of the pits, scarps and slopes of the heavily eroded carbon dioxide ice of the South Polar Residual Cap (SPRC). These features had not been observed before, and have not been observed since. We report on the results of an observational study of these halos, and spectral modeling of the SPRC surface at the time of their appearance. Image analysis was performed using data from MRO’s Context Camera (CTX), and High Resolution Imaging Science Experiment (HiRISE), as well as images from Mars Global Surveyor’s (MGS) Mars Orbiter Camera (MOC). Data from MRO’s Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) were used for the spectral analysis of the SPRC ice at the time of the halos. These data were compared with a Hapke reflectance model of the surface to constrain their formation mechanism. We find that the unique appearance of the halos is intimately linked to a near-perihelion global dust storm that occurred shortly before they were observed. The combination of vigorous summertime sublimation of carbon dioxide ice from sloped surfaces on the SPRC and simultaneous settling of dust from the global storm, resulted in a sublimation wind that deflected settling dust particles away from the edges of these slopes, keeping these areas relatively free of dust compared to the rest of the cap. The fact that the halos were not exhumed in subsequent years indicates a positive mass-balance for flat portions of the SPRC in those years. A net accumulation mass-balance on flat surfaces of the SPRC is required to preserve the cap, as it is constantly being eroded by the expansion of the pits and scarps that populate its surface.

Rising river flows and glacial mass balance in central Karakoram

Summary Field observations and geodetic measurements suggest that in the central part of the Karakoram Mountains, glaciers are either stable or expanding since 1990, in sharp contrast to glacier retreats that are prevalent in the Himalayas and adjoining high-altitude terrains of central Asia. Decreased discharge in the rivers originating from Karakoram is cited as a supporting evidence of this somewhat anomalous phenomenon. Here, in contrary to those citations, we show that, river discharge during the melting season of the glaciers in central Karakoram has increased from 1985 to 2010. We have implemented a sophisticated statistical procedure involving non-parametric tests combined with a benchmark smoothing technique that has proven to be a powerful method for empirical trend analysis to remove the stochastic component from the trend component in a time series. Melt water from winter snows is the dominant constituent of June and July flows. Glacial melts predominantly contribute to August and partially to September flows, which are controlled by monsoonal snowfall too at elevations approximately >3500 m. For all four summer months, flows increased from 1985 to 2000. August flows, which actually reflect the states of the glaciers, have continued to rise steadily after 2000 at the same rate as those did for the period 1985–2000. However, the rising trends of June and July flows changed to slightly declining trends from 2000 to 2010. These trends most likely indicate drop in winter snowfall over Karakoram and do not provide direct indications about the states of the glaciers. The rising trend of August discharge is due to change in glacial storage at a steadily decreasing rate of approximately 0.04–0.05 mm/day/year for the period 1985–2010. This rate is nearly equal to the rate of increase in precipitation during the summer months over Karakoram Mountains in recent decades as determined from the ERA-40 and GPCP precipitation datasets. Thereby, this is most plausibly, why the glacial mass balance in central Karakoram is nearly neutral. The rising river flows accompanying non-negative glacier mass balance are consistent with predicted future river flows derived from hydrologic modeling coupled with a climate projection suggesting increasing temperature and precipitation with unchanged glacier covers. This investigation reconciles two apparently contradictory observations namely rising river flows and either zero or slightly positive mass balance of central Karakoram glaciers.