Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Limitations of suction‐controlled triaxial tests in the characterization of unsaturated soils

Behavior of unsaturated soils is influenced by many factors, and the influences of these factors are usually coupled together. Suction‐controlled triaxial (SCTX) tests are considered to allow researchers to investigate influences of individual variables on unsaturated soils under specified stress path with controls of stresses, pore water, and air pressures. In the past 50 years, SCTX testing method has been established as a standard approach to characterize constitutive behavior of unsaturated soils. Most important concepts for modern unsaturated soil mechanics were developed upon results from the SCTX tests. Among these, one of the most important contributions in the constitutive modeling of elasto‐plastic behavior for unsaturated soils is the Barcelona basic model (BBM) proposed by Alonso et al. in 1990. The BBM successfully explained many features of unsaturated soils and received extensive acceptance. However, the SCTX tests are designed based upon the divide‐and‐conquer approach in which an implicit assumption is used: soil behavior is stress‐path independent. However, it is well‐established that unsaturated soil behavior is elasto‐plastic and stress‐path dependent. It is found that the SCTX tests in fact cannot control the stress path of an unsaturated soil during loading. This incapability, in combination with complicated loading/collapse behavior of unsaturated soils, makes the SCTX tests for characterizing unsaturated soil questionable. This paper discusses the limitations of the SCTX tests in the characterization of unsaturated soils. A possible solution to the problem was proposed based on a newly developed modified state surface approach. The discussions are limited for isotropic conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Characterisation of hydroclimatological trends and variability in the Lake Naivasha basin,Kenya

Recent hydro‐climatological trends and variability characteristics were investigated for the Lake Naivasha basin with the aim of understanding the changes in water balance components and their evolution over the past 50 years. Using a Bayesian change point analysis and modified Mann–Kendall tests, time series of annual mean, maximum, minimum, and seasonal precipitation and flow, as well as annual mean lake volumes, were analysed for the period 1960–2010 to uncover possible abrupt shifts and gradual trends. Double cumulative curve analysis was used to investigate the changes in hydrological response attributable to either human influence or climatic variability. The results indicate a significant decline in lake volumes at a mean rate of 9.35 × 10⁶ m³ year^?1. Most of the river gauging stations showed no evidence of trends in the annual mean and maximum flows as well as seasonal flows. Annual minimum flows, however, showed abrupt shifts and significant (upward/downward) trends at the main outlet stations. Precipitation in the basin showed no evidence of abrupt shifts, but a few stations showed gradual decline. The observed changes in precipitation could not explain the decline in both minimum flows and lake volumes. The findings show no evidence of any impact of climate change for the Lake Naivasha basin over the past 50 years. This implies that other factors, such as changes in land cover and infrastructure development, have been responsible for the observed changes in streamflow and lake volumes. Copyright © 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.     

Discharge‐sediment processes of the Zhadang glacier on the Tibetan Plateau measured with a high frequency data acquisition system

In high elevation cold regions of the Tibetan Plateau, suspended sediment transfer from glacier meltwater erosion is one of the important hydrological components. The Zhadang glacier is a typical valley‐type glacier in the Nyainqentanglha Mountains on the Tibetan Plateau. To make frequent and long period records of meltwater runoff and sediment processes in the very high elevation and isolated regions, an automatic system was installed near the glacier snout (5400 m a.s.l) in August 2013, to measure the transient discharge and sediment processes at 5‐min interval, which is shorter than the time span for the water flow to traverse the catchment from the farthest end to the watershed outlet. Diurnal variations of discharge, and suspended sediment concentration (SSC) were recorded at high frequency for the Zhadang glacier, before suspended sediment load (SSL) was computed. Hourly SSC varied from the range of 0.2 kg/m³ to 0.5 kg/m³ (at 8:00–9:00) to the range of 2.0 kg/m³ to 4.0 kg/m³ (at 17:00–18:00). The daily SSL was 32.24 t during the intense ablation period. Hourly SSC was linearly correlated with discharge (r = 0.885**, n = 18, p < 0.01). A digit‐eight hysteresis loop was observed for the sediment transport in the glacier area. Air temperature fluctuations influence discharge, and then result in the sediment variations. The results of this study provide insight into the responses of suspended sediment delivery processes with a high frequency data in the high elevation cold regions. Copyright © 2016 John Wiley & Sons, Ltd.     

A probabilistic approach for computing water retention of particulate systems from statistics of grain size and tessellated pore network

The paper offers an analytical determination of the hydraulic properties of an unsaturated soil with reference to its retention curve, which describes the relationship between the volumetric water content and capillarity through matric suction. The analysis combines a particulate approach focused on the physics at the pore scale, including microstructural aspects, with a probabilistic approach where the void space and grain size are considered as random variables. In the end, the soil water characteristic curve of an unsaturated granular medium along a drying path can be derived analytically based on the sole information of particle size distribution. The analysis hinges on the tessellation of a wet granular system into an assemblage of tetrahedral unit cells revealing a pore network upon which capillary physics are computed with respect to pore throat invasion by a non-wetting fluid with evolving pendular capillary bridges. The crux of the paper is to pass from particle size probability distribution to a matching void space distribution to eventually reveal key information such as void cell and solid volume statistics. Making reasonable statistically based assumptions to render calculations tractable, the water retention curve can be readily constructed. Model predictions compare quite favourably with experimental data available for actual soils, especially in the high saturation range. Having a sound scientific basis, the model can be made amenable to address a variety of soils with a wider range of particle sizes.     

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.