Grain size analyses in tin-lead glazes based on 2d-sections

A problem in the archaeometric classification of Catalan Renaissance pottery is the fact, thatthe clay supply of the pottery workshops was centrally organized by guilds, and thereforeusually all potters of a single production centre produced chemically similar ceramics.However, analysing the glazes of the ware usually a large number of inclusions in the glaze isfound, which reveal technological differences between single workshops. These inclusionshave been used by the potters in order to opacify the transparent glaze and to achieve a whitebackground for further decoration.In order to distinguish different technological preparation procedures of the single workshops,at a Scanning Electron Microscope the chemical composition of those inclusions as well astheir size in the two-dimensional cut is recorded. Based on the latter, a frequency distributionof the apparent diameters is estimated for each sample and type of inclusion.Following an approach by S.D. Wicksell (1925), it is principally possible to transform thedistributions of the apparent 2D-diameters back to those of the true three-dimensional bodies.The applicability of this approach and its practical problems are examined using differentways of kernel density estimation and Monte-Carlo tests of the methodology. Finally, it istested in how far the obtained frequency distributions can be used to classify the pottery

The chemical variability at the surface of Mars: implications for exogenic processes

The chemical composition of sediments and rocks, as well as their distribution at theMartian surface, represent a long term archive of processes, which have formed theplanetary surface. A survey of chemical compositions by means of Compositional DataAnalysis represents a valuable tool to extract direct evidence for weathering processesand allows to quantify weathering and sedimentation rates. clr-biplot techniques areapplied for visualization of chemical relationships across the surface (“chemical maps”).The variability among individual suites of data is further analyzed by means of clr-PCA,in order to extract chemical alteration vectors between fresh rocks and their crusts andfor an assessment of different source reservoirs accessible to soil formation. Bothtechniques are applied to elucidate the influence of remote weathering by combinedanalysis of several soil forming branches. Vector analysis in the Simplex provides theopportunity to study atmosphere surface interactions, including the role andcomposition of volcanic gases

Statistical evidences of cyclic changes in volcanic gas chemistry composition by inverse modelling

The identification of compositional changes in fumarolic gases of active and quiescent volcanoes is one of the mostimportant targets in monitoring programs. From a general point of view, many systematic (often cyclic) and randomprocesses control the chemistry of gas discharges, making difficult to produce a convincing mathematical-statisticalmodelling.Changes in the chemical composition of volcanic gases sampled at Vulcano Island (Aeolian Arc, Sicily, Italy) fromeight different fumaroles located in the northern sector of the summit crater (La Fossa) have been analysed byconsidering their dependence from time in the period 2000-2007. Each intermediate chemical composition has beenconsidered as potentially derived from the contribution of the two temporal extremes represented by the 2000 and 2007samples, respectively, by using inverse modelling methodologies for compositional data. Data pertaining to fumarolesF5 and F27, located on the rim and in the inner part of La Fossa crater, respectively, have been used to achieve theproposed aim. The statistical approach has allowed us to highlight the presence of random and not random fluctuations,features useful to understand how the volcanic system works, opening new perspectives in sampling strategies and inthe evaluation of the natural risk related to a quiescent volcano

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.