Application of ICP-MS trace element analysis in study of ancient Chinese ceramics

Thirty-nine trace elements of the Song-Yuan period (960-1368 AD) porcelain bodies from Cizhou, Jizhou and Longquanwu kilns were analyzed with ICP-MS, a technique rarely used in Chinese archaeometry, to investigate its potential application in such studies. Trace element compositions clearly reflect the distinctive raw materials and their mineralogy at the three kilns and allow their products to be distinguished. Significant chemical variations are also observed between Yuan and Song-Jing dynasties samples from Cizhou as well as fine and coarse porcelain bodies from Longquanwu. In Cizhou, porcelains of better quality which imitate the famous Ding kiln have trace element features distinctive from ordinary Cizhou products, that indicates geochemically distinctive raw materials were used and which possibly also underwent extra refining prior to use. The distinct trace element features of different kilns and the various types of porcelains from an individual kiln can be interpreted from a geochemical perspective. ICP-MS can provide a large amount of valuable information about ancient Chinese ceramics as it is capable of analyzing >40 elements with a typical of precision < 2%.

ICP-MS Trace Element Analysis of Song Dynasty Porcelains from Ding, Jiexiu and Guantai kilns, North China

Bodies of Ding kiln white porcelains and their imitations from Guantai and Jiexiu kilns of the Chinese Song dynasty (960-1279 AD) were analysed for 40 trace elements by inductively coupled plasma mass spectrometry (ICP-MS). Numerous trace element compositions and ratios allow these visually similar products to be distinguished, and a Ding-style shard of uncertain origin is identified as a likely genuine Ding product. In Jiexiu kiln, Ding-style products have trace element features distinctive from blackwares of an inferior quality intended for the lower end market. Based on geochemical behaviour of these trace elements, we propose that geochemically distinctive raw materials were used for Ding-style products of a higher quality, which possibly also underwent purification by levigation prior to use. Capable of analysing over 40 elements with a typical long term precision of < 2%, this high precision ICP-MS method proves to be very powerful for grouping and characterising archaeological ceramics. Combined with geochemical interpretation, it can provide insights into the raw materials and techniques used by ancient potters. (C) 2004 Elsevier Ltd. All rights reserved.

Characterisation of Chinese Tang sancai from Gongxian and Yaozhou kilns using ICP-MS trace element and TIMS Sr–Nd isotopic analysis

Tang sancai is one of the most important types of Chinese ceramics. To determine the provenance of Tang sancai is important to study ancient trade and other issues. In this paper we compare ICP-MS trace elements and TIMS Sr-Nd isotopes of visibly similar Tang sancai from two major production centres Gongxian and Yaozhou. The variation in contents/ratios of many of > 40 trace elements is small for samples from Gongxian, yet is considerably bigger for that from Yaozhou. However, the variation in Sr-87/Sr-86 and Nd-143/Nd-144 isotopic ratios is very small for samples from both places. Gongxian and Yaozhou samples have distinctive Sr-Nd isotopic and trace element features despite their similarity in major elements, and these analysis data can be interpreted with geochemistry, indicating that Sr and Nd isotopes have great potential in ceramic provenance studies. The distinct characterisation of these samples provides valuable criteria for identifying provenance of Tang sancai of uncertain origin. Two modern fakes are also analysed, and they can as well be distinguished from antique Tang sancai using above criteria. (c) 2005 Published by Elsevier Ltd.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element analysis of mechanical properties in discontinuously reinforced metal matrix composites with ultrafine microstructure

This investigation focused on the finite element analyses of elastic and plastic properties of aluminium/alumina composite materials with ultrafine microstructure. The commonly used unit cell model was used to predict the elastic properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the plastic properties of the composites and the associated strengthening mechanism within the metal matrix material were investigated. The grain size of the matrix material was found to be an important factor influencing the mechanical properties of the composites studied. (C) 1997 Elsevier Science S.A.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.

Geological and trace element evidence for a marine sedimentary environment of deposition and biogenicity of 3.45 Ga stromatolitic carbonates in the Pilbara Craton, and support for a reducing Archaean ocean

Bedded carbonate rocks from the 3.45 Ga Warrawoona Group, Pilbara Craton, contain structures that have been regarded either as the oldest known stromatolites or as abiotic hydrothermal deposits. We present new field and petrological observations and high-precision REE + Y data from the carbonates in order to test the origin of the deposits. Trace element geochemistry from a number of laminated stromatolitic dolomite samples of the c. 3.40 Ga Strelley Pool Chert conclusively shows that they precipitated from anoxic seawater, probably in a very shallow environment consistent with previous sedimentological observations. Edge-wise conglomerates in troughs between stromatolites and widespread cross-stratification provide additional evidence of stromatolite construction, at least partly, from layers of particulate sediment, rather than solely from rigid crusts. Accumulation of particulate sediment on steep stromatolite sides in a high-energy environment suggests organic binding of the surface. Relative and absolute REE + Y contents are exactly comparable with Late Archaean microbial carbonates of widely agreed biological origin. Ankerite from a unit of bedded ankerite–chert couplets from near the top of the stratigraphically older (3.49 Ga) Dresser Formation, which immediately underlies wrinkly stromatolites with small, broad, low-amplitude domes, also precipitated from anoxic seawater. The REE + Y data of carbonates from the Strelley Pool Chert and Dresser Formation contrast strongly with those from siderite layers in a jasper–siderite–Fe-chlorite banded iron-formation from the base of the Panorama Formation (3.45 Ga), which is clearly hydrothermal in origin. The geochemical results, together with sedimentological data, strongly support: (1) deposition of Dresser Formation and Strelley Pool Chert carbonates from Archaean seawater, in part as particulate carbonate sediment; (2) biogenicity of the stromatolitic carbonates; (3) a reducing Archaean atmosphere; (4) ongoing extensive terrestrial erosion prior to ∼3.45 Ga.

Characterisation of early Archaean chemical sediments by trace element signatures

Rare earth element (REE) plus yttrium (Y) patterns of modem seawater have characteristic features that can be used as chemical fingerprints. Reliable proxies for marine REE + Y chemistry have been demonstrated from a large geological time span, including Archaean banded iron formation (BIF), stromatolitic limestone, Phanerozoic reef carbonate and Holocene microbialite. Here we present new REE + Y data for two distinct suites of early Archaean (ca. 3.7-3.8 Ga) metamorphosed rocks from southern West Greenland, whose interrelationships, if any, have been much debated in recent literature. The first suite comprises mangetite-quartz BIF, magnetite-carbonate BIF and banded magnetite-rich quartz rock, mostly from the Isua Greenstone Belt (IGB). The REE + Y patterns, particularly diagnostic anomalies (Ce/Ce*, Pr/Pr*), are closely related to those of published seawater proxies. The second suite includes banded quartz-pyroxene-amphibole +/- garnet rocks with minor magnetite from the so-called Akilia Association enclaves (in early Archaean granitoid gneisses) of the coastal region, some 150 km southwest of the IGB. Rocks of this type from one much publicised and highly debated locality (the island of Akilia) have been identified by some workers [Nature 384 (1996) 55; Geochim. Cosmochim. Acta 61 (1997) 2475] as BIF-facies, and their C-13-depleted signature in trace graphite interpreted as a proxy for earliest life on Earth. However, REE + Y patterns of the Akilia Association suite (except for one probably genuine magnetite-rich BIF from Ugpik) are inconsistent with a seawater origin. We agree with published geological and geochemical (including REE) work [Science 296 (2002) 1448] that most of the analysed Akilia rocks are not chemical sediments, and that C-isotopes in such rocks therefore cannot be used as biological proxies. Application of the REE + Y discriminant for the above two rock suites has been facilitated in this study by the use of MC-ICP technique which yields a more complete and precise REE + Y spectrum than was available in many previous studies. (C) 2004 Elsevier B.V. All rights reserved.

Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Potential of Sr isotopic analysis in ceramic provenance studies: Characterisation of Chinese stonewares

We compare the trace element and Sr isotopic compositions of stoneware bodies made in Yaozhou and Jizhou to characterise these Chinese archaeological ceramics and examine the potential of Sr isotopes in provenance studies. Element concentrations determined by ICP-MS achieve distinct characterisation for Jizhou samples due to their restricted variation, yet had limited success with Yaozhou wares because of their large variability. In contrast, Sr-87/Sr-86 ratios in Yaozhou samples have a very small variation and are all significantly lower than those of Jizhou samples, which show a large variation and cannot be well characterised with Sr isotopes. Geochemical interpretation reveals that Sr-87/Sr-86 ratios will have greater potential to characterise ceramics made of low Rb/Sr materials such as kaolin clay, yet will show larger variations in ceramics made of high Rb/Sr materials such as porcelain stone. (c) 2005 Elsevier B.V. All rights reserved.

A trace element study of siderite-jasper banded iron formation in the 3.45 Ga Warrawoona Group, Pilbara Craton - Formation from hydrothermal fluids and shallow seawater

Shale-normalised rare earth element and yttrium (REE + Y) patterns for siderite-jasper couples in a banded iron formation of the 3.45 Ga Panorama Formation, Warrawoona Group, eastern Pilbara Craton, display distinct positive Y and Eu anomalies and weak positive La and Gd anomalies, combined with depleted light REE relative to middle and heavy REE. Ambient seawater and hydrothermal fluids are identified as major sources of REE + Y for the BIF. In the case of siderites, strong correlations between incompatible trace elements and trace element ratios diagnostic of seawater indicate variable input from a terrigenous source (e.g. volcanic ash). We propose a volcanic caldera setting as a likely depositional environment where jasper and siderite precipitated as alternating bands in response to episodic changes in ambient water chemistry. The episodicity was either driven by fluctuations in the intensity of hydrothermal activity or changes in magma chamber activity, which in turn controlled relative sea level. In this context, precipitation of jasper probably reflects background conditions during which seawater was saturated in silica due to evaporative conditions, while siderites were deposited most likely during intermittent periods of enhanced volcanic activity when seawater was more acidic due to the release of exhalative phases (e.g. CO2). © 2005 Elsevier B.V. All rights reserved.