首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 38 毫秒
1.
Declining water quality on the south coast of Western Australia has been linked to current agricultural practices. Riparian buffers were identified as a tool available to farmers and catchment managers to achieve water quality improvements. This study compares 10 m wide regenerating grass and Eucalyptus globulus buffer performance. Surface and subsurface water quality were monitored over a 3‐year period. Nutrient and sediment transport were both dominated by subsurface flow, in particular through the B‐horizon, and this may seriously limit the surface‐runoff‐related functions of the riparian buffers. Riparian buffer trapping efficiencies were variable on an event basis and annual basis. The grass buffer reduced total phosphorus, filterable reactive phosphorus, total nitrogen and suspended sediment loads from surface runoff by 50 to 60%. The E. globulus buffer was not as effective, and total load reductions in surface runoff ranged between 10 and 40%. A key difference between the grass and E. globulus buffers was the seasonality of sediment and nutrient transport. Surface runoff, and therefore sediment and nutrient transport, occurred throughout the year in the E. globulus buffer, but only during the winter in the grass buffer. As a consequence of high summer nutrient and sediment concentrations, half the annual loads moving via surface runoff pathways through the E. globulus buffer were transported during intense summer storms. This study demonstrates that grass and E. globulus riparian buffers receiving runoff from pasture under natural rainfall can reduce sediment and nutrient loads from surface runoff. However, in this environment the B‐horizon subsurface flow is the dominant flowpath for nutrient transport through the riparian buffers, and this subsurface flow pathway carries contaminant loads at least three times greater than surface runoff. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

2.
Three main reservoirs were identified that contribute to the shallow subsurface flow regime of a valley drained by a fourth‐order stream in Brittany (western France). (i) An upland flow that supplied a wetland area, mainly during the high‐water period. It has high N‐NO3? and average Cl? concentrations. (ii) A deep confined aquifer characterized by low nitrate and low chloride concentrations that supplied the floodplain via flow upwelling. (iii) An unconfined aquifer under the riparian zone with high Cl? and low N‐NO3? concentrations where biological processes removed groundwater nitrate. This aquifer collected the upland flow and supplied a relict channel that controlled drainage from the whole riparian zone. Patterns of N‐NO3? and Cl? concentrations along riparian transects, together with calculated high nitrate removal, indicate that removal occurred mainly at the hillslope–riparian zone interface (i.e. first few metres of wetland), whereas dilution occurred in lower parts of the transects, especially during low‐water periods and at the beginning of recharge periods. Stream flow was modelled as a mixture of water from the three reservoirs. An estimation of these contributions revealed that the deep aquifer contribution to stream flow averaged 37% throughout the study period, while the contribution of the unconfined reservoir below the riparian zone and hillslope flow was more variable (from ca 6 to 85%) relative to rainfall events and the level of the riparian water table. At the entire riparian zone scale, NO3? removal (probably from denitrification) appeared most effective in winter, despite higher estimated upland NO3? fluxes entering the riparian zone during this period. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

3.
The riparian zone is in intimate contact with the river and, as such, is a critical zone for understanding hydrological problems. This paper presents a general modelling methodology for the assessment of riparian hydrological processes. It is applicable to a wide range of riparian spaces and incorporates current expertise in numerical methods. A core part of the modelling methodology is the random walk particle method (RWPM). We develop an RWPM as part of the ESTEL2D subsurface flow model, test it against analytical solutions and apply it to the simulation of parcels of water as they move through the riparian zone. The modelling methodology provides a new opportunity to assess fundamental hydrological process issues such as the proportioning of pre‐event and event water storm runoff, and reversals of flow in floodplains. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

4.
5.
Although riparian zones are well known to reduce nitrogen (N) and phosphorus (P) runoff to streams, they also have the potential to affect greenhouse gas (CO2, N2O, and CH4) fluxes to the atmosphere. Following large storms, soil biogeochemical conditions often become more reduced, especially in oxbow depressions and side channels, which can lead to hot moments of greenhouse gas production. Here, we investigate the impact of the remnants of Hurricane Irene and Tropical Storm Lee on riparian zone hydrology (water table: WT), and biogeochemistry (oxidation‐reduction potential [ORP], dissolved oxygen [DO], NO3?, PO43?, CO2, N2O, CH4). Results indicate that large storms have the potential to reset WT levels for weeks to months. Overbank flooding at our site following Irene and Lee led to the infiltration of well‐oxygenated water at depth (higher DO and ORP) while promoting the development of anoxic conditions within soil aggregates near the soil surface (increased N2O and CH4 fluxes). A short‐term increase in CO2 emission was observed following Irene at our site where aerobic respiration was water‐limited. Over a 2‐year period, an oxbow depression exhibited higher WT, higher N2O and CH4 fluxes (hot moment), higher CO2 fluxes (seasonal), and lower NO3? concentrations (seasonal) than the rest of the riparian zone. However, neither Irene, nor Lee, nor the oxbow depression significantly impacted PO43?. Dissolved organic carbon, ORP, and DO data illustrate the time‐lag (>20 years) between the creation of an oxbow depression and the development of reducing conditions despite clear differences in riparian zone and oxbow WT dynamics.  相似文献   

6.
Little is known about how active stream network expansion during rainstorms influences the ability of riparian buffers to improve water quality. We used aerial photographs to quantify stream network expansion during the wet winter season in five agricultural catchments in western Oregon, USA. Winter stream drainage densities were nearly two orders of magnitude greater than summer stream densities, and agricultural land use was much more abundant along transient portions (e.g. swales, road ditches) of stream networks. Water moving from agricultural fields into expanded stream networks during large hydrologic events has the opportunity to bypass downstream riparian buffers along perennial streams and contribute nonpoint‐source pollutants directly into perennial stream channels. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

7.
Riparian buffer zones in agriculture dominated watersheds play important roles in reducing nonpoint source pollution into aquatic ecosystems and are widely used as a Best Management Practice. Assessment of the effectiveness of riparian buffer zones by modeling method is widely used for watershed management as field measurement‐based assessment is difficult and expensive. The integration of Riparian Ecosystem Management Model (REMM) and Soil and Water Assessment Tool (SWAT) has been developed to simulate the effect of nonpoint source pollution reduction by riparian buffer zones at subbasin scale. However, there are problems in using the integrated model at subbasin scale, as the size of subbasin partition could affect the pollutant reduction rate by riparian buffers. In this study, we partitioned a large watershed with size of 1331 ha into sub‐watersheds with sizes of 666, 333, 166, 83, 51, and 29 ha, and then compared the different simulation results. We found that the modeling could yield more convergent results when the sub‐watersheds were partitioned into suitable size. In the studied area, the suitable sub‐watershed size was less than about 166 ha for runoff and nitrogen and 83 ha for sediment and phosphorus. Among the eight sub‐watersheds (partitioned based on the size of 166 ha), results showed that the effects of riparian buffers on runoff and nutrient loading varied drastically. The reduction rate varied from 0.26% to 30.13% for runoff, 29.4% to 74.07% for sediment, 9.61% to 57.85% for nitrogen, and 18.61% to 68.12% for phosphorus, respectively.  相似文献   

8.
Philippe Vidon 《水文研究》2012,26(21):3207-3215
Determining how riparian zone hydrological conditions may change in response to precipitation in various geomorphic settings is critical to determine the occurrence of hot moments of biogeochemical transformations for phosphorus, nitrogen, sulfate, mercury and greenhouse gases in these systems. The author investigate water table response to precipitation at a high temporal resolution (15 min) in a riparian zone located in a deeply incised glacial till valley (20 m) with approximately 2 m of alluvium over a confining layer, in Indiana, USA. During storms, larger water table fluctuations (approximately 100 cm) occurred near the stream than near the toe slope (10–25 cm). A quick rise in water table near the stream occurred for all storms, with partial flow reversals occurring for three of seven storms. The quick rise of the water table near the stream was associated with a decrease in hillslope water contributions to the stream during storms and the development of a water table down valley gradient for most storms. Water table fluctuations, groundwater flow velocities and electrical conductivity data indicated that riparian zone water table response to precipitation was primarily regulated by pressure wave processes. Regardless of the storm, high water tables persisted for at least 2 days after the cessation of precipitation. Although this suggests that high‐resolution precipitation data may be useful to quantify hot moments of biogeochemical transformation associated with high water tables in riparian zones, precipitation data alone are not sufficient to correctly estimate the magnitude of riparian water table level changes during storms. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

9.
The hydrology and nitrogen biogeochemistry of a riparian zone were compared before and after the construction of beaver dams along an agricultural stream in southern Ontario, Canada. The beaver dams increased surface flooding and raised the riparian water table by up to 1·0 m. Increased hydraulic gradients inland from the stream limited the entry of oxic nitrate‐rich subsurface water from adjacent cropland. Permeable riparian sediments overlying dense till remained saturated during the summer and autumn months, whereas before dam construction a large area of the riparian zone was unsaturated in these seasons each year. Beaver dam construction produced significant changes in riparian groundwater chemistry. Median dissolved oxygen concentrations were lower in riparian groundwater after dam construction (0·9–2·1 mg L?1) than in the pre‐dam period (2·3–3·9 mg L?1). Median NO3‐N concentrations in autumn and spring were also lower in the post‐dam (0·03–0·07 mg L?1) versus the pre‐dam period (0·1–0·3 mg L?1). In contrast, median NH4‐N concentrations in autumn and spring months were higher after dam construction (0·3–0·4 mg L?1) than before construction (0·13–0·14 mg L?1). Results suggest that beaver dams can increase stream inflow to riparian areas that limit water table declines and increase depths of saturated riparian soils which become more anaerobic. These changes in subsurface hydrology and chemistry have the potential to affect the transport and transformation of nitrate fluxes from adjacent cropland in agricultural landscapes. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

10.
Sugarcane is an annual crop with a dynamic canopy that changes over time mainly because of genetic adaptation. There is uncertainty about the temporal trends of throughfall (TF) in this important commercial crop. In the present paper, we used troughs to measure TF in a third and fourth ratoon and subsequently in a fourth and fifth ratoon. Additional measurements were carried out in an adjacent riparian forest. There were no significant differences between cycles of sugarcane, growth phases and riparian forest. The TF results for ratoon crop and riparian forest in 2011/2012 were 76% and 79.5% of gross rainfall, respectively, while in 2012/2013, they were 79% and 78%, respectively. However, TF was remarkably lower in the riparian forest relative to ratoon from the second half of the culm formation and elongation phase (280 days after harvest) until harvest. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
Abstract

The effects of acidic precipitation on stream chemistry were measured on an east-central Pennsylvania basin. When combined with flow and chemical mass balances, the data can help quantify hydrological source areas and their contributions to acidic storm hydrographs. For small storms on the well-buffered agricultural basin, small volumes of acidic precipitation falling directly on the stream surface react with more alkaline inflows from subsurface flow and surface runoff components to reduce streamflow pH temporarily by approximately one-half unit. During larger storms, the pH of surface runoff approaches that of precipitation, causing a relatively large acidic loading to the stream. However, this large input is buffered by a correspondingly larger subsurface flow component which results in stream pH reductions similar to those observed during the smaller events. Hydrological interpretations derived from a pH based mass balance are reinforced by a mass balance based on electrical conductivity and are consistent with the variable source area concept of basin hydrology.  相似文献   

12.
We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

13.
The management of the riparian zone has been suggested as a technique for controlling the amounts of phosphorus (P) entering watercourses draining pasture catchments. A study was therefore made of P entering a stream from various sources (rainfall, surface and subsurface derived runoff, direct fallout from aerial topdressing), with the object of providing a rational basis for the design of effective riparian management schemes. P entrained in surface runoff could account for virtually all of the P entering the stream during storms. Approximately 20 per cent of the annual P export from the catchment could be accounted for by direct aerial input of P to the stream during autumn fertilizer topdressing. More than 85 per cent of the P was exported from the catchment as particulate P. Stream sediment had higher P sorption capacities, and were enriched with P relative to the soils from which they were derived. There was a pronounced seasonal variation in sediment enrichment which could be predicted (r2 = 0.92) by the logarithm of the rainfall since fertilizer topdressing (LNFERT) and flood intensity. The amount of P lost in streamflow during any flood event was predicted (r2 = 0.94) by peak flow, seven day antecedent peak flow and LNFERT. Approximately 40 per cent of the 1.3 kg P ha?1 exported during 1981 occurred in four storms with recurrence intervals of more than three months. From a P budget compiled from nine events it was hypothesized that the stream acted as a net sink for P at baseflow and low-medium intensity floods but was a source of P at higher flood intensities. It was concluded that P losses from hill pasture catchments could be reduced by avoidance of direct application of P fertilizer to the stream channel, and by fencing out stock from seasonally saturated areas during periods of saturation. The ultimate success of the latter technique would depend on whether buffer vegetation could retain accumulated P during extreme storm events.  相似文献   

14.
The influence of riparian woodland on stream temperature, micro‐climate and energy exchange was investigated over seven calendar years. Continuous data were collected from two reaches of the Girnock Burn (a tributary of the Aberdeenshire Dee, Scotland) with contrasting land use characteristics: (1) semi‐natural riparian forest and (2) open moorland. In the moorland reach, wind speed and energy fluxes (especially net radiation, latent heat and sensible heat) varied considerably between years because of variable riparian micro‐climate coupled strongly to prevailing meteorological conditions. In the forested reach, riparian vegetation sheltered the stream from meteorological conditions that produced a moderated micro‐climate and thus energy exchange conditions, which were relatively stable between years. Net energy gains (losses) in spring and summer (autumn and winter) were typically greater in the moorland than the forest. However, when particularly high latent heat loss or low net radiation gain occurred in the moorland, net energy gain (loss) was less than that in the forest during the spring and summer (autumn and winter) months. Spring and summer water temperature was typically cooler in the forest and characterised by less inter‐annual variability due to reduced, more inter‐annually stable energy gain in the forested reach. The effect of riparian vegetation on autumn and winter water temperature dynamics was less clear because of the confounding effects of reach‐scale inflows of thermally stable groundwater in the moorland reach, which strongly influenced the local heat budget. These findings provide new insights as to the hydrometeorological conditions under which semi‐natural riparian forest may be effective in mitigating river thermal variability, notably peaks, under present and future climates. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.  相似文献   

15.
The Xiaolangdi Dam, completed in 2000, is second in scale in China to the Three Gorges Project. It has generated remarkable economic and social benefits but with profound impacts to the riverine and regional environments. This paper reports field monitoring of riparian groundwater in the Kouma section of the Yellow River to illustrate the interactions between dam‐regulated river flow and riparian groundwater. The results show that the hydrological condition in riparian zones downstream from the dam has changed from a typical wet–dry cycle to a condition of semi‐permanent dryness, resulting in degradation of the typical attributes and functions of the wetland ecosystem. Hydrological processes in the riparian zone have changed from a complex multiple flooding regime to a simple regime of dominant groundwater drainage towards the river, which only reverses temporarily during the water and sediment regulation period of the dam. Data on groundwater level and groundwater quality show that there are two key points, at ca 200 and 400 m from the river bank, which distinguish zones with different sensitivity to changes of river flow and indicate different interactions between river water and groundwater. The shallow groundwater quality also is negatively affected by the intensive agricultural development that has occurred since the dam was completed. Ecological restoration needs to be carried out to construct a protective natural riparian zone within ca 200 m from the river, this being an ecotone, which is key to the protection of both riparian groundwater and the river. The riparian zone from 200 to 400 m also should be treated as a transitional zone. In addition, ecologically sensitive agriculture and ecotourism organized by local communities would be beneficial in the area beyond 400 m. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

16.
Water flow in the soil–root–stem system was studied in a flooded riparian hardwood forest in the upper Rhine floodplain. The study was undertaken to identify the vertical distribution of water uptake by trees in a system where the groundwater is at a depth of less than 1 m. The three dominant ligneous species (Quercus robur, Fraxinus excelsior and Populus alba) were investigated for root structure (vertical extension of root systems), leaf and soil water potential (Ψm), isotopic signal (18O) of soil water and xylem sap. The root density of oak and poplar was maximal at a depth of 20 to 60 cm, whereas the roots of the ash explored the surface horizon between 0 and 30 cm, which suggests a complementary tree root distribution in the hardwood forest. The flow density of oak and poplar was much lower than that of the ash. However, in the three cases the depth of soil explored by the roots reached 1·2 m, i.e. just above a bed of gravel. The oak roots had a large lateral distribution up to a distance of 15 m from the trunk. The water potential of the soil measured at 1 m from the trunk showed a zone of strong water potential between 20 and 60 cm deep. The vertical profile of soil water content varied from 0·40 to 0·50 cm3 cm?3 close to the water table, and 0·20 to 0·30 cm3 cm?3 in the rooting zone. The isotopic signal of stem water was constant over the whole 24‐h cycle, which suggested that the uptake of water by trees occurred at a relatively constant depth. By comparing the isotopic composition of water between soil and plant, it was concluded that the water uptake occurred at a depth of 20 to 60 cm, which was in good agreement with the root and soil water potential distributions. The riparian forest therefore did not take water directly from the water table but from the unsaturated zone through the effect of capillarity. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

17.
Streambank erosion is a pathway for sediment and nutrient loading to streams, but insufficient data exist on the magnitude of this source. Riparian protection can significantly decrease streambank erosion in some locations, but estimates of actual sediment load reductions are limited. The objective of this research was to quantify watershed‐scale streambank erosion and estimate the benefits of riparian protection. The research focused on Spavinaw Creek within the Eucha‐Spavinaw watershed in eastern Oklahoma, where composite streambanks consist of a small cohesive topsoil layer underlain by non‐cohesive gravel. Fine sediment erosion from 2003 to 2013 was derived using aerial photography and processed in ArcMap to quantify eroded area. ArcMap was also utilized in determining the bank retreat rate at various locations in relation to the riparian vegetation buffer width. Box and whisker plots clearly showed that sites with riparian vegetation had on average three times less bank retreat than unprotected banks, statistically significant based on non‐parametric t‐tests. The total soil mass eroded from 2003 to 2013 was estimated at 7.27 × 107 kg yr.?1, and the average bank retreat was 2.5 m yr.?1. Many current erosion models assume that fluvial erosion is the dominant stream erosion process. Bank retreat was positively correlated with stream discharge and/or stream power, but with considerable variability, suggesting that mass wasting plays an important role in streambank erosion within this watershed. Finally, watershed monitoring programs commonly characterize erosion at only a few sites and may scale results to the entire watershed. Selection of random sites and scaling to the watershed scale greatly underestimated the actual erosion and loading rates. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

18.
To determine how soil frost changes flowpaths of runoff water along a hillslope, a transect consisting of four soil profiles directed towards a small stream in a mature forest stand was investigated at Svartberget, near Vindeln in northern Sweden. Soil temperature, unfrozen water content, groundwater level and snow depth were investigated along the transect, which started at the riparian peat, and extended 30 m upslope into mineral soils. The two, more organic‐rich profiles closest to the stream had higher water retention and wetter autumn conditions than the sandy mineral soils further upslope. The organic content of the soil influenced the variation in frost along the transect. The first winter (1995–96) had abnormally low snow precipitation, which gave a deep frost down to 40–80 cm, whereas the two following winters had frost depths of 5–20 cm. During winter 1995–96, the two organic profiles close to the stream had a shallower frost depth than the mineral soil profile higher upslope, but a considerably larger amount of frozen water. The fraction of water that did not freeze despite several minus degrees in the soil was 5–7 vol.% in the mineral soil and 10–15 vol.% in the organic soil. From the measurements there were no signs of perched water tables during any of the three snowmelt periods, which would have been strong evidence for changed water flowpaths due to soil frost. When shallow soil layers became saturated during snowmelt, especially in 1997 and 1998, it was because of rising groundwater levels. Several rain on frozen ground events during spring 1996 resulted in little runoff, since most of the rain either froze in the soil or filled up the soil water storage. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

19.
We examined how and why dominant peak-flow runoff-generation mechanisms differ among neighbouring headwater catchments. We monitored runoff and groundwater levels and performed terrain analyses in a granitic second-order catchment and its four neighbouring subcatchments in the Kiryu Experimental Watershed in Japan. Our analysis of lag times from peak rainfall to peak runoff suggests differences in the dominant peak-flow runoff-generation mechanisms among the five catchments. For two of the three zero-order catchments, with few perennial groundwater bodies, subsurface flow from hillslopes was the dominant mechanism at some events. However, the dominant mechanisms were channel precipitation and riparian runoff at almost all events in first- and second-order catchments and in the third zero-order catchment, which has a large perennial groundwater body over a bedrock depression in the riparian zone. In this zero-order catchment, the quick-flow ratio was the smallest of the five catchments because subsurface flow from the hillslope was buffered at the riparian zone. These facts suggest that the channel length, riparian buffering, and hillslope connectivity were the factors governing the different dominant peak-flow runoff-generation mechanisms among the catchments. Riparian buffering was affected, not only by surface topography, but also by bedrock topography and bedrock groundwater (BGW) dynamics. Our findings indicate that both of BGW dynamics and topography are important for catchment classification, and the relative importance of topography increases with the change from baseflow to stormflow. Furthermore, mismatching between a geographic source and a flow path resulted in different catchment classifications depending on the approach. Therefore, multiple approaches during both baseflow and stormflow periods are necessary for catchment classification to apply information obtained from one headwater catchment to other headwater catchments within the same region.  相似文献   

20.
The processes of stormflow generation were studied in a hill pasture catchment near Hamilton, New Zealand. Although rainfall was relatively evenly distributed throughout the year, stormflow was highly seasonal and over 65 per cent occurred during the winter. Three main processes contributing to stormflow were identified which could be related to soil type and physiographic position. On gleyed soils derived from rhyolitic colluvium, saturation overland flow was the dominant process. Hydrographs from ‘Whipkey’ throughflow troughs also indicated that there was a subsurface response (saturated wedge) from this soil type. On steeper convex slopes, more permeable soils were derived from weathered greywacke. The presence of ephemeral springs on the hillslopes and direct observation during storm events indicated that storm runoff was generated as return flow from this soil. It was noted that nitrate concentrations from subsurface sources were 5–10 times higher than surface runoff. This difference in concentration was utilized in a chemical mixing equation which partitioned stormflow sources. This was compared with the stormflow predicted from rain falling on to saturated areas. There was good agreement between the two models for winter-spring events with respect to the volumes of surface runoff predicted, however the saturated areas model underestimated total stormflow. The results of the study are briefly discussed in terms of the potential for water quality management.  相似文献   

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

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