The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution          下载免费PDF全文

A. A. Ameli  J. J. McDonnell  K. Bishop 《水文研究》2016,30(14):2438-2450

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Analysis of rectangular piles subjected to lateral loads in nonhomogeneous soil using a Winkler model approach          下载免费PDF全文

Hiroyoshi Hirai 《国际地质力学数值与分析法杂志》2015,39(9):937-968

An analytical approach using a Winkler model based on two lateral soil displacement components in a three‐dimensional soil is investigated to provide analytical solutions of horizontal response of a rectangular pile subjected to lateral loads in nonhomogeneous soil. The two lateral displacement components of a soil surrounding the rectangular pile are represented by the Fourier series of displacement potential functions in the elastic three‐dimensional analysis. The lateral stiffness coefficient of the rectangular pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient taking into account rocking rotation of a rigid pile shaft. The relationship between horizontal displacement, rotation, moment, and shear force for the rectangular pile subjected to horizontal loads in nonhomogeneous soil is obtainable in the form of the recurrence equation. The formulation of lateral displacement and rotation for a rectangular pile subjected to lateral loads on the pile base in nonhomogeneous soil is proposed by taking into account Mindlin's equation and the equivalent thickness for soil layers in the equivalent elastic method. The difference of lateral behavior between square and circular piles subjected to lateral loads is insignificant. The effect of aspect ratio of the rectangular pile on the lateral behavior is great for the lower stiffness ratio between pile and soil and the larger length–equivalent diameter ratio. The effect of the value of Poisson's ratio of soil on lateral stiffness coefficient is relatively small except Poisson's ratio close to 0.5. The comparison of the results calculated by the current method for a rectangular pile subjected to lateral loads in nonhomogeneous soil has shown good agreement with those obtained from the analytical methods and the finite element method. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment          下载免费PDF全文

Christian Birkel  Chris Soulsby 《水文研究》2016,30(24):4477-4493

Assessing catchment runoff response remains a key research frontier because of limitations in current observational techniques to fully characterize water source areas and transit times in diverse geographical environments. Here, we report a study that combines empirical data with modelling to identify dominant runoff processes in a sparsely monitored humid tropical catchment. The analysis integrated isotope tracers into conceptual rainfall–runoff models of varying complexity (from 5 to 11 calibrated parameters) that are able to simulate discharge and tracer concentrations and track the evolving age of stream water exiting the catchment. The model structures can be seen as competing hypotheses of catchment functioning and were simultaneously calibrated against uncertain streamflow gaugings and a 2‐year daily isotope rainfall–runoff record. Comparison of the models was facilitated using global parameter sensitivity analysis and the resulting effect on calibration. We show that a variety of tested model structures reproduced water and tracer dynamics in stream, but the simpler models failed to adequately reproduce both. The resulting water age distributions of the tested models varied significantly with little similarity between the stream water age and stored water age distributions. The sensitivity analysis revealed that only some of the more complex models (from eight parameters) could be better constrained to infer more plausible water age distributions and catchment storage estimates. These models indicated that the age of water stored in the catchment is generally older compared with the age of water fluxes, with evapotranspiration age being younger compared with streamflow. However, the water age distributions followed a similar temporal behaviour dominated by climatic seasonality. Stream water ages increased during the dry season (greater than 1 year) and decreased with increased streamflow (a few weeks old) during the wet season. We further show that the ratios of the streamwater age to stored water age distribution and the water age distribution of actual evapotranspiration to the stored water age distribution from constrained models could potentially serve as useful hydrological indicators of catchment functioning. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

A two‐dimensional coupled hydromechanical discontinuum model for simulating rock hydraulic fracturing          下载免费PDF全文

Yu‐Yong Jiao  Huan‐Qiang Zhang  Xiu‐Li Zhang  Hai‐Bo Li  Qing‐Hui Jiang 《国际地质力学数值与分析法杂志》2015,39(5):457-481

Within the framework of our discontinuous deformation analysis for rock failure algorithm, this paper presents a two‐dimensional coupled hydromechanical discontinuum model for simulating the rock hydraulic fracturing process. In the proposed approach, based on the generated joint network, the calculation of fluid mechanics is performed first to obtain the seepage pressure near the tips of existing cracks, and then the fluid pressure is treated as linearly distributed loads on corresponding block boundaries. The contribution of the hydraulic pressure to the initiation/propagation of the cracks is considered by adding the components of these blocks into the force matrix of the global equilibrium equation. Finally, failure criteria are applied at the crack tips to determine the occurrence of cracking events. Several verification examples are simulated, and the results show that this newly proposed numerical model can simulate the hydraulic fracturing process correctly and effectively. Although the numerical and experimental verifications focus on one unique preexisting crack, because of the capability of discontinuous deformation analysis in simulating block‐like structures, the proposed approach is capable of modeling rock hydraulic fracturing processes. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Hydrologic and water isotope characterization of a regulated Canadian Shield river basin     

Nancy England  April L. James  Krystopher J. Chutko  Richard S. Pyrce  Huaxia Yao 《水文研究》2019,33(6):905-919

Stable water isotope surveys have increasingly been integrated into river basins studies, but fewer have used them to evaluate impact of hydropower regulation. This study applies hydrologic and water isotope survey approaches to a Canadian Shield river basin with both regulated and natural flows. Historical streamflow records were used to evaluate the influence of three hydroelectric reservoirs and unregulated portions of the basin on downstream flows and changes in water level management implemented after an extreme flood year (1979). In 2013, water isotope surveys of surface and source waters (e.g., rainfall, groundwater, snowmelt) were conducted to examine spatial and temporal variation in contributions to river flow. Seasonal changes in relative groundwater contribution were assessed using a water‐isotope mass balance approach. Within the basin, two regulated reservoirs exhibited inverted hydrographs with augmented winter flows, whereas a third exhibited a hydrograph dominated by spring snowmelt. In 2013, spatial variation in rain‐on‐snow and air temperatures resulted in a critical lag in snowmelt initiation in the southern and northern portions of the basin resulting in a dispersed, double peak spring hydrograph, contrasting with 1979 when a combination of rain‐on‐snow and coincident snowmelt led to the highest flood on record. Although eastern basin reservoirs become seasonally enriched in δ¹⁸O and δ²H values, unregulated western basin flows remain less variable due to groundwater driven baseflow with increasing influence downstream. Combined analysis of historical streamflow (e.g., flood of 1979, drought of 2010) and the 2013 water isotope surveys illustrate extreme meteorological conditions that current management activities are unable to prevent. In this study, the influence of evaporative fractionation on large surface water reservoirs provides important evidence of streamflow partitioning, illustrating the value of stable water isotope tracers for study of larger catchments.  相似文献   

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils     

Panagiotis Sitarenios  Michael Kavvadas 《国际地质力学数值与分析法杂志》2020,44(4):435-454

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.  相似文献   

How preferential flow delivers pre‐event groundwater rapidly to streams     

Stephen R.H. Worthington 《水文研究》2019,33(17):2373-2380

Groundwater often accounts for a substantial fraction of flood hydrographs, but the processes responsible for this have been unclear. However, many aquifers have preferential flow and this explains how aquifers can be so responsive. In bedrock aquifers, weathering enhances the connectivity and apertures along the most efficient flow paths and hence enhances the permeability. This results in celerities and velocities of the preferential flow in these dual‐porosity aquifers that are two to three orders of magnitude higher than if the aquifers behaved as single‐porosity media. The celerities have been determined from artificial and natural flood pulses, from tidal lags, and from pumping tests. Preferential‐flow velocities have been calculated from tests using applied tracers. Celerities in bedrock aquifers are typically one to two orders of magnitude faster than velocities. The ubiquitous preferential flow in aquifers provides an additional explanation, besides groundwater ridging, for the rapid release of groundwater to streams during storm events.  相似文献   

Modelling atmospheric and hydrologic processes for assessment of meadow restoration impact on flow and sediment in a sparsely gauged California watershed          下载免费PDF全文

N. Ohara  M. L. Kavvas  Z. Q. Chen  L. Liang  M. Anderson  J. Wilcox  L. Mink 《水文研究》2014,28(7):3053-3066

The restoration of meadowland using the pond and plug technique of gully elimination was performed in a 9‐mile segment along Last Chance Creek, Feather River Basin, California, in order to rehabilitate floodplain functions such as mitigating floods, retaining groundwater, and reducing sediment yield associated with bank erosion and to significantly alter the hydrologic regime. However, because the atmospheric and hydrological conditions have evolved over the restoration period, it was difficult to obtain a comprehensible evaluation of the impact of restoration activities by means of field measurements. In this paper, a new use of physically based models for environmental assessment is described. The atmospheric conditions over the sparsely gauged Last Chance Creek watershed (which does not have any precipitation or weather stations) during the combined historical critical dry and wet period (1982–1993) were reconstructed over the whole watershed using the atmospheric fifth‐generation mesoscale model driven with the US National Center for Atmospheric Research and US National Center for Environmental Prediction reanalysis data. Using the downscaled atmospheric data as its input, the watershed environmental hydrology (WEHY) model was applied to this watershed. All physical parameters of the WEHY model were derived from the existing geographic information system and satellite‐driven data sets. By comparing the prerestoration and postrestoration simulation results under the identical atmospheric conditions, a more complete environmental assessment of the restoration project was made. Model results indicate that the flood peak may be reduced by 10–20% during the wet year and the baseflow may be enhanced by 10–20% during the following dry seasons (summer to fall) in the postrestoration condition. The model results also showed that the hydrologic impact of the land management associated with the restoration mitigates bank erosion and sediment discharge during winter storm events. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

ISA model: A constitutive model for soils with yield surface in the intergranular strain space          下载免费PDF全文

W. Fuentes  Th. Triantafyllidis 《国际地质力学数值与分析法杂志》2015,39(11):1235-1254

In this article, a new constitutive model for soils is proposed. It is formulated by means of plasticity, but in contrast to the precedent works, it presents a yield function describing a surface within the intergranular strain space. This latter is a state variable providing information of the recent strain history. An expression for the plastic strain rate has been proposed to guarantee the stress rate continuity. Under the application of medium or large strain amplitudes, the constitutive equation becomes independent of the intergranular strain and delivers a mathematical structure similar to some Karlsruhe hypoplastic models. Some simulations of monotonic and cyclic triaxial test are provided to evaluate and analyze the model performance. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Assessing the reliability of probabilistic flood inundation model predictions          下载免费PDF全文

Elisabeth Stephens  Paul Bates 《水文研究》2015,29(19):4264-4283

An ability to quantify the reliability of probabilistic flood inundation predictions is a requirement not only for guiding model development but also for their successful application. Probabilistic flood inundation predictions are usually produced by choosing a method of weighting the model parameter space, but previous study suggests that this choice leads to clear differences in inundation probabilities. This study aims to address the evaluation of the reliability of these probabilistic predictions. However, a lack of an adequate number of observations of flood inundation for a catchment limits the application of conventional methods of evaluating predictive reliability. Consequently, attempts have been made to assess the reliability of probabilistic predictions using multiple observations from a single flood event. Here, a LISFLOOD‐FP hydraulic model of an extreme (>1 in 1000 years) flood event in Cockermouth, UK, is constructed and calibrated using multiple performance measures from both peak flood wrack mark data and aerial photography captured post‐peak. These measures are used in weighting the parameter space to produce multiple probabilistic predictions for the event. Two methods of assessing the reliability of these probabilistic predictions using limited observations are utilized; an existing method assessing the binary pattern of flooding, and a method developed in this paper to assess predictions of water surface elevation. This study finds that the water surface elevation method has both a better diagnostic and discriminatory ability, but this result is likely to be sensitive to the unknown uncertainties in the upstream boundary condition. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献