Theoretical Developments in Electromagnetic Induction Geophysics with Selected Applications in the Near Surface

Near-surface applied electromagnetic geophysics is experiencing an explosive period of growth with many innovative techniques and applications presently emergent and others certain to be forthcoming. An attempt is made here to bring together and describe some of the most notable advances. This is a difficult task since papers describing electromagnetic induction methods are widely dispersed throughout the scientific literature. The traditional topics discussed herein include modeling, inversion, heterogeneity, anisotropy, target recognition, logging, and airborne electromagnetics (EM). Several new or emerging techniques are introduced including landmine detection, biogeophysics, interferometry, shallow-water electromagnetics, radiomagnetotellurics, and airborne unexploded ordnance (UXO) discrimination. Representative case histories that illustrate the range of exciting new geoscience that has been enabled by the developing techniques are presented from important application areas such as hydrogeology, contamination, UXO and landmines, soils and agriculture, archeology, and hazards and climate.     

A sedimentary record of the rise and fall of the metal industry in Bergslagen,south central Sweden

Detailed stratigraphic analyses of sediments deposited in Lake Botjärnen, a small boreal forest lake in the shield terrain of central Sweden, clearly reflect progressively increasing human impact on terrestrial and aquatic ecosystems following settlement and establishment of an iron industry in the 17th century. Rising frequencies of pollen and spores from light-demanding plants provide evidence of extensive forest clearance for charcoal and timber production, which peaked in the early 20th century. An associated increase in catchment erosion is reflected by changing carbon and nitrogen elemental content and carbon–isotope composition of sediment organic matter and by increased magnetic susceptibility of the sediments. Records of air-borne pollutants (lead, zinc and sulphur) can be correlated to the development of local and regional mining and metal industry as inferred from historical accounts. Rapid recession of the iron industry led to re-forestation and recovery of the aquatic nutrient status to pre-industrial conditions over the past 100 years. The chronology of the sediment succession, which is based on ²¹⁰Pb and ¹³⁷Cs radionuclide data in combination with radiocarbon dating, is confirmed by historical lead pollution trends established for the region.     

High temperature gas–solid reactions in calc–silicate Cu–Au skarn formation; Ertsberg,Papua Province,Indonesia

The 2.7–3 Ma Ertsberg East Skarn System (Indonesia), adjacent to the giant Grasberg Porphyry Copper deposit, is part of the world’s largest system of Cu–Au skarn deposits. Published fluid inclusion and stable isotope data show that it formed through the flux of magma-derived fluid through contact metamorphosed carbonate rock sequences at temperatures well above 600° C and pressures of less than 50 MPa. Under these conditions, the fluid has very low density and the properties of a gas. Combining a range of micro-analytical techniques, high-resolution QEMSCAN mineral mapping and computer-assisted X-ray micro-tomography, an array of coupled gas–solid reactions may be identified that controlled reactive mass transfer through the ~ 1 km³ hydrothermal skarn system. Vacancy-driven mineral chemisorption reactions are identified as a new type of reactive transport process for high-temperature skarn alteration. These gas–solid reactions are maintained by the interaction of unsatisfied bonds on mineral surfaces and dipolar gas-phase reactants such as SO₂ and HCl that are continuously supplied through open fractures and intergranular diffusion. Principal reactions are (a) incongruent dissolution of almandine-grossular to andradite and anorthite (an alteration mineral not previously recognized at Ertsberg), and (b) sulfation of anorthite to anhydrite. These sulfation reactions also generate reduced sulfur with consequent co-deposition of metal sulfides. Diopside undergoes similar reactions with deposition of Fe-enriched pyroxene in crypto-veins and vein selvedges. The loss of calcium from contact metamorphic garnet to form vein anhydrite necessarily results in Fe-enrichment of wallrock, and does not require Fe-addition from a vein fluid as is commonly assumed.     

Assessment of drought over the past two millennia using near-shore sediment cores from a Canadian boreal lake

Lack of long-term studies on drought in the boreal region of northwest Ontario limits our ability to assess the vulnerability of this region to climate change. We investigated changes in diatoms, scaled chrysophytes, and sedimentary pigments in two near-shore cores from Gall Lake to infer limnological and water-level changes over the last two millennia. The two coring locations, at modern water depths of 7.5 and 11.5 m, were selected to contrast inferences for past lake level based on distance from the modern water-depth boundary between predominantly benthic and planktonic diatom assemblages in surface sediments (6.0 m). Diatom-inferred depth inferences were more variable in the core from 7.5-m water depth, consistent with the hypothesis that the greatest changes occurred at water depths closest to the benthic:planktonic boundary. Both sites revealed a pronounced drought from ~AD 950 to 1300, synchronous with the medieval climate anomaly (MCA). This finding suggests a northeast expansion of the arid MCA into northwest Ontario, extending the known spatial extent of this megadrought. Scaled chrysophytes and sedimentary pigments suggest a recent increase in thermal stratification. Our findings indicate this region is more susceptible to changes in moisture than was previously suspected based on the instrumental record for the past century.     

Lake ecosystem dynamics and links to climate change inferred from a stable isotope and organic palaeorecord from a mountain lake in southwestern China (ca. 22.6–10.5 cal ka BP)

A detailed understanding of long-term climatic and environmental change in southwestern China is hampered by a lack of long-term regional palaeorecords. Organic analysis (%TOC, %TN, C/N ratios and δ¹³C values) of a sediment sequence from Lake Shudu, Yunnan Province (ca. 22.6–10.5 cal ka BP) indicates generally low aquatic palaeoproductivity rates over millennial timescales in response to cold, dry climatic conditions. However, the record is punctuated by two marked phases of increased aquatic productivity from ca. 17.7 to 17.1 cal ka BP and from ca. 11.9 to 10.5 cal ka BP. We hypothesise that these shifts reflect a marked, stepwise lacustrine response to Asian summer monsoon strengthening during the last deglaciation.     

The strengthening East Australian Current, its eddies and biological effects — an introduction and overview

The poleward flowing East Australian Current (EAC) is characterised by its separation from the coast, 100-200 nautical miles north of Sydney, to form the eastward flowing Tasman Front and a southward flowing eddy field. The separation zone greatly influences coastal ecosystems for the relatively narrow continental shelf (only 15-50 km wide), particularly between 32-34°S. In this region the continental shelf has a marked shift in the seasonal temperature-salinity relationship and elevated surface nitrate concentrations. This current parallels the portion of the coast where Australia’s population is concentrated and has a long history of scientific research. However, understanding of physical and biological processes driven by the EAC, particularly in linking circulation to ecosystems, is limited. In this special issue of 16 papers on the EAC, we examine the effects of climatic wind-stress forced ocean dynamics on EAC transport variability and coastal sea level, from ENSO to multi-decadal time scales; eddy formation and structure; fine scale connectivity and larval retention. Comparisons with the poleward-flowing Leeuwin Current on Australia’s west coast show differences in ecosystem productivity that can be attributed to the underlying physics in each region. On average there is double the chlorophyll a concentration on the east coast than the west. In comparison to the Leeuwin, the EAC may have less local retention of larvae and act as a partial barrier to onshore transport, which may also be related to the local spawning and early life history of small pelagic fish on each coast. Inter-annual variations in the EAC transport produce a detectable sea-level signal in Sydney Harbour, which could provide a useful fisheries index as does the Fremantle sea level and Leeuwin Current relationship. The EAC’s eddy structure and formation by the EAC are examined. A particular cold-core eddy is shown to have a “tilt” towards the coast, and that during a rotation the flow of particles may rise up to the euphotic zone and then down beneath. In a warm-core eddy, surface flooding is shown to produce a new shallower surface mixed layer and promote algal growth. An assessment of plankton data from 1938-1942 showed that the local, synoptic conditions had to be incorporated before any comparison with the present. There are useful relationships of water mass characteristics in the Tasman Sea and separation zone with larval fish diversity and abundance, as well as with long-line fisheries. These fisheries-pelagic habitat relationships are invaluable for fisheries management, as well as for climate change assessments.There is further need to examine the EAC influence on rainfall, storm activity, dust deposition, and on the movements by fish, sharks and whales. The Australian Integrated Marine Observing System (IMOS) has provided new infrastructure to determine the changing behaviour of the EAC and its bio-physical interaction with the coasts and estuaries. The forecasting and hindcasting capability developed under the Bluelink project has provided a new tool for data synthesis and dynamical analysis. The impact of a strengthening EAC and how it influences the livelihoods of over half the Australian population, from Brisbane to Sydney, Hobart and Melbourne, is just being realised.     

Genesis and solvus relations of submicroscopically intergrown paragonite and phengite in a blueschist from northern California

Electron microbeam techniques have been used to examine submicroscopically intergrown paragonite, phengite and chlorite from the South Fork Mountain Schist of the Franciscan Terrane of northern California, which was subjected to blueschist facies metamorphism. The sample also contains quartz, albite, lawsonite, and rutile. The subassemblage albite-lawsonite-rutile requires metamorphic conditions on the low-temperature side of the equilibrium albite+lawsonite+rutile=paragonite+sphene+quartz+H₂O (T<200° C and P<7.4 kbars based on thermodynamic data of Holland and Powell 1990). The white micas appear to be optically homogeneous, but back-scattered electron images can distinguish two different micas by their slight difference in contrast. Electron microprobe analyses (EMPA) of micas show Na/(Na+K) ranging from 0.2 to 0.8. The two micas are resolved by transmission electron microscopy (TEM) as packets of phengite and paragonite that range from 20 to several hundred nm in thickness. The compositions, determined by analytical electron microscopy (AEM), constrain the limbs of the phengite-paragonite solvus to values of Na/(Na+K)=<0.02 and 0.97, representing less mutual solid solution than ever reported by EMPA. The textural relations imply that the sheet silicates were derived from reactions between fluids and detrital clays and that they are in an intermediate stage of textural development. We caution that microprobe analyses of apparently homogeneous sheet silicates may yield erroneous data and lead to faulty conclusions using phengite barometry and paragonite-muscovite thermometry, especially in fine-grained rocks that formed at relatively low temperatures. Contribution no. 473 from the Mineralogical Laboratory, Department of Geological Sciences, The University of Michigan, Ann Arbor, Michigan, USA     

A Common Reference Material for Cadmium Isotope Studies – NIST SRM 3108

Research into natural mass‐dependent stable isotope fractionation of cadmium has rapidly expanded in the past few years. Methodologies are diverse with MC‐ICP‐MS favoured by all but one laboratory, which uses thermal ionisation mass spectrometry (TIMS). To quantify the isotope fractionation and correct for instrumental mass bias, double‐spike techniques, sample‐calibrator bracketing or element doping has been used. However, easy comparison between data sets has been hampered by the multitude of in‐house Cd solutions used as zero‐delta reference in different laboratories. The lack of a suitable isotopic reference material for Cd is detrimental for progress in the long term. We have conducted a comprehensive round‐robin assay of NIST SRM 3108 and the Cd isotope offsets to commonly used in‐house reference materials. Here, we advocate NIST SRM 3108 both as an isotope standard and the isotopic reference point for Cd and encourage its use as ‘zero‐delta’ in future studies. The purity of NIST SRM 3108 was evaluated regarding isobaric and polyatomic molecular interferences, and the levels of Zn, Pd and Sn found were not significant. The isotope ratio ¹¹⁴Cd/¹¹⁰Cd for NIST SRM 3108 lies within ～ 10 ppm Da^?1 of best estimates for the Bulk Silicate Earth and is validated for all measurement technologies currently in use.     

Rapid transport of anthropogenic lead through soils in southeast Missouri

To investigate Pb transport and cycling, soils from the forest floor and cores from White Oaks (Quercus alba L.) were collected near a Pb smelter in SE Missouri at varying depths from the surface and varying distances. Lead concentrations in soil samples at the surface drop dramatically with distance from approximately 1500 mg/kg at less than 2 km from the smelter to around 100 mg/kg at localities greater than 2 km from the smelter. Lead contents in tree rings are below 0.5 mg/kg in samples dated prior to 1970, and rapidly increase in 1975–1990 samples. Isotopic compositions of soils and tree rings exhibit systematic variations of Pb isotopic compositions with depth and tree ring age. Distinguishable isotopic signatures for Pb sources allowed quantification of the contribution of smelter Pb to the soils. At depths where Pb concentrations decreased and approached constant values (10–25 cm, 10–30 mg/kg), 50–90%, 40–50% and 10–50% of the Pb could be derived from the smelter for the samples at locations less than 2, 2–4 and over 4 km from the smelter, respectively. The remaining portion was attributable to automobile emission and bedrock sources. Because the smelter operated from 1963 to 2003 and samples were collected in 1999, it is estimated that smelter Pb infiltrates at rates of 1 cm/yr (30 cm in 30 yr). At distances less than 1.5 km from the smelter, even though Pb concentrations become asymptotic at a depth of 30 cm, isotopic evidence suggests that Pb has migrated below this depth, presumably through exchange with naturally occurring Pb in the soil matrix. This implies that soils heavily polluted by Pb can exceed their Pb carrying capacity, which could have potential impacts on shallow groundwater systems and risk further exposure to human and ecological receptors.     

Laboratory Demonstration of Hydrocarbon Migration in the Vadose Zone: Effectiveness of the U-Tube Design for Underground Storage Tank Leak Detection Monitoring

The experimental results derived from a laboratory model of gasoline migration from a simulated leaking underground storage tank (UST) demonstrate that the direction of migration through the simulated fill material is essentially straight down. Consequently, the U-tube design for leak detection monitoring of USTs has limited utility because that design relies on the assumption that the path of migration of the leaked fluids will follow the surface of the tank to the bottom and be intercepted by the U-tube below.