首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The drainage networks of catchment areas burned by wildfire were analysed at several scales. The smallest scale (1–1000 m2) representative of hillslopes, and the small scale (1000 m2 to 1 km2), representative of small catchments, were characterized by the analysis of field measurements. The large scale (1–1000 km2), representative of perennial stream networks, was derived from a 30‐m digital elevation model and analysed by computer analysis. Scaling laws used to describe large‐scale drainage networks could be extrapolated to the small scale but could not describe the smallest scale of drainage structures observed in the hillslope region. The hillslope drainage network appears to have a second‐order effect that reduces the number of order 1 and order 2 streams predicted by the large‐scale channel structure. This network comprises two spatial patterns of rills with width‐to‐depth ratios typically less than 10. One pattern is parallel rills draining nearly planar hillslope surfaces, and the other pattern is three to six converging rills draining the critical source area uphill from an order 1 channel head. The magnitude of this critical area depends on infiltration, hillslope roughness and critical shear stress for erosion of sediment, all of which can be substantially altered by wildfire. Order 1 and 2 streams were found to constitute the interface region, which is altered by a disturbance, like wildfire, from subtle unchannelized drainages in unburned catchments to incised drainages. These drainages are characterized by gullies also with width‐to‐depth ratios typically less than 10 in burned catchments. The regions (hillslope, interface and channel) had different drainage network structures to collect and transfer water and sediment. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

2.
High and moderate severity wildfires should increase sediment production from unpaved roads due to the increased surface runoff from upslope, and increase road–stream connectivity due to the decrease in downslope surface roughness as well as the increase in surface runoff and erosion. Because no study has documented these effects, we surveyed road surface erosion features and quantified road–stream connectivity as a function of fire severity and road segment characteristics. The data were collected one year after the High Park wildfire from 141 hydrologically distinct road segments along 6.8 km of an unpaved road west of Fort Collins, Colorado. Road segments below areas burned at high and moderate severity had significantly more rills than road segments below areas that burned at low severity. Road segment slope was an important control on the proportion of segment length with rills, and the strength of the relationship between road segment slope and the amount of rilling increased with burn severity. Flatter road segments tended to capture the sediment eroded from upslope burned areas. In areas burned at high and moderate severity all of the road segments had drainage features extending to a stream, and 78% of the segments in areas burned at low severity also were connected. These exceptionally high rates of road–stream connectivity are attributed to the increased runoff from upslope, the segment‐scale collection and funneling of hillslope and road surface runoff to a single drainage point, and the reduced infiltration and trapping capacity of the burned area below the road. The results show the need to either outslope the roads or increase the frequency of constructed drainage features after wildfires, particularly for steeper road segments in areas burned at high or moderate severity. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

3.
Different hydraulic gradients, especially due to seepage or drainage, at different locations on a hillslope profile may have a profound effect on the dominant erosion processes. A laboratory study was designed to simulate hillslope processes and quantify effects of surface hydraulic gradients on erosion for a Glynwood clay loam soil (fine, illitic, mesic Aquic Hapludalf). A 5 m long, 1·2 m wide soil pan was used at 5 and 10 per cent slopes with an external watering tube to vary the soil bed's hydrological conditions. Different combinations of slope steepness with seepage or drainage gradients were used to simulate the hydrologic conditions on a 5 m segment of a hillslope profile. Runoff samples were taken during rainfall-only and rainfall with added inflow. Results showed that, under drainage conditions, interrill processes dominated and rilling was limited. The surface contained scattered crescent-shaped pits after the run. Under seepage conditions, rilling processes dominated and the inflow introduced at the top of the soil pan further accelerated the headward erosion of the rills. Erosion rates increased by as much as 60 times under seepage conditions representative of the lower backslope when compared to drainage conditions that generally occur at the upper backslope. This indicated that rills and gullies on backslopes and footslopes may be catalysed or enhanced by seepage conditions rather than form from flow hydraulic shear stress alone. An understanding of spatial and temporal changes that affect both hillslope hydrology and erosional processes is needed to develop accurate process-based erosion prediction models. This knowledge may lead to different management practices on landscape positions where seepage occurs. © 1998 John Wiley & Sons, Ltd.  相似文献   

4.
This article reports the results of a field investigation aimed to characterize the morphology of both rills monitored at Sparacia experimental area and two ephemeral gullies (EGs) located in the Tremamargi basin, Sicily, Italy. At first, the available literature data together with the measurements carried out in this investigation were used to show that the EG length is a key parameter for the estimation of the eroded volume. Then, the comparison among the pairs length and volume corresponding to measured rills, EGs and gullies showed that the exponent of the power relationship is independent of the channelized erosion type (rill, EG and gully), while a different scale factor has to be used for each erosion process. Finally, a single relationship applicable to all channelized erosion processes was deduced applying the dimensional analysis and the self‐similarity theory. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

5.
The effects of basin hydrology on hydraulic geometry of channels variability for incised streams were investigated using available field data sets and models of watershed hydrology and channel hydraulics for the Yazoo River basin,USA.The study presents the hydraulic geometry relations of bankfull discharge,channel width,mean depth,cross-sectional area,longitudinal slope,unit stream power,and mean velocity at bankfull discharge as a function of drainage area using simple linear regression.The hydraulic geometry relations were developed for 61 streams,20 of them are classified as channel evolution model(CEM) Types Ⅳ and Ⅴ and 41 of them are CEM streams Types Ⅱ and Ⅲ.These relationships are invaluable to hydraulic and water resources engineers,hydrologists,and geomorphologists involved in stream restoration and protection.These relations can be used to assist in field identification of bankfull stage and stream dimension in un-gauged watersheds as well as estimation of the comparative stability of a stream channel.A set of hydraulic geometry relations are presented in this study,these empirical relations describe physical correlations for stable and incised channels.Cross-sectional area,which combines the effects of channel width and mean channel depth,was found to be highly responsive to changes in drainage area and bankfull discharge.Analyses of cross-sectional area,channel width,mean channel depth,and mean velocity in conjunction with changes in drainage area and bankfull discharge indicated that the channel width is much more responsive to changes in both drainage area and bankfull discharge than are mean channel depth or mean velocity.  相似文献   

6.
Wildfire increases the potential connectivity of runoff and sediment throughout watersheds due to greater bare soil, runoff and erosion as compared to pre-fire conditions. This research examines the connectivity of post-fire runoff and sediment from hillslopes (< 1.5 ha; n = 31) and catchments (< 1000 ha; n = 10) within two watersheds (< 1500 ha) burned by the 2012 High Park Fire in northcentral Colorado, USA. Our objectives were to: (1) identify sources and quantify magnitudes of post-fire runoff and erosion at nested hillslopes and watersheds for two rain storms with varied duration, intensity and antecedent precipitation; and (2) assess the factors affecting the magnitude and connectivity of runoff and sediment across spatial scales for these two rain storms. The two summer storms that are the focus of this research occurred during the third summer after burning. The first storm had low intensity rainfall over 11 hours (return interval <1–2 years), whereas the second event had high intensity rainfall over 1 hour (return interval <1–10 years). The lower intensity storm was preceded by high antecedent rainfall and led to low hillslope sediment yields and channel incision at most locations, whereas the high intensity storm led to infiltration-excess overland flow, high sediment yields, in-stream sediment deposition and channel substrate fining. For both storms, hillslope-to-stream sediment delivery ratios and area-normalised cross-sectional channel change increased with the percent of catchment that burned at high severity. For the high intensity storm, hillslope-to-stream sediment delivery ratios decreased with unconfined channel length (%). The findings quantify post-fire connectivity and sediment delivery from hillslopes and streams, and highlight how different types of storms can cause varying magnitues and spatial patterns of sediment transport and deposition from hillslopes through stream channel networks.  相似文献   

7.
Loess gullies are the most active and changeable landform unit on the Loess Plateau of China. Under the influence of inhomogeneous internal and external forces, various gully morphologies have been identified as specific forms of asymmetrical loess gullies in the northeastern Loess Plateau. Thus, the formation mechanisms of asymmetrical gullies should be examined to better understand the gully evolution processes in this area. In this study, a typical asymmetrical gully area and its geological background in the northeastern Loess Plateau are investigated. Then, the asymmetrical gullies are extracted and ordered under different watershed hierarchies using 5 m horizontal resolution digital elevation models. The asymmetrical gullies are characterized using the gully deviation index and gully asymmetrical coefficient to quantitatively and qualitatively describe the gully formation from the perspective of gully morphology. Subsequently, environmental factors, such as the bedrock, climate, vegetation and interactions with neighbouring watersheds, are combined to achieve an in-depth understanding of the mechanisms of asymmetrical gully formation. The results show that most watersheds shift to the right side of the watershed geometric centre line, thereby forming a specific asymmetrical gully morphology. The phenomenon in which the asymmetrical degree characteristics decrease with the increase in drainage area suggests evident morphological differences on both sides of the main channel on a small scale, and relatively weak morphological differences on both sides of the main channel on a large scale. The degree of loess gully asymmetry appears higher in the area where only the windward slope is covered by loess than in areas where all slopes are covered by loess. The interaction between adjacent watersheds also influences the formation of asymmetrical gullies. These results support the understanding of asymmetrical gully formation in relation to the underlying bedrock structure and gully reorganization, thereby contributing to the development of process-based gully evolution models.  相似文献   

8.
Wildfires are landscape scale disturbances that can significantly affect hydrologic processes such as runoff generation and sediment and nutrient transport to streams. In Fall 2016, multiple large drought-related wildfires burned forests across the southern Appalachian Mountains. Immediately after the fires, we identified and instrumented eight 28.4–344 ha watersheds (four burned and four unburned) to measure vegetation, soil, water quantity, and water quality responses over the following two years. Within burned watersheds, plots varied in burn severity with up to 100% tree mortality and soil O-horizon loss. Watershed scale high burn severity extent ranged from 5% to 65% of total watershed area. Water quantity and quality responses among burned watersheds were closely related to the high burn severity extent. Total water yield (Q) was up to 39% greater in burned watersheds than unburned reference watersheds. Total suspended solids (TSS) concentration during storm events were up to 168 times greater in samples collected from the most severely burned watershed than from a corresponding unburned reference watershed, suggesting that there was elevated risk of localized erosion and sedimentation of streams. NO3-N concentration, export, and concentration dependence on streamflow were greater in burned watersheds and increased with increasing high burn severity extent. Mean NO3-N concentration in the most severely burned watershed increased from 0.087 mg L−1 in the first year to 0.363 mg L−1 (+317%) in the second year. These results suggest that the 2016 wildfires degraded forest condition, increased Q, and had negative effects on water quality particularly during storm events.  相似文献   

9.
This paper describes an analysis of natural and anthropogenic factors controlling the evolution of gullies in a rural basin in the basaltic upland in the State of Rio Grande do Sul, Southern Brazil. In this region of deep ferrallitic soils with more than 60% clay, runoff and erosion are of increasing concern. In the TaboAo drainage basin (100 km^2), gully erosion was studied in a field survey that measured rills and gullies. Eighty-four gullies were identified. They had an average length of 136 m, were 10 m wide, and 3 m deep and had a volume of 15.458 m3. Each gully was characterised in terms of factors that included slope, geological structure, presence of piping, drainage, soil use, and the presence of surface and subsurface flow. On average, the main channels had knickpoints varying from 2 m to 7 m, and their evolution in the vertical plane increased until bed-rock basalt material was reached, after which gullies increase in width and length. Gully development was also monitored from 1991 to 2003. Subsurface flow appears to be the principal agent controlling their development. Results show that both natural (slope, surface curvature, geological structure and rainfall) and anthropogenic (soil use, road construction) factors are important in gully development. The change in cultural practices throughout the drainage basin from conventional to direct seeding has led to increased subsurface flow, which was more important than surface runoff in causing erosion. However, the higher rainfall during E1 Nifio Southern Oscillation (ENSO) events and the consequently higher subsurface flow were the dominant factors. From 1991 to 2003 a total land loss of 1,013 m3 was observed in one gully, with 236 m^3 lost during the 1992 ENSO and 702 m3 during the 1997 ENSO; 95% of the total volume lost occurred during ENSO periods.  相似文献   

10.
A study was carried out on a rural catchment located in northwest Spain to examine the sediment yield from the catchment by measuring suspended sediments during rainfall events. Within the catchment regular surveys were conducted to obtain data on the suspended sediment sources. Important variations in sediment load were detected at event scale (0·3–21·0 Mg); some of these can be explained in terms of event size, antecedent conditions, rainfall distribution and soil surface erosion. To study the variables controlling suspended sediment yield during the events in the catchment, several event and pre‐event variables were calculated for all events. The sediment load is strongly influenced by discharge variables. During the events discharge–suspended sediments were also analysed. When the soil surface was unprotected, the formation of rills and ephemeral gullies on agricultural land at the catchment head was an important source of suspended sediments in the catchment. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

11.
Factors affecting rill erosion of unpaved loess roads in China   总被引:1,自引:0,他引:1       下载免费PDF全文
This study explores factors that affect road surface erosion in a small watershed on the Loess Plateau. Global positioning system (GPS)‐assisted field surveys and geographical information system methods were applied. The results show that road surface rills in the watershed are more easily formed on main roads, which are disturbed by intensive human activities. Secondary unpaved road networks occupied the largest road surface area and contributed 49% of the total road surface rill volumes. Spatial analysis reveals that roads near residential areas or leading to other human‐disturbed land‐use types are at high risk of soil loss. In each road segment, slope gradient, road segment length and drainage area have impacts on surface rill formation and development. Among these factors, slope gradients have been verified as a controlling factor of rill erosion intensification. Both road segment length (R = 0.83, N = 82) and drainage area (R = 0.72 for road segment and 0.76 for upslope drainage areas, N = 82) significantly influence total road surface rill volumes. The interaction variable of road segment length multiplied by slope is more closely correlated with road segment soil loss than that of the independent variables alone. Linear equations composed of slope gradient, road segment length and upslope drainage area are proposed. The new equation performs much better at predicting surface soil loss from secondary road segments compared with the previous models, which have not considered upslope drainage areas. The relationships and equations from this study will be helpful for road erosion evaluation in a small watershed of the study area. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

12.
A computer model has been used to estimate soil loss and sediment yield from irregular field-size units of small watersheds. Input to the model includes spring data (i.e. relating to February through May) for the independent variables of the Universal Soil Loss Equation, and for factors such as surface roughness, an index of overland runoff, and proximity to the stream. Output from the model includes maps of seasonal estimates of potential soil losses, field sediment delivery ratios, and expected sediment yields. On the basis of selected erosion and sediment yield tolerances, the output information has been analysed to identify watershed areas which (1) exhibit both erosion and sediment yield problems; (2) exhibit only erosion problems; (3) exhibit only sediment yield problems; and (4) exhibit neither erosion nor sediment yield problems. The percentage of the watershed area in each category and the percentage of the watershed soil loss and sediment loads contributed by each category are also identified. Application of the procedure for planning remedial control programs for five watersheds is discussed.  相似文献   

13.
流域水系自动提取的方法和应用   总被引:31,自引:3,他引:28  
李昌峰  冯学智  赵锐 《湖泊科学》2003,15(3):205-212
讨论由栅格数字高程模型(DEM)自动提取流域水系的原理和方法,并以西苕溪中上游流域为研究区进行了河网生成实验. 研究表明:在山地丘陵区和平均地形坡度不小于3°的区域,所生成河网具有很高的可靠性. 为了解决在平均地形坡度小于3°的平坦区域河网生成中产生的虚拟河网与自然水系偏差较大的问题,提出了利用主干河道和平原水系数字化作为约束条件的生成河网的方法,取得了与实际情况比较接近的结果,从而使水文要素的模拟更具有实际意义.  相似文献   

14.
1 INTRODUCTION Erosion caused by ephemeral flows is a frequent phenomenon in nature and contributes to the shape of the landscape. This type of erosion may cause great soil losses in agricultural areas, which are quickly transferred to the watershed outlets through the rill and gully network (Bennett et al., 2000; Poesen et al., 2003). Concentrated flow erosion is controlled by the erodibility of surface materials, climate, soil use and management, and watershed topography. Several metho…  相似文献   

15.
An extensive survey and topographic analysis of five watersheds draining the Luquillo Mountains in north‐eastern Puerto Rico was conducted to decouple the relative influences of lithologic and hydraulic forces in shaping the morphology of tropical montane stream channels. The Luquillo Mountains are a steep landscape composed of volcaniclastic and igneous rocks that exert a localized lithologic influence on the stream channels. However, the stream channels also experience strong hydraulic forcing due to high unit discharge in the humid rainforest environment. GIS‐based topographic analysis was used to examine channel profiles, and survey data were used to analyze downstream changes in channel geometry, grain sizes, stream power, and shear stresses. Results indicate that the longitudinal profiles are generally well graded but have concavities that reflect the influence of multiple rock types and colluvial‐alluvial transitions. Non‐fluvial processes, such as landslides, deliver coarse boulder‐sized sediment to the channels and may locally determine channel gradient and geometry. Median grain size is strongly related to drainage area and slope, and coarsens in the headwaters before fining in the downstream reaches; a pattern associated with a mid‐basin transition between colluvial and fluvial processes. Downstream hydraulic geometry relationships between discharge, width and velocity (although not depth) are well developed for all watersheds. Stream power displays a mid‐basin maximum in all basins, although the ratio of stream power to coarse grain size (indicative of hydraulic forcing) increases downstream. Excess dimensionless shear stress at bankfull flow wavers around the threshold for sediment mobility of the median grain size, and does not vary systematically with bankfull discharge; a common characteristic in self‐forming ‘threshold’ alluvial channels. The results suggest that although there is apparent bedrock and lithologic control on local reach‐scale channel morphology, strong fluvial forces acting over time have been sufficient to override boundary resistance and give rise to systematic basin‐scale patterns. Copyright © 2010 John Wiley and Sons, Ltd.  相似文献   

16.
Sediment yield is a complex function of many environmental factors including climate,hydrology,vegetation,basin topography,soil types,and land cover.We present a new semi-physical watershed sediment yield model for the estimation of suspended sediment in loess region.This model is composed by three modules in slope,gully,and stream phases.For slope sediment yield,a balance equation is established based on the concept of hydraulic erosion capacity and soil erosion resistance capacity.According to the statistical analysis of watershed characteristics,we use an exponential curve to approximately describe the spatial variability of watershed soil erosion resistance capacity.In gully phase,the relationship between gully sediment concentration and flow velocity is established based on the Bagnold'stream power function.In the stream phase,we assume a linear dependence of the sediment volume in the reach on the weighted sediment input and output.The proposed sediment yield model is operated in conjunction with a conceptual hydrologic model,and is tested over 16 regions including testing grounds,and small,medium and large watersheds in the loess plateau region in the mid-reach of Yellow River.Our results indicate that the model is reasonable in structure and is able to provide a good simulation of sediment generation and transportation processes at both flood event scale and inter-annual time scale.The proposed model is generally applicable to the watersheds with soil texture similar to that of the loess plateau region in the Yellow River basin in China.  相似文献   

17.
The formation of erosion rills and gullies is a critical step in land surface development, but possibilities to study initial unaffected surface development under natural conditions and with well‐defined initial and boundary conditions are rare. The objective of this study was to characterize rill network development from ’point zero’ in the artificially‐created catchment ‘Hühnerwasser’. To ensure unaffected development, the study was largely restricted to the analysis of remotely‐sensed data. We analyzed a series of photogrammetry‐based digital elevation models (DEMs) for 10 points in time, over a period of five years and beginning with the initial state. The evolving erosion rill network was quantitatively described based on mapping from aerial photographs. DEMs and rill network maps were combined to specifically analyze the development of morphometry for different parts of the network and to characterize energy dissipation and connectivity. The restriction to remote‐sensing data did not allow for analyzing specific processes governing rill network development, nevertheless, two major development phases could be characterized. We observed a phase of growth of the rill network along with variations in drainage patterns during the first two years of development and a subsequent phase of reduction of its area along with comparably stable patterns. Region‐specific analysis of morphometry indicates that, besides effects of changing sediment characteristics and vegetation cover development, locally evolving hydro‐geomorphic feedback cycles influenced this development. Results show an increasing similarity of overall statistical characteristics (e.g. drainage density) for two parts of the catchment, but a persistent influence of initial conditions on specific rill geometry. The observed development towards higher orderliness and increased connectivity is consistent with experiments and concepts on drainage network evolution across scales; however, we did not observe major influences of rill piracy and cross grading or a reduction of energy dissipation with network development. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

18.
Large-scale flood modelling approaches designed for regional to continental scales usually rely on relatively simple assumptions to represent the potentially highly complex river bathymetry at the watershed scale based on digital elevation models (DEMs) with a resolution in the range of 25–30 m. Here, high-resolution (1 m) LiDAR DEMs are employed to present a novel large-scale methodology using a more realistic estimation of bathymetry based on hydrogeomorphological GIS tools to extract water surface slope. The large-scale 1D/2D flood model LISFLOOD-FP is applied to validate the simulated flood levels using detailed water level data in four different watersheds in Quebec (Canada), including continuous profiles over extensive distances measured with the HydroBall technology. A GIS-automated procedure allows to obtain the average width required to run LISFLOOD-FP. The GIS-automated procedure to estimate bathymetry from LiDAR water surface data uses a hydraulic inverse problem based on discharge at the time of acquisition of LiDAR data. A tiling approach, allowing several small independent hydraulic simulations to cover an entire watershed, greatly improves processing time to simulate large watersheds with a 10-m resampled LiDAR DEM. Results show significant improvements to large-scale flood modelling at the watershed scale with standard deviation in the range of 0.30 m and an average fit of around 90%. The main advantage of the proposed approach is to avoid the need to collect expensive bathymetry data to efficiently and accurately simulate flood levels over extensive areas.  相似文献   

19.
Extreme rainfall events (>50 mm day?1) falling on cultivated land which is relatively bare of vegetation cover, typically give rise to networks of rills and gullies with ephemeral gullies in depressions and valley bottoms. Farming practices such as the removal of field boundaries, the presence of wheelings and rolled surfaces encourage runoff. The coincidence of vulnerable crops such as maize, potatoes and sugar beet with erodible soils and sloping sites may lead to high rates of erosion associated with single events or wet seasons. Not all extreme rainfall events lead to runoff and erosion, this depends on timing with respect to the growing crop. Rates of erosion associated with extreme events may be high but when placed in a long‐term temporal context, they tend to be quite low. Extreme events frequently lead to off‐site impacts most notably muddy flooding of properties and the pollution of watercourses. Landscapes may be protected from extreme events by standard soil conservation techniques; off‐site impacts may similarly be alleviated by flood‐protection measures. In both cases, the challenge is to put in place adequate economic incentives, social pressures and governmental policy frameworks to incentivise effective action. Predicted rainfall changes in the future include wetter winters and increases in rain per rain‐day. In this case, the risk of erosion on cultivated land will increase. However, erosion mitigation strategies should still address the issue of the incidence of high‐risk crops on vulnerable sites. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

20.
Headwater streams expand, contract, and disconnect in response to seasonal moisture conditions or those related to individual precipitation events. The fluctuation of the surface flow extent, or active drainage network, reflects catchment storage characteristics and has important impacts on stream ecology; however, the hydrological mechanisms that drive this phenomenon are still uncertain. Here, we present field surveys of the active drainage networks of four headwater streams in Central Idaho's Frank Church‐River of No Return Wilderness (7–21 km2) spanning the spring and summer months of 2014. We report the total length of the active drainage networks, which varied as a power law function with stream discharge with an average exponent of 0.11 ± 0.03 (range of 0.05–0.20). Generally, these active drainage networks were less responsive to changes in discharge than many streams in past studies. We observed that the locations where surface flow originates, or flowheads, were often stable, and an average of 64% of the change in active drainage network length was explained by downstream discontinuities. Analysis of geologic and geomorphic characteristics of individual watersheds and flowheads suggests that most flowheads below approximately 2200 m are supported by stable flowpaths controlled by bedrock structure. At higher elevations, small accumulation areas and saturation of shallow and conductive soil and colluvium after snowmelt result in more mobile flowhead locations. The dynamics of active drainage networks can help illuminate the spatiotemporal structure of flowpaths supporting surface flow. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

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

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