Elemental analysis of New Zealand Nephrite Jade (pounamu) for forensic purposes

Pounamu (NZ jade), or nephrite, is a protected mineral in its natural form following thetransfer of ownership back to Ngai Tahu under the Ngai Tahu (Pounamu Vesting) Act 1997.Any theft of nephrite is prosecutable under the Crimes Act 1961. Scientific evidence isessential in cases where origin is disputed. A robust method for discrimination of thismaterial through the use of elemental analysis and compositional data analysis is required.Initial studies have characterised the variability within a given nephrite source. This hasincluded investigation of both in situ outcrops and alluvial material. Methods for thediscrimination of two geographically close nephrite sources are being developed.Key Words: forensic, jade, nephrite, laser ablation, inductively coupled plasma massspectrometry, multivariate analysis, elemental analysis, compositional data analysis

A Factor analysis of hidrochemical composition of Llobregat river basin

Hydrogeological research usually includes some statistical studies devised to elucidate mean background state, characterise relationships among different hydrochemical parameters, and show the influence of human activities. These goals are achieved either by means of a statistical approach or by mixing modelsbetween end-members. Compositional data analysis has proved to be effective with the first approach, but there is no commonly accepted solution to the end-member problem in a compositional framework.We present here a possible solution based on factor analysis of compositions illustrated with a case study.We find two factors on the compositional bi-plot fitting two non-centered orthogonal axes to the most representative variables. Each one of these axes defines a subcomposition, grouping those variables thatlay nearest to it. With each subcomposition a log-contrast is computed and rewritten as an equilibrium equation. These two factors can be interpreted as the isometric log-ratio coordinates (ilr) of three hiddencomponents, that can be plotted in a ternary diagram. These hidden components might be interpreted as end-members.We have analysed 14 molarities in 31 sampling stations all along the Llobregat River and its tributaries, with a monthly measure during two years. We have obtained a bi-plot with a 57% of explained totalvariance, from which we have extracted two factors: factor G, reflecting geological background enhanced by potash mining; and factor A, essentially controlled by urban and/or farming wastewater. Graphicalrepresentation of these two factors allows us to identify three extreme samples, corresponding to pristine waters, potash mining influence and urban sewage influence. To confirm this, we have available analysisof diffused and widespread point sources identified in the area: springs, potash mining lixiviates, sewage, and fertilisers. Each one of these sources shows a clear link with one of the extreme samples, exceptfertilisers due to the heterogeneity of their composition.This approach is a useful tool to distinguish end-members, and characterise them, an issue generally difficult to solve. It is worth note that the end-member composition cannot be fully estimated but only characterised through log-ratio relationships among components. Moreover, the influence of each endmember in a given sample must be evaluated in relative terms of the other samples. These limitations areintrinsic to the relative nature of compositional data

Quantifying rock fabrics: a test of autocorrelation of the spatial distribution of cristals

A novel test of spatial independence of the distribution of crystals or phases in rocksbased on compositional statistics is introduced. It improves and generalizes the commonjoins-count statistics known from map analysis in geographic information systems.Assigning phases independently to objects in RD is modelled by a single-trial multinomialrandom function Z(x), where the probabilities of phases add to one and areexplicitly modelled as compositions in the K-part simplex SK. Thus, apparent inconsistenciesof the tests based on the conventional joins{count statistics and their possiblycontradictory interpretations are avoided. In practical applications we assume that theprobabilities of phases do not depend on the location but are identical everywhere inthe domain of de nition. Thus, the model involves the sum of r independent identicalmultinomial distributed 1-trial random variables which is an r-trial multinomialdistributed random variable. The probabilities of the distribution of the r counts canbe considered as a composition in the Q-part simplex SQ. They span the so calledHardy-Weinberg manifold H that is proved to be a K-1-affine subspace of SQ. This isa generalisation of the well-known Hardy-Weinberg law of genetics. If the assignmentof phases accounts for some kind of spatial dependence, then the r-trial probabilitiesdo not remain on H. This suggests the use of the Aitchison distance between observedprobabilities to H to test dependence. Moreover, when there is a spatial uctuation ofthe multinomial probabilities, the observed r-trial probabilities move on H. This shiftcan be used as to check for these uctuations. A practical procedure and an algorithmto perform the test have been developed. Some cases applied to simulated and realdata are presented.Key words: Spatial distribution of crystals in rocks, spatial distribution of phases,joins-count statistics, multinomial distribution, Hardy-Weinberg law, Hardy-Weinbergmanifold, Aitchison geometry

When zero doesn't mean it and other geomathematical mischief

There is almost not a case in exploration geology, where the studied data doesn’tincludes below detection limits and/or zero values, and since most of the geological dataresponds to lognormal distributions, these “zero data” represent a mathematicalchallenge for the interpretation.We need to start by recognizing that there are zero values in geology. For example theamount of quartz in a foyaite (nepheline syenite) is zero, since quartz cannot co-existswith nepheline. Another common essential zero is a North azimuth, however we canalways change that zero for the value of 360°. These are known as “Essential zeros”, butwhat can we do with “Rounded zeros” that are the result of below the detection limit ofthe equipment?Amalgamation, e.g. adding Na2O and K2O, as total alkalis is a solution, but sometimeswe need to differentiate between a sodic and a potassic alteration. Pre-classification intogroups requires a good knowledge of the distribution of the data and the geochemicalcharacteristics of the groups which is not always available. Considering the zero valuesequal to the limit of detection of the used equipment will generate spuriousdistributions, especially in ternary diagrams. Same situation will occur if we replace thezero values by a small amount using non-parametric or parametric techniques(imputation).The method that we are proposing takes into consideration the well known relationshipsbetween some elements. For example, in copper porphyry deposits, there is always agood direct correlation between the copper values and the molybdenum ones, but whilecopper will always be above the limit of detection, many of the molybdenum values willbe “rounded zeros”. So, we will take the lower quartile of the real molybdenum valuesand establish a regression equation with copper, and then we will estimate the“rounded” zero values of molybdenum by their corresponding copper values.The method could be applied to any type of data, provided we establish first theircorrelation dependency.One of the main advantages of this method is that we do not obtain a fixed value for the“rounded zeros”, but one that depends on the value of the other variable.Key words: compositional data analysis, treatment of zeros, essential zeros, roundedzeros, correlation dependency

Modelling changes in stable isotope compositions of minerals during net transfer reactions in a contact aureole: Wollastonite growth at the northern Hunter Mountain Batholith (Death Valley National Park, USA)

One of the world's largest wollastonite deposits was formed at the contact of the northern Hunter Mountain Batholith (California, USA) in Paleozoic sediments. Wollastonite occurs as zones of variable thickness surrounding layers or nodules of quartzite in limestones. A minimum formation temperature of 650 degrees C is estimated from isolated periclase-bearing lenses in that area. Contact metamorphism of siliceous carbonates has produced mineral assemblages that are consistent with heterogeneous, and partly limited infiltration of water-rich fluids, compatible with O-18/O-16 and C-13/C-12 isotopic patterns recorded in carbonates. Oxygen isotope compositions of wollastonites in the study area may also not require infiltration of large quantities of externally-derived fluids that were out of equilibrium with the rocks. 8180 values of wollastonite are high (14.8 parts per thousand to 25.0 parts per thousand; median: 19.7 parts per thousand) and close to those of the host limestone (19.7 parts per thousand to 28 parts per thousand; median: 24.9 parts per thousand) and quartz (18.0 parts per thousand. to 29.1 parts per thousand; median: 22.6 parts per thousand). Isotopic disequilibrium exists at quartz/wollastonite and wollastonite/calcite boundaries. Therefore, classical batch/Rayleigh fractionation models based on reactant and product equilibrium are not applicable to the wollastonite rims. An approach that relies on local instantaneous mass balance for the reactants, based on the wollastonite-forming reaction is suggested as an alternative way to model wollastonite reaction rims. This model reproduces many of the measured delta O-18 values of wollastonite reaction rims of the current study to within +/- 1 parts per thousand, even though the wollastonite compositions vary by almost 10 parts per thousand. (C) 2011 Elsevier B.V. All rights reserved.

Indicator taxa of benthic diatom communities: a case study in Mediterranean streams

A key issue in the implementation of the Water Framework Directive is the classification of streams and rivers using biological quality parameters and type-specific reference conditions. Four groups of stream types were defined in NE Spain on the basis of 152 diatom samples by means of detrended correspondence analysis and classification techniques. Diatom analysis was restricted to epilithic taxa, and the sites included gradients ranging from near-natural streams to sites with poor ecological quality. The main gradient shows a clear separation of sites in relation to the degree of human influence: polluted streams (mainly located in the lowlands) differ from streams in mountainous areas and in the Pyrenees. A second gradient is related to physiographical features. Headwater streams can be distinguished by their catchment geology. The type-specific diatom taxa for the stream types studied were determined by using indicator species analysis (IndVal). The type-specific taxa from near-natural streams are coincident with the indicator taxa for high ecological status. Human impact reduced the typological heterogeneity of the diatom community composition. Overall, the diatom communities in NE Spain exhibit a regional distribution pattern that closely corresponds with that observed in river systems elsewhere. Physiographical differences are only evident in undisturbed sites, while nutrient enrichment and other human disturbances may mask the regional differences in the distribution of diatom communities

Reconstructing paleoelevation in eroded orogens

Hydrogen isotope and Ar-40/Ar-39 geochronological data are presented from muscovite within a crustal-scale extensional detachment of the Shuswap Metamorphic Complex, North American Cordillera. The hydrogen isotope compositions (deltaD(ms)) of precisely dated muscovite attain values as low as -156parts per thousand in the detachment mylonite, whereas footwall quartzite has a deltaD(ms) value of -81parts per thousand. The very low deltaD(ms) values in the detachment are best explained by infiltration of meteoric water, with maximum deltaD values of -135parts per thousand +/- 3parts per thousand, during extensional unroofing of the orogen at 49.0-47.9 Ma. On the basis of the empirically determined relationship between elevation and isotopic composition of precipitation, the reconstructed early Eocene paleoelevations of the orogen are 4060 +/- 250 m to 4320 +/- 250 m, at least 1000 m higher than the highest present-day peaks. We propose that the isotopic composition of surface-derived waters in extensional detachments represents a newly recognized method to estimate maximum paleoelevations attained immediately preceding extensional orogenic collapse.

Microstructural, chemical and textural records during growth of snowball garnet

The growth history of two populations of snowball garnet from the Lukmanier Pass area (central Swiss Alps) was examined through a detailed analysis of three-dimensional geometry, chemical zoning and crystallographic orientation. The first population, collected in the hinge of a chevron-type fold, shows an apparent rotation of 360 degrees. The first 270 degrees are characterized by spiral-shaped inclusion trails, gradual and concentric Mn zoning and a single crystallographic orientation, whereas in the last 90 degrees, crenulated inclusion trails and secondary Mn maxima centred on distinct crystallographic garnet domains are observed. Microstructural, geochemical and textural data indicate a radical change in growth regime between the two growth sequences. In the first 270 degrees, growth occurred under rotational non-coaxial flow, whereas in the last 90 degrees, garnet grew under a non-rotational shortening regime. The second population, collected in the limb of the same chevron-type fold structure, is characterized by a spiral geometry that does not exceed 270 degrees of apparent rotation. These garnet microstructures do not record any evidence for a modification of the stress field during garnet growth. Concentric Mn zoning as well as a single crystallographic orientation are observed for the entire spiral. Electron backscatter diffraction data indicate that nearly all central domains in the snowball garnet are characterized by one [001] axis oriented (sub-)parallel to the symmetry axis and by another [001] axis oriented (sub-)parallel to the orientation of the internal foliation. These features suggest that the crystallographic orientation across the garnet spiral is not random and that a relation exists among the symmetry axis, the internal foliation and the crystallographic orientation.

Geologia del Baix Ter

Article que analitza el valor i lasingularitat de les zones humides del Baix Ter, tenint en compte la descripció de les característiques geològiques de la Plana del Ter i dels diferents ambients morfodinàmics que s' hi troben

Three-dimensional sequence stratigraphy and subtle stratigraphic traps associated with systems tracts: West Cameron region, offshore Louisiana, Gulf of Mexico

Three-dimensional sequence stratigraphy is a potent exploration and development tool for the discovery of subtle stratigraphic traps. Reservoir morphology, heterogeneity and subtle stratigraphic trapping mechanisms can be better understood through systematic horizontal identification of sedimentary facies of systems tracts provided by three-dimensional attribute maps used as an important complement to the sequential analysis on the two-dimensional seismic lines and the well log data. On new prospects as well as on already-producing fields, the additional input of sequential analysis on three-dimensional data enables the identification, location and precise delimitation of new potentially productive zones. The first part of this paper presents four typical horizontal seismic facies assigned to the successive systems tracts of a third- or fourth-order sequence deposited in inner to outer neritic conditions on a elastic shelf. The construction of this synthetic representative sequence is based on the observed reproducibility of the horizontal seismic facies response to cyclic eustatic events on more than 35 sequences registered in the Gulf coast Plio-Pleistocene and Late Miocene, offshore Louisiana in the West Cameron region of the Gulf of Mexico. The second part shows how three-dimensional sequence stratigraphy can contribute in localizing and understanding sedimentary facies associated with productive zones. A case study in the early Middle Miocene Cibicides opima sands shows multiple stacked gas accumulations in the top slope fan, prograding wedge and basal transgressive systems tract of the third-order sequence between SB15.5 and SB 13.8 Ma.

Bitumens in the late Variscan hydrothermal vein-type uranium deposit of Pribram, Czech Republic: Sources, radiation-induced alteration, and relation to mineralization

The late Variscan (275-278 Ma) Pribram uranium deposit is one of the largest known accumulations of uraniferous bitumens in hydrothermal veins. The deposit extends along the northwestern boundary of the Central Bohemian pluton (345-335 Ma) with low-grade metamorphosed Late Proterozoic and unmetamorphosed Cambrian rocks. From a net uranium production of 41,742 metric tons (t), more than 6,000 t were extracted from bitumen-uraninite ores during 43 years of exploration and mining. Three morphological varieties of solid bitumen are recognized: globular, asphaltlike, and cokelike. While the globular bitumen is uranium free, the other two types are uraniferous. The amount of bitumen in ore veins gradually decreases toward the contact with the plutonic body and increases with depth. Two types of bitumen microtextures are recognized using high-resolution transmission electron microscopy: amorphous and microporous, the former being less common in uraniferous samples. A lower Raman peak area ratio (1,360/1,575 cm(-1)) in mineralized bitumens (0.9) compared with uranium-free samples (2.0) indicates a lower degree of microtextural organization in the latter The H/C and O/C atomic ratios in uranium-free bitumens (0.9-1.1 and 0.09, respectively) are higher than those in mineralized samples (H/C = 0.3-0.8, O/C = 0.03-0.09). The chloroform extractable matter yield is Very low in uranium-free bitumens (0.30-0.35% of the total organic carbon,TOC) and decreases with uranium content increase. The extracted solid uraniferous bitumen infrared spectra show depletion in aliphatic CH2 and CH3 groups compared to uranium-free samples. The concentration of oxygen-bearing functional groups relative to aromatic bonds in the IR spectra of uranium-free and mineralized bitumen, however, do not differ significantly. C-13 NMR confirmed than the aromaticity of a uraniferous sample is higher (F-ar = 0.61) than in the uranium-free bitumen (F-ar = 0.51). Pyrolysates from uraniferous and nonuraniferous bitumens do not differ significantly, being predominantly cresol, alkylphenols, alkylbenzenes, and alkylnaphthalenes. The liquid pyrolysate yield decreases significantly with increasing uranium content. The delta(13)C Values of bulk uranium-free bitumens and low-grade uraniferous, asphaltlike bitumens range from -43.6 to 52.3 per mil. High-grade, cokelike, uraniferous bitumens are more C-13 depleted (54.5 to -58.4 parts per thousand). In contrast to the very light isotopic ratios of the high-grade uraniferous cokelike bitumen bulk carbon, the individual n-alkanes and isoprenoids (pristane and phytane) extracted from the same sample are significantly C-13 enriched. The isotopic composition of the C13-24 n-alkanes extracted from the high-grade uraniferous sample (delta(13)C = -28.0 to 32.6 parts per thousand) are heavier compared with the same compounds in a uranium-free sample (delta(13)C = 31.9 to 33.8 parts per thousand). It is proposed that the bitumen source was the isotopically light (delta(13)C = 35.8 to 30.2 parts per thousand) organic matter of the Upper Proterozoic host rocks that were pyrolyzed during intrusion of the Central Bohemian pluton. The C-13- depleted pyrolysates were mobilized from the innermost part of the contact-metamorphic aureole, accumulated in structural traps in less thermally influenced parts of the sedimentary complex and were later extracted by hydrothermal fluids. Bitumens at the Pribram deposit are younger than the main part of the uranium mineralization and were formed through water-washing and radiation-induced polymerization of both the gaseous and liquid pyrolysates. Direct evidence for pyrolysate reduction of uranium in the hydrothermal system is difficult to obtain as the chemical composition of the original organic fluid phase was modified during water-washing and radiolytic alteration. However, indirect evidence-e.g., higher O/C atomic ratios in uranium-free bitumens (0.1) relative to the Upper Proterozoic source rocks (0.02-0.05), isotopically very light carbon in associated whewellite (delta(13)C = 31.7 to -28.4 parts per thousand), and the striking absence of bitumens in the pre-uranium, hematite stage of the mineralization-indicates that oxidation of organic fluids may have contributed to lowering of aO(2) and uraninite precipitation.

Oxygen isotope compositions of iron oxides from high-grade BIF-hosted iron ore deposits of the Central Hamersley Province, Western Australia: Constraints on the evolution of hydrothermal fluids

The Hamersley province of northwest Australia is one of the world's premier iron ore regions with high-grade martite-microplaty hematite iron ore deposits mostly hosted within banded iron formation (BIF) sequences of the Brockman Iron Formations of the Hamersley Group. These high-grade iron ores contain between 60 and 68 wt percent Fe, and formed by the multistage interaction of hydrothermal fluids with the host BIF formation. The oxygen isotope compositions of magnetite and hematite from BIF, hydrothermal alteration assemblages, and high-grade iron Ore were analyzed from the Mount Tom Price, Paraburdoo, and Charmar iron ore deposits. The delta(18)O values of magnetite and hematite from hydrothermal alteration assemblages and high-grade iron ore range from -9.0 to -2.9 per mil, a depletion of 5 to 15 per mil relative to the host BIF. The delta(18)O values are spatially controlled by faults within the deposits, a response to higher fluid flux and larger influence the isotopic compositions by the hydrothermal fluids. The oxygen isotope composition of hydrothermal fluids (delta(18)O(fluid)) indicates that the decrease in the (18)O content of iron oxides was due to the interaction of both basinal brines and meteoric fluids with the original BIF. Late-stage talc-bearing ore at the Mount Tom Price deposit formed in the presence of a pulse of delta(18)O-enriched basinal brine, indicating that hydrothermal fluids may have repeatedly interacted with the BIFs during the Paleoproterozoic.

Anachronistic facies from a drowned Lower Triassic carbonate platform: Lower member of the Alwa Formation (Ba'id Exotic), Oman Mountains

- The lower member of the Alwa Formation (Lower Olenekian), found within the Ba'id Exotic in the Oman Mountains (Sultanate of Oman), consists of ammonoid-bearing, pelagic limestones that were deposited on an isolated, drowned carbonate platform on the Neotethyan Gondwana margin. The strata contain a variety of unusual carbonate textures and features, including thrombolites, Frutexites-bearing microbialites that contain synsedimentary cements, matrix-free breccias surrounded by isopachous calcite cement, and fissures and cavities filled with large botryoidal cements. Thrombolites are found throughout the study interval, and occur as 0.5-1.0 m thick lenses or beds that contain laterally laterally-linked stromatactis cavities. The Frutexites-bearing microbialites occur less frequently, and also form lenses or beds, up to 30 cm thick; the microbialites may be laminated, and often developed on hardgrounds. In addition, the Frutexites-bearing microbialites also contain synsedimentary calcite cement crusts and botryoids (typically <1 cm thick) that harbour layers or pockets of what appear to be bacterial sheaths and coccoids, and are indicative of biologically mediated precipitation of the cement bodies. Slumping following lithification led to fracturing of the limestone and the precipitation of large, botryoidal aragonite cements in fissures that cut across the primary fabric. Environmental conditions, specifically palaeoxygenation and the degree of calcium carbonate supersaturation, likely controlled whether the thrombolites (high level of calcium carbonate supersaturation associated with vertical mixing of water masses and dysoxic conditions) or Frutexites-bearing microbialites (low level of calcium carbonate supersaturation associated with anoxic conditions and deposition below a stable chemocline) formed. The results of this study point to continued environmental stress in the region during the Early Triassic that likely contributed to the uneven recovery from the Permian-Triassic mass extinction.

Preservation of an extreme transient geotherm in the Raft River detachment shear zone

Extensional detachment systems separate hot footwalls from cool hanging walls, but the degree to which this thermal gradient is the product of ductile or brittle deformation or a preserved original transient geotherm is unclear. Oxygen isotope thermometry using recrystallized quartz-muscovite pairs indicates a smooth thermal gradient (140 degrees C/100 m) across the gently dipping, quartzite-dominated detachment zone that bounds the Raft River core complex in northwest Utah (United States). Hydrogen isotope values of muscovite (delta D-Ms similar to-100 parts per thousand) and fluid inclusions in quartz (delta D-Fluid similar to-85 parts per thousand) indicate the presence of meteoric fluids during detachment dynamics. Recrystallized grain-shape fabrics and quartz c-axis fabric patterns reveal a large component of coaxial strain (pure shear), consistent with thinning of the detachment section. Therefore, the high thermal gradient preserved in the Raft River detachment reflects the transient geotherm that developed owing to shearing, thinning, and the potentially prominent role of convective flow of surface fluids.

The Fry method applied to an augen orthogneiss: Problems and results

The application of the Fry method to measure strain in deformed porphyritic granites is discussed. This method requires that the distribution of markers has to satisfy at least two conditions. It has to be homogeneous and isotropic. Statistics on point distribution with the help of a Morishita diagram can easily test homogeneity. Isotropy can be checked with a cumulative histogram of angles between points. Application of these tests to undeformed (Mte Capanne granite, Elba) and to deformed (Randa orthogneiss, Alps of Switzerland) porphyritic granite reveals that their K-feldspars phenocrysts both satisfy these conditions and can be used as strain markers with the Fry method. Other problems are also examined. One is the possible distribution of deformation on discrete shear-bands. Providing several tests are met, we conclude that the Fry method can be used to estimate strain in deformed porphyritic granites. (c) 2006 Elsevier Ltd. All rights reserved.