Дисертації з теми "Electromagnetic induction geophysic"

Numerical simulations were used to study spatial averaging in low-induction-number frequency-domain electromagnetic induction (LIN FEM) instruments. Local ( LS) and cumulative (CS) sensitivity were used to analyze three different aspects of LIN FEM spatial sensitivity. LS is the variation in a measured property given a small change at a given location of the property of interest. CS contours are derived from LS and reveal the shape and the fraction of total instrument sensitivity enclosed within the contours. The first study re-evaluated the asymptotic approach to LIN FEM spatial sensitivity. Using this approach, LIN FEM measurements have often been assumed to represent electrical conductivity (sigma) at discreet depths that do not vary with the sigma of the ground. This assumption was tested using simulations of electromagnetic fields in environments with homogeneous and layered sigma distributions. When the induction number was greater than 0.01, the 1-D vertical CS distribution and the depth of investigation varied up to 20% over the range of sigma simulated. As sigma increased, CS contours and depth of investigation decreased in depth. In the second study a small perturbation approach was used to calculate CS distributions so that each distribution is unique to a given LS distribution. CS was summed from regions of high to low LS, and retained information on the magnitude and location of LS. As sigma increased, CS became focused around the highest LS values. The maximum reduction in depth of investigation was about 40% at the highest sigma investigated. In the final study, a series of small, electrically conductive perturbations was simulated in a three-dimensional, homogeneous environment. Three-dimensional LS varied markedly with a large difference between horizontal (HMD) and vertical (VMD) orientations of the transmitter and receiver dipoles. In some regions, the calculated magnetic field intensity with the perturbation was less than that calculated for the host without the perturbation. This occurred for both VMD and HMD orientations of the transmitter. CS contours were highly complex. One dimensional, vertical LS curves extracted from the three-dimensional data were very different from curves from infinite layer simulations.

The solution of the secondary coupled-vector potential formulation of Maxwell??s equations governing the controlled-source electromagnetic (CSEM) response of an arbitrary, threedimensionalconductivitymodelmust be calculatednumerically.The finite elementmethod is attractive, because it allows the model to be discretized into an unstructured mesh, permitting the specification of realistic irregular conductor geometries, and permitting the mesh to be refined locally, where finer resolution is needed. The calculated results for a series ofsimple test problems, ranging from one-dimensionalscalar differentialequations to three-dimensional coupled vector equations match the known analytic solutions well, with error values several orders of magnitude smaller than the calculated values. The electromagnetic fields of a fully three-dimensional CSEM model, recovered from the potentials using the moving least squares interpolation numerical differentiation algorithm, compares well with published numerical modeling results, particularly when local refinement is applied. Multiple buried conductors in a conductive host interact via mutual induction and current flow through the host due to the dissipation of charge accumulated on the conductor boundary. The effect of this interaction varies with host conductivity, transmitter frequency, and conductor geometry, orientation, and conductivity. For three test models containingtwo highly conductive plate-like targets, oriented in various geometries (parallel, perpendicular, and horizontal), mutual coupling ranges as high as twenty times the total magnetic field. The effect of varying host conductivity is significant, especially at high frequencies. Numerical modeling also shows that the vorticity of the currents density induced in a vertically oriented plate-like conductor rotates from vertical at high frequencies, to horizontal at low frequencies, a phenomenon confirmed by comparison with time domain field data collected in Brazos County, Texas. Furthermore, the effect of the presence of a simple horst on the CSEM response of a homogeneous conductive earth is significant, even when the height of the horst is only a fraction of the skin depth of the model. When the transmitter is placedon topofthe horst, the currents inducedtherein account for nearly all of the total magnetic field of the model, indicating that topography, like mutual coupling must be accounted for when interpreting CSEM data.

Electromagnetic induction (EM) methods have been utilised in a recent surge of archaeological applications across continental Europe, Ireland and Scandinavia. Development of multi-exploration depth instruments and improvements to instrument stability have improved its reputation as an effective method for mapping archaeological remains. Despite these advances, EM methods are comparatively lacking in rigour when for British sites. Through a structured scheme of experimental analysis and fieldwork, this thesis develops an understanding of the responses of EM instruments over a range of British archaeology, including earthworks, field systems, burials, modern remains, and a Cistercian abbey; the results of which demonstrate its effective over a diversity of environments. The impact of instrument-based issues on the collected measurements was quantified through a scheme of experiments targeting instrument drift, calibration and elevation. Dedicated instrument operation and processing workflows were developed based on the collective field and experimental results, which recommend best practice guidelines for improving the quality and accuracy of collected data. The link between instrument measurements and buried archaeology was further developed through a structured analysis of the EM datasets with complementary earth resistance and magnetic results. The integration of the EM, earth resistance and magnetic datasets was utilised to develop an enhanced archaeological characterisation of subsurface features. While the earth resistance and magnetic methods generally responded to different aspects of the buried archaeology, the EM surveys were able to detect a range of responses evident in the results of the former methods. Therefore, the role of EM methods within this characterisation are shown to “bridge the gap” between the earth resistance and magnetic methods, while providing a comprehensive characterisation of the remains in their own right.

An integrated electromagnetic (EM) and seismic geophysical study was performed to evaluate non-invasive approaches to estimate depth to shallow groundwater in arid environments with elevated soil salinity where the installation of piezometers would be impractical or prohibited. Both methods were tested in two study areas (semi-arid and arid respectively), one in Palmyra, Utah, USA near the shore of Utah Lake where groundwater is shallow and unconfined in relatively homogeneous lacustrine sediments. The other area is Carson Slough, Nevada, USA near Ash Meadows National Wildlife Refuge in Amargosa Valley. The area is underlain by valley fill, with generally variable shallow depths to water in an ephemeral braided stream environment. The methods used include frequency domain electromagnetic induction allowing for multiple antenna-receiver spacings. High resolution compressional P-wave seismic profiles using a short (0.305 m) geophone spacing for common depth-point reflection stacking and first arrival modeling were also acquired. Both methods were deployed over several profiles where shallow piezometer control was present. The semi-arid Palmyra site with its simpler geohydrology serves as an independent calibration to be compared to the Carson Slough Site. EM results at both sites show that water surfaces correspond with a drop in conductivity. This is due to elevated concentrations of evaporative salts in the vadose zone immediately above the water table. EM and seismic profiles at the Palmyra site were readily correlated to depth to groundwater in monitoring wells demonstrating that the method is ideal under laterally homogeneous conditions. Interpreting the EM and seismic profiles at Carson Slough was challenging due to the laterally and vertically variable soil types, segmented perched water surfaces, and strong salinity variations. The high-resolution images and models provided by the seismic profiles confirm the simple soil and hydrological structure at the Palmyra site as well as the laterally complex structure at Carson Slough. The EM and seismic results indicate that an integrated geophysical approach is necessary for an area like Carson Slough, where continued leaching of salts combined with braided stream deposition has created a geophysically complex soil and groundwater system.

Geology
M.S.
Vacant lots in cities and surrounding urban areas can potentially be used for stormwater management because they are pervious. However, the extent to which vacant lots provide pervious cover to increase infiltration and reduce stormflow is poorly understood. The goal of this study was to develop faster methods for monitoring stormwater infiltration to improve characterization of heterogeneous urban systems. Geophysical techniques are capable of mapping and characterizing subsurface materials, but are often limited by time and sensitivity constraints. In this study, the infiltration characteristics of a vacant lot created by the demolition of a house was characterized using a series of modeling, field and lab experiments. Site characterization under background conditions with an EM Profiler was used to map zones of different fill materials. Three zones were identified in the study site: grass area, driveway area, and a former house area. Transient soil moisture conditions were monitored during irrigation tests using two geophysical methods (electrical resistivity tomography [ERT] and electromagnetic induction [EM]) to evaluate method sensitivity and differences between the three zones. ERT proved more sensitive than EM profiling at detecting changes in the three zones. Soil moisture changes in the driveway area were particularly difficult to detect using EM. The EM Profiler showed a reduction rather than increase in conductivity at the start of irrigation and storms, which was attributed to flushing of high conductivity pore fluids by dilute irrigation or rain water. This explanation was supported using Archie’s Law to model the response of apparent conductivity under highly conductive pore fluid conditions. The EM Profiler was also used under natural precipitation conditions to quickly monitor areas too large for the ERT to reasonably survey. The results suggested that EM instrument drift needs to be corrected to make the method more sensitive. It was difficult to detect differences in hydrologic characterization between areas of the vacant lot using traditional soil point measurements because of the inherent spatial variability. The most useful point measurement was soil moisture loggers. Data from soil moisture loggers was used to parameterize the model; in addition, the soil moisture loggers showed a slow drying period. By combining the EM Profiler method with soil moisture data and applying corrections for drift, some improvement in sensitivity might be achieved. Quantitative characterization of fill material was shown by ERT, which detected more heterogeneous infiltration in the area of the former house than in the grass area.
Temple University--Theses

Runion is a protohistoric Native American village located on the floodplain of the Nolichucky River in western Washington County. Previous archaeological excavations and radiocarbon dates suggest that the village was occupied during the mid-16th to mid-17th century. The Nolichucky River, in contrast, has been flowing through the area for millennia. Geophysical surveys are used to image the subsurface non-invasively, without disturbing protected land and/or organisms. Preliminary geophysical data collected at Runion include ground penetrating radar (GPR), electromagnetic induction (EMI), and magnetometry. These data show a linear feature surrounding the protohistoric village. Given its placement around the margins of the village, the feature could be interpreted as a fortification ditch, which is often paired with a palisade wall to defend a village from attack. The feature is also consistent with typical meandering floodplain stratigraphy, where sections of channel are often abandoned to form oxbow lakes. Over time these abandoned channels fill in and are called paleochannels. Each geophysical method measures the properties and characteristics of the linear feature, a presumed paleochannel. GPR sends electromagnetic radar waves into the ground, which reflect off different subsurface layers and are recorded as radargrams. Magnetometry measures subtle changes in earth magnetism, including the magnetization of rocks, soils, and/or ferrous objects. EMI systems transmit low frequency electromagnetic waves to measure both electrical conductivity (EC) and magnetic susceptibility (MS). Each of these instruments are used to collect data in transects and then processed to produce profiles, maps and, in the case of GPR, three-dimensional datasets of the subsurface. It is anticipated that GPR will reveal details about the stratigraphy of the linear feature. Magnetic, EC, and MS measurements will further help to interpret the GPR data by distinguishing between different types of sediments. These data may show if the feature is a paleochannel or a ditch excavated into older stratigraphic layers by village inhabitants for fortification. Ultimately, the feature will be tested with soil cores to study the sediments directly. At this preliminary stage the feature is interpreted to be a paleochannel. The stratigraphic layers revealed by GPR show a broad depression with stratigraphic layers characteristic of a paleochannel. In addition, magnetic readings are anomalously low on the eastern margin (closer to the modern river channel) and high on the western margin. This could indicate paired point bar sands and paleochannel fill, respectively. This interpretation is still tentative, however, because we have not yet integrated the EMI data, extracted soil cores, or dated the feature. Radiocarbon dates might help determine the relative age of the feature if organic carbon is present. In conclusion, preliminary data currently suggests that the structure is geological rather than archaeological. In the coming months we will collect more GPR data with different frequency GPR antennas, integrate the EMI data, and test the findings by extracting soil cores and reconstructing the stratigraphy.

This dissertation focuses on two different types of responses of Earth; that is, seismic and electromagnetic, and aims to better understand Earth processes at a wider range of scales than those conventional approaches offer. Electromagnetic responses resulting from the subsurface diffusion of applied electromagnetic fields through heterogeneous geoelectrical structures are utilized to characterize the underlying geology. Geology exhibits multiscale hierarchical structure which brought about by almost all geological processes operating across multiple length scales and the relationship between multiscale electrical properties of underlying geology and the observed electromagnetic response has not yet been fully understood. To quantify this relationship, the electromagnetic responses of textured and spatially correlated, stochastic geologic media are herein presented. The modelling results demonstrate that the resulting electromagnetic responses present a power law distribution, rather than a smooth response polluted with random, incoherent noise as commonly assumed; moreover, they are examples of fractional Brownian motion. Furthermore, the results indicate that the fractal behavior of electromagnetic responses is correlated with the degree of the spatial correlation, the contrasts in ground electrical conductivity, and the preferred orientation of small-scale heterogeneity. In addition, these inferences are also supported by the observed electromagnetic responses from a fault zone comprising different lithological units and varying wavelengths of geologic heterogeneity. Seismic signals generated by aftershocks are generally recorded by local aftershock networks consisted of insufficient number of stations which result in strongly spatially-aliased aftershock data. This limits aftershock detections and locations at smaller magnitudes. Following the 23 August 2011 Mineral, Virginia earthquake, to drastically reduce spatial aliasing, a temporary dense array (AIDA) consisting of ~200 stations at 200-400 m spacing was deployed near the epicenter to record the 12 days of the aftershocks. The backprojection imaging method is applied to the entire AIDA dataset to detect and locate aftershocks. The method takes advantage of staking of many seismograms and improves the signal-to-noise ratio for detection. The catalog obtained from the co-deployed, unusually large temporal traditional network of 36 stations enabled a quantitative comparison. The aftershock catalog derived from the dense AIDA array and the backprojection indicates event detection an order of magnitude smaller including events as small as M–1.8. The catalog is complete to magnitude –1.0 while the traditional network catalog was complete to M–0.27 for the same time period. The AIDA backprojection catalog indicate the same major patterns of seismicity in the epicentral region, but additional details are revealed indicating a more complex fault zone and a new shallow cluster. The b-value or the temporal decay constant were not changed by inclusion of the small events; however, they are different for two completeness periods and are different at shallow depth than greater depth.
Ph. D.

The need to consider groundwater and surface water as a single resource has fostered the interest of the scientific community on the interactions between surface water and groundwater. The region below and alongside rivers where surface hydrology and subsurface hydrology concur is the hyporheic zone. This is the region where water exchange determines many biogeochemical and ecological processes of great impact on the functioning of rivers. However, the complex processes taking place in the hyporheic zone require a multidisciplinary approach. The combination of innovative point and distributed techniques originally developed in separated disciplines is of great advantage for the indirect identification of water exchange in the hyporheic zone. Distributed techniques using temperature as a tracer such as fiber-optic distributed temperature sensing can identify the different components of groundwater-surface water interactions based on their spatial and temporal thermal patterns at the sediment-water interface. In particular, groundwater, interflow discharge and local hyporheic exchange flows can be differentiated based on the distinct size, duration and sign of the temperature anomalies. The scale range and resolution of fiber-optic distributed temperature sensing are well complemented by geophysics providing subsurface structures with a similar resolution and scale. Thus, the use of fiber-optic distributed temperature sensing to trace flux patterns supported by the exploration of subsurface structures with geophysics enables spatial and temporal investigation of groundwater-surface water interactions with an unprecedented level of accuracy and resolution. In contrast to the aforementioned methods that can be used for pattern identification at the interface, other methods such as point techniques are required to quantify hyporheic exchange fluxes. In the present PhD thesis, point methods based on hydraulic gradients and thermal profiles are used to quantify hyporheic exchange flows. However, both methods are one-dimensional methods and assume that only vertical flow occurs while the reality is much more complex. The study evaluates the accuracy of the available methods and the factors that impact their reliability. The applied methods allow not only to quantify hyporheic exchange flows but they are also the basis for an interpretation of the sediment layering in the hyporheic zone. For upscaling of the previous results three-dimensional modelling of flow and heat transport in the hyporheic zone combines pattern identification and quantification of fluxes into a single framework. Modelling can evaluate the influence of factors governing groundwater-surface water interactions as well as assess the impact of multiple aspects of model design and calibration of high impact on the reliability of the simulations. But more importantly, this modelling approach enables accurate estimation of water exchange at any location of the domain with unparalleled resolution. Despite the challenges in 3D modelling of the hyporheic zone and in the integration of point and distributed data in models, the benefits should encourage the hyporheic community to adopt an integrative approach comprising from the measurement to the upscaling of hyporheic processes.

This dissertation consists of three projects utilizing electric and electromagnetic (EM) methods to better understand critical-zone hydrogeologic conditions in select Florida wetlands and waters. First, a time-lapse electrical resistivity (ER) survey was conducted in section of mangrove forest on a barrier island in southeast Florida to image changes in pore-water salinity in the root zone. ER data show the most variability in the root zone over a 24-hour period, and, generally, the ground is more resistive during the day than overnight. Second, a suite of three-dimensional forward models, based on varying lateral boundaries and conductivities typical of a coastal wetland, were run to simulate the EM response of a commerical electromagnetic induction instrument crossing over said boundaries. Normalized profiles show the transition is sharper in a hypersaline regime than one where freshwater and clay are present. Furthermore, enough variability exists in hypersaline regimes to justify collecting profile measurements in multiple coil configurations to constrain the nature of a lateral boundary. Also, under certain circumstances, there are kinks in the EMI response even across abrupt boundaries due to concentrated current density at a layer's edge. Lastly, geophysical surveys were conducted at six wetlands in west-central Florida to characterize potential hydrostratigraphic units and compare/contrast them to the current conceptual model for cypress dome wetlands. ER was used to image the geometry of the top of limestone; ground penetrating radar (GPR) was used to image stratigraphy beneath and surrounding wetlands. These wetlands can be grouped into two models. Topographic highs surrounding wetlands are controlled by the undulating top of limestone at sites where the region is characterized by limestone ridges. In contrast, topographic highs are controlled by thick sand packages at sites regionally characterized by sand dunes over scoured limestone.

La méthode géophysique électromagnétique inductive dans le domaine temporel (TDEM) est employée dans des domaines tels que l’hydrologie ou l’exploration minière et permet l’établissement de modèles de résistivité électrique du sous-sol. L’application de cette méthode avec des dispositifs d’acquisition réduits, avec des boucles de quelques mètres de côté, rend la mesure plus sensible aux autres paramètres électromagnétiques : la permittivité diélectrique et la viscosité magnétique. L’objet du travail de recherche présenté ici est de réaliser des sondages TDEM avec des dispositifs d’acquisition réduits et d’en extraire non seulement un modèle de résistivité électrique du sous-sol, mais aussi de permittivité diélectrique et de viscosité magnétique. L’acquisition d’un sondage TDEM avec des dispositifs à petites boucles implique également une sensibilité à la réponse du système sur lui-même, plus importante. Des simulations de cette réponse ont été menées afin de pouvoir l’estimer le plus précisément possible et la prendre en compte lors du traitement des données. Une géométrie d’acquisition des mesures est ensuite proposée pour rendre l’interprétation multi-paramètres de la mesure plus accessible. Les effets des propriétés électromagnétiques se superposent sur la mesure en configuration centrale. L’utilisation de géométries avec une sensibilité plus faible à la permittivité diélectrique et la viscosité magnétique permettra de séparer les contributions des différents paramètres. Ces développements méthodologiques ont enfin été testés lors d’acquisitions sur le site test de Garchy (Nièvre)
The geophysical electromagnetic inductive method in the time domain (TDEM) is used in fields such as hydrology or in mining exploration. It allows the establishment of electrical resistivity models of the subsurface. Using this methods with reduced acquisition loops with a width of only a few meters, makes the measurement more sensitive to additional parameters : the dielectric permittivity and the magnetic permeability. The present research work aims to acquire a small-loop TDEM sounding and to extract not only an electical resistivity model but also an estimation of the dielctric permittivity and of the magnetic viscosity. The acquisition of a small-loop TDEM sounding implies a increased sensitivity to the response from the acquisition system itself, in part because of the reduced distance between the transmission and reception loops. Simulations aiming to reproduce this system response have been carried out to be able to take it into account during the data processing step. A specific acquisition geometry is then proposed to ease the multi-parameter interpretation of the TDEM data. The effects from the electric, dielectric and magnetic properties of the ground are overlapping on the measurements in the central configuration. Using different geometry with lower sensivities to dielectric permittivity and magnetic viscosity should make the contribution of each parameter easier to isolate. These methological developments were then tested with field measurement on the test site of Garchy (Nièvre, France)

The Pimpama estuarine plain in subtropical southeast Queensland is comprised of Quaternary sediments infilling older bedrock. These multilayered unconsolidated sediments have various depositional origins, and are highly heterogeneous. The plain is low-lying and the surface drainage is controlled by flood mitigation measures including tidal gates and channelised streams. The control of surface drainage potentially affects the shallow water table. This modification of hydrology has implications for future viability of agriculture and also the environmental health of waterways. Increased landscape modification and water management is likely in the coming years. The combination of sediment heterogeneity, low hydraulic gradients, and artificial drainage modification result in the plain being hydrogeologically complex. In order to understand hydrologic processes in this setting, a multi-disciplinary research programme was conducted which included a drilling program, overland electromagnetic induction and other geophysical surveys (downhole gamma log, electromagnetic induction and magnetic susceptibility) to initially establish the geologic framework. These surveys were followed by hydrogeochemical testing which includes for major and minor ions and also stable isotopes, and mineralogical analysis of drillhole material. Underlying basement rock occurs at up to 60 m depth. Unconsolidated gravel and sand deposits occur within incised paleo-valleys and are overlain by predominantly low-permeability fluvial sandy clays and estuarine and lagoonal muds. Fine-grained delta sands occur in the top 15 m of the sub-surface. Within the unconsolidated sediments, hydrodynamic trends clearly discriminated between upper unconfined and lower semi-confined aquifer systems. A comparison of surface water and shallow groundwater levels indicate limited interaction of groundwater and surface water. Hydrogeochemical analysis effectively distinguished between groundwater bodies, and also distinguished saline groundwater from seawater. Trends in major ion chemistry in the semi-confined system (particularly Na/Cl and Ca/Cl ratios) showed ion exchange accompanying saline intrusion. However, due to factors such as mineral dissolution, major ion chemistry does not clearly identify solute flux trends in the shallow aquifer system. Water stable isotope analysis (δ18O and δ2H) indicated the provenance of fresh and saline groundwater and also the relative importance of the principal hydrologic processes, i.e. evaporation and water uptake by plants. Groundwater exhibited a wide range in salinity, from very fresh to hypersaline. The formation of hypersaline groundwater was attributed largely to uptake of water by mangrove forests. Since mangrove forests were more extensive at the time of the Holocene maximum sea level (approximately 6,000 years ago) than at present, some of this groundwater may represent relict salinity from this earlier time. The relationship of relict salinity to low permeability sediments, particularly at intermediate depths, and their depositional history was examined. Vertical salinity gradients and hydrogeochemistry within these sediments varied according to position within the plain, suggesting deposition under various hydrological and sea level regimes. A preliminary investigation using analysis of stable sulfate isotopes (δ34S and δ18OSO4) was made. This study shows substantial potential for the application of this technique for quantification of solute flux and sulfur chemical transformations within settings such as this coastal plain. To establish shallow groundwater flow processes, a MODFLOW-based numerical model was used to inversely estimate aquifer parameters under various recharge scenarios. The model was designed to examine the relative importance of evapotranspiration and discharge to surface waters. However, largely due to the complexity of the drainage network and non-uniform surface water flows, the quantification of surface water- groundwater interaction by consideration of hydrodynamics is problematic. Therefore, the chemistry of groundwater and surface water was compared. While the estimated contribution of rainfall to groundwater level fluctuations was significant (46%), high evapotranspiration rates reduced net recharge and it was concluded that baseflow to drains and creeks during dry periods was insignificant, and groundwater velocities in the shallow aquifer are low. The study illustrates the value of both hydrodynamic and hydrogeochemical analyses in estuarine settings where relict salinity and groundwater-aquifer interactions impact significantly on water quality. Saline groundwater is chemically distinct from theoretical mixtures of seawater and freshwater. The study also demonstrates the value of particular chemical parameters, e.g. Na/Cl and SO4/Cl ratios and stable water isotopes, for identifying hydrologic processes in this setting.

The Pimpama estuarine plain in subtropical southeast Queensland is comprised of Quaternary sediments infilling older bedrock. These multilayered unconsolidated sediments have various depositional origins, and are highly heterogeneous. The plain is low-lying and the surface drainage is controlled by flood mitigation measures including tidal gates and channelised streams. The control of surface drainage potentially affects the shallow water table. This modification of hydrology has implications for future viability of agriculture and also the environmental health of waterways. Increased landscape modification and water management is likely in the coming years. The combination of sediment heterogeneity, low hydraulic gradients, and artificial drainage modification result in the plain being hydrogeologically complex. In order to understand hydrologic processes in this setting, a multi-disciplinary research programme was conducted which included a drilling program, overland electromagnetic induction and other geophysical surveys (downhole gamma log, electromagnetic induction and magnetic susceptibility) to initially establish the geologic framework. These surveys were followed by hydrogeochemical testing which includes for major and minor ions and also stable isotopes, and mineralogical analysis of drillhole material. Underlying basement rock occurs at up to 60 m depth. Unconsolidated gravel and sand deposits occur within incised paleo-valleys and are overlain by predominantly low-permeability fluvial sandy clays and estuarine and lagoonal muds. Fine-grained delta sands occur in the top 15 m of the sub-surface. Within the unconsolidated sediments, hydrodynamic trends clearly discriminated between upper unconfined and lower semi-confined aquifer systems. A comparison of surface water and shallow groundwater levels indicate limited interaction of groundwater and surface water. Hydrogeochemical analysis effectively distinguished between groundwater bodies, and also distinguished saline groundwater from seawater. Trends in major ion chemistry in the semi-confined system (particularly Na/Cl and Ca/Cl ratios) showed ion exchange accompanying saline intrusion. However, due to factors such as mineral dissolution, major ion chemistry does not clearly identify solute flux trends in the shallow aquifer system. Water stable isotope analysis (δ18O and δ2H) indicated the provenance of fresh and saline groundwater and also the relative importance of the principal hydrologic processes, i.e. evaporation and water uptake by plants. Groundwater exhibited a wide range in salinity, from very fresh to hypersaline. The formation of hypersaline groundwater was attributed largely to uptake of water by mangrove forests. Since mangrove forests were more extensive at the time of the Holocene maximum sea level (approximately 6,000 years ago) than at present, some of this groundwater may represent relict salinity from this earlier time. The relationship of relict salinity to low permeability sediments, particularly at intermediate depths, and their depositional history was examined. Vertical salinity gradients and hydrogeochemistry within these sediments varied according to position within the plain, suggesting deposition under various hydrological and sea level regimes. A preliminary investigation using analysis of stable sulfate isotopes (δ34S and δ18OSO4) was made. This study shows substantial potential for the application of this technique for quantification of solute flux and sulfur chemical transformations within settings such as this coastal plain. To establish shallow groundwater flow processes, a MODFLOW-based numerical model was used to inversely estimate aquifer parameters under various recharge scenarios. The model was designed to examine the relative importance of evapotranspiration and discharge to surface waters. However, largely due to the complexity of the drainage network and non-uniform surface water flows, the quantification of surface water- groundwater interaction by consideration of hydrodynamics is problematic. Therefore, the chemistry of groundwater and surface water was compared. While the estimated contribution of rainfall to groundwater level fluctuations was significant (46%), high evapotranspiration rates reduced net recharge and it was concluded that baseflow to drains and creeks during dry periods was insignificant, and groundwater velocities in the shallow aquifer are low. The study illustrates the value of both hydrodynamic and hydrogeochemical analyses in estuarine settings where relict salinity and groundwater-aquifer interactions impact significantly on water quality. Saline groundwater is chemically distinct from theoretical mixtures of seawater and freshwater. The study also demonstrates the value of particular chemical parameters, e.g. Na/Cl and SO4/Cl ratios and stable water isotopes, for identifying hydrologic processes in this setting.

The need to consider groundwater and surface water as a single resource has fostered the interest of the scientific community on the interactions between surface water and groundwater. The region below and alongside rivers where surface hydrology and subsurface hydrology concur is the hyporheic zone. This is the region where water exchange determines many biogeochemical and ecological processes of great impact on the functioning of rivers. However, the complex processes taking place in the hyporheic zone require a multidisciplinary approach. The combination of innovative point and distributed techniques originally developed in separated disciplines is of great advantage for the indirect identification of water exchange in the hyporheic zone. Distributed techniques using temperature as a tracer such as fiber-optic distributed temperature sensing can identify the different components of groundwater-surface water interactions based on their spatial and temporal thermal patterns at the sediment-water interface. In particular, groundwater, interflow discharge and local hyporheic exchange flows can be differentiated based on the distinct size, duration and sign of the temperature anomalies. The scale range and resolution of fiber-optic distributed temperature sensing are well complemented by geophysics providing subsurface structures with a similar resolution and scale. Thus, the use of fiber-optic distributed temperature sensing to trace flux patterns supported by the exploration of subsurface structures with geophysics enables spatial and temporal investigation of groundwater-surface water interactions with an unprecedented level of accuracy and resolution. In contrast to the aforementioned methods that can be used for pattern identification at the interface, other methods such as point techniques are required to quantify hyporheic exchange fluxes. In the present PhD thesis, point methods based on hydraulic gradients and thermal profiles are used to quantify hyporheic exchange flows. However, both methods are one-dimensional methods and assume that only vertical flow occurs while the reality is much more complex. The study evaluates the accuracy of the available methods and the factors that impact their reliability. The applied methods allow not only to quantify hyporheic exchange flows but they are also the basis for an interpretation of the sediment layering in the hyporheic zone. For upscaling of the previous results three-dimensional modelling of flow and heat transport in the hyporheic zone combines pattern identification and quantification of fluxes into a single framework. Modelling can evaluate the influence of factors governing groundwater-surface water interactions as well as assess the impact of multiple aspects of model design and calibration of high impact on the reliability of the simulations. But more importantly, this modelling approach enables accurate estimation of water exchange at any location of the domain with unparalleled resolution. Despite the challenges in 3D modelling of the hyporheic zone and in the integration of point and distributed data in models, the benefits should encourage the hyporheic community to adopt an integrative approach comprising from the measurement to the upscaling of hyporheic processes.

A significant problem with exploring for electrically conductive mineral deposits with airborne electromagnetic (AEM) methods is that many of the most valuable sulphide deposits are too conductive to be detected with conventional systems. High-grade sulphide deposits with bulk electrical conductivities on the order of 100,000 S/m can appear as “perfect conductors” to most EM systems because the decay of secondary fields (the “time constant” of the deposit) generated in the target by the system transmitter takes much longer than the short measuring time of EM systems. Their EM response is essentially undetectable with off-time measurements. One solution is to make measurements during the transmitter on-time when the secondary field of the target produced by magnetic flux exclusion is large. The difficulty is that the secondary field must be measured in the presence of a primary field which is orders of magnitude larger. The goal of this thesis is to advance the methodology of making AEM measurements during transmitter on-time by analysing experimental data from three different AEM systems. The first system analysed is a very large separation, two helicopter system where geometry is measured using GPS sensors. In order to calculate the primary field at the receiver with sufficient accuracy, the very large (nominally 400 m) separation requires geometry to be known to better than 1 m. Using the measured geometry to estimate and remove the primary field, I show that a very conductive target can be detected at depths of 200m using the total secondary field. I then used fluxgate magnetometers to correct for receiver rotation which allowed the component of the secondary field to be determined. The second system I examined was a large separation fixed-wing AEM system. Using a towed receiver bird with a smaller (˜ 135m) separation, the geometry must be known much more accurately. In the absence of direct measurement of this geometry, I used a least-squares prediction approach using measurements of aircraft manoeuvres which allowed primary field contamination to be estimated. Subtracting this estimate from the recorded signal increased the maximum time constant observed in a field survey for conductive targets by a factor of seven. Finally, a study of a nominally rigid helicopter EM system employing a bucking coil to cancel primary field showed that system geometry (specifically, the position of the receiver coil relative to the transmitter and bucking coils) must be known to better than 0.01 mm to detect deep targets. Again, direct measurements of system geometry were not available. A least-squares prediction filter using helicopter manoeuvre and system pitch and roll measurements was applied, but was not able to estimate primary field well enough to provide an accurate secondary on-time response. Direct measurements of relative motion of the system components might solve this problem.

Remediation of unexploded ordnance (UXO) and prioritization of excavation procedures for archaeological artifacts using electromagnetic (EM) induction are studied in this dissertation. Lowering of the false alarm rates that require excavation and artifact excavation prioritization can reduce the costs associated with unnecessary procedures. Data were taken over 5 areas at the San Jacinto Battleground near Houston, Texas, using an EM-63 metal detection instrument. The areas were selected using the archaeological concepts of cultural and natural formation processes applied to what is thought to be areas that were involved in the 1836 Battle of San Jacinto. Innovative use of a Statistical Point Pattern Analysis (PPA) is employed to identify clustering of EM anomalies. The K-function uses point {x,y} data to look for possible clusters in relation to other points in the data set. The clusters once identified using K-function will be further examined for classification and prioritization using the Weighted K-function. The Weighted K-function uses a third variable such as millivolt values or time decay to aid in segregation and prioritization of anomalies present. Once the anomalies of interest are identified, their locations are determined using the Gi-Statistics Technique. The Gi*-Statistic uses the individual Cartesian{x, y} points as origin locations to establish a range of distances to other cluster points in the data set. The segregation and location of anomalies supplied by this analysis will have several benefits. Prioritization of excavations will narrow down what areas should be excavated first. Anomalies of interest can be located to guide excavation procedures within the areas surveyed. Knowing what anomalies are of greater importance than others will help to lower false alarm rates for UXO remediation or for archaeological artifact selection. Knowing significant anomaly location will reduce the number of excavations which will subsequently save time and money. The procedures and analyses presented here are an interdisciplinary compilation of geophysics, archaeology and statistical analysis brought together for the first time to examine problems associated with UXO remediation as well as archaeological artifact selection at historic battlegrounds using electromagnetic data.

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

Dissertation

No
Searching for and mapping the physical extent of unmarked graves using geophysical techniques has proven difficult in many cases. The success of individual geophysical techniques for detecting graves depends on a site-by-site basis. Significantly, detection of graves often results from measured contrasts that are linked to the background soils rather than the type of archaeological feature associated with the grave. It is evident that investigation of buried remains should be considered within a 3D space as the variation in burial environment can be extremely varied through the grave. Within this paper, we demonstrate the need for a multi-method survey strategy to investigate unmarked graves, as applied at a "planned" but unmarked pauper's cemetery. The outcome from this case study provides new insights into the strategy that is required at such sites. Perhaps the most significant conclusion is that unmarked graves are best understood in terms of characterization rather than identification. In this paper, we argue for a methodological approach that, while following the current trends to use multiple techniques, is fundamentally dependent on a structured approach to the analysis of the data. The ramifications of this case study illustrate the necessity of an integrated strategy to provide a more holistic understanding of unmarked graves that may help aid in management of these unseen but important aspects of our heritage. It is concluded that the search for graves is still a current debate and one that will be solved by methodological rather than technique-based arguments.

No
In this article we assess the abilities of a new electromagnetic (EM) system, the CMD Mini-Explorer, for prospecting of archaeological features in Ireland and the UK. The Mini-Explorer is an EM probe which is primarily aimed at the environmental/geological prospecting market for the detection of pipes and geology. It has long been evident from the use of other EM devices that such an instrument might be suitable for shallow soil studies and applicable for archaeological prospecting. Of particular interest for the archaeological surveyor is the fact that the Mini-Explorer simultaneously obtains both quadrature (conductivity') and in-phase (relative to magnetic susceptibility') data from three depth levels. As the maximum depth range is probably about 1.5m, a comprehensive analysis of the subsoil within that range is possible. As with all EM devices the measurements require no contact with the ground, thereby negating the problem of high contact resistance that often besets earth resistance data during dry spells. The use of the CMD Mini-Explorer at a number of sites has demonstrated that it has the potential to detect a range of archaeological features and produces high-quality data that are comparable in quality to those obtained from standard earth resistance and magnetometer techniques. In theory the ability to measure two phenomena at three depths suggests that this type of instrument could reduce the number of poor outcomes that are the result of single measurement surveys. The high success rate reported here in the identification of buried archaeology using a multi-depth device that responds to the two most commonly mapped geophysical phenomena has implications for evaluation style surveys. Copyright (c) 2013 John Wiley & Sons, Ltd.

Yes
Archaeological research at Stonehenge (UK) is increasingly aimed at understanding the dynamic of the wider archaeological landscape. Through the application of state-of-the-art geophysical techniques, unprecedented insight is being gathered into the buried archaeological features of the area. However, applied survey techniques have rarely targeted natural soil variation, and the detailed knowledge of the palaeotopography is consequently less complete. In addition, metallic topsoil debris, scattered over different parts of the Stonehenge landscape, often impacts the interpretation of geophysical datasets. The research presented here demonstrates how a single multi-receiver electromagnetic induction (EMI) survey, conducted over a 22 ha area within the Stonehenge landscape, offers detailed insight into natural and anthropogenic soil variation at Stonehenge. The soil variations that were detected through recording the electrical and magnetic soil variability, shed light on the genesis of the landscape, and allow for a better definition of potential palaeoenvironmental and archaeological sampling locations. Based on the multi-layered dataset, a procedure was developed to remove the influence of topsoil metal from the survey data, which enabled a more straightforward identification of the detected archaeology. The results provide a robust basis for further geoarchaeological research, while potential to differentiate between modern soil disturbances and the underlying sub-surface variations can help in solving conservation and management issues. Through expanding this approach over the wider area, we aim at a fuller understanding of the human–landscape interactions that have shaped the Stonehenge landscape.