Search NASA⌕ Search

SEARCH · Search NASA

Results for “tracer flow & transport”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Utilizing Environmental Tracers to Reduce Groundwater Flow and Transport Model Parameter Uncertainties

Abstract Non‐uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium ( 3 H), chlorofluorocarbon‐12 (CFC‐12), and sulfur hexafluoride (SF 6 ) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady‐state 2‐D synthetic aquifer and a transient 3‐D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post‐calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration‐constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

HdsaLib

HdsaLib performs hyper-differential sensitivity analysis in optimization problems that are constrained by partial differential equation problems. The software aids researchers in understanding the impact of uncertainties on their optimization models. The tool has been applied for use in ice sheet modeling, fluid flow modeling, atmospheric tracer transport modeling, and subsurface flow modeling. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

van Bloemen Waanders, Bart [Sandia National Lab. (↗
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗
DOE OSTI2024

Sub-Diurnal Methane Variations on Mars Driven by Barometric Pumping and Planetary Boundary Layer Evolution

In recent years, the Tunable Laser Spectrometer within the Sample Analysis at Mars (TLS-SAM) instrument on board the Mars Science Laboratory (MSL) Curiosity rover has detected methane variations in the atmosphere at Gale crater. Methane concentrations appear to fluctuate seasonally as well as sub-diurnally, which is difficult to reconcile with an as-yet-unknown transport mechanism delivering the gas from underground to the atmosphere. To potentially explain the fluctuations, we consider barometrically induced transport of methane from an underground source to the surface, modulated by temperature-dependent adsorption. The subsurface fractured-rock seepage model is coupled to a simplified 1-D atmospheric mixing model to provide insights on the pattern of atmospheric methane concentrations in response to transient surface methane emissions, as well as to predict sub-diurnal variation in methane abundance for the northern summer period, which is a candidate time frame for a MSL Curiosity sampling campaign. Our analysis suggests that there is a lower limit to the subsurface fracture density that can produce the observed methane patterns, below which the atmospheric methane variations would be out of phase with the observations. The best-performing model scenarios indicate a significant, short-lived methane pulse just prior to sunrise, the detection of which by TLS-SAM would be a potential indicator of the contribution of barometric pumping to Mars' atmospheric methane variations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling framework for evaluating the impacts of hydrodynamic pressure on hydrologic exchange fluxes and residence time for a large-scale river section over a long-term period

Quantifying hydrologic exchange fluxes (HEF) at the river and subsurface interface and their residence times (RT) in subsurface are important for managing the water quality and ecosystem health in dynamic river corridor systems. In this study, a modeling framework is developed for coupling the three-dimensional (3D) multi-phase surface, subsurface flow transport, and numerical tracer model for RT in a large-scale river section over a long period. The framework is utilized to evaluate the impacts of hydrodynamic pressure on HEFs and RT for a 30 km section of the Columbia River in Washington State over a three-year period. Based on comparisons between model simulations with and without considering hydrodynamic pressure, we found that hydrodynamic pressure increases the net HEFs by 7% with river gaining water from the subsurface domain, and leads to slight reduction of RT.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characterization of flow and transport in a fracture network at the EGS Collab field experiment through stochastic modeling of tracer recovery

Energy extraction from subsurface reservoirs is important for addressing the increasing energy demand and environmental concerns such as global warming. However, the characterization of subsurface reservoirs, particularly reservoirs dominated by fracture networks remains a challenge due to the lack of means to directly observe subsurface processes. This study explores the feasibility and efficacy of characterizing fracture flow and transport processes in an enhanced geothermal system (EGS) testbed through stochastic tracer modeling. There are two enabling factors that allow application of stochastic modeling to characterize a subsurface reservoir. First, an abundance of geological and geophysical measurements enables the development of a high-fidelity and well-constrained fracture network model. Second, high-performance computing (HPC) allows running massive realizations efficiently. Six conservative tracer tests were stochastically modeled and produced satisfactory realizations that successfully reproduce field tracer recovery data from each tracer test. The evolution of flow and transport processes in the fracture network was then analyzed from these satisfactory realizations. This work demonstrates that stochastic tracer modeling on a high-fidelity fracture network model is feasible and can provide important insights regarding flow and transport characteristics in subsurface fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

Understanding the relative importance of vertical and horizontal flow in ice-wedge polygons

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Surface-to-Interior Transport Timescales and Ventilation Patterns in a Time-Dependent Circulation Driven by Sustained Climate Warming

Abstract The effect of climate warming in response to rising atmospheric CO 2 on the ventilation of the ocean remains uncertain. Here we make theoretical advances in elucidating the relationship between ideal age and transit time distribution (TTD) in a time-dependent flow. Subsequently, we develop an offline tracer-transport model to characterize the ventilation patterns and time scales in the time-evolving circulation for the 1850–2300 period as simulated with the Community Earth System Model version 1 (CESMv1) under a business-as-usual warming scenario. We found that by 2300 2.1% less water originates from the high-latitude deep water formation regions (both hemispheres) compared to 1850. In compensation, there is an increase in the water originating from the subantarctic. We also found that slowing meridional overturning circulation causes a gradual increase in mean age during the 1850–2300 period, with a globally averaged mean-age increase of ∼110 years in 2300. Where and when the water will be re-exposed to the atmosphere depends on the post-2300 circulation. For example, if we assume that the circulation persists in its year-2300 state (scenario 1), the mean interior-to-surface transit time in year 1850 is ∼1140 years. In contrast, if we assume that the circulation abruptly recovers to its year-1850 state (scenario 2), the mean interior-to-surface transit time in 1850 is only ∼740 years. By 2300, these differences become even larger; in scenario 1, the mean interior-to-surface transit time increases by ∼200 years, whereas scenario 2 decreases by ∼80 years. The dependence of interior-to-surface transit time on the future ocean circulation produces an additional unavoidable uncertainty in the long-term durability of marine carbon dioxide removal strategies. Significance Statement The ocean’s circulation, when altered by climate warming, can affect its capacity to absorb heat and CO 2 , which are crucial for the global climate. In our study, we investigated how global warming, caused by rising CO 2 levels, might impact the ocean circulation—the way water moves from deep ocean to the surface and vice versa. We discovered that by 2300, if we continue on our current warming trajectory, the origins of water within the ocean will shift, with less coming from deep, cold zones near the poles and more from subantarctic regions. As a result, deep water will take longer time before it resurfaces than shallow water. How quickly this water travels from deep regions to the surface could change, depending on the state of future ocean circulation. If the circulation remains as predicted in 2300, this journey will take longer. Conversely, if it reverts to the pattern in 1850, the process will be quicker. This variability introduces added uncertainty to strategies aimed at mitigating climate change by storing CO 2 in the ocean. Our work highlights the intricate ways in which climate change can influence our oceans, potentially affecting our plans to mitigate global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electrical monitoring of saline tracers to reveal subsurface flow pathways in a flat ditch-drained field

Time-lapse electrical resistivity imaging (ERI) combined with salt tracers can enhance our understanding of subsurface flow and transport in ditch-drained agroecosystems. In a nearly level field adjacent to a drainage ditch, we installed 192 electrodes in a 72-m2 plot along with five hydrometric stations to enable time-lapse ERI monitoring during storms. Immediately prior to a wintertime rainstorm, we applied 5 kg of KBr to a shallow trench located roughly 11 m upgradient from the ditch. Antecedent conditions prior to the storm were wet, assuring a hydrologically connected state that favored tracer dissolution and movement in the subsurface. We monitored the evolution of the tracer plume with ERI over a seven-day period. Four-dimensional (4D) delineation of precipitation and tracer transport using time-lapse ERI showed that the tracer plume moved preferentially downward through a discontinuity in the argillic horizon within the first 57 hrs. Over the following 79 hrs, the tracer plume exhibited steady, lateral movement towards the drainage ditch within a more permeable sandy horizon situated just below the argillic. A spatial moment analysis of logarithmic changes in conductivity in the sand layer revealed seepage velocities ranging from 0.0027 to 0.0482 m d-1; these values were positively correlated with hydraulic gradients, suggesting laminar flow processes typical of sandy aquifers. Increasing groundwater velocities with hydraulic gradients suggest that solute transfers from groundwater to ditch are likely to be faster in storms, especially during wet periods. Even so, predominantly slow lateral transport within the sandy horizon, even during storms, indicates that critical source areas of pollutant loss are likely within a few meters of the ditch.

Robinson, Judith L.↗
NASA NTRS1980

Reactive mixing in swirling flows

The effect of swirl on the mixing of two gas streams of the same density and axial velocity in a coaxial annular geometry has been examined. The flow is studied via a fluorescent tracer, 2-3 biacetyl. Through direct excitation, collisional excitation, and collisional de-excitation of the tracer, the turbulent transport, molecularly mixed and unmixed regions in the flow are visualized. The visualization results indicate that there are interesting three dimensional inviscid large scale structures. Mixing in the molecular level is found to be mainly localized in these structures.

Cheng, W. K.↗
DOE OSTI2022

Tracer Gas Model Development and Verification in PFLOTRAN

Tracer gases, whether they are chemical or isotopic in nature, are useful tools in examining the flow and transport of gaseous or volatile species in the underground. One application is using detection of short-lived argon and xenon radionuclides to monitor for underground nuclear explosions. However, even chemically inert species, such as the noble gases, have bene observed to exhibit non-conservative behavior when flowing through porous media containing certain materials, such as zeolites, due to gas adsorption processes. This report details the model developed, implemented, and tested in the open source and massively parallel subsurface flow and transport simulator PFLOTRAN for future use in modeling the transport of adsorbing tracer gases.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS2019

Large-Scale Transport into the Arctic: The Roles of the Midlatitude Jet and the Hadley Cell

Transport from the Northern Hemisphere (NH) midlatitudes to the Arctic plays a crucial role in determining the abundance of trace gases and aerosols that are important to Arctic climate via impacts on radiation and chemistry. Here we examine this transport using an idealized tracer with a fixed lifetime and predominantly midlatitude land-based sources in models participating in the Chemistry Climate Model Initiative (CCMI). We show that there is a 25%-45% difference in the Arctic concentrations of this tracer among the models. This spread is correlated with the spread in the location of the Pacific jet, as well as the spread in the location of the Hadley Cell (HC) edge, which varies consistently with jet latitude. Our results suggest that it is likely that the HC-related zonal-mean meridional transport rather than the jet-related eddy mixing is the major contributor to the inter-model spread in the transport of land-based tracers into the Arctic. Specifically, in models with a more northern jet, the HC generally extends further north and the tracer source region is mostly covered by surface southward flow associated with the lower branch of the HC, resulting in less efficient transport poleward to the Arctic. During boreal summer, there are poleward biases in jet location in free-running models, and these models likely underestimate the rate of transport into the Arctic. Models using specified dynamics do not have biases in the jet location, but do have biases in the surface meridional flow, which may result in differences in transport into the Arctic. In addition to the land-based tracer, the midlatitude-to-Arctic transport is further examined by another idealized tracer with zonally uniform sources. With equal sources from both land and ocean, the inter-model spread of this zonally uniform tracer is more related to variations in parameterized convection over oceans rather than variations in HC extent, particularly during boreal winter. This suggests that transport of land-based and oceanic tracers or aerosols towards the Arctic differs in pathways and therefore their corresponding inter-model variabilities result from different physical processes.

chemistry climate models↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1999

Synoptic Mapping of Chemical Composition, Thermal Structure, and Air Motion from UARS Observations

An operational mapping algorithm was developed to process measurements of individual species observed by different satellite instruments on board UARS. Based on Fast Fourier Synoptic Mapping (FFSM), the algorithm accounts for the precessing orbit of UARS, the error of individual instruments, and gaps associated with instrument duty cycle and the satellite yaw maneuver. It provides synoptic structure and evolution on periods as short as 1 day, derived collectively from all observations of an individual species. The algorithm was applied to synoptically map temperature, thickness, and several chemical constituents observed by the instruments: MLS, CLAES, ISAMS, and HALOE. Each field variable observed by these instruments was mapped twice-daily in continuous global time series several months long. Mapped behavior produced via FFSM was compared against standard archived products generated via Kalman filtering. The standard map products reveal similar behavior, but are limited to features of larger scale. Tracer structure that develops through flow deformation and attending transport is therefore represented in those products more coarsely, eventually being sheared down to scales that are no longer represented properly. The synoptic time series also reveal a diurnal cycle for several of the constituents, one that emerges clearly in their space-time spectra. The zonal and meridional structure of diurnal variability was mapped by filtering the space-time spectrum, an intermediate product of FFSM, to those scales resolved by the UARS sampling. Geographical variations of the diurnal cycle have also been evaluated in its seasonal-mean structure by compositing distributions at individual local times. Both were compared to diurnal variations in chemical models and in dynamical models that account for transport by the diurnal tide. Diurnal variations were found to introduce spurious behavior into the archived products generated via Kalman filtering, behavior that is misrepresented in those synoptic time series. Aliasing by under-sampled diurnal variations introduces a bias into the time-mean structure of archived products. That systematic error is manifested as anomalous wave structure in the time-mean fields. By contrast, synoptic time series generated via FFSM faithfully represent such behavior as varying diurnally, so it does not corrupt time-mean structure. Tracer structure mapped via the above algorithm was being applied to determine the Brewer-Dobson circulation of the stratosphere through Lagrangian considerations. Preliminary results are encouraging: They reveal diabatic ascent in the tropics, balanced by diabatic descent over the winter hemisphere. This Lagrangian description of transport derived from UARS tracer measurements would be potentially valuable for diagnosing interannual changes of transport, specifically, in relation to interannual changes of total ozone. Regrettably, support was reduced to 50% of the approved funding level, leaving this component incomplete.

Salby, Murry L.↗
NASA NTRS1994

Tracer evolution in winds generated by a global spectral mechanistic model

The lower boundary of a spectral mechanistic model is prescribed with 100 hPa geopotentials, and its performance during a November 1989 through March 1990 integration is compared with National Meteorological Center observations. Although the stratopause temperatures quickly become biased near the pole in both hemispheres, the model develops a residual mean circulation which shows significant descent over the winter pole and ascent in the tropics and over the summer pole at pressures less than 10 hPa. The daily correspondence of observed to modeled features in the upper stratosphere and mesosphere degrades after one month. However, the long-term variability qualitatively follows the observations. The results of off-line transport experiments are also described. A passive tracer is instantaneously injected into the flow over the poles and evolves in a manner which is consistent with the residual mean circulation. It demonstrates a significant cross-equatorial flux in the mesosphere near solstice, and air which originates in the southern hemisphere polar mesosphere can be found descending deep into the nothern polar stratosphere at the end of the integration. Nitrous oxide is also transported, and its ability to act as a dynamical tracer is evaluated by comparison to the evolution of the passive tracer.

Nielsen, J. E.↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗