PHT3D‐UZF: A Reactive Transport Model for Variably‐Saturated Porous Media

A modified version of the MODFLOW/MT3DMS‐based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably‐saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D‐UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated‐zone flow (UZF1) package. A volume‐averaged approach similar to the method used in UZF‐MT3DMS was adopted. The PHREEQC‐based computation of chemical processes within PHT3D‐UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional‐scale applications, UZF1 simulates downward‐only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably‐saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.     

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

Fluid pressure responses for a Devil's Slide‐like system: problem formulation and simulation

This study employs a hydrogeologic simulation approach to investigate subsurface fluid pressures for a landslide‐prone section of the central California, USA, coast known as Devil's Slide. Understanding the relative changes in subsurface fluid pressures is important for systems, such as Devil's Slide, where slope creep can be interrupted by episodic slip events. Surface mapping, exploratory core, tunnel excavation records, and dip meter data were leveraged to conceptualize the parameter space for three‐dimensional (3D) Devil's Slide‐like simulations. Field observations (i.e. seepage meter, water retention, and infiltration experiments; well records; and piezometric data) and groundwater flow simulation (i.e. one‐dimensional vertical, transient, and variably saturated) were used to design the boundary conditions for 3D Devil's Slide‐like problems. Twenty‐four simulations of steady‐state saturated subsurface flow were conducted in a concept‐development mode. Recharge, heterogeneity, and anisotropy are shown to increase fluid pressures for failure‐prone locations by up to 18.1, 4.5, and 1.8% respectively. Previous estimates of slope stability, driven by simple water balances, are significantly improved upon with the fluid pressures reported here. The results, for a Devil's Slide‐like system, provide a foundation for future investigations. Copyright © 2014 John Wiley & Sons, Ltd.     

Understanding the water balance paradox in the Athabasca River Basin,Canada

This study demonstrates the importance of the including and appropriately parameterizing peatlands and forestlands for basin‐scale integrated surface–subsurface models in the northern boreal forest, with particular emphasis on the Athabasca River Basin (ARB). With a long‐term water balance approach to the ARB, we investigate reasons why downstream mean annual stream flow rates are consistently higher than upstream, despite the subhumid water deficit conditions in the downstream regimes. A high‐resolution 3D variably saturated subsurface and surface water flow and evapotranspiration model of the ARB is constructed based on the bedrock and surficial geology and the spatial distribution of peatlands and their corresponding eco‐regions. Historical climate data were used to drive the model for calibration against 40‐year long‐term average surface flow and groundwater observations during the historic instrumental period. The simulation results demonstrate that at the basin‐scale, peatlands and forestlands can have a strong influence on the surface–subsurface hydrologic systems. In particular, peatlands in the midstream and downstream regimes of the ARB increase the water availability to the surface–subsurface water systems by reducing water loss through evapotranspiration. Based on the comparison of forestland evapotranspiration between observation and simulation, the overall spatial average evapotranspiration in downstream forestlands is larger than that in peatlands and thus the water contribution to the stream flow in downstream areas is relatively minor. Therefore, appropriate representation of peatlands and forestlands within the basin‐scale hydrologic model is critical to reproduce the water balance of the ARB.     

A kinetic approach for simulating redox‐controlled fringe and core biodegradation processes in groundwater: model development and application to a landfill site in Piedmont,Italy

A three‐dimensional model for predicting redox controlled, multi‐species reactive transport processes in groundwater systems is presented. The model equations were fully integrated within a MODFLOW‐family reactive transport code, RT3D. The model can simulate organic compound biodegradation coupled to different terminal electron acceptor processes. A computational approach, which uses the spatial and temporal distribution of the rates of different redox reactions, is proposed to map redox zones. The method allows one to quantify and visualize the biological degradation reactions occurring in three distinct patterns involving fringe, pseudo‐core and core processes. The capabilities of the numerical model are demonstrated using two hypothetical examples: a batch problem and a simplified two‐dimensional reactive transport problem. The model is then applied to an unconfined aquifer underlying a leaking landfill located near the city of Turin, in Piedmont (Italy). At this site, high organic load from the landfill leachate activates different biogeochemical processes, including aerobic degradation, denitrification, manganese reduction, iron reduction, sulfate reduction and methanogenesis. The model was able to describe and quantify these complex biogeochemical processes. The proposed model offers a rational framework for simulating coupled reactive transport processes occurring beneath a landfill site. Copyright © 2008 John Wiley & Sons, Ltd.     

VS2DRTI: Simulating Heat and Reactive Solute Transport in Variably Saturated Porous Media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software.     

Simulated effect of a forest road on near‐surface hydrologic response and slope stability

Forest management practices often result in significant changes to hydrologic and geomorphic responses at or near the earth's surface. A well‐known, but not fully tested, hypothesis in hillslope hydrology[sol ]geomorphology is that a near‐surface permeability contrast, caused by the surface compaction associated with forest roads, can result in diverted subsurface flow paths that produce increased up‐slope pore pressures and slope failure. The forest road focused on in this study is located in a steep forested, zero‐order catchment within the H. J. Andrews Experimental Forest (Oregon). A three‐phase modelling effort was employed to test the aforementioned hypothesis: (i) two‐dimensional (vertical slice), steady‐state, heterogeneous, saturated subsurface flow simulations at the watershed scale for establishing the boundary conditions for the catchment‐scale boundary‐value problem in (ii); (ii) two‐dimensional (vertical slice), transient, heterogeneous, variably saturated subsurface flow simulations at the catchment scale for estimating near‐surface hydrologic response and pore pressure distributions; and (iii) slope stability analyses, using the infinite slope approach, driven by the pore pressure distributions simulated in (ii), for assessing the impact of the forest road. Both observed and hypothetical rainfall events are used to drive the catchment‐scale simulations. The results reported here support the hypothesis that a forest road can have an effect on slope stability. The permeability contrast associated with the forest road in this study led to a simulated altering of slope‐parallel subsurface flow with increased pore pressures up‐slope of the road and, for a large rainfall event, a slope failure prediction. Copyright © 2005 John Wiley & Sons, Ltd.     

A two‐dimensional model of hydraulic performance of stormwater infiltration systems

Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd.     

Aquifer Vulnerability Assessment Based on Sequence Stratigraphic and 39Ar Transport Modeling

A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of ³⁹Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.     

Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

In this work we study mixed finite element approximations of Richards’ equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61–78]), the Raviart–Thomas mixed finite element method (RT₀) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565–581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373–394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869–1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148–167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701–724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi-Douglas-Marini element (BDM₁) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM₁ approach to the RT₀ one is observed, which however is less significant for the accompanying solute transport.     

Investigation of the thermal regime and subsurface properties of a tidally affected,variably saturated streambed

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.     

A coupled model for simulating surface water and groundwater interactions in coastal wetlands

Coastal wetlands are characterized by strong, dynamic interactions between surface water and groundwater. This paper presents a coupled model that simulates interacting surface water and groundwater flow and solute transport processes in these wetlands. The coupled model is based on two existing (sub) models for surface water and groundwater, respectively: ELCIRC (a three‐dimensional (3‐D) finite‐volume/finite‐difference model for simulating shallow water flow and solute transport in rivers, estuaries and coastal seas) and SUTRA (a 3‐D finite‐element/finite‐difference model for simulating variably saturated, variable‐density fluid flow and solute transport in porous media). Both submodels, using compatible unstructured meshes, are coupled spatially at the common interface between the surface water and groundwater bodies. The surface water level and solute concentrations computed by the ELCIRC model are used to determine the boundary conditions of the SUTRA‐based groundwater model at the interface. In turn, the groundwater model provides water and solute fluxes as inputs for the continuity equations of surface water flow and solute transport to account for the mass exchange across the interface. Additionally, flux from the seepage face was routed instantaneously to the nearest surface water cell according to the local sediment surface slope. With an external coupling approach, these two submodels run in parallel using time steps of different sizes. The time step (Δt_g) for the groundwater model is set to be larger than that (Δt_s) used by the surface water model for computational efficiency: Δt_g = M × Δt_s where M is an integer greater than 1. Data exchange takes place between the two submodels through a common database at synchronized times (e.g. end of each Δt_g). The coupled model was validated against two previously reported experiments on surface water and groundwater interactions in coastal lagoons. The results suggest that the model represents well the interacting surface water and groundwater flow and solute transport processes in the lagoons. Copyright © 2011 John Wiley & Sons, Ltd.     

Transport and fate of nonaqueous phase liquid (NAPL) in variably saturated porous media with evolving scales of heterogeneity

 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.     

Evaluation of model complexity and space–time resolution on the prediction of long‐term soil salinity dynamics,western San Joaquin Valley,California

The numerical simulation of long‐term large‐scale (field to regional) variably saturated subsurface flow and transport remains a computational challenge, even with today's computing power. Therefore, it is appropriate to develop and use simplified models that focus on the main processes operating at the pertinent time and space scales, as long as the error introduced by the simpler model is small relative to the uncertainties associated with the spatial and temporal variation of boundary conditions and parameter values. This study investigates the effects of various model simplifications on the prediction of long‐term soil salinity and salt transport in irrigated soils. Average root‐zone salinity and cumulative annual drainage salt load were predicted for a 10‐year period using a one‐dimensional numerical flow and transport model (i.e. UNSATCHEM) that accounts for solute advection, dispersion and diffusion, and complex salt chemistry. The model uses daily values for rainfall, irrigation, and potential evapotranspiration rates. Model simulations consist of benchmark scenarios for different hypothetical cases that include shallow and deep water tables, different leaching fractions and soil gypsum content, and shallow groundwater salinity, with and without soil chemical reactions. These hypothetical benchmark simulations are compared with the results of various model simplifications that considered (i) annual average boundary conditions, (ii) coarser spatial discretization, and (iii) reducing the complexity of the salt‐soil reaction system. Based on the 10‐year simulation results, we conclude that salt transport modelling does not require daily boundary conditions, a fine spatial resolution, or complex salt chemistry. Instead, if the focus is on long‐term salinity, then a simplified modelling approach can be used, using annually averaged boundary conditions, a coarse spatial discretization, and inclusion of soil chemistry that only accounts for cation exchange and gypsum dissolution–precipitation. We also demonstrate that prediction errors due to these model simplifications may be small, when compared with effects of parameter uncertainty on model predictions. The proposed model simplifications lead to larger time steps and reduced computer simulation times by a factor of 1000. Copyright © 2006 John Wiley & Sons, Ltd.     

Phast4Windows: A 3D Graphical User Interface for the Reactive‐Transport Simulator PHAST

Phast4Windows is a Windows® program for developing and running groundwater‐flow and reactive‐transport models with the PHAST simulator. This graphical user interface allows definition of grid‐independent spatial distributions of model properties—the porous media properties, the initial head and chemistry conditions, boundary conditions, and locations of wells, rivers, drains, and accounting zones—and other parameters necessary for a simulation. Spatial data can be defined without reference to a grid by drawing, by point‐by‐point definitions, or by importing files, including ArcInfo® shape and raster files. All definitions can be inspected, edited, deleted, moved, copied, and switched from hidden to visible through the data tree of the interface. Model features are visualized in the main panel of the interface, so that it is possible to zoom, pan, and rotate features in three dimensions (3D). PHAST simulates single phase, constant density, saturated groundwater flow under confined or unconfined conditions. Reactions among multiple solutes include mineral equilibria, cation exchange, surface complexation, solid solutions, and general kinetic reactions. The interface can be used to develop and run simple or complex models, and is ideal for use in the classroom, for analysis of laboratory column experiments, and for development of field‐scale simulations of geochemical processes and contaminant transport.     

Complexities in Hindcasting Models—When Should We Say Enough Is Enough?

Groundwater models are routinely used in hindcasting applications to predict the past concentration levels in contaminated aquifers. These predictions are used in risk assessment and epidemiological studies, which are often completed either for resolving a court case or for developing a public‐policy solution. Hindcast groundwater modeling studies utilize a variety of computer tools with complexity levels ranging from simple analytical models to detailed three‐dimensional, multiphase, multispecies, reactive transport models. The aim of this study is to explore the value of using complex reactive transport models in hindcasting studies that have limited historic data. I review a chlorinated solvent exposure problem that occurred at a U.S. Marine Corp Base in Camp Lejeune, North Carolina and use it as an example to discuss the limits of hindcasting modeling exercises. The lessons learned from the study are used to reflect upon the following questions related to model complexity: How should we decide how much is enough? Who should decide when enough is enough?     

Modelling long‐term contaminant migration in a catchment at fine spatial and temporal scales using the UP system

A method has been developed to simulate the long‐term migration of radionuclides in the near‐surface of a river catchment, following their release from a deep underground repository for radioactive waste. Previous (30‐year) simulations, conducted using the SHETRAN physically based modelling system, showed that long‐term (many decades) simulations are required to allow the system to reach steady state. Physically based, distributed models, such as SHETRAN, tend to be too computationally expensive for this task. Traditional lumped catchment‐scale models, on the other hand, do not give sufficiently detailed spatially distributed results. An intermediate approach to modelling has therefore been developed which allows flow and transport processes to be simulated with the spatial resolution normally associated with distributed models, whilst being computationally efficient.The approach involves constructing a lumped model in which the catchment is represented by a number of conceptual water storage compartments. The flow rates to and from these compartments are prescribed by functions that summarize the results from physically based distributed models run for a range of characteristic flow regimes. The physically based models used were, SHETRAN for the subsurface compartments, a particle tracking model for overland flow and an analytical model for channel routing. One important advantage of the method used in constructing the lumped model is that it makes down scaling possible, in the sense that fine‐scale information on the distributed hydrological regime, as simulated by the physically based distributed models, can be inferred from the variables in the lumped model that describe the hydrology at the catchment scale. A 250‐year flow simulation has been run and the down scaling process used to infer a 250‐year time‐series of three‐dimensional velocity fields for the subsurface of the catchment. This series was then used to drive a particle tracking simulation of contaminant migration. The concentration and spatial distribution of contaminants simulated by this model for the first 30 years were in close agreement with SHETRAN results. The remaining 220 years highlighted the fact that some of the most important transport pathways to the surface carry contaminants only very slowly so both the magnitude and spatial distribution of concentration in surface soils are not apparent over the shorter SHETRAN simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

How Does Subsurface Characterization Affect Simulations of Hyporheic Exchange?

We investigated the role of increasingly well‐constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time‐lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10^?4 vs. 10^?7 m²/s, respectively). Median residence times for uncalibrated and calibrated models are 10³ and 10⁰ h, respectively. Increasingly well‐resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.