The role of baseflow in dissolved solids delivery to streams in the Upper Colorado River Basin

Salinity has a major effect on water users in the Colorado River Basin, estimated to cause almost $300 million per year in economic damages. The Colorado River Basin Salinity Control Program implements and manages projects to reduce salinity loads, investing millions of dollars per year in irrigation upgrades, canal projects, and other mitigation strategies. To inform and improve mitigation efforts, there is a need to better understand sources of salinity to streams and how salinity has changed over time. This study explores salinity in the baseflow fraction of streamflow, assessing whether groundwater is a significant contributor of dissolved solids to streams in the Upper Colorado River Basin (UCRB). Chemical hydrograph separation was used to estimate baseflow discharge and baseflow dissolved solids loads at stream gages (n = 69) across the UCRB. On average, it is estimated that 89% of dissolved solids loads originate from the baseflow fraction of streamflow, indicating that subsurface transport processes play a dominant role in delivering dissolved solids to streams in the UCRB. A statistical trend analysis using weighted regressions on time, discharge, and season was used to evaluate changes in baseflow dissolved solids loads in streams (n = 27) from 1986 to 2011. Decreasing trends in baseflow dissolved solids loads were observed at 63% of streams. At the three most downstream sites, Green River at Green River, UT, Colorado River at Cisco, UT, and the San Juan River near Bluff, UT, baseflow dissolved solids loads decreased by a combined 823,000 metric tons (mT), which is approximately 69% of projected basin‐scale decreases in total dissolved solids loads as a result of salinity control efforts. Decreasing trends in baseflow dissolved solids loads suggest that salinity mitigation projects, landscape changes, and/or climate are reducing dissolved solids transported to streams through the subsurface. Notably, the pace and extent of decreases in baseflow dissolved solids loads declined during the most recent decade; average decreasing loads during the 2000s (28,200 mT) were only 54% of average decreasing loads in the 1990s (51,700 mT).     

Application of Distributed Temperature Sensing for coupled mapping of sedimentation processes and spatio‐temporal variability of groundwater discharge in soft‐bedded streams

The delineation of groundwater discharge areas based on Distributed Temperature Sensing (DTS) data of the streambed can be difficult in soft‐bedded streams where sedimentation and scouring processes constantly change the position of the fibre optic cable relative to the streambed. Deposition‐induced temperature anomalies resemble the signal of groundwater discharge while scouring will cause the cable to float in the water column and measure stream water temperatures. DTS applied in a looped layout with nine fibre optic cable rows in a 70 × 5 m section of a soft‐bedded stream made it possible to detect variability in streambed temperatures between October 2011 and January 2012. Detailed monthly streambed elevation surveys were carried out to monitor the position of the fibre optic cable relative to the streambed and to quantify the effect of sedimentation processes on streambed temperatures. Based on the simultaneous interpretation of streambed temperature and elevation data, a method is proposed to delineate potential high‐groundwater discharge areas and identify deposition‐induced temperature anomalies in soft‐bedded streams. Potential high‐discharge sites were detected using as metrics the daily minimum, maximum and mean streambed temperatures as well as the daily amplitude and standard deviation of temperatures. The identified potential high‐discharge areas were mostly located near the channel banks, also showing temporal variability because of the scouring and redistribution of streambed sediments, leading to the relocation of pool‐riffle sequences. This study also shows that sediment deposits of 0.1 m thickness already resulted in an increase in daily minimum streambed temperatures and decrease in daily amplitude and standard deviation. Scouring sites showed lower daily minimum streambed temperatures and higher daily amplitude and standard deviation compared with areas without sedimentation and scouring. As a limitation of the approach, groundwater discharge occurring at depositional and scouring areas cannot be identified by the metrics applied. Copyright © 2015 John Wiley & Sons, Ltd.     

A long‐term study of stable isotopes as tracers of processes governing water flow and quality in a lowland river basin: the upper Thames,UK

A long‐term study of O, H and C stable isotopes has been undertaken on river waters across the 7000‐km² upper Thames lowland river basin in the southern UK. During the period, flow conditions ranged from drought to flood. A 10‐year monthly record (2003–2012) of the main River Thames showed a maximum variation of 3‰ (δ¹⁸O) and 20‰ (δ²H), although interannual average values varied little around a mean of –6.5‰ (δ¹⁸O) and –44‰ (δ²H). A δ²H/δ¹⁸O slope of 5.3 suggested a degree of evaporative enrichment, consistent with derivation from local rainfall with a weighted mean of –7.2‰ (δ¹⁸O) and –48‰ (δ²H) for the period. A tendency towards isotopic depletion of the river with increasing flow rate was noted, but at very high flows (>100 m³/s), a reversion to the mean was interpreted as the displacement of bank storage by rising groundwater levels (corroborated by measurements of specific electrical conductivity). A shorter quarterly study (October 2011–April 2013) of isotope variations in 15 tributaries with varying geology revealed different responses to evaporation, with a well‐correlated inverse relationship between Δ¹⁸O and baseflow index for most of the rivers. A comparison with aquifer waters in the basin showed that even at low flow, rivers rarely consist solely of isotopically unmodified groundwater. Long‐term monitoring (2003–2007) of carbon stable isotopes in dissolved inorganic carbon (DIC) in the Thames revealed a complex interplay between respiration, photosynthesis and evasion, but with a mean interannual δ¹³C‐DIC value of –14.8 ± 0.5‰, exchange with atmospheric carbon could be ruled out. Quarterly monitoring of the tributaries (October 2011–April 2013) indicated that in addition to the aforementioned factors, river flow variations and catchment characteristics were likely to affect δ¹³C‐DIC. Comparison with basin groundwaters of different alkalinity and δ¹³C‐DIC values showed that the origin of river baseflow is usually obscured. The findings show that long‐term monitoring of environmental tracers can help to improve the understanding of how lowland river catchments function. Copyright © NERC 2015. Hydrological Processes © 2015 John Wiley & Sons, Ltd.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

Using hydrogeochemical signatures of stream water to assess pathways for rainfall events: towards a predictive model

Thanks to its simple division into agricultural and forestry land use, the Corbeira catchment (Galicia, Spain) is used as a case study to build a predictive model using hydrogeochemical signatures. Stream data acquired under recessional flow conditions over a one year period were obtained from a sampling station near the downstream end of the catchment, and using principal component analysis, it is shown that some of the analytical parameters are covariant, and some are negatively correlated. These findings support inferences about the pathways of rainfall in the catchment. Specific signatures may be associated with the dominant hydrological source, either surface runoff or subsurface waters: additionally, the dominant land use in that part of the catchment, where the flow originated, can also be predicted. The dominant runoff shows a strong covariance between suspended solids (SS) and particulate phosphorus (PP), with a clear negative correlation with pH. Dissolved organic carbon (DOC) data are associated with this covariant set when these compounds are available in the soils in question. Dissolved phosphorus, total organic nitrogen and dissolved nitrates are also associated with the same covariant set when the runoff flows through areas of extensive agricultural use. The SS ? PP covariance is less significant at lower flows. Typical base flow regimes show a significant covariance between salinity and pH, with a marked negative correlation with SS ? PP set, confirming the dominance of subsurface waters in the baseflow, as expected. Seasonally divergent DOC ? SS behaviour proves to be a useful tracer for rainfall regimes. The DOC trend shows a sinusoidal annual variation in amplitude, determined by the rainfall regime. As a result, flow from the catchment is dominated by surface water whenever there is synchronicity between the peaks of DOC and SS. Copyright © 2013 John Wiley & Sons, Ltd.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Submarine groundwater discharge revealed by aerial thermal infrared imagery: a case study on Jeju Island,Korea

Submarine groundwater discharge (SGD) is a global phenomenon that carries large volumes of groundwater and dissolved chemical species such as nutrient, metals, and organic compounds to coastal zones. We report the influence of SGD on the coastal waters of Jeju Island, Korea, using high‐resolution aerial thermal infrared (TIR) mapping techniques and field investigations. An aircraft‐based system was implemented using a cost‐effective TIR camera for aerial TIR mapping. Ground‐based calibrations and system integration with GPS/IMU (global positioning system/inertial measurement unit) were performed for the aerial systems. The aerial surveys showed distinct low‐temperature signatures of SGD along the coasts of Jeju Island, revealing large groundwater inputs from the coastal aquifers to the ocean. Multiple aerial surveys over a range of seasons and tidal stages revealed that SGD rates dynamically affect the sea surface temperature (SST) of the coastal zone. The in‐situ measurements supported that SGD has a substantial influence on the coastal water chemistry as well as SST. Our observations highlight the extent to which aerial‐based TIR mapping can serve as a powerful tool for studying SGD and other coastal processes. Copyright © 2016 John Wiley & Sons, Ltd.