首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
Formation and development of alternate bars in an engineered mountainous reach of the Arc River, France, is studied using photo analysis, 1D modelling and by applying theoretical and empirical models for alternate bar systems. Alternate bars already existed in the 80s in the form of a stable confined wandering system. In 1994, the river bed was flattened after engineering works. However, aerial photographs and cross-sectional profiles show that bars rapidly recovered within a few years. The alternate bar system evolved rapidly with a reduction of the number of bars and so an increasing bar length. The width-to-depth ratio, the slope change, the bend upstream of the reach, and the sediment supplies are the main controls of bar formation and evolution. The system appears to lead to force bars due to the bend but also due to a bridge in the downstream part of the reach. Nevertheless some free mobile bars are still observed in the middle of the reach. A discussion on the alternate bar formation is provided using empirical and analytical models. Finally, impacts of low flows and vegetation seem to be significant in the stabilization of the system toward a confined wandering system as observed before the engineering works.  相似文献   

2.
This paper analyses types and rates of change in river meander morphology and the links between mechanisms of change and emergent behaviour of planform morphology. It uses evidence of four dates of aerial photography combined with annual field mapping and ground photography to examine the morphological changes and mechanisms of change in a series of bends on an active meandering river, the River Dane in NW England, over a 25 year period. This unique data set allows insight into the spatial and temporal variability of bank line movement and component processes. Bank lines were mapped photogrametrically from air photos of 1984, 1996, 2001 and 2007 and the digitised courses compared in ArcGIS to produce calculations of erosional and depositional areas and rates. Most bends exhibit morphological change that largely follows the autogenic sequence, identified in qualitative models of meander development, from low sinuosity curves through simple symmetric and asymmetric bends to compound forms with lobe development in the apex region. Rates of erosion and bankline movement increase through this sequence until the compound phase. Relationships of amounts of movement to various curvature measures of bend morphology are complex. Several new loops, distinct from compound bend behaviour, have developed during the study period in formerly straight sections. Mechanisms of morphological change are illustrated for four types of bends: new, rapid growth bend; sharp‐angled bend with mid‐channel bar development; symmetric migrating bend; and simple to compound bend development. The changes take place in phases that are not simply related to discharge but to inherent sequences and feedbacks in development of bars and bend morphology and timescales for these are identified. Overall, emergent behaviour of systematic planform change, moderated by channel confinement and boundary features, is produced from spatially and temporally varied channel processes. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

3.
A sustained dynamic inflow perturbation and bar–floodplain conversion are considered crucial to dynamic meandering. Past experiments, one-dimensional modelling and linear theory have demonstrated that the initiation and persistence of dynamic meandering require a periodic transverse motion of the inflow. However, it remains unknown whether the period of the inflow perturbation affects self-formed meander dynamics. Here, we numerically study the effect of the inflow perturbation period on the development and meander dynamics of a chute-cutoff-dominated river, which requires two-dimensional modelling with vegetation forming floodplain on bars. We extended the morphodynamic model Nays2D with growth and mortality rules of vegetation to allow for meandering. We tested the effect of a transversely migrating inflow boundary by varying the perturbation period between runs over an order of magnitude around typical modelled meander periods. Following the cutoff cascade after initial meander formation from a straight channel, all runs with sufficient vegetation show series of growing meanders terminated by chute cutoffs. This generates an intricate channel belt topography with point bar complexes truncated by chutes, oxbow lakes, and scroll-bar-related vegetation age patterns. The sinuosity, braiding index and meander period, which emerge from the inherent biomorphological feedback loops, are unrelated to the inflow perturbation period, although the spin-up to dynamic equilibrium takes a longer time and distance for weak and absent inflow perturbations. This explains why, in previous experimental studies, dynamic meandering was only accomplished with a sustained upstream perturbation in flumes that were short relative to the meander wavelength. Our modelling of self-formed meander patterns is evidence that scroll-bar-dominated and chute-cutoff-dominated meanders develop from downstream convecting instabilities. This insight extends to many more fluvial, estuarine and coastal systems in morphological models and experiments, which require sustained dynamic perturbations to form complex patterns and develop natural dynamics. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.  相似文献   

4.
Scroll bars across a 65-km stretch of the Trinity River in Texas, USA were studied using LiDAR data as well as with a series of 11 trenches spread out across the survey area. We conclude that scroll bars are levees that are deposited along the inner banks of these meandering river bends. Scroll bar crests were found to have similar elevations to those of outer bank levee crests, implying that they are constructional features that create positive topographic relief above the elevation of the floodplain. Trenches reveal that scroll bars are built from reworked suspended sediment, with common ripple-scale cross stratification, planar laminations and muddy bioturbated layers – characteristics often associated with levee sedimentation in other systems. LiDAR observation of the erosion of scroll bars by bed material transport during flood implies that scroll bar spacing is an imperfect proxy for estimating overall channel migration rates. In addition, interspersed lenses of coarser sediment with dune-scale cross stratification represent the stratigraphic record of these erosional events and suggest that erosion of the channel-ward edge of the scroll bar is not uncommon. Preservation of scroll bars is unlikely, given that they are responsible for an average of only the uppermost 12% of the total inner bank relief. We suggest that misidentification of point bar lateral accretion surfaces as scroll bars is common and can lead to issues with reconstructing channel properties due to systematic differences between point bar and scroll bar planform geometries. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.  相似文献   

5.
Channel bars and banks strongly affect the morphology of both braided and meandering rivers. Accordingly, bar formation and bank erosion processes have been greatly explored. There is, however, a lack of investigations addressing the interactions between bed and bank morphodynamics, especially over short timescales. One major implication of this gap is that the processes leading to the repeated accretion of mid‐channel bars and associated widenings remain unsolved. In a restored section of the Drau River, a gravel‐bed river in Austria, mid‐channel bars have developed in a widening channel. During mean flow conditions, the bars divert the flow towards the banks. One channel section exhibited both an actively retreating bank and an expanding mid‐channel bar, and was selected to investigate the morphodynamic processes involved in bar accretion and channel widening at the intra‐event timescale. We repeatedly surveyed riverbed and riverbank topography, monitored riverbank hydrology and mounted a time‐lapse camera for continuous observation of riverbank erosion processes during four flow events. The mid‐channel bar was shown to accrete when it was submerged during flood events, which at the subsequent flow diversion during lower discharges narrowed the branch along the bank and increased the water surface elevation upstream from the riffle, which constituted the inlet into the branch. These changes of bed topography accelerated the flow along the bank and triggered bank failures up to 20 days after the flood events. Four analysed flow events exhibited a total bar expansion from initially 126 m2 to 295 m2, while bank retreat was 6 m at the apex of the branch. The results revealed the forcing role of bar accretion in channel widening and highlighted the importance of intra‐event scale bed morphodynamics for bank erosion, which were summarized in a conceptual model of the observed bar–bank interactions. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
7.
In this paper we use multiple field surveys spanning several decades to systematically evaluate the geomorphic consequences of a change in flow hydraulics from uniform flow to backwater flow for the lower Trinity River in east Texas, USA. Spatial changes in lateral migration rate, channel geometry, and point bar size correspond to two distinct geomorphic zones. Within the upstream uniform flow reach, the river channel is defined by fully developed point bars and a high rate of lateral channel migration. This zone transitions where the median channel bottom elevation drops below sea level. At this point flow is affected by the backwater influence of the Trinity Bay water surface elevation, as opposed to being bed slope control dominated. The change in hydraulics within the backwater zone is reflected in the channel morphology, which is characterized by smaller point bars, narrower and more symmetrical cross-sectional channel geometry, lower channel migration rates, and little to no bend deformation or cutoffs. Studying the connection between channel geometry, river bend kinematics, sediment transport, and fluid mechanics in each zone provides a deeper understanding of the relationship between channel shape and river mechanics. © 2019 John Wiley & Sons, Ltd.  相似文献   

8.
River meandering has been extensively investigated. Two fundamental features to be explored in order to make further progress are nonlinearity and unsteadiness. Linear steady models have played an important role in the development of the subject but suffer from a number of limits. Moreover, rivers are not steady systems; rather their states respond to hydrologic forcing subject to seasonal oscillations, punctuated by the occurrence of flood events. We first derive a classification of river bends based on a systematic assessment of the various physical mechanisms affecting their morphodynamic equilibrium and their evolution in response to variations of hydrodynamic forcing. Using the database by Lagasse et al. ( 2004 ) we also show that natural meanders are typically mildly curved and long, i.e. such that both the centrifugal and the topographic secondary flows are weak, but they are almost invariably nonlinear. We then review some recent developments which allow us to treat analytically the flow and bed topography of mildly curved and long nonlinear bends subject to steady forcing, taking advantage of the fact that flow and bed topography in mildly curved long bends are slowly varying. Results show that nonlinearity has a number of consequences: most notably damping of the morphodynamic response and upstream shifting of the location of the nonlinear peak of the flow speed. Next we extend the latter model to the case of unsteady forcing. Results are found to depend crucially on the ratio between the flood duration and a morphodynamic timescale. It turns out that, in a channel subject to a repeated sequence of floods, the system reaches a dynamic equilibrium. We conclude the paper discussing how the present assessment relates to the debate on meander modelling of the late 1980s and suggesting what we see as promising lines of future developments.  相似文献   

9.
Along the lower reaches of the Waipaoa River, New Zealand, cross‐section survey data indicate there was a 23 per cent decrease in bankfull width and a 22 per cent reduction in channel cross‐section area between 1948 and 2000, as the channel responded to increased inputs of fine (suspended) sediment following deforestation of the headwaters in late C19 and early C20. We determined the bankfull discharge within a ~39 km long reach by routing known discharges through the one‐dimensional MIKE 11 flow model. The model runs suggest that the bankfull discharge varies between ~800 and ~2300 m3 s?1 and that the average recurrence interval is 4 ± 2 years on the annual maximum series; by contrast, the effective flow (360 m3 s?1) is equaled or exceeded three times a year. The variability in bankfull discharge arises because the banks tend to be lower in places where flood flows are constricted than in reaches where overbank flow is dispersed over a wide area, and because scour has counteracted aggradation in some locations. There is no downstream variation in Shields stress, or in relative shear stress, within the study reach. Bankfull shear stress is, on average, five times greater than the shear stress required to initiate motion. At the effective discharge it is more than twice the threshold value. The effective discharge probably has more relevance than the bankfull discharge to the overall picture of sediment movement in the lower reaches of the Waipaoa River but, because width is constrained by the stability and resistance of the bank material to erosion during high flows that also scour the bed, the overall channel geometry is likely determined by discharges at or near bankfull. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

10.
Bar colonization by vegetation and subsequent island formation is a key bio-geomorphological process in fluvial landscape evolution. Here we investigate morphological and ecological evolution of river islands over timescales from single floods to decades, focusing on islands initiated by deposited trees that sprout to form vegetated patches. On a braided reach of the high-energy Tagliamento River, Italy, we monitored 30 pioneer islands of 1 to 17 years age in comparison with unvegetated bar surfaces, open areas between islands, and established island surfaces. We integrated morphological, surface sediment and vegetation properties of islands initiated by different flood events, combining evidence from remotely-sensed and ground observations, flow and climate time series. At a decadal timescale, pioneer islands aggrade rapidly to the elevation of the mean annual flood, showing a steady increase in vegetation canopy height, fining of surface sediments from predominantly gravel to silty-sand with a notable clay and organic fraction. The standing vegetation included over 130 species, with the largest number on island surfaces of intermediate elevation and flood disturbance. As islands age, standing vegetation becomes comprised mainly of competitor species with transient seed banks and typical of woodland, scrub, pasture and wetland habitats, whereas the winter seedbank is dominated on all surfaces by ruderal species with persistent seedbanks, mainly associated with aquatic, wetland, pasture, arable and wasteland habitats. At shorter timescales, the biogeomorphological trajectory of pioneer islands is initiated by large flood events that control the elevation of deposited trees, and subsequent flows that control tree survival and establishment. Island morphological evolution depends on the frequency-magnitude of sediment and seed delivery and redistribution by flood and possibly wind events, whereas island ability to retain sediments reflects the degree of vegetation establishment, which in the short-term may vary with seasonal to annual moisture supply, substrate characteristics and climatic growth conditions. © 2018 John Wiley & Sons, Ltd.  相似文献   

11.
Digital elevation models and topographic pro?les of a beach with intertidal bar and trough (ridge‐and‐runnel) morphology in Merlimont, northern France, were analysed in order to assess patterns of cross‐shore and longshore intertidal bar mobility. The beach exhibited a pronounced dual bar–trough system that showed cross‐shore stationarity. The bars and troughs were, however, characterized by signi?cant longshore advection of sand under the in?uence of suspension by waves and transport by strong tide‐ and wind‐driven longshore currents. Pro?le changes were due in part to the longshore migration of medium‐sized bedforms. The potential for cross‐shore bar migration appears to be mitigated by the large size of the two bars relative to incident wave energy, which is modulated by high vertical tidal excursion rates on this beach due to the large tidal range (mean spring tidal range = 8·3 m). Cross‐shore bar migration is also probably hindered by the well‐entrenched troughs which are maintained by channelled high‐energy intertidal ?ows generated by swash bores and by tidal discharge and drainage. The longshore migration of intertidal bars affecting Merlimont beach is embedded in a regional coastal sand transport pathway involving tidal and wind‐forced northward residual ?ows affecting the rectilinear northern French coast in the eastern English Channel. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

12.
Sedimentary deposits provide records of environmental change quantifying erosion fluxes conditioned by natural and anthropogenic disturbances. These fluxes are lagged by internal storage, particularly within floodplains, complicating reconstruction of environmental changes. The time sediment remains in storage underpins the interpretation of sedimentary records and accurate monitoring of pollutant fluxes. Turnover time is a measure of the timeframe to erode every floodplain surface. CAESAR-Lisflood is used to simulate fluvial evolution at reach scale, providing a basis for quantifying environmental changes on the timescales of sediment storage. We evaluate the accuracy of CAESAR-Lisflood simulations of channel changes and turnover times for alluvial floodplains using historical channel changes reconstructed for 10 reaches in northern England to quantify model accuracy in replicating mean annual erosion, deposition and channel lateral migration rates, alongside planform morphology. Here, a split-sample testing approach is adopted, whereby five of the reaches were calibrated and the resulting parameter values were applied to the other reaches to evaluate the transferability of parameter settings. The lowest overall integrated error identified the best-fit simulations and showed that modelled process rates were within ~25–50% of rates from historical reconstructions, generally. Calibrated parameters for some reaches are widely transferable, producing accurate geomorphic changes for some uncalibrated sites. However, large errors along some reaches indicate that reach-specific parameterization is recommended. Turnover times are underpinned by the assumption that areas of floodplain previously unvisited by the channel are reworked. This assumption has been challenged by studies that show floodplain (re)occupation rates vary spatially. However, this limitation is less important for the short-duration simulations presented here. The simulations reconstruct floodplain turnover times estimated by mapped rates mostly successfully, demonstrating the potential applicability of calibrated parameters over much longer timescales. Errors in the form of under-predicted erosion rates propagated, resulting in over-predicted turnover times by even greater magnitudes. © 2020 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.  相似文献   

13.
Riparian vegetation influences hydraulic and morphodynamic river processes and may contribute to sediment stabilization. In turn, vegetation recruitment and growth on non‐cohesive fluvial deposits strongly depends on river hydrology and the ability of roots to develop and to anchor efficiently to resist flow erosion. In this paper, we examine the above‐ground and the below‐ground seasonal growth dynamics of Salix cuttings in relation to local river hydrodynamics and morphodynamics, on the basis of a detailed and unique data set. During the two season‐long campaigns in 2009 and 2010, 1188 and 1152 cuttings, respectively, were organized in square plots and planted on a gravel island of the restored reach of the River Thur (Neunforn, Thurgau, Switzerland). Each year, all cuttings were monitored almost regularly from the beginning until the end of the growing season (April–September). Root development statistics were also obtained from high‐resolution scanner analysis of carefully uprooted samples from selected plots. Our results show how cutting survival and the nature and strength of correlations between island topography and cutting growth statistics depend on river hydrology. An empirical functional form that links root development based on the measured main stem length is then proposed for predictive purposes. Cutting mortality following flood events is shown to depend nonlinearly on both erosion and deposition processes, whereas it appears more linearly related to the magnitude of the bed shear stress distribution generated by the maximum seasonal flood. This analysis allows an identification of an important threshold for plant survival within different erosion and deposition regimes, which explains the spatial and temporal distribution of the surviving cuttings within the plots. These results have practical implications, for instance, for evaluating, planning and managing the use of riparian trees in restoration projects. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
15.
A digital elevation model (DEM) of a fluvial environment represented landform surface variability well and provided a medium for monitoring morphological change over time. Elevation was measured above an arbitrary datum using a ground‐based three‐dimensional tacheometric survey in two reaches of the River Nent, UK, in July 1998, October 1998 (after flood conditions) and June 1999. A detailed geostatistical analysis of the elevation data was used to model the spatial variation of elevation and to produce DEMs in each reach and for each survey period. Maps of the difference in elevation were produced and volumetric change was calculated for each reach and each survey period. The parameters of variogram models were used to describe the morphological character of each reach and to elucidate the linkages between process and the form of channel change operating at different spatial and temporal scales. The analysis of channel change on the River Nent shows the potential of geostatistics for investigating the magnitude and frequency of geomorphic work in other rivers. A flood modified the channel features, but low magnitude and high frequency flows rationalized the morphology. In spite of relatively small amounts of net flux the channel features changed as a consequence of the reworking of existing material. The blocking of chute entrances and redirection of the channel had a considerable effect on the behaviour of the channel. Such small changes suggested that the distributary system was sensitive to variation in sediment regime. Plots of the kriging variances against sampling intervals were used to quantify the temporal variation in sampling redundancy (ranging between ?11 per cent and +93 per cent). These curves illustrated the importance of bespoke sampling designs to reduce sampling effort by incorporating anisotropic variation in space and geomorphic information on flow regime. Variation in the nugget parameter of the variogram models was interpreted as sampling inaccuracy caused by variability in particle size and is believed to be important for future work on surface roughness. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

16.
The Rhine bifurcation at Pannerden forms the major distribution point for water supply in the Netherlands, distributing not only water and sediment but also flooding risks and navigability. Its morphological stability has been a concern for centuries. We present experiences from more than two decades of numerical morphological modelling of this bifurcation with a gravel–sand bed and a meandering planform. Successive computations have shown the importance of upstream approach conditions, the necessity to include physical mechanisms for grain sorting and alluvial roughness, and the need to assume a thicker active layer of the river bed than is suggested by laboratory flume experiments using a constant discharge. The active layer must be thicker in the model to account for river bed variations due to higher‐frequency discharge variations that are filtered out in morphological modelling. We discuss limitations in calibration and verification, but argue that, notwithstanding these limitations, 2D and 3D morphological models are valuable tools, not only for pragmatic applications to engineering problems, but also for revealing the limitations of established knowledge and understanding of the relevant physical processes. The application of numerical models to the Pannerden bifurcation appeared to reveal shortcomings in established model formulations that do not pose particular problems in other cases. This application is therefore particularly useful for setting the agenda for further research. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

17.
Channelization of the severely polluted Odra and Vistula Rivers in Poland induced intensive accumulation of fine‐grained deposits rich in organic matter and heavy metals. These sediments have been identified in vertical profiles in a narrow zone along river banks both in groyne‐created basins and on the floodplain. Grain size, organic matter, zinc (Zn), lead (Pb), copper (Cu) content and cesium‐137 (137Cs) was used for sediment dating and, stratigraphy and chemistry have been diagnostic features for these deposits, named industrial alluvium. In the most polluted river reaches stabilized by bank reinforcements and groynes, 2‐m‐thick slack water groyne deposits are composed of uniform strata of polluted silts with organic matter content over 10%, Zn content over 1000 mg/kg and average Cu and Pb over 100 mg/kg. The average rate of sediment accretion in groynes is higher than on the floodplain and reaches 5 cm/yr. Stratification which appears at higher levels in the groyne fields and on the levees reflects a change from in‐channel to overbank deposition and is typified by dark layers separated by bright, sandy, and less polluted strata. Stratified, 4‐m‐thick, sediment sequences have been found in groyne fields of incised river reaches. The average rate of sediment accretion in these reaches is of the order of 5 cm/yr. In stable and relatively less polluted river reaches, vertical‐accretion organic deposits are finely laminated and the average rate of deposition amounts to a few millimeters per year. Investigations indicate that groyne construction favors conditions for long‐term storage of sediments at channel banks. For this reason, groynes should be considered as structures that efficiently limit sudden release of sediment‐associated heavy metals stored in channels and in floodplains of the historically polluted rivers. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

18.
Alluvial channels arise through the interaction between morphology, hydraulics, and sediment transport, known as the ‘fluvial trinity’. Over relatively short timescales where climate and geology are fixed but discharge and sediment supply may vary, this process facilitates adjustments towards steady state, where the system oscillates around a mean condition. The relationship between changes in conditions and geomorphic response may be highly complex and nonlinear, especially in systems with multiple modes of adjustment. This study examines the adjustment of an erodible channel with fixed banks and a widely graded sediment mixture to successive increases in discharge. With each increase in discharge, components of the fluvial trinity adjusted towards a steady state. Particularly at relatively low discharges, adjustments were controlled by intrinsic thresholds and highlighted important morphodynamic processes. Notably, there was a strong interplay between channel morphology and sediment transport, and an effect whereby larger-than-average grains controlled channel deformation. These two processes occurred at the bar scale and were highly spatialised, which has two important implications: (1) reach-averaged representations of process provide only partial insight into morphodynamics; and (2) models of rivers that suppress these process feedbacks and size-dependent transport may not replicate morphodynamics that typically occur in field conditions. The experiments provide quantitative evidence for conceptual models describing exponential approaches towards steady state and the potential for transiency if disturbance frequency exceeds the recovery time. They also highlight how in natural rivers, particularly those with greater degrees of freedom for adjustment (notably, lateral adjustment and meandering), continuous changes in discharge may lead to nonlinear rather than steady-state behaviour. In these settings, more holistic analytical frameworks that embrace different aspects of the system are critical in understanding the direction, magnitude and timing of channel adjustments.  相似文献   

19.
20.
Mountain rivers can be subject to strong constraints imposed by changes in gradient and grain size supplied by processes such as glaciation and rockfall. Nonetheless, adjustments in the channel geometry and hydraulics of mountain rivers at the reach scale can produce discernible patterns analogous to those in fully alluvial rivers. Mountain rivers can differ in that imposed reach‐scale gradient is an especially important control on reach‐scale channel characteristics, as indicated by examination of North St Vrain Creek in Colorado. North St Vrain Creek drains 250 km2 of the Rocky Mountains. We used 25 study reaches within the basin to examine controls on reach‐scale channel geometry. Variables measured included channel geometry, large woody debris, grain size, and mean velocity. Drainage area at the study reaches ranged from 2·2 to 245 km2, and gradient from 0·013 to 0·147 m m?1. We examined correlations among (1) potential reach‐scale response variables describing channel bankfull dimension and shape, hydraulics, bedform wavelength and amplitude, grain size, ?ow resistance, standard deviation of hydraulic radius, and volume of large woody debris, and (2) potential control variables that change progressively downstream (drainage area, discharge) or that are likely to re?ect a reach‐speci?c control (bed gradient). We tested the hypothesis that response variables correlate most strongly with local bed gradient because of the segmented nature of mountain channels. Results from simple linear regression analyses indicate that most response variables correlate best with gradient, although channel width and width/depth ratio correlate best with discharge. Multiple regression analyses using Mallow's Cp selection criterion and log‐transformation of all variables produced similar results in that most response variables correlate strongly with gradient. These results suggest that the hypothesis is partially supported: channel bed gradient is likely to be a good predictor for many reach‐scale response variables along mountain rivers, but discharge is also an important predictor for some response variables. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

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

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