Disconnected Surface Water and Groundwater: From Theory to Practice

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923) , it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.     

U-series disequilibria in Kick’em Jenny submarine volcano lavas: A new view of time-scales of magmatism in convergent margins

We present data for U and its decay series nuclides ²³⁰Th, ²²⁶Ra, ²³¹Pa, and ²¹⁰Po for 14 lavas from Kick’em Jenny (KEJ) submarine volcano to constrain the time-scales and processes of magmatism in the Southern Lesser Antilles, the arc having the globally lowest plate convergence rate. Although these samples are thought to have been erupted in the last century, most have (²²⁶Ra)/(²¹⁰Po) within ±15% of unity. Ten out of 14 samples have significant ²²⁶Ra excesses over ²³⁰Th, with (²²⁶Ra)/(²³⁰Th) up to 2.97, while four samples are in ²²⁶Ra-²³⁰Th equilibrium within error. All KEJ samples have high (²³¹Pa)/(²³⁵U), ranging from 1.56 to 2.64 and high ²³⁸U excesses (up to 43%), providing a global end-member of high ²³⁸U and high ²³¹Pa excesses. Negative correlations between Sr, sensitive to plagioclase fractionation, and Ho/Sm, sensitive to amphibole fractionation, or K/Rb, sensitive to open system behavior, indicate that differentiation at KEJ lavas was dominated by amphibole fractionation and open-system assimilation. While (²³¹Pa)/(²³⁵U) does not correlate with differentiation indices such as Ho/Sm, (²³⁰Th)/(²³⁸U) shows a slight negative correlation, likely due to assimilation of materials with slightly higher (²³⁰Th)/(²³⁸U). Samples with ²²⁶Ra excess have higher Sr/Th and Ba/Th than those in ²²⁶Ra-²³⁰Th equilibrium, forming rough positive correlations of (²²⁶Ra)/(²³⁰Th) with Sr/Th and Ba/Th similar to those observed in many arc settings. We interpret these correlations to reflect a time-dependent magma differentiation process at shallow crustal levels and not the process of recent fluid addition at the slab-wedge interface.The high ²³¹Pa excesses require an in-growth melting process operating at low melting rates and small residual porosity; such a model will also produce significant ²³⁸U-²³⁰Th and ²²⁶Ra-²³⁰Th disequilibrium in erupted lavas, meaning that signatures of recent fluid addition from the slab are unlikely to be preserved in KEJ lavas. We instead propose that most of the ²³⁸U-²³⁰Th, ²²⁶Ra-²³⁰Th, and ²³⁵U-²³¹Pa disequilibria in erupted KEJ lavas reflect the in-growth melting process in the mantle wedge (reflecting variations in U/Th, daughter-parent ratios, fO₂, and thermal structure), followed by modification by magma differentiation at crustal depths. Such a conclusion reconciles the different temporal implications from different U-series parent-daughter pairs and relaxes the time constraint on mass transfer from slab to eruption occurring in less than a few thousand years imposed by models whereby ²²⁶Ra excess is derived from the slab.     

Understanding the fate of iron in a modern temperate estuary: Leirárvogur,Iceland

Fluvial dissolved Fe concentrations decrease upon mixing with seawater, resulting in the formation of Fe-floccules. However, a clear understanding of the fate of these floccules has yet to be established. Assessing how tidal processes affect the formation of Fe-colloids in the Leirárvogur estuary, SW Iceland, is an important step in understanding the formation and potential deposition of estuarine Fe-rich minerals within this estuarine system. The Leirárvogur estuary drains predominately Fe-rich basalt, increasing the likelihood of detecting changes in Fe-phases. Fluvial waters and local lake waters that drain into the estuary were compared and the effects of seasonal changes were considered, in an attempt to understand how varying end-members and external factors play a role in Fe-rich mineral formation. Aqueous and colloidal Fe concentrations were found to be greater towards the head of the Leirárvogur estuary, suggesting that potential Fe-rich minerals and complexes are forming at sites of fluvial input. Increasing suspended colloidal Fe towards the estuary mouth suggests that Fe-colloids are readily transported seaward.     

The effects of climate and landscape position on chemical denudation and mineral transformation in the Santa Catalina mountain critical zone observatory

Understanding the interactions of climate, physical erosion, chemical weathering and pedogenic processes is essential when considering the evolution of critical zone systems. Interactions among these components are particularly important to predicting how semiarid landscapes will respond to forecasted changes in precipitation and temperature under future climate change. The primary goal of this study was to understand how climate and landscape structure interact to control chemical denudation and mineral transformation across a range of semiarid ecosystems in southern Arizona. The research was conducted along the steep environmental gradient encompassed by the Santa Catalina Mountains Critical Zone Observatory (SCM-CZO). The gradient is dominated by granitic parent materials and spans significant range in both mean annual temperature (>10 °C) and precipitation (>50 cm a^?1), with concomitant shift in vegetation communities from desert scrub to mixed conifer forest. Regolith profiles were sampled from divergent and convergent landscape positions in five different ecosystems to quantify how climate-landscape position interactions control regolith development. Regolith development was quantified as depth to paralithic contact and degree of chemical weathering and mineral transformation using a combination of quantitative and semi-quantitative X-ray diffraction (XRD) analyses of bulk soils and specific particle size classes. Depth to paralithic contact was found to increase systematically with elevation for divergent positions at approximately 28 cm per 1000 m elevation, but varied inconsistently for convergent positions. The relative differences in depth between convergent and divergent landscape positions was greatest at the low and high elevation sites and is hypothesized to be a product of changes in physical erosion rates across the gradient. Quartz/Plagioclase (Q/P) ratios were used as a general proxy for bulk regolith chemical denudation. Q/P was generally higher in divergent landscape positions compared to the adjacent convergent hollows. Convergent landscape positions appear to be collecting solute-rich soil–waters from divergent positions thereby inhibiting chemical denudation. Clay mineral assemblage of the low elevation sites was dominated by smectite and partially dehydrated halloysite whereas vermiculite and kaolinite were predominant in the high elevation sites. The increased depth to paralithic contact, chemical denudation and mineral transformation are likely functions of greater water availability and increased primary productivity. Landscape position within a given ecosystem exerts strong control on chemical denudation as a result of the redistribution of water and solutes across the landscape surface. The combined data from this research demonstrates a strong interactive control of climate, landscape position and erosion on the development of soil and regolith.     

A 50-year record of NO<Subscript>x</Subscript> and SO<Subscript>2</Subscript> sources in precipitation in the Northern Rocky Mountains,USA

Ice-core samples from Upper Fremont Glacier (UFG), Wyoming, were used as proxy records for the chemical composition of atmospheric deposition. Results of analysis of the ice-core samples for stable isotopes of nitrogen (δ¹⁵N, ) and sulfur (δ³⁴S, ), as well as and deposition rates from the late-1940s thru the early-1990s, were used to enhance and extend existing National Atmospheric Deposition Program/National Trends Network (NADP/NTN) data in western Wyoming. The most enriched δ³⁴S value in the UFG ice-core samples coincided with snow deposited during the 1980 eruption of Mt. St. Helens, Washington. The remaining δ³⁴S values were similar to the isotopic composition of coal from southern Wyoming. The δ¹⁵N values in ice-core samples representing a similar period of snow deposition were negative, ranging from -5.9 to -3.2 ‰ and all fall within the δ¹⁵N values expected from vehicle emissions. Ice-core nitrate and sulfate deposition data reflect the sharply increasing U.S. emissions data from 1950 to the mid-1970s.     

Detection of deep water benthic macroalgae using image-based classification techniques on multibeam backscatter at Cashes Ledge,Gulf of Maine,USA

Benthic macroalgae form an important part of temperate marine ecosystems, exhibiting a complex three-dimensional character which represents a vital foraging and spawning ground for many juvenile fish species. In this research, image-based techniques for classification of multibeam backscatter are explored for the detection of benthic macroalgae at Cashes Ledge in the Gulf of Maine, USA. Two classifications were performed using QTC-Multiview, differentiated by application of a threshold filter, and macroalgal signatures were independently extracted from the raw sonar datagrams in Matlab. All classifications were validated by comparison with video ground-truth data. The unfiltered classification shows a high degree of complexity in the shallowest areas within the study site; the filtered demonstrates markedly less variation by depth. The unfiltered classification shows a positive agreement with the video ground-truth data; 82.6% of observations recording Laminaria sp., 39.1% of Agarum cribrosum and 100.0% (n = 3) of mixed macroalgae occur within the same acoustically distinct group of classes. These are discrete from the 8.1% recorded agreement with absences and nulls (>40 m) of macrophytes (n = 32) from a total of 86 ground-truth locations. The results of the water column data extraction (WCDE) show similar success, accurately predicting 78.3% of Laminaria sp. and 30.4% of A. cribrosum observations.     

Response of juvenile lemon sharks, Negaprion brevirostris, to a magnetic barrier simulating a beach net

Beach nets are preventative devices used to minimize interactions between potentially harmful sharks and unsuspecting swimmers. Quantitative studies demonstrated that beach nets drastically reduced local elasmobranch populations, as well as caused considerable bycatch mortality. For this experiment, a beach net-like device was constructed and the behaviors of six juvenile lemon sharks (Negaprion brevirostris) were analyzed. Induced by olfactory and gustatory cues, sharks were given the choice to swim through a magnetic or control opening in the net. In the first trial, all six sharks avoided the magnetic region and significantly preferred to swim through the control region of the fence. The magnetic stimulus no longer affected the swimming behavior of three sharks retested after resting 24 h. Results from the retested sharks were correlated with those from repeated tonic immobility trials, which demonstrated a linear decrease in sensitivity to repeated magnetic stimulation. This study serves as a baseline experiment demonstrating that permanent magnets may substantially decrease elasmobranch mortality within beach nets.     

Discretizing the fracture-matrix interface to simulate solute transport

This article examines the required spatial discretization perpendicular to the fracture-matrix interface (FMI) for numerical simulation of solute transport in discretely fractured porous media. The discrete-fracture, finite-element model HydroGeoSphere ( Therrien et al. 2005 ) and a discrete-fracture implementation of MT3DMS ( Zheng 1990 ) were used to model solute transport in a single fracture, and the results were compared to the analytical solution of Tang et al. (1981) . To match analytical results on the relatively short timescales simulated in this study, very fine grid spacing perpendicular to the FMI of the scale of the fracture aperture is necessary if advection and/or dispersion in the fracture is high compared to diffusion in the matrix. The requirement of such extremely fine spatial discretization has not been previously reported in the literature. In cases of high matrix diffusion, matching the analytical results is achieved with larger grid spacing at the FMI. Cases where matrix diffusion is lower can employ a larger grid multiplier moving away from the FMI. The very fine spatial discretization identified in this study for cases of low matrix diffusion may limit the applicability of numerical discrete-fracture models in such cases.