Holocene tephras highlight complexity of volcanic signals in Greenland ice cores

Acidity peaks in Greenland ice cores have been used as critical reference horizons for synchronizing ice core records, aiding the construction of a single Greenland Ice Core Chronology (GICC05) for the Holocene. Guided by GICC05, we examined sub-sections of three Greenland cores in the search for tephra from specific eruptions that might facilitate the linkage of ice core records, the dating of prehistoric tephras and the understanding of the eruptions. Here we report the identification of 14 horizons with tephra particles, including 11 that have not previously been reported from the North Atlantic region and that have the potential to be valuable isochrons. The positions of tephras whose major element data are consistent with ash from the Katmai AD 1912 and Öraefajökull AD 1362 eruptions confirm the annually resolved ice core chronology for the last 700 years. We provide a more refined date for the so-called “AD860B” tephra, a widespread isochron found across NW Europe, and present new evidence relating to the 17th century BC Thera/Aniakchak debate that shows N. American eruptions likely contributed to the acid signals at this time. Our results emphasize the variable spatial and temporal distributions of volcanic products in Greenland ice that call for a more cautious approach in the attribution of acid signals to specific eruptive events.

Delineation of a quick clay zone at Smørgrav, Norway, with electromagnetic methods under geotechnical constraints

In many coastal areas of North America and Scandinavia, post-glacial clay sediments have emerged above sea level due to iso-static uplift. These clays are often destabilised by fresh water leaching and transformed to so-called quick clays as at the investigated area at Smørgrav, Norway. Slight mechanical disturbances of these materials may trigger landslides. Since the leaching increases the electrical resistivity of quick clay as compared to normal marine clay, the application of electromagnetic (EM) methods is of particular interest in the study of quick clay structures.

For the first time, single and joint inversions of direct-current resistivity (DCR), radiomagnetotelluric (RMT) and controlled-source audiomagnetotelluric (CSAMT) data were applied to delineate a zone of quick clay. The resulting 2-D models of electrical resistivity correlate excellently with previously published data from a ground conductivity metre and resistivity logs from two resistivity cone penetration tests (RCPT) into marine clay and quick clay. The RCPT log into the central part of the quick clay identifies the electrical resistivity of the quick clay structure to lie between 10 and 80 O m. In combination with the 2-D inversion models, it becomes possible to delineate the vertical and horizontal extent of the quick clay zone. As compared to the inversions of single data sets, the joint inversion model exhibits sharper resistivity contrasts and its resistivity values are more characteristic of the expected geology. In our preferred joint inversion model, there is a clear demarcation between dry soil, marine clay, quick clay and bedrock, which consists of alum shale and limestone.

Hydrobiogeochemical Interactions along “Flow Tubes” Controls: Watershed Scale Contaminant Flow and Transformation at the Oak Ridge IFRC

Geochemical,spectrographic, microbiological and hydrogeologic studies at the ORIFRC site indicate that groundwater transport in structured media may behave as a system of parallel flow tubes. These tubes are preferred flowpaths that enable contaminant transport parallel to bedding planes (strike) over distances of 1000s of meters. A significant flux of groundwater is focused within an interval defined by the interface between the competent bedrock and overlying highly-weathered saprolite, commonly referred to as the"transition zone." Characteristics of this transition zone are dense fractures and the relative absence of weathering products (e.g. clays)results in a significantly higher permeability compared to both the overlying clay-saprolite and underlying bedrock. Several stratabound low seismic velocity zones located below the transition zone were identified during geophysics studies and were also determined to be fractured high permeability preferred contaminant transport pathways during subsequent drilling activities. XANES analysis of precipitates collected from these deeper flow zones indicate 95% or more of the U deposited is U(VI). Linear combination fitting of the EXAFS data shows that precipitates are ~51±5% U(VI)-carbonate-like phase (e.g., liebigite) and ~49±5% U(VI) associated with an iron oxide phase; inclusion of a third component in the fit suggests that up to 15% of the U(VI) may be associated with a phosphate phase or OH- phase (e.g.,schoepite). Although precipitates with similar U(VI)-carbonate and/or phosphate associations were identified in the transition zone pathways,there were also U(VI) complexes adsorbed to mineral surfaces that would tend to be more readily mobilized. Groundwater in the different flow tubes has been determined to consist of different water quality types that vary with the solid phase encountered (e.g., clays, carbonates, clastics) as contaminants migrate along the flow paths. This lateral and vertical variability in geochemistry, particularly pH, has a significant impact on microbiological community composition and activity. Ribosomal RNA gene analyses coupled with physiological and genomic analyses suggest that bacteria from the genus Rhodanobacter(a diverse population of denitrifiers that are moderately acid tolerant) have a high relative abundance in the acidic source zone at the ORIFRC site.Watershed-scale analysis across different flow paths/tubes revealed strong negative correlation between pH and the absolute and relative abundance of Rhodanobacter. Recent studies also confirmed that the ORIFRC site hosts a diverse fungal community, with significant differences observed between acidic (pH <5) and circumneutral (>5) wells. The lack of nitrous oxide reduction capability in fungi, and the detection of denitrification potential in slurry microcosms suggest that fungi may have aheretofore under appreciated role in biogeochemical transformations, with implications forsite remediation and greenhouse gas emissions. Further research is needed to determine if these organisms can influence U(VI) mobility either directly through immobilization or indirectly through the depletion of nitrate.In conclusion, additional studies are required to quantify the processes (e.g., solid phase reactions, recharge, diffusion, microbial interactions) that are occurring along the groundwater flow tubes identified at the ORIFRC so predictive models can be parameterized and used to assess long-term contaminant fate and transport and remedial options.

An experimental and numerical investigation into the seismic performance of a multi-storey concentrically braced plan irregular structure

In this investigation, the seismic torsional response of a multi-storey concentrically braced frame (CBF) plan irregular structure is evaluated numerically and experimentally through a series of hybrid tests. CBF structures have become popular in seismic design because they are one of the most efficient types of steel structures to resist earthquake loading. However, their response under plan irregular conditions has received little focus mostly in part
due to their complex behaviour under seismic loading conditions. The majority of research on the seismic response of plan irregular structures is based purely on numerical investigations. This paper provides much needed experimental investigation of the seismic response of a CBF plan irregular structure with the aim of characterising the response of this class of structure. The effectiveness of the Eurocode 8 torsional effects provision as a method of designing for
low levels of mass eccentricity is evaluated. Results indicate that some of the observations made by purely numerical models are valid in that; torsionally stiff structures perform well and the stiff side of the structure is subjected to a greater ductility demand compared to the flexible side of the structure. The Eurocode 8 torsional effects provision is shown to be adequate in terms of ductility and interstorey drift however the structure performs poorly
in terms of floor rotation. Importantly, stiffness eccentricity occurs when the provision is applied to the structure when no mass eccentricity exists and results in a significant increase in floor rotations.

Ash from Changbaishan Millennium eruption recorded in Greenland ice: Implications for determining the eruption's timing and impact

Major volcanic eruptions can impact on global climate by injecting large quantities of aerosols and ash into the atmosphere that alter the radiative balance and chemical equilibrium of the stratosphere. The Millennium eruption of Tianchi (Paektu), China/North Korea, was one of the largest Late Holocene eruptions. Uncertainty about the precise timing of the eruption has hindered the recognition of its climate impact in palaeoclimate and historical records. Here we report the compelling identification of the eruption's volcanic signal in Greenland ice cores through the association of geochemically-characterized volcanic glass, represented in by bimodal populations that compare with proximal material from the source eruption. The eruption most probably occurred in the AD 940?s, seven years after the Eldgjá eruption on Iceland. We examine the eruption's potential for climate forcing using the sulfate records from the ice-cores and conclude that it was unlikely to have had a global or extra-regional impact.

Geophysical and hydrogeological characterisation of the impacts of on-site wastewater treatment discharge to groundwater in a poorly productive bedrock acquifer

Contaminants discharging from on-site wastewater treatment systems (OSWTSs) can impact groundwater quality, threatening human health and surface water ecosystems. Risk of negative impacts becomes elevated in areas of extreme vulnerability with high water tables, where thin unsaturated intervals limit vadose zone attenuation. A combined geophysical/hydrogeological investigation into the effects of an OSWTS, located over a poorly productive aquifer (PPA) with thin subsoil cover, aimed to characterise effluent impacts on groundwater. Groundwater, sampled from piezometers down-gradient of the OSWTS percolation area displayed spatially erratic, yet temporally consistent, contaminant distributions. Electrical resistivity tomography identified an area of gross groundwater contamination close to the percolation area and, when combined with seismic refraction and water quality data, indicated that infiltrating effluent reaching the water table discharged to a deeper more permeable zone of weathered shale resting on more competent bedrock. Subsurface structure, defined by geophysics, indicated that elevated chemical and microbiological contaminant levels encountered in groundwater samples collected from piezometers, down-gradient of sampling points with lower contaminant levels, corresponded to those locations where piezometers were screened close to the weathered shale/competent rock interface; those immediately up-gradient were too shallow to intercept this interval, and thus the more impacted zone of the contaminant plume. Intermittent occurrence of faecal indicator bacteria more than 100 m down gradient of the percolation area suggested relatively short travel times. Study findings highlight the utility of geophysics as part of multidisciplinary investigations for OSWTS contaminant plume characterisation, while also demonstrating the capacity of effluent discharging to PPAs to impact groundwater quality at distance. Comparable geophysical responses observed in similar settings across Ireland suggest the phenomena observed in this study are more widespread than previously suspected.

Spatial statistics to estimate peat thickness using airborne radiometric data

Soil carbon stores are a major component of the annual returns required by EU governments to the Intergovernmental Panel on Climate Change. Peat has a high proportion of soil carbon due to the relatively high carbon density of peat and organic-rich soils. For this reason it has become increasingly important to measure and model soil carbon stores and changes in peat stocks to facilitate the management of carbon changes over time. The approach investigated in this research evaluates the use of airborne geophysical (radiometric) data to estimate peat thickness using the attenuation of bedrock geology radioactivity by superficial peat cover. Remotely sensed radiometric data are validated with ground peat depth measurements combined with non-invasive geophysical surveys. Two field-based case studies exemplify and validate the results. Variography and kriging are used to predict peat thickness from point measurements of peat depth and airborne radiometric data and provide an estimate of uncertainty in the predictions. Cokriging, by assessing the degree of spatial correlation between recent remote sensed geophysical monitoring and previous peat depth models, is used to examine changes in peat stocks over time. The significance of the coregionalisation is that the spatial cross correlation between the remote and ground based data can be used to update the model of peat depth. The result is that by integrating remotely sensed data with ground geophysics, the need is reduced for extensive ground-based monitoring and invasive peat depth measurements. The overall goal is to provide robust estimates of peat thickness to improve estimates of carbon stocks. The implications from the research have a broader significance that promotes a reduction in the need for damaging onsite peat thickness measurement and an increase in the use of remote sensed data for carbon stock estimations.

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.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Rapid Petrogenesis of Dacites at Fonualei Volcano

Fonualei is unusual amongst subaerial volcanoes in the Tonga arc because it has erupted dacitic vesicular lavas, tuffs and phreomagmatic deposits for the last 165 years. The total volume of dacite may approach 5 km(3) and overlies basal basaltic andesite and andesite lavas that are constrained to be less than a few millennia in age. All of the products are crystal-poor and formed from relatively low-viscosity magmas inferred to have had temperatures of 1100-1000 degrees C, 2-4 wt % H2O and oxygen fugacities 1-2 log units above the quartz-fayalite-magnetite buffer. Major and trace element data, along with Sr-Nd-Pb and U-Th-Ra isotope data, are used to assess competing models for the origin of the dacites. Positive correlations between Sc and Zr and Sr rule out evolution of the within-dacite compositional array by closed-system crystal fractionation of a single magma batch. An origin by partial melting of lower crustal amphibolites cannot reproduce these data trends or, arguably, any of the dacites either. Instead, we develop a model in which the dacites reflect mixing between two dacitic magmas, each the product of fractional crystallization of basaltic andesite magmas formed by different degrees of partial melting. Mixing was efficient because the two magmas had similar temperatures and viscosities. This is inferred to have occurred at shallow (2-6 km) depths beneath the volcano. U-Th-Ra disequilibria in the basaltic andesite and andesite indicate that the parental magmas had fluids added to their mantle source regions less than 8 kyr ago and that fractionation to the dacitic compositions took less than a few millennia. The 165 year eruption period for the dacites implies that mixing occurred on a similar timescale, possibly during ascent in conduits. The composition of the dacites renders them unsuitable candidates as contributors to average continental crust.

Rapid magmatic processes accompany arc-continent collision: the Western Bismarck arc, Papua New Guinea

New U-Th-Ra, major and trace element, and Sr-Nd-Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc-continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc-continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated Pb-208/Pb-204. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. U-238 and Ra-226 excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.