N-Map: High-resolution groundwater N-retention mapping and modelling by integration of geophysical, geological, geochemical, and hydrological data.

A key aspect of protecting aquatic ecosystems from agricultural nitrogen (N) is to locate (i) farmlands where nitrate leaches from the bottom of the root zone and (ii) denitrifying zones in the aquifers where nitrate is removed before entering the surface water (N-retention). N-retention affects the choice of field mitigation measures to reduce delivered N to surface water. Farmland parcels associated with high N-retention gives the lowest impact of the targeted field measures and vice versa. In Denmark, a targeted N-regulation approach is currently implemented on small catchment scale (approx. 15 km2). Although this regulatory scale is much more detailed than what has been used previously, it is still so large that regulation for most individual fields will be either over- or under-regulated due to large spatial variation in the N-retention. The potential cost reduction for farmers is of up to 20-30% from detailed retention mapping at the field scale compared to the current small catchment scale. In this study, we present a mapping framework (N-Map) for differentiating farmland according to their N-retention, which can be used for improving the effectiveness of targeted N-regulation. The framework currently only includes N-retention in the groundwater. The framework benefits from the incorporation of innovative geophysics in hydrogeological and geochemical mapping and modelling. To capture and describe relevant uncertainties a large number of equally probable realizations are created through Multiple Point Statistical (MPS) methods. This allows relevant descriptions of uncertainties of parts of the model structure and includes other relevant uncertainty measures that affects the obtained N-retention. The output is data-driven high-resolution groundwater N-retention maps, to be used by the individual farmers to manage their cropping systems due to the given regulatory boundary conditions. The detailed mapping allows farmers to use this information in the farm planning in order to optimize the use of field measures to reduce delivered agricultural N to the surface water and thereby lower the costs of the field measures. From farmer interviews, however, it is clear that not all farms will have an economic gain from the detailed mapping as the mapping costs will exceed the potential economic gains for the farmers. The costs of N-Map is here estimated to 5-7 €/ha/year plus implementation costs at the farm. At the society level, the N-retention maps allow authorities to point out opportunities for a more targeted implementation of field measures to efficiently reduce the delivered N-load to surface waters.

Geophysical imaging of the Yellowstone hydrothermal plumbing system.

The nature of Yellowstone National Park's plumbing system linking deep thermal fluids to its legendary thermal features is virtually unknown. The prevailing concepts of Yellowstone hydrology and chemistry are that fluids reside in reservoirs with unknown geometries, flow laterally from distal sources and emerge at the edges of lava flows1-4. Here we present a high-resolution synoptic view of pathways of the Yellowstone hydrothermal system derived from electrical resistivity and magnetic susceptibility models of airborne geophysical data5,6. Groundwater and thermal fluids containing appreciable total dissolved solids significantly reduce resistivities of porous volcanic rocks and are differentiated by their resistivity signatures7. Clay sequences mapped in thermal areas8,9 and boreholes10 typically form at depths of less than 1,000 metres over fault-controlled thermal fluid and/or gas conduits11-14. We show that most thermal features are located above high-flux conduits along buried faults capped with clay that has low resistivity and low susceptibility. Shallow subhorizontal pathways feed groundwater into basins that mixes with thermal fluids from vertical conduits. These mixed fluids emerge at the surface, controlled by surficial permeability, and flow outwards along deeper brecciated layers. These outflows, continuing between the geyser basins, mix with local groundwater and thermal fluids to produce the observed geochemical signatures. Our high-fidelity images inform geochemical and groundwater models for hydrothermal systems worldwide.

Rapid Mapping of Hydrological Systems in Tanzania Using a Towed Transient Electromagnetic System.

Limited knowledge of local groundwater systems often results in the failure of boreholes to yield water of the required quantity and quality. This is particularly problematic in the developing world, where financial resources are often limited, and failed wells represent a significant financial burden. To enhance understanding of local hydrological systems, noninvasive geophysical methods can aid the understanding of hydrogeological structures and identification of groundwater sources needed to optimize siting of wells. Here, we highlight the utility of a relatively new towed-transient electromagnetic system, called tTEM. This system is a rapidly deployable mobile geophysical method well-suited to cost-efficient characterization of local-to-regional groundwater systems. Results from tTEM surveys conducted in two refugee camps and several host communities in western Tanzania demonstrate the capability of the method to characterize shallow aquifer systems with high lateral and vertical resolution, with data collection typically exceeding 15 to 20 line-kilometers (km) per day. This work focuses on tTEM's ability to provide semiquantitative insights into regional hydrogeological settings when supporting data required for more rigorous interpretation/modeling is lacking. The system provided useful data within communities with low density of electrification and near buildings with metal roofs and walls. tTEM-derived resistivity profiles were correlated with limited local borehole lithologic information to develop conceptual models of the local groundwater systems. These models were used to successfully guide the siting of a production well and to identify future drilling targets in the refugee camps and surrounding communities.

Characterizing the diverse hydrogeology underlying rivers and estuaries using new floating transient electromagnetic methodology.

The hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

Cross-borehole tomography with full-decay spectral time-domain induced polarization for mapping of potential contaminant flow-paths.

Soil contamination from industrial activities is a large problem in urban areas worldwide. Understanding the spreading of contamination to underlying aquifers is crucial to make adequate risk assessments and for designing remediation actions. A large part of the northern hemisphere has quaternary deposits consisting of glacial clayey till. The till often has a complex hydrogeological structure consisting of networks of fractures, sand stringers and sand lenses that each contribute to a transport network for water, free phase and dissolved contaminants. Thus, to determine the possible flow-paths of contaminants, the geology must be described in great detail. Normally, multiple boreholes would be drilled in order to describe the geology, but boreholes alone do not provide the needed resolution to map such sand lenses and their connectivity. Cross-borehole full-decay time-domain induced polarization (TDIP) is a new tool that allows for quantitatively mapping not only contrasts in bulk resistivity, but also contrasts in spectral IP parameters. We present a feasibility study with synthetic tests and a field application on a clayey moraine environment with embedded sand lenses, with hitherto unseen ground-truth verification. Indeed, the investigated area was above the water table, which allowed for digging out the entire area after the investigation for an unprecedented description of the lens interconnectivity. The TDIP data were acquired with a full-waveform acquisition at high sampling rate, signal-processed by harmonic denoising, background removal, and de-spiking, and subsequently the full-waveform data were stacked in log-increasing tapered gates (with 7 gates per decade). The resulting TDIP decays, with usable time-gates as early as two milliseconds, were inverted in terms of a re-parameterization of the Cole-Cole model. The inverted models of the field data show a remarkable delineation of the sand lenses/layers at the site, with structure in both the resistivity and the IP parameters matching the results from the ground-truthing. The synthetic examples show that in models both below and above the groundwater table, sand-lenses with thicknesses comparable to the vertical electrode spacing can be well resolved. This suggests that full-decay cross-borehole TDIP is an ideal tool for high-resolution sand-lens imaging.

Large Scale Mapping of Fractures and Groundwater Pathways in Crystalline Hardrock By AEM.

In hardrocks that cover about 20% of the Earth's surface, it is difficult to locate steady sources for groundwater due to inadequate understanding of the fracture networks. A comprehensive knowledge of fracture distribution at the regional scale is necessary to delineate sustainable aquifers and manage them efficiently. The resistivity maps derived from the airborne electromagnetic (AEM) survey over the Ankasandra watershed in Karnataka, India, reveal sharp and deep zones of low formation resistivity, which indicate groundwater-bearing zones. It is found that some of these zones are hydrogeologically connected through fracture networks resulting in augmented yield. AEM results in combination with an in-depth understanding of the geological structures successfully map these groundwater-saturated fracture networks (or hydrogeological lineaments) that we term as 'Hydrolins'. As groundwater occurrence is generally associated with lineaments, we analyzed the drilling and geophysical logs from 21 wells within a 380 sq.km area to study the relationships of various lineaments with 'Hydrolins', particularly in respect of their groundwater potential. AEM results, though calibrated and correlated with a limited number of well data, revealed a threshold groundwater horizon (TGWH), found to be at 80 m depth for Ankasandra watershed, beyond which a strong correlation exists between the depth of a well and its yield. While the TGWH may differ for different watersheds, the approach presented here can be readily adopted to map sustainable groundwater sources in hardrocks worldwide.

Nitrate reduction in geologically heterogeneous catchments--a framework for assessing the scale of predictive capability of hydrological models.

In order to fulfil the requirements of the EU Water Framework Directive nitrate load from agricultural areas to surface water in Denmark needs to be reduced by about 40%. The regulations imposed until now have been uniform, i.e. the same restrictions for all areas independent of the subsurface conditions. Studies have shown that on a national basis about 2/3 of the nitrate leaching from the root zone is reduced naturally, through denitrification, in the subsurface before reaching the streams. Therefore, it is more cost-effective to identify robust areas, where nitrate leaching through the root zone is reduced in the saturated zone before reaching the streams, and vulnerable areas, where no subsurface reduction takes place, and then only impose regulations/restrictions on the vulnerable areas. Distributed hydrological models can make predictions at grid scale, i.e. at much smaller scale than the entire catchment. However, as distributed models often do not include local scale hydrogeological heterogeneities, they are typically not able to make accurate predictions at scales smaller than they are calibrated. We present a framework for assessing nitrate reduction in the subsurface and for assessing at which spatial scales modelling tools have predictive capabilities. A new instrument has been developed for airborne geophysical measurements, Mini-SkyTEM, dedicated to identifying geological structures and heterogeneities with horizontal and lateral resolutions of 30-50 m and 2m, respectively, in the upper 30 m. The geological heterogeneity and uncertainty are further analysed by use of the geostatistical software TProGS by generating stochastic geological realisations that are soft conditioned against the geophysical data. Finally, the flow paths within the catchment are simulated by use of the MIKE SHE hydrological modelling system for each of the geological models generated by TProGS and the prediction uncertainty is characterised by the variance between the predictions of the different models.

Calibrating a salt water intrusion model with time-domain electromagnetic data.

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground-based electromagnetic surveys, electrical resistivity models can be obtained to provide high-resolution three-dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion-State (CHI-S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time-domain electromagnetic (TDEM) dataset was collected. For this location, a simple two-dimensional cross-sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root-mean-square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.