ALBERT

Description: The Chemistry CATT-BRAMS model (CCATT-BRAMS 4.5): a regional atmospheric model system for integrated air quality and weather forecasting and research Geoscientific Model Development, 6, 1389-1405, 2013 Author(s): K. M. Longo, S. R. Freitas, M. Pirre, V. Marécal, L. F. Rodrigues, J. Panetta, M. F. Alonso, N. E. Rosário, D. S. Moreira, M. S. Gácita, J. Arteta, R. Fonseca, R. Stockler, D. M. Katsurayama, A. Fazenda, and M. Bela Coupled Chemistry Aerosol-Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CCATT-BRAMS, version 4.5) is an on-line regional chemical transport model designed for local and regional studies of atmospheric chemistry from the surface to the lower stratosphere suitable both for operational and research purposes. It includes gaseous/aqueous chemistry, photochemistry, scavenging and dry deposition. The CCATT-BRAMS model takes advantage of BRAMS-specific development for the tropics/subtropics as well as the recent availability of preprocessing tools for chemical mechanisms and fast codes for photolysis rates. BRAMS includes state-of-the-art physical parameterizations and dynamic formulations to simulate atmospheric circulations down to the meter. This on-line coupling of meteorology and chemistry allows the system to be used for simultaneous weather and chemical composition forecasts as well as potential feedback between the two. The entire system is made of three preprocessing software tools for user-defined chemical mechanisms, aerosol and trace gas emissions fields and the interpolation of initial and boundary conditions for meteorology and chemistry. In this paper, the model description is provided along with the evaluations performed by using observational data obtained from ground-based stations, instruments aboard aircrafts and retrieval from space remote sensing. The evaluation accounts for model applications at different scales from megacities and the Amazon Basin up to the intercontinental region of the Southern Hemisphere.

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

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: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 1: Implementation and verification Geoscientific Model Development, 6, 1481-1491, 2013 Author(s): D. M. Roche A new 18 O stable water isotope scheme is developed for three components of the i LOVECLIM coupled climate model: atmospheric, oceanic and land surface. The equations required to reproduce the fractionation of stable water isotopes in the simplified atmospheric model ECBilt are developed consistently with the moisture scheme. Simplifications in the processes are made to account for the simplified vertical structure including only one moist layer. Implementation of these equations together with a passive tracer scheme for the ocean and a equilibrium fractionation scheme for the land surface leads to the closure of the (isotopic-) water budget in our climate system. Following the implementation, verification of the existence of usual δ 18 O to climatic relationships are performed for the Rayleigh distillation, the Dansgaard relationship and the δ 18 O –salinity relationship. Advantages and caveats of the approach taken are outlined. The isotopic fields simulated are shown to reproduce most expected oxygen-18–climate relationships with the notable exception of the isotopic composition in Antarctica.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 3: A palaeo-perspective based on present-day data–model comparison for oxygen stable isotopes in carbonates Geoscientific Model Development, 6, 1505-1516, 2013 Author(s): T. Caley and D. M. Roche Oxygen stable isotopes (δ 18 O) are among the most useful tools in palaeoclimatology/palaeoceanography. Simulation of oxygen stable isotopes allows testing how the past variability of these isotopes in water can be interpreted. By modelling the proxy directly in the model, the results can also be directly compared with the data. Water isotopes have been implemented in the global three-dimensional model of intermediate complexity i LOVECLIM, allowing fully coupled atmosphere–ocean simulations. In this study, we present the validation of the model results for present-day climate against the global database for oxygen stable isotopes in carbonates. The limitation of the model together with the processes operating in the natural environment reveal the complexity of use the continental calcite-δ 18 O signal of speleothems for a global quantitative data–model comparison exercise. On the contrary, the reconstructed surface ocean calcite-δ 18 O signal in i LOVECLIM does show a very good agreement with the late Holocene database (foraminifers) at the global and regional scales. Our results indicate that temperature and the isotopic composition of the seawater are the main control on the fossil-δ 18 O signal recorded in foraminifer shells when all species are grouped together. Depth habitat, seasonality and other ecological effects play a more significant role when individual species are considered. We argue that a data–model comparison for surface ocean calcite δ 18 O in past climates, such as the Last Glacial Maximum (≈ 21 000 yr), could constitute an interesting tool for mapping the potential shifts of the frontal systems and circulation changes throughout time. Similarly, the potential changes in intermediate oceanic circulation systems in the past could be documented by a data (benthic foraminifers)-model comparison exercise whereas future investigations are necessary in order to quantitatively compare the results with data for the deep ocean.

Description: Assimilating water column and satellite data for marine export production estimation Geoscientific Model Development, 6, 1575-1590, 2013 Author(s): X. Yao and R. Schlitzer Recent advances in satellite retrieval methodology now allow for estimation of particular organic carbon (POC) concentration in ocean surface waters directly from satellite-based optical data. Because of the good coverage, these data reveal small-scale spatial and temporal concentration gradients and document the evolution of surface water POC as well as the underlying driving biogeochemical processes throughout the seasons. Water column nutrient data also reveal biogeochemical activity. However, because of the scarcity of data, the deduction of temporal changes of particle production and export is not possible in most parts of the ocean. Here we present first results from a new study combining both data streams, thereby exploiting the high spatio-temporal resolution of surface POC concentrations from satellite optical sensors with water column nutrient data having sparser coverage but providing information throughout the entire water column. We use a medium-resolution global model with steady-state 3-D circulation that has been optimized by fitting to a large number of hydrographic parameters and tracers, including CFCs and natural radiocarbon. Production and export of POC is allowed to vary monthly, and the magnitudes of the monthly export fluxes are determined by fitting the model to satellite POC data as well as water column nutrient data using the adjoint method. Two cases have been investigated: (1) the production rate of POC is set to be proportional to export production (EP) and the seasonal changes are assumed sinusoidal (meridionally varying amplitude and phase), and (2) the POC production rate is linked to primary production rates (literature). Both cases were run with the same initial state and model settings, and show total cost function decreases of 12 and 95%, respectively. The POC misfit term alone decreased by 75 and 99.8%. The integrated annual global POC exports of the two cases are 9.9 and 12.3 Gt C yr −1 , respectively. Overall, the remaining POC and phosphate misfits of both solutions are considered too large, and the difference fields still exhibit significant systematic geographical patterns. This indicates that the present model runs are too simplistic and do not fully explain the data. Further, more refined model setups are needed.

Description: Snow on the Ross Ice Shelf: comparison of reanalyses and observations from automatic weather stations The Cryosphere, 7, 1399-1410, 2013 Author(s): L. Cohen and S. Dean Snow accumulation measurements from automatic weather stations (AWS) around the Ross Ice Shelf (RIS), Antarctica, are used to provide a new set of ground-based observations which are compared to precipitation from the ECMWF ERA-Interim and NCEP/NCAR Reanalysis-2 datasets. The high temporal resolution of the AWS snow accumulation measurements allow for an event-based comparison of reanalyses precipitation to the in situ observations. Snow accumulation records from nine AWS provide multiple years of accumulation data between 2008 and 2012 over a relatively large, homogeneous region of Antarctica, and also provide the basis for a statistical evaluation of accumulation and precipitation events. The complex effects of wind on snow accumulation (which can both limit and enhance accumulation) complicate the use of the accumulation measurements, but this analysis shows that they can provide a valuable source of ground-based observations for comparisons to modelled precipitation on synoptic timescales. The analysis shows that ERA-Interim reproduces more precipitation events than NCEP-2, and these events correspond to an average 8.2% more precipitation. Significant correlations between reanalyses and AWS event sizes are seen at several stations and show that ERA-Interim consistently produces larger precipitation events than NCEP-2.

Description: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity The Cryosphere, 7, 1473-1480, 2013 Author(s): H. Löwe, F. Riche, and M. Schneebeli Finding relevant microstructural parameters beyond density is a longstanding problem which hinders the formulation of accurate parameterizations of physical properties of snow. Towards a remedy, we address the effective thermal conductivity tensor of snow via anisotropic, second-order bounds. The bound provides an explicit expression for the thermal conductivity and predicts the relevance of a microstructural anisotropy parameter Q , which is given by an integral over the two-point correlation function and unambiguously defined for arbitrary snow structures. For validation we compiled a comprehensive data set of 167 snow samples. The set comprises individual samples of various snow types and entire time series of metamorphism experiments under isothermal and temperature gradient conditions. All samples were digitally reconstructed by micro-computed tomography to perform microstructure-based simulations of heat transport. The incorporation of anisotropy via Q considerably reduces the root mean square error over the usual density-based parameterization. The systematic quantification of anisotropy via the two-point correlation function suggests a generalizable route to incorporate microstructure into snowpack models. We indicate the inter-relation of the conductivity to other properties and outline a potential impact of Q on dielectric constant, permeability and adsorption rate of diffusing species in the pore space.

Description: Influence of high-order mechanics on simulation of glacier response to climate change: insights from Haig Glacier, Canadian Rocky Mountains The Cryosphere, 7, 1527-1541, 2013 Author(s): S. Adhikari and S. J. Marshall Evolution of glaciers in response to climate change has mostly been simulated using simplified dynamical models. Because these models do not account for the influence of high-order physics, corresponding results may exhibit some biases. For Haig Glacier in the Canadian Rocky Mountains, we test this hypothesis by comparing simulation results obtained from 3-D numerical models that deal with different assumptions concerning physics, ranging from simple shear deformation to comprehensive Stokes flow. In glacier retreat scenarios, we find a minimal role of high-order mechanics in glacier evolution, as geometric effects at our site (the presence of an overdeepened bed) result in limited horizontal movement of ice (flow speed on the order of a few meters per year). Consequently, high-order and reduced models all predict that Haig Glacier ceases to exist by ca. 2080 under ongoing climate warming. The influence of high-order mechanics is evident, however, in glacier advance scenarios, where ice speeds are greater and ice dynamical effects become more important. Although similar studies on other glaciers are essential to generalize such findings, we advise that high-order mechanics are important and therefore should be considered while modeling the evolution of active glaciers. Reduced model predictions may be adequate for other glaciologic and topographic settings, particularly where flow speeds are low and where mass balance changes dominate over ice dynamics in determining glacier geometry.

Description: Meteorological drivers of ablation processes on a cold glacier in the semi-arid Andes of Chile The Cryosphere, 7, 1513-1526, 2013 Author(s): S. MacDonell, C. Kinnard, T. Mölg, L. Nicholson, and J. Abermann Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s −1 ), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m −2 ; summer hourly maximum 1354 W m −2 ) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Description: Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates The Cryosphere, 7, 1499-1512, 2013 Author(s): I. Sasgen, H. Konrad, E. R. Ivins, M. R. Van den Broeke, J. L. Bamber, Z. Martinec, and V. Klemann We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 ± 18 Gt yr −1 . Therefore, our AIS mass balance of −114 ± 23 Gt yr −1 is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (−26 ± 3 Gt yr −1 and −127 ± 7 Gt yr −1 , respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 ± 13 Gt yr −1 ), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 ± 7 Gt yr −1 (ca. 0.4 mm yr −1 ) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.

Description: Enhancing the representation of subgrid land surface characteristics in land surface models Geoscientific Model Development, 6, 1609-1622, 2013 Author(s): Y. Ke, L. R. Leung, M. Huang, and H. Li Land surface heterogeneity has long been recognized as important to represent in the land surface models. In most existing land surface models, the spatial variability of surface cover is represented as subgrid composition of multiple surface cover types, although subgrid topography also has major controls on surface processes. In this study, we developed a new subgrid classification method (SGC) that accounts for variability of both topography and vegetation cover. Each model grid cell was represented with a variable number of elevation classes and each elevation class was further described by a variable number of vegetation types optimized for each model grid given a predetermined total number of land response units (LRUs). The subgrid structure of the Community Land Model (CLM) was used to illustrate the newly developed method in this study. Although the new method increases the computational burden in the model simulation compared to the CLM subgrid vegetation representation, it greatly reduced the variations of elevation within each subgrid class and is able to explain at least 80% of the total subgrid plant functional types (PFTs). The new method was also evaluated against two other subgrid methods (SGC1 and SGC2) that assigned fixed numbers of elevation and vegetation classes for each model grid (SGC1: M elevation bands– N PFTs method; SGC2: N PFTs– M elevation bands method). Implemented at five model resolutions (0.1°, 0.25°, 0.5°, 1.0°and 2.0°) with three maximum-allowed total number of LRUs (i.e., N LRU of 24, 18 and 12) over North America (NA), the new method yielded more computationally efficient subgrid representation compared to SGC1 and SGC2, particularly at coarser model resolutions and moderate computational intensity ( N LRU = 18). It also explained the most PFTs and elevation variability that is more homogeneously distributed spatially. The SGC method will be implemented in CLM over the NA continent to assess its impacts on simulating land surface processes.

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: An efficient method to generate a perturbed parameter ensemble of a fully coupled AOGCM without flux-adjustment Geoscientific Model Development, 6, 1447-1462, 2013 Author(s): P. J. Irvine, L. J. Gregoire, D. J. Lunt, and P. J. Valdes We present a simple method to generate a perturbed parameter ensemble (PPE) of a fully-coupled atmosphere-ocean general circulation model (AOGCM), HadCM3, without requiring flux-adjustment. The aim was to produce an ensemble that samples parametric uncertainty in some key variables and gives a plausible representation of the climate. Six atmospheric parameters, a sea-ice parameter and an ocean parameter were jointly perturbed within a reasonable range to generate an initial group of 200 members. To screen out implausible ensemble members, 20 yr pre-industrial control simulations were run and members whose temperature responses to the parameter perturbations were projected to be outside the range of 13.6 ± 2 °C, i.e. near to the observed pre-industrial global mean, were discarded. Twenty-one members, including the standard unperturbed model, were accepted, covering almost the entire span of the eight parameters, challenging the argument that without flux-adjustment parameter ranges would be unduly restricted. This ensemble was used in 2 experiments; an 800 yr pre-industrial and a 150 yr quadrupled CO 2 simulation. The behaviour of the PPE for the pre-industrial control compared well to ERA-40 reanalysis data and the CMIP3 ensemble for a number of surface and atmospheric column variables with the exception of a few members in the Tropics. However, we find that members of the PPE with low values of the entrainment rate coefficient show very large increases in upper tropospheric and stratospheric water vapour concentrations in response to elevated CO 2 and one member showed an implausible nonlinear climate response, and as such will be excluded from future experiments with this ensemble. The outcome of this study is a PPE of a fully-coupled AOGCM which samples parametric uncertainty and a simple methodology which would be applicable to other GCMs.

Description: GAPPARD: a computationally efficient method of approximating gap-scale disturbance in vegetation models Geoscientific Model Development, 6, 1517-1542, 2013 Author(s): M. Scherstjanoi, J. O. Kaplan, E. Thürig, and H. Lischke Models of vegetation dynamics that are designed for application at spatial scales larger than individual forest gaps suffer from several limitations. Typically, either a population average approximation is used that results in unrealistic tree allometry and forest stand structure, or models have a high computational demand because they need to simulate both a series of age-based cohorts and a number of replicate patches to account for stochastic gap-scale disturbances. The detail required by the latter method increases the number of calculations by two to three orders of magnitude compared to the less realistic population average approach. In an effort to increase the efficiency of dynamic vegetation models without sacrificing realism, we developed a new method for simulating stand-replacing disturbances that is both accurate and faster than approaches that use replicate patches. The GAPPARD (approximating GAP model results with a Probabilistic Approach to account for stand Replacing Disturbances) method works by postprocessing the output of deterministic, undisturbed simulations of a cohort-based vegetation model by deriving the distribution of patch ages at any point in time on the basis of a disturbance probability. With this distribution, the expected value of any output variable can be calculated from the output values of the deterministic undisturbed run at the time corresponding to the patch age. To account for temporal changes in model forcing (e.g., as a result of climate change), GAPPARD performs a series of deterministic simulations and interpolates between the results in the postprocessing step. We integrated the GAPPARD method in the vegetation model LPJ-GUESS, and evaluated it in a series of simulations along an altitudinal transect of an inner-Alpine valley. We obtained results very similar to the output of the original LPJ-GUESS model that uses 100 replicate patches, but simulation time was reduced by approximately the factor 10. Our new method is therefore highly suited for rapidly approximating LPJ-GUESS results, and provides the opportunity for future studies over large spatial domains, allows easier parameterization of tree species, faster identification of areas of interesting simulation results, and comparisons with large-scale datasets and results of other forest models.

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: Evidence of meltwater retention within the Greenland ice sheet The Cryosphere, 7, 1433-1445, 2013 Author(s): A. K. Rennermalm, L. C. Smith, V. W. Chu, J. E. Box, R. R. Forster, M. R. Van den Broeke, D. Van As, and S. E. Moustafa Greenland ice sheet mass losses have increased in recent decades with more than half of these attributed to surface meltwater runoff. However, the magnitudes of englacial storage, firn retention, internal refreezing and other hydrologic processes that delay or reduce true water export to the global ocean remain less understood, partly due to a scarcity of in situ measurements. Here, ice sheet surface meltwater runoff and proglacial river discharge between 2008 and 2010 near Kangerlussuaq, southwestern Greenland were used to establish sub- and englacial meltwater storage for a small ice sheet watershed (36–64 km 2 ). This watershed lacks significant potential meltwater storage in firn, surface lakes on the ice sheet and in the proglacial area, and receives limited proglacial precipitation. Thus, ice sheet surface runoff not accounted for by river discharge can reasonably be attributed to retention in sub- and englacial storage. Evidence for meltwater storage within the ice sheet includes (1) characteristic dampened daily river discharge amplitudes relative to ice sheet runoff; (2) three cold-season river discharge anomalies at times with limited ice sheet surface melt, demonstrating that meltwater may be retained up to 1–6 months; (3) annual ice sheet watershed runoff is not balanced by river discharge, and while near water budget closure is possible as much as 54% of melting season ice sheet runoff may not escape to downstream rivers; (4) even the large meltwater retention estimate (54%) is equivalent to less than 1% of the ice sheet volume, which suggests that storage in en- and subglacial cavities and till is plausible. While this study is the first to provide evidence for meltwater retention and delayed release within the Greenland ice sheet, more information is needed to establish how widespread this is along the Greenland ice sheet perimeter.

Description: High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records The Cryosphere, 7, 1481-1498, 2013 Author(s): S. Kutuzov, M. Shahgedanova, V. Mikhalenko, P. Ginot, I. Lavrentiev, and S. Kemp The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

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: Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls The Cryosphere, 7, 1455-1471, 2013 Author(s): A. Rabatel, A. Letréguilly, J.-P. Dedieu, and N. Eckert We present time series of equilibrium-line altitude (ELA) measured from the end-of-summer snow line altitude computed using satellite images, for 43 glaciers in the western Alps over the 1984–2010 period. More than 120 satellite images acquired by Landsat, SPOT and ASTER were used. In parallel, changes in climate variables, summer cumulative positive degree days (CPDD) and winter precipitation, were analyzed over the same time period using 22 weather stations located inside and around the study area. Assuming a continuous linear trend over the study period: (1) the average ELA of the 43 glaciers increased by about 170 m; (2) summer CPDD increased by about 150 PDD at 3000 m a.s.l.; and (3) winter precipitation remained rather stationary. Summer CPDD showed homogeneous spatial and temporal variability; winter precipitation showed homogeneous temporal variability, but some stations showed a slightly different spatial pattern. Regarding ELAs, temporal variability between the 43 glaciers was also homogeneous, but spatially, glaciers in the southern part of the study area differed from glaciers in the northern part, mainly due to a different precipitation pattern. A sensitivity analysis of the ELAs to climate and morpho-topographic variables (elevation, aspect, latitude) highlighted the following: (1) the average ELA over the study period of each glacier is strongly controlled by morpho-topographic variables; and (2) the interannual variability of the ELA is strongly controlled by climate variables, with the observed increasing trend mainly driven by increasing temperatures, even if significant nonlinear, low-frequency fluctuations appear to be driven by winter precipitation anomalies. Finally, we used an expansion of Lliboutry's approach to reconstruct fluctuations in the ELA of any glacier of the study area with respect to morpho-topographic and climate variables, by quantifying their respective weight and the related uncertainties in a consistent manner within a hierarchical Bayesian framework. This method was tested and validated using the ELA measured on the satellite images.

Description: A note on the water budget of temperate glaciers The Cryosphere, 7, 1557-1564, 2013 Author(s): J. Oerlemans In this note, the total dissipative melting in temperate glaciers is studied. The analysis is based on the notion that the dissipation is determined by the loss of potential energy due to the downward motion of mass (ice, snow, meltwater and rain). A mathematical formulation of the dissipation is developed and applied to a simple glacier geometry. In the next step, meltwater production resulting from enhanced ice motion during a glacier surge is calculated. The amount of melt energy available follows directly from the lowering of the centre of gravity of the glacier. To illustrate the concept, schematic calculations are presented for a number of glaciers with different geometric characteristics. Typical dissipative melt rates, expressed as water-layer depth averaged over the glacier, range from a few centimetres per year for smaller glaciers to half a metre per year for Franz Josef Glacier, one of the most active glaciers in the world (in terms of mass turnover). The total generation of meltwater during a surge is typically half a metre. For Variegated Glacier a value of 70 cm is found, for Kongsvegen 20 cm. These values refer to water layer depth averaged over the entire glacier. The melt \textit{rate} depends on the duration of the surge. It is generally an order of magnitude greater than water production by `normal' dissipation. On the other hand, the additional basal melt rate during a surge is comparable in magnitude with the water input from meltwater and precipitation. This suggests that enhanced melting during a surge does not grossly change the total water budget of a glacier. Basal water generated by enhanced sliding is an important ingredient in many theories of glacier surges. It provides a positive feedback mechanism that actually makes the surge happen. The results found here suggest that this can only work if water generated by enhanced sliding accumulates in a part of the glacier base where surface meltwater and rain have no or very limited access. This finding seems compatible with the fact that, on many glaciers, surges are initiated in the lower accumulation zone.

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.

Description: Global glacier changes: a revised assessment of committed mass losses and sampling uncertainties The Cryosphere, 7, 1565-1577, 2013 Author(s): S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp Most glaciers and ice caps (GIC) are out of balance with the current climate. To return to equilibrium, GIC must thin and retreat, losing additional mass and raising sea level. Because glacier observations are sparse and geographically biased, there is an undersampling problem common to all global assessments. Here, we further develop an assessment approach based on accumulation-area ratios (AAR) to estimate committed mass losses and analyze the undersampling problem. We compiled all available AAR observations for 144 GIC from 1971 to 2010, and found that most glaciers and ice caps are farther from balance than previously believed. Accounting for regional and global undersampling errors, our model suggests that GIC are committed to additional losses of 32 ± 12% of their area and 38 ± 16% of their volume if the future climate resembles the climate of the past decade. These losses imply global mean sea-level rise of 163 ± 69 mm, assuming total glacier volume of 430 mm sea-level equivalent. To reduce the large uncertainties in these projections, more long-term glacier measurements are needed in poorly sampled regions.

Description: The physical basis for gas transport through polar firn: a case study at Summit, Greenland The Cryosphere Discussions, 7, 2455-2487, 2013 Author(s): A. C. Adolph and M. R. Albert Compared to other natural porous materials, relatively little is known about the physical nature of polar firn. This intricate network of ice and pore space that comprises the top 60–100 m of the polar ice sheets is the framework that forms the natural archive of past climate information. Despite the many implications for ice core interpretation, direct measurements of physical properties throughout the firn column are limited. Models of gas transport through firn are used to interpret in-situ chemical data which is retrieved to analyze past atmospheric composition. These traditional models treat the firn as a "black box," with gas transport parameters tuned to match gas concentrations with depth to known atmospheric histories. Though this method has been largely successful and provided very useful insights, there are still many questions and uncertainties to be addressed. This work seeks to understand the impact of firn structure on gas transport in firn from a first principles standpoint through direct measurements of permeability, gas diffusivity and microstructure. The relationships between gas transport properties and microstructure will be characterized and compared to existing relationships for general porous media. Direct measurements of gas diffusivity are compared to diffusivities deduced from models based on firn air chemical sampling. Our comparison illuminates the primary importance of including microstructural parameters, beyond just porosity or density, in mass transport modeling, and it provides insights about the nature of gas transport throughout the firn column. Guidance is provided for development of next-generation firn air transport models.

Description: Recent changes in spring snowmelt timing in the Yukon River basin detected by passive microwave satellite data The Cryosphere, 7, 905-916, 2013 Author(s): K. A. Semmens and J. M. Ramage Spring melt is a significant feature of high latitude snowmelt dominated drainage basins influencing hydrological and ecological processes such as snowmelt runoff and green-up. Melt duration, defined as the transition period from snowmelt onset until the end of the melt refreeze, is characterized by high diurnal amplitude variations (DAV) where the snowpack is melting during the day and refreezing at night, after which the snowpack melts constantly until depletion. Determining trends for this critical period is necessary for understanding how the Arctic is changing with rising temperatures and provides a baseline from which to assess future change. To study this dynamic period, brightness temperature ( T b ) data from the Special Sensor Microwave Imager (SSM/I) 37 V-GHz frequency from 1988 to 2010 were used to assess snowmelt timing trends for the Yukon River basin, Alaska/Canada. Annual T b and DAV for 1434 Equal-Area Scalable Earth (EASE)-Grid pixels (25 km resolution) were processed to determine melt onset and melt refreeze dates from T b and DAV thresholds previously established in the region. Temporal and spatial trends in the timing of melt onset and melt refreeze, and the duration of melt were analyzed for the 13 sub-basins of the Yukon River basin with three different time interval approaches. Results show a lengthening of the melt period for the majority of the sub-basins with a significant trend toward later end of melt refreeze after which the snowpack melts day and night leading to snow clearance, peak discharge, and green-up. Earlier melt onset trends were also found in the higher elevations and northernmost sub-basins (Porcupine, Chandalar, and Koyukuk rivers). Latitude and elevation displayed the dominant controls on melt timing variability and spring solar flux was highly correlated with melt timing in middle (∼600–1600 m) elevations.

Description: Solving Richards Equation for snow improves snowpack meltwater runoff estimations The Cryosphere Discussions, 7, 2373-2412, 2013 Author(s): N. Wever, C. Fierz, C. Mitterer, H. Hirashima, and M. Lehning The runoff from the snow cover during spring snow melt or rain-on-snow events is an important factor in the hydrological cycle. In this study, water transport schemes for a 1-dimensional physical based snowpack model are compared to 14 yr of lysimeter measurements at a high alpine site. The schemes include a simple bucket-type approach, an approximation of Richards Equation (RE), and the full RE. The results show that daily sums of runoff are strongly related to a positive energy balance of the snow cover and therefore, all water transport schemes show very similar performance in terms of Nash–Sutcliffe efficiency (NSE) coefficients (around 0.59) and r 2 values (around 0.77). Timing of the arrival of meltwater in spring at the bottom of the snowpack showed differences between the schemes, where especially in the bucket-type and approximated RE approach, meltwater release is slower than in the measurements. Overall, solving RE for the snow cover yields the best agreement between modelled and measured runoff. On sub-daily time scales, the water transport schemes behave very differently. Also here, solving RE provides the highest agreement between modelled and measured runoff in terms of NSE coefficient (0.48), where other water transport schemes loose any predictive power. This appears to be mainly due to bad timing of meltwater release during the day. Accordingly, solving RE for the snow cover improves several aspects of modelling snow cover runoff. The additional computational cost was found to be in the order of a factor of 1.5.

Description: An upper-bound estimate for the accuracy of volume-area scaling The Cryosphere Discussions, 7, 2293-2331, 2013 Author(s): D. Farinotti and M. Huss Volume-area scaling is the most popular method for estimating the ice volume of large glacier samples. Here, a series of resampling experiments based on different sets of synthetic data are presented in order to derive an upper-bound estimate (i.e. a level achieved only with ideal conditions) for the accuracy of its application. We also quantify the maximum accuracy expected when scaling is used for determining the glacier volume change, and area change of a given glacier population. A comprehensive set of measured glacier areas, volumes, area and volume changes is evaluated to investigate the impact of real-world data quality on the so assessed accuracies. For populations larger than a few thousand glaciers, the total ice volume can be recovered within 30% if all measurements available worldwide are used for estimating the scaling coefficients. Assuming no systematic biases in ice volume measurements, their uncertainty is of secondary importance. Knowing the individual areas of a glacier sample for two points in time allows recovering the corresponding ice volume change within 40% for populations larger than a few hundred glaciers, both for steady-state and transient geometries. If ice volume changes can be estimated without bias, glacier area changes derived from volume-area scaling show similar uncertainties as for the volume changes. This paper does not aim at making a final judgement about the suitability of volume-area scaling, but provides the means for assessing the accuracy expected from its application.

Description: Evaluation of the snow regime in dynamic vegetation land surface models using field measurements The Cryosphere Discussions, 7, 2333-2372, 2013 Author(s): E. Kantzas, M. Lomas, S. Quegan, and E. Zakharova An increasing number of studies have demonstrated the significant climatic and ecological changes occurring in the northern latitudes over the past decades. As coupled, earth-system models attempt to describe and simulate the dynamics and complex feedbacks of the Arctic environment, it is important to reduce their uncertainties in short-term predictions by improving the description of both the systems processes and its initial state. This study focuses on snow-related variables and extensively utilizes a historical data set (1966–1996) of field snow measurements acquired across the extend of the Former Soviet Union (FSU) to evaluate a range of simulated snow metrics produced by a variety of land surface models, most of them embedded in IPCC-standard climate models. We reveal model-specific issues in simulating snow dynamics such as magnitude and timings of SWE as well as evolution of snow density. We further employ the field snow measurements alongside novel and model-independent methodologies to extract for the first time (i) a fresh snow density value (57–117 kg m –3 ) for the region and (ii) mean monthly snowpack sublimation estimates across a grassland-dominated western (November–February) [9.2, 6.1, 9.15, 15.25] mm and forested eastern sub-sector (November–March) [1.53, 1.52, 3.05, 3.80, 12.20] mm; we subsequently use the retrieved values to assess relevant model outputs. The discussion session consists of two parts. The first describes a sensitivity study where field data of snow depth and snow density are forced directly into the surface heat exchange formulation of a land surface model to evaluate how inaccuracies in simulating snow metrics affect important modeled variables and carbon fluxes such as soil temperature, thaw depth and soil carbon decomposition. The second part showcases how the field data can be assimilated with ready-available optimization techniques to pinpoint model issues and improve their performance.

Description: Seasonal controls on snow distribution and aerial ablation at the snow-patch and landscape scales, McMurdo Dry Valleys, Antarctica The Cryosphere, 7, 917-931, 2013 Author(s): J. W. Eveland, M. N. Gooseff, D. J. Lampkin, J. E. Barrett, and C. D. Takacs-Vesbach Accumulated snow in the McMurdo Dry Valleys, while limited, has great ecological significance to subnivian soil environments. Though sublimation dominates the ablation process in this region, measurable increases in soil moisture and insulation from temperature extremes provide more favorable conditions with respect to subnivian soil communities. While precipitation is not substantial, significant amounts of snow can accumulate, via wind transport, in topographic lees along the valley bottoms, forming thousands of discontinuous snow patches. These patches have the potential to act as significant sources of local meltwater, controlling biogeochemical cycling and the landscape distribution of microbial communities. Therefore, determining the spatial and temporal dynamics of snow at multiple scales is imperative to understanding the broader ecological role of snow in this region. High-resolution satellite imagery acquired during the 2009–2010 and 2010–2011 austral summers was used to quantify the distribution of snow across Taylor and Wright valleys. Extracted snow-covered area from the imagery was used as the basis for assessing inter-annual variability and seasonal controls on accumulation and ablation of snow at multiple scales. In addition to landscape analyses, fifteen 1 km 2 plots (3 in each of 5 study regions) were selected to assess the prevalence of snow cover at finer spatial scales, referred to herein as the snow-patch scale. Results confirm that snow patches tend to form in the same locations each year with some minor deviations observed. At the snow-patch scale, neighboring patches often exhibit considerable differences in aerial ablation rates, and particular snow patches do not reflect trends for snow-covered area observed at the landscape scale. These differences are presumably related to microtopographic influences acting on individual snow patches, such as wind sheltering and differences in snow depth due to the underlying topography. This highlights the importance of both the landscape and snow-patch scales in assessing the effects of snow cover on biogeochemical cycling and microbial communities.

Description: On the parallelization of atmospheric inversions of CO 2 surface fluxes within a variational framework Geoscientific Model Development, 6, 783-790, 2013 Author(s): F. Chevallier The variational formulation of Bayes' theorem allows inferring CO 2 sources and sinks from atmospheric concentrations at much higher time–space resolution than the ensemble or analytical approaches. However, it usually exhibits limited scalable parallelism. This limitation hinders global atmospheric inversions operated on decadal time scales and regional ones with kilometric spatial scales because of the computational cost of the underlying transport model that has to be run at each iteration of the variational minimization. Here, we introduce a physical parallelization (PP) of variational atmospheric inversions. In the PP, the inversion still manages a single physically and statistically consistent window, but the transport model is run in parallel overlapping sub-segments in order to massively reduce the computation wall-clock time of the inversion. For global inversions, a simplification of transport modelling is described to connect the output of all segments. We demonstrate the performance of the approach on a global inversion for CO 2 with a 32 yr inversion window (1979–2010) with atmospheric measurements from 81 sites of the NOAA global cooperative air sampling network. In this case, we show that the duration of the inversion is reduced by a seven-fold factor (from months to days), while still processing the three decades consistently and with improved numerical stability.

Description: Spatial debris-cover effect on the maritime glaciers of Mount Gongga, south-eastern Tibetan Plateau The Cryosphere Discussions, 7, 2413-2453, 2013 Author(s): Y. Zhang, Y. Hirabayashi, K. Fujita, S. Liu, and Q. Liu The Tibetan Plateau and surroundings contain a large number of debris-covered glaciers, on which debris cover affects glacier response to climate change by altering ice melting rates and spatial patterns of mass loss. Insufficient spatial distribution of debris thickness data makes it difficult to analyze regional debris-cover effects. Mount Gongga glaciers, maritime glaciers in the south-eastern Tibetan Plateau, are characterized by a substantial reduction in glacier length and ice mass in recent decades. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)-derived thermal property of the debris layer reveals that 68% of the glaciers have extensive mantles of supraglacial debris in their ablation zones, in which the proportion of debris cover to total glacier area varies from 1.74% to 53.0%. Using a surface energy-mass balance model accounting for the debris-cover effect applied at a regional scale, we find that although the presence of supraglacial debris has a significant insulating effect on heavily debris-covered glaciers, it accelerates ice melting on ~ 10.2% of the total ablation area and produces rapid wastage of ~ 25% of the debris-covered glaciers, resulting in the similar mass losses between debris-covered and debris-free glaciers. Widespread debris cover also facilitates the development of active terminus regions. Regional differences in the debris-cover effect are apparent, highlighting the importance of debris cover for understanding glacier status and hydrology in both the Tibetan Plateau and other mountain ranges around the world.

Description: Decadal changes from a multi-temporal glacier inventory of Svalbard The Cryosphere Discussions, 7, 2489-2532, 2013 Author(s): C. Nuth, J. Kohler, M. König, A. von Deschwanden, J. O. Hagen, A. Kääb, G. Moholdt, and R. Pettersson We present a multi-temporal digital inventory of Svalbard glaciers with the most recent from the late 2000s containing 33 775 km 2 of glaciers, or 57% of the total land area of the archipelago. At present, 68% of the glaciated area of Svalbard drains through tidewater glaciers that have a summed terminus width of ~ 740 km. The glaciated area over the entire archipelago has decreased by an average of 80 km 2 a −1 over the past ~ 30 yr, representing a reduction of 7%. For a sample of ~ 400 glaciers (10 000 km 2 ) in the south and west of Spitsbergen, three digital inventories are available from 1930/60s, 1990 and 2007 from which we calculate average changes during 2 epochs. In the more recent epoch, the terminus retreat was larger than in the earlier epoch while area shrinkage was smaller. The contrasting pattern may be explained by the decreased lateral wastage of the glacier tongues. Temporal retreat rates for individual glaciers show a mix of accelerating and decelerating trends, reflecting the large spatial variability of glacier types and climatic/dynamic response times in Svalbard. Last, retreat rates estimated by dividing glacier area changes by the tongue width are larger than centerline retreat due to a more encompassing frontal change estimate with inclusion of lateral area loss.

Description: Inclusion of ash and SO 2 emissions from volcanic eruptions in WRF-Chem: development and some applications Geoscientific Model Development, 6, 457-468, 2013 Author(s): M. Stuefer, S. R. Freitas, G. Grell, P. Webley, S. Peckham, S. A. McKeen, and S. D. Egan We describe a new functionality within the Weather Research and Forecasting (WRF) model with coupled Chemistry (WRF-Chem) that allows simulating emission, transport, dispersion, transformation and sedimentation of pollutants released during volcanic activities. Emissions from both an explosive eruption case and a relatively calm degassing situation are considered using the most recent volcanic emission databases. A preprocessor tool provides emission fields and additional information needed to establish the initial three-dimensional cloud umbrella/vertical distribution within the transport model grid, as well as the timing and duration of an eruption. From this source condition, the transport, dispersion and sedimentation of the ash cloud can be realistically simulated by WRF-Chem using its own dynamics and physical parameterization as well as data assimilation. Examples of model applications include a comparison of tephra fall deposits from the 1989 eruption of Mount Redoubt (Alaska) and the dispersion of ash from the 2010 Eyjafjallajökull eruption in Iceland. Both model applications show good coincidence between WRF-Chem and observations.

Description: PORT, a CESM tool for the diagnosis of radiative forcing Geoscientific Model Development, 6, 469-476, 2013 Author(s): A. J. Conley, J.-F. Lamarque, F. Vitt, W. D. Collins, and J. Kiehl The Parallel Offline Radiative Transfer (PORT) model is a stand-alone tool, driven by model-generated datasets, that can be used for any radiation calculation that the underlying radiative transfer schemes can perform, such as diagnosing radiative forcing. In its present distribution, PORT isolates the radiation code from the Community Atmosphere Model (CAM4) in the Community Earth System Model (CESM1). The current configuration focuses on CAM4 radiation with the constituents as represented in present-day conditions in CESM1, along with their optical properties. PORT includes an implementation of stratospheric temperature adjustment under the assumption of fixed dynamical heating, which is necessary to compute radiative forcing in addition to the more straightforward instantaneous radiative forcing. PORT can be extended to use radiative constituent distributions from other models or model simulations. Ultimately, PORT can be used with various radiative transfer models. As illustrations of the use of PORT, we perform the computation of radiative forcing from doubling of carbon dioxide, from the change of tropospheric ozone concentration from the year 1850 to 2000, and from present-day aerosols. The radiative forcing from tropospheric ozone (with respect to 1850) generated by a collection of model simulations under the Atmospheric Chemistry and Climate Model Intercomparison Project is found to be 0.34 (with an intermodel standard deviation of 0.07) W m −2 . Present-day aerosol direct forcing (relative to no aerosols) is found to be −1.3 W m −2 .

Description: The effect of black carbon on reflectance of snow in the accumulation area of glaciers in the Baspa basin, Himachal Pradesh, India The Cryosphere Discussions, 7, 1359-1382, 2013 Author(s): A. V. Kulkarni, G. Vinay Kumar, H. S. Negi, J. Srinivasan, and S. K. Satheesh Himalayan glaciers are being extensively debated in scientific and public forums, as changes in their distribution can significantly affect the availability of water in many rivers originating in the region. The distribution of glaciers can be influenced by mass balance, and most of the glaciers located in the Pir Panjal and Greater Himalayan mountain ranges are losing mass at the rate of almost a meter per year. The Equilibrium Line Altitude (ELA) has also shifted upward by 400 m in the last two decades. This upward migration of ELA and the loss in mass could have been influenced by changes in temperature, precipitation and by the deposition of black carbon in the accumulation area of glaciers. The deposition of black carbon can reduce the albedo of snow in the accumulation area leading to faster melting of snow and causing more negative mass balance. In this investigation, a change in reflectance in the accumulation area of the Baspa basin is analysed for the year 2009, as the region has experienced extensive forest fires along with northern Indian biomass burning. The investigation has shown that: (1) The number of forest fires in the summer of 2009 was substantially higher than in any other year between 2001 and 2010; (2) the drop in reflectance in the visible region from April to May in the accumulation area was significantly higher in the year 2009 than in any other year from 2000 to 2012; (3) the temperature of the region was substantially lower than the freezing point during the active fire period of 2009, indicating the small influence of liquid water and grain size; (4) the drop in reflectance was observed only in the visible region, indicating role of contamination; (5) in the visible region, a mean drop in reflectance of 21± 5% was observed during the active fire period in the accumulation area. At some places, the drop was as high as 50 ± 5%. This can only be explained by the deposition of black carbon. The study suggests that a change in snow albedo in the accumulation area due to the deposition of black carbon from anthropogenic and natural causes can influence the mass balance of the glaciers in the Baspa basin, Himachal Pradesh, India.

Description: Ice-shelf buttressing and the stability of marine ice sheets The Cryosphere, 7, 647-655, 2013 Author(s): G. H. Gudmundsson Ice-shelf buttressing and the stability of marine-type ice sheets are investigated numerically. Buttressing effects are analysed for a situation where a stable grounding line is located on a bed sloping upwards in the direction of flow. Such grounding-line positions are known to be unconditionally unstable in the absence of transverse flow variations. It is shown that ice-shelf buttressing can restore stability under these conditions. Ice flux at the grounding line is, in general, not a monotonically increasing function of ice thickness. This, possibly at first somewhat counterintuitive result, is found to be fully consistent with recent theoretical work. Grounding lines on retrograde slopes are conditionally stable, and the stability regime is a non-trivial function of bed and ice-shelf geometry. The stability of grounding lines cannot be assessed from considerations of local bed slope only.

Description: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data The Cryosphere, 7, 615-630, 2013 Author(s): M. Tedesco, X. Fettweis, T. Mote, J. Wahr, P. Alexander, J. E. Box, and B. Wouters A combined analysis of remote sensing observations, regional climate model (RCM) outputs and reanalysis data over the Greenland ice sheet provides evidence that multiple records were set during summer 2012. Melt extent was the largest in the satellite era (extending up to ∼97% of the ice sheet) and melting lasted up to ∼2 months longer than the 1979–2011 mean. Model results indicate that near surface temperature was ∼3 standard deviations (σ) above the 1958–2011 mean, while surface mass balance (SMB) was ∼3σ below the mean and runoff was 3.9σ above the mean over the same period. Albedo, exposure of bare ice and surface mass balance also set new records, as did the total mass balance with summer and annual mass changes of, respectively, −627 Gt and −574 Gt, 2σ below the 2003–2012 mean. We identify persistent anticyclonic conditions over Greenland associated with anomalies in the North Atlantic Oscillation (NAO), changes in surface conditions (e.g., albedo, surface temperature) and preconditioning of surface properties from recent extreme melting as major driving mechanisms for the 2012 records. Less positive if not increasingly negative SMB will likely occur should these characteristics persist.

Description: Balanced conditions or slight mass gain of glaciers in the Lahaul and Spiti region (northern India, Himalaya) during the nineties preceded recent mass loss The Cryosphere, 7, 569-582, 2013 Author(s): C. Vincent, Al. Ramanathan, P. Wagnon, D. P. Dobhal, A. Linda, E. Berthier, P. Sharma, Y. Arnaud, M. F. Azam, P. G. Jose, and J. Gardelle The volume change of the Chhota Shigri Glacier (India, 32° 20 N, 77° 30' E) between 1988 and 2010 has been determined using in situ geodetic measurements. This glacier has experienced only a slight mass loss between 1988 and 2010 (−3.8 ± 2.0 m w.e. (water equivalent) corresponding to −0.17 ± 0.09 m w.e. yr −1 ). Using satellite digital elevation models (DEM) differencing and field measurements, we measure a negative mass balance (MB) between 1999 and 2010 (−4.8 ± 1.8 m w.e. corresponding to −0.44 ± 0.16 m w.e. yr −1 ). Thus, we deduce a slightly positive or near-zero MB between 1988 and 1999 (+1.0 ± 2.7 m w.e. corresponding to +0.09 ± 0.24 m w.e. yr −1 ). Furthermore, satellite DEM differencing reveals that the MB of the Chhota Shigri Glacier (−0.39 ± 0.15 m w.e. yr −1 ) has been only slightly less negative than the MB of a 2110 km 2 glaciarized area in the Lahaul and Spiti region (−0.44 ± 0.09 m w.e. yr −1 ) during 1999−2011. Hence, we conclude that the ice wastage is probably moderate in this region over the last 22 yr, with near equilibrium conditions during the nineties, and an ice mass loss after. The turning point from balanced to negative mass budget is not known but lies probably in the late nineties and at the latest in 1999. This positive or near-zero MB for Chhota Shigri Glacier (and probably for the surrounding glaciers of the Lahaul and Spiti region) during at least part of the 1990s contrasts with a recent compilation of MB data in the Himalayan range that indicated ice wastage since 1975. However, in agreement with this compilation, we confirm more negative balances since the beginning of the 21st century.

Description: Supercooled interfacial water in fine grained soils probed by dielectric spectroscopy The Cryosphere Discussions, 7, 1441-1493, 2013 Author(s): A. Lorek and N. Wagner Water as thermodynamic state parameter affects nearly all physical, chemical and biological processes on the earth. Recent Mars observations as well as laboratory investigations suggest that water is also a key factor of current physical and chemical processes on the martian surface, e.g. rheological phenomena. Therefore it is of particular interest to get information about the liquid like state of water on martian analog soils in the temperature range below 0 °C. In this context, a parallel plate capacitor has been developed to obtain isothermal dielectric spectra of fine grained soils in the frequency range from 10 Hz to 1.1 MHz at martian like temperatures down to −70 °C. Two martian analogue soils have been investigated: a Ca-Bentonite (specific surface of 237 m 2 g −1 , up to 9.4% w/w gravimetric water content) and JSC Mars 1, a volcanic ash (specific surface of 146 m 2 g −1 , up to 7.4% w/w ). Three soil-specific relaxation processes are observed in the investigated frequency-temperature range: two weak high frequency processes (bound or hydrated water as well as ice) and a strong low frequency process due to counter ion relaxation and the Maxwell–Wagner effect. To characterize the dielectric relaxation behavior, a generalized fractional dielectric relaxation model is applied assuming three active relaxation processes with relaxation time of the i th process according to an Eyring equation. The real part of effective complex soil permittivity at 350 kHz was used to determine ice and liquid like water content by means of the Birchak or CRIM equation. There are evidence that Bentonite down to −70 °C has a liquid like water content of 1.17 monolayers and JSC Mars 1 a liquid like water content of 1.96 mono layers.

Description: The geomorphological effect of cornice fall avalanches in the Longyeardalen valley, Svalbard The Cryosphere, 7, 1361-1374, 2013 Author(s): M. Eckerstorfer, H. H. Christiansen, L. Rubensdotter, and S. Vogel The study of snow avalanches and their geomorphological effect in the periglacial parts of the cryosphere is important for enhanced geomorphological process understanding as well as hazard-related studies. Only a few field studies, and particularly few in the High Arctic, have quantified avalanche sedimentation. Snow avalanches are traditionally ranked behind rockfall in terms of their significance for mass-wasting processes of rockslopes. Cornice fall avalanches are at present the most dominant snow avalanche type at two slope systems, called Nybyen and Larsbreen, in the valley Longyeardalen in central Svalbard. Both slope systems are on northwest-facing lee slopes underneath a large summit plateau, with annual cornices forming on the top. High-frequency and magnitude cornice fall avalanching is observed by daily automatic time-lapse photography. In addition, rock debris sedimentation by cornice fall avalanches was measured directly in permanent sediment traps or by snow inventories. The results from a maximum of seven years of measurements in a total of 13 catchments show maximum mean rock debris sedimentation rates ranging from 8.2 to 38.7 kg m −2 at Nybyen, and from 0.8 to 55.4 kg m −2 at Larsbreen. Correspondingly, avalanche fan surfaces accreted from 2.6 to 8.8 mm yr −1 at Nybyen, and from 0.2 to 13.9 mm yr −1 at Larsbreen. This comparably efficient rockslope mass wasting is due to collapsing cornices producing cornice fall avalanches containing large amounts of rock debris throughout the entire winter. The rock debris of different origin stems from the plateau crests, the adjacent free rock face and the transport pathway, accumulating distinct avalanche fans at both slope systems. Cornice fall avalanche sedimentation also contributed to the development of a rock glacier at the Larsbreen site during the Holocene. We have recorded present maximum rockwall retreat rates of 0.9 mm yr −1 at Nybyen, but as much as 6.7 mm yr −1 at Larsbreen, while average Holocene rockwall retreat rates of 1.1 mm yr −1 at Nybyen have been determined earlier. As cornice fall avalanches are the dominant type of avalanche in central Svalbard, the related geomorphological effect is assumed to be of significance at periglacial landscape scale. A climate-induced shift in prevailing winter wind direction could change the rockslope sedimentation effectively by changing the snow avalanche activity.

Description: Downscaling a global climate model to simulate climate change over the US and the implication on regional and urban air quality Geoscientific Model Development, 6, 1429-1445, 2013 Author(s): M. Trail, A. P. Tsimpidi, P. Liu, K. Tsigaridis, Y. Hu, A. Nenes, and A. G. Russell Climate change can exacerbate future regional air pollution events by making conditions more favorable to form high levels of ozone. In this study, we use spectral nudging with the Weather Research and Forecasting (WRF) model to downscale NASA earth system GISS modelE2 results during the years 2006 to 2010 and 2048 to 2052 over the contiguous United States in order to compare the resulting meteorological fields from the air quality perspective during the four seasons of five-year historic and future climatological periods. GISS results are used as initial and boundary conditions by the WRF regional climate model (RCM) to produce hourly meteorological fields. The downscaling technique and choice of physics parameterizations used are evaluated by comparing them with in situ observations. This study investigates changes of similar regional climate conditions down to a 12 km by 12 km resolution, as well as the effect of evolving climate conditions on the air quality at major US cities. The high-resolution simulations produce somewhat different results than the coarse-resolution simulations in some regions. Also, through the analysis of the meteorological variables that most strongly influence air quality, we find consistent changes in regional climate that would enhance ozone levels in four regions of the US during fall (western US, Texas, northeastern, and southeastern US), one region during summer (Texas), and one region where changes potentially would lead to better air quality during spring (Northeast). Changes in regional climate that would enhance ozone levels are increased temperatures and stagnation along with decreased precipitation and ventilation. We also find that daily peak temperatures tend to increase in most major cities in the US, which would increase the risk of health problems associated with heat stress. Future work will address a more comprehensive assessment of emissions and chemistry involved in the formation and removal of air pollutants.

Description: The SOCOL version 3.0 chemistry–climate model: description, evaluation, and implications from an advanced transport algorithm Geoscientific Model Development, 6, 1407-1427, 2013 Author(s): A. Stenke, M. Schraner, E. Rozanov, T. Egorova, B. Luo, and T. Peter We present the third generation of the coupled chemistry–climate model (CCM) SOCOL (modeling tools for studies of SOlar Climate Ozone Links). The most notable modifications compared to the previous model version are (1) the dynamical core has been updated with the fifth generation of the middle-atmosphere general circulation model MA-ECHAM (European Centre/HAMburg climate model), and (2) the advection of the chemical species is now calculated by a mass-conserving and shape-preserving flux-form transport scheme instead of the previously used hybrid advection scheme. The whole chemistry code has been rewritten according to the ECHAM5 infrastructure and transferred to Fortran95. In contrast to its predecessors, SOCOLvs3 is now fully parallelized. The performance of the new SOCOL version is evaluated on the basis of transient model simulations (1975–2004) with different horizontal (T31 and T42) resolutions, following the approach of the CCMVal-1 model validation activity. The advanced advection scheme significantly reduces the artificial loss and accumulation of tracer mass in regions with strong gradients that was observed in previous model versions. Compared to its predecessors, SOCOLvs3 generally shows more realistic distributions of chemical trace species, especially of total inorganic chlorine, in terms of the mean state, but also of the annual and interannual variability. Advancements with respect to model dynamics are for example a better representation of the stratospheric mean state in spring, especially in the Southern Hemisphere, and a slowdown of the upward propagation in the tropical lower stratosphere. Despite a large number of improvements model deficiencies still remain. Examples include a too-fast vertical ascent and/or horizontal mixing in the tropical stratosphere, the cold temperature bias in the lowermost polar stratosphere, and the overestimation of polar total ozone loss during Antarctic springtime.

Description: Stable water isotopes in the coupled atmosphere–land surface model ECHAM5-JSBACH Geoscientific Model Development, 6, 1463-1480, 2013 Author(s): B. Haese, M. Werner, and G. Lohmann In this study we present first results of a new model development, ECHAM5-JSBACH-wiso, where we have incorporated the stable water isotopes H 2 18 O and HDO as tracers in the hydrological cycle of the coupled atmosphere–land surface model ECHAM5-JSBACH. The ECHAM5-JSBACH-wiso model was run under present-day climate conditions at two different resolutions (T31L19, T63L31). A comparison between ECHAM5-JSBACH-wiso and ECHAM5-wiso shows that the coupling has a strong impact on the simulated temperature and soil wetness. Caused by these changes of temperature and the hydrological cycle, the δ 18 O in precipitation also shows variations from −4‰ up to 4‰. One of the strongest anomalies is shown over northeast Asia where, due to an increase of temperature, the δ 18 O in precipitation increases as well. In order to analyze the sensitivity of the fractionation processes over land, we compare a set of simulations with various implementations of these processes over the land surface. The simulations allow us to distinguish between no fractionation, fractionation included in the evaporation flux (from bare soil) and also fractionation included in both evaporation and transpiration (from water transport through plants) fluxes. While the isotopic composition of the soil water may change for δ 18 O by up to +8&permil:, the simulated δ 18 O in precipitation shows only slight differences on the order of ±1‰. The simulated isotopic composition of precipitation fits well with the available observations from the GNIP (Global Network of Isotopes in Precipitation) database.

Description: δ 18 O water isotope in the i LOVECLIM model (version 1.0) – Part 2: Evaluation of model results against observed δ 18 O in water samples Geoscientific Model Development, 6, 1493-1504, 2013 Author(s): D. M. Roche and T. Caley The H 2 18 O stable isotope was previously introduced in the three coupled components of the earth system model i LOVECLIM: atmosphere, ocean and vegetation. The results of a long (5000 yr) pre-industrial equilibrium simulation are presented and evaluated against measurement of H 2 18 O abundance in present-day water for the atmospheric and oceanic components. For the atmosphere, it is found that the model reproduces the observed spatial distribution and relationships to climate variables with some merit, though limitations following our approach are highlighted. Indeed, we obtain the main gradients with a robust representation of the Rayleigh distillation but caveats appear in Antarctica and around the Mediterranean region due to model limitation. For the oceanic component, the agreement between the modelled and observed distribution of water δ 18 O is found to be very good. Mean ocean surface latitudinal gradients are faithfully reproduced as well as the mark of the main intermediate and deep water masses. This opens large prospects for the applications in palaeoclimatic context.

Description: The Rock Geochemical Model (RokGeM) v0.9 Geoscientific Model Development, 6, 1543-1573, 2013 Author(s): G. Colbourn, A. Ridgwell, and T. M. Lenton A new model of terrestrial rock weathering – the Rock Geochemical Model (RokGeM) – was developed for incorporation into the GENIE Earth System modelling framework. In this paper we describe the model. We consider a range of previously devised parameterizations, ranging from simple dependencies on global mean temperature following Berner et al. (1983), to spatially explicit dependencies on run-off and temperature (GKWM, Bluth and Kump, 1994; GEM-CO2, Amiotte-Suchet et al., 2003) – fields provided by the energy-moisture balance atmosphere model component in GENIE. Using long-term carbon cycle perturbation experiments, we test the effects of a wide range of model parameters, including whether or not the atmosphere was "short-circuited" in the carbon cycle; the sensitivity and feedback strength of temperature and run-off on carbonate and silicate weathering; different river-routing schemes; 0-D (global average) vs. 2-D (spatially explicit) weathering schemes; and the lithology dependence of weathering. Included are details of how to run the model and visualize the results.

Description: Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram The Cryosphere, 7, 1385-1398, 2013 Author(s): R. Bhambri, T. Bolch, P. Kawishwar, D. P. Dobhal, D. Srivastava, and B. Pratap Glaciers in the Karakoram show long-term irregular behaviour with comparatively frequent and sudden advances. A glacier inventory of the upper Shyok valley situated in northeast Karakoram has been generated for the year 2002 using Landsat ETM+ and SRTM3 DEM as baseline data for the investigations and subsequent change analysis. The upper Shyok valley contained 2123 glaciers (larger than 0.02 km 2 in size) with an area of 2977.9 ± 95.3 km 2 in 2002. Out of these, 18 glaciers with an area of 1004.1 ± 32.1 km 2 showed surge-type behaviour. Change analysis based on Hexagon KH-9 (years 1973 and 1974) and Landsat TM/ETM+ (years 1989, 2002 and 2011) images had to be restricted to a subset of 136 glaciers (covering an area of 1609.7 ± 51.5 km 2 in 2002) due to adverse snow conditions. The area of the investigated glaciers, including the 18 surge-type glaciers identified, showed no significant changes during all studied periods. However, the analysis provides a hint that the overall glacier area slightly decreased until about 1989 (area 1973: 1613.6 ± 43.6 km 2 ; area 1989: 1602.0 ± 33.6 km 2 ) followed by an increase (area 2002: 1609.7 ± 51.5; area 2011: 1615.8 ± 35.5 km 2 ). Although the overall change in area is insignificant, advances in glacier tongues since the end of the 1980s are clearly visible. Detailed estimations of length changes for individual glaciers since the 1970s and for Central Rimo Glacier since the 1930s confirm the irregular retreat and advance.

Description: Influence of snow depth distribution on surface roughness in alpine terrain: a multi-scale approach The Cryosphere Discussions, 7, 4633-4680, 2013 Author(s): J. Veitinger, B. Sovilla, and R. S. Purves In alpine terrain, the snow covered winter surface deviates from its underlying summer terrain due to the progressive smoothing caused by snow accumulation. Terrain smoothing is believed to be an important factor in avalanche formation, avalanche dynamics and affects surface heat transfer, energy balance as well as snow depth distribution. To characterize the effect of snow on terrain we use the concept of roughness. Roughness is calculated for several snow surfaces and its corresponding underlying terrain for three alpine basins in the Swiss Alps characterized by low medium and high terrain roughness. To this end, elevation models of winter and summer terrain are derived from high-resolution (1 m) measurements performed by airborne and terrestrial LIDAR. We showed that on basin scale terrain smoothing not only depends on mean snow depth in the basin but also on its variability. Terrain smoothing can be modelled in function of mean snow depth and its standard deviation using a power law. However, a relationship between terrain smoothing and snow depth does not exist on a pixel scale. Further we demonstrated the high persistence of snow surface roughness even in between winter seasons. Those persistent patterns might be very useful to improve the representation of a winter terrain without modelling of the snow cover distribution. This can potentially improve avalanche release area definition and in the long term natural hazard management strategies.

Description: An approach to computing direction relations between separated object groups Geoscientific Model Development, 6, 1591-1599, 2013 Author(s): H. Yan, Z. Wang, and J. Li Direction relations between object groups play an important role in qualitative spatial reasoning, spatial computation and spatial recognition. However, none of existing models can be used to compute direction relations between object groups. To fill this gap, an approach to computing direction relations between separated object groups is proposed in this paper, which is theoretically based on gestalt principles and the idea of multi-directions. The approach firstly triangulates the two object groups, and then it constructs the Voronoi diagram between the two groups using the triangular network. After this, the normal of each Voronoi edge is calculated, and the quantitative expression of the direction relations is constructed. Finally, the quantitative direction relations are transformed into qualitative ones. The psychological experiments show that the proposed approach can obtain direction relations both between two single objects and between two object groups, and the results are correct from the point of view of spatial cognition.

Description: Weak layer fracture: facets and depth hoar The Cryosphere, 7, 1447-1453, 2013 Author(s): I. Reiweger and J. Schweizer Understanding failure initiation within weak snow layers is essential for modeling and predicting dry-snow slab avalanches. We therefore performed laboratory experiments with snow samples containing a weak layer consisting of either faceted crystals or depth hoar. During these experiments the samples were loaded with different loading rates and at various tilt angles until fracture. The strength of the samples decreased with increasing loading rate and increasing tilt angle. Additionally, we took pictures of the side of four samples with a high-speed video camera and calculated the displacement using a particle image velocimetry (PIV) algorithm. The fracture process within the weak layer could thus be observed in detail. Catastrophic failure started due to a shear fracture just above the interface between the depth hoar layer and the underlying crust.

Description: Scatter of mass changes estimates at basin scale for Greenland and Antarctica The Cryosphere, 7, 1411-1432, 2013 Author(s): V. R. Barletta, L. S. Sørensen, and R. Forsberg During the last decade, the GRACE mission has provided valuable data for determining the mass changes of the Greenland and Antarctic ice sheets. Yet, discrepancies still exist in the published mass balance results, and comprehensive analyses on the sources of errors and discrepancies are lacking. Here, we present monthly mass changes together with trends derived from GRACE data at basin scale for both the Greenland and Antarctic ice sheets, and we assess the variability and errors for each of the possible sources of discrepancies, and we do this in an unprecedented systematic way, taking into account mass inference methods, data sets and background models. We find a very good agreement between the monthly mass change results derived from two independent methods, which represents a cross validation. For the monthly solutions, we find that most of the scatter is caused by the use of the two different data sets rather than the two different methods applied. Besides the well-known GIA trend uncertainty, we find that the geocenter motion and the recent de-aliasing corrections significantly impact the trends, with contributions of +13.2 Gt yr −1 and −20 Gt yr −1 , respectively, for Antarctica, which is more affected by these than Greenland. We show differences between the use of release RL04 and the new RL05 and confirm a lower noise content in the new release. The overall scatter of the solutions well exceeds the uncertainties propagated from the data errors and the leakage (as done in the past); hence we calculate new sound total errors for the monthly solutions and the trends. We find that the scatter in the monthly solutions caused by applying different estimates of geocenter motion time series (degree-1 corrections) is significant – contributing with up to 40% of the total error. For the whole GRACE period (2003–2011) our trend estimate for Greenland is −234 ± 20 Gt yr −1 and −83 ± 36 Gt yr −1 for Antarctica (−111 ± 15 Gt yr −1 in the western part). We also find a clear (with respect to our errors) increase of mass loss in the last four years.

Description: A method to represent ozone response to large changes in precursor emissions using high-order sensitivity analysis in photochemical models Geoscientific Model Development, 6, 1601-1608, 2013 Author(s): G. Yarwood, C. Emery, J. Jung, U. Nopmongcol, and T. Sakulyanontvittaya Photochemical grid models (PGMs) are used to simulate tropospheric ozone and quantify its response to emission changes. PGMs are often applied for annual simulations to provide both maximum concentrations for assessing compliance with air quality standards and frequency distributions for assessing human exposure. Efficient methods for computing ozone at different emission levels can improve the quality of ozone air quality management efforts. This study demonstrates the feasibility of using the decoupled direct method (DDM) to calculate first- and second-order sensitivity of ozone to anthropogenic NO x and VOC emissions in annual PGM simulations at continental scale. Algebraic models are developed that use Taylor series to produce complete annual frequency distributions of hourly ozone at any location and any anthropogenic emission level between zero and 100%, adjusted independently for NO x and VOC. We recommend computing the sensitivity coefficients at the midpoint of the emissions range over which they are intended to be applied, in this case with 50% anthropogenic emissions. The algebraic model predictions can be improved by combining sensitivity coefficients computed at 10 and 50% anthropogenic emissions. Compared to brute force simulations, algebraic model predictions tend to be more accurate in summer than winter, at rural than urban locations, and with 100% than zero anthropogenic emissions. Equations developed to combine sensitivity coefficients computed with 10 and 50% anthropogenic emissions are able to reproduce brute force simulation results with zero and 100% anthropogenic emissions with a mean bias of less than 2 ppb and mean error of less than 3 ppb averaged over 22 US cities.

Description: A Lagrangian advection scheme with shape matrix (LASM) for solving advection problems Geoscientific Model Development, 7, 2951-2968, 2014 Author(s): L. Dong, B. Wang, and L. Liu A new Lagrangian advection scheme with shape matrix (LASM) is proposed to take advantage of the extreme low numerical diffusion of the Lagrangian methods. The tracer is discretized into finite parcels, which move along the downstream trajectories. Different from other Lagrangian schemes, the parcel shape is simulated explicitly by a linear transformation matrix. By doing so, the aliasing error in the Lagrangian schemes is largely reduced without introducing substantial interparcel mixing in the pure advection stage, because the flow information will be respected when remapping tracer density onto the fixed model grids. An adaptive interparcel mixing algorithm is constructed to ensure the validity of the linear approximation of the parcel shape, where the mixing is only triggered when it is necessary and resembles the physical mixing. The total tracer mass on the parcels is conserved exactly. The new scheme is validated by using several test cases.

Description: Air temperature variability over three glaciers in the Ortles-Cevedale (Italian Alps): effects of glacier disintegration, intercomparison of calculation methods, and impacts on mass balance modeling The Cryosphere Discussions, 8, 6147-6192, 2014 Author(s): L. Carturan, F. Cazorzi, F. De Blasi, and G. Dalla Fontana Glacier mass balance models rely on accurate spatial calculation of input data, in particular air temperature. Lower temperatures (the so-called glacier cooling effect), and lower temperature variability (the so-called glacier damping effect) generally occur over glaciers, compared to ambient conditions. These effects, which depend on the geometric characteristics of glaciers and display a high spatial and temporal variability, have been mostly investigated on medium- to large-size glaciers so far, while observations on smaller ice bodies are scarce. Using a dataset from 8 on-glacier and 4 off-glacier weather stations, collected in summer 2010 and 2011, we analyzed the air temperature variability and wind regime over three different glaciers in the Ortles-Cevedale. The magnitude of the cooling effect and the occurrence of katabatic boundary layer (KBL) processes showed remarkable differences among the three ice bodies, suggesting the likely existence of important reinforcing mechanisms during glacier decay and disintegration. None of the methods proposed in the literature for calculating on-glacier temperature from off-glacier data fully reproduced our observations. Among them, the more physically-based procedure of Greuell and Böhm (1998) provided the best overall results where the KBL prevail, but it was not effective elsewhere (i.e. on smaller ice bodies and close to the glacier margins). The accuracy of air temperature estimations strongly impacted the results from a mass balance model which was applied to the three investigated glaciers. Most importantly, even small temperature deviations caused distortions in parameter calibration, thus compromising the model generalizability.

Description: Sensitivity of the Mediterranean sea level to atmospheric pressure and free surface elevation numerical formulation in NEMO Geoscientific Model Development, 7, 3001-3015, 2014 Author(s): P. Oddo, A. Bonaduce, N. Pinardi, and A. Guarnieri The sensitivity of the dynamics of the Mediterranean Sea to atmospheric pressure and free surface elevation formulation using NEMO (Nucleus for European Modelling of the Ocean) was evaluated. Four different experiments were carried out in the Mediterranean Sea using filtered or explicit free surface numerical schemes and accounting for the effect of atmospheric pressure in addition to wind and buoyancy fluxes. Model results were evaluated by coherency and power spectrum analysis with tide gauge data. We found that atmospheric pressure plays an important role for periods shorter than 100 days. The free surface formulation is important to obtain the correct ocean response for periods shorter than 30 days. At frequencies higher than 15 days −1 the Mediterranean basin's response to atmospheric pressure was not coherent and the performance of the model strongly depended on the specific area considered. A large-amplitude seasonal oscillation observed in the experiments using a filtered free surface was not evident in the corresponding explicit free surface formulation case, which was due to a phase shift between mass fluxes in the Gibraltar Strait and at the surface. The configuration with time splitting and atmospheric pressure always performed best; the differences were enhanced at very high frequencies.

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Strategy of valid 14 C dates choice in syngenetic permafrost The Cryosphere Discussions, 8, 5589-5621, 2014 Author(s): Y. K. Vasil'chuk and A. C. Vasil'chuk The main problem of radiocarbon dating within permafrost is the uncertain reliability of the 14 C dates. Syngenetic sediments contain allochthonous organic deposit that originated at a distance from its present position. Due to the very good preservation of organic materials in permafrost conditions and numerous re-burials of the fossils from ancient deposits into younger ones the dates could be both younger and older than the true age of dated material. The strategy for the most authentic radiocarbon date selection for dating of syncryogenic sediments is considered taking into account the fluvial origin of the syngenetic sediments. The re-deposition of organic material is discussed in terms of cyclic syncryogenic sedimentation and also the possible re-deposition of organic material in subaerial-subaqueous conditions. The advantages and the complications of dating organic micro-inclusions from ice wedges by the accelerator mass spectrometry (AMS) method are discussed applying to true age of dated material search. Radiocarbon dates of different organic materials from the same samples are compared. The younger age of the yedoma from cross-sections of Duvanny Yar in Kolyma River and Mamontova Khayata in the mouth of Lena River is substantiated due to the principle of the choice of the youngest 14 C date from the set.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Simulating the Antarctic ice sheet in the Late-Pliocene warm period: PLISMIP-ANT, an ice-sheet model intercomparison project The Cryosphere Discussions, 8, 5539-5588, 2014 Author(s): B. de Boer, A. M. Dolan, J. Bernales, E. Gasson, H. Goelzer, N. R. Golledge, J. Sutter, P. Huybrechts, G. Lohmann, I. Rogozhina, A. Abe-Ouchi, F. Saito, and R. S. W. van de Wal In the context of future climate change, understanding the nature and behaviour of ice sheets during warm intervals in Earth history is of fundamental importance. The Late-Pliocene warm period (also known as the PRISM interval: 3.264 to 3.025 million years before present) can serve as a potential analogue for projected future climates. Although Pliocene ice locations and extents are still poorly constrained, a significant contribution to sea-level rise should be expected from both the Greenland ice sheet and the West and East Antarctic ice sheets based on palaeo sea-level reconstructions. Here, we present results from simulations of the Antarctic ice sheet by means of an international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP-ANT). For the experiments, ice-sheet models including the shallow ice and shelf approximations have been used to simulate the complete Antarctic domain (including grounded and floating ice). We compare the performance of six existing numerical ice-sheet models in simulating modern control and Pliocene ice sheets by a suite of four sensitivity experiments. Ice-sheet model forcing fields are taken from the HadCM3 atmosphere–ocean climate model runs for the pre-industrial and the Pliocene. We include an overview of the different ice-sheet models used and how specific model configurations influence the resulting Pliocene Antarctic ice sheet. The six ice-sheet models simulate a comparable present-day ice sheet, although the models are setup with their own parameter settings. For the Pliocene simulations using the Bedmap1 bedrock topography, some models show a small retreat of the East Antarctic ice sheet, which is thought to have happened during the Pliocene for the Wilkes and Aurora basins. This can be ascribed to either the surface mass balance, as the HadCM3 Pliocene climate shows a significant increase over the Wilkes and Aurora basin, or the initial bedrock topography. For the latter, our simulations with the recently published Bedmap2 bedrock topography indicate a significantly larger contribution to Pliocene sea-level rise from the East Antarctic ice sheet for all six models relative to the simulations with Bedmap1.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Description: Snow mass decrease in the Northern Hemisphere (1979/80–2010/11) The Cryosphere Discussions, 8, 5623-5644, 2014 Author(s): Z. Li, J. Liu, L. Huang, N. Wang, B. Tian, J. Zhou, Q. Chen, and P. Zhang Snow cover has a key effect on climate change and hydrological cycling, as well as water supply to a sixth of the world's population across the Northern Hemisphere. However, reliable data on trends in snow cover in the Northern Hemisphere is lacking. Snow water equivalent (SWE) is a common measure of the amount of equivalent water of the snow pack. Here we verify the accuracy of three existing global SWE products and merge the most accurate aspects of them to generate a new SWE product covering the last 32 years (1979/80–2010/11). Using this new SWE product, we show that there has been a significant decreasing trend in the total mass of snow in the Northern Hemisphere. The most notable changes in total snow mass are −16.45 ± 6.68 and −13.55 ± 7.80 Gt year −1 in January and February, respectively. These are followed by March and December, which have trends of −12.58 ± 6.88 and −10.70 ± 5.62 Gt year −1 , respectively, from 1979/80 to 2010/11. During the same period, the temperature in the study area raised 0.17 °C decade −1 , which is thought to be the main reason of SWE decline.

Description: Factors controlling Slope Environmental Lapse Rate (SELR) of temperature in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya The Cryosphere Discussions, 8, 5645-5686, 2014 Author(s): R. J. Thayyen and A. P. Dimri Moisture, temperature and precipitation interplay forced through the orographic processes sustains the Himalayan cryospheric system. However, factors controlling the Slope Environmental Lapse Rate (SELR) of temperature along the higher Himalayan mountain slopes across various glacio-hydrologic regimes remain as a key knowledge gap. Present study dwells on the orographic processes driving the moisture–temperature interplay in the monsoon and cold-arid glacio-hydrological regimes of the Himalaya. Systematic data collection at three altitudes between 2540 and 3763 m a.s.l. in the Garhwal Himalaya (hereafter called monsoon regime) and between 3500 and 5600 m a.s.l. in the Ladakh Himalaya (herefater called cold-arid regime) revealed moistrue control on temperature distribution at temporal and spatial scales. Observed daily SELR of temperature ranges between 9.0 to 1.9 °C km −1 and 17.0 to 2.8 °C km −1 in the monsoon and cold-arid regimes respectively highlighting strong regional variability. Moisture influx to the region, either from Indian summer monsoon (ISM) or from Indian winter monsoon (IWM) forced lowering of SELR. This phenophena of "monsoon lowering" of SELR is due to the release latent heat of condensation from orographically focred lifted air parcel. Seasonal response of SELR in the monsoon regime is found to be closly linked with the variations in the local lifting condensation levels (LCL). Contrary to this, cold-arid system is characterised by the extremely high values of daily SELR upto 17 °C km −1 signifying the extremely arid conditions prevailing in summer. Distinctly lower SELR devoid of monsoon lowering at higher altitude sections of monsoon and cold-arid regimes suggests sustained wetter high altitude regimes. We have proposed a SELR model for both glacio-hydrological regimes demostrating with two sections each using a derivative of the Clausius–Clapeyron relationship by deriving monthly SELR indices. It has been proposed that the manifestations of presence or absence of moisture is the single most important factor determining the temperature distribution along the higher Himalayan slopes driven by the orographic forcings. This work also suggests that the arbitary use of temperature lapse rate to extrapolate temperature to the higher Himalaya is extremely untenable.

Description: Gaseous chemistry and aerosol mechanism developments for version 3.5.1 of the online regional model, WRF-Chem Geoscientific Model Development, 7, 2557-2579, 2014 Author(s): S. Archer-Nicholls, D. Lowe, S. Utembe, J. Allan, R. A. Zaveri, J. D. Fast, Ø. Hodnebrog, H. Denier van der Gon, and G. McFiggans We have made a number of developments to the Weather, Research and Forecasting model coupled with Chemistry (WRF-Chem), with the aim of improving model prediction of trace atmospheric gas-phase chemical and aerosol composition, and of interactions between air quality and weather. A reduced form of the Common Reactive Intermediates gas-phase chemical mechanism (CRIv2-R5) has been added, using the Kinetic Pre-Processor (KPP) interface, to enable more explicit simulation of VOC degradation. N 2 O 5 heterogeneous chemistry has been added to the existing sectional MOSAIC aerosol module, and coupled to both the CRIv2-R5 and existing CBM-Z gas-phase schemes. Modifications have also been made to the sea-spray aerosol emission representation, allowing the inclusion of primary organic material in sea-spray aerosol. We have worked on the European domain, with a particular focus on making the model suitable for the study of nighttime chemistry and oxidation by the nitrate radical in the UK atmosphere. Driven by appropriate emissions, wind fields and chemical boundary conditions, implementation of the different developments are illustrated, using a modified version of WRF-Chem 3.4.1, in order to demonstrate the impact that these changes have in the Northwest European domain. These developments are publicly available in WRF-Chem from version 3.5.1 onwards.

Description: A model study of Abrahamsenbreen, a surging glacier in northern Spitsbergen The Cryosphere Discussions, 8, 5687-5726, 2014 Author(s): J. Oerlemans and W. J. J. van Pelt The climate sensitivity of Abrahamsenbreen, a 20 km long surge-type glacier in northern Spitsbergen, is studied with a simple glacier model. A scheme to describe the surges is included, which makes it possible to account for the effect of surges on the total mass budget of the glacier. A climate reconstruction back to AD 1300, based on ice-core data from Lomonosovfonna and climate records from Longyearbyen, is used to drive the model. The model is calibrated by requesting that it produces the correct Little Ice Age maximum glacier length and simulates the observed magnitude of the 1978-surge. Abrahamsenbreen is strongly out of balance with the current climate. If climatic conditions will remain as they were for the period 1989–2010, the glacier will ultimately shrink to a length of about 4 km (but this will take hundreds of years). For a climate change scenario involving a 2 m yr −1 rise of the equilibrium line from now onwards, we predict that in the year 2100 Abrahamsenbreen will be about 12 km long. The main effect of a surge is to lower the mean surface elevation and to increase the ablation area, thereby causing a negative perturbation of the mass budget. We found that the occurrence of surges leads to a somewhat stronger retreat of the glacier in a warming climate. Because of the very small bed slope, Abrahamsenbreen is sensitive to small perturbations in the equilibrium-line altitude E . For a decrease of E of only 160 m, the glacier would steadily grow into the Woodfjorddalen until after 2000 years it would reach the Woodfjord and calving could slow down the advance.

Description: A strategy for GIS-based 3-D slope stability modelling over large areas Geoscientific Model Development, 7, 2969-2982, 2014 Author(s): M. Mergili, I. Marchesini, M. Alvioli, M. Metz, B. Schneider-Muntau, M. Rossi, and F. Guzzetti GIS-based deterministic models may be used for landslide susceptibility mapping over large areas. However, such efforts require specific strategies to (i) keep computing time at an acceptable level, and (ii) parameterize the geotechnical data. We test and optimize the performance of the GIS-based, 3-D slope stability model r.slope.stability in terms of computing time and model results. The model was developed as a C- and Python-based raster module of the open source software GRASS GIS and considers the 3-D geometry of the sliding surface. It calculates the factor of safety (FoS) and the probability of slope failure ( P f ) for a number of randomly selected potential slip surfaces, ellipsoidal or truncated in shape. Model input consists of a digital elevation model (DEM), ranges of geotechnical parameter values derived from laboratory tests, and a range of possible soil depths estimated in the field. Probability density functions are exploited to assign P f to each ellipsoid. The model calculates for each pixel multiple values of FoS and P f corresponding to different sliding surfaces. The minimum value of FoS and the maximum value of P f for each pixel give an estimate of the landslide susceptibility in the study area. Optionally, r.slope.stability is able to split the study area into a defined number of tiles, allowing parallel processing of the model on the given area. Focusing on shallow landslides, we show how multi-core processing makes it possible to reduce computing times by a factor larger than 20 in the study area. We further demonstrate how the number of random slip surfaces and the sampling of parameters influence the average value of P f and the capacity of r.slope.stability to predict the observed patterns of shallow landslides in the 89.5 km 2 Collazzone area in Umbria, central Italy.

Description: Brief communication "Snow profile associated measurements (SPAM) – a new instrument for quick snow profile measurements" The Cryosphere Discussions, 5, 1737-1748, 2011 Author(s): P. Lahtinen A new instrument concept (SPAM) for snow profile associated measurements is presented. The potential of the concept is demonstrated by presenting preliminary results obtained with the prototype instrument. With this concept it is possible to retrieve rapid snow profiles of e.g. light extinction, reflectance, temperature and snow layer structure with high vertical resolution. As a side-product, also snow depth is retrieved.

Description: Modeling the impact of wintertime rain events on the thermal regime of permafrost The Cryosphere Discussions, 5, 1697-1736, 2011 Author(s): S. Westermann, J. Boike, M. Langer, T. V. Schuler, and B. Etzelmüller In this study, we present field measurements and numerical process modeling from Western Svalbard showing that the ground surface temperature below the snow is impacted by strong wintertime rain events. During such events, rain water percolates to the bottom of the snow pack, where it freezes and releases latent heat. In the winter season 2005/2006, on the order of 20 to 50 % of the wintertime precipitation fell as rain, thus confining the surface temperature to close to 0 °C for several weeks. The measured average ground surface temperature during the snow-covered period is −0.6 °C, despite of a snow surface temperature of on average −8.5 °C. For the considered period, the temperature threshold below which permafrost is sustainable on long timescales is exceeded. We present a simplified model of rain water infiltration in the snow coupled to a transient permafrost model. While small amounts of rain have only minor impact on the ground surface temperature, strong rain events have a long-lasting impact. We show that consecutively applying the conditions encountered in the winter season 2005/2006 results in the formation of an unfrozen zone in the soil after three to five years, depending on the prescribed soil properties. If water infiltration in the snow is disabled in the model, more time is required for the permafrost to reach a similar state of degradation.

Description: Semi-Lagrangian methods in air pollution models Geoscientific Model Development, 4, 511-541, 2011 Author(s): A. B. Hansen, J. Brandt, J. H. Christensen, and E. Kaas Various semi-Lagrangian methods are tested with respect to advection in air pollution modeling. The aim is to find a method fulfilling as many of the desirable properties by Rasch andWilliamson (1990) and Machenhauer et al. (2008) as possible. The focus in this study is on accuracy and local mass conservation. The methods tested are, first, classical semi-Lagrangian cubic interpolation, see e.g. Durran (1999), second, semi-Lagrangian cubic cascade interpolation, by Nair et al. (2002), third, semi-Lagrangian cubic interpolation with the modified interpolation weights, Locally Mass Conserving Semi-Lagrangian (LMCSL), by Kaas (2008), and last, semi-Lagrangian cubic interpolation with a locally mass conserving monotonic filter by Kaas and Nielsen (2010). Semi-Lagrangian (SL) interpolation is a classical method for atmospheric modeling, cascade interpolation is more efficient computationally, modified interpolation weights assure mass conservation and the locally mass conserving monotonic filter imposes monotonicity. All schemes are tested with advection alone or with advection and chemistry together under both typical rural and urban conditions using different temporal and spatial resolution. The methods are compared with a current state-of-the-art scheme, Accurate Space Derivatives (ASD), see Frohn et al. (2002), presently used at the National Environmental Research Institute (NERI) in Denmark. To enable a consistent comparison only non-divergent flow configurations are tested. The test cases are based either on the traditional slotted cylinder or the rotating cone, where the schemes' ability to model both steep gradients and slopes are challenged. The tests showed that the locally mass conserving monotonic filter improved the results significantly for some of the test cases, however, not for all. It was found that the semi-Lagrangian schemes, in almost every case, were not able to outperform the current ASD scheme used in DEHM with respect to accuracy.

Description: GEOCLIM reloaded (v 1.0): a new coupled earth system model for past climate change Geoscientific Model Development, 4, 451-481, 2011 Author(s): S. Arndt, P. Regnier, Y. Goddéris, and Y. Donnadieu We present a new version of the coupled Earth system model GEOCLIM. The new release, GEOCLIM reloaded (v 1.0), links the existing atmosphere and weathering modules to a novel, temporally and spatially resolved model of the global ocean circulation, which provides a physical framework for a mechanistic description of the marine biogeochemical dynamics of carbon, nitrogen, phosphorus and oxygen. The ocean model is also coupled to a fully formulated, vertically resolved diagenetic model. GEOCLIM reloaded is thus a unique tool to investigate the short- and long-term feedbacks between climatic conditions, continental inputs, ocean biogeochemical dynamics and diagenesis. A complete and detailed description of the resulting Earth system model and its new features is first provided. The performance of GEOCLIM reloaded is then evaluated by comparing steady-state simulation under present-day conditions with a comprehensive set of oceanic data and existing global estimates of bio-element cycling in the pelagic and benthic compartments.

Description: The surface energy balance of a polygonal tundra site in northern Siberia – Part 2: Winter The Cryosphere, 5, 509-524, 2011 Author(s): M. Langer, S. Westermann, S. Muster, K. Piel, and J. Boike In this study, we present the winter time surface energy balance at a polygonal tundra site in northern Siberia based on independent measurements of the net radiation, the sensible heat flux and the ground heat flux from two winter seasons. The latent heat flux is inferred from measurements of the atmospheric turbulence characteristics and a model approach. The long-wave radiation is found to be the dominant factor in the surface energy balance. The radiative losses are balanced to about 60 % by the ground heat flux and almost 40 % by the sensible heat fluxes, whereas the contribution of the latent heat flux is small. The main controlling factors of the surface energy budget are the snow cover, the cloudiness and the soil temperature gradient. Large spatial differences in the surface energy balance are observed between tundra soils and a small pond. The ground heat flux released at a freezing pond is by a factor of two higher compared to the freezing soil, whereas large differences in net radiation between the pond and soil are only observed at the end of the winter period. Differences in the surface energy balance between the two winter seasons are found to be related to differences in snow depth and cloud cover which strongly affect the temperature evolution and the freeze-up at the investigated pond.

Description: The CSIRO Mk3L climate system model version 1.0 – Part 1: Description and evaluation Geoscientific Model Development, 4, 483-509, 2011 Author(s): S. J. Phipps, L. D. Rotstayn, H. B. Gordon, J. L. Roberts, A. C. Hirst, and W. F. Budd The CSIRO Mk3L climate system model is a coupled general circulation model, designed primarily for millennial-scale climate simulations and palaeoclimate research. Mk3L includes components which describe the atmosphere, ocean, sea ice and land surface, and combines computational efficiency with a stable and realistic control climatology. This paper describes the model physics and software, analyses the control climatology, and evaluates the ability of the model to simulate the modern climate. Mk3L incorporates a spectral atmospheric general circulation model, a z -coordinate ocean general circulation model, a dynamic-thermodynamic sea ice model and a land surface scheme with static vegetation. The source code is highly portable, and has no dependence upon proprietary software. The model distribution is freely available to the research community. A 1000-yr climate simulation can be completed in around one-and-a-half months on a typical desktop computer, with greater throughput being possible on high-performance computing facilities. Mk3L produces realistic simulations of the larger-scale features of the modern climate, although with some biases on the regional scale. The model also produces reasonable representations of the leading modes of internal climate variability in both the tropics and extratropics. The control state of the model exhibits a high degree of stability, with only a weak cooling trend on millennial timescales. Ongoing development work aims to improve the model climatology and transform Mk3L into a comprehensive earth system model.

Description: A full Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution The Cryosphere, 5, 495-508, 2011 Author(s): H. Seddik, R. Greve, T. Zwinger, and L. Placidi A three-dimensional, thermo-mechanically coupled ice flow model with induced anisotropy has been applied to a ~200 × 200 km domain around the Dome Fuji drill site, Antarctica. The model ("Elmer/Ice") is based on the open-source multi-physics package Elmer ( http://www.csc.fi/elmer/ ) and solves the full Stokes equations. Flow-induced anisotropy in ice is accounted for by an implementation of the C ontinuum-mechanical, A nisotropic F low model, based on an anisotropic F low E nhancement factor ("CAFFE model"). Steady-state simulations for present-day climate conditions are conducted. The main findings are: (i) the flow regime at Dome Fuji is a complex superposition of vertical compression, horizontal extension and bed-parallel shear; (ii) for an assumed geothermal heat flux of 60 mW m −2 the basal temperature at Dome Fuji reaches the pressure melting point and the basal melting rate is ~0.35 mm a −1 ; (iii) in agreement with observational data, the fabric shows a strong single maximum at Dome Fuji, and the age of the ice is decreased compared to an isotropic scenario; (iv) as a consequence of spatially variable basal melting conditions, the basal age tends to be smaller where the ice is thicker and larger where the ice is thinner. The latter result is of great relevance for the consideration of a future drill site in the area.

Description: Medusa-1.0: a new intermediate complexity plankton ecosystem model for the global domain Geoscientific Model Development, 4, 381-417, 2011 Author(s): A. Yool, E. E. Popova, and T. R. Anderson The ongoing, anthropogenically-driven changes to the global ocean are expected to have significant consequences for plankton ecosystems in the future. Because of the role that plankton play in the ocean's "biological pump", changes in abundance, distribution and productivity will likely have additional consequences for the wider carbon cycle. Just as in the terrestrial biosphere, marine ecosystems exhibit marked diversity in species and functional types of organisms. Predicting potential change in plankton ecosystems therefore requires the use of models that are suited to this diversity, but whose parameterisation also permits robust and realistic functional behaviour. In the past decade, advances in model sophistication have attempted to address diversity, but have been criticised for doing so inaccurately or ahead of a requisite understanding of underlying processes. Here we introduce MEDUSA-1.0 ( M odel of E cosystem D ynamics, nutrient U tilisation, S equestration and A cidification), a new "intermediate complexity" plankton ecosystem model that expands on traditional nutrient-phytoplankton-zooplankton-detritus (NPZD) models, and remains amenable to global-scale evaluation. MEDUSA-1.0 includes the biogeochemical cycles of nitrogen, silicon and iron, broadly structured into "small" and "large" plankton size classes, of which the "large" phytoplankton class is representative of a key phytoplankton group, the diatoms. A full description of MEDUSA-1.0's state variables, differential equations, functional forms and parameter values is included, with particular attention focused on the submodel describing the export of organic carbon from the surface to the deep ocean. MEDUSA-1.0 is used here in a multi-decadal hindcast simulation, and its biogeochemical performance evaluated at the global scale.

Description: PREP-CHEM-SRC – 1.0: a preprocessor of trace gas and aerosol emission fields for regional and global atmospheric chemistry models Geoscientific Model Development, 4, 419-433, 2011 Author(s): S. R. Freitas, K. M. Longo, M. F. Alonso, M. Pirre, V. Marecal, G. Grell, R. Stockler, R. F. Mello, and M. Sánchez Gácita The preprocessor PREP-CHEM-SRC presented in the paper is a comprehensive tool aiming at preparing emission fields of trace gases and aerosols for use in atmospheric-chemistry transport models. The considered emissions are from the most recent databases of urban/industrial, biogenic, biomass burning, volcanic, biofuel use and burning from agricultural waste sources. For biomass burning, emissions can be also estimated directly from satellite fire detections using a fire emission model included in the tool. The preprocessor provides emission fields interpolated onto the transport model grid. Several map projections can be chosen. The inclusion of these emissions in transport models is also presented. The preprocessor is coded using Fortran90 and C and is driven by a namelist allowing the user to choose the type of emissions and the databases.

Description: Investigating changes in basal conditions of Variegated Glacier prior and during its 1982–1983 surge The Cryosphere Discussions, 5, 1461-1494, 2011 Author(s): M. Jay-Allemand, F. Gillet-Chaulet, O. Gagliardini, and M. Nodet The Variegated Glacier (Alaska) is known to surge periodically after a sufficient amount of cumulative mass balance is reached, but this observation is difficult to link with changes in the basal conditions. Here, using a 10-year dataset, consisting in surface topography and surface velocity observations along a flow line for 25 dates, we have reconstructed the evolution of the basal conditions prior and during the 1982–1983 surge. The model solves the full-Stokes problem along the central flow line using the finite element method. For the 25 dates of the dataset, the basal friction parameter distribution is inferred using the inverse method proposed by Arthern and Gudmundson (2010). This method is here slightly modified by incorporating a regularisation term in the cost function to avoid short wave length changes in the friction parameter. Our results indicate that dramatic changes in the basal conditions occurred between 1973 to 1983. Prior to the surge, periodical changes can be observed between winter and summer, with a regular increase of the sliding from 1973 to 1982. During the surge, the basal friction decreased dramatically and an area of very low friction moved from the upper part of the glacier to its terminus. Using a more complex friction law, these changes in basal sliding are then interpreted in terms of basal water pressure. It confirms that dramatic changes took place in the subglacial drainage system of Variegated Glacier, moving from a relatively efficient drainage system prior to the surge to an inefficient one during the surge. By reconstructing the water pressure evolution at the base of the glacier it is possible to infer realistic scenarios for the hydrological history leading to the occurrence of a surge.

Description: Scale-dependent measurement and analysis of ground surface temperature variability in alpine terrain The Cryosphere, 5, 431-443, 2011 Author(s): S. Gubler, J. Fiddes, M. Keller, and S. Gruber Measurements of environmental variables are often used to validate and calibrate physically-based models. Depending on their application, the models are used at different scales, ranging from few meters to tens of kilometers. Environmental variables can vary strongly within the grid cells of these models. Validating a model with a single measurement is therefore delicate and susceptible to induce bias in further model applications. To address the question of uncertainty associated with scale in permafrost models, we present data of 390 spatially-distributed ground surface temperature measurements recorded in terrain of high topographic variability in the Swiss Alps. We illustrate a way to program, deploy and refind a large number of measurement devices efficiently, and present a strategy to reduce data loss reported in earlier studies. Data after the first year of deployment is presented. The measurements represent the variability of ground surface temperatures at two different scales ranging from few meters to some kilometers. On the coarser scale, the dependence of mean annual ground surface temperature on elevation, slope, aspect and ground cover type is modelled with a multiple linear regression model. Sampled mean annual ground surface temperatures vary from −4 °C to 5 °C within an area of approximately 16 km 2 subject to elevational differences of approximately 1000 m. The measurements also indicate that mean annual ground surface temperatures vary up to 6 °C (i.e., from −2 °C to 4 °C) even within an elevational band of 300 m. Furthermore, fine-scale variations can be high (up to 2.5 °C) at distances of less than 14 m in homogeneous terrain. The effect of this high variability of an environmental variable on model validation and applications in alpine regions is discussed.

Description: Tropical troposphere to stratosphere transport of carbon monoxide and long-lived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development, 7, 2895-2916, 2014 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, F. Ploeger, B. Vogel, M. Tao, C. M. Hoppe, G. Günther, N. Spelten, L. Hoffmann, H.-C. Pumphrey, S. Viciani, F. D'Amato, C. M. Volk, P. Hoor, H. Schlager, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact of these transport fluxes on the composition of the tropical lower stratosphere. Anomaly patterns of carbon monoxide (CO) and long-lived tracers in the lower tropical stratosphere allow conclusions about the rate and the variability of tropical upwelling to be drawn. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F (CFC-11), CCl 2 F 2 (CFC-12), and CO 2 ) in the lower tropical stratosphere. For the long-lived trace substances, the boundary conditions at the surface are prescribed based on ground-based measurements in the lowest model level. The boundary condition for CO in the lower troposphere (below about 4 km) is deduced from MOPITT measurements. Due to the lack of a specific representation of mixing and convective uplift in the troposphere in this model version, enhanced CO values, in particular those resulting from convective outflow are underestimated. However, in the tropical tropopause layer and the lower tropical stratosphere, there is relatively good agreement of simulated CO with in situ measurements (with the exception of the TROCCINOX campaign, where CO in the simulation is biased low ≈10–15 ppbv). Further, the model results (and therefore also the ERA-Interim winds, on which the transport in the model is based) are of sufficient quality to describe large scale anomaly patterns of CO in the lower stratosphere. In particular, the zonally averaged tropical CO anomaly patterns (the so called "tape recorder" patterns) simulated by this model version of CLaMS are in good agreement with observations, although the simulations show a too rapid upwelling compared to observations as a consequence of the overestimated vertical velocities in the ERA-Interim reanalysis data set. Moreover, the simulated tropical anomaly patterns of N 2 O are in good agreement with observations. In the simulations, anomaly patterns of CH 4 and CFC-11 were found to be very similar to those of N 2 O; for all long-lived tracers, positive anomalies are simulated because of the enhanced tropical upwelling in the easterly shear phase of the quasi-biennial oscillation.

Description: Comparison of aeolian snow transport events and snow mass fluxes between observations and simulations made by the regional climate model MAR in Adélie Land, East Antarctica The Cryosphere Discussions, 8, 6007-6032, 2014 Author(s): A. Trouvilliez, H. Gallée, F. Naaim-Bouvet, C. Genthon, C. Amory, V. Favier, C. Agosta, L. Piard, and H. Bellot The regional climate model MAR including a coupled snow pack/aeolian snow transport parameterisation is compared with aeolian snow mass fluxes at a fine spatial resolution (5 km horizontally and 2 m vertically) and at a fine temporal resolution (30 min) over 1 month in Antarctica. Numerous feedbacks are taken into account in the MAR including the drag partitioning caused by the roughness elements. Wind speed is correctly simulated with a positive value of the Nash test (0.60 and 0.37) but the wind speeds above 10 m s −1 are underestimated. The aeolian snow transport events are correctly reproduced with a good temporal resolution except for the aeolian snow transport events with a particles' maximum height below 1 m. The simulated threshold friction velocity, calculated without snowfall, is overestimated. The simulated aeolian snow mass fluxes between 0 to 2 m have the same variations but are underestimated compared to the second-generation FlowCapt values and so is the simulated relative humidity at 2 m. This underestimation is not entirely due to the underestimation of the simulated wind speed. The MAR underestimates the aeolian snow quantity that pass through the first two meters by a factor ten compared to the second-generation FlowCapt value (13 990 kg m −1 and 151 509 kg m −1 respectively). It will conduct the MAR, with this parametrisation, to underestimate the effect of the aeolian snow transport on the Antarctic surface mass balance.

Description: Influence of weak layer heterogeneity and slab properties on slab tensile failure propensity and avalanche release area The Cryosphere Discussions, 8, 6033-6057, 2014 Author(s): J. Gaume, G. Chambon, N. Eckert, M. Naaim, and J. Schweizer Dry-snow slab avalanches are generally caused by a sequence of fracture processes including failure initiation in a weak snow layer underlying a cohesive slab followed by crack propagation within the weak layer (WL) and tensile fracture through the slab. During past decades, theoretical and experimental work has gradually improved our knowledge of the fracture process in snow. However, our limited understanding of crack propagation and fracture arrest propensity prevents the evaluation of avalanche release sizes and thus impedes hazard assessment. To address this issue, slab tensile failure propensity is examined using a mechanically-based statistical model of the slab–WL system based on the finite element method. This model accounts for WL heterogeneity, stress redistribution by elasticity of the slab and the slab possible tensile failure. Two types of avalanche release are distinguished in the simulations: (1) full-slope release if the heterogeneity is not sufficient to stop crack propagation and to trigger a tensile failure within the slab, (2) partial-slope release if fracture arrest and slab tensile failure occurs due to the WL heterogeneity. The probability of these two release types is presented as a function of the characteristics of WL heterogeneity and of the slab. One of the main outcomes is that, for realistic values of the parameters, the tensile failure propensity is mainly influenced by slab properties. Hard and thick snow slabs are more prone to wide-scale crack propagation and thus lead to larger avalanches (full-slope release). In this case, the avalanche size is mainly influenced by topographical and morphological features such as rocks, trees, slope curvature and the spatial variability of the snow depth as it is often claimed in the literature.

Description: A new approach to mapping permafrost and change incorporating uncertainties in ground conditions and climate projections The Cryosphere, 8, 2177-2194, 2014 Author(s): Y. Zhang, I. Olthof, R. Fraser, and S. A. Wolfe Spatially detailed information on permafrost distribution and change with climate is important for land use planning, infrastructure development, and environmental assessments. However, the required soil and surficial geology maps in the North are coarse, and projected climate scenarios vary widely. Considering these uncertainties, we propose a new approach to mapping permafrost distribution and change by integrating remote sensing data, field measurements, and a process-based model. Land cover types from satellite imagery are used to capture the general land conditions and to improve the resolution of existing permafrost maps. For each land cover type, field observations are used to estimate the probabilities of different ground conditions. A process-based model is used to quantify the evolution of permafrost for each ground condition under three representative climate scenarios (low, medium, and high warming). From the model results, the probability of permafrost occurrence and the most likely permafrost conditions are determined. We apply this approach at 20 m resolution to a large area in Northwest Territories, Canada. Mapped permafrost conditions are in agreement with field observations and other studies. The data requirements, model robustness, and computation time are reasonable, and this approach may serve as a practical means to mapping permafrost and changes at high resolution in other regions.

Description: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice The Cryosphere, 8, 2219-2233, 2014 Author(s): S. Arndt and M. Nicolaus Arctic sea ice has not only decreased in volume during the last decades, but has also changed in its physical properties towards a thinner and more seasonal ice cover. These changes strongly impact the energy budget, and might affect the ice-associated ecosystems. In this study, we quantify solar shortwave fluxes through sea ice for the entire Arctic during all seasons. To focus on sea-ice-related processes, we exclude fluxes through open water, scaling linearly with sea ice concentration. We present a new parameterization of light transmittance through sea ice for all seasons as a function of variable sea ice properties. The maximum monthly mean solar heat flux under the ice of 30 × 10 5 Jm −2 occurs in June, enough heat to melt 0.3 m of sea ice. Furthermore, our results suggest that 96% of the annual solar heat input through sea ice occurs during only a 4-month period from May to August. Applying the new parameterization to remote sensing and reanalysis data from 1979 to 2011, we find an increase in transmitted light of 1.5% yr −1 for all regions. This corresponds to an increase in potential sea ice bottom melt of 63% over the 33-year study period. Sensitivity studies reveal that the results depend strongly on the timing of melt onset and the correct classification of ice types. Assuming 2 weeks earlier melt onset, the annual transmitted solar radiation to the upper ocean increases by 20%. Continuing the observed transition from a mixed multi-year/first-year sea ice cover to a seasonal ice cover results in an increase in light transmittance by an additional 18%.

Description: Climatic signals from 76 shallow firn cores in Dronning Maud Land, East Antarctica The Cryosphere Discussions, 8, 5961-6005, 2014 Author(s): S. Altnau, E. Schlosser, E. Isaksson, and D. Divine The spatial and temporal distribution of surface mass balance (SMB) and δ 18 O were investigated in the first comprehensive study of a set of 76 firn cores retrieved by various expeditions during the past three decades in Dronning Maud Land, East Antarctica. The large number of cores was used to calculate stacked records of SMB and δ 18 O, which considerably increased the signal-to-noise ratio compared to earlier studies and facilitated the detection of climatic signals. Considerable differences between cores from the interior plateau and the coastal cores were found. The δ 18 O of both the plateau and the ice shelf cores exhibit a slight positive trend over the second half of the 20th century. In the corresponding period, the SMB has a negative trend in the ice shelf cores, but increases on the plateau. Comparison with meteorological data from Neumayer Station revealed that for the ice shelf regions atmospheric dynamic effects are more important than thermodynamics, while on the plateau, the temporal variations of SMB and δ 18 O occur mostly in parallel, thus can be explained by thermodynamic effects. The Southern Annular Mode (SAM) exhibits a positive trend since the mid-1960s, which is assumed to lead to a cooling of East Antarctica. This is not confirmed by the firn core data in our data set. Changes in the atmospheric circulation that result in a changed seasonal distribution of precipitation/accumulation could partly explain the observed features in the ice shelf cores.

Description: Post-LIA glacier changes along a latitudinal transect in the Central Italian Alps The Cryosphere, 8, 2235-2252, 2014 Author(s): R. Scotti, F. Brardinoni, and G. B. Crosta The variability of glacier response to atmospheric temperature rise in different topo-climatic settings is still a matter of debate. To address this question in the Central Italian Alps, we compile a post-LIA (Little Ice Age) multitemporal glacier inventory (1860–1954–1990–2003–2007) along a latitudinal transect that originates north of the continental divide in the Livigno Mountains and extends south through the Disgrazia and Orobie ranges, encompassing continental-to-maritime climatic settings. In these sub-regions, we examine the area change of 111 glaciers. Overall, the total glacierized area has declined from 34.1 to 10.1 km 2 , with a substantial increase in the number of small glaciers due to fragmentation. The average annual decrease (AAD) in glacier area has risen by about 1 order of magnitude from 1860–1990 (Livigno: 0.45; Orobie: 0.42; and Disgrazia: 0.39 % a −1 ) to 1990–2007 (Livigno: 3.08; Orobie: 2.44; and Disgrazia: 2.27 % a −1 ). This ranking changes when considering glaciers smaller than 0.5 km 2 only (i.e., we remove the confounding caused by large glaciers in Disgrazia), so that post-1990 AAD follows the latitudinal gradient and Orobie glaciers stand out (Livigno: 4.07; Disgrazia: 3.57; and Orobie: 2.47 % a −1 ). More recent (2007–2013) field-based mass balances in three selected small glaciers confirm post-1990 trends showing the consistently highest retreat in continental Livigno and minimal area loss in maritime Orobie, with Disgrazia displaying transitional behavior. We argue that the recent resilience of glaciers in Orobie is a consequence of their decoupling from synoptic atmospheric temperature trends, a decoupling that arises from the combination of local topographic configuration (i.e., deep, north-facing cirques) and high winter precipitation, which ensures high snow-avalanche supply, as well as high summer shading and sheltering. Our hypothesis is further supported by the lack of correlations between glacier change and glacier attributes in Orobie, as well as by the higher variability in ELA,sub〉0 positioning, post-LIA glacier change, and interannual mass balances, as we move southward along the transect.

Description: Evaluation of the ECHAM family radiation codes performance in the representation of the solar signal Geoscientific Model Development, 7, 2859-2866, 2014 Author(s): T. Sukhodolov, E. Rozanov, A. I. Shapiro, J. Anet, C. Cagnazzo, T. Peter, and W. Schmutz Solar radiation is the main source of energy for the Earth's atmosphere and in many respects defines its composition, photochemistry, temperature profile and dynamics. The magnitude of the solar irradiance variability strongly depends on the wavelength, making difficult its representation in climate models. Due to some deficiencies in the applied radiation codes, several models fail to show a clear response in middle stratospheric heating rates to solar spectral irradiance variability; therefore, it is important to evaluate model performance in this respect before doing multiple runs. In this work we evaluate the performance of three generations of ECHAM (4, 5 and 6) solar radiation schemes by a comparison with the reference high-resolution libRadtran code. We found that all original ECHAM radiation codes miss almost all solar signals in the heating rates in the mesosphere. In the stratosphere the two-band ECHAM4 code (E4) has an almost negligible radiative response to solar irradiance changes and the six-band ECHAM5 code (E5c) reproduces only about half of the reference signal, while representation in the ECHAM6 code (E6) is better – it misses a maximum of about 15% in the upper stratosphere. On the basis of the comparison results we suggest necessary improvements to the ECHAM family codes by the inclusion of available parameterizations of the heating rate due to absorption by oxygen (O 2 ) and ozone (O 3 ). Improvement is presented for E5c and E6, and both codes, with the introduced parameterizations, represent the heating rate response to the spectral solar irradiance variability simulated with libRadtran much better without a substantial increase in computer time. The suggested parameterizations are recommended to be applied in the middle-atmosphere version of the ECHAM-5 and 6 models for the study of the solar irradiance influence on climate.

Description: GEOtop 2.0: simulating the combined energy and water balance at and below the land surface accounting for soil freezing, snow cover and terrain effects Geoscientific Model Development, 7, 2831-2857, 2014 Author(s): S. Endrizzi, S. Gruber, M. Dall'Amico, and R. Rigon GEOtop is a fine-scale grid-based simulator that represents the heat and water budgets at and below the soil surface. It describes the three-dimensional water flow in the soil and the energy exchange with the atmosphere, considering the radiative and turbulent fluxes. Furthermore, it reproduces the highly non-linear interactions between the water and energy balance during soil freezing and thawing, and simulates the temporal evolution of the water and energy budgets in the snow cover and their effect on soil temperature. Here, we present the core components of GEOtop 2.0 and demonstrate its functioning. Based on a synthetic simulation, we show that the interaction of processes represented in GEOtop 2.0 can result in phenomena that are significant and relevant for applications involving permafrost and seasonally frozen soils, both in high altitude and latitude regions.

Description: The North American Carbon Program Multi-scale Synthesis and Terrestrial Model Intercomparison Project – Part 2: Environmental driver data Geoscientific Model Development, 7, 2875-2893, 2014 Author(s): Y. Wei, S. Liu, D. N. Huntzinger, A. M. Michalak, N. Viovy, W. M. Post, C. R. Schwalm, K. Schaefer, A. R. Jacobson, C. Lu, H. Tian, D. M. Ricciuto, R. B. Cook, J. Mao, and X. Shi Ecosystems are important and dynamic components of the global carbon cycle, and terrestrial biospheric models (TBMs) are crucial tools in further understanding of how terrestrial carbon is stored and exchanged with the atmosphere across a variety of spatial and temporal scales. Improving TBM skills, and quantifying and reducing their estimation uncertainties, pose significant challenges. The Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal multi-scale and multi-model intercomparison effort set up to tackle these challenges. The MsTMIP protocol prescribes standardized environmental driver data that are shared among model teams to facilitate model–model and model–observation comparisons. This paper describes the global and North American environmental driver data sets prepared for the MsTMIP activity to both support their use in MsTMIP and make these data, along with the processes used in selecting/processing these data, accessible to a broader audience. Based on project needs and lessons learned from past model intercomparison activities, we compiled climate, atmospheric CO 2 concentrations, nitrogen deposition, land use and land cover change (LULCC), C3 / C4 grasses fractions, major crops, phenology and soil data into a standard format for global (0.5° × 0.5° resolution) and regional (North American: 0.25° × 0.25° resolution) simulations. In order to meet the needs of MsTMIP, improvements were made to several of the original environmental data sets, by improving the quality, and/or changing their spatial and temporal coverage, and resolution. The resulting standardized model driver data sets are being used by over 20 different models participating in MsTMIP. The data are archived at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov ) to provide long-term data management and distribution.

Description: Melting trends over the Greenland ice sheet (1958–2009) from spaceborne microwave data and regional climate models The Cryosphere, 5, 359-375, 2011 Author(s): X. Fettweis, M. Tedesco, M. van den Broeke, and J. Ettema To study near-surface melt changes over the Greenland ice sheet (GrIS) since 1979, melt extent estimates from two regional climate models were compared with those obtained from spaceborne microwave brightness temperatures using two different remote sensing algorithms. The results from the two models were consistent with those obtained with the remote sensing algorithms at both daily and yearly time scales, encouraging the use of the models for analyzing melting trends before the satellite era (1958–1979), when forcing data is available. Differences between satellite-derived and model-simulated results still occur and are used here to identify (i) biases in the snow models (notably in the albedo parametrization, in the thickness of a snow layer, in the maximum liquid water content within the snowpack and in the snowfall impacting the bare ice appearance in summer) and (ii) limitations in the use of passive microwave data for snowmelt detection at the edge of the ice sheet due to mixed pixel effect (e.g., tundra or rock nearby the ice sheet). The results from models and spaceborne microwave sensors confirm a significant (p-value = 0.01) increase in GrIS surface melting since 1979. The melt extent recorded over the last years (1998, 2003, 2005 and 2007) is unprecedented in the last 50 yr with the cumulated melt area in the 2000's being, on the average, twice that of the 1980's.

Description: Warming of waters in an East Greenland fjord prior to glacier retreat: mechanisms and connection to large-scale atmospheric conditions The Cryosphere Discussions, 5, 1335-1364, 2011 Author(s): P. Christoffersen, R. I. Mugford, K. J. Heywood, I. Joughin, J. A. Dowdeswell, J. P. M. Syvitski, A. Luckman, and T. J. Benham Hydrographic data acquired in Kangerlugssuaq Fjord and adjacent seas in 1993 and 2004 are used together with ocean reanalysis to elucidate water mass change and ice-ocean-atmosphere interactions in East Greenland. The hydrographic data show substantial warming of fjord waters between 1993 and 2004 and warm subsurface conditions coincide with the rapid retreat of Kangerlugssuaq Glacier in 2004–2005. The ocean reanalysis shows that the warm properties of fjord waters in 2004 are related to a major peak in oceanic shoreward heat flux into a cross-shelf trough on the outer continental shelf. The heat flux into this trough varies according to seasonal exchanges with the atmosphere as well as from deep seasonal intrusions of subtropical waters. Both mechanisms contribute to high (low) shoreward heat flux when winds from the northeast are weak (strong). The combined effect of surface heating and inflow of subtropical waters is seen in the hydrographic data, which were collected after periods when along-shore coastal winds from the north were strong (1993) and weak (2004). We show that coastal winds vary according to the pressure gradient defined by a semi-permanent atmospheric pressure system over Greenland and a persistent atmospheric low situated near Iceland. The magnitude of this pressure gradient is controlled by longitudinal variability in the position of the Icelandic Low.

Description: Recent wind driven high sea ice export in the Fram Strait contributes to Arctic sea ice decline The Cryosphere Discussions, 5, 1311-1334, 2011 Author(s): L. H. Smedsrud, A. Sirevaag, K. Kloster, A. Sorteberg, and S. Sandven Arctic sea ice area decrease has been visible for two decades, and continues at a steady rate. Apart from melting, the southward drift through Fram Strait is the main loss. We present high resolution sea ice drift across 79° N from 2004 to 2010. The ice drift is based on radar satellite data and correspond well with variability in local geostrophic wind. The underlying current contributes with a constant southward speed close to 5 cm s −1 , and drives about 33 % of the ice export. We use geostrophic winds derived from reanalysis data to calculate the Fram Strait ice area export back to 1957, finding that the sea ice area export recently is about 25 % larger than during the 1960's. The increase in ice export occurred mostly during winter and is directly connected to higher southward ice drift velocities, due to stronger geostrophic winds. The increase in ice drift is large enough to counteract a decrease in ice concentration of the exported sea ice. Using storm tracking we link changes in geostrophic winds to more intense Nordic Sea low pressure systems. Annual sea ice export likely has a significant influence on the summer sea ice variability and we find low values in the 60's, the late 80's and 90's, and particularly high values during 2005–2008. The study highlight the possible role of variability in ice export as an explanatory factor for understanding the dramatic loss of Arctic sea ice the last decades.

Description: The atmospheric chemistry box model CAABA/MECCA-3.0 Geoscientific Model Development, 4, 373-380, 2011 Author(s): R. Sander, A. Baumgaertner, S. Gromov, H. Harder, P. Jöckel, A. Kerkweg, D. Kubistin, E. Regelin, H. Riede, A. Sandu, D. Taraborrelli, H. Tost, and Z.-Q. Xie We present version 3.0 of the atmospheric chemistry box model CAABA/MECCA. In addition to a complete update of the rate coefficients to the most recent recommendations, a number of new features have been added: chemistry in multiple aerosol size bins; automatic multiple simulations reaching steady-state conditions; Monte-Carlo simulations with randomly varied rate coefficients within their experimental uncertainties; calculations along Lagrangian trajectories; mercury chemistry; more detailed isoprene chemistry; tagging of isotopically labeled species. Further changes have been implemented to make the code more user-friendly and to facilitate the analysis of the model results. Like earlier versions, CAABA/MECCA-3.0 is a community model published under the GNU General Public License.

Description: Spatio-temporal measurements and analysis of snow depth in a rock face The Cryosphere Discussions, 5, 1383-1418, 2011 Author(s): V. Wirz, M. Schirmer, S. Gruber, and M. Lehning Snow in rock faces plays a key role in the alpine environment for permafrost distribution, snow water storage or run off in spring. However, a detailed assessment of snow depths in steep rock walls has never been attempted. To understand snow distribution in rock walls a high-resolution terrestrial laser scanner (TLS), including a digital camera, was used to obtain snow depth (HS) data with a resolution of one metre. The mean HS, the snow covered area and their evolution in the rock face were compared to a neighbouring smoother catchment and a flat field station at similar elevation. Further we analyzed the patterns of HS distribution in the rock face after different periods and investigated the main factors contributing to them. In a first step we could show that with TLS reliable information on surface data of a steep rocky surface can be obtained. In comparison to the flatter sites in the vicinity, mean HS in the rock face was lower during the entire winter, but trends of snow depth changes were similar. We observed repeating accumulation and ablation patterns in the rock face, while maximum snow depth loss always occurred at those places with maximum snow depth gain. Further analysis of the main factors contributing to the snow depth distribution in the rock face revealed terrain-wind-interaction processes to be dominant. Processes related to slope angle seem to play a role, but no linear function of slope angle and snow depth was found. Further analyses should involve measurements in rock faces with other characteristics and higher temporal resolutions to be able to distinguish individual processes better. Additionally the relation of spatial and temporal distribution of snow depth to terrain-wind interactions should be tested.

Description: Utility of late summer transient snowline migration rate on Taku Glacier, Alaska The Cryosphere Discussions, 5, 1365-1382, 2011 Author(s): M. Pelto On Taku Glacier, Alaska a combination of field observations of snow water equivalent (SWE) from snowpits and probing in the vicinity of the transient snowline (TSL) are used to quantify the mass balance gradient. The balance gradient is determined from the difference in elevation and SWE from the TSL to snowpits at 1000 m from 1998–2010 and ranges from 2.6–3.8 mm m −1 . Probing transects from 950 m–1100 m directly measure SWE and yield a slightly higher balance gradient of 3.3–3.8 mm m −1 . TSL is identified in MODIS and Landsat 4 and 7 Thematic Mapper imagery for 31 dates during the 2004–2010 period on Taku Glacier to assess the consistency of its rate of rise and usefulness in assessing mass balance. In 2010, the TSL rose from 750 m on 28 July, 800 m on 5 August, 875 m on 14 August, 925 m on 30 August, and to 975 m on 20 September. The mean observed probing balance gradient was 3.3 mm m −1 and TSL rise was 3.7 m day −1 , yielding an ablation rate of 12.2 mm day −1 on Taku Glacier from mid-July to mid-September. A comparison of the TSL rise in the region from 750–1100 m on Taku Glacier during eleven different periods of more than 14 days during the ablation season with repeat imagery indicates a mean TSL rise of 3.7 m day −1 on Taku Glacier, the rate of rise is relatively consistent ranging from 3.0 to 4.8 m day −1 . This is useful for ascertaining the final ELA if imagery or observations are not available within a week or two of the end of the ablation season. From mid-July-mid-September the mean ablation from 750–1100 m determined from the TSL rise and the observed balance gradient varied from 11 to 18 mm day −1 on Taku Glacier during the 2004–2010 period.

Description: Isotope hydrological studies of the perennial ice deposit of Saarhalle, Mammuthöhle, Dachstein Mts, Austria The Cryosphere, 5, 291-298, 2011 Author(s): Z. Kern, I. Fórizs, R. Pavuza, M. Molnár, and B. Nagy A 5.28 m-long ice core was extracted from a major cave ice body in the Mammuthöhle cave system. The upper ~1.2 m of ice most likely originate from precipitation fallen before the 1960s (based on

Description: Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change The Cryosphere, 5, 271-290, 2011 Author(s): C. Nuth and A. Kääb There are an increasing number of digital elevation models (DEMs) available worldwide for deriving elevation differences over time, including vertical changes on glaciers. Most of these DEMs are heavily post-processed or merged, so that physical error modelling becomes difficult and statistical error modelling is required instead. We propose a three-step methodological framework for assessing and correcting DEMs to quantify glacier elevation changes: (i) remove DEM shifts, (ii) check for elevation-dependent biases, and (iii) check for higher-order, sensor-specific biases. A simple, analytic and robust method to co-register elevation data is presented in regions where stable terrain is either plentiful (case study New Zealand) or limited (case study Svalbard). The method is demonstrated using the three global elevation data sets available to date, SRTM, ICESat and the ASTER GDEM, and with automatically generated DEMs from satellite stereo instruments of ASTER and SPOT5-HRS. After 3-D co-registration, significant biases related to elevation were found in some of the stereoscopic DEMs. Biases related to the satellite acquisition geometry (along/cross track) were detected at two frequencies in the automatically generated ASTER DEMs. The higher frequency bias seems to be related to satellite jitter , most apparent in the back-looking pass of the satellite. The origins of the more significant lower frequency bias is uncertain. ICESat-derived elevations are found to be the most consistent globally available elevation data set available so far. Before performing regional-scale glacier elevation change studies or mosaicking DEMs from multiple individual tiles (e.g. ASTER GDEM), we recommend to co-register all elevation data to ICESat as a global vertical reference system.