首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In semi‐arid areas, high‐intensity rainfall events are often held responsible for the main part of soil erosion. Long‐term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event‐based soil erosion models are better suited for this purpose but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event‐based erosion model was applied to the medium‐sized (～50 km2) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium‐sized catchments, (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that because of large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in Ksat values. Scenario calculations show that although ～50% of soil erosion occurs as a result of high frequency, low‐intensity rainfall events, large‐magnitude, low‐frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude–frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

2.
Distributed erosion models, which simulate the physical processes of water flow and soil erosion, are effective for predicting soil erosion in forested catchments. Although subsurface flow through multiple pathways is dominant for runoff generation in forested headwater catchments, the process-based erosion model, Geo-spatial interface for Water Erosion Prediction Project(Geo WEPP), does not have an adequate subsurface component for the simulation of hillslope water flow. In the current study, t...  相似文献   

3.
Accurate runoff and soil erosion modeling is constrained by data availability, particularly for physically based models such as OpenLISEM that are data demanding, as the processes are calculated on a cell‐by‐cell basis. The first decision when using such models is to select mapping units that best reflect the spatial variability of the soil and hydraulic properties in the catchment. In environments with limited data, available maps are usually generic, with large units that may lump together the values of the soil properties, affecting the spatial patterns of the predictions and output values in the outlet. Conversely, the output results may be equally acceptable, following the principle of equifinality. To studyhow the mapping method selected affects the model outputs, four types of input maps with different degrees of complexity were created: average values allocated to general soil map units (ASG1), average values allocated to detailed map units (ASG2), values interpolated by ordinary kriging (OK) and interpolated by kriging with external drift (KED). The study area was Ribeira Seca, a 90 km2 catchment located in Santiago Island, Cape Verde (West Africa), a semi‐arid country subject to scarce but extreme rainfall during the short tropical summer monsoon. To evaluate the influence of rainfall on runoff and erosion, two storm events with different intensity and duration were considered. OK and KED inputs produced similar results, with the latter being closer to the observed hydrographs. The highest soil losses were obtained with KED (43 ton ha? 1 for the strongest event). To improve the results of soil loss predictions, higher accurate spatial information on the processes is needed; however, spatial information of input soil properties alone is not enough in complex landscapes. The results demonstrate the importance of selecting the appropriate mapping strategy to obtain reliable runoff and erosion estimates. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

4.
A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.  相似文献   

5.
This article investigates the soil moisture dynamics within two catchments (Stanley and Krui) in the Goulburn River in NSW during a 3‐year period (2005–2007) using the HYDRUS‐1D soil water model. Sensitivity analyses indicated that soil type, and leaf area index were the key parameters affecting model performance. The model was satisfactorily calibrated on the Stanley microcatchment sites with a single point rainfall record from this microcatchment for both surface 30 cm and full‐profile soil moisture measurements. Good correlations were obtained between observed and simulated soil water storage when calibrations for one site were applied to the other sites. We extended the predictions of soil moisture to a larger spatial scale using the calibrated soil and vegetation parameters to the sites in the Krui catchment where soil moisture measurement sites were up to 30 km distant from Stanley. Similarly good results show that it is possible to use a calibrated soil moisture model with measurements at a single site to extrapolate the soil moisture to other sites for a catchment with an area of up to 1000 km2 given similar soils and vegetation and local rainfall data. Site predictions were effectively improved by our simple data assimilation method using only a few sample data collected from the site. This article demonstrates the potential usefulness of continuous time, point‐scale soil moisture data (typical of that measured by permanently installed TDR probes) and simulations for predicting the soil wetness status over a catchment of significant size (up to 1000 km2). Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

6.
Concentrated erosion, a major feature of land degradation, represents a serious problem for soil and water resources management and a threat to ecosystems. Understanding the internal mechanisms (de-)coupling sediment pathways can improve the management and resilience of catchments. In this study, concentrated erosion and deposition forms were mapped accurately through field and aerial unmanned aerial vehicle (UAV) campaigns, in order to assess the evolution of connectivity pathways over a series of three contrasted and consecutive flood events occurring between October 2016 and January 2017 (return period ranging from 0.5 to 25 years) in a small Mediterranean agricultural catchment (Can Revull, Mallorca, Spain; 1.4 km2). In addition, a morphometric index of connectivity (IC) was used to identify the potential trajectories of different concentrated erosion forms and deposition areas. IC predictions were calibrated by identifying the optimal critical thresholds, i.e. those most consistent with field observations after each of the events studied. The results found that the index performed well in predicting the occurrence and the length/area of the different type of landforms, giving kappa (κ) coefficients of variation ranging between 0.21 and 0.92 and linear correlations R2 between 0.33 and 0.72. The type of landform affected the correspondence of IC predictions and field observations, with lower thresholds the greater the magnitude of their associated geomorphic processes. Rainfall magnitude proved to be a very important factor controlling the development of erosion and deposition landforms, with large differences in length/area between the contrasted events. The evolution of the observed trajectories revealed feedback dynamics between the structural and functional connectivity of the catchment, in which morphological changes determined the spatial distribution of the processes’ activity in the successive events and vice versa. © 2020 John Wiley & Sons, Ltd.  相似文献   

7.
An efficient calibration with remotely sensed (RS) data is important for accurate predictions at ungauged catchments. This study investigates the advantages of streamflow-sensitive regionalization on calibration with RS evapotranspiration (ET). Regionalization experiments are performed at 28 catchments in Australia. The catchments are classified into three groups based on annual rainfall and runoff coefficients. Streamflow, RS ET, and a multi-objective RS ET-streamflow calibration are performed using the DiffeRential Evolution Adaptive Metropolis algorithm in each catchment. Simplified Australian Water Resource Assessment-Landscape model is calibrated for a selection of five parameters. Posterior probability distributions of parameters from three calibrations performed at donor catchments in each group are inspected to find the parameter for regionalization in the individual group. In group 1 of wetter catchments, regionalization of parameter FsoilEmax (soil evaporation scaling factor) helps to simplify the calibration without any deterioration in ET, soil moisture (SM) and streamflow predictions. Regionalization of parameter Beta (coefficient describing rate of hydraulic conductivity increase with water content) in group 2 assists to improve the streamflow predictions with no decrement in ET and SM predictions. However, regionalization is not able to provide satisfactory results in group 3. Group 3 includes low-yielding catchments, with average annual rainfall below 1000 mm/year and runoff coefficient less than 0.1, where traditional streamflow calibration also fails to produce accurate results. This study concludes that streamflow-sensitive regionalization is effective for improving the efficacy of RS ET calibration in wetter catchments.  相似文献   

8.
In this paper, we present MHYDAS‐Erosion, a dynamic and distributed single‐storm water erosion model developed as a module of the existing hydrological MHYDAS model. As with many catchment erosion models, MHYDAS‐Erosion is able to simulate sediment transport, erosion and deposition by rill and interrill processes. Its originality stems from its capacity to integrate the impact of land management practices (LMP) as key elements controlling the sedimentological connectivity in agricultural catchments. To this end, the water‐sediment pathways are first determined by a specific process‐oriented procedure defined and controlled by the user, which makes the integration of LMP easier. The LMP dynamic behaviours are then integrated into the model as a time‐dependent function of hydrological variables and LMP characteristics. The first version of the model was implemented for vegetative filters and tested using water and sediment discharge measurements at three nested scales of a densely instrumented catchment (Roujan, OMERE Observatory, southern France). The results of discharge and soil loss for simulated rainfall events have been found to acceptably compare with available data. The average R2 values for water and sediment discharge are 0·82 and 0·83, respectively. The sensitivity of the model to changes in the proportion of LMP was assessed for a single rain event by considering three scenarios of the Roujan catchment management with vegetative filters: 0% (Scenario 1), 18% (Scenario 2, real case) and 100% (Scenario 3). Compared to Scenario 2 (real case), soil losses decreased for Scenario 3 by 65% on the agricultural plot scale, 62% on the sub‐catchment scale and 45% at the outlet of the catchment and increased for Scenario 1 by 0% on the plot scale, 26% on the sub‐catchment scale and 18% at the outlet of the catchment. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
In this study, we investigate the surface flow time of rise in response to rainfall and snowmelt events at different spatial scales and the main sources originating channel runoff and spring water in a steep nested headwater catchment (Rio Vauz, Italian Dolomites), characterized by a marked elevation gradient. We monitored precipitation at different elevations and measured water stage/streamflow at the outlet of two rocky subcatchments of the same size, representative of the upper part of the catchment dominated by outcropping bedrock, at the outlet of a soil‐mantled and vegetated subcatchment of similar size but different morphology, and at the outlet of the main catchment. Hydrometric data are coupled with stable isotopes and electrical conductivity sampled from different water sources during five years, and used as tracers in end‐member mixing analysis, application of two component mixing models and analysis of the slope of the dual‐isotope regression line. Results reveal that times of rise are slightly shorter for the two rocky subcatchments, particularly for snowmelt and mixed rainfall/snowmelt events, compared to the soil‐mantled catchment and the entire Rio Vauz Catchment. The highly‐variable tracer signature of the different water sources reflects the geomorphological and geological complexity of the study area. The principal end‐members for channel runoff and spring water are identified in rainfall and snowmelt, which are the dominant water sources in the rocky upper part of the study catchment, and soil water and shallow groundwater, which play a relevant role in originating baseflow and spring water in the soil‐mantled and vegetated lower part of the catchment. Particularly, snowmelt contributes up to 64 ± 8% to spring water in the concave upper parts of the catchment and up to 62 ± 11% to channel runoff in the lower part of the catchment. These results offer new experimental evidences on how Dolomitic catchments capture and store rain water and meltwater, releasing it through a complex network of surface and subsurface flow pathways, and allow for the construction of a preliminary conceptual model on water transmission in snowmelt‐dominated catchments featuring marked elevation gradients.  相似文献   

10.
1 INTRODUCTION Soil erosion in the foothills of the Hindu Kush-Himalayas (HKH) is considered to be a hot topic in land degradation research in the region (Scherr and Yadav, 1996). The land degradation research has mainly addressed the issue of topsoil los…  相似文献   

11.
A general trend of decreasing soil loss rates with increasing vegetation cover fraction is widely accepted. Field observations and experimental work, however, show that the form of the cover‐erosion function can vary considerably, in particular for low cover conditions that prevail on arid and semiarid hillslopes. In this paper the structured spatial distribution of the vegetation cover and associated soil attributes is proposed as one of the possible causes of variation in cover–erosion relationships, in particular in dryland environments where patchy vegetation covers are common. A simulation approach was used to test the hypothesis that hillslope discharge and soil loss could be affected by variation in the spatial correlation structure of coupled vegetation cover and soil patterns alone. The Limburg Soil Erosion Model (LISEM) was parameterized and verified for a small catchment with discontinuous vegetation cover at Rambla Honda, SE Spain. Using the same parameter sets LISEM was subsequently used to simulate water and sediment fluxes on 1 ha hypothetical hillslopes with simulated spatial distributions of vegetation and soil parameters. Storms of constant rainfall intensity in the range of 30–70 mm h?1 and 10–30 min duration were applied. To quantify the effect of the spatial correlation structure of the vegetation and soil patterns, predicted discharge and soil loss rates from hillslopes with spatially structured distributions of vegetation and soil parameters were compared with those from hillslopes with spatially uniform distributions. The results showed that the spatial organization of bare and vegetated surfaces alone can have a substantial impact on predicted storm discharge and erosion. In general, water and sediment yields from hillslopes with spatially structured distributions of vegetation and soil parameters were greater than from identical hillslopes with spatially uniform distributions. Within a storm the effect of spatially structured vegetation and soil patterns was observed to be highly dynamic, and to depend on rainfall intensity and slope gradient. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

12.
Accurate estimations of water retention and detention are needed to simulate surface runoff and soil erosion following a rainfall event in a catchment. Several equations to estimate the amount of surface depressional storage, the fraction of the soil surface covered by water and the amount of rainfall excess needed to start surface runoff have been developed by Onstad (1984). The random roughness and slope gradient are needed for those estimations. Surface micro-elevation data have been gathered by a photographic method. The random roughness was determined from those elevation measurements. Several factors which have an impact on the soil surface roughness were taken into account. The main sources of influence are the type of land use, the crop stage within the growing period and tillage direction. Analyses of variance indicated that the variation in the RR-index could be explained mainly by type of land use, orientation and field type. The temporal variation was relatively small. Gradient data have been determined from a digital elevation model, constructed by digitizing contours. Combining the random roughness and the steepness of slope, the amounts of surface water retention and detention could be estimated. Knowledge of water retention and detention will improve the estimations of runoff and soil erosion modelling in catchments, such as those made with the LISEM model. The agricultural systems examined in this study have similar random roughness values in summer. Different soil erosion rates for several types of land use can not therefore be explained by the random roughness.  相似文献   

13.
This paper presents an erosion model, ARMOUR, which simulates time‐varying runoff, erosion, deposition and surface armour evolution down a hillslope either as a result of a single erosion event or as the cumulative impact of many events over periods up to decades. ARMOUR simulates sediment transport for both cohesive and non‐cohesive soil and dynamically differentiates between ‘transport‐limited’ and ‘source‐limited’ processes. A variety of feasible processes for entrainment of different size classes can be modelled and evaluated against data. The generalized likelihood of uncertainty estimation (GLUE) technique was used to calibrate and validate ARMOUR using data collected during rainfall simulator experiments at two contrasting sites: (1) non‐cohesive stony sediments at Ranger Uranium Mine, Northern Territory, Australia; and (2) cohesive silty sediments at Northparkes Gold Mine, NSW, Australia. The spatial and temporal variations of model predictions within the individual runoff events showed that some entrainment processes could not model the spikes in concentration and subsequent depletion, while the hiding model of Andrews and Parker best simulated the concentration trends for both calibrated and independent runoff events. ARMOUR also successfully captured the coarsening of the surface material, though small, over the duration of the rainfall simulator trials. This was driven by the depletion of the finest size class of the soil. For a constant discharge, ARMOUR simulated higher sediment flux at the start of the storm with the sediment flux and concentration diminishing with time. For natural rainfall a power law relationship between sediment flux and discharge was observed. The calibration exercise showed that sediment concentration and discharge alone are insufficient to calibrate all aspects of the physics, in particular the armour depth. This appears to be because the armouring during the short duration events is driven by depletion of the finest classes of the sediments (diameters less then 62·5 mm), which are not normally measured. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

14.
Modelling mean annual sediment yield using a distributed approach   总被引：3，自引：0，他引：3  
In this paper a spatially distributed model for the calculation of sediment delivery to river channels is presented (SEDEM: SEdiment DElivery Model). The model consists of two components: (1) the calculation of a spatial pattern of mean annual soil erosion rates in the catchment using a RUSLE (Revised Soil Erosion Equation) approach; and (2) the routing of the eroded sediment to the river channel network taking into account the transport capacity of each spatial unit. If the amount of routed sediment exceeds the local transport capacity, sediment deposition occurs. An existing dataset on sediment yield for 24 catchments in central Belgium was used to calibrate the transport capacity parameters of the model. A validation of the model results shows that the sediment yield for small and medium sized catchments (10–5000 ha) can be predicted with an average accuracy of 41 per cent. The predicted sediment yield values with SEDEM are significantly more accurate than the predictions using a lumped regression model. Moreover a spatially distributed approach allows simulation of the effect of different land use scenarios and soil conservation techniques. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

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

16.
ABSTRACT

Poorly monitored catchments could pose a challenge in the provision of accurate flood predictions by hydrological models, especially in urbanized areas subject to heavy rainfall events. Data assimilation techniques have been widely used in hydraulic and hydrological models for model updating (typically updating model states) to provide a more reliable prediction. However, in the case of nonlinear systems, such procedures are quite complex and time-consuming, making them unsuitable for real-time forecasting. In this study, we present a data assimilation procedure, which corrects the uncertain inputs (rainfall), rather than states, of an urban catchment model by assimilating water-level data. Five rainfall correction methods are proposed and their effectiveness is explored under different scenarios for assimilating data from one or multiple sensors. The methodology is adopted in the city of São Carlos, Brazil. The results show a significant improvement in the simulation accuracy.  相似文献   

17.
This paper presents preliminary results from the application of a transfer‐function rainfall–runoff model to ephemeral streams in Mediterranean Spain. Flow simulations have been conducted for two small catchments (Carraixet and Poyo basins), located in close proximity to one another yet with significantly different geological characteristics. Analysis of flow simulations for a number of high‐flow events has revealed the dominant influence of the rainfall on the catchment response, particularly for high‐rainfall events. Particular success has been attained modelling the highest magnitude events in both catchments and for all events in the faster responding (Poyo) catchment. In order to investigate the viability of the model for forecasting floods in ungauged catchments, additional investigations have been conducted by calibrating the model for one catchment (donor catchment) and then applying it to another (receptor catchment). The results indicate that this can be successful when either the donor catchment is a fast response catchment or when the model is calibrated using a high‐magnitude event in the donor catchment, providing that the modelled receptor catchment event is of a lower magnitude. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

18.
Abstract

Soil erosion vulnerability and extreme rainfall characteristics over the Mediterranean semi-arid region of Tunisia are crucial input for estimation of siltation rate in artificial reservoirs. A comprehensive high-resolution database on erosive rainfall, together with siltation records for 28 small reservoirs, were analysed for this region, the Tunisian Dorsal (the easternmost part of the Atlas Mountains). The general life-span of these reservoirs is only about 14 years. Depending on the soil degradation in the different catchments, the corresponding reservoirs display a wide range of soil erosion rates. The average soil loss was 14.5 t ha?1 year?1 but some catchments display values of up to 36.4 t ha?1 year?1. The maximum 15-min duration rainfall intensity was used to determine the spatial distribution of rainfall erosivity. The northwestern parts of the Tunisian Dorsal display the most extreme rainfall erosivity. Spatial erosion patterns are to some extent similar; however, they vary greatly according to their location in the “soil degradation cycle”. This cycle determines the soil particle delivery potential of the catchment. In general, the northwestern parts of the Dorsal display modest soil erosion patterns due to the already severely degraded soil structure. Here, the soil surface is often the original bedrock. However, the greatest soil erosion occurs in the mid-eastern parts of the Dorsal, which represents the “degradation front”. The latter corresponds to the area with highest erosion, which is continuously progressing westward in the Dorsal. The large variation between the erosive rainfall events and the annual soil loss rates was explained by two important factors. The first relates to the soil degradation cycle. The second factor corresponds to the degradation front with the highest soil loss rates. At present this front is located at 300 m altitude and appears to be moving along an 80-km westward path starting from the east coast. A better understanding of the above can be used to better manage soils and soil covers in the Tunisian Dorsal area and, eventually, to decrease the soil erosion and reservoir siltation risk.

Citation Jebari, S., Berndtsson, R., Bahri, A. & Boufaroua, M. (2010) Spatial soil loss risk and reservoir siltation in semi-arid Tunisia. Hydrol. Sci. J. 55(1), 121–137.  相似文献   

19.
This paper explores a scale‐adapted erosion mapping method which aims at a rapid assessment of field erosion and sediment transport pathways in catchments up to several square kilometres and compares the results with the output of a well‐known erosion model (LISEM). The mapping method is based on an event‐defined classification scheme of erosion intensity (zero, weak, moderate and strong) that is applied to arable fields, in combination with incision measurements of erosion features for each erosion intensity class on a small sample of fields. Sediment deposition is classified on the basis of quantity indicators and abundance. In addition, relevant conditions and erosion factors are determined for each field. The method was applied to an agricultural catchment (4·2 km2) in the Sundgau (Alsace), after a short but violent thunderstorm in May 2001, to illustrate its potential use and its limitations. The rainfall event led to strong erosion on the arable fields and a muddy flow that caused significant damage in the built‐up area. On the basis of the analyses of the incision measurements in combination with the mapping of erosion intensity classes, total erosion for the catchment was estimated as 15 000 t (an average of about 36 t[sol ]ha). Sediment deposition was found to occur in three major locations: (1) in thalwegs at the interface between maize and downslope winter wheat fields, (2) in downslope headlands where the flow direction suddenly changed due to oriented tillage structures in the perpendicular direction, and (3) the lowest corners of fields which collect all the runoff from the field. Preliminary data analyses suggest that erosion intensity is related to field size and[sol ]or tillage direction and to slope morphology. Model output (LISEM) appeared to depend more strongly on slope gradient than the results obtained with the mapping method. The method yields a database, which can be used as a foundation for conservation strategies in small regions with similar land use and geomorphology. The mapping and modelling methods are compared, and their complementary aspects are highlighted. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
This paper focuses on the problem of quantifying real world catchment response using a distributed model and discusses the ability of the model to capture that response. The rainfall–runoff responses of seven small agricultural catchments in the eastern wheatbelt region of south-western Australia are examined. The variability in runoff generation and the factors that contribute to that variability (i.e. rainfall intensity, soil properties and topography) are investigated to determine if their influence can be captured in a mathematical model. The spatially distributed rainfall–runoff model used in this study is based on the TOPMODEL concepts of Beven and Kirkby (1979), and simulates runoff generation by both the infiltration excess and saturation excess mechanisms. Simulations with the model revealed the highly complex nature of catchment response to rainfall events. Runoff generation was highly heterogeneous in both space and time, with the runoff response being governed by the spatial variability of soil properties and topography, and by the temporal variation in rainfall intensity. Although the model proved capable of simulating catchment response for many events, the investigation has demonstrated that not all aspects of the variability associated with agricultural catchments (particularly the effects of land management) can be captured using this relatively simple model. © 1997 by John Wiley & Sons, Ltd  相似文献   

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

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