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.     

A 53 year seasonally resolved oxygen and carbon isotope record from a modern Gibraltar speleothem: Reconstructed drip water and relationship to local precipitation

The response of a climate proxy against measured temperature, rainfall and atmospheric circulation patterns at sub-annual resolution is the ultimate test of proxy fidelity but very few data exist showing the level of correspondence between speleothem climate proxies and the instrumental climate record. Cave sites on the Gibraltar peninsula provide a unique opportunity to calibrate speleothem climate proxies with the longest known available precipitation isotopes and instrumental records. An actively growing speleothem sampled from New St. Michaels Cave in 2004 is composed of paired laminae consisting of light columnar calcite and a darker microsparitic calcite. Stable isotope analysis of samples micromilled in 100 μm steps at the equivalent of bi-monthly intervals reveals fabric-correlated annual cycles in carbon isotopes, oxygen isotopes and trace elements responding to seasonal changes in cave microclimate, hydrology and ventilation patterns. Calcite δ¹³C values reach a minimum in the light columnar fabric and evidence from trace element behaviour and cave monitoring indicates that this grows under cave ‘winter’ conditions of highest pCO₂, whereas the dark microsparitic calcite, characterised by elevated δ¹³C and δ¹⁸O values grows under low ‘summer’ pCO₂ conditions. Drip water δ¹³C_DIC reaches a minimum in March–April, at which time the attenuated δ¹⁸O signal becomes most representative of winter precipitation. An age model based on cycle counting and the position of the ¹⁴C bomb carbon spike yields a precisely dated winter oxygen isotope proxy of cave seepage water for comparison with the GNIP and instrumental climate record for Gibraltar. The δ¹⁸O characteristics of calcite deposited from drip water representing winter precipitation for each year can be derived from the seasonally resolved record and allows reconstruction of the δ¹⁸O drip water representing winter precipitation for each year from 1951–2004. These data show an encouraging level of correspondence (r² = 0.47) with the δ¹⁸O of rainfall falling each year between October and March and on a decadal scale the δ¹⁸O of reconstructed winter drip water mirrors secular change in mean winter temperatures.     

Epikarst hydrology and implications for stalagmite capture of climate changes at Grotta di Ernesto (NE Italy): results from long‐term monitoring

Grotta di Ernesto is a cave site well suited for palaeoclimate studies because it contains annually laminated stalagmites and was monitored from 1995 to the end of 2008 for microclimate, hydrology and hydrochemistry. Long‐term monitoring highlighted that cave drips show three different hydrological responses to rainfall and infiltration: (1) fast seasonal drips in the upper part of the cave, which are mostly fed by fractures, (2) slow seasonal drips, located at mid‐depth in the cave characterized by mixed feeding and (3) slow drips, mostly located in the deeper gallery, which are fed by seepage flow from bulk porosity with a minor fracture‐fed component. The slow drips display daily cycles during spring thaw. Monitoring also indicated that drip waters are only slightly modified by degassing within the soil zone and aquifer and by prior calcite precipitation. Hydrochemical studies show a clear seasonality in calcite saturation index, which results in most cave calcite precipitation occurring during late autumn and winter with similar amounts of precipitated calcite on most stalagmites, regardless of drip rate (discharge) differences. Drip rate, and drip rate variability, therefore, has a minor role in modulating the amount of annual calcite formation. In contrast, drip rate, when associated with moderate reduction in calcite saturation index, clearly influences stalagmite morphology. Increasing drip rate yields a passage from candle‐, to cone‐ to dome‐shaped stalagmites. Very high drip rates feed speleothems with flowstone morphology. In summary, monitoring provides information about the karst aquifer and how hydrology influences those physical and chemical characteristics of speleothems which are commonly used as climate proxies. Copyright © 2010 John Wiley & Sons, Ltd.     

Temporal and spatial variations in the discharge and dissolved organic carbon of drip waters in Beijing Shihua Cave,China

To detect the causal relationship between cave drip waters and stalagmite laminae, which have been used as a climate change proxy, three drip sites in Beijing Shihua Cave were monitored for discharge and dissolved organic carbon (DOC). Drip discharges and DOC were determined at 0 to 14‐day intervals over the period 2004–2006. Drip discharges show two types of response to surface precipitation variations: (1) a rapid response; and (2) a time‐lagged response. Intra‐annual variability in drip discharge is significantly higher than inter‐annual variability. The content of DOC in all drip waters varies inter‐ and intra‐annually and has good correlation with drip water discharge at the rapid response sites. High DOC was observed in July and August in the three years observed. The flushing of soil organic matter is dependent upon the intensity of rain events. The DOC content of drip water increases sharply above a threshold rainfall intensity (>50 mm d^?1) and shows several pulses corresponding with intense rain events (>25 mm d^?1). The DOC content was lower and less variable during the dry period than during the rainy period. The shape of DOC peak also varies from year to year as it is influenced by the intensity and frequency of rainfall. The different drip sites show marked differences in DOC response, which are dominated by hydrological behaviour linked to the recharge of the soil and karst micro‐fissure/porosity network. The results explain why not all stalagmite laminae are consistent with climate changes and suggest that the structure of the rainy season events could be preserved in speleothems. Copyright © 2008 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.     

Drying increases organic colloidal mobilization in the karst vadose zone: Evidence from a 15-year cave-monitoring study

Understanding the hydrological processes of colloids within the karst vadose zone is vital to the security of karst groundwater and providing appropriate paleohydrological explanations of colloid-facilitated metals in speleothem. This study addresses the mobilization mechanisms driving colloidal organic matter (COM) transport in the karst vadose zone using a 15-year long monthly monitoring dataset from a cave drip point (HS4) in Heshang Cave, Qingjiang Valley, China. Variations in COM concentrations were reported as the fluorescence difference values of raw and filtered (<0.22 μm) samples at an excitation wavelength of 320 nm and emission wavelength of ~400 nm. A fluorescence humification index (HIX) lower than 0.8 and an autochthonous index (BIX) higher than 1.2 indicated that the origin of COM was mainly from the karst vadose zone, rather than the soil zone. The COM concentration varied from 0.001 to 0.038 Raman Unit (RU), with evident seasonal fluctuations. Rising limbs for COM values occurred prior to rising limbs within a dripwater hydrograph; moreover, the COM peak values corresponding to the beginning of the increasing hydrograph generally suggested that the mobilization of COM reflected the movement of the air–water interface (AWI) in the karst vadose zone rather than rainfall intensity or flow velocity. COM peak values were positively correlated with the antecedent drying duration and negatively correlated with HIX values. These phenomena may be explained by the increased amount of organic matter that was aggregated and absorbed on the surface of carbonate in the karst vadose zone during a longer drying duration. Moreover, the longer drying duration was also beneficial to autochthonous biological activity, which subsequently decreased the HIX value of the organic matter in the karst vadose zone. The movement of AWI and the drying duration are both controlled by the outside weather conditions. This study is therefore conducive to evaluating the security of karst groundwater in response to climate change, and challenges prevailing paleoclimate interpretations of colloid-facilitated metal abundance timeseries reported from speleothems.     

Modelling hydrological processes influenced by soil,rock and vegetation in a small karst basin of southwest China

Hydrological processes in karst basins are controlled by permeable multimedia, consisting of soil pores, epikarst fractures, and underground conduits. Distributed modelling of hydrological dynamics in such heterogeneous hydrogeological conditions is a challenging task. Basing on the multilayer structure of the distributed hydrology‐soil‐vegetation model (DHSVM), a distributed hydrological model for a karst basin was developed by integrating mathematical routings of porous Darcy flow, fissure flow and underground channel flow. Specifically, infiltration and saturated flow movement within epikarst fractures are expressed by the ‘cubic law’ equation which is associated with fractural width, direction, and spacing. A small karst basin located in Guizhou province of southwest China was selected for this hydrological simulation. The model parameters were determined on the basis of field measurement and calibrated against the observed soil moisture contents, vegetation interception, surface runoff, and underground flow discharges from the basin outlet. The results show that due to high permeability of the epikarst zone, a significant amount of surface runoff is only generated after heavy rainfall events during the wet season. Rock exposure and the epikarst zone significantly increase flood discharge and decrease evapotranspiration (ET) loss; the peak flood discharge is directly proportional to the size of the aperture. Distribution of soil moisture content (SMC) primarily depends on topographic variations just after a heavy rainfall, while SMC and actual ET are dominated by land cover after a period of consecutive non‐rainfall days. The new model was able to capture the sharp increase and decrease of the underground streamflow hydrograph, and as such can be used to investigate hydrological effects in such rock features and land covers. Copyright © 2011 John Wiley & Sons, Ltd.     

广西岩溶洞穴沉积物与地震记录初探

广西是西南地区的一个岩溶大省,其岩溶面积占土地面积的42%.广西岩溶洞穴十分发育,至少有6万个洞穴.洞穴景观资源丰富,为广两的旅游创造了大量的财富.岩溶洞穴沉积物不仅具有旅游开发价值,而且蕴含了丰富的古气候环境信息,同时,也是地震信息的存储器,是研究地震遗迹的良好材料.地震活动不仅导致岩溶洞穴及洞穴沉积物发生破裂、坍塌...     

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.     

Advances in interpretation of subsurface processes with time‐lapse electrical imaging

Electrical geophysical methods, including electrical resistivity, time‐domain induced polarization, and complex resistivity, have become commonly used to image the near subsurface. Here, we outline their utility for time‐lapse imaging of hydrological, geochemical, and biogeochemical processes, focusing on new instrumentation, processing, and analysis techniques specific to monitoring. We review data collection procedures, parameters measured, and petrophysical relationships and then outline the state of the science with respect to inversion methodologies, including coupled inversion. We conclude by highlighting recent research focused on innovative applications of time‐lapse imaging in hydrology, biology, ecology, and geochemistry, among other areas of interest. Copyright © 2014 John Wiley & Sons, Ltd.     

A multi‐species dendroclimatic reconstruction of Chilko River streamflow,British Columbia,Canada

Dendroclimatological data were used to reconstruct the discharge history of Chilko River, which drains a glacierized watershed in the Coast Mountains of British Columbia. We correlated ring‐width records from Engelmann spruce (ES) (Picea engelmanni) and mountain hemlock (MH) (Tsuga mertensiana) trees to historical hydroclimate data. Over the period of record, spruce and hemlock radial growth correlates significantly with temperature and snow depth, respectively. We found that a multi‐species approach provided a better model fit and reconstructive power. Using these relationships, we developed generalized linear models for mean June, July, and June‐July discharge. The proxy records provide insights into streamflow variability of a typical Coast Mountains river over the past 240 years and confirm the long‐term influence of the Pacific Decadal Oscillation (PDO) on hydroclimatic regimes in the region. A relationship also exists between the reconstructed June‐July discharge record and the North Pacific (NP) Index, suggesting that winter atmospheric patterns over the North Pacific influence the hydrology of coastal British Columbia. Copyright © 2010 John Wiley & Sons, Ltd.     

Groundwater levels,climate and anthropogenic factors affect the hydrology and water quality of an intermittent and a regulated subtropical stream

Stream hydrology and water quality are highly interconnected and impacted by climate, land use and geology. We examined this connection using monitoring data from 2000 to 2019 for two streams with contrasting hydrological regimes—intermittent and regulated perennial—in subtropical Queensland, Australia. Our main objective was to evaluate relationships between groundwater levels, climate and flow regulation on the hydrology and water quality of an intermittent and a regulated subtropical stream. In intermittently flowing Lockyer Creek, flow was highly dependent on groundwater levels and occurred when the aquifer was recharged to elevations exceeding the upper 90-percentile value. With 9.4% of the catchment area in irrigated horticulture, flow in Lockyer Creek was also likely to be reduced by drawdown of the aquifer for irrigation, with no flow for 30% to 81% of days over the observation period for stations in Lockyer Creek. In contrast, flow in the mid-Brisbane River was continuous, regulated by discharge from a large upstream dam. Nutrient and suspended sediment concentrations in Lockyer Creek were generally higher than in the mid-Brisbane River, likely associated with runoff from agricultural areas adjacent to the stream, while the upstream dam likely reduced the concentration and variability of nutrients and suspended sediment in the mid-Brisbane River. During periods of low flow in the mid-Brisbane River, longitudinal changes in nutrient and suspended sediment concentrations occurred, notably a significant decrease in total and dissolved inorganic nitrogen concentrations downstream (p < 0.05), indicating a possible effect of in-stream algal uptake and denitrification. This study highlights the impact of human modifications on stream hydrology and water quality in the face of climate change. The findings can inform decision-making on groundwater irrigation or dam release control for water security.     

Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

Tree‐ring records unveil long‐term influence of the Pacific Decadal Oscillation on snowpack dynamics in the Stikine River basin,northern British Columbia

Long hydroclimate records are essential elements for the assessment and management of changing freshwater resources. These records are especially important in transboundary watersheds where international cooperation is required in the joint planning and management process of shared basins. Dendrochronological techniques were used to develop a multicentury record of April 1 snow water equivalent (SWE) for the Stikine River basin in northern British Columbia, Canada, from moisture‐sensitive white spruce (Picea glauca) tree rings. Explaining 43% of the instrumental SWE variability, to our knowledge, this research represents the first attempt to develop long‐term snowpack reconstructions in northern British Columbia. The results indicated that 15 extreme low April 1 SWE events occurred from 1789 to the beginning of the instrumental record in 1974. The reconstruction record also shows that the occurrence of hydrological extremes in the Stikine River basin is characterized by persistent below‐average periods in SWE consistent with phase shifts of the Pacific Decadal Oscillation (PDO). Spectral analyses indicate a very distinct in‐phase (positive) relationship between the multidecadal frequencies of variability (~40 years) extracted from the SWE tree‐ring reconstruction and other reconstructed winter and spring PDO indices. Comparison of the reconstructed SWE record with other tree‐ring‐derived PDO proxy records shows coherence at multidecadal frequencies of variability. The research has significant implications for regional watershed management by highlighting the hydrological response of the Stikine River basin to prior climate changes.     

Using 2D inversion of magnetic resonance soundings to locate a water-filled karst conduit

A new methodology for magnetic resonance sounding (MRS) data acquisition and interpretation was developed for locating water-filled karst cavities. This methodology was used to investigate the Ouysse karst system in the Poumeyssens shaft in the Causse de Gramat (France). A new 2D numerical MRS response model was designed for improved accuracy over the previous 1D MRS approach. A special survey performed by cave divers confirmed the accuracy of the MRS results. Field results demonstrated that in favourable conditions (a low EM noise environment and a relatively shallow, large target) the MRS method, used with a coincident transmitter/receiver loop, can be an effective tool for locating a water-filled karst conduit. It was shown numerically that because an a priori orientation of the MRS profile with the karst conduit is used in the inversion scheme (perpendicular for instance), any error in this assumption introduces an additional error in locating the karst. However, the resulting error is within acceptable limits when the deviation is less than 30°. The MRS results were compared with an electrical resistivity tomography (ERT) survey. It was found that in Poumeyssens, ERT is not able to locate the water-filled karst. On the other hand, ERT provides additional information about heterogeneities in the limestone.     

Time‐lapse electrical resistivity imaging of solute transport in a karst conduit

The use of electrical resistivity surveys to locate karst conduits has shown mixed success. However, time‐lapse electrical resistivity imaging combined with salt injection improves conduit detection and can yield valuable insight into solute transport behaviour. We present a proof of concept above a known karst conduit in the Kentucky Horse Park (Lexington, Kentucky). A salt tracer solution was injected into a karst window over a 45‐min interval, and repeat resistivity surveys were collected every 20 min along a 125‐m transect near a monitoring well approximately 750 m downgradient from the injection site. In situ fluid conductivity measurements in the well peaked at approximately 25% of the initial value about 3 h after salt injection. Time‐lapse electrical resistivity inversions show two broad zones at the approximate conduit depth where resistivity decreased and then recovered in general agreement with in situ measurements. Combined salt injection and electrical resistivity imaging are a promising tool for locating karst conduits. The method is also useful for gaining insight into conduit geometry and could be expanded to include multiple electrical resistivity transects. Copyright © 2015 John Wiley & Sons, Ltd.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulating effect of anthropogenic activities and climate variation on Liulin Springs discharge depletion by using the ARIMAX model

Based on the groundwater development process, and regional economic and social developing history, we divided the spring hydrological process of the Liulin Springs Basin into two periods: pre‐1973 and post‐1974. In the first period (i.e. 1957–1973), the spring discharge was affected by climate variation alone, and in the second period (i.e. 1974–2009), the spring discharge charge was influenced by both climate variation and human activities. A piecewise analysis strategy was used to differentiate the contribution of anthropogenic activities from climate variation on karst spring discharge depletion in the second period. Then, the ARIMAX model was applied to spring flow time series of the first period to develop a model for the effects of climate variation only. Using this model, we estimated the spring discharge in the second period solely under the influence of climate variation. Based on the water budget, we subtracted observed spring discharge from the estimated spring discharge and acquired the contribution of human activities on spring discharge depletion for the second period. The results of the analysis indicated that the contribution of climate variation to the spring discharge depletion is?0.20 m³/s from 1970s to 2000s. The contribution of anthropogenic activities to the spring flow depletion was ?2.56 m³/s in 2000s, which was about 13 times more than that of climate variation. Our analysis further indicates that groundwater exploitation only accounts for 29% of the spring flow depletion due to the effects of human activities. The remaining 71% of the depletion is likely to be caused by other human activities, including dam building, dewatering during coal mining, and deforestation. Copyright © 2012 John Wiley & Sons, Ltd.     

Using reanalysis‐driven regional climate model outputs for hydrology modelling

Bias correction methods are usually applied to climate model outputs before using these outputs for hydrological climate change impact studies. However, the use of a bias correction procedure is debatable, due to the lack of physical basis and the bias nonstationarity of climate model outputs between future and historical periods. The direct use of climate model outputs for impact studies has therefore been recommended in a few studies. This study investigates the possibility of using reanalysis‐driven regional climate model (RCM) outputs directly for hydrological modelling by comparing the performance of bias‐corrected and nonbias‐corrected climate simulations in hydrological simulations over 246 watersheds in the Province of Québec, Canada. When using RCM outputs directly, the hydrological model is specifically calibrated using RCM simulations. Two evaluation metrics (Nash–Sutcliffe efficiency [NSE] and transformed root mean square error [TRMSE]) and three hydrological indicators (mean, high, and low flows) are used as criteria for this comparison. Two reanalysis‐driven RCMs with resolutions of 45 km and 15 km are used to investigate the scale effect of climate model simulations and bias correction approaches on hydrology modelling. The results show that nonbias‐corrected simulations perform better than bias‐corrected simulations for the reproduction of the observed streamflows when using NSE and TRMSE as criteria. The nonbias‐corrected simulations are also better than or comparable with the bias‐corrected simulations in terms of reproducing the three hydrological indicators. These results imply that the raw RCM outputs driven by reanalysis can be used directly for hydrological modelling with a specific calibration of hydrological models using these datasets when gauged observations are scarce or unavailable. The nonbias‐corrected simulations (at a minimum) should be provided to end users, along with the bias‐corrected ones, especially for studying the uncertainty of hydrological climate change impacts. This is especially true when using an RCM with a high resolution, since the scale effect is observed when the RCM resolution increases from a 45‐km to a 15‐km scale.