Regulating the ultimate sink: managing the risks of geologic CO2 storage

The geologic storage (GS) of carbon dioxide (CO2) is emerging as an important tool for managing carbon. While this Journal recently published an excellent review of GS technology (Bruant, R. G.; Guswa, A. J.; Celia, M. A.; Peters, C. A. Environ. Sci. Technol. 2002, 36, 240A-245A), few studies have explored the regulatory environment for GS or have compared it with current underground injection experience. We review the risks and regulatory history of deep underground injection on the U.S. mainland and surrounding continental shelf. Our treatment is selective, focusing on the technical and regulatory aspects that are most likely to be important in assessing and managing the risks of GS. We also describe current underground injection activities and explore how these are now regulated.

Geomycology

Fungi have long been known to have capabilities for reduction and alkylation of arsenate and selenate but their general capabilities for solubilizing and accumulating metallic substances have been given serious attention only in recent years. Common members of the Aspergillaceae cultured on boron, copper, lead and other metals or oxides can solubilize and concentrate the elements or their compounds. To account for bio-solubilization of the metals, we have set up a model study, incubating selected metals, e.g. mercury, in solutions of various metabolites including L-lysine and citric acid. Results of 100-300 days incubation showed that many metals can in fact be readily solubilized, and in some cases more effectively at pH 6-7 than at pH 1.5-2.

Phylogeography of the endangered darkling beetle species of Pimelia endemic to Gran Canaria (Canary Islands)

Phylogenetic and geographical nested clade analysis (NCA) methods were applied to mitochondrial DNA sequences of Pimelia darkling beetles (Coleoptera, Tenebrionidae) endemic to Gran Canaria, an island in the Canary archipelago. The three species P. granulicollis, P. estevezi and P. sparsa occur on the island, the latter with three recognized subspecies. Another species, P. fernandezlopezi (endemic to the island of La Gomera) is a close relative of P. granulicollis based on partial Cytochrome Oxidase I mtDNA sequences obtained in a previous study. Some of these beetles are endangered, so phylogeographical structure within species and populations can help to define conservation priorities. A total of about 700 bp of Cytochrome Oxidase II were examined in 18 populations and up to 75 individuals excluding outgroups. Among them, 22 haplotypes were exclusive to P. granulicollis and P. estevezi and 31 were from P. sparsa. Phylogenetic analysis points to the paraphyly of Gran Canarian Pimelia, as the La Gomera P. fernandezlopezi haplotypes are included in them, and reciprocal monophyly of two species groups: one constituted by P. granulicollis, P. estevezi and P. fernandezlopezi (subgenus Aphanaspis), and the other by P. sparsa'sensu lato'. The two species groups show a remarkably high mtDNA divergence. Within P. sparsa, different analyses all reveal a common result, i.e. conflict between current subspecific taxonomic designations and evolutionary units, while P. estevezi and P. fernandezlopezi are very close to P. granulicollis measured at the mtDNA level. Geographical NCA identifies several cases of nonrandom associations between haplotypes and geography that may be caused by allopatric fragmentation of populations with some cases of restriction of gene flow or range expansion. Analyses of molecular variance and geographical NCA allow definition of evolutionary units for conservation purposes in both species-groups and suggest scenarios in which vicariance caused by geological history of the island may have shaped the pattern of the mitochondrial genetic diversity of these beetles.

Isotopes as tracers in a contaminated fractured chalk aquitard

Clusters of industrial plants often generate contaminant plumes with several potential sources. Prevention of further pollution and designing suitable remedial measures require identification of the contributing source among all potential ones and the sorting of currently active sources from historical ones. In the study area, an industrial complex in the Negev desert, Israel, contaminants could not serve as indicators for the contamination sources because of their extensive spatial distribution across the site. However, stable isotopes of oxygen, hydrogen and sulfur, as well as tritium, proved to be efficient tools for this task. The isotopic characterization of the potential end members provided the criteria for constraining a contaminating source when several alternative sources appeared viable. The isotopic fractionation of oxygen and hydrogen isotopes could be tied to the various disposal phases of the industrial wastewater. The three case studies presented here confirm the important role of isotopes as tracers in contaminated sites.

Hydrogeological and hydrochemical investigation for below-sea-level quarrying at cement raw material site (Kocaeli-Darica, Turkey)

A research has been carried out to investigate the effects of below sea level mining on the cement raw material quality of a limestone quarry located adjacent to the shoreline near Darica-Kocaeli-Turkey. Field studies involved rock mass characterization through discontinuity surveys performed at the working benches of the quarry as well as on the core samples, monitoring of groundwater levels, performance of water pressure tests, and in-situ hydrochemical measurements. Hydrogeological data suggest that the carbonate sequence forms a poor unconfined aquifer having hydraulic conductivity values ranging between 10(-6) m/s and 10(-8) m/s. In the quarry, water seepages can only be observed at the shear zones. Electrical conductivity profiles taken from the boreholes located at various distances from the shore line indicated the present position of the salt water wedge.

Modeling of flow and contaminant transport in coupled stream-aquifer systems

This paper describes an integrated approach for modeling flow and contaminant transport in hydraulically connected stream-aquifer systems. The code, FTSTREAM, extended the capabilities of the ground-water model, FTWORK, to incorporate chemical fate and transport in streams. Flow in the stream network is modeled as an unsteady, spatially varying flow, while transport modeling is based on a one-dimensional advection-dispersion equation. In addition to sorption and decay during transport in ground water, the model incorporates volatilization, settling and decay during transport in surface water. The interaction between surface water and ground water is accommodated by a leakage term and is implemented in the model using an iterative Picard-type procedure to ensure mass conservation across the interface between the two systems. The modeling approach is used to simulate contaminant transport in the Mad River, Ohio, which is hydraulically connected to a buried valley aquifer of sand and gravel outwash. The river is a receiving stream in the upstream part of the modeled area. Downstream, heavy pumping from a municipal well field causes the river to become a loosing stream. Induced infiltration from the river is responsible for a considerable portion of the well yield. The flow and transport model, developed for this study, simulates coupling between flow in the aquifer and the river. Hypothetical sources of contamination are introduced at selected locations in the upstream portion of the aquifer. The model is then used to simulate the expected transport in both the aquifer and the stream. A series of simulations elucidates the role of the river in facilitating the transport of the hypothetical contaminants in ground water and surface water. Effect of sorption, retardation and volatilization on contaminant transport is also examined for the case of the volatile organic compounds.

Environmental protection problems in the vicinity of the Zelazny most flotation wastes depository in Poland

The Zelazny Most depository of wastes from copper-ore processing, located in southwest Poland, is the largest mineral wastes repository in Europe. Moreover, it is located in a seismically active area. The seismicity is induced and is connected with mining works in the nearby underground copper mines. Any release of the contents of the repository to the environment could have devastating and even catastrophic consequences. For this reason, geophysical methods are used for continuous monitoring the state of the repository containment dams. The article presents examples of the application of geoelectric methods for detecting sites of leakage of contaminated water and a sketch of the seismic hazard analysis, which was used to predict future seismic vibrations of the repository dams.

Biogeochemistry of mercury in the Amazonian environment

In this paper, the processes that affect mercury (Hg) cycling in the Amazonian environment were reviewed, criticized and new directions of research are proposed. The discussion of the origin of the mercury contamination, whether natural or anthropogenic is marked by a lack of fundamented arguments from both sides. Undoubtedly mercury inputs from gold mining have locally increased environmental concentrations, but in the whole Amazon, these loads would be insignificant, considering the high concentrations observed by some authors in remote soils. A reasonable process that should explain these elevated concentrations in soil is that B horizons function as a mercury "sponge" that have been accumulating mercury over a geological time scale, releasing it back to cycling during erosion and forest fires. The environmental degradation of the Amazonian forest due to human activities is probably enhancing the release of that mercury to the cycle. Mercury transformations in reduced, anoxic environments--sediments and waters--are also a key problem for the understanding of the environmental methylation. The studies that have been carried out in the Amazonian environment are too restricted and results permit only circumstantial conclusions. Large efforts must be directed to monitoring programs considering time and space variability.

Hot rocks or no hot rocks: overnight retreat availability and selection by a diurnal lizard

I used radio telemetry to determine the effects of substrate size and composition on overnight retreat site selection by western fence lizards (Sceloporus occidentalis). In watersheds of northern California (USA), these lizards occupy two habitat types differing in substrate characteristics: rocky cobble bars found in the dry, active channels of rivers and grassy upland meadows. Rocky substrates, found almost exclusively on cobble bars, provided warmer potential retreat sites than all available retreat sites on meadows during the first 5 h of inactivity. Only cobble and sand substrates provided retreats with temperatures within the preferred daily active range (32-36 degrees C) during the inactive period for these lizards (1900-0900 hours). Females on a cobble bar used rocks as retreats on >90% of nights during the breeding season whereas females on a meadow used wood (>70% of nights) and burrows (>25% of nights). In contrast to females, cobble bar males used rocks significantly less frequently (<70%) and slept in the open air significantly more frequently (25% vs. <1%). Cobble bar females further, showed a significant preference for cobbles 15 cm thick, whereas the rocks used by males did not differ significantly in thickness from those measured in randomly placed transects. Rocks 15 cm thick were the warmest retreats commonly available on this habitat type. Thus, thermal microenvironments available to and chosen by gravid female lizards differ considerably between river and non-river habitats.

Biogeochemical assessment of resilient humus formations from virgin and cultivated northern Botswana soils

Humic matter in virgin and cultivated Vertisols and Arenosols from Pandamatenga (northern Botswana) has been studied by chemical fractionation and visible and infrared derivative spectroscopies. The in vitro soil respiration was also determined. Soils contained <15 g of total C kg(-1) and displayed scant mineralization activity. In Vertisols, cultivation has led to scarce significant changes in humus characteristics, pointing to a noteworthy resilience of the organic matter. Humic acids showed a very dark color, indicating a large concentration of aromatic structures and stable free radicals. Infrared spectra were featureless and alike. This strong structural stability is also suggested by Curie-point pyrolysis of humic acids, which failed to yield substantial amounts of diagnostic products, mainly in the case of Vertisols. Some unexpected similarities between Vertisols and Arenosols indicate that the influence of external factors on the humic acid formation processes prevails on that of the geological substrate. In Pandamatenga soils only small amounts of recalcitrant C and N forms are sequestered; they represent stable pools relatively independent from short- or medium-term climatic changes or management practices.

Mercury in sediments of Ulhas estuary

Hg levels in water, suspended particulate matter and sediment of the Ulhas estuary are under considerable environmental stress due to the indiscriminate release of effluents from a variety of industries including chlor-alkali plants. Concentration ranges of dissolved (0.04-0.61 micro gl(-1)) and particulate (1.13-6.43 micro gg(-1)) Hg reveal a definite enhancement of levels in the estuary. The Hg burden in sediment upstream of the weir that limits the tidal influence is low (0.08-0.19 micro gg(-1)) with low C(org) content (1.8-2.9%). The high Hg content of the sediment just below the weir varies seasonally (highest concentration recorded being 38.45 micro gg(-1)) due to incremental accretion of sediment as the fresh water flow over the weir progressively decreases. The 30 km segment of the estuary sustains markedly high levels of Hg in the sediment with an exponential decrease in the seaward direction from the weir. Higher concentrations than the expected background prevail in all the estuarine cores up to the bottom, though the overall concentration decreases from about 20 micro gg(-1) in core 7 (inner estuary) to 1 micro gg(-1) in core 31 (outer estuary). The Hg in sediment is associated with C(org), while its correlation with Al, Fe and Mn is poor. The Hg profiles in cores from the Arabian Sea (stations 34, 35 and 37) have a distinct horizon of enhanced concentration in the 5-60 cm segment. Based on 210Pb dating of core 37, the sediment at the bottom of this core is inferred to have been deposited in the year 1949, roughly two year prior to the establishment of the first chlor-alkali plant and represents the background (0.06-0.10 micro gg(-1)). The Hg profiles in the offshore cores indicate a marked increase in transfer of Hg to sediment subsequent to 1980, with a peak around 1990-1992. Based on the index of geoaccumulation it is considered that the estuarine segment between stations 4 and 23 is extremely polluted, while the sediment from the open coast is moderately polluted in the top 25-30 cm with respect to Hg. The enrichment factor of Hg in the sediment is 350-700 for core 4 and decreases to 0-7 for the open-shore cores.

Characterizing benthic substrates of Santa Monica Bay with seafloor photography and multibeam sonar imagery

Seafloor photography from three cruises is combined with multibeam sonar imagery to characterize benthic substrates and associated fauna of Santa Monica Bay, California. The multibeam EM1000 imagery was collected in 1996. Two sampling cruises (in 1998 and 1999) provided photographs at 142 sites throughout the Bay; a final cruise (in 2000) collected still photographs and continuous video along nine transects on the mainland shelf from Pt. Dume to the Palos Verdes peninsula. Muddy substrates (typically low backscatter) were the predominant habitat throughout the Santa Monica Bay, from the 20 m isobath to the adjacent Santa Monica basin floor (780 m). Bioturbation was pervasive as evidenced by abundant open burrows, mounds, and faunal tracks and trails. Sandy substrates (typically intermediate to high backscatter) were restricted to the innermost mainland shelf and a narrow outer shelf band north of Santa Monica Canyon. Cobble and gravel substrates (high backscatter) were restricted to the innermost shelf south of El Segundo and limited parts of the shelf edge. Rocky substrates (high backscatter) with interspersed patches of sand and gravel occurred on the high-relief marginal plateau and along parts of the shelf break offshore of Malibu.

Vicariance and dispersal across Baja California in disjunct marine fish populations

Population disjunctions, as a first step toward complete allopatry, present an interesting situation to study incipient speciation. The geological formation of the Baja California Peninsula currently divides 19 species of fish into disjunct populations that are found on its Pacific Coast and in the northern part of the Gulf of California (also called the Sea of Cortez), but are absent from the Cape (Cabo San Lucas) region. We studied the genetic makeup of disjunct populations for 12 of these 19 fish species. Phylogeographic patterns for the 12 species can be separated into two major classes: a first group (eight species) showed reciprocal monophyly and high genetic divergence between disjunct populations. A second group (four species) displayed what appeared to be panmictic populations. Population structure between Pacific Coast populations, across the Punta Eugenia biogeographic boundary, was also evaluated. While dispersal potential (inferred by pelagic larval duration) was a poor predictor of population structure between Gulf of California and Pacific populations, we found that population genetic subdivision along the Pacific Coast at Punta Eugenia was always positively correlated with differentiation between Pacific and Gulf of California populations. Vicariant events, ongoing gene flow, and ecological characteristics played essential roles in shaping the population structures observed in this study.

Speciation in North American black basses, Micropterus (Actinopterygii: Centrarchidae)

The Pleistocene Epoch has been frequently cited as a period of intense speciation for a significant portion of temperate continental biotas. To critically assess the role of Pleistocene glaciations on the evolution of the freshwater fish clade Micropterus, we use a phylogenetic analysis of complete gene sequences from two mitochondrial genes (cytochrome b and ND2), and a fossil calibration of the molecular clock to estimate ages of speciation events and rates of diversification. The absence of substantial morphological and ecological divergence together with endemism of five of the eight species in North American tributaries of the Gulf of Mexico may be interpreted as the result of a recent Pleistocene origin for these species. Speciation dates in Micropterus range from 1.01 +/- 0.32 to 11.17 +/- 1.02 million years ago. Only one speciation event is dated to the Pleistocene, and rates of diversification are not significantly variable in Micropterus. The premise that the Pleistocene was an exceptional period of speciation in Micropterus is not supported. Instead, a Gulf Coast allopatric speciation model is proposed, and predicts periods of dynamic speciation driven by sea level fluctuations in the Late Miocene and Pliocene. The Pleistocene, however, was a period of significant intraspecific mitochondrial lineage diversification. The application of the Gulf Coast allopatric speciation model to the remaining aquatic fauna of the Gulf of Mexico coast in North America will rely on robust phylogenetic hypotheses and accurate age estimations of speciation events.

Microbial respiration and diffusive transport of O2, 16O2, and 18O15O in unsaturated soils and geologic sediments

Molecular oxygen (O2) in unsaturated geologic sediments plays an important role in soil respiration, biodegradation of organic contaminants, metal oxidation, and global oxygen and carbon cycling, yet little is known about oxygen isotope fractionation during the consumption and transport of O2 in unsaturated zones. We used a laboratory kinetic cell technique to quantify isotope fractionation due to respiration and a numerical model to quantify both consumptive and diffusive fractionation of O2 isotopes at a field site comprised of unsaturated lacustrine sandy materials. The combined use of laboratory-based kinetic cell experiments and field-based isotope transport modeling provided an effective tool to characterize microbial respiration in unsaturated media. Based on results from the closed-system kinetic cells, O2 consumption and isotope fractionation were attributed to the alternative cyanide-resistant respiration pathway. At the field site, the modeled depth profiles for O2 and delta18O matched the measured in situ data and confirmed that the consumption of O2 was via the alternative respiration pathway. If the cyanide-resistant respiration pathway is indeed widespread in soils, its high oxygen isotope enrichment factor could help to explain the discrepancy between the predicted present-day Dole effect (+20.8/1000) and the observed Dole effect (+23.5/1000). Thus, further soil O2 isotope studies are needed to better characterize and model the fractionation of oxygen isotopes during subsurface respiration and the potential impact on the isotopic content of atmospheric O2.

Main results and experience obtained on Mir space station and experiment program for Russian segment of ISS

This article presents main scientific and practical results obtained in course of scientific and applied research and experiments on Mir space station. Based on Mir experience, processes of research program formation for the Russian Segment of the ISS are briefly described. The major trends of activities planned in the frames of these programs as well as preliminary results of increment research programs implementation in the ISS' first missions are also presented.

Application of harmonic analysis of water levels to determine vertical hydraulic conductivities in clay-rich aquitards

A harmonic analysis method was used to determine vertical hydraulic conductivities (Kv) in geologic media between vertically separated piezometers using water level measurements. In this method, each water level time series was filtered and then decomposed using harmonic analysis into a sum of trigonometric components. The phase and amplitude of each harmonic function were calculated. These data were used to estimate Kv values between vertically separated data sets assuming one-dimensional transient flow. The method was applied to water level data collected from nested piezometers at two thick clay-rich till aquitards in Saskatchewan, Canada. At one site, routine water levels were measured in 12 piezometers (installed between 1 and 29 m below ground surface) since installation (1995). At the other site, water levels were measured in seven piezometers (installed between 4 and 53 m below ground surface) since installation (1998-1999). The Kv calculated using harmonic analysis decreased with depth below the water table at both sites, approaching matrix estimates of hydraulic conductivity between 10 and 11 m and between 21 and 43 m below ground surface. These depths reflected the depth of extensive vertical fracturing at the sites and showed that the depth of fracturing may be site specific.

Effect of randomly methylated beta-cyclodextrin on physical properties of soils

The application of cyclodextrins in several soil remediation technologies has been increasingly studied, but little is known about their effects on soil physical properties. One of the popular soil remediation additives, randomly methylated beta-cyclodextrin (RAMEB), was found to significantly alter surface and pore properties of soil clay minerals. Therefore, in this paper we studied the effect of various RAMEB doses on physical properties of selected soils, representing a wide range of clay content (3-49%). The results showed that soil physical properties were greatly modified by RAMEB treatment. Analysis of water vapor adsorption isotherms revealed that RAMEB increased water adsorption and surface area in sandy soils and decreased them in clayey soils. An increase in adsorption energy of water in RAMEB-treated soils indicated that desorption of nonpolar pollutants can be enhanced. Water vapor desorption isotherms showed that the volumes and radii of micropores (nanometers range) increased above 1% RAMEB concentration. The micropores became more rough and complex after RAMEB treatment as deduced from an increase in fractal dimensions. The volume of soil mesopores measured by mercury intrusion porosimetry (micrometers range) gradually decreased in most soils with an increase in RAMEB concentration whereas the average mesopore radius increased, indicating that finer mesopores were blocked by RAMEB. Measurements of the granulometric composition of soils by sedimentation analysis showed that the amount of coarse-size soil fractions increased on the expense of finer fractions due to aggregation of smaller particles. Behavior of the studied soils after RAMEB treatment depended on their clay content and the dose of cyclodextrin. In clay-rich soils, strong interactions of cyclodextrins with the soil solid phase governed the resulting soil properties. In clay-poor soils, the cyclodextrin excess (not interacted with clays) played a dominant role. Modification of surface, pore, and aggregation properties of soils by RAMEB can have a significant effect in soil remediation technologies.

Measurement and analysis of non-Fickian dispersion in heterogeneous porous media

Contaminant breakthrough behavior in a variety of heterogeneous porous media was measured in laboratory experiments, and evaluated in terms of both the classical advection-dispersion equation (ADE) and the continuous time random walk (CTRW) framework. Heterogeneity can give rise to non-Fickian transport patterns, which are distinguished by "anomalous" early arrival and late time tails in breakthrough curves. Experiments were conducted in two mid-scale laboratory flow cells packed with clean, sieved sand of specified grain sizes. Three sets of experiments were performed, using a "homogeneous" packing, a randomly heterogeneous packing using sand of two grain sizes, and an exponentially correlated structure using sand of three grain sizes. Concentrations of sodium chloride tracer were monitored at the inflow reservoir and measured at the outflow reservoir. Breakthrough curves were then analyzed by comparison to fitted solutions from the ADE and CTRW formulations. In all three systems, including the "homogeneous" one, subtle yet measurable differences between Fickian and non-Fickian transport were observed. Quantitative analysis demonstrated that the CTRW theory characterized the full shape of the breakthrough curves far more effectively than the ADE.

Spatial analysis of grain size in Santa Monica Bay

Maps are useful scientific tools for presenting environmental information, but the statistical techniques necessary to prepare scientifically rigorous maps have primarily focused on terrestrial habitats. This study compares three popular techniques (triangulation, kriging, and co-kriging) to map sediment grain size in Santa Monica Bay, California. Two grain size data sets, one collected in 1994 (79 sites) and one collected in 1997 and 1998 (149 sites) were used for model development. A bathymetric data set collected in 1997 was used as a model covariate. A third grain size data set (40 sites) collected in 1996 from independent sites was used for model evaluation. Predictions were compared to validation data by average difference, prediction mean square error (PMSE), and a goodness-of-prediction measure, G. The average difference between prediction and truth was similar for all methods, but the PMSE for triangulation was more than twice that for kriging or co-kriging, which were similar. The G measure also shows triangulation to be a far worse predictor than kriging and co-kriging. Small-scale differences were observed between kriging and co-kriging at steep depth contours, where co-kriging predicted values commensurate with the expected depth-defined grain size.

Magnitude and variability of Holocene sediment accumulation in Santa Monica Bay, California

The spatial variability of Holocene (past 10,000 years) sediment accumulation in Santa Monica Bay (California) was examined to identify controls sediment trapping in a bathymetrically complex coastal embayment and to provide geologic context for the post-industrial sedimentary record and associated pollution gradients. Sediment chronologies based on downcore AMS 14C dates were used to quantify long-term (millennia) accumulation rates in an effort to elucidate particle-transport pathways and sinks. Sediment accumulation rates for the full range of bayfloor environments (50-630 m water depths) range from 22 to 102 mg/cm(2)/year (15-88 mm/100 year), have an overall mean of 51+/-21 mg/cm(2)/year (1 sigma, n=11), and are comparable to rates reported for adjacent borderland basins. Maximal accumulation rates on the Malibu shelf and within a reentrant to Redondo canyon are interpreted to reflect (1) proximity to sediment sources and (2) localized oceanographic and topographic conditions conducive to sediment trapping and deposition. The 14C-derived accumulation rates are 2-10 times lower than rates determined through (210)Pb geochronology for the same sites in a related study, revealing that Holocene sediment accumulation has been non-steady-state. Santa Monica Bay is an important sink for suspended matter; averaged over the past several millennia a mass of sediment equivalent to 10-80% of the modern annual river supply is sequestered yearly. Net influx of suspended matter derived from the adjacent Palos Verdes shelf is evinced by a concentration gradient of p,p'-DDE in bayfloor sediments, whereas the distribution of anthropogenic silver suggests transport from Santa Monica shelf to the southeastern boundary of the bay. The results of this study provide new insight to the long-term fates of particulate matter in Los Angeles coastal waters.

Simulating bioremediation of uranium-contaminated aquifers; uncertainty assessment of model parameters

Bioremediation of trace metals and radionuclides in groundwater may require the manipulation of redox conditions via the injection of a carbon source. For example, after nitrate has been reduced, soluble U(VI) can be reduced simultaneously with other electron acceptors such as Fe(III) or sulfate to U(IV), which may precipitate as a solid (uraninite). To simulate the numerous biogeochemical processes that will occur during the bioremediation of trace-metal-contaminated aquifers, a time-dependent one-dimensional reactive transport model has been developed. The model consists of a set of coupled mass balance equations, accounting for advection, hydrodynamic dispersion, and a kinetic formulation of the biological or chemical transformations affecting an organic substrate, electron acceptors, corresponding reduced species, and trace metal contaminants of interest, uranium in this study. This set of equations is solved numerically, using a finite difference approximation. The redox conditions of the domain are characterized by estimating the pE, based on the concentration of the dominant terminal electron acceptor and its corresponding reduced species. This pE and the concentrations of relevant species are then used by a modified version of MINTEQA2, which calculates the speciation/sorption and precipitation/dissolution of the species of interest under equilibrium conditions. Kinetics of precipitation/dissolution processes are described as being proportional to the difference between the actual and calculated equilibrium concentration. A global uncertainty assessment, determined by Random Sampling High Dimensional Model Representation (RS-HDMR), was performed to attain a phenomenological understanding of the origins of output variability and to suggest input parameter refinements as well as to provide guidance for field experiments to improve the quality of the model predictions. By decomposing the model output variance into its different input contributions, RS-HDMR can identify the model inputs with the most influence on various model outputs, as well as their behavior pattern on the model output. Simulations are performed to illustrate the effect of biostimulation on the fate of uranium in a saturated aquifer, and to identify the key processes that need to be characterized with the highest accuracy prior to designing a uranium bioremediation scheme.