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     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

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.     

Biological and Geochemical Record of Anthropogenic Impacts in Recent Sediments from Lake Pupuke, Auckland City, New Zealand

Frozen sediment cores from Lake Pupuke in Auckland City, New Zealand, contain a high resolution decadal to annual scale record of changing lake paleoenvironments and geochemistry that reflects changing landuse and hydrology in the catchment over the past c. 190 years. A reliable chronology is available from AMS ¹⁴C and ²¹⁰Pb dating of the sediments, with the timing of the older part of the record supported by the first appearance of pollen of introduced Pinus and Plantago lanceolata associated with European settlement of Auckland in the early 1840s. Diatom stratigraphy, sediment elemental and carbon isotope geochemistry reflect changes in sediment sources and lake conditions commensurate with European development of the Lake Pupuke catchment, in particular enhanced algal productivity controlled by the influx of nutrients after c. 1920 AD. Attempts to prevent nuisance algal blooms in 1933, 1934 and 1939 using CuSO₄ addition produced Cu spikes in the sediment that allowed confirmation of the accuracy of the ²¹⁰Pb chronology. Hence, the elemental and isotopic composition of the Lake Pupuke sediments reflect the timing of significant anthropogenic activities, rather than climatic variations, that have occurred within the watershed over the past c. 190 years. The comparison of records of land use change in the catchment with the multi-proxy record of changes in the sediments usually allowed unambiguous identification of the signatures of change and their causes.