Towards a better understanding of the nature and contribution of the interflow to catchment dynamics.

A detailed understanding of flow and contaminant transfer along each of the key hydrological pathways within a catchment is critical for designing and implementing cost effective Programmes of Measures under the Water
Framework Directive.
The Contaminant Movement along Pathways Project (’The Pathways Project’) is an Irish, EPA STRIVE funded, large multi-disciplinary project which is focussed on understanding and modelling flow and attenuation along each of these pathways for the purposes of developing a catchment management tool. The tool will be used by EPA and RBD catchment managers to assess and manage the impacts of diffuse contamination on stream aquatic ecology. Four main contaminants of interest — nitrogen, phosphorus, sediment and pathogens — are being
investigated in four instrumented test catchments. In addition to the usual hydrological and water chemistry/quality parameters typically captured in catchment studies, field measurements at the test catchments include ecological
sampling, sediment dynamics, soil moisture dynamics, and groundwater levels and chemistry/quality, both during and between significant rainfall events. Spatial and temporal sampling of waters directly from the pathways of
interest is also being carried out.
Sixty-five percent of Ireland is underlain by poorly productive aquifers. In these hydrogeological settings, the main pathways delivering flow to streams are overland flow, interflow and shallow bedrock flow. Little is
known about the interflow pathway and its relative importance in delivery of flow and contaminants to the streams. Interflow can occur in both the topsoil and subsoil, and may include unsaturated matrix flow, bypass or macropore
flow, saturated flow in locally perched water tables and artificial field drainage.
Results to date from the test catchment experiments show that artificial field drains play an important role in the delivery of interflow to these streams, during and between rainfall events when antecedent conditions are
favourable. Hydrochemical mixing models, using silica and SAC254 (the absorbance of UV light at a wavelength of 254 nm which is a proxy for dissolved organic matter) as tracers, show that drain flow is an important end
member contributing to the stream and that proportionally, its contribution is relatively high.
Results from the study also demonstrate that waters originating from one pathway often mix with the waters from another, and are subsequently delivered to the stream at rates, and with chemical/quality characteristics,
that are not typical of either pathway. For example, pre-event shallow groundwater not far from the catchment divide comes up to the surface as rejected recharge during rainfall events and is rapidly delivered to the stream
via overland flow and/or artificial land drainage, bringing with it higher nitrate than would often be expected from a quickflow pathway contribution. This is contrary to the assumption often made in catchment studies that the
deeper hydrological pathways have slower response times in stream hydrographs during a rainfall event, and it emphasizes that it is critical to have a strong three-dimensional conceptual model as the basis for the interpretation
of catchment data.

Pathways Project Final Report Volume 1: Field Investigation and Catchment Conceptual Models

viii
Executive Summary
The Pathways Project field studies were targeted at improving the understanding of contaminant transport along different hydrological pathways in Irish catchments, including their associated impacts on water quality and river ecology. The contaminants of interest were phosphorus, nitrogen and sediment. The working Pathways conceptual model included overland flow, interflow, shallow groundwater flow, and deep groundwater flow. This research informed the development of a set of Catchment Management Support Tools (CMSTs) comprising an Exploratory Tool, Catchment Characterization Tool (CCT) and Catchment Modelling Tool (CMT) as outlined in Pathways Project Final Reports Volumes 3 and 4.
In order to inform the CMST, four suitable study catchments were selected following an extensive selection process, namely the Mattock catchment, Co. Louth/Meath; Gortinlieve catchment, Co. Donegal; Nuenna catchment, Co. Kilkenny and the Glen Burn catchment, Co. Down. The Nuenna catchment is well drained as it is underlain by a regionally important karstified limestone aquifer with permeable limestone tills and gravels, while the other three catchments are underlain by poorly productive aquifers and low permeability clayey tills, and are poorly drained.
All catchments were instrumented, and groundwater, surface and near-surface water and aquatic ecology were monitored for a period of two years. Intensive water quality sampling during rainfall events was used to investigate the pathways delivering nutrients. The proportion of flow along each pathway was determined using chemical and physical hydrograph separation techniques, supported by numerical modelling.
The outcome of the field studies broadly supported the use of the initial four-pathway conceptual model used in the Pathways CMT (time-variant model). The artificial drainage network was found to be a significant contributing pathway in the poorly drained catchments, at low flows and during peak flows in wet antecedent conditions. The transition zone (TZ), i.e. the broken up weathered zone at the top of the bedrock, was also found to be an important pathway. It was observed to operate in two contrasting hydrogeological scenarios: in groundwater discharge zones the TZ can be regarded as being part of the shallow groundwater pathway, whereas in groundwater recharge zones it behaves more like interflow.
In the catchments overlying poorly productive aquifers, only a few fractures or fracture zones were found to be hydraulically active and the TZ, where present, was the main groundwater pathway. In the karstified Nuenna catchment, the springs, which are linked to conduits as well as to a diffuse fracture network, delivered the majority of the flow. These findings confirm the two-component groundwater contribution from bedrock but suggest that the size and nature of the hydraulically active fractures and the nature of the TZ are the dominant factors at the scale of a stream flow event.
Diffuse sources of nitrate were found to be typically delivered via the subsurface pathways, especially in the TZ and land drains in the poorly productive aquifer catchments, and via the bedrock groundwater in the Nuenna. Phosphorus was primarily transported via overland flow in both particulate and soluble forms. Where preferential flow paths existed in the soil and subsoil, soluble P, and to a lesser extent particulate P, were also transported via the TZ and in drains and ditches. Arable land was found to be the most important land use for
ix
the delivery of sediment, although channel bank and in-stream sources were the most significant in the Glen Burn catchment. Overland flow was found to be the predominant transport sediment pathway in the poorly productive catchments. These findings informed the development of the transport and attenuation equations used in the CCT and CMT. From an assessment of the relationship between physico-chemical and biological conditions, it is suggested that in the Nuenna, Glen Burn and Gortinlieve catchments, a relationship may exist between biological water quality and nitrogen concentrations, as well as with P. In the Nuenna, there was also a relationship between macroinvertebrate status and alkalinity.
Further research is recommended on the transport and delivery of phosphorus in groundwater, the transport and attenuation dynamics in the TZ in different hydrogeological settings and the relationship between macroinvertebrates and co-limiting factors. High resolution temporal and spatial sampling was found to be important for constraining the conceptual understanding of nutrient and sediment dynamics which should also be considered in future studies.

The influence of bedrock hydrogeology on catchment-scale nitrate fate and transport in fractured aquifers

Characterising catchment scale biogeochemical processes controlling nitrate fate in groundwater constitutes a fundamental consideration when applying programmes of measures to reduce risks posed by diffuse agricultural pollutants to water quality. Combining hydrochemical analyses with nitrate isotopic data and physical hydrogeological measurements permitted characterisation of biogeochemical processes influencing nitrogen fate and transport in the groundwater in two fractured bedrock aquifers with contrasting hydrogeology but comparable nutrient loads. Hydrochemical and isotopic analyses of groundwater samples collected from moderately fractured, diffusely karstified limestone indicated nitrification controlled dissolved nitrogen fate and delivery to aquatic receptors. By contrast nitrate concentrations in groundwater were considerably lower in a low transmissivity highly lithified sandstone and pyrite-bearing shale unit with patchy subsoil cover. Geophysical and hydrochemical investigations showed shallower intervals contained hydraulically active fractures where denitrification was reflected through lower nitrogen levels and an isotopic enrichment ratio of 1.7 between δ15N and δ18O. Study findings highlight the influence of bedrock hydrogeological conditions on aqueous nitrogen mobility. Investigation results demonstrate that bedrock conditions need to be considered when implementing catchment management plans to reduce the impact of agricultural practices on the quality of groundwater and baseflow in receiving rivers.

Dynamics of Finnish starlings in 1951-2005 : from monitoring to population modelling

Farmland bird species have been declining in Europe. Many declines have coincided with general intensification of farming practices. In Finland, replacement of mixed farming, including rotational pastures, with specialized cultivation has been one of the most drastic changes from the 1960s to the 1990s. This kind of habitat deterioration limits the persistence of populations, as has been previously indicated from local populations. Integrated population monitoring, which gathers species-specific information of population size and demography, can be used to assess the response of a population to environment changes also at a large spatial scale. I targeted my analysis at the Finnish starling (Sturnus vulgaris). Starlings are common breeders in farmland habitats, but severe declines of local populations have been reported from Finland in the 1970s and 1980s and later from other parts of Europe. Habitat deterioration (replacement of pasture and grassland habitats with specialized cultivation areas) limits reproductive success of the species. I analysed regional population data in order to exemplify the importance of agricultural change to bird population dynamics. I used nestling ringing and nest-card data from 1951 to 2005 in order to quantify population trends and per capita reproductive success within several geographical regions (south/north and west/east aspects). I used matrix modelling, acknowledging age-specific survival and fecundity parameters and density-dependence, to model population dynamics. Finnish starlings declined by 80% from the end of the 1960s up to the end of the 1980s. The observed patterns and the model indicated that the population decline was due to the decline of the carrying capacity of farmland habitats. The decline was most severe in north Finland where populations largely become extinct. However, habitat deterioration was most severe in the southern breeding areas. The deteriorations in habitat quality decreased reproduction, which finally caused the decline. I suggest that poorly-productive northern populations have been partly maintained by immigration from the highly-productive southern populations. As the southern populations declined, ceasing emigration caused the population extinction in north. This phenomenon was explained with source sink population dynamics, which I structured and verified on the basis of a spatially explicit simulation model. I found that southern Finnish starling population exhibits ten-year cyclic regularity, a phenomenon that can be explained with delayed density-dependence in reproduction.

Brackish-water fishery resources

Ceylon has about 300,000 acres of coastal brackish-water areas of which about 100,000 acres constitute shallow lagoons, tidal flats, mangrove swamps and saline marshes, and the rest deeper lagoons and estuaries. While the former represent a vast potential resource with regard to fish farming, the latter are the sites of important fisheries. W. H. Schuster (1951) estimated the average natural production of Ceylon’s brackish waters to be less than 20 lbs. per acre per annum. Since then estimates have been made by the author for a rich lagoon, the Negombo lagoon, a poorly productive lagoon, the Ratgama lake (Dodanduwa) and studies are in progress of some of the other lagoons. The natural production of the Ratgama lake was estimated in 1959 to be 18.5 lbs. per acre per annum while that of Negombo lagoon was estimated in 1960 to be 65 lbs. per acre per annum. It is reasonable to estimate the average production of Ceylon's brackish-waters to be 25 lbs. per acre per annum. Thus the total production is about 3,350 tons per annum. Considering the fact that the island's present total production is 90,000 tons per annum, the brackish-waters contribute 3.7% of it. Schuster (1951) further states that the natural production in the brackish-waters of other countries is around 80 lbs. per acre per annum. In order to increase our average natural production to this value it would seem necessary to consider the nature, biology and fish resources of the brackish-waters and draw some conclusions with regard to their proper exploitation.

Hydraulic Subsurface Measurements and Hydrodynamic Modeling as Indicators for Groundwater Flow Systems in the Rotondo Granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.

Detection of Groundwater flow zones in low productivity boreholes (Submitted)

Conventional methods of detecting groundwater flow zones in open boreholes installed in fractured bedrock aquifers rely on either contrasts in water quality or on significant rates of vertical flow. In low productivity boreholes these methods have proced of limited value. Tracer tests completed in six low productivity bedrock boreholes installed into fresh Precambrian crystalline bedrock revealed measurable dilution, despite very low tranmissivities. Fluorescent tracer profiles generated during discharge pumping permitted identification of the principle zones contributing to flow. High resolution acoustic televiewer strike and dip measurements of fractures corresponding to these zones revealed a consistent pattern of regional lineament trends and suggested a strongly anisotropic flow pattern in bedrock.

Simulation of Ground-Water Flow in the Cedar River Alluvial Aquifer Flow System, Cedar Rapids, Iowa Scientific

The Cedar River alluvial aquifer is the primary source of municipal water in the Cedar Rapids, Iowa, area. Since 1992, the U.S. Geological Survey, in cooperation with the City of Cedar Rapids, has investigated the hydrogeology and water quality of the Cedar River alluvial aquifer. This report describes a detailed analysis of the ground-water flow system in the alluvial aquifer, particularly near well field areas. The ground-water flow system in the Cedar Rapids area consists of two main components, the unconsolidated Quaternary deposits and the underlying carbonate bedrock that has a variable fracture density. Quaternary deposits consist of eolian sand, loess, alluvium, and glacial till. Devonian and Silurian bedrock aquifers overlie the Maquoketa Shale (Formation) of Ordovician age, a regional confining unit. Ground-water and surface-water data were collected during the study to better define the hydrogeology of the Cedar River alluvial aquifer and Devonian and Silurian aquifers. Stream stage and discharge, ground-water levels, and estimates of aquifer hydraulic properties were used to develop a conceptual ground-water flow model and to construct and calibrate a model of the flow system. This model was used to quantify the movement of water between the various components of the alluvial aquifer flow system and provide an improved understanding of the hydrology of the alluvial aquifer.

Constraining age distributions of groundwater from public supply wells in diverse hydrogeological settings in Scania, Sweden

Twenty-five public supply wells throughout the hydrogeologically diverse region of Scania, southern Sweden are subjected to environmental tracer analysis (³H–³He,⁴He, CFCs, SF₆ and for one well only also ⁸⁵Kr and ³⁹Ar) to study well and aquifer vulnerability and evaluate possibilities of groundwater age distribution assessment. We find CFC and SF₆ concentrations well above solubility equilibrium with modern atmosphere, indicating local contamination, as well as indications of CFC degradation. The tracer-specific complications considerably constrain possibilities for sound quantitative regional ground- water age distribution assessment and demonstrate the importance of initial qualitative assessment of tracer-specific reliability, as well a need for additional, complementary tracers (e.g. ⁸⁵Kr,³⁹Ar and potentially also ¹⁴C). Lumped parameter modelling yields credible age distribution assessments for representative wells in four type aquifers. Pollution vulnerability of the aquifer types was based on the selected LPM models and qualitative age characterisation. Most vulnerable are unconfined dual porosity and fractured bedrock aquifers, due to a large component of very young groundwater. Unconfined sedimentary aquifers are vulnerable due to young groundwater and a small pre-modern component. Less vulnerable are semi-confined sedimentary or dual-porosity aquifers, due to older age of the modern component and a larger pre-modern component. Confined aquifers appear least vulnerable, due an entirely pre-modern groundwater age distribution (recharged before 1963). Tracer complications aside, environmental tracer analyses and lumped parameter modelling aid in vulnerability assessment and protection of regional groundwater resources.

A comparison of nitogen fate and transport in catchments underlain by high and low permeability aquifers using chemical and isotopic tools.

Variability in nitrogen fate and transport in different catchments types is often not considered. This research considers the importance of such nitrogen processes within groundwater pathways in two agricultural catchments in Ireland; a well drained catchment, underlain by karstified Carboniferous limestone, and a poorly drained catchment, underlain by Silurian greywacke.
Depth specific low-flow groundwater sampling was used to evaluate the hydrochemical stratification in groundwater. Groundwater samples, as well as surface water samples, along river courses were analysed for nitrogen species (NO3, NH4 and NO2) and nitrate isotopes (d15N and d18O) as well as field parameters and major ions
.
The dominant nitrate (NO3) groundwater pathway in the poorly drained greywacke catchment is through the shallow weathered bedrock, as indicated by transmissivity values and the ionic and isotopic signatures, and a clear reduction in NO3 concentration is observed with depth. A similar chloride trend would suggest dilution is a major factor, however d15N and d18O isotopic values producing an enrichment ratio of 1.8 indicate that denitrification is also an important process involved in the fate of the NO3 within the groundwater flow system. This consistent trend with depth is in contrast to the stratification pattern observed in the karstified catchment. NO3 was not detected in the shallow groundwater pathway; the dominant groundwater pathway is in the deeper groundwater where there is little change in the nitrate isotope values with depth (d15N values range between 4.1 and 4.6 ‰). This deeper groundwater contributes the dominant proportion of the river flow through a number of springs. As a result, the deeper groundwater, springs and river have a similar ionic signature and NO3 concentration range (23 ± 3 mg/l). Despite this pattern, the NO3 isotopes show a distinct difference in isotopic values between the deeper groundwater in the diffuse karst and the springs indicating some denitrification is occurring during groundwater discharge into the river. Furthermore the isotopes give an indication of the variability of the spatial extent of the springs and the complexities of the fissures through which they are fed. The results of this study clearly show the importance of the geology in the fate and transport of NO3 in agricultural catchments.

Assessment of interactions between alluvial, Great Artesian Basin (GAB) and volcanic aquifers using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia

A detailed 3D lithological model framework was developed using GOCAD software to understand interactions between alluvial, volcanic and GAB aquifers and the spatial and temporal distribution of groundwater recharge to the alluvium of the Lockyer Valley. Groundwater chemistry, isotope data (H20-δ2H and δ18O , 87Sr/86Sr, 3H and 14C) and groundwater level time-series data from approximately 550 observation wells were integrated into the catchment-wide 3D model to assess the recharge processes involved. This approach enabled the identification of zones where recharge to the alluvium primarily occurs from stream water during episodic flood events. Importantly, the study also demonstrates that in some sections of the alluvium recharge is also from storm rainfall and seepage discharge from the underlying GAB aquifers. These other sources of recharge are indicated by (a) the absence of a response of groundwater levels to flooding in some areas, (b) old radiocarbon ages, and (c) distinct bedrock water chemistry and δ2H and δ18O signatures in alluvial groundwater at these locations. Integration of isotopes, water chemistry and time-series displays of groundwater levels before and after the 2010/2011 flood into the 3D model suggest that the spatial variations in the alluvial groundwater response are mostly controlled by valley morphology and lithological (i.e. permeability) variations within the alluvium. Examination of the groundwater level variations in the 3D model also enabled quantification of the volumetric change of groundwater stored in the unconfined alluvial aquifer prior to and post-flood events.

Mechanics of sediment transport and bedrock erosion in steep landscapes

Erosion is concentrated in steep landscapes such that, despite accounting for only a small fraction of Earth’s total surface area, these areas regulate the flux of sediment to downstream basins, and their rugged morphology records transient changes (or lack thereof) in geologic and climatic forcing. Steep landscapes are geomorphically active; large sediment fluxes and rapid landscape evolution rates can create or destroy habitat for humans and wildlife alike, and landslides, debris flows, and floods common in mountainous areas represent a persistent natural and structural hazard. Despite the central role that steep landscapes play in the geosciences and in landscape management, the processes controlling their evolution have been poorly studied compared to lower-gradient areas. This thesis focuses on the basic mechanics of sediment transport and bedrock incision in steep landscapes, as these are the fundamental processes which set the pace and style of landscape evolution. Chapter 1 examines the spatial distribution of slow-moving landslides; these landslides can dominate sediment fluxes to river networks, but the controls on their occurrence are poorly understood. Using a case-study along the San Andreas Fault, California, I show that slow-moving landslides preferentially occur near the fault, suggesting a rock-strength control on landslide distribution. Chapter 2 provides the first field-measurements of incipient sediment motion in streams steeper than 14% and shows a large influence of slope-dependent flow hydraulics and grain-scale roughness on particle motion. Chapter 3 presents experimental evidence for bedrock erosion by suspended sediment, suggesting that, in contrast to prevailing theoretical predictions, suspension-regime transport in steep streams can be the dominant erosion agent. Steep streams are often characterized by the presence of waterfalls and bedrock steps which can have locally high rates of erosion; Chapters 4 and 5 present newly developed, experimentally validated theory on sediment transport through and bedrock erosion in waterfall plunge pools. Finally, Chapter 6 explores the formation of a bedrock slot canyon where interactions between sediment transport and bedrock incision lead to the formation of upstream-propagating bedrock step-pools and waterfalls.

The typology of Irish hard rock aquifers based on an integrated hydrogeological and geophysical approach

Groundwater flow in hard-rock aquifers is strongly controlled by the characteristics and distribution of structural heterogeneity. A methodology for catchment-scale characterisation is presented, based on the integration of complementary, multi-scale hydrogeological, geophysical and geological approaches. This was applied to three contrasting catchments underlain by metamorphic rocks in the northern parts of Ireland (Republic of Ireland and Northern Ireland, UK). Cross-validated surface and borehole geophysical investigations confirm the discontinuous overburden, lithological compartmentalisation of the bedrock and important spatial variations of the weathered bedrock profiles at macro-scale. Fracture analysis suggests that the recent (Alpine) tectonic fabric exerts strong control on the internal aquifer structure at meso-scale, which is likely to impact on the anisotropy of aquifer properties. The combination of the interpretation of depth-specific hydraulic-test data with the structural information provided by geophysical tests allows characterisation of the hydrodynamic properties of the identified aquifer units. Regionally, the distribution of hydraulic conductivities can be described by inverse power laws specific to the aquifer litho-type. Observed groundwater flow directions reflect this multi-scale structure. The proposed integrated approach applies widely available investigative tools to identify key dominant structures controlling groundwater flow, characterising the aquifer type for each catchment and resolving the spatial distribution of relevant aquifer units and associated hydrodynamic parameters.

Codeposition of organic carbon and arsenic in Bengal Delta aquifers

We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p <0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As-OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments.

Co-deposition of organic carbon and arsenic in Bengal Delta aquifers.

We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p < 0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As/OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments. We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p < 0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As?OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments.