Effects of nonlinear model response on allocation of streamflow depletion: exemplified by the case of Beaver Creek,USA

Anomalies found when apportioning responsibility for streamflow depletion are examined. The anomalies arise when responsibility is assigned to the two states that contribute to depletion of Beaver Creek in the Republican River Basin in the United States. The apportioning procedure for this basin presumes that the sum of streamflow depletions, computed by comparing simulation model runs with and without groundwater pumping from individual states, approximates the streamflow depletion when both states are pumping. In the case study presented here, this presumed superposition fails dramatically. The stream drying and aquifer-storage depletion, as represented in the simulation model used for allocation, are examined in detail to understand the hydrologic and numerical basis for the severe nonlinear response. Users of apportioning procedures that rely on superposition should be aware of the presence and likely magnitude of nonlinear responses in modeling tools.     

Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria,Austria

An approach to establish the recharge component of managed aquifer recharge (MAR) has recently been proposed that uses small-diameter shallow wells installed using relatively inexpensive drilling methods such as direct push. As part of further development of that approach, a generalized procedure is presented for a technical and economic assessment of the approach’s potential in comparison to other systems. Following this procedure, the use of small-diameter wells was evaluated both experimentally and numerically for a site located in southern Styria, Austria. MAR is currently done at the site using a horizontal pipe infiltration system, and system expansion has been proposed with a target rate of 12 l/s using small-diameter wells as one possible option. A short-duration single-well field recharge experiment (recharge rate 1.3–3.5 l/s) was performed (recharge by gravity only). Numerical modeling of the injection test was used to estimate hydraulic conductivity (K). Quasi-steady-state, single-well recharge simulations for different locations, as well as a long-term transient simulation, were performed using the K value calibrated from the field injection test. Results indicate that a recharge capacity of 4.1 l/s was achievable with a maximum head rise of 0.2 m at the injection well. Finally, simulations were performed for three different well fields (4, 6 and 8 wells, respectively) designed to infiltrate a target rate of 12 l/s. The experimental and numerical assessments, supported by a cost analysis of the small-diameter wells, indicate that the small-diameter wells are a viable, cost-effective recharge approach at this and other similar sites.     

A multi-tracer approach to delineate groundwater dynamics in the Rio Actopan Basin,Veracruz State,Mexico

Geochemistry and environmental tracers were used to understand groundwater resources, recharge processes, and potential sources of contamination in the Rio Actopan Basin, Veracruz State, Mexico. Total dissolved solids are lower in wells and springs located in the basin uplands compared with those closer to the coast, likely associated with rock/water interaction. Geochemical results also indicate some saltwater intrusion near the coast and increased nitrate near urban centers. Stable isotopes show that precipitation is the source of recharge to the groundwater system. Interestingly, some high-elevation springs are more isotopically enriched than average annual precipitation at higher elevations, indicating preferential recharge during the drier but cooler winter months when evapotranspiration is reduced. In contrast, groundwater below 1,200 m elevation is more isotopically depleted than average precipitation, indicating recharge occurring at much higher elevation than the sampling site. Relatively cool recharge temperatures, derived from noble gas measurements at four sites (11–20 °C), also suggest higher elevation recharge. Environmental tracers indicate that groundwater residence time in the basin ranges from 12,000 years to modern. While this large range shows varying groundwater flowpaths and travel times, ages using different tracer methods (¹⁴C, ³H/³He, CFCs) were generally consistent. Comparing multiple tracers such as CFC-12 with CFC-113 indicates piston-flow to some discharge points, yet binary mixing of young and older groundwater at other points. In summary, groundwater within the Rio Actopan Basin watershed is relatively young (Holocene) and the majority of recharge occurs in the basin uplands and moves towards the coast.     

Hydrologic testing during drilling: application of the flowing fluid electrical conductivity (FFEC) logging method to drilling of a deep borehole

Drilling of a deep borehole does not normally allow for hydrologic testing during the drilling period. It is only done when drilling experiences a large loss (or high return) of drilling fluid due to penetration of a large-transmissivity zone. The paper proposes the possibility of conducting flowing fluid electrical conductivity (FFEC) logging during the drilling period, with negligible impact on the drilling schedule, yet providing important information on depth locations of both high- and low-transmissivity zones and their hydraulic properties. The information can be used to guide downhole fluid sampling and post-drilling detailed testing of the borehole. The method has been applied to the drilling of a 2,500-m borehole at Åre, central Sweden, firstly when the drilling reached 1,600 m, and then when the drilling reached the target depth of 2,500 m. Results unveil eight hydraulically active zones from 300 m down to borehole bottom, with depths determined to within the order of a meter. Further, the first set of data allows the estimation of hydraulic transmissivity values of the six hydraulically conductive zones found from 300 to 1,600 m, which are very low and range over one order of magnitude.     

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.     

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.     

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.