Tracer tests and the structure of permeability in the Corallian limestone aquifer of northern England, UK

The Corallian limestone of northern England (UK) is widely exploited for water supplies and exhibits the karstic phenomena of sinking rivers, conduit development and groundwater velocities of several kilometres per day. To test a number of model-derived source protection zones and elucidate contaminant transport mechanisms in the aquifer, three tracer tests were conducted from a set of swallow-holes draining the River Derwent toward public water supply wells in the eastern part of the aquifer. Tracers used included: Enterobacter cloacae (bacteriophage), Photine C (optical brightener), sodium fluorescein (fluorescent dye) and sulphur hexafluoride (dissolved gas), the varying properties of which make them suitable analogues for different types of potential contaminant. Observed tracer transport times and arrival patterns indicate that tracer transport occurs through karstic channels embedded in a network of primary fissures which exert control over tracer concentrations once initial tracer plumes have passed. A dipole flow system is observed between the swallow-holes and the closest abstraction well, whilst previously modelled source protection zones do not accurately reflect either groundwater velocity or those areas of the aquifer supplying the wells. These findings imply that managing such aquifers for potential contamination should rely upon empirical tracer evidence for source-protection zone modelling.     

Hydrogeophysical Methods for Analyzing Aquifer Storage and Recovery Systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.     

Design and evaluation of a geovisual analytics system for uncovering patterns in spatio-temporal event data

ABSTRACT

It remains difficult to develop a clear understanding of geo-located events and their relationships to one another, particularly when it comes to identifying patterns of events in less-structured textual sources, such as news feeds and social media streams. Here we present a geovisualization tool that can leverage computational methods, such as T-pattern analysis, for extracting patterns of interest from event data streams. Our system, STempo, includes coordinated-view geovisualization components designed to support visual exploration and analysis of event data, and patterns extracted from those data, in terms of time, geography, and content. Through a user evaluation, we explore the usability and utility of STempo for understanding patterns of recent political, social, economic, and military events in Syria.     

Investigating snowpack across scale in the northern Great Lakes–St. Lawrence forest region of Central Ontario,Canada

This study investigates scaling issues by evaluating snow processes and quantifying bias in snowpack properties across scale in a northern Great Lakes–St. Lawrence forest. Snow depth and density were measured along transects stratified by land cover over the 2015/2016 and 2016/2017 winters. Daily snow depth was measured using a time‐lapse (TL) camera at each transect. Semivariogram analysis of the transect data was conducted, and no autocorrelation was found, indicating little spatial structure along the transects. Pairwise differences in snow depth and snow water equivalent (SWE) between land covers were calculated and compared across scales. Differences in snowpack between forested sites at the TL points corresponded to differences in canopy cover, but this relationship was not evident at the transect scale, indicating a difference in observed process across scale. TL and transect estimates had substantial bias, but consistency in error was observed, which indicates that scaling coefficients may be derived to improve point scale estimates. TL and transect measurements were upscaled to estimate grid scale means. Upscaled estimates were compared and found to be consistent, indicating that appropriately stratified point scale measurements can be used to approximate a grid scale mean when transect data are not available. These findings are important in remote regions such as the study area, where frequent transect data may be difficult to obtain. TL, transect, and upscaled means were compared with modelled depth and SWE. Model comparisons with TL and transect data indicated that bias was dependent on land cover, measurement scale, and seasonality. Modelled means compared well with upscaled estimates, but model SWE was underestimated during spring melt. These findings highlight the importance of understanding the spatial representativeness of in situ measurements and the processes those measurements represent when validating gridded snow products or assimilating data into models.     

Rietveld analysis of X‐ray powder diffraction patterns as a potential tool for the identification of impact‐deformed carbonate rocks

Abstract— Previous X‐ray powder diffraction (XRD) studies revealed that shock deformed carbonates and quartz have broader XRD patterns than those of unshocked samples. Entire XRD patterns, single peak profiles and Rietveld refined parameters of carbonate samples from the Sierra Madera impact crater, west Texas, unshocked equivalent samples from 95 miles north of the crater and the Mission Canyon Formation of southwest Montana and western Wyoming were used to evaluate the use of X‐ray powder diffraction as a potential tool for distinguishing impact deformed rocks from unshocked and tectonically deformed rocks. At Sierra Madera dolostone and limestone samples were collected from the crater rim (lower shock intensity) and the central uplift (higher shock intensity). Unshocked equivalent dolostone samples were collected from well cores drilled outside of the impact crater. Carbonate rocks of the Mission Canyon Formation were sampled along a transect across the tectonic front of the Sevier and Laramide orogenic belts. Whereas calcite subjected to significant shock intensities at the Sierra Madera impact crater can be differentiated from tectonically deformed calcite from the Mission Canyon Formation using Rietveld refined peak profiles, weakly shocked calcite from the crater rim appears to be indistinguishable from the tectonically deformed calcite. In contrast, Rietveld analysis readily distinguishes shocked Sierra Madera dolomite from unshocked equivalent dolostone samples from outside the crater and tectonically deformed Mission Canyon Formation dolomite.     

Hydrogen-isotopic variability in lipids from Santa Barbara Basin sediments

We conducted an extensive survey of hydrogen-isotopic compositions (D/H ratios) of diverse sedimentary lipids from the Santa Barbara Basin (SBB), offshore southern California. The main goal of this survey was to assess the diversity of D/H ratios in lipids from marine sediments, in order to provide a more detailed understanding of relevant biological and geochemical factors impacting lipid isotopic variability. A total of 1182 individual δD values are reported from two stations in SBB, one located in the suboxic basin depocenter and the other on the fully oxic flank of the basin. Sediments collected from the basin depocenter span a depth of ∼2.5 m and reach the methanogenic zone. Lipids that were analyzed include n-alkanes, n-alkanols and alkenols, short- and long-chain fatty acids, linear isoprenoids, steroids, and hopanoids, and exhibit several systematic patterns. First, there are no significant differences in δD values between the two sampling locations, nor with increasing depth for most lipids, indicating that degradation does not influence sedimentary lipid δD values. Second, relatively large differences in δD values among differing molecular structures are observed in all samples. n-Alkyl lipids of probable marine origin have typical δD values between −150 and −200‰, those from terrestrial leaf waxes and aquatic plants range from −80 to −170‰, while petroleum n-alkanes are typically −90 to −150‰. Third, lipids inferred to derive from bacteria (branched fatty acids and hopanols) living at the sediment surface or in the water column tend to be D-enriched relative to similar algal products by 30‰ or more. At the same time, several other lipids have δD values that decrease strongly with depth, presumably as a result of in situ production by anaerobic bacteria. This dichotomy in isotopic compositions of bacterial lipids is inconsistent with a nearly constant D/H fractionation during lipid biosynthesis, and likely reflects significant variations associated with metabolism.     

Geometry and architecture of faults in a syn-rift normal fault array: The Nukhul half-graben, Suez rift, Egypt

The geometry and architecture of a well exposed syn-rift normal fault array in the Suez rift is examined. At pre-rift level, the Nukhul fault consists of a single zone of intense deformation up to 10 m wide, with a significant monocline in the hanging wall and much more limited folding in the footwall. At syn-rift level, the fault zone is characterised by a single discrete fault zone less than 2 m wide, with damage zone faults up to approximately 200 m into the hanging wall, and with no significant monocline developed. The evolution of the fault from a buried structure with associated fault-propagation folding, to a surface-breaking structure with associated surface faulting, has led to enhanced bedding-parallel slip at lower levels that is absent at higher levels. Strain is enhanced at breached relay ramps and bends inherited from pre-existing structures that were reactivated during rifting. Damage zone faults observed within the pre-rift show ramp-flat geometries associated with contrast in competency of the layers cut and commonly contain zones of scaly shale or clay smear. Damage zone faults within the syn-rift are commonly very straight, and may be discrete fault planes with no visible fault rock at the scale of observation, or contain relatively thin and simple zones of scaly shale or gouge. The geometric and architectural evolution of the fault array is interpreted to be the result of (i) the evolution from distributed trishear deformation during upward propagation of buried fault tips to surface faulting after faults breach the surface; (ii) differences in deformation response between lithified pre-rift units that display high competence contrasts during deformation, and unlithified syn-rift units that display low competence contrasts during deformation, and; (iii) the history of segmentation, growth and linkage of the faults that make up the fault array. This has important implications for fluid flow in fault zones.