首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Modelling the hydrology of North American Prairie watersheds is complicated because of the existence of numerous landscape depressions that vary in storage capacity. The Soil and Water Assessment Tool (SWAT) is a widely applied model for long‐term hydrological simulations in watersheds dominated by agricultural land uses. However, several studies show that the SWAT model has had limited success in handling prairie watersheds. In past works using SWAT, landscape depression storage heterogeneity has largely been neglected or lumped. In this study, a probability distributed model of depression storage is introduced into the SWAT model to better handle landscape storage heterogeneity. The work utilizes a probability density function to describe the spatial heterogeneity of the landscape depression storages that was developed from topographic characteristics. The integrated SWAT–PDLD model is tested using datasets for two prairie depression dominated watersheds in Canada: the Moose Jaw River watershed, Saskatchewan; and the Assiniboine River watershed, Saskatchewan. Simulation results were compared to observed streamflow using graphical and multiple statistical criterions. Representation of landscape depressions within SWAT using a probability distribution (SWAT–PDLD) provides improved estimations of streamflow for large prairie watersheds in comparison to results using a lumped, single storage approach. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
The Soil and Water Assessment Tool (SWAT) is a physically‐based hydrologic model developed for agricultural watersheds, which has been infrequently validated for forested watersheds, particularly those with deep overwinter snow accumulation and abundant lakes and wetlands. The goal of this study was to determine the applicability of SWAT for modelling streamflow in two watersheds of the Ontonagon River basin of northern Michigan which differ in proportion of wetland and lake area. The forest‐dominated East Branch watershed contains 17% wetland and lake area, whereas the wetland/lake‐dominated Middle Branch watershed contains 26% wetland and lake area. The specific objectives were to: (1) calibrate and validate SWAT models for the East Branch and Middle Branch watersheds to simulate monthly stream flow, and (2) compare the effects of wetland and lake abundance on the magnitude and timing of streamflow. Model calibration and validation was satisfactory, as determined by deviation of discharge D and Nash and Sutcliffe coefficient values E that compared simulated monthly mean discharge versus measured monthly mean discharge. Streamflow simulation discrepancies occurred during summer and fall months and dry years. Several snow melting parameters were found to be critical for the SWAT simulation: TIMP (snow temperature lag factor) and SMFMX and SMFMN (melting factors). Snow melting parameters were not transferable between adjacent watersheds. Differences in seasonal pattern of long‐term monthly streamflow were found, with the forest‐dominated watershed having a higher peak flow during April but a lower flow during the remainder of the year in comparison to the wetland and lake‐dominated watershed. The results suggested that a greater proportion of wetland and lake area increases the capacity of a watershed to impound surface runoff and to delay storm and snow melting events. Representation of wetlands and lakes in a watershed model is required to simulate monthly stream flow in a wetland/lake‐dominated watershed. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

3.
This paper describes the results of benchmark testing of land use change impact on direct runoff using Soil Conservation Service-Curve Number (SCS-CN) model in two ungauged neighbouring urban watersheds (Çınar and Kadıyakuplu) in Istanbul, Turkey. To examine this impact, the model was applied to daily rainfall data using three different dated (1982, 1996 and 2012) hydrological soil groups and land use of the two ungauged urban watersheds. Finally, the impact of land use change and model performance were evaluated with the rainfall-runoff regression, the coefficient of determination and the NSE test using benchmark runoff data based on 1982 land use conditions. The results of the analysis indicate that the changing of land use types from natural surfaces to impervious surfaces has a significant impact on surface runoff. Additionally, remarkable spatial variations of the land use changes and their impact on the runoff in 1996 and 2012 were more detected in the Çınar watershed compared with the Kadıyakuplu watershed. The planning decision on land use of the watersheds, has vital role in these differences. The results of this research also reveal that change to intensive land use in urban watersheds has a significantly larger impact on runoff generation than those rainfall.  相似文献   

4.
Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

5.
  相似文献   

6.
Approaches to modeling the continuous hydrologic response of ungauged basins use observable physical characteristics of watersheds to either directly infer values for the parameters of hydrologic models, or to establish regression relationships between watershed structure and model parameters. Both these approaches still have widely discussed limitations, including impacts of model structural uncertainty. In this paper we introduce an alternative, model independent, approach to streamflow prediction in ungauged basins based on empirical evidence of relationships between watershed structure, climate and watershed response behavior. Instead of directly estimating values for model parameters, different hydrologic response behaviors of the watershed, quantified through model independent streamflow indices, are estimated and subsequently regionalized in an uncertainty framework. This results in expected ranges of streamflow indices in ungauged watersheds. A pilot study using 30 UK watersheds shows how this regionalized information can be used to constrain ensemble predictions of any model at ungauged sites. Dominant controlling characteristics were found to be climate (wetness index), watershed topography (slope), and hydrogeology. Main streamflow indices were high pulse count, runoff ratio, and the slope of the flow duration curve. This new approach provided sharp and reliable predictions of continuous streamflow at the ungauged sites tested.  相似文献   

7.
Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
Physiography and land cover determine the hydrologic response of watersheds to climatic events. However, vast differences in climate regimes and variation of landscape attributes among watersheds (including size) have prevented the establishment of general relationships between land cover and runoff patterns across broad scales. This paper addresses these difficulties by using power spectral analysis to characterize area‐normalized runoff patterns and then compare these patterns with landscape features among watersheds within the same physiographic region. We assembled long‐term precipitation and runoff data for 87 watersheds (first to seventh order) within the eastern Piedmont (USA) that contained a wide variety of land cover types, collected environmental data for each watershed, and compared the datasets using a variety of statistical measures. The effect of land cover on runoff patterns was confirmed. Urban‐dominated watersheds were flashier and had less hydrologic memory compared with forest‐dominated watersheds, whereas watersheds with high wetland coverage had greater hydrologic memory. We also detected a 10–15% urban threshold above which urban coverage became the dominant control on runoff patterns. When spectral properties of runoff were compared across stream orders, a threshold after the third order was detected at which watershed processes became dominant over precipitation regime in determining runoff patterns. Finally, we present a matrix that characterizes the hydrologic signatures of rivers based on precipitation versus landscape effects and low‐frequency versus high‐frequency events. The concepts and methods presented can be generally applied to all river systems to characterize multiscale patterns of watershed runoff. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

9.
Most semi‐distributed watershed water quality models divide the watershed into hydrologic response units (HRU) with no flow among them. This is problematic when watersheds are delineated to include variable source areas (VSAs) because it is the lateral flows from upslope areas to downslope areas that generate VSAs. Although hydrologic modellers have often successfully calibrated these types of models, there can still be considerable uncertainty in model results. In this paper, a topographic‐index‐based method is described and tested to distribute effective soil water holding capacity among HRUs, which can be subsequently adjusted using the watershed baseflow coefficient. The method is tested using a version of the Soil and Water Assessment Tool (SWAT) model that simulates VSA runoff and is applied to two watersheds: a New York State (NYS) watershed, and one in the head waters of the Blue Nile Basin (BNB) in Ethiopia. Daily streamflow predicted using effective soil water storage capacities based only on the topographic index were reassuringly accurate in both the NYS watershed (daily Nash Sutcliffe (E) = 0·73) and in the BNB (E = 0·70). Using the baseflow coefficient to adjust the effective soil water storage capacity only slightly improved streamflow predictions in NYS (E = 0·75) but substantially improved the BNB predictions (E = 0·80). By comparison, the standard SWAT model, which uses the traditional look‐up tables to determine a runoff curve number, performed considerably less accurately in un‐calibrated form (E = 0·51 for NYS and E = 0·45 for BNB), but improved substantially when explicitly calibrated to streamflow measurements (E = 0·76 for NYS and E = 0·67 for the BNB). The calibration method presented here provides a parsimonious, systematic approach to using established models in VSA watersheds that reduces the ambiguity inherent in model calibration. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
Precipitation is a key control on watershed hydrologic modelling output, with errors in rainfall propagating through subsequent stages of water quantity and quality analysis. Most watershed models incorporate precipitation data from rain gauges; higher‐resolution data sources are available, but they are associated with greater computational requirements and expertise. Here, we investigate whether the Multisensor Precipitation Estimator (MPE or Stage IV Next‐Generation Radar) data improve the accuracy of streamflow simulations using the Soil and Water Assessment Tool (SWAT), compared with rain gauge data. Simulated flows from 2002 to 2010 at five timesteps were compared with observed flows for four nested subwatersheds of the Neuse River basin in North Carolina (21‐, 203‐, 2979‐, and 10 100‐km2 watershed area), using a multi‐objective function, informal likelihood‐weighted calibration approach. Across watersheds and timesteps, total gauge precipitation was greater than radar precipitation, but radar data showed a conditional bias of higher rainfall estimates during large events (>25–50 mm/day). Model parameterization differed between calibrations with the two datasets, despite the fact that all watershed characteristics were the same across simulation scenarios. This underscores the importance of linking calibration parameters to realistic processes. SWAT simulations with both datasets underestimated median and low flows, whereas radar‐based simulations were more accurate than gauge‐based simulations for high flows. At coarser timesteps, differences were less pronounced. Our results suggest that modelling efforts in watersheds with poor rain gauge coverage can be improved with MPE radar data, especially at short timesteps. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.  相似文献   

11.
Design flood estimation in ungauged catchments is of great importance in hydrologic practice especially where there is no available data about streamflow. Except the watershed of Anseghmir who is equipped with a gauge station, all the other watersheds are ungauged catchments. The use of frequency analysis of series of rainfall and streamflow is very important for the characterization of the hydrologic resources of the Upper Moulouya. The region has a semiarid climate that requires a good knowledge of the watershed's potential water to assist policy makers in forecasting extreme events, managing water resources and decision making. The frequency analysis was used to determine the design flood of different return periods. The results obtained are used in Gradex method to estimate the hydrologic variables of each subcatchment of the Upper Moulouya. Once the hydrologic study is completed, a principal components analysis was made to highlight the affinities between the different subcatchments and to deduce the hydrologic and hydrographic parameters that better characterize them. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

12.
Using a mass balance algorithm, this study develops an extension module that can be embedded in the commonly used Soil and Water Assessment Tool (SWAT). This module makes it possible to assess effects of riparian wetlands on runoff and sediment yields at a watershed scale, which is very important for aquatic ecosystem management but rarely documented in the literature. In addition to delineating boundaries of a watershed and its subwatersheds, the module groups riparian wetlands within a subwatershed into an equivalent wetland for modelling purposes. Further, the module has functions to compute upland drainage area and other parameters (e.g. maximum volume) for the equivalent wetland based on digital elevation model, stream network, land use, soil and wetland distribution GIS datasets. SWAT is used to estimate and route runoff and sediment generated from upland drainage area. The lateral exchange processes between riparian wetlands and their hydraulically connected streams are simulated by the extension module. The developed module is empirically applied to the 53 km2 Upper Canagagigue Creek watershed located in Southern Ontario of Canada. The simulation results indicate that the module can make SWAT more reasonably predict flow and sediment loads at the outlet of the watershed and better represent the hydrologic processes within it. The simulation is sensitive to errors of wetland parameters and channel geometry. The approach of embedding the module into SWAT enables simulation of hydrologic processes in riparian wetlands, evaluation of wetland effects on regulating stream flow and sediment loading and assessment of various wetland restoration scenarios. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

13.
Efficiency of hydrological models mostly depends on the quality of the calibration performed prior to use. In this paper, an automatic calibration framework for the distributed hydrological model HYDROTEL is proposed. The calibration procedure was performed for three watersheds characterized with different hydroclimatological conditions: the Sassandra located in Ivory Coast, Africa, and the Montmorency and Beaurivage watersheds located in Quebec (Canada). Results of one‐a‐time (OAT) sensitivity analysis showed that the order of the most sensitive parameters differs for each watershed. Thus, the sensitivity depends on the hydroclimatic and physiographic characteristics of the watersheds. Co‐linearity indices showed that all model parameters were identifiable, that is, none of the studied parameters could be explained by a combination of the other parameters. Following these findings, an automatic calibration was run. Results indicated there was good agreement between simulated and measured streamflows at the outlet of each watershed; Nash–Sutcliffe efficiency (NSE) ranging between 0.77 and 0.92 and R2 ranging from 0.87 to 0.97. When comparing NSE and R2 values obtained using a process‐oriented, multiple‐objective, manual calibration strategy, a slight increase in model efficiency was reached with the automatic calibration procedure (4.15% for NSE and 2.95% for R2) improving predictions of peak flows for the Montmorency and Beaurivage watersheds (temperate climate conditions) and flows beyond the rainfall season in the Sassandra watershed. The proposed automatic calibration procedure introduced in this paper may be applied to other distributed hydrological model. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
This research was conducted to develop relationships among evapotranspiration (ET), percolation (PERC), groundwater discharge to the stream (GWQ), and water table fluctuations through a modeling approach. The Soil and Water Assessment Tool (SWAT) hydrologic and crop models were applied in the Big Sunflower River watershed (BSRW; 7660 km2) within the Yazoo River Basin of the Lower Mississippi River alluvial plain. Results of this study showed good to very good model performances with the coefficient of determination (R2) and Nash‐Sutcliffe efficiency (NSE) index from 0.4 to 0.9, respectively, during both hydrologic and crop model calibration and validation. An empirical relationship between ET, PERC, GWQ, and water table fluctuations was able to predict 64% of the water table variation of the alluvial plain in this study. Thematic maps were developed to identify areas with overuse of groundwater, which can help watershed managers to develop water resource programs.  相似文献   

15.
The curve number (CN) method is widely used for rainfall–runoff modelling in continuous hydrologic simulation models. A sound continuous soil moisture accounting procedure is necessary for models using the CN method. For shallow soils and soils with low storage, the existing methods have limitations in their ability to reproduce the observed runoff. Therefore, a simple one‐parameter model based on the Soil Conservation Society CN procedure is developed for use in continuous hydrologic simulation. The sensitivity of the model parameter to runoff predictions was also analysed. In addition, the behaviour of the procedure developed and the existing continuous soil moisture accounting procedure used in hydrologic models, in combination with Penman–Monteith and Hargreaves evapotranspiration (ET) methods was also analysed. The new CN methodology, its behaviour and the sensitivity of the depletion coefficient (model parameter) were tested in four United States Geological Survey defined eight‐digit watersheds in different water resources regions of the USA using the SWAT model. In addition to easy parameterization for calibration, the one‐parameter model developed performed adequately in predicting runoff. When tested for shallow soils, the parameter is found to be very sensitive to surface runoff and subsurface flow and less sensitive to ET. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

16.
With the availability of spatially distributed data, distributed hydrologic models are increasingly used for simulation of spatially varied hydrologic processes to understand and manage natural and human activities that affect watershed systems. Multi‐objective optimization methods have been applied to calibrate distributed hydrologic models using observed data from multiple sites. As the time consumed by running these complex models is increasing substantially, selecting efficient and effective multi‐objective optimization algorithms is becoming a nontrivial issue. In this study, we evaluated a multi‐algorithm, genetically adaptive multi‐objective method (AMALGAM) for multi‐site calibration of a distributed hydrologic model—Soil and Water Assessment Tool (SWAT), and compared its performance with two widely used evolutionary multi‐objective optimization (EMO) algorithms (i.e. Strength Pareto Evolutionary Algorithm 2 (SPEA2) and Non‐dominated Sorted Genetic Algorithm II (NSGA‐II)). In order to provide insights into each method's overall performance, these three methods were tested in four watersheds with various characteristics. The test results indicate that the AMALGAM can consistently provide competitive or superior results compared with the other two methods. The multi‐method search framework of AMALGAM, which can flexibly and adaptively utilize multiple optimization algorithms, makes it a promising tool for multi‐site calibration of the distributed SWAT. For practical use of AMALGAM, it is suggested to implement this method in multiple trials with relatively small number of model runs rather than run it once with long iterations. In addition, incorporating different multi‐objective optimization algorithms and multi‐mode search operators into AMALGAM deserves further research. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

17.
With the popularity of complex hydrologic models, the time taken to run these models is increasing substantially. Comparing and evaluating the efficacy of different optimization algorithms for calibrating computationally intensive hydrologic models is becoming a nontrivial issue. In this study, five global optimization algorithms (genetic algorithms, shuffled complex evolution, particle swarm optimization, differential evolution, and artificial immune system) were tested for automatic parameter calibration of a complex hydrologic model, Soil and Water Assessment Tool (SWAT), in four watersheds. The results show that genetic algorithms (GA) outperform the other four algorithms given model evaluation numbers larger than 2000, while particle swarm optimization (PSO) can obtain better parameter solutions than other algorithms given fewer number of model runs (less than 2000). Given limited computational time, the PSO algorithm is preferred, while GA should be chosen given plenty of computational resources. When applying GA and PSO for parameter optimization of SWAT, small population size should be chosen. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
This study was conducted under the USDA‐Conservation Effects Assessment Project (CEAP) in the Cheney Lake watershed in south‐central Kansas. The Cheney Lake watershed has been identified as ‘impaired waters’ under Section 303(d) of the Federal Clean Water Act for sediments and total phosphorus. The USDA‐CEAP seeks to quantify environmental benefits of conservation programmes on water quality by monitoring and modelling. Two of the most widely used USDA watershed‐scale models are Annualized AGricultural Non‐Point Source (AnnAGNPS) and Soil and Water Assessment Tool (SWAT). The objectives of this study were to compare hydrology, sediment, and total phosphorus simulation results from AnnAGNPS and SWAT in separate calibration and validation watersheds. Models were calibrated in Red Rock Creek watershed and validated in Goose Creek watershed, both sub‐watersheds of the Cheney Lake watershed. Forty‐five months (January 1997 to September 2000) of monthly measured flow and water quality data were used to evaluate the two models. Both models generally provided from fair to very good correlation and model efficiency for simulating surface runoff and sediment yield during calibration and validation (correlation coefficient; R2, from 0·50 to 0·89, Nash Sutcliffe efficiency index, E, from 0·47 to 0·73, root mean square error, RMSE, from 0·25 to 0·45 m3 s?1 for flow, from 158 to 312 Mg for sediment yield). Total phosphorus predictions from calibration and validation of SWAT indicated good correlation and model efficiency (R2 from 0·60 to 0·70, E from 0·63 to 0·68) while total phosphorus predictions from validation of AnnAGNPS were from unsatisfactory to very good (R2 from 0·60 to 0·77, E from ? 2·38 to 0·32). The root mean square error–observations standard deviation ratio (RSR) was estimated as excellent (from 0·08 to 0·25) for the all model simulated parameters during the calibration and validation study. The percentage bias (PBIAS) of the model simulated parameters varied from unsatisfactory to excellent (from 128 to 3). This study determined SWAT to be the most appropriate model for this watershed based on calibration and validation results. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

19.
The present effect of watershed subdivision on simulated water balance components using the thoroughly tested Soil and Water Assessment Tool (SWAT) model has been evaluated for the Nagwan watershed in eastern India. Observed meteorological and hydrological data (daily rainfall, temperature, relative humidity and runoff) for the years 1995 to 1998 were collected and used. The watershed and sub‐watershed boundaries, slope and soil texture maps were generated using a geographical information system. A supervised classification method was used for land‐use/cover classification from satellite imagery of 1996. In order to study the effect of watershed subdivision, the watershed was spatially defined into three decomposition schemes, namely a single watershed, and 12 and 22 sub‐watersheds. The simulation using the SWAT model was done for a period of 4 years (1995 to 1998). Results of the study showed a perfect water balance for the Nagwan watershed under all of the decomposition schemes. Results also revealed that the number and size of sub‐watersheds do not appreciably affect surface runoff. Except for runoff, there was a marked variation in the individual components of the water balance under the three decomposition schemes. Though the runoff component of the water balance showed negligible variation among the three cases, variations were noticed in the other components: evapotranspiration (5 to 48%), percolation (2 to 26%) and soil water content (0·30 to 22%). Thus, based on this study, it is concluded that watershed subdivision has a significant effect on the water balance components. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
The current study focuses on understanding key factors controlling geochemical export in eight diverse coastal watersheds at seasonal and annual time scales. Geochemical, atmospheric and hydrologic data across a range of hydro‐climatic regimes and varying land uses were investigated and relationships analysed. A hyperbolic dilution model was fitted for each watershed system to evaluate discharge–concentration relationships. Nitrate concentration effects were observed in watersheds exposed to high atmospheric deposition rates as well as agricultural watersheds, whereas urban watersheds showed nitrate dilution effects. Dilution patterns were observed for calcium, magnesium and sulfate for almost all watersheds. Seasonal loads for almost all constituents were noted to be mainly driven by hydrologic seasonality, but are also dependent on inputs (atmospheric deposition and land use sources). Understanding the primary controls on hydro‐chemical interactions is critical for developing and refining predictive water quality models, especially in coastal watersheds where sensitive downstream ecosystems act as receiving waters for upstream pollutant loads. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

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

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