A comparison of groundwater fluxes computed with MODFLOW and a mixing model using deuterium: Application to the eastern Nevada Test Site and vicinity

被引:17
作者
Carroll, Rosemary W. H. [1 ]
Pohll, Greg M. [1 ]
Earman, Sam [1 ]
Hershey, Ronald L. [1 ]
机构
[1] Desert Res Inst, Div Hydrol Sci, Reno, NV 89512 USA
基金
美国能源部;
关键词
MODFLOW; Isotope mixing model; Nevada Test Site; Deuterium;
D O I
10.1016/j.jhydrol.2008.08.005
中图分类号
TU [建筑科学];
学科分类号
0813 ;
摘要
The primary objective of this study was to verify groundwater flows in the vicinity of the eastern Nevada Test Site (NTS) computed with a hydraulically defined flow model against groundwater delta D values computed by a steady-state mixing model. The United States Geological Survey's Death Valley regional flow model (DVRFM) is a transient, three-dimensional, groundwater model that uses the public domain, finite-difference code MODFLOW. The mixing model (Discrete-State Compartment Model-Shuffled Complex Evolution, or DSCM-SCE) is a recently developed code that is able to autocalibrate groundwater fluxes (both magnitude and direction) to best match observed tracer concentrations in groundwater. To compare modeling approaches, DVRFM boundary conditions and cell-to-cell interactions were implemented into a previously developed 15-cell DSCM-SCE delta D model of the eastern NTS. Analysis of delta D and delta(18) data conducted throughout the model domain suggests recharge and mixing are the dominant mechanisms for groundwater isotopic enrichment in the downgradient direction. Therefore, evaporation, at least at the regional scale, was ignored. Model results showed that DVRFM boundary fluxes and cell outflow volumes reproduced observed groundwater delta D values in nine of the 11 hydrographic basins that contained deuterium data. The remaining four modeled basins did not have any deuterium data, however, modeling results closely matched estimated delta D for two of these basins. The isotope mixing model. independently verified that most of the hydraulically defined groundwater flows simulated by MODFLOW are reasonable. Optimization of the DSCM-SCE was then done by adjusting groundwater fluxes between cells to improve predicted groundwater delta D values. This significantly towered the weighted root-mean-squared error and helped to identify several basins where DVRFM-simulated boundary conditions and groundwater fluxes should be reexamined. Overall, the two modeling approaches complemented each other, such that the isotope mixing model verified the hydraulically based groundwater flow model and refined groundwater flux estimates to better understand the physical system. (C) 2008 Elsevier B.V. All rights reserved.
引用
收藏
页码:371 / 385
页数:15
相关论文
共 46 条
  • [1] ANDERMAN ER, 2003, 03347 US GEOL SURV
  • [2] [Anonymous], THESIS U NEVADA RENO
  • [3] [Anonymous], 2000, US GEOL SURV OPEN FI
  • [4] [Anonymous], 994079 US GEOL SURV
  • [5] [Anonymous], TECHNICAL REPORT SER
  • [6] [Anonymous], THESIS U NEVADA RENO
  • [7] [Anonymous], 2000, MODFLOW 2000 US GEOL
  • [8] [Anonymous], THESIS U ARIZONA TUC
  • [9] [Anonymous], 20045205 US GEOL SUR
  • [10] BELCHER WR, 2004, 20045205 US GEOL SUR, P3