首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 12 毫秒
1.
This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km2) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

2.
Reliable assessment of the spatial distribution of soil erosion is important for making land management decisions, but it has not been thoroughly evaluated in karst geo‐environments. The objective of this study was to modify a physically based, spatially distributed erosion model, the revised Morgan, Morgan and Finney (RMMF) model, to estimate the superficial (as opposed to subsurface creep) soil erosion rates and their spatial patterns in a 1022 ha karst catchment in northwest Guangxi, China. Model parameters were calculated using local data in a raster geographic information system (GIS) framework. The cumulative runoff on each grid cell, as an input to the RMMF model for erosion computations, was computed using a combined flow algorithm that allowed for flow into multiple cells with a transfer grid considering infiltration and runoff seepage to the subsurface. The predicted spatial distributions of soil erosion rates were analyzed relative to land uses and slope zones. Results showed that the simulated effective runoff and annual soil erosion rates of hillslopes agreed well with the field observations and previous quantified redistribution rates with caesium‐137 (137Cs). The estimated average effective runoff and annual erosion rate on hillslopes of the study catchment were 18 mm and 0.27 Mg ha?1 yr?1 during 2006–2007. Human disturbances played an important role in accelerating soil erosion rates with the average values ranged from 0.1 to 3.02 Mg ha?1 yr?1 for different land uses. The study indicated that the modified model was effective to predict superficial soil erosion rates in karst regions and the spatial distribution results could provide useful information for developing local soil and water conservation plans. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
A mathematical model was developed for simulating runoff generation and soil erosion on hillslopes. The model is comprised of three modules: one for overland flow, one for soil infiltration, and one for soil erosion including rill erosion and interrill erosion. Rainfall and slope characteristics affecting soil erosion on hillslopes were analysed. The model results show that the slope length and gradient, time distribution rainfall, and distribution of rills have varying influence on soil erosion. Erosion rate increases nonlinearly with increase in the slope length; a long slope length leads to more serious erosion. The effect of the slope gradient on soil erosion can be both positive and negative. Thus, there exists a critical slope gradient for soil erosion, which is about 45° for the rate of erosion at the end of the slope and about 25° for the accumulated erosion. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

4.
Eight runoff plots, located within a small catchment within the Walnut Gulch Experimental Watershed, southern Arizona, were constructed to test the argument that sediment yield (kg m?2) decreases as plot length increases. The plots ranged in length from 2 m to 27·78 m. Runoff and sediment loss from these plots were obtained for ten natural storm events. The pattern of sediment yield from these plots conforms to the case in which sediment yield first increases as plot length increases, but then subsequently decreases. Data from the present experiment indicate that maximum sediment yield would occur from a plot 7 m long. Analysis of both runoff and sediment yield from the plots indicates that the relationship of sediment yield to plot length derives both from the limited travel distance of individual entrained particles and from a decline in runoff coefficient as plot length increases. Particle‐size analysis of eroded sediment confirms the role of travel distance in controlling sediment yield. Whether in response to the finite travel distance of entrained particles or the relationship of runoff coefficient to plot length, the experiment clearly demonstrates that the erosion rates for hillslopes and catchments cannot be simply extrapolated from plot measurements, and that alternative methods for estimating large‐area erosion rates are required. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

5.
Runoff generation and soil loss from slopes have been studied for decades, but the relationships among runoff, soil loss and rill development are still not well understood. In this paper, rainfall simulation experiments were conducted in two neighbouring plots (scale: 1 m by 5 m) with four varying slopes (17.6%, 26.8%, 36.4% and 46.6%) and two rainfall intensities (90 and 120 mm h?1) using two loess soils. Data on rill development were extracted from the digital elevation models by means of photogrammetry. The effects of rainfall intensity and slope gradient on runoff, soil loss and rill development were different for the two soils. The runoff and soil loss from the Anthrosol surface were generally higher than those from the Calcaric Cambisol surface. Higher rainfall intensity produced less runoff and more sediment for almost each treatment. With increasing slope gradient, the values of cumulative runoff and soil loss peaked, except for the treatments with 90 mm h?1 rainfall on the slopes with Anthrosol. With rainfall duration, runoff discharge decreased for Anthrosol and increased for Calcaric Cambisol for almost all the treatments. For both soils, sediment concentration was very high at the onset of rainfall and decreased quickly. Almost all the sediment concentrations increased on the 17.6% and 26.8% slopes and peaked on the 36.4% and 46.6% slopes. Sediment concentrations were higher on the Anthrosol slopes than on the Calcaric Cambisol slopes. At 90 mm h?1 rainfall intensity, increasingly denser rills appeared on the Anthrosol slope as the slope gradient increased, while only steep slopes (36.4% and 46.6%) developed rills for the Calcaric Cambisol soil. The contributions of rill erosion ranged from 36% to 62% of the cumulative soil losses for Anthrosol, while the maximum contribution of rill erosion to the cumulative soil loss was only 37.9% for Calcaric Cambisol. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

6.
In most regions of the world overgrazing plays a major role in land degradation and thus creates a major threat to natural ecosystems. Several feedbacks exist between overgrazing, vegetation, soil infiltration by water and soil erosion that need to be better understood. In this study of a sub‐humid overgrazed rangeland in South Africa, the main objective was to evaluate the impact of grass cover on soil infiltration by water and soil detachment. Artificial rains of 30 and 60 mm h?1 were applied for 30 min on 1 m2 micro‐plots showing similar sandy‐loam Acrisols with different proportions of soil surface coverage by grass (Class A: 75–100%; B: 75–50%; C: 50–25%; D: 25–5%; E: 5–0% with an outcropping A horizon; F: 0% with an outcropping B horizon) to evaluate pre‐runoff rainfall (Pr), steady state water infiltration (I), sediment concentration (SC) and soil losses (SL). Whatever the class of vegetal cover and the rainfall intensity, with the exception of two plots probably affected by biological activity, I decreased regularly to a steady rate <2 mm h?1 after 15 min rain. There was no significant correlation between I and Pr with vegetal cover. The average SC computed from the two rains increased from 0·16 g L?1 (class A) to 48·5 g L?1 (class F) while SL was varied between 4 g m?2 h?1 for A and 1883 g m?2 h?1 for F. SL increased significantly with decreasing vegetal cover with an exponential increase while the removal of the A horizon increased SC and SL by a factor of 4. The results support the belief that soil vegetation cover and overgrazing plays a major role in soil infiltration by water but also suggest that the interrill erosion process is self‐increasing. Abandoned cultivated lands and animal preferred pathways are more vulnerable to erosive processes than simply overgrazed rangelands. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

7.
Changes of soil surface roughness under water erosion process   总被引:5,自引:0,他引:5       下载免费PDF全文
The objective of this study was to determine the changing characteristics of soil surface roughness under different rainfall intensities and examine the interaction between soil surface roughness and different water erosion processes. Four artificial management practices (raking cropland, artificial hoeing, artificial digging, and contour tillage) were used according to the local agriculture customs of the Loess Plateau of China to simulate different types of soil surface roughness, using an additional smooth slope for comparison purposes. A total of 20 rainfall simulation experiments were conducted in five 1 m by 2 m boxes under two rainfall intensities (0.68 and 1.50 mm min?1) on a 15° slope. During splash erosion, soil surface roughness decreased in all treatments except raking cropland and smooth baseline under rainfall intensity of 0.68 mm min?1, while increasing for all treatments except smooth baseline under rainfall intensity of 1.50 mm min?1. During sheet erosion, soil surface roughness decreased for all treatments except hoeing cropland under rainfall intensity of 0.68 mm min?1. However, soil surface roughness increased for the artificial hoeing and raking cropland under rainfall intensity of 1.50 mm min?1. Soil surface roughness has a control effect on sheet erosion for different treatments under two rainfall intensities. For rill erosion, soil surface roughness increased for raking cropland and artificial hoeing treatments, and soil surface roughness decreased for artificial digging and the contour tillage treatments under two rainfall intensities. Under rainfall intensity of 0.68 mm min?1, the critical soil surface roughness was 0.706 cm for the resistance control of runoff and sediment yield. Under rainfall intensity of 1.50 mm min?1, the critical soil surface roughness was 1.633 cm for the resistance control of runoff, while the critical soil surface roughness was 0.706 cm for the resistance control of sediment yield. These findings have important implications for clarifying the erosive nature of soil surface roughness and harnessing sloped farmland. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
In this paper, we analyse how the performance and calibration of a distributed event‐based soil erosion model at the hillslope scale is affected by different simplifications on the parameterizations used to compute the production of suspended sediment by rainfall and runoff. Six modelling scenarios of different complexity are used to evaluate the temporal variability of the sedimentograph at the outlet of a 60 m long cultivated hillslope. The six scenarios are calibrated within the generalized likelihood uncertainty estimation framework in order to account for parameter uncertainty, and their performance is evaluated against experimental data registered during five storm events. The Nash–Sutcliffe efficiency, percent bias and coverage performance ratios show that the sedimentary response of the hillslope in terms of mass flux of eroded soil can be efficiently captured by a model structure including only two soil erodibility parameters, which control the rainfall and runoff production of suspended sediment. Increasing the number of parameters makes the calibration process more complex without increasing in a noticeable manner the predictive capability of the model. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
The hydrology and contrasting erosional responses of two snowmelt events on arable farmland in Fife, Scotland, are compared. Snowmelt-generated runoff in January 1993 caused widespread soil erosion across eastern Scotland. Gullying was exemplified by three sites in Fife, where thaw of a drifted snowpack was augmented by rainfall to produce a larger erosive response than meteorological data alone would have predicted. Up to 127 m3 of soil was lost from individual gullies in fields sown to winter cereals. In February 1996 snowfall of comparable depth again covered the field area, but a more uniform snowpack, slower thaw, greater crop cover and lower rainfall during the thaw phase combined to lessen the impact of erosion. These case studies demonstrate the complexity of the erosion/runoff relationship for rain on snow events, in which erosional severity depends not just on snow depth but on snow distribution, thaw rate and the amount and timing of rainfall during the thaw phase. © 1998 John Wiley & Sons, Ltd.  相似文献   

10.
11.
On the basis of detailed rill surveys carried out on bare plots of different lengths at slopes of 12 per cent, basic rill parameters were derived. Rill width and maximum depth increased with plot length, whereas rill amount and cross‐sectional area, expressed per unit length, remained similar. On smaller plots, all rills were connected in a continuous transport system reaching the plot outlet, whilst on larger plots (10 and 20 m long) part of the rills ended with a deposition areas inside the plots. Amounts of erosion, calculated from rill volume and soil bulk density, were compared with soil loss measured at the plot outlets. On plots 10 and 20 m long, erosion estimated from volume of all rills was larger than measured soil loss. The latter was larger than erosion estimated from volume of contributing rills. To identify contributing soil loss area on these plots, two methods were applied: (i) ratio of total soil loss to maximum soil loss per unit area, and (ii) partition of plot area according to the ratio of contributing to total rill volume. Both methods resulted in similar areas of 21·8–23·5 m2 for the plot 10 m long and 31·2 m2 for the plot 20 m long. Identification of contributing areas enabled rill (5·9 kg m?2) and interrill (2·6 kg m?2) erosion rate to be calculated, the latter being very close to the value predicted from the Universal Soil Loss Equation. Although rill and interrill rates seemed to be similar on all plots, their ratio increased slightly with plot length. Application of this ratio to compute slope length factor of the Revised Universal Soil Loss Equation resulted in similar values to those predicted with the model. The achieved balance of soil loss suggested that all the sediment measured at the plot outlet originated from contributing rills and associated contributing rill areas. The results confirmed the utility of different plot lengths as a research tool for analysing the dynamic response of soil to rainfall–runoff. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

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

13.
Soil erosion on hillslopes occurs by processes of soil splash from raindrop impacts and sediment entrainment by surface water flows. This study investigates the process of soil erosion by surface water flow on a stony soil in a semiarid environment. A field experimental method was developed whereby erosion by concentrated flow could be measured in predefined flow areas without disturbing the soil surface. The method allowed for measurements in this study of flow erosion at a much wider range of slopes (2·6 to 30·1 per cent) and unit discharge rates (0·0007 to 0·007 m2 s−1) than have been previously feasible. Flow velocities were correlated to discharge and hydraulic radius, but not to slope. The lack of correlation between velocity and slope might have been due to the greater rock cover on the steeper slopes which caused the surface to be hydraulically rougher and thus counteract the expected effect of slope on flow velocity. The detachment data illustrated limitations in applying a linear hydraulic shear stress model over the entire range of the data collected. Flow detachment rates were better correlated to a power function of either shear stress (r2 = 0·51) or stream power (r2 = 0·59). Published in 1999 by John Wiley & Sons, Ltd.  相似文献   

14.
Sediment, nutrients and pollutants discharged from sheet and interrill erosion areas by rain‐impacted flows may influence water quality in streams and rivers. The depth of water on the soil surface influences the capacity of raindrop impacts to detach soil material underlying rain‐impacted flows, and a number of so‐called process‐based and mechanistic models erroneously use equations on the basis of the effect of water depth on splash erosion to account for this effect. Also, a number of these models require complex mathematical solutions to make them operate and can only predict sediment composition and discharges well if many of their parameters are calibrated specifically to the situations where they are being applied. Experiments with rain‐impacted flows, where flow depth and velocity over eroding surfaces have been controlled, have been reported in the literature and provide more appropriate equations to account for the drop size – flow depth interactions that affect detachment and transport of particles in rain‐impacted flows. There is a need to develop modeling approaches that rely on relevant data obtained under well‐controlled flow conditions where flow depths and velocities are known. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

15.
Concentrated flow erosion is the dominant form of winter erosion in northern France. This study correlates the ephemeral rill and gully volumes measured in 20 cultivated catchments (4–95 ha) for three consecutive winters with the size of the potential runoff-contributing areas. These areas were identified by characterizing soil surface state through crust development stage, importance of surface wheel tracks and roughness grade. A single and significant relationship was found between the size of runoff-contributing areas, estimated by this criterion, and the rill and gully volumes. This identified the proportion of the catchment area occupied by fields with a degraded surface structure as the main factor controlling the variability of erosion in a context of concentrated flow erosion on cultivated land. The extension of degraded areas was shown to be controlled by dynamic interactions between weather, land occupation and soil physical properties. This criterion accounts for the uneven distribution of rainfall in space and time. Morphological factors, such as talweg length and slope, are believed to determine part of the residual variability.  相似文献   

16.
Although the protective role of leaf litter cover against soil erosion is known for a long time, little research has been conducted on the processes involved. Moreover, the impact of soil meso‐ and macrofauna within the litter layer on erosion control is not clear. To investigate how leaf litter cover and diversity as well as meso‐ and macrofauna influence sediment discharge in subtropical forest ecosystems, a field experiment has been carried out in Southeast China. A full‐factorial random design with 96 micro‐scale runoff plots and 7 domestic leaf species was established and erosion was triggered by a rainfall simulator. Our results demonstrate that leaf litter cover protects soil from erosion (?82 % sediment discharge on leaf covered plots) by rainfall and this protection is removed as litter decomposes. The protective effect is influenced by the presence or absence of soil meso‐ and macrofauna. Fauna presence increases soil erosion rates significantly by 58 %, while leaf species diversity shows a non‐significant negative trend. We assume that the faunal effect arises from arthropods slackening and processing the soil surface as well as fragmenting and decomposing the protecting leaf litter covers. Even though the diversity level did not show a significant influence, single leaf species in monocultures show rather different impacts on sediment discharge and thus, erosion control. In our experiment, runoff plots with leaf litter from Machilus thunbergii showed the highest sediment discharge (68.0 g m?2) whereas plots with Cyclobalanopsis glauca showed the smallest rates (7.9 g m?2). Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

17.
Empirical prediction of soil erosion has both scientific and practical importance. This investigation tested USLE and USLE‐based procedures to predict bare plot soil loss at the Sparacia area, in Sicily. Event soil loss per unit area, Ae, did not vary appreciably with plot length, λ, because the decrease in runoff with λ was offset by an increase in sediment concentration. Slope steepness, s, had a positive effective on Ae, and this result was associated with a runoff coefficient that did not vary appreciably with s and a sediment concentration generally increasing with s. Plot steepness did not have a statistically detectable effect on the calculations of the soil erodibility factor of both the USLE, K, and the USLE‐M, KUM, models, but a soil‐independent relationship between KUM and K was not found. The erosivity index of the USLE‐MM model performed better than the erosivity index of the Central and Southern Italy model. In conclusion, the importance of an approach allowing soil loss predictions that do not necessarily increase with λ was confirmed together with the usability of already established and largely applied relationships to predict steepness effects. Soil erodibility has to be determined with reference to the specific mathematical scheme and conversion between different schemes seems to need taking into account the soil characteristics. The USLE‐MM shows promise for further developments. The evolutionary concept applied in the development of the USLE should probably be rediscovered to improve development of soil erosion prediction tools. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

18.
The enrichment of organic matter in interrill sediment is well documented; however, the respective roles of soil organic matter (SOM) and interrill erosion processes for the enrichment are unclear. In this study, organic matter content of sediment generated on two silts with almost identical textures, but different organic matter contents and aggregations, was tested. Artificial rainfall was applied to the soils in wet, dry and crusted initial conditions to determine the effects of soil moisture and rainfall and drying history on organic matter enrichment in interrill sediment. While erosional response of the soils varied significantly, organic matter enrichment of sediment was not sensitive to initial soil conditions. However, enrichment was higher on the silt with a lower organic matter content and lower interrill erodibility. The results show that enrichment of organic matter in interrill sediment is not directly related to either SOM content or soil interrill erodibility, but is dominated by interrill erosion processes. As a consequence of the complex interaction between soil, organic matter and interrill erosion processes, erodibility of organic matter should be treated as a separate variable in erosion models. Further research on aggregate breakdown, in particular the content and fate of the organic matter in the soil fragments, is required. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

19.
Soil erosion hinders the recovery and development of ecosystems in semiarid regions. Rainstorms, coupled with the absence of vegetation and improper land management, are important causes of soil erosion in such areas. Greater effort should be made to quantify the initial erosion processes and try to find better solutions for soil and water conservation. In this research, 54 rainfall simulations were performed to assess the impacts of vegetation patterns on soil erosion in a semiarid area of the Loess Plateau, China. Three rainfall intensities (15 mm h‐1, 30 mm h‐1 and 60 mm h‐1) and six vegetation patterns (arbors‐shrubs‐grass ‐A‐S‐G‐, arbors‐grass‐shrubs ‐A‐G‐S‐, shrubs‐arbors‐grass ‐S‐A‐G‐, shrubs‐grass‐arbors ‐S‐G‐A‐, grass‐shrubs‐arbors ‐G‐S‐A‐ and grass‐arbors‐shrubs ‐G‐A‐S‐) were examined at different slope positions (summits, backslopes and footslopes) in the plots (33.3%, 33.3%, 33.3%), respectively. Results showed that the response of soil erosion to rainfall intensity differed under different vegetation patterns. On average, increasing rainfall intensity by 2 to 4 times induced increases of 3.1 to 12.5 times in total runoff and 6.9 to 46.4 times in total sediment yield, respectively. Moreover, if total biomass was held constant across the slope, the patterns of A‐G‐S and A‐S‐G (planting arbor at the summit position) had the highest runoff (18.34 L m‐2 h‐1) and soil losses (197.98 g m‐2 h‐1), while S‐A‐G had the lowest runoff (5.51 L m‐2 h‐1) and soil loss (21.77 g m‐2 h‐1). As indicated by redundancy analysis (RDA) and Pearson correlation results, a greater volume of vegetation located on the back‐ and footslopes acted as effective buffers to prevent runoff generation and sediment yield. Our findings indicated that adjusting vegetation position along slopes can be a crucial tool to control water erosion and benefit ecosystem restoration on the Loess Plateau and other similar regions of the world. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

20.
Crop residues in conservation tillage systems are known to cause both a reduction in the erosive runoff power and an increase in the topsoil erosion resistance. In this study, the relative importance of both mechanisms in reducing soil loss by concentrated flow erosion is examined. Therefore, a method to calculate the effective flow shear stress responsible for soil detachment in the presence of a residue cover is applied. The determination of effective flow shear stress is based on the recalculation of the hydraulic radius for residue treatments. The method was tested in a laboratory flume by comparing soil detachment rates of identical pairs of soil samples that only differ in the presence or absence of crop residues. This shear stress partitioning approach and a soil detachment correction were then applied to a dataset of soil detachment measurements on undisturbed topsoil samples from a no‐till field plot on a loess‐derived soil, sampled during one growing season. Results indicate that only a small fraction (10% on average) of the difference in soil detachment rate between conventional and conservation tillage can be attributed to the dissipation of shear forces on the residues. The remaining decrease in soil detachment during concentrated runoff after a two‐year application of conservation tillage can be explained by the increased dry bulk density and root and crop residue content in the topsoil that reduces soil erodibility. After correcting for the presence of residues, the temporal variability in soil detachment rates (Dr) during concentrated flow for a given flow shear stress (τ) for both treatments can be predicted fairly well (R2 = 0·87) from dry soil bulk density (DBD, representing consolidation effects), soil moisture content (SMC, representing antecedent rainfall conditions), the dry mass of organic material (OM, representing root growth and residue decomposition) and saturated soil shear strength σs, sat using an equation of the form: This study is the first to show that the effect of conservation tillage on soil detachment rates is a result of soil property modifications affecting soil erodibility, rather than a result of the surface residue decreasing flow erosivity. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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