首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
With the recent development of distributed hydrological models, the use of multi‐site observed data to evaluate model performance is becoming more common. Distributed hydrological model have many advantages, and at the same time, it also faces the challenge to calibrate over‐do parameters. As a typical distributed hydrological model, problems also exist in Soil and Water Assessment Tool (SWAT) parameter calibration. In the paper, four different uncertainty approaches – Particle Swarm Optimization (PSO) techniques, Generalized Likelihood Uncertainty Estimation (GLUE), Sequential Uncertainty Fitting algorithm (SUFI‐2) and Parameter Solution (PARASOL) – are taken to a comparative study with the SWAT model applied in Peace River Basin, central Florida. In our study, the observed river discharge data used in SWAT model calibration were collected from the three gauging stations at the main tributary of the Peace River. Behind these approaches, there is a shared philosophy; all methods seek out many parameter set to fit the uncertainties due to the non‐uniqueness in model parameter evaluation. On the basis of the statistical results of four uncertainty methods, difficulty level of each method, the number of runs and theoretical basis, the reasons that affected the accuracy of simulation were analysed and compared. Furthermore, for the four uncertainty method with SWAT model in the study area, the pairwise correlation between parameters and the distributions of model fit summary statistics computed from the sampling over the behavioural parameter and the entire model calibration parameter feasible spaces were identified and examined. It provided additional insight into the relative identifiability of the four uncertainty methods Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
Lake Tana Basin is of significant importance to Ethiopia concerning water resources aspects and the ecological balance of the area. Many years of mismanagement, wetland losses due to urban encroachment and population growth, and droughts are causing its rapid deterioration. The main objective of this study was to assess the performance and applicability of the soil water assessment tool (SWAT) model for prediction of streamflow in the Lake Tana Basin, so that the influence of topography, land use, soil and climatic condition on the hydrology of Lake Tana Basin can be well examined. The physically based SWAT model was calibrated and validated for four tributaries of Lake Tana. Sequential uncertainty fitting (SUFI‐2), parameter solution (ParaSol) and generalized likelihood uncertainty estimation (GLUE) calibration and uncertainty analysis methods were compared and used for the set‐up of the SWAT model. The model evaluation statistics for streamflows prediction shows that there is a good agreement between the measured and simulated flows that was verified by coefficients of determination and Nash Sutcliffe efficiency greater than 0·5. The hydrological water balance analysis of the basin indicated that baseflow is an important component of the total discharge within the study area that contributes more than the surface runoff. More than 60% of losses in the watershed are through evapotranspiration. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

3.
Streams play an important role in linking the land with lakes. Nutrients released from agricultural or urban sources flow via streams to lakes, causing water quality deterioration and eutrophication. Therefore, accurate simulation of streamflow is helpful for water quality improvement in lake basins. Lake Dianchi has been listed in the ‘Three Important Lakes Restoration Act’ in China, and the degradation of its water quality has been of great concern since the 1980s. To assist environmental decision making, it is important to assess and predict hydrological processes at the basin scale. This study evaluated the performance of the soil and water assessment tool (SWAT) and the feasibility of using this model as a decision support tool for predicting streamflow in the Lake Dianchi Basin. The model was calibrated and validated using monthly observed streamflow values at three flow stations within the Lake Dianchi Basin through application of the sequential uncertainty fitting algorithm (SUFI‐2). The results of the autocalibration method for calibrating and the prediction uncertainty from different sources were also examined. Together, the p‐factor (the percentage of measured data bracketed by 95% prediction of uncertainty, or 95PPU) and the r‐factor (the average thickness of the 95PPU band divided by the standard deviation of the measured data) indicated the strength of the calibration and uncertainty analysis. The results showed that the SUFI‐2 algorithm performed better than the autocalibration method. Comparison of the SUFI‐2 algorithm and autocalibration results showed that some snowmelt factors were sensitive to model output upstream at the Panlongjiang flow station. The 95PPU captured more than 70% of the observed streamflow at the three flow stations. The corresponding p‐factors and r‐factors suggested that some flow stations had relatively large uncertainty, especially in the prediction of some peak flows. Although uncertainty existed, statistical criteria including R2 and Nash–Sutcliffe efficiency were reasonably determined. The model produced a useful result and can be used for further applications. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

4.
Abstract

The SWAT model was tested to simulate the streamflow of two small Mediterranean catchments (the Vène and the Pallas) in southern France. Model calibration and prediction uncertainty were assessed simultaneously by using three different techniques (SUFI-2, GLUE and ParaSol). Initially, a sensitivity analysis was conducted using the LH-OAT method. Subsequent sensitive parameter calibration and SWAT prediction uncertainty were analysed by considering, firstly, deterministic discharge data (assuming no uncertainty in discharge data) and secondly, uncertainty in discharge data through the development of a methodology that accounts explicitly for error in the rating curve (the stage?discharge relationship). To efficiently compare the different uncertainty methods and the effect of the uncertainty of the rating curve on model prediction uncertainty, common criteria were set for the likelihood function, the threshold value and the number of simulations. The results show that model prediction uncertainty is not only case-study specific, but also depends on the selected uncertainty analysis technique. It was also found that the 95% model prediction uncertainty interval is wider and more successful at encompassing the observations when uncertainty in the discharge data is considered explicitly. The latter source of uncertainty adds additional uncertainty to the total model prediction uncertainty.
Editor D. Koutsoyiannis; Associate editor D. Gerten

Citation Sellami, H., La Jeunesse, I., Benabdallah, S., and Vanclooster, M., 2013. Parameter and rating curve uncertainty propagation analysis of the SWAT model for two small Mediterranean watersheds. Hydrological Sciences Journal, 58 (8), 1635?1657.  相似文献   

5.
Abandoned underground mines (AUM) have caused dramatic environmental effects that are closely linked to regional sustainability. This paper explores the potential hydrological impact of AUM in the Monday Creek Watershed, a typically mined area in Appalachian region, using the Soil and Water Assessment Tool (SWAT 2005) model and Sequential Uncertainty Fitting (SUFI‐2), calibrated at both the global and local scales. The locally calibrated model better incorporates those key parameters relevant to AUM for specific sub‐basins and hydrologic response units. Data from the years 2003–2004 were used for calibration and 2005–2006 for validation. The results were quite satisfactory; both the coefficient of determination (R2) and the Nash–Sutcliffe efficiency statistic were over 0.80. The potential influences of AUM were assessed by modelling an alternative scenario assuming no AUM for the period 2003–2009. Results show that the hydrological process of lateral subsurface flow plays a dominant role in linking AUM to overall watershed hydrology. The potential hydrological impact of AUM is an increased annual lateral flow of 82.1%, and a decrease in annual surface flow by 15%, leading to an increase of 16.9% in annual water yield for the Monday Creek Watershed. The seasonal fluctuation of water yield has a similar trend to lateral flow, decreasing from March to August and increasing from August to January. Higher volume, higher flow peaks and higher recession constants characterized the hydrograph of daily streamflow from AUM. The results indicate that more emphasis should be put on lateral flow for further study of acid mine drainage and flooding control in those watersheds with AUM. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

6.
Understanding hydrological processes at catchment scale through the use of hydrological model parameters is essential for enhancing water resource management. Given the difficulty of using lump parameters to calibrate distributed catchment hydrological models in spatially heterogeneous catchments, a multiple calibration technique was adopted to enhance model calibration in this study. Different calibration techniques were used to calibrate the Soil and Water Assessment Tool (SWAT) model at different locations along the Logone river channel. These were: single-site calibration (SSC); sequential calibration (SC); and simultaneous multi-site calibration (SMSC). Results indicate that it is possible to reveal differences in hydrological behavior between the upstream and downstream parts of the catchment using different parameter values. Using all calibration techniques, model performance indicators were mostly above the minimum threshold of 0.60 and 0.65 for Nash Sutcliff Efficiency (NSE) and coefficient of determination (R 2) respectively, at both daily and monthly time-steps. Model uncertainty analysis showed that more than 60% of observed streamflow values were bracketed within the 95% prediction uncertainty (95PPU) band after calibration and validation. Furthermore, results indicated that the SC technique out-performed the other two methods (SSC and SMSC). It was also observed that although the SMSC technique uses streamflow data from all gauging stations during calibration and validation, thereby taking into account the catchment spatial variability, the choice of each calibration method will depend on the application and spatial scale of implementation of the modelling results in the catchment.  相似文献   

7.
ABSTRACT

Climate models and hydrological parameter uncertainties were quantified and compared while assessing climate change impacts on monthly runoff and daily flow duration curve (FDC) in a Mediterranean catchment. Simulations of the Soil and Water Assessment Tool (SWAT) model using an ensemble of behavioural parameter sets derived from the Generalized Likelihood Uncertainty Estimation (GLUE) method were approximated by feed-forward artificial neural networks (FF-NN). Then, outputs of climate models were used as inputs to the FF-NN models. Subsequently, projected changes in runoff and FDC were calculated and their associated uncertainty was partitioned into climate model and hydrological parameter uncertainties. Runoff and daily discharge of the Chiba catchment were expected to decrease in response to drier and warmer climatic conditions in the 2050s. For both hydrological indicators, uncertainty magnitude increased when moving from dry to wet periods. The decomposition of uncertainty demonstrated that climate model uncertainty dominated hydrological parameter uncertainty in wet periods, whereas in dry periods hydrological parametric uncertainty became more important.
Editor M.C. Acreman; Associate editor S. Kanae  相似文献   

8.
9.
This paper addresses the application of a data‐based mechanistic (DBM) modelling approach using transfer function models (TFMs) with non‐linear rainfall filtering to predict runoff generation from a semi‐arid catchment (795 km2) in Tanzania. With DBM modelling, time series of rainfall and streamflow were allowed to suggest an appropriate model structure compatible with the data available. The model structures were evaluated by looking at how well the model fitted the data, and how well the parameters of the model were estimated. The results indicated that a parallel model structure is appropriate with a proportion of the runoff being routed through a fast flow pathway and the remainder through a slow flow pathway. Finally, the study employed a Generalized Likelihood Uncertainty Estimation (GLUE) methodology to evaluate the parameter sensitivity and predictive uncertainty based on the feasible parameter ranges chosen from the initial analysis of recession curves and calibration of the TFM. Results showed that parameters that control the slow flow pathway are relatively more sensitive than those that control the fast flow pathway of the hydrograph. Within the GLUE framework, it was found that multiple acceptable parameter sets give a range of predictions. This was found to be an advantage, since it allows the possibility of assessing the uncertainty in predictions as conditioned on the calibration data and then using that uncertainty as part of the decision‐making process arising from any rainfall‐runoff modelling project. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

10.
Many methods developed for calibration and validation of physically based distributed hydrological models are time consuming and computationally intensive. Only a small set of input parameters can be optimized, and the optimization often results in unrealistic values. In this study we adopted a multi‐variable and multi‐site approach to calibration and validation of the Soil Water Assessment Tool (SWAT) model for the Motueka catchment, making use of extensive field measurements. Not only were a number of hydrological processes (model components) in a catchment evaluated, but also a number of subcatchments were used in the calibration. The internal variables used were PET, annual water yield, daily streamflow, baseflow, and soil moisture. The study was conducted using an 11‐year historical flow record (1990–2000); 1990–94 was used for calibration and 1995–2000 for validation. SWAT generally predicted well the PET, water yield and daily streamflow. The predicted daily streamflow matched the observed values, with a Nash–Sutcliffe coefficient of 0·78 during calibration and 0·72 during validation. However, values for subcatchments ranged from 0·31 to 0·67 during calibration, and 0·36 to 0·52 during validation. The predicted soil moisture remained wet compared with the measurement. About 50% of the extra soil water storage predicted by the model can be ascribed to overprediction of precipitation; the remaining 50% discrepancy was likely to be a result of poor representation of soil properties. Hydrological compensations in the modelling results are derived from water balances in the various pathways and storage (evaporation, streamflow, surface runoff, soil moisture and groundwater) and the contributions to streamflow from different geographic areas (hill slopes, variable source areas, sub‐basins, and subcatchments). The use of an integrated multi‐variable and multi‐site method improved the model calibration and validation and highlighted the areas and hydrological processes requiring greater calibration effort. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

11.
Soils affect the distribution of hydrological processes by partitioning precipitation into different components of the water balance. Therefore, understanding soil-water dynamics at a catchment scale remains imperative to future water resource management. In this study, the value of hydropedological insights was examined to calibrate a processes-based model. Soil morphology was used as soft data to assist in the calibration of the Soil Water Assessment Tool (SWAT+) model at five different catchment scales (48, 56, 174, 674, and 2421 km2) in the Sabie River catchment, South Africa. The aim of this study was to calibrate the SWAT+ model to accurately simulate long-term monthly streamflow predictions as well as to reflect internal soil hydrological processes using a procedure focusing on hydropedology as a calibration tool in a multigauge system. Results indicated that calibration improved streamflow predictions where R2 improved by 2%–8%. Nash-Sutcliffe Efficiency (NSE) improved from negative correlations to values exceeding 0.5 at four of the five catchment scales compared to the uncalibrated model. Results confirm that soil mapping units can be calibrated individually within SWAT+ to improve the representation of hydrological processes. Particularly, the spatial linkage between hydropedology and hydrological processes, which is captured within the soil map of the catchment, can be adequately reflected within the model simulations after calibration. This research will lead to an improved understanding of hydropedology as soft data to improve hydrological modelling accuracy.  相似文献   

12.
Provision of reliable scientific support to socio‐economic development and eco‐environmental conservation is challenged by complexities of irregular nonlinearities, data uncertainties, and multivariate dependencies of hydrological systems in the Three Gorges Reservoir (TGR) region, China. Among them, the irregular nonlinearities mainly represent unreliability of regular functions for robust simulation of highly complicated relationships between variables. Based on the proposed discrete principal‐monotonicity inference (DPMI) approach, streamflow generation in the Xingshan Watershed, a representative watershed in this region, is examined. Based on system characterization, predictor identification, and streamflow distribution transformation, DPMI parameters are calibrated through a two‐stage strategy. Results indicate that the modelling efficiency of DPMI is satisfactory for streamflow simulation under these complexities. The distribution transformation method and the two‐stage calibration strategy can deal with non‐normality of streamflow and temporally unstable accuracy of hydrological models, respectively. The DPMI process and results reveal that both streamflow uncertainty and its rising tendency increase with flow levels. The dominant driving forces of streamflow generation are daily lowest temperature and daily cumulative precipitation in consideration of performances in global and local scales. The temporal heterogeneity of local significances to streamflow is insignificant for meteorological conditions. There is significant nonlinearity between meteorological conditions and streamflow and dependencies among meteorological conditions. The generation mechanism of low flows is more complicated than medium flows and high flows. The DPMI approach can facilitate improving robustness of hydro‐system analysis studies in the Xingshan Watershed or the TGR region. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

13.
Accurate analysis of water flow pathways from rainfall to streams is critical for simulating water use, climate change impact, and contaminants transport. In this study, we developed a new scheme to simultaneously calibrate surface flow (SF) and baseflow (BF) simulations of soil and water assessment tool (SWAT) by combing evolutionary multi‐objective optimization (EMO) and BF separation techniques. The application of this scheme demonstrated pronounced trade‐off of SWAT's performance on SF and BF simulations. The simulated major water fluxes and storages variables (e.g. soil moisture, evapotranspiration, and groundwater) using the multiple parameters from EMO span wide ranges. Uncertainty analysis was conducted by Bayesian model averaging of the Pareto optimal solutions. The 90% confidence interval (CI) estimated using all streamflows substantially overestimate the uncertainty of low flows on BF days while underestimating the uncertainty of high flows on SF days. Despite using statistical criteria calculated based on streamflow for model selection, it is important to conduct diagnostic analysis of the agreement of SWAT behaviour and actual watershed dynamics. The new calibration technique can serve as a useful tool to explore the trade‐off between SF and BF simulations and provide candidates for further diagnostic assessment and model identification. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

14.
Keith Beven  Andrew Binley 《水文研究》2014,28(24):5897-5918
This paper reviews the use of the Generalized Likelihood Uncertainty Estimation (GLUE) methodology in the 20 years since the paper by Beven and Binley in Hydrological Processes in (1992), which is now one of the most highly cited papers in hydrology. The original conception, the on‐going controversy it has generated, the nature of different sources of uncertainty and the meaning of the GLUE prediction uncertainty bounds are discussed. The hydrological, rather than statistical, arguments about the nature of model and data errors and uncertainties that are the basis for GLUE are emphasized. The application of the Institute of Hydrology distributed model to the Gwy catchment at Plynlimon presented in the original paper is revisited, using a larger sample of models, a wider range of likelihood evaluations and new visualization techniques. It is concluded that there are good reasons to reject this model for that data set. This is a positive result in a research environment in that it requires improved models or data to be made available. In practice, there may be ethical issues of using outputs from models for which there is evidence for model rejection in decision making. Finally, some suggestions for what is needed in the next 20 years are provided. © 2013 The Authors. Hydrological Processes published by John Wiley & Sons, Ltd.  相似文献   

15.
ABSTRACT

In this study, the distributed catchment-scale model, DiCaSM, was applied on five catchments across the UK. Given its importance, river flow was selected to study the uncertainty in streamflow prediction using the Generalized Likelihood Uncertainty Estimation (GLUE) methodology at different timescales (daily, monthly, seasonal and annual). The uncertainty analysis showed that the observed river flows were within the predicted bounds/envelope of 5% and 95% percentiles. These predicted river flow bounds contained most of the observed river flows, as expressed by the high containment ratio, CR. In addition to CR, other uncertainty indices – bandwidth B, relative bandwidth RB, degrees of asymmetry S and T, deviation amplitude D, relative deviation amplitude RD and the R factor – also indicated that the predicted river flows have acceptable uncertainty levels. The results show lower uncertainty in predicted river flows when increasing the timescale from daily to monthly to seasonal, with the lowest uncertainty associated with annual flows.  相似文献   

16.
Agricultural pollutant runoff is a major source of water contamination in California's Sacramento River watershed where 8500 km2 of agricultural land influences water quality. The Soil and Water Assessment Tool (SWAT) hydrology, sediment, nitrate and pesticide transport components were assessed for the Sacramento River watershed. To represent flood conveyance in the area, the model was improved by implementing a flood routing algorithm. Sensitivity/uncertainty analyses and multi‐objective calibration were incorporated into the model application for predicting streamflow, sediment, nitrate and pesticides (chlorpyrifos and diazinon) at multiple watershed sites from 1992 to 2008. Most of the observed data were within the 95% uncertainty interval, indicating that the SWAT simulations were capturing the uncertainties that existed, such as model simplification, observed data errors and lack of agricultural management data. The monthly Nash–Sutcliffe coefficients at the watershed outlet ranged from 0.48 to 0.82, indicating that the model was able to successfully predict streamflow and agricultural pollutant transport after calibration. Predicted sediment loads were highly correlated to streamflow, whereas nitrate, chlorpyrifos and diazinon were moderately correlated to streamflow. This indicates that timing of agricultural management operations plays a role in agricultural pollutant runoff. Best management practices, such as pesticide use limits during wet seasons, could improve water quality in the Sacramento River watershed. The calibrated model establishes a modelling framework for further studies of hydrology, water quality and ecosystem protection in the study area. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

17.
This paper proposes a new orientation to address the problem of hydrological model calibration in ungauged basin. Satellite radar altimetric observations of river water level at basin outlet are used to calibrate the model, as a surrogate of streamflow data. To shift the calibration objective, the hydrological model is coupled with a hydraulic model describing the relation between streamflow and water stage. The methodology is illustrated by a case study in the Upper Mississippi Basin using TOPEX/Poseidon (T/P) satellite data. The generalized likelihood uncertainty estimation (GLUE) is employed for model calibration and uncertainty analysis. We found that even without any streamflow information for regulating model behavior, the calibrated hydrological model can make fairly reasonable streamflow estimation. In order to illustrate the degree of additional uncertainty associated with shifting calibration objective and identifying its sources, the posterior distributions of hydrological parameters derived from calibration based on T/P data, streamflow data and T/P data with fixed hydraulic parameters are compared. The results show that the main source is the model parameter uncertainty. And the contribution of remote sensing data uncertainty is minor. Furthermore, the influence of removing high error satellite observations on streamflow estimation is also examined. Under the precondition of sufficient temporal coverage of calibration data, such data screening can eliminate some unrealistic parameter sets from the behavioral group. The study contributes to improve streamflow estimation in ungauged basin and evaluate the value of remote sensing in hydrological modeling. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

18.
19.
A. Montenegro  R. Ragab 《水文研究》2010,24(19):2705-2723
Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R2 of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

20.
Understanding the impacts of land‐use changes on hydrology at the watershed scale can facilitate development of sustainable water resource strategies. This paper investigates the hydrological effects of land‐use change in Zanjanrood basin, Iran. The water balance was simulated using the Soil and Water Assessment Tool (AVSWAT2000). Model calibration and uncertainty analysis were performed with sequential uncertainty fitting (SUFI‐2). Simulation results from January 1998 to December 2002 were used for parameter calibration, and then the model was validated for the period of January 2003 to December 2004. The predicted monthly streamflow matched the observed values: during calibration the correlation coefficient was 0·86 and the Nash–Sutcliffe coefficient 0·79, compared with 0·80 and 0·79, respectively, during validation. The model was used to simulate the main components of the hydrological cycle, in order to study the effects of land‐use changes in 1967, 1994 and 2007. The study reveals that during 1967 a 34·5% decrease of grassland with concurrent increases of shrubland (13·9%), rain‐fed agriculture (12·1%), bare ground (5·5%) irrigated agriculture (2·2%), and urban area (0·7%) led to a 33% increase in the amount of surface runoff and a 22% decrease in the groundwater recharge. Furthermore, the area of sub‐basins that was influenced by high runoff (14–28 mm) increased. The results indicate that the hydrological response to overgrazing and the replacing of rangelands (grassland and shrubland) with rain‐fed agriculture and bare ground (badlands) is nonlinear and exhibits a threshold effect. The runoff rises dramatically when more than 60% of the rangeland is removed. For groundwater this threshold lies at an 80% decrease in rangeland. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

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

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