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.     

Lander radioscience for obtaining the rotation and orientation of Mars

The paper presents the concept, the objectives, the approach used, and the expected performances and accuracies of a radioscience experiment based on a radio link between the Earth and the surface of Mars. This experiment involves radioscience equipment installed on a lander at the surface of Mars. The experiment with the generic name lander radioscience (LaRa) consists of an X-band transponder that has been designed to obtain, over at least one Martian year, two-way Doppler measurements from the radio link between the ExoMars lander and the Earth (ExoMars is an ESA mission to Mars due to launch in 2013). These Doppler measurements will be used to obtain Mars’ orientation in space and rotation (precession and nutations, and length-of-day variations). More specifically, the relative position of the lander on the surface of Mars with respect to the Earth ground stations allows reconstructing Mars’ time varying orientation and rotation in space.Precession will be determined with an accuracy better by a factor of 4 (better than the 0.1% level) with respect to the present-day accuracy after only a few months at the Martian surface. This precession determination will, in turn, improve the determination of the moment of inertia of the whole planet (mantle plus core) and the radius of the core: for a specific interior composition or even for a range of possible compositions, the core radius is expected to be determined with a precision decreasing to a few tens of kilometers.A fairly precise measurement of variations in the orientation of Mars’ spin axis will enable, in addition to the determination of the moment of inertia of the core, an even better determination of the size of the core via the core resonance in the nutation amplitudes. When the core is liquid, the free core nutation (FCN) resonance induces a change in the nutation amplitudes, with respect to their values for a solid planet, at the percent level in the large semi-annual prograde nutation amplitude and even more (a few percent, a few tens of percent or more, depending on the FCN period) for the retrograde ter-annual nutation amplitude. The resonance amplification depends on the size, moment of inertia, and flattening of the core. For a large core, the amplification can be very large, ensuring the detection of the FCN, and determination of the core moment of inertia.The measurement of variations in Mars’ rotation also determines variations of the angular momentum due to seasonal mass transfer between the atmosphere and ice caps. Observations even for a short period of 180 days at the surface of Mars will decrease the uncertainty by a factor of two with respect to the present knowledge of these quantities (at the 10% level).The ultimate objectives of the proposed experiment are to obtain information on Mars’ interior and on the sublimation/condensation of CO₂ in Mars’ atmosphere. Improved knowledge of the interior will help us to better understand the formation and evolution of Mars. Improved knowledge of the CO₂ sublimation/condensation cycle will enable better understanding of the circulation and dynamics of Mars’ atmosphere.     

The catastrophe structure of thermohaline convection in a two-dimensional fluid model and a comparison with low-order box models

Abstract

We impose a surface forcing on the 2D, Boussinesq, thermohaline equations in a rectangular domain, in the form of equatorially symmetric cosine distributions of salinity flux and temperature. This system may be seen as an idealization of the ocean thermohaline circulation on the global scale over intervals of centuries or millenia. Multiple steady states are found numerically. They reflect the competition between the opposite signs of the temperature and salinity-driven equatorially symmetric circulations. There are also pole-to-pole, equatorially asymmetric circulations. In the control space of the temperature and salinity-flux forcing amplitudes, these equilibria form two cusp catastrophes, and transitions between stable equilibria occur through several distinct bifurcations. These catastrophes can be reproduced in simple box models connecting stirred reservoirs through capillary pipes. This steady-state analysis may provide a framework for a better understanding of climatic transitions between different stable regimes of the ocean-atmosphere system.     

Sulfur isotope measurement of sulfate and sulfide by high-resolution MC-ICP-MS

We have developed a technique for the accurate and precise determination of ³⁴S/³²S isotope ratios (δ³⁴S) in sulfur-bearing minerals using solution and laser ablation multiple-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). We have examined and determined rigorous corrections for analytical difficulties such as instrumental mass bias, unresolved isobaric interferences, blanks, and laser ablation- and matrix-induced isotopic fractionation. Use of high resolution sector-field mass spectrometry removes major isobaric interferences from O₂⁺. Standard-sample bracketing is used to correct for the instrumental mass bias of unknown samples. Background on sulfur masses arising from memory effects and residual oxygen-tailing are typically minor (< 0.2‰, within analytical error), and are mathematically removed by on-peak zero subtraction and by bracketing of samples with standards determined at the same signal intensity (within 20%). Matrix effects are significant (up to 0.7‰) for matrix compositions relevant to many natural sulfur-bearing minerals. For solution analysis, sulfur isotope compositions are best determined using purified (matrix-clean) sulfur standards and sample solutions using the chemical purification protocol we present. For in situ analysis, where the complex matrix cannot be removed prior to analysis, appropriately matrix-matching standards and samples removes matrix artifacts and yields sulfur isotope ratios consistent with conventional techniques using matrix-clean analytes. Our method enables solid samples to be calibrated against aqueous standards; a consideration that is important when certified, isotopically-homogeneous and appropriately matrix-matched solid standards do not exist. Further, bulk and in situ analyses can be performed interchangeably in a single analytical session because the instrumental setup is identical for both. We validated the robustness of our analytical method through multiple isotope analyses of a range of reference materials and have compared these with isotope ratios determined using independent techniques. Long-term reproducibility of S isotope compositions is typically 0.20‰ and 0.45‰ (2σ) for solution and laser analysis, respectively. Our method affords the opportunity to make accurate and relatively precise S isotope measurement for a wide range of sulfur-bearing materials, and is particularly appropriate for geologic samples with complex matrix and for which high-resolution in situ analysis is critical.     

A thermodynamic model for Fe–Mg dioctahedral K white micas using data from phase-equilibrium experiments and natural pelitic assemblages

The purpose of this study is to derive a solid-solution model for potassic white micas (KWM) encountered in rocks of various bulk compositions, over a wide range of P-T conditions. A compilation of phengite compositions lead us to propose a seven-thermodynamic-component (muscovite, Fe²⁺-Al-celadonite, Mg-Al-celadonite, annite, phlogopite, pyrophyllite and paragonite) ionic solid-solution model which accounts for the Tschermak, Fe-Mg, di/trioctahedral, pyrophyllitic and paragonitic substitutions observed in nature. A four-site mixing model with symmetric Margules parameters to model the Tschermak substitutions, asymmetric Margules parameters to model the other substitutions, and ideal intersite interaction has been adopted. In contrast to previous models, the relevant thermodynamic data and solid-solution properties are calibrated with independent sets of published experiments conducted for the KMASH, KFASH, KFMASH, and KNASH systems, as well as about 200 natural data involving KWM assemblages. The constraints span a wide range of pressure and temperature conditions (150 to 750 °C, 0.5 to 30 kbar), so that our model does not need to be extrapolated outside the calibration range to be used for P-T thermobarometric purposes. The calculated thermodynamic data are interconsistent with the TWQ thermodynamic database and solid-solution models, including that recently published for chlorites.     

Zircon geochronology and geochemistry to constrain the youngest eruption events and magma evolution of the Mid-Miocene ignimbrite flare-up in the Pannonian Basin,eastern central Europe

A silicic ignimbrite flare-up episode occurred in the Pannonian Basin during the Miocene, coeval with the syn-extensional period in the region. It produced important correlation horizons in the regional stratigraphy; however, they lacked precise and accurate geochronology. Here, we used U–Pb (LA-ICP-MS and ID-TIMS) and (U–Th)/He dating of zircons to determine the eruption ages of the youngest stage of this volcanic activity and constrain the longevity of the magma storage in crustal reservoirs. Reliability of the U–Pb data is supported by (U–Th)/He zircon dating and magnetostratigraphic constraints. We distinguish four eruptive phases from 15.9 ± 0.3 to 14.1 ± 0.3 Ma, each of which possibly includes multiple eruptive events. Among these, at least two large volume eruptions (>10 km³) occurred at 14.8 ± 0.3 Ma (Demjén ignimbrite) and 14.1 ± 0.3 Ma (Harsány ignimbrite). The in situ U–Pb zircon dating shows wide age ranges (up to 700 kyr) in most of the crystal-poor pyroclastic units, containing few to no xenocrysts, which implies efficient recycling of antecrysts. We propose that long-lived silicic magma reservoirs, mostly kept as high-crystallinity mushes, have existed in the Pannonian Basin during the 16–14 Ma period. Small but significant differences in zircon, bulk rock and glass shard composition among units suggest the presence of spatially separated reservoirs, sometimes existing contemporaneously. Our results also better constrain the time frame of the main tectonic events that occurred in the Northern Pannonian Basin: We refined the upper temporal boundary (15 Ma) of the youngest counterclockwise block rotation and the beginning of a new deformation phase, which structurally characterized the onset of the youngest volcanic and sedimentary phase.