Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers

The validation of variable-density flow models simulating seawater intrusion in coastal aquifers requires information about concentration distribution in groundwater. Electrical resistivity tomography (ERT) provides relevant data for this purpose. However, inverse modeling is not accurate because of the non-uniqueness of solutions. Such difficulties in evaluating seawater intrusion can be overcome by coupling geophysical data and groundwater modeling. First, the resistivity distribution obtained by inverse geo-electrical modeling is established. Second, a 3-D variable-density flow hydrogeological model is developed. Third, using Archie's Law, the electrical resistivity model deduced from salt concentration is compared to the formerly interpreted electrical model. Finally, aside from that usual comparison-validation, the theoretical geophysical response of concentrations simulated with the groundwater model can be compared to field-measured resistivity data. This constitutes a cross-validation of both the inverse geo-electrical model and the groundwater model.
[Comte, J.-C., and O. Banton (2007), Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers, Geophys. Res. Lett., 34, L10402, doi:10.1029/2007GL029981.]

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.

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.

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.

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.