Engineering topics

Gains-Germain, Leslie

Publications and source records attributed to Gains-Germain, Leslie.

9 records

DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

A Preliminary Radiological Risk Assessment Model for Disposition of Remote-Handled Transuranic Wastes at Los Alamos National Laboratory Area G - 20116

The U.S. Department of Energy (DOE) operates a low-level radioactive waste (LLW) disposal site at Material Disposal Area G, in Los Alamos, New Mexico, USA. Area G has been the primary LLW disposal site for Los Alamos National Laboratory (LANL) since the 1960's. In addition to LLW, Area G is host to a variety of other wastes, the disposition of which must be determined before closure of the site. A probabilistic Radiological Risk Assessment (RRA) for Area G is used in order to support decision making regarding some wastes that are not addressed in the extant Area G Performance Assessment (PA) and Composite Analysis (CA). Between 1979 and 1987, 33 special shafts were augered into the Bandelier Tuff at Area G. This volcanic tuff is present across Pajarito Plateau on the eastern slopes of the Jemez Mountains, and varies widely in its consistency, from weakly indurated non-welded layers to welded layers that uphold the mesa cliffs of the Plateau. These mesas are home to LANL, Area G, and the townsites of Los Alamos and White Rock, with residences about 1400 m from Area G. The 33 Shafts were lined with steel casing, and contain remote-handled (RH) transuranic wastes (TRU) resulting from experiments and analysis performed in special glove boxes at the Chemistry and Metallurgy Research (CMR) facility at LANL. Some of these wastes originated as used nuclear fuel. The purpose of the Area G RRA is to evaluate the potential future risk to humans and the environment from the RH TRU in the 33 Shafts in the context of the risk associated with the surrounding wastes at Area G. The analysis is responsive to expectations outlined in DOE Order 458.1, Radiation Protection of the Public and the Environment, and is informed by the Manual and Guidance accompanying DOE O 435.1, Radioactive Waste Management. Because the waste meets the definition of TRU, the regulatory context necessarily takes into consideration the regulation governing the disposal of TRU from the U.S. Environmental Protection Agency (EPA): 40 CFR 191, Environmental Radiation Protection Standards for Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Radioactive Wastes. Given the broader regulatory context for the RRA, the analysis is subject to different assumptions from those made in the existing DOE O 435.1 PA and CA, such as allowing for future occupation of the site. The analysis begins with a comprehensive evaluation of features, events, processes, and exposure scenarios (FEPS) for Area G and the wastes it contains. These FEPSs are screened to eliminate from further consideration those of extremely low probability and/or consequence, and a conceptual site model (CSM) is subsequently developed. The scope and structure of the Area G RRA Model is informed by this CSM, and the Area G RRA Model is developed using the GoldSim systems analysis modeling platform. This paper presents the initial version of a defensible, transparent, and reasonably realistic model, which is based on the state of knowledge of the wastes, the site, and the FEPSs that govern contaminant transport from wastes into the environment and subsequent exposures to humans and other biota. Probabilistic model input distributions represent uncertainties inherent in the real and modeled systems. The results of the Area G RRA Model inform decisions regarding the disposition of the RH TRU in the 33 Shafts. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Radionuclide Risk Coefficient Uncertainty Distributions for Use in Performance Assessments - 19156

Traditionally, radiation dose has been used as the regulatory criterion for decision making in performance assessments (PAs), despite the concept of cancer risk arguably being a more easily understood endpoint than dose for the public, and a preferable metric for decision analysis. However, the United States Environmental Protection Agency (EPA) has developed risk parameters for use in risk assessments and PAs in the form of radionuclide-specific risk coefficients, or RCs, which are published in various Federal Guidance Reports (FGRs). RCs are employed as scalars to convert environmental radionuclide exposure (via ingestion, inhalation, and external exposure) estimates to incremental lifetime cancer risk (ILCR). ILCR is the additional cancer risk averaged over a lifetime, over-and-above all-cause cancer risk, associated with carcinogen exposure. Uncertainties associated with RCs are discussed in FGR 13 and in a joint report published by EPA and Oak Ridge National Laboratory, but statistical distributions for the RCs were not directly developed. The uncertainty assessments in the EPA/ONRL report were based on dosimetric modeling and elicitation of nine experts. From these sources, five uncertainty categories were defined based on ratios of the 5. and 95. quantiles of uncertainty. Thus, the representations of uncertainty defined by EPA/ORNL are not distributions in a traditional sense. Furthermore, the authors only addressed uncertainty in RCs for inhalation and ingestion pathways. We develop a methodology for assigning statistical distributions to EPA's published mortality RC values (i.e. to estimate the probability of dying from radiation-induced cancers) for the inhalation (slow, medium, and fast absorption types), ingestion (water, food, and soil), and external exposure pathways associated with all cancer types. We employed the same dosimetric concepts and assumptions as EPA/ORNL, including cancer risk (dose-response) distributional information from a previous expert elicitation, which was also employed by EPA/ORNL in their work. We assumed that uncertainty in RCs can be described by lognormal distributions based on the modeling conducted previously by EPA/ORNL, and that point estimates tabulated in FGR 13 represent median values of the distributions. We also extended EPA/ORNL's uncertainty ratios to statistical distributions for external exposure RCs. The results are a set of radionuclide-specific distributions that are applicable in PAs that focus on cancer mortality. These distributions are currently and will subsequently be employed in PAs at various sites, in addition to the required regulatory dose limits. Quantification of ILCRs and associated uncertainty in PAs will provide regulators and the public with additional information for good decision making. Future directions include further refinement of the distributions and development of statistical distributions for morbidity RCs, which are used to estimate the probability of developing, but not necessarily dying as a result of, radiation-induced cancer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

A Preliminary Probabilistic Performance Assessment Model for Radiological and Chemical Contamination at the West Valley Site, New York - 19295

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States at the WNYNSC. NFS built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (640 metric tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations halted in 1972 and never restarted, leaving behind radioactive and chemical wastes. The reprocessing wastes include high-level waste (HLW), so U.S. Department of Energy (DOE) was compelled to complete certain waste management activities under the West Valley Demonstration Project (WVDP) Act of 1980, including decommissioning of WVDP facilities. As collaborating agencies, NYSERDA and the DOE are tasked with making decisions about decommissioning and risk reduction for the West Valley Site. Neptune and Company, Inc. (Neptune) was contracted to develop a probabilistic performance assessment (PPA) model to assist the agencies in their decision making process for decommissioning the WVDP and WNYNSC. One important tool that is needed in order to inform the decision-making process is a science-based model of the contamination at the West Valley Site and its potential implications for human health and the environment. The West Valley PPA Model (the Model), developed using the GoldSim modeling software, is a tool intended to provide support for decision making with uncertainty, in a manner that is transparent, defensible, and robust. The PPA Model includes contaminant transport and health effects components, and is organized around geographically-organized contaminated facilities. These include the waste disposal areas licensed by the U.S. Nuclear Regulatory Commission and the State of New York, a waste tank farm for storage of HLW resulting from reprocessing operations, and several areas contaminated with radioactive and chemical constituents. The Model evaluates contaminant transport from these sources to points of exposure across the site and into receiving surface waters and sediments. Hypothetical people and wildlife are exposed to contamination at these locations, and the effects of these exposures are evaluated. Contaminant transport processes to be evaluated in the Model include groundwater and surface water transport, contaminant translocation by plants and animals, diffusion in subsurface air and water, and erosion of the Site. The evaluation of exposures to people in this preliminary model is limited to a resident farmer scenario, and ecological assessment is performed at the level of a screening analysis. The results of these preliminary evaluations inform future model developments. Model results are subjected to sensitivity analysis in order to determine those pathways and parameters that are most significant in influencing the results. This information allows analysts and decision makers to focus on those aspects of Site behavior and processes that contribute most to overall uncertainty. With this information, the Supplemental Environmental Impact Statement (SEIS) can inform defensible decisions regarding remediation of the Site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Development of Waste Inventory for Probabilistic Performance Assessment Modeling - 19361

Estimation of radioactive waste inventory is a key initial step in the development of an environmental fate and transport model for performance assessment (PA) of radioactive waste sites. PA model output is often sensitive to parameters that represent inventory as a source term, yet estimation is a difficult problem. Inventory estimation for a probabilistic performance assessment (PPA) is related to, but distinct from, that performed for waste excavation, characterization, and offsite disposal. The latter appropriately focuses on obtaining the most detailed information possible on specific hazards and locations; however, these questions cannot be fully answered without intrusive site characterization efforts. This type of characterization is often not feasible when assembling information regarding a source term in a PPA model. Development of waste inventory for a PPA is one piece of a larger hierarchical model. The physical process model represented by Goldsim computer code represents one level of the hierarchy, and the second level of the hierarchy involves estimation of inputs such as inventory, transport, and risk parameters. The estimation of inputs is not a straightforward problem. The parameters of a computer model often do not have an equivalence to the physical quantities they describe, and estimation of the physical quantities themselves can be difficult due to limitations of available information. Inventory information for the West Valley Site in upstate New York is based on data from historical shipping and disposal forms. In the 1990's, the West Valley Demonstration Project (WVDP) compiled information from existing databases and historical waste records into one complete modern database. A conceptual model was then applied and used with recorded information to characterize the waste mass and radionuclide composition for each waste burial. Information was ultimately aggregated over all waste records to estimate the mass of each radionuclide within each trench, hole, or caisson. Three decades later, we revisit the master database, make updates with newly-identified information, and incorporate the conceptual model previously used for inventory estimation in the hierarchical modeling framework. We characterize uncertainty using probability distributions for waste mass and radionuclide composition at the individual shipment or individual burial level. Uncertainty incorporated at the lowest level of the hierarchy is ultimately carried through to uncertainty in PPA model endpoints. Sensitivity analysis of early model results is then applied to identify influential inputs and guide subsequent efforts towards parameters that drive risk. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Explicitly Accounting for Multiple Spatial Scales in Data Informing Kd Distribution Development - 18414

This work provides an example of incorporating information about spatial scale into input parameter distribution development for probabilistic performance assessment (PA) models and highlights the importance of having an interdisciplinary team working together in the distribution development process. The novel approach is demonstrated for contaminant partitioning coefficients for water/rock (Kds), providing a specific context for discussing the approach and challenges, but the strategy is general and applies to many other PA model inputs. In this example, a PA model of a radioactive waste disposal site is assumed to have large cells representing spatial scales of up to thousands of square meters in area and several to tens of meters in depth. Spatial variability in Kds is expected to occur at scales much smaller than the units in the model domain over which values are applied, and therefore values are intended to represent possible averages over the modeled volume (i.e., spatial scale represented by a PA model cell). Information about possible Kd values for each material category comes from a variety of literature sources, and the reported results used to inform the distribution are available at a wide range of scales (e.g., individual batch experiment results from small-scale samples, averages based on samples collected over site-scales, averages based on samples from large regions, etc.). Effectively addressing the mismatch in spatial scale between available information and the PA Model is a challenge within PA distribution development and it takes effort from database managers, subject matter experts, modelers, and statisticians. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Capturing the 'Representative Person' Using Exposure Equations and Input Variable Distributions - 18416

Radiological performance assessments (PAs) for radioactive waste employ the concept of future human receptors, or persons potentially exposed to radiation releases. The International Commission on Radiological Protection (ICRP) has defined the 'representative person' as a receptor, to replace the earlier concept of an 'average member of the critical group', for use in individual dose assessments and comparisons to regulatory criteria. A representative person can be assessed in PAs by defining a critical group, which is a relatively highly exposed subpopulation; and then estimating average exposures, doses, and risks for that group. Probabilistic PA models employ probability distributions to characterize uncertainty in exposure characteristics/variables at a scale aligned with contaminant transport and other aspects of the PA model. The concept of a representative person from a critical group, in the context of the first stage of an iterative approach and a probabilistic PA model, is illustrated using a receptor scenario and a set of exposure equations that incorporate probability distributions for input variables. Different scenarios may require different equations and distributions, but the approach taken is broadly generalizable. PA practitioners and regulators may benefit from such a generalizable method, as it may simplify approaches as well as facilitate reviews and comparisons between different PAs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗