This paper discusses a method for overcoming the problem of weak sinks representing wells that result from spatial discretization effects when using MODPATH, the particle-tracking postprocessor for the ground water flow model MODFLOW. Weak sink cells are model cells that represent a well that does not discharge at a sufficiently large rate to capture all of the flow entering the cell; therefore, flowpaths within these cells cannot be uniquely defined because it is impossible to know whether a given water particle discharges to the well or passes through the cell. Creating a submodel of the well cell by using the nested rediscretization method can eliminate this ambiguity by converting the weak sink cell into a strong sink cell. The method is designed to be run manually for each well and for steady-state conditions. Other advantages, disadvantages, technical considerations, and limitations of the method are presented. Software created for the method consists of five Fortran programs that are operated using a set of instructions. A practical application of the method is presented by using an example wellhead-protection problem that demonstrates that nested rediscretization can provide more accurate particle-tracking results than those obtained by using a coarsely discretized model alone.

The Regional Ocean Modeling System (ROMS) has been coupled with the Water Quality Analysis Simulation Program (WASP) to be used in a comprehensive analysis of water quality in Barnegat Bay, New Jersey. The coupler can spatially aggregate hydrodynamic information in ROMS cells into larger WASP segments. It can also be used to resample ROMS output at a finer temporal scale to meet WASP time-stepping requirements. The coupler aggregates flow components, temperature, and salinity in ROMS output for input to WASP via a hydrodynamic linkage file. The coupler was tested initially with idealized cases designed to verify the water mass balance and conservation of constituent mass using oneto-one and one-to-many connectivity options between segments. A realistic example from the Toms River embayment, a subdomain of Barnegat Bay, was used to demonstrate the functionality of the coupling. A WASP eutrophication model accounting for dissolved oxygen (DO), nitrogen, and constant phytoplankton concentrations was applied to explore the distribution and trends in DO and nitrogen in the embayment for the period of July-August 2012. Results of DO modeling indicate satisfactory agreement with measurements collected at in-bay stations and also indicate that this coupled approach, despite substantial differences in spatiotemporal discretization between the models, provides adequate predictive capabilities. ADDITIONAL INDEX WORDS:ROMS, WASP, hydrodynamic linkage, coupler.

No abstract

No abstract

No abstract