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.     

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.     

Improved regional groundwater flow modeling using drainage features: a case study of the central northern karst aquifer system of Puerto Rico (USA)

In northern Puerto Rico (USA), subsurface conduit networks with unknown characteristics, and surface features such as springs, rivers, lagoons and wetlands, drain the coastal karst aquifers. In this study, drain lines connecting sinkholes and springs are used to improve the developed regional model by simulating the drainage effects of conduit networks. Implemented in an equivalent porous media (EPM) approach, the model with drains is able to roughly reproduce the spring discharge hydrographs in response to rainfall. Hydraulic conductivities are found to be scale dependent and significantly increase with higher test radius, indicating scale dependency of the EPM approach. Similar to other karst regions in the world, hydraulic gradients are steeper where the transmissivity is lower approaching the coastline. This study enhances current understanding of the complex flow patterns in karst aquifers and suggests that using a drainage feature improves modeling results where available data on conduit characteristics are minimal.     

Identification of hydrochemical facies in the Roswell Artesian Basin,New Mexico (USA), using graphical and statistical methods

Analysis of groundwater chemistry can yield important insights about subsurface conditions, and provide an alternative and complementary method for characterizing basin hydrogeology, especially in areas where hydraulic data are limited. More specifically, hydrochemical facies have been used for decades to help understand basin flow and transport, and a set of facies were developed for the Roswell Artesian Basin (RAB) in a semi-arid part of New Mexico, USA. The RAB is an important agricultural water source, and is an excellent example of a rechargeable artesian system. However, substantial uncertainties about the RAB hydrogeology and groundwater chemistry exist. The RAB was a great opportunity to explore hydrochemcial facies definition. A set of facies, derived from fingerprint diagrams (graphical approach), existed as a basis for testing and for comparison to principal components, factor analysis, and cluster analyses (statistical approaches). Geochemical data from over 300 RAB wells in the central basin were examined. The statistical testing of fingerprint-diagram-based facies was useful in terms of quantitatively evaluating differences between facies, and for understanding potential controls on basin groundwater chemistry. This study suggests the presence of three hydrochemical facies in the shallower part of the RAB (mostly unconfined conditions) and three in the deeper artesian system of the RAB. These facies reflect significant spatial differences in chemistry in the basin that are associated with specific stratigraphic intervals as well as structural features. Substantial chemical variability across faults and within fault blocks was also observed.     

The removal kinetics of dissolved organic matter and the optical clarity of groundwater

Concentrations of dissolved organic matter (DOM) and ultraviolet/visible light absorbance decrease systematically as groundwater moves through the unsaturated zones overlying aquifers and along flowpaths within aquifers. These changes occur over distances of tens of meters (m) implying rapid removal kinetics of the chromophoric DOM that imparts color to groundwater. A one-compartment input-output model was used to derive a differential equation describing the removal of DOM from the dissolved phase due to the combined effects of biodegradation and sorption. The general solution to the equation was parameterized using a 2-year record of dissolved organic carbon (DOC) concentration changes in groundwater at a long-term observation well. Estimated rates of DOC loss were rapid and ranged from 0.093 to 0.21 micromoles per liter per day (μM d^?1), and rate constants for DOC removal ranged from 0.0021 to 0.011 per day (d^?1). Applying these removal rate constants to an advective-dispersion model illustrates substantial depletion of DOC over flow-path distances of 200 m or less and in timeframes of 2 years or less. These results explain the low to moderate DOC concentrations (20–75 μM; 0.26–1 mg L^?1) and ultraviolet absorption coefficient values (a ₂₅₄?<?5 m^?1) observed in groundwater produced from 59 wells tapping eight different aquifer systems of the United States. The nearly uniform optical clarity of groundwater, therefore, results from similarly rapid DOM-removal kinetics exhibited by geologically and hydrologically dissimilar aquifers.     

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.     

Groundwater flow in a transboundary fault-dominated aquifer and the importance of regional modeling: the case of the city of Querétaro,Mexico

The city of Querétaro, located near the political boundary of the Mexican states of Querétaro and Guanajuato, relies on groundwater as it sole water supply. Groundwater extraction in the city increased from 21?×?10⁶ m³/yr in 1970 to 104?×?10⁶ m³/yr in 2010, with an associated drawdown of 100 m in some parts of the aquifer. A three-dimensional numerical groundwater-flow model has been developed that represents the historical evolution of the aquifer’s potentiometric levels and is used to simulate the effect of two scenarios: (1) a 40 % reduction in the extraction rate from public water supply wells in early 2011 (thus reducing the extraction to 62?×?10⁶ m³/yr), and (2) a further reduction in 2021 to 1?×?10⁶ m³/yr. The modeling results project a temporary recovery of the potentiometric levels after the 40 % reduction of early 2011, but a return to 2010 levels by 2020. If scenario 2 is implemented in 2021, the aquifer will take nearly 30 years to recover to the simulated levels of 1995. The model also shows that the wells located in the city of Querétaro started to extract water from part of the aquifer beneath the State of Guanajuato in the late 1970s, thus showing that the administrative boundaries used in Mexico to study and develop water resources are inappropriate, and consideration should be given to physical boundaries instead. A regional approach to studying aquifers is needed in order to adequately understand groundwater flow dynamics.     

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.