Evaluating the source and seasonality of submarine groundwater discharge using a radon-222 pore water transport model

Pore water radon (²²²Rn) distributions from Indian River Lagoon, Florida, are characterized by three zones: a lower zone where pore water ²²²Rn and sediment-bound radium (²²⁶Ra) are in equilibrium and concentration gradients are vertical; a middle zone where ²²²Rn is in excess of sediment-bound ²²⁶Ra and concentration gradients are concave-downward; and an upper zone where ²²²Rn concentration gradients are nearly vertical. These ²²²Rn data are simulated in a one-dimensional numerical model including advection, diffusion, and non-local exchange to estimate magnitudes of submarine groundwater discharge components (fresh or marine). The numerical model estimates three parameters, fresh groundwater seepage velocity, irrigation intensity, and irrigation attenuation, using two Monte Carlo (MC) simulations that (1) ensure the minimization algorithm converges on a global minimum of the merit function and the parameter estimates are consistent within this global minimum, and (2) provide 90% confidence intervals on the parameter estimates using the measured ²²²Rn activity variance. Model estimates of seepage velocities and discharge agree with previous estimates obtained from numerical groundwater flow models and seepage meter measurements and show the fresh water component decreases offshore and varies seasonally by a factor of nine or less. Comparison between the discharge estimates and precipitation patterns suggests a mean residence time in unsaturated and saturated zones on the order of 5 to 7 months. Irrigation rates generally decrease offshore for all sampling periods. The mean irrigation rate is approximately three times greater than the mean seepage velocity although the ranges of irrigation rates and seepage velocities are the same. Possible mechanisms for irrigation include density-driven convection, wave pumping, and bio-irrigation. Simulation of both advection and irrigation allows the separation of submarine groundwater discharge into fresh groundwater and (re)circulated lagoon water.     

Numerical determination of residence time and age in a partially mixed estuary using three-dimensional hydrodynamic model

This paper documents a numerical modeling study to calculate the residence time and age of dissolved substances in a partially mixed estuary. A three-dimensional, time-dependent hydrodynamic model was established and applied to the Danshuei River estuarine system and adjacent coastal sea in Taiwan. The model showed good agreement with observations of surface elevation, tidal currents and salinity made in 2002. The model was then applied to calculate the residence time and age distribution response to different freshwater discharges with and without density-induced circulations in the Danshuei River estuarine system. Regression analysis of model results reveals that an exponential equation can be used to correlate the residence time to change of freshwater input. The simulated results show it takes approximately 10, 4.5, and 3 days, respectively, for a water parcel that has entered the headwaters of the estuary to be transported out of the estuary under low, mean, and high flow conditions with density-induced circulation. The calculated age with density-induced circulation is less than that without density-induced circulation. The age of the surface layer is less than that at the bottom layer. Overall the study shows that freshwater discharges are the important factors in controlling the transport of dissolved substances in the Danshuei River estuarine system.     

Decadal link between longitudinal morphological changes in branching channels of Yangtze estuary and movement of the offshore depo-center

In estuaries, the morphology of inland and offshore areas usually evolves synergistically. This study examines the decadal link between longitudinal changes in morphology of branching channels and movement of the offshore depo-center (where sediment deposition rate is maximum) of the Yangtze River estuary, under intense human interference. Integrated data analysis is provided on morphology, runoff discharge, and ebb partition ratio from 1950 to 2017. Channel-volume reductions and change rates between isobaths in branching channels reflect the impact of estuarine engineering projects. Ebb partition ratio and duration of discharge ≥ 60 000 m³ s^-1 act as proxies for the water excavating force in branching channels and runoff intensity. It is found that deposition occurs in the lower/upper sub-reaches (or further downstream/upstream channels) of the inland north/south branching channels, and the offshore depo-center moves southward or southeastward, as runoff intensity grows; the reverse occurs as runoff intensity declines. This is because the horizontal circumfluence in the Yangtze estuary rotates clockwise as ebb partition ratios of the north/south branching channels increase/decrease for increasing runoff, and conversely rotates anticlockwise for decreasing runoff. Land reclamation activities, the Deepwater Channel Project, and the Qingcaosha Reservoir have impacted greatly on longitudinal changes of morphology in the North Branch and the South Passage and on ebb partition ratio variations in the North/South Channel and the North/South Passage. Dam-induced runoff flattening has enhanced deposition in the upper/lower sub-reaches of the north/south branching channels and caused northward movement of the offshore depo-center, except in areas affected by estuarine engineering projects. Dam-induced longitudinal evolution of branching channel morphology and offshore depo-center movement will likely persist in the future, given the ongoing construction of large cascade dams in the upper Yangtze and the completion of major projects in the Yangtze estuary. © 2020 John Wiley & Sons, Ltd.     

Determination of groundwater discharge rates and water residence time of groundwater‐fed lakes by stable isotopes of water (18O, 2H) and radon (222Rn) mass balances

Lacustrine groundwater discharge (LGD) and the related water residence time are crucial parameters for quantifying lake matter budgets and assessing its vulnerability to contaminant input. Our approach utilizes the stable isotopes of water (δ¹⁸O, δ²H) and the radioisotope radon (²²²Rn) for determining long‐term average and short‐term snapshots in LGD. We conducted isotope balances for the 0.5‐km² Lake Ammelshainer See (Germany) based on measurements of lake isotope inventories and groundwater composition accompanied by good quality and comprehensive long‐term meteorological and isotopic data (precipitation) from nearby monitoring stations. The results from the steady‐state annual isotope balances that rely on only two sampling campaigns are consistent for both δ¹⁸O and δ²H and suggested an overall long‐term average LGD rate that was used to infer the water residence time of the lake. These findings were supported by the good agreement of the simulated LGD‐driven annual cycles of δ¹⁸O and δ²H lake inventories with the observed lake isotope inventories. However, radon mass balances revealed lower values that might be the result of seasonal LGD variability. For obtaining further insights into possible seasonal variability of groundwater–lake interaction, stable water isotope and radon mass balances could be conducted more frequently (e.g., monthly) in order to use the derived groundwater discharge rates as input for time‐variant isotope balances.     

Influence of seasonal variability in benthic invertebrate community structure on the use of biotic indices to assess the ecological status of a Portuguese estuary

The present study focused on the use of benthic invertebrate communities to assess the ecological quality of a Portuguese estuary characterized by strong seasonal changes and with eutrophication problems. Seasonal benthic samples were collected during a flood year and the methodology proposed by the WFD Portuguese group was used to classify benthic assemblages into five different quality classes. Factor analysis was applied to classify stations based on their physical-chemical status. Different classifications were obtained with different indices and among seasons and there was low agreement between indices and index-season interactions. Diversity indices were better correlated to eutrophication related variables than AMBI and ABC method. Predictable responses of benthic indices to anthropogenic stress symptoms were stronger during the dry period.     

Travel time controls the magnitude of nitrate discharge in groundwater bypassing the riparian zone to a stream on Virginia's coastal plain

Groundwater that bypasses the riparian zone by travelling along deep flow paths may deliver high concentrations of fertilizer‐derived NO₃^? to streams, or it may be impacted by the NO₃^? removal process of denitrification in streambed sediments. In a study of a small agricultural catchment on the Atlantic coastal plain of Virginia's eastern shore, we used seepage meters deployed in the streambed to measure specific discharge of groundwater and its solute concentrations for various locations and dates. We used values of Cl^? concentration to discriminate between bypass water recharged distal to the stream and that contained high NO₃^? but low Cl^? concentrations and riparian‐influenced water recharged proximal to the stream that contained low NO₃^? and high Cl^? concentrations. The travel time required for bypass water to transit the 30‐cm‐thick, microbially active denitrifying zone in the streambed determined the extent of NO₃^? removal, and hydraulic conductivity determined travel time through the streambed sediments. At all travel times greater than 2 days, NO₃^? removal was virtually complete. Comparison of the timescales for reaction and transport through the streambed sediments in this system confirmed that the predominant control on nitrate flux was travel time rather than denitrification rate coefficients. We conclude that extensive denitrification can occur in groundwater that bypasses the riparian zone, but a residence time in biologically active streambed sediments sufficient to remove a large fraction of the NO₃^? is only achieved in relatively low‐conductivity porous media. Instead of viewing them as separate, the streambed and riparian zone should be considered an integrated NO₃^? removal unit. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Long‐term monitoring of a mercury contaminated estuary (Ria de Aveiro,Portugal): the effect of weather events and management in mercury transport

The main aim of this research was to assess the mercury transport from an estuarine basin with a background of anthropogenic contamination during a spring tidal cycle (year 2009) and compare it with two previous tidal cycles (years 1994 and 1999), as part of a long‐term monitoring program. Results showed that effective mercury transport occurs both in the dissolved and particulate fractions (0.18 and 0.20 kg per tidal cycle, respectively), and despite an overall decrease in environmental contamination, results more than double previous findings on particulate transport in the system. These findings result essentially from changes in the tidal prism (net export of 2 million m³ of water), given that both dissolved and particulate concentrations did not increase over time. Hydrodynamic simulations were performed to evaluate the effect of physical disturbance (dredging) and weather events (increased freshwater flow) in these processes, and results suggest the increased freshwater flow into the system as the main forcing function for the mercury transport increment. These results highlight the importance of long‐term monitoring programs, since despite an overall improvement in local contamination levels, the enhancement of transport processes through hydrological changes increases environmental pressure away from the contamination source. Copyright © 2012 John Wiley & Sons, Ltd.     

Delineating coastal groundwater discharge processes in a wetland area by means of electrical resistivity imaging, 224Ra and 222Rn

Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and ²²²Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (²²⁴Ra, ²²²Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.     

Semi‐analytical solution for a slug test in partially penetrating wells including the effect of finite‐thickness skin

This paper presents a new semi‐analytical solution for a slug test in a well partially penetrating a confined aquifer, accounting for the skin effect. This solution is developed based on the solution for a constant‐flux pumping test and a formula given by Peres and co‐workers in 1989. The solution agrees with that of Cooper and co‐workers and the KGS model when the well is fully penetrating. The present solution can be applied to simulate the temporal and spatial head distributions in both the skin and formation zones. It can also be used to demonstrate the influences of skin type or skin thickness on the well water level and to estimate the hydraulic parameters of the skin and formation zones using a least‐squares approach. The results of this study indicate that the determination of hydraulic conductivity using a conventional slug‐test data analysis that neglects the presence of a skin zone will give an incorrect result if the aquifer has a skin zone. Copyright © 2008 John Wiley & Sons, Ltd.     

Real‐time monitoring and estimation of the discharge of flash floods in a steep mountain catchment

Synchronously and accurately estimating the flood discharges and dynamic changes in the fluid density is essential for hydraulic analysis and forecasting of flash floods, as well as for risk assessment. However, such information is rare for steep mountain catchments, especially in regions that are hotspots for earthquakes. Therefore, six hydrological monitoring sites were established in the main stream and tributaries of the 78.3‐km² Longxi River catchment, an affected region of the Wenchuan earthquake region in China. Direct real‐time monitoring equipment was installed to measure the flow depths, velocities, and fluid total pressures of the flood hydrographs. On the basis of field measurements, real‐time mean cross‐sectional velocities during the flood hydrographs could be derived from easily obtainable parameters: cross‐sectional maximum velocities and the calibrated dimensionless parameter K_h . Real‐time discharges were determined on the basis of a noncontact method to establish the effective rating curves of this mountainous stream, ranging from 1.46 to 386.34 m³/s with the root mean square errors of ≤10.22 m³/s. Compared with the traditional point‐velocity method and empirical Manning's formula, the proposed noncontact method was reliable and safe for monitoring whole flood hydrographs. Additionally, the real‐time fluid density during the flood hydrographs was calculated on the basis of the direct monitoring parameters for fluid total pressures and water depths. During the flood hydrograph, transient flow behaviour with higher fluid density generally occurred downstream during the flood peak periods when the flow was in the supercritical flow regime. The observed behaviour greatly increased the threat of damage to infrastructure and human life near the river. Thus, it is important to accurately estimate flood discharge and identify for fluid densities so that people at risk from an impending flash flood are given reliable, advanced warning.     

Role of bedrock groundwater in the rainfall–runoff process in a small headwater catchment underlain by volcanic rock

The role of bedrock groundwater in rainfall–runoff processes is poorly understood. Hydrometric, tracer and subsurface water potential observations were conducted to study the role of bedrock groundwater and subsurface flow in the rainfall–runoff process in a small headwater catchment in Shiranui, Kumamoto prefecture, south‐west Japan. The catchment bedrock consists of a strongly weathered, fractured andesite layer and a relatively fresh continuous layer. Major chemical constituents and stable isotopic ratios of δ¹⁸O and δD were analysed for spring water, rainwater, soil water and bedrock groundwater. Temporal and spatial variation in SiO₂ showed that stream flow under the base flow condition was maintained by bedrock groundwater. Time series of three components of the rainstorm hydrograph (rainwater, soil water and bedrock groundwater) separated by end member mixing analysis showed that each component fluctuated during rainstorm, and their patterns and magnitudes differed between events. During a typical mid‐magnitude storm event, a delayed secondary runoff peak with 1·0 l s⁻¹ was caused by increase in the bedrock groundwater component, whereas during a large rainstorm event the bedrock groundwater component increased to ≈ 2·5 l s⁻¹. This research shows that the contribution of bedrock groundwater and soil water depends strongly on the location of the groundwater table, i.e. whether or not it rises above the soil–bedrock interface. Copyright © 2010 John Wiley & Sons, Ltd.