Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

Volcanic edifice stability during cryptodome intrusion

Limit equilibrium analyses were applied to the 1980 Mount St. Helens and 1956 Bezymianny failures in order to examine the influence on stability of structural deformation produced by cryptodome emplacement. Weakening structures associated with the cryptodome include outward-dipping normal faults bounding a summit graben and a flat shear zone at the base of the bulged flank generated by lateral push of the magma. Together with the head of the magmatic body itself, these structures serve directly to localize failure along a critical surface with low stability deep within the interior of the edifice. This critical surface, with the safety coefficient reduced by 25-30%, is then very sensitive to stability condition variation, in particular to the pore-pressure ratio (r_u) and seismicity coefficient (n). For r_u=0.3, or n=0.2, the deep surface suffers catastrophic failure, removing a large volume of the edifice flank. In the case of Mount St. Helens, failure occurred within a material with angle of friction ~40°, cohesion in the range 10⁵-10⁶ Pa, and probably significant water pore pressure. On 18 May 1980, detachment of slide block I occurred along a newly formed rupture surface passing through the crest of the bulge. Although sliding of block I may have been helped by the basal shear zone, significant pore pressure and a triggering earthquake were required (r_u=0.3 and n=0.2). Detachment of the second block was guided by the summit normal fault, the front of the cryptodome, and the basal shear zone. This occurred along a deep critical surface, which was on the verge of failure even before the 18 May 1980 earthquake. The stability of equivalent surfaces at Bezymianny Volcano appears significantly higher. Thus, although magma had already reached the surface, weaker materials, or higher pore pressure and/or seismic conditions were probably required to reach the rupture threshold. From our analysis, we find that deep-seated sector collapses formed by removing the edifice summit cannot generally result from a single slide. Cryptodome-induced deformation does, however, provide a deep potential slip surface. As previously thought, it may assist deep-seated sector collapse because it favors multiple retrogressive slides. This leads to explosive depressurization of the magmatic and hydrothermal systems, which undermines the edifice summit and produces secondary collapses and explosive blasts.     

Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

Depositional environment of a Kimmeridgian carbonate 'black band' (Akkuyu Formation, south-western Turkey)

The Upper Jurassic Akkuyu Formation of south-western Turkey consists of pale-coloured carbonate facies, but includes a peculiar ‘black band’ at its base. In order to characterize the depositional environment of this level, a sedimentological, organic and inorganic geochemical study was performed on 30 samples. Three units are distinguished. Unit 1 consists of black marlstones, very rich in organic carbon (22–30%) with HI (hydrogen index) values around 600. The organic matter is an amorphous homogeneous material even at nanoscopic scale. This unit is also characterized by a very high concentration of V ( >1000 p.p.m.) and a δ¹³C_org value around ?26·5‰. Unit 2 shows thin-bedded black limestones with chert layers. Its organic content varies from 2% to 6%. As in unit 1, the HI values are uniformly high (600), and gas chromatography shows a high quantity of hopanes. The V concentration is also high (≈400 p.p.m.). The unit has a high Sr/Ca ratio, an organic carbon isotopic ratio close to ?28·7‰ and diffuse, nanoscopically amorphous organic material that is closely associated with pyrite. Unit 3 corresponds to the white limestones overlying the ‘black band’. This calcitic facies is almost devoid of organic carbon (<0·1%). The δ¹³C_org values are similar to those of unit 1 (?25·5‰ to ?26‰). Deposition occurred in a basinal environment surrounded by a carbonate platform. Primary productivity was high but not exceptional; terrestrial organic inputs were scarce. Bottom waters were probably euxinic during the deposition of unit 1. The presence of both hopane-rich and sulphur-rich nanoscopically amorphous organic matter may be the result of strong sulphate-reducing bacterial activity. Coupled with low amounts of iron, this led to intense vulcanization of the organic matter and thus good preservation. With the onset of unit 2 deposition, environmental conditions became less extreme, perhaps because of cooler water temperatures. Finally, the depositional conditions became oxic, leading to the accumulation of pale-coloured organic-poor limestones of unit 3.     

Assessment of sandstone deterioration at Ta Keo temple (Angkor): first results and future prospects

A first application of geomorphological methods to the assessment of sandstone deterioration at Angkor is presented. Damage diagnosis was carried out on the first eastern tier of the central pyramid of the 1,000 year-old Ta Keo temple. Methods combine field observations and measurements at 230 sampling points, high-resolution lasergrammetry and stereophotogrammetry on a 2-m² test zone, and SEM observations. The first results indicate that decay operates through a synergistic combination of weathering phenomena dominated by scaling and solution, and exhibits a high spatial variability. Percentages of deteriorated surfaces vary from 17.6 to 93.8%, and average stone recession values from 0.00 to 2.71 cm (minimum) and 0.34 to 5.49 cm (maximum). On the test-zone, stereophotogrammetry and 3D-mapping of the present and reconstructed initial states using lasergrammetry indicate that erosion scars up to 6 cm deep have formed since 1963. On the whole, the amount of deteriorated surfaces more than tripled between 1963 and 2008. The degree of implication of salts in stone decay remains unclear for most efflorescences are composed of calcite (CaCO3), with secondary importance of barite (BaSO4) and gypsum (CaSO₄·2H₂O). Future prospects aim to evaluate the impact on stone decay of the clearing out of the temple from the forest in the 1920s.     

Conflicting structural and geochronological data from the Ibituruna quartz-syenite (SE Brazil): Effect of protracted “hot” orogeny and slow cooling rate?

The Ibituruna quartz-syenite was emplaced as a sill in the Ribeira-Araçuaí Neoproterozoic belt (Southeastern Brazil) during the last stages of the Gondwana supercontinent amalgamation. We have measured the Anisotropy of Magnetic Susceptibility (AMS) in samples from the Ibituruna sill to unravel its magnetic fabric that is regarded as a proxy for its magmatic fabric. A large magnetic anisotropy, dominantly due to magnetite, and a consistent magnetic fabric have been determined over the entire Ibituruna massif. The magmatic foliation and lineation are strikingly parallel to the solid-state mylonitic foliation and lineation measured in the country-rock. Altogether, these observations suggest that the Ibituruna sill was emplaced during the high temperature (~ 750 °C) regional deformation and was deformed before full solidification coherently with its country-rock. Unexpectedly, geochronological data suggest a rather different conclusion. LA-ICP-MS and SHRIMP ages of zircons from the Ibituruna quartz-syenite are in the range 530–535 Ma and LA-ICP-MS ages of zircons and monazites from synkinematic leucocratic veins in the country-rocks suggest a crystallization at ~ 570–580 Ma, i.e., an HT deformation > 35My older than the emplacement of the Ibituruna quartz-syenite. Conclusions from the structural and the geochronological studies are therefore conflicting. A possible explanation arises from ⁴⁰Ar–³⁹Ar thermochronology. We have dated amphiboles from the quartz-syenite, and amphiboles and biotites from the country-rock. Together with the ages of monazites and zircons in the country-rock, ⁴⁰Ar–³⁹Ar mineral ages suggest a very low cooling rate: < 3 °C/My between 570 and ~ 500 Ma and ~ 5 °C/My between 500 and 460 Ma. Assuming a protracted regional deformation consistent over tens of My, under such stable thermal conditions the fabric and microstructure of deformed rocks may remain almost unchanged even if they underwent and recorded strain pulses separated by long periods of time. This may be a characteristic of slow cooling “hot orogens” that rocks deformed at significantly different periods during the orogeny, but under roughly unchanged temperature conditions, may display almost indiscernible microstructure and fabric.     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Sumixam maximus gen. and sp. nov., a New Pennsylvanian Species of Archaeorthoptera (Neoptera; Indiana, USA)

A new Mazon-Creek-type fossil specimen is described as Sumixam maximus gen. and sp. nov., assigned to the taxon Archaeorthoptera Béthoux and Nel, 2002 (composed of Orthoptera and some fossil relatives): it exhibits one of the diagnostic character states of this taxon, which is the fusion of CuA (emerging from M + CuA) with CuP, or a branch of it. A more precise taxonomic assignment is out of reach. The new taxon exhibits a fusion of the anterior branch of MA with RP, and a point of divergence of MA and MP located near the point of divergence of CuA (from M + CuA), but lack a branching of CuPa. This unique combination of character states is probably derived with respect to a number of contemporaneous species. Sumixam maximus gen. and sp. nov. is most likely the closest known relative of the panorthopterans, which include all the recent Archaeorthoptera.