Sedimentologic and magnetic data of sediment cores on the western slope of the Mid-Atlantic Ridge

The relative paleointensity (RPI) method assumes that the intensity of post depositional remanent magnetization (PDRM) depends exclusively on the magnetic field strength and the concentration of the magnetic carriers. Sedimentary remanence is regarded as an equilibrium state between aligning geomagnetic and randomizing interparticle forces. Just how strong these mechanical and electrostatic forces are, depends on many petrophysical factors related to mineralogy, particle size and shape of the matrix constituents. We therefore test the hypothesis that variations in sediment lithology modulate RPI records. For 90 selected Late Quaternary sediment samples from the subtropical and subantarctic South Atlantic Ocean a combined paleomagnetic and sedimentological dataset was established. Misleading alterations of the magnetic mineral fraction were detected by a routine Fe/kappa test (Funk, J., von Dobeneck, T., Reitz, A., 2004. Integrated rock magnetic and geochemical quantification of redoxomorphic iron mineral diagenesis in Late Quaternary sediments from the Equatorial Atlantic. In: Wefer, G., Mulitza, S., Ratmeyer, V. (Eds.), The South Atlantic in the Late Quaternary: reconstruction of material budgets and current systems. Springer-Verlag, Berlin/Heidelberg/New York/Tokyo, pp. 239-262). Samples with any indication of suboxic magnetite dissolution were excluded from the dataset. The parameters under study include carbonate, opal and terrigenous content, grain size distribution and clay mineral composition. Their bi- and multivariate correlations with the RPI signal were statistically investigated using standard techniques and criteria. While several of the parameters did not yield significant results, clay grain size and chlorite correlate weakly and opal, illite and kaolinite correlate moderately to the NRM/ARM signal used here as a RPI measure. The most influential single sedimentological factor is the kaolinite/illite ratio with a Pearson's coefficient of 0.51 and 99.9% significance. A three-member regression model suggests that matrix effects can make up over 50% of the observed RPI dynamics.

Geochemical parameters of bottom sediments from the Tajura Rift, Gulf of Aden

Results of a geochemical study of bottom sediments from the Tadjura rift zone are reported. The sediments were analyzed for CaCO3, Si, Al, Ti, Fe, Mn, Cu, Zn, Ni, Co, Cr, V, Zr, Ga, Yb and Y. It was found that formation of chemical composition of the sediments was controlled by factors being appropriate for a near-continental area of the arid climatic zone (aeolian supply of terrigenous material and high biological productivity), as well as by hydrothermal activity in the rift valley. It was shown that high Mn contents were typical for the sediments in study while maximal contents of Fe were found near supposed hydrothermal sources. Total flux of Mn into sediments was been calculated. Diagenetic redistribution gives the main contribution of Mn in surface layer sediments. Speciations of Fe, Mn, Cu, Zn, Ni, Co, and Al were studied. In the surface layer sediments iron and manganese were in hydroxides. Model calculations of contents of chemical elements in sediments of the area in study are given.

(Table T1) Geochemical composition of ODP Site 177-1089 sediments

A primary objective of Ocean Drilling Program Leg 177 was to document changes in circulation and biogeochemical cycling on glacial/interglacial time scales across a wide latitudinal range of the south Atlantic Ocean. One of the more northerly sites drilled, Site 1089 (41°S, 10°E), is located within the present-day Subantarctic Zone, south of the Subtropical Front. The drilling site itself is located in the southern Cape Basin at a water depth of 4620 m. Pleistocene sediments at this site are dominated by interbedded carbonate and opal oozes. Initial shipboard stratigraphy identified the opal-rich sediments as deposited during glacial intervals and the carbonate-rich sediments as deposited during interglacial intervals (Gersonde, Hodell, Blum, et al., 1999, doi:10.2973/odp.proc.ir.177.1999). Postcruise isotopic stratigraphy, however, verified that this site displayed a Pacific Pleistocene sedimentation pattern with glacial intervals marked by high carbonate content (Hodell and Charles, 1999). To assess changes in biological productivity and terrigenous inputs at this site, a number of geochemical indicators were determined. Phosphorus concentrations and P/metal ratios were determined to assess changes in export production on glacial/interglacial time scales. Metal concentrations, along with elemental ratios, were used to assess terrigenous inputs. Sediment geochemistry allows us to identify changes in the lithologic component using elemental data based on Fe, Al, and Ti concentrations. Records of concentrations and ratios of biologically related elements identify changes in export production. The P and metal results are important to assess the glacial/interglacial changes in P burial and the relationships between a major nutrient such as P with metals (and possibly trace nutrients) like Fe.

Data-driven fuzzy analysis in quantitative mineral resource assessment

The integration of geo-information from multiple sources and of diverse nature in developing mineral favourability indexes (MFIs) is a well-known problem in mineral exploration and mineral resource assessment. Fuzzy set theory provides a convenient framework to combine and analyse qualitative and quantitative data independently of their source or characteristics. A novel, data-driven formulation for calculating MFIs based on fuzzy analysis is developed in this paper. Different geo-variables are considered fuzzy sets and their appropriate membership functions are defined and modelled. A new weighted average-type aggregation operator is then introduced to generate a new fuzzy set representing mineral favourability. The membership grades of the new fuzzy set are considered as the MFI. The weights for the aggregation operation combine the individual membership functions of the geo-variables, and are derived using information from training areas and L, regression. The technique is demonstrated in a case study of skarn tin deposits and is used to integrate geological, geochemical and magnetic data. The study area covers a total of 22.5 km(2) and is divided into 349 cells, which include nine control cells. Nine geo-variables are considered in this study. Depending on the nature of the various geo-variables, four different types of membership functions are used to model the fuzzy membership of the geo-variables involved. (C) 2002 Elsevier Science Ltd. All rights reserved.

Clay mineralogical, geochemical and isotopic tracing of the evolution of the Woodleigh impact structure, Southern Carnarvon Basin, Western Australia

Chaotically structured diamictite from the inner ring syncline surrounding the central uplift of the Woodleigh impact structure contains shocked metamorphic and impact melt-rock fragments, largely derived from Ordovician and Devonian target sandstones. Coarse illite fractions (< 2 mu m) from the sandstones containing no K-feldspar yield K-Ar ages of around 400 Ma, whereas the K-Ar ages of authigenic clays of > 0.2 mu m fractions from the diamictite without smectite and K-feldspar cluster around 360 Ma, consistent with Rb-Sr data. Crystallisation of newly formed illite in the impact melt rock clasts and recrystallisation of earlier formed illite in the sandstone clasts preserved in the diamictite, are attributed to impact-induced hydrothermal processes in the Late Devonian. The illitic clays from the diamictite and from the sandstones have very similar trace element compositions, with significantly enriched incompatible lithophile elements, which increase in concentrations correlatively with those of the compatible ferromagnesian elements. The unusual trace element associations in the clays may be due to the involvement of hot gravity-driven basinal fluids that interacted with rocks of the Precambrian craton to the east of the study area, or with such material transported and reworked in the studied sedimentary succession.

Geochemical and Sr, Nd, Pb isotope investigation of the New Caledonia ophiolite

The New Caledonia ophiolite hosts one of the largest obducted mantle section in the world, hence providing a unique insight for the study of upper mantle processes. These mantle rocks belong to an “atypical” ophiolitic sequence, which is dominated by refractory harzburgites but it also includes minor spinel and plagioclase lherzolites. Upper crust is notably absent in the ophiolite, with the exception of some mafic-ultramafic cumulates cropping out in the southern part of the island. Although the New Caledonia ophiolite has been under investigation for decades, its ultra-depleted nature has made its characterization an analytical challenge, so that few trace element data are available, while isotopic data are completely missing. In this thesis a comprehensive geochemical study (major, trace element and Sr-Nd-Pb isotopes) of the peridotites and the associated intrusive mafic rocks from the New Caledonia ophiolite has been carried out. The peridotites are low-strain tectonites showing porphyroclastic textures. Spinel lherzolites are undepleted lithotypes, as attested by the presence of 7-8 vol% of Na2O and Al2O3-rich clinopyroxene (up to 0.5 wt% Na2O; 6.5 wt% Al2O3), Fo content of olivine (88.5-90.0 mol%) and low Cr# of spinel (13-17). Conversely, harzburgites display a refractory nature, proven by the remarkable absence of primary clinopyroxene, very high Fo content in olivine (90.9-92.9 mol%), high Mg# in orthopyroxene (89.8-94.2) and Cr# in spinel (39-71). REE contents show abyssal-type patterns for spinel lherzolites, while harzburgites display U-shaped patterns, typical of fore-arc settings. Spinel lherzolites REE compositions are consistent with relatively low degree (8-9%) of fractional melting of a DMM source, starting in the garnet stability field. Conversely, REE models for harzburgites indicate high melting degrees (20-25%) of a DMM mantle source under spinel faies conditions, consistent with hydrous melting in forearc setting. Plagioclase lherzolites exhibit melt impregnation microtextures, Cr- and TiO2-enriched spinels and REE, Ti, Y, Zr progressive increase with respect to spinel lherzolites. Impregnation models indicate that plagioclase lherzolites may derive from spinel lherzolites by entrapment of highly depleted MORB melts in the shallow oceanic lithosphere. Mafic intrusives are olivine gabbronorites with a very refractory composition, as attested by high Fo content of olivine (87.3-88.9 mol.%), very high Mg# of clinopyroxene (87.7-92.2) and extreme anorthitic content of plagioclase (An = 90-96 mol%). The high Mg#, low TiO2 concentrations in pyroxenes and the anorthitic composition of plagioclase point out an origin from ultra-depleted primitive magmas in a convergent setting. Geochemical trace element models show that the parental melts of gabbronorites are primitive magmas with striking depleted compositions, bearing only in part similarities with the primitive boninitic melts of Bonin Islands. The first Sr, Nd and Pb isotope data obtained for the New Caledonia ophiolite highlight the presence of DM mantle source variably modified by different processes. Nd-Sr-Pb isotopic ratios for the lherzolites (+6.98≤epsilon Ndi≤+10.97) indicate a DM source that suffered low-temperature hydrothermal reactions. Harzburgites are characterized by a wide variation of Sr, Nd and Pb isotopic values, extending from DM-type to EM2 compositions (-0.82≤ epsilon Ndi≤+17.55), suggesting that harzburgite source was strongly affected by subduction-related processes. Conversely, combined trace element and Sr-Nd-Pb isotopic data for gabbronorites indicate a derivation from a source with composition similar to Indian-type mantle, but affected by fluid input in subduction environment. These geochemical features point out an evolution in a pre-Eocenic marginal basin setting, possibly in the proximity of a transform fault, for the lherzolites. Conversely, the harzburgites acquired their main geochemical and isotopic fingerprint in subduction zone setting.

Data visualisation and exploration with prior knowledge

Visualising data for exploratory analysis is a major challenge in many applications. Visualisation allows scientists to gain insight into the structure and distribution of the data, for example finding common patterns and relationships between samples as well as variables. Typically, visualisation methods like principal component analysis and multi-dimensional scaling are employed. These methods are favoured because of their simplicity, but they cannot cope with missing data and it is difficult to incorporate prior knowledge about properties of the variable space into the analysis; this is particularly important in the high-dimensional, sparse datasets typical in geochemistry. In this paper we show how to utilise a block-structured correlation matrix using a modification of a well known non-linear probabilistic visualisation model, the Generative Topographic Mapping (GTM), which can cope with missing data. The block structure supports direct modelling of strongly correlated variables. We show that including prior structural information it is possible to improve both the data visualisation and the model fit. These benefits are demonstrated on artificial data as well as a real geochemical dataset used for oil exploration, where the proposed modifications improved the missing data imputation results by 3 to 13%.

Exploratory data analysis with non-linear and missing data in geochemistry

Exploratory analysis of data seeks to find common patterns to gain insights into the structure and distribution of the data. In geochemistry it is a valuable means to gain insights into the complicated processes making up a petroleum system. Typically linear visualisation methods like principal components analysis, linked plots, or brushing are used. These methods can not directly be employed when dealing with missing data and they struggle to capture global non-linear structures in the data, however they can do so locally. This thesis discusses a complementary approach based on a non-linear probabilistic model. The generative topographic mapping (GTM) enables the visualisation of the effects of very many variables on a single plot, which is able to incorporate more structure than a two dimensional principal components plot. The model can deal with uncertainty, missing data and allows for the exploration of the non-linear structure in the data. In this thesis a novel approach to initialise the GTM with arbitrary projections is developed. This makes it possible to combine GTM with algorithms like Isomap and fit complex non-linear structure like the Swiss-roll. Another novel extension is the incorporation of prior knowledge about the structure of the covariance matrix. This extension greatly enhances the modelling capabilities of the algorithm resulting in better fit to the data and better imputation capabilities for missing data. Additionally an extensive benchmark study of the missing data imputation capabilities of GTM is performed. Further a novel approach, based on missing data, will be introduced to benchmark the fit of probabilistic visualisation algorithms on unlabelled data. Finally the work is complemented by evaluating the algorithms on real-life datasets from geochemical projects.

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.