Spatial and temporal constraints on regional-scale groundwater flow in the Pampa del Tamarugal Basin,Atacama Desert,Chile

Aquifers within the Pampa del Tamarugal Basin (Atacama Desert, northern Chile) are the sole source of water for the coastal city of Iquique and the economically important mining industry. Despite this, the regional groundwater system remains poorly understood. Although it is widely accepted that aquifer recharge originates as precipitation in the Altiplano and Andean Cordillera to the east, there remains debate on whether recharge is driven primarily by near-surface groundwater flow in response to periodic flood events or by basal groundwater flux through deep-seated basin fractures. In addressing this debate, the present study quantifies spatial and temporal variability in regional-scale groundwater flow paths at 20.5°S latitude by combining a two-dimensional model of groundwater and heat flow with field observations and δ¹⁸O isotope values in surface water and groundwater. Results suggest that both previously proposed aquifer recharge mechanisms are likely influencing aquifers within the Pampa del Tamarugal Basin; however, each mechanism is operating on different spatial and temporal scales. Storm-driven flood events in the Altiplano readily transmit groundwater to the eastern Pampa del Tamarugal Basin through near-surface groundwater flow on short time scales, e.g., 10⁰–10¹ years, but these effects are likely isolated to aquifers in the eastern third of the basin. In addition, this study illustrates a physical mechanism for groundwater originating in the eastern highlands to recharge aquifers and salars in the western Pampa del Tamarugal Basin over timescales of 10⁴–10⁵ years.     

Testing the suitability of geologic frameworks for extrapolating hydraulic properties across regional scales

The suitability of geologic frameworks for extrapolating hydraulic conductivity (K) to length scales commensurate with hydraulic data is difficult to assess. A novel method is presented for evaluating assumed relations between K and geologic interpretations for regional-scale groundwater modeling. The approach relies on simultaneous interpretation of multiple aquifer tests using alternative geologic frameworks of variable complexity, where each framework is incorporated as prior information that assumes homogeneous K within each model unit. This approach is tested at Pahute Mesa within the Nevada National Security Site (USA), where observed drawdowns from eight aquifer tests in complex, highly faulted volcanic rocks provide the necessary hydraulic constraints. The investigated volume encompasses 40 mi³ (167 km³) where drawdowns traversed major fault structures and were detected more than 2 mi (3.2 km) from pumping wells. Complexity of the five frameworks assessed ranges from an undifferentiated mass of rock with a single unit to 14 distinct geologic units. Results show that only four geologic units can be justified as hydraulically unique for this location. The approach qualitatively evaluates the consistency of hydraulic property estimates within extents of investigation and effects of geologic frameworks on extrapolation. Distributions of transmissivity are similar within the investigated extents irrespective of the geologic framework. In contrast, the extrapolation of hydraulic properties beyond the volume investigated with interfering aquifer tests is strongly affected by the complexity of a given framework. Testing at Pahute Mesa illustrates how this method can be employed to determine the appropriate level of geologic complexity for large-scale groundwater modeling.     

Double trouble: subsidence and CO<Subscript>2</Subscript> respiration due to 1,000 years of Dutch coastal peatlands cultivation

Coastal plains are amongst the most densely populated areas in the world. Many coastal peatlands are drained to create arable land. This is not without consequences; physical compaction of peat and its degradation by oxidation lead to subsidence, and oxidation also leads to emissions of carbon dioxide (CO₂). This study complements existing studies by quantifying total land subsidence and associated CO₂ respiration over the past millennium in the Dutch coastal peatlands, to gain insight into the consequences of cultivating coastal peatlands over longer timescales. Results show that the peat volume loss was 19.8 km³, which lowered the Dutch coastal plain by 1.9 m on average, bringing most of it below sea level. At least 66 % of the volume reduction is the result of drainage, and 34 % was caused by the excavation and subsequent combustion of peat. The associated CO₂ respiration is equivalent to a global atmospheric CO₂ concentration increase of ~0.39 ppmv. Cultivation of coastal peatlands can turn a carbon sink into a carbon source. If the path taken by the Dutch would be followed worldwide, there will be double trouble: globally significant carbon emissions and increased flood risk in a globally important human habitat. The effects would be larger than the historic ones because most of the cumulative Dutch subsidence and peat loss was accomplished with much less efficient techniques than those available now.     

Using a Thermal Proxy to Examine Sediment–Water Exchange in Mid-Continental Shelf Sandy Sediments

Fluid exchange across the sediment–water interface in a sandy open continental shelf setting was studied using heat as a tracer. Summertime tidal oscillation of cross-shelf thermal fronts on the South Atlantic Bight provided a sufficient signal at the sediment–water interface to trace the advective and conductive transport of heat into and out of the seabed, indicating rapid flushing of ocean water through the upper 10–40 cm of the sandy seafloor. A newly developed transport model was applied to the in situ temperature data set to estimate the extent to which heat was transported by advection rather than conduction. Heat transported by shallow 3-D porewater flow processes was accounted for in the model by using a dispersion term, the depth and intensity of which reflected the depth and intensity of shallow flushing. Similar to the results of past studies in shallower and more energetic nearshore settings, transport of heat was greater when higher near-bed velocities and shear stresses occurred over a rippled bed. However, boundary layer processes by themselves were insufficient to promote non-conductive heat transport. Advective heat transport only occurred when both larger boundary layer stresses and thermal instabilities within the porespace were present. The latter process is dependent on shelf-scale heating and cooling of bottom water associated with upwelling events that are not coupled to local-scale boundary layer processes.     

A multiobjective steepest descent method with applications to optimal well control

Multiobjective optimization deals with mathematical optimization problems where two or more objective functions (cost functions) are to be optimized (maximized or minimized) simultaneously. In most cases of interest, the objective functions are in conflict, i.e., there does not exist a decision (design) vector (vector of optimization variables) at which every objective function takes on its optimal value. The solution of a multiobjective problem is commonly defined as a Pareto front, and any decision vector which maps to a point on the Pareto front is said to be Pareto optimal. We present an original derivation of an analytical expression for the steepest descent direction for multiobjective optimization for the case of two objectives. This leads to an algorithm which can be applied to obtain Pareto optimal points or, equivalently, points on the Pareto front when the problem is the minimization of two conflicting objectives. The method is in effect a generalization of the steepest descent algorithm for minimizing a single objective function. The steepest-descent multiobjective optimization algorithm is applied to obtain optimal well controls for two example problems where the two conflicting objectives are the maximization of the life-cycle (long-term) net-present-value (NPV) and the maximization of the short-term NPV. The results strongly suggest the multiobjective steepest-descent (MOSD) algorithm is more efficient than competing multiobjective optimization algorithms.     

Heat tracing to examine seasonal groundwater flow beneath a low-gradient stream in rural central Illinois,USA

The thermal profile of a streambed is affected by a number of factors including: temperatures of stream water and groundwater, hydraulic conductivity, thermal conductivity, heat capacity of the streambed, and the geometry of hyporheic flow paths. Changes in these parameters over time cause changes in thermal profiles. In this study, temperature data were collected at depths of 30, 60, 90 and 150 cm at six streambed wells 5 m apart along the thalweg of Little Kickapoo Creek, in rural central Illinois, USA. This is a third-order low-gradient baseflow-fed stream. A positive temperature gradient with inflection at 90-cm depth was observed during the summer period. A negative temperature gradient with inflection at 30 cm was observed during the winter period, which suggests greater influence of stream-water temperatures in the substrate during the summer. Thermal models of the streambed were built using VS2DHI to simulate the thermal profiles observed in the field. Comparison of the parameters along with analysis of temperature envelopes and Peclet numbers suggested greater upwelling and stability in temperatures during the winter than during the summer. Upwelling was more pronounced in the downstream reach of the pool in the riffle and pool sequence.     

The Influence of Bioturbation on Iron and Sulphur Cycling in Marine Sediments: A Model Analysis

The geochemical cycles of iron and sulphur in marine sediments are strongly intertwined and give rise to a complex network of redox and precipitation reactions. Bioturbation refers to all modes of transport of particles and solutes induced by larger organisms, and in the present-day seafloor, bioturbation is one of the most important factors controlling the biogeochemical cycling of iron and sulphur. To better understand how bioturbation controls Fe and S cycling, we developed reactive transport model of a coastal sediment impacted by faunal activity. Subsequently, we performed a model sensitivity analysis, separately investigating the two different transport modes of bioturbation, i.e. bio-mixing (solid particle transport) and bio-irrigation (enhanced solute transport). This analysis reveals that bio-mixing and bio-irrigation have distinct—and largely opposing effects on both the iron and sulphur cycles. Bio-mixing enhances transport between the oxic and suboxic zones, thus promoting the reduction of oxidised species (e.g. iron oxyhydroxides) and the oxidation of reduced species (e.g. iron sulphides). Through the re-oxidation of iron sulphides, bio-mixing strongly enhances the recycling of Fe and S between their reduced and oxidised states. Bio-irrigation on the other hand removes reduced solutes, i.e. ferrous iron and free sulphide, from the sediment pore water. These reduced species are then reoxidised in the overlying water and not recycled within the sediment column, which leads to a decrease in Fe and S recycling. Overall, our results demonstrate that the ecology of the macrofauna (inducing bio-mixing or bio-irrigation, or both) matters when assessing their impact on sediment geochemistry. This finding seems particularly relevant for sedimentary cycling across Cambrian transition, when benthic fauna started colonizing and reworking the seafloor.     

Assessment of Submarine Groundwater Discharge (SGD) as a Source of Dissolved Radium and Nutrients to Moorea (French Polynesia) Coastal Waters

Previous work has documented large fluxes of freshwater and nutrients from submarine groundwater discharge (SGD) into the coastal waters of a few volcanic oceanic islands. However, on the majority of such islands, including Moorea (French Polynesia), SGD has not been studied. In this study, we used radium (Ra) isotopes and salinity to investigate SGD and associated nutrient inputs at five coastal sites and Paopao Bay on the north shore of Moorea. Ra activities were highest in coastal groundwater, intermediate in coastal ocean surface water, and lowest in offshore surface water, indicating that high-Ra groundwater was discharging into the coastal ocean. On average, groundwater nitrate and nitrite (N + N), phosphate, ammonium, and silica concentrations were 12, 21, 29, and 33 times greater, respectively, than those in coastal ocean surface water, suggesting that groundwater discharge could be an important source of nutrients to the coastal ocean. Ra and salinity mass balances indicated that most or all SGD at these sites was saline and likely originated from a deeper, unsampled layer of Ra-enriched recirculated seawater. This high-salinity SGD may be less affected by terrestrial nutrient sources, such as fertilizer, sewage, and animal waste, compared to meteoric groundwater; however, nutrient-salinity trends indicate it may still have much higher concentrations of nitrate and phosphate than coastal receiving waters. Coastal ocean nutrient concentrations were virtually identical to those measured offshore, suggesting that nutrient subsidies from SGD are efficiently utilized.