Fast switch from extensional exhumation to thrusting of the Ronda Peridotites (South Spain)

The Alboran Domain, situated at the western end of the Mediterranean subduction system, is characterized by the Ronda Peridotites, one of the world's largest exposures of sub‐continental mantle. Using U–Pb (LA‐ICP‐MS) and Ar–Ar dating, we precisely dated two tectonic events associated with the Tertiary exhumation of the Ronda Peridotites. First, shearing along the Crust–Mantle Extensional Shear Zone caused, at ca. 22.5 Ma, mantle exhumation, local partial melting in the deep crust and coeval cooling in the upper crust. Second, the Ronda Peridotites Thrust triggered the final emplacement of the peridotites onto the continental crust at c. 21 Ma, as testified by granitic intrusions in the thrust hangingwall. The tectonic evolution of the western Alboran Domain is therefore characterized by a fast switch from continental lithospheric extension in a backarc setting, with sub‐continental mantle exhumation, to a rift inversion by thrusting driven by shortening of the upper plate.     

Annual fossil organic carbon delivery due to mechanical and chemical weathering of marly badlands areas

A key issue in the study of the carbon cycle is constraining the stocks and fluxes in and between C‐reservoirs. Among these, the role and importance of fossil organic carbon (FOC) release by weathering of outcropping sedimentary rocks on continental surfaces is still debated and remains poorly constrained. Our work focuses on FOC fluxes due to chemical and mechanical weathering of marls in two experimental watersheds with typical badlands geomorphology (Draix watersheds, Laval and Moulin, Alpes de Haute Provence, France). Organic matter from bedrock, soil litter and riverine particles are characterized by Rock‐Eval 6 pyrolysis. FOC fluxes due to mechanical weathering are then estimated by monitoring the annual particulate solid exports at the outlets of the watersheds (1985–2005 period). FOC fluxes from chemical weathering were calculated using Ca²⁺ concentrations in dissolved loads (year 2002) to assess the amount of FOC released by the dissolution of the carbonate matrix. Results show that FOC delivery is mainly driven by mechanical weathering, with a yield ranging from 30 to 59 t km^‐2 yr^‐1 in the Moulin (0.08 km²) and Laval (0.86 km²) catchments, respectively, (1985–2005 average). The release of FOC attributed to chemical weathering was 2.2 to 4.2 t km^‐2 for the year 2002. These high FOC fluxes from badlands are similar to those observed in tectonically active mountain catchments. At a regional scale, badland outcropping within the Durance watershed does not exceed 0.25% in area of the Rhône catchment, but could annually deliver 12 000 t yr^‐1 of FOC. This flux could correspond to 27% of the total particulate organic carbon (POC) load exported by the Rhône River to the Mediterranean Sea. At a global scale, our findings suggest that erosion of badlands may contribute significantly to the transfer of FOC from continental surfaces to depositional environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Hierarchy of non-extensive mechanical processes in fracturing sea ice

The fracture-induced oscillations in sea ice were detected by seismographs and seismic tiltmeters established on the Arctic ice pack. Field observations were supplemented with a laboratory experiment. The energy distributions in elastic waves generated during: (i) large-scale ice pack fragmentation over area of about 10⁵ km², (ii) local crack propagation in ice floe, and (iii) laboratory ice crashing were constructed and analyzed using principles of the Tsallis statistics. The energy release regimes at different stages of fracturing were characterized by the parameter of nonextensivity q. In terms of the non-extensive statistical mechanics, q > 1 evidences the correlated (non-extensive) dynamics of the process in nonequilibrium system, q = 1 responds to the additivity of events occurring in equilibrium system, and q < 1 takes place when the energy release is additive and limited by an upper cut-off. All these scenarios were revealed in fracture processes occurring at three hierarchic levels. The variation of the q-value demonstrates high thermodynamic changeability of the fracture process driven by irregular external source. The role of energy conservation in fracturing sea ice is discussed in connection with the observed reversible transitions between extensive and non-extensive modes of fracture.     

Simulation of masonry out‐of‐plane failure modes by multi‐body dynamics

The seismic assessment of the local failure modes in existing masonry buildings is currently based on the identification of the so‐called local mechanisms, often associated with the out‐of‐plane wall behavior, whose stability is evaluated by static force‐based approaches and, more recently, by some displacement‐based proposals. Local mechanisms consist of kinematic chains of masonry portions, often regarded as rigid bodies, with geometric nonlinearity and concentrated nonlinearity in predefined contact regions (unilateral no‐tension behavior, possible sliding with friction). In this work, the dynamic behavior of local mechanisms is simulated through multi‐body dynamics, to obtain the nonlinear response with efficient time history analyses that directly take into account the characteristics of the ground motion. The amplification/filtering effects of the structure are considered within the input motion. The proposed approach is validated with experimental results of two full‐scale shaking‐table tests on stone masonry buildings: a sacco‐stone masonry façade tested at Laboratório Nacional de Engenharia Civil and a two‐storey double‐leaf masonry building tested at European Centre for Training and Research in Earthquake Engineering (EUCENTRE). Copyright © 2015 John Wiley & Sons, Ltd.     

Titan's highly dynamic magnetic environment: A systematic survey of Cassini magnetometer observations from flybys TA-T62

We analyze the variability of the ambient magnetic field near Titan during Cassini encounters TA-T62 (October 2004-October 2009). Cassini magnetometer (MAG) data show that the moon's magnetic environment is strongly affected by its proximity to Saturn's warped and highly dynamic magnetodisk. In the nightside sector of Saturn's magnetosphere, the magnetic field near Titan is controlled by intense vertical flapping motions of the magnetodisk current sheet, alternately exposing the moon to radially stretched lobe-type fields and to more dipolar, but highly distorted current sheet fields. In southern summer, when most of the Cassini encounters took place, the magnetodisk current sheet was on average located above Titan's orbital plane. However, around equinox in August 2009, the distortions of Titan's magnetic environment due to the rapidly moving current sheet reached a maximum, thus suggesting that the equilibrium position of the sheet at that time was significantly closer to the moon's orbital plane. In the dayside magnetosphere, the formation of the magnetodisk lobes is partially suppressed due to the proximity of the magnetopause. Therefore, during most encounters that took place near noon, Titan was embedded in highly distorted current sheet fields. Within the framework of this study, we not only provide a systematic classification of all Titan flybys between October 2004 and October 2009 as lobe-type or current sheet scenarios, but we also calculate the magnetospheric background field near Titan's orbit whenever possible. Our results show that so far, there is not a single Cassini flyby that matches the frequently applied picture of Titan's plasma interaction from the pre-Cassini era (background field homogeneous, stationary and perpendicular to the moon's orbital plane). The time scales upon which the ambient magnetospheric field close to Titan undergoes significant changes range between only a few minutes and up to several hours. The implications for the development of numerical models for Titan's local plasma interaction are discussed as well.     

The Transfiguration continental red-bed Cu–Pb–Zn–Ag deposit,Quebec Appalachians,Canada

The Transfiguration Cu–Pb–Zn–Ag deposit, enclosed within reduced grey sandstone, is associated with continental red beds of the Lower Silurian Robitaille Formation in the Quebec Appalachians, Canada. The Robitaille Formation rests unconformably on foliated Cambro-Ordovician rocks. The unconformity is locally cut by barite veins. The basal unit of the Robitaille Formation comprises green wacke and pebble conglomerate, which locally contain calcite nodules. The latter have microstructures characteristic of alpha-type calcretes, such as “floating” fabrics, calcite-filled fractures (crystallaria) and circumgranular cracks. Massive, grey sandstone overlies the basal green wacke and pebble conglomerate unit, which is overlain, in turn, by red, fine-grained sandstone. Mineralisation occurred underneath the red sandstone unit, chiefly in the grey sandstone unit, as disseminated and veinlet sulphides. Chalcopyrite, the most abundant Cu sulphide, replaced early pyrite. Calcrete, disseminated carbonate and vein carbonate have stable isotope ratios varying from −7.5‰ to −1.1‰ δ¹³C and from 14.7‰ to 21.3‰ δ¹⁸O. The negative δ¹³C values indicate the oxidation of organic matter in a continental environment. Sulphur isotope ratios for pyrite, chalcopyrite and galena vary from −19‰ to 25‰ δ³⁴S, as measured on mineral concentrates by a conventional SO₂ technique. Laser-assisted microanalyses (by fluorination) of S isotopes in pyrite show an analogous range in δ³⁴S values, from −21‰ to 25‰. Negative and positive δ³⁴S values are compatible with bacterial sulphate reduction (BSR) in systems open and closed with respect to sulphate. We interpret similarly high δ³⁴S values for sulphide concentrates (25.1‰) and for vein barite (26.2‰) to result from rapid and complete thermochemical reduction of pore-water sulphate. Two early to late diagenetic stages of mineralisation best explain the origin of the Transfiguration deposit. The first stage was characterised by the ponding of groundwater over the Taconian unconformity, recorded by calcrete and early pyrite formation via BSR in grey sandstone. Early pyrite contains up to 2 wt.% Pb, which is consistent with Pb fixation by sulphate-reducing bacteria. The second stage (II) is defined by the replacement of early pyrite by chalcopyrite, as well as by sulphide precipitation via either BSR or thermochemical sulphate reduction (TSR) in grey sandstone. This event resulted from the synsedimentary fault-controlled percolation and mixing of (1) an oxidising, sulphate-bearing cupriferous fluid migrating per descensum from the red-bed sequence and (2) a hydrocarbon-bearing fluid migrating per ascensum from the Cambro-Ordovician basement. Mixing between the two fluids led to sulphate reduction, causing Cu sulphide precipitation. The positive correlation between Cu and Fe³⁺/Fe²⁺ bulk rock values suggests that Fe acted as a redox agent during sulphate reduction. Stage II diagenetic fluid migration is tentatively attributed to the Late Silurian Salinic extensional event.     

On elevated fluoride and boron concentrations in groundwaters associated with the Lake Saint-Martin impact structure,Manitoba

Hydrogeological investigations conducted by the Geological Survey of Canada in the Lake Saint-Martin region of Manitoba have confirmed earlier reports of naturally elevated F⁻ and B concentrations in local groundwaters. Fluoride and B concentrations are highly correlated (r² = 0.905) and reach 15.1 mg/L and 8.5 mg/L, respectively. Virtually all groundwaters with F⁻ concentrations greater than the drinking water limit of 1.5 mg/L are from wells within the Lake Saint-Martin impact structure, a 208 Ma complex crater 23 km in diameter underlying a large part of the study area. The high-F⁻ groundwaters can be classified into two groups according to their anionic and isotopic compositions. Group I samples consist of Na-mixed anion groundwaters, with Cl greater than 100 mg/L and highly depleted ¹⁸O compositions indicative of recharge under much cooler climatic conditions than at present. Samples belonging to this group exhibit a striking relationship to crater morphology, and are found in an arcuate belt within the southern rim of the impact structure. Group II high-F⁻ samples consist of Na–HCO₃–SO₄ groundwaters, with little Cl, and less depleted ¹⁸O compositions. Samples belonging to this group are associated with groundwaters recharged locally, on a low ridge within the impact structure. This paper traces the probable source of high-F⁻ groundwaters to phosphatic pellets in shales of the Winnipeg Formation, a regional basal clastic unit which sub-crops at shallow depth beneath the crater rim as a result of more than 200 m of structural uplift associated with the impact event. This extensive aquifer is known elsewhere in southern Manitoba for its naturally-softened groundwaters and locally elevated F⁻ concentrations. Group I groundwaters are interpreted as discharge from the Winnipeg Formation where it abuts against crater-fill deposits. Group II high-F⁻ groundwaters are interpreted as modern recharge from within the impact structure, displacing Group I groundwaters. Thus, elevated F⁻ and B concentrations observed in groundwaters of the Lake Saint-Martin area represent the geochemical signature of upwelling from a deep regional aquifer. The previously unsuspected discharge zone occurs within an isolated sub-crop of the aquifer formed as a result of structural uplift caused by the impact event.     

Block Simulation of Multiple Correlated Variables

Numerical representations of multivariate natural phenomena, including characteristics of mineral deposits, petroleum reservoirs and geo-environmental attributes, need to consider and reproduce the spatial relationships between correlated attributes of interest. There are, however, only a few methods that can practically jointly simulate large size multivariate fields. This paper presents a method for the conditional simulation of a non-Gaussian vector random field directly on block support. The method is derived from the group sequential simulation paradigm and the direct block simulation algorithm which leads to the efficient joint simulation of large multivariate datasets jointly and directly on the block support. This method is a multistage process. First, a vector random function is orthogonalized with minimum/maximum autocorrelation factors (MAF). Blocks are then simulated by performing LU simulation on their discretized points, which are later back-rotated and averaged to yield the block value. The internal points are then discarded and only the block value is stored in memory to be used for further conditioning through a joint LU, resulting in the reduction of memory requirements. The method is termed direct block simulation with MAF or DBMAFSIM. A proof of the concept using an exhaustive data set demonstrates the intricacies and the performance of the proposed approach.     

Internal structure of Piton de la Fournaise volcano from seismic wave propagation and earthquake distribution

Arrival times of seismic waves from local earthquakes are inverted for both locating the source and defining the 3-D velocity heterogeneity of Piton de la Fournaise.The lateral heterogeneity of the 2632 m high edifice is resolved as a high-velocity plug, 1.5 km in diameter, surrounded by a low-velocity ring, which may be interpreted as due to the construction of Fournaise on the flank of the older volcano Piton des Neiges. Wave mode conversion detected on three-component seismograms provides evidence for boundaries of contrasted velocities.Pre-eruptive swarm earthquakes cluster in the high-velocity zone, under the Dolomieu summit crater. Low strength and cohesion of the surrounding material account for the lack of seismicity for the final 1–3 km radial flow of magma to the vents in Enclos Fouqué.Beneath the high-velocity plug the existence of a body with low velocity for P, and even for S, waves is well constrained. However, the walls and base are poorly defined because of the lack of deep earthquakes for sampling. The few earthquakes that are located in this depth region usually occur at a depth of around 1.5 km below sea level in the region of the cone. This can be considered providing the upper constraint on the lower limit of the aseismic part of the low-velocity body. The coincidence in time of their occurrence with the swarms above sea level and the eruptions suggests magmatic activation of the low-velocity aseismic volume 1.5 km below sea level under the high-velocity plug of the cone. Further down, the concentration of seismicity in two swarms, between 2 and 4 km, under the eastern flank does not allow the structure to be sampled effectively.