Travel Time Distributions Derived from Particle Tracking in Models Containing Weak Sinks

A travel time distribution based on a particle-tracking analysis in a ground water model containing weak sinks is often uncertain because whether a particle is discharged or allowed to pass through a weak sink is unresolved by particle-tracking theory. We present a probability-based method to derive an objective travel time distribution in models containing weak sinks. The method discharges a fraction of the particle at the weak sink and allows the remaining fraction to pass through the weak sink. The weight of the discharged fraction depends on the ratio of the sink flux to the influx into the weak sink cell. We tested this approach on a coarse (100 × 100 m) and a fine (25 × 25 m) horizontal resolution regional scale ground water model (34.5 × 24 km). We compared the travel time distributions in a small subcatchment derived from particle-tracking analysis with one derived from a transport model. We found that the particle-tracking analysis with the coarse model underestimated the travel time distribution of the catchment compared to the transport solution or a particle-tracking analysis with the fine model. The underestimation of travel times with the coarse model was a result of a large area covered by sink cells in this model and the more accurate flow patterns simulated by the fine model. The probability-based method presented here compares favorably with a solute transport solution and provides an accurate travel time distribution when used with a fine-resolution ground water model.     

Impact of Sea-Level Rise on Sea Water Intrusion in Coastal Aquifers

Despite its purported importance, previous studies of the influence of sea-level rise on coastal aquifers have focused on specific sites, and a generalized systematic analysis of the general case of the sea water intrusion response to sea-level rise has not been reported. In this study, a simple conceptual framework is used to provide a first-order assessment of sea water intrusion changes in coastal unconfined aquifers in response to sea-level rise. Two conceptual models are tested: (1) flux-controlled systems, in which ground water discharge to the sea is persistent despite changes in sea level, and (2) head-controlled systems, whereby ground water abstractions or surface features maintain the head condition in the aquifer despite sea-level changes. The conceptualization assumes steady-state conditions, a sharp interface sea water-fresh water transition zone, homogeneous and isotropic aquifer properties, and constant recharge. In the case of constant flux conditions, the upper limit for sea water intrusion due to sea-level rise (up to 1.5 m is tested) is no greater than 50 m for typical values of recharge, hydraulic conductivity, and aquifer depth. This is in striking contrast to the constant head cases, in which the magnitude of salt water toe migration is on the order of hundreds of meters to several kilometers for the same sea-level rise. This study has highlighted the importance of inland boundary conditions on the sea-level rise impact. It identifies combinations of hydrogeologic parameters that control whether large or small salt water toe migration will occur for any given change in a hydrogeologic variable.     

Ground Water Development—The Time to Full Capture Problem

Ground water systems can be categorized with respect to quantity into two groups: (1) those that will ultimately reach a new equilibrium state where pumping can be continued indefinitely and (2) those in which the stress is so large that a new equilibrium is impossible; hence, the system has a finite life. Large ground water systems, where a new equilibrium can be reached and in which the pumping is a long distance from boundaries where capture can occur, take long times to reach a new equilibrium. Some systems are so large that the new equilibrium will take a millennium or more to reach a new steady-state condition. These large systems pose a challenge to the water manager, especially when the water manager is committed to attempting to reach a new equilibrium state in which water levels will stabilize and the system can be maintained indefinitely.     

Evaluation of Noble Gas Recharge Temperatures in a Shallow Unconfined Aquifer

Water table temperatures inferred from dissolved noble gas concentrations (noble gas temperatures, NGT) are useful as a quantitative proxy for air temperature change since the last glacial maximum. Despite their importance in paleoclimate research, few studies have investigated the relationship between NGT and actual recharge temperatures in field settings. This study presents dissolved noble gas data from a shallow unconfined aquifer heavily impacted by agriculture. Considering samples unaffected by degassing, NGT calculated from common physically based interpretive gas dissolution models that correct measured noble gas concentrations for "excess air" agreed with measured water table temperatures (WTT). The ability to fit data to multiple interpretive models indicates that model goodness-of-fit does not necessarily mean that the model reflects actual gas dissolution processes. Although NGT are useful in that they reflect WTT, caution is recommended when using these interpretive models. There was no measurable difference in excess air characteristics (amount and degree of fractionation) between two recharge regimes studied (higher flux recharge primarily during spring and summer vs. continuous, low flux recharge). Approximately 20% of samples had dissolved gas concentrations below equilibrium concentration with respect to atmospheric pressure, indicating degassing. Geochemical and dissolved gas data indicate that saturated zone denitrification caused degassing by gas stripping. Modeling indicates that minor degassing (<10% ΔNe) may cause underestimation of ground water recharge temperature by up to 2°C. Such errors are problematic because degassing may not be apparent and degassed samples may be fit by a model with a high degree of certainty.     

Global Optimal Design of Ground Water Monitoring Network Using Embedded Kriging

We present a methodology for global optimal design of ground water quality monitoring networks using a linear mixed-integer formulation. The proposed methodology incorporates ordinary kriging (OK) within the decision model formulation for spatial estimation of contaminant concentration values. Different monitoring network design models incorporating concentration estimation error, variance estimation error, mass estimation error, error in locating plume centroid, and spatial coverage of the designed network are developed. A big-M technique is used for reformulating the monitoring network design model to a linear decision model while incorporating different objectives and OK equations. Global optimality of the solutions obtained for the monitoring network design can be ensured due to the linear mixed-integer programming formulations proposed. Performances of the proposed models are evaluated for both field and hypothetical illustrative systems. Evaluation results indicate that the proposed methodology performs satisfactorily. These performance evaluation results demonstrate the potential applicability of the proposed methodology for optimal ground water contaminant monitoring network design.     

Submarine Ground Water Discharge Driven by Tidal Pumping in a Heterogeneous Aquifer

Submarine ground water discharge (SGD) is now recognized as an important water pathway between land and sea. It is difficult to quantitatively predict SGD owing to its significant spatial and temporal variability. This study focuses on quantitative estimation of SGD caused by tidally induced sea water recirculation and a terrestrial hydraulic gradient. A two-dimensional hydrogeological model was developed to simulate SGD from a coastal unconfined aquifer in the northeastern Gulf of Mexico, where previous SGD studies were performed. A density-variable numerical code, SEAWAT2000, was applied to simulate SGD. To accurately predict discharge, various influencing factors such as heterogeneity in conductivity, uncertain boundary conditions, and tidal pumping were systematically assessed. The tidally influenced sea water recirculation zone and the fresh water–salt water mixing zone under various tidal patterns, tidal ranges, and water table heights were also investigated. The model was calibrated and validated from long-term, intensive measurements at the study site. The percentage of fresh SGD relative to total SGD ranged from 4% to 50% under normal conditions. Based on simulations of two field measurements in summer and spring, respectively, the fresh water ratios were 9% and 15%, respectively. These results support the hypothesis that the SGD induced by tidally driven sea water recirculation is much larger than terrestrial fresh ground water discharge at this site. The estimates of total and fresh SGD are at the low and high ends, respectively, of the estimation ranges obtained from geochemical tracers (e.g., ²²²Rn).     

Influences of Aquifer Properties on Flow Dimensions in Dolomites

The paper focuses on analyses and correlations of flow dimensions in different dolomite aquifers in Slovenia. Flow dimensions are obtained through the reinterpretation of 72 pumping tests with the generalized radial flow model, based on the fractional flow dimension. The average value of flow dimensions is 2.16 for all dolomites. A study of flow dimensions in individual aquifers categorized according to their lithological properties shows that higher dimensions occur in massive late-diagenetic Cordevolian and Anisian dolomites compared with bedded Main, Bača, and especially Lower Triassic dolomites, which contain a greater proportion of noncarbonate minerals. Partially penetrating wells have higher flow dimensions than fully penetrating wells. Flow dimensions are poorly correlated with hydraulic conductivities of fractures. When comparing the quantities of major dissolved minerals, obtained by hydrogeochemical inverse modeling, with the values of flow dimensions, the Cordevolian and Anisian dolomites are found to exhibit the highest values of both dissolved dolomite and flow dimensions, indicating that greater dissolution occurs at higher flow dimensions. For other aquifers, data points are more scattered and the correlation is mostly poor. When compared with three-dimensional fractal dimensions of fracture networks, there is no correlation with flow dimensions. However, almost all the values of flow dimensions are lower than the corresponding fractal dimensions in dolomites (average D = 2.77), possibly indicating the channeling of flow within the available space of the fracture networks, consequently reducing the flow dimensions.     

Determining Hydraulic Conductivity Using Pumping Data from Low-Flow Sampling

Hydraulic conductivity values computed using the steady-state discharge and drawdown attained while low-flow sampling were evaluated to determine if they were equivalent to those determined from slug testing. Based on testing 12 wells, it was found that the results were statistically equivalent. Conductivity values computed using low-flow sampling parameters were also evaluated as to their reproducibility in actual practice by analyzing consultant data for three wells sampled over three quarterly monitoring periods by four field technicians. The results were found to be reproducible within about a factor of 2 or better. Since the method is based on only one pair of parameters, diligence is required in attaining steady state and in accurately measuring the flow rate and drawdown. Conductivity values computed using this approach can enhance the use of low-flow data gathered in water quality sampling, avoid the need for slug testing in a subsequent phase of investigation, and help reduce the cost of characterizing sites when multilevel samplers are used. Given the practical range of discharge in low-flow sampling, the method was found to be applicable at conductivity values somewhat greater than 10⁻⁶ cm/s. Given the typical accuracy of water level meters and pressure transducers and a maximum discharge of 1 L/min, as mandated by regulatory guidance, the method has a calculated upper conductivity limit in the range of 10⁻³ to 10⁻² cm/s.     

Influence of an oxygen minimum zone and macroalgal enrichment on benthic megafaunal community composition in a NE Pacific submarine canyon

Megafaunal diversity in the deep sea shows a parabolic pattern with depth. It can be affected by factors such as low oxygen concentration, which suppresses diversity, or the presence of submarine canyons, which enhances it. Barkley Canyon, located off the west coast of British Columbia, Canada, is a submarine canyon that extends from the continental margin (200 m) into the deep ocean (2,000 m). This canyon receives drift kelp from shoreline kelp forests and contains an oxygen minimum zone (OMZ) at 500 to 1,500 m depth. Our study investigated the abundance and diversity of epibenthic megafauna over a range of depths (200–2,000 m) and oxygen concentrations (0.5–5.0 ml/L) within Barkley Canyon, as well as changes in abundance near detrital kelp. Video was collected using the remotely operated vehicle ROPOS along seven 1‐km cross‐canyon (i.e., across the axis of the canyon) transects and three 40‐m perpendicular cross‐transects over kelp. Taxonomic groups were associated with depth, temperature, and the presence of pebbles. The OMZ restricted pennatulids, and edge effects along OMZ boundaries were observed for ophiuroids. The geomorphology of the sea floor affected the distribution of taxa across the canyon, with Porifera mainly found along the walls and Echinoidea within the canyon axis. Expected richness exhibited a bimodal pattern, peaking at 300 and 2,000 m, possibly due to the combined effect of the OMZ and the submarine canyon. Echinoidea aggregated near drift kelp at 200 and 300 m. We found that faunal communities in Barkley Canyon were influenced by several confounded factors including depth, oxygen and substrate. Understanding faunal patterns is paramount with increased exploitation and a changing climate.