Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography

We present an adaptive unstructured triangular grid finite element approach for effectively simulating plane-wave diffusive electromagnetic fields in 2-D conductivity structures.
The most striking advantage of irregular grids is their potential to incorporate arbitrary geometries including surface and seafloor topography. Adaptive mesh refinement strategies using an a posteriori error estimator yield most efficient numerical solutions since meshes are only refined where required.
We demonstrate the robustness of this approach by comparison with analytical solutions and previously published numerical simulations. Maximum errors may systematically be reduced to, for example, 0.8 per cent for the apparent resistivity and 0.2° in the phase.
An additional accuracy study of the thickness of the air layer in E-polarization suggests to keep a minimum thickness depending on lateral conductivity contrasts within the earth.
Furthermore, we point out the new quality and flexibility of our simulation technique by addressing two marine magnetotelluric applications. In the first case, we discuss topographic effects associated with a synthetic sinusoidal sea bottom model and in the second case, we show a close-to-reality scenario using real bathymetry data from the East Pacific Rise at 17°S.     

87Sr/86Sr ratio: a natural tracer to monitor groundwater flow paths during artificial recharge in the Bangkok area, Thailand

The aquifer system in the Thon Buri sedimentary basin below the deltaic flood plain of the Chao Phraya River, central Thailand, has been exploited for public water supply for the capital Bangkok since the early 1920s. Groundwater withdrawal, currently 1.4 million m³/d, has resulted in a maximum decline in hydraulic head of up to 40 m. This has induced land subsidence of as much as 1.7 m (1940–1992) in the eastern suburbs of the metropolis. Artificial injection of purified water within an area-wide network of recharge wells could constitute a remedy to slow the water level depression within the sedimentary basin, and thus the subsidence. This requires a prior shutdown of water withdrawal. The flow paths of the injected water can be traced by changes in the ⁸⁷Sr/⁸⁶Sr ratio of the groundwater and injected water mixture within the three main aquifers in the basin that are used for public supply. The ratios, monitored at five monitoring stations within the cone of depression, have been constant over 3 years. Injection of the calculated cone volume of 5.2?×?10⁹ m³ would take at least 10 years, depending on the injection pressure and the number and position of wells.     

Rare earth element fractionation in magmatic Ca-rich garnets

Igneous garnets have the potential to strongly fractionate rare earth elements (REE). Yet informations on partition coefficients are very scant, and criteria for distinguishing between hydrothermal and magmatic garnets are ambiguous. To fill this gap, we present trace element and isotopic data for two types of Ca-rich garnets from phonolites (Mt. Somma-Vesuvius). Both Ca-garnet populations are different in their style and dynamics of fractionation: one population is progressively strongly depleted in HREE from core to rim, reflecting REE fractionation in the host phonolite via earlier-crystallized garnets. Such examples for extreme changes in HREE in garnets are only known for hydrothermal grandites by REE-bearing fluids. The second garnet population is homogeneous and formed in a closed system. Near-flat patterns between Sm and Lu confirm experimental data indicating lower D(Sm)/D(Lu) for Ca-rich garnets than for e.g. pyrope-rich garnets. It follows: D ^Grt/PhMelt for La = 0.5, Sm = 48 and Yb = 110.     

Niobium-tantalum Fractionation During Slab Subduction:Implications for the Formation of Continental Crust

Compared with the oceanic crust, knowledge about the formation of the continental crust (CC) is relatively poor. Although melting of subducted slabs in the early history of the Earth has been considered as the major way that shaped the chemical characteristics of the CC by most geologists, as the CC shares many characteristics with modern adakites, some geologists argued that Archean TTG was formed in the same way as modern arcs rather than slab melting, whereas others proposed that melting at the bottom of the thickened oceanic crust was more important. Recently,the debate is mainly focused on the unique subchondritic Nb/Ta value of the CC, and particularly, how Nb and Ta fractionated from each other and consequently how, in detail, the CC was built.     

A new relation between the central spectral solar H I Lyman α irradiance and the line irradiance measured by SUMER/SOHO during the cycle 23

The spectral irradiance at the center of the solar H I Lyman α (, referred to as Lyα in this paper) line profile is the main excitation source responsible for the atomic hydrogen resonant scattering of cool material in our Solar System. It has therefore to be known with the best possible accuracy in order to model the various Lyα emissions taking place in planetary, cometary, and interplanetary environments. Since the only permanently monitored solar irradiance is the total one (i.e. integrated over the whole Lyα line profile), Vidal-Madjar [1975. Evolution of the solar Lyman alpha flux during four consecutive years. Solar Phys. 40, 69-86] using Orbiting Solar Observatory 5 (OSO-5) satellite Lyα data, established a semi-empirical formula allowing him to deduce the central spectral Lyα irradiance from the total one. This relation has been extensively used for three decades. But, at the low altitude of the OSO-5 orbit, the central part of the solar line profile was deeply absorbed by a large column of exospheric atomic hydrogen. Consequently, the spectral irradiance at the center of the line was obtained by a complex procedure confronting the observations with simulations of both the geocoronal absorption and the self-reversed shape of the solar Lyα profile. The SUMER spectrometer onboard SOHO positioned well outside the hydrogen geocorona, provided full-Sun Lyα profiles, not affected by such an absorption [Lemaire et al., 1998. Solar H I Lyman α full disk profile obtained with the SUMER/SOHO spectrometer. Astron. Astrophys. 334, 1095-1098; 2002. Variation of the full Sun Hydrogen Lyman α and β profiles with the activity cycle. Proc. SOHO 11 Symposium, ESA SP-508, 219-222; 2004. Variation of the full Sun Hydrogen Lyman profiles through solar cycle 23. COSPAR 2004 Meeting], making it—for the first time—possible to measure the spectral and total Lyα solar irradiances directly and simultaneously. A new relation between these two quantities is derived in an expression that is formally similar to the previous one, but with significantly different parameters. After having discussed the potential causes for such differences, it is suggested that the new relation should replace the old one for any future modeling of the numerous Lyα absorptions and emissions observed in the Solar System.     

Significance of strain localization and fracturing in relation to hydrothermal mineralization at Mount Isa, Australia

The rheology of layered meta-sedimentary rocks, and their orientation and position relative to major fault systems were the key controls on Proterozoic hydrothermal copper mineralization at Mount Isa, Australia. Compositional layering in the host rock partitioned mechanical behavior and strain, leading to selective permeability generation and the focusing of fluid flow. Shale layers preferentially failed by plastic shearing, whereas meta-siltstones remained elastic or failed in tension depending on magnitude of deformation and fluid pressure. Numerical simulations support the hypothesis that the orientation of layering and the proximity to major fault systems controlled fracturing and permeability increase in the Urquhart shale. The dilating shale provided a pathway for an upward-flowing, reduced basement fluid, from which quartz was precipitated during cooling. During a later event, the reactivation of steep structures provided access to surface derived oxidized metal-bearing brine, causing the precipitation of dolomite followed by chalcopyrite ore in the brecciated silicified shale.     

SIZE-DISTRIBUTIONS OF CHONDRULE TYPES IN THE INMAN AND ALLAN HILLS A77011 L3 CHONDRITES

A petrographc study of 9 thin sections of Inman (L3) and 18 thin sections of ALHA77011 (L3) served to determine the size-distributions of different chondrule textural types. Inman chondrules are significantly larger than those in ALHA77011, but in each chondrite, there is no statistically significant difference between the size-distributions of barred olivine and radial pyroxene plus cryptocrystalline chondrules. In ALHA77011, barred olivine chondrules outnumber radial pyroxene plus cryptocrystalline chondrules, whereas in Inman, the reverse is true. Because compound and cratered chondrules were formed by the collision of similarly-sized objects, the dustball precursors of chondrules must have been size-sorted prior to chondrule formation. The region of dustball size-sorting in the solar nebula must have been very large, similarly affecting the physically-separated precursors of different chondrule types. Size-sorting was probably accomplished by aerodynamic particle-gas interactions. Zones of dustball melting (i.e., chondrule formation) were relatively small, generally affecting only dustballs of one compositional type and relatively uniform size. Different chondrule types were then mixed together in somewhat variable ratios. Within the region where chondrites of a particular compositional group agglomerated, there were sub-reservoirs that contained (roughly) uniformly large or uniformly small chondrules with different mixtures of textural types.