Xiaolong Geng

Numerical study of solute transport in shallow beach aquifers subjected to waves and tides

A numerical study was conducted to investigate the fate of solute in a laboratory beach in response to waves and tides. A new temporal upscaling approach labeled “net inflow” was introduced to address impacts of waves on solute transport within beaches. Numerical simulations using a computational fluid dynamic model were used as boundary conditions for the two-dimensional variably saturated flow and solute transport model MARUN. The modeling approach was validated against experimental data of solute transport due to waves and tides. Exchange fluxes across the beach face and subsurface solute transport (e.g., trajectory, movement speed, and residence time) were quantified. Simulation results revealed that waves increased the exchange fluxes, and engendered a wider exchange flux zone along the beach surface. Compared to tide-only forcing, waves superimposed on tide caused the plume to be deeper into the beach, and to migrate more seaward. The infiltration into the beach was found to be directly proportional to the general hydraulic gradient in the beach and inversely proportional to the matrix retention (or capillary) capacity. The simulations showed that a higher inland water table would attenuate wave-caused seawater infiltration, which might impact beach geochemical processes (e.g., nutrient recycle and redox condition), especially at low tide zone. The concept of biochemical residence time maps (BRTM) was introduced to account for the net effect of limiting concentration of chemicals on biochemical reactions. It was found that waves shifted the BRTMs downward and seaward in the beach, and subsequently they engendered different biochemical conditions within the beach.

Numerical study of wave effects on groundwater flow and solute transport in a laboratory beach

A numerical study was undertaken to investigate the effects of waves on groundwater flow and associated inland-released solute transport based on tracer experiments in a laboratory beach. The MARUN model was used to simulate the density-dependent groundwater flow and subsurface solute transport in the saturated and unsaturated regions of the beach subjected to waves. The Computational Fluid Dynamics (CFD) software, Fluent, was used to simulate waves, which were the seaward boundary condition for MARUN. A no-wave case was also simulated for comparison. Simulation results matched the observed water table and concentration at numerous locations. The results revealed that waves generated seawater-groundwater circulations in the swash and surf zones of the beach, which induced a large seawater-groundwater exchange across the beach face. In comparison to the no-wave case, waves significantly increased the residence time and spreading of inland-applied solutes in the beach. Waves also altered solute pathways and shifted the solute discharge zone further seaward. Residence Time Maps (RTM) revealed that the wave-induced residence time of the inland-applied solutes was largest near the solute exit zone to the sea. Sensitivity analyses suggested that the change in the permeability in the beach altered solute transport properties in a nonlinear way. Due to the slow movement of solutes in the unsaturated zone, the mass of the solute in the unsaturated zone, which reached up to 10% of the total mass in some cases, constituted a continuous slow release of solutes to the saturated zone of the beach. This means of control was not addressed in prior studies.

Simulation of the Landfall of the Deepwater Horizon Oil on the Shorelines of the Gulf of Mexico

We conducted simulations of oil transport from the footprint of the Macondo Well on the water surface throughout the Gulf of Mexico, including deposition on the shorelines. We used the U.S. National Oceanic Atmospheric Administration (NOAA) model General NOAA Operational Modeling Environment (GNOME) and the same parameter values and input adopted by NOAA following the Deepwater Horizon (DWH) blowout. We found that the disappearance rate of oil off the water surface was most likely around 20% per day based on satellite-based observations of the disappearance rate of oil detected on the sea surface after the DWH wellhead was capped. The simulations and oil mass estimates suggest that the mass of oil that reached the shorelines was between 10,000 and 30,000 tons, with an expected value of 22,000 tons. More than 90% of the oil deposition occurred on the Louisiana shorelines, and it occurred in two batches. Simulations revealed that capping the well after 2 weeks would have resulted in only 30% of the total oil depositing on the shorelines, while capping after 3 weeks would have resulted in 60% deposition. Additional delay in capping after 3 weeks would have averted little additional shoreline oiling over the ensuing 4 weeks.

Biodegradation of subsurface oil in a tidally influenced sand beach: Impact of hydraulics and interaction with pore water chemistry

The aerobic biodegradation of oil in tidally influenced beaches was investigated numerically in this work using realistic beach and tide conditions. A numerical model BIOMARUN, coupling a multiple-Monod kinetic model BIOB to a density-dependent variably saturated groundwater flow model 2-D MARUN, was used to simulate the biodegradation of low-solubility hydrocarbon and transport processes of associated solute species (i.e., oxygen and nitrogen) in a tidally influenced beach environment. It was found that different limiting factors affect different portions of the beach. In the upper intertidal zone, where the inland incoming nutrient concentration was large (1.2 mg N/L), oil biodegradation occurred deeper in the beach (i.e., 0.3 m below the surface). In the midintertidal zone, a reversal was noted where the biodegradation was fast at shallow locations (i.e., 0.1 m below the surface), and it was due to the decrease of oxygen with depth due to consumption, which made oxygen the limiting factor for biodegradation. Oxygen concentration in the midintertidal zone exhibited two peaks as a function of time. One peak was associated with the high tide, when dissolved oxygen laden seawater filled the beach and a second oxygen peak was observed during low tides, and it was due to pore oxygen replenishment from the atmosphere. The effect of the capillary fringe (CF) height was investigated, and it was found that there is an optimal CF for the maximum biodegradation of oil in the beach. Too large a CF (i.e., very fine material) would attenuate oxygen replenishment (either from seawater or the atmosphere), while too small a CF (i.e., very coarse material) would reduce the interaction between microorganisms and oil in the upper intertidal zone due to rapid reduction in the soil moisture at low tide.

Evidence of salt accumulation in beach intertidal zone due to evaporation.

In coastal environments, evaporation is an important driver of subsurface salinity gradients in marsh systems. However, it has not been addressed in the intertidal zone of sandy beaches. Here, we used field data on an estuarine beach foreshore with numerical simulations to show that evaporation causes upper intertidal zone pore-water salinity to be double that of seawater. We found the increase in pore-water salinity mainly depends on air temperature and relative humidity, and tide and wave actions dilute a fraction of the high salinity plume, resulting in a complex process. This is in contrast to previous studies that consider seawater as the most saline source to a coastal aquifer system, thereby concluding that seawater infiltration always increases pore-water salinity by seawater-groundwater mixing dynamics. Our results demonstrate the combined effects of evaporation and tide and waves on subsurface salinity distribution on a beach face. We anticipate our quantitative investigation will shed light on the studies of salt-affected biological activities in the intertidal zone. It also impacts our understanding of the impact of global warming; in particular, the increase in temperature does not only shift the saltwater landward, but creates a different salinity distribution that would have implications on intertidal biological zonation.

BIOB: a mathematical model for the biodegradation of low solubility hydrocarbons.

Modeling oil biodegradation is an important step in predicting the long term fate of oil on beaches. Unfortunately, existing models do not account mechanistically for environmental factors, such as pore water nutrient concentration, affecting oil biodegradation, rather in an empirical way. We present herein a numerical model, BIOB, to simulate the biodegradation of insoluble attached hydrocarbon. The model was used to simulate an experimental oil spill on a sand beach. The biodegradation kinetic parameters were estimated by fitting the model to the experimental data of alkanes and aromatics. It was found that parameter values are comparable to their counterparts for the biodegradation of dissolved organic matter. The biodegradation of aromatics was highly affected by the decay of aromatic biomass, probably due to its low growth rate. Numerical simulations revealed that the biodegradation rate increases by 3-4 folds when the nutrient concentration is increased from 0.2 to 2.0 mg N/L.

A-DROP: A predictive model for the formation of oil particle aggregates (OPAs).

Oil-particle interactions play a major role in removal of free oil from the water column. We present a new conceptual-numerical model, A-DROP, to predict oil amount trapped in oil-particle aggregates. A new conceptual formulation of oil-particle coagulation efficiency is introduced to account for the effects of oil stabilization by particles, particle hydrophobicity, and oil-particle size ratio on OPA formation. A-DROP was able to closely reproduce the oil trapping efficiency reported in experimental studies. The model was then used to simulate the OPA formation in a typical nearshore environment. Modeling results indicate that the increase of particle concentration in the swash zone would speed up the oil-particle interaction process; but the oil amount trapped in OPAs did not correspond to the increase of particle concentration. The developed A-DROP model could become an important tool in understanding the natural removal of oil and developing oil spill countermeasures by means of oil-particle aggregation.

Impacts of evaporation on subsurface flow and salt accumulation in a tidally influenced beach

We coupled a mass-transfer formulation of evaporation to the density-dependent variably saturated finite element model MARUN to predict the salinity under realistic conditions in a beach in the Gulf of Mexico. The results showed that evaporation almost doubled the pore water salt concentration in the intertidal zone in comparison to the case with no evaporation. The results also showed that for the relatively wet atmospheric conditions, the maximum evaporation did not occur from the supratidal zone, rather from the intertidal zone, and it was due to the difference in moisture between pore water and atmosphere. The highest salinity occurred during the spring tide period, and extended into the neap tide period. When the atmospheric condition was assumed to be constant and conducive to large evaporation (high temperature and low air humidity), the maximum salinity region shifted to the supratidal zone and persisted there due to low mixing. The maximum salinity value in the intertidal zone for this case was not larger than that of the variable atmospheric condition case. However, the maximum was more landward than the (more realistic) variable atmospheric condition case, which suggests that accurate characterization of pore water salinity in tidally influenced beaches requires accurate knowledge of atmospheric conditions in addition to beach and hydraulic properties. The constant atmospheric condition simulations suggest that antecedent extreme conditions greatly affect the salinity in the supratidal zone but they are dampened in the intertidal zone.

Oil droplets transport due to irregular waves: Development of large-scale spreading coefficients

The movement of oil droplets due to waves and buoyancy was investigated by assuming an irregular sea state following a JONSWAP spectrum and four buoyancy values. A technique known as Wheeler stretching was used to model the movement of particles under the moving water surface. In each simulation, 500 particles were released and were tracked for a real time of 4.0 h. A Monte Carlo approach was used to obtain ensemble properties. It was found that small eddy diffusivities that decrease rapidly with depth generated the largest horizontal spreading of the plume. It was also found that large eddy diffusivities that decrease slowly with depth generated the smallest horizontal spreading coefficient of the plume. The increase in buoyancy resulted in a decrease in the horizontal spreading coefficient, which suggests that two-dimensional (horizontal) models that predict the transport of surface oil could be overestimating the spreading of oil.

Simulation of oil bioremediation in a tidally influenced beach: Spatiotemporal evolution of nutrient and dissolved oxygen

Numerical experiments of oil bioremediation of tidally influenced beach were simulated using the model BIOMARUN. Nutrient and dissolved oxygen were assumed present in a solution applied on the exposed beach face, and the concentration of these amendments was tracked throughout the beach for up to 6 months. It was found that, in comparison to natural attenuation, bioremediation increased the removal efficiency by 76% and 65% for alkanes and aromatics, respectively. Increasing the nutrient concentration in the applied solution did not always enhance biodegradation as oxygen became limiting even when the beach was originally oxygen-rich. Therefore, replenishment of oxygen to oil-contaminated zone was also essential. Stimulation of oil biodegradation was more evident in the upper and midintertidal zone of the beach, and less in the lower intertidal zone. This was due to reduced nutrient and oxygen replenishment, as very little of the amendment solution reached that zone. It was found that under continual application, most of the oil biodegraded within 2 months, while it persisted for 6 months under natural conditions. While the difference in duration suggests minimal long-term effects, there are situations where the beach would need to be cleaned for major ecological functions, such as temporary nesting or feeding for migratory birds. Biochemical retention time map (BRTM) showed that the duration of solution application was dependent upon the stimulated oil biodegradation rate. By contrast, the application rate of the amendment solution was dependent upon the subsurface extent of the oil-contaminated zone. Delivery of nutrient and oxygen into coastal beach involved complex interaction among amendment solution, groundwater, and seawater. Therefore, approaches that ignore the hydrodynamics due to tide are unlikely to provide the optimal solutions for shoreline bioremediation.