首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Spatial heterogeneity of soil has great impacts on dynamic processes of the hydrological systems. However, it is challenging and expensive to obtain spatial distribution of soil hydraulic properties, which often requires extensive soil sampling and observations and intensive laboratory analyses, especially in high elevation, hard to access mountainous areas. This study evaluates the impacts of soil heterogeneity on hydrological process in a high elevation, topographically complex watershed in Northwest China. Two approaches were used to derive the spatial heterogeneity of soil properties in the study watershed: (1) the spatial clustering method, Full‐Order‐CLK was used to determine five soil heterogeneous clusters (configurations 97, 80, 60, 40 and 20) through large number of soil sampling and in situ observations, and (2) the average values of soil hydraulic properties for each soil type were derived from the coarse provincial soil data sets (Gansu Soil Handbook at 1 : 1 000 000 scale). Subsequently, Soil and Water Assessment Tool model was used to quantify the impact of the spatial heterogeneity of soil hydraulic properties on hydrological process in the study watershed. Results show the simulations by Soil and Water Assessment Tool with the spatially clustered soil hydraulic information from the field sampling data had much better representation of the soil heterogeneity and had more accurate performance than the model using the average soil property values for each soil type derived from the coarse soil data sets. Thus, incorporating detailed field sampling, soil heterogeneity data greatly improve performance in hydrological modelling. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

2.
Many researchers have examined the impact of detailed soil spatial information on hydrological modelling due to the fact that such information serves as important input to hydrological modelling, yet is difficult and expensive to obtain. Most research has focused on the effects at single scales; however, the effects in the context of spatial aggregation across different scales are largely missing. This paper examines such effects by comparing the simulated runoffs across scales from watershed models based on two different levels of soil spatial information: the 10‐m‐resolution soil data derived from the Soil‐Land Inference Model (SoLIM) and the 1:24000 scale Soil Survey Geographic (SSURGO) database in the United States. The study was conducted at three different spatial scales: two at different watershed size levels (referred to as full watershed and sub‐basin, respectively) and one at the model minimum simulation unit level. A fully distributed hydrologic model (WetSpa) and a semi‐distributed model (SWAT) were used to assess the effects. The results show that at the minimum simulation unit level the differences in simulated runoff are large, but the differences gradually decrease as the spatial scale of the simulation units increases. For sub‐basins larger than 10 km2 in the study area, stream flows simulated by spatially detailed SoLIM soil data do not significantly vary from those by SSURGO. The effects of spatial scale are shown to correlate with aggregation effect of the watershed routing process. The unique findings of this paper provide an important and unified perspective on the different views reported in the literature concerning how spatial detail of soil data affects watershed modelling. Different views result from different scales at which those studies were conducted. In addition, the findings offer a potentially useful basis for selecting details of soil spatial information appropriate for watershed modelling at a given scale. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

3.
A spatially distributed representation of basin hydrology and transport processes in hydrologic models facilitates the identification of critical source areas and the placement of management and conservation measures. Floodplains are critical landscape features that differ from neighbouring uplands in terms of their hydrological processes and functions. Accordingly, an important step in watershed modelling is the representation of floodplain and upland areas within a watershed. The aim of this study is (1) to evaluate four floodplain–upland delineation methods that use readily available topographic data (topographic wetness index, slope position, uniform flood stage, and variable flood stage) with regard to their suitability for hydrological models and (2) to introduce an evaluation scheme for the delineated landscape units. The methods are tested in three U.S. watersheds ranging in size from 334 to 629 km2 with different climatic, hydrological, and geomorphological characteristics. Evaluation of the landscape delineation methods includes visual comparisons, error matrices (i.e. cross‐tabulations of delineated vs reference data), and geometric accuracy metrics. Reference data were obtained from the Soil Survey Geographic (SSURGO) database and Federal Emergency Management Agency (FEMA) flood maps. Results suggest that the slope position and the variable flood stage method work very well in all three watersheds. Overall percentages of floodplain and upland areas allocated correctly were obtained by comparing delineated and reference data. Values range from 83 to 93% for the slope position and from 80 to 95% for the variable flood stage method. Future studies will incorporate these two floodplain–upland delineation methods into the subwatershed‐based hydrologic model Soil and Water Assessment Tool (SWAT) to improve the representation of hydrological processes within floodplain and upland areas. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

4.
Two hydrological models with different structures and spatial capabilities are selected to simulate the runoff and actual evapotranspiration (AET) in Yingluoxia watershed, the upper reaches of Heihe River basin in northwest of China, to validate their performances in simulating hydrological processes. They are calibrated against the observed runoff at the watershed outlet (Yingluoxia station) for the period from 1990 to 1996 and validated for the period from 1997 to 2000. Results show that in terms of the simulated hydrograph against observations and the two selected objective functions, the conceptual, lumped Water And Snow balance MODeling system (WASMOD) with simple model structure could give the same, even better results than the semi‐distributed Soil and Water Assessment Tool (SWAT) with complex structure. Compared with other model applications to the watershed, simulation for monthly runoff made in this study seems better. With regard to AET, results calculated from both models are comparable as well. Both WASMOD and SWAT are proved to be suitable and satisfactory tools in simulating hydrological processes in the study area, although both of them have strengths and limitations in applications. WASMOD model may be one of the promising alternatives in hydrological modelling. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

5.
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.  相似文献   

6.
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.  相似文献   

7.
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.  相似文献   

8.
Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.  相似文献   

9.
《水文科学杂志》2013,58(6):953-970
Abstract

The 5000 km2 topographically closed Estancia basin in central New Mexico has been the focus of several palaeoclimatic studies based on changes in the level of late Pleistocene Lake Estancia. A large, unknown volume of surface runoff and groundwater from adjacent mountains contributed to the hydrological balance during highstands and lowstands. The US Department of Agriculture hydrological model SWAT (Soil and Water Assessment Tool) and the US Geological Survey groundwater flow model MODFLOW, with the LAK2 package, were used in this study to estimate runoff and water balance under present climate. A Geographic Information Systems (GIS) interface was used for SWAT, digitized data were applied for soils and vegetation, and limited streamflow data were used to obtain an approximate calibration for the model. Simulated streamflow is generally within 30% of observed values, and simulated runoff for the entire basin is about 8% of the annual inflow volume needed to support lowstands of the former Lake Estancia. Results from the combined models suggest application to other palaeoclimate investigations in semiarid lake basins.  相似文献   

10.
Simulation of watershed scale hydrologic and water quality processes is important for watershed assessments. Proper characterization of the accuracy of these simulations, particularly in cases with limited observed data, is critical. The Soil & Water Assessment Tool (SWAT) is frequently used for watershed scale simulation. The accuracy of the model was assessed by extrapolating calibration results from a well studied Coastal Plain watershed in Southwest Georgia, USA, to watersheds within the same geographic region without further calibration. SWAT was calibrated and validated on a 16.7‐km2 subwatershed within the Little River Experimental Watershed by varying six model parameters. The optimized parameter set was then applied to a watershed of similar land use and soils, a smaller watershed with different land use and soils and three larger watersheds within the same drainage system without further calibration. Simulation results with percent bias (PB) ±15% ≤ PB < ±25% and Nash–Sutcliffe efficiency (NSE) 0.50 < NSE ≤ 0.65 were considered to be satisfactory, whereas those with PB < ±10% and 0.75 < NSE ≤ 1.00 were considered very good. With these criteria, simulation results for the five non‐calibration watersheds were satisfactory to very good. Differences across watersheds were attributed to differences in soils, land use, and surficial aquifer characteristics. These results indicate that SWAT can be a useful tool for predicting streamflow for ungauged watersheds with similar physical characteristics to the calibration watershed studied here and provide an indication of the accuracy of hydrologic simulations for ungauged watersheds. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
ABSTRACT

Surface runoff generation capacity can be modified by land-use and climate changes. Annual runoff volumes have been evaluated in a small watershed of tropical forest (Brazil), using the Soil and Water Assessment Tool (SWAT) model. Firstly, the accuracy of SWAT in runoff predictions has been assessed by default input parameters and improved by automatic calibration, using 20-year observations. Then, the hydrological response under land uses (cropland, pasture and deforested soil) alternative to tropical forest and climate change scenarios has been simulated. SWAT application has showed that, if forest was replaced by crops or pasture, the watershed’s hydrological response would not significantly be affected. Conversely, a complete deforestation would slightly increase its runoff generation capacity. Under forecasted climate scenarios, the runoff generation capacity of the watershed will tend to decrease and will not be noticeably different among the representative concentration pathways. Pasture and bare soil will give the lowest and highest runoff coefficients, respectively.  相似文献   

12.
基于改进型SIMTOP参数化径流方案和新安江模型的三层土壤水量平衡计算方法,本文构建了一个输入数据和率定参数较少、同时具有地形指数尺度转换机制、较好描述二维水文过程的简单高效的大尺度水文模型TOPX,并将其与区域环境系统集成模式RIEMS紧密耦合,以增强区域气候模式对大尺度流域径流量的定量数值模拟能力.TOPX模型在酉水河流域和泾河流域的离线测试表明:该模型对小尺度流域的径流量模拟精度较高,能够较好地描述流域水文变化过程;同时,该模型在大尺度上具有较强的分布式模拟能力,能够捕捉陆面水文过程的主要特征和时空演变特点.TOPX与RIEMS的耦合模式在泾河流域进行了在线测试,借助TOPX模型中的地形指数降尺度转换和水文过程产汇流机制,耦合模式实现了利用区域气候模式模拟的气象资料来驱动水文模型进行大尺度流域日径流量的模拟.进一步分析还表明:区域气候模式RIEMS模拟的降水时空分布数据的精度是影响耦合模式对径流量模拟效果的关键因素.  相似文献   

13.
The hydrological sensitivities to long-term climate change of a watershed in Eastern Canada were analysed using a deterministic watershed runoff model developed to simulate watershed acidification. This model was modified to study atmospheric change effects in the watershed. Water balance modelling techniques, modified for assessing climate effects, were developed and tested for a watershed using atmospheric change scenarios from both state of the art general circulation models and a series of hypothetical scenarios. The model computed daily surface, inter- and groundwater flows from the watershed. The moisture, infiltration and recharge rate are also computed in the soil reservoirs. The thirty years of simulated data can be used to evaluate the effects of climatic change on soil moisture, recharge rate and surface and subsurface flow systems. The interaction between surface and subsurface water is discussed in relation to climate change. These hydrological results raise the possibility of major environmental and socioeconomic difficulties and have significant implications for future water resource planning and management. © 1997 John Wiley & Sons, Ltd.  相似文献   

14.
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.  相似文献   

15.
Watershed areal rainfall estimation, which is one of the most important and fundamental aspects in hydrological forecasting and various kinds of catchment‐scale hydrological models, is widely used in the analysis of hydrological regime change, and its precision has a direct influence on the accuracy of hydrological forecasting and hydrological simulation. In China, it is difficult to obtain the watershed areal rainfall estimate with reliable precision and avoid the phenomenon of ‘the same effect of different parameters’ because of the low density of the rain gauge network. Therefore, a watershed rainfall data recovery approach of improving the precision of watershed areal rainfall estimation is proposed here. This approach is to build new observatories, establish the time–space relations of rainfall between newly built observatories and previously built observatories in a relatively short interval and then recover the rainfall data of newly built observatories prior to their construction through simulating the relations over a longer time. As a result, watershed rainfall information could be elaborated to improve the precision of watershed areal rainfall estimate and avoid the phenomenon of ‘the same effect of different parameters’ to a certain degree in the process of hydrological simulation. The approach is used in the hydrological simulation of Hali River catchment. In combination with the Soil Water Assessment Tool model, a better result can be obtained in the hydrological simulation. Therefore, the approach can be used in other similar catchments. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

16.
Polders are one of the most common artificial hydrological entities in the plain river network regions of China. Due to enclosed dikes, manual drainage, and irrigation intake operations, polders have had a significant impact on the hydrological processes of these areas. Distributed hydrological models are effective tools to understand and reproduce the hydrological processes of a watershed. To date, however, few models are able to simulate the drainage and irrigation intake interactions of polders at a watershed scale. This study develops a modified version of the Soil and Water Assessment Tool (SWAT) model, which is designed to better represent polders (SWATpld). The SWATpld model simulates drainage and irrigation intake processes by calculating the excess‐water storage in the inner rivers and irrigation schedule for paddy rice in the polder. Both SWAT and SWATpld models were tested for the Liyang watershed. SWATpld outperformed SWAT in simulating the daily discharge and intake of the experimental polder and predicting the monthly peak flow at the outlet of the Liyang watershed, which suggests that the modified model simulates the hydrological responses of the study watershed with polder operations more realistically than the original SWAT model does. Further evaluation at various locations and in various climate conditions would increase the confidence of this model.  相似文献   

17.
Problems related to scale continue to be at the forefront of research in hydrology. Past research into issues of scale has focused mainly on digital elevation model grid size, the appropriate number and size of sub‐areas for subdividing a watershed, parameter transferability between watersheds and appropriate scales for linking hydrological and general circulation models. Much less attention has been given to the effects of scale on the representation of land cover and hydrological model response. Recent studies with respect to changes in land cover and hydrologic response have tended to focus on the issue of land cover maturity and the conversion of land through agricultural and forestry practices. The focus of this study is to examine the impact of the level of detail at which land cover is represented in modelling the hydrological response of Wolf Creek Basin in northwest Canada. A grid‐based land cover map with a spatial resolution of 30 m is coarsened or smoothed using several common grid‐based methods of aggregating categorical data, including: pixel thinning, modal smoothing and modal aggregation. A majority rule method based on polygons is also applied to the 30 m base cover. The SLURP hydrologic model is calibrated for the base cover and used as a reference for comparing simulations for the coarsened or ‘generalized’ land cover maps. Results of the simulations are compared to examine the sensitivity of hydrologic response to generalized land cover information. Comparisons of the SLURP model runs for Wolf Creek suggest that reducing the level of detail of land cover information generally has a limited effect on hydrologic response at the outlet. However, results for averages of water balance components across the basin suggest that the local variability of hydrologic response is affected in general. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

18.
19.
The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.  相似文献   

20.
Soil moisture plays an important role in hydrology. Understanding factors (such as topography, vegetation, and meteorological conditions) that influence spatio‐temporal variability in soil moisture, and how this influence is manifested, is important for understanding hydrological processes. A number of distributed (quasi‐)physical hydrological models have been developed to investigate this subject. Previous studies have shown that the spatial differences in the distribution of soil types (residual and colluvial soils) dominantly reflect spatio‐temporal fluctuations in soil moisture and runoff. We present a methodology for assessing the spatial distribution of residual and colluvial soils, which differ with respect to their physical characteristics, in a 0·88 km2 forested catchment with complex topography and a complex land‐use history. Our method is based on penetration resistance profile data; in this data set, each data point represents soil physical characteristics within an area of about 25 m2. If the spatial distribution of soils under similar meteorological, geological, historical land use, and other conditions could be characterized on the basis of similarity in topographic features, then the spatial distribution of soil could be predicted based on relationships between various topographic indices (e.g. topographic index and local slope). We tested whether our model correctly assessed the reference data. The model's results were 90·5% correct for residual soils and 87·3% correct for colluvial soils. Further studies will quantify the relationships between topographic features of land covered by residual and colluvial soils and changes in spatio‐temporal variations in the catchment (e.g. vegetation and land use) as a function of geology or meteorology. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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