Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

Fractional crystallization models and B–Be–Li systematics at Mt Somma-Vesuvius volcano (Southern Italy)

A compilation of B–Be–Li data on rocks that cover the entire eruptive history of Somma-Vesuvius is presented and interpreted in the light of evolution models for the Somma-Vesuvius rocks. Using major and trace element data, fractional crystalllization models are presented for different geochemical units. These data were used to constrain the source mineralogy of the Somma-Vesuvius rocks (ol-opx-cpx-gar-amp of 0.4-0.3-0.1-0.1-0.1), the amount of sediment added (5–10%) and the melt fraction from batch partial melting computations (0.05–0.1). From the B–Li data it is inferred that the main process responsible for the B isotopic signature is sediment recycling. However, the B–Li data show a major variation in Li abundances respect to B which is explained with Li dehydration before the fluid enriched the mantle wedge that produced the arc magmas. The Somma-Vesuvius B isotope composition is intermediate between that of the Campi Flegrei and the broad field of the Eolian Island arc. A low Be isotopes in the recent volcanic rocks can be explained as: (a) the top 1–22 m of the incoming sediment is accreted, (b) large amounts of sediment erosion, (c) a slow rate of subduction which have provoked a long magmatic history for the Vesuvius magma, (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Italy.     

A dynamic box model of bioactive elements in the southern Taiwan Strait

A dynamic box model was applied to study the characteristics of biogeochemical cycling of PO_4-P,NO_3-N,AOU,POC and PON in the southern Taiwan Strait region based on field data of the"Minnan Taiwan Bank Fishing Ground Upwelling Ecosystem Study" during the period of Dec.1987-Nov.1988.According to the unique hydrological and topographical features of the region,six boxesand three layers were considered in the model.The variation rates and fluxes of elements induced by hor-izontal current,upwelling,by diffusion,sinking of particles and biogeochemical processes were estimatedrespectively.Results further confirmed that upwellings had important effects in this region.Thenearshore upwelling areas had net input fluxes of nutrients brought by upwelling water,also had high de-pletion rates of nutrients and production rates of particulate organic matter and dissolved oxygen.Theabnormal net production of nutrients in the middle layer(10-30m) indicated the important role of bacte-ria in this high production region.Th     

ANALYSIS OF A CASE OF EAST ASIA COASTAL CYCLOGENESIS*

The evolution of an explosive cyclone off the East Asia coast in March 1979 is described.A shortwave trough in the southern branch of upper-level westerlies initiated the incipient cyclone.Later,a polar trough in the north amplified and became in phase with the southern shortwave to form a major trough.This major trough was responsible for the rapid intensification of the surface cyclone.In the early development stage,warm and moist air was transported northward to the developing area by a strong low-level jet.The ageostrophic wind associated with the low-level jet contributed to the frontogenesis,creating a favorable low-level environment for the rapid deepening.A low-level positive potential vorticity anomaly was created prior to the onset of rapid deepening.It was a result of frontal cloud condensation.The cyclone intensified rapidly when stratospheric air with high potential vorticity penetrated to the mid-troposphere.The rapid deepening took place at a location under the left-exit region of an amplifying jet streak behind the major trough and the right-entrance region of another anticyclonically-curved subtropical jet streak in a quasi-stationary ridge overJapan.Due to the blocking effect of the Tibetan Plateau,two shortwave disturbances were observed in the upper-level westerlies on the north and south sides of the Plateau.The southern disturbance initiated the incipient surface cyclone,while the amplifying northern disturbance was responsible for the rapid deepening.Thus,the evolution of the explosive cyclone in this case can be regarded as consecutive Petterssen's "type-B" cyclogenesis in two separate stages.     

GGR Critical Review of Analytical Developments in 2004–2005

In 2005 Geostandards and Geoanalytical Research embarked upon a new initiative for its readers. Key researchers in various fields of geoanalytical technique development and their application were identified and invited to provide reviews pertinent to their expertise. As noted in the first of these publications "…instead of revisiting the historical context or decades of development in each analytical technique, the goal here has been to capture a snapshot of "hot topics" across a range of fields as represented in the… literature" (Hergt et al . 2005). Rather than prepare an annual review, a decision was taken earlier this year to provide a biennial summary of progress and accomplishments, in this case for the years 2004–2005. The principal techniques employed in Earth and environmental sciences are covered here, and include laser ablation and multicollector ICP-MS, ICP-AES, thermal ionisation and secondary ion mass spectrometry, as well as neutron activation analysis, X-ray fluorescence and atomic absorption spectrometry. A comprehensive review of the development of reference materials, often essential to these techniques, is also provided. The contributions assembled serve both to keep readers informed of advances they may be unfamiliar with, but also as a means of showcasing examples of the breadth and depth of work being conducted in these fields.     

Reference Materials in Geoanalytical Research -Review for 2004 and 2005

This review gives an overview of the use and development of reference materials of geochemical and environmental interest in the literature of the years 2004 and 2005. In these years the performance of existing methods has been improved and new geochemical applications using new techniques have been developed. Accordingly, there was an increasing need for new reference materials, especially for in situ microanalysis and for precise stable isotope measurements. In addition, there was a notable trend for further characterisation of existing reference materials, mainly for the platinum-group elements. This review focuses on five topics: reference materials for platinum-group elements, reference glasses for in situ microanalysis, zircon reference materials, isotopic reference materials, and the development and certification of reference materials.     

长白山锦北雨养泥炭剖面元素富集规律分析

为了揭示长白山锦北泥炭沼泽区泥炭沉积环境和重金属元素的历史变化及其对环境的指示意义,利用荷兰W ardenaar泥炭采样器采取泥炭样品并对样品进行高分辨率切割。用原子吸收分光光度计和等离子体发射光谱仪测定不同深度剖面Ca,Mg,Na,K,A l,Ti,Cu,Zn,Fe,Mn,Pb,Cr,N i和Co的浓度。研究分析了Ca,Mg,Na,K,A l,Ti,Mn,Zn,Cu,Pb,Cr,N i和Co在泥炭剖面的分布趋势并对其进行了相关分析和富集分析,得出Mn,Zn,Cu,Pb,Cr,N i和Co呈显著相关且在泥炭表层明显累积,揭示出泥炭表层重金属元素的大气沉降来源,表明近年来人类重金属消耗量及排放量的增加是大气沉降重金属元素的主要释放源,而且环境中重金属污染有加重的趋势。     

The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

Monazite behaviour and age significance in poly-metamorphic high-grade terrains: A case study from the western Musgrave Block, central Australia

Integrated, in situ textural, chemical and electron microprobe age analysis of monazite grains in a migmatitic metapelitic gneiss from the western Musgrave Block, central Australia has identified evidence for multiple events of growth and recrystallisation during poly-metamorphism in the Mesoproterozoic. Garnet + sillimanite-bearing metapelite underwent partial melting and segregation to palaeosome and leucosome during metamorphism between 1330 and 1296 Ma, with monazite grains in leucosome recording crystallisation at 1300 Ma. Monazite breakdown during melting is inferred to have occurred in the palaeosome. During a subsequent granulite facies event at 1200 Ma, deformation and metamorphism of leucosome and palaeosome resulted in partial disturbance of ages and potential minor growth on 1300 Ma monazite in leucosome. Growth of new, high-Y (+HREE) monazite in palaeosome domains occurred during garnet breakdown in the presence of sillimanite to cordierite and spinel, as a result of post-peak isothermal decompression. Diffusive enrichment of resorbed garnet rims in Y + HREE suggests garnet breakdown occurred slower than volume diffusion of REE. Monazite in both palaeosome and leucosome were subsequently partially to penetratively recrystallised during a retrogression event that is suggested to have occurred at 1150–1130 Ma. The intensity of recrystallisation and disturbance of ages appears linked to proximity to retrogressed garnet porphyroblasts and their occurrence in the relatively reactive or ‘fertile’ local environments provided by the palaeosome/mesosome volumes, which caused localised changes in retrogressive fluids towards compositions more aggressive to monazite. Like reaction textures, it is apparent that domainal equilibrium and reaction may control or at least strongly influence monazite REE and U–Th–Pb chemistry and hence ages.