ALBERT

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

Description: Nature Geoscience 6, 828 (2013). doi:10.1038/ngeo1931 Authors: Shang-Ping Xie, Bo Lu & Baoqiang Xiang Spatial variations in ocean warming have been linked to regional changes in tropical cyclones, precipitation and monsoons. But development of reliable regional climate projections for climate change mitigation and adaptation remains challenging. The presence of anthropogenic aerosols, which are highly variable in space and time, is thought to induce spatial patterns of climate response that are distinct from those of well-mixed greenhouse gases. Using CMIP5 climate simulations that consider aerosols and greenhouse gases separately, we show that regional responses to changes in greenhouse gases and aerosols are similar over the ocean, as reflected in similar spatial patterns of ocean temperature and precipitation. This similarity suggests that the climate response to radiative changes is relatively insensitive to the spatial distribution of these changes. Although anthropogenic aerosols are largely confined to the Northern Hemisphere, simulations that include aerosol forcing predict decreases in temperature and westerly wind speed that reach the pristine Southern Hemisphere oceans. Over land, the climate response to aerosol forcing is more localized, but larger scale spatial patterns are also evident. We suggest that the climate responses induced by greenhouse gases and aerosols share key ocean–atmosphere feedbacks, leading to a qualitative resemblance in spatial distribution.

Description: Nature Geoscience 6, 856 (2013). doi:10.1038/ngeo1928 Authors: Chloé Michaut, Yanick Ricard, David Bercovici & R. Steve J. Sparks Eruptions at active silicic volcanoes are often cyclical. For example, at the Soufrière Hills volcano in Montserrat, Mount Pinatubo in the Philippines, and Sakurajima in Japan, episodes of intense activity alternate with repose intervals over periods between several hours and a day. Abrupt changes in eruption rates have been explained with the motion of a plug of magma that alternatively sticks or slides along the wall of the volcanic conduit. However, it is unclear how the static friction that prevents the plug from sliding is periodically overcome. Here we use two-phase flow equations to model a gas-rich, viscous magma ascending through a volcanic conduit. Our analyses indicate that magma compaction yields ascending waves comprised of low- and high-porosity bands. However, magma ascent to lower pressures also causes gas expansion. We find that the competition between magma compaction and gas expansion naturally selects pressurized gas waves with specific periods. At the surface, these waves can induce cyclical eruptive behaviour with periods between 1 and 100 hours, which compares well to the observations from Soufrière Hills, Mount Pinatubo and Sakurajima. We find that the period is insensitive to volcano structure, but increases weakly with magma viscosity. We conclude that observations of a shift to a longer eruption cycle imply an increase in magma viscosity and thereby enhanced volcanic hazard.

Description: Nature Geoscience 6, 866 (2013). doi:10.1038/ngeo1933 Authors: Suzanne M. Carbotte, Milena Marjanović, Helene Carton, John C. Mutter, Juan Pablo Canales, Mladen R. Nedimović, Shuoshuo Han & Michael R. Perfit The global mid-ocean ridge is segmented in its seafloor morphology and magmatic systems, but the origin of and relationships between this tectonic and magmatic segmentation are poorly understood. At fast-spreading ridges, tectonic segmentation is observed on a fine scale, but it is unclear whether this partitioning also occurs in the magmatic system. Fine-scale tectonic segmentation could have a deep origin, arising from the distribution of upwelling mantle melt, or a shallow origin, linked to offset intruding dikes from long, more continuous crustal reservoirs. Here we use seismic reflection data from the fast-spreading East Pacific Rise, between 8° 20′ N and 10° 10′ N, which includes a unique area where two documented volcanic eruptions have occurred, to image the crustal magma bodies in high resolution. We find that the magma reservoirs form 5- to 15-km-long segments that coincide with the fine-scale tectonic segmentation at the seafloor and that three lens segments fed the recent eruptions. Transitions in composition, volume and morphology of erupted lavas coincide with disruptions in the lens that define magmatic segments. We conclude that eruptions at the East Pacific Rise are associated with the vertical ascent of magma from lenses that are mostly physically isolated, leading to the eruption of distinct lavas at the surface that coincide with fine-scale tectonic segmentation.

Description: Nature Geoscience 6, 871 (2013). doi:10.1038/ngeo1911 Authors: Judith A. Coggon, Ambre Luguet, Geoffrey M. Nowell & Peter W. U. Appel Partial melting of the Earth’s mantle is a key process in the generation of crustal material and the formation of continents. Crustal samples record the generation of crust up to 4.4 billion years (Gyr) ago, yet the complementary record in the mantle extends to only 3.3 Gyr ago, with sparse evidence for differentiation occurring 3.9–4.1 Gyr ago. Here we use the Pt–Os isotope chronometer to show that a Hadean record of mantle depletion is preserved in Earth’s oldest known ultramafic rocks, the Ujaragssuit Nunât intrusion of southwest Greenland. We identify two distinct age populations at approximately 4.1 and 2.9 Gyr. We suggest that the younger age population records a regional metamorphic event and the older one records mantle depletion. We also identify individual sample ages of up to 4.36 Gyr old, thus extending the record of large mantle-melting events into the Hadean. Furthermore, the preservation of Hadean model ages in Os-rich mantle-derived rocks supports the theory that re-enrichment of Os in the mantle during the Late Heavy Bombardment—after expected partitioning into the Earth’s core—occurred at least 0.2 Gyr earlier than previously thought. This also implies that the Earth could have been habitable by 4.1 Gyr ago.

Description: Nature Geoscience 6, 875 (2013). doi:10.1038/ngeo1903 Authors: Jun Wu & Peter R. Buseck A significant fraction of Earth’s carbon resides in the mantle, but the mode of carbon storage presents a long-standing problem. The mantle contains fluids rich in carbon dioxide and methane, carbonate-bearing melts, carbonate minerals, graphite, diamond and carbides, as well as dissolved carbon atoms in metals. However, it is uncertain whether these can sufficiently account for the total amount of carbon thought to be stored in the mantle and the volume of carbon degassed from the mantle at volcanoes. Moreover, such carbon hosts should significantly affect the physical and chemical behaviour of the mantle, including its melting temperature, electrical conductivity and oxidation state. Here we use in situ transmission electron microscopy to measure the storage of carbon within common mantle mineral analogues—nickel-doped lanthanum chromate perovskite and titanium dioxide—in laboratory experiments at high pressure and temperature. We detect elevated carbon concentrations at defect sites in the nanocrystals, maintained at high pressures within annealed carbon nanocages. Specifically, our experiments show that small stacking faults within the mantle analogue materials are effective carbon sinks at mantle conditions, potentially providing an efficient mechanism for carbon storage in the mantle. Furthermore, this carbon can be readily released under lower pressure conditions, and may therefore help to explain carbon release in volcanic eruptions.

Description: Nature Geoscience 6, 885 (2013). doi:10.1038/ngeo1922 Authors: Lars Möller, Todd Sowers, Michael Bock, Renato Spahni, Melanie Behrens, Jochen Schmitt, Heinrich Miller & Hubertus Fischer

Description: Nature Geoscience 6, 847 (2013). doi:10.1038/ngeo1905 Authors: Marcello Campione & Gian Carlo Capitani Earthquakes generated in subduction zones are caused by unstable movements along faults. This fault-slip instability is determined by frictional forces that depend on the temperature, pressure, morphology and deformation state of the fault rocks. Fault friction may also be influenced by preferred mineral orientations. Over-thrusting of rocks at the interface between a subducting slab and the overlying mantle wedge generates shear deformation that causes minerals to align, and this preferred mineral orientation affects the propagation of shear seismic waves. Here we use laboratory experiments to simulate fault slip in antigorite, the most abundant hydrous mineral phase within Earth’s upper mantle. Using atomic force microscopy, we show that antigorite single crystals possess strong frictional anisotropy on their basal slip surface and that preferred mineral alignment extends this property to a regional scale. Depending on the alignment, fault movements can occur along a high-friction direction, creating stick-slip behaviour that generates earthquakes. In contrast, if movements occur along a low-friction direction, the mantle wedge will deform aseismically. Our results imply that mantle rocks in subduction-zone thrust faults can exhibit two opposite frictional behaviours, seismic and aseismic.

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

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

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

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

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

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

Description: Nature Geoscience 6, 801 (2013). doi:10.1038/ngeo1983 The latest report on global warming brings yet another rise in confidence that human actions are altering the Earth's climate. But in contrast to its 2007 predecessor, it is unlikely to cause a stir.

Description: Nature Geoscience 6, 802 (2013). doi:10.1038/ngeo1971 Authors: Thomas Shea, Julia Hammer & Emily First

Description: Nature Geoscience 6, 803 (2013). doi:10.1038/ngeo1970 Authors: Melissa D. Rotella, Colin J. N. Wilson, Simon J. Barker & Ian C. Wright

Description: Nature Geoscience 6, 805 (2013). doi:10.1038/ngeo1981 Author: Amy Whitchurch

Description: Nature Geoscience 6, 807 (2013). doi:10.1038/ngeo1960 Author: Julie Brigham-Grette During the Last Glacial Maximum, ice sheets in Eurasia terminated at the edge of the Laptev Sea. Seismic data now suggest that a separate ice sheet was repeatedly centred further east, in the East Siberian Sea, during previous glacial periods.

Description: Nature Geoscience 6, 805 (2013). doi:10.1038/ngeo1979 Author: Tamara Goldin

Description: Nature Geoscience 6, 811 (2013). doi:10.1038/ngeo1951 Author: W. Roger Buck Intrusions of magma into the crust help accommodate the divergence between tectonic plates. A magnetotelluric survey of the crust and mantle beneath Afar, Ethiopia, has identified enough magma to accommodate plate separation there for about 50,000 years.

Description: Nature Geoscience 6, 805 (2013). doi:10.1038/ngeo1978 Author: Anna Armstrong

Description: Nature Geoscience 6, 813 (2013). doi:10.1038/ngeo1955 Authors: Stefanie Kirschke, Philippe Bousquet, Philippe Ciais, Marielle Saunois, Josep G. Canadell, Edward J. Dlugokencky, Peter Bergamaschi, Daniel Bergmann, Donald R. Blake, Lori Bruhwiler, Philip Cameron-Smith, Simona Castaldi, Frédéric Chevallier, Liang Feng, Annemarie Fraser, Martin Heimann, Elke L. Hodson, Sander Houweling, Béatrice Josse, Paul J. Fraser, Paul B. Krummel, Jean-François Lamarque, Ray L. Langenfelds, Corinne Le Quéré, Vaishali Naik, Simon O'Doherty, Paul I. Palmer, Isabelle Pison, David Plummer, Benjamin Poulter, Ronald G. Prinn, Matt Rigby, Bruno Ringeval, Monia Santini, Martina Schmidt, Drew T. Shindell, Isobel J. Simpson, Renato Spahni, L. Paul Steele, Sarah A. Strode, Kengo Sudo, Sophie Szopa, Guido R. van der Werf, Apostolos Voulgarakis, Michiel van Weele, Ray F. Weiss, Jason E. Williams & Guang Zeng

Description: Nature Geoscience 6, 805 (2013). doi:10.1038/ngeo1980 Author: Alicia Newton

Description: Nature Geoscience 6, 809 (2013). doi:10.1038/ngeo1953 Author: Klaus Hasselmann

Description: Nature Geoscience 6, 804 (2013). doi:10.1038/ngeo1965 Author: Emily Lakdawalla

Description: Nature Geoscience 6, 833 (2013). doi:10.1038/ngeo1887 Authors: J. N. Bassis & S. Jacobs Iceberg calving has been implicated in the retreat and acceleration of glaciers and ice shelves along the margins of the Greenland and Antarctic ice sheets. Accurate projections of sea-level rise therefore require an understanding of how and why calving occurs. Unfortunately, calving is a complex process and previous models of the phenomenon have not reproduced the diverse patterns of iceberg calving observed in nature. Here we present a numerical model that simulates the disparate calving regimes observed, including the detachment of large tabular bergs from floating ice tongues, the disintegration of ice shelves and the capsizing of smaller bergs from grounded glaciers that terminate in deep water. Our model treats glacier ice as a granular material made of interacting boulders of ice that are bonded together. Simulations suggest that different calving regimes are controlled by glacier geometry, which controls the stress state within the glacier. We also find that calving is a two-stage process that requires both ice fracture and transport of detached icebergs away from the calving front. We suggest that, as a result, rapid iceberg discharge is possible in regions where highly crevassed glaciers are grounded deep beneath sea level, indicating portions of Greenland and Antarctica that may be vulnerable to rapid ice loss through catastrophic disintegration.

Description: Nature Geoscience 6, 824 (2013). doi:10.1038/ngeo1923 Authors: C. T. Adcock, E. M. Hausrath & P. M. Forster If the chemistry essential to life was present in water-containing environments on Mars, the processes that led to life on Earth may have also occurred on the red planet. Phosphate is one of the chemical nutrients thought to be essential for life and is also considered critical to reactions that may have led to life on Earth. However, low prebiotic availability of phosphate may have been a complicating factor in terrestrial abiogenesis, suggesting that a similar hurdle may have confronted the development of life on Mars. Phosphate available for biological reactions can be introduced into aqueous environments through dissolution of primary phosphate minerals during water–rock interactions, but little is known about the dissolution of the dominant phosphate minerals found in martian meteorites and presumably on Mars. Here we present dissolution rates, phosphate release rates and solubilities of phosphate minerals found in martian rocks as determined from laboratory measurements. Our experimental findings predict phosphate release rates during water–rock interactions on Mars that are as much as 45 times higher than on Earth and phosphate concentrations of early wet martian environments more than twice those of Earth. We suggest that available phosphate may have mitigated one of the hurdles to abiogenesis on Mars.

Description: Nature Geoscience 6, 837 (2013). doi:10.1038/ngeo1924 Authors: María Fernanda Sánchez Goñi, Edouard Bard, Amaelle Landais, Linda Rossignol & Francesco d’Errico A period of continental ice growth between about 80,000 and 70,000 years ago was controlled by a decrease in summer insolation, and was among the four largest ice expansions of the past 250,000 years. The moisture source for this ice sheet expansion, known as the Marine Isotope Stage (MIS) 5a/4 transition, has been proposed to be the warm subpolar and northern subtropical Atlantic Ocean. However, the mechanism by which glaciers kept growing through three suborbital cooling events within this period, which were associated with iceberg discharge in the North Atlantic and cooling over Greenland, is unclear. Here we reconstruct parallel records of sea surface and air temperatures from marine microfossil and pollen data, respectively, from two sediment cores collected within the northern subtropical gyre. The thermal gradient between the cold air and warmer sea increased throughout the MIS5a/4 transition, and was marked by three intervals of even more pronounced thermal gradients associated with the C20, C19 and C18’ cold events. We argue that the warm ocean surface along the western European margin provided a source of moisture that was transported, through northward-tracking storms, to feed ice sheets in colder Greenland, northern Europe and the Arctic.

Description: Nature Geoscience 6, 842 (2013). doi:10.1038/ngeo1904 Authors: Frank Niessen, Jong Kuk Hong, Anne Hegewald, Jens Matthiessen, Rüdiger Stein, Hyoungjun Kim, Sookwan Kim, Laura Jensen, Wilfried Jokat, Seung-Il Nam & Sung-Ho Kang During the Pleistocene glaciations, Arctic ice sheets on western Eurasia, Greenland and North America terminated at their continental margins. In contrast, the exposed continental shelves in the Beringian region of Siberia are thought to have been covered by a tundra landscape. Evidence of grounded ice on seafloor ridges and plateaux off the coast of the Beringian margin, at depths of up to 1,000 m, have generally been attributed to ice shelves or giant icebergs that spread oceanwards during glacial maxima. Here we identify marine glaciogenic landforms visible in seismic profiles and detailed bathymetric maps along the East Siberian continental margin. We interpret these features, which occur in present water depths of up to 1,200 m, as traces from grounding events of ice sheets and ice shelves. We conclude that the Siberian Shelf edge and parts of the Arctic Ocean were covered by ice sheets of about 1 km in thickness during several Pleistocene glaciations before the most recent glacial period, which must have had a significant influence on albedo and oceanic and atmospheric circulation.

Description: Nature Geoscience 6, 861 (2013). doi:10.1038/ngeo1925 Authors: M. Desissa, N. E. Johnson, K. A. Whaler, S. Hautot, S. Fisseha & G. J. K. Dawes Shallow magma reservoirs exist in the crust beneath volcanoes and mid-ocean ridges, yet there are no reports of extensive magma bodies within the uppermost mantle. Indeed the buoyancy of magma should cause it to intrude into the crust, preventing it from ponding in the mantle below. The Dabbahu magmatic segment in Afar, Ethiopia, marks the late stages of continental rifting. This segment has been active since 2005 and has experienced repeated magma intrusions. Here we use magnetotelluric data to image magma bodies beneath it. We identify a 30-km-wide region of very high electrical conductivity that reaches down to about 35 km depth. We interpret this region as a large volume of magma of at least 500 km3 that extends well into the mantle and contains about 13% melt fraction. The magma volume is orders of magnitude larger than that intruded during a typical rifting episode, implying that the magma reservoir persists for several tens of thousands of years. This is in marked contrast to the situation beneath mid-ocean ridges, where melt supply is thought to be episodic. Large magma reservoirs within the mantle may therefore be responsible for the localization of strain that accompanies the final stages of continental break-up.

Description: Nature Geoscience 6, 852 (2013). doi:10.1038/ngeo1927 Authors: Pascal Audet & Susan Y. Schwartz Subduction zones can exhibit variable seismic behaviour, ranging from great earthquakes to slow slip. This variability may be linked to fault frictional properties, and the rheology and structure of the upper plate. The subduction zone beneath the Nicoya Peninsula, Costa Rica, is characterized by strong variations in fault-slip behaviour and a lateral change in the origin of the subducting plate. In the northwest, the plate interface is locked, and experiences large, infrequent earthquakes, and the subducting plate is formed at the East Pacific Rise. In contrast, in the southeast, slow-slip events occur frequently and the subducting plate is formed at the Cocos–Nazca spreading centre. Here we use seismic receiver-function data to analyse the structure of the subduction zone beneath the Nicoya Peninsula. We find extremely high P–S seismic-velocity ratios within the entire subducting oceanic crust that we interpret as high pore-fluid pressure. Velocity ratios in the overriding continental crust, however, change from lower values in the northwest to higher ones in the southeast, indicating a disparity in fluid accumulation. We infer that this disparity is caused by a higher supply of fluid from the subducting slab in the southeast, owing to the permeability structure of oceanic crust formed at the Cocos–Nazca spreading centre. We suggest that the spatial gradient in fluid content influences upper-plate strength and controls the segmentation of seismogenic behaviour in this subduction zone.

Description: Nature Geoscience 6, 806 (2013). doi:10.1038/ngeo1929 Author: Matthew Pasek Phosphorus is an important element for biogeochemical development. According to a set of experiments, martian phosphate minerals dissolve more quickly than terrestrial ones, possibly providing nutrients in aqueous environments for early martian life.

Description: Nature Geoscience 6, 810 (2013). doi:10.1038/ngeo1967 Author: Daniel Cossa The neurotoxin methylmercury can accumulate in marine food webs, contaminating seafood. An analysis of the isotopic composition of fish in the North Pacific suggests that much of the mercury that enters the marine food web originates from low-oxygen subsurface waters.

Description: Nature Geoscience 6, 879 (2013). doi:10.1038/ngeo1918 Authors: Joel D. Blum, Brian N. Popp, Jeffrey C. Drazen, C. Anela Choy & Marcus W. Johnson

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: Nature Geoscience 7, 860 (2014). doi:10.1038/ngeo2301 Author: Jay J. Ague Aqueous subduction-zone fluids contain CO2 and methane. New calculations indicate that these fluids also host a wide array of organic carbon species, in concentrations sufficient to influence the deep carbon cycle.

Description: Nature Geoscience 7, 925 (2014). doi:10.1038/ngeo2278 Authors: Devin McPhillips, Paul R. Bierman & Dylan H. Rood

Description: Nature Geoscience 7, 920 (2014). doi:10.1038/ngeo2284 Authors: Peter B. Reich, Sarah E. Hobbie & Tali D. Lee

Description: Nature Geoscience 7, 931 (2014). doi:10.1038/ngeo2296 Author: Colin Goldblatt, Tyler D. Robinson, Kevin J. Zahnle & David Crisp

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

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

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

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

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

Description: Tweets on Earth Nature Geoscience 4, 209 (2011). doi:10.1038/ngeo1137 Nature Geoscience has joined Twitter. We share our take on exciting developments in the Earth and planetary sciences as they happen.

Description: Closed ranks in oceanography Nature Geoscience 4, 211 (2011). doi:10.1038/ngeo1113 Authors: LuAnne Thompson, Renellys C. Perez & Amelia E. Shevenell

Description: Not the earthquake's fault Nature Geoscience 4, 210 (2011). doi:10.1038/ngeo1116 Authors: Brian G. McAdoo & Lisabeth Paravisini-Gebert

Description: Coal's true cost Nature Geoscience 4, 209 (2011). doi:10.1038/ngeo1136 The deaths of birds have become a rallying point against the proliferation of wind farms. Yet the loss of human life in mines is rarely linked with coal as an energy source.

Description: Our choice from the recent literature Nature Geoscience 4, 215 (2011). doi:10.1038/ngeo1131

Description: Science-based insurance Nature Geoscience 4, 213 (2011). doi:10.1038/ngeo1117 Authors: Molly E. Brown, Daniel E. Osgood & Miguel A. Carriquiry Crops are at risk in a changing climate. Farmers in the developing world will be able to insure against harvest failure if robust insurance packages, based on a geophysical index rather than individual loss, become widely available.

Description: Earth's core: A new twist on inner-core spin Nature Geoscience 4, 216 (2011). doi:10.1038/ngeo1091 Author: Mathieu Dumberry An observed hemispheric structure in the Earth's inner core has been hard to reconcile with evidence that it rotates faster than the mantle. Detection of a shift of the hemisphere boundary that occurred over geological timescales removes the contradiction.

Description: Structural geology: Natural fault lubricants Nature Geoscience 4, 217 (2011). doi:10.1038/ngeo1126 Author: Ben van der Pluijm Motion along faults can occur in sudden earthquakes or through steady, aseismic creep. Rock samples retrieved by drilling deep into a creeping section of the San Andreas Fault show that clay minerals in fault rock promote creep behaviour.

Description: Transformation of tectonic and climatic signals from source to sedimentary archive Nature Geoscience 4, 231 (2011). doi:10.1038/ngeo1087 Authors: John J. Armitage, Robert A. Duller, Alex C. Whittaker & Philip A. Allen The Earth’s sedimentary successions are an archive of past tectonic and climate events. The physical characteristics of the sediment record are controlled by three main factors: the sediment supply from the eroding source region, the grain size distribution of that sediment supply, and the area available for sediment accumulation in the downstream regions. The interplay of these factors can make the interpretation of a climatic or tectonic signal complex, particularly as these processes are nonlinear. Here we assess the evolution of a tectonically active landscape undergoing erosional and depositional processes, using a numerical model that incorporates variations in grain size and the volumetric sediment budget. Our simulations indicate that changes in precipitation and tectonic uplift both generate characteristic patterns of grain size and stratigraphic geometry. An increase in catchment precipitation results in the deposition of a laterally extensive sheet of coarse gravel. The responses to a changing tectonic regime are more diverse: a large increase in uplift rate results first in the deposition of sediments of larger grain size at proximal sites, followed by a reduction in grain size at distal locations. We conclude that the stratigraphic record is strongly controlled by the grain size of sediments released from catchments undergoing tectonic or climatic change.

Description: Early Solar System: Gathering dust Nature Geoscience 4, 219 (2011). doi:10.1038/ngeo1112 Author: Jeff Cuzzi Chondritic meteorites are remnants of the ancient Solar System. Analysis of the dust rims often found on their constituent particles shows that the rims were swept up while the particles wafted about and collided in a weakly turbulent protoplanetary nebula.

Description: Oceanography: Chemical twins, separated Nature Geoscience 4, 220 (2011). doi:10.1038/ngeo1125 Author: Martin Frank How the chemical composition of sea water changes on its journey through the world's oceans is poorly understood. Systematic measurements of dissolved trace metals across the Pacific Ocean suggest that these metals may help track sources and mixing of water masses.

Description: Weakening of the equatorial Atlantic cold tongue over the past six decades Nature Geoscience 4, 222 (2011). doi:10.1038/ngeo1078 Authors: Hiroki Tokinaga & Shang-Ping Xie Seasonal and interannual variations of the equatorial cold tongue are defining features of the tropical Atlantic Ocean, with significant climatic and biogeochemical effects. However, its long-term changes are poorly understood owing to biases in observations and climate models. Here we use a suite of bias-corrected observations, and find that cold-tongue variability has weakened during the past six decades. We find that sea surface temperature has increased across the basin, with a local enhancement over the eastern equatorial Atlantic. This warming pattern of the sea surface is most pronounced during boreal summer, reducing the annual cycle through a positive ocean–atmosphere feedback. Specifically, the eastward-intensified warming leads to enhanced atmospheric convection in the equatorial eastern Atlantic region, as well as to less vigorous trade winds. These in turn deepen the thermocline in the east, and reinforce the sea surface warming pattern. The flattened thermocline and reduced thermocline feedback weaken interannual variability of equatorial sea surface temperatures and Guinea coast precipitation associated with the Atlantic Niño. We suggest that the observed changes could be associated with cooling by anthropogenic aerosols, an effect that is stronger in the Northern than in the Southern Hemisphere. If the aerosol emissions decrease in the next decades, the tropical Atlantic may experience yet another shift as the greenhouse gas forcing increases.

Description: Strong elemental fractionation of Zr–Hf and Nb–Ta across the Pacific Ocean Nature Geoscience 4, 227 (2011). doi:10.1038/ngeo1114 Authors: M. Lutfi Firdaus, Tomoharu Minami, Kazuhiro Norisuye & Yoshiki Sohrin Understanding the circulation of water masses in the world’s oceans is critical to our knowledge of the Earth’s climate system. Trace elements and their isotopes have been explored as tracers for the movement of water masses . One type of candidate elements are the high-field-strength elements zirconium (Zr), hafnium (Hf), niobium (Nb) and tantalum (Ta). Here we measure the distributions of dissolved Zr, Hf, Nb and Ta along two meridional sections in the Pacific Ocean that extend from 65° to 10° S and from 10° to 50° N. We find that all four elements tend to be depleted in surface water. In the deep oceans, their concentrations rise along our transects from the Southern Ocean to the North Pacific Ocean, and show strong correlations with the concentration of silicate. These results indicate that terrigenous sources are important to the budget of Zr, Hf, Nb and Ta in sea water, compared with hydrothermal input. Unexpectedly, the weight ratios for Zr/Hf fall between 45 and 350 and those for Nb/Ta between 14 and 85 in Pacific sea water, higher than the ratios observed in fresh water, in the silicate Earth or in chondritic meteorites. We conclude that the fractionation of Zr/Hf and Nb/Ta ratios will be useful for tracing water masses in the ocean.

Description: Earliest rock fabric formed in the Solar System preserved in a chondrule rim Nature Geoscience 4, 244 (2011). doi:10.1038/ngeo1120 Authors: Philip A. Bland, Lauren E. Howard, David J. Prior, John Wheeler, Robert M. Hough & Kathryn A. Dyl Rock fabrics—the preferred orientation of grains—provide a window into the history of rock formation, deformation and compaction. Chondritic meteorites are among the oldest materials in the Solar System and their fabrics should record a range of processes occurring in the nebula and in asteroids. However, owing to abundant fine-grained material, chondrites have largely resisted traditional in situ fabric analysis. Here we use high-resolution electron backscatter diffraction to map the orientation of submicrometre grains in the Allende CV carbonaceous chondrite. We look at the fine-grained rims surrounding the chondrules—spherical grains cooled from molten droplets before accretion in the meteorite—as well as the matrix material between the chondrules. Although the matrix exhibits a bulk uniaxial fabric indicative of a compressive event in the parent asteroid, we find that the chondrule rims preserve a spherically symmetric fabric centred on the chondrule. We define a method to quantitatively relate fabric intensity to net compression, and reconstruct an initial rim porosity of 70–80%. Our calculations provide meteoritic evidence that the first solids formed in the Solar System accreted with high porosity, similar to modelling and laboratory estimates. We conclude that the chondrule rim textures formed in a nebula setting and may therefore represent the first rock fabric in the Solar System.

Description: The protracted development of the continent–ocean transition in Afar Nature Geoscience 4, 248 (2011). doi:10.1038/ngeo1095 Authors: Ian D. Bastow & Derek Keir Continental breakup and the transition to seafloor spreading is characterized by extensional faulting, thinning of the lithosphere and, at magmatic margins, voluminous intrusive and extrusive magmatism. It is difficult to discriminate between different mechanisms of extension and magmatism at ancient continental margins because the continent–ocean transition is buried beneath thick layers of volcanic and sedimentary rocks and the tectonic activity that characterized breakup has ceased. Instead, the timing of these mechanisms is inferred from theoretical models or from the geological record preserved at the fully developed, ancient rifted margins. Ongoing rifting in Ethiopia offers a unique opportunity to address these problems because it exposes subaerially the transition between continental rifting towards the south and seafloor spreading further northward. Here we synthesize constraints on the spatial and temporal evolution of magmatism and extension in Ethiopia. We show that although intrusion of magma maintains crustal thickness during the early stages of the continent–ocean transition, subsidence of the margin below sea level, and eruption of voluminous basalt flows, is initiated by late-stage thinning of the heavily intruded, weakened plate just before the onset of seafloor spreading. We thus conclude that faulting, stretching and magma intrusion are each important, but at different times during breakup.

Description: Explosive eruptions at mid-ocean ridges driven by CO2-rich magmas Nature Geoscience 4, 260 (2011). doi:10.1038/ngeo1104 Authors: Christoph Helo, Marc-Antoine Longpré, Nobumichi Shimizu, David A. Clague & John Stix The abundance of volatile compounds, and particularly CO2, in the upper oceanic mantle affects the style of volcanic eruptions. At mid-ocean ridges, eruptions are generally dominated by the gentle effusion of basaltic lavas with a low volatile content. However explosive volcanism has been documented at some ocean spreading centres, indicative of abundant volatile compounds. Estimates of the initial CO2 concentration of primary magmas can be used to constrain the CO2 content of the upper oceanic mantle, but these estimates vary greatly. Here we present ion microprobe measurements of the CO2 content of basaltic melt trapped in plagioclase crystals. The crystals are derived from volcanic ash deposits erupted explosively at Axial Seamount, Juan de Fuca Ridge, in the northeast Pacific Ocean. We report unusually high CO2 concentrations of up to 9,160 ppm, which indicate that the upper oceanic mantle is more enriched in carbon than previously thought. We furthermore suggest that CO2 fluxes along mid-ocean ridges vary significantly. Our results demonstrate that elevated fluxes of CO2 from the upper oceanic mantle can drive explosive eruptions at mid-ocean ridges.

Description: Fan-delta uplift and mountain subsidence during the Haiti 2010 earthquake Nature Geoscience 4, 255 (2011). doi:10.1038/ngeo1115 Authors: Manabu Hashimoto, Yo Fukushima & Yukitoshi Fukahata The relative motion between the Caribbean and North American plates is accommodated by several active faults around Hispaniola Island. The Enriquillo–Plantain Garden fault in southern Haiti is one of these structures. Strain equivalent to a magnitude 7.2 earthquake is estimated to have accumulated along this fault since its last significant activity. The Haiti earthquake of 12 January 2010 was initially reported to have occurred along this fault, but more recent studies proposed slips on previously unrecognized, neighbouring faults. Here we use interferometric synthetic aperture radar data to show that surface deformation caused by the earthquake does not correspond to the present topography. Alluvial fan deltas were uplifted on the north side of the Enriquillo–Plantain Garden fault, whereas mountains located on the south side of the fault subsided, implying that faults other than the Enriquillo–Plantain Garden fault were responsible for the deformation. To determine fault structure, we fit the satellite surface deformation data to a fault model. We show that slip occurred on a fault dipping northward at 42°, with large thrust components. The maximum displacement on the fault was about 4 m at 10–20 km depth, offshore from the Tiburon peninsula. We confirm that the earthquake ruptured a blind thrust fault and show that the fault could not be identified from large-scale present-day topography.

Description: Pleistocene water intrusions from the Mediterranean and Caspian seas into the Black Sea Nature Geoscience 4, 236 (2011). doi:10.1038/ngeo1106 Authors: S. Badertscher, D. Fleitmann, H. Cheng, R. L. Edwards, O. M. Göktürk, A. Zumbühl, M. Leuenberger & O. Tüysüz The hydrological balance of the Black Sea is governed by riverine input and by the exchange with the Mediterranean Sea through the shallow Bosporus Strait. These sources have distinctly different oxygen isotope (δ18O) signatures. Therefore, the δ18O of Black Sea water directly reflects the presence or absence of a connection with the Mediterranean Sea, as well as hydrological changes in the vast watersheds of the Black and Caspian seas. However, the timing of late to middle Pleistocene water intrusions to the Black Sea is poorly constrained in sedimentary sequences. Here we present a stacked speleothem δ18O record from Sofular Cave in northern Turkey that tracks the isotopic signature of Black Sea surface water, and thus allows a reconstruction of the precise timing of hydrological shifts of the Black Sea. Our record, which extends discontinuously over the last 670,000 years, suggests that the connection between the Black Sea and Mediterranean Sea has been open for a significant period at least twelve times since 670,000 yr ago, more often than previously suggested. Distinct minima in the Sofular δ18O record indicate at least seven intervals when isotopically depleted freshwater from the Caspian Sea entered the Black Sea. Our data provide precisely dated evidence for a highly dynamic hydrological history of the Black Sea.

Description: Haematite pseudomicrofossils present in the 3.5-billion-year-old Apex Chert Nature Geoscience 4, 240 (2011). doi:10.1038/ngeo1084 Authors: Craig P. Marshall, Julienne R. Emry & Alison Olcott Marshall Microstructures in the ∼3.5 Gyr Apex Chert Formation were initially described as the oldest bacterial fossils on Earth over 20 years ago. However, the identification of the structures (which resemble cyanobacteria) as biological in origin remains controversial. Here we determine the petrology and geochemistry of similar structures from the original Apex Chert locality using thin sections and Raman spectroscopy. Based on the microscopic examination of thin sections, we identify features similar to those previously identified as microfossils as a series of quartz and haematite-filled fractures. Raman spectroscopy of the fractures shows that carbonaceous material is not, as previously reported, associated with the structures, but is instead disseminated in the surrounding quartz matrix. We suggest that although the microstuctures analysed are not microfossils, the presence of carbonaceous material in the surrounding matrix is consistent with the existence of microbial life at this time, and with evidence of early Archaean life found at other sites. Furthermore, we caution against identifying microstructures as biological in origin without a full morphological and geochemical assessment.

Description: Weakness of the San Andreas Fault revealed by samples from the active fault zone Nature Geoscience 4, 251 (2011). doi:10.1038/ngeo1089 Authors: B. M. Carpenter, C. Marone & D. M. Saffer Understanding the strength and slip behaviour of tectonic faults is a central problem in earthquake physics and seismic-hazard assessment. Many major faults, including the San Andreas Fault, are weak compared with the surrounding rock, but the cause of this weakness is debated. Previous measurements of the frictional strength of San Andreas Fault rocks are too high to explain the observed weakness. However, these measurements relied on samples taken at a distance from the active fault or from weathered surface samples. Recent drilling into the San Andreas Fault has provided material from the actively slipping fault at seismogenic depths. Here we present systematic measurements of the frictional properties and composition of the San Andreas Fault at 2.7 km depth, including the wall rock and active fault. We find that the fault is weak relative to the surrounding rock and that the fault rock exhibits stable sliding friction behaviour. The fault zone contains the weak mineral smectite and exhibits no frictional healing—bonds in the material do not heal after rupture. Taken together, the low inherent strength and lack of healing of the fault-zone material could explain why the San Andreas Fault slips by aseismic creep and small earthquakes in central California, rather than by large, destructive earthquakes.

Description: Core surprise Nature Geoscience 4, 268 (2011). doi:10.1038/ngeo1128 Brett Carpenter and colleagues got more than they bargained for when cleaning the mud off a metre-long piece of core from the San Andreas Fault.

Description: Reconciling the hemispherical structure of Earth’s inner core with its super-rotation Nature Geoscience 4, 264 (2011). doi:10.1038/ngeo1083 Authors: Lauren Waszek, Jessica Irving & Arwen Deuss Earth’s solid inner core grows through solidification of material from the fluid outer core onto its surface at rates of about 1 mm per year, freezing in core properties over time and generating an age–depth relation for the inner core. A hemispherical structure of the inner core is well-documented: an isotropic eastern hemisphere with fast seismic velocities contrasts with a slower, anisotropic western hemisphere. Independently, the inner core is reported to super-rotate at rates of up to 1° per year. Considering the slow growth, steady rotation rates of this magnitude would erase ’frozen-in’ regional variation and cannot coexist with hemispherical structure. Here, we exploit the age–depth relation, using the largest available PKIKP–PKiKP seismic travel time data set, to confirm hemispherical structure in the uppermost inner core, and to constrain the locations of the hemisphere boundaries. We find consistent eastward displacement of these boundaries with depth, from which we infer extremely slow steady inner core super-rotation of 0.1°–1° per million years. Our estimate of long-term super-rotation reconciles inner core rotation with hemispherical structure, two properties previously thought incompatible. It is in excellent agreement with geodynamo simulations, while not excluding the possibility that the much larger rotation rates inferred earlier correspond to fluctuations in inner core rotation on shorter timescales.

Description: Spatial and temporal variability of snow accumulation in Dronning Maud Land, East Antarctica, including two deep ice coring sites at Dome Fuji and EPICA DML The Cryosphere Discussions, 5, 2061-2114, 2011 Author(s): S. Fujita, P. Holmlund, I. Andersson, I. Brown, H. Enomoto, Y. Fujii, K. Fujita, K. Fukui, T. Furukawa, M. Hansson, K. Hara, Y. Hoshina, M. Igarashi, Y. Iizuka, S. Imura, S. Ingvander, T. Karlin, H. Motoyama, F. Nakazawa, H. Oerter, L. E. Sjöberg, S. Sugiyama, S. Surdyk, J. Ström, R. Uemura, and F. Wilhelms To better understand the spatio-temporal variability of the glaciological environment in Dronning Maud Land (DML), East Antarctica, investigations were carried out along the 2800-km-long Japanese-Swedish IPY 2007/2008 traverse. The route covers ice sheet ridges and two deep ice coring sites at Dome Fuji and EPICA DML. The surface mass balance (SMB) distribution was derived based on analysis of isochrones within snow pits, firn cores and subsurface radar signals. The SMB averaged over various time scales in the Holocene was determined. This was then compared with various glaciological data. We find that the large-scale distribution of the SMB depends on the surface elevation, continentality and interactions between ice sheet ridges and the prevailing counterclockwise windfield in DML. A different SMB is found for the windward and leeward sides of the ridges. Local-scale variability in the SMB is essentially governed by bedrock topography which determines the local surface topography. In the eastern part of DML, the accumulation rate in the second half of the 20th century is found to be higher by 15 % compared to averages over longer periods of 722 a or 7.9 ka before AD 2008. A similar trend has been reported for many inland plateau sites in East Antarctica.

Description: Find faults and funds Nature Geoscience 4, 345 (2011). doi:10.1038/ngeo1180 Seismic risk is poorly known in many places on Earth. To save lives it is necessary — but by no means sufficient — to map the faults that pose a threat more accurately.

Description: Not just family matters Nature Geoscience 4, 346 (2011). doi:10.1038/ngeo1163 Author: Helena L. Filipsson

Description: Uncharted seismic risk Nature Geoscience 4, 348 (2011). doi:10.1038/ngeo1168 Authors: Philip England & James Jackson The potential for earthquakes along the plate boundaries has been mapped with reasonable success. Our attention should now focus on the threat posed by unanticipated quakes located in the continental interiors.

Description: Ocean science: If Gaia could talk Nature Geoscience 4, 351 (2011). doi:10.1038/ngeo1175 Author: Maurice Levasseur Sulphur aerosols contribute to the regulation of the Earth's climate. Biogenic sulphur emissions from the Southern Ocean have been underestimated, and could be further enhanced owing to climate warming, according to observations and model simulations.

Description: Our choice from the recent literature Nature Geoscience 4, 350 (2011). doi:10.1038/ngeo1178

Description: Atmospheric science: Asia under a high-level brown cloud Nature Geoscience 4, 352 (2011). doi:10.1038/ngeo1166 Author: Mark G. Lawrence Gaseous pollutants such as ozone and carbon monoxide from Asia are lifted to altitudes of more than 10 km during the summer monsoon season. Satellite observations show that aerosol particles, too, can rise high and spread across thousands of kilometres.

Description: Biogenicity of Apex Chert microstructures Nature Geoscience 4, 346 (2011). doi:10.1038/ngeo1172 Authors: J. William Schopf & Anatoliy B. Kudryavtsev

Description: Geography: Strabo's Greece Nature Geoscience 4, 355 (2011). doi:10.1038/ngeo1177 Author: Alicia Newton

Description: Planetary science: Titan's nitrogenesis Nature Geoscience 4, 356 (2011). doi:10.1038/ngeo1162 Author: Catherine Neish Observations from the Cassini–Huygens mission have produced potentially contradictory constraints on the origin of Titan's atmosphere. Experiments and a simple model demonstrate that a new mechanism for late formation is plausible.

Description: Tectonics: Unexpected Andean earthquakes Nature Geoscience 4, 357 (2011). doi:10.1038/ngeo1167 Author: Nadine McQuarrie Great earthquakes along the western, subduction zone boundary of the Andes Mountains in South America are expected. Measurements of surface motion along the eastern boundary highlight the potential for equally large earthquakes in the east.