ALBERT

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: 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: Halocline water modification and along slope advection at the Laptev Sea continental margin Ocean Science Discussions, 10, 1581-1617, 2013 Author(s): D. Bauch, S. Torres-Valdes, I. Polyakov, A. Novikhin, I. Dmitrenko, J. McKay, and A. Mix A general pattern in water mass distribution and potential shelf-basin exchanges is revealed at the Laptev Sea continental slope based on hydrochemical and stable oxygen isotope data from summers 2005–2009. Despite considerable interannual variations, a frontal system can be inferred between shelf, continental slope and central Eurasian Basin waters in the upper 100 m of the water column along the continental slope. Net sea-ice melt is consistently found at the continental slope; however the sea-ice meltwater signal is independent from the local retreat of the ice cover and appears to be advected from upwind locations. In addition to the along-slope frontal system at the continental shelf break a strong gradient is identified on the Laptev Sea shelf between 122 and 126° E with an eastward increase of riverine and sea-ice related brine water contents. These waters cross the shelf break at ~ 140° E and feed the Low Salinity Halocline Water (LSHW, salinity S 〈 33) in the upper 50 m of the water column. Extremely high silicate concentrations in Laptev Sea bottom waters may lead to speculation on a link to the local silicate maximum found within the salinity range of ~ 33 to 34.5, typical for the Lower Halocline Water (LHW) at the continental slope. But brine signatures and nutrient ratios from the central Laptev Sea differ from those at the continental slope. Thus a significant contribution of Laptev Sea bottom waters to the LHW at the continental slope can be excluded. The silicate maximum within the LHW at the continental slope may be formed locally or at the outer Laptev Sea shelf. Similar to the advection of the sea-ice melt signal along the Laptev Sea continental slope the nutrient signal at 50–70 m water depth within the LHW might also be fed by advection parallel to the slope. Thus, our analyses suggest that advective processes from upwind locations play a significant role in the halocline formation in the northern Laptev Sea.

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: 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: 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: 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: 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: The transient distributions of nuclear weapon-generated tritium and its decay product 3 He in the Mediterranean Sea, 1952–2011, and their oceanographic potential Ocean Science Discussions, 10, 649-690, 2013 Author(s): W. Roether, P. Jean-Baptiste, E. Fourré, and J. Sültenfuß We present a comprehensive account of tritium and 3 He in the Mediterranean Sea since the appearance of the tritium generated by the atmospheric nuclear-weapon testing in the 1950's and early 1960's, based on essentially all available observations. Tritium in surface waters rose to 20–30 TU in 1964 (TU = 10 18 · [ 3 H]/[H]), a factor of about 100 above the natural level, and thereafter declined 30-fold up to 2011. The decline was largely due to radioactive tritium decay, which produced significant amounts of its stable daughter 3 He. We present the scheme by which we separate the tritiugenic part of 3 He and the part due to release from the sea floor (terrigenic part). We show that the tritiugenic component can be quantified throughout the Mediterranean waters, typically to a ±0.15 TU equivalent, mostly because the terrigenic part is low in 3 He. This fact makes the Mediterranean unique in offering a potential for the use of tritiugenic 3 He as a tracer. The transient distributions of the two tracers are illustrated by a number of sections spanning the entire sea and relevant features of their distributions are noted. By 2011, the 3 He concentrations in the top few hundred meters had become low, in response to the decreasing tritium concentrations combined with a flushing out by the general westward drift of these waters. Tritium- 3 He ages in Levantine Intermediate Water (LIW) were obtained repeated in time at different locations, defining transit times from the LIW source region east of Rhodes. The ages show an upward trend with the time elapsed since the surface-water tritium maximum, which arises because the repeated observations represent increasingly slower moving parts of the full transit time spectrum of LIW. The transit time dispersion found by this new application of tritium- 3 He dating is considerable. We find mean transit times of 12 ± 2 yr up to the Strait of Sicily, 18 ± 3 yr up to the Tyrrhenian Sea, and 22 ± 4 yr up into the Western Mediterranean. We furthermore present full Eastern Mediterranean sections of terrigenic 3 He and tritium- 3 He age in 1987, the latter one similarly showing an effect of the transit time dispersion. We conclude that the available tritium and 3 He data, in particular if combined with other tracer data, are useful for constraining the subsurface circulation and mixing of the Mediterranean Sea.

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: 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: The effects of climate change on persistent organic pollutants (POPs) in the North Sea Ocean Science Discussions, 10, 1525-1557, 2013 Author(s): K. O'Driscoll, B. Mayer, J. Su, and M. Mathis The fate and cycling of two selected legacy persistent organic pollutants (POPs), PCB 153 and γ-HCH, in the North Sea in the 21st century have been modelled with combined hydrodynamic and fate and transport ocean models. To investigate the impact of climate variability on POPs in the North Sea in the 21st century, future scenario model runs for three 10 yr periods to the year 2100 using plausible levels of both in situ concentrations and atmospheric, river and open boundary inputs are performed. Since estimates of future concentration levels of POPs in the atmosphere, oceans and rivers are not available, our approach was to reutilise 2005 values in the atmosphere, rivers and at the open ocean boundaries for every year of the simulations. In this way, we attribute differences between the three 10 yr simulations to climate change only. For the HAMSOM and atmospheric forcing, results of the IPCC A1B (SRES) 21st century scenario are utilised, where surface forcing is provided by the REMO downscaling of the ECHAM5 global atmospheric model, and open boundary conditions are provided by the MPIOM global ocean model. Dry gas deposition and volatilisation of γ-HCH increase in the future relative to the present. In the water column, total mass of γ-HCH and PCB 153 remain fairly steady in all three runs. In sediment, γ-HCH increases in the future runs, relative to the present, while PCB 153 in sediment decreases exponentially in all three runs, but even faster in the future, both of which are the result of climate change. Annual net sinks exceed sources at the ends of all periods.

Description: Observations of phytoplankton spring bloom onset triggered by a density front in NW Mediterranean Ocean Science Discussions, 10, 1559-1580, 2013 Author(s): A. Olita, S. Sparnocchia, S. Cusí, L. Fazioli, R. Sorgente, J. Tintoré, and A. Ribotti Phytoplankton bloom in NW Mediterranan sea is a seasonal event that mainly occurrs in a limited area (Gulf of Lyon and Provençal basin) where this phenomenon is promoted by a cyclonic circulation, strong wind-driven mixing and subsequent spring restratification. At the southern boundary of this area a density front (North Balearic Front) separating denser waters from the lighter Modified Atlantic Waters reservoir at south is suspected to trigger weaker and earlier (late winter) blooms by (a) enhanced pumping of nutrients into the euphotic layer and (b) promoting an early restratification of the water column (by frontal instabilities). A multisensor glider round trip, equipped with CTD and fluorimeter, crossing the frontal area in February–March 2013, allowed to observe the bloom triggering after the decrease of intense wind-driven turbulent convection and mixing. Satellite imagery supports and confirms in-situ observations. It was shown that frontal activity has a relevant role in the promotion and acceleration of the dynamical restratification, with a consequent biological response in terms of primary production. Restratification is necessary preconditioning factor for bloom triggering in frontal area, net of other involved mechanism promoting the bloom as the enhanced biological pump. So, like for high-latitude fronts (Taylor and Ferrari, 2011a), also for this mid-latitude oligotrophic region front seems to promote new production by dynamically enahnced restratification inhibiting mixing. Finally, we argued that Sverdrup's Critical Depth criterion seems to apply in the northern well-mixed area, where the zeroing of heat fluxes (and related turbulent convection) does not correspond to a prompt onset of the bloom (which appeared 1 month later).

Description: Observed decline of the Atlantic Meridional Overturning Circulation 2004 to 2012 Ocean Science Discussions, 10, 1619-1645, 2013 Author(s): D. A. Smeed, G. McCarthy, S. A. Cunningham, E. Frajka-Williams, D. Rayner, W. E. Johns, C. S. Meinen, M. O. Baringer, B. I. Moat, A. Duchez, and H. L. Bryden The Atlantic Meridional Overturning Circulation (AMOC) has been observed continuously at 26° N since April 2004. The AMOC and its component parts are monitored by combining a transatlantic array of moored instruments with submarine-cable based measurements of the Gulf Stream and satellite derived Ekman transport. The time series has recently been extended to October 2012 and the results show a downward trend since 2004. From April~2008 to March 2012 the AMOC was an average of 2.7 Sv weaker than in the first four years of observation (95% confidence that the reduction is 0.3 Sv or more). Ekman transport reduced by about 0.2 Sv and the Gulf Stream by 0.5 Sv but most of the change (2.0 Sv) is due to the mid-ocean geostrophic flow. The change of the mid-ocean geostrophic flow represents a strengthening of the subtropical gyre above the thermocline. The increased southward flow of warm waters is balanced by a decrease in the southward flow of Lower North Atlantic Deep Water below 3000 m. The transport of Lower North Atlantic Deep Water slowed by 7% per year (95% confidence that the rate of slowing is greater than 2.5% per year).

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: 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: Global representation of tropical cyclone-induced ocean thermal changes using Argo data – Part 2: Estimating air–sea heat fluxes and ocean heat content changes Ocean Science Discussions, 11, 2907-2937, 2014 Author(s): L. Cheng, J. Zhu, and R. L. Sriver We use Argo temperature data to examine changes in ocean heat content (OHC) and air–sea heat fluxes induced by tropical cyclones (TC)s on a global scale. A footprint technique that analyzes the vertical structure of cross-track thermal responses along all storm tracks during the period 2004–2012 is utilized (see part I). We find that TCs are responsible for 1.87 PW (11.05 W m −2 when averaging over the global ocean basin) of heat transfer annually from the global ocean to the atmosphere during storm passage (0–3 days) on a global scale. Of this total, 1.05 ± 0.20 PW (4.80 ± 0.85 W m −2 ) is caused by Tropical storms/Tropical depressions (TS/TD) and 0.82 ± 0.21 PW (6.25 ± 1.5 W m −2 ) is caused by hurricanes. Our findings indicate that ocean heat loss by TCs may be a substantial missing piece of the global ocean heat budget. Net changes in OHC after storm passage is estimated by analyzing the temperature anomalies during wake recovery following storm events (4–20 days after storm passage) relative to pre-storm conditions. Results indicate the global ocean experiences a 0.75 ± 0.25 PW (5.98 ± 2.1W m −2 ) net heat gain annually for hurricanes. In contrast, under TS/TD conditions, ocean experiences 0.41 ± 0.21 PW (1.90 ± 0.96 W m −2 ) net ocean heat loss, suggesting the overall oceanic thermal response is particularly sensitive to the intensity of the event. The net ocean heat uptake caused by all storms is 0.34 PW.

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: 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: 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: Dynamics of turbulent western boundary currents at low latitude in a shallow water model Ocean Science Discussions, 11, 2461-2493, 2014 Author(s): C. Q. C. Akuetevi and A. Wirth The dynamics of low latitude turbulent western boundary currents crossing the equator is considered using numerical results from integrations of a reduced gravity shallow-water model. For viscosity values of 1000 m 2 s −1 and more, the boundary layer dynamics compares well to the analytical Munk-layer solution. When the viscosity is reduced, the boundary layer becomes turbulent and coherent structures in form of anticyclonic eddies, bursts (violent detachments of the viscous sub-layer) and dipoles appear. Three distinct boundary layers emerge, the viscous sub-layer, the advective boundary layer and the extended boundary layer. The first is characterized by a dominant vorticity balance between the viscous transport and the advective transport of vorticity. The second by a balance between the advection of planetary vorticity and the advective transport of relative vorticity. The extended boundary layer is the area to which turbulent motion from the boundary extends. The scaling of the three boundary layer thicknesses with viscosity is evaluated. Characteristic scales of the dynamics and dissipation are determined. A pragmatic approach to determine the eddy viscosity diagnostically for coarse resolution numerical models is proposed.

Description: Calculating the water and heat balances of the Eastern Mediterranean basin using ocean modelling and available meteorological, hydrological, and ocean data Ocean Science Discussions, 8, 1301-1338, 2011 Author(s): M. Shaltout and A. Omstedt This paper analyses the Eastern Mediterranean water and heat balances over a 52-yr period. The modelling uses a process-oriented approach resolving the one-dimensional equations of momentum, heat, and salt conservation, with turbulence modelled using a two-equation model. The exchange through the Sicily Channel connecting the Eastern and Western basins is calculated from satellite altimeter data. The results illustrates that calculated surface temperature and salinity follow the reanalysed data well and with biases of −0.4 °C and −0.004, respectively. Monthly and yearly temperature and salinity cycles are also satisfactory simulated. Reanalysed data and calculated water mass structure and heat balance components are in good agreement, indicating that the air-sea interaction and the turbulent mixing are realistically simulated. The study illustrates that the water balance in the Eastern Mediterranean basin is controlled by the difference between inflows/outflows through the Sicily Channel and by the net precipitation rates. The heat balance is controlled by the heat loss from the water surface, sun radiation into the sea, and heat flow through the Sicily Channel, the first two displaying both climate trends. An annual net heat loss of approximately 9 W m −2 was balanced by net heat in flow through the Sicily Channel.

Description: Extraction of spatial-temporal rules from mesoscale eddies in the South China Sea Based on rough set theory Ocean Science Discussions, 8, 1261-1300, 2011 Author(s): Y. Du, X. Fan, Z. He, F. Su, C. Zhou, H. Mao, and D. Wang In this paper, a rough set theory is introduced to represent spatial-temporal relationships and extract the corresponding rules from typical mesoscale-eddy states in the South China Sea (SCS). Three decision attributes are adopted in this study, which make the approach flexible in retrieving spatial-temporal rules with different features. Spatial-temporal rules of typical states in the SCS are extracted as three decision attributes, which then are confirmed by the previous works. The results demonstrate that this approach is effective in extracting spatial-temporal rules from typical mesoscale-eddy states, and therefore provides a powerful approach to forecasts in the future. Spatial-temporal rules in the SCS indicate that warm eddies following the rules are generally in the southeastern and central SCS around 2000 m isobaths in winter. Their intensity and vorticity are weaker than those of cold eddies. They usually move a shorter distance. By contrast, cold eddies are in 2000 m-deeper regions of the southwestern and northeastern SCS in spring and fall. Their intensity and vorticity are strong. Usually they move a long distance. In winter, a few rules are followed by cold eddies in the northern tip of the basin and southwest of Taiwan Island rather than warm eddies, indicating cold eddies may be well-regulated in the region. Several warm-eddy rules are achieved west of Luzon Island, indicating warm eddies may be well-regulated in the region as well. Otherwise, warm and cold eddies are distributed not only in the jet flow off southern Vietnam induced by intraseasonal wind stress in summer-fall, but also in the northern shallow water, which should be a focus of future study.

Description: Tracer distribution in the Pacific Ocean following a release off Japan – what does an oceanic general circulation model tell us? Ocean Science Discussions, 8, 1441-1466, 2011 Author(s): H. Dietze and I. Kriest In the aftermath of an earthquake and tsunami on 11 March 2011 considerable amounts of radioactive materials were accidentally released into the sea off Fukushima-Daiichi, Japan. This study uses a three-dimensional eddy-resolving oceanic general circulation model to explore potential pathways of a tracer, similar to 137 Cs, from the coast to the open ocean. Results indicate that enhanced concentrations meet a receding spring bloom offshore and that the area of enhanced concentrations offshore is strongly determined by surface mixed layer dynamics. However, huge uncertainties remain. Among them are the realism of the simulated cross-shelf transport and apparently inconsistent estimates of the particle reactivity of 137 Cs which are discussed in a brief literature review. We argue that a comprehensive set of 137 Cs measurements, including sites offshore, could be a unique opportunity to both evaluate and advance the evaluation of oceanic general circulation models.

Description: Usefulness of high resolution coastal models for operational oil spill forecast: the Full City accident Ocean Science Discussions, 8, 1467-1504, 2011 Author(s): G. Broström, A. Carrasco, L. R. Hole, S. Dick, F. Janssen, J. Mattsson, and S. Berger Oil spill modeling is considered to be an important decision support system (DeSS) useful for remedial action in case of accidents, as well as for designing the environmental monitoring system that is frequently set up after major accidents. Many accidents take place in coastal areas implying that low resolution basin scale ocean models is of limited use for predicting the trajectories of an oil spill. In this study, we target the oil spill in connection with the Full City accident on the Norwegian south coast and compare three different oil spill models for the area. The result of the analysis is that all models do a satisfactory job. The "standard" operational model for the area is shown to have severe flaws but including an analysis based on a higher resolution model (1.5 km resolution) for the area the model system show results that compare well with observations. The study also shows that an ensemble using three different models is useful when predicting/analyzing oil spill in coastal areas.

Description: Impact of combining GRACE and GOCE gravity data on ocean circulation estimates Ocean Science Discussions, 8, 1535-1573, 2011 Author(s): T. Janjić, J. Schröter, R. Savcenko, W. Bosch, A. Albertella, R. Rummel, and O. Klatt In this work we examine the impact of assimilation of multi-mission-altimeter data and the GRACE/GOCE gravity fields into the finite element ocean model (FEOM), with the focus on the Southern Ocean circulation. In order to do so, we use the geodetic approach for obtaining the dynamical ocean topography (DOT), that combines the multi-mission-altimeter data and the GRACE/GOCE gravity fields, and requires that both fields be spectrally consistent. The spectral consistency is achieved by filtering of the sea surface height and the geoid using profile approach. Combining the GRACE and GOCE data, a considerably shorter filter length resolving more DOT details can be used. In order to specify the spectrally consistent geodetic DOT we applied the Jekeli-Wahr filter corresponding to 241 km, 121 km, 97 km and 81 km halfwidths for the GRACE/GOCE based gravity field model GOCO01S and to the sea surface. More realistic features of the ocean assimilation were obtained in the Weddel gyre area due to increased resolution of the data fields, particularly for temperature field at the 800 m depth compared to Argo data.

Description: Wind forcing of salinity anomalies in the Denmark Strait overflow Ocean Science Discussions, 8, 1403-1440, 2011 Author(s): S. Hall, S. R. Dye, K. J. Heywood, and M. R. Wadley The overflow of dense water from the Nordic Seas to the North Atlantic through Denmark Strait is an important part of the global thermohaline circulation. The salinity of the overflow plume has been measured by an array of current meters across the continental slope off the coast of Angmagssalik, southeast Greenland since September 1998. During 2004 the salinity of the overflow plume changed dramatically, with the entire width of the array (70 km) freshening between January 2004 and July 2004, with a significant negative salinity anomaly of about 0.06 in May. The event in May represents a fresh anomaly of over 3 standard deviations from the mean since recording began in 1998. We show that the OCCAM 1/12° Ocean General Circulation Model not only reproduces the 2004 freshening event ( r =0.96, p

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: Mesoscale variability of water masses in the Arabian Sea as revealed by ARGO floats Ocean Science Discussions, 8, 1369-1402, 2011 Author(s): X. Carton and P. L'Hegaret By analysing ARGO float data over the last four years, some aspects of the mesoscale variability of water masses in the Arabian Sea are described. The Red Sea Water outflow is strong in the Southwestern Gulf of Aden, in particular when a cyclonic gyre predominates in this region. Salinities of 36.5 and temperatures of 16 °C are found there between 600 and 1000 m depths. The Red Sea Water is more dilute in the eastern part of the Gulf, and fragments of this water mass can be advected offshore across the gulf or towards its northern coast by the regional gyres. The Red Sea Water outflow is also detected along the northeastern coast of Socotra, and fragments of RSW are found between one and three degrees of latitude north of this island. In the whole Gulf of Aden, the correlation between the deep motions of the floats and the SSH measured by altimetry is strong, at regional scale. The finer scale details of the float trajectories are more often related to the anomalous water masses that they encounter. The Persian Gulf Water (PGW) is found in the float profiles near Ras ash Sharbatat (near 57° E, 18° N), again with 36.5 in salinity and about 18–19 °C in temperature. These observations were achieved in winter when the southwestward monsoon currents can advect PGW along the South Arabian coast. Fragments of PGW are found in the Arabian Sea between 18 and 20° N and 63 and 65° E, showing that this water mass can escape the Gulf of Oman southeastward, in particular during summer.

Description: Using dissolved oxygen concentrations to determine mixed layer depths in the Bellingshausen Sea Ocean Science Discussions, 8, 1505-1533, 2011 Author(s): K. Castro-Morales and J. Kaiser Concentrations of oxygen (O 2 ) and other dissolved gases in the oceanic mixed layer are often used to calculate air-sea gas exchange fluxes; for example, in the context of net and gross biological production estimates. The mixed layer depth ( z mix ) may be defined using criteria based on temperature or density differences to a reference depth near the ocean surface. However, temperature criteria fail in regions with strong haloclines such as the Southern Ocean where heat, freshwater and momentum fluxes interact to establish mixed layers. Moreover, the time scales of air-sea exchange differ for gases and heat, so that z mix defined using O 2 may be different to z mix defined using temperature or density. Here, we propose to define an O 2 -based mixed layer depth, z mix (O 2 ), as the depth where the relative difference between the O 2 concentration and a reference value at a depth equivalent to 10 dbar equals 0.5 %. This definition was established by numerical analysis of O 2 profiles in coastal areas of the Southern Ocean and corroborated by visual inspection. Comparisons of z mix (O 2 ) with z mix based on potential temperature differences, i.e. z mix (Δ θ = 0.2 °C) and z mix (Δ θ = 0.5 °C), and potential density differences, i.e. z mix (Δ σ θ = 0.03 kg m −3 ) and z mix (Δ σ θ = 0.125 kg m −3 ), showed that z mix (O 2 ) closely follows z mix (Δ σ θ = 0.03 kg m −3 ). Further comparisons with published z mix climatologies and z mix derived from World Ocean Atlas 2005 data were also performed. To establish z mix for use with biological production estimates in the absence of O 2 profiles, we suggest using z mix (Δ σ θ = 0.03 kg m −3 ), which is also the basis for the climatology by de Boyer Montégut et al. (2004).

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: N/P ratio of nutrient uptake in the Baltic Sea Ocean Science Discussions, 8, 1233-1259, 2011 Author(s): Z. Wan, L. Jonasson, and H. Bi The N/P ratio of nutrient uptake, i.e., the ratio of dissolved inorganic nitrogen (DIN) to dissolved inorganic phosphorus (DIP) taken by primary producers, varies in different basins and in different seasons in the Baltic Sea. The N/P ratio of nutrient alteration fore and after spring blooms is not same as the N/P ratio of nutrient uptake, but the former can be regarded as an indicator for the later in the Baltic Sea. Based on the observed N/P ratio of nutrient alteration, we hypothesize a non-Redfield N/P ratio of nutrient uptake. The 3D-ecosystem model ERGOM coupled with the circulation model DMI-BSHcmod was used to test the hypothesis. When the Redfield ratio was used in the model, the DIP surplus after spring blooms was too high and resulted in the overly growth of cyanobacteria and too much nitrogen fixation. When the non-Redfield ratio was used in the model, the corresponding problem tended to disappear. In summary, we show that: (1) the Redfield N/P ratio of nutrient uptake in the Baltic Sea tends to be too high; (2) a lower N/P ratio 10:1 appears to work better than the Redfield value; and (3) the N/P ratio of nutrient uptake in the Baltic Proper during spring blooms is around 6:1.

Description: Spectrophotometric high-precision seawater pH determination for use in underway measuring systems Ocean Science Discussions, 8, 1339-1367, 2011 Author(s): S. Aßmann, C. Frank, and A. Kötzinger Autonomous sensors are required for a comprehensive documentation of the changes in the marine carbon system and thus to differentiate between its natural variability and anthropogenic impacts. Spectrophotometric determination of pH – a key variable of the seawater carbon system – is particularly suited to achieve precise and drift-free measurements. However, available spectrophotometric instruments are not suitable for integration into automated measurement systems (e.g. FerryBox) since they do not meet the major requirements of reliability, stability, robustness and moderate cost. Here we report on the development and testing of a new indicator-based pH sensor that meets all of these requirements. This sensor can withstand the rough conditions during long-term deployments on ships of opportunities and is applicable on the open ocean as well as in coastal waters with complex background and highly variable conditions. The sensor uses a high resolution CCD spectrometer as detector connected via optical fibers to a custom-made cuvette designed to reduce the impact of air bubbles. The sample temperature can be precisely adjusted (25 °C ± 0.006 °C) using computer-controlled power supplies and Peltier elements thus avoiding the widely used water bath. The overall setup achieves a measurement frequency of 1 min −1 with a precision of ± 0.0007 pH units and an average offset of +0.0018 pH units to a pH reference during shipboard operation. Application of this sensor allows monitoring of seawater pH in autonomous underway systems, providing a key variable for characterization and understanding the marine carbon system.

Description: A pre-operational 3-D variational data assimilation system in the North/Baltic Sea Ocean Science Discussions, 8, 1131-1160, 2011 Author(s): S. Y. Zhuang, W. W. Fu, and J. She This paper describes the implementation and evaluation of a pre-operational three dimensional variational (3DVAR) data assimilation system for the North/Baltic Sea. The univariate analysis for temperature and salinity is applied in a 3DVAR scheme in which the horizontal component of the background error covariance is modeled by an isotropic recursive filter (IRF) and the vertical component is represented by dominant Empirical Orthogonal Functions (EOFs) of the background error. Observations of temperature and salinity ( T / S ) profiles in the North/Baltic Sea are assimilated in the year of 2005. Effectiveness of the data assimilation scheme is assessed by comparison with the control run that no assimilation is done. The statistical analysis indicates that the model simulation is significantly improved with the 3DVAR scheme. On average, the root mean square error (RMSE) of temperature and salinity is reduced by 0.2 °C and 0.25 psu in the North/Baltic Sea. In addition, the bias of temperature and salinity is also decreased by 0.1 °C and 0.2 psu, respectively. Starting from an analyzed initial state, one month simulation without assimilation is carried out with the aim of examining the persistence of the initial impact. It is shown that the assimilated initial state can impact the model simulation for nearly two weeks. The influence on salinity is more pronounced than temperature.

Description: Numerical simulation and decomposition of kinetic energies in the Central Mediterranean Sea: insight on mesoscale circulation and energy conversion Ocean Science Discussions, 8, 1161-1214, 2011 Author(s): R. Sorgente, A. Olita, P. Oddo, L. Fazioli, and A. Ribotti The spatial and temporal variability of eddy and mean kinetic energy of the Central Mediterranean Sea has been investigated, from January 2008 to December 2010, by mean of a numerical simulation mainly to quantify the mesoscale dynamics and their relationships with physical forcing. In order to understand the energy redistribution processes, the baroclinic energy conversion has been analysed, suggesting hypotheses about the drivers of the mesoscale activity in this area. The ocean model used is based on the Princeton Ocean Model implemented at 1/32° horizontal resolution. Surface momentum and buoyancy fluxes are interactively computed by mean of standard bulk formulae using predicted model Sea Surface Temperature and atmospheric variables provided by the European Centre for Medium Range Weather Forecast operational analyses. At its lateral boundaries the model is one-way nested within the Mediterranean Forecasting System operational products. The model domain has been subdivided in four sub-regions: Sardinia channel and southern Tyrrhenian Sea, Sicily channel, eastern Tunisian shelf and Libyan Sea. Temporal evolution of eddy and mean kinetic energy has been analysed, on each of the four sub-regions composing the model domain, showing different behaviours. On annual scales and within the first 5 m depth, the eddy kinetic energy represents approximately the 60 % of the total kinetic energy over the whole domain, confirming the strong mesoscale nature of the surface current flows in this area. The analyses show that the model well reproduces the path and the temporal behaviour of the main known sub-basin circulation features. New mesoscale structures have been also identified, from numerical results and direct observations, for the first time as the Pantelleria Vortex and the Medina Gyre. The classical the kinetic energy decomposition (eddy and mean) allowed to depict and to quantify the stable and fluctuating parts of the circulation in the region, and to differentiate the four sub-regions as function of relative and absolute strength of the mesoscale activity. Furthermore the Baroclinic Energy Conversion term shows that in the Sardinia Channel the mesoscale activity, due to baroclinic instabilities, is significantly larger than in the other sub-regions, while a negative sign of the energy conversion, meaning a transfer of energy from the Eddy Kinetic Energy to the Eddy Available Potential Energy, has been recorded only for the surface layers of the Sicily Channel during summer.

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: Modelling of the anthropogenic tritium transient and its decay product helium-3 in the Mediterranean Sea using a high-resolution regional model Ocean Science Discussions, 11, 2691-2732, 2014 Author(s): M. Ayache, J.-C. Dutay, P. Jean-Baptiste, K. Beranger, T. Arsouze, J. Beuvier, J. Palmieri, B. Le-vu, and W. Roether This numerical study provides the first simulation of the anthropogenic tritium invasion and its decay product helium-3 ( 3 He) in the Mediterranean Sea. The simulation covers the entire tritium ( 3 H) transient generated by the atmospheric nuclear-weapon tests performed in the 1950s and early 1960s and run till 2011. Tritium, helium-3 and their derived age estimates are particularly suitable for studying intermediate and deep-water ventilation and spreading of water masses at intermediate/deep levels. The simulation is made using a high resolution regional model NEMO-MED12 forced at the surface with prescribed tritium evolution derived from observations. The simulation is compared to measurements of tritium and helium-3 performed along large-scale transects in the Mediterranean Sea during the last few decades on cruises of Meteor M5/6, M31/1, M44/4, M51/2, M84/3, and Poseidon 234. The results show that the input function used for the tritium, generates a realistic distribution of the main hydrographic features of the Mediterranean Sea circulation. In the eastern basin, the results highlight the weak formation of Adriatic Deep Water in the model, which explains its weak contribution to the Eastern Mediterranean Deep Water in the Ionian sub-basin. It produces a realistic representation of the Eastern Mediterranean Transient signal, simulating a deep-water formation in the Aegean sub-basin at the beginning of the 1993, with a realistic timing of deep-water renewal in the eastern basin. In the western basin, the unusual intense deep convection event of winter 2005 in the Gulf of Lions during the Western Mediterranean Transition is simulated. However the spreading of the recently ventilated deep water toward the South is too weak. The ventilation and spreading of the Levantine Intermediate Water from the eastern basin toward the western basin is simulated with realistic tracer-age distribution compared to observation-based estimates.

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: 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: 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: The Aegean Sea marine security decision support system Ocean Science Discussions, 8, 1025-1053, 2011 Author(s): L. Perivoliotis, G. Krokos, K. Nittis, and G. Korres As part of the integrated ECOOP (European Coastal Sea Operational observing and Forecasting System) project, HCMR upgraded the already existing standalone Oil Spill Forecasting System for the Aegean Sea, initially developed for the Greek Operational Oceanography System (POSEIDON), into an active element of the European Decision Support System (EuroDeSS). The system is accessible through a user friendly web interface where the case scenarios can be fed into the oil spill drift model component, while the synthetic output contains detailed information about the distribution of oil spill particles and the oil spill budget and it is provided both in text based ECOOP common output format and as a series of sequential graphics. The main development steps that were necessary for this transition were the modification of the forcing input data module in order to allow the import of other system products which are usually provided in standard formats such as NetCDF and the transformation of the model's calculation routines to allow use of current, density and diffusivities data in z instead of sigma coordinates. During the implementation of the Aegean DeSS, the system was used in operational mode in order support the Greek marine authorities in handling a real accident that took place in North Aegean area. Furthermore, the introduction of common input and output files by all the partners of EuroDeSS extended the system's interoperability thus facilitating data exchanges and comparison experiments.

Description: How well can we derive Global Ocean Indicators from Argo data? Ocean Science Discussions, 8, 999-1024, 2011 Author(s): K. von Schuckmann and P.-Y. Le Traon Argo deployments began in the year 2000 and by November 2007 the array was 100 % complete, covering the global ocean from the surface down to 2000 m depth. In this study, Argo temperature and salinity measurements during the period 2005 to 2010 are used to develop a revised estimation of Global Ocean Indicators (GOIs) such as heat content variability, freshwater content and steric height. These revised indices are based on a simple box averaging scheme using a weighted mean. They include a proper estimation of the errors due to data handling methods and climatology uncertainties. A global ocean heat content change (OHC) trend of 0.55 ± 0.1 W m −2 is estimated over the time period 2005–2010. Similarly, a global steric sea level (GSSL) rise of 0.69 ± 0.14 mm yr −1 is observed. The global ocean freshwater content (OFC) trend is barely significant. Results show that there is significant interannual variability at global scale, especially for global OFC. Annual mean GOIs from the today's Argo samling can be derived with an accuracy of ±0.10 cm for GSSL, ±0.21 × 10 8 J m −2 for global OHC, and ±700 km 3 for global OFC. Long-term trends (15 yr) of GOIs based on the complete Argo sampling (10–1500 m depth) can be performed with an accuracy of about ±0.03 mm yr −1 for steric rise, ±0.02 W m −2 for ocean warming and ±20 km 3 yr −1 for global OFC trends – under the assumption that no systematic errors remain in the observing system.

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.

Description: Present dynamics and future prognosis of a slowly surging glacier The Cryosphere, 5, 299-313, 2011 Author(s): G. E. Flowers, N. Roux, S. Pimentel, and C. G. Schoof Glacier surges are a well-known example of an internal dynamic oscillation whose occurrence is not a direct response to the external climate forcing, but whose character (i.e. period, amplitude, mechanism) may depend on the glacier's environmental or climate setting. We examine the dynamics of a small (∼5 km 2 ) valley glacier in Yukon, Canada, where two previous surges have been photographically documented and an unusually slow surge is currently underway. To characterize the dynamics of the present surge, and to speculate on the future of this glacier, we employ a higher-order flowband model of ice dynamics with a regularized Coulomb-friction sliding law in both diagnostic and prognostic simulations. Diagnostic (force balance) calculations capture the measured ice-surface velocity profile only when non-zero basal water pressures are prescribed over the central region of the glacier, coincident with where evidence of the surge has been identified. This leads to sliding accounting for 50–100% of the total surface motion in this region. Prognostic simulations, where the glacier geometry evolves in response to a prescribed surface mass balance, reveal a significant role played by a bedrock ridge beneath the current equilibrium line of the glacier. Ice thickening occurs above the ridge in our simulations, until the net mass balance reaches sufficiently negative values. We suggest that the bedrock ridge may contribute to the propensity for surges in this glacier by promoting the development of the reservoir area during quiescence, and may permit surges to occur under more negative balance conditions than would otherwise be possible. Collectively, these results corroborate our interpretation of the current glacier flow regime as indicative of a slow surge that has been ongoing for some time, and support a relationship between surge incidence or character and the net mass balance. Our results also highlight the importance of glacier bed topography in controlling ice dynamics, as observed in many other glacier systems.

Description: Silicon pool dynamics and biogenic silica export in the Southern Ocean, inferred from Si-isotopes Ocean Science Discussions, 8, 639-674, 2011 Author(s): F. Fripiat, A.-J. Cavagna, F. Dehairs, S. Speich, L. André, and D. Cardinal Water column silicon isotopic signatures (δ 30 Si) of silicic acid (Si(OH) 4 ) in the Southern Ocean were measured along a meridional transect from South Africa (Subtropical Zone) down to 57° S (northern Weddell Gyre). These data are the first reported for a summer transect across the whole Antarctic Circumpolar Current (ACC). δ 30 Si variations are large in the upper 1000 m, reflecting the effect of the silica pump superimposed upon meridional transfer across the ACC: the transport of Antarctic surface waters northward by a net Ekman drift and their convergence and mixing with warmer upper-ocean Si-depleted waters to the north. Using Si isotopic signatures, we determined different mixing interfaces between ACC water masses: the Antarctic Surface Water (AASW), the Antarctic Intermediate Water (AAIW), and the thermoclines in the low latitude areas. The residual silicic acid concentrations of end-members control the δ 30 Si alteration of the mixing products. With the exception of AASW, all mixing interfaces have a highly Si-depleted mixed layer end-member. These processes deplete the silicic acid AASW concentration across the different interfaces northward without significantly changing the AASW δ 30 Si. By comparing our new results with a previous study in the Australian sector we show that during the circumpolar transport of the ACC eastward, there is a slight but significant Si-isotopic lightening of the silicic acid pools from the Atlantic to the Australian sectors. This results either from the dissolution of biogenic silica in the deeper layers and/or from an isopycnal mixing with the deep water masses in the different oceanic basins: North Atlantic Deep Water in the Atlantic, and Indian Ocean deep water in the Indo-Australian sector. This eastward lightening is further transmitted to the subsurface waters, representing mixing interfaces between the surface and deeper layers. Using the Si-isotopic constraint, we estimate for the Greenwich Meridian a net biogenic silica production which should be representative of the annual export, at 4.5 ± 1.1 and 1.5 ± 0.4 mol Si m −2 for the Antarctic Zone and Polar Front Zone, respectively, in agreement with previous estimations. The summertime Si-supply into the mixed layer via vertical mixing was also assessed at 1.5 ± 0.4 and 0.1 ± 0.5 mol Si m −2 , respectively.

Description: Quality control of automated hyperspectral remote sensing measurements from a seaborne platform Ocean Science Discussions, 8, 613-638, 2011 Author(s): S. P. Garaba, M. R. Wernand, and O. Zielinski In this study four data quality flags are presented for automated and unmanned above-water hyperspectral optical measurements collected underway in the North Sea, The Minch, Irish Sea and Celtic Sea in April/May 2009. Coincident to these optical measurements a DualDome D12 (Mobotix, Germany) camera system was used to capture sea surface and sky images. The first three flags are based on meteorological conditions, to select erroneous incoming solar irradiance ( E S ) taken during dusk, dawn, before significant incoming solar radiation could be detected or under rainfall. Furthermore, the relative azimuthal angle of the optical sensors to the sun is used to identify possible sunglint free sea surface zones. A total of 629 spectra remained after applying the meteorological masks (first three flags). Based on this dataset, a fourth flag for sunglint was generated by analysing and evaluating water leaving radiance ( L W ) and remote sensing reflectance ( R RS ) spectral behaviour in the presence and absence of sunglint salient in the simultaneously available sea surface images. Spectra conditions satisfying "mean L W (700–950 nm) 〈 2 mW m −2 nm −1 Sr −1 " or alternatively "minimum R RS (700–950 nm) 〈 0.010 Sr −1 ", mask the most measurements affected by sunglint, providing efficient flagging of sunglint in automated quality control. It is confirmed that valid optical measurements can be performed 0° ≤ Φ ≤ 360° although 90° ≤ Φ ≤ 135° is recommended.

Description: Numerical modelling of POC yearly dynamics in the southern Baltic under variable scenarios of nutrients, light and temperature Ocean Science Discussions, 8, 675-700, 2011 Author(s): L. Dzierzbicka-Glowacka, K. Kuliński, A. Maciejewska, J. Jakacki, and J. Pempkowiak This paper presents various scenarios of the particulate organic carbon (POC) in the southern Baltic Sea. The study is based on a one-dimensional Particulate Organic Carbon model (1-D POC). Mathematically, the pelagic variables of 1-D POC model are described by a second-order partial differential equations of the diffusion type with biogeochemical sources and sinks. The POC concentration is determined as the sum of phytoplankton, zooplankton and dead organic matter (detritus) concentrations. The temporal changes in the phytoplankton biomass are caused by primary production, mortality, grazing by zooplankton and sinking. The zooplankton biomass is affected by ingestion, excretion, faecal production, mortality, and carnivorous grazing. The changes in the pelagic detritus concentration are determined by input of: dead phytoplankton and zooplankton, natural mortality of predators, faecal pellets, and sinks: sedimentation, zooplankton grazing and biochemical decomposition. The 1-D POC model was used to simulate temporal dynamics of POC in the southern Baltic Sea (Gdansk Deep, Bornholm Deep and Gotland Deep) under scenarios characterized by different temperature, nutrients and light. Daily, monthly, seasonal and annual variabilities of POC in the upper water layer are presented for the different scenarios. The starting-point of the numerical simulations was assumed as average values of the investigated pelagic variables for 1965–1998 period. Two- to three-fold increases of POC concentrations in late spring were revealed as well as the shift towards postponed maximum POC concentration. It is speculated that, due to POC increase, oxygenation of under-halocline water layer will decrease, while supply of food to organisms from higher trophic level should increase.

Description: A computational method for determining XBT depths Ocean Science Discussions, 8, 1777-1802, 2011 Author(s): J. Stark, J. Gorman, M. Hennessey, F. Reseghetti, J. Willis, J. Lyman, and J. Abraham A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method combines a forward-stepping calculation which incorporates all of the forces on the XBT devices during their descent. Of particular note are drag forces which are calculated using a new drag coefficient expression. That expression, obtained entirely from computational fluid dynamic modeling, accounts for local variations in the ocean environment. Consequently, the method allows for accurate determination of depths for any local temperature environment. The results, which are entirely based on numerical simulation, are compared with an experimental descent of an LM-Sippican T-5 XBT. It is found that the calculated depths differ by less than 3 % from depth estimates using the industry standard FRE. Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid patterns along the outer surface of the probe as well as in the interior channel. The simulations take account of complex flow phenomena such as laminar-turbulent transition and flow separation.

Description: In-situ multispectral and bathymetric measurements over a supraglacial lake in western Greenland using a remotely controlled watercraft The Cryosphere, 5, 445-452, 2011 Author(s): M. Tedesco and N. Steiner Supraglacial lakes form from meltwater on the Greenland ice sheet in topographic depressions on the surface, affecting both surface and sub-glacial processes. As the reflectance in the visible and near-infrared regions of a column of water is modulated by its height, retrieval techniques using spaceborne remote sensing data (e.g. Landsat, MODIS) have been proposed in the literature for the detection of lakes and estimation of their volume. These techniques require basic assumptions on the spectral properties of the water as well as the bottom of the lake, among other things. In this study, we report results obtained from the analysis of concurrent in-situ multi-spectral and depth measurements collected over a supraglacial lake during early July 2010 in West Greenland (Lake Olivia, 69°36'35" N, 49°29'40" W) and aim to assess some of the underlying hypotheses in remote sensing based bathymetric approaches. In particular, we focus our attention on the analysis of the lake bottom albedo and of the water attenuation coefficient. The analysis of in-situ data (collected by means of a remotely controlled boat equipped with a GPS, a sonar and a spectrometer) highlights the exponential trend of the water-leaving reflectance with lake depth. The values of the attenuation factor obtained from in-situ data are compared with those computed using approaches proposed in the literature. Also, the values of the lake bottom albedo from in-situ measurements are compared with those obtained from the analysis of reflectance of shallow waters. Finally, we quantify the error between in-situ measured and satellite-estimated lake depth values for the lake under study.

Description: A vertical-mode decomposition to investigate low-frequency internal motion across the Atlantic at 26° N Ocean Science Discussions, 8, 2047-2100, 2011 Author(s): Z. B. Szuts, J. R. Blundell, M. P. Chidichimo, and J. Marotzke Hydrographic data from full-depth moorings maintained by the RAPID/MOCHA project that span the Atlantic at 26° N are decomposed into vertical modes, in order to give a dynamical framework for interpreting the observed fluctuations. Vertical modes at each mooring are fit to pressure perturbations using a Gauss-Markov inversion. Away from boundaries, the vertical structure is almost entirely described by the first baroclinic mode, as confirmed by high correlation between the original signal and reconstructions using only the first baroclinic mode. These first baroclinic motions are also highly coherent with altimetric sea surface height (SSH). On both the western and eastern boundaries, however, the decomposition contains significant variance at higher modes, and there is a corresponding decrease in the agreement between SSH and either the original signal or the first baroclinic mode reconstruction. At the boundaries, the transport fluctuations described by the first baroclinic mode represent less than 10% of the variance of the full transport signal. At the eastern boundary, a linear combination of many baroclinic modes is required to explain the observed vertical density profile of the seasonal cycle, a result that is consistent with the oceanic response to wind-forcing not propagating far from the eastern boundary.

Description: Temperature variability and offset in steep alpine rock and ice faces The Cryosphere, 5, 977-988, 2011 Author(s): A. Hasler, S. Gruber, and W. Haeberli The thermal condition of high-alpine mountain flanks can be an important determinant of climate change impact on slope stability and correspondingly down-slope hazard regimes. In this study we analyze time-series from 17 shallow temperature-depth profiles at two field sites in steep bedrock and ice. Extending earlier studies that revealed the topographic variations in temperatures, we demonstrate considerable differences of annual mean temperatures for variable surface characteristics and depths within the measured profiles. This implies that measurements and model related to compact and near-vertical bedrock temperatures may deviate considerably from conditions in the majority of bedrock slopes in mountain ranges that are usually non-vertical and fractured. For radiation-exposed faces mean annual temperatures at depth are up to 3 °C lower and permafrost is likely to exist at lower elevations than reflected by estimates based on near-vertical homogeneous cases. Retention of a thin snow cover and ventilation effects in open clefts are most likely responsible for this cooling. The measurements presented or similar data could be used in the future to support the development and testing of models related to the thermal effect of snow-cover and fractures in steep bedrock.

Description: On the influence of model physics on simulations of Arctic and Antarctic sea ice The Cryosphere, 5, 687-699, 2011 Author(s): F. Massonnet, T. Fichefet, H. Goosse, M. Vancoppenolle, P. Mathiot, and C. König Beatty Two hindcast (1983–2007) simulations are performed with the global, ocean-sea ice models NEMO-LIM2 and NEMO-LIM3 driven by atmospheric reanalyses and climatologies. The two simulations differ only in their sea ice component, while all other elements of experimental design (resolution, initial conditions, atmospheric forcing) are kept identical. The main differences in the sea ice models lie in the formulation of the subgrid-scale ice thickness distribution, of the thermodynamic processes, of the sea ice salinity and of the sea ice rheology. To assess the differences in model skill over the period of investigation, we develop a set of metrics for both hemispheres, comparing the main sea ice variables (concentration, thickness and drift) to available observations and focusing on both mean state and seasonal to interannual variability. Based upon these metrics, we discuss the physical processes potentially responsible for the differences in model skill. In particular, we suggest that (i) a detailed representation of the ice thickness distribution increases the seasonal to interannual variability of ice extent, with spectacular improvement for the simulation of the recent observed summer Arctic sea ice retreats, (ii) the elastic-viscous-plastic rheology enhances the response of ice to wind stress, compared to the classical viscous-plastic approach, (iii) the grid formulation and the air-sea ice drag coefficient affect the simulated ice export through Fram Strait and the ice accumulation along the Canadian Archipelago, and (iv) both models show less skill in the Southern Ocean, probably due to the low quality of the reanalyses in this region and to the absence of important small-scale oceanic processes at the models' resolution (~1°).

Description: Arctic Ocean circulation and variability – advection and external forcing encounter constraints and local processes Ocean Science Discussions, 8, 2313-2376, 2011 Author(s): B. Rudels The first hydrographic data from the Arctic Ocean, the section from the Laptev Sea to the passage between Greenland and Svalbard obtained by Nansen on the drift by Fram 1893–1896, aptly illustrate the main features of Arctic Ocean oceanography and indicate possible processes active in transforming the water masses in the Arctic Ocean. Many, perhaps most, of these processes were identified already by Nansen, who put his mark on almost all subsequent research in the Arctic Ocean. Here we shall revisit some key questions and follow how our understanding has evolved from the early 20th century to present. What questions, if any, can now be regarded as solved and which remain still open? Five different but connected topics will be discussed: (1) The low salinity surface layer and the storage and export of freshwater. (2) The vertical heat transfer from the Atlantic water to sea ice and to the atmosphere. (3) The circulation and mixing of the two Atlantic inflow branches. (4) The formation and circulation of deep and bottom waters in the Arctic Ocean. (5) The exchanges through Fram Strait. Foci will be on the potential effects of increased freshwater input and reduced sea ice export on the freshwater storage and residence time in the Arctic Ocean, on the deep waters of the Makarov Basin and on the circulation and relative importance of the two inflows, over the Barents Sea and through Fram Strait, for the distribution of heat in the intermediate layers of the Arctic Ocean.

Description: Glacier contribution to streamflow in two headwaters of the Huasco River, Dry Andes of Chile The Cryosphere, 5, 1099-1113, 2011 Author(s): S. Gascoin, C. Kinnard, R. Ponce, S. Lhermitte, S. MacDonell, and A. Rabatel Quantitative assessment of glacier contribution to present-day streamflow is a prerequisite to the anticipation of climate change impact on water resources in the Dry Andes. In this paper we focus on two glaciated headwater catchments of the Huasco Basin (Chile, 29° S). The combination of glacier monitoring data for five glaciers (Toro 1, Toro 2, Esperanza, Guanaco, Estrecho and Ortigas) with five automatic streamflow records at sites with glacier coverage of 0.4 to 11 % allows the estimation of the mean annual glacier contribution to discharge between 2003/2004 and 2007/2008 hydrological years. In addition, direct manual measurements of glacier runoff were conducted in summer at the snouts of four glaciers, which provide the instantaneous contribution of glacier meltwater to stream runoff during summer. The results show that the mean annual glacier contribution to streamflow ranges between 3.3 and 23 %, which is greater than the glaciated fraction of the catchments. We argue that glacier contribution is partly enhanced by the effect of snowdrift from the non-glacier area to the glacier surface. Glacier mass loss is evident over the study period, with a mean of −0.84 m w.e. yr −1 for the period 2003/2004–2007/2008, and also contributes to increase glacier runoff. An El Niño episode in 2002 resulted in high snow accumulation, modifying the hydrological regime and probably reducing the glacier contribution in favor of seasonal snowmelt during the subsequent 2002/2003 hydrological year. At the hourly timescale, summertime glacier contributions are highly variable in space and time, revealing large differences in effective melting rates between glaciers and glacierets (from 1 mm w.e. h −1 to 6 mm w.e. h −1 ).

Description: Coupling of climate models and ice sheet models by surface mass balance gradients: application to the Greenland Ice Sheet The Cryosphere, 6, 255-272, 2012 Author(s): M. M. Helsen, R. S. W. van de Wal, M. R. van den Broeke, W. J. van de Berg, and J. Oerlemans It is notoriously difficult to couple surface mass balance (SMB) results from climate models to the changing geometry of an ice sheet model. This problem is traditionally avoided by using only accumulation from a climate model, and parameterizing the meltwater run-off as a function of temperature, which is often related to surface elevation ( H s ). In this study, we propose a new strategy to calculate SMB, to allow a direct adjustment of SMB to a change in ice sheet topography and/or a change in climate forcing. This method is based on elevational gradients in the SMB field as computed by a regional climate model. Separate linear relations are derived for ablation and accumulation, using pairs of H s and SMB within a minimum search radius. The continuously adjusting SMB forcing is consistent with climate model forcing fields, also for initially non-glaciated areas in the peripheral areas of an ice sheet. When applied to an asynchronous coupled ice sheet – climate model setup, this method circumvents traditional temperature lapse rate assumptions. Here we apply it to the Greenland Ice Sheet (GrIS). Experiments using both steady-state forcing and glacial-interglacial forcing result in realistic ice sheet reconstructions.

Description: Estimating ice phenology on large northern lakes from AMSR-E: algorithm development and application to Great Bear Lake and Great Slave Lake, Canada The Cryosphere, 6, 235-254, 2012 Author(s): K.-K. Kang, C. R. Duguay, and S. E. L. Howell Time series of brightness temperatures ( T B ) from the Advanced Microwave Scanning Radiometer–Earth Observing System (AMSR-E) are examined to determine ice phenology variables on the two largest lakes of northern Canada: Great Bear Lake (GBL) and Great Slave Lake (GSL). T B measurements from the 18.7, 23.8, 36.5, and 89.0 GHz channels (H- and V- polarization) are compared to assess their potential for detecting freeze-onset/melt-onset and ice-on/ice-off dates on both lakes. The 18.7 GHz (H-pol) channel is found to be the most suitable for estimating these ice dates as well as the duration of the ice cover and ice-free seasons. A new algorithm is proposed using this channel and applied to map all ice phenology variables on GBL and GSL over seven ice seasons (2002–2009). Analysis of the spatio-temporal patterns of each variable at the pixel level reveals that: (1) both freeze-onset and ice-on dates occur on average about one week earlier on GBL than on GSL (Day of Year (DY) 318 and 333 for GBL; DY 328 and 343 for GSL); (2) the freeze-up process or freeze duration (freeze-onset to ice-on) takes a slightly longer amount of time on GBL than on GSL (about 1 week on average); (3) melt-onset and ice-off dates occur on average one week and approximately four weeks later, respectively, on GBL (DY 143 and 183 for GBL; DY 135 and 157 for GSL); (4) the break-up process or melt duration (melt-onset to ice-off) lasts on average about three weeks longer on GBL; and (5) ice cover duration estimated from each individual pixel is on average about three weeks longer on GBL compared to its more southern counterpart, GSL. A comparison of dates for several ice phenology variables derived from other satellite remote sensing products (e.g. NOAA Interactive Multisensor Snow and Ice Mapping System (IMS), QuikSCAT, and Canadian Ice Service Database) show that, despite its relatively coarse spatial resolution, AMSR-E 18.7 GHz provides a viable means for monitoring of ice phenology on large northern lakes.

Description: Assessment of a physical-biogeochemical coupled model system for operational service in the Baltic Sea Ocean Science Discussions, 9, 835-876, 2012 Author(s): Z. Wan, J. She, M. Maar, L. Jonasson, and J. Baasch-Larsen Thanks to the abundant observation data, we are able to deploy the traditional point-to-point comparison and statistical measures in combination with a comprehensive model validation scheme to assess the skills of the biogeochemical model ERGOM in providing an operational service for the Baltic Sea. The model assessment concludes that the operational products can resolve the main observed seasonal features for phytoplankton biomass, dissolved inorganic nitrogen, dissolved inorganic phosphorus and dissolved oxygen in euphotic layers, as well as their vertical profiles. This assessment reflects that the model errors of the operational system at the current stage are mainly caused by insufficient light penetration, excessive organic particle export downward, insufficient regional adaptation and some of improper initialization. This study highlights the importance of applying multiple schemes in order to assess model skills rigidly and identify main causes for major model errors.

Description: Numerical tools to estimate the flux of a gas across the air-water interface and assess the heterogeny of its forcing functions Ocean Science Discussions, 9, 909-975, 2012 Author(s): V. M. N. de C. da S. Vieira A numerical tool was developed for the estimation of gas fluxes across the air water interface. The primary objective is to use it to estimate CO 2 fluxes. Nevertheless application to other gases is easily accomplished by changing the values of the parameters related to the physical properties of the gases. A user friendly software was developed allowing to build upon a standard kernel a custom made gas flux model with the preferred parametrizations. These include single or double layer models; several numerical schemes for the effects of wind in the air-side and water-side transfer velocities; the effect of turbulence from current drag with the bottom; and the effects on solubility of water temperature, salinity, air temperature and pressure. It was also developed an analysis which decomposes the difference between the fluxes in a reference situation and in alternative situations into its several forcing functions. This analysis relies on the Taylor expansion of the gas flux model, requiring the numerical estimation of partial derivatives by a multivariate version of the collocation polynomial. Both the flux model and the difference decomposition analysis were tested with data taken from surveys done in the lagoonary system of Ria Formosa, south Portugal, in which the CO 2 fluxes were estimated using the IRGA and floating chamber method whereas the CO 2 concentrations were estimated using the IRGA and degasification chamber. Observations and estimations show a remarkable fit.

Description: Arctic surface temperatures from Metop AVHRR compared to in situ ocean and land data Ocean Science Discussions, 9, 1009-1043, 2012 Author(s): G. Dybkjær, R. Tonboe, and J. Høyer The ice surface temperature (IST) is an important boundary condition for both atmospheric and ocean and sea ice models and for coupled systems. An operational ice surface temperature product using satellite Metop AVHRR infra-red data was developed for MyOcean. The IST can be mapped in clear sky regions using a split window algorithm specially tuned for sea ice. Clear sky conditions are prevailing during spring in the Arctic while persistent cloud cover limits data coverage during summer. The cloud covered regions are detected using the EUMETSAT cloud mask. The Metop IST compares to 2 m temperature at the Greenland ice cap Summit within STD error of 3.14 °C and to Arctic drifting buoy temperature data within STD error of 3.69 °C. A case study reveal that the in situ radiometer data versus satellite IST STD error can be much lower (0.73 °C) and that the different in situ measures complicates the validation. Differences and variability between Metop IST and in situ data are analysed and discussed. An inter-comparison of Metop IST, numerical weather prediction temperatures and in situ observation indicates large biases between the different quantities. Because of the scarcity of conventional surface temperature or surface air temperature data in the Arctic the satellite IST data with its relatively good coverage can potentially add valuable information to model analysis for the Arctic atmosphere.

Description: Towards high resolution mapping of 3-D mesoscale dynamics from observations: preliminary comparison of retrieval techniques and models within MESCLA project Ocean Science Discussions, 9, 1045-1083, 2012 Author(s): B. Buongiorno Nardelli, S. Guinehut, A. Pascual, Y. Drillet, S. Ruiz, and S. Mulet Within the MyOcean R&D project MESCLA (MEsoSCale dynamical Analysis through combined model, satellite and in situ data), different estimates of the vertical velocities derived from observations have been compared. Two main approaches have been considered, one based on the retrieval of 3-D fields from the observations alone and one based on the analyses provided by MyOcean MERCATOR models. The motivation for this double approach is that, while data assimilation in numerical models is crucial to obtain more accurate analyses and forecasts, its results might be significantly influenced by specific model configurations (e.g. forcing, parameterization of smaller scale processes and spatial resolution). On the other hand, the purely observation-based approach is limited by the underlying assumptions of simplified dynamical models and by the relatively low resolution of present products. MESCLA tested innovative methods for the high resolution mapping of 3-D mesoscale dynamics from observations, developing new products that might be used to gradually build the next generations of operational observation-based products.

Description: On the Shelf Resonances of the Gulf of Carpentaria and the Arafura Sea Ocean Science Discussions, 9, 443-497, 2012 Author(s): D. J. Webb A numerical model is used to investigate the resonances of the Gulf of Carpentaria and the Arafura Sea. The model is forced at the shelf edge, first with physically realistic real values of angular velocity. The response functions at points within the region show maxima and other behaviour which imply that resonances are involved but it is difficult to be more specific. The study is then extended to complex angular velocities and the results then show a clear pattern of gravity wave and Rossby wave like resonances. The properties of the resonances are investigated and used to reinterpret the responses at real values of angular velocity. It is found that in some regions the response is dominated by modes trapped between the shelf edge and the coast or between opposing coastlines. In other regions the resonances show cooperative behaviour, possibly indicating the importance of other physical processes.

Description: Large surface meltwater discharge from the Kangerlussuaq sector of the Greenland ice sheet during the record-warm year 2010 explained by detailed energy balance observations The Cryosphere, 6, 199-209, 2012 Author(s): D. van As, A. L. Hubbard, B. Hasholt, A. B. Mikkelsen, M. R. van den Broeke, and R. S. Fausto This study uses data from six on-ice weather stations, calibrated MODIS-derived albedo and proglacial river gauging measurements to drive and validate an energy balance model. We aim to quantify the record-setting positive temperature anomaly in 2010 and its effect on mass balance and runoff from the Kangerlussuaq sector of the Greenland ice sheet. In 2010, the average temperature was 4.9 °C (2.7 standard deviations) above the 1974–2010 average in Kangerlussuaq. High temperatures were also observed over the ice sheet, with the magnitude of the positive anomaly increasing with altitude, particularly in August. Simultaneously, surface albedo was anomalously low in 2010, predominantly in the upper ablation zone. The low albedo was caused by high ablation, which in turn profited from high temperatures and low winter snowfall. Surface energy balance calculations show that the largest melt excess (∼170%) occurred in the upper ablation zone (above 1000 m), where higher temperatures and lower albedo contributed equally to the melt anomaly. At lower elevations the melt excess can be attributed to high atmospheric temperatures alone. In total, we calculate that 6.6 ± 1.0 km 3 of surface meltwater ran off the ice sheet in the Kangerlussuaq catchment in 2010, exceeding the reference year 2009 (based on atmospheric temperature measurements) by ∼150%. During future warm episodes we can expect a melt response of at least the same magnitude, unless a larger wintertime snow accumulation delays and moderates the melt-albedo feedback. Due to the hypsometry of the ice sheet, yielding an increasing surface area with elevation, meltwater runoff will be further amplified by increases in melt forcings such as atmospheric heat.

Description: NEMO on the shelf: assessment of the Iberia-Biscay-Ireland configuration Ocean Science Discussions, 10, 83-151, 2013 Author(s): C. Maraldi, J. Chanut, B. Levier, N. Ayoub, P. De Mey, G. Reffray, F. Lyard, S. Cailleau, M. Drévillon, E. A. Fanjul, M. G. Sotillo, and P. Marsaleix The Iberia-Biscay-Ireland (IBI) system serves one of the 7 MyOcean "Monitoring and Forecasting Centres". A high resolution simulation covering the IBI region is set-up over July 2007–February 2009. The NEMO (Nucleus for European Modelling of the Ocean) model is used with a 1/36° horizontal resolution and 50 z-levels in the vertical. New developments have been incorporated in NEMO to make it suitable to open- as well as coastal-ocean modelling. In this paper, we pursue three main objectives: (1) give an overview of the model configuration used for the simulations; (2) give a broad-brush account of one particular aspect of this work, namely consistency verification; this type of validation is conducted upstream of the implementation of the system before it is used for production and routinely validated; it is meant to guide model development in identifying gross deficiencies in the modelling of several key physical processes; (3) show that such a regional modelling system has potential as a complement to patchy observations (an integrated approach) to give information on non-observed physical quantities and to provide links between observations by identifying broader-scale patterns and processes. We concentrate on the year 2008. We first provide domain-wide consistency verification results in terms of barotropic tides, transports, sea surface temperature and stratification. We then focus on two dynamical sub-regions: the Celtic shelves and the Bay of Biscay slope and deep regions. The model-data consistency is checked for variables and processes such as tidal currents, tidal fronts, internal tides, residual elevation. We also examine the representation in the model of a seasonal pattern of the Bay of Biscay circulation: the warm extension of the Iberian Poleward Current along the northern Spanish coast (Navidad event) in winter 2007–2008.

Description: Climatic drivers of seasonal glacier mass balances: an analysis of 6 decades at Glacier de Sarennes (French Alps) The Cryosphere, 7, 47-66, 2013 Author(s): E. Thibert, N. Eckert, and C. Vincent Refined temporal signals extracted from a winter and summer mass balance series recorded at Glacier de Sarennes (French Alps) using variance decomposition are related to local meteorological data and large-scale North Atlantic Oscillation (NAO) anomalies in terms of interannual variability, trends of the low-frequency signals, and breaks in the time series. The winter balance has increased by +23% since 1976 due to more precipitation in early and late winter. The summer balance has decreased since 1982 due to a 43% increase in snow and ice melt. A 24-day lengthening of the ablation period – mainly due to longer ice ablation – is the main component in the overall increase in ablation. In addition, the last 25 yr have seen increases in ablation rates of 14 and 10% for snow and ice, respectively. A simple degree-day analysis can account for both the snow/ice melt rate rise and the lengthening of the ablation period as a function of higher air temperatures. From the same analysis, the equilibrium-line altitude of this 45° N latitude south-facing glacier has a sensitivity to temperature of +93 m °C −1 around its mean elevation of 3100 m a.s.l. over 6 decades. The sensitivity of summer balance to temperature is −0.62 m w.e. yr −1 °C −1 for a typical 125-day-long ablation season. Finally, the correlation of winter and summer mass balance terms with NAO anomalies is investigated. Singularly, highest values are obtained between winter NAO anomalies and summer balance. Winter NAO anomalies and winter balance and precipitation are almost disconnected. However, these results strongly depend on how the NAO signal is smoothed, so that the link between Sarennes mass balance seasonal terms and NAO signal remains tenuous and hard to interpret.

Description: A synthesis of the Antarctic surface mass balance during the last 800 yr The Cryosphere, 7, 303-319, 2013 Author(s): M. Frezzotti, C. Scarchilli, S. Becagli, M. Proposito, and S. Urbini Global climate models suggest that Antarctic snowfall should increase in a warming climate and mitigate rises in the sea level. Several processes affect surface mass balance (SMB), introducing large uncertainties in past, present and future ice sheet mass balance. To provide an extended perspective on the past SMB of Antarctica, we used 67 firn/ice core records to reconstruct the temporal variability in the SMB over the past 800 yr and, in greater detail, over the last 200 yr. Our SMB reconstructions indicate that the SMB changes over most of Antarctica are statistically negligible and that the current SMB is not exceptionally high compared to the last 800 yr. High-accumulation periods have occurred in the past, specifically during the 1370s and 1610s. However, a clear increase in accumulation of more than 10% has occurred in high SMB coastal regions and over the highest part of the East Antarctic ice divide since the 1960s. To explain the differences in behaviour between the coastal/ice divide sites and the rest of Antarctica, we suggest that a higher frequency of blocking anticyclones increases the precipitation at coastal sites, leading to the advection of moist air in the highest areas, whereas blowing snow and/or erosion have significant negative impacts on the SMB at windy sites. Eight hundred years of stacked records of the SMB mimic the total solar irradiance during the 13th and 18th centuries. The link between those two variables is probably indirect and linked to a teleconnection in atmospheric circulation that forces complex feedback between the tropical Pacific and Antarctica via the generation and propagation of a large-scale atmospheric wave train.

Description: Impact of the Indonesian throughflow on Agulhas leakage Ocean Science Discussions, 10, 353-391, 2013 Author(s): D. Le Bars, H. A. Dijkstra, and W. P. M. De Ruijter Using ocean models of different complexity we show that opening the Indonesian Passage between the Pacific and the Indian Ocean increases the input of Indian Ocean water into the South Atlantic via the Agulhas leakage. In a strongly eddying global ocean model this response results from an increased Agulhas Current transport and a constant proportion of Agulhas retroflection south of Africa. The leakage increases through an increased frequency of ring shedding events. In an idealized two-layer and flat-bottom eddy resolving model, the proportion of the Agulhas Current transport that retroflects is (for a wide range of wind stress forcing) not affected by an opening of the Indonesian Passage. A linear ocean model is not able to explain this behavior which reveals the importance of mixed barotropic/baroclinic instabilities in controlling the Agulhas leakage.

Description: On the shelf resonances of the English Channel and Irish Sea Ocean Science Discussions, 10, 393-433, 2013 Author(s): D. J. Webb The resonances of the English Channel and Irish Sea are investigated using the methods of Webb (2012) together with an Arakawa C-grid model of the region under study. In the semi-diurnal tidal band, the high tides of the Bristol Channel and Gulf of St. Malo are shown to be due to two shelf resonances which strongly couple the two regions. In the diurnal band, the response is complicated by the presence of continental shelf waves.

Description: Future Arctic marine access: analysis and evaluation of observations, models, and projections of sea ice The Cryosphere, 7, 321-332, 2013 Author(s): T. S. Rogers, J. E. Walsh, T. S. Rupp, L. W. Brigham, and M. Sfraga There is an emerging need for regional applications of sea ice projections to provide more accuracy and greater detail to scientists, national, state and local planners, and other stakeholders. The present study offers a prototype for a comprehensive, interdisciplinary study to bridge observational data, climate model simulations, and user needs. The study's first component is an observationally based evaluation of Arctic sea ice trends during 1980–2008, with an emphasis on seasonal and regional differences relative to the overall pan-Arctic trend. Regional sea ice loss has varied, with a significantly larger decline of winter maximum (January–March) extent in the Atlantic region than in other sectors. A lead–lag regression analysis of Atlantic sea ice extent and ocean temperatures indicates that reduced sea ice extent is associated with increased Atlantic Ocean temperatures. Correlations between the two variables are greater when ocean temperatures lag rather than lead sea ice. The performance of 13 global climate models is evaluated using three metrics to compare sea ice simulations with the observed record. We rank models over the pan-Arctic domain and regional quadrants and synthesize model performance across several different studies. The best performing models project reduced ice cover across key access routes in the Arctic through 2100, with a lengthening of seasons for marine operations by 1–3 months. This assessment suggests that the Northwest and Northeast Passages hold potential for enhanced marine access to the Arctic in the future, including shipping and resource development opportunities.

Description: Greenland ice sheet surface mass balance: evaluating simulations and making projections with regional climate models The Cryosphere, 6, 1275-1294, 2012 Author(s): J. G. L. Rae, G. Aðalgeirsdóttir, T. L. Edwards, X. Fettweis, J. M. Gregory, H. T. Hewitt, J. A. Lowe, P. Lucas-Picher, R. H. Mottram, A. J. Payne, J. K. Ridley, S. R. Shannon, W. J. van de Berg, R. S. W. van de Wal, and M. R. van den Broeke Four high-resolution regional climate models (RCMs) have been set up for the area of Greenland, with the aim of providing future projections of Greenland ice sheet surface mass balance (SMB), and its contribution to sea level rise, with greater accuracy than is possible from coarser-resolution general circulation models (GCMs). This is the first time an intercomparison has been carried out of RCM results for Greenland climate and SMB. Output from RCM simulations for the recent past with the four RCMs is evaluated against available observations. The evaluation highlights the importance of using a detailed snow physics scheme, especially regarding the representations of albedo and meltwater refreezing. Simulations with three of the RCMs for the 21st century using SRES scenario A1B from two GCMs produce trends of between −5.5 and −1.1 Gt yr −2 in SMB (equivalent to +0.015 and +0.003 mm sea level equivalent yr −2 ), with trends of smaller magnitude for scenario E1, in which emissions are mitigated. Results from one of the RCMs whose present-day simulation is most realistic indicate that an annual mean near-surface air temperature increase over Greenland of ~ 2°C would be required for the mass loss to increase such that it exceeds accumulation, thereby causing the SMB to become negative, which has been suggested as a threshold beyond which the ice sheet would eventually be eliminated.

Description: Simulating snow maps for Norway: description and statistical evaluation of the seNorge snow model The Cryosphere, 6, 1323-1337, 2012 Author(s): T. M. Saloranta Daily maps of snow conditions have been produced in Norway with the seNorge snow model since 2004. The seNorge snow model operates with 1 × 1 km resolution, uses gridded observations of daily temperature and precipitation as its input forcing, and simulates, among others, snow water equivalent (SWE), snow depth (SD), and the snow bulk density (ρ). In this paper the set of equations contained in the seNorge model code is described and a thorough spatiotemporal statistical evaluation of the model performance from 1957–2011 is made using the two major sets of extensive in situ snow measurements that exist for Norway. The evaluation results show that the seNorge model generally overestimates both SWE and ρ, and that the overestimation of SWE increases with elevation throughout the snow season. However, the R 2 -values for model fit are 0.60 for (log-transformed) SWE and 0.45 for ρ, indicating that after removal of the detected systematic model biases (e.g. by recalibrating the model or expressing snow conditions in relative units) the model performs rather well. The seNorge model provides a relatively simple, not very data-demanding, yet nonetheless process-based method to construct snow maps of high spatiotemporal resolution. It is an especially well suited alternative for operational snow mapping in regions with rugged topography and large spatiotemporal variability in snow conditions, as is the case in the mountainous Norway.

Description: Ground penetrating radar detection of subsnow slush on ice-covered lakes in interior Alaska The Cryosphere, 6, 1435-1443, 2012 Author(s): A. Gusmeroli and G. Grosse Lakes are abundant throughout the pan-Arctic region. For many of these lakes ice cover lasts for up to two thirds of the year. The frozen cover allows human access to these lakes, which are therefore used for many subsistence and recreational activities, including water harvesting, fishing, and skiing. Safe traveling condition onto lakes may be compromised, however, when, after significant snowfall, the weight of the snow acts on the ice and causes liquid water to spill through weak spots and overflow at the snow-ice interface. Since visual detection of subsnow slush is almost impossible our understanding on overflow processes is still very limited and geophysical methods that allow water and slush detection are desirable. In this study we demonstrate that a commercially available, lightweight 1 GHz, ground penetrating radar system can detect and map extent and intensity of overflow. The strength of radar reflections from wet snow-ice interfaces are at least twice as much in strength than returns from dry snow-ice interface. The presence of overflow also affects the quality of radar returns from the base of the lake ice. During dry conditions we were able to profile ice thickness of up to 1 m, conversely, we did not retrieve any ice-water returns in areas affected by overflow.

Description: The footprint of Asian monsoon dynamics in the mass and energy balance of a Tibetan glacier The Cryosphere, 6, 1445-1461, 2012 Author(s): T. Mölg, F. Maussion, W. Yang, and D. Scherer Determinations of glacier-wide mass and energy balance are still scarce for the remote mountains of the Tibetan Plateau, where field measurements are challenging. Here we run and evaluate a physical, distributed mass balance model for Zhadang Glacier (central Tibet, 30° N) based on in-situ measurements over 2009–2011 and an uncertainty estimate by Monte Carlo and ensemble strategies. The model application aims to provide the first quantification of how the Indian Summer Monsoon (ISM) impacts an entire glacier over the various stages of the monsoon's annual cycle. We find a strong and systematic ISM footprint on the interannual scale. Early (late) monsoon onset causes higher (lower) accumulation, and reduces (increases) the available energy for ablation primarily through changes in absorbed shortwave radiation. By contrast, only a weak footprint exists in the ISM cessation phase. Most striking though is the core monsoon season: local mass and energy balance variability is fully decoupled from the active/break cycle that defines large-scale atmospheric variability during the ISM. Our results demonstrate quantitatively that monsoon onset strongly affects the ablation season of glaciers in Tibet. However, we find no direct ISM impact on the glacier in the main monsoon season, which has not been acknowledged so far. This result also adds cryospheric evidence that, once the monsoon is in full swing, regional atmospheric variability prevails on the Tibetan Plateau in summer.

Description: On the use of the Strouhal/Stokes number to explain the dynamics and water column structure on shelf seas Ocean Science Discussions, 9, 3723-3738, 2012 Author(s): A. J. Souza In recent years coastal oceanographers have suggested the use of the "Strouhal" number or it's inverse the "Stokes" number, which have been defined as the ratios of the frictional depth (δ) to the water column depth ( h ) or vice versa, to describe the effect of bottom boundary layer turbulence on the vertical structure of both density and currents. Although they have mention that the effects of rotation should be important, they have tended to omit it. This omission may be important when talking about tidal currents as the frictional depth from a fully cyclonic to a fully anticyclonic tidal ellipse can vary up to an order of magnitude in the mid latitudes; so that the stokes number might appear smaller (larger) than it is resulting in frictional effects being underestimated (overestimated). Here a way to calculate a Stokes number, in which the effect of the Earth's rotation is taken into account, is suggested. Then the standard Stokes and the rotational Stokes numbers are used as predictors for the position of the tidal mixing fronts in the Irish Sea. Results show that the rotational number improves prediction of the front in shallow cyclonic areas of the eastern Irish Sea. This suggest that the effect of rotation on the water column structure will be more important in shallow shelf seas and estuaries with strong rotational currents.

Description: The first complete inventory of the local \newline glaciers and ice caps on Greenland The Cryosphere, 6, 1483-1495, 2012 Author(s): P. Rastner, T. Bolch, N. Mölg, H. Machguth, R. Le Bris, and F. Paul Glacier inventories provide essential baseline information for the determination of water resources, glacier-specific changes in area and volume, climate change impacts as well as past, potential and future contribution of glaciers to sea-level rise. Although Greenland is heavily glacierised and thus highly relevant for all of the above points, a complete inventory of its glaciers was not available so far. Here we present the results and details of a new and complete inventory that has been compiled from more than 70 Landsat scenes (mostly acquired between 1999 and 2002) using semi-automated glacier mapping techniques. A digital elevation model (DEM) was used to derive drainage divides from watershed analysis and topographic attributes for each glacier entity. To serve the needs of different user communities, we assigned to each glacier one of three connectivity levels with the ice sheet (CL0, CL1, CL2; i.e. no, weak, and strong connection) to clearly, but still flexibly, distinguish the local glaciers and ice caps (GIC) from the ice sheet and its outlet glaciers. In total, we mapped ~ 20 300 glaciers larger than 0.05 km 2 (of which ~ 900 are marine terminating), covering an area of 130 076 ± 4032 km 2 , or 89 720 ± 2781 km 2 without the CL2 GIC. The latter value is about 50% higher than the mean value of more recent previous estimates. Glaciers smaller than 0.5 km 2 contribute only 1.5% to the total area but more than 50% (11 000) to the total number. In contrast, the 25 largest GIC (〉 500 km 2 ) contribute 28% to the total area, but only 0.1% to the total number. The mean elevation of the GIC is 1700 m in the eastern sector and around 1000 m otherwise. The median elevation increases with distance from the coast, but has only a weak dependence on mean glacier aspect.

Description: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica The Cryosphere, 7, 375-393, 2013 Author(s): P. Fretwell, H. D. Pritchard, D. G. Vaughan, J. L. Bamber, N. E. Barrand, R. Bell, C. Bianchi, R. G. Bingham, D. D. Blankenship, G. Casassa, G. Catania, D. Callens, H. Conway, A. J. Cook, H. F. J. Corr, D. Damaske, V. Damm, F. Ferraccioli, R. Forsberg, S. Fujita, Y. Gim, P. Gogineni, J. A. Griggs, R. C. A. Hindmarsh, P. Holmlund, J. W. Holt, R. W. Jacobel, A. Jenkins, W. Jokat, T. Jordan, E. C. King, J. Kohler, W. Krabill, M. Riger-Kusk, K. A. Langley, G. Leitchenkov, C. Leuschen, B. P. Luyendyk, K. Matsuoka, J. Mouginot, F. O. Nitsche, Y. Nogi, O. A. Nost, S. V. Popov, E. Rignot, D. M. Rippin, A. Rivera, J. Roberts, N. Ross, M. J. Siegert, A. M. Smith, D. Steinhage, M. Studinger, B. Sun, B. K. Tinto, B. C. Welch, D. Wilson, D. A. Young, C. Xiangbin, and A. Zirizzotti We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km 3 ) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.

Description: Variability and trends in Laptev Sea ice outflow between 1992–2011 The Cryosphere, 7, 349-363, 2013 Author(s): T. Krumpen, M. Janout, K. I. Hodges, R. Gerdes, F. Girard-Ardhuin, J. A. Hölemann, and S. Willmes Variability and trends in seasonal and interannual ice area export out of the Laptev Sea between 1992 and 2011 are investigated using satellite-based sea ice drift and concentration data. We found an average total winter (October to May) ice area transport across the northern and eastern Laptev Sea boundaries (NB and EB) of 3.48 × 10 5 km 2 . The average transport across the NB (2.87 × 10 5 km 2 ) is thereby higher than across the EB (0.61 × 10 5 km 2 ), with a less pronounced seasonal cycle. The total Laptev Sea ice area flux significantly increased over the last decades (0.85 × 10 5 km 2 decade −1 , p 〉 0.95), dominated by increasing export through the EB (0.55 × 10 5 km 2 decade −1 , p 〉 0.90), while the increase in export across the NB is smaller (0.3 × 10 5 km 2 decade −1 ) and statistically not significant. The strong coupling between across-boundary SLP gradient and ice drift velocity indicates that monthly variations in ice area flux are primarily controlled by changes in geostrophic wind velocities, although the Laptev Sea ice circulation shows no clear relationship with large-scale atmospheric indices. Also there is no evidence of increasing wind velocities that could explain the overall positive trends in ice export. The increased transport rates are rather the consequence of a changing ice cover such as thinning and/or a decrease in concentration. The use of a back-propagation method revealed that most of the ice that is incorporated into the Transpolar Drift is formed during freeze-up and originates from the central and western part of the Laptev Sea, while the exchange with the East Siberian Sea is dominated by ice coming from the central and southeastern Laptev Sea. Furthermore, our results imply that years of high ice export in late winter (February to May) have a thinning effect on the ice cover, which in turn preconditions the occurence of negative sea ice extent anomalies in summer.

Description: Retention and radiative forcing of black carbon in eastern Sierra Nevada snow The Cryosphere, 7, 365-374, 2013 Author(s): K. M. Sterle, J. R. McConnell, J. Dozier, R. Edwards, and M. G. Flanner When contaminated by absorbing particles, such as refractory black carbon (rBC) and continental dust, snow's albedo decreases and thus its absorption of solar radiation increases, thereby hastening snowmelt. For this reason, an understanding of rBC's affect on snow albedo, melt processes, and radiation balance is critical for water management, especially in a changing climate. Measurements of rBC in a sequence of snow pits and surface snow samples in the eastern Sierra Nevada of California during the snow accumulation and ablation seasons of 2009 show that concentrations of rBC were enhanced sevenfold in surface snow (~25 ng g –1 ) compared to bulk values in the snowpack (~3 ng g –1 ). Unlike major ions, which were preferentially released during the initial melt, rBC and continental dust were retained in the snow, enhancing concentrations well into late spring, until a final flush occurred during the ablation period. We estimate a combined rBC and continental dust surface radiative forcing of 20 to 40 W m −2 during April and May, with dust likely contributing a greater share of the forcing.

Description: Surface undulations of Antarctic ice streams tightly controlled by bedrock topography The Cryosphere, 7, 407-417, 2013 Author(s): J. De Rydt, G. H. Gudmundsson, H. F. J. Corr, and P. Christoffersen Full Stokes flow-line models predict that fast-flowing ice streams transmit information about their bedrock topography most efficiently to the surface for basal undulations with length scales between 1 and 20 times the mean ice thickness. This typical behaviour is independent of the precise values of the flow law and sliding law exponents, and should be universally observable. However, no experimental evidence for this important theoretical prediction has been obtained so far, hence ignoring an important test for the physical validity of current-day ice flow models. In our work we use recently acquired airborne radar data for the Rutford Ice Stream and Evans Ice Stream, and we show that the surface response of fast-flowing ice is highly sensitive to bedrock irregularities with wavelengths of several ice thicknesses. The sensitivity depends on the slip ratio, i.e. the ratio between mean basal sliding velocity and mean deformational velocity. We find that higher values of the slip ratio generally lead to a more efficient transfer, whereas the transfer is significantly dampened for ice that attains most of its surface velocity by creep. Our findings underline the importance of bedrock topography for ice stream dynamics on spatial scales up to 20 times the mean ice thickness. Our results also suggest that local variations in the flow regime and surface topography at this spatial scale cannot be explained by variations in basal slipperiness.

Description: Grounding line transient response in marine ice sheet models The Cryosphere, 7, 395-406, 2013 Author(s): A. S. Drouet, D. Docquier, G. Durand, R. Hindmarsh, F. Pattyn, O. Gagliardini, and T. Zwinger Marine ice-sheet stability is mostly controlled by the dynamics of the grounding line, i.e. the junction between the grounded ice sheet and the floating ice shelf. Grounding line migration has been investigated within the framework of MISMIP (Marine Ice Sheet Model Intercomparison Project), which mainly aimed at investigating steady state solutions. Here we focus on transient behaviour, executing short-term simulations (200 yr) of a steady ice sheet perturbed by the release of the buttressing restraint exerted by the ice shelf on the grounded ice upstream. The transient grounding line behaviour of four different flowline ice-sheet models has been compared. The models differ in the physics implemented (full Stokes and shallow shelf approximation), the numerical approach, as well as the grounding line treatment. Their overall response to the loss of buttressing is found to be broadly consistent in terms of grounding line position, rate of surface elevation change and surface velocity. However, still small differences appear for these latter variables, and they can lead to large discrepancies (〉 100%) observed in terms of ice sheet contribution to sea level when cumulated over time. Despite the recent important improvements of marine ice-sheet models in their ability to compute steady state configurations, our results question the capacity of these models to compute short-term reliable sea-level rise projections.

Description: Glacier changes from 1966–2009 in the Gongga Mountains, on the south-eastern margin of the Qinghai-Tibetan Plateau and their climatic forcing The Cryosphere, 6, 1087-1101, 2012 Author(s): B. T. Pan, G. L. Zhang, J. Wang, B. Cao, H. P. Geng, J. Wang, C. Zhang, and Y. P. Ji In order to monitor the changes of the glaciers in the Gongga Mountain region on the south-eastern margin of the Qinghai-Tibetan Plateau, 74 monsoonal temperate glaciers were investigated by comparing the Chinese Glacier Inventory (CGI), recorded in the 1960s, with Landsat MSS in 1974, Landsat TM in 1989, 1994, 2005, and ASTER data in 2009. The remote sensing data have been applied to map the glacier outline by threshold ratio images (TM4/TM5). Moreover, the glacier outlines were verified by GPS survey on four large glaciers (Hailuogou (HLG), Mozigou (MZG), Yanzigou (YZG), and Dagongba (DGB)) in 2009. The results show that the area dominated by the 74 glaciers has shrunk by 11.3% (29.2 km 2 ) from 1966 to 2009. Glacier area on the eastern and western slopes of the Gongga Mountains decreased by 9.8% and 14.6% since 1966, respectively. The loss in glacier area and length is, respectively, 0.8 km 2 and 1146.4 m for the HLG Glacier, 2.1 km 2 and 501.8 m for the MZG Glacier, 0.8 km 2 and 724.8 m for the YZG Glacier, and 2.4 km 2 and 1002.3 m for the DGB Glacier. Decades of climate records obtained from three meteorological stations in the Gongga Mountains were analyzed to evaluate the impact of the temperature and precipitation on glacier retreat. The mean annual temperatures over the eastern and western slopes of the Gongga Mountains have been increasing by 0.34 K decade −1 and 0.24 K decade −1 (1988–2009), respectively. Moreover, mean annual precipitation has only increased by 1% in the past 50 yr. The increasing amount of precipitation could not compensate for the glacier mass loss due to the temperature increase in the Gongga Mountains. This suggests that the warming of the climate is probably also responsible for the glacier retreat in the study region. At the region scale, glacier changes were also controlled by local topographical factors.

Description: Chaotic variability of the meridional overturning circulation on subannual to interannual timescales Ocean Science Discussions, 9, 3191-3238, 2012 Author(s): J. J.-M. Hirschi, A. T. Blaker, B. Sinha, A. Coward, B. de Cuevas, S. Alderson, and G. Madec Observations and numerical simulations have shown that the meridional overturning circulation (MOC) exhibits substantial variability on sub- to interannual timescales. This variability is not fully understood. In particular it is not known what fraction of the MOC variability is caused by processes such as mesoscale ocean eddies and internal waves which are ubiquitous in the ocean. Here we analyse twin experiments performed with a global ocean model at eddying (1/4°) and non-eddying (1°) resolutions. The twin experiments are forced with the same surface fluxes for the 1958 to 2001 period but start from different initial conditions. Our results show that on subannual to interannual timescales a large fraction of MOC variability directly reflects variability in the surface forcing. Nevertheless, in the eddy-permitting case there is an initial condition dependent MOC variability (hereinafter referred to as "chaotic" variability) of several Sv (1 Sv = 10 6 m 3 s −1 ) in the Atlantic and the Indo-Pacific. In the Atlantic the chaotic MOC variability represents up to 30% of the total variability at the depths where the maximum MOC occurs. In comparison the chaotic MOC variability is only 5–10% in the non-eddying case. The surface forcing being identical in the twin experiments suggests that mesoscale ocean eddies are the most likely cause for the increased chaotic MOC variability in the eddying case. The exact formation time of eddies is determined by the initial conditions which are different in the two accordance with and as a consequence the mesoscale eddy field is decorrelated in the twin experiments. In regions where eddy activity is high in the eddy-permitting model, the correlation of sea surface height variability in the twin runs is close to zero. In the non-eddying case in contrast, we find high correlations (0.9 or higher) over most regions. Looking at the sub- and interannual MOC components separately reveals that despite the amplitude of the chaotic variability being larger on subannual than on interannual timescales, the ratio of chaotic to total MOC variability is larger on interannual than on subannual timescales.

Description: Accelerated contributions of Canada's Baffin and Bylot Island glaciers to sea level rise over the past half century The Cryosphere, 6, 1103-1125, 2012 Author(s): A. Gardner, G. Moholdt, A. Arendt, and B. Wouters Canadian Arctic glaciers have recently contributed large volumes of meltwater to the world's oceans. To place recently observed glacier wastage into a historical perspective and to determine the region's longer-term (~50 years) contribution to sea level, we estimate mass and volume changes for the glaciers of Baffin and Bylot Islands using digital elevation models generated from airborne and satellite stereoscopic imagery and elevation postings from repeat airborne and satellite laser altimetry. In addition, we update existing glacier mass change records from GRACE satellite gravimetry to cover the period from 2003 to 2011. Using this integrated approach, we find that the rate of mass loss from the region's glaciers increased from 11.1 ± 3.4 Gt a −1 (271 ± 84 kg m −2 a −1 ) for the period 1963–2006 to 23.8 ± 6.1 Gt a −1 (581 ± 149 kg m −2 a −1 ) for the period 2003–2011. The doubling of the rate of mass loss is attributed to higher temperatures in summer with little change in annual precipitation. Through both direct and indirect effects, changes in summer temperatures accounted for 70–98% of the variance in the rate of mass loss, to which the Barnes Ice Cap was found to be 1.7 times more sensitive than either the Penny Ice Cap or the region's glaciers as a whole. This heightened sensitivity is the result of a glacier hypsometry that is skewed to lower elevations, which are shown to have a higher mass change sensitive to temperature compared to glacier surfaces at higher elevations. Between 2003 and 2011 the glaciers of Baffin and Bylot Islands contributed 0.07 ± 0.02 mm a −1 to sea level rise accounting for 16% of the total contribution from glaciers outside of Greenland and Antarctica, a rate much higher than the longer-term average of 0.03 ± 0.01 mm a −1 (1963 to 2006).

Description: Ice velocity changes in the Ross and Ronne sectors observed using satellite radar data from 1997 and 2009 The Cryosphere, 6, 1019-1030, 2012 Author(s): B. Scheuchl, J. Mouginot, and E. Rignot We report changes in ice velocity of a 6.5 million km 2 region around South Pole encompassing the Filchner-Ronne and Ross Ice Shelves and a significant portion of the ice streams and glaciers that constitute their catchment areas. Using the first full interferometric synthetic aperture radar (InSAR) coverage of the region completed in 2009 and partial coverage acquired in 1997, we processed the data to assemble a comprehensive map of ice speed changes between those two years. On the Ross Ice Shelf, our results confirm a continued deceleration of Mercer and Whillans Ice Streams with a 12-yr velocity difference of −50 m yr −1 (−16.7%) and −100 m yr −1 (−25.3%) at their grounding lines. The deceleration spreads 450 km upstream of the grounding line and more than 500 km onto the shelf, beyond what was previously known. Ross and Filchner Ice Shelves exhibit signs of pre-calving events, representing the largest observed changes, with an increase in speed in excess of +100 m yr −1 in 12 yr. Other changes in the Ross Ice Shelf region are less significant. The observed changes in glacier speed extend on the Ross Ice Shelf along the ice streams' flow lines. Most tributaries of the Filchner-Ronne Ice Shelf show a modest deceleration or no change between 1997 and 2009. Slessor Glacier shows a small deceleration over a large sector. No change is detected on the Bailey, Rutford, and Institute Ice Streams. On the Filchner Ice Shelf itself, ice decelerated rather uniformly with a 12-yr difference in speed of −50 m yr −1 , or −5% of its ice front speed, which we attribute to a 12 km advance in its ice front position. Our results show that dynamic changes are present in the region. They highlight the need for continued observation of the area with a primary focus on the Siple Coast. The dynamic changes in Central Antarctica between 1997 and 2009 are generally second-order effects in comparison to losses on glaciers in the Bellingshausen and Amundsen Seas region and on the Antarctic Peninsula. We therefore conclude that the dynamic changes shown here do not have a strong impact on the mass budget of the Antarctic continent.

Description: Uncertainties in the global temperature change caused by carbon release from permafrost thawing The Cryosphere, 6, 1063-1076, 2012 Author(s): E. J. Burke, I. P. Hartley, and C. D. Jones Under climate change thawing permafrost will cause old carbon which is currently frozen and inert to become vulnerable to decomposition and release into the climate system. This paper develops a simple framework for estimating the impact of this permafrost carbon release on the global mean temperature (P-GMT). The analysis is based on simulations made with the Hadley Centre climate model (HadGEM2-ES) for a range of representative CO 2 concentration pathways. Results using the high concentration pathway (RCP 8.5) suggest that by 2100 the annual methane (CH 4 ) emission rate is 2–59 Tg CH 4 yr −1 and 50–270 Pg C has been released as CO 2 with an associated P-GMT of 0.08–0.36 °C (all 5th–95th percentile ranges). P-GMT is considerably lower – between 0.02 and 0.11 °C – for the low concentration pathway (RCP2.6). The uncertainty in climate model scenario causes about 50% of the spread in P-GMT by the end of the 21st century. The distribution of soil carbon, in particular how it varies with depth, contributes to about half of the remaining spread, with quality of soil carbon and decomposition processes contributing a further quarter each. These latter uncertainties could be reduced through additional observations. Over the next 20–30 yr, whilst scenario uncertainty is small, improving our knowledge of the quality of soil carbon will contribute significantly to reducing the spread in the, albeit relatively small, P-GMT.

Description: Co-existence of wind seas and swells along the west coast of India during non-monsoon season Ocean Science Discussions, 9, 3097-3125, 2012 Author(s): R. Rashmi, V. M. Aboobacker, P. Vethamony, and M. P. John Wave data collected along the west coast of India (off Goa, Ratnagiri and Dwarka) during non-monsoon season have been analysed to study the co-existence of wind seas and swells. Diurnal variation in wind and wave parameters is noticeable along the central west coast of India (off Goa and Ratnagiri), and this is not present along the northwest coast of India (off Dwarka). Swells are predominantly mature (91%) and old (88%) during late pre-monsoon and post-monsoon seasons, respectively. Sea Swell Energy Ratio quantifies wind sea, swell and mixed seas prevailing in the regions during non-monsoon season. Intermodal Distance (ID) between the energy peaks is moderately separated during non-monsoon season, whereas, during the shamal events, energy peaks are very close to each other (ID ~ 0). However, pure wind seas (ID ~ 1) are found to co-exist with the swells during non-monsoon season. Wind seas are growing, when wind and wind seas are opposite to swell direction. Wind seas have minimum angular spreads in multimodal state. Under low winds, the interaction between wind sea and swell dominates and thereby the multimodal state reduces to unimodal state. The fetch available for the evolution of the wind sea spectrum has been estimated, and it is found to be less than 150 km. For the fetch limited condition, a non-dimensional empirical relation has been derived relating the significant wind sea height in terms of wind speed and peak wind sea period, and this relation fits for the west coast of India.

Description: Refreezing on the Greenland ice sheet: a comparison of parameterizations The Cryosphere, 6, 743-762, 2012 Author(s): C. H. Reijmer, M. R. van den Broeke, X. Fettweis, J. Ettema, and L. B. Stap Retention and refreezing of meltwater are acknowledged to be important processes for the mass budget of polar glaciers and ice sheets. Several parameterizations of these processes exist for use in energy and mass balance models. Due to a lack of direct observations, validation of these parameterizations is difficult. In this study we compare a set of 6 refreezing parameterizations against output of two Regional Climate Models (RCMs) coupled to an energy balance snow model, the Regional Atmospheric Climate Model (RACMO2) and the Modèle Atmosphérique Régional (MAR), applied to the Greenland ice sheet. In both RCMs, refreezing is explicitly calculated in a snow model that calculates vertical profiles of temperature, density and liquid water content. Between RACMO2 and MAR, the ice sheet-integrated amount of refreezing differs by only 4.9 mm w.e yr −1 (4.5 %), and the temporal and spatial variability are very similar. For consistency, the parameterizations are forced with output (surface temperature, precipitation and melt) of the RCMs. For the ice sheet-integrated amount of refreezing and its inter-annual variations, all parameterizations give similar results, especially after some tuning. However, the spatial distributions differ significantly and the spatial correspondence between the RCMs is better than with any of the parameterizations. Results are especially sensitive to the choice of the depth of the thermally active layer, which determines the cold content of the snow in most parameterizations. These results are independent of which RCM is used to force the parameterizations.

Description: Vapor flux and recrystallization during dry snow metamorphism under a steady temperature gradient as observed by time-lapse micro-tomography The Cryosphere, 6, 1141-1155, 2012 Author(s): B. R. Pinzer, M. Schneebeli, and T. U. Kaempfer Dry snow metamorphism under an external temperature gradient is the most common type of recrystallization of snow on the ground. The changes in snow microstructure modify the physical properties of snow, and therefore an understanding of this process is essential for many disciplines, from modeling the effects of snow on climate to assessing avalanche risk. We directly imaged the microstructural changes in snow during temperature gradient metamorphism (TGM) under a constant gradient of 50 K m −1 , using in situ time-lapse X-ray micro-tomography. This novel and non-destructive technique directly reveals the amount of ice that sublimates and is deposited during metamorphism, in addition to the exact locations of these phase changes. We calculated the average time that an ice volume stayed in place before it sublimated and found a characteristic residence time of 2–3 days. This means that most of the ice changes its phase from solid to vapor and back many times in a seasonal snowpack where similar temperature conditions can be found. Consistent with such a short timescale, we observed a mass turnover of up to 60% of the total ice mass per day. The concept of hand-to-hand transport for the water vapor flux describes the observed changes very well. However, we did not find evidence for a macroscopic vapor diffusion enhancement. The picture of {temperature gradient metamorphism} produced by directly observing the changing microstructure sheds light on the micro-physical processes and could help to improve models that predict the physical properties of snow.