Ground Water Replenishment with Reclaimed Water in the City of Los Angeles

A two-year pilot study involving the recharge of a ground water basin with reclaimed water was completed in the city of Los Angeles. The city's Department of Water and Power is planning to initiate several ground water recharge projects using reclaimed water in the near future. One such project is the Headworks Recharge Project, the focus of this paper, Approximately 1 cfs of the Los Angeles (LA) River water comprised of 70% tertiary treated reclaimed water was recharged on a two-day wet and five-day dry cycle. The recharge water was then extracted from the basin approximately 1000 feet downgradient. Results showed greater than 4-log removal of coliform bacteria, up to 87% reduction in TOC, and compliance of the product water with federal and state drinking water standards. Model simulation showed after 15 years of recharging 3000 acre-feet per year of the LA River water and extracting about 10,000 acre-feet from the basin, the product water would contain from 5% to 15% reclaimed water. This is well below the maximum allowable limit of 20% stipulated by the California regulation.     

The effects of uranium speciation on the rate of U(VI) reduction by Shewanella oneidensis MR-1

We measured the kinetics of U(VI) reduction by Shewanella oneidensis MR-1 under anaerobic conditions in the presence of variable concentrations of either EDTA or dissolved Ca. We measured both total dissolved U and U(VI) concentrations in solution as a function of time. In separate experiments, we also measured the extent of U(VI) adsorption onto S. oneidensis in order to quantify the thermodynamic stabilities of the important U(VI)-bacterial surface complexes. In the EDTA experiments, the rate of U(IV) production increased with increasing EDTA concentration. However, the total dissolved U concentrations remained constant and identical to the initial U concentrations during the course of the experiments for all EDTA-bearing systems. Additionally, the U(VI) reduction rate in the EDTA experiments exhibited a strong correlation to the concentration of the aqueous U⁴⁺-EDTA complex. We conclude that the U(VI) reduction rate increases with increasing EDTA concentration, likely due to U⁴⁺-EDTA aqueous complexation which removes U(IV) from the cell surface and prevents UO₂ precipitation.In the Ca experiments, the U(VI) reduction rate decreased as Ca concentration increased. Our thermodynamic modeling results based on the U(VI) adsorption data demonstrate that U(VI) was adsorbed onto the bacterial surface in the form of a Ca-uranyl-carbonate complex in addition to a number of other Ca-free uranyl complexes. The observed U(VI) reduction rates in the presence of Ca exhibit a strong negative correlation to the concentration of the Ca-uranyl-carbonate bacterial surface complex, but a strong positive correlation to the total concentration of all the other Ca-free uranyl surface complexes. Thus, the concentration of these Ca-free uranyl surface complexes appears to control the rate of U(VI) reduction by S. oneidensis in the presence of dissolved Ca. Our results demonstrate that U speciation, both of U(VI) before reduction and of U(IV) after reduction, affects the reduction kinetics, and that thermodynamic modeling of the U speciation may be useful in the prediction of reduction kinetics in realistic geologic settings.     

The deglaciation of the western sector of the Irish Ice Sheet from the inner continental shelf to its terrestrial margin

This paper provides a new deglacial chronology for retreat of the Irish Ice Sheet from the continental shelf of western Ireland to the adjoining coastline, a region where the timing and drivers of ice recession have never been fully constrained. Previous work suggests maximum ice-sheet extent on the outer western continental shelf occurred at ~26–24 cal. ka BP with the initial retreat of the ice marked by the production of grounding-zone wedges between 23–21.1 cal. ka BP. However, the timing and rate of ice-sheet retreat from the inner continental shelf to the present coast are largely unknown. This paper reports 31 new terrestrial cosmogenic nuclide (TCN) ages from erratics and ice-moulded bedrock and three new optically stimulated luminescence (OSL) ages on deglacial outwash. The TCN data constrain deglaciation of the near coast (Aran Islands) to ~19.5–18.5 ka. This infers ice retreated rapidly from the mid-shelf after 21 ka, but the combined effects of bathymetric shallowing and pinning acted to stabilize the ice at the Aran Islands. However, marginal stability was short-lived, with multiple coastal sites along the Connemara/Galway coasts demonstrating ice recession under terrestrial conditions by 18.2–17. ka. This pattern of retreat continued as ice retreated eastward through inner Galway Bay by 16.5 ka. South of Galway, the Kilkee–Kilrush Moraine Complex and Scattery Island moraines point to late stage re-advances of the ice sheet into southern County Clare ~14.1–13.3 ka, but the large errors associated with the OSL ages make correlation with other regional re-advances difficult. It seems more likely that these moraines are the product of regional ice lobes adjusting to internal ice-sheet dynamics during deglaciation in the time window 17–16 ka.     

The Role of Ocean Tides on Groundwater-Surface Water Exchange in a Mangrove-Dominated Estuary: Shark River Slough,Florida Coastal Everglades,USA

Low-relief environments like the Florida Coastal Everglades (FCE) have complicated hydrologic systems where surface water and groundwater processes are intimately linked yet hard to separate. Fluid exchange within these low-hydraulic-gradient systems can occur across broad spatial and temporal scales, with variable contributions to material transport and transformation. Identifying and assessing the scales at which these processes operate is essential for accurate evaluations of how these systems contribute to global biogeochemical cycles. The distribution of ²²²Rn and ^223,224,226Ra have complex spatial patterns along the Shark River Slough estuary (SRSE), Everglades, FL. High-resolution time-series measurements of ²²²Rn activity, salinity, and water level were used to quantify processes affecting radon fluxes out of the mangrove forest over a tidal cycle. Based on field data, tidal pumping through an extensive network of crab burrows in the lower FCE provides the best explanation for the high radon and fluid fluxes. Burrows are irrigated during rising tides when radon and other dissolved constituents are released from the mangrove soil. Flushing efficiency of the burrows—defined as the tidal volume divided by the volume of burrows—estimated for the creek drainage area vary seasonally from 25 (wet season) to 100 % (dry season) in this study. The tidal pumping of the mangrove forest soil acts as a significant vector for exchange between the forest and the estuary. Processes that enhance exchange of O₂ and other materials across the sediment-water interface could have a profound impact on the environmental response to larger scale processes such as sea level rise and climate change. Compounding the material budgets of the SRSE are additional inputs from groundwater from the Biscayne Aquifer, which were identified using radium isotopes. Quantification of the deep groundwater component is not obtainable, but isotopic data suggest a more prevalent signal in the dry season. These findings highlight the important role that both tidal- and seasonal-scale forcings play on groundwater movement in low-gradient hydrologic systems.     

Spatial variability assessment of soil nutrients in an intense agricultural area,a case study of Rugao County in Yangtze River Delta Region,China

Topsoil samples (0–20 cm) (n = 237) were collected from Rugao County, China. Geostatistical variogram analysis, sequential Gaussian simulation (SGS), and principal component (PC) analysis were applied to assess spatial variability of soil nutrients, identify the possible areas of nutrient deficiency, and explore spatial scale of variability of soil nutrients in the county. High variability of soil nutrient such as soil organic matter (SOM), total nitrogen (TN), available P, K, Fe, Mn, Cu, Zn, and B concentrations were observed. Soil nutrient properties displayed significant differences in their spatial structures, with available Cu having strong spatial dependence, SOM and available P having weak spatial dependence, and other nutrient properties having moderate spatial dependence. The soil nutrient deficiency, defined here as measured nutrient concentrations which do not meet the advisory threshold values specific to the county for dominant crops, namely rice, wheat, and rape seeds, was observed in available K and Zn, and the deficient areas covered 38 and 11%, respectively. The first three PCs of the nine soil nutrient properties explained 62.40% of the total variance. TN and SOM with higher loadings on PC1 are closely related to soil texture derived from different parent materials. The PC2 combined intermediate response variables such as available Zn and P that are likely to be controlled by land use and soil pH. Available B has the highest loading on PC3 and its variability of concentrations may be primarily ascribed to localized anthropogenic influence. The amelioration of soil physical properties (i.e. soil texture) and soil pH may improve the availability of soil nutrients and the sustainability of the agricultural system of Rugao County.     

Formation of syngenetic and early diagenetic iron minerals in the late Archean Mt. McRae Shale, Hamersley Basin, Australia: New insights on the patterns, controls and paleoenvironmental implications of authigenic mineral formation

A section through the late Archean Mt. McRae Shale comprising, in ascending order, a lower shale interval (LSI), a banded iron formation (BIF), an upper shale (USI) and a carbonate (C1) has been analyzed for total Fe and Al contents and authigenic Fe present as carbonate, oxide, sulfide and silicate phases. The authigenic mineralogy is controlled by the episodic addition of Fe from hydrothermal activity and removal of Fe by sulfide, relative to rates of clastic sedimentation. The LSI and BIF have mean Fe_T/Al values of 2 and 5, respectively, that record iron enrichment from hydrothermal sources. Iron was precipitated primarily as siderite accompanied by Fe-rich chlorite from anoxic bottom waters rich in dissolved Fe. Pyrite formation was probably limited by the availability of sulfate, which was present at low concentrations and became rapidly depleted. The USI has generally lower Fe_T/Al values (0.6-1.3), similar to those found in Paleozoic shales, with the exception of one interval where enrichment may reflect either a weak hydrothermal source or the operation of an iron shuttle. This interval contains authigenic Fe predominantly as pyrite, where high values for DOP (>0.8) indicate the existence of a water column that became rich in dissolved sulfide (euxinic) when sulfate concentrations increased due to a transient or secular increase in ocean/atmosphere oxygenation. High concentrations of dissolved sulfide maintained low concentrations of dissolved Fe, which allowed only minor amounts of Fe to be precipitated as carbonates and silicates. The USI also has elevated concentrations of organic matter that most probably reflect increased productivity and likely limited euxinia to midportions of the water column on the basin margin. The carbonate C1 represents a basinal chemistry where sulfide has been removed and Fe_T/Al values are ∼1 indicating that hydrothermal activity again produced dissolved Fe-rich bottom waters. Detailed iron speciation of the Mt. McRae Shale can be used to recognize spatial and temporal variations in iron and sulfur inputs to the late Archean Hamersley Basin, just prior to the Paleoproterozoic rise in atmospheric oxygenation, and our refined methods have relevance to all Fe-rich deposits.     

Assessment of the UV camera sulfur dioxide retrieval for point source plumes

Digital cameras, sensitive to specific regions of the ultra-violet (UV) spectrum, have been employed for quantifying sulfur dioxide (SO₂) emissions in recent years. The instruments make use of the selective absorption of UV light by SO₂ molecules to determine pathlength concentration. Many monitoring advantages are gained by using this technique, but the accuracy and limitations have not been thoroughly investigated. The effect of some user-controlled parameters, including image exposure duration, the diameter of the lens aperture, the frequency of calibration cell imaging, and the use of the single or paired bandpass filters, have not yet been addressed. In order to clarify methodological consequences and quantify accuracy, laboratory and field experiments were conducted. Images were collected of calibration cells under varying observational conditions, and our conclusions provide guidance for enhanced image collection. Results indicate that the calibration cell response is reliably linear below 1500 ppm m, but that the response is significantly affected by changing light conditions. Exposure durations that produced maximum image digital numbers above 32 500 counts can reduce noise in plume images. Sulfur dioxide retrieval results from a coal-fired power plant plume were compared to direct sampling measurements and the results indicate that the accuracy of the UV camera retrieval method is within the range of current spectrometric methods.     

Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (V_p and V_s, respectively) decrease with increasing porosity and that V_p decreases approximately twice as fast as V_s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model.     

The dynamics of mountain erosion: Cirque growth slows as landscapes age

Glacial cirques are widely used palaeoenvironmental indicators, and are key to understanding the role of glaciers in shaping mountain topography. However, notable uncertainty persists regarding the rate and timing of cirque erosion. In order to address this uncertainty, we analyse the dimensions of 2208 cirques in Britain and Ireland and model ice accumulation to investigate the degree of coupling between glacier occupation times and cirque growth. Results indicate that during the last ~120 ka, cirques were glacier-free for an average of 52.0 ± 21.2 ka (43 ± 18%); occupied by small (largely cirque-confined) glaciers for 16.2 ± 9.9 ka (14 ± 8%); and occupied by large glaciers, including ice sheets, for 51.8 ± 18.6 ka (43 ± 16%). Over the entire Quaternary (i.e. 2.6 Ma), we estimate that cirques were glacier-free for 1.1 ± 0.5 Ma; occupied by small glaciers for 0.3 ± 0.2 Ma; and occupied by large glaciers for 1.1 ± 0.4 Ma. Comparing occupation times to cirque depths, and calculating required erosion rates, reveals that continuous cirque growth during glacier occupation is unlikely. Instead, we propose that cirques attained much of their size during the first occupation of a non-glacially sculpted landscape (perhaps during the timeframe of a single glacial cycle). During subsequent glacier occupations, cirque growth may have slowed considerably, with the highest rates of subglacial erosion focused during periods of marginal (small glacier) glaciation. We propose comparatively slow rates of growth following initial cirque development because a ‘least resistance’ shape is formed, and as cirques deepen, sediment becomes trapped subglacially, partly protecting the bedrock from subsequent erosion. In support of the idea of rapid cirque growth, we present evidence from northern British Columbia, where cirques of comparable size to those in Britain and Ireland developed in less than 140 ka. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Sensitivity of Simulated Hyporheic Exchange to River Bathymetry: The Steinlach River Test Site

This study determines the aspects of river bathymetry that have the greatest influence on the predictive biases when simulating hyporheic exchange. To investigate this, we build a highly parameterized HydroGeoSphere model of the Steinlach River Test Site in southwest Germany as a reference. This model is then modified with simpler bathymetries, evaluating the changes to hyporheic exchange fluxes and transit time distributions. Results indicate that simulating hyporheic exchange with a high-resolution detailed bathymetry using a three-dimensional fully coupled model leads to nested multiscale hyporheic exchange systems. A poorly resolved bathymetry will underestimate the small-scale hyporheic exchange, biasing the simulated hyporheic exchange towards larger scales, thus leading to overestimates of hyporheic exchange residence times. This can lead to gross biases in the estimation of a catchment's capacity to attenuate pollutants when extrapolated to account for all meanders along an entire river within a watershed. The detailed river slope alone is not enough to accurately simulate the locations and magnitudes of losing and gaining river reaches. Thus, local bedforms in terms of bathymetric highs and lows within the river are required. Bathymetry surveying campaigns can be more effective by prioritizing bathymetry measurements along the thalweg and gegenweg of a meandering channel. We define the gegenweg as the line that connects the shallowest points in successive cross-sections along a river opposite to the thalweg under average flow conditions. Incorporating local bedforms will likely capture the nested nature of hyporheic exchange, leading to more physically meaningful simulations of hyporheic exchange fluxes and transit times.