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.     

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.     

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.     

Submergence and Herbivory as Divergent Causes of Marsh Loss in Long Island Sound

Tidal marsh degradation has been attributed to a number of different causes, but few studies have examined multiple potential factors at the same sites. Differentiating the diverse drivers of marsh loss is critical to prescribing successful interventions for conservation and restoration of this important habitat. We evaluated two hypotheses for vegetation loss at two marshes in Long Island Sound (LIS): (1) marsh submergence, caused by an imbalance between sea-level rise and marsh accretion, and (2) defoliation associated with herbivory by the purple marsh crab, Sesarma reticulatum. At our western LIS site, we found no evidence of herbivory: crabs were scarce, and crab-exclusion cages provided no benefit. We attribute degradation at that site to submergence, a conclusion supported by topographic and hydrologic data showing that loss of vegetation occurred only in wetter parts of the marsh. In contrast, at our central LIS site, our observations were consistent with herbivory as a driving force: There were substantial populations of Sesarma, crab-exclusion cages allowed plants to thrive, and vegetation loss took place across a variety of elevations. We also analyzed soil conditions at both sites, in order to determine the signatures of different degradation processes and assess the potential for restoration. At the submergence site, unvegetated soils exhibited high bulk density, low organic content, and low soil strength, posing significant biogeochemical challenges to re-colonization by vegetation. At the herbivory site, unvegetated soils had a characteristic “riddled-peat” appearance, resulting from expansion and erosion of Sesarma burrow networks. The high redox potential and organic content of those soils suggested that revegetation at the herbivory site would be likely if Sesarma populations could be controlled before erosion leads to elevation loss.     

Destructive near-fault strong ground motion from the 2016 Kumamoto prefecture,Japan, M7.3 earthquake

The sequence of the 2016 Kumamoto, Japan, earthquake, which included an initial M6.5 foreshock on April 14, followed by a larger M7.3 mainshock on April 16, and subsequently occurred high aftershock activity, caused significant damage in Kumamoto and neighboring regions. The near-field strong motion record by strong motion network (K-NET and KiK-net) and the intensity meter network demonstrated clearly the characteristics of the strong ground motion developed by the shallow (H = 12 km), inland earthquake comprising short-time duration (<15–20 s) but large (>1G) ground accelerations. The velocity response spectra of the near-fault motion at Mashiki and Nishihara showed large levels (>300–550 cm/s) in the short-period range (T = 1–2 s), several times larger than that of the near-field record of the destructive 1995 Kobe earthquake (M7.3) and that of the 2004 Mid-Niigata earthquake (M6.8). This period corresponds to the collapse vulnerability of Japanese wooden-frame houses, and is the major cause of severe damage during the Kumamoto earthquake. The response spectra also showed extremely large levels (>240–340 cm/s) in the long-period (T > 3 s) band, which is potentially disastrous for high-rise buildings, large oil storage tanks, etc. to have longer resonant period. Such long-period motion was, for the most parts, developed by the static displacement of the fault movement rather than by the seismic waves radiating from the source fault. Thus, the extreme near-fault long-period motion was hazardous only close to the fault but it attenuated very rapidly away from the fault.