A multi‐storage groundwater concept for the SWAT model to emphasize nonlinear groundwater dynamics in lowland catchments

Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.     

Flow dynamics at the continental scale: Streamflow correlation and hydrological similarity

Streamflow variability in space and time critically affects anthropic water uses and ecosystem services. Unfortunately, spatiotemporal patterns of flow regimes are often unknown, as discharge measurements are usually recorded at a limited number of hydrometric stations unevenly distributed along river networks. Advances in understanding the physical processes that control the spatial patterns of river flows are therefore necessary to predict water availability at ungauged locations or to extrapolate pointwise streamflow observations. This work explores the use of the spatial correlation of river flows as a metric to quantify the similarity between hydrological responses of two catchments. Following a stochastic framework, 340,000 cross‐correlations between pairs of daily streamflows time series are predicted at a seasonal timescale across the contiguous United States using 413 catchments of the MOPEX dataset. Model predictions of streamflow correlation obtained in absence of run‐off information are successfully used to identify catchment outlets sharing similar discharge dynamics and flow regimes across a broad range of geomorphoclimatic conditions, without relying on calibration. The selection of reference streamgauges based on predicted streamflow correlation generally outperforms the selection based on spatial proximity, especially as the density of available gauged sections decreases. Interestingly, correlated outlets share a broad spectrum of hydrological signatures (mean discharge, flow variability, and recession properties), suggesting that catchments forced by analogous frequency and intensity of effective rainfall events might exhibit common geomorphoecological traits leading to similar hydrological responses. The proposed framework provides a physical basis to assist the regionalization of flow dynamics and to interpret the spatial variability of flow regimes along stream networks.     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

Estimation of discharges of water flows and debris floods in a small watershed

Floods in small mountainous watersheds cover a wide spectrum of flow. They can range from clear water flows and hyperconcentrated flows to debris floods and debris flows, and calculation of the peak discharge is crucial for predicting and mitigating such hazards. To determine the optimal approach for discharge estimation, this study compared water flow monitoring hydrographs to investigate the performance of five hydrological models that incorporate different runoff yields and influx calculation methods. Two of the models performed well in simulating the peak discharge, peak time, and total flow volume of the water flood. The ratio (γ) of the monitored debris flood discharge (Q_d) to the simulated water flow discharge (Q_w) was investigated. Qualitatively, γ initially increased with Q_w but then decreased when Q_w exceeded a certain threshold, which corresponded to rainfall of 95 and 120 mm in a 6- and 24-h event with a normal distribution of precipitation, respectively. The decrease might be attributable to a threshold of sediment availability being reached, beyond which increased flow rate is not matched by increased sediment input in the large watershed. Uncertainty of hydrological calculation was evaluated by dividing the catchment into sub-basins and adopting different rainfall time steps as input. The efficiency of using a distributed simulation exhibited marginal improvement potential compared with a lumped simulation. Conversely, the rainfall time step input significantly affected the simulation results by delaying the peak time and decreasing the peak discharge. This research demonstrates the applicability of a discharge estimation method that combines a hydrological water flow simulation and an estimation of γ. The results were verified on the basis of monitored flow densities and videos obtained in two watersheds with areas of 2.34 and 32.4 km².     

Interpretation of spring recession curves

Recession curves contain information on storage properties and different types of media such as porous, fractured, cracked lithologies and karst. Recession curve analysis provides a function that quantitatively describes the temporal discharge decay and expresses the drained volume between specific time limits (Hall 1968). This analysis also allows estimating the hydrological significance of the discharge function parameters and the hydrological properties of the aquifer. In this study, we analyze data from perennial springs in the Judean Mountains and from others in the Galilee Mountains, northern Israel. All the springs drain perched carbonate aquifers. Eight of the studied springs discharge from a karst dolomite sequence, whereas one flows out from a fractured, slumped block of chalk. We show that all the recession curves can be well fitted by a function that consists of two exponential terms with exponential coefficients alpha1 and alpha2. These coefficients are approximately constant for each spring, reflecting the hydraulic conductivity of different media through which the ground water flows to the spring. The highest coefficient represents the fast flow, probably through cracks, or quickflow, whereas the lower one reflects the slow flow through the porous medium, or baseflow. The comparison of recession curves from different springs and different years leads to the conclusion that the main factors that affect the recession curve exponential coefficients are the aquifer lithology and the geometry of the water conduits therein. In normal years of rainy winter and dry summer, alpha1 is constant in time. However, when the dry period is longer than usual because of a dry winter, alpha1 slightly decreases with time.     

Adjustment of self-formed bankfull channel geometry of meandering rivers: modelling study

We consider the evolution of the hydraulic geometry of sand-bed meandering rivers. We study the difference between the timescale of longitudinal river profile adjustment and that of channel width and depth adjustment. We also study the effect of hydrological regime alteration on the evolution of bankfull channel geometry. To achieve this, a previously developed model for the spatiotemporal co-evolution of bankfull channel characteristics, including bankfull discharge, bankfull width, bankfull depth and down-channel bed slope, is used. In our modelling framework, flow variability is considered in terms of a specified flow duration curve. Taking advantage of this unique feature, we identify the flow range responsible for long-term bankfull channel change within the specified flow duration curve. That is, the relative importance of extremely high short-duration flows compared to moderately high longer duration flows is examined. The Minnesota River, MN, USA, an actively meandering sand-bed stream, is selected for a case study. The longitudinal profile of the study reach has been in adjustment toward equilibrium since the end of the last glaciation, while its bankfull cross-section is rapidly widening due to hydrological regime change in the last several decades. We use the model to demonstrate that the timescale for longitudinal channel profile adjustment is much greater than the timescale for cross-sectional profile adjustment due to a lateral channel shift. We also show that hydrological regime shift is responsible for the recent rapid widening of the Minnesota River. Our analysis suggests that increases in the 5–25% exceedance flows play a more significant role in recent bankfull channel enlargement of the Minnesota River than increase in either the 0.1% exceedance flow or the 90% exceedance flow. © 2020 John Wiley & Sons, Ltd.     

Effect of climate change on overland flow generation: a case study in central Germany

The impact of global climate change on runoff components, especially on the type of overland flow, is of utmost significance. High‐resolution temporal rainfall plays an important role in determining the hydrological response of quick runoff components. However, hydrological climate change scenario analyses with high temporal resolution are rare. This study investigates the impact of climate change on discharge peak events generated by rainfall, snowmelt, and soil‐frost induced runoff using high‐resolution hydrological modelling. The study area is Schäfertal catchment (1.44 km²) in the lower Harz Mountains in central Germany. The WaSiM‐ETH hydrological model is used to investigate the rainfall response of runoff components under near future (2021–2050) and far‐distant future (2071–2100) climatic conditions. Disaggregated daily climate variables of WETTREG2010 SRES scenario A1B are used on a temporal resolution of 10 min. Hydrological model parameter optimization and uncertainty analysis was conducted using the Differential Evolution Adaptive Metropolis (DREAM_(ZS)) uncertainty tool. The scenario results show that total runoff and interflow will increase by 3.8% and 3.5% in the near future and decrease by 32.85% and 31% in the far‐distant future compared to the baseline scenario. In contrast, overland flow and the number and size of peak runoff will decrease moderately for the near future and drastically for the far‐distant future compared to the baseline scenario. We found the strongest decrease for soil‐frost induced discharge peaks at 79.6% in the near future and at 98.2% in the far‐distant future scenario. It can be concluded that high‐resolution hydrological modelling can provide detailed predictions of future hydrological regimes and discharge peak events of the catchment. Copyright © 2014 John Wiley & Sons, Ltd.     

Climate impacts on river flow: projections for the Medway catchment,UK, with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.     

A framework to assess the cumulative impacts of dams on hydrological regime: A case study of the Yangtze River

Comparisons of flow time series between preimpact and postimpact periods have been widely used to determine hydrological alterations caused by reservoir operation. However, preimpact and postimpact periods might also be characterized by different climatological properties, a problem that has not been well addressed. In this study, we propose a framework to assess the cumulative impact of dams on hydrological regime over time. The impacts of the Three Gorges Dam (TGD) on the flow regime of the Yangtze River were investigated using this framework. We reconstructed the unregulated flow series to compare with the regulated flow series during the same period (2010 to 2015). Eco‐surplus and eco‐deficit and the Indicators of Hydrologic Alteration (IHA) parameters were used to examine hydrological regime change. Among 32 IHA parameters, Wilcoxon signed‐rank test and principal component analysis identified the October median flow, 1‐ and 3‐day maximum flows, 1‐day minimum flow, and rise rate as representative indicators of hydrological alterations. Eco‐surplus and eco‐deficit showed that the reservoir also changed the seasonal regime of the flows by reducing autumn flow and increasing winter flow. Changes in annual extreme flows and October flows lead to negative ecological implications downstream of the TGD. Ecological considerations should be taken into account during operation of the TGD in order to mitigate the negative effects on the fluvial ecosystem in the middle reach of Yangtze River. The framework proposed here could be a robust method to assess the cumulative impacts of reservoir operation over time.     

On characterizing the temporal dominance patterns of model parameters and processes

Diagnostic analyses of hydrological models intend to improve the understanding of how processes and their dynamics are represented in models. Temporal patterns of parameter dominance could be precisely characterized with a temporally resolved parameter sensitivity analysis. In this way, the discharge conditions are characterized, that lead to a parameter dominance in the model. To achieve this, the analysis of temporal dynamics in parameter sensitivity is enhanced by including additional information in a three‐tiered framework on different aggregation levels. Firstly, temporal dynamics of parameter sensitivity provide daily time series of their sensitivities to detect variations in the dominance of model parameters. Secondly, the daily sensitivities are related to the flow duration curve (FDC) to emphasize high sensitivities of model parameters in relation to specific discharge magnitudes. Thirdly, parameter sensitivities are monthly averaged separately for five segments of the FDC to detect typical patterns of parameter dominances for different discharge magnitudes. The three methodical steps are applied on two contrasting catchments (upland and lowland catchment) to demonstrate how the temporal patterns of parameter dynamics represent different hydrological regimes. The discharge dynamic in the lowland catchment is controlled by groundwater parameters for all discharge magnitudes. In contrast, different processes are relevant in the upland catchment, because the dominances of parameters from fast and slow runoff components in the upland catchment are changing over the year for the different discharge magnitudes. The joined interpretation of these three diagnostic steps provides deeper insights of how model parameters represent hydrological dynamics in models for different discharge magnitudes. Thus, this diagnostic framework leads to a better characterization of model parameters and their temporal dynamics and helps to understand the process behaviour in hydrological models. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of spring flows from a karst aquifer with an artificial neural network

In China, 9·5% of the landmass is karst terrain and of that 47,000 km² is located in semiarid regions. In these regions the karst aquifers feed many large karst springs within basins of thousands of square kilometres. Spring discharges reflect the fluctuation of ground water level and variability of ground water storage in the basins. However, karst aquifers are highly heterogeneous and monitoring data are sparse in these regions. Therefore, for sustainable utilization and conservation of karst ground water it is necessary to simulate the spring flows to acquire better understanding of karst hydrological processes. The purpose of this study is to develop a parsimonious model that accurately simulates spring discharges using an artificial neural network (ANN) model. The karst spring aquifer was treated as a non‐linear input/output system to simulate the response of karst spring flow to precipitation and applied the model to the Niangziguan Springs, located in the east of Shanxi Province, China and a representative of karst springs in a semiarid area. Moreover, the ANN model was compared with a previous time‐lag linear model and it was found that the ANN model performed better. Copyright © 2007 John Wiley & Sons, Ltd.     

A revised range of variability approach considering the periodicity of hydrological indicators

Comprehensively assessing the hydrological alteration of river flows is a prerequisite for the scientific management, protection, and restoration of rivers. The range of variability approach (RVA), which is based on indicators of hydrological alteration (IHAs), is a widely used method to evaluate hydrological alteration. However, the RVA only considers the frequency of each IHA, neglecting the equally important temporal order of these IHAs. The order of IHA event can be reflected by its periodicity. On the basis of the RVA, in this study, we propose a revised RVA that considers both the frequency and periodicity of IHAs. In the revised RVA, first, the periodic time of each IHA is identified; next, the periodicity alteration (P) of river flow is calculated by comparing the periodic times of the pre‐impact‐period and post‐impact‐period IHAs; finally, P and the frequency alteration(D) in traditional RVA are incorporated into a single index (H) to reflect the overall hydrological alteration. A case study of the Xi Dayang (XDY) Reservoir and rearranged flow suggests that the traditional RVA underestimates hydrological alterations because it neglects the alteration of periodicity. Compared with the traditional RVA and its alternatives, the revised RVA could give a more comprehensive representation of hydrological alteration caused by human and nature impacts. Thus, better protection of an ecosystem could be obtained by applying this method in the evaluation and management of water resources. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydrologic and geomorphic flow thresholds in the Lower Brazos River,Texas, USA

Abstract

River science and management often require a design or reference discharge. The common (and sometimes unavoidable) use of such discharges may, however, obscure the fact that the magnitude and frequency of critical flows can differ due to various hydrological, geomorphological, and ecological criteria. Threshold stages and discharges were identified for six lower Brazos River, Texas gaging stations corresponding to thalweg connectivity, bed inundation, high sub-banktop flows, channel–floodplain connectivity (CFC), and overbank flooding. Critical flows were also identified for estimated thresholds for sandy bedform and medium gravel mobility, critical specific stream power for potential channel modifications, and cohesive-bank channel erosion. These thresholds have variable relationships to mean, median, and maximum flows. For four of the six stations, daily recurrence probabilities for all but flood flows are at least 1%, and as high as 11%. All stations achieve channel–floodplain connectivity at stages less than banktop. Estimated threshold flows for sediment mobility and channel erosion occur relatively frequently, with daily probabilities of 2–77%. Critical flows for bank erosion occur least often, and for sandy bedform and gravel mobility most often. Thalweg connectivity is always maintained at all sites, while bed inundation flows have a daily probability of about 80% or more. Overall, results suggest that no single flow level is dominant in hydrological or geomorphic dynamics, and that the frequency of a given threshold varies considerably even along a single river. The results support the idea that multiple flow levels and ranges are necessary to create and maintain the hydrological, geomorphological, and ecological characteristics of rivers, and that no single flow level is a reliable determinant of fluvial state.

Editor Z.W. Kundzewicz; Associate editor Q. Zhang     

Modelling the effects of permafrost loss on discharge from a wetland‐dominated,discontinuous permafrost basin

Permafrost degradation in the peat‐rich southern fringe of the discontinuous permafrost zone is catalysing substantial changes to land cover with expansion of permafrost‐free wetlands (bogs and fens) and shrinkage of forest‐dominated permafrost peat plateaux. Predicting discharge from headwater basins in this region depends upon understanding and numerically representing the interactions between storage and discharge within and between the major land cover types and how these interactions are changing. To better understand the implications of advanced permafrost thaw‐induced land cover change on wetland discharge, with all landscape features capable of contributing to drainage networks, the hydrological behaviour of a channel fen sub‐basin in the headwaters of Scotty Creek, Northwest Territories, Canada, dominated by peat plateau–bog complexes, was modelled using the Cold Regions Hydrological Modelling platform for the period of 2009 to 2015. The model construction was based on field water balance observations, and performance was deemed adequate when evaluated against measured water balance components. A sensitivity analysis was conducted to assess the impact of progressive permafrost loss on discharge from the sub‐basin, in which all units of the sub‐basin have the potential to contribute to the drainage network, by incrementally reducing the ratio of wetland to plateau in the modelled sub‐basin. Simulated reductions in permafrost extent decreased total annual discharge from the channel fen by 2.5% for every 10% decrease in permafrost area due to increased surface storage capacity, reduced run‐off efficiency, and increased landscape evapotranspiration. Runoff ratios for the fen hydrological response unit dropped from 0.54 to 0.48 after the simulated 50% permafrost area loss with a substantial reduction of 0.47 to 0.31 during the snowmelt season. The reduction in peat plateau area resulted in decreased seasonal variability in discharge due to changes in the flow path routing, with amplified low flows associated with small increases in subsurface discharge, and decreased peak discharge with large reductions in surface run‐off.     

Hydrology of cave drip waters at varying bedrock depths from a karst system in southeastern Australia

The study of how cave drip‐water discharge responds to recharge events is fundamental to evaluating the potential of actively forming speleothems as high‐resolution climate archives. Most previous research has focused on caves of the Northern Hemisphere middle latitudes, where recharge is strongly seasonal. Few studies have explored drip‐water behaviour from regions where the expected seasonal rainfall pattern is significantly perturbed on an irregular basis by changing regional atmospheric circulation patterns. Here, we report the results of a 4‐year study of cave drip‐water–climate relationships from two caves in eastern Australia. The discharge of 10 drip sites located beneath bedrock thicknesses of 12, 22 and 45 m was monitored either continuously (using automated infrared sensors) or at discrete approximately monthly intervals and compared with local rainfall and water balance data. The study period traversed two major droughts, including the severe 2002–2003 El Niño. Drips at 12 and 22 m depths responded almost simultaneously to individual recharge events, although the time lag between individual events varied according to the volume of recharge and pre‐event storage. Overall, a steady decline in discharge is evident through the moisture‐deficit period, with increased flows through phases of positive water balance. Speleothems growing at these and similar shallow‐chamber sites have potential for reconstructing palaeo‐rainfall trends at high‐resolution, although the highly variable nature of year‐to‐year recharge would make it difficult to obtain data on a calendrical time‐scale. Drips at 45 m depth did not respond consistently to individual recharge events and displayed hydrological behaviour markedly dissimilar to one another and to the near‐surface drip sites, indicating great complexity in karst architecture and the absence of fissure flow. Although speleothems at this depth may well preserve information on longer‐term rainfall trends, their potential to encode a palaeo‐rainfall variability signal at interannual resolution is poor. Copyright © 2007 John Wiley & Sons, Ltd.     

Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone

Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km² catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimating minimum streamflow from measurements at ungauged sites in regions with streamflow‐gauging networks

Estimation of low flows in rivers continues to be a vexing problem despite advances in statistical and process‐based hydrological models. We develop a method to estimate minimum streamflow at seasonal to annual timescales from measured streamflow based on regional similarity in the deviations of daily streamflow from minimum streamflow for a period of interest. The method is applied to 1,019 gauged sites in the Western United States for June to December 2015. The gauges were clustered into six regions with distinct timing and magnitude of low flows. A gamma distribution was fit each day to the deviations in specific discharge (daily streamflow divided by drainage area) from minimum specific discharge for gauges in each region. The Kolmogorov–Smirnov test identified days when the gamma distribution was adequate to represent the distribution of deviations in a region. The performance of the gamma distribution was evaluated at gauges by comparing daily estimates of minimum streamflow with estimates from area‐based regression relations for minimum streamflow. Each region had at least 8 days during the period when streamflow measurements would provide better estimates than the regional regression equation, but the number of such days varied by region depending on aridity and homogeneity of streamflow within the region. Synoptic streamflow measurements at ungauged sites have value for estimating minimum streamflow and improving the spatial resolution of hydrological model in regions with streamflow‐gauging networks.     

Influence of bedrock groundwater on streamflow characteristics in a volcanic catchment

Groundwater movements in volcanic mountains and their effects on streamflow discharge and representative elementary area (REA) have remained largely unclear. We surveyed the discharge and chemical composition of spring and stream water in two catchments: the Hontani river (NR) catchment (6.6 km²) and the Hosotani river (SR) catchment (4.0 km²) at the southern part of Daisen volcano, Japan. Daisen volcano is a young volcano (17 × 10³ years) at an early stage of erosion. Our study indicated that deep groundwater that moved through thick lava and pyroclastic flows and that could not be explained by shallow movements controlled by surface topography contributed dominantly to streamflow at larger catchment areas. At the NR catchment, the deep groundwater contribution clearly increased at a catchment boundary defined by an area of 3.0 km² and an elevation of 800 m. At the SR catchment, the contribution deep groundwater to the stream also increased suddenly at a boundary threshold of 2.0 and 700 m. Beyond these thresholds, the contributions of deep bedrock groundwater remained constant, indicating that the REA is between 2 and 3 km² at the observed area. These results indicate that the hydrological conditions of base flow were controlled mainly by the deep bedrock groundwater that moved through thick lava and pyroclastic flows in the undissected volcanic body of the upper part of the catchment. Our study demonstrates that deep and long groundwater movements via a deep bedrock layer including thick deposits of volcanic materials at the two catchments on Daisen volcano strongly determined streamflow discharge instead of the mixing of small‐scale hydrological conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.