Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

Coordinate mapping of analytical contaminant transport solutions to non-uniform flow fields

Existing analytical solutions to 2D and 3D contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. An approximate method is developed herein for coordinate mapping of 2D (vertically-averaged) transport solutions to non-uniform steady-state irrotational and divergence-free flow fields in single-layer aquifers. The method enables existing analytical transport solutions to be applied to aquifer systems with wells, non-uniform saturated thickness, surface water features, and (to a limited degree) heterogeneous hydraulic conductivity and recharge. This mass-conservative coordinate mapping approach is inexact in its approximation of the dispersion process but is still sufficiently accurate for many simple flow systems. The degree of model error is directly proportional to the variation of velocity magnitude within the domain. These mapped analytical solutions are compared to numerical simulation results and the coordinate mapping errors are investigated. The methods described herein may be used in the traditional capacity of analytical transport models, i.e., screening and preliminary site assessment, without sacrificing accuracy by assuming locally uniform flow conditions or applying an ad-hoc coordinate transformation. The solutions benefit from the traditional advantages of analytical methods, particularly the removal of artifacts due to spatial and temporal discretization: no time-stepping or numerical discretization is required.     

Triggered earthquakes and deep well activities

Earthquakes can be triggered by any significant perturbation of the hydrologic regime. In areas where potentially active faults are already close to failure, the increased pore pressure resulting from fluid injection, or, alternatively, the massive extraction of fluid or gas, can induce sufficient stress and/or strain changes that, with time, can lead to sudden catastrophic failure in a major earthquake. Injection-induced earthquakes typically result from the reduction in frictional strength along preexisting, nearby faults caused by the increased formation fluid pressure. Earthquakes associated with production appear to respond to more complex mechanisms of subsidence, crustal unloading, and poroelastic changes in response to applied strains induced by the massive withdrawal of subsurface material. As each of these different types of triggered events can occur up to several years after well activities have begun (or even several years after all well activities have stopped), this suggests that the actual triggering process may be a very complex combination of effects, particularly if both fluid extraction and injection have taken place locally. To date, more than thirty cases of earthquakes triggered by well activities can be documented throughout the United States and Canada. Based on these case histories, it is evident that, owing to preexisting stress conditions in the upper crust, certain areas tend to have higher probabilities of exhibiting such induced seismicity.     

Selenium and mercury concentrations in harbor seals (Phoca vitulina) from central California: Health implications in an urbanized estuary

We measured total selenium and total mercury concentrations ([TSe] and [THg]) in hair (n = 138) and blood (n = 73) of harbor seals (Phoca vitulina) from California to assess variation by geography and sex, and inferred feeding relationships based on carbon, nitrogen, and sulfur stable isotopes. Harbor seals from Hg-contaminated sites had significantly greater [THg], and lesser [TSe] and TSe:THg molar ratios than seals from a relatively uncontaminated site. Males had significantly greater [THg] than females at all locations. Sulfur stable isotope values explained approximately 25% of the variability in [THg], indicating increased Hg exposure for seals with a greater use of estuarine prey species. Decreased [TSe] in harbor seals from Hg-contaminated regions may indicate a relative Se deficiency to mitigate the toxic effects of Hg. Further investigation into the Se status and the potential negative impact of Hg on harbor seals from Hg-contaminated sites is warranted.     

Geophysical Expression of Low Sulphidation Epithermal Au-Ag Deposits and Exploration Implications –Examples from the Hokusatsu Region of SW Kyushu, Japan–

Abstract: Low sulphidation epithermal Au-Ag deposits constitute a major exploration target for many mining companies worldwide. These deposits are commonly characterized by quartz vein systems which formed at temperatures often in the range 150–250°C giving rise to major alteration halos. Accompanying mineralization and associated alteration is the destruction and redistribution of magnetite in the host rocks and the formation of various K-bearing minerals such as adularia, illite and alunite. These alteration halos can often cover an area of at least several square kilometres and as such have the potential to be detected during reconnaissance geological surveys utilizing high resolution airborne aeromagnetic and radiometric sensors.
Analysis of airborne magnetic and radiometric data recently acquired over the Hokusatsu Region of SW Kyushu indicates that major low sulphidation epithermal Au-Ag deposits are characterized by several distinct geophysical features varying with the degree of post-mineralization cover and level of erosion. Several deposits coincide with extensive magnetic "flat" and high radiometric K–counts. In addition, several workers have shown that a very close association exists between Bouguer gravity highs and epithermal Au mineralization in the Hokusatsu Region. The combined use of areomagnetics, radiometrics and gravity data therefore provides an exceptionally effective approach to exploration for low sulphidation epithermal Au–Ag deposits in the Hokusatsu Region of Japan.     

Evaluating legacy contaminants and emerging chemicals in marine environments using adverse outcome pathways and biological effects-directed analysis

Natural and synthetic chemicals are essential to our daily lives, food supplies, health care, industries and safe sanitation. At the same time protecting marine ecosystems and seafood resources from the adverse effects of chemical contaminants remains an important issue. Since the 1970s, monitoring of persistent, bioaccumulative and toxic (PBT) chemicals using analytical chemistry has provided important spatial and temporal trend data in three important contexts; relating to human health protection from seafood contamination, addressing threats to marine top predators and finally providing essential evidence to better protect the biodiversity of commercial and non-commercial marine species. A number of regional conventions have led to controls on certain PBT chemicals over several years (termed ‘legacy contaminants’; e.g. cadmium, lindane, polycyclic aromatic hydrocarbons [PAHs] and polychlorinated biphenyls [PCBs]). Analytical chemistry plays a key role in evaluating to what extent such regulatory steps have been effective in leading to reduced emissions of these legacy contaminants into marine environments. In parallel, the application of biomarkers (e.g. DNA adducts, CYP1A-EROD, vitellogenin) and bioassays integrated with analytical chemistry has strengthened the evidence base to support an ecosystem approach to manage marine pollution problems. In recent years, however, the increased sensitivity of analytical chemistry, toxicity alerts and wider environmental awareness has led to a focus on emerging chemical contaminants (defined as chemicals that have been detected in the environment, but which are currently not included in regulatory monitoring programmes and whose fate and biological impacts are poorly understood). It is also known that natural chemicals (e.g. algal biotoxins) may also pose a threat to marine species and seafood quality. Hence complex mixtures of legacy contaminants, emerging chemicals and natural biotoxins in marine ecosystems represent important scientific, economic and health challenges. In order to meet these challenges and pursue cost-effective scientific approaches that can provide evidence necessary to support policy needs (e.g. the European Marine Strategy Framework Directive), it is widely recognised that there is a need to (i) provide marine exposure assessments for priority contaminants using a range of validated models, passive samplers and biomarkers; (ii) integrate chemical monitoring data with biological effects data across spatial and temporal scales (including quality controls); and (iii) strengthen the evidence base to understand the relationship between exposure to complex chemical mixtures, biological and ecological impacts through integrated approaches and molecular data (e.g. genomics, proteomics and metabolomics). Additionally, we support the widely held view that (iv) that rather than increasing the analytical chemistry monitoring of large number of emerging contaminants, it will be important to target analytical chemistry towards key groups of chemicals of concern using effects-directed analysis. It is also important to evaluate to what extent existing biomarkers and bioassays can address various classes of emerging chemicals using the adverse outcome pathway (AOP) approach now being developed by the Organization for Economic Cooperation and Development (OECD) with respect to human toxicology and ecotoxicology.     

Transport of ionizing radiation in terrestrial-like exoplanet atmospheres

The propagation of ionizing radiation through model atmospheres of terrestrial-like exoplanets is studied for a large range of column densities and incident photon energies using a Monte Carlo code we have developed to treat Compton scattering and photoabsorption. Incident spectra from parent star flares, supernovae, and gamma-ray bursts are modeled and compared to energetic particles in importance. Large irradiation events with fluences of 10⁶-10⁹ erg cm⁻² at the conventional habitable zone can occur at a rate from many per day (flares from young low-mass parent stars) to ∼100 per Gyr (supernovae and gamma-ray bursts). We find that terrestrial-like exoplanets with atmospheres thinner than about 100 g cm⁻² block nearly all X-rays, but transmit and reprocess a significant fraction of incident γ-rays, producing a characteristic, flat surficial spectrum. Thick atmospheres (?100 g cm⁻²) efficiently block even γ-rays, but nearly all the incident energy is redistributed into diffuse UV and visible aurora-like emission, increasing the effective atmospheric transmission by many orders of magnitude. Depending on the presence of molecular UV absorbers and atmospheric thickness, up to 10% of the incident energy can reach the surface as UV reemission. For the Earth, between 2×10⁻³ and 4×10⁻² of the incident flux reaches the ground in the biologically effective 200-320 nm range, depending on O₂/O₃ shielding. For atmospheres thicker than ∼50 g cm⁻² in the case of pure Rayleigh scattering and ∼100 g cm⁻² in the case of O₂/O₃ absorption, the UV reemission exceeds the surficial transmitted ionizing radiation. We also discuss the effects of angle of incidence and derive a modified two-stream approximation solution for the UV transfer. Finally, we suggest that transient atmospheric ionization layers can be frequently created at altitudes lower than the equilibrium layers that result from steady irradiation and winds from the parent star. We suggest that these events can produce frequent fluctuations in atmospheric ionization levels and surficial UV fluxes on terrestrial-like planets.     

Assessment of channel dynamics,in-stream structures and post-project channel adjustments in North Carolina and its implications to effective stream restoration

Site assessment and monitoring data were analyzed for 26 stream restoration projects in North Carolina where the channel was reconfigured. Post-project changes in channel capacity were highly variable from site to site, but more than 60% of the projects underwent, on average at a given site, at least a 20% change in channel capacity. An analysis of site and basin geomorphology revealed that large post-construction adjustments were associated with highly dynamic stream channels characterized by a combination of high sediment transport capacity, large sediment supply, and/or easily eroded bank materials. In-stream structures along dynamic, reconfigured channels also exhibited high incidences of damage. Thus, the design and construction of channels in a state of equilibrium, which do not exhibit excessive erosion or deposition along highly dynamic rivers is currently problematic. In light of these findings, a conceptual framework based on geomorphic parameters is put forth to assess the likelihood of project success early in the design process to (1) eliminate high risk sites from consideration of channel reconfiguration and (2) improve upon the implemented management strategies that are ultimately used. It is also argued that where space permits an enhanced natural channel, adjustment approach is likely to be more effective than projects based on natural channel design.     

GEOSECS AtlanticSi profiles

Measurements of five cosmogenic³²Si vertical profiles in Atlantic waters (27°N to 60°S) are presented. The amounts of dissolved SiO₂ extracted range from 2 to 54 g; the amounts of water from which SiO₂ was extracted range between 540 kg and 270, 000 kg. In additon, SiO₂ recovered from four surface particulate composites (64°N to 61°S) were also analyzed for³²Si.³²Si measurements were made by milking and counting the daughter activity, ³²P. The net³²P activities range from 0.7 to 6.8 cph; typical errors in measurements of the³²P activities are 20–30%.The³²Si concentrations vary from 0.6 dpm/10⁶ kg of water in the North Atlantic surface waters to 235 dpm/10⁶ kg at 400 m depth in the circumpolar waters. The vertical profiles of³²Si at the five Atlantic stations approximately follow the Si profiles but the depth gradients are different. This would be expected also considering the in-situ release mechanisms due to dissolution and advection/diffusion from the bottom waters. Except for the circumpolar station 89, where the Si and³²Si profiles show the effect of marked vertical mixing (nearly depth independent profiles), the profiles show the following features: (1) specific activities of³²Si (³²Si/SiO₂ ratios) are lowest at intermediate depths, and (2) on an average the surface specific activities are higher, by 2–4 times, than the bottom water values. These data are consistent with generation of the highest specific activity³²Si waters at the surface, where Si concentrations are lowest and precipitation adds cosmogenic³²Si scavenged from the troposphere. Rapid removal of biogenic silica to the water-sediment interface, without much dissolution during transit, leads to the second regime of high³²Si specific activities.The³²Si inventories in the water column in the latitude belt 27°N-27°S are in the range (1–1.4) × 10⁻² dpm³²Si/cm², which is consistent with the expected fallout of cosmogenic³²Si. However, the³²Si column inventories south of 40°S are higher by a factor of 5–7, whereas the corresponding Si inventories increase by only a factor of 3. This excess³²Si in the Southern Ocean cannot be explained by direct fallout from the stratosphere or by melting of Antarctic snow and ice. Instead, this excess is maintained primarily by the southward deep-water transport of³²Si dissolved from sinking particulates.     

Rank of coal beds of the Narragansett basin,Massachusetts and Rhode Island

Coal of the Narragansett basin generally has been considered to be anthracite and/or meta-anthracite. However, no single reliable method has been used to distinguish these two ranks in this basin. Three methods — chemical, X-ray, and petrographic — have been used with some degree of success on coal of the Narragansett basin, but too often the results are in conflict. Chemical methods have been limited by inadequate sampling on a coal-bed-by-coal-bed basis and by a lack of analyses made according to (American Society for Testing and Materials, 1974) standard specifications.In addition, when corrections are made by using the Parr formulas, as required by the ASTM (1974) procedures, the generally high to very high ash content of coal from the Narragansett basin causes the fixed-carbon content to appear higher than it actually is. X-ray methods using the degree of graphitization as a measure of rank are not reliable because some of the graphite is related to shearing and brecciation associated with folding and faulting. Petrographic methods using reflectance on vitrinite give results that are generally consistent with results from chemical determinations. However, it is not clear whether the mean maximum reflectance or mean bireflectance is a better indicator of similar rank of such high-rank coals that have been structurally deformed.Coal from the Cranston Mine, RI, is probably meta-anthracite and coal from the Portsmouth Mine is probably anthracite. These ranks are based on chemical,X-ray, and petrographic data and are supported by associated metamorphic mineral assemblages that indicate that the Cranston Mine is in a higher metamorphic zone than the zone containing the Porthmouth Mine. Interpretation of the rank of Mansfield, MA, coal on the basis of extant chemical data is difficult because it is an impure coal with an ash content of 33 to 50%. Reflectance data indicate that the Mansfield, Foxborough, and Plainville coals in the northern part of the Narragansett basin are meta-anthracite but this is in disagreement with the rank suggested by the low degree of metamorphism of the associated rocks.