ALBERT

Description: Decadal trends in the Antarctic sea ice extent ultimately controlled by ice-ocean feedback The Cryosphere Discussions, 7, 4585-4632, 2013 Author(s): H. Goosse and V. Zunz The large natural variability of the Antarctic sea ice is a key characteristic of the system that might be responsible for the small positive trend in sea ice extent observed since 1979. In order to gain insight in the processes responsible for this variability, we have analysed in a control simulation performed with a coupled climate model a strong positive ice-ocean feedback that amplifies sea ice variations. When sea ice concentration increases in a region, in particular close to the ice edge, the mixed layer depth tends to decrease. This can be caused by a net inflow of ice and thus of freshwater that stabilizes the water column. Another stabilizing mechanism at interannual time scales that appears more widespread in our simulation is associated with the downward salt transport due to the seasonal cycle of ice formation: brine is released in winter when ice is formed and mixed over a deep layer while the freshwater flux caused by ice melting is included in a shallow layer, resulting in a net vertical transport of salt. Because of this stronger stratification due to the presence of sea ice, more heat is stored at depth in the ocean and the vertical oceanic heat flux is reduced, which contributes to maintain a higher ice extent. This positive feedback is not associated with a particular spatial pattern. Consequently, the spatial distribution of the trend in ice concentration is largely imposed by the wind changes that can provide the initial perturbation. A positive freshwater flux could alternatively be the initial trigger but the amplitude of the final response of the sea ice extent is finally set up by the amplification related to ice-ocean feedback. Initial conditions have also an influence as the chance to have a large increase in ice extent is higher if starting from a state characterized by a low value.

Description: A data set of world-wide glacier length fluctuations The Cryosphere Discussions, 7, 4775-4811, 2013 Author(s): P. W. Leclercq, J. Oerlemans, H. J. Basagic, I. Bushueva, A. J. Cook, and R. Le Bris Glacier fluctuations contribute to variations in sea level and historical glacier length fluctuations are natural indicators of climate change. To study these subjects, long-term information of glacier change is needed. In this paper we present a~data set of global long-term glacier length fluctuations. The data set is a compilation of available information on changes in glacier length world-wide, including both measured and reconstructed glacier length fluctuations. All 471 length series start before 1950 and cover at least four decades. The longest record starts in 1534, but the majority of time series start after 1850. The number of available records decreases again after 1962. The data set has global coverage including records from all continents. However, the Canadian Arctic is not represented in the data set. The glacier length series show relatively small fluctuations until the mid-19th century followed by a global retreat that was strongest in the first half of the 20th century, although large variability in the length change of the different glaciers is observed. During the 20th century, calving glaciers retreated more than land terminating glaciers, but their relative length change was approximately equal. Besides calving, the glacier slope is the most important glacier property determining length change: steep glaciers have retreated less than glaciers with a gentle slope.

Description: Implementation and evaluation of prognostic representations of the optical diameter of snow in the detailed snowpack model SURFEX/ISBA-Crocus The Cryosphere Discussions, 7, 4443-4500, 2013 Author(s): C. M. Carmagnola, S. Morin, M. Lafaysse, F. Domine, B. Lesaffre, Y. Lejeune, G. Picard, and L. Arnaud In the SURFEX/ISBA-Crocus multi-layer snowpack model, the snow microstructure was up to now characterized by the grain size and by semi-empirical shape variables which cannot be measured easily in the field or linked to other relevant snow properties. In this work we introduce a new formulation of snow metamorphism directly based on equations describing the rate of change of the optical diameter ( d opt ). This variable is considered here to be equal to the equivalent sphere optical diameter, which is inversely proportional to the specific surface area (SSA). d opt thus represents quantitatively some of the geometric characteristics of a porous medium. Different prognostic rate equations of d opt , including a re-formulation of the original Crocus scheme and the parametrizations from Taillandier et al. (2007) and Flanner and Zender (2006), were evaluated by comparing their predictions to field measurements carried out at Summit Camp (Greenland) in May and June 2011 and at Col de Porte (French Alps) during the 2009/10 and 2011/12 winter seasons. We focused especially on results in terms of SSA. In addition, we tested the impact of the different formulations on the simulated density profile, the total snow height, the snow water equivalent (SWE) and the surface albedo. Results indicate that all formulations perform well, with median values of the RMSD between measured and simulated SSA lower than 10 m 2 kg −1 . Incorporating the optical diameter as a fully-fledged prognostic variable is an important step forward in the quantitative description of the snow microstructure within snowpack models, because it opens the way to data assimilation of various electromagnetic observations.

Description: Distribution and recent variations of supraglacial lakes on dendritic-type glaciers in the Khan Tengri-Tomur Mountains, Central Asia The Cryosphere Discussions, 7, 4545-4584, 2013 Author(s): Q. Liu, C. Mayer, and S. Liu Supraglacial lakes are widely distributed on glaciers in the Tomur-Khan Tengri Tianshan Mountains, Central Asia. The existence and development of supraglacial lakes play an important role in the ice melting processes and also in the storage and release of glacial melt water. Here we mapped the supraglacial lakes of eight typical debris-covered dendritic-type glaciers around the Tomur-Khan Tengri peaks based on 9 Landsat TM/ETM+ images acquired in the summers of 1990 until 2011. With a lower area limit of 3600 m 2 for a conservative identification of glacial lakes, we mapped 775 supraglacial lakes and 38 marginal glacial lakes in total. Our results indicate that supraglacial lakes (area 〉 3600 m 2 ) in the study region never develop beyond an elevation of about 3850 m a.s.l., 800 m lower than the maximum upper boundary of debris cover (4650 m a.s.l.). The area-elevation distribution shows that lakes are predominantly occured close to the altitude of 3250 m a.s.l., where the clean ice simultaneously disappears. The majority of the supraglacial lakes are found on the Tomur Glacier and the South Inylchek Glacier, two strongly debris-covered dendritic-type glaciers in the region. As for the multi-year variation of lake area, the summer total and mean areas of supraglacial lakes show some variability from 1990 and 2005 but increased noticeably between 2005 and 2011. The mean area of the mapped lakes reached a maximum in 2010. We found that the area of supraglacial lakes is positively correlated to the total precipitation in summer (July to September) but negatively correlated to the mean spring air temperature (April to June). Pre-summer air temperature fluctuations likely have a stronger impact on the different evolution processes of glacial drainage, evolving from unconnected to connected systems, which may lead to the drainage of larger supraglacial lakes and results in shrinkage of the total and mean lake area during the summer.

Description: Fracture-induced softening for large-scale ice dynamics The Cryosphere Discussions, 7, 4501-4544, 2013 Author(s): T. Albrecht and A. Levermann Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by fracture-coupled dynamic processes within the protective ice shelves leading to accelerated ice flow and hence to a sea-level contribution. In order to better understand the role of fractures in ice dynamics we apply a large-scale continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM). To this end we introduce a higher-order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of a non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in damage threshold. As a result of prognostic flow simulations, flow patterns with realistically large across-flow velocity gradients in fracture-weakened regions as seen in observations are reproduced. This model framework is expandable to grounded ice streams and accounts for climate-induced effects on fracturing and hence on the ice-flow dynamics. It further allows for an enhanced fracture-based calving parameterization.

Description: Sea ice detection with space-based LIDAR The Cryosphere Discussions, 7, 4681-4701, 2013 Author(s): S. Rodier, Y. Hu, and M. Vaughan Monitoring long-term climate change in the Polar Regions relies on accurate, detailed and repeatable measurements of geophysical processes and states. These regions are among the Earth's most vulnerable ecosystems, and measurements there have shown rapid changes in the seasonality and the extent of snow and sea ice coverage. The authors have recently developed a promising new technique that uses lidar surface measurements from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission to infer ocean surface ice-water phase. CALIPSO's 532 nm depolarization ratio measurements of the ocean surface are uniquely capable of providing information about the ever-changing sea surface state within the Polar Regions. With the finer resolution of the CALIPSO footprint (90 m diameter, spaced 335 m apart) and its ability to acquire measurements during both daytime and nighttime orbit segments and in the presence of clouds, the CALIPSO sea ice product provides fine-scale information on mixed phase scenes and can be used to assess/validate the estimates of sea-ice concentration currently provided by passive sensors. This paper describes the fundamentals of the CALIPSO sea-ice detection and classification technique. We present retrieval results from a six-year study, which are compared to existing data sets obtained by satellite-based passive remote sensors.

Description: Simulating the role of gravel on the dynamics of permafrost on the Qinghai-Tibetan Plateau The Cryosphere Discussions, 7, 4703-4740, 2013 Author(s): S. Yi, J. Chen, Q. Wu, and Y. Ding Gravel (particle size ≥ 2 mm) is common in soil profiles of the Qinghai-Tibetan Plateau (QTP). It has different thermal and hydrological properties than other fine mineral soils (particle size 〈 2 mm), which may have significant impacts on the thermal and hydrological processes of soil. However, few models have considered gravel. In this study, we implemented the thermal and hydraulic properties of gravel into the Dynamic Organic Soil-Terrestrial Ecosystem Model to develop new schemes to simulate the dynamics of permafrost on the QTP. Results showed that: (1) the widely used Farouki thermal scheme always simulated higher thermal conductivity of frozen soils than unfrozen soils with the same soil water content; therefore it tends to overestimate permafrost thickness strongly; (2) there exists a soil moisture threshold, below which the new set of schemes with gravel simulated smaller thermal conductivity of frozen soils than unfrozen soils; (3) soil with gravel has higher hydraulic conductivity and poorer water retention capability; and simulations with gravel were usually drier than those without gravel; and (4) the new schemes simulated faster upward degradation than downward degradation; and the simulated permafrost thicknesses were sensitive to the fraction of gravel, the gravel size, the thickness of soil with gravel, and the subsurface drainage. To reduce the uncertainties in the projection of permafrost degradation on the QTP, more effort should be made to: (1) developing robust relationships between soil thermal and hydraulic properties and gravel characteristics based on laboratory work; and (2) compiling spatial datasets of the vertical distribution of gravel content based on measurements during drilling or the digging of soil pits.

Description: Ice volume estimates for the Himalaya–Karakoram region: evaluating different methods The Cryosphere Discussions, 7, 4813-4854, 2013 Author(s): H. Frey, H. Machguth, M. Huss, C. Huggel, S. Bajracharya, T. Bolch, A. Kulkarni, A. Linsbauer, N. Salzmann, and M. Stoffel Ice volume estimates are crucial for assessing water reserves stored in glaciers. A variety of different methodologies exist but there is a lack of systematic comparative analysis thereof. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan-Karakoram (HK) region. Here, three volume–area (V–A) relations, a slope-dependent estimation method, and two ice-thickness distribution models are applied to a complete glacier inventory of the HK region. An uncertainty and sensitivity assessment is performed to investigate the influence of the input glacier areas, and model approaches and parameters on the resulting total ice volumes. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 km 3 to 6455 km 3 , depending on the approach. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses while V–A relations show stronger deviations. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation. More ice-thickness measurements are needed to improve models and results in the future.

Description: Black carbon concentrations from a Tibetan Plateau ice core spanning 1843–1982: recent increases due to emissions and glacier melt The Cryosphere Discussions, 7, 4855-4880, 2013 Author(s): M. Jenkins, S. Kaspari, S. Kang, B. Grigholm, and P. A. Mayewski Black carbon (BC) deposited on snow and glacier surfaces can reduce albedo and lead to accelerated melt. An ice core recovered from Guoqu glacier on Mt. Geladaindong and analyzed using a Single Particle Soot Photometer provides the first long-term (1843–1982) record of BC concentrations from the Central Tibetan Plateau. The highest concentrations are observed from 1975–1982, which corresponds to a 2.0-fold and 2.4-fold increase in average and median values, respectively, relative to 1843–1940. BC concentrations post-1940 are also elevated relative to the earlier portion of the record. Causes for the higher BC concentrations include increased regional BC emissions and subsequent deposition, and melt induced enrichment of BC, with the melt potentially accelerated due to the presence of BC at the glacier surface. A qualitative comparison of the BC and Fe (used as a dust proxy) records suggests that if changes in the concentrations of absorbing impurities at the glacier surface have influenced recent glacial melt, the melt may be due to the presence of BC rather than dust. Guoqu glacier has received no net ice accumulation since the 1980s, and is a potential example of a glacier where an increase in the equilibrium line altitude is exposing buried high impurity layers. That BC concentrations in the uppermost layers of the Geladaindong ice core are not substantially higher relative to deeper in the ice core suggests that some of the BC that must have been deposited on Guoqu glacier via wet or dry deposition between 1983 and 2005 has been removed from the surface of the glacier, potentially via supraglacial or englacial meltwater.

Description: 3D-VAR multilayer assimilation of X-band SAR data into a detailed snowpack model The Cryosphere Discussions, 7, 4881-4912, 2013 Author(s): X. V. Phan, L. Ferro-Famil, M. Gay, Y. Durand, M. Dumont, S. Morin, S. Allain, G. D'Urso, and A. Girard We introduce a variational data assimilation scheme to assimilate X-band Synthetic Aperture Radar (SAR) data into a snowpack evolution model. The structure properties of a snowpack, such as snow density and grain optical diameter of each layer, are simulated over a period of time by the snow metamorphism model Crocus, fed by the local reanalysis SAFRAN at a French alpine location. These parameters are used as inputs of an Electromagnetic Backscattering Model (EBM) based on Dense Media Radiative Transfer (DMRT) theory, which calculates the simulated total backscattering coefficient. Next, 3D-VAR data assimilation is implemented in order to minimize the discrepancies between model simulations and observations obtained from SAR acquisitions, by modifying the parameters of a multilayer snowpack calculated by Crocus. The algorithm then reinitializes Crocus with the optimized snowpack structure properties, and therefore allows it to continue the simulation of snowpack evolution where adjustments based on remote sensing data has been taken into account. Results obtained using TerraSAR-X acquisitions on Argentière Glacier (Mont-Blanc massif, French Alps) show the high potential of this method for improving snow cover simulation.