Data‐Driven Approach for Analyzing Hydrogeology and Groundwater Quality Across Multiple Scales

Recent trends of assimilating water well records into statewide databases provide a new opportunity for evaluating spatial dynamics of groundwater quality and quantity. However, these datasets are scarcely rigorously analyzed to address larger scientific problems because they are of lower quality and massive. We develop an approach for utilizing well databases to analyze physical and geochemical aspects of groundwater systems, and apply it to a multiscale investigation of the sources and dynamics of chloride (Cl^?) in the near‐surface groundwater of the Lower Peninsula of Michigan. Nearly 500,000 static water levels (SWLs) were critically evaluated, extracted, and analyzed to delineate long‐term, average groundwater flow patterns using a nonstationary kriging technique at the basin‐scale (i.e., across the entire peninsula). Two regions identified as major basin‐scale discharge zones—the Michigan and Saginaw Lowlands—were further analyzed with regional‐ and local‐scale SWL models. Groundwater valleys (“discharge” zones) and mounds (“recharge” zones) were identified for all models, and the proportions of wells with elevated Cl^? concentrations in each zone were calculated, visualized, and compared. Concentrations in discharge zones, where groundwater is expected to flow primarily upwards, are consistently and significantly higher than those in recharge zones. A synoptic sampling campaign in the Michigan Lowlands revealed concentrations generally increase with depth, a trend noted in previous studies of the Saginaw Lowlands. These strong, consistent SWL and Cl^? distribution patterns across multiple scales suggest that a deep source (i.e., Michigan brines) is the primary cause for the elevated chloride concentrations observed in discharge areas across the peninsula.     

Contribution of baseflow nitrate export to non-point source pollution

As a common pollutant of nitrogen in groundwater, nitrate contamination has become a major concern worldwide. Baseflow, one of the dominant hydrological pathways for nitrate migration to streamflow, has been confirmed as a leading nitrate source for stream water where groundwater or subsurface flow contaminated heavily by nitrate. That is, sufficient improvements of water quality may not be attained without proper management for baseflow, even if non-point sources (NPS) pollutants discharged through surface runoff are being well managed. This article reviews the primary nitrate sources, the main factors affecting its transport, and the methodologies for baseflow nitrate estimation, to give some recommendations for future works, including: (1) giving sufficient consideration for the effects of climatological, morphological, and geological factors on baseflow recessions to obtain more reliable and accurate baseflow separation; (2) trying to solve calibration and validation problems for baseflow loads determining in storm flow period; (3) developing a simple and convenient algorithm with certain physics that can be used to separate baseflow NPS pollution from the total directly in different regions, for a reliable estimation of baseflow NPS pollution at larger scale (e.g., national scale); (4) improving groundwater quality simulation module of existing NPS pollution models to have a better simulation for biogeochemical processes in shallow aquifers; (5) taking integrated measures of “source control”, “process interception” and “end remediation” to prevent and control NPS nitrate pollution effectively, not just only the strict control of nutrients loss from surface runoff.     

Convergence of Groundwater Discharge through the Hyporheic Zone of Streams

Significant attention has been given to hyporheic water fluxes induced by hydromorphologic processes in streambeds and the effects they have on stream ecology. However, the impact of hyporheic fluxes on regional groundwater flow discharge zones as well as the interaction of these flows are much less investigated. The groundwater-hyporheic interactive flow not only governs solute mass and heat transport in streams but also controls the retention of solute and contamination following the discharge of deep groundwater, such as naturally occurring solutes and leakage from geological waste disposal facilities. Here, we applied a physically based modeling approach combined with extensive hydrologic, geologic and geographical data to investigate the effect of hyporheic flow on groundwater discharge in the Krycklan catchment, located in a boreal landscape in Sweden. Regional groundwater modeling was conducted using COMSOL Multiphysics by considering geologic heterogeneity and infiltration constraint of the groundwater circulation intensity. Moreover, the hyporheic flow was analyzed using an exact spectral solution accounting for the fluctuating streambed topography and superimposed with the regional groundwater flow. By comparing the discharge flow fields with and without consideration of hyporheic flows, we found that the divergence of the discharge was substantially enhanced and the distribution of the travel times of groundwater was significantly shifted toward shorter times due to the presence of hyporheic flow. Particularly important is that the groundwater flow paths contract near the streambed interface due to the hyporheic flow, which leads to a phenomenon that we name “fragmentation” of coherent areas of groundwater upwelling in pinhole-shaped stream tubes.     

Understanding fen hydrology across multiple scales

Fens, which are among the most biodiverse of wetland types in the USA, typically occur in glacial landscapes characterized by geo‐morphologic variability at multiple spatial scales. As a result, the hydrologic systems that sustain fens are complex and not well understood. Traditional approaches for characterizing such systems use simplifying assumptions that cannot adequately capture the impact of variability in geology and topography. In this study, a hierarchical, multi‐scale groundwater modelling approach coupled with a geologic model is used to understand the hydrology of a fen in Michigan. This approach uses high‐resolution data to simulate the multi‐scale topographic and hydrologic framework and lithologic data from more than 8500 boreholes in a statewide water well database to capture the complex geology. A hierarchy of dynamically linked models is developed that simulates groundwater flow at all scales of interest and to delineate the areas that contribute groundwater to the fen. The results show the fen receiving groundwater from multiple sources: an adjacent wetland, local recharge, a nearby lake and a regional groundwater mound. Water from the regional mound flows to an intermediate source before reaching the fen, forming a ‘cascading’ connection, while other sources provide water through ‘direct’ connections. The regional mound is also the source of water to other fens, streams and lakes in this area, thus creating a large, interconnected hydrologic system that sustains the entire ecosystem. In order to sustainably manage such systems, conservation efforts must include both site‐based protection and management, as well as regional protection and management of groundwater source areas. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of Flow in a Complex Aquifer System Subjected to Long-Term Well Network Growth

In west-central Lower Peninsula of Michigan, population growth and expanded agricultural activities over recent decades have resulted in significant increases in distributed groundwater withdrawals. The growth of the extensive well network and anecdotes of water shortages (dry wells) have raised concerns over the region's groundwater sustainability. We developed an unsteady, three-dimensional (3D) groundwater flow model to describe system dynamics over the last 50 years and evaluate long-term impacts of groundwater use. Simulating this large aquifer system was challenging; the site is characterized by strong, spatially distributed, and statistically nonstationary heterogeneity, making it difficult to avoid over-parameterization using traditional approaches for conceptualizing and calibrating a flow model. Moreover, traditional pumping and water level data were lacking and prohibitively expensive to collect given the large-scale and long-term nature of this study. An integrated, stochastic-deterministic approach was developed to characterize the system and calibrate the flow model through innovative use of high-density water well datasets. This approached allowed (1) implementation of a “zone-based,” nonstationary stochastic approach to conceptualize complex spatial variability using a small set of geologic material types; (2) modeling the spatiotemporal evolution of many water well withdrawals across several decades using sector-based parameterization; and (3) critical analysis of long-term water level changes at different locations in the aquifer system for characterizing the system dynamics and calibrating the model. Results show the approach is reasonably successful in calibrating a complex model for a highly complex site in a way that honors complex distributed heterogeneity and stress configurations.     

Isotopically Enriched Geogenic δ81Br and δ37Cl: Primary Evidence for the Ascending Brine Model

Mass balance calculations and hydrodynamics of groundwater flow suggest that the solutes in brines of the coastal sabkha aquifer from the Emirate of Abu Dhabi are derived largely from ascending geologic brines into the sabkha from the underlying formations. Solute interpretation for the ascending brine model (ABM) was based on two independent but secondary lines of evidence (solute ratios and solute fluxes). In the current study, direct primary evidence for this ABM was provided through analyses of δ⁸¹Br, δ³⁷Cl, and ⁸⁷Sr/⁸⁶Sr. Different solute histories of geologic brine and sea water provide an “isotopic fingerprint” that can uniquely distinguish between the two possible sources. Samples from the coastal sabkha aquifer of Abu Dhabi were determined to have a mean δ⁸¹Br of 1.17‰ that is statistically equal, at the 95% confidence level, to the mean of 1.11‰ observed in the underlying geologic brine and statistically different than sea water. Similarly, the δ³⁷Cl in sabkha brine has a mean of 0.25‰ and is statistically equal to a mean of 0.21‰ in the underlying geologic brines at the 95% confidence level and statistically different from sea water. Also, dissolved strontium isotope data are consistent with the ABM and even with the complex set of processes in the sabkha, the variance in strontium isotope results is similar to the geologic brine. These observations provide primary direct evidence consistent that the major source of these solutes (and presumably others in the aquifer) is from discharging geologic brines, not from adjacent sea water.     

Potential Contaminant Pathways from Hydraulically Fractured Shale to Aquifers

Hydraulic fracturing of deep shale beds to develop natural gas has caused concern regarding the potential for various forms of water pollution. Two potential pathways—advective transport through bulk media and preferential flow through fractures—could allow the transport of contaminants from the fractured shale to aquifers. There is substantial geologic evidence that natural vertical flow drives contaminants, mostly brine, to near the surface from deep evaporite sources. Interpretative modeling shows that advective transport could require up to tens of thousands of years to move contaminants to the surface, but also that fracking the shale could reduce that transport time to tens or hundreds of years. Conductive faults or fracture zones, as found throughout the Marcellus shale region, could reduce the travel time further. Injection of up to 15,000,000 L of fluid into the shale generates high pressure at the well, which decreases with distance from the well and with time after injection as the fluid advects through the shale. The advection displaces native fluids, mostly brine, and fractures the bulk media widening existing fractures. Simulated pressure returns to pre‐injection levels in about 300 d. The overall system requires from 3 to 6 years to reach a new equilibrium reflecting the significant changes caused by fracking the shale, which could allow advective transport to aquifers in less than 10 years. The rapid expansion of hydraulic fracturing requires that monitoring systems be employed to track the movement of contaminants and that gas wells have a reasonable offset from faults.     

Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution,spreadsheet tool,and field applications

Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction‐dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson ( 2004 ), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater‐surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux‐LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi‐decadal surface signals have disturbed subsurface thermal regimes.     

Groundwater Modeling in Integrated Water Resources Management—Visions for 2020

Groundwater modeling is undergoing a change from traditional stand-alone studies toward being an integrated part of holistic water resources management procedures. This is illustrated by the development in Denmark, where comprehensive national databases for geologic borehole data, groundwater-related geophysical data, geologic models, as well as a national groundwater-surface water model have been established and integrated to support water management. This has enhanced the benefits of using groundwater models. Based on insight gained from this Danish experience, a scientifically realistic scenario for the use of groundwater modeling in 2020 has been developed, in which groundwater models will be a part of sophisticated databases and modeling systems. The databases and numerical models will be seamlessly integrated, and the tasks of monitoring and modeling will be merged. Numerical models for atmospheric, surface water, and groundwater processes will be coupled in one integrated modeling system that can operate at a wide range of spatial scales. Furthermore, the management systems will be constructed with a focus on building credibility of model and data use among all stakeholders and on facilitating a learning process whereby data and models, as well as stakeholders' understanding of the system, are updated to currently available information. The key scientific challenges for achieving this are (1) developing new methodologies for integration of statistical and qualitative uncertainty; (2) mapping geological heterogeneity and developing scaling methodologies; (3) developing coupled model codes; and (4) developing integrated information systems, including quality assurance and uncertainty information that facilitate active stakeholder involvement and learning.     

Hydrological impact of roadside ditches in an agricultural watershed in Central New York: implications for non‐point source pollutant transport

Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint‐source (NPS) pollution by acting as hydrological links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field‐ and watershed‐scale hydrology, and pollutant transport in a 38 km² agricultural watershed in south‐central NY. Salient findings include the following: (i) 94% of road ditches discharged to natural streams, effectively doubling the drainage density; (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches; and (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.     

Tracking Groundwater Discharge to a Large River using Tracers and Geophysics

Few studies have investigated large reaches of rivers in which multiple sources of groundwater are responsible for maintaining baseflow. This paper builds upon previous work undertaken along the Fitzroy River, one of the largest perennial river systems in north‐western Australia. Synoptic regional‐scale sampling of both river water and groundwater for a suite of environmental tracers (⁴He, ⁸⁷Sr/⁸⁶Sr, ²²²Rn and major ions), and subsequent modeling of tracer behavior in the river, has enabled definition and quantification of groundwater input from at least three different sources. We show unambiguous evidence of both shallow “local” groundwater, possibly recharged to alluvial aquifers beneath the adjacent floodplain during recent high‐flow events, and old “regional” groundwater introduced via artesian flow from deep confined aquifers. We also invoke hyporheic exchange and either bank return flow or parafluvial flow to account for background ²²²Rn activities and anomalous chloride trends along river reaches where there is no evidence of the local or regional groundwater inputs. Vertical conductivity sections acquired through an airborne electromagnetic (AEM) survey provide insights to the architecture of the aquifers associated with these sources and general groundwater quality characteristics. These data indicate fresh groundwater from about 300 m below ground preferentially discharging to the river, at locations consistent with those inferred from tracer data. The results demonstrate how sampling rivers for multiple environmental tracers of different types—including stable and radioactive isotopes, dissolved gases and major ions—can significantly improve conceptualization of groundwater—surface water interaction processes, particularly when coupled with geophysical techniques in complex hydrogeological settings.     

Effects of shield brine on the safe disposal of waste in deep geologic environments

The salinity of groundwater increases with depth in the Canadian Shield (up to 1.3 kg/L of density). The existence of brine can be critically important for the safe geologic disposal of radioactive wastes, as dense brine can significantly retard the upward migration of radionuclides released from repositories. Static and flushing conditions of the deep brine are analyzed using a U-tube analogy model. Velocity reduction due to the presence of dense brine is derived under flushing conditions. A set of illustrative numerical simulations in a two-dimensional cross section is presented to demonstrate that dense brine can significantly influence regional groundwater flow patterns in a shield environment. It is implied from the results that (1) the existence of Shield brine can be an indicator of a hydrogeologically stable environment, (2) activities near ground surface may not perturb the stable groundwater environment in the deep brine region, and thus, (3) the deep brine region can be considered as a candidate geologic site for the safe disposal of waste. In addition to brine, other issues associated with long-term waste disposal, such as geological, glacial and seismic events, may need to be considered for the safe storage of spent nuclear fuel in a shield environment.     

Wetland-Scale Mapping of Preferential Fresh Groundwater Discharge to the Colorado River

Quantitative evaluation of groundwater/surface water exchange dynamics is universally challenging in large river systems, because existing methodology often does not yield spatially-distributed data and is difficult to apply in deeper water. Here we apply a combined near-surface geophysical and direct groundwater chemical toolkit to refine fresh groundwater discharge estimates to the Colorado River through a 4-km² wetland that borders the town of Moab, Utah, USA. Preliminary characterization of raw electromagnetic imaging (EMI) data, collected by kayak and by walking, was used to guide additional direct-contact electrical measurements and installation of new monitoring wells. Chemical data from the wells strongly supported the EMI spatial characterization of preferential fresh groundwater discharge embedded in natural brine groundwaters and weighted to the southern wetland section. Inversion of the EMI data revealed sub-meter scale detail regarding bulk electrical conductivity zonation across approximately 15.5 km of transects, collected in only 3 days. This electrical detail indicates processes such as salinization of the unsaturated zone and direct discharge through the Colorado River sediments and a tributary creek bed. Overall, the study contributed to a substantial reduction in fresh groundwater discharge estimates previously made using sparse existing well data and a simplified assumption of diffuse fresh groundwater discharge below the entire wetland. EMI will likely become a widely used tool in systems with natural electrical contrast as groundwater/surface water hydrogeologists continue to recognize the prevalence of preferential groundwater discharge processes.     

Groundwater potential zones using a combination of geospatial technology and geophysical approach: case study in Dehradun,India

ABSTRACT

A combination of geospatial, geophysical and statistical models using satellite data, the weighted index overlay (WIO) method and two-dimensional electrical resistivity tomography (2D-ERT) is applied to generate the highest potential groundwater area and to further explore the groundwater in Dehradun, India. The results show that of 19.7 km² total basin area, 0.26% falls under the “poor” category as a prospect zone for groundwater, 4.3% is “moderate”, 10.10% “moderately good”, 4.9% “good” and 0.17% “very good”. In addition, the demonstration of the geophysical survey is presented, in which Purkal Youth Society Division (PYSD) site is categorized as a shallow aquifer zone and the Guru Nanak Fifth Centenary School (GNFCS) site is a deeper aquifer zone. Our study emphasizes remote sensing and geographic information system integrated with a geophysical survey to support prospecting the most probable area and confirm the existence of groundwater.     

Use of the MODFLOW 6 Water Mover Package to Represent Natural and Managed Hydrologic Connections

The latest release of MODFLOW 6, the current core version of the MODFLOW groundwater modeling software, debuted a new package dubbed the “mover” (MVR). Using a generalized approach, MVR facilitates the transfer of water among any arbitrary combination of simulated features (i.e., pumping wells, stream, drains, lakes, etc.) within a MODFLOW 6 simulation. Four “rules” controlling the amount of water transferred from a providing feature to a receiving feature are currently available. In this way, MVR can represent natural connections between features, for example, streams entering or exiting lakes, and perhaps more interestingly, it also can transfer water among simulated features to more accurately simulate water management. An example model representative of an agricultural setting demonstrates some of the available MVR connections. For example, an irrigation event that transfers surface water from an irrigation delivery ditch to multiple cropped areas demonstrates a “one-to-many” connection that is possible within MVR. Conversely, irrigation or precipitation runoff from multiple fields may be routed to a particular stream segment using “many-to-one” MVR connections. MVR supports many additional connection types, several of which are demonstrated by the included example problem.     

Origin and dynamics of groundwater salinity in the alluvial plains of western Delhi and adjacent territories of Haryana State,India

Groundwater salinity is a widespread problem and a challenge to water resources management. It is an increasing concern in the alluvial plains of Delhi and neighbouring Haryana state as well as a risk for agricultural production water supply and sustainable development. This study aims to identify potential sources of dissolved salts and the driving mechanisms of salinity ingress in the shallow aquifer. It combines a comprehensive review of environmental conditions and the analysis of groundwater samples from 25 sampling points. Major ions are analysed to describe the composition and distribution of saline groundwater and dissolution/precipitation dynamics. Density stratification and local upconing of saline waters were identified by multilevel monitoring and temperature logging. Bromide–chloride ratios hold information on the formation of saline waters, and nitrate is used as an indicator for anthropogenic influences. In addition, stable isotope analysis helps to identify evaporation and to better understand recharge processes and mixing dynamics in the study region. The results lead to the conclusion that surface water and groundwater influx into the poorly drained semiarid basin naturally results in the accumulation of salts in soil, sediments and groundwater. Human‐induced changes of environmental conditions, especially the implementation of traditional canal and modern groundwater irrigation, have augmented evapotranspiration and led to waterlogging in large areas. In addition, water‐level fluctuations and perturbation of the natural hydraulic equilibrium favour the mobilisation of salts from salt stores in the unsaturated zone and deeper aquifer sections. The holistic approach of this study demonstrates the importance of various salinity mechanisms and provides new insights into the interference of natural and anthropogenic influences. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling non-stationary water ages in a tropical rainforest: A preliminary spatially distributed assessment

Pristine tropical forests play a critical role in regional and global climate systems. For a better understanding of the eco-hydrology of tropical “evergreen” vegetation, it is essential to know the partitioning of water into transpiration and evaporation, runoff and associated water ages. For this purpose, we evaluated how topography and vegetation influence water flux and age dynamics at high temporal (hourly) and spatial (10 m) resolution using the Spatially Distributed Tracer-Aided Rainfall-Runoff model for the tropics (STARRtropics). The model was applied in a tropical rainforest catchment (3.2 km²) where data were collected biweekly to monthly and during intensive monitoring campaigns from January 2013 to July 2018. The STARRtropics model was further developed, incorporating an isotope mass balance for evapotranspiration partitioning into transpiration and evaporation. Results exhibited a rapid streamflow response to rainfall inputs (water and isotopes) with limited mixing and a largely time-invariant baseflow isotope composition. Simulated soil water storage showed a transient response to rainfall inputs with a seasonal component directly resembling the streamflow dynamics which was independently evaluated using soil water content measurements. High transpiration fluxes (max 7 mm/day) were linked to lower slope gradients, deeper soils and greater leaf area index. Overall water partitioning resulted in 65% of the actual evapotranspiration being driven by vegetation with high transpiration rates over the drier months compared to the wet season. Time scales of water age were highly variable, ranging from hours to a few years. Stream water ages were conceptualized as a mixture of younger soil water and slightly older, deeper soil water and shallow groundwater with a maximum age of roughly 2 years during drought conditions (722 days). The simulated soil water ages ranged from hours to 162 days and for shallow groundwater up to 1,200 days. Despite the model assumptions, experimental challenges and data limitation, this preliminary spatially distributed model study enhances knowledge about the water ages and overall young water dominance in a tropical rainforest with little influence of deeper and older groundwater.     

Spatial patterns of groundwater‐lake exchange – implications for acid neutralization processes in an acid mine lake

Exchange of groundwater and lake water with typically quite different chemical composition is an important driver for biogeochemical processes at the groundwater‐lake interface, which can affect the water quality of lakes. This is of particular relevance in mine lakes where anoxic and slightly acidic groundwater mixes with oxic and acidic lake water (pH < 3). To identify links between groundwater‐lake exchange rates and acid neutralization processes in the sediments, exchange rates were quantified and related to pore‐water pH, sulfate and iron concentrations as well as sulfate reduction rates within the sediment. Seepage rates measured with seepage meters (?2.5 to 5.8 L m^‐2 d^‐1) were in reasonable agreement with rates inverted from modeled chloride profiles (?1.8 to 8.1 L m^‐2 d^‐1). Large‐scale exchange patterns were defined by the (hydro)geologic setting but superimposed by smaller scale variations caused by variability in sediment texture. Sites characterized by groundwater upwelling (flow into the lake) and sites where flow alternated between upwelling and downwelling were identified. Observed chloride profiles at the alternating sites reflected the transient flow regime. Seepage direction, as well as seepage rate, were found to influence pH, sulfate and iron profiles and the associated sulfate reduction rates. Under alternating conditions proton‐consuming processes, for example, sulfate reduction, were slowed. In the uppermost layer of the sediment (max. 5 cm), sulfate reduction rates were significantly higher at upwelling (>330 nmol g^‐1 d^‐1) compared to alternating sites (<220 nmol g^‐1 d^‐1). Although differences in sulfate reduction rates could not be explained solely by different flux rates, they were clearly related to the prevailing groundwater‐lake exchange patterns and the associated pH conditions. Our findings strongly suggest that groundwater‐lake exchange has significant effects on the biogeochemical processes that are coupled to sulfate reduction such as acidity retention and precipitation of iron sulfides. Copyright © 2012 John Wiley & Sons, Ltd.     

Modeling fresh water lens damage and recovery on atolls after storm-wave washover

The principal natural source of fresh water on scattered coral atolls throughout the tropical Pacific Ocean is thin unconfined groundwater lenses within islet substrates. Although there are many threats to the viability of atoll fresh water lenses, salinization caused by large storm waves washing over individual atoll islets is poorly understood. In this study, a mathematical modeling approach is used to examine the immediate responses, longer-term behavior, and subsequent (partial) recovery of a Pacific atoll fresh water lens after saline damage caused by cyclone-generated wave washover under different scenarios. Important findings include: (1) the saline plume formed by a washover event mostly migrates downward first through the top coral sand and gravel substrate, but then exits the aquifer to the ocean laterally through the more permeable basement limestone; (2) a lower water table position before the washover event, rather than a longer duration of storm washover, causes more severe damage to the fresh water lens; (3) relatively fresher water can possibly be found as a preserved horizon in the deeper part of an aquifer after disturbance, especially if the fresh water lens extends into the limestone under normal conditions; (4) post-cyclone accumulation of sea water in the central depression (swamp) of an atoll islet prolongs the later stage of fresh water lens recovery.     

Pollution Source Investigation and Water Quality Management in the Carp Lake Watershed,Taiwan

In this study, a full survey of pollutant sources and water quality was conducted, followed by the application of a water quality model (Water Quality Analysis Simulation Program, WASP) to establish strategies of water quality control in Carp Lake, Taiwan. Results of the field investigation show that both point and non‐point source (NPS) pollutants were responsible for the poor water quality. The contributions of biochemical oxygen demand (BOD) from point source and NPS pollution were 45.9 and 55.1%, respectively. About 80% of total phosphorus (TP) were contributed by NPS. Additionally, point source and NPS pollution discharged 55.5 and 44.5% of NH₃–N load, respectively. The Carlson's Trophic State Index ranged from 61.9 to 69.2 showing serious eutrophic problems in Carp Lake. The calculated BOD, NH₃–N, and TP carrying capacity were approximately 2.8, 0.42, and 0.15 kg per day, respectively. However, the current pollutant loadings are approximately 3.0–5.5 times the calculated carrying capacity. With the help of the calibrated WASP model, remedial strategies for the lake water from short‐term to long‐term were developed. The completion of the small local sewer system to remove 80% of the point source pollution can serve as a short‐term goal while 40–60% of NPS removal by natural treatment systems may serve as a mid‐term goal. Furthermore, 80% of both source point and NPS pollution removal can be considered as a long‐term strategy. Results of heavy metal analysis show that the enriched sediment would be safe for agricultural applications.