Contaminant Spread and Flushing in Fractured Rocks Near Oak Ridge, Tennessee

Liquid wastes, including metals dissolved in nitric acid, were discharged into the S-3 Ponds from 1951 to 1983. During this period, contaminants in ground water spread along shallow fracture flow paths toward nearby streams. Also, a high concentration of nitrate in one well at a depth of 110 to 240 in shows that some contaminants may have moved downdip because of differences in fluid density. Neutralization of the ponds in June 1983 caused a dramatic decrease in the contaminant concentrations of Bear Creek about three months later. Since then, the contaminant concentrations of Bear Creek have decreased at a first-order exponential rate. This average rate, which is the same for both more reactive and less reactive constituents, can be interpreted to show that the contaminant reservoir consists of the unfractured rock matrix.
Flushing caused by the natural recharge and discharge of ground water is occurring at all locations, but contaminant concentrations are controlled by the relative rates of molecular diffusion from the rock matrix and advection along the fracture flow paths. Hushing has thus been most effective near the water table. If the exponential decrease in contaminant concentrations continues, water in Bear Creek will meet drinking water standards by 2012: regardless of any remedial action, contaminants will remain in the rocks for many years.     

Absolute Calibration of the TOPEX/POSEIDON Altimeters using UK Tide Gauges, GPS, and Precise, Local Geoid-Differences

This article describes an "absolute" calibration of TOPEX/POSEIDON (T/P) altimeter bias using UK tide gauges equipped with Global Positioning System (GPS) receivers. The method is an extension of earlier work using the Newhaven tide gauge in the English Channel. However, the present analysis extends the research to a number of gauges around the UK and incorporates several improvements. The time-averaged TOPEX and POSEIDON biases are obtained with a precision of approximately 2 and 3 cm, respectively. The research complements work on bias determination by other groups in the T/P Science Working Team and can, in principle, be applied at other locations for which precise, local geoid-differences are available. The relatively sparse POSEIDON data set has been used as a test of our ability to perform an absolute calibration of upcoming Jason-1 altimetry as soon as possible after launch.     

Significance of metallothioneins and lysosomes in cadmium toxicity and homeostasis in the digestive gland cells of mussels exposed to the metal in presence or absence of phenanthrene

Metallothioneins and lysosomes are known to be involved in cellular detoxication and sequestration of certain metals^1–3 and both have been identified in this role in elimination of copper from marine mussels (Mytilus edulis/galloprovincialis).³ Cadmium (Cd), however, has been shown to persist in the cells of the digestive gland for long periods with only minimal elimination. An experiment was designed to test the effects of Cd on the fragility of lysosomal membranes in the digestive cells as a measure of cellular injury,^4,5 metallothionein content of the digestive gland and cadmium concentration in this organ. Phenanthrene was used also to destabilise lysosomal membranes⁶ in order to test if increased lysosomal fragility interfered with cadmium metabolism and detoxication. The results demonstrated that Cd induced metallothionein synthesis and that elimination of Cd was minimal after 28 days in clear seawater. Lysosomal fragility was initially increased but this effect was soon reversed, even with continued exposure to Cd. The lysosomal destabiliser, phenanthrene, did not appear to affect accumulation of Cd or levels of metallothionein.     

Ra and Ra in the mixing zones of the Pee Dee River-Winyah Bay, Yangtze River and Delaware Bay Estuaries

²²⁶Ra and ²²⁸Ra have non-conservative excess concentrations in the mixing zones of the Pee Dee River-Winyah Bay estuary, the Yangtze River estuary, and the Delaware Bay estuary. Laboratory experiments, using Pee Dee River sediment, indicate desorption of ²²⁶Ra to increase with increasing salinities up to 20‰. In Winyah Bay desorption from river-borne sediments could contribute almost all of the increases for both isotopes. Desorption adds only a portion of the excess ²²⁸Ra measured in the Yangtse River and adjacent Shelf waters and Delaware Bay. In the Yangtze River the mixing zone extends over a considerable portion of the Continental Shelf where ²²⁸Ra is added to the water column by diffusion from bottom sediments, while ²²⁶Ra concentrations decrease from dilution. Diffusion of ²²⁸Ra from bottom sediments in Delaware Bay primarily occurs in the upper part of the bay (< 22‰ water) where fine grained sediments predominate. A diffusive flux for ²²⁸Ra of 0·33 dpm cm⁻² year was determined for Delaware Bay.     

Submarine groundwater discharge: A large, previously unrecognized source of dissolved iron to the South Atlantic Ocean

This paper reports the initial results of a study of groundwater and coastal waters of southern Brazil adjacent to a 240 km barrier spit separating the Patos Lagoon, the largest coastal lagoon in South America, from the South Atlantic Ocean. The objective of this research is to assess the chemical alteration of freshwater and freshwater–seawater mixtures advecting through coastal permeable sands, and the influence of the submarine discharge of these fluids (SGD) on the chemistry of coastal waters. Here we focus on dissolved iron in this system and use radium isotopic tracers to quantify SGD and cross-shelf fluxes. Iron concentrations in groundwaters vary between 0.6 and 180 μM. The influence of the submarine discharge of these fluids into the surf zone produces dissolved Fe concentrations as high as several micromolar in coastal surface waters. The offshore gradient of dissolved Fe, coupled with results for Ra isotopes, is used to quantify the SGD flux of dissolved Fe from this coastline. We estimate the SGD flux to be 2 × 10⁶ mol day^− 1 and the cross-shelf flux to be 3.2 × 10⁵ mol day^− 1. This latter flux is equal to about 10% of the soluble atmospheric Fe flux to the entire South Atlantic Ocean. We speculate on the importance of this previously unrecognized iron input to regional ocean production and on the potential significance of this source to understanding variations in glacial–interglacial ocean production.     

Assessing the potential for CO2 adsorption in a subbituminous coal,Huntly Coalfield,New Zealand,using small angle scattering techniques

Small angle scattering techniques (SAXS and SANS) have been used to investigate the microstructural properties of the subbituminous coals (R_max 0.42–0.45%) from the Huntly Coalfield, New Zealand. Samples were collected from the two thick (> 5 m) coal seams in the coalfield and have been analysed for methane and carbon dioxide sorption capacity, petrography, pore size distribution, specific surface area and porosity.Specific surface area (SSA) available for carbon dioxide adsorption, extrapolated to a probe size of 4 Å, ranged from 1.25 × 10⁶ cm^? 1 to 4.26 × 10⁶ cm^? 1 with total porosity varying from 16% to 25%. Porosity was found to be predominantly composed of microporosity, which contributed the majority of the available SSA. Although considerable variation was seen between samples, the results fit well with published rank trends.Gas holding capacity at the reservoir pressure (approximately 4 MPa) ranged from 2.63 to 4.18 m³/t for methane on a dry, ash-free basis (daf) and from 22.00 to 23.72 m³/t daf for carbon dioxide. The resulting ratio of CO₂:CH₄ ranged from 5.7 to 8.6, with an average of 6.7:1.Holding capacities for both methane and carbon dioxide on a dry ash free basis (daf) were found to be correlated with sample microporosity. However, holding capacities for the two gases on an as analysed (aa) basis (that is including mineral matter and moisture), showed no such correlation. Carbon dioxide (aa) does show a negative correlation with both specific surface area and microporosity. As the coals have low inorganic matter content, the reversal is thought to be related to moisture which is likely concentrated in the pore size range 12.5–125 Å. Methane holding capacity, both daf and aa, correlates with macroporosity, thus suggesting that the holding capacity of micropores is diminished by the presence of moisture in the pores.     

Uncertainty of gas saturation estimates in a subbituminous coal seam

To assess the commercial viability of a coalbed methane prospect two of the key geological parameters measured are gas content (desorbed gas) and gas holding capacity (adsorption capacity). These two measures, together with reservoir pressure, give an estimate of the gas saturation of the reservoir. Typically gas saturation has been assessed by collecting one adsorption isotherm sample and assuming it is representative of the whole seam reservoir conditions. This study addresses that assumption.To understand the level of variation, and thus the inherent uncertainty in saturation, one core (Jasper-1) from the Huntly coalfield in New Zealand was analysed in detail. Ten samples (representing the whole coal seam) were placed into gas desorption canisters and desorbed for ten days and then analysed for adsorption capacity. Desorption analyses for total measured gas content (average in-situ basis) ranged from 2.32 to 2.89 m³/t (standard deviation (sd) = 0.18) and gas adsorptive capacity at 4 MPa (average in-situ basis) from 2.11 to 3.51 m³/t (sd = 0.38) resulting in saturations ranging from 66% to 120% (sd = 15).Determination of how many samples are required to make a realistic assessment of average reservoir properties requires a consideration of: (i) the level of accuracy desired, (ii) the limit of accuracy possible, which is governed by the magnitude of experimental error, and (iii) the innate variability of the seam. It was found that a minimum of five samples each for adsorption and desorption were required in order to significantly decrease the uncertainty in gas saturation estimates for a subbituminous coal.     

Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA

Karst aquifer studies often focus on allogenic water inputs and large conduit flow. However, diffuse recharge can be significant, particularly in unconfined eogenetic karst aquifers that retain high matrix permeability. This study examines an unconfined region of the upper Floridan aquifer (USA) that hosts a sinking stream, its resurgence, and a large conduit system. Daily diffuse recharge was approximated using a water-budget method and ranged from 17% of precipitation during a low precipitation year to >53% during the highest precipitation year, illustrating the highly variable nature of diffuse recharge in this region. The total allogenic input via the sinking stream over the 5 years of the study was significantly larger than the volume of diffuse recharge. However, only about 2% of the allogenic recharge flows from the conduit into the surrounding aquifer. That flow is restricted to storm events when hydraulic heads in the conduits exceed those in the surrounding aquifer. The estimated volume of dissolution is similar for allogenic recharge and diffuse recharge to the unconfined region surrounding the conduits, but dissolution from the diffuse recharge is distributed over a larger area than dissolution from allogenic recharge. These results exemplify how recharge type impacts flow and water–rock interactions in eogenetic karst aquifers.     

Toward petascale earthquake simulations

Earthquakes are among the most complex terrestrial phenomena, and modeling of earthquake dynamics is one of the most challenging computational problems in science. Computational capabilities have advanced to a state where we can perform wavefield simulations for realistic three-dimensional earth models, and gain more insights into the earthquakes that threaten California and many areas of the world. The Southern California Earthquake Center initiated a major earthquake research program called TeraShake to perform physics-based numerical simulations of earthquake processes for large geographical regions, at high resolution, and for high frequencies. For a large scale simulation such as TeraShake, optimization problems tend to emerge that are not significant in smaller scale simulations. This involves both large parallel computation and also massive data management and visualization coordination. In this paper, we describe how we performed single-processor optimization of the TeraShake AWM application, optimization of the I/O handling, and optimization of initialization. We also look at the challenges presented by run-time data archive management and visualization. The improvements made to the TeraShake AWM code enabled execution on the 40k IBM Blue Gene processors and have created a community code that can be used by seismologists to perform petascale earthquake simulations.