首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
For 27 hillslope profiles located within four first-order drainage basins on the dip slope of the South Downs in East Sussex, contour curvature (P), measured in degrees per 100 m, is highly variable within any one hillslope. Hillslopes may be described as wholly concave in plan where P is less than ?40, wholly convex in plan where P is greater than +40 and convexo-concave in plan where P lies between ?40 and +40. One measure of P taken at the profile sampling line is found to give a better estimate of P than measuring P at the steepest point on the hillslope. P is found to be significantly negatively correlated with slope in 38 per cent of cases, correlated with distance from the divide in 55 per cent of cases (positively where P <0 and negatively where P <0) and positively correlated with profile curvature in 70 per cent of cases. Implications of these correlations for hillslope process-response models are discussed.  相似文献   

2.
The areally weighted surface erosion for Niwot Ridge, an alpine interfluve in the Colorado Front Range, is 10−1 mm/y. This may be subdivided into rates for three generalized cover types: tundra meadow, 10−2 mm/y; dry tundra, 10−1 mm/y; late-lying snow patches, 10° mm/y. Tundra meadow (about 50 per cent of the interfluve area) yields about 5 per cent of the eroded material; dry tundra (35 per cent of the area) contributes slightly less than 50 per cent of the eroded material; while nivation hollows occupied by late-lying snow patches occupy only about 3 per cent of the area they contribute 50 per cent of the eroded material. The bulk of the surficial erosion is accomplished between June and September, primarily by rainsplash, except where snowmelt is important. The overall estimated surface lowering rate presented here is substantially higher than those reported previously.  相似文献   

3.
Karst solution processes are investigated on Oligocene limestones in the Waitomo district, west central North Island, New Zealand. Estimates of the inputs, throughputs and outputs of water and dissolved calcium and magnesium in two drainage basins were used to establish the rate of limestone solution by autogenic waters. The best estimate for solution loss from the basins during the study year is 69 m3/km2. The potential measurement errors inherent in each parameter used in the erosion rate computations were assessed and the probable maximum and minimum erosion rates were estimated to be 88 and 61 m3/km2. In both basins approximately 67 per cent of the annual solute load is transported by flows greater than the mean annual discharge, over 15 per cent being transported by flood flows that are exceeded only 5 per cent of the time. Almost half of the annual load is transported during the three winter months (June-August), but no one month accounts for more than 18 per cent or less than 2.7 per cent of the annual total. Approximately 37 per cent of solution takes place within the soil profile, and most of the remainder is concentrated in 5–10 m of weathered bedrock (the subcutaneous zone) beneath this. Thus, it is likely that at least 85 per cent of the total solutional erosion contributes to the surface lowering of soil and bedrock.  相似文献   

4.
Predicting sediment yield at the catchment scale is one of the main challenges in geomorphologic research. The application of both physics‐based models and regression models has until now not provided very satisfying results for prediction of sediment yield for medium to large sized catchments (c. >50 km2). The explanation for this lies in a combination of the large data requirements of most models and a lack of knowledge to describe all processes and process interactions at the catchment scale. In particular, point sources of sediment (e.g. gullies, mass movements), connectivity and sediment transport remain difficult to describe in most models. From reservoir sedimentation data of 44 Italian catchments, it appeared that there was a (non‐significant) positive relation between catchment area and sediment yield. This is in contrast to what is generally expected from the theory of decreasing sediment delivery rates with increasing catchment area. Furthermore, this positive relation suggests that processes other than upland erosion are responsible for catchment sediment yield. Here we explore the potential of the Factorial Scoring Model (FSM) and the Pacific Southwest Interagency Committee (PSIAC) model to predict sediment yield, and indicate the most important sediment sources. In these models different factors are used to characterize a drainage basin in terms of sensitivity to erosion and connectivity. In both models an index is calculated that is related to sediment yield. The FSM explained between 36 and 61 per cent of the variation in sediment yield, and the PSIAC model between 57 and 62 per cent, depending on the factors used to characterize the catchments. The FSM model performed best based on a factor to describe gullies, lithology, landslides, catchment shape and vegetation. Topography and catchment area did not explain additional variance. In particular, the addition of the landslide factor resulted in a significantly increased model performance. The FSM and PSIAC model both performed better than a spatially distributed model describing water erosion and sediment transport, which was applied to the same catchments but explained only between 20 and 51 per cent of the variation in sediment yield. Model results confirmed the hypothesis that processes other than upland erosion are probably responsible for sediment yield in the Italian catchments. A promising future development of the models is by the use of detailed spatially distributed data to determine the scores, decrease model subjectivity and provide spatially distributed output. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

5.
A sediment budget was developed for the 1.7 km2 Maluna Creek drainage basin located in the Hunter Valley, New South Wales, Australia, for the period 1971-86. the impact of viticulture, which commenced at Maluna in 1971, was studied using erosion plots, with caesium-137 as an indicator of both soil erosion and sedimentation. Two methods were used to estimate vineyard soil losses from caesium-137 measurements. Sediment output from the catchment was measured for three years, and extrapolated from readings taken at a nearby long-term stream flow gauging station for the remaining 13 years. Relative amounts of soil loss from forest (60 per cent basin area), grazing land (30 per cent) and vineyards (10 per cent) were calculated. Soil losses by rain splash detachment were ten times greater from bare/cultivated sufaces than from the forest. Erosion plots of area 2 m2 showed no significant differences in soil loss between forest and grassland but, under bare soil, losses were 100 times greater. the 137Cs method was employed to calculate net soil loss from all vineyard blocks using both a previously established calibration curve and a proportional model. the latter method gave estimates of soil loss which were 3-9 times greater than by the calibration curve, and indicated that average soil losses from the vineyard were equivalent to 62 t ha?1 y?1 (1971-86). It was estimated that the forest contributed 1-8 per cent, the grazing land 1.6 per cent, and the vineyard 96.6 per cent of the total soil loss during that period. Sediment storages within the fluvial system adjacent to the vineyard ws 9460 t for the period, whereas sediment output was equivalent to 215 t km?1 y?1. Independent measurements of soil erosion, storage, and output showed that 56 per cent of the eroded sediment remained in the catchment, and 34 per cent was transported out by Maluna Creek. the budget was able to be balanced to within 10 per cent.  相似文献   

6.
The influence of woody debris on channel morphology and aquatic habitat has been recognized for many years. Unlike sediment, however, little is known about how wood moves through river systems. We examined some dynamics of wood transport in streams through a series of flume experiments and observed three distinct wood transport regimes: uncongested, congested and semi-congested. During uncongested transport, logs move without piece-to-piece interactions and generally occupy less than 10 per cent of the channel area. In congested transport, the logs move together as a single mass and occupy more than 33 per cent of the channel area. Semi-congested transport is intermediate between these two transport regimes. The type of transport regime was most sensitive to changes in a dimensionless input rate, defined as the ratio of log volume delivered to the channel per second (Qlog) to discharge (QW); this ratio varied between 0·015 for uncongested transport and 0·20 for congested transport. Depositional fabrics within stable log jams varied by transport type, with deposits derived from uncongested and semi-congested transport regimes having a higher proportion of pieces orientated normal to flow than those from congested transport. Because wood input rates are higher and channel dimensions decrease relative to piece size in low-order channels, we expect congested transport will be more common in low-order streams while uncongested transport will dominate higher-order streams. Single flotation models can be used to model the stability of individual pieces, especially in higher-order channels, but are insufficient for modelling the more complex intractions that occur in lower-order streams. © 1997 John Wiley & Sons, Ltd.  相似文献   

7.
A one-day field investigation on an unvegetated backbeach documents the importance of surface sediment drying to aeolian transport. Surface sediments were well sorted fine sand. Moisture content of samples taken in the moist areas on the backbeach varied from 2·9 to 9·2 per cent. Lack of dry sediment inhibited transport prior to 08:50. By 09:10 conspicuous streamers of dry sand moved across the moist surface. Barchan-shaped bedforms, 30 to 40 mm high and composed of dry sand (moisture content <0·10 per cent), formed where sand streamers converged. The surface composed of dry sand increased from 5 per cent of the area of the backbeach at 09:50 to 90 per cent by 12:50 Mean wind speeds were beetween 5·6 and 8·6 m s−1 at 6 m above the backbeach. Corresponding shear velocities were always above the entrainment threshold for dry sand and below the threshold for the moist sand on the backbeach. Measured rates of sand trapped (by vertical cylindrical traps) increased during the day relative to calculated rates. The measured rate of sand trapped on the moist foreshore was higher than the rate trapped on the backbeach during the same interval, indicating that the moist foreshore (moisture content 18 per cent) was an efficient transport surface for sediment delivered from the dry portion of the beach upwind. Measured rates of sand trapped show no clear relationship to shear velocities unless time-dependent surface moisture content is considered. Results document conditions that describe transport across moist surfaces in terms of four stages including: (1) entrainment of moist sediment from a moist surface; (2) in situ drying of surface grains from a moist surface followed by transport across the surface; (3) entrainment and transport of dry sediment from bedforms that have accumulated on the moist surface; and (4) entrainment of sand from a dry upwind source and transport across a moist downwind surface. © 1997 John Wiley & Sons, Ltd.  相似文献   

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

9.
Thirty‐six runoff plot experiments provide data on flow depths, speeds, and Darcy–Weisbach friction coefficients (f) on bare soil surfaces, and surfaces to which were added sufficient extra plant litter or surface stones to provide projected cover of 5, 10 and 20 per cent. Precision flow depth data were derived with a computer‐controlled gantry and needle gauge for two different discharges for each plot treatment. Taking a fixed flow intensity (Reynolds number, Re = 150) for purposes of comparison shows means of f = 17·7 for bare soil surfaces, f = 11·4 for added stone treatments, and f = 23·8 for added litter treatments. Many individual values of f for stone treatments are lower than for the bare soil surface, but all litter treatments show increases in fcompared to bare soil. The lowering of f in stone treatments relates to the submerged volume that the stones occupied, and the associated concentration of flow onto a smaller part of the plot surface. This leads to locally higher flow intensities and lower frictional drag along threads of flow that the obstacles create. Litter causes higher frictional drag because the particles are smaller, and, for the same cover fraction, are 100 times more numerous and provide 20 times the edge or perimeter length. Along these edges, which in total exceed 2·5 m g?1 (equivalent to 500 m m?2 for a loading of 2 t ha?1), surface tension draws up water from between the litter particles. This reduces flow depth there, and as a consequence of the lower flow intensity, frictional drag rises. Furthermore, no clear passage remains for the establishment of flow threads. These findings apply to shallow interrill flows in which litter is largely immobile. The key new result from these experiments is that under these conditions, a 20 per cent cover of organic litter can generate interrill frictional retardation that exceeds by nearly 41 per cent that of a bare soil surface, and twice that contributed by the same cover fraction of surface stones. Even greater dominance by litter can be anticipated at the many dryland sites where litter covers exceed those tested here. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

10.
Rates of splash detachment from a sandy soil of the Cottenham Series, subjected to a five-minute design storm of 50 mm/h with a kinetic energy of 127 J/m2 and a median volume drop size of 3·2 mm supplied from a rotating-disc rainfall simulator, are determined without a plant cover and with the cover of a single Brussels sprouts plant. Measurements are made at regular intervals throughout the growing season. Plant canopies of 10 to 25 per cent result in reductions of 10 to 25 per cent in rainfall volume and 10 to 81 per cent in rainfall energy. The volume and energy of the rain beneath the plant are significantly (P < 0·05) correlated with its number of leaves (r = ?0·84 and ?0·92 respectively for n = 49). No reduction was observed in the rate of splash detachment which averages 1·2 kg/m2 for the design storm with and without the plant. The detachment rate was found rather surprisingly to be inversely related to the energy of the rainfall under the plant and positively related to the number of leaves. As the number of leaves increases, so does the detachment rate per unit of rainfall energy.  相似文献   

11.
Some 300 measurements under a wide range of flow conditions were used to construct regression models to predict concentrations of six major ions in a moorland stream draining intermediate-basic lavas in Scotland. Ca, Mg, and HCO3 were very closely related to discharge in log-log regressions (r2 > 90 percent). Na and Cl required inclusion of seasonal functions for best prediction (r2 = 76 per cent in both cases). K concentrations were predicted from discharge, direction of stage change, and a seasonal function. Annual losses of Ca, Mg, and HCO3 were 252, 36, and 871 kg ha ?1, substantially larger than those for upland catchments on acid rocks. Losses of K, Na, and Cl were of a similar order to other upland sites.  相似文献   

12.
Lee-side windspeed and sediment transport were measured over a small (1·2 m) transverse ridge in the Silver Peak dunefield, west-central Nevada, USA, using an intensive array of 25 cup anemometers and seven total flux traps. During crest-transverse and transporting flow conditions (u0·3crest ≈ 8·4 m s−1), windspeed near the surface of the lee slope averaged half (48 per cent) that of crest speeds. Dimensionless speeds in the separation zone ranged from 0·2 to 0·8 that of the outer flow (u12). Along the boundary of the separation cell, windspeed increased by 10 per cent of the crest speed before separation. Equilibrium of upper and lower wake regions was not observed by the documented eight dune heights, suggesting that wake recovery may not occur over closely spaced dunes. Sediment transport measured directly on both the lee slope and interdune surfaces averaged approximately 15 per cent of crest inputs. This suggests that a significant amount (c. 70–95 per cent) of sediment transported over the crest moved as fallout. For this data set, flux was approximately proportional to the cube of the near-surface windspeed (u0·3) and in general there was an order of magnitude difference between flux measured at the crest and that measured within the separation zone. Transport direction in the separation zone was acutely oblique to the incident direction owing to secondary flow deflection. Beyond the interdune, transport direction progressed from oblique to crest-transverse. This indicates that an appreciable amount of sediment may move laterally along the lee slope and interdune corridor under crest-transverse flows. Regarding the grain size and sorting properties of transported sediment, there was no significant difference in mean grain size over the dune, although in general particles were finer and more poorly sorted in the lee. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

13.
The effect of turbulent flow structures on saltation sand transport was studied during two convective storms in Niger, West Africa. Continuous, synchronous measurements of saltation fluxes and turbulent velocity fluctuations were made with a sampling frequency of 1 Hz. The shear stress production was determined from the vertical and streamwise velocity fluctuations. The greatest stress-bearing events were classified as turbulent structures, with sweep, ejection, inward interaction, and outward interaction described according to the quadrant technique. The classified turbulent structures accounted for 63·5 per cent of the average shear stress during the first storm, and 56·0 per cent during the second storm. The percentage of active time was only 20·6 per cent and 15·8 per cent, respectively. High saltation fluxes were associated with sweeps and outward interactions. These two structures contribute positively (sweeps) and negatively (outward interactions) to the shear stress, but have in common that the streamwise velocity component is higher than average. Therefore, the horizontal drag force seems primarily responsible for saltation sand transport, and not the shear stress. This was also reflected by the low correlation coefficients (r) between shear stress and saltation flux (0·12 and 0·14, respectively), while the correlation coefficients between the streamwise velocity component and saltation flux were much higher (0·65 and 0·57, respectively). © 1998 John Wiley & Sons, Ltd.  相似文献   

14.
Sediment transport capacity, Tc, defined as the maximum amount of sediment that a flow can carry, is the basic concept in determining detachment and deposition processes in current process-based erosion models. Although defined conceptually and used extensively in modelling erosion, Tc was rarely measured. Recently, a series of laboratory studies designed to quantify effects of surface hydrologic conditions on erosion processes produced data sets feasible to evaluate the concept of Tc. A dual-box system, consisting of 1·8 m long sediment feeder box and a 5 m long test box, was used. Depending on the relative magnitudes of sediment delivery from feeder and test boxes, five scenarios are proposed ranging from deposition-dominated to transport-dominated sediment regimes. Results showed that at 5 per cent slope under seepage or 10 per cent slope under drainage conditions, the runoff from the feeder box caused in the additional sediment transport in the test box, indicating a transport-dominated sediment regime. At 5 per cent slope under drainage conditions, deposition occurred at low rainfall intensities. Increases in slope steepness, rainfall intensity and soil erodibility shifted the dominant erosion process from deposition to transport. Erosion process concepts from the Meyer–Wishmeier, Foster–Meyer and Rose models were compared with the experimental data, and the Rose model was found to best describe processes occurring during rain. A process-based erosion model needs to have components that can represent surface conditions and physical processes and their dynamic interactions. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

15.
Four techniques for soil erosion assessment were compared over two consecutive seasons for bare-fallow plots and a maize-cowpea sequence in 1985 at IITA, Ibadan, Nigeria. The techniques used were: tracer (aluminium paint), nails (16 and 25), the rill method, and the Universal Soil Loss Equation (USLE). Soil loss estimated by these techniques was compared with that determined using the runoff plot technique. There was significantly more soil loss (P < 0·01) in bare-fallow than in plots under maize (Zea mays) or cowpea (Vigna unguiculata). In the first season, soil loss from plots sown to maize was 40·2 Mg ha?1 compared with 153·3 Mg ha?1 from bare-fallow plots. In the second season, bare-fallow plots lost 87·5 Mg ha?1 against 39·4 Mg ha?1 lost from plots growing cowpea. The techniques used for assessing erosion had no influence on the magnitude of soil erosion and did not interfere with the processes of erosion. There was no significant difference (P < 0·05) between soil erosion determined by the nails and the runoff plot technique. Soil loss determined on six plots (three under maize, three bare-fallow) by the rill technique, at the end of the season, was significantly lower (P < 0·05) than that determined by the runoff plot technique. The soil loss estimated by the rill method was 143·2, 108·8 and 121·9 Mg ha?1 for 11, 11, and 8 per cent slopes respectively, in comparison with 201·5, 162·0, and 166·4 Mg ha?1 measured by the runoff plot method. Soil loss measured on three bare-fallow plots on 10 different dates by the rill technique was also significantly lower (P < 0·01) than that measured by the runoff plot. In the first season the USLE significantly underestimated soil loss. On 11, 11, and 8 per cent slopes, respectively, soil loss determined by the USLE was 77, 92, and 63 per cent of that measured by the runoff plot. However, in the second season there was no significant difference between soil loss determined by the USLE and that determined by the conventional runoff plot technique.  相似文献   

16.
The wave energy loss as waves break on nearshore bars is an important geomorphological process as it controls the stability of the adjacent beaches and dunes. The effect is particularly important during storms as it serves to guard against severe shoreline erosion. Field measurements from Florida and Ireland of waves crossing both single and multiple bar bedforms, indicate that where bar-breaking occurs, between 78 and 99 per cent of wave energy may be dissipated from individual waves. However, a further 20 per cent energy may reach the shore if wave frequencies increase during the reformation process. This latter effect was not noted at the multiple bar site where spilling mode breakers dominate. The level of energy dissipation is related to the time taken to cross that section of the submerged bar below the critical depth, d < 1.28H. In the field this may be approximated by the ratio of wave speed to the width of the bar breaker zone.  相似文献   

17.
This paper presents an evaluation of the feasibility and the reliability of a visual characterization technique for gravel–cobble river bed surface substrate. Based on principal axis regressions, using phi scale (ϕ), comparisons of visual estimation and grid sampling techniques show that useful predictive relations (R2 = 0·78–0·88) exist between visual estimates of the surface d16, d50 and d84 and estimates obtained for the same percentiles with the grid sampling technique. Comparisons of visual estimation and the surface‐bulk sampling technique also indicate a predictive relation (R2 = 0·70) between the d50 of the two methods. Trained operators can visually estimate gravel–cobble bed surface d16 to uncertainties of 41 per cent, d50 to 15 per cent and d84 to 11 per cent (for example, there is a 5·5 mm error on a d84 size of 50 mm). Furthermore, evidence shows that if operators are properly trained, a calibration relation for each percentile can be applied independently of operators. This visual characterization allows effective detailed mapping of spatial patterns in substrate size distribution along extensive reaches of gravel‐bed rivers. The technique can be very useful in creating terrain models for various geomorphological, hydrological and biological applications such as the determination of entrainment thresholds, hydraulic roughness and substrate suitability for benthic insects or salmonid habitat. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

18.
ABSTRACT

The evaporation losses from the container of an old galvanized Hellmann gauge, 7·1 cm2 aperture area, were five times greater than the losses from the container of a new gauge, 1·8 cm2 aperture area, of almost the same grey colour. The maximum evaporation from the old gauge amounted to 0·75 mm per day. The same evaporation losses in the Baye of Montreux basin over a period of 10 years from April to September amounted to 0·09 mm per measurement or to 0·7 per cent of the precipitation catch. There is a relationship between the monthly percentage evaporation losses and the ratio of evaporation time and rainfall duration (r = 0·803).  相似文献   

19.
Bagnold's theory of bedload stresses and bedforms requires that the tangential stress exerted at the base of a moving bedload dispersion in high stage plane beds is carried entirely by the grain stress T. Derivation of T from experimental data in flow channnels shows, however, that the fluid contribution τ′ is still dominant (≈ 1 to 2 × T). Particle concentrations in the experimental bedload dispersions approached 10 per cent and it is inferred that no appreciable concentration gradient exists. Grain collisions are thought to be most probable whenever more than 1 former stationary grain layer is in forward motion. High stage plane beds may be unstable when grain collisions cease. Current ripples evolve from the amplification and migration of chance defects due to unsteady fluid impulses generated by injection, or burst, phenomena in the viscous sublayer and not through any inherent static shear resistance deficit exerted by stationary bed grains.  相似文献   

20.
Fragments of weathered granitic rocks from the Kunlun Shan, Qinghai Plateau (China) were investigated to elucidate the in?uence of biotic crusts on the breakdown of granitic rocks in an alpine environment. Scanning electron microscopy with energy dispersive system and X‐ray diffractometry were used to describe the nature and properties of mineral accumulations on the rock surface. Results showed that organic salts such as calcium oxalate and calcium formate are associated with Aspicilia caesiocinera (Nyl.ex Malbr.) Arnold, Caloplaca sp., Xanthoria elegans (Link) Th.Fr., and Lecidea plana (Lahm) Nyl. Secondary accumulations of 2 : 1 clays minerals are found in A. caesiocinera while oxides of manganese are associated with X. elegans. Coatings of goethite (iron oxides) are believed to form from biological activity associated with the presence of hyphae and rodlet structures on the ?akes. Calcium oxalate crystallizes into several morphologies such as druse, hexagonal plates, and lenticular containing between 20 and 48 per cent calcium by weight. Calcium formate and iron oxide (goethite) occur together in the form of ‘red’ desert varnish. Observed ‘black’ coatings contain as much as 37 per cent manganese and 22 per cent iron. Clay accumulations have plate‐like morphology and contain c. 2 : 1 silicon to aluminium contents. We argue that organic acids from the activities of biotic crusts contribute to the breakdown of granitic rocks. Fungi accelerate the breakdown of granitic rocks through the growth of fungal hyphae along the 001 cleavage planes in primary chloritic minerals. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

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

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