Lithophysal Cavities Research Articles

Word to the Wise: Spherulites and Lithophysae

SPHERULITIC OBSIDIAN AND LITHOPHYSAL CAVITIES that host secondary minerals can be prized finds in the collecting community. Both occur in glassy rocks associated with obsidian flows and domes, in r...

Relations between thermal history and secondary structures of ignimbrites exclusive of rheomorphism

Newly calculated model density profiles of ignimbrites also provide model thermal histories, which serve as a framework to help analyze the development of cooling joints, devitrification textures, and lithophysal cavities. Only those parts of sections of Rattlesnake Tuff in central Oregon that resided >600 °C for at least two years show devitrification. Low tensile strength means that joints form by thermoelastic contraction after cooling by as little as 25 °C. This means that columnar joints formed in as little time as a few weeks after deposition of the Aravaipa Tuff in SE Arizona. Compaction is rapid as well, so both columnar jointing and compaction are complete before the onset of devitrification in deposits <40 m thick. A section of Rattlesnake Tuff shows stratabound occurrences of devitrified spots and cavities; devitrification appears to have begun at scattered spots. Most spots are bounded by crescentic cavities in formerly ductile shards. The equation for conductive cooling of a spherical heat source shows that for rock volumes larger than ∼15 cm diameter, the rate of cooling is insufficient to prevent heating by latent heat by as much as 12 °C. Inflation therefore appears to have been caused by vapor released by devitrification and its slight adiabatic expansion. Eventual wholesale devitrification of matrix resulted in more-devitrified spots scattered in less devitrified matrix. Lithophysal cavities in the Rattlesnake Tuff and in the Peach Springs Tuff in NW Arizona are more abundant in lower density horizons, showing that cavity growth is favored in permeable zones between impermeable horizons. The Peach Springs Tuff and other deposits described in the literature have horizontal joints that formed during cooling and whose origin has yet to be deterministically modeled. The longest joints are most closely spaced in the zone of cavities and formed after columnar joints. It is inferred that after formation, each vertical column responds independently to evolving stresses; at this time, asperities together with gas pressures accompanying cavity formation result in horizontal joints. Columnar joints do not completely relieve growing gas pressures during devitrification, likely because of sealing by wholesale inflation of the rock mass and by mineral deposition. Some devitrified horizons of Rattlesnake Tuff and Peach Springs Tuff are densely fractured, resulting in a rubbly surface. These small fractures show zigzag traces, bifurcations, abrupt terminations, and pinch-and-swell walls, similar to ductile fractures described by [Eichhubl (2004)][1]. Their origin is interpreted to be the result of tensile stress due to porosity increase during later stages of devitrification. Closely adjacent sections of Rattlesnake Tuff at one site differ by 18% in thickness, reflecting buried paleotopography. Each section is devitrified. The thicker profile has a zone of lithophysal cavities near its base; the thinner profile has few cavities. The greater thickness translates to only 0.18 MPa additional lithostatic pressure at the base of the thicker profile, which serves to emphasize the subtle differences in initial conditions that can lead to formation of cavities. [1]: #ref-17

Limited hydrologic response to Pleistocene climate change in deep vadose zones — Yucca Mountain, Nevada

Understanding the movement of water through thick vadose zones, especially on time scales encompassing long-term climate change, is increasingly important as societies utilize semi-arid environments for both water resources and sites viewed as favorable for long-term disposal or storage of hazardous waste. Hydrologic responses to Pleistocene climate change within a deep vadose zone in the eastern Mojave Desert at Yucca Mountain, Nevada, were evaluated by uranium-series dating of finely layered hyalitic opal using secondary ion mass spectrometry. Opal is present within cm-thick secondary hydrogenic mineral crusts coating floors of lithophysal cavities in fractured volcanic rocks at depths of 200 to 300 m below land surface. Uranium concentrations in opal fluctuate systematically between 5 and 550 μg/g. Age-calibrated profiles of uranium concentration correlate with regional climate records over the last 300,000 years and produce time-series spectral peaks that have distinct periodicities of 100- and 41-ka, consistent with planetary orbital parameters. These results indicate that the chemical compositions of percolating solutions varied in response to near-surface, climate-driven processes. However, slow (micrometers per thousand years), relatively uniform growth rates of secondary opal and calcite deposition spanning several glacial–interglacial climate cycles imply that water fluxes in the deep vadose zone remained low and generally buffered from the large fluctuations in available surface moisture during different climates.

Methodology for bounding calculations of nuclear criticality of fissile material accumulations external to a waste container at Yucca Mountain, Nevada

The possibility of nuclear criticality, however remote, in the vicinity of the proposed repository at Yucca Mountain, Nevada generates justified concerns and may impact the performance of the repository. A heuristic approach is presented here for determining the amount, spatial distribution and other characteristics of fissile material accumulation in the rock beneath a waste package that could contribute to such an event. This study is concerned primarily with waste packages containing special spent fuel from the Department of Energy and high-level nuclear waste glass. Mixing with less alkaline waters and the subsequent drop in pH is the mechanism that is most efficient for precipitating fissile material from the waste package internal leachate, in contrast to natural deposits in which redox changes are the main precipitation driver. External accumulation size is determined by (1) computing the chemical composition of the leachate leaving a package as its internal materials degrade (with the batch geochemical code EQ3/6), (2) determining precipitation of fissile material into mineral phases (using the 1D geochemical code PHREEQC) as the effluent mixes with percolation water, and (3) heuristically scaling results to a 3D volume and computing the criticality coefficient (using the code MCNP). Loci for accumulation are the multiple lithophysal cavities and the fracture system. A bounding conservative approach is used by necessity in Step 3. Nuclear criticality is sensitive to small variations in the distribution of fissile material and parameters of natural systems vary by orders of magnitude. Because the most likely combinations of parameters are not conducive to nuclear criticality, this study focuses on extreme values of parameter probabilistic distributions, such as limited flow into the package associated with a large percolation rate, combinations of material degradation rates favoring actinide release, and very high host-rock porosity values. By considering these combinations, most favorable to criticality but unlikely, it was concluded that external nuclear criticality is not a concern at the proposed repository.

Infiltration into fractured bedrock

One potential consequence of global climate change and rapid changes in land use is an increased risk of flooding. Proper understanding of floodwater infiltration thus becomes a crucial component of our preparedness to meet the environmental challenges of projected climate change. In this paper, we present the results of a long‐term infiltration experiment performed on fractured ash flow tuff. Water was released from a 3 × 4 m2 infiltration plot (divided into 12 square subplots) with a head of ∼0.04 m, over a period of ∼800 days. This experiment revealed peculiar infiltration patterns not amenable to current infiltration models, which were originally developed for infiltration into soils over a short duration. In particular, we observed that in part of the infiltration plot, the infiltration rate abruptly increased a few weeks into the infiltration tests. We suggest that these anomalies result from increases in fracture permeability during infiltration, which may be caused by swelling of clay fillings and/or erosion of infill debris. Interaction of the infiltration water with subsurface natural cavities (lithophysal cavities) could also contribute to such anomalies. This paper provides a conceptual model that partly describes the observed infiltration patterns in fractured rock and highlights some of the pitfalls associated with direct extension of soil infiltration models to fractured rock over a long period.

Mechanical degradation of emplacement drifts at Yucca Mountain—A modeling case study: Part II: Lithophysal rock

This paper outlines rock mechanics investigations associated with mechanical degradation of planned emplacement drifts at Yucca Mountain, which is the designated site for the proposed US high-level nuclear waste repository. The welded tuff emplacement horizon consists of two groups of rock with distinct engineering properties: nonlithophysal and lithophysal units, based on the relative proportion of lithophysal cavities. The term ‘lithophysal’ refers to hollow, bubble like cavities in volcanic rock that are surrounded by a porous rim formed by fine-grained alkali feldspar, quartz, and other minerals. Lithophysae are typically a few centimeters to a few decimeters in diameter. Part I of the paper concentrated on the degradation behavior of the generally hard, strong, and fractured nonlithophysal rock. Part II concentrates on the host rock in the lithophysal units. Lithophysal rock is characterized by lithophysal cavities interconnected by fracturing. Fracture sets are not as clearly defined as in the nonlithophysal rock. The rock mass porosity in the lithophysal units has been shown to be the primary physical factor that governs elastic and strength properties. The degradation behavior of the tunnels in the lithophysal rock is controlled by the spalling of the surrounding rock mass. A discontinuum material model has been developed to represent the drift scale rock mass properties in which slip and separation of contacting rock blocks can be estimated. Two-dimensional discontinuum analyses were developed with the consideration of in situ, thermal, and seismic loads. Time-dependent degradation was also considered. In this study, field and laboratory data and numerical analyses are well integrated to provide a solution for the unique problem of modeling drift degradation.

Mechanical degradation of emplacement drifts at Yucca Mountain—A modeling case study—Part I: Nonlithophysal rock

This paper outlines rock mechanics investigations associated with mechanical degradation of planned emplacement drifts at Yucca Mountain, which is the designated site for the proposed US high-level nuclear waste repository. The factors leading to drift degradation include stresses from the overburden, stresses induced by the heat released from the emplaced waste, stresses due to seismically related ground motions, and time-dependent strength degradation. The welded tuff emplacement horizon consists of two groups of rock with distinct engineering properties: nonlithophysal units and lithophysal units, based on the relative proportion of lithophysal cavities. The term ‘lithophysal’ refers to hollow, bubble like cavities in volcanic rock that are surrounded by a porous rim formed by fine-grained alkali feldspar, quartz, and other minerals. Lithophysae are typically a few centimeters to a few decimeters in diameter. Part I of the paper concentrates on the generally hard, strong, and fractured nonlithophysal rock. The degradation behavior of the tunnels in the nonlithophysal rock is controlled by the occurrence of keyblocks. A statistically equivalent fracture model was generated based on extensive underground fracture mapping data from the Exploratory Studies Facility at Yucca Mountain. Three-dimensional distinct block analyses, generated with the fracture patterns randomly selected from the fracture model, were developed with the consideration of in situ, thermal, and seismic loads. In this study, field data, laboratory data, and numerical analyses are well integrated to provide a solution for the unique problem of modeling drift degradation.

Fault-zone deformation in welded tuffs at Yucca Mountain, Nevada, USA

Field and microstructural analyses of faults with millimeters to hundreds of meters of displacement within welded tuffs at Yucca Mountain, Nevada, has led us to identify four different fault zone architectures. We designate these fault zones as Classes A–D. We conclude that Classes A, C, and D are genetically related and that observed differences in their morphology and deformation mechanisms are related to differences in displacement magnitude and changes in the rheology of faulted rocks. We show that Class A, C, and D fault zones were formed by progressive cataclasis, leading to the development of a foliated gouge core and wide damage zone at the highest displacement magnitudes. By inference, displacement variations across individual high displacement faults are expected to produce significant lateral changes in fault zone architecture. The broad range of fault rock textures and architectures found on the Class A, C, and D continuum may produce important corresponding variability in hydrologic properties in ignimbrites. By contrast, Class B fault zones have unique mineralogies and microtextures relative to other Yucca Mountain fault zones. Most notable is the occurrence of distinctive jigsaw puzzle texture within fault core breccias, which consists of angular wall rock fragments floating in a coarsely crystalline secondary calcite cement matrix. The calcite cement constitutes more than 65% of the Class B fault zones. Comparisons with other occurrences of secondary minerals within open lithophysal cavities, which appear to have formed in unsaturated conditions, indicate that the secondary mineralization history at Yucca Mountain was complex and likely polygenetic. The age and origin of the secondary calcite minerals in the Class B fault zones remain undetermined.

Reply to “Commentary: Assessment of past infiltration fluxes through Yucca Mountain on the basis of the secondary mineral record—is it a viable methodology?”, by Y.V. Dublyansky and S.Z. Smirnov

Xu et al. (2003) [Xu, T., Sonnenthal, E., Bodvarsson, G., 2003. A reaction-transport model for calcite precipitation and evaluation of infiltration–percolation fluxes in unsaturated fractured rock. J. Contam. Hydrol., 64, 113–127.] presented results of a reaction-transport model for calcite deposition in the unsaturated zone at Yucca Mountain, and compared the model results to measured abundances in core from a surface-based borehole. Marshall et al. (2003) [Marshall, B.D., Neymark, L.A., Peterman, Z.E., 2003. Estimation of past seepage volumes from calcite distribution in the Topopah Spring Tuff, Yucca Mountain, Nevada. J. Contam. Hydrol., 62–63, 237–247.] used the calcite distribution in the Topopah Spring Tuff to estimate past seepage into lithophysal cavities as an analog for seepage into the potential repository waste emplacement drifts at Yucca Mountain in southern Nevada (USA). Dublyansky and Smirnov (2005) [Dublyansky, Y.V., Smirnov, S.Z., 2005. Commentary: assessment of past infiltration fluxes through Yucca mountain on the basis of the secondary mineral record—is it a viable methodology? J. Contam. Hydrol. (this issue).] wrote a commentary paper to Marshall et al. (2003) [Marshall, B.D., Neymark, L.A., Peterman, Z.E., 2003. Estimation of past seepage volumes from calcite distribution in the Topopah Spring Tuff, Yucca Mountain, Nevada. J. Contam. Hydrol., 62–63, 237–247.] and Xu et al. (2003) [Xu, T., Sonnenthal, E., Bodvarsson, G., 2003. A reaction-transport model for calcite precipitation and evaluation of infiltration–percolation fluxes in unsaturated fractured rock. J. Contam. Hydrol., 64, 113–127.], containing two points: (1) questionable phenomenological model for the secondary mineral deposits and (2) inappropriate thermal boundary conditions. In this reply we address primarily the modeling approach by showing results of a sensitivity simulation regarding the effect of an elevated temperature history that approximates the temperature history inferred from fluid inclusions by Wilson et al. (2003) [Wilson, N.S.F., Cline, J.S., Amelin, Y.V., 2003. Origin, timing, and temperature of secondary calcite–silica mineral formation at Yucca Mountain, Nevada. Geochimica et Cosmochimica Acta, 67 (6), 1145–1184.]. Modeled calcite abundances using the time-varying temperature history are similar to the results for the steady-state ambient temperature profile ( Xu, T., Sonnenthal, E., Bodvarsson, G., 2003. A reaction-transport model for calcite precipitation and evaluation of infiltration–percolation fluxes in unsaturated fractured rock. J. Contam. Hydrol., 64, 113–127), and are still consistent with the measured abundances at the proposed repository horizon.

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture–matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault–matrix and fracture–matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture–matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.

Porosity dependence of the elastic modulus of lithophysae-rich tuff: numerical and experimental investigations

Abstract The influence of porosity on the mechanical properties of rock has received much research attention. Portions of the Yucca Mountain high-level nuclear waste repository may be placed in tuff units containing lithophysal cavities, which are large, generally noninterconnected cavities that can be considered a form of macroscopic porosity. This paper presents the results of numerical modeling and uniaxial compression testing of analog models and tuff rock, in order to assess the relationships between elastic modulus and porosity. The first part of the paper presents numerical simulation of uniaxial compression testing to calculate the elastic modulus of two-dimensional models containing randomly distributed circular holes in plane strain. The range of porosities investigated is approximately 5–40%. In the second part, the elastic modulus determined from the uniaxial compression testing of analog models and tuff specimens is presented. The analog specimens were made of plaster of Paris containing varying amounts of spherical shaped Styrofoam® inclusions to simulate a cavity structure similar to tuff. The results from the numerical analysis and analog material testing show an exponential decrease in elastic modulus with increasing porosity, whereas the elastic moduli of tuff show a linear decrease. The difference in the two behaviors can be attributed to the nonuniform cavity shapes in the tuff specimens.

Measurement system for systematic hydrological characterization of unsaturated fractured welded tuff in a mined underground tunnel.

A field investigation of unsaturated flow through a lithophysal unit of fractured welded tuff containing lithophysal cavities has been initiated. To characterize flow in this spatially heterogeneous medium, a systematic approach has been developed to perform tests in boreholes drilled at regular intervals in an underground tunnel (drift). The purpose of the testing is to quantify the amounts of water seeping into the drift versus the amount of water moving around the drift when released into boreholes at many equidistant locations along the drift. In this paper, we describe the test equipment system that has been built for this purpose. Because the field-scale measurements--of liquid flow in the unsaturated, fractured rocks--require continuous testing for periods of days to weeks, the control of test equipment has been fully automated, allowing operation with no human presence at the field site. Preliminary results from the first set of tests indicate that, while the effects of evaporation on characterization of hydrological properties of the rock can be significant, these effects can be controlled and quantified. These tests give insight into the role of the cavities as potential storage during the initial transient flow prior to the breakthrough of water at the drift crown, as well as the role of connected fractures that provide the subsequent quasi-steady flow. In addition to the stated purpose of realizing the flow partitioning, the results yield values for the effective porosity in the pathways for liquid flow in the regions tested thus far.

A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock

The percolation flux in the unsaturated zone (UZ) is an important parameter addressed in site characterization and flow and transport modeling of the potential nuclear-waste repository at Yucca Mountain, NV, USA. The US Geological Survey (USGS) has documented hydrogenic calcite abundances in fractures and lithophysal cavities at Yucca Mountain to provide constraints on percolation fluxes in the UZ. The purpose of this study was to investigate the relationship between percolation flux and measured calcite abundances using reactive transport modeling. Our model considers the following essential factors affecting calcite precipitation: (1) infiltration, (2) the ambient geothermal gradient, (3) gaseous CO 2 diffusive transport and partitioning in liquid and gas phases, (4) fracture–matrix interaction for water flow and chemical constituents, and (5) water–rock interaction. Over a bounding range of 2–20 mm/year infiltration rate, the simulated calcite distributions capture the trend in calcite abundances measured in a deep borehole (WT-24) by the USGS. The calcite is found predominantly in fractures in the welded tuffs, which is also captured by the model simulations. Simulations showed that from about 2 to 6 mm/year, the amount of calcite precipitated in the welded Topopah Spring tuff is sensitive to the infiltration rate. This dependence decreases at higher infiltration rates owing to a modification of the geothermal gradient from the increased percolation flux. The model also confirms the conceptual model for higher percolation fluxes in the fractures compared to the matrix in the welded units, and the significant contribution of Ca from water–rock interaction. This study indicates that reactive transport modeling of calcite deposition can yield important constraints on the unsaturated zone infiltration-percolation flux and provide useful insight into processes such as fracture–matrix interaction as well as conditions and parameters controlling calcite deposition.

Estimation of past seepage volumes from calcite distribution in the Topopah Spring Tuff, Yucca Mountain, Nevada

Low-temperature calcite and opal record the past seepage of water into open fractures and lithophysal cavities in the unsaturated zone at Yucca Mountain, Nevada, site of a proposed high-level radioactive waste repository. Systematic measurements of calcite and opal coatings in the Exploratory Studies Facility (ESF) tunnel at the proposed repository horizon are used to estimate the volume of calcite at each site of calcite and/or opal deposition. By estimating the volume of water required to precipitate the measured volumes of calcite in the unsaturated zone, seepage rates of 0.005 to 5 liters/year (l/year) are calculated at the median and 95th percentile of the measured volumes, respectively. These seepage rates are at the low end of the range of seepage rates from recent performance assessment (PA) calculations, confirming the conservative nature of the performance assessment. However, the distribution of the calcite and opal coatings indicate that a much larger fraction of the potential waste packages would be contacted by this seepage than is calculated in the performance assessment.

Origin, timing, and temperature of secondary calcite–silica mineral formation at Yucca Mountain, Nevada

The origin of secondary calcite–silica minerals in primary and secondary porosity of the host Miocene tuffs at Yucca Mountain has been hotly debated during the last decade. Proponents of a high-level nuclear waste repository beneath Yucca Mountain have interpreted the secondary minerals to have formed from cool, descending meteoric fluids in the vadose zone; critics, citing the presence of two-phase fluid inclusions, argued that the minerals could only have formed in the phreatic zone from ascending hydrothermal fluids. Understanding the origin, temperature, and timing of these minerals is critical in characterizing geologically recent fluid flux at the site, and has significant implications to whether waste should be stored at Yucca Mountain. Petrographic and paragenetic studies of 155 samples collected from the Exploratory Studies Facility (ESF) and repository block cross drift (ECRB) tunnels indicate that heterogeneously distributed calcite with lesser chalcedony, quartz, opal, and fluorite comprise the oldest secondary minerals. These are typically overgrown by intermediate-aged calcite, often exhibiting distinctive bladed habits. The youngest event recorded across the site is the deposition of Mg-enriched (up to ∼1 wt%) and depleted, growth-zoned calcite intergrown with U-enriched opal. The cyclical variation in Mg enrichment and depletion is probably related to climate changes that have occurred during the last few million years. The distribution of secondary minerals is consistent with precipitation in the vadose zone. Fluid inclusion petrography of sections from the 155 samples determined that 96% of the fluid inclusion assemblages (FIAs) contained liquid-only inclusions that formed at ambient temperatures (<35°C). However, 50% of the samples (n = 78) contained relatively rare FIA that contain both liquid-only and liquid plus vapor inclusions (herein termed two-phase FIAs) that formed at temperatures above 35°C. Virtually all of these two-phase FIAs occur in paragenetically old calcite; rare two-phase inclusion assemblages were also observed in early fluorite and quartz, and early–intermediate calcite. Homogenization temperatures (≡ trapping temperatures) across Yucca Mountain are generally 45 to 60°C, but higher temperatures reaching 83°C were recorded in calcite from the north portal and ramp of the ESF. Cooler temperatures of ∼35 to 45°C were recorded in the intensely fractured zone. Multiple populations of two-phase FIAs from lithophysal cavities in the ESF and ECRB cross drift indicate early fluid cooling with time from temperatures >45°C in early calcite, to <35 to 45°C in paragenetically younger calcite. Freezing point depressions range from −0.2 to −1.6°C, indicating trapping of a low salinity fluid. The majority of intermediate calcite and all outermost Mg-enriched calcite contains rare all-liquid inclusions and formed from ambient temperature (<35°C) fluids. Carbon and oxygen isotope data reveal a consistent trend of decreasing δ 13C (from 9.5 to −8.5‰) and increasing δ 18O (from 5.2 to 22.1‰) values from paragenetically early calcite to Mg-enriched growth-zoned calcite. Depleted δD values (−131 to −90‰) of inclusion fluids from intermediate and the youngest Mg-enriched calcite indicate derivation from surface meteoric fluids. Recalculation of δ 18O H 2 O values of −12 to −10‰ is consistent with derivation from paleometeoric fluids. Results of integrated U-Pb dating (opal and chalcedony) and fluid inclusion microthermometry indicate that two-phase FIAs that trapped fluids of >50°C are older than 6.29 ± 0.30 Ma. Two-phase FIAs in paragenetically later calcite, which formed from fluids of 35 to 45°C, are older than 5.32 ± 0.02 Ma. There is no evidence for trapping of fluids with elevated temperatures during the past 5.32 my. The youngest Mg-enriched calcite intergrown with opal began to precipitate between about 1.9 to 2.9 Ma and has continued to precipitate within the past half million years. The presence of liquid-only inclusions and the consistent occurrence of Mg-enriched calcite and opal as the youngest event indicate a minor, but chemically distinct, ambient temperature (<35°C) fluid flux during the past 2 to 3 my.

Physical and stable-isotope evidence for formation of secondary calcite and silica in the unsaturated zone, Yucca Mountain, Nevada

Calcite and silica form coatings on fracture footwalls and cavity floors in the welded tuffs at Yucca Mountain, the potential site of a high-level radioactive waste repository. These secondary mineral deposits are heterogeneously distributed in the unsaturated zone (UZ) with fewer than 10% of possible depositional sites mineralized. The paragenetic sequence, compiled from deposits throughout the UZ, consists of an early-stage assemblage of calcite±fluorite±zeolites that is frequently capped by chalcedony±quartz. Intermediate- and late-stage deposits consist largely of calcite, commonly with opal on buried growth layers or outermost crystal faces of the calcite. Coatings on steep-dipping fractures usually are thin (⩽3 mm) with low-relief outer surfaces whereas shallow-dipping fractures and lithophysal cavities typically contain thicker, more coarsely crystalline deposits characterized by unusual thin, tabular calcite blades up to several cms in length. These blades may be capped with knobby or corniced overgrowths of late-stage calcite intergrown with opal. The observed textures in the fracture and cavity deposits are consistent with deposition from films of water fingering down fracture footwalls or drawn up faces of growing crystals by surface tension and evaporated at the crystal tips. Fluid inclusion studies have shown that most early-stage and some intermediate-stage calcite formed at temperatures of 35 to 85 °C. Calcite deposition during the past several million years appears to have been at temperatures <30 °C. The elevated temperatures indicated by the fluid inclusions are consistent with temperatures estimated from calcite δ 18O values. Although others have interpreted the elevated temperatures as evidence of hydrothermal activity and flooding of the tuffs of the potential repository, the authors conclude that the temperatures and fluid-inclusion assemblages are consistent with deposition in a UZ environment that experienced prolonged heat input from gradual cooling of nearby plutons. The physical restriction of the deposits (and, therefore, fluid flow) to fracture footwalls and cavity floors and the heterogeneous and limited distribution of the deposits provides compelling evidence that they do not reflect flooding of the thick UZ at Yucca Mountain. The textures and isotopic and chemical compositions of these mineral deposits are consistent with deposition in a UZ setting from meteoric waters percolating downward along fracture flow paths.

U-Pb ages of secondary silica at Yucca Mountain, Nevada: implications for the paleohydrology of the unsaturated zone

Uranium, Th and Pb isotopes were analyzed in layers of opal and chalcedony from individual mm- to cm-thick calcite and silica coatings at Yucca Mountain, Nevada, USA, a site that is being evaluated for a potential high-level nuclear waste repository. These calcite and silica coatings on fractures and in lithophysal cavities in Miocene-age tuffs in the unsaturated zone (UZ) precipitated from descending water and record a long history of percolation through the UZ. Opal and chalcedony have high concentrations of U (10 to 780 ppm) and low concentrations of common Pb as indicated by large values of 206Pb/ 204Pb (up to 53,806), thus making them suitable for U-Pb age determinations. Interpretations of U-Pb isotope systems in opal samples at Yucca Mountain are complicated by the incorporation of excess 234U at the time of mineral formation, resulting in reverse discordance of U-Pb ages. However, the 207Pb/ 235U ages are much less affected by deviation from initial secular equilibrium and provide reliable ages of most silica deposits between 0.6 and 9.8 Ma. For chalcedony subsamples showing normal age discordance, these ages may represent minimum times of deposition. Typically, 207Pb/ 235U ages are consistent with the microstratigraphy in the mineral coating samples, such that the youngest ages are for subsamples from outer layers, intermediate ages are from inner layers, and oldest ages are from innermost layers. 234U and 230Th in most silica layers deeper in the coatings are in secular equilibrium with 238U, which is consistent with their old age and closed system behavior during the past ∼0.5 Ma. The ages for subsamples of silica layers from different microstratigraphic positions in individual calcite and silica coating samples collected from lithophysal cavities in the welded part of the Topopah Spring Tuff yield slow long-term average growth rates of 1 to 5 mm/Ma. These data imply that the deeper parts of the UZ at Yucca Mountain maintained long-term hydrologic stability over the past 10 Ma. despite significant climate variations. U-Pb ages for subsamples of silica layers from different microstratigraphic positions in individual calcite and silica coating samples collected from fractures in the shallower part of the UZ (welded part of the overlying Tiva Canyon Tuff) indicate larger long-term average growth rates up to 23 mm/Ma and an absence of recently deposited materials (ages of outermost layers are 3–5 Ma.). These differences between the characteristics of the coatings for samples from the shallower and deeper parts of the UZ may indicate that the nonwelded tuffs (PTn), located between the welded parts of the Tiva Canyon and Topopah Spring Tuffs, play an important role in moderating UZ flow.

Scientific Serials

Scientific Serials