首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 718 毫秒
1.
Physically based soil erosion simulation models require input parameters of soil detachment and sediment transport owing to the action and interactions of both raindrops and overland flow. A simple interrill soil water transport model is applied to a laboratory catchment to investigate the application of raindrop detachment and transport in interrill areas explicitly. A controlled laboratory rainfall simulation study with slope length simulation by flow addition was used to assess the raindrop detachment and transport of detached soil by overland flow in interrill areas. Artificial rainfall of moderate to high intensity was used to simulate intense rain storms. However, experiments were restricted to conditions where rilling and channelling did not occur and where overland flow covered most of the surface. A simple equation with a rainfall intensity term for raindrop detachment, and a simple sediment transport equation with unit discharge and a slope term were found to be applicable to the situation where clear water is added at the upper end of a small plot to simulate increased slope length. The proposed generic relationships can be used to predict raindrop detachment and the sediment transport capacity of interrill flow and can therefore contribute to the development of physically‐based erosion models. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

2.
Road survey and field rainfall simulation experiments have shown that the erodibility of a road surface is dynamic. In the absence of extreme runoff events, dynamic erodibility results from the generation and removal of easily entrained surface material by human road surface maintenance activities, vehicular detachment and overland flow events. Maintenance activities introduce easily transportable material to the road surface where it can be entrained by overland flow. Traffic in dry conditions detaches material that is quickly removed during subsequent overland flow events. The pre‐storm erodibility of a road is therefore largely a function of maintenance and vehicle activity since the last overland flow event. During rainstorms, vehicle passes increase sediment production by detaching/redistributing surface material and creating efficient overland flow pathways for sediment transport. However, if incision of tracks by overland flow does not occur, post‐pass sediment transport quickly returns to pre‐pass rates. Field rainfall simulation data suggest that sediment transport resulting from during‐storm vehicle passes is greatly influenced by the presence of existing loose material, which again is a function of prior road usage and maintenance activities. Incorporation of vehicular effects into physically based road erosion models may be possible by parameterizing both during‐storm and inter‐storm changes in the supply of loose surface material as changes in surface erodibility. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

3.
This study examines the size characteristics of sediment removed from a semiarid hillslope by interrill overland flow. Rainfall simulation experiments were conducted on a runoff plot 18 m wide and 35 m long established on a piedmont hillslope in southern Arizona. The top of the plot coincided with the hillslope divide, and its outlet was located within a shallow rill. Samples of runoff were obtained from two cross-sections located in the interrill portion of the plot upslope of the rill and from a calibrated flume through which was directed interrill overland flow reaching the bottom of the plot. Analyses of sediment contained in these samples showed that sediment in interrill flow is finer than the matrix soil. The fineness of the interrill sediment compared to the matrix soil appears to be due to the inability of interrill overland flow to transport the coarser fraction of the sediment supplied to it by raindrop detachment. This finding implies that the rate of soil erosion in interrill areas is not. as is commonly supposed, limited by the rate at which raindrops can detach sediment but by the rate at which they detach sediment of a size that the overland flow is competent to transport. The relative fineness of sediment eroded from this hillslope is consistent with other evidence for the recent evolution of shrub-covered hillslopes in southern Arizona.  相似文献   

4.
The event‐ and physics‐based KINEROS2 runoff/erosion model for predicting overland flow generation and sediment production was applied to unpaved mountain roads. Field rainfall simulations conducted in northern Thailand provided independent data for model calibration and validation. Validation shows that KINEROS2 can be parameterized to simulate total discharge, sediment transport and sediment concentration on small‐scale road plots, for a range of slopes, during simulated rainfall events. The KINEROS2 model, however, did not accurately predict time‐dependent changes in sediment output and concentration. In particular, early flush peaks and the temporal decay in sediment output were not predicted, owing to the inability of KINEROS2 to model removal of a surface sediment layer of finite depth. After 15–20 min, sediment transport declines as the supply of loose superficial material becomes depleted. Modelled erosion response was improved by allowing road erodibility to vary during an event. Changing the model values of erosion detachment parameters in response to changes in surface sediment availability improved model accuracy of predicted sediment transport by 30–40%. A predictive relationship between road erodibility ‘states’ and road surface sediment depth is presented. This relationship allows implementation of the dynamic erodibility (DE) method to events where pre‐storm sediment depth can be estimated (e.g., from traffic usage variables). Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
Effects of rainfall patterns on runoff and rainfall-induced erosion   总被引：3，自引：0，他引：3  
Rainfall-induced erosion involves the detachment of soil particles by raindrop impact and their transport by the combined action of the shallow surface runoff and raindrop impact.Although temporal variation in rainfall intensity(pattern)during natural rainstorms is a common phenomenon,the available information is inadequate to understand its effects on runoff and rainfall-induced erosion processes.To address this issue,four simulated rainfall patterns(constant,increasing,decreasing,and increasing-decreasing)with the same total kinetic energy were designed.Two soil types(sandy and sandy loam)were subjected to simulated rainfall using 15 cm×30 cm long detachment trays under infiltration conditions.For each simulation,runoff and sediment concentration were sampled at regular intervals.No obvious difference was observed in runoff across the two soil types,but there were significant differences in soil losses among the different rainfall patterns and stages.For varying-intensity rainfall patterns,the dominant sediment transport mechanism was not only influenced by raindrop detachment but also was affected by raindrop-induced shallow flow transport.Moreover,the efficiency of equations that predict the interrill erosion rate increased when the integrated raindrop impact and surface runoff rate were applied.Although the processes of interrill erosion are complex,the findings in this study may provide useful insight for developing models that predict the effects of rainfall pattern on runoff and erosion.  相似文献   

6.
The development and application of the physically-based and spatially-distributed mathematical model CTSS8-SED is presented. The model simulates hydrologic-hydraulic processes produced by storm events and related soil erosion and sediment transport processes at basin scale in lowland areas. The model simulates (i) storm runoff,(ii) soil detachment by raindrop impact and overland flow (gross sediment yield),(iii) sediment transport by overland flow and associated erosion-deposition processes and (iv) sediment transport by stream flow and riverbed erosion-deposition processes. A quasi two-dimensional representation of water flow and sediment transport routing is made by means of interconnected cells approach. The model is applied to simulate two flooding events in the Luduea Creek basin (Santa Fe, Argentina) occurred in April 1994 and March 2007 due to extraordinary rainfalls.  相似文献   

7.
P. I. A. Kinnell 《水文研究》2005,19(14):2815-2844
Raindrop‐impact‐induced erosion is initiated when detachment of soil particles from the surface of the soil results from an expenditure of raindrop energy. Once detachment by raindrop impact has taken place, particles are transported away from the site of the impact by one or more of the following transport processes: drop splash, raindrop‐induced flow transport, or transport by flow without stimulation by drop impact. These transport processes exhibit varying efficiencies. Particles that fall back to the surface as a result of gravity produce a layer of pre‐detached particles that provides a degree of protection against the detachment of particles from the underlying soil. This, in turn, influences the erodibility of the eroding surface. Good understanding of rainfall erosion processes is necessary if the results of erosion experiments are to be properly interpreted. Current process‐based erosion prediction models do not deal with the issue of temporal variations in erodibility during a rainfall event or variabilities in erodibility associated with spatial changes in dominance of the transport processes that follow detachment by drop impact. Although more complex erosion models may deal with issues like this, their complexity and high data requirement may make them unsuitable for use as general prediction tools. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

8.
For interrill erosion, raindrop‐induced detachment and transport of sediment by rainfall‐disturbed sheet flow are the predominant processes, while detachment by sheet flow and transport by raindrop impact are negligible. In general, interrill subprocesses are inter‐actively affected by rainfall, soil and surface properties. The objective of this work was to study the relationships among interrill runoff and sediment loss and some selected para‐meters, for cultivated soils in central Greece, and also the development of a formula for predicting single storm sediment delivery. Runoff and soil loss measurement field experiments have been conducted for a 3·5‐year period, under natural storms. The soils studied were developed on Tertiary calcareous materials and Quaternary alluvial deposits and were textured from sandy loam to clay. The second group of soils showed greater susceptibility to sealing and erosion than the first group. Single storm sediment loss was mainly affected by rain and runoff erosivity, being significantly correlated with rain kinetic energy (r = 0·64***), its maximum 30‐minute intensity (r = 0·64***) and runoff amount (r = 0·56***). Runoff had the greatest correlation with rain kinetic energy (r = 0·64***). A complementary effect on soil loss was detected between rain kinetic energy and its maximum 30‐minute intensity. The same was true for rain kinetic energy and topsoil aggregate instability, on surface seal formation and thus on infiltration characteristics and overland flow rate. Empirical analysis showed that the following formula can be used for the successful prediction of sediment delivery (Di): Di = 0·638βEI30tan(θ) (R2 = 0·893***), where β is a topsoil aggregate instability index, E the rain kinetic energy, I30 the maximum 30‐minute rain intensity and θ the slope angle. It describes soil erodibility using a topsoil aggregate instability index, which can be determined easily by a simple laboratory technique, and runoff through the product of this index and rain kinetic energy. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

9.
Monthly runoff and soil loss data of three fallow experimental plots are presented, comprising a summer and following winter season. The fallow plots were only tilled once, at the end of April. Summer runoff appeared to be controlled by rainfall intensity and conforms to the Horton model of overland flow generation. Winter runoff was primarily controlled by rainfall amount and conforms to the saturation or storage control model of runoff generation. Summer runoff volume was one fourth of winter runoff volume. Summer soil loss was twice as high as winter soil loss and was caused by high intensity, high energy rainfall. Winter soil loss was due to detachment limited erosion, caused by low intensity, low energy rainfall. Mean sediment concentration of winter runoff was one seventh of that of summer runoff. Implications for runoff and erosion of climatic change, involving increased rainfall amounts or intensities in summer or winter, are given.  相似文献   

10.
IINTRODUCTIONTheinterrillerosiononafieldplotisaffectedbythekineticenergyoftherainfall,wind,topographyfactors,propertiesofsoilandthecanopy.Theinterrillerosionoccursasthefirstdropimpactsthehillslopes.Theinterrillerosionoccursinallkindsofrainfallandtheamountofthesplasherosion,whichisthemainpartofinterrillerosion,canaccountforagreatpanofthetotalerosionamountinaheavystorm(Baner1990,Glymph1957,QianandWan1986,Zhou1981).Therefore,itisveryimportanttorevealthemechanismtoestimatetheamountofinterri…  相似文献   

11.
Distributed physically based erosion and phosphorus (P) transport models, run by the overland flow model described in Taskinen and Bruen (2006. Hydrological Processes 20 : this issue), are described. In the erosion model, the additional components to the basic model were the outflow of the particles by infiltration and a new model component, i.e. deposition when rainfall stops. Two ways of calculating the shielding factor due to the flow depth were compared. The P transport model had both dissolved P (DP) and particulate P (PP) components. The processes included in the DP model were desorption from the soil surface, advection, storage in the overland flow and infiltration. The PP model accounted for advection, storage in the flow, infiltration, detachment from the soil surface by flow and rainfall and deposition both when transport capacity of suspended solids (SS) is exceeded and when rainfall ceases. When the models were developed and validated in small agricultural fields of cohesive soil types in southern Finland, comparisons were made between corresponding processes and the significance of added components were estimated in order to find out whether increased model complexity improves the model performance. The sedigraphs were found to follow the dynamics of rainfall, emphasizing the importance of the rainfall splash component. The basic model was too slow to react to changes in rainfall and flow rates, but infiltration and deposition that acts during the cessation in rainfall improved the model significantly by enabling the modelled SS to fall sharply enough. The shielding effect of flow depth from the splash detachment was found to play a significant role. Transport capacity should also be included in erosion models when they are applied to cohesive soils. In this study, the Yalin method worked well. A strong correlation was obvious between the measured SS and total P concentrations, indicating that the main form of P in runoff is PP. This emphasizes the importance of a good sediment transport model in P transport modelling. The submodel used for DP desorption from the soil surface produced plausible results without any calibration. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
Regolith surface characteristics and response were examined over a three‐year period in a badland area in a Mediterranean middle‐mountain zone near Vallcebre (Eastern Pyrenees). Preliminary work carried out in this area indicated clear seasonal patterns of regolith properties driven by frost heaving in winter and crusting and erosion in the rest of the year. Rainfall simulations were performed with a small portable nozzle simulator in order to study seasonal changes in runoff generation, erosion rates and raindrop effect on bulk density changes. The results showed large seasonal variations in runoff and erosion responses. In?ltration rates after runoff start were correlated with precipitation depth before runoff start; runoff generation was therefore related to regolith saturation only to a very limited extent. Erosion rates were more controlled by runoff rates than by the weakness of regolith against raindrop splash, and sediment grain size increased with concentration. The combined role of antecedent regolith moisture and bulk density explained most of the seasonal variability in in?ltration, bulk density changes during rainfall and erosion rates, but some seasonal differences in sediment detachability were not explained by these variables and may be attributed to changes in roughness. Overall, runoff and erosion responses were relatively stable during spring and autumn, whereas wide variations in in?ltration rates and sediment detachment occurred in winter and summer respectively. Experiments conducted in a single season would have produced poorly representative, if not erroneous, results. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

13.
Physically based models are useful frameworks for testing intervention strategies designed to reduce elevated sediment loads in agricultural catchments. Evaluating the success of these strategies depends on model accuracy, generally established by a calibration and evaluation process. In this contribution, the physically based SHETRAN model was assessed in two similar U.K. agricultural catchments. The model was calibrated on the Blackwater catchment (18 km2) and evaluated in the adjacent Kit Brook catchment (22 km2) using 4 years of 15 min discharge and suspended sediment flux data. Model sensitivity to changes in single and multiple combinations of parameters and sensitivity to changes in digital elevation model resolution were assessed. Model flow performance was reasonably accurate with a Nash–Sutcliffe efficiency coefficient of 0.78 in Blackwater and 0.60 in Kit Brook. In terms of event prediction, the mean of the absolute percentage of difference (μAbsdiff) between measured and simulated flow volume (Qv), peak discharge (Qp), sediment yield (Sy), and peak sediment flux (Sp) showed larger values in Kit Brook (48% [Qv], 66% [Qp], 298% [Sy], and 438% [Sp]) compared with the Blackwater catchment (30% [Qv], 41% [Qp], 106% [Sy], and 86% [Sp]). Results indicate that SHETRAN can produce reasonable flow prediction but performs less well in estimation of sediment flux, despite reasonably similar hydrosedimentary behaviour between catchments. The sensitivity index showed flow volume sensitive to saturated hydraulic conductivity and peak discharge to the Strickler coefficient; sediment yield was sensitive to the overland flow erodibility coefficient and peak sediment flux to raindrop/leaf soil erodibility coefficient. The multiparameter sensitivity analysis showed that different combinations of parameters produced similar model responses. Model sensitivity to grid resolution presented similar flow volumes for different digital elevation model resolutions, whereas event peak and duration (for both flow and sediment flux) were highly sensitive to changes in grid size.  相似文献   

14.
  相似文献   

15.
Changing fire regimes and prescribed‐fire use in invasive species management on rangelands require improved understanding of fire effects on runoff and erosion from steeply sloping sagebrush‐steppe. Small (0·5 m2) and large (32·5 m2) plot rainfall simulations (85 mm h–1, 1 h) and concentrated flow methodologies were employed immediately following burning and 1 and 2 years post‐fire to investigate infiltration, runoff and erosion from interrill (rainsplash, sheetwash) and rill (concentrated flow) processes on unburned and burned areas of a steeply sloped sagebrush site on coarse‐textured soils. Soil water repellency and vegetation were assessed to infer relationships in soil and vegetation factors that influence runoff and erosion. Runoff and erosion from rainfall simulations and concentrated flow experiments increased immediately following burning. Runoff returned to near pre‐burn levels and sediment yield was greatly reduced with ground cover recovery to 40 per cent 1 year post‐fire. Erosion remained above pre‐burn levels on large rainfall simulation and concentrated flow plots until ground cover reached 60 per cent two growing seasons post‐fire. The greatest impact of the fire was the threefold reduction of ground cover. Removal of vegetation and ground cover and the influence of pre‐existing strong soil‐water repellency increased the spatial continuity of overland flow, reduced runoff and sediment filtering effects of vegetation and ground cover, and facilitated increased velocity and transport capacity of overland flow. Small plot rainfall simulations suggest ground cover recovery to 40 per cent probably protected the site from low‐return‐interval storms, large plot rainfall and concentrated flow experiments indicate the site remained susceptible to elevated erosion rates during high‐intensity or long duration events until ground cover levels reached 60 per cent. The data demonstrate that the persistence of fire effects on steeply‐sloped, sandy sagebrush sites depends on the time period required for ground cover to recover to near 60 per cent and on the strength and persistence of ‘background’ or fire‐induced soil water repellency. Published in 2009 by John Wiley & Sons, Ltd.  相似文献   

16.
Water is a major limiting factor in arid and semi‐arid agriculture. In the Sahelian zone of Africa, it is not always the limited amount of annual rainfall that constrains crop production, but rather the proportion of rainfall that enters the root zone and becomes plant‐available soil moisture. Maximizing the rain‐use efficiency and therefore limiting overland flow is an important issue for farmers. The objectives of this research were to model the processes of infiltration, runoff and subsequent erosion in a Sahelian environment and to study the spatial distribution of overland flow and soil erosion. The wide variety of existing water erosion models are not developed for the Sahel and so do not include the unique Sahelian processes. The topography of the Sahelian agricultural lands in northern Burkina Faso is such that field slopes are generally low (0–5°) and overland flow mostly occurs in the form of sheet flow, which may transport large amounts of fine, nutrient‐rich particles despite its low sediment transport capacity. Furthermore, pool formation in a field limits overland flow and causes resettlement of sediment resulting in the development of a surface crust. The EUROSEM model was rewritten in the dynamic modelling code of PCRaster and extended to account for the pool formation and crust development. The modelling results were calibrated with field data from the 2001 rainy season in the Katacheri catchment in northern Burkina Faso. It is concluded that the modified version of EUROSEM for the Sahel is a fully dynamic erosion model, able to simulate infiltration, runoff routing, pool formation, sediment transport, and erosion and deposition by inter‐rill processes over the land surface in individual storms at the scale of both runoff plots and fields. A good agreement is obtained between simulated and measured amounts of runoff and sediment discharge. Incorporating crust development during the event may enhance model performance, since the process has a large influence on infiltration capacity and sediment detachment in the Sahel. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

17.
Soil moisture dynamics have a significant effect on overland flow generation. Catchment aspect is one of the major controlling factors of overland flow and soil moisture behaviour. A few experimental studies have been carried out in the uneven topography of the Himalayas. This study presents plot‐scale experiments using portable rainfall simulator at an altitude of 1,230 m above mean sea level and modelling of overland flow using observed datasets. Two plots were selected in 2 different aspects of Aglar watershed of Lesser Himalaya; the agro‐forested (AF) plot was positioned at the north aspect whereas the degraded (DE) plot was located at the south aspect of the hillslope. HS flumes and rain gauges were installed to measure the runoff at the outlet of the plot and the rainfall depth during rainfall simulation experiments. Moreover, 10 soil moisture sensors were installed at upslope and downslope locations of both the plots at 5, 15, 25, 35, and 45 cm depth from ground level to capture the soil moisture dynamics. The tests were conducted at intensities of 79.8 and 75 mm/hr in AF plot and 82.2 and 72 mm/hr in the DE plot during Test 1 and Test 2, respectively. The observed data indicate the presence of reinfiltration process only in the AF plot. The high water holding capacity and the presence of reinfiltration process results in less runoff volume in the AF plot compared with the DE plot. The Hortonian overland flow mechanism was found to be the dominant overland flow mechanism as only a few layers of top soil get saturated during all of the rainfall–runoff experiments. The runoff, rainfall, and soil moisture data were subsequently used to calibrate the parameters of HYDRUS‐2D overland flow module to simulate the runoff hydrograph and soil moisture. The components of hydrograph were evaluated in terms of peak discharge, runoff volume and time of concentration, the results were found to be within the satisfactory range. The goodness of fit of simulated hydrographs were more than 0.85 and 0.95 for AF and DE plot, respectively. The model produced satisfactory simulation results of soil moisture for all of the rainfall–runoff experiments. The HYDRUS‐2D overland flow module was found promising to simulate the runoff hydrograph and soil moisture in plot‐scale research.  相似文献   

18.
The rangeland hydrology and erosion model (RHEM) is a new process‐based model developed by the USDA Agricultural Research Service. RHEM was initially developed for functionally intact rangelands where concentrated flow erosion is minimal and most soil loss occurs by rain splash and sheet flow erosion processes. Disturbance such as fire or woody plant encroachment can amplify overland flow erosion by increasing the likelihood of concentrated flow formation. In this study, we enhanced RHEM applications on disturbed rangelands by using a new approach for the prediction and parameterization of concentrated flow erosion. The new approach was conceptualized based on observations and results of experimental studies on rangelands disturbed by fire and/or by tree encroachment. The sediment detachment rate for concentrated flow was calculated using soil erodibility and hydraulic (flow width and stream power) parameters. Concentrated flow width was calculated based on flow discharge and slope using an equation developed specifically for disturbed rangelands. Soil detachment was assumed to begin with concentrated flow initiation. A dynamic erodibility concept was applied where concentrated flow erodibility was set to decrease exponentially during a run‐off event because of declining sediment availability. Erodibility was estimated using an empirical parameterization equation as a function of vegetation cover and surface soil texture. A dynamic partial differential sediment continuity equation was used to model the total detachment rate of concentrated flow and rain splash and sheet flow. The enhanced version of the model was evaluated against rainfall simulation data for three different sites that exhibit some degree of disturbance by fire and/or by tree encroachment. The coefficient of determination (R2) and Nash–Sutcliffe efficiency were 0.78 and 0.71, respectively, which indicates the capability of the model using the new approach for predicting soil loss on disturbed rangeland. By using the new concentrated flow modelling approach, the model was enhanced to be a practical tool that utilizes readily available vegetation and soil data for quantifying erosion and assessing erosion risk following rangeland disturbance. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
Over a two-year period, rainfall, runoff and sediment output were measured in six small agricultural catchments (3–10 ha) in south Limburg (The Netherlands). These measurements were needed for validation of an erosion model for South Limburg (LISEM). In this paper, results of the measurements are presented and processes that determine surface runoff and sediment yield during winter and summer rainfall are identified. Before the start of the measurement programme, surface slaking and crust formation on the erodible loess soils were thought to be the main cause of overland flow and soil erosion in South Limburg. This was the starting point for soil conservation measures in the area. The measurement results discussed in this paper show that in some catchments much runoff occurred in winter and that soil moisture storage capacity may be just as important for runoff generation as infiltration capacity. Therefore, when modelling soil erosion and optimizing erosion control measures for South Limburg, runoff generation through Hortonian as well as through saturation overland flow must be considered.  相似文献   

20.
Experiments were undertaken to determine the feasibility of tracing sediment movement in interrill overland flow. Crushed magnetite was introduced as a source-line 10 cm wide by 8 m long on a runoff plot 18 m wide by 29 m long located in southern Arizona. Initial magnetic susceptibilities along this source line, and along three transects located 0·25, 2·95 and 5 m downslope of the source-line, were measured. Movement of the magnetite in response to three rainfall simulation experiments was monitored. During the first two experiments, overland flow discharge was sampled at miniature flumes located along two cross sections on the plot downslope of the source-line, and at a supercritical flume at the plot outlet. Magnetic susceptibilities along the source-line and transects were measured after all three experiments. Results show that the magnetite moves very early in the experiments and that it reaches one of the flumes 2 m downslope of the source-line in 3 min. Most of the tracer moves a very short distance: 29·7 per cent is deposited within 25 cm of the source-line and only 2·2 per cent is deposited 2·95 m away. The deposition rate appears to decrease exponentially away from the source-line. Very little magnetite is recorded in the flow through the miniature flumes: in general it makes up less than 1 per cent of the total sediment load. No temporal pattern in these percentages is observed. Magnetite appears to be an effective tracer of sediment movement in interrill overland flow, though its higher density than natural soil may affect its detachment and transport.  相似文献   

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

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