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.

Soil Carbon stocks and schanges in the Republic of Ireland

The soil carbon (C) stock of the Republic of Ireland is estimated to have been 2048 Mt in 1990 and 2021 Mt in 2000. Peat holds around 53% of the soil C stock, but on 17% of the land area. The C density of soils (t C ha-1) is mapped at 2 km*2 km resolution. The greatest soil C densities occur where deep raised bogs are the dominant soil; in these grid squares C density can reach 3000 t C ha-1. Most of the loss of soil C between 1990 and 2000-up to 23 Mt C (1% of 1990 soil C stock)-was through industrial peat extraction. The average annual change in soil C stocks from 1990 to 2000 due to land use change was estimated at around 0.02% of the 1990 stock. Considering uncertainties in the data used to calculate soil C stocks and changes, the small average annual 'loss' could be regarded as 'no change'.

Deep weathering of calcareous sedimentary rock and the redistribution of iron and manganese in soil and saprolite

Iron and Mn redistribute in soil and saprolite during weathering. The geological weathering fronts ofcalcareous sedimentary rock were investigated by examining the bulk density, porosity, and distribution ofCa, Fe, and Mn. Core samples were taken ofsoil, saprolite, and bedrock material from both summit (HHMS-4B) and sideslope (HHMS-5A) positions on an interbedded Nolichucky shale and Maryville limestone landform in Solid Waste Storage Area 6 (SWSA-6). This is a low-level radioactive solids waste disposal site on the Dept. ofEnergy (DOE) Oak Ridge Reservation in Roane County Tennessee. This work was initiated because data about the properties of highly weathered sedimentary rock on this site were limited. The core samples were analyzed for pH, calcium carbonate equivalence (CCE), hydroxylamine-extractable (HA) Mn, and dithionite-citrate (CBD)-extractable Fe and Mn. Low pH values occurred from the soil surface down to the depth of the oxidized and leached saprolite in both cores. The CCE and HA-extractable Mn results were also influenced by the weathering that has occurred in these zones. Extractable Mn oxide was higher at a lower depth in the oxidized and leached saprolite compared with the Fe oxide, which was higher in the overlying soil solum. Amounts of Mn oxides were higher in the sideslope core (HHMS-5A) than in the summit core (HHMS-4B). Iron was more abundant in the deeper weathered summit core, but the highest value, 39.4 g kg-1, was found at 1.8 to 2.4 m in the sideslope core. The zone encompassing the oxidized and partially leached saprolite down to the unoxidized and unleached bedrock had higher densities and larger quantities of CaCO3 than the soil solum and oxidized and leached saprolite. The overlying soil and oxidized and leached saprolite had lower pH and CCE values and were higher in Fe and Mn oxides than the oxidized and unleached saprolite. The distribution of Fe and Mn is important when evaluating soil and saprolite for hazardous waste disposal site assessment.

High density mangrove plantation enhances surface accretion, surface elevation change, and tree survival in coastal areas susceptible to sea-level rise

Survival, growth, above ground biomass accumulation, soil surface elevation dynamics and nitrogen accumulation in accreted sediments were studied in experimental treatments planted with four different densities (6.96, 3.26, 1.93 and 0.95 seedlings m-2) of the mangrove Rhizophora mucronata in Puttalam Lagoon, Sri Lanka. Measurements were taken over a period of 1171 days and were compared with those from unplanted controls. Trees at the lowest density showed significantly reduced survival, whilst measures of individual tree growth did not differ significantly among treatments. Rates of surface sediment accretion (means ± S.E.) were 13.0 (±1.3), 10.5 (±0.9), 8.4 (±0.3), 6.9 (±0.5) and 5.7 (±0.3) mm yr-1 at planting densities of 6.96, 3.26, 1.93, 0.95, and 0 (unplanted control) seedlings m-2, respectively, showing highly significant differences among treatments. Mean (± S.E.) rates of surface elevation change were much lower than rates of accretion at 2.8 (±0.2), 1.6 (±0.1), 1.1 (±0.2), 0.6 (±0.2) and -0.3 (±0.1) mm yr-1 for 6.96, 3.26, 1.93, 0.95, and 0 seedlings m-2, respectively. All planted treatments appeared to accumulate greater nitrogen concentrations in the sediment compared to the unplanted control, and suggests one potential causal mechanism for the facilitatory effects observed; high densities of plants potentially contribute to the accretion of greater amounts of nutrient rich sediment. While this potential process needs further study, this study demonstrated how higher densities of mangroves enhance rates of sediment accretion and surface elevation, processes that may be crucial in mangrove ecosystem adaptation to sea level rise. There was no evidence that increasing plant density evoked a trade-off with growth and survival of the planted trees. Rather facilitatory effects enhanced survival at high densities, suggesting that local land managers may be able to take advantage of plantation densities to help mitigate sea-level rise effects by encouraging positive soil surface elevation increment, and perhaps even greater nutrient retention to promote mangrove growth and ameliorate nearshore eutrophication in tropical island environments.

Detection of soil compaction using seismic surface waves

Seismic geophysical methods have rarely been used in precision agriculture, predominantly due to the perception that they are slow and results require a complex evaluation. This paper explores the possibility of using a recently developed surface wave seismic geophysical approach, the multichannel analysis of surface waves (MASW) method, for assessment of agricultural compaction. This approach has the advantage of being non-intrusive, rapid and is able to produce 2D ground models with a relatively high density of spatial sampling points. The method, which was tested on a research site in Oakpark, Ireland, detected a significant difference in shear wave velocity between a heavily compacted headland and an uncompacted location. The results from this approach compared favourably with those obtained
from measurements of bulk density and penetrometer resistance and demonstrate that the MASW approach can distinguish between the extreme states of heavily compacted and uncompacted soil.

Stable isotope probing analysis of the influence of liming on root exudate utilization by soil microorganisms

Rhizosphere microorganisms play an important role in soil carbon flow, through turnover of root exudates, but there is little information on which organisms are actively involved or on the influence of environmental conditions on active communities. In this study, a (CO2)-C-13 pulse labelling field experiment was performed in an upland grassland soil, followed by RNA-stable isotope probing (SIP) analysis, to determine the effect of liming on the structure of the rhizosphere microbial community metabolizing root exudates. The lower limit of detection for SIP was determined in soil samples inoculated with a range of concentrations of C-13-labelled Pseudomonas fluorescens and was found to lie between 10(5) and 10(6) cells per gram of soil. The technique was capable of detecting microbial communities actively assimilating root exudates derived from recent photo-assimilate in the field. Denaturing gradient gel electrophoresis (DGGE) profiles of bacteria, archaea and fungi derived from fractions obtained from caesium trifluoroacetate (CsTFA) density gradient ultracentrifugation indicated that active communities in limed soils were more complex than those in unlimed soils and were more active in utilization of recently exuded C-13 compounds. In limed soils, the majority of the community detected by standard RNA-DGGE analysis appeared to be utilizing root exudates. In unlimed soils, DGGE profiles from C-12 and C-13 RNA fractions differed, suggesting that a proportion of the active community was utilizing other sources of organic carbon. These differences may reflect differences in the amount of root exudation under the different conditions.

Direct, positive feedbacks produce instability in models of interrelationships among soil structure, plants and arbuscular mycorrhizal fungi

Current conceptual models of reciprocal interactions linking soil structure, plants and arbuscular mycorrhizal fungi emphasise positive feedbacks among the components of the system. However, dynamical systems with high dimensionality and several positive feedbacks (i.e. mutualism) are prone to instability. Further, organisms such as arbuscular mycorrhizal fungi (AMF) are obligate biotrophs of plants and are considered major biological agents in soil aggregate stabilization. With these considerations in mind, we developed dynamical models of soil ecosystems that reflect the main features of current conceptual models and empirical data, especially positive feedbacks and linear interactions among plants, AMF and the component of soil structure dependent on aggregates. We found that systems become increasingly unstable the more positive effects with Type I functional response (i.e., the growth rate of a mutualist is modified by the density of its partner through linear proportionality) are added to the model, to the point that increasing the realism of models by adding linear effects produces the most unstable systems. The present theoretical analysis thus offers a framework for modelling and suggests new directions for experimental studies on the interrelationship between soil structure, plants and AMF. Non-linearity in functional responses, spatial and temporal heterogeneity, and indirect effects can be invoked on a theoretical basis and experimentally tested in laboratory and field experiments in order to account for and buffer the local instability of the simplest of current scenarios. This first model presented here may generate interest in more explicitly representing the role of biota in soil physical structure, a phenomenon that is typically viewed in a more process- and management-focused context. (C) 2011 Elsevier Ltd. All rights reserved.

Assessing abundance and diversity patterns of soil microarthropod assemblages in northern Victoria Land (Antarctica)

This study aimed to: (1) assess differences between two quantitative sampling methods of soil microarthropods (visual census vs. stone washing) in ice-free areas located along a latitudinal gradient (from 72 degrees 37'S to 74 degrees 42'S) in northern Victoria Land (Antarctica); (2) furnish preliminary results on the abundance and diversity of mites and springtails in the studied areas. Visual census yielded reliable density estimates for adult collembolans and larger prostigmatic mites but did not detect small species. The study updates the distribution of several mites, including the southernmost record of an Oribatida species at global scale. Species composition was correlated with latitude but the uneven abundance distribution and local high beta-diversity probably reflect habitat fragmentation and population isolation. Under this circumstance nested sampling design should be usefully employed. Priorities and suitable methods for studying terrestrial microarthropod communities in continental Antarctica are discussed.

Using geostatistical Bayesian Updating to integrate airborne radiometrics and soil geochemistry to improve mapping for mineral exploration

Mineral exploration programmes around the world use data from remote sensing, geophysics and direct sampling. On a regional scale, the combination of airborne geophysics and ground-based geochemical sampling can aid geological mapping and economic minerals exploration. The fact that airborne geophysical and traditional soil-sampling data are generated at different spatial resolutions means that they are not immediately comparable due to their different sampling density. Several geostatistical techniques, including indicator cokriging and collocated cokriging, can be used to integrate different types of data into a geostatistical model. With increasing numbers of variables the inference of the cross-covariance model required for cokriging can be demanding in terms of effort and computational time. In this paper a Gaussian-based Bayesian updating approach is applied to integrate airborne radiometric data and ground-sampled geochemical soil data to maximise information generated from the soil survey, to enable more accurate geological interpretation for the exploration and development of natural resources. The Bayesian updating technique decomposes the collocated estimate into a production of two models: prior and likelihood models. The prior model is built from primary information and the likelihood model is built from secondary information. The prior model is then updated with the likelihood model to build the final model. The approach allows multiple secondary variables to be simultaneously integrated into the mapping of the primary variable. The Bayesian updating approach is demonstrated using a case study from Northern Ireland where the history of mineral prospecting for precious and base metals dates from the 18th century. Vein-hosted, strata-bound and volcanogenic occurrences of mineralisation are found. The geostatistical technique was used to improve the resolution of soil geochemistry, collected one sample per 2 km2, by integrating more closely measured airborne geophysical data from the GSNI Tellus Survey, measured over a footprint of 65 x 200 m. The directly measured geochemistry data were considered as primary data in the Bayesian approach and the airborne radiometric data were used as secondary data. The approach produced more detailed updated maps and in particular maximized information on mapped estimates of zinc, copper and lead. Greater delineation of an elongated northwest/southeast trending zone in the updated maps strengthened the potential to investigate stratabound base metal deposits.

Spatial and Temporal Changes in Soil Organic Carbon Assessment using an Integrated Geostatistical and Compositional Data Analysis Approach

The Irish and UK governments, along with other countries, have made a commitment to limit the concentrations of greenhouse gases in the atmosphere by reducing emissions from the burning of fossil fuels. This can be achieved (in part) through increasing the sequestration of CO2 from the atmosphere including monitoring the amount stored in vegetation and soils. A large proportion of soil carbon is held within peat due to the relatively high carbon density of peat and organic-rich soils. This is particularly important for a country such as Ireland, where some 16% of the land surface is covered by peat. For Northern Ireland, it has been estimated that the total amount of carbon stored in vegetation is 4.4Mt compared to 386Mt stored within peat and soils. As a result it has become increasingly important to measure and monitor changes in stores of carbon in soils. The conservation and restoration of peat covered areas, although ongoing for many years, has become increasingly important. This is summed up in current EU policy outlined by the European Commission (2012) which seeks to assess the relative contributions of the different inputs and outputs of organic carbon and organic matter to and from soil. Results are presented from the EU-funded Tellus Border Soil Carbon Project (2011 to 2013) which aimed to improve current estimates of carbon in soil and peat across Northern Ireland and the bordering counties of the Republic of Ireland.
Historical reports and previous surveys provide baseline data. To monitor change in peat depth and soil organic carbon, these historical data are integrated with more recently acquired airborne geophysical (radiometric) data and ground-based geochemical data generated by two surveys, the Tellus Project (2004-2007: covering Northern Ireland) and the EU-funded Tellus Border project (2011-2013) covering the six bordering counties of the Republic of Ireland, Donegal, Sligo, Leitrim, Cavan, Monaghan and Louth. The concept being applied is that saturated organic-rich soil and peat attenuate gamma-radiation from underlying soils and rocks. This research uses the degree of spatial correlation (coregionalization) between peat depth, soil organic carbon (SOC) and the attenuation of the radiometric signal to update a limited sampling regime of ground-based measurements with remotely acquired data. To comply with the compositional nature of the SOC data (perturbations of loss on ignition [LOI] data), a compositional data analysis approach is investigated. Contemporaneous ground-based measurements allow corroboration for the updated mapped outputs. This provides a methodology that can be used to improve estimates of soil carbon with minimal impact to sensitive habitats (like peat bogs), but with maximum output of data and knowledge.

Diversity and abundance of soil arthropods in urban and suburban holm oak stands

We investigated the soil arthropod communities of urban and suburban holm oak (Quercus ilex L.) stands in a small (Siena) and a large Italian city (Naples) and tested whether the abundance and diversity of higher arthropod taxa are affected by the biotic and abiotic conditions of urban forest soils, including pollution. Acarina and Collembola were the dominant taxa in both cities. In Siena the total number of arthropod individuals collected in the samples was over 1/3 greater than in Naples, but all diversity indices scored higher in Naples than in Siena, probably in response to the higher heterogeneity of microclimatic and pedological conditions found in Naples study area. Oribatids resulted twice more abundant in Siena and so were the total mites with respect to Collembola. While “taxonomic richness” per site increased with distance from road traffic, entropy and evenness indices scored higher at the two ends of the impact gradient in both cities. The overall variation in basic pedological and microbiological soil parameters positively correlated with the total abundance of arthropods, and negatively correlated with their taxonomic richness. At the resolution employed, no significant relation emerged between anthropogenic factors, such as traffic load and soil pollution, and the arthropod fauna density and variety. These results are consistent with conclusions drawn from a previous study on the enchytraeid fauna examined at species level, which is remarkable considering the different taxonomic resolutions of the two studies. CCA results suggest that the higher abundance of Oribatid mites, Protura and Thysanura and the lower abundance of Diplopoda and Symphyla in Siena could depend on a higher fungi/bacteria ratio. This observation can be interpreted in terms of differences in fungi and bacteria between the two cities: Siena is shifted towards the fungal decomposition channel, which supports taxa such as oribatid mites, while Naples is shifted towards the bacterial channel, which supports chiefly detritivorous groups, such as diplopods.

Toxic effects of cadmium and zinc on ectomycorrhizal colonization of scots pine (Pinus sylvestris L.) from soil inoculum

Scots pine seedlings colonized by ectomycorrhizal (ECM) fungi from natural soil inoculum were exposed to a range of Cd or Zn concentrations to investigate the effects of metals on ECM fungi-Scots pine associations in a realistic soil environment. Experiments focused on the relationship between the sensitivity of ECM fungi and their host plants, the influence of metals on ECM community dynamics on Scots pine roots, and the effects of metal exposure on ECM colonization from soil-borne propagules. Ectomycorrhizal colonization was inhibited by Cd and Zn, with a decrease in the proportion of ECM-colonized root tips. Shoot and root biomass, total root length, and total root-tip density, however, were unaffected by Cd or Zn. A decrease in the diversity of ECM morphotypes also occurred, which could have a negative effect on tree vigor. Overall, colonization by ECM fungi was more sensitive than seedling growth to Cd and Zn, and this could have serious implications for successful tree establishment on metal-contaminated soils.

Catchment and in-stream influences on metal concentration and ochre deposit density in upland streams, Northern Ireland

Metal concentrations from stream waters in two geological blocks in Northern Ireland were compared to determine the contributions of catchment characteristics and in-stream conditions. One block is composed of metamorphosed schist and unconsolidated glacial drift with peat or peaty podzol (mainly humic) soils, while the other block consists of tertiary basalt with brown earth and gley soils. Water samples were collected from 52 stream sites and analysed for Fe, Mn and Al as well as a range of other chemical determinands known to affect metal solubility. Densities of metal-rich ochre deposit were determined for stream bed stone samples. Higher conductivities and concentrations of bicarbonate, alkalinity, Ca and Mg occurred on basalt than on schist. Despite higher Fe and Mn oxide concentrations in basalt-derived non-humic soils, stream water concentrations were much lower and ochre deposit densities only one third of those on schist overlain by humic soils. Neither rock nor soil type predicted Al concentrations, but pH and dissolved oxygen did. Peat-generated acidity and the limited acid neutralising capacity of base-poor metamorphosed schist have resulted in elevated concentrations of metals and ochre deposit in surface waters.