Cryptic crustal contamination of MORB primitive melts recorded in olivine-hosted glass and mineral inclusions

The origin of compositional heterogeneities among the magmas parental to mid-ocean ridge basalts (MORB) was investigated using a single rock piece of the olivine-phyric basalt from 43°N, Mid-Atlantic Ridge (AII D11-177). The exceptional feature of this sample is presence of very primitive olivine crystals (90–91 mol% Fo) that are significantly variable in terms of CaO (0.15–0.35 wt%). A population of low-Ca olivine (0.15–0.25 wt% CaO) is also notably distinct from high-Ca olivine population in AII D11-177, and primitive MORB olivine in general, in having unusual assemblage of trapped mineral and glass inclusions. Mineral inclusions are represented by high-magnesian (Mg# 90.7–91.1 mol%) orthopyroxene and Cr-spinel, distinctly enriched in TiO₂ (up to 5 wt%, c.f. <1 wt% in common MORB spinel). Glass inclusions associated with orthopyroxene and high-Ti Cr-spinel have andesitic compositions (53–58 wt% SiO₂). Compared to the pillow-rim glass and “normal” MORB inclusions, the Si-rich glass inclusions in low-Ca olivine have strongly reduced Ca and elevated concentrations of Ti, Na, K, P, Cl, and highly incompatible trace elements. Strong variability is recorded among glass inclusions within a single olivine phenocrysts. We argue that the observed compositional anomalies are mineralogically controlled, and thus may arise from the interaction between hot MORB magmas and crystal cumulates in the oceanic crust or magma chamber.     

<Superscript>57</Superscript>Fe Mössbauer spectroscopy,X-ray single-crystal diffractometry,and electronic structure calculations on natural alexandrite

Natural alexandrite Al₂BeO₄:Cr from Malyshevo near Terem Tschanka, Sverdlovsk, Ural, Russia, has been characterized by ⁵⁷Fe Mössbauer spectroscopy, electron microprobe, X-ray single-crystal diffractometry and by electronic structure calculations in order to determine oxidation state and location of iron. The sample contains 0.3 wt% of total iron oxide. The ⁵⁷Fe Mössbauer spectrum can be resolved into three doublets. Two of them with hyperfine parameters typical for octahedrally coordinated high-spin Fe³⁺ and Fe²⁺, respectively, are assigned to iron substituting for Al in the octahedral M2-site. The third doublet is attributed to Fe³⁺ in hematite. Electronic structure calculations in the local spin density approximation are in reasonable agreement with experimental data provided that expansion and/or distortion of the coordination octahedra are presumed upon iron substitution. The calculated hyperfine parameters of Fe³⁺ are almost identical for the M1 and M2 positions, but the calculated ligand-field splitting is by far too large for high-spin Fe³⁺ on M1.     

Changing poverty distribution in Bolivia: the role of rural–urban migration and urban services

Bolivia is a country with high levels of poverty and inequality among its peoples and regions. For the nation and its urban and rural areas, trends in the social and spatial distribution of poverty (and extreme poverty) are identified from 1976 to 2003 using UBN data with minor support where appropriate from poverty lines. The main survey between 1992 and 2001 uses composite and selected UBN to track detailed poverty change for the country’s nine departments, its ten largest cities and a selection of other smaller urban and rural municipalities. Because of rising background increases in population in the various surveyed administrative units, many instances of relative reductions in poverty are accompanied by rising absolute increases. Marked spatial variations in poverty and development in the country over the last several decades are identified as the main driver for the country’s quickening pace of rural–urban migration. As a result, the paper concludes by assessing two different but closely related views. One investigation tests the notion that because more poor people have been living in Bolivia’s cities than in its rural areas since the mid to late 1990s, rapid rural–urban migration has simply shifted the locus of poverty from the countryside to the cities in a process called, the ‘urbanisation of poverty.’ A second, more challenging, investigation assesses the view that the flow of poor rural people to the better serviced urban areas of Bolivia has actually acted to alleviate national poverty levels.     

Earthquakes in Switzerland and surrounding regions during 2006

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2006. During this period, 572 earthquakes and 91 quarry blasts were detected and located in the region under consideration. Of these earthquakes, two occurred in conjunction with the construction of the new Gotthard railway tunnel and 165 were induced artificially by the stimulation of a proposed geothermal reservoir beneath the city of Basel. With 20 events with M _L ≥ 2.5, five of which were artificially induced, the seismic activity in the year 2006 was far below the average over the previous 31 years. Nevertheless, six events were felt by the public, most prominently the strongest of the induced Basel events (M _L 3.4), which caused some non-structural building damage. Noteworthy are also the two earthquakes near Cortaillod (M _L 3.2), on the shore of Lake Neuchatel, and in Val Mora (M _L 3.5), between the Engadin and Val Müstair, as well as the 42 aftershocks of the M _L 4.9 Vallorcine earthquake, between Martigny and Chamonix, of September 2005. Editorial handling: Stefan Bucher     

Limitations of real-time models for forecasting river flooding from monsoon rainfall

Very intense rainfall during the southwest and northeast monsoons causes severe river flooding in India. Some traditional techniques used for real-time forecasting of flooding involve the relationship between effective rainfall and direct surface runoff, which simplifies the complex interactions between rainfall and runoff processes. There are, however, serious problems in deducing these variables in real time, so it is highly desirable to have a real-time flood forecasting model that would directly relate the observed discharge hydrograph to the observed rainfall. The storage routing model described by Baba and Hoshi (1997), Tanaka et al. (1997), and Baba et al. (2000), and a simplified version of this model, have been used to compute observed river discharge directly from observed hourly rainfall. This method has been used to study rainfall–runoff data of the Ajay River Basin in eastern India. Five intense rainfall events of this basin were studied. Our results showed that the Nash–Sutcliffe efficiency of discharge prediction for these five events was 98.6%, 94.3%, 86.9%, 85.6%, and 67%. The hindcast for the first two events is regarded as completely satisfactory whereas for the next two events it is deemed reasonable and for the fifth it is unsatisfactory. It seems the models will yield accurate hindcast if the rainfall is uniform over the drainage basin. When the rainfall is not uniform the performance of the model is unsatisfactory. In future this problem can, in principle, be corrected by using a weighted amount if rainfall is based upon multiple rain-gauge observations over the drainage basin. This would provide some measure of the dispersion in the rainfall. The model also seems unable to simulate flooding events with multiple peaks.     

Reliability analysis of drilled shaft behavior using finite difference method and <Emphasis Type="Italic">Monte Carlo</Emphasis> simulation

Load displacement analysis of drilled shafts can be accomplished by utilizing the “t-z” method, which models soil resistance along the length and tip of the drilled shaft as a series of springs. For non-linear soil springs, the governing differential equation that describes the soil-structure interaction may be discretized into a set of algebraic equations based upon finite difference methods. This system of algebraic equations may be solved to determine the load–displacement behavior of the drilled shaft when subjected to compression or pullout. By combining the finite difference method with Monte Carlo simulation techniques, a probabilistic load–displacement analysis can be conducted. The probabilistic analysis is advantageous compared to standard factor of safety design because uncertainties with the shaft–soil interface and tip properties can be independently quantified. This paper presents a reliability analysis of drilled shaft behavior by combining the finite difference technique for analyzing non-linear load–displacement behavior with Monte Carlo simulation method. As a result we develop probabilistic relationships for drilled shaft design for both total stress (undrained) and effective stress (drained) parameters. The results are presented in the form of factor of safety or resistance factors suitable for serviceability design of drilled shafts.     

Discovering the ancient lakelet of the archeological site of pasargadae in Iran using a multi-layer technique

Pasargadae complex had been the main center of Persian history during the Achaemenian period (560–330 bc), registered as one of the world heritage sites by UNESCO in 2004. This research is an attempt to recognize the natural setting and landscape of Pasargadae during the Achaemenian period, in Iran. The method is based on the application of a multi-layer technique. Collected data by fieldwork, especially in geology, hydrology, and hydrogeology, along with the analysis of aerial data and satellite information provided the baseline information for preparing historic and environmental layers and have made it possible to recognize a natural lakelet in the this area. The results show that the geometry of lakelet together with the ancient river of Pasargadae had been the main elements that formed the character of the site; they were also important in the locations of buildings in the Pasargadae complex.     

Coupling of oceanic and continental crust during Eocene eclogite-facies metamorphism: evidence from the Monte Rosa nappe,western Alps

High precision U–Pb geochronology of rutile from quartz–carbonate–white mica–rutile veins that are hosted within eclogite and schist of the Monte Rosa nappe, western Alps, Italy, indicate that the Monte Rosa nappe was at eclogite-facies metamorphic conditions at 42.6 ± 0.6 Ma. The sample area [Indren glacier, Furgg zone; Dal Piaz (2001) Geology of the Monte Rosa massif: historical review and personal comments. SMPM] consists of eclogite boudins that are exposed inside a south-plunging overturned synform within micaceous schist. Associated with the eclogite and schist are quartz–carbonate–white mica–rutile veins that formed in tension cracks in the eclogite and along the contact between eclogite and surrounding schist. Intrusion of the veins at about 42.6 Ma occurred at eclogite-facies metamorphic conditions (480–570°C, >1.3–1.4 GPa) based on textural relations, oxygen isotope thermometry, and geothermobarometry. The timing of eclogite-facies metamorphism in the Monte Rosa nappe determined in this study is identical to that of the Gran Paradiso nappe [Meffan-Main et al. (2004) J Metamorphic Geol 22:261–281], confirming that these two units have shared the same Alpine metamorphic history. Furthermore, the Gran Paradiso and Monte Rosa nappes underwent eclogite-facies metamorphism within the same time interval as the structurally overlying Zermatt-Saas ophiolite [∼50–40 Ma; e.g., Amato et al. (1999) Earth Planet Sci Lett 171:425–438; Mayer et al. (1999) Eur Union Geosci 10:809 (abstract); Lapen et al. (2003) Earth Planet Sci Lett 215:57–72]. The nearly identical P–T–t histories of the Gran Paradiso, Monte Rosa, and Zermatt-Saas units suggest that these units shared a common Alpine tectonic and metamorphic history. The close spatial and temporal associations between high pressure (HP) ophiolite and continental crust during Alpine orogeny indicates that the HP internal basement nappes in the western Alps may have played a key role in exhumation and preservation of the ophiolitic rocks through buoyancy-driven uplift. Coupling of oceanic and continental crust may therefore be critical in preventing permanent loss of oceanic crust to the mantle.     

Roughness control on hydraulic conductivity in fractured rocks

The influence of joint roughness on the typologies of fluid flow inside fractures is well known and, thanks to experiences in the field of hydraulics, it has been studied from both a physical and mathematical point of view. Nevertheless, the formulations adopted by traditional hydraulic models are hardly applicable in the geological field, because of the difficulty encountered in the roughness parameter estimation. Normally this parameter can be estimated using the joint roughness coefficient (JRC), which considers both the asperity height and its regularity and directional trend. The main advantage in using the JRC arises from the fact that it can easily be obtained from geological-technical surveys and from comparison with the standard Barton profiles. Some relationships have been built up that allow for the estimation of the hydraulic conductivity tensor (an essential parameter for understanding water flow in fractured rock masses), not only as a function of traditional parameters like aperture, spacing, dip and dip direction, etc., but also of joint roughness, precisely expressed in terms of the roughness coefficient. These relationships have been studied initially from a theoretical point of view and then practically, through laboratory investigations.

---

Resumen  Se conoce muy bien la influencia de la rugosidad de las grietas en las tipologías del flujo de fluidos a lo interior de las fracturas y gracias a las experiencias en el campo de hidráulica ha sido posible estudiarla desde puntos de vista matemáticos y físicos. Sin embargo, las formulaciones adoptadas por los modelos hidráulicos tradicionales tienen poca aplicabilidad en el campo geológico debido a la dificultad relacionada con la estimación del parámetro de rugosidad. Normalmente este parámetro puede estimarse usando el coeficiente de rugosidad de grieta (JRC) el cual considera tanto la altura de la aspereza como su regularidad y tendencia direccional. La principal ventaja de utilizar el JRC se deriva del hecho que puede obtenerse fácilmente de levantamientos técnico-geológicos y de la comparación con los perfiles Standard Barton. Se han construido algunas relaciones que permiten la estimación del tensor de conductividad hidráulica (un parámetro esencial para el entendimiento del flujo de agua en masas de roca fracturadas), no solo en función de parámetros tradicionales como apertura, espaciado, buzamiento y dirección de buzamiento, etc., sino también en función de la rugosidad de la grieta estimada con precisión en términos del coeficiente de rugosidad. Estas relaciones se han estudiado inicialmente desde un punto de vista teórico y luego de modo práctico a través de investigaciones de laboratorio.

---

Résumé  L’influence de la rugosité des joints sur les types d’écoulement de fluide dans les fractures est bien connue et a été étudiée aussi bien du point de vue physique que mathématique grace à des expériences menées dans le domaine de l’hydraulique. Cependant les formulations adoptées dans les modèles hydrauliques traditionnels sont difficilement applicables dans le domaine de la géologie à cause de la difficulté rencontrée pour estimer la rugosité. Ce paramètre peut normalement être apprécié grace au coefficient de rugosité du joint (JRC), lequel prend en compte à la fois la hauteur de l’aspérité ainsi que sa régularité et sa direction. Le principal avantage dans l’utilisation du JRC réside dans le fait qu’il peut facilement être obtenu à partir d’études techniques-géologiques et par comparaison avec la classification de Barton. Des relations qui permettent une estimation du tenseur de conductivité hydraulique (un paramètre essentiel pour comprendre l’écoulement de l’eau dans les masses rocheuses fracturées) ont été élaborées, pas seulement en fonction de paramètres traditionnels tels que l’ouverture, l’espacement, l’inclinaison et la direction d’inclinaison, etc , mais aussi en prenant en compte la rugosité des joints à travers le coefficient de rugosité. Ces relations ont initialement été étudiées d’un point de vue théorique puis expérimentalement à travers des recherches en laboratoire.

     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.