Compositional data-driven alternatives to single geochemical component maps

Conventional practice in Regional Geochemistry includes as a final step of any geochemical campaign the generation of a series of maps, to show the spatial distribution of each of the components considered. Such maps, though necessary, do not comply with the compositional, relative nature of the data, which unfortunately make any conclusion based on them sensitive
to spurious correlation problems. This is one of the reasons why these maps are never interpreted isolated. This contribution aims at gathering a series of statistical methods to produce individual maps of multiplicative combinations of components (logcontrasts), much in the flavor of equilibrium constants, which are designed on purpose to capture certain aspects of the data.
We distinguish between supervised and unsupervised methods, where the first require an external, non-compositional variable (besides the compositional geochemical information) available in an analogous training set. This external variable can be a quantity (soil density, collocated magnetics, collocated ratio of Th/U spectral gamma counts, proportion of clay particle fraction, etc) or a category (rock type, land use type, etc). In the supervised methods, a regression-like model between the external variable and the geochemical composition is derived in the training set, and then this model is mapped on the whole region. This case is illustrated with the Tellus dataset, covering Northern Ireland at a density of 1 soil sample per 2 square km, where we map the presence of blanket peat and the underlying geology. The unsupervised methods considered include principal components and principal balances
(Pawlowsky-Glahn et al., CoDaWork2013), i.e. logcontrasts of the data that are devised to capture very large variability or else be quasi-constant. Using the Tellus dataset again, it is found that geological features are highlighted by the quasi-constant ratios Hf/Nb and their ratio against SiO2; Rb/K2O and Zr/Na2O and the balance between these two groups of two variables; the balance of Al2O3 and TiO2 vs. MgO; or the balance of Cr, Ni and Co vs. V and Fe2O3. The largest variability appears to be related to the presence/absence of peat.

What defines a geochemical map? Implications for environmental assessment.

The environmental quality of land is often assessed by the calculation of threshold values which aim to differentiate between concentrations of elements based on whether the soils are in residential or industrial sites. In Europe, for example, soil guideline values exist for agricultural and grazing land. A threshold is often set to differentiate between concentrations of the element that naturally occur in the soil and concentrations that result from diffuse anthropogenic sources. Regional geochemistry and, in particular, single component geochemical maps are increasingly being used to determine these baseline environmental assessments. The key question raised in this paper is whether the geochemical map can provide an accurate interpretation on its own. Implicit is the thought that single component geochemical maps represent absolute abundances. However,because of the compositional (closed) nature of the data univariate geochemical maps cannot be compared directly with one another.. As a result, any interpretation based on them is vulnerable to spurious correlation problems. What does this mean for soil geochemistry mapping, baseline quality documentation, soil resource assessment or risk evaluation? Despite the limitation of relative abundances, individual raw geochemical maps are deemed fundamental to several applications of geochemical maps including environmental assessments. However, element toxicity is related to its bioavailable concentration, which is lowered if its source is mixed with another source. Elements interact, for example under reducing conditions with iron oxides, its solid state is lost and arsenic becomes soluble and mobile. Both of these matters may be more adequately dealt with if a single component map is not interpreted in isolation to determine baseline and threshold assessments. A range of alternative compositionally compliant representations based on log-ratio and log-contrast approaches are explored to supplement the classical single component maps for environmental assessment. Case study examples are shown based on the Tellus soil geochemical dataset, covering Northern Ireland and the results of in vitro oral bioaccessibility testing carried out on a sub-set of archived Tellus Survey shallow soils following the Unified BARGE (Bioaccessibility Research Group of Europe).

Geochemical study of the rehabilitation of a shallow anoxic reservoir, in volcanic terrains, under semi-arid climate, in Cape Verde Archipelago

Poilão dam reservoir (Cape Verde Archipelago) is in critical conditions, owing the excessive silting up, the high concentration of nitrates and the pronounced anoxia all over the lake. Considering that the most suitable remediation strategy is the removal of the bottom sediments where nutrients are preferentially concentrated, we have done a geochemical study, in order of evaluating their suitability to agricultural use. Analyses indicate that sediments are rich in a few key nutrients, when compared with parent soils. Thus, adding suitable sediments to nearby degraded soils can improve food crops for smallholder farmers living in close proximity to this system.

Geochemical analysis of sediments from a polluted waterway in the vicinity of a metallurgical plant, Minas Gerais, Brazil

Mining and metallurgical industries are the most responsible for heavy metal contamination. These contaminants are often associated to poor management strategies and the lack of suitable containment areas. This is the case of one of the largest metal alloy producers in the São Francisco river basin, in Minas Gerais state, Brazil. During operation, the tailings of this activity have been accumulated in the surrounding area, leading to heavy metal contamination (Mn, Cu, Zn, Ni, Pb, Cd, Cr and As), much above the critical levels, which accumulate in the sediments of the nearby waterways, in the soils of the drainage area, and in the water column. In this work we are going to discuss both the contamination level and its extent in the Consciência river, a small tributary of the São Francisco river, the longest river entirely Brazilian, by analyzing the geochemistry and mineralogy of the accumulated sediments. It will be also discuss the reclamation strategies of this site and the impact on the environment.

Petrologic and geochemical characterization of São Jorge island volcanism, Azores

The island of São Jorge (38º 45’ 24’’ N - 28º 20’ 44’’W and 38º 33’ 00’’ N - 27º 44’ 32’’ W) is one of the nine islands of the Azores Archipelago that is rooted in the Azores Plateau, a wide and complex region which encompasses the triple junction between the American, Eurasia and Nubia plates. São Jorge Island has grown by fissural volcanic activity along fractures with the regional WNW-ESE trend, unveiling the importance of the regional tectonics during volcanic activity. The combination of the volcanostratigraphy (Forjaz & Fernandes, 1975; and Madeira, 1998) with geochronological data evidences that the island developed during two main volcanic phases. The first subaerial phase that occurred between 1.32 and 1.21 Ma ago (Hildenbrand et al. 2008) is recorded on the lava sequence forming the cliff at Fajã de São João, while the second phase started at 757 ka ago, is still active, and edified the rest of the island. This second phase edified the east side of the island that corresponds to Topo Volcanic Complex, in the period between 757 and 543 ka ago, while the west side named Rosais Volcanic Complex, started at 368 ka ago (Hildenbrand et al. 2008) and was still active at 117 ka ago. After the onset of Rosais, volcanic activity migrates to the center of São Jorge edifying Manadas Volcanic Complex. The volcanism on São Jorge is dominantly alkaline, with a narrow lithological composition ranging between the basanites/tefrites through the basaltic trachyandesites, in spite of this the two volcanic phases show distinct mineralogical, petrographic and geochemical characteristics that should be related with different petrogenetic conditions and growth rates of the island. Abstract viii During the first volcanic phase, growth rates are faster (≈3.4 m/ka), the lavas are slightly less alkaline and plagioclase-richer, pointing to the existence of a relative shallow and dynamic magma chamber where fractional crystallization associated with gravitational segregation and accumulation processes, produced the lavas of Fajã de São João sequence. The average growth rates during the second volcanic phase are lower (≈1.9 m/ka) and the lavas are mainly alkaline sodic, with a mineralogy composed by olivine, pyroxene, plagioclase and oxide phenocrysts, in a crystalline groundmass. The lavas are characterized by enrichment in incompatible trace element and light REE, but show differences for close-spaced lavas that unveil, in some cases, slight different degrees of fertilization of the mantle source along the island. These differences might also result from higher degrees of partial melting, as observed in the early stages of Topo and Rosais volcanic complexes, of a mantle source with residual garnet and amphibole, and/or from changing melting conditions of the mantle source as pressure. The subtle geochemical differences of the lavas contrast with the isotopic signatures, obtained from Sr-Nd-Pb-Hf isotopes, that São Jorge Island volcanism exhibit along its volcanic complexes. The lavas from Topo Volcanic Complex and from the submarine flank, i.e. the lavas located east of Ribeira Seca Fault, sample a mantle source with similar isotopic signature that, in terms of lead, overlaps Terceira Island. The lavas from Rosais and Manadas volcanic complexes, the western lavas, sample a mantle source that becomes progressively more distinct towards the west end of the island and that, in terms of lead isotopes, trends towards the isotopic composition of Faial Island. The two isotopic signatures of São Jorge, observed from the combination of lead isotopes with the other three systems, seem to result from the mixing of three distinct end-members. These end-members are (1) the common component related with the Azores Plateau and the MAR, (2) the eastern component with a FOZO signature and possibly related with the Azores plume located beneath Terceira, and (3) the western component, similar to Faial, where the lithosphere could have been entrained by an ancient magmatic liquid, isolated for a period longer than 2Ga. The two trends observed in the island reinforce the idea of small-scale mantle heterogeneities beneath the Azores region, as it has been proposed to explain the isotopic diversity observed in the Archipelago.

Integration of geochemical and biological modifications in mining areas: waters, sediments and diatoms

The development of mining activities over thousands of years in the region of Aljustrel is nowadays visible as a vast area of ore tailings, slag and host rocks of sulphides mineralization. The generation of acidic waters by the alteration of pyritic minerals - Acid Mine Drainage (AMD) - causes a significant impact on the river system both in the south of the village (Rib ª. Água Forte) and in the north of it (Rib ª. Água Azeda and Barranco do Farrobo), which is reflected in extremely low pH values (< 3) and high concentrations of As, Cd, Cu, Fe, Mn, Pb, Zn and sulphates. This study aimed to assess the environmental impacts extent, integrating geochemical (surface waters and stream sediments) and biological (diatoms) parameters. Three groups of sites were defined, based on sediments and water analysis, which integration with diatom data showed the same association of groups: Group 1- impacted, with acidic pH (1.9-5.1), high metal contents (0.4-1975 mg L-1) and Fe-Mg-sulphate waters, being metals more bioavailable in waters in cationic form (Me2+); mineralogically the sediments were characterized by phyllosilicates and sulphates/oxy-hydroxysulphate phases, easily solubilized, retaining a high amount of metals when precipitated; dominant taxon was Pinnularia aljustrelica (a new species); Group 2- slightly impacted, weak acid to neutral pH (5.0-6.8), metal contents not so high (0.2-25 mg L-1) and Fe-Mg-sulphate to Mg-chloride waters; dominant taxa were Brachysira neglectissima and Achnanthidium minutissimum; Group 3- unimpacted, alkaline pH (7.0-8.4), low metal contents (0-7 mg L-1) with Mg-chloride waters. In this group, metals were associated to the primary phases (e.g. sulphides), not so easily available; the existence of high chloride contents explained the presence of typical taxa of brackish/marine (e.g. Entomoneis paludosa) waters. Taxonomical aspects of the diatoms were studied (discovery of a new species: Pinnularia aljustrelica Luis, Almeida et Ector sp. nov.), as well as morphometric (size decrease of diatoms valves, as well as the appearance of deformed valves of Eunotia exigua in Group 1 and A. minutissimum in Group 2) and physiological (effective to assess the effects of metals/acidity in the photosynthetic efficiency through PAM Fluorometry) aspects. A study was carried out in an artificial river system (microcosm) that aimed to mimic Aljustrel’s extreme conditions in controlled laboratory conditions. The chronic effects of Fe, SO42- and acidity in field biofilms, inoculated in the artificial rivers, were evaluated as well as their contribution to the communities’ tolerance to metal toxicity, through acute tests with two metals (Cu and Zn). In general, the effects caused by low pH values and high concentrations of Fe and SO42- were reflected at the community level by the decrease in diversity, the predominance of acidophilic species, the decrease in photosynthetic efficiency and the increase of enzymatic (e.g. catalase, superoxide dismutase) and non-enzymatic activities (e.g. total glutathione and total phytochelatins). However, it was possible to verify that acidity performed a protective effect in the communities, upon Cu and Zn addition. A comparative study between Aljustrel mining area and New Brunswick mining area was carried out, both with similar mining and geological conditions, reflected in similar diatom communities in both mines, but in very different geographic and climatic areas.

Comparing Zostera and Spartina environments in relation to carbon burial: a sedimentary and geochemical approach from Ria Formosa

Dissertação de mestrado, Biologia Marinha, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

Impacts from aquaculture facilities on intertidal ecosystems using quadrant surveys and antibiotic tracers

Senior thesis written for Oceanography 444

Accuracy assessment of depth estimates from variation in pre- and post-processing of multi-beam SONAR surveys in a fjord environment

Senior thesis written for Oceanography 445

pH Limitations in environmental geochemical monitoring

Senior thesis written for Oceanography 445

Conditional probability applied to the definition of limestone geochemical facies: A new Methodological Approach

A new method, based on linear correlation and phase diagrams was successfully developed for processes like the sedimentary process, where the deposition phase can have different time duration - represented by repeated values in a series - and where the erosion can play an important rule deleting values of a series. The sampling process itself can be the cause of repeated values - large strata twice sampled - or deleted values: tiny strata fitted between two consecutive samples. What we developed was a mathematical procedure which, based upon the depth chemical composition evolution, allows the establishment of frontiers as well as the periodicity of different sedimentary environments. The basic tool isn't more than a linear correlation analysis which allow us to detect the existence of eventual evolution rules, connected with cyclical phenomena within time series (considering the space assimilated to time), with the final objective of prevision. A very interesting discovery was the phenomenon of repeated sliding windows that represent quasi-cycles of a series of quasi-periods. An accurate forecast can be obtained if we are inside a quasi-cycle (it is possible to predict the other elements of the cycle with the probability related with the number of repeated and deleted points). We deal with an innovator methodology, reason why it's efficiency is being tested in some case studies, with remarkable results that shows it's efficacy. Keywords: sedimentary environments, sequence stratigraphy, data analysis, time-series, conditional probability.

Optimization of River Water Quality Surveys by Multivariate Analysis of Physicochemical, Bacteriological and Ecotoxicological Data

This study aims to optimize the water quality monitoring of a polluted watercourse (Leça River, Portugal) through the principal component analysis (PCA) and cluster analysis (CA). These statistical methodologies were applied to physicochemical, bacteriological and ecotoxicological data (with the marine bacterium Vibrio fischeri and the green alga Chlorella vulgaris) obtained with the analysis of water samples monthly collected at seven monitoring sites and during five campaigns (February, May, June, August, and September 2006). The results of some variables were assigned to water quality classes according to national guidelines. Chemical and bacteriological quality data led to classify Leça River water quality as “bad” or “very bad”. PCA and CA identified monitoring sites with similar pollution pattern, giving to site 1 (located in the upstream stretch of the river) a distinct feature from all other sampling sites downstream. Ecotoxicity results corroborated this classification thus revealing differences in space and time. The present study includes not only physical, chemical and bacteriological but also ecotoxicological parameters, which broadens new perspectives in river water characterization. Moreover, the application of PCA and CA is very useful to optimize water quality monitoring networks, defining the minimum number of sites and their location. Thus, these tools can support appropriate management decisions.