首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Abstract

As watershed models become increasingly sophisticated and useful, there is a need to extend their applicability to locations where they cannot be calibrated or validated. A new methodology for the regionalization of a watershed model is introduced and evaluated. The approach involves calibration of a watershed model to many sites in a region, concurrently. Previous research that has sought to relate the parameters of monthly water balance models to physical drainage basin characteristics in a region has met with limited success. Previous studies have taken the two-step approach: (a) estimation of watershed model parameters at each site, followed by (b) attempts to relate model parameters to drainage basin characteristics. Instead of treating these two steps as independent, both steps are implemented concurrently. All watershed models in a region are calibrated simultaneously, with the dual objective of reproducing the behaviour of observed monthly streamflows and, additionally, to obtain good relationships between watershed model parameters and basin characteristics. The approach is evaluated using 33 basins in the southeastern region of the United States by comparing simulations using the regional models for three catchments which were not used to develop the regional regression equations. Although the regional calibration approach led to nearly perfect regional relationships between watershed model parameters and basin characteristics, these “improved” regional relationships did not result in improvements in the ability to model streamflow at ungauged sites. This experiment reveals that improvements in regional relationships between watershed model parameters and basin characteristics will not necessarily lead to improvements in the ability to calibrate a watershed model at an ungauged site.  相似文献   

2.
The basic aim of this paper is to formulate rigorous conservation equations for mass, momentum, energy and entropy for a watershed organized around the channel network. The approach adopted is based on the subdivision of the whole watershed into smaller discrete units, called representative elementary watersheds (REW), and the formulation of conservation equations for these REWs. The REW as a spatial domain is divided into five different subregions: (1) unsaturated zone; (2) saturated zone; (3) concentrated overland flow; (4) saturated overland flow; and (5) channel reach. These subregions all occupy separate volumina. Within the REW, the subregions interact with each other, with the atmosphere on top and with the groundwater or impermeable strata at the bottom, and are characterized by typical flow time scales.The balance equations are derived for water, solid and air phases in the unsaturated zone, water and solid phases in the saturated zone and only the water phase in the two overland flow zones and the channel. In this way REW-scale balance equations, and respective exchange terms for mass, momentum, energy and entropy between neighbouring subregions and phases, are obtained. Averaging of the balance equations over time allows to keep the theory general such that the hydrologic system can be studied over a range of time scales. Finally, the entropy inequality for the entire watershed as an ensemble of subregions is derived as constraint-type relationship for the development of constitutive relationships, which are necessary for the closure of the problem. The exploitation of the second law and the derivation of constitutive equations for specific types of watersheds will be the subject of a subsequent paper.  相似文献   

3.
This study presents a Geographic Information System (GIS)‐based distributed rainfall‐runoff model for simulating surface flows in small to large watersheds during isolated storm events. The model takes into account the amount of interception storage to be filled using a modified Merriam ( 1960 ) approach before estimating infiltration by the Smith and Parlange ( 1978 ) method. The mechanics of overland and channel flow are modelled by the kinematic wave approximation of the Saint Venant equations which are then numerically solved by the weighted four‐point implicit finite difference method. In this modelling the watershed was discretized into overland planes and channels using the algorithms proposed by Garbrecht and Martz ( 1999 ). The model code was first validated by comparing the model output with an analytical solution for a hypothetical plane. Then the model was tested in a medium‐sized semi‐forested watershed of Pathri Rao located in the Shivalik ranges of the Garhwal Himalayas, India. Initially, a local sensitivity analysis was performed to identify the parameters to which the model outputs like runoff volume, peak flow and time to peak flow are sensitive. Before going for model validation, calibration was performed using the Ordered‐Physics‐based Parameter Adjustment (OPPA) method. The proposed Physically Based Distributed (PBD) model was then evaluated both at the watershed outlet as well as at the internal gauging station, making this study a first of its kind in Indian watersheds. The results of performance evaluation indicate that the model has simulated the runoff hydrographs reasonably well within the watershed as well as at the watershed outlet with the same set of calibrated parameters. The model also simulates, realistically, the temporal variation of the spatial distribution of runoff over the watershed and the same has been illustrated graphically. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

5.
Recognizing that simple watershed conceptual models such as the Nash cascade ofn equal linear reservoirs continue to be reasonable means to approximate the Instantaneous Unit Hydrograph (IUH), it is natural to accept that random errors generated by climatological variability of data used in fitting an imprecise conceptual model will produce an IUH which is random itself. It is desirable to define the random properties of the IUH in a watershed in order to have a more realistic hydrologic application of this important function. Since in this case the IUH results from a series of differential equations where one or more of the uncertain parameters is treated in stochastic terms, then the statistical properties of the IUH are best described by the solution of the corresponding Stochastic Differential Equations (SDE's). This article attempts to present a methodology to derive the IUH in a small watershed by combining a classical conceptual model with the theory of SDE's. The procedure is illustrated with the application to the Middle Thames River, Ontario, Canada, and the model is verified by the comparison of the simulated statistical measures of the IUH with the corresponding observed ones with good agreement.  相似文献   

6.
Recognizing that simple watershed conceptual models such as the Nash cascade ofn equal linear reservoirs continue to be reasonable means to approximate the Instantaneous Unit Hydrograph (IUH), it is natural to accept that random errors generated by climatological variability of data used in fitting an imprecise conceptual model will produce an IUH which is random itself. It is desirable to define the random properties of the IUH in a watershed in order to have a more realistic hydrologic application of this important function. Since in this case the IUH results from a series of differential equations where one or more of the uncertain parameters is treated in stochastic terms, then the statistical properties of the IUH are best described by the solution of the corresponding Stochastic Differential Equations (SDE's). This article attempts to present a methodology to derive the IUH in a small watershed by combining a classical conceptual model with the theory of SDE's. The procedure is illustrated with the application to the Middle Thames River, Ontario, Canada, and the model is verified by the comparison of the simulated statistical measures of the IUH with the corresponding observed ones with good agreement.  相似文献   

7.
Representation of agricultural conservation practices with SWAT   总被引：5，自引：0，他引：5  
Results of modelling studies for the evaluation of water quality impacts of agricultural conservation practices depend heavily on the numerical procedure used to represent the practices. Herein, a method for the representation of several agricultural conservation practices with the Soil and Water Assessment Tool (SWAT) is developed and evaluated. The representation procedure entails identifying hydrologic and water quality processes that are affected by practice implementation, selecting SWAT parameters that represent the affected processes, performing a sensitivity analysis to ascertain the sensitivity of model outputs to selected parameters, adjusting the selected parameters based on the function of conservation practices, and verifying the reasonableness of the SWAT results. This representation procedure is demonstrated for a case study of a small agricultural watershed in Indiana in the Midwestern USA. The methods developed in the present work can be applied with other watershed models that employ similar underlying equations to represent hydrologic and water quality processes. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

8.
The quantification of the various components of hydrological processes in a watershed remains a challenging topic as the hydrological system is altered by internal and external drivers. Watershed models have become essential tools to understand the behaviour of a catchment under dynamic processes. In this study, a physically based watershed model called Soil Water Assessment Tool was used to understand the hydrologic behaviour of the Upper Tiber River Basin, Central Italy. The model was successfully calibrated and validated using observed weather and flow data for the period of 1963–1970 and 1971–1978, respectively. Eighteen parameters were evaluated, and the model showed high relative sensitivity to groundwater flow parameters than the surface flow parameters. An analysis of annual hydrological water balance was performed for the entire upper Tiber watershed and selected subbasins. The overall behaviour of the watershed was represented by three categories of parameters governing surface flow, subsurface flow and whole basin response. The base flow contribution has shown that 60% of the streamflow is from shallow aquifer in the subbasins. The model evaluation statistics that evaluate the agreement between the simulated and the observed streamflow at the outlet of a watershed and other three different subbasins has shown a coefficient of determination (R2) from 0.68 to 0.81 and a Nash–Sutcliffe efficiency (ENS) between 0.51 and 0.8 for the validation period. The components of the hydrologic cycle showed variation for dry and wet periods within the watershed for the same parameter sets. On the basis of the calibrated parameters, the model can be used for the prediction of the impact of climate and land use changes and water resources planning and management. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

9.
Model calibration is important for streamflow simulations using distributed hydrological models, especially in highland and cold areas of northwest China with scarce data. Quantitative analysis of water balance based on the accurate simulation is also essential for reasonably planning and managing water resources in these river basins facing a severe water shortage. In this study, a comprehensive method was proposed to calibrate the Soil and Water Assessment Tool (SWAT) model in the Yingluoxia watershed, upstream area of the Heihe River basin; it was based on multi-temporal, multi-variable and multi-site integrated drainage characteristics. Meanwhile a fresh approach of the parameter transferability and model validation was used by applying the set of calibrated parameters in its tributary to other area of the watershed. The results indicated that the method was effective and feasible; the values of Nash–Sutcliffe Efficiency (NSE) and Coefficient of Determination (r2) were greater than 0.81 and as high as 0.94 and the absolute values of the Percent Bias (PBIAS) were less than 2. Based the output of model the water balance in the Yingluoxia watershed was analyzed, that the mean annual precipitation, evapotranspiration, and discharge of the watershed from 1990 to 2000 were 491.8 mm, 334 mm, and 157.8 mm, respectively. The comprehensive calibration method based on multi-temporal, multi-variable and multi-site integrated drainage characteristics can better portray the hydrological processes of watershed and improve the model simulation; and the output of the model then provide a reliable reference for assessing and managing water resource of the watershed.  相似文献   

10.
— We propose an algorithm for local evaluation of weak anisotropy (WA) parameters from measurements of slowness vector components and/or of particle motions of q P waves at individual receivers in a borehole in a multi-azimuthal multiple-source offset VSP experiment. As a byproduct the algorithm yields approximate angular variation of q P-wave phase velocity. The formulae are derived under assumption of weak but arbitrary anisotropy and lateral inhomogeneity of the medium. The algorithm is thus independent of structural complexities between the source and the receiver. If complete slowness vector is determinable from observed data, then the information about polarization can be used as an independent additional constraint. If only the component of the slowness along the borehole can be determined from observations (which is mostly the case), the inversion without information about polarization is impossible. We present several systems of equations which can be used when different numbers of components of the slowness vector are available. The SVD algorithm is used to solve an overdetermined system of linear equations for WA parameters for two test examples of synthetic multi-azimuthal multiple-source offset VSP data. The system of equations results from approximate first-order perturbation equations for the slowness and polarization vectors of the q P wave. Analysis of singular values and of variances of WA parameters is used for the estimation of chances to recover the sought parameters. Effects of varying number of profiles with sources and of noise added to “observed” data are illustrated. An important observation is that although, due to insufficient data, we often cannot recover all individual WA parameters with sufficient accuracy, angular phase velocity variation can be recovered rather well.  相似文献   

11.
This paper presents an ecological–economic model for a lake and its watershed systems. We describe the linkage between the watershed system and the lake aquatic ecosystem and the modeling process. The lake–watershed system was divided into six subsystems: social system, economic system, terrestrial ecosystem, lake water system, pollutant system, and lake aquatic ecosystem. The model equations were constructed based on five main assumptions. The Lake Qionghai watershed in southwestern China, which is undergoing rapid eutrophication, was used as a case study. The targeted goals for total phosphorus (TP) and chlorophyll a (Chl a) concentrations in the lake in 2015 are 0.025 and 10.0 mg m−3, respectively. We present two scenarios from 2004 to 2015 based on the ecological–economic model. In both scenarios, the TP and Chl a concentrations in the lake are predicted to increase under the effects of watershed pressures and the targeted goals cannot be met. The application of techniques to reduce pollutants loading and the corresponding pollutants reductions are reflected again in the constructed model. The model predicts that TP and Chl a concentrations will decrease to 0.024 and 7.71 mg m−3, respectively, which meet the targeted thresholds. The model results provide directions for local government management of watersheds and lake aquatic ecosystem restoration.  相似文献   

12.
Mineral weathering rates and a forest macronutrient uptake stoichiometry were determined for the forested, metabasaltic Hauver Branch watershed in north‐central Maryland, USA. Previous studies of Hauver Branch have had an insufficient number of analytes to permit determination of rates of all the minerals involved in chemical weathering, including biomass. More equations in the mass‐balance matrix were added using existing mineralogic information. The stoichiometry of a deciduous biomass term was determined using multi‐year weekly to biweekly stream‐water chemistry for a nearby watershed, which drains relatively unreactive quartzite bedrock. At Hauver Branch, calcite hosts ~38 mol% of the calcium ion (Ca2+) contained in weathering minerals, but its weathering provides ~90% of the stream water Ca2+. This occurs in a landscape with a regolith residence time of more than several Ka (kiloannum). Previous studies indicate that such old regolith does not typically contain dissolving calcite that affects stream Ca2+/Na+ ratios. The relatively high calcite dissolution rate likely reflects dissolution of calcite in fractures of the deep critical zone. Of the carbon dioxide (CO2) consumed by mineral weathering, calcite is responsible for approximately 27%, with the silicate weathering consumption rate far exceeding that of the global average. The chemical weathering of mafic terrains in decaying orogens thus may be capable of influencing global geochemical cycles, and therefore, climate, on geological timescales. Based on carbon‐balance calculations, atmospheric‐derived sulfuric acid is responsible for approximately 22% of the mineral weathering occurring in the watershed. Our results suggest that rising air temperatures, driven by global warming and resulting in higher precipitation, will cause the rate of chemical weathering in the Hauver Branch watershed to increase until a threshold temperature is reached. Beyond the threshold temperature, increased recharge would produce a shallower groundwater table and reduced chemical weathering rates. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

13.
The time of concentration (Tc) is a fundamental parameter in the design of hydrological projects for watersheds. In this study a graphical methodology is described for estimating Tc in a watershed, and this is applied to 17 rainfall–runoff events from a rural watershed located near the capital city of Mato Grosso do Sul State, in the Brazilian Cerrado. The Tc values obtained through the graphical method were compared to Tc values estimated using 20 equations from various references. The equations were selected by considering those that were not developed using data for urban watersheds, and the results of the graphical method were compared to those derived by applying the equations to sub-basin data. The graphical method was reliable in determining Tc, and Ventura’s equation was found to present the best performance for a rural watershed in a tropical climate region.  相似文献   

14.
Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤Mw≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (Mw=6.8) and 1999 Chamoli earthquake (Mw=6.4) from Garhwal Himalaya region.  相似文献   

15.
Gully morphology characteristics can be used effectively to describe the status of gully development. The Chabagou watershed, located in the hilly‐gully region of the Loess Plateau in China, was selected to investigate gully morphological characteristics using a 3D laser scanning technique (LIDAR). Thirty‐one representative gullies located at different watershed locations and gully orders were chosen to quantitatively describe gully morphology and establish empirical equations for estimating gully volume based on gully length and gully surface area. Images and point cloud data for the 31 gullies were collected, and digital elevation models (DEMs) with 10‐cm resolution were generated. ArcGIS 10.1 was then used to extract fundamental gully morphological parameters covering gully length (L), gully width (WT) and gully depth (D), and some derivative morphological parameters, including gully head curvature (C), gully width–depth ratio (w/d), gully bottom‐to‐top width ratio (WB/WT), gully surface area (Ag) and gully volume (Vg). The results indicated that gullies in the upper watershed and the second order were more developed based on their high values of gully head curvature. The potential for gully development increased from the second order to the fourth order. Within the same gully orders, gullies in the lower watershed were more active with more development potential. A method for differentiating between gully head and gully sidewalls based on the gully head curvature value was proposed with a mean relative error of 8.77%. U‐shaped cross‐sections were widely distributed in the upper watershed and upper positions of a gully, while V‐shaped cross‐sections were widely distributed in the lower watershed and lower positions of a gully. V–L and V–Ag empirical equations with acceptable accuracy were established and can be used to estimate gully erosion in the Loess hilly‐gully region. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

16.
The hydrology of Quebec, Canada, boreal fens is poorly documented. Many peatlands are located in watersheds with impounded rivers. In such cases, their presence influences reservoir inflows. In recent years, some fens have been subjected to an increase of their wet area, a sign that they may be evolving towards an aquatic ecosystem. This dynamic process is called aqualysis. This article presents the seasonal and monthly hydrological budgets of a small watershed including a highly aqualysed fen (James Bay region). The monitoring of precipitation (P), runoff (Q) and groundwater levels (WL) was conducted during the ice‐free season. Three semiempirical equations (Thornthwaite, Priestley–Taylor and Penman–Monteith) were used and compared to calculate potential evapotranspiration. The first two equations, having fewer parameters, estimate higher potential evapotranspiration values than the third equation. The use of pressure‐level gauges installed in wells, for the calculation of peatland water storage, is inconclusive. Swelling of peat, peat decomposition and plant composition could be responsible for nonnegligible amounts of absorbed water, which are not entirely accounted for by well levels. The estimation of peat matrix water storage is potentially the largest source of error and the limiting factor to calculate water balances in this environment. The results show that the groundwater level and the water storage vary depending on the season and especially after a heavy rainfall. Finally, the results illustrate the complexity of water routing through the site and thus raise several questions to be resolved in the future. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

17.
在地震子波非因果、混合相位的假设下,本文应用自回归滑动平均(ARMA)模型对地震子波进行参数化建模,并提出利用线性(矩阵方程法)和非线性(ARMA拟合方法)相结合的参数估计方式对该模型进行参数估计.在利用矩阵方程法确定模型参数范围的基础上,利用累积量拟合法精确估计参数.理论分析和仿真结果表明,该方式有较好的适应性：一方面提高了子波估计精度,避免单独使用矩阵方程法在短数据地震记录情况下可能带来的估计误差;另一方面提高了子波提取运算效率,降低了ARMA模型拟合方法参数范围确定的复杂性,避免了单纯使用滑动平均（MA）模型拟合法估计过多参数所导致的运算规模过大问题.初步应用结果表明该方法是有效可行的.  相似文献   

18.
The balance equations for mass and momentum, averaged over the scale of a watershed entity, need to be supplemented with constitutive equations relating flow velocities, pressure potential differences, as well as mass and force exchanges within and across the boundaries of a watershed. In this paper, the procedure for the derivation of such constitutive relationships is described in detail. This procedure is based on the method pioneered by Coleman and Noll through exploitation of the second law of thermodynamics acting as a constraint-type relationship. The method is illustrated by its application to some common situations occurring in real world watersheds. Thermodynamically admissible and physically consistent constitutive relationships for mass exchange terms among the subregions constituting the watershed (subsurface zones, overland flow regions, channel) are proposed. These constitutive equations are subsequently combined with equations of mass balance for the subregions. In addition, constitutive relationships for forces exchanged amongst the subregions are also derived within the same thermodynamic framework. It is shown that, after linearisation of the latter constitutive relations in terms of the velocity, a watershed-scale Darcy's law governing flow in the unsaturated and saturated zones can be obtained. For the overland flow, a second order constitutive relationship with respect to velocity is proposed for the momentum exchange terms, leading to a watershed-scale Chezy formula. For the channel network REW-scale Saint–Venant equations are derived. Thus, within the framework of this approach new relationships governing exchange terms for mass and momentum are obtained and, moreover, some well-known experimental results are derived in a rigorous manner.  相似文献   

19.
Generalized inversion of the S-wave amplitude spectra from the strong-motion network data in the East-Central Iran has been used to estimate simultaneously source parameters, site response and the S-wave attenuation (Qs). In this regard, 190 three-component records were used corresponded to 40 earthquakes with the magnitudes M3.5–M7.3. These earthquakes were recorded at 42 stations in the hypocentral distance range from 9 to 200 km. The inverse problem was solved in 20 logarithmically equally spaced points in the frequency band from 0.4 to 15 Hz. The frequency-dependent site amplification was imposed, as a constraint, on two reference site responses in order to remove the undetermined degree of freedom in the inversion and obtain a unique inverse solution. Also, a geometrical spreading factor was assumed for removing the trade-off between geometrical spreading and anelastic attenuation. Different source parameters, such as seismic moment (M0), seismic energy (Es), corner frequency (fc) and Brune stress drop (Δσ), were estimated for each event by fitting an ω2 model to the spectra obtained from the inversion. The stress drop values of earthquakes, obtained in this research, are in good agreement with those of other studies. Also average site response values were correlated to the average shear wave velocities in the uppermost 30 m, in high and low frequency bands. The peak frequencies of site amplifications, estimated by the generalized inversion method, where in good agreement with those of horizontal to vertical (H/V) spectral ratios for the S-wave portion of records. However, no perfect matching in amplitude was obtained due to the deficiencies of the H/V ratio technique. By supposing a free shape for Q factor, a frequency dependent function was found, the logarithm of which could be approximated by a linear function, Q(f)=151f0.75. The uncertainties of model parameters have been evaluated by covariance matrix of least-square fit. The residuals were also analyzed in order to assess the validity of the model. The analysis of residuals with respect to magnitude and distance indicates that they are distributed normally with approximately zero mean. The robustness of the results has been studied concerning their sensitivities to the omission of different datasets, selected randomly from original database. The results obtained here can be used in predicting ground-motion parameters applying stochastic methods.  相似文献   

20.
To recognize the geographical characteristics of the landslide areas will be helpful for the watershed management in the reservoir watershed.According to the quantitative analysis,we'll take different scores and weighting for the potential parameters of the landslide areas in the Tsengwen reservoir watershed,and in the meanwhile,we'll extract the different factors,including the slope,aspect,altitude,soil and geological textures etc.,and the results shown as maximum one-day rainfall,ratio of forests and average relief is the most affecting parameters on the potential risk map of landslide areas.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号