首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 687 毫秒
1.
2.
Interrill erosion processes on gentle slopes are affected by mechanisms of raindrop impact, overland flow and their interaction. However, limited experimental work has been conducted to understand how important each of the mechanisms are and how they interact, in particular for peat soil. Laboratory simulation experiments were conducted on peat blocks under two slopes (2.5° and 7.5°) and three treatments: Rainfall, where rainfall with an intensity of 12 mm h?1 was simulated; Inflow, where upslope overland flow at a rate of 12 mm h?1 was applied; and Rainfall + Inflow which combined both Rainfall and Inflow. Overland flow, sediment loss and overland flow velocity data were collected and splash cups were used to measure the mass of sediment detached by raindrops. Raindrop impact was found to reduce overland flow by 10 to 13%, due to increased infiltration, and reduce erosion by 47% on average for both slope gradients. Raindrop impact also reduced flow velocity (80–92%) and increased roughness (72–78%). The interaction between rainfall and flow was found to significantly reduce sediment concentrations (73–85%). Slope gradient had only a minor effect on overland flow and sediment yield. Significantly higher flow velocities and sediment yields were observed under the Rainfall + Inflow treatment compared to the Rainfall treatment. On average, upslope inflow was found to increase erosion by 36%. These results indicate that overland flow and erosion processes on peat hillslopes are affected by upslope inflow. There was no significant relationship between interrill erosion and overland flow, whereas stream power had a strong relationship with erosion. These findings help improve our understanding of the importance of interrill erosion processes on peat. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

3.
Although numerous studies have acknowledged that vegetation can reduce erosion, few process-based studies have examined how vegetation cover affect runoff hydraulics and erosion processes. We present field observations of overland flow hydraulics using rainfall simulations in a typical semiarid area in China. Field plots (5 × 2 m2) were constructed on a loess hillslope (25°), including bare soil plot as control and three plots with planted forage species as treatments—Astragalus adsurgens, Medicago sativa and Cosmos bipinnatus. Both simulated rainfall and simulated rainfall + inflow were applied. Forages reduced soil loss by 55–85% and decreased overland flow rate by 12–37%. Forages significantly increased flow hydraulic resistance expressed by Darcy–Weisbach friction factor by 188–202% and expressed by Manning's friction factor by 66–75%; and decreased overland flow velocity by 28–30%. The upslope inflow significantly increased overland flow velocity by 67% and stream power by 449%, resulting in increased sediment yield rate by 108%. Erosion rate exhibited a significant linear relationship with stream power. M. sativa exhibited the best in reducing soil loss which probably resulted from its role in reducing stream power. Forages on the downslope performed better at reducing sediment yield than upslope due to decreased rill formation and stream power. The findings contribute to an improved understanding of using vegetation to control water and soil loss and land degradation in semiarid environments.  相似文献   

4.
Bank erosion rates and processes across a range of spatial scales are poorly understood in most environments, especially in the seasonally wet tropics of northern Australia where sediment yields are among global minima. A total of 177 erosion pins was installed at 45 sites on four sand‐bed streams (Tributaries North and Central, East Tributary and Ngarradj) in the Ngarradj catchment in the Alligator Rivers Region. Bank erosion was measured for up to 3·5 years (start of 1998/99 wet season to end of 2001/02 wet season) at three spatial scales, namely a discontinuous gully (0·6 km2) that was initiated by erosion of a grass swale between 1975 and 1981, a small continuous channel (2·5 km2) on an alluvial fan that was formed by incision of a formerly discontinuous channel between 1964 and 1978, and three medium‐sized, continuous channels (8·5–43·6 km2) with riparian vegetation. The bank erosion measurements during a period of average to above‐average rainfall established that substantial bank erosion occurred during the wet season on the two smaller channels by rapid lateral migration (Tributary Central) and by erosion of gully sidewalls due to a combination of within‐gully flows and overland flow plunging over the sidewalls (Tributary North). Minor bank erosion also occurred during the dry season by faunal activity, by desiccation and loss of cohesion of the sandy bank sediments and by dry flow processes. The larger channels with riparian vegetation (East Tributary and Ngarradj) did not generate significant amounts of sediment by bank erosion. Deposition (i.e. negative pin values) was locally significant at all scales. Bank profile form and channel planform exert a strong control on erosion rates during the wet season but not during the dry season. Copyright © 2006 Commonwealth Government of Australia.  相似文献   

5.
The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

6.
Four runs of experimental landform development, with the same uplift rate, different rainfall intensity, and the same material of different permeability adjusted by the degree of compaction, showed complicated effects of rainfall and mound-forming material. In the run with more rainfall on less permeable material, low separated ridges developed in the uplifted area, because abundant overland flow promoted valley erosion and slope processes from early stages. In the run with less rainfall on less permeable material, valley incision proceeded mostly in major valleys where surface water converges. Canyons developed during early stages and later a high massive mountain emerged. The effect of rainfall difference, however, appeared completely opposite on more permeable material accompanied by lower shear strength. In the run with more rainfall on more permeable material, a massive mountain similar to that with less rainfall on less permeable material appeared, and low separated ridges appeared in the run with less rainfall on more permeable material as in the run with more rainfall on less permeable material. In the former case, similar amount of water available for Hortonian overland flow in early stages estimated from rainfall rate and permeability can explain the development of similar landforms. In the latter case, while abundant surface water with more rainfall on less permeable material made fluvial erosion active from early stages, the deficiency in surface water with less rainfall on more permeable material apparently attenuated fluvial erosion but possibly accentuated slope processes and slope failures by seepage water flow through more permeable material of low shear strength. The active erosion from early stages apparently resulted in the development of enduring similar low landforms later in the dynamic equilibrium stage. These experimental results indicate that similar landforms can emerge from different environmental and lithologic controls, and that process does not necessarily follow from form.  相似文献   

7.
ABSTRACT

Many studies have focused on soil erosion in unmanaged Japanese cypress plantations because the sparse understory vegetation and litter covering the forest ground enhance soil erosion. In this study, soil erosion, litter, and overland flow measurements were conducted over 14 months to identify the spatio-temporal variation and examine the optimal sample size. Fifteen traps (each 0.25 m wide) were installed in line along the bottom of a 15-m-wide slope. Soil erosion and overland flow had large spatial variations as compared to litter. The temporal coefficient of variation of soil erosion and overland flow was highest during dry seasons, while smaller during wet seasons. The random sampling analysis showed that the rate of decrease in spatio-temporal variation became moderate as the sample size increased beyond six. This result indicated that the optimal sample size was five, the total width of which was equivalent to about 8% of the monitored slope width.  相似文献   

8.
A new hydrological and soil erosion model has been developed and tested: LISEM, the Limburg soil erosion model. The model uses physically based equations to describe interception, infiltration and soil water transport, storage in surface depressions, splash and flow detachment, transport capacity and overland and channel flow. From the validation results it is clear that, although the model has several advantages over other models, the results of LISEM 1.0 are far from perfect. Based on the sensitivity analysis and field observations, the main reasons for these differences seems to be the spatial and temporal variability of the soil hydraulic conductivity and the initial pressure head at the basin scale. Another reason for the differences between measured and simulated results is our lack or understanding of the theory of hydrological and soil erosion processes.  相似文献   

9.
This work addresses the temporal dynamics of riparian vegetation in large braided rivers, exploring the relationship between vegetation erosion and flood magnitude. In particular, it investigates the existence of a threshold discharge, or a range of discharges, above which erosion of vegetated patches within the channel occurs. The research was conducted on a 14 km long reach of the Tagliamento River, a braided river in north‐eastern Italy. Ten sets of aerial photographs were used to investigate vegetation dynamics in the period 1954–2011. By using different geographic information system (GIS) procedures, three aspects of geomorphic‐vegetation dynamics and interactions were addressed: (i) long‐term (1954–2011) channel evolution and vegetation dynamics; (ii) the relationship between vegetation erosion/establishment and flow regime; (iii) vegetation turnover, in the period 1986–2011. Results show that vegetation turnover is remarkably rapid in the study reach with 50% of in‐channel vegetation persisting for less than 5–6 years and only 10% of vegetation persisting for more than 18–19 years. The analysis shows that significant vegetation erosion is determined by relatively frequent floods, i.e. floods with a recurrence interval of c. 1–2.5 years, although some differences exist between sub‐reaches with different densities of vegetation cover. These findings suggest that the erosion of riparian vegetation in braided rivers may not be controlled solely by very large floods, as is the case for lower energy gravel‐bed rivers. Besides flow regime, other factors seem to play a significant role for in‐channel vegetation cover over long time spans. In particular, erosion of marginal vegetation, which supplies large wood elements to the channel, increased notably over the study period and was an important factor for in‐channel vegetation trends. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

10.
The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

12.
MODELING EPHEMERAL GULLY EROSION FOR CONSERVATION PLANNING   总被引:9,自引:0,他引:9  
1INTRODUCTIONEphemeral gully erosion,which is caused by concentrated flow within cultivated farm fields,is distinct from rill erosion.Ephemeral gully erosion is also distinct from gully erosion in permanent,deep,incised channels,formed by headcuts moving upstream.Ephemeral gully erosion is often overlooked.It is not estimated with rill-interrill erosion prediction technology such as the Revised Universal Soil Loss Equation(Renard et al.,1997),and it is often not measured in field survey…  相似文献   

13.
This study uses evidence for the long-term (35 years) pattern of soil redistribution within two agricultural fields in the UK to identify the relative importance of tillage and overland flow erosion. Spatially distributed long-term total soil redistribution data for the fields (Dalicott Farm and Rufford Forest Farm) were obtained using the caesium-137 (137Cs) technique. These data were compared with predicted patterns of soil redistribution. Recent studies have demonstrated that the redistribution of soil by tillage may be described as a diffusive process. A two-component model was, therefore, developed which accounts for soil redistribution by both overland flow and diffusive processes. Comparison of the predicted patterns of overland flow erosion alone with the observed (137Cs-derived) data indicated a poor agreement (r2 = 0.17 and 0.11). In contrast, a good agreement exists between the predicted pattern of diffusive redistribution and the observed data (r2 = 0.43 and 0.41). These results give a clear indication that diffusive processes are dominant in soil redistribution within these fields. Possible diffusive processes include splash erosion, soil creep and tillage. However, the magnitude of the diffusion coefficients for the optimum predicted pattern (c. 350–400 kg m−1 a−1) demonstrates that tillage is the only process capable of explaining the very significant soil redistribution which is indicated by the 137Cs data. Consideration is given to the implications of these results for both soil erosion prediction and landscape interpretation.  相似文献   

14.
Intensive agricultural land use in the 18th to early 20th centuries on the southeastern Piedmont resulted in substantial soil erosion and gully development. Today, many historically farmed areas have been abandoned and afforested, and such landscapes are an opportunity to study channel network recovery from disturbance by gullying. Channel initiation mapping, watershed area–slope relationships, and field monitoring of flow generation processes are used to identify channel network extent and place it in hydrologic, historical and landscape evolution context. In six study areas in the North Carolina Piedmont, 100 channel heads were mapped in fully‐forested watersheds, revealing a channel initiation relationship of 380 = AS1.27, where A is contributing area (m2) and S is local slope (m/m). Flow in these channels is generated by subsurface and overland flow. The measured relative slope exponent is lower than expected based on literature values of ~2 for forested watersheds with subsurface and overland flow, suggesting that the channel network extent may reflect a former hydrological regime. However, geomorphic evidence of recovery in channel heads within fully forested watersheds is greater than those with present day pasture. Present day channel heads lie within hollows or downslope of unchanneled valleys, which may be remnants of historical gullies, and area–slope relationships provide evidence of colluvial aggradation within the valleys. Channel network extent appears to be sensitive to land use change, with recovery beginning within decades of afforestation. Channel initiation mapping and area–slope relationships are shown to be useful tools for interpreting geomorphic effects of land use change. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

15.
The dynamics of vegetation‐driven spatial heterogeneity (VDSH) and its function in structuring runoff and sediment fluxes have received increased attention from both geomorphological and ecological perspectives, particularly in arid regions with sparse vegetation cover. This paper reviews the recent findings in this area obtained from field evidence and numerical simulation experiments, and outlines their implications for soil erosion assessment. VDSH is often observed at two scales, individual plant clumps and stands of clumps. At the patch scale, the local outcomes of vegetated patches on soil erodibility and hydraulic soil properties are well established. They involve greater water storage capacity as well as increased organic carbon and nutrient inputs. These effects operate together with an enhanced capacity for the interception of water and windborne resources, and an increased biological activity that accelerates breakdown of plant litter and nutrient turnover rates. This suite of relationships, which often involve positive feedback mechanisms, creates vegetated patches that are increasingly different from nearby bare ground areas. By this way a mosaic builds up with bare ground and vegetated patches coupled together, respectively, as sources and sinks of water, sediments and nutrients. At the stand scale within‐storm temporal variability of rainfall intensity controls reinfiltration of overland flow and its decay with slope length. At moderate rainfall intensity, this factor interacts with the spatial structure of VDSH and the mechanism of overland flow generation. Reinfiltration is greater in small‐grained VDSH and topsoil saturation excess overland flow. Available information shows that VDSH structures of sources and sinks of water and sediments evolve dynamically with hillslope fluxes and tune their spatial configurations to them. Rainfall simulation experiments in large plots show that coarsening VDSH leads to significantly greater erosion rates even under heavy rainfall intensity because of the flow concentration and its velocity increase. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

16.
In Mediterranean areas the dynamics of gully development act as an important indicator of desertification. However, little is known about the influence of climate and land‐use changes, and almost no field data exist to assess the sensitivity of a landscape to gully erosion. Two important components of gully erosion studies are the prediction of where gullies begin and where they end. To address some of these issues, topographical thresholds for gully initiation and sedimentation in six different Mediterranean study areas were established. Field measurements of local soil surface slope (S) and drainage‐basin area (A) at the point of initiation of ephemeral gullies in intensively cultivated fields (five datasets) and permanent gullies in rangelands (three datasets) were carried out. A negative power relationship of the form S = aAb was fitted through all datasets, and defined as the mean topographical threshold for gullying in the respective area. Topographically controlled slopes of sedimentation at the gully bottom were also measured. Compared to theoretical relationships for channel initiation by overland flow, relatively low values for b are obtained, suggesting a dominance of overland flow and an influence of subsurface flow. The influence of landsliding at steeper slopes appeared from the flattening of the overall negative trend in the higher slope range (S > 0·30) of the integrated dataset. Comparing the threshold lines of our datasets to the average trend lines through data found in literature revealed that vegetation type and cover could better explain differences in topographical thresholds level than climatic conditions. In cultivated fields, soil structure and moisture conditions, as determined by the rainfall distribution, are critical factors influencing topographical thresholds rather than daily rainfall amounts of the gully‐initiating events. In rangelands, vegetation cover at the time of incision appears to be the most important factor differentiating between topographical thresholds, overruling the effect of average annual rainfall amounts. Soil texture and rock fragment cover contributed little to the explanation of the relative threshold levels. Differences in regression slopes (b) between the S–A relationships found in this study have been attributed to the soil characteristics in the different study areas, determining the relative importance of subsurface flow and Hortonian overland flow. Sedimentation slopes where both ephemeral and permanent gullies end were generally high because of the high rock fragment content of the transported sediment. A positive relationship was found between the rock fragment content at the apex of the sedimentation fan and the slope of the soil surface at this location. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

17.
Little attention has been paid to the role of sidewall processes in gully development. Simple estimates of linear incision/sidewall erosion ratios at two localities in New South Wales, Australia suggest that sidewall erosion is responsible for more than half the gully volume. Studies at three localities where overland flow and throughflow are limited to only one sidewall indicate that these processes are responsible for 10-30 per cent of gully volume. These observations have important implications for gully management.  相似文献   

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

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

20.
A replicated field study using rainfall simulation and overland flow application was conducted in central Oahu, Hawaii, on a clay‐dominated Oxisol with a 9% slope. Three main treatment groups were examined: a bare treatment, a group of four rolled erosion control systems (RECSs) with open weave designs, and a group of five randomly oriented fibre RECSs. A total of 1122 measurements of runoff and erosion were made to examine treatment differences and to explore temporal patterns in runoff and sediment flux. All erosion control systems significantly delayed the time required to generate plot runoff under both simulated rainfall (35 mm h?1) and the more intense trickle flow application (114 mm h?1). Once runoff was generated during the rainfall application phase, the bare treatment runoff coefficients were significantly lower than those from the two groups of RECSs, as surface seal disruption by rilling is inferred to have enhanced infiltration in the bare treatments. During the more intense phase of overland flow application, the reverse pattern was observed. Interrill contributing‐area roughness was reduced on the bare treatment, facilitating increased runoff to well‐developed rill networks. Meanwhile, the form roughness associated with the RECSs delayed interrill flow to the poorly organized rills that formed under some of the RECSs. Regardless of runoff variations between treatments, sediment output was significantly lower from all surfaces covered by RECSs. The median cumulative sediment output from the bare surfaces was 6·9 kg, compared with 1·2 kg from the open‐weave RECSs and 0·2 kg from the random‐fibre RECSs. The random‐fibre systems were particularly effective under the more stressful overland flow application phase, with 63 times less sediment eroded than the bare treatments and 12 times less than that from the open‐weave systems. Architectural design differences between the two groups of RECSs are discussed in light of their relation to erosion process dynamics and shear stress partitioning. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

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

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