https://www.stt.msu.edu/users/mcubed/frade.pdf

Finite difference approximations for fractional advection-dispersion flow equations

Among the current problems that hydrogeologists face, perhaps there is none as challenging as the characterization of fractured rock (Faybishenko and Benson 2000). This paper discusses issues associated with the quantification of flow and transport through fractured rocks on scales not exceeding those typically associated with single- and multi-well pressure (or flow) and tracer tests. As much of the corresponding literature has focused on fractured crystalline rocks and hard sedimentary rocks such as sandstones, limestones (karst is excluded) and chalk, so by default does this paper. Direct quantification of flow and transport in such rocks is commonly done on the basis of fracture geometric data coupled with pressure (or flow) and tracer tests, which therefore form the main focus. Geological, geophysical and geochemical (including isotope) data are critical for the qualitative conceptualization of flow and transport in fractured rocks, and are being gradually incorporated in quantitative flow and transport models, in ways that this paper unfortunately cannot describe but in passing. The hydrogeology of fractured aquifers and other earth science aspects of fractured rock hydrology merit separate treatments. All evidence suggests that rarely can one model flow and transport in a fractured rock consistently by treating it as a uniform or mildly nonuniform isotropic continuum. Instead, one must generally account for the highly erratic heterogeneity, directional dependence, dual or multicomponent nature and multiscale behavior of fractured rocks. One way is to depict the rock as a network of discrete fractures (with permeable or impermeable matrix blocks) and another as a nonuniform (single, dual or multiple) continuum. A third way is to combine these into a hybrid model of a nonuniform continuum containing a relatively small number of discrete dominant features. In either case the description can be deterministic or stochastic. The paper contains a brief assessment of these trends in light of recent experimental and theoretical findings, ending with a short list of prospects and challenges for the future.

Trends, prospects and challenges in quantifying flow and transport through fractured rocks

[1] Two large-scale cross-hole pumping tests were conducted at depths of 191–226 m and 662–706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (Ss) tomograms, as well as their uncertainties in three dimensions The equivalent K and Ss estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort Results show several, distinct, high K and low Ss zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and Ss heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2007WR006715

Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan

https://digital.library.unt.edu/ark:/67531/metadc897878/m2/1/high_res_d/929003.pdf

Field-scale effective matrix diffusion coefficient for fractured rock : Results from literature survey

Enhancement of bedrock permeability by weathering

[1] A cross-hole tracer test involving the simultaneous injection of two nonsorbing solute tracers with different diffusion coefficients (bromide and pentafluorobenzoate) and a weakly sorbing solute tracer (lithium ion) was conducted in a fractured granite near an underground nuclear test cavity in central Nevada. The test was conducted to (1) test a conceptual radionuclide transport model for the site and (2) obtain transport parameter estimates for predictive modeling. The differences between the responses of the two nonsorbing tracers (when normalized to injection masses) are consistent with a dual-porosity transport system in which matrix diffusion is occurring. The large concentration attenuation of the sorbing tracer relative to the nonsorbing tracers suggests that diffusion occurs primarily into matrix pores, not simply into stagnant water within the fractures. The relative responses of the tracers at late times suggest that the diffusion-accessible matrix pore volume is possibly limited to only half the total volume of the flow system, implying that the effective retardation factor due to matrix diffusion may be as small as 1.5 for nonsorbing solutes in the system. The lower end of the range of possible sorption Kd values deduced from the lithium response is greater than the upper 95% confidence bound of Kd values measured in laboratory sorption tests using crushed granite from the site. This result suggests that the practice of using laboratory sorption data in field-scale transport predictions of cation-exchanging radionuclides, such as 137Cs+ and 90Sr++, should be conservative for the SHOAL site.

Testing and parameterizing a conceptual model for solute transport in a fractured granite using multiple tracers in a forced-gradient test

Multiple isotope forensics of nitrate in a wild horse poisoning incident.

/pdf/evaluation-of-chemical-data-from-samples-collected-february-1tg9el21up.pdf

Evaluation of Chemical Data from Samples Collected February 2008 at the Cactus Flat Main Lake Depression and Surrounding Area, Nevada Test and Training Range

Groundwater samples in the Yucca Mountain area were collected for chemical and isotopic analyses and measurements of water temperature, pH, specific conductivity, and alkalinity were obtained at the well or spring at the time of sampling. For this project, groundwater samples were analyzed for major-ion chemistry, deuterium, oxygen-18, and carbon isotopes of dissolved inorganic carbon (DIC) and dissolved organic carbon (DOC). The U.S. Geological Survey (USGS) performed all the fieldwork on this project including measurement of water chemistry field parameters and sample collection. The major ions dissolved in the groundwater, deuterium, oxygen-18, and carbon isotopes of dissolved inorganic carbon (DIC) were analyzed by the USGS. All preparation and processing of samples for DOC carbon isotopic analyses and geochemical modeling were performed by the Desert Research Institute (DRI). Analysis of the DOC carbon dioxide gas produced at DRI to obtain carbon-13 and carbon-14 values was conducted at the University of Arizona Accelerator Facility (a NSHE Yucca Mountain project QA qualified contract facility). The major-ion chemistry, deuterium, oxygen-18, and carbon isotopes of DIC were used in geochemical modeling (NETPATH) to determine groundwater sources, flow paths, mixing, and ages. The carbon isotopes of DOC were used to calculate groundwater ages that are independent of DIC model corrected carbon-14 ages. The DIC model corrected carbon-14 calculated ages were used to evaluate groundwater travel times for mixtures of water including water beneath Yucca Mountain. When possible, groundwater travel times were calculated for groundwater flow from beneath Yucca Mountain to down gradient sample sites. DOC carbon-14 groundwater ages were also calculated for groundwaters in the Yucca Mountain area. When possible, groundwater travel times were estimated for groundwater flow from beneath Yucca Mountain to down gradient groundwater sample sites using the DOC calculated groundwater ages. The DIC calculated groundwater ages were compared with DOC calculated groundwater ages and both of these ages were compared to travel times developed in ground-water flow and transport models. If nuclear waste is stored in Yucca Mountain, the saturated zone is the final barrier against the release of radionuclides to the environment. The most recent rendition of the TSPA takes little credit for the presence of the saturated zone and is a testament to the inadequate understanding of this important barrier. If radionuclides reach the saturated zone beneath Yucca Mountain, then there is a travel time before they would leave the Yucca Mountain area and flow down gradient to the Amargosa Valley area. Knowing how long it takes groundwater in the saturated zone to flow from beneath Yucca Mountain to down gradient areas is critical information for potential radionuclide transport. Radionuclide transport in groundwater may be the quickest pathway for radionuclides in the proposed Yucca Mountain repository to reach land surface by way of groundwater pumped in Amargosa Valley. An alternative approach to ground-water flow and transport models to determine the travel time of radionuclides from beneath Yucca Mountain to down gradient areas in the saturated zone is by carbon-14 dating of both inorganic and organic carbon dissolved in the groundwater. A standard method of determining ground-water ages is to measure the carbon-13 and carbon-14 of DIC in the groundwater and then correct the measured carbon-14 along a flow path for geochemical reactions that involve carbon containing phases. These geochemical reactions are constrained by carbon-13 and isotopic fractionations. Without correcting for geochemical reactions, the ground-water ages calculated from only the differences in carbon-14 measured along a flow path (assuming the decrease in carbon-14 is due strictly to radioactive decay) could be tens of thousands of years too old. The computer program NETPATH, developed by the USGS, is the best geochemical program for correcting carbon-14 activities for geochemical reactions. The DIC carbon-14 corrected ages can be further constrained by measuring the carbon isotopes of DOC. Because the only source of organic carbon in aquifers is almost always greater than 40,000 years old, any organic carbon that may be added to the groundwater would contain no carbon-14. Thus, ground-water ages determined by carbon isotopes of DOC should be maximum ages that can be used to constrain DIC corrected ages.

/pdf/yucca-mountain-area-saturated-zone-dissolved-organic-carbon-578mgm4zsj.pdf

Yucca Mountain Area Saturated Zone Dissolved Organic Carbon Isotopic Data

The purpose of this study is to quantify the flow and transport model uncertainty for the Central Nevada Test Area (CNTA). Six parameters were identified as uncertain, including the specified head boundary conditions used in the flow model, the spatial distribution of the underlying welded tuff unit, effective porosity, sorption coefficients, matrix diffusion coefficient, and the geochemical release function which describes nuclear glass dissolution. The parameter uncertainty was described by assigning prior statistical distributions for each of these parameters. Standard Monte Carlo techniques were used to sample from the parameter distributions to determine the full prediction uncertainty. Additional analysis is performed to determine the most cost-beneficial characterization activities. The maximum radius of the tritium and strontium-90 contaminant boundary was used as the output metric for evaluation of prediction uncertainty. The results indicate that combining all of the uncertainty in the parameters listed above propagates to a prediction uncertainty in the maximum radius of the contaminant boundary of 234 to 308 m and 234 to 302 m, for tritium and strontium-90, respectively. Although the uncertainty in the input parameters is large, the prediction uncertainty in the contaminant boundary is relatively small. The relatively small prediction uncertainty is primarily due to the small transport velocities such that large changes in the uncertain input parameters causes small changes in the contaminant boundary. This suggests that the model is suitable in terms of predictive capability for the contaminant boundary delineation.

/pdf/assessing-groundwater-model-uncertainty-for-the-central-t8jhuee2tx.pdf

Todd Mihevc

Papers

Testing and parameterizing a conceptual model for solute transport in a fractured granite using multiple tracers in a forced-gradient test

Multiple isotope forensics of nitrate in a wild horse poisoning incident.

Evaluation of Chemical Data from Samples Collected February 2008 at the Cactus Flat Main Lake Depression and Surrounding Area, Nevada Test and Training Range

Yucca Mountain Area Saturated Zone Dissolved Organic Carbon Isotopic Data

Assessing Groundwater Model Uncertainty for the Central Nevada Test Area