Geology Reports⌕ Search

Geology topics

C.D. Langevin

Publications and source records attributed to C.D. Langevin.

13 recordsLinked to original sources
USGS2013

Seawater circulation in sediments driven by interactions between seabed topography and fluid density

Measurements of submarine groundwater discharge (SGD) in coastal areas often show that the saltwater discharge component is substantially greater than the freshwater discharge. Several mechanisms have been proposed to explain these high saltwater discharge values, including saltwater circulation driven by wave and tidal pumping, wave and tidal setup in intertidal areas, currents over bedforms, and density gradients resulting from mixing along the freshwater-saltwater interface. In this study, a new mechanism for saltwater circulation and discharge is proposed and evaluated. The process results from interaction between bedform topography and buoyancy forces, even without flow or current over the bedform. In this mechanism, an inverted salinity (and density) profile in the presence of both a bedform on the seafloor and an upward flow of fresher groundwater from depth induces a downward flow of saline pore water under the troughs and upward flow under the adjacent crest of the bedform. The magnitude and occurrence of the mechanism were tested using numerical methods. The results indicate that this mechanism could drive seawater circulation under a limited range of conditions and contribute 20%–30% of local SGD when and where the process is operative. Bedform shape, hydraulic conductivity, hydraulic head, and salinity at depth in the porous media, aquifer thickness, effective porosity, and hydrodynamic dispersion are among the factors that control the occurrence and magnitude of the circulation of seawater by this mechanism.

Water Resources Research↗
USGS2012

Numerical simulation of saltwater intrusion in response to sea-level rise

A two dimensional numerical model of variable-density groundwater flow and dispersive solute transport was used to predict the extent, rate, and lag time of saltwater intrusion in response to various sea-level rise scenarios. Three simulations were performed with varying rates of sea-level rise. For the first simulation, sea-level rise was specified at a rate of 0.9 mm/yr, which is the slowest rate of sea-level rise estimated by the Intergovernmental Panel on Climate Change (IPCC). After 100 years, the 250 mg/L chloride isochlor moved inland by about 40 m, and required an additional 8 years for the system to reach equilibrium. For the next simulation, sea-level rise was specified at 4.8 mm/yr, which is the central value of the IPCC estimate. For this moderate rate of sea-level rise, the 250 mg/L isochlor moved inland by about 740 m after 100 years, and required an additional 10 years for the system to reach equilibrium. For the fastest rate of sea level rise estimated by IPCC (8.8 mm/yr), the 250 mg/L isochlor moved inland by about 1800 m after 100 years, and required more than 50 years to reach equilibrium. Copyright ASCE 2005.

Conference Paper↗
USGS2011

Use of upscaled elevation and surface roughness data in two-dimensional surface water models

In this paper, we present an approach that uses a combination of cell-block- and cell-face-averaging of high-resolution cell elevation and roughness data to upscale hydraulic parameters and accurately simulate surface water flow in relatively low-resolution numerical models. The method developed allows channelized features that preferentially connect large-scale grid cells at cell interfaces to be represented in models where these features are significantly smaller than the selected grid size. The developed upscaling approach has been implemented in a two-dimensional finite difference model that solves a diffusive wave approximation of the depth-integrated shallow surface water equations using preconditioned Newton–Krylov methods. Computational results are presented to show the effectiveness of the mixed cell-block and cell-face averaging upscaling approach in maintaining model accuracy, reducing model run-times, and how decreased grid resolution affects errors. Application examples demonstrate that sub-grid roughness coefficient variations have a larger effect on simulated error than sub-grid elevation variations.

Advances in Water Resources↗
USGS2010

Quantifying data worth toward reducing predictive uncertainty

The present study demonstrates a methodology for optimization of environmental data acquisition. Based on the premise that the worth of data increases in proportion to its ability to reduce the uncertainty of key model predictions, the methodology can be used to compare the worth of different data types, gathered at different locations within study areas of arbitrary complexity. The method is applied to a hypothetical nonlinear, variable density numerical model of salt and heat transport. The relative utilities of temperature and concentration measurements at different locations within the model domain are assessed in terms of their ability to reduce the uncertainty associated with predictions of movement of the salt water interface in response to a decrease in fresh water recharge. In order to test the sensitivity of the method to nonlinear model behavior, analyses were repeated for multiple realizations of system properties. Rankings of observation worth were similar for all realizations, indicating robust performance of the methodology when employed in conjunction with a highly nonlinear model. The analysis showed that while concentration and temperature measurements can both aid in the prediction of interface movement, concentration measurements, especially when taken in proximity to the interface at locations where the interface is expected to move, are of greater worth than temperature measurements. Nevertheless, it was also demonstrated that pairs of temperature measurements, taken in strategic locations with respect to the interface, can also lead to more precise predictions of interface movement. Journal compilation ?? 2010 National Ground Water Association.

Ground Water↗
USGS2008

Modeling axisymmetric flow and transport

Unmodified versions of common computer programs such as MODFLOW, MT3DMS, and SEAWAT that use Cartesian geometry can accurately simulate axially symmetric ground water flow and solute transport. Axisymmetric flow and transport are simulated by adjusting several input parameters to account for the increase in flow area with radial distance from the injection or extraction well. Logarithmic weighting of interblock transmissivity, a standard option in MODFLOW, can be used for axisymmetric models to represent the linear change in hydraulic conductance within a single finite-difference cell. Results from three test problems (ground water extraction, an aquifer push-pull test, and upconing of saline water into an extraction well) show good agreement with analytical solutions or with results from other numerical models designed specifically to simulate the axisymmetric geometry. Axisymmetric models are not commonly used but can offer an efficient alternative to full three-dimensional models, provided the assumption of axial symmetry can be justified. For the upconing problem, the axisymmetric model was more than 1000 times faster than an equivalent three-dimensional model. Computational gains with the axisymmetric models may be useful for quickly determining appropriate levels of grid resolution for three-dimensional models and for estimating aquifer parameters from field tests.

Ground Water↗
USGS2008

Utilizing spectral analysis of coastal discharge computed by a numerical model to determine boundary influence

In the present study, a spectral analysis was applied to field data and a numerical model of southeastern Everglades and northeastern Florida Bay that involved computing and comparing the power spectrum of simulated and measured flows at the primary coastal outflow creek. Four dominant power frequencies, corresponding to the S1, S2, M2, and 01 tidal periods, were apparent in the measured outflows. The model seemed to reproduce the magnitudes of the S1 and S2 components better than those of the M2 and 01 components. To determine the cause of the relatively poor representation of the M2 and 01 components, we created a steady-base version of the model by setting the time-varying forcing functions - rainfall, evapotranspiration, wind, and inland and tidal boundary conditions - to averaged values. The steady-base model was then modified to produce multiple simulations with only one time-varying forcing function for each model run. These experimental simulations approximated the individual effects of each forcing function on the system. The spectral analysis of the experimental simulations indicated that temporal fluctuations in rainfall, evapotranspiration, and inland water level and discharge boundaries have negligible effects on coastal creek flow fluctuations with periods of less than 48 hours. The tidal boundary seems to be the only forcing function inducing the M2 and 01 frequency flow fluctuations in the creek. An analytical formulation was developed, relating the errors induced by the tidal water-level gauge resolution to the errors in the simulated discharge fluctuations at the coastal creek. This formulation yielded a discharge-fluctuation error similar in magnitude to the errors observed when comparing the spectrum of the simulated and measured discharge. The dominant source of error in the simulation of discharge fluctuation magnitude is most likely the resolution of the water-level gauges used to create the model boundary.

Journal of Coastal Research↗
USGS2007

Simulation of submarine groundwater discharge salinity and temperature variations: Implications for remote detection

A hydrological analysis using a numerical simulation was done to identify the transient response of the salinity and temperature of submarine groundwater discharge (SGD) and utilize the results to guide data collection. Results indicate that the amount of SGD fluctuates depending on the ocean stage and geology, with the greatest amount of SGD delivered at low tide when the aquifer is in direct hydraulic contact with the ocean. The salinity of SGD remains lower than the ocean throughout the year; however, the salinity difference between the aquifer and ocean is inversely proportional to the ocean stage. The temperature difference between the ocean and SGD fluctuates seasonally, with the greatest temperature differences occurring in summer and winter. The outcome of this research reveals that numerical modelling could potentially be used to guide data collection including aerial surveys using electromagnetic (EM) resistivity and thermal imagery.

Conference Paper↗
USGS2006

Addition of simultaneous heat and solute transport and variable fluid viscosity to SEAWAT

SEAWAT is a finite-difference computer code designed to simulate coupled variable-density ground water flow and solute transport. This paper describes a new version of SEAWAT that adds the ability to simultaneously model energy and solute transport. This is necessary for simulating the transport of heat and salinity in coastal aquifers for example. This work extends the equation of state for fluid density to vary as a function of temperature and/or solute concentration. The program has also been modified to represent the effects of variable fluid viscosity as a function of temperature and/or concentration. The viscosity mechanism is verified against an analytical solution, and a test of temperature-dependent viscosity is provided. Finally, the classic Henry-Hilleke problem is solved with the new code. ?? 2006 Elsevier Ltd. All rights reserved.

Computers & Geosciences↗
USGS2006

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software↗
USGS2005

Stochastic Ground Water Flow Simulation with a Fracture Zone Continuum Model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.

Ground Water↗
USGS2005

Simulation of Submarine Ground Water Discharge to a Marine Estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only ???10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.

Ground Water↗
USGS2004

The rotating movement of three immiscible fluids - A benchmark problem

A benchmark problem involving the rotating movement of three immiscible fluids is proposed for verifying the density-dependent flow component of groundwater flow codes. The problem consists of a two-dimensional strip in the vertical plane filled with three fluids of different densities separated by interfaces. Initially, the interfaces between the fluids make a 45??angle with the horizontal. Over time, the fluids rotate to the stable position whereby the interfaces are horizontal; all flow is caused by density differences. Two cases of the problem are presented, one resulting in a symmetric flow field and one resulting in an asymmetric flow field. An exact analytical solution for the initial flow field is presented by application of the vortex theory and complex variables. Numerical results are obtained using three variable-density groundwater flow codes (SWI, MOCDENS3D, and SEAWAT). Initial horizontal velocities of the interfaces, as simulated by the three codes, compare well with the exact solution. The three codes are used to simulate the positions of the interfaces at two times; the three codes produce nearly identical results. The agreement between the results is evidence that the specific rotational behavior predicted by the models is correct. It also shows that the proposed problem may be used to benchmark variable-density codes. It is concluded that the three models can be used to model accurately the movement of interfaces between immiscible fluids, and have little or no numerical dispersion. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗