Cave and karst development in the Craven Basin, UK

The Craven Basin, which was an area of crustal subsidence during much of the Carboniferous, contains carbonate strata that pre-date the deposition of the better known carbonates of the Askrigg Block to the north. Later basin-fill deposits consist mainly of clastic sedimentary rocks. However, interbedded carbonate horizons and Waulsortian mud-mounds are also present. Karst development has taken place locally on the exposed carbonate strata. including establishment of the only example of a sulphurous cave stream currently known in the British Isles. The caves so far examined in the area are described, alongside an account of their geological setting

Numerical simulation of spring hydrograph recession curves for evaluating behavior of the East Yorkshire chalk aquifer

The Cretaceous Chalk aquifer is the most important in the UK for the provision of water to public supply and agriculture. The Chalk has both matrix and fracture porosity and is thus best considered as a dual porosity aquifer system. Although the matrix porosity is large, typically around 0.35 in the study area of East Yorkshire, UK (ESI, 2010), pore diameters are typically very small, and the water contained in them is virtually immobile. The high permeability fracture network is responsible for the ability of water to drain; spatial variations in fracture network properties mean conventional approaches to aquifer characterization such as borehole pumping tests are of limited utility. Hence this study attempts to better understand the flow system and characterise aquifer properties from the recession response seen at springs during the spring/summer period when recharge is minimal. This approach has the advantage that spring hydrographs represent the sum of the response from entire catchments. This paper reports numerical modeling for simulating aquifer and spring responses during hydrological recession. Firstly, available geological and hydrogeological information for the study area was used to develop hydrogeological conceptual models. Four different numerical models have been constructed representing four possible scenarios that could represent the aquifer in the selected area. These are: single reservoir aquifer, double reservoir aquifer, single reservoir aquifer with highly permeable vertical zone intersecting the spring location and single reservoir aquifer containing tunnel shaped highly permeable zone at the spring elevation respectively. The sensitivity of spring recession response to various external and internal parameter values was investigated, to understand relations between spring recession, hydrological inputs (recharge) and aquifer structure. Spring hydrographs from the real aquifer were compared with the hydrographs generated from models, in order to estimate aquifer properties. The work aims to identify the utility of spring hydrographs in eliciting aquifer permeability structure, as well as identifying the conceptual scenario which best represents the Chalk Aquifer in East Yorkshire, UK

Oxygen isotopes in tree rings show good coherence between species and sites in Bolivia

A tree ring oxygen isotope (δ18OTR) chronology developed from one species (Cedrela odorata) growing in a single site has been shown to be a sensitive proxy for rainfall over the Amazon Basin, thus allowing reconstructions of precipitation in a region where meteorological records are short and scarce. Although these results suggest there should be large-scale (> 100 km) spatial coherence of δ18OTR records in the Amazon, this has not been tested. Furthermore, it is of interest to investigate whether other, possibly longer -lived, species similarly record interannual variation of Amazon precipitation, and can be used to develop climate sensitive isotope chronologies. In this study, we measured δ18O in tree rings from seven lowland and one highland tree species from Bolivia. We found that cross-dating with δ18OTR gave more accurate tree ring dates than using ring width. Our “isotope cross-dating approach” is confirmed with radiocarbon “bomb-peak” dates, and has the potential to greatly facilitate development of δ18OTR records in the tropics, identify dating errors, and check annual ring formation in tropical trees. Six of the seven lowland species correlated significantly with C. odorata, showing that variation in δ18OTR has a coherent imprint across very different species, most likely arising from a dominant influence of source water δ18O on δ18OTR. In addition we show that δ18OTR series cohere over large distances, within and between species. Comparison of two C. odorata δ18OTR chronologies from sites several hundreds of kilometres apart showed a very strong correlation (r = 0.80, p < 0.001, 1901–2001), and a significant (but weaker) relationship was found between lowland C. odorata trees and a Polylepis tarapacana tree growing in the distant Altiplano (r = 0.39, p < 0.01, 1931–2001). This large-scale coherence of δ18OTR records is probably triggered by a strong spatial coherence in precipitation δ18O due to large-scale controls. These results highlight the strength of δ18OTR as a precipitation proxy, and open the way for temporal and spatial expansion of precipitation reconstructions in South America

Earliest Triassic microbialites in the South China Block and other areas; controls on their growth and distribution

Earliest Triassic microbialites (ETMs) and inorganic carbonate crystal fans formed after the end-Permian mass extinction (ca. 251.4 Ma) within the basal Triassic Hindeodus parvus conodont zone. ETMs are distinguished from rarer, and more regional, subsequent Triassic microbialites. Large differences in ETMs between northern and southern areas of the South China block suggest geographic provinces, and ETMs are most abundant throughout the equatorial Tethys Ocean with further geographic variation. ETMs occur in shallow-marine shelves in a superanoxic stratified ocean and form the only widespread Phanerozoic microbialites with structures similar to those of the Cambro-Ordovician, and briefly after the latest Ordovician, Late Silurian and Late Devonian extinctions. ETMs disappeared long before the mid-Triassic biotic recovery, but it is not clear why, if they are interpreted as disaster taxa. In general, ETM occurrence suggests that microbially mediated calcification occurred where upwelled carbonate-rich anoxic waters mixed with warm aerated surface waters, forming regional dysoxia, so that extreme carbonate supersaturation and dysoxic conditions were both required for their growth. Long-term oceanic and atmospheric changes may have contributed to a trigger for ETM formation. In equatorial western Pangea, the earliest microbialites are late Early Triassic, but it is possible that ETMs could exist in western Pangea, if well-preserved earliest Triassic facies are discovered in future work

A qualitative and quantitative model for climate-driven lake formation on carbonate platforms based on examples from the Bahamian archipelago

Lakes on carbonate platform islands such as the Bahamas display wide variability in morphometry, chemistry, and fauna. These parameters are ultimately driven by climate, sea level, and carbonate accumulation and dissolution. The authors propose a model that integrates climatological, geomorphological, and stratigraphic frameworks to understand processes of carbonate-hosted lake formation and limnological characteristics in modern day environments, with applications to carbonate lake sedimentary records. Fifty-two lakes from San Salvador Island and Eleuthera, Bahamas, were examined for water chemistry, basin morphology, conduit development, conductivity, and major ions. Using non-metric, multi-dimensional scaling ordination methods, the authors derived a model dividing lakes into either constructional or destructional formational modes. Constructional lakes were further divided into pre-highstand and highstand types based on whether their formation occurred during a marine regressive or transgressive phase. Destructional lakes are created continually by dissolution of bedrock at fresh/saline water interfaces and their formation is therefore related to changing climate and sea level. This model shows that lake formation is influenced by the hydrologic balance associated with climatic conditions that drives karst dissolution as well as the deposition of aeolian dune ridges that isolate basins due to sea-level fluctuations. It allows for testing and examining the climatic and hydrologic regime as related to carbonate accumulation and dissolution through time, and for an improved understanding of lake sensitivity and response to climate as preserved in the lacustrine sedimentary record

Function of a deltaic silt deposit as a repository and long-term source of sulfate and related weathering products in a glaciofluvial aquifer derived from organic-rich shale (North Dakota, USA)

A shallow unconfined glaciofluvial aquifer in North Dakota (USA) has largest groundwater sulfate concentrations near the bottom boundary. A deltaic silt layer underlying the aquifer, at >16 m, is the modern proximate sulfate source for the aquifer. The original sulfate source was pyrite in the organic-rich shale component of the aquifer and silt grain matrix. An oxidizing event occurred during which grain-matrix pyrite sulfur was oxidized to sulfate. Thereafter the silt served as a "conserving" layer, slowly feeding sulfate into the lower part of the aquifer and the underlying till. A method was developed for estimating the approximate initial sulfate concentration in the source layer and the redistribution time since the oxidizing event, using a semi-generic convection-dispersion model. The convection-dispersion model and a model for the evolution of modern sulfate δ 34S in silt-layer pore water from the initial grain-matrix pyrite δ 34S, both estimated that the oxidizing event occurred several thousand years ago, and was likely related to the dry conditions of the Hypsithermal Interval. The silt layer also serves as an arsenic source. Results indicate that deltaic silts derived from organic-rich shale parent materials in a glacial environment can provide long-term sources for sulfate and arsenic and possibly other related oxidative weathering products

Carbon and sulphur geochemistry and clay mineralogy of the West Runton Freshwater Bed

Weight % organic carbon and weight % total sulphur (C/S) in the Cromerian West Runton Freshwater Bed (WRFWB) of north Norfolk are not simply related to depositional environment as in the classic Berner and Raiswell, (1984) model. This is because the sediment is organic matter-rich, with total sulphur representing a variable mixture of pyrite sulphurs and organic sulphur. In addition, remobilisation of sulphur associated with post-depositional groundwater flow has modified depositional values. In the upper 20 cm of the bed, pyrite sulphur has been oxidised and largely removed. This sulphur was subsequently reprecipitated as later diagenic pyrite at about 40-60 cm in the bed under reducing conditions. Below 60 cm C/S ratios are probably close to depositional values and can be used as palaeosalinity indicators. As expected, C/pyrite-S ratios in this lower part of the bed indicate a freshwater depositional environment. The groundwater which oxidised pyrite in the top of the bed also destroyed organic carbon and probably dissolved aragonitic shell material. However, preservation of aragonitic shells below 60 cm in the WRFWB defines the extent to which groundwater penetrated the bed. Smectite is one of the dominant clay minerals in the WRFWB and potassium saturation suggests that the smectite may have had a volcanic orgin. The Eifel region of Germany is a possible source of fine-grained ash or dust at this tim

Characterisation of fractured carbonate aquifers using ambient borehole dilution tests

Fractured carbonate aquifers derive their transmissivity essentially from a well-developed network of solutionally-enhanced fractures and conduits that can lead to high groundwater velocities and high vulnerability to contamination of water quality. Characterisation of the variation of hydraulic properties with depth is important for delineating source protection areas, characterising contaminant fate and transport, determination of the effectiveness of aquifer remediation, and parameter estimation for models. In this work, ambient open borehole uniform and point injection dilution tests were conducted on observation boreholes in the unconfined Cretaceous Chalk aquifer of East Yorkshire, UK, and interpreted in conjunction with other data via the implementation of a new work flow. This resulted in the characterisation of flow in these boreholes and the inference of properties such as groundwater flow patterns and velocities in the surrounding aquifer formation. Our workflow allowed sections of open boreholes showing horizontal versus vertical flow to be distinguished, and the magnitude of such flows and exchanges with the aquifer to be determined. Flow within boreholes were then used to characterise: i) presence and direction of vertical hydraulic gradients; ii) nature and depth distribution of flowing features; iii) depth interval porosity and permeability estimation of the flowing features from overall borehole transmissivity and geophysical image or caliper logs; iv) groundwater velocity estimation in the surrounding aquifer. Discrete flowing features were distributed across the range of depths sampled by the observation boreholes (typically up to 45–60 mbgl), but the majority were located in the zone of water table fluctuation marked by solutionally enlarged flow features. Quantitative interpretation of both uniform injection (tracer distributed throughout the open borehole section) and point injection (slug of tracer introduced at targeted depth) yielded vertical velocities within the borehole water column in broad agreement with those measured by flow logging. Depth specific fracture kinematic porosities inferred from the ambient dilution data combined with long-interval pump test and geophysical log data ranged between 3.7 × 10−4–4.1 × 10−3 with an average of 2.1 × 10−3; these values were in excellent agreement with those from other methods applied to the same aquifer such as larger scale pumping tests. A new approach to estimation of groundwater velocities from the dilution test data using externally measured hydraulic gradients gave inferred horizontal groundwater velocities ranging between 60 and 850 m/day, in full agreement with those from previously conducted borehole-to-borehole tracer tests. These results confirm that the studied aquifer is karstic, with rapid preferential pathways which have implication for flow and transport modelling, and pollution vulnerability. Our study results indicate that ambient single-borehole dilution approaches can provide an inexpensive and reliable approach for the characterisation of fractured and karstic aquifers

Sulphur springs of the Craven Basin, NW England: indicators of natural methane leakage?

Sulphidic springs are a common feature of the Craven Basin and have long histories, indicating that they are natural phenomena, and not a result of recent anthropogenic contamination of groundwater. This study identified several extant sulphur springs and these were sampled and analysed for hydrochemistry and isotopic composition (δS of sulphate and sulphide, δO of sulphate and δC of total dissolved inorganic carbon (TDIC)). The springs are associated with limestone units and are located in anticlinal crests that are sometimes faulted. Sulphide concentrations at these 'sulphur springs' range from <0.5 to 96 mg l and light sulphide isotopic compositions indicate that sulphide originates by microbial sulphate reduction. One site, 'Stinky Bottoms', with high sulphide concentration (46 mg l) also has very low δC-TDIC (-25.3‰ V-PDB) characteristic of a significant component of TDIC generation via methane oxidation. Other sites have δC-TDIC more typical of shallow groundwaters (-14.2 to -16.3‰ V-PDB) and two sites with highest Cl concentration and elevated Sr:Ca have heavier δCTDIC (-11.3 and -6.8‰ V-PDB), indicative of a more evolved, long residence-time brine component. At Stinky Bottoms a strong case can be made that generation of sulphidic groundwater is related to subsurface methane. At other sites, sulphide generation may also be related to hydrocarbon or methane seeps (in some cases associated with a brine component) but the δC-TDIC values are more equivocal