Automated image analysis for experimental investigations of salt water intrusion in coastal aquifers

A novel methodology has been developed to quantify important saltwater intrusion parameters in a sandbox style experiment using image analysis. Existing methods found in the literature are based mainly on visual observations, which are subjective, labour intensive and limits the temporal and spatial resolutions that can be analysed. A robust error analysis was undertaken to determine the optimum methodology to convert image light intensity to concentration. Results showed that defining a relationship on a pixel-wise basis provided the most accurate image to concentration conversion and allowed quantification of the width of mixing zone between the saltwater and freshwater. A large image sample rate was used to investigate the transient dynamics of saltwater intrusion, which rendered analysis by visual observation unsuitable. This paper presents the methodologies developed to minimise human input and promote autonomy, provide high resolution image to concentration conversion and allow the quantification of intrusion parameters under transient conditions.

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.

Experimental Saltwater Intrusion in Coastal aquifers Using Automated Image Analysis: Applications to Homogeneous Aquifers.

This paper presents the applications of a novel methodology to quantify saltwater intrusion parameters in laboratory-scale experiments. The methodology uses an automated image analysis procedure, minimizing manual inputs and the subsequent systematic errors that can be introduced. This allowed the quantification of the width of the mixing zone which is difficult to measure in experimental methods that are based on visual observations. Glass beads of different grain sizes were tested for both steady-state and transient conditions. The transient results showed good correlation between experimental and numerical intrusion rates. The experimental intrusion rates revealed that the saltwater wedge reached a steady state condition sooner while receding than advancing. The hydrodynamics of the experimental mixing zone exhibited similar
traits; a greater increase in the width of the mixing zone was observed in the receding saltwater wedge, which indicates faster fluid velocities and higher dispersion. The angle of intrusion analysis revealed the formation of a volume of diluted saltwater at the toe position when the saltwater wedge is prompted to recede. In addition, results of different physical repeats of the experiment produced an average coefficient of variation less than 0.18 of the measured toe length and width of the mixing zone.

Investigating Fractured Flow from Orthopaedics to Rehabilitation

Patients often spend time as inpatients in hospitals outside their home area because of the need to access specialist services. If there is a need for ongoing care after the need for specialist care has passed, patients are transferred (or ‘repatriated’) back to the inpatient care of a hospital in their local Health and Social Care Trust. The need for bed space in specialist units means that there is pressure for this transition to occur in a timely way. We investigated the flow of patients through a trauma and orthopaedics unit using the 6M Design® framework and Vitals Charts® in order to investigate concerns about delayed repatriation. We found that repatriation was part of a complex system that had interdependent components. There was considerable variation in the number of discharges (to any destination) by day of week, with a reduction on Saturdays and Sundays. Understanding that the pressure for quicker repatriation was really due to high work-in-progress led us to model the effects of strategies to address the high work-in-progress. We found that, because only a small proportion of patients require repatriation, expediting the repatriation process by one day for each patient would only reduce WIP by an average of 1.6 patients. Reducing the average length of stay for all trauma and orthopaedics inpatients by one day would reduce the WIP by 10 patients, which would make a much greater impact on the problem of high bed occupancy. Though the smooth and timely repatriation of patients to rehabilitation units is desirable, it is unlikely that efforts to achieve this will have a substantial impact on the problem of high WIP, so other strategies will be required. We will model the effects of strategies to reduce variation in daily discharges by the day of week in a future essay.

Analysis of diffusion process in fractured reservoirs using fractional derivative approach

The fractal geometry is used to model of a naturally fractured reservoir and the concept of fractional derivative is applied to the diffusion equation to incorporate the history of fluid flow in naturally fractured reservoirs. The resulting fractally fractional diffusion (FFD) equation is solved analytically in the Laplace space for three outer boundary conditions. The analytical solutions are used to analyze the response of a naturally fractured reservoir considering the anomalous behavior of oil production. Several synthetic examples are provided to illustrate the methodology proposed in this work and to explain the diffusion process in fractally fractured systems.

Hydrogeophysical data integration at larger scales: Application of Bayesian sequential simulation for the characterization of heterogeneous alluvial aquifers

Knowledge of the spatial distribution of hydraulic conductivity (K) within an aquifer is critical for reliable predictions of solute transport and the development of effective groundwater management and/or remediation strategies. While core analyses and hydraulic logging can provide highly detailed information, such information is inherently localized around boreholes that tend to be sparsely distributed throughout the aquifer volume. Conversely, larger-scale hydraulic experiments like pumping and tracer tests provide relatively low-resolution estimates of K in the investigated subsurface region. As a result, traditional hydrogeological measurement techniques contain a gap in terms of spatial resolution and coverage, and they are often alone inadequate for characterizing heterogeneous aquifers. Geophysical methods have the potential to bridge this gap. The recent increased interest in the application of geophysical methods to hydrogeological problems is clearly evidenced by the formation and rapid growth of the domain of hydrogeophysics over the past decade (e.g., Rubin and Hubbard, 2005).

The surficial geology, sedimentology and geochemistry of the late glacial sediments and Paleozoic bedrock in the Campbellford area, Ontario, with special reference to the Dummer Complex /

The Dummer Complex extends 180 km along the Precambrian - Paleozoic contact from Tamworth to Lake Simcoe. It is composed of coarse, angular Paleozoic clasts in discontinuous, pitted, hummocky deposits. Deposits are usually separated by bare or boulder strewn bedrock, but have been found in the southern drumlinized till sheet. Dummer Complex deposits show rough alignment with ice-flow. Eskers cross-cut many of the deposits. Dummer sediment subfacies are defined on the basis of dominant coarse grain size and lithology, which relate directly to the underlying Paleozoic formation. Three subglacial tills are identified based on the degree of comminution and distance of transport; the immature facies of the Dummer Complex; the mature facies of the drumlinized till sheet and; the submature facies which is transitional. Carbonate geochemistry was used for till-bedrock correlation in various grain sizes. Of the 3 Paleozoic formations underlying the Dummer Complex, the Gull River Fm. is geochemically distinctive from the Bobcaygeon and Verulam Formations using Ca, Mg, Sr, Cu, Mn, Fe and Na. The Bobcaygeon Fm. and Verulam Fm. can be differentiated using Ca and the Sr/Ca ratio. The immature facies from 1.0 phi and finer is dominated by the non-carbonate, long distance transported component which decreases slightly downice. The submature till facies contains more long distance material than the immature facies. Sr and Mn can be used to correlate the Gull River immature till facies to the underlying bedrock the other subfacies could not be distinguished from each other or their respective source formation. This method proved to be ineffective for sediments with greater than 35% non-carbonate component, due to leaching of elements by the dissolving acid.The Dummer Complex is produced subglacially , as the compressional ice encounters the permeable Paleozoic carbonates. The increased shear strength of the ice and pore pressures in the carbonates results in the basal ice zones becoming debris ladden. Cleaner ice overrides the basal debris . laden dead ice which then acts as the glacier bed. During retreat, the Simcoe lobe stagnates as flow is cut-off by the Algonquin Highlands.

A hypothesis for the geochemical anomaly in the North Creek Watershed /

A regional geochemical reconnaissance by bottom stream sediment sampling, has delineated an area of high metal content in the north central sector of the North Creek Watershed. Development of a geochemical model, relating to the relative chemical concentrations derived from the chemical analyses of bottom sediments, suspended sediments, stream waters and well waters collected from the north central sector, was designed to discover the source of the anomaly. Samples of each type of material were analysed by the A.R.L. Direct Reading Multi-element Emission Spectrograph Q.A. 137 for elements: Na, K, Ca, Sr, Si, As, Pb, Zn, Cd, Ni, Ti, Ag, Mo, Be, Fe, AI, Mn, Cu, Cr, P and Y. Anomalous results led to the discovery of a spring, the waters of which carried high concentrations of Zn, Cd, Pb, As, Ni, Ti, Ag, Sr and Si. In addition, the spring waters had high concentrations of Na, Ca, Mg, 504 , alkalinity, N03' and low concentrations of K, Cl and NH3. Increased specific conductivity (up to 2500 ~mho/cm.) was noted in the spring waters as well as increased calculated total dissolved solids (up to 2047 mg/l) and increased ionic strength (up to 0.06). On the other hand, decreases were noted in water temperature (8°C), pH (pH 7.2) and Eh (+.154 volts). Piezometer nests were installed in the anomalous north central sector of the watershed. In accordance with the slope of the piezometric surface from wells cased down to the till/bedrock interface, groundwater flow is directed from the recharge area (northwest of the anomaly) towards the artesian spring via the highly fractured dolostone aquifer of the Upper Eramosa Member. The bedrock aquifer is confined by the overlying Halton till and the underlying Lower Eramosa Member (Vinemount Shale). The oxidation of sphalerite and galena and the dissolution of gypsum, celestite, calcite, and dolomite within the Eramosa Member, contributed its highly, dissolved constituents to the circulating groundwaters, the age of which is greater than 20 years as determined by tritium dating. Groundwater is assumed to flow along the Vinemount Shale and discharge as an artesian spring where the shale unit becomes discontinuous. The anomaly is located on a topographic low where bedrock is close to the surface. Thermodynamic evaluation of the major ion speciation from the anomalous spring and surface waters, showed gypsum to be supersaturated in these spring waters. Downstream from the spring, the loss of carbon dioxide from the spring waters resulted in the supersaturation with respect to calcite, aragonite, magnesite and dolomite. This corresponded with increases in Eh (+.304 volts) and pH (pH 8.5) in the anomalous surface waters. In conclusion, the interaction of groundwaters within the highly, mineralized carbonate source (Eramosa Member) resulted in the characteristic Ca*Mg*HC03*S04 spring water at the anomalous site, which appeared to be the principle effect upon controlling the anomalous surface water chemistry.

Hydraulic and fluvial geomorphological models for a bedrock channel reach of the Twenty Mile Creek /

Bedrock channels have been considered challenging geomorphic settings for the application of numerical models. Bedrock fluvial systems exhibit boundaries that are typically less mobile than alluvial systems, yet they are still dynamic systems with a high degree of spatial and temporal variability. To understand the variability of fluvial systems, numerical models have been developed to quantify flow magnitudes and patterns as the driving force for geomorphic change. Two types of numerical model were assessed for their efficacy in examining the bedrock channel system consisting of a high gradient portion of the Twenty Mile Creek in the Niagara Region of Ontario, Canada. A one-dimensional (1-D) flow model that utilizes energy equations, HEC RAS, was used to determine velocity distributions through the study reach for the mean annual flood (MAF), the 100-year return flood and the 1,000-year return flood. A two-dimensional (2-D) flow model that makes use of Navier-Stokes equations, RMA2, was created with the same objectives. The 2-D modeling effort was not successful due to the spatial complexity of the system (high slope and high variance). The successful 1 -D model runs were further extended using very high resolution geospatial interpolations inherent to the HEC RAS extension, HEC geoRAS. The modeled velocity data then formed the basis for the creation of a geomorphological analysis that focused upon large particles (boulders) and the forces needed to mobilize them. Several existing boulders were examined by collecting detailed measurements to derive three-dimensional physical models for the application of fluid and solid mechanics to predict movement in the study reach. An imaginary unit cuboid (1 metre by 1 metre by 1 metre) boulder was also envisioned to determine the general propensity for the movement of such a boulder through the bedrock system. The efforts and findings of this study provide a standardized means for the assessment of large particle movement in a bedrock fluvial system. Further efforts may expand upon this standardization by modeling differing boulder configurations (platy boulders, etc.) at a high level of resolution.

The erring archive in Anne Carson

L’archive erronée dans l’œuvre d’Anne Carson enquête sur les effets que peuvent entraîner l’archive classique sur la poésie d’Anne Carson et révèle que le travail de cette dernière est issu de l’espace situé entre la critique et la créativité, ce qui génère ce qu’on appellera une « poétique de l’erreur ». La poésie de Carson se démarque par sa prédilection pour les accidents, les imperfections et les impondérables de la transmission. La présente dissertation émerge des attitudes critiques ambivalentes face à la dualité de l’identité de Carson, autant poète qu’universitaire, et leur offrira une réponse. Alors que l’objectif traditionnel du philologue classique est de reconstruire le sens du texte « original », l’approche poétique de Carson sape en douce les prétentions universitaires d’exactitude, de précision et de totalisation. La rencontre de Carson avec l’archive classique embrasse plutôt les bourdes, les mauvaises lectures et les erreurs de traduction inhérentes à la transmission et à la réception de traductions classiques. La poésie de Carson est ludique, sexuée et politique. Sa manière de jouer avec l’épave du passé classique torpille la patri-archive, telle que critiquée par Derrida dans Mal d’Archive ; c’est-à-dire cette archive considérée comme un point d’origine stable grâce auquel s’orienter. De plus, en remettant en question la notion de l’archive classique en tant qu’origine de la civilisation occidentale, Carson offre simultanément une critique de l’humanisme, en particulier au plan de la stabilité, du caractère mesurable et de l’autonomie de « l’homme ». L’archive, pour Carson, est ouverte, en cours et incomplète ; les manques linguistiques, chronologiques et affectifs de l’archive classique représentent ainsi des sources d’inspiration poétique. La présente dissertation étudie quatre dimensions de l’archive classique : la critique, la saphique, l’élégiaque et l’érotique. Grâce à ces coordonnées, on y établit le statut fragmentaire et fissuré du passé classique, tel que conçu par Carson. Si le fondement classique sur lequel la culture occidentale a été conçue est fissuré, qu’en est-il de la stabilité, des frontières et des catégories que sont le genre, la langue et le texte ? L’ouverture de l’archive critique de manière implicite les désirs de totalité associés au corps du texte, à la narration, à la traduction et à l’érotisme. En offrant une recension exhaustive de sa poétique, L’archive erronée dans l’œuvre d’Anne Carson tente d’analyser l’accueil hostile qu’elle a subi, contribue à renforcer la documentation sans cesse croissante dont elle fait l’objet et anticipe sa transmutation actuelle de médium et de genre, sa migration de la page à la scène.

An Integrated Hydro geological Study of the Muvattupuzha River Basin, Kerala, India

The present investigation on the Muvattupuzha river basin is an integrated approach based on hydrogeological, geophysical, hydrogeochemical parameters and the results are interpreted using satellite data. GIS also been used to combine the various spatial and non-spatial data. The salient finding of the present study are accounted below to provide a holistic picture on the groundwaters of the Muvattupuzha river basin. In the Muvattupuzha river basin the groundwaters are drawn from the weathered and fractured zones. The groundwater level fluctuations of the basin from 1992 to 2001 reveal that the water level varies between a minimum of 0.003 m and a maximum of 3.45 m. The groundwater fluctuation is affected by rainfall. Various aquifer parameters like transmissivity, storage coefficient, optimum yield, time for full recovery and specific capacity indices are analyzed. The depth to the bedrock of the basin varies widely from 1.5 to 17 mbgl. A ground water prospective map of phreatic aquifer has been prepared based on thickness of the weathered zone and low resistivity values (<500 ohm-m) and accordingly the basin is classified in three phreatic potential zones as good, moderate and poor. The groundwater of the Muvattupuzha river basin, the pH value ranges from 5.5 to 8.1, in acidic nature. Hydrochemical facies diagram reveals that most of the samples in both the seasons fall in mixing and dissolution facies and a few in static and dynamic natures. Further study is needed on impact of dykes on the occurrence and movement of groundwater, impact of seapages from irrigation canals on the groundwater quality and resources of this basin, and influence of inter-basin transfer of surface water on groundwater.

A comparison of cross-hole electrical and seismic data in fractured rock

Cross-hole anisotropic electrical and seismic tomograms of fractured metamorphic rock have been obtained at a test site where extensive hydrological data were available. A strong correlation between electrical resistivity anisotropy and seismic compressional-wave velocity anisotropy has been observed. Analysis of core samples from the site reveal that the shale-rich rocks have fabric-related average velocity anisotropy of between 10% and 30%. The cross-hole seismic data are consistent with these values, indicating that observed anisotropy might be principally due to the inherent rock fabric rather than to the aligned sets of open fractures. One region with velocity anisotropy greater than 30% has been modelled as aligned open fractures within an anisotropic rock matrix and this model is consistent with available fracture density and hydraulic transmissivity data from the boreholes and the cross-hole resistivity tomography data. However, in general the study highlights the uncertainties that can arise, due to the relative influence of rock fabric and fluid-filled fractures, when using geophysical techniques for hydrological investigations.