Distinguishing volcanic debris avalanche deposits from their reworked products: the Perrier sequence (French Massif Central)

Debris avalanches associated with volcanic sector collapse are usually high-volume high-mobility phenomena. Debris avalanche deposit remobilisation by cohesive debris flows and landslides is common, so they can share textural characteristics such as hummocks and jigsaw cracks. Distinguishing original deposits from reworked products is critical for geological understanding and hazard assessment because of their different origin, frequency and environmental impact. We present a methodology based on field evidence to differentiate such epiclastic breccias. Basal contact mapping constrained by accurate altitude and location data allows the reconstruction of deposit stratigraphy and geometry. Lithological analysis helps to distinguish the different units. Incorporation structures, kinematic indicators and component mingling textures are used to characterise erosion and transport mechanisms. We apply this method to the enigmatic sequence at Perrier (French Massif Central), where four units (U1–U4) have been interpreted either as debris flow or debris avalanche deposits. The sequence results from activity on the Monts Dore Volcano about 2 Ma ago. The epiclastic units are matrix supported with an almost flat top. U2 and U3 have clear debris flow deposit affinities such as rounded clasts and intact blocks (no jigsaw cracks). U1 and U4 have jigsaw cracked blocks with matrix injection and stretched sediment blocks. U1 lacks large blocks (>10 m wide) and has a homogenous matrix with an upward increase of trapped air vesicle content and size. This unit is interpreted as a cohesive debris flow deposit spawned from a debris avalanche upstream. In contrast, U4 has large mega-blocks (up to 40 m wide), sharp contacts between mixed facies zones with different colours and numerous jigsaw fit blocks (open jigsaw cracks filled by monogenic intra-clast matrix). Mega-blocks are concentrated near the deposit base and are spatially associated with major substratum erosion. This deposit has a debris avalanche distal facies with local debris flow affinities due to partial water saturation. We also identify two landslide deposits (L1 and L2) resulting from recent reworking that has produced a similar facies to U1 and U4. These are distinguishable from the original deposits, as they contain blocks of mixed U1/U4 facies, a distinctly less consolidated and more porous matrix and a fresh hummocky topography. This work shows how to differentiate epiclastic deposits with similar characteristics, but different origins. In doing so, we improve understanding of present and past instability of the Monts Dore and identify present landslide hazards at Perrier.     

Spectral density regression for bivariate extremes

We introduce a density regression model for the spectral density of a bivariate extreme value distribution, that allows us to assess how extremal dependence can change over a covariate. Inference is performed through a double kernel estimator, which can be seen as an extension of the Nadaraya–Watson estimator where the usual scalar responses are replaced by mean constrained densities on the unit interval. Numerical experiments with the methods illustrate their resilience in a variety of contexts of practical interest. An extreme temperature dataset is used to illustrate our methods.     

A Snapshot of the Sun Near Solar Minimum: The Whole Heliosphere Interval

We present an overview of the data and models collected for the Whole Heliosphere Interval, an international campaign to study the three-dimensional solar?Cheliospheric?Cplanetary connected system near solar minimum. The data and models correspond to solar Carrington Rotation 2068 (20 March??C?16 April 2008) extending from below the solar photosphere, through interplanetary space, and down to Earth??s mesosphere. Nearly 200 people participated in aspects of WHI studies, analyzing and interpreting data from nearly 100 instruments and models in order to elucidate the physics of fundamental heliophysical processes. The solar and inner heliospheric data showed structure consistent with the declining phase of the solar cycle. A closely spaced cluster of low-latitude active regions was responsible for an increased level of magnetic activity, while a highly warped current sheet dominated heliospheric structure. The geospace data revealed an unusually high level of activity, driven primarily by the periodic impingement of high-speed streams. The WHI studies traced the solar activity and structure into the heliosphere and geospace, and provided new insight into the nature of the interconnected heliophysical system near solar minimum.     

Potential for small scale added value of RCM’s downscaled climate change signal

In recent decades, the need of future climate information at local scales have pushed the climate modelling community to perform increasingly higher resolution simulations and to develop alternative approaches to obtain fine-scale climatic information. In this article, various nested regional climate model (RCM) simulations have been used to try to identify regions across North America where high-resolution downscaling generates fine-scale details in the climate projection derived using the “delta method”. Two necessary conditions were identified for an RCM to produce added value (AV) over lower resolution atmosphere-ocean general circulation models in the fine-scale component of the climate change (CC) signal. First, the RCM-derived CC signal must contain some non-negligible fine-scale information—independently of the RCM ability to produce AV in the present climate. Second, the uncertainty related with the estimation of this fine-scale information should be relatively small compared with the information itself in order to suggest that RCMs are able to simulate robust fine-scale features in the CC signal. Clearly, considering necessary (but not sufficient) conditions means that we are studying the “potential” of RCMs to add value instead of the AV, which preempts and avoids any discussion of the actual skill and hence the need for hindcast comparisons. The analysis concentrates on the CC signal obtained from the seasonal-averaged temperature and precipitation fields and shows that the fine-scale variability of the CC signal is generally small compared to its large-scale component, suggesting that little AV can be expected for the time-averaged fields. For the temperature variable, the largest potential for fine-scale added value appears in coastal regions mainly related with differential warming in land and oceanic surfaces. Fine-scale features can account for nearly 60 % of the total CC signal in some coastal regions although for most regions the fine scale contributions to the total CC signal are of around ～5 %. For the precipitation variable, fine scales contribute to a change of generally less than 15 % of the seasonal-averaged precipitation in present climate with a continental North American average of ～5 % in both summer and winter seasons. In the case of precipitation, uncertainty due to sampling issues may further dilute the information present in the downscaled fine scales. These results suggest that users of RCM simulations for climate change studies in a delta method framework have little high-resolution information to gain from RCMs at least if they limit themselves to the study of first-order statistical moments. Other possible benefits arising from the use of RCMs—such as in the large scale of the downscaled fields– were not explored in this research.     

The effect of As, Co, and Ni impurities on pyrite oxidation kinetics: Batch and flow-through reactor experiments with synthetic pyrite

Pyrite samples synthesized with As, Co, or Ni impurities and without added impurities were oxidized in batch and mixed flow-through reactors in the presence of 1 mM ferric iron, at pH 2. Six samples from each dopant population were used to provide a statistically robust comparison; two natural samples from Leadville, CO (major impurities Pb, As, Bi, Ag, Zn) and Elba, Italy (Co, As) were also included. In each experiment, three reaction progress variables were monitored: ferric iron, ferrous iron, and sulfate. The pyrite samples with impurities have average oxidation rates that are faster than the undoped samples, with As- and Co-doped pyrite having the highest rates. As, Co, and Ni were released to solution in accordance with their concentrations in the solid samples. As concentrations in the batch reactor experiments tended to remain constant, in contrast to Co and Ni, which increased over time. Initial rates, calculated from the batch reactor experiments, were faster than the steady-state rates calculated from the mixed flow-through reactor experiments. Apparent rates calculated using sulfate were faster than apparent rates calculated using ferric and ferrous iron, reflecting oxidation of ferrous iron in solution by dissolved oxygen. The results imply that impurities in pyrite do contribute to its reactivity, in agreement with studies using electrochemical methods. Oxidation rate differences among pyrite samples with different impurities are probably too small to warrant explicit consideration in environmental modeling applications, but are important to understanding pyrite oxidation mechanisms and semiconducting properties.     

Ice and climate modeling: An editorial essay

The growth and decay of ice sheets are driven by forces affecting the seasonal cycles of snowfall and snowmelt. The external forces are likely to be variations in the earth's orbit which cause differences in the solar radiation received. Radiational control of snowmelt is modulated by the seasonal cycles of snow albedo and cloud cover. The effects of orbital changes can be magnified by feedbacks involving atmospheric CO₂ content, ocean temperatures and desert areas. Climate modeling of the causes of the Pleistocene ice ages involves modeling the interactions of all components of the climate system; snow, sea ice, glacier ice, the ocean, the atmosphere, and the solid earth. Such modeling is also necessary for interpreting oxygen isotope records from ice and ocean as paleoclimatic evidence.     

Finite element implementation of a homogenized constitutive law for stone column-reinforced foundation soils,with application to the design of structures

Making use of the homogenization method for periodic media, developed in the context of elastoplasticity, a simplified constitutive law is proposed for a stone column-reinforced soil, regarded as a homogeneous but anisotropic medium. The closed-form expressions derived for such a constitutive law allow for its implementation into a f.e.m-based numerical procedure. The computational code so obtained is then applied to simulating the response of a foundation soil reinforced by a group of floating columns, expressed in terms of load–settlement curves drawn up to the ultimate bearing capacity.     

Late Pleistocene glaciation of the Mt Giluwe volcano,Papua New Guinea

The Mt Giluwe shield volcano was the largest area glaciated in Papua New Guinea during the Pleistocene. Despite minimal cooling of the sea surface during the last glacial maximum, glaciers reached elevations as low as 3200 m. To investigate changes in the extent of ice through time we have re-mapped evidence for glaciation on the southwest flank of Mt Giluwe. We find that an ice cap has formed on the flanks of the mountain on at least three, and probably four, separate occasions. To constrain the ages of these glaciations we present 39 new cosmogenic ³⁶Cl exposure ages complemented by new radiocarbon dates. Direct dating of the moraines identifies that the maximum extent of glaciation on the mountain was not during the last glacial maximum as previously thought. In conjunction with existing potassium/argon and radiocarbon dating, we recognise four distinct glacial periods between 293–306 ka (Gogon Glaciation), 136–158 ka (Mengane Glaciation), centred at 62 ka (Komia Glaciation) and from >20.3–11.5 ka (Tongo Glaciation). The temperature difference relative to the present during the Tongo Glaciation is likely to be of the order of at least 5 °C which is a minimum difference for the previous glaciations. During the Tongo Glaciation, ice was briefly at its maximum for less than 1000 years, but stayed near maximum levels for nearly 4000 years, until about 15.4 ka. Over the next 4000 years there was more rapid retreat with ice free conditions by the early Holocene.