Monte Carlo simulations of metasomatic enrichment in the lithosphere and implications for the source of alkaline basalts

One hypothesis for the origin of alkaline lavas erupted on oceanic islands and in intracontinental settings is that they represent the melts of amphibole-rich veins in the lithosphere (or melts of their dehydrated equivalents if metasomatized lithosphere is recycled into the convecting mantle). Amphibole-rich veins are interpreted as cumulates produced by crystallization of low-degree melts of the underlying asthenosphere as they ascend through the lithosphere. We present the results of trace-element modelling of the formation and melting of veins formed in this way with the goal of testing this hypothesis and for predicting how variability in the formation and subsequent melting of such cumulates (and adjacent cryptically and modally metasomatized lithospheric peridotite) would be manifested in magmas generated by such a process. Because the high-pressure phase equilibria of hydrous near-solidus melts of garnet lherzolite are poorly constrained and given the likely high variability of the hypothesized accumulation and remelting processes, we used Monte Carlo techniques to estimate how uncertainties in the model parameters (e.g. the compositions of the asthenospheric sources, their trace-element contents, and their degree of melting; the modal proportions of crystallizing phases, including accessory phases, as the asthenospheric partial melts ascend and crystallize in the lithosphere; the amount of metasomatism of the peridotitic country rock; the degree of melting of the cumulates and the amount of melt derived from the metasomatized country rock) propagate through the process and manifest themselves as variability in the trace-element contents and radiogenic isotopic ratios of model vein compositions and erupted alkaline magma compositions. We then compare the results of the models with amphibole observed in lithospheric veins and with oceanic and continental alkaline magmas. While the trace-element patterns of the near-solidus peridotite melts, the initial anhydrous cumulate assemblage (clinopyroxene +/- garnet +/- olivine +/- orthopyroxene), and the modelled coexisting liquids do not match the patterns observed in alkaline lavas, our calculations show that with further crystallization and the appearance of amphibole (and accessory minerals such as rutile, ilmenite, apatite, etc.) the calculated cumulate assemblages have trace-element patterns that closely match those observed in the veins and lavas. These calculated hydrous cumulate assemblages are highly enriched in incompatible trace elements and share many similarities with the trace-element patterns of alkaline basalts observed in oceanic or continental setting such as positive Nb/La, negative Ce/Pb, and similiar slopes of the rare earth elements. By varying the proportions of trapped liquid and thus simulating the cryptic and modal metasomatism observed in peridotite that surrounds these veins, we can model the variations in Ba/Nb, Ce/Pb, and Nb/U ratios that are observed in alkaline basalts. If the isotopic compositions of the initial low-degree peridotite melts are similar to the range observed in mid-ocean ridge basalt, our model calculations produce cumulates that would have isotopic compositions similar to those observed in most alkaline ocean island basalt (OIB) and continental magmas after similar to 0 center dot 15 Gyr. However, to produce alkaline basalts with HIMU isotopic compositions requires much longer residence times (i.e. 1-2 Gyr), consistent with subduction and recycling of metasomatized lithosphere through the mantle. such as a heterogeneous asthenosphere. These modelling results support the interpretation proposed by various researchers that amphibole-bearing veins represent cumulates formed during the differentiation of a volatile-bearing low-degree peridotite melt and that these cumulates are significant components of the sources of alkaline OIB and continental magmas. The results of the forward models provide the potential for detailed tests of this class of hypotheses for the origin of alkaline magmas worldwide and for interpreting major and minor aspects of the geochemical variability of these magmas.

Bentho-planktonic evidence from the Austrian Alps for a decline in sea-surface carbonate production at the end of the Triassic

A high-resolution micropalaeontological study, combined with geochemical and sedimentological analyses was performed on the Tiefengraben, Schlossgraben and Eiberg sections (Austrian Alps) in order to characterize sea-surface carbonate production during the end-Triassic crisis. At the end-Rhaetian, the dominant calcareous nannofossil Prinsiosphaera triassica shows a decrease in abundance and size and this is correlated with a increase in delta O-18 and a gradual decline in delta C-13(carb) values. Simultaneously, benthic foraminiferal assemblages show a decrease in diversity and abundance of calcareous taxa and a dominance of infaunal agglutinated taxa. The smaller size of calcareous nannofossils disturbed the vertical export balance of the biological carbon pump towards the sea-bottom, resulting in changes in feeding strategies within the benthic foraminiferal assemblages from deposit feeders to detritus feeders and bacterial scavengers. These micropalaeontological data combined with geochemical proxies suggest that changes in seawater chemistry and/or cooling episodes might have occurred in the latest Triassic, leading to a marked decrease of carbonate production. This in turn culminated in the quasi-absence of calcareous nannofossils and benthic foraminifers in the latest Triassic. The aftermath (latest Triassic earliest Jurassic) was characterised by abundance peaks of ``disaster'' epifaunal agglutinated foraminifera Trochammina on the sea-floor. Central Atlantic Magmatic Province (CAMP) paroxysmal activity, superimposed on a major worldwide regressive phase, is assumed to be responsible for a deterioration in marine palaeoenvironments. CAMP sulfuric emissions might have been the trigger for cooling episodes and seawater acidification leading to disturbance of the surface carbonate production at the very end-Triassic.

U-Pb and Ar-Ar geochronological data from the Pelagonian basement in Evia (Greece): geodynamic implications for the evolution of Paleotethys

High precision U-Pb zircon and Ar-40/Ar-39 mica geochronological data on metagranodiorites, metagranites and mica schists from north and central Evia island (Greece) are presented in this study. U-Pb zircon ages range from 308 to 1912 Ma, and indicate a prolonged magmatic activity in Late Carboniferous. Proterozoic ages represent inherited cores within younger crystals. Muscovite Ar-40/Ar-39 plateau ages of 288 to 297 Ma are interpreted as cooling ages of the magmatic bodies and metamorphic host rocks in upper greenschist to epidote-amphibolite metamorphic conditions. The multistage magmatism had a duration between 308 and 319 hla but some older intrusions, as well as metamorphic events, cannot be excluded. Geochemical analyses and zircon typology indicate calc-alkaline affinities for the granites of central Evia and alkaline to calc-alkaline characteristics for the metagranodiorites from the northern part of the island. The new data point towards the SE continuation, in Evia and the Cyclades, of a Variscan continental crust already recognised in northern Greece (Pelagonian basement). The Late Carboniferous magmatism is viewed as a result of northward subduction of the Paleotethys under the Eurasian margin.

La gravimétrie expéditive appliquée à la recherche d'aquifères en zone aride. Cas de la nappe alluviale du Teloua (Agadez, Niger)

This paper illustrates the practicality and efficiency of gravimetry for aquifer prospecting in arid zones. Known for the long and tedious data-processing it requires, this method becomes expeditious when simplified as presented here. Its use is then fully justified in a survey of this kind. During the study of the Teloua alluvial aquifer (Agadez, Niger), several ancient channels were clearly and rapidly located. Comparison of the results obtained here with those from previous studies demonstrates anew that for comprehensive prospecting, several complementary geophysical methods should always be employed.

Exploration of geological variability and possible processes through the use of compositional data analysis: an example using scottish metamorphosed

Developments in the statistical analysis of compositional data over the last twodecades have made possible a much deeper exploration of the nature of variability,and the possible processes associated with compositional data sets from manydisciplines. In this paper we concentrate on geochemical data sets. First we explainhow hypotheses of compositional variability may be formulated within the naturalsample space, the unit simplex, including useful hypotheses of subcompositionaldiscrimination and specific perturbational change. Then we develop through standardmethodology, such as generalised likelihood ratio tests, statistical tools to allow thesystematic investigation of a complete lattice of such hypotheses. Some of these tests are simple adaptations of existing multivariate tests but others require specialconstruction. We comment on the use of graphical methods in compositional dataanalysis and on the ordination of specimens. The recent development of the conceptof compositional processes is then explained together with the necessary tools for astaying- in-the-simplex approach, namely compositional singular value decompositions. All these statistical techniques are illustrated for a substantial compositional data set, consisting of 209 major-oxide and rare-element compositions of metamorphosed limestones from the Northeast and Central Highlands of Scotland.Finally we point out a number of unresolved problems in the statistical analysis ofcompositional processes

New insights on river water chemistry by using non-centred simplicial principal component analysis: a case study

The use of perturbation and power transformation operations permits the investigation of linear processes in the simplex as in a vectorial space. When the investigated geochemical processes can be constrained by the use of well-known starting point, the eigenvectors of the covariance matrix of a non-centred principalcomponent analysis allow to model compositional changes compared with a reference point.The results obtained for the chemistry of water collected in River Arno (central-northern Italy) have open new perspectives for considering relative changes of the analysed variables and to hypothesise the relative effect of different acting physical-chemical processes, thus posing the basis for a quantitative modelling

Compositional analysis of archaeological glasses

At CoDaWork'03 we presented work on the analysis of archaeological glass composi-tional data. Such data typically consist of geochemical compositions involving 10-12variables and approximates completely compositional data if the main component, sil-ica, is included. We suggested that what has been termed `crude' principal componentanalysis (PCA) of standardized data often identi ed interpretable pattern in the datamore readily than analyses based on log-ratio transformed data (LRA). The funda-mental problem is that, in LRA, minor oxides with high relative variation, that maynot be structure carrying, can dominate an analysis and obscure pattern associatedwith variables present at higher absolute levels. We investigate this further using sub-compositional data relating to archaeological glasses found on Israeli sites. A simplemodel for glass-making is that it is based on a `recipe' consisting of two `ingredients',sand and a source of soda. Our analysis focuses on the sub-composition of componentsassociated with the sand source. A `crude' PCA of standardized data shows two clearcompositional groups that can be interpreted in terms of di erent recipes being used atdi erent periods, reected in absolute di erences in the composition. LRA analysis canbe undertaken either by normalizing the data or de ning a `residual'. In either case,after some `tuning', these groups are recovered. The results from the normalized LRAare di erently interpreted as showing that the source of sand used to make the glassdi ered. These results are complementary. One relates to the recipe used. The otherrelates to the composition (and presumed sources) of one of the ingredients. It seemsto be axiomatic in some expositions of LRA that statistical analysis of compositionaldata should focus on relative variation via the use of ratios. Our analysis suggests thatabsolute di erences can also be informative

Magma generation at the easternmost section of the Hellenic arc: Hf, Nd, Pb and Sr isotope geochemistry of Nisyros and Yali volcanoes (Greece)

Geochemical and petrographical studies of lavas and ignimbrites from the Quaternary Nisyros-Yali volcanic system in the easternmost part of the Hellenic arc (Greece) reveal insight into magma generating processes. A compositional gap between 61 and 68 wt.% SiO2 is recognized that coincides with the stratigraphic distinction between pre-caldera and postcaldera volcanic units. Trace element systematics support the subdivision of Nisyros and Yali volcanic units into two distinct suites of rocks. The variation of Nd and Hf present day isotope data and the fact that they are distinct from the isotope compositions of MORB rule out an origin by pure differentiation and require assimilation of a crustal component. Lead isotope ratios of Nisyros and Yali volcanic rocks support mixing of mantle material with a lower crust equivalent. However, Sr-87/Sr-86 ratios of 0.7036-0.7048 are incompatible with a simple binary mixing scenario and give low depleted mantle extraction ages (< 0.1 Ga), in contrast with Pb model ages of 0.3 Ga and Hf and Nd model ages of ca. 0.8 Ga. The budget of fluid-mobile elements Sr and Pb is likely to be dominated by abundant hydrous fluids characterised by mantle-like Sr isotope ratios. Late stage fluids probably were enriched in CO2, needed to explain the high Th concentrations. The occurrence of hydrated minerals (e.g., amphibole) in the first post-caldera unit with the lowermost Sr-87/Sr-86 ratio of 0.7036 +/- 2 can be interpreted as the result of the increased water activity in the source. The presence of two different plagioclase phenocryst generations in the first lava subsequent to the caldera-causing event is indicative for a longer storage time of this magma at a shallower level. A model capable of explaining these observations involves three evolutionary stages. First stage, assimilation of lower crustal material by a primitive magma of mantle origin (as modelled by Nd-Hf isotope systematics). This stage ended by an interruption in replenishment that led to an increase of crystallization and, hence, an increase in viscosity, suppressing eruption. During this time gap, differentiation by fractional crystallization led to enrichment of incompatible species, especially aqueous fluids, to silica depolymerisation and to a decrease in viscosity, finally enabling eruption again in the third stage. (c) 2005 Elsevier B.V. All rights reserved.

Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of beingcompositional data and, for a given level of significance, are likely to be beyond thecapabilities of laboratories to distinguish between minute concentrations and completeabsence, thus preventing laboratories from reporting extremely low concentrations of theanalyte. Instead, what is reported is the detection limit, which is the minimumconcentration that conclusively differentiates between presence and absence of theelement. A spatially distributed exhaustive sample is employed in this study to generateunbiased sub-samples, which are further censored to observe the effect that differentdetection limits and sample sizes have on the inference of population distributionsstarting from geochemical analyses having specimens below detection limit (nondetects).The isometric logratio transformation is used to convert the compositional data in thesimplex to samples in real space, thus allowing the practitioner to properly borrow fromthe large source of statistical techniques valid only in real space. The bootstrap method isused to numerically investigate the reliability of inferring several distributionalparameters employing different forms of imputation for the censored data. The casestudy illustrates that, in general, best results are obtained when imputations are madeusing the distribution best fitting the readings above detection limit and exposes theproblems of other more widely used practices. When the sample is spatially correlated, itis necessary to combine the bootstrap with stochastic simulation

A compositional approach to stable isotope data analysis

Isotopic data are currently becoming an important source of information regardingsources, evolution and mixing processes of water in hydrogeologic systems. However, itis not clear how to treat with statistics the geochemical data and the isotopic datatogether. We propose to introduce the isotopic information as new parts, and applycompositional data analysis with the resulting increased composition. Results areequivalent to downscale the classical isotopic delta variables, because they are alreadyrelative (as needed in the compositional framework) and isotopic variations are almostalways very small. This methodology is illustrated and tested with the study of theLlobregat River Basin (Barcelona, NE Spain), where it is shown that, though verysmall, isotopic variations comp lement geochemical principal components, and help inthe better identification of pollution sources

Assessing sedimentary evolution by means of Sr-isotope ratios : 3 case studies on the Caribbean plate : (Cretaceous: Nicoya Peninsula, Costa Rica, Tertiary: Hess Rise, and La Désirade, Guadeloupe, France)

The understanding of sedimentary evolution is intimately related to the knowledge of the exact ages of the sediments. When working on carbonate sediments, age dating is commonly based on paleontological observations and established biozonations, which may prove to be relatively imprecise. Dating by means of strontium isotope ratios in marine bioclasts is the probably best method in order to precisely date carbonate successions, provided that the sample reflects original marine geochemical characteristics. This requires a precise study of the samples including its petrography, SEM and cathodoluminescence observations, stable carbon and oxygen isotope geochemistry and finally the strontium isotope measurement itself. On the Nicoya Peninsula (Northwestern Costa Rica) sediments from the Piedras Blancas Formation, Nambi Formation and Quebrada Pavas Formation were dated by the means of strontium isotope ratios measured in Upper Cretaceous Inoceramus shell fragments. Results have shown average 87Sr/86Sr values of 0.707654 (middle late Campanian) for the Piedras Blancas Formation, 0.707322 (Turonian-Coniacian) for the Nambi Formation and 0.707721 (late Campanian-Maastrichtian) for the Quebrada Pavas Formation. Abundant detrital components in the studied formations constitute a difficulty to strontium isotope dating. In fact, the fossil bearing sediments can easily contaminate the target fossil with strontium mobilized form basalts during diagenesis and thus the obtained strontium isotope ratios may be influenced significantly and so will the obtained ages. The new and more precise age assignments allow for more precision in the chronostratigraphic chart of the sedimentary and tectonic evolution of the Nicoya Peninsula, providing a better insight on the evolution of this region. Meteor Cruise M81 dredged shallow water carbonates from the Hess Rise and Hess Escarpment during March 2010. Several of these shallow water carbonates contain abundant Larger Foraminifera that indicates an Eocene-Oligocene age. In this study the strontium isotope values ranging from 0.707847 to 0.708238 can be interpreted as a Rupelian to Chattian age of these sediments. These platform sediments are placed on seamounts, now located at depths reaching 1600 m. Observation of sedimentologic characteristics of these sediments has helped to resolve apparent discrepancies between fossil and strontium isotope ages. Hence, it is possible to show that the subsidence was active during early Miocene times. On La Désirade (Guadeloupe France), the Neogene to Quaternary carbonate cover has been dated by microfossils and some U/Th-ages. Disagreements subsisted in the paleontological ages of the formations. Strontium isotope ratios ranging from 0.709047 to 0.709076 showed the Limestone Table of La Désirade to range from an Early Pliocene to Late Pliocene/early Pleistocene age. A very late Miocene age (87Sr/86Sr =0.709013) can be determined to the Detrital Offshore Limestone. The flat volcanic basement had to be eroded by wave-action during a long-term stable relative sea-level. Sediments of the Table Limestone on La Désirade show both low-stand and high-stand facies that encroach on the igneous basement, implying deposition during a major phase of subsidence creating accommodation space. Subsidence is followed by tectonic uplift documented by fringing reefs and beach rocks that young from the top of the Table Limestone (180 m) towards the present coastline. Strontium isotope ratios from two different fringing reefs (0.707172 and 0.709145) and from a beach rock (0.709163) allow tentative dating, (125ky, ~ 400ky, 945ky) and indicate an uplift rate of about 5cm/ky for this time period of La Désirade Island. The documented subsidence and uplift history calls for a new model of tectonic evolution of the area.

Sr, C and O isotope systematics in the Pucara Basin, central Peru - Comparison between Mississippi Valley-type deposits and barren areas

A combined Sr, O and C isotope study has been carried out in the Pucara basin, central Peru, to compare local isotopic trends of the San Vicente and Shalipayco Zn-Pb Mississippi Valley-type (MVT) deposits with regional geochemical patterns of the sedimentary host basin. Gypsum, limestone and regional replacement dolomite yield Sr-87/Sr-86 ratios that fall within or slightly below the published range of seawater Sr-87/Sr-86 values for the Lower Jurassic and the Upper Triassic. Our data indicate that the Sr isotopic composition of seawater between the Hettangian and the Toarcian may extend to lower Sr-87/Sr-86 ratios than previously published values. An Sr-87-enrichment is noted in (1) carbonate rocks from the lowermost part of the Pucara basin, and (2) different carbonate generations at the MVT deposits. This indicates that host rocks at MVT deposits and in the lowermost part of the carbonate sequence interacted with Sr-87-enriched fluids. The fluids acquired their radiogenic nature by interaction with lithologies underlying the carbonate rocks of the Pucara basin. The San Ramon granite, similar Permo-Triassic intrusions and their elastic derivatives in the Mitu Group are likely sources of radiogenic Sr-87. The Brazilian shield and its erosion products are an additional potential source of radiogenic Sr-87. Volcanic rocks of the Mitu Group are not a significant source for radiogenic Sr-87; however, molasse-type sedimentary rocks and volcaniclastic rocks cannot be ruled out as a possible source of radiogenic Sr-87. The marked enrichment in Sr-87 of carbonates toward the lower part of the Pucara Group is accompanied by only a slight decrease in delta(18)O values and essentially no change in delta(13)C values, whereas replacement dolomite and sparry carbonates at the MVT deposits display a coherent trend of progressive Sr-87-enrichment, and O-18- and C-13-depletion. The depletion in O-18 in carbonates from the MVT deposits are likely related to a temperature increase, possibly coupled with a O-18-enrichment of the ore-forming fluids. Progressively lower delta(13)C values throughout the paragenetic sequence at the MVT deposits are interpreted as a gradually more important contribution from organically derived carbon. Quantitative calculations show that a single fluid-rock interaction model satisfactorily reproduces the marked Sr-87-enrichment and the slight decrease in delta(18)O values in carbonate rocks from the lower part of the Pucara Group. By contrast, the isotopic covariation trends of the MVT deposits are better reproduced by a model combining fluid mixing and fluid-rock interaction. The modelled ore-bearing fluids have a range of compositions between a hot, saline, radiogenic brine that had interacted with lithologies underlying the Pucara sequence and cooler, dilute brines possibly representing local fluids within the Pucara sequence. The composition of the local fluids varies according to the nature of the lithologies present in the neighborhood of the different MVT deposits. The proportion of the radiogenic fluid in the modelled fluid mixtures interacting with the carbonate host rocks at the MVT deposits decreases as one moves up in the stratigraphic sequence of the Pucara Group.

Sulfur and strontium isotope composition of the Llobregat River (NE Spain): Tracers of natural and anthropogenic chemicals in stream waters

The use of sulfur and strontium isotopes as tracers for the source/s of water contaminants have been applied to the water of the Llobregat River system (NE Spain). Surface water samples from June 1997 were collected from the Llobregat River and its main tributaries and creeks. The chemistry of most stream waters are controlled mainly by the weathering of Tertiary chemical sediments within the drainage basin. The largest variation in delta(34)S values were found in the small creeks with values ranging from -9.9 to 15parts per thousand, whilst in the main river channels values ranged from 6.3 to 12.4parts per thousand. The Sr-87/Sr-86 ratio for dissolved strontium ranged from 0.70795 for a non-polluted site to 0.70882 for a polluted one. Most of the waters with high NO3 and low Ca/Na ratio converge to the same Sr-87/Sr-86 value, pointing to dominant pollutant end member contribution or a mixing of pollutants with an isotopic composition around 0.7083-0.7085. Although the concentration of the natural inputs in the river for sulfate and strontium are high, as a result of the sulfate outcrops within the geology of the basin, their isotopic characteristics suggest that they can be used as a discriminating device in water pollution problems. However to establish the detailed characteristics of the isotopes as geochemical tools, specific high-resolution case studies are necessary in small areas, where the inputs are well known.

Multi-isotopic and compositional exploration of factors controlling nitrate pollution

In Catalonia, according to the nitrate directive (91/676/EU), nine areas have been declared as vulnerable to nitrate pollution from agricultural sources (Decret 283/1998 and Decret 479/2004). Five of these areas have been studied coupling hydro chemical data with a multi-isotopic approach (Vitòria et al. 2005, Otero et al. 2007, Puig et al. 2007), in an ongoing research project looking for an integrated application of classical hydrochemistry data, with a comprehensive isotopic characterisation (δ15N and δ18O of dissolved nitrate, δ34S and δ18O of dissolved sulphate, δ13C of dissolved inorganic carbon, and δD and δ18O of water). Within this general frame, the contribution presented explores compositional ways of: (i) distinguish agrochemicals and manure N pollution, (ii) quantify natural attenuation of nitrate (denitrification), and identify possible controlling factors.To achieve this two-fold goal, the following techniques have been used. Separate biplots of each suite of data show that each studied region has a distinct δ34S and pH signatures, but they are homogeneous with regard to NO3- related variables. Also, the geochemical variables were projected onto the compositional directions associated with the possible denitrification reactions in each region. The resulting balances can be plot together with some isotopes, to assess their likelihood of occurrence

Studies in stoichiometry with compositional data

Geochemical data that is derived from the whole or partial analysis of various geologic materialsrepresent a composition of mineralogies or solute species. Minerals are composed of structuredrelationships between cations and anions which, through atomic and molecular forces, keep the elementsbound in specific configurations. The chemical compositions of minerals have specific relationships thatare governed by these molecular controls. In the case of olivine, there is a well-defined relationshipbetween Mn-Fe-Mg with Si. Balances between the principal elements defining olivine composition andother significant constituents in the composition (Al, Ti) have been defined, resulting in a near-linearrelationship between the logarithmic relative proportion of Si versus (MgMnFe) and Mg versus (MnFe),which is typically described but poorly illustrated in the simplex.The present contribution corresponds to ongoing research, which attempts to relate stoichiometry andgeochemical data using compositional geometry. We describe here the approach by which stoichiometricrelationships based on mineralogical constraints can be accounted for in the space of simplicialcoordinates using olivines as an example. Further examples for other mineral types (plagioclases andmore complex minerals such as clays) are needed. Issues that remain to be dealt with include thereduction of a bulk chemical composition of a rock comprised of several minerals from which appropriatebalances can be used to describe the composition in a realistic mineralogical framework. The overallobjective of our research is to answer the question: In the cases where the mineralogy is unknown, arethere suitable proxies that can be substituted?Kew words: Aitchison geometry, balances, mineral composition, oxides