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: Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions Annales Geophysicae, 31, 1559-1567, 2013 Author(s): E. K. J. Kilpua, H. Hietala, H. E. J. Koskinen, D. Fontaine, and L. Turc Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure ( P dyn ) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P dyn ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma.

Description: Cluster observation of few-hour-scale evolution of structured plasma in the inner magnetosphere Annales Geophysicae, 31, 1569-1578, 2013 Author(s): M. Yamauchi, I. Dandouras, H. Rème, R. Lundin, and L. M. Kistler Using Cluster Ion Spectrometry (CIS) data from the spacecraft-4 perigee traversals during the 2001–2006 period (nearly 500 traversals after removing those that are highly contaminated by radiation belt particles), we statistically examined the local time distribution of structured trapped ions at sub- to few-keV range as well as inbound–outbound differences of these ion signatures in intensities and energy–latitude dispersion directions. Since the Cluster orbit during this period was almost constant and approximately north–south symmetric at nearly constant local time near the perigee, inbound–outbound differences are attributed to temporal developments in a 1–2 h timescale. Three types of structured ions at sub- to few keV range that are commonly found in the inner magnetosphere are examined: – Energy–latitude dispersed structured ions at less than a few keV, – Short-lived dispersionless ion stripes at wide energy range extending 0.1–10 keV, – Short-lived low-energy ion bursts at less than a few hundred eV. The statistics revealed that the wedge-like dispersed ions are most often observed in the dawn sector (60% of traversals), and a large portion of them show significant enhancement during the traversals at all local times. The short-lived ion stripes are predominantly found near midnight, where most stripes are significantly enhanced during the traversals and are associated with substorm activities with geomagnetic AL 〈 −300 nT. The low-energy bursts are observed at all local times and under all geomagnetic conditions, with moderate peak of the occurrence rate in the afternoon sector. A large portion of them again show significant enhancement or decay during the traversals.

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: Influence of the terrestrial magnetic field geometry on the cutoff rigidity of cosmic ray particles Annales Geophysicae, 31, 1637-1643, 2013 Author(s): K. Herbst, A. Kopp, and B. Heber Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Størmer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity R C ; that is, the minimum momentum per charge a particle must have in order to reach a certain geographical location. Employing Monte Carlo simulations with the PLANETOCOSMICS code we investigate the correlation between the geomagnetic field structure and the cutoff rigidity. We show that the geometry of the magnetic field has a considerable influence on the resulting cutoff rigidity distribution. Furthermore, we will present a simple geometry-based parameter, δ B , which is able to reflect the location-dependent cutoff rigidity. We show that this correlation is also visible in the temporal evolution of the Earth's magnetic field, at least over the last 100 yr. Using latitude scans with neutron monitors, changes of the relative counting rates at different positions are calculated, showing small variations for, e.g., Kiel and Moscow, while large ones occur at Mexico City as well as on the British Virgin Islands.

Description: Electron acceleration at Jupiter: input from cyclotron-resonant interaction with whistler-mode chorus waves Annales Geophysicae, 31, 1619-1630, 2013 Author(s): E. E. Woodfield, R. B. Horne, S. A. Glauert, J. D. Menietti, and Y. Y. Shprits Jupiter has the most intense radiation belts of all the outer planets. It is not yet known how electrons can be accelerated to energies of 10 MeV or more. It has been suggested that cyclotron-resonant wave-particle interactions by chorus waves could accelerate electrons to a few MeV near the orbit of Io. Here we use the chorus wave intensities observed by the Galileo spacecraft to calculate the changes in electron flux as a result of pitch angle and energy diffusion. We show that, when the bandwidth of the waves and its variation with L are taken into account, pitch angle and energy diffusion due to chorus waves is a factor of 8 larger at L -shells greater than 10 than previously shown. We have used the latitudinal wave intensity profile from Galileo data to model the time evolution of the electron flux using the British Antarctic Survey Radiation Belt (BAS) model. This profile confines intense chorus waves near the magnetic equator with a peak intensity at ∼5° latitude. Electron fluxes in the BAS model increase by an order of magnitude for energies around 3 MeV. Extending our results to L = 14 shows that cyclotron-resonant interactions with chorus waves are equally important for electron acceleration beyond L = 10. These results suggest that there is significant electron acceleration by cyclotron-resonant interactions at Jupiter contributing to the creation of Jupiter's radiation belts and also increasing the range of L -shells over which this mechanism should be considered.

Description: The coupling between the solar wind and proton fluxes at GEO Annales Geophysicae, 31, 1631-1636, 2013 Author(s): R. J. Boynton, S. A. Billings, O. A. Amariutei, and I. Moiseenko The relationship between the solar wind and the proton flux at geosynchronous Earth orbit (GEO) is investigated using the error reduction ratio (ERR) analysis. The ERR analysis is able to search for the most appropriate inputs that control the evolution of the system. This approach is a black box method and is able to derive a mathematical model of a system from input-output data. This method is used to analyse eight energy ranges of the proton flux at GEO from 80 keV to 14.5 MeV. The inputs to the algorithm were solar wind velocity, density and pressure; the Dst index; the solar energetic proton (SEP) flux; and a function of the interplanetary magnetic field (IMF) tangential magnitude and clock angle. The results show that for lowest five energy channels (80 to 800 keV) the GEO proton fluxes are controlled by the solar wind velocity with a lag of two to three days. However, above 350 keV, the SEP fluxes, accounts for a significant portion of the GEO proton flux variance. For the highest three energy channels (0.74 to 14.5 MeV), the SEPs account for the majority of the ERR. The results also show an anisotropy of protons with gyrocenters inside GEO and outside GEO, where the protons inside GEO are controlled partly by the Dst index and also an IMF-clock angle function.

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: A statistical study on O + flux in the dayside magnetosheath Annales Geophysicae, 31, 1005-1010, 2013 Author(s): R. Slapak, H. Nilsson, and L. G. Westerberg Studies on terrestrial oxygen ion (O + ) escape into the interplanetary space have considered a number of different escape paths. Recent observations however suggest a yet insufficiently investigated additional escape route for hot O + : along open magnetic field lines in the high altitude cusp and mantle. Here we present a statistical study on O + flux in the high-latitude dayside magnetosheath. The O + is generally seen relatively close to the magnetopause, consistent with observations of O + flowing primarily tangentially to the magnetopause. We estimate the total escape flux in this region to be ~ 7 × 10 24 s −1 , implying this escape route to significantly contribute to the overall total O + loss into interplanetary space.

Description: Simulation of the influence of historical land cover changes on the global climate Annales Geophysicae, 31, 995-1004, 2013 Author(s): Y. Wang, X. Yan, and Z. Wang In order to estimate biogeophysical effects of historical land cover change on climate during last three centuries, a set of experiments with a climate system model of intermediate complexity (MPM-2) is performed. In response to historical deforestation, the model simulates a decrease in annual mean global temperature in the range of 0.07–0.14 °C based on different grassland albedos. The effect of land cover changes is most pronounced in the middle northern latitudes with maximum cooling reaching approximately 0.6 °C during northern summer. The cooling reaches 0.57 °C during northern spring owing to the large effects of land surface albedo. These results suggest that land cover forcing is important for study on historical climate change and that more research is necessary in the assessment of land management options for climate change mitigation.

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: IMF effect on the polar cap contraction and expansion during a period of substorms Annales Geophysicae, 31, 1021-1034, 2013 Author(s): A. T. Aikio, T. Pitkänen, I. Honkonen, M. Palmroth, and O. Amm The polar cap boundary (PCB) location and motion in the nightside ionosphere has been studied by using measurements from the EISCAT radars and the MIRACLE magnetometers during a period of four substorms on 18 February 2004. The OMNI database has been used for observations of the solar wind and the Geotail satellite for magnetospheric measurements. In addition, the event was modelled by the GUMICS-4 MHD simulation. The simulation of the PCB location was in a rather good agreement with the experimental estimates at the EISCAT longitude. During the first three substorm expansion phases, neither the local observations nor the global simulation showed any poleward motions of the PCB, even though the electrojets intensified. Rapid poleward motions of the PCB took place only in the early recovery phases of the substorms. Hence, in these cases the nightside reconnection rate was locally higher in the recovery phase than in the expansion phase. In addition, we suggest that the IMF B z component correlated with the nightside tail inclination angle and the PCB location with about a 17-min delay from the bow shock. By taking the delay into account, the IMF northward turnings were associated with dipolarizations of the magnetotail and poleward motions of the PCB in the recovery phase. The mechanism behind this effect should be studied further.

Description: A comparison of bow shock models with Cluster observations during low Alfvén Mach number magnetic clouds Annales Geophysicae, 31, 1011-1019, 2013 Author(s): L. Turc, D. Fontaine, P. Savoini, H. Hietala, and E. K. J. Kilpua Magnetic clouds (MCs) are very geoeffective solar wind structures. Their properties in the interplanetary medium have been extensively studied, yet little is known about their characteristics in the Earth's magnetosheath. The Cluster spacecraft offer the opportunity to observe MCs in the magnetosheath, but before MCs reach the magnetosphere, their structure is altered when they interact with the terrestrial bow shock (BS). The physics taking place at the BS strongly depends on Θ Bn , the angle between the shock normal and the interplanetary magnetic field. However, in situ observations of the BS during an MC's crossing are seldom available. In order to relate magnetosheath observations to solar wind conditions, we need to rely on a model to determine the shock's position and normal direction. Yet during MCs, the models tend to be less accurate, because the Alfvén Mach number ( M A ) is often significantly lower than in regular solar wind. On the contrary, the models are generally optimised for high M A conditions. In this study, we compare the predictions of four widely used models available in the literature (Wu et al., 2000; Chapman and Cairns, 2003; Jeřáb et al., 2005; Měrka et al., 2005b) to Cluster's dayside BS crossings observed during five MC events. Our analysis shows that the Θ Bn angle is well predicted by all four models. On the other hand, the Jeřáb et al. (2005) model yields the best estimates of the BS position during low M A MCs. The other models locate the BS either too far from or too close to Earth. The results of this paper can be directly used to estimate the BS parameters in all studies of MC interaction with Earth's magnetosphere.

Description: Investigation of radiative effects of the optically thick dust layer over the Indian tropical region Annales Geophysicae, 31, 647-663, 2013 Author(s): S. K. Das, J.-P. Chen, M. Venkat Ratnam, and A. Jayaraman Optical and physical properties of aerosols derived from multi-satellite observations (MODIS-Aqua, OMI-Aura, MISR-Terra, CALIOP-CALIPSO) have been used to estimate radiative effects of the dust layer over southern India. The vertical distribution of aerosol radiative forcing and heating rates are calculated with 100 m resolution in the lower atmosphere, using temperature and relative humidity data from balloon-borne radiosonde observations. The present study investigates the optically thick dust layer of optical thickness 0.18 ± 0.06 at an altitude of 2.5 ± 0.7 km over Gadanki, transported from the Thar Desert, producing radiative forcing and heating rate of 11.5 ± 3.3 W m −2 and 0.6 ± 0.26 K day −1 , respectively, with a forcing efficiency of 43 W m −2 and an effective heating rate of 4 K day −1 per unit dust optical depth. Presence of the dust layer increases radiative forcing by 60% and heating rate by 60 times at that altitude compared to non-dusty cloud-free days. Calculation shows that the radiative effects of the dust layer strongly depend on the boundary layer aerosol type and mass loading. An increase of 25% of heating by the dust layer is found over relatively cleaner regions than urban regions in southern India and further 15% of heating increases over the marine region. Such heating differences in free troposphere may have significant consequences in the atmospheric circulation and hydrological cycle over the tropical Indian region.

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: Low level jet intensification by mineral dust aerosols Annales Geophysicae, 31, 625-632, 2013 Author(s): O. Alizadeh Choobari, P. Zawar-Reza, and A. Sturman Modification of the intensity of a low level jet (LLJ) and near-surface wind speed by mineral dust is important as it has implications for dust emission and its long-range transport. Using the Weather Research and Forecasting with Chemistry (WRF/Chem) regional model, it is shown that direct radiative forcing by mineral dust reduces temperature in the lower atmosphere, but increases it in the layers aloft. The surface cooling is shown to be associated with a reduction of turbulent kinetic energy (TKE) and hence vertical mixing of horizontal momentum. Changes in the vertical profile of temperature over the regions that are under the influence of a LLJ are shown to result in an intensification of the LLJ and near-surface wind speed, but a decrease of winds aloft. These changes in the wind speed profile differ from results of previous research which suggested a decrease of wind speed in the lower atmosphere and its increase in the upper boundary layer.

Description: Parametric validations of analytical lifetime estimates for radiation belt electron diffusion by whistler waves Annales Geophysicae, 31, 599-624, 2013 Author(s): A. V. Artemyev, D. Mourenas, O. V. Agapitov, and V. V. Krasnoselskikh The lifetimes of electrons trapped in Earth's radiation belts can be calculated from quasi-linear pitch-angle diffusion by whistler-mode waves, provided that their frequency spectrum is broad enough and/or their average amplitude is not too large. Extensive comparisons between improved analytical lifetime estimates and full numerical calculations have been performed in a broad parameter range representative of a large part of the magnetosphere from L ~ 2 to 6. The effects of observed very oblique whistler waves are taken into account in both numerical and analytical calculations. Analytical lifetimes (and pitch-angle diffusion coefficients) are found to be in good agreement with full numerical calculations based on CRRES and Cluster hiss and lightning-generated wave measurements inside the plasmasphere and Cluster lower-band chorus waves measurements in the outer belt for electron energies ranging from 100 keV to 5 MeV. Comparisons with lifetimes recently obtained from electron flux measurements on SAMPEX, SCATHA, SAC-C and DEMETER also show reasonable agreement.

Description: Coherent amplitude modulation of electron-beam-driven Langmuir waves Annales Geophysicae, 31, 633-638, 2013 Author(s): K. Baumgärtel A linear approach to the phenomenon of irregular amplitude modulation of beam-driven Langmuir waves, developed in a previous paper, is extended to explain periodic modulation as well. It comes about by beating of the fastest growing mode of the instability with beam-aligned plasma oscillations. They are naturally generated in a uniform domain of beam–plasma interaction prior to the onset of the instability. Particle-in-cell (PIC) simulations support the results of the linear analysis.

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: Seasonal variation of the ion upflow in the topside ionosphere during SAPS (subauroral polarization stream) periods Annales Geophysicae, 31, 1521-1534, 2013 Author(s): H. Wang and H. Lühr A statistical study has been performed by using two years of DMSP (Defense Meteorological Satellite Program) plasma observations to investigate the seasonal effect of SAPS (subauroral polarization stream) on the ion upflow in the duskside ionosphere of the Northern Hemisphere. There are obvious upflows occurring in the topside ionosphere around the SAPS region, exceeding 200 m s −1 at winter solstice, indicating an important relationship between SAPS and the local plasma upward motion. Both SAPS and ion upward velocities show similar seasonal variations, largest in winter and smallest in summer, irrespective of geomagnetic activity. A good correlation is found and a linear relationship is derived between SAPS and the ion upflow velocities. During December solstice the average upflow flux can reach about 2 × 10 8 cm −2 s −1 for more disturbed periods, which is comparable to the typical upflow flux in the dayside cusp region. The depression of the ion temperatures around the peak SAPS region can be understood in terms of the adiabatic cooling. The hot ion cools down when expanding into the low ion concentration region. The electron temperature elevates around the SAPS region because of the reduced Coulomb cooling in the low ion density region. Both the changes of ion and electron temperatures are larger in winter than in summer, however, for Kp 〈 4 the electron temperatures are almost seasonably independent. The present work highlights the important role of the SAPS-related frictional heating at mid-latitudes on the local formation of the strong upward flow, which might provide a direct ionospheric ion source for the ring current and plasmasphere in the duskside sector.

Description: Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability Annales Geophysicae, 31, 1535-1541, 2013 Author(s): K.-I. Nishikawa, P. Hardee, B. Zhang, I. Duţan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, and D. H. Hartmann We have investigated the generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin–Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of the kinetic Kelvin–Helmholtz instability (KKHI) of our jet-sheath configuration is slightly different, even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field E z , perpendicular to the flow boundary, and the magnetic field B y , transverse to the flow direction. After the B y component is excited, an induced electric field E x , parallel to the flow direction, becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios m i / m e = 1836 and m i / m e = 20 are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case (γ j = 1.5) is larger than for a relativistic jet case (γ j = 15).

Description: Ionospheric response to total solar eclipse of 22 July 2009 in different Indian regions Annales Geophysicae, 31, 1549-1558, 2013 Author(s): S. Kumar, A. K. Singh, and R. P. Singh The variability of ionospheric response to the total solar eclipse of 22 July 2009 has been studied analyzing the GPS data recorded at the four Indian low-latitude stations Varanasi (100% obscuration), Kanpur (95% obscuration), Hyderabad (84% obscuration) and Bangalore (72% obscuration). The retrieved ionospheric vertical total electron content (VTEC) shows a significant reduction (reflected by all PRNs (satellites) at all stations) with a maximum of 48% at Varanasi (PRN 14), which decreases to 30% at Bangalore (PRN 14). Data from PRN 31 show a maximum of 54% at Kanpur and 26% at Hyderabad. The maximum decrement in VTEC occurs some time (2–15 min) after the maximum obscuration. The reduction in VTEC compared to the quiet mean VTEC depends on latitude as well as longitude, which also depends on the location of the satellite with respect to the solar eclipse path. The amount of reduction in VTEC decreases as the present obscuration decreases, which is directly related to the electron production by the photoionization process. The analysis of electron density height profile derived from the COSMIC (Constellation Observing System for Meteorology, Ionosphere & Climate) satellite over the Indian region shows significant reduction from 100 km altitude up to 800 km altitude with a maximum of 48% at 360 km altitude. The oscillatory nature in total electron content data at all stations is observed with different wave periods lying between 40 and 120 min, which are attributed to gravity wave effects generated in the lower atmosphere during the total solar eclipse.

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: Application of nonlinear autoregressive moving average exogenous input models to geospace: advances in understanding and space weather forecasts Annales Geophysicae, 31, 1579-1589, 2013 Author(s): R. J. Boynton, M. A. Balikhin, S. A. Billings, and O. A. Amariutei The nonlinear autoregressive moving average with exogenous inputs (NARMAX) system identification technique is applied to various aspects of the magnetospheres dynamics. It is shown, from an example system, how the inputs to a system can be found from the error reduction ratio (ERR) analysis, a key concept of the NARMAX approach. The application of the NARMAX approach to the Dst (disturbance storm time) index and the electron fluxes at geostationary Earth orbit (GEO) are reviewed, revealing new insight into the physics of the system. The review of studies into the Dst index illustrate how the NARMAX approach is able to find a coupling function for the Dst index from data, which was then analytically justified from first principles. While the review of the electron flux demonstrates how NARMAX is able to reveal new insight into the physics of the acceleration and loss processes within the radiation belt.

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: Support vector machines for TEC seismo-ionospheric anomalies detection Annales Geophysicae, 31, 173-186, 2013 Author(s): M. Akhoondzadeh Using time series prediction methods, it is possible to pursue the behaviors of earthquake precursors in the future and to announce early warnings when the differences between the predicted value and the observed value exceed the predefined threshold value. Support Vector Machines (SVMs) are widely used due to their many advantages for classification and regression tasks. This study is concerned with investigating the Total Electron Content (TEC) time series by using a SVM to detect seismo-ionospheric anomalous variations induced by the three powerful earthquakes of Tohoku (11 March 2011), Haiti (12 January 2010) and Samoa (29 September 2009). The duration of TEC time series dataset is 49, 46 and 71 days, for Tohoku, Haiti and Samoa earthquakes, respectively, with each at time resolution of 2 h. In the case of Tohoku earthquake, the results show that the difference between the predicted value obtained from the SVM method and the observed value reaches the maximum value (i.e., 129.31 TECU) at earthquake time in a period of high geomagnetic activities. The SVM method detected a considerable number of anomalous occurrences 1 and 2 days prior to the Haiti earthquake and also 1 and 5 days before the Samoa earthquake in a period of low geomagnetic activities. In order to show that the method is acting sensibly with regard to the results extracted during nonevent and event TEC data, i.e., to perform some null-hypothesis tests in which the methods would also be calibrated, the same period of data from the previous year of the Samoa earthquake date has been taken into the account. Further to this, in this study, the detected TEC anomalies using the SVM method were compared to the previous results (Akhoondzadeh and Saradjian, 2011; Akhoondzadeh, 2012) obtained from the mean, median, wavelet and Kalman filter methods. The SVM detected anomalies are similar to those detected using the previous methods. It can be concluded that SVM can be a suitable learning method to detect the novelty changes of a nonlinear time series such as variations of earthquake precursors.

Description: Magnetosheath dynamic pressure enhancements: occurrence and typical properties Annales Geophysicae, 31, 319-331, 2013 Author(s): M. O. Archer and T. S. Horbury The first comprehensive statistical study of large-amplitude (〉 100%) transient enhancements of the magnetosheath dynamic pressure reveals events of up to ~ 15 times the ambient dynamic pressure with durations up to 3 min and an average duration of around 30 s, predominantly downstream of the quasi-parallel shock. The dynamic pressure transients are most often dominated by velocity increases along with a small fractional increase in the density, though the velocity is generally only deflected by a few degrees. Superposed wavelet transforms of the magnetic field show that, whilst most enhancements exhibit changes in the magnetosheath magnetic field, the majority are not associated with changes in the Interplanetary Magnetic Field (IMF). However, there is a minority of enhancements that do appear to be associated with solar wind discontinuities which cannot be explained simply by random events. In general, it is found that during periods of magnetosheath dynamic pressure enhancements the IMF is steadier than usual. This suggests that a stable foreshock and hence foreshock structures or processes may be important in the generation of the majority of magnetosheath dynamic pressure enhancements.

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: Response of polar mesosphere summer echoes to geomagnetic disturbances in the Southern and Northern Hemispheres: the importance of nitric oxide Annales Geophysicae, 31, 333-347, 2013 Author(s): S. Kirkwood, E. Belova, P. Dalin, M. Mihalikova, D. Mikhaylova, D. Murtagh, H. Nilsson, K. Satheesan, J. Urban, and I. Wolf The relationship between polar mesosphere summer echoes (PMSE) and geomagnetic disturbances (represented by magnetic K indices) is examined. Calibrated PMSE reflectivities for the period May 2006–February 2012 are used from two 52.0/54.5 MHz radars located in Arctic Sweden (68° N, geomagnetic latitude 65°) and at two different sites in Queen Maud Land, Antarctica (73°/72° S, geomagnetic latitudes 62°/63°). In both the Northern Hemisphere (NH) and the Southern Hemisphere (SH) there is a strong increase in mean PMSE reflectivity between quiet and disturbed geomagnetic conditions. Mean volume reflectivities are slightly lower at the SH locations compared to the NH, but the position of the peak in the lognormal distribution of PMSE reflectivities is close to the same at both NH and SH locations, and varies only slightly with magnetic disturbance level. Differences between the sites, and between geomagnetic disturbance levels, are primarily due to differences in the high-reflectivity tail of the distribution. PMSE occurrence rates are essentially the same at both NH and SH locations during most of the PMSE season when a sufficiently low detection threshold is used so that the peak in the lognormal distribution is included. When the local-time dependence of the PMSE response to geomagnetic disturbance level is considered, the response in the NH is found to be immediate at most local times, but delayed by several hours in the afternoon sector and absent in the early evening. At the SH sites, at lower magnetic latitude, there is a delayed response (by several hours) at almost all local times. At the NH (auroral zone) site, the dependence on magnetic disturbance is highest during evening-to-morning hours. At the SH (sub-auroral) sites the response to magnetic disturbance is weaker but persists throughout the day. While the immediate response to magnetic activity can be qualitatively explained by changes in electron density resulting from energetic particle precipitation, the delayed response can largely be explained by changes in nitric oxide concentrations. Observations of nitric oxide concentration at PMSE heights by the Odin satellite support this hypothesis. Sensitivity to geomagnetic disturbances, including nitric oxide produced during these disturbances, can explain previously reported differences between sites in the auroral zone and those at higher or lower magnetic latitudes. The several-day lifetime of nitric oxide can also explain earlier reported discrepancies between high correlations for average conditions (year-by-year PMSE reflectivities and K indices) and low correlations for minute-to-day timescales.

Description: Observations of poleward-propagating large-scale traveling ionospheric disturbances in southern China Annales Geophysicae, 31, 377-385, 2013 Author(s): F. Ding, W. Wan, B. Ning, B. Zhao, Q. Li, Y. Wang, L. Hu, R. Zhang, and B. Xiong We report here on two cases of poleward-propagating large-scale traveling ionospheric disturbances (LSTIDs) in China during a medium-scale storm between 27 May and 1 June 2011. The observations were conducted by making use of the Global Positioning System network and ionosondes in China and Southeast Asia. One northeastward-propagating LSTID occurred on the morning of 30 May, while the other was observed during the nighttime of 1 June. Both poleward-traveling LSTIDs occurred during the storm's recovery phase in southern China's low-latitude region (geomagnetic latitude ~ 7.3–24° N) and experienced severe dissipation during their propagation from south to north. Although the initial relative amplitude of the nighttime LSTID was ~ 60% larger than that of the morning event, the nighttime event dissipated more quickly than the morning event because of a strong nighttime enhancement in background total electronic content (TEC) during storm time, which led to strong ion-drag dissipation during the evening. The poleward-propagating LSTIDs exhibit a narrower latitudinal range, a smaller amplitude, and a slightly higher elevation compared with the equatorward-moving LSTIDs observed in the same region. Given these features, the poleward-propagating LSTIDs were likely excited by some local source near southern China. Excitation of secondary LSTIDs during the dissipation of some primary medium-scale disturbances from the lower atmosphere is a possible mechanism.

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: A new method for solving the MHD equations in the magnetosheath Annales Geophysicae, 31, 419-437, 2013 Author(s): C. Nabert, K.-H. Glassmeier, and F. Plaschke We present a new analytical method to derive steady-state magnetohydrodynamic (MHD) solutions of the magnetosheath in different levels of approximation. With this method, we calculate the magnetosheath's density, velocity, and magnetic field distribution as well as its geometry. Thereby, the solution depends on the geomagnetic dipole moment and solar wind conditions only. To simplify the representation, we restrict our model to northward IMF with the solar wind flow along the stagnation streamline. The sheath's geometry, with its boundaries, bow shock and magnetopause, is determined self-consistently. Our model is stationary and time relaxation has not to be considered as in global MHD simulations. Our method uses series expansion to transfer the MHD equations into a new set of ordinary differential equations. The number of equations is related to the level of approximation considered including different physical processes. These equations can be solved numerically; however, an analytical approach for the lowest-order approximation is also presented. This yields explicit expressions, not only for the flow and field variations but also for the magnetosheath thickness, depending on the solar wind parameters. Results are compared to THEMIS data and offer a detailed explanation of, e.g., the pile-up process and the corresponding plasma depletion layer, the bow shock and magnetopause geometry, the magnetosheath thickness, and the flow deceleration.

Description: Investigation of a mesospheric bore event over northern China Annales Geophysicae, 31, 409-418, 2013 Author(s): Q. Li, J. Xu, J. Yue, X. Liu, W. Yuan, B. Ning, S. Guan, and J. P. Younger A mesospheric bore event was observed using an OH all-sky airglow imager (ASAI) at Xinglong (40.2° N, 117.4° E), in northern China, on the night of 8–9 January 2011. Simultaneous observations by a Doppler meteor radar, a broadband sodium lidar, and TIMED/SABER OH intensity and temperature measurements are used to investigate the characteristics and environment of the bore propagation and the possible relations with the Na density perturbations. The bore propagated from northeast to southwest and divided the sky into bright and dark halves. The calculations show that the bore has an average phase velocity of 68 m s −1 . The crests following the bore have a horizontal wavelength of ~ 22 km. These parameters are consistent with the hydraulic jump theory proposed by Dewan and Picard, as well as the previous bore reports. Simultaneous wind measurements from the Doppler meteor radar at Shisanling (40.3° N, 116.2° E) and temperature data from SABER on board the TIMED satellite are used to characterize the propagating environment of the bore. The result shows that a thermal-Doppler duct exists near the OH layer that supports the horizontal propagation of the bore. Simultaneous Na lidar observations at Yanqing (40.4° N, 116.0° E) suggest that there is a downward displacement of Na density during the passage of the mesospheric bore event.

Description: Thermal conductivity of snow measured by three independent methods and anisotropy considerations The Cryosphere, 7, 217-227, 2013 Author(s): F. Riche and M. Schneebeli The thermal conductivity of snow determines the temperature gradient, and by this, it has a direct effect on the rate of snow metamorphism. It is therefore a key property of snow. However, thermal conductivities measured with the transient needle probe and the steady-state, heat flux plate differ. In addition, the anisotropy of thermal conductivity plays an important role in the accuracy of thermal conductivity measurements. In this study, we investigated three independent methods to measure snow thermal conductivity and its anisotropy: a needle probe with a long heating time, a guarded heat flux plate, and direct numerical simulation at the microstructural level of the pore and ice structure. The three methods were applied to identical snow samples. We analyzed the consistency and the difference between these methods. As already shown in former studies, we observed a distinct difference between the anisotropy of thermal conductivity in small rounded grains and in depth hoar. Indeed, the anisotropy between vertical and horizontal thermal conductivity components ranges between 0.5–2. This can cause a difference in thermal conductivity measurements carried out with needle probes of up to –25 % to +25 % if the thermal conductivity is calculated only from a horizontally inserted needle probe. Based on our measurements and the comparison of the three methods studied here, the direct numerical simulation is the most reliable method, as the tensorial components of the thermal conductivity can be calculated and the corresponding microstructure is precisely known.

Description: The ultra-fast Kelvin waves in the equatorial ionosphere: observations and modeling Annales Geophysicae, 31, 209-215, 2013 Author(s): A. N. Onohara, I. S. Batista, and H. Takahashi The main purpose of this study is to investigate the vertical coupling between the mesosphere and lower thermosphere (MLT) region and the ionosphere through ultra-fast Kelvin (UFK) waves in the equatorial atmosphere. The effect of UFK waves on the ionospheric parameters was estimated using an ionospheric model which calculates electrostatic potential in the E-region and solves coupled electrodynamics of the equatorial ionosphere in the E- and F-regions. The UFK wave was observed in the South American equatorial region during February–March 2005. The MLT wind data obtained by meteor radar at São João do Cariri (7.5° S, 37.5° W) and ionospheric F-layer bottom height ( h 'F) observed by ionosonde at Fortaleza (3.9° S; 38.4° W) were used in order to calculate the wave characteristics and amplitude of oscillation. The simulation results showed that the combined electrodynamical effect of tides and UFK waves in the MLT region could explain the oscillations observed in the ionospheric parameters.

Description: Payload charging events in the mesosphere and their impact on Langmuir type electric probes Annales Geophysicae, 31, 187-196, 2013 Author(s): T. A. Bekkeng, A. Barjatya, U.-P. Hoppe, A. Pedersen, J. I. Moen, M. Friedrich, and M. Rapp Three sounding rockets were launched from Andøya Rocket Range in the ECOMA campaign in December 2010. The aim was to study the evolution of meteoric smoke particles during a major meteor shower. Of the various instruments onboard the rocket payload, this paper presents the data from a multi-Needle Langmuir Probe (m-NLP) and a charged dust detector. The payload floating potential, as observed using the m-NLP instrument, shows charging events on two of the three flights. These charging events cannot be explained using a simple charging model, and have implications towards the use of fixed bias Langmuir probes on sounding rockets investigating mesospheric altitudes. We show that for a reliable use of a single fixed bias Langmuir probe as a high spatial resolution relative density measurement, each payload should also carry an additional instrument to measure payload floating potential, and an instrument that is immune to spacecraft charging and measures absolute plasma density.

Description: Brief communication "Important role of the mid-tropospheric atmospheric circulation in the recent surface melt increase over the Greenland ice sheet" The Cryosphere, 7, 241-248, 2013 Author(s): X. Fettweis, E. Hanna, C. Lang, A. Belleflamme, M. Erpicum, and H. Gallée Since 2007, there has been a series of surface melt records over the Greenland ice sheet (GrIS), continuing the trend towards increased melt observed since the end of the 1990's. The last two decades are characterized by an increase of negative phases of the North Atlantic Oscillation (NAO) favouring warmer and drier summers than normal over GrIS. In this context, we use a circulation type classification based on daily 500 hPa geopotential height to evaluate the role of atmospheric dynamics in this surface melt acceleration for the last two decades. Due to the lack of direct observations, the interannual melt variability is gauged here by the summer (June–July–August) mean temperature from reanalyses at 700 hPa over Greenland; analogous atmospheric circulations in the past show that ~70% of the 1993–2012 warming at 700 hPa over Greenland has been driven by changes in the atmospheric flow frequencies. Indeed, the occurrence of anticyclones centred over the GrIS at the surface and at 500 hPa has doubled since the end of 1990's, which induces more frequent southerly warm air advection along the western Greenland coast and over the neighbouring Canadian Arctic Archipelago (CAA). These changes in the NAO modes explain also why no significant warming has been observed these last summers over Svalbard, where northerly atmospheric flows are twice as frequent as before. Therefore, the recent warmer summers over GrIS and CAA cannot be considered as a long-term climate warming but are more a consequence of NAO variability affecting atmospheric heat transport. Although no global model from the CMIP5 database projects subsequent significant changes in NAO through this century, we cannot exclude the possibility that the observed NAO changes are due to global warming.

Description: Restoring mass conservation to shallow ice flow models over complex terrain The Cryosphere, 7, 229-240, 2013 Author(s): A. H. Jarosch, C. G. Schoof, and F. S. Anslow Numerical simulation of glacier dynamics in mountainous regions using zero-order, shallow ice models is desirable for computational efficiency so as to allow broad coverage. However, these models present several difficulties when applied to complex terrain. One such problem arises where steep terrain can spuriously lead to large ice fluxes that remove more mass from a grid cell than it originally contains, leading to mass conservation being violated. This paper describes a vertically integrated, shallow ice model using a second-order flux-limiting spatial discretization scheme that enforces mass conservation. An exact solution to ice flow over a bedrock step is derived for a given mass balance forcing as a benchmark to evaluate the model performance in such a difficult setting. This benchmark should serve as a useful test for modellers interested in simulating glaciers over complex terrain.

Description: Borehole temperatures reveal a changed energy budget at Mill Island, East Antarctica, over recent decades The Cryosphere, 7, 263-273, 2013 Author(s): J. L. Roberts, A. D. Moy, T. D. van Ommen, M. A. J. Curran, A. P. Worby, I. D. Goodwin, and M. Inoue A borehole temperature record from the Mill Island (East Antarctica) icecap reveals a large surface warming signal manifested as a 0.75 K temperature difference over the approximate 100 m depth in the zone of zero annual amplitude below the seasonally varying zone. The temperature profile shows a break in gradient around 49 m depth, which we model with inverse numerical simulations, indicating that surface warming started around the austral summer of 1980/81 AD ±5 yr. This warming of approximately 0.37 K per decade is consistent with trends seen in both instrumental and other reconstructions for Antarctica and, therefore, suggests that regional- rather than local-scale processes are largely responsible. Alteration of the surface energy budget arising from changes in radiation balances due to local cloud, the amount of liquid deposition and local air temperatures associated with altered air/sea exchanges also potentially plays a role at this location due to the proximity of the Shackleton Ice Shelf and sea-ice zone.

Description: Analysis of ice phenology of lakes on the Tibetan Plateau from MODIS data The Cryosphere, 7, 287-301, 2013 Author(s): J. Kropáček, F. Maussion, F. Chen, S. Hoerz, and V. Hochschild The Tibetan Plateau includes a large system of endorheic (closed basin) lakes. Lake ice phenology, i.e. the timing of freeze-up and break-up and the duration of the ice cover may provide valuable information about climate variations in this region. The ice phenology of 59 large lakes on the Tibetan Plateau was derived from Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day composite data for the period from 2001 to 2010. Ice cover duration appears to have a high variability in the studied region due to both climatic and local factors. Mean values for the duration of ice cover were calculated for three groups of lakes defined by clustering, resulting in relatively compact geographic regions. In each group several lakes showed anomalies in ice cover duration in the studied period. Possible reasons for such anomalous behaviour are discussed. Furthermore, many lakes do not freeze up completely during some seasons. This was confirmed by inspection of high resolution optical data. Mild winter seasons, large water volume and/or high salinity are the most likely explanations. Trends in the ice cover duration derived by linear regression for all the studied lakes show a high variation in space. A correlation of ice phenology variables with parameters describing climatic and local conditions showed a high thermal dependency of the ice regime. It appears that the freeze-up tends to be more thermally determined than break-up for the studied lakes.

Description: Electric potential differences across auroral generator interfaces Annales Geophysicae, 31, 251-261, 2013 Author(s): J. De Keyser and M. Echim Strong localized high-altitude auroral electric fields, such as those observed by Cluster, are often associated with magnetospheric interfaces. The type of high-altitude electric field profile (monopolar, bipolar, or more complicated) depends on the properties of the plasmas on either side of the interface, as well as on the total electric potential difference across the structure. The present paper explores the role of this cross-field electric potential difference in the situation where the interface is a tangential discontinuity. A self-consistent Vlasov description is used to determine the equilibrium configuration for different values of the transverse potential difference. A major observation is that there exist limits to the potential difference, beyond which no equilibrium configuration of the interface can be sustained. It is further demonstrated how the plasma densities and temperatures affect the type of electric field profile in the transition, with monopolar electric fields appearing primarily when the temperature contrast is large. These findings strongly support the observed association of monopolar fields with the plasma sheet boundary. The role of shear flow tangent to the interface is also examined.

Description: Variability of ionospheric TEC during solar and geomagnetic minima (2008 and 2009): external high speed stream drivers Annales Geophysicae, 31, 263-276, 2013 Author(s): O. P. Verkhoglyadova, B. T. Tsurutani, A. J. Mannucci, M. G. Mlynczak, L. A. Hunt, and T. Runge We study solar wind–ionosphere coupling through the late declining phase/solar minimum and geomagnetic minimum phases during the last solar cycle (SC23) – 2008 and 2009. This interval was characterized by sequences of high-speed solar wind streams (HSSs). The concomitant geomagnetic response was moderate geomagnetic storms and high-intensity, long-duration continuous auroral activity (HILDCAA) events. The JPL Global Ionospheric Map (GIM) software and the GPS total electron content (TEC) database were used to calculate the vertical TEC (VTEC) and estimate daily averaged values in separate latitude and local time ranges. Our results show distinct low- and mid-latitude VTEC responses to HSSs during this interval, with the low-latitude daytime daily averaged values increasing by up to 33 TECU (annual average of ~20 TECU) near local noon (12:00 to 14:00 LT) in 2008. In 2009 during the minimum geomagnetic activity (MGA) interval, the response to HSSs was a maximum of ~30 TECU increases with a slightly lower average value than in 2008. There was a weak nighttime ionospheric response to the HSSs. A well-studied solar cycle declining phase interval, 10–22 October 2003, was analyzed for comparative purposes, with daytime low-latitude VTEC peak values of up to ~58 TECU (event average of ~55 TECU). The ionospheric VTEC changes during 2008–2009 were similar but ~60% less intense on average. There is an evidence of correlations of filtered daily averaged VTEC data with Ap index and solar wind speed. We use the infrared NO and CO 2 emission data obtained with SABER on TIMED as a proxy for the radiation balance of the thermosphere. It is shown that infrared emissions increase during HSS events possibly due to increased energy input into the auroral region associated with HILDCAAs. The 2008–2009 HSS intervals were ~85% less intense than the 2003 early declining phase event, with annual averages of daily infrared NO emission power of ~ 3.3 × 10 10 W and 2.7 × 10 10 W in 2008 and 2009, respectively. The roles of disturbance dynamos caused by high-latitude winds (due to particle precipitation and Joule heating in the auroral zones) and of prompt penetrating electric fields (PPEFs) in the solar wind–ionosphere coupling during these intervals are discussed. A correlation between geoeffective interplanetary electric field components and HSS intervals is shown. Both PPEF and disturbance dynamo mechanisms could play important roles in solar wind–ionosphere coupling during prolonged (up to days) external driving within HILDCAA intervals.

Description: Paleo ice flow and subglacial meltwater dynamics in Pine Island Bay, West Antarctica The Cryosphere, 7, 249-262, 2013 Author(s): F. O. Nitsche, K. Gohl, R. D. Larter, C.-D. Hillenbrand, G. Kuhn, J. A. Smith, S. Jacobs, J. B. Anderson, and M. Jakobsson Increasing evidence for an elaborate subglacial drainage network underneath modern Antarctic ice sheets suggests that basal meltwater has an important influence on ice stream flow. Swath bathymetry surveys from previously glaciated continental margins display morphological features indicative of subglacial meltwater flow in inner shelf areas of some paleo ice stream troughs. Over the last few years several expeditions to the eastern Amundsen Sea embayment (West Antarctica) have investigated the paleo ice streams that extended from the Pine Island and Thwaites glaciers. A compilation of high-resolution swath bathymetry data from inner Pine Island Bay reveals details of a rough seabed topography including several deep channels that connect a series of basins. This complex basin and channel network is indicative of meltwater flow beneath the paleo-Pine Island and Thwaites ice streams, along with substantial subglacial water inflow from the east. This meltwater could have enhanced ice flow over the rough bedrock topography. Meltwater features diminish with the onset of linear features north of the basins. Similar features have previously been observed in several other areas, including the Dotson-Getz Trough (western Amundsen Sea embayment) and Marguerite Bay (SW Antarctic Peninsula), suggesting that these features may be widespread around the Antarctic margin and that subglacial meltwater drainage played a major role in past ice-sheet dynamics.

Description: Case study of stratospheric gravity waves of convective origin over Arctic Scandinavia – VHF radar observations and numerical modelling Annales Geophysicae, 31, 239-250, 2013 Author(s): A. Réchou, J. Arnault, P. Dalin, and S. Kirkwood Orography is a well-known source of gravity and inertia-gravity waves in the atmosphere. Other sources, such as convection, are also known to be potentially important but the large amplitude of orographic waves over Scandinavia has generally precluded the possibility to study such other sources experimentally in this region. In order to better understand the origin of stratospheric gravity waves observed by the VHF radar ESRAD (Esrange MST radar) over Kiruna, in Arctic Sweden (67.88° N, 21.10° E), observations have been compared to simulations made using the Weather Research and Forecasting model (WRF) with and without the effects of orography and clouds. This case study concerns gravity waves observed from 00:00 UTC on 18 February to 12:00 UTC on 20 February 2007. We focus on the wave signatures in the static stability field and vertical wind deduced from the simulations and from the observations as these are the parameters which are provided by the observations with the best height coverage. As is common at this site, orographic gravity waves were produced over the Scandinavian mountains and observed by the radar. However, at the same time, southward propagation of fronts in the Barents Sea created short-period waves which propagated into the stratosphere and were transported, embedded in the cyclonic winds, over the radar site.

Description: Long-term trends in the ionospheric F2 region with different solar activity indices Annales Geophysicae, 31, 291-303, 2013 Author(s): J. Mielich and J. Bremer A new comprehensive data collection by Damboldt and Suessmann (2012a) with monthly fo F2 and M(3000)F2 median values is an excellent basis for the derivation of long-term trends in the ionospheric F2 region. Ionospheric trends have been derived only for stations with data series of at least 22 years (124 stations with fo F2 data and 113 stations with M(3000)F2 data) using a twofold regression analysis depending on solar and geomagnetic activity. Three main results have been derived: Firstly, it could be shown that the solar 10.7 cm radio flux F10.7 is a better index for the description of the solar activity than the relative solar sunspot number R as well as the solar EUV proxy E10.7. Secondly, the global mean fo F2 and hm F2 trends derived for the interval between 1948 and 2006 are in surprisingly good agreement with model calculations of an increasing atmospheric greenhouse effect (Rishbeth and Roble, 1992). Thirdly, during the years 2007 until 2009, the hm F2 values and to a smaller amount the fo F2 values strongly decrease. The reason for this effect is a reduction of the thermospheric density and ionization due to a markedly reduced solar EUV irradiation and extremely small geomagnetic activity during the solar cycle 23/24 minimum.

Description: Simultaneous measurement of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow at Kolhapur (17° N), India Annales Geophysicae, 31, 197-208, 2013 Author(s): N. Parihar, A. Taori, S. Gurubaran, and G. K. Mukherjee Near-simultaneous measurements of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow were carried out at Kolhapur (17° N), India during February–March 2007. Atmospheric temperatures around 87 and 94 km were derived from the knowledge of intensity measurements of spectral features OH (6, 2) Meinel Band and O 2 Atmospheric Band, respectively. An account of the behaviour of derived temperatures has been presented. The nocturnal behaviour of OH and O 2 temperatures is governed by the waves of tidal origin, whereas the signatures of planetary wave-like oscillations is noted in the night-to-night variation of two temperatures. This is probably the first report of planetary waves observed in nightglow temperatures from the Indian subcontinent.

Description: Refilling process in the plasmasphere: a 3-D statistical characterization based on Cluster density observations Annales Geophysicae, 31, 217-237, 2013 Author(s): G. Lointier, F. Darrouzet, P. M. E. Décréau, X. Vallières, S. Kougblénou, J. G. Trotignon, and J.-L. Rauch The Cluster mission offers an excellent opportunity to investigate the evolution of the plasma population in a large part of the inner magnetosphere, explored near its orbit's perigee, over a complete solar cycle. The WHISPER sounder, on board each satellite of the mission, is particularly suitable to study the electron density in this region, between 0.2 and 80 cm −3 . Compiling WHISPER observations during 1339 perigee passes distributed over more than three years of the Cluster mission, we present first results of a statistical analysis dedicated to the study of the electron density morphology and dynamics along and across magnetic field lines between L = 2 and L = 10. In this study, we examine a specific topic: the refilling of the plasmasphere and trough regions during extended periods of quiet magnetic conditions. To do so, we survey the evolution of the ap index during the days preceding each perigee crossing and sort out electron density profiles along the orbit according to three classes, namely after respectively less than 2 days, between 2 and 4 days, and more than 4 days of quiet magnetic conditions (ap ≤ 15 nT) following an active episode (ap 〉 15 nT). This leads to three independent data subsets. Comparisons between density distributions in the 3-D plasmasphere and trough regions at the three stages of quiet magnetosphere provide novel views about the distribution of matter inside the inner magnetosphere during several days of low activity. Clear signatures of a refilling process inside an expended plasmasphere in formation are noted. A plasmapause-like boundary, at L ~ 6 for all MLT sectors, is formed after 3 to 4 days and expends somewhat further after that. In the outer part of the plasmasphere ( L ~ 8), latitudinal profiles of median density values vary essentially according to the MLT sector considered rather than according to the refilling duration. The shape of these density profiles indicates that magnetic flux tubes are not fully replenished after 6 days of quiet conditions. In addition, the outer plasmasphere in the night and dawn sectors (22:00 to 10:00 MLT range) maintains an overall clear deficit of ionospheric population, when compared to the situation in the noon and dusk sectors (10:00 to 22:00 MLT range).

Description: A recent tipping point in the Arctic sea-ice cover: abrupt and persistent increase in the seasonal cycle since 2007 The Cryosphere, 7, 275-286, 2013 Author(s): V. N. Livina and T. M. Lenton There is ongoing debate over whether Arctic sea ice has already passed a "tipping point", or whether it will do so in the future. Several recent studies argue that the loss of summer sea ice does not involve an irreversible bifurcation, because it is highly reversible in models. However, a broader definition of a "tipping point" also includes other abrupt, non-linear changes that are neither bifurcations nor necessarily irreversible. Examination of satellite data for Arctic sea-ice area reveals an abrupt increase in the amplitude of seasonal variability in 2007 that has persisted since then. We identified this abrupt transition using recently developed methods that can detect multi-modality in time-series data and sometimes forewarn of bifurcations. When removing the mean seasonal cycle (up to 2008) from the satellite data, the residual sea-ice fluctuations switch from uni-modal to multi-modal behaviour around 2007. We originally interpreted this as a bifurcation in which a new lower ice cover attractor appears in deseasonalised fluctuations and is sampled in every summer–autumn from 2007 onwards. However, this interpretation is clearly sensitive to how the seasonal cycle is removed from the raw data, and to the presence of continental land masses restricting winter–spring ice fluctuations. Furthermore, there was no robust early warning signal of critical slowing down prior to the hypothesized bifurcation. Early warning indicators do however show destabilization of the summer–autumn sea-ice cover since 2007. Thus, the bifurcation hypothesis lacks consistent support, but there was an abrupt and persistent increase in the amplitude of the seasonal cycle of Arctic sea-ice cover in 2007, which we describe as a (non-bifurcation) "tipping point". Our statistical methods detect this "tipping point" and its time of onset. We discuss potential geophysical mechanisms behind it, which should be the subject of further work with process-based models.

Description: Spatiotemporal variability analysis of diffuse radiation in China during 1981–2010 Annales Geophysicae, 31, 277-289, 2013 Author(s): X. L. Ren, H. L. He, L. Zhang, L. Zhou, G. R. Yu, and J. W. Fan Solar radiation is the primary driver of terrestrial plant photosynthesis and the diffuse component can enhance canopy light use efficiency (LUE), which in turn influences the carbon balance of terrestrial ecosystems. In this study we calculated the spatial data of diffuse radiation in China from 1981 to 2010, using a radiation decomposition model and spatial interpolation method based on observational data. Furthermore, we explored the spatiotemporal characteristics of diffuse radiation using GIS and trend analysis techniques. The results show the following: (1) The spatial patterns of perennial average of annual diffuse radiation during 1981–2010 are complex and inhomogeneous in China, generally lower in the north and higher in the south and west. The perennial average ranges from 1730.20 to 3064.41 MJ m −2 yr −1 across the whole country. (2) There is an increasing trend of annual diffuse radiation in China from 1981 to 2010 on the whole, with mean increasing amplitude of 7.03 MJ m −2 yr −1 per decade. Whereas a significant downtrend was observed in the first 10 years, distinct anomalies in 1982, 1983, 1991 and 1992 occurred due to the eruptions of El Chinchon and Pinatubo. (3) The spatial distribution of the temporal variability of diffuse radiation showed significant regional heterogeneity in addition to the seasonal differences. Northwestern China has the most evident downtrend, with highest decreasing rate of 6% per decade, while the Tibetan Plateau has the most evident uptrend, with highest increasing rate of up to 9% per decade. Such quantitative spatiotemporal characteristics of diffuse radiation are essential in regional scale modeling of terrestrial carbon dynamics.

Description: Analysis of the snow-atmosphere energy balance during wet-snow instabilities and implications for avalanche prediction The Cryosphere, 7, 205-216, 2013 Author(s): C. Mitterer and J. Schweizer Wet-snow avalanches are notoriously difficult to predict; their formation mechanism is poorly understood since in situ measurements representing the thermal and mechanical evolution are difficult to perform. Instead, air temperature is commonly used as a predictor variable for days with high wet-snow avalanche danger – often with limited success. As melt water is a major driver of wet-snow instability and snow melt depends on the energy input into the snow cover, we computed the energy balance for predicting periods with high wet-snow avalanche activity. The energy balance was partly measured and partly modelled for virtual slopes at different elevations for the aspects south and north using the 1-D snow cover model SNOWPACK. We used measured meteorological variables and computed energy balance and its components to compare wet-snow avalanche days to non-avalanche days for four consecutive winter seasons in the surroundings of Davos, Switzerland. Air temperature, the net shortwave radiation and the energy input integrated over 3 or 5 days showed best results in discriminating event from non-event days. Multivariate statistics, however, revealed that for better predicting avalanche days, information on the cold content of the snowpack is necessary. Wet-snow avalanche activity was closely related to periods when large parts of the snowpack reached an isothermal state (0 °C) and energy input exceeded a maximum value of 200 kJ m −2 in one day, or the 3-day sum of positive energy input was larger than 1.2 MJ m −2 . Prediction accuracy with measured meteorological variables was as good as with computed energy balance parameters, but simulated energy balance variables accounted better for different aspects, slopes and elevations than meteorological data.

Description: Effect of higher-order stress gradients on the centennial mass evolution of the Greenland ice sheet The Cryosphere, 7, 183-199, 2013 Author(s): J. J. Fürst, H. Goelzer, and P. Huybrechts We use a three-dimensional thermo-mechanically coupled model of the Greenland ice sheet to assess the effects of marginal perturbations on volume changes on centennial timescales. The model is designed to allow for five ice dynamic formulations using different approximations to the force balance. The standard model is based on the shallow ice approximation for both ice deformation and basal sliding. A second model version relies on a higher-order Blatter/Pattyn type of core that resolves effects from gradients in longitudinal stresses and transverse horizontal shearing, i.e. membrane-like stresses. Together with three intermediate model versions, these five versions allow for gradually more dynamic feedbacks from membrane stresses. Idealised experiments are conducted on various resolutions to compare the time-dependent response to imposed accelerations at the marine ice front. If such marginal accelerations are to have an appreciable effect on total mass loss on a century timescale, a fast mechanism to transmit such perturbations inland is required. While the forcing is independent of the model version, inclusion of direct horizontal coupling allows the initial speed-up to reach several tens of kilometres inland. Within one century, effects from gradients in membrane stress alter the inland signal propagation and transmit additional dynamic thinning to the ice sheet interior. But the centennial overall volume loss differs only by some percents from the standard model, as the dominant response is a diffusive inland propagation of geometric changes. For the experiments considered, this volume response is even attenuated by direct horizontal coupling. The reason is a faster adjustment of the sliding regime by instant stress transmission in models that account for the effect of membrane stresses. Ultimately, horizontal coupling decreases the reaction time to perturbations at the ice sheet margin. These findings suggest that for modelling the mass evolution of a large-scale ice sheet, effects from diffusive geometric adjustments dominate effects from successively more complete dynamic approaches.

Description: Environmental controls on the thermal structure of alpine glaciers The Cryosphere, 7, 167-182, 2013 Author(s): N. J. Wilson and G. E. Flowers Water entrapped in glacier accumulation zones represents a significant latent heat contribution to the development of thermal structure. It also provides a direct link between glacier environments and thermal regimes. We apply a two-dimensional mechanically-coupled model of heat flow to synthetic glacier geometries in order to explore the environmental controls on flowband thermal structure. We use this model to test the sensitivity of thermal structure to physical and environmental variables and to explore glacier thermal response to environmental changes. In different conditions consistent with a warming climate, mean glacier temperature and the volume of temperate ice may either increase or decrease, depending on the competing effects of elevated meltwater production, reduced accumulation zone extent and thinning firn. For two model reference states that exhibit commonly-observed thermal structures, the fraction of temperate ice is shown to decline with warming air temperatures. Mass balance and aquifer sensitivities play an important role in determining how the englacial thermal regimes of alpine glaciers will adjust in the future.

Description: Comment on "Storming the Bastille: the effect of electric fields on the ionospheric F-layer" by Rishbeth et al. (2010) Annales Geophysicae, 31, 145-150, 2013 Author(s): B. T. Tsurutani, A. J. Mannuccci, O. P. Verkhoglyadova, and G. S. Lakhina No abstract available.

Description: Brief Communication "Expansion of meltwater lakes on the Greenland Ice Sheet" The Cryosphere, 7, 201-204, 2013 Author(s): I. M. Howat, S. de la Peña, J. H. van Angelen, J. T. M. Lenaerts, and M. R. van den Broeke Forty years of satellite imagery reveal that meltwater lakes on the margin of the Greenland Ice Sheet have expanded substantially inland to higher elevations with warming. These lakes are important because they provide a mechanism for bringing water to the ice bed, causing sliding. Inland expansion of lakes could accelerate ice flow by bringing water to previously frozen bed, potentially increasing future rates of mass loss. Increasing lake elevations closely follow the rise of the mass balance equilibrium line over much of the ice sheet, suggesting no physical limit on lake expansion. Data are not yet available to detect a corresponding change in ice flow, and the potential effects of lake expansion on ice sheet dynamics are not included in ice sheet models.

Description: Interferometric swath processing of Cryosat data for glacial ice topography The Cryosphere, 7, 1857-1867, 2013 Author(s): L. Gray, D. Burgess, L. Copland, R. Cullen, N. Galin, R. Hawley, and V. Helm We have derived digital elevation models (DEMs) over the western part of the Devon Ice Cap in Nunavut, Canada, using "swath processing" of interferometric data collected by Cryosat between February 2011 and January 2012. With the standard ESA (European Space Agency) SARIn (synthetic aperture radar interferometry) level 2 (L2) data product, the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time-delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~ 0.5 to ~ 2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna so that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped, allowing the creation of DEMs from a sequence of L1b SARIn Cryosat data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order 1 m with a standard deviation of about 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.

Description: Statistical study of foreshock cavitons Annales Geophysicae, 31, 2163-2178, 2013 Author(s): P. Kajdič, X. Blanco-Cano, N. Omidi, K. Meziane, C. T. Russell, J.-A. Sauvaud, I. Dandouras, and B. Lavraud In this work we perform a statistical analysis of 92 foreshock cavitons observed with the Cluster spacecraft 1 during the period 2001–2006. We analyze time intervals during which the spacecraft was located in the Earth's foreshock with durations longer than 10 min. Together these amount to ~ 50 days. The cavitons are transient structures in the Earth's foreshock. Their main signatures in the data include simultaneous depletions of the magnetic field intensity and plasma density, which are surrounded by a rim of enhanced values of these two quantities. Cavitons form due to nonlinear interaction of transverse and compressive ultra-low frequency (ULF) waves and are therefore always surrounded by intense compressive ULF fluctuations. They are carried by the solar wind towards the bow shock. This work represents the first systematic study of a large sample of foreshock cavitons. We find that cavitons appear for a wide range of solar wind and interplanetary magnetic field conditions and are therefore a common feature upstream of Earth's quasi-parallel bow shock with an average occurrence rate of ~ 2 events per day. We also discuss their observational properties in the context of other known upstream phenomena and show that the cavitons are a distinct structure in the foreshock.

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

Description: Theoretical constraints on the cross-tail width of bursty bulk flows Annales Geophysicae, 31, 2179-2192, 2013 Author(s): C. X. Chen The characteristic cross-tail width of bursty bulk flows (BBFs) in earth's plasma sheet was investigated at two stages of its life, one at its onset, the other when it is fully developed. Equilibrium domains with gradient of magnetic field are constructed. Interchange instability analysis of such domains yields the most unstable mode with the half wave length comparable with the observed cross-tail width of a flow burst and the inverse of growth rate comparable with its duration. The thickness of the plasma sheet for the most unstable mode is also comparable to the width of BBFs in the north–south direction. We found that viscosity, the dimension of the unstable domain, the thickness of the plasma sheet and gradient of the magnetic field together determine the most unstable mode. The ion Larmor radius plays an important role in viscosity as half effective mean free path. For a fully developed flow, however, velocity-caused pressure difference between the leading and trailing sides of a flow burst also plays a role. The equatorial cross section of flow is reshaped and its cross-tail width is changed as well. Representing the surrounding medium with empirical magnetic field and plasma models, the force balance of the fast flow is analyzed. The cross-section area of flow burst is estimated to be one to several square earth radii, and the cross-tail width of fast flow is estimated to be 1 to 3 earth radii, which is consistent with observations of BBFs.

Description: Ray tracing model of the auroral kilometric radiation generation in the 3-D plasma cavity Annales Geophysicae, 31, 1645-1652, 2013 Author(s): T. M. Burinskaya Propagation and amplification of the auroral kilometric radiation (AKR) in a three-dimensional plasma cavity is investigated using the approximation of the geometrical optics, and taking into account both the slightly relativistic electrons propagating inside a cavity and the background cold electrons. It is shown that the global magnetic field inhomogeneity plays a key role in a wave escape from a thin plasma cavity. The main contribution to the AKR spectrum is made by waves initially generated with the component of group velocity directed to the Earth and with the optimum relationship between the wave vector components, controlling the value of the linear grow rate and duration of the ray lifetime inside a source.

Description: Empirical regional models for the short-term forecast of M3000F2 during not quiet geomagnetic conditions over Europe Annales Geophysicae, 31, 1653-1671, 2013 Author(s): M. Pietrella Twelve empirical local models have been developed for the long-term prediction of the ionospheric characteristic M3000F2 , and then used as starting point for the development of a short-term forecasting empirical regional model of M3000F2 under not quiet geomagnetic conditions. Under the assumption that the monthly median measurements of M3000F2 are linearly correlated to the solar activity, a set of regression coefficients were calculated over 12 months and 24 h for each of 12 ionospheric observatories located in the European area, and then used for the long-term prediction of M3000F2 at each station under consideration. Based on the 12 long-term prediction empirical local models of M3000F2 , an empirical regional model for the prediction of the monthly median field of M3000F2 over Europe (indicated as RM_M3000F2 ) was developed. Thanks to the IFELM_foF2 models, which are able to provide short-term forecasts of the critical frequency of the F2 layer ( fo F2 STF ) up to three hours in advance, it was possible to considerer the Brudley–Dudeney algorithm as a function of fo F2 STF to correct RM_M3000F2 and thus obtain an empirical regional model for the short-term forecasting of M3000F2 (indicated as RM_M3000F2_BD ) up to three hours in advance under not quiet geomagnetic conditions. From the long-term predictions of M3000F2 provided by the IRI model, an empirical regional model for the forecast of the monthly median field of M3000F2 over Europe (indicated as IRI_RM_M3000F2 ) was derived. IRI_RM_M3000F2 predictions were modified with the Bradley–Dudeney correction factor, and another empirical regional model for the short-term forecasting of M3000F2 (indicated as IRI_RM_M3000F2_BD ) up to three hours ahead under not quiet geomagnetic conditions was obtained. The main results achieved comparing the performance of RM_M3000F2 , RM_M3000F2_BD , IRI_RM_M3000F2 , and IRI_RM_M3000F2_BD are (1) in the case of moderate geomagnetic activity, the Bradley–Dudeney correction factor does not improve significantly the predictions; (2) under disturbed geomagnetic conditions, the Bradley–Dudeney formula improves the predictions of RM_M3000F2 in the entire European area; (3) in the case of very disturbed geomagnetic conditions, the Bradley–Dudeney algorithm is very effective in improving the performance of IRI_RM_M3000F2 ; (4) under moderate geomagnetic conditions, the long-term prediction maps of M3000F2 generated by RM_M3000F2 can be considered as short-term forecasting maps providing very satisfactory results because quiet geomagnetic conditions are not so diverse from moderate geomagnetic conditions; (5) the forecasting maps originated by RM_M3000F2 , RM_M3000F2_BD , and IRI_RM_M3000F2_BD show some regions where the forecasts are not satisfactory, but also wide sectors where the M3000F2 forecasts quite faithfully match the M3000F2 observations, and therefore RM_M3000F2 , RM_M3000F2_BD , and IRI_RM_M3000F2_BD could be exploited to produce short-term forecasting maps of M3000F2 over Europe up to 3 h in advance.

Description: The sensitivity of flowline models of tidewater glaciers to parameter uncertainty The Cryosphere, 7, 1579-1590, 2013 Author(s): E. M. Enderlin, I. M. Howat, and A. Vieli Depth-integrated (1-D) flowline models have been widely used to simulate fast-flowing tidewater glaciers and predict change because the continuous grounding line tracking, high horizontal resolution, and physically based calving criterion that are essential to realistic modeling of tidewater glaciers can easily be incorporated into the models while maintaining high computational efficiency. As with all models, the values for parameters describing ice rheology and basal friction must be assumed and/or tuned based on observations. For prognostic studies, these parameters are typically tuned so that the glacier matches observed thickness and speeds at an initial state, to which a perturbation is applied. While it is well know that ice flow models are sensitive to these parameters, the sensitivity of tidewater glacier models has not been systematically investigated. Here we investigate the sensitivity of such flowline models of outlet glacier dynamics to uncertainty in three key parameters that influence a glacier's resistive stress components. We find that, within typical observational uncertainty, similar initial (i.e., steady-state) glacier configurations can be produced with substantially different combinations of parameter values, leading to differing transient responses after a perturbation is applied. In cases where the glacier is initially grounded near flotation across a basal over-deepening, as typically observed for rapidly changing glaciers, these differences can be dramatic owing to the threshold of stability imposed by the flotation criterion. The simulated transient response is particularly sensitive to the parameterization of ice rheology: differences in ice temperature of ~ 2 °C can determine whether the glaciers thin to flotation and retreat unstably or remain grounded on a marine shoal. Due to the highly non-linear dependence of tidewater glaciers on model parameters, we recommend that their predictions are accompanied by sensitivity tests that take parameter uncertainty into account.

Description: On the supply of heavy planetary material to the magnetotail of Mercury Annales Geophysicae, 31, 1673-1679, 2013 Author(s): D. C. Delcourt We examine the transport of low-energy heavy ions of planetary origin (O + , Na + , Ca + ) in the magnetosphere of Mercury. We show that, in contrast to Earth, these ions are abruptly energized after ejection into the magnetosphere due to enhanced curvature-related parallel acceleration. Regardless of their mass-to-charge ratio, the parallel speed of these ions is rapidly raised up to ~ 2 V E × B (denoting by V E × B the magnitude of the local E × B drift speed), in a like manner to Fermi-type acceleration by a moving magnetic mirror. This parallel energization is such that ions with very low initial energies (a few tenths of eVs) can overcome gravity and, regardless of species or convection rate, are transported over comparable distances into the nightside magnetosphere. The region of space where these ions reach the magnetotail is found to extend over altitudes similar to those where enhanced densities are noticeable in the MESSENGER data, viz., from ~ 1000 km up to ~ 6000 km in the pre-midnight sector. The observed density enhancements may thus follow from E × B related focusing of planetary material of dayside origin into the magnetotail. Due to the planetary magnetic field offset, an asymmetry is found between drift paths anchored in the Northern and Southern hemispheres, which puts forward a predominant role of heavy material originating in the Northern Hemisphere in populating the innermost region of Mercury's magnetotail.

Description: Evaluating the applicability of the finite element method for modelling of geoelectric fields Annales Geophysicae, 31, 1689-1698, 2013 Author(s): B. Dong, D. W. Danskin, R. J. Pirjola, D. H. Boteler, and Z. Z. Wang Geomagnetically induced currents in power systems are due to space weather events which create geomagnetic disturbances accompanied by electric fields at the surface of the Earth. The purpose of this paper is to evaluate the use of the finite element method (FEM) to calculate the magnetic and electric fields to which long transmission lines of power systems on the Earth are exposed. The well-known technique of FEM is used for the first time to simulate magnetic and electric fields applicable to power systems. Several test cases are modelled and compared with known solutions. It is shown that FEM is an effective modelling technique that can be applied to determine the electric fields which affect power systems. FEM enables an increased capability beyond the traditional methods for modelling electric and magnetic fields with layered earth conductivity structures, as spatially more complex structures can be considered using FEM. As an example results are presented for induction, due to a line current source, in adjacent regions with different layered conductivity structures. The results show the electric field away from the interface is the same as calculated for a single region; however near the interface the electric field is influenced by both regions.

Description: Stereoscopic determination of all-sky altitude map of aurora using two ground-based Nikon DSLR cameras Annales Geophysicae, 31, 1543-1548, 2013 Author(s): R. Kataoka, Y. Miyoshi, K. Shigematsu, D. Hampton, Y. Mori, T. Kubo, A. Yamashita, M. Tanaka, T. Takahei, T. Nakai, H. Miyahara, and K. Shiokawa A new stereoscopic measurement technique is developed to obtain an all-sky altitude map of aurora using two ground-based digital single-lens reflex (DSLR) cameras. Two identical full-color all-sky cameras were set with an 8 km separation across the Chatanika area in Alaska (Poker Flat Research Range and Aurora Borealis Lodge) to find localized emission height with the maximum correlation of the apparent patterns in the localized pixels applying a method of the geographical coordinate transform. It is found that a typical ray structure of discrete aurora shows the broad altitude distribution above 100 km, while a typical patchy structure of pulsating aurora shows the narrow altitude distribution of less than 100 km. Because of its portability and low cost of the DSLR camera systems, the new technique may open a unique opportunity not only for scientists but also for night-sky photographers to complementarily attend the aurora science to potentially form a dense observation network.

Description: A particle based simulation model for glacier dynamics The Cryosphere, 7, 1591-1602, 2013 Author(s): J. A. Åström, T. I. Riikilä, T. Tallinen, T. Zwinger, D. Benn, J. C. Moore, and J. Timonen A particle-based computer simulation model was developed for investigating the dynamics of glaciers. In the model, large ice bodies are made of discrete elastic particles which are bound together by massless elastic beams. These beams can break, which induces brittle behaviour. At loads below fracture, beams may also break and reform with small probabilities to incorporate slowly deforming viscous behaviour in the model. This model has the advantage that it can simulate important physical processes such as ice calving and fracturing in a more realistic way than traditional continuum models. For benchmarking purposes the deformation of an ice block on a slip-free surface was compared to that of a similar block simulated with a Finite Element full-Stokes continuum model. Two simulations were performed: (1) calving of an ice block partially supported in water, similar to a grounded marine glacier terminus, and (2) fracturing of an ice block on an inclined plane of varying basal friction, which could represent transition to fast flow or surging. Despite several approximations, including restriction to two-dimensions and simplified water-ice interaction, the model was able to reproduce the size distributions of the debris observed in calving, which may be approximated by universal scaling laws. On a moderate slope, a large ice block was stable and quiescent as long as there was enough of friction against the substrate. For a critical length of frictional contact, global sliding began, and the model block disintegrated in a manner suggestive of a surging glacier. In this case the fragment size distribution produced was typical of a grinding process.

Description: Pine Island glacier ice shelf melt distributed at kilometre scales The Cryosphere, 7, 1543-1555, 2013 Author(s): P. Dutrieux, D. G. Vaughan, H. F. J. Corr, A. Jenkins, P. R. Holland, I. Joughin, and A. H. Fleming By thinning and accelerating, West Antarctic ice streams are contributing about 10% of the observed global sea level rise. Much of this ice loss is from Pine Island Glacier, which has thinned since at least 1992, driven by changes in ocean heat transport beneath its ice shelf and retreat of the grounding line. Details of the processes driving this change, however, remain largely elusive, hampering our ability to predict the future behaviour of this and similar systems. Here, a Lagrangian methodology is developed to measure oceanic melting of such rapidly advecting ice. High-resolution satellite and airborne observations of ice surface velocity and elevation are used to quantify patterns of basal melt under the Pine Island Glacier ice shelf and the associated adjustments to ice flow. At the broad scale, melt rates of up to 100 m yr −1 occur near the grounding line, reducing to 30 m yr −1 just 20 km downstream. Between 2008 and 2011, basal melting was largely compensated by ice advection, allowing us to estimate an average loss of ice to the ocean of 87 km 3 yr −1 , in close agreement with 2009 oceanographically constrained estimates. At smaller scales, a network of basal channels typically 500 m to 3 km wide is sculpted by concentrated melt, with kilometre-scale anomalies reaching 50% of the broad-scale basal melt. Basal melting enlarges the channels close to the grounding line, but farther downstream melting tends to diminish them. Kilometre-scale variations in melt are a key component of the complex ice–ocean interaction beneath the ice shelf, implying that greater understanding of their effect, or very high resolution models, are required to predict the sea-level contribution of the region.

Description: Cluster as current sheet surveyor in the magnetotail Annales Geophysicae, 31, 1605-1610, 2013 Author(s): Y. Narita, R. Nakamura, and W. Baumjohann A novel analysis technique is presented to estimate the current sheet thickness unambiguously and directly, without associating time series data with spatial structure. The technique is a combination of eigenvalue analysis and minimum variance estimator adapted to Harris current sheet geometry, and needs one-time, four-point magnetic field data as provided by the Cluster spacecraft. Two current sheet parameters, thickness and distance to the spacecraft, can be determined at each time step of the magnetic field measurements. An example is shown from a Cluster magnetotail crossing under quiet magnetospheric conditions, yielding the result that the current sheet thickness is on the scale of the proton gyroradius. The analysis technique can also be used to track the dynamical evolution of the current sheet structure in three dimensions.

Description: Changes in the ultra-low frequency wave field during the precursor phase to the Sichuan earthquake: DEMETER observations Annales Geophysicae, 31, 1597-1603, 2013 Author(s): S. N. Walker, V. Kadirkamanathan, and O. A. Pokhotelov Electromagnetic phenomena observed in association with increases in seismic activity have been studied for several decades. These phenomena are generated during the precursory phases of an earthquake as well as during the main event. Their occurrence during the precursory phases may be used in short-term prediction of a large earthquake. In this paper, we examine ultra-low frequency (ULF) electric field data from the DEMETER satellite during the period leading up to the Sichuan earthquake. It is shown that there is an increase in ULF wave activity observed as DEMETER passes in the vicinity of the earthquake epicentre. This increase is most obvious at lower frequencies. Examination of the ULF spectra shows the possible occurrence of geomagnetic pearl pulsations, resulting from the passage of atmospheric gravity waves generated in the vicinity of the earthquake epicentre.

Description: Determination of wave vectors using the phase differencing method Annales Geophysicae, 31, 1611-1617, 2013 Author(s): S. N. Walker and I. Moiseenko Due to the collisionless nature of space plasmas, plasma waves play an important role in the redistribution of energy between the various particle populations in many regions of geospace. In order to fully comprehend such mechanisms it is necessary to characterise the nature of the waves present. This involves the determination of properties such as wave vector k . There are a number of methods used to calculate k based on the multipoint measurements that are now available. These methods rely on the fact that the same wave packet is simultaneously observed at two or more locations whose separation is small in comparison to the correlation length of the wave packet. This limitation restricts the analysis to low frequency (MHD) waves. In this paper we propose an extension to the phase differencing method to enable the correlation of measurements that were not made simultaneously but differ temporally by a number of wave periods. The method is illustrated using measurements of magnetosonic waves from the Cluster STAFF search coil magnetometer. It is shown that it is possible to identify wave packets whose coherence length is much less than the separation between the measurement locations. The resulting dispersion is found to agree with theoretical results.

Description: Effective CO 2 lifetime and future CO 2 levels based on fit function Annales Geophysicae, 31, 1591-1596, 2013 Author(s): G. R. Sonnemann and M. Grygalashvyly The estimated global CO 2 emission rates and the measured atmospheric CO 2 concentrations show that only a certain share of the emitted CO 2 accumulates in the atmosphere. For given atmospheric emissions of CO 2 , the effective lifetime determines its accumulation in the atmosphere and, consequently, its impact on the future global warming. We found that on average the inferred effective lifetime of CO 2 decreases as its atmospheric concentration increases, reducing the rate of its accumulation in the atmosphere. We derived a power function that fits the varying lifetimes. Based on this fitting function, we calculated the increase of CO 2 for different scenarios of future global emission rates.

Description: Enhanced EISCAT UHF backscatter during high-energy auroral electron precipitation Annales Geophysicae, 31, 1681-1687, 2013 Author(s): N. M. Schlatter, N. Ivchenko, T. Sergienko, B. Gustavsson, and B. U. E. Brändström Natural enhancements in the backscattered power of incoherent scatter radars up to 5 orders of magnitudes above the thermal backscatter are sometimes observed at high latitudes. Recently observations of enhancements in the backscattered power including a feature at zero Doppler shift have been reported. These enhancements are limited in altitude to tens of kilometers. The zero Doppler shift feature has been interpreted as a signature of electron density cavitation. Enhanced plasma lines during these observations have also been reported. We report on the first EISCAT UHF observations of enhanced backscattered radar power including a zero Doppler shift feature. The enhancements originated from two distinct and intermittent layers at about 200 km altitude. The altitude extent of the enhancements, observed during auroral high-energy electron precipitation, was 〈 2 km.

Description: Seasonal dependence of the longitudinal variations of nighttime ionospheric electron density and equivalent winds at southern midlatitudes Annales Geophysicae, 31, 1699-1708, 2013 Author(s): X. Luan and X. Dou It has been indicated that the observed Weddell Sea anomaly (WSA) appeared to be an extreme manifestation of the longitudinal variations in the Southern Hemisphere, since the WSA is characterized by greater evening electron density than the daytime density in the region near the Weddell Sea. In the present study, the longitudinal variations of the nighttime F2-layer peak electron density at southern midlatitudes are analyzed using the observations of the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites between 2006 and 2008. It is found that significant longitudinal difference (〉 150%) relative to the minimum density at each local time prevails in all seasons, although the WSA phenomenon is only evident in summer under this solar minimum condition. Another interesting feature is that in summer, the maximum longitudinal differences occur around midnight (~ 23:00–00:00 LT) rather than in the evening (19:00–21:00 LT) in the evening, when the most prominent electron density enhancement occurs for the WSA phenomenon. Thus the seasonal–local time patterns of the electron density longitudinal variations during nighttime at southern midlatitudes cannot be simply explained in terms of the WSA. Meanwhile, the variations of the geomagnetic configuration and the equivalent magnetic meridional winds/upward plasma drifts are analyzed to explore their contributions to the longitudinal variations of the nighttime electron density. The maximum longitudinal differences are associated with the strongest wind-induced vertical plasma drifts after 21:00 LT in the Western Hemisphere. Besides the magnetic declination–zonal wind effects, the geographic meridional winds and the magnetic inclination also have significant effects on the upward plasma drifts and the resultant electron density.

Description: Extracting gravity wave parameters during the September 2002 Southern Hemisphere major sudden stratospheric warming using a SANAE imaging riometer Annales Geophysicae, 31, 1709-1719, 2013 Author(s): N. Mbatha, V. Sivakumar, H. Bencherif, and S. Malinga Using absorption data measured by imaging riometer for ionospheric studies (IRIS) located at the South Africa National Antarctic Expedition (SANAE), Antarctica (72° S, 3° W), we extracted the parameters of gravity waves (GW) of periods between 40 and 50 min during late winter/spring of the year 2002, a period of the unprecedented major sudden stratospheric warming (SSW) in the Southern Hemisphere middle atmosphere. During this period, an unprecedented substantial increase of temperature by about 25–30 K throughout the stratosphere was observed. During the period of the occurrence of the major stratospheric warming, there was a reduction of both the GW horizontal phase speeds and the horizontal wavelengths at 90 km. The GW phase speeds and horizontal wavelengths were observed to reach minimum values of about 7 m s −1 and 19 km, respectively, while during the quiet period the average value of the phase speed and horizontal wavelength was approximately 23 m s −1 and 62 km, respectively. The observed event is discussed in terms of momentum flux and also a potential interaction of gravity waves, planetary waves and mean circulation.

Description: GNSS troposphere tomography based on two-step reconstructions using GPS observations and COSMIC profiles Annales Geophysicae, 31, 1805-1815, 2013 Author(s): P. Xia, C. Cai, and Z. Liu Traditionally, balloon-based radiosonde soundings are used to study the spatial distribution of atmospheric water vapour. However, this approach cannot be frequently employed due to its high cost. In contrast, GPS tomography technique can obtain water vapour in a high temporal resolution. In the tomography technique, an iterative or non-iterative reconstruction algorithm is usually utilised to overcome rank deficiency of observation equations for water vapour inversion. However, the single iterative or non-iterative reconstruction algorithm has their limitations. For instance, the iterative reconstruction algorithm requires accurate initial values of water vapour while the non-iterative reconstruction algorithm needs proper constraint conditions. To overcome these drawbacks, we present a combined iterative and non-iterative reconstruction approach for the three-dimensional (3-D) water vapour inversion using GPS observations and COSMIC profiles. In this approach, the non-iterative reconstruction algorithm is first used to estimate water vapour density based on a priori water vapour information derived from COSMIC radio occultation data. The estimates are then employed as initial values in the iterative reconstruction algorithm. The largest advantage of this approach is that precise initial values of water vapour density that are essential in the iterative reconstruction algorithm can be obtained. This combined reconstruction algorithm (CRA) is evaluated using 10-day GPS observations in Hong Kong and COSMIC profiles. The test results indicate that the water vapor accuracy from CRA is 16 and 14% higher than that of iterative and non-iterative reconstruction approaches, respectively. In addition, the tomography results obtained from the CRA are further validated using radiosonde data. Results indicate that water vapour densities derived from the CRA agree with radiosonde results very well at altitudes above 2.5 km. The average RMS value of their differences above 2.5 km is 0.44 g m −3 .

Description: Corrigendum to "Estimates of aerosol absorption over India using multi-satellite retrieval" published in Ann. Geophys., 31, 1773–1778, 2013 Annales Geophysicae, 31, 1791-1791, 2013 Author(s): D. Narasimhan and S. K. Satheesh No abstract available.

Description: Frequency variations of gravity waves interacting with a time-varying tide Annales Geophysicae, 31, 1731-1743, 2013 Author(s): C. M. Huang, S. D. Zhang, F. Yi, K. M. Huang, Y. H. Zhang, Q. Gan, and Y. Gong Using a nonlinear, 2-D time-dependent numerical model, we simulate the propagation of gravity waves (GWs) in a time-varying tide. Our simulations show that when a GW packet propagates in a time-varying tidal-wind environment, not only its intrinsic frequency but also its ground-based frequency would change significantly. The tidal horizontal-wind acceleration dominates the GW frequency variation. Positive (negative) accelerations induce frequency increases (decreases) with time. More interestingly, tidal-wind acceleration near the critical layers always causes the GW frequency to increase, which may partially explain the observations that high-frequency GW components are more dominant in the middle and upper atmosphere than in the lower atmosphere. The combination of the increased ground-based frequency of propagating GWs in a time-varying tidal-wind field and the transient nature of the critical layer induced by a time-varying tidal zonal wind creates favorable conditions for GWs to penetrate their originally expected critical layers. Consequently, GWs have an impact on the background atmosphere at much higher altitudes than expected, which indicates that the dynamical effects of tidal–GW interactions are more complicated than usually taken into account by GW parameterizations in global models.

Description: A ten-year record of supraglacial lake evolution and rapid drainage in West Greenland using an automated processing algorithm for multispectral imagery The Cryosphere, 7, 1869-1877, 2013 Author(s): B. F. Morriss, R. L. Hawley, J. W. Chipman, L. C. Andrews, G. A. Catania, M. J. Hoffman, M. P. Lüthi, and T. A. Neumann The rapid drainage of supraglacial lakes introduces large pulses of meltwater to the subglacial environment and creates moulins, surface-to-bed conduits for future melt. Introduction of water to the subglacial system has been shown to affect ice flow, and modeling suggests that variability in water supply and delivery to the subsurface play an important role in the development of the subglacial hydrologic system and its ability to enhance or mitigate ice flow. We developed a fully automated method for tracking meltwater and rapid drainages in large (〉 0.125 km 2 ) perennial lakes and applied it to a 10 yr time series of ETM+ and MODIS imagery of an outlet glacier flow band in West Greenland. Results indicate interannual variability in maximum coverage and spatial evolution of total lake area. We identify 238 rapid drainage events, occurring most often at low ( 〈 900 m) and middle (900–1200 m) elevations during periods of net filling or peak lake coverage. We observe a general progression of both lake filling and draining from lower to higher elevations but note that the timing of filling onset, peak coverage, and dissipation are also variable. Lake coverage is sensitive to air temperature, and warm years exhibit greater variability in both coverage evolution and rapid drainage. Mid-elevation drainages in 2011 coincide with large surface velocity increases at nearby GPS sites, though the relationships between ice-shed-scale dynamics and meltwater input are still unclear.

Description: Diurnal variation in gravity wave activity at low and middle latitudes Annales Geophysicae, 31, 2123-2135, 2013 Author(s): V. F. Andrioli, D. C. Fritts, P. P. Batista, B. R. Clemesha, and D. Janches We employ a modified composite day extension of the Hocking (2005) analysis method to study gravity wave (GW) activity in the mesosphere and lower thermosphere using 4 meteor radars spanning latitudes from 7° S to 53.6° S. Diurnal and semidiurnal modulations were observed in GW variances over all sites. Semidiurnal modulation with downward phase propagation was observed at lower latitudes mainly near the equinoxes. Diurnal modulations occur mainly near solstice and, except for the zonal component at Cariri (7° S), do not exhibit downward phase propagation. At a higher latitude (SAAMER, 53.6° S) these modulations are only observed in the meridional component where we can observe diurnal variation from March to May, and semidiurnal, during January, February, October (above 88 km) and November. Some of these modulations with downward phase progression correlate well with wind shear. When the wind shear is well correlated with the maximum of the variances the diurnal tide has its largest amplitudes, i.e., near equinox. Correlations exhibiting variations with tidal phases suggest significant GW-tidal interactions that have different characters depending on the tidal components and possible mean wind shears. Modulations that do not exhibit phase variations could be indicative of diurnal variations in GW sources.

Description: Influence of grain shape on light penetration in snow The Cryosphere, 7, 1803-1818, 2013 Author(s): Q. Libois, G. Picard, J. L. France, L. Arnaud, M. Dumont, C. M. Carmagnola, and M. D. King The energy budget and the photochemistry of a snowpack depend greatly on the penetration of solar radiation in snow. Below the snow surface, spectral irradiance decreases exponentially with depth with a decay constant called the asymptotic flux extinction coefficient. As with the albedo of the snowpack, the asymptotic flux extinction coefficient depends on snow grain shape. While representing snow by a collection of spherical particles has been successful in the numerical computation of albedo, such a description poorly explains the decrease of irradiance in snow with depth. Here we explore the limits of the spherical representation. Under the assumption of geometric optics and weak absorption by snow, the grain shape can be simply described by two parameters: the absorption enhancement parameter B and the geometric asymmetry factor g G . Theoretical calculations show that the albedo depends on the ratio B /(1- g G ) and the asymptotic flux extinction coefficient depends on the product B (1- g G ). To understand the influence of grain shape, the values of B and g G are calculated for a variety of simple geometric shapes using ray tracing simulations. The results show that B and (1- g G ) generally covary so that the asymptotic flux extinction coefficient exhibits larger sensitivity to the grain shape than albedo. In particular it is found that spherical grains propagate light deeper than any other investigated shape. In a second step, we developed a method to estimate B from optical measurements in snow. A multi-layer, two-stream, radiative transfer model, with explicit grain shape dependence, is used to retrieve values of the B parameter of snow by comparing the model to joint measurements of reflectance and irradiance profiles. Such measurements were performed in Antarctica and in the Alps yielding estimates of B between 0.8 and 2.0. In addition, values of B were estimated from various measurements found in the literature, leading to a wider range of values (1.0–9.9) which may be partially explained by the limited accuracy of the data. This work highlights the large variety of snow microstructure and experimentally demonstrates that spherical grains, with B = 1.25, are inappropriate to model irradiance profiles in snow, an important result that should be considered in further studies dedicated to subsurface absorption of short-wave radiation and snow photochemistry.

Description: Influence of ice-sheet geometry and supraglacial lakes on seasonal ice-flow variability The Cryosphere, 7, 1185-1192, 2013 Author(s): I. Joughin, S. B. Das, G. E. Flowers, M. D. Behn, R. B. Alley, M. A. King, B. E. Smith, J. L. Bamber, M. R. van den Broeke, and J. H. van Angelen Supraglacial lakes play an important role in establishing hydrological connections that allow lubricating seasonal meltwater to reach the base of the Greenland Ice Sheet. Here we use new surface velocity observations to examine the influence of supraglacial lake drainages and surface melt rate on ice flow. We find large, spatially extensive speedups concurrent with times of lake drainage, showing that lakes play a key role in modulating regional ice flow. While surface meltwater is supplied to the bed via a geographically sparse network of moulins, the observed ice-flow enhancement suggests that this meltwater spreads widely over the ice-sheet bed. We also find that the complex spatial pattern of speedup is strongly determined by the combined influence of bed and surface topography on subglacial water flow. Thus, modeling of ice-sheet basal hydrology likely will require knowledge of bed topography resolved at scales (sub-kilometer) far finer than existing data (several km).

Description: The effects of additional black carbon on the albedo of Arctic sea ice: variation with sea ice type and snow cover The Cryosphere, 7, 1193-1204, 2013 Author(s): A. A. Marks and M. D. King The response of the albedo of bare sea ice and snow-covered sea ice to the addition of black carbon is calculated. Visible light absorption and light-scattering cross-sections are derived for a typical first-year and multi-year sea ice with both "dry" and "wet" snow types. The cross-sections are derived using data from a 1970s field study that recorded both reflectivity and light penetration in Arctic sea ice and snow overlying sea ice. The variation of absorption cross-section over the visible wavelengths suggests black carbon is the dominating light-absorbing impurity. The response of first-year and multi-year sea ice albedo to increasing black carbon, from 1 to 1024 ng g −1 , in a top 5 cm layer of a 155 cm-thick sea ice was calculated using a radiative-transfer model. The albedo of the first-year sea ice is more sensitive to additional loadings of black carbon than the multi-year sea ice. An addition of 8 ng g −1 of black carbon causes a decrease to 98.7% of the original albedo for first-year sea ice compared to a decrease to 99.7% for the albedo of multi-year sea ice, at a wavelength of 500 nm. The albedo of sea ice is surprisingly unresponsive to additional black carbon up to 100 ng g −1 . Snow layers on sea ice may mitigate the effects of black carbon in sea ice. Wet and dry snow layers of 0.5, 1, 2, 5 and 10 cm depth were added onto the sea ice surface. The albedo of the snow surface was calculated whilst the black carbon in the underlying sea ice was increased. A layer of snow 0.5 cm thick greatly diminishes the effect of black carbon in sea ice on the surface albedo. The albedo of a 2–5 cm snow layer (less than the e -folding depth of snow) is still influenced by the underlying sea ice, but the effect of additional black carbon in the sea ice is masked.

Description: A study of solar and interplanetary parameters of CMEs causing major geomagnetic storms during SC 23 Annales Geophysicae, 31, 1285-1295, 2013 Author(s): C. Oprea, M. Mierla, D. Beşliu-Ionescu, O. Stere, and G. Mariş Muntean In this paper we analyse 25 Earth-directed and strongly geoeffective interplanetary coronal mass ejections (ICMEs) which occurred during solar cycle 23, using data provided by instruments on SOHO (Solar and Heliospheric Observatory), ACE (Advanced Composition Explorer) and geomagnetic stations. We also examine the in situ parameters, the energy transfer into magnetosphere, and the geomagnetic indexes. We compare observed travel times with those calculated by observed speeds projected into the plane of the sky and de-projected by a simple model. The best fit was found with the projected speeds. No correlation was found between the importance of a flare and the geomagnetic Dst (disturbance storm time) index. By comparing the in situ parameters with the Dst index we find a strong connection between some of these parameters (such as Bz, Bs · V and the energy transfer into the magnetosphere) with the strength of the geomagnetic storm. No correlation was found with proton density and plasma temperature. A superposed epoch analysis revealed a strong dependence of the Dst index on the southward component of interplanetary magnetic field, Bz, and to the Akasofu coupling function, which evaluates the energy transfer between the ICME and the magnetosphere. The analysis also showed that the geomagnetic field at higher latitudes is disturbed before the field around the Earth's equator.