The origins of Antarctic rock glaciers: periglacial or glacial features?

In extensively glaciarized permafrost areas such as Northern Victoria Land, rock glaciers are quite common and are considered postglacial cryotic landforms. This paper reveals that two rock glaciers in Northern Victoria Land (at Adélie Cove and Strandline) that are located close to the Italian Antarctic Station (Mario Zucchelli Station) should have the same origin, although they were previously mapped as Holocene periglacial landforms and subsequently considered ice‐cored and ice‐cemented rock glaciers, respectively. In fact, by integrating different geophysical investigations and borehole stratigraphy, we show that both landforms have similar internal structures and cores of buried glacier ice. Therefore, this kind of rock glacier is possibly related to the long‐term creep of buried ice rather than to permafrost creep alone. This interpretation can be extended to the larger part of the features mapped as rock glaciers in Antarctica. In addition, a high‐reflective horizon sub‐parallel to the topographic surface was detected in Ground Probing Radar (GPR) data over a large part of the study area. Combining all the available information, we conclude that it cannot be straightforwardly interpreted as the base of the active layer but rather represents the top of a cryo‐lithological unit characterized by ice lenses within sediments that could be interpreted as the transition zone between the active layer and the long‐term permafrost table. More generally, knowledge of the subsurface ice content and, in particular, the occurrence of massive ice and its depth is crucial to make realistic and affordable forecasts regarding thermokarst development and related feedbacks involving GHG emissions, especially in the case of cryosoils rich in carbon content. Copyright © 2017 John Wiley & Sons, Ltd.     

A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.     

Thermal expansion and phase transitions in åkermanite and gehlenite

Thermal expansion has been measured by laboratory and synchrotron X-ray powder diffraction for end-member åkermanite (ak, Ca₂MgSi₂O₇) and gehlenite (ge, Ca₂Al₂SiO₇) in the range 20–1,500 K. In ak in the range 340–390 K, there is a negative linear thermal expansion in [001] direction. This is related to the phase transition from an incommensurate modulated structure (IC) to a normal one (N). The volumetric mean thermal expansion coefficients for ak and ge, obtained with a linear fit of the experimental data in the temperature range 298–1,400 K, are respectively 32.1×10^–6 and 28.3×10^–6 K^–1 . The variation of the c/a ratio with temperature, due to different thermal expansion along the crystallographic axes, can be related to the different behaviour of the tetrahedral layers in the N and IC phases. Analysis of the variation of the superstructure peaks intensity across the phase transition confirms the tricritical behaviour of the IC/N transition in ak.     

A new view of the He-Ar-CO2 degassing at mid-ocean ridges: Homogeneous composition of magmas from the upper mantle

Deep-sea exploration is rapidly improving our understanding of volatiles geochemistry in mid-ocean-ridge igneous products. It is also placing greater constraints on degassing processes of the Earth’s mantle, with the result that degassing models based on vapour-melt equilibrium are no longer able to explain the increasing number of data. In fact, such models force to postulate an upper mantle strongly heterogeneous at any scale, and cannot account for the widespread carbon supersaturation of the recovered igneous products. Here we review the global He-Ar-CO₂ dataset of fluid inclusions in mid-ocean-ridge glasses using the framework of advanced modelling of multicomponent bubble growth in magmas. We display that non-equilibrium fractionations among He, Ar and CO₂, driven by their different diffusivities in silicate melts, are common in most of the natural conditions of magma decompression and their signature strongly depends on pressure of degassing. Due to the comparable Ar and CO₂ diffusivity, magma degassing at low pressure fractionates both the He/Ar and He/CO₂ ratio by a similar extent, while the slower CO₂ diffusion at high pressure causes early kinetic effects on Ar/CO₂ ratio and dramatically changes the degassing path. On this ground, the very different geochemical signatures among suites of data coming from different ridge segments mainly depend on the depth of the magma chamber where the melt was stored. Besides, the variations inside a single suite highlight variable ascent speed and cooling rate of the emplaced lava. The large variations in both the He/CO₂ and Ar/CO₂ ratios at almost constant He/Ar, displayed in glasses coming from the Mid-Atlantic Ridge 24-30°N segment and the Rodriguez Triple Junction, are therefore interpreted as a high-pressure signature. In contrast, the simultaneous increase in both He/CO₂ and He/Ar of the East Pacific Rise, Pito Seamount and South-East Indian Ridge data sets suggests the dominance of low-pressure fractionation, implying that the shallow magma chambers are at a lower depth than those of the Mid-Atlantic Ridge 24-30°N and Rodriguez Triple Junction. Our conclusions support the presence of a relationship between spreading rate and depth of high-temperature zones below ridges, and are consistent with the depth of magma chambers as suggested from seismic studies. Non-equilibrium degassing explains the volatile systematics of mid-ocean-ridge basalts by starting from a single mantle-derived magma, dispensing with the supposed need for heterogeneities in abundance ratios of volatiles in the mantle below oceanic ridges.     

Interpretation of the solar-like pulsational behaviour of η Bootis

HR 5235, better known as η Bootis, is a bright and well-known star for which very accurate observations have recently enabled Kjeldsen et al. (2003) and Carrier, Bouchy, and Eggenberger (2003) not only to confirm the presence of solar-like oscillations, but also to identify the excitation in the oscillation spectrum of several p-mode frequencies with harmonic degrees l = 0 – 2. Here we show how such observational success, through the calculation and the investigation of theoretical structure models and the comparison of the observed oscillation spectra with the predicted p-mode frequencies of oscillations, permits one to draw conclusions about the actual evolutionary state of this star and on the physical properties of its internal structure. The computation of the structure models is based on the use of updated global parameters and includes overshooting from the convective core. In particular, we consider the effect on the stellar structure, and hence on the theoretical frequencies, of employing different equations of state and different formalisms to describe the convective energy transport.     

Spitzer IRAC infrared colours of submillimetre-bright galaxies

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the Infra Red Array Camera (IRAC) colour–colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (≲100 Myr old), extreme starbursts at z > 2 also occupy the same general area as AGNs in the IRAC colour–colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (≳1 mag bluer in the observed [4.5]–[8.0] colour). Only about 20 per cent of the objects overlap in the colour–colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red infrared colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of infrared colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Multiband Imaging Photometer for Spitzer (MIPS) 24 μm data.     

Plan‐form evolution of ancient meandering rivers reconstructed from longitudinal outcrop sections

The mode of channel‐bend transformation (i.e. expansion, translation, rotation or a combination thereof) has a direct bearing on the dimensions, shape, bedding architecture and connectivity of point‐bar sandstone bodies within a fluvial meander belt, but is generally difficult to recognize in vertical outcrops. This study demonstrates how the bend transformation mode and relative rate of channel‐floor aggradation can be deciphered from longitudinal outcrop sections aligned parallel to the meander‐belt axis, as a crucial methodological aid to the reconstruction of ancient fluvial systems and the development of outcrop analogue models for fluvial petroleum reservoirs. The study focuses on single‐storey and multi‐storey fluvial meander‐belt sandstone bodies in the Palaeogene piggyback Boyabat Basin of north‐central Turkey. The sandstone bodies are several hundred metres wide, 5 to 40 m thick and encased in muddy floodplain deposits. The individual channel‐belt storeys are 5 to 9 m thick and their transverse sections show lateral‐accretion bed packages representing point bars. Point bars in longitudinal sections are recognizable as broad mounds whose parts with downstream‐inclined, subhorizontal and upstream‐inclined bedding represent, respectively, the bar downstream, central and upstream parts. The inter‐bar channel thalweg is recognizable as the transition zone between adjacent point‐bar bedsets with opposing dip directions into or out of the outcrop section. The diverging or converging adjacent thalweg trajectories, or a trajectory migrating in up‐valley direction, indicate point‐bar broadening and hence channel‐bend expansion. A concurrent down‐valley migration of adjacent trajectories indicates channel‐bend translation. Bend rotation is recognizable from the replacement of a depositional riffle by an erosional pool zone or vice versa along the thalweg trajectory. The steepness of the thalweg trajectory reflects the relative rate of channel‐floor aggradation. This study discusses further how the late‐stage foreland tectonics, with its alternating pulses of uplift and subsidence and a progressive narrowing of the basin, has forced aggradation of fluvial channels and caused vertical stacking of meander belts.     

Last glacial maximum environments in northwestern Patagonia revealed by fossil small mammals

Comparisons of historical and modern assemblages of mammals can yield important insights into patterns and processes of environmental change. Here, we present the first analyses of small mammal assemblages present in northern Patagonia during the last glacial maximum (LGM). Using remains obtained from owl pellets excavated from an archeological cave site (Arroyo Corral I, levels VII–V, carbon dates of 22,400–21,530 cal yr BP), we generate estimates of the minimum number of individuals for all species detected; these estimates, in turn are used to determine relative species abundances. Comparisons of these data with similar analyses of small mammal remains obtained from a second archeological site (ACoII, levels IV–V, carbon dates of 10,010–9220 cal yr BP) as well as from modern owl pellets reveal pronounced changes in relative species abundance since the LGM. In particular, Euneomys chinchilloides and Ctenomys sociabilis – the predominant species during the LGM – declined markedly, suggesting a change from open, bare habitat punctuated by patches of wet meadows and shrubs to the more densely vegetated mosaic of ecotone habitats found in this region today. These data provide important new insights into the environmental changes that have occurred in northern Patagonia over the last 20,000 years.