首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Effects of large organic material on channel form and fluvial processes   总被引:1,自引:0,他引:1  
Stream channel development in forested areas is profoundly influenced by large organic debris (logs, limbs and rootwads greater than 10 cm in diameter) in the channels. In low gradient meandering streams large organic debris enters the channel through bank erosion, mass wasting, blowdown, and collapse of trees due to ice loading. In small streams large organic debris may locally influence channel morphology and sediment transport processes because the stream may not have the competency to redistribute the debris. In larger streams flowing water may move large organic debris, concentrating it into distinct accumulations (debris jams). Organic debris may greatly affect channel form and process by: increasing or decreasing stability of stream banks; influencing development of midchannel bars and short braided reaches; and facilitating, with other favourable circumstances, development of meander cutoffs. In steep gradient mountain streams organic debris may enter the channel by all the processes mentioned for low gradient streams. In addition, considerable debris may also enter the channel by way of debris avalanches or debris torrents. In small to intermediate size mountain streams with steep valley walls and little or no floodplain or flat valley floor, the effects of large organic debris on the fluvial processes and channel form may be very significant. Debris jams may locally accelerate or retard channel bed and bank erosion and/or deposition; create sites for significant sediment storage; and produce a stepped channel profile, herein referred to as ‘organic stepping’, which provides for variable channel morphology and flow conditions. The effect of live or dead trees anchored by rootwads into the stream bank may not only greatly retard bank erosion but also influence channel width and the development of small scour holes along the channel beneath tree roots. Once trees fall into the stream, their influence on the channel form and process may be quite different than when they were defending the banks, and, depending on the size of the debris, size of the stream, and many other factors, their effects range from insignificant to very important.  相似文献   

2.
Streambed hydraulic conductivity is one of the main factors controlling variability in surface water‐groundwater interactions, but only few studies aim at quantifying its spatial and temporal variability in different stream morphologies. Streambed horizontal hydraulic conductivities (Kh) were therefore determined from in‐stream slug tests, vertical hydraulic conductivities (Kv) were calculated with in‐stream permeameter tests and hydraulic heads were measured to obtain vertical head gradients at eight transects, each comprising five test locations, in a groundwater‐dominated stream. Seasonal small‐scale measurements were taken in December 2011 and August 2012, both in a straight stream channel with homogeneous elevation and downstream of a channel meander with heterogeneous elevation. All streambed attributes showed large spatial variability. Kh values were the highest at the depositional inner bend of the stream, whereas high Kv values were observed at the erosional outer bend and near the middle of the channel. Calculated Kv values were related to the thickness of the organic streambed sediment layer and also showed higher temporal variability than Kh because of sedimentation and scouring processes affecting the upper layers of the streambed. Test locations at the channel bend showed a more heterogeneous distribution of streambed properties than test locations in the straight channel, whereas within the channel bend, higher spatial variability in streambed attributes was observed across the stream than along the stream channel. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
Field observations on hydraulics and sediment dynamics during extreme floods in two mountain torrents show the influence of man-made constructions such as bridges and check dams, in addition to the sediment supplied naturally by the basin and the channel network, on the formation of hyperconcentrated flows. In the Pyrenean Arás basin, hyperconcentrated flow occurred after collapse of a bridge, which in turn mobilized large volumes of sediment from the stream channel and, subsequently, destroyed a series of check dams. Boulders up to several metres in size were transported in a mixture of sand and fine material. A minimum of 100000 tonnes of sediment were deposited on the alluvial fan during the event. Prior to bridge destruction, mean bedload transport rates had reached 0.4t m−1 s−1 upstream. In the alpine Lainbach basin, the flood was characterized by transportation of large amounts of slope material, including debris flows. Along its main tributary an intensive hyperconcentrated flow occurred during the rising stage, whereas in the main valley smaller flows occurred after failure of check dams. The depth of coarse material deposited reached 80 cm. The effectiveness of the Aràs and Lainbach floods was attained due to exceptional rates of energy expediture. Flood power ranged from 20000 W m−2 to 40000 W m−2 on average. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

4.
Streambed hydraulic conductivity (K) and vertical K (K v) are key controls on groundwater and surface water exchange and biogeochemical fluxes through the hyporheic zone, but drivers of transient hydraulic properties in different fluvial environments are poorly understood. This study combines hydrogeology, geophysics, and sedimentology to reveal mechanisms of K and K v transience in the upper 0.5 m of a sandy streambed during low discharge. Hydraulic tests (44 slug tests, 130 falling-head permeameter tests) and 130 grain-size analyses were repeated three times over 8 weeks on a 1,200 m2 grid spanning: (a) a channel with continuously flowing water and mobile bed load; (b) an adjacent mid-channel bar that was stationary and infrequently submerged. Aerial photographs and ground-penetrating radar show scour and complete reworking of fluvial sediments in the channel. Bar sediments below the water table remained immobile, but infrequent flows of moderate discharge reworked the uppermost few centimetres of the bar top. Despite differences in sediment mobility and stream flow characteristics across environments, K and K v exhibited order-of-magnitude differences in spatial heterogeneity and temporal variability in both the channel and bar. Mean K and K v values in the channel were comparatively stable over time. In the immobile bar, mean K declined 20% and K v declined 26% after increased discharge temporarily inundated the bar. Grain-size distributions were steady across both environments over time, but repeat geophysical surveys of the bar show a decrease in electrical conductivity, likely from porosity reduction. These findings suggest that sediment dynamics and stream flow characteristics in different streambed environments are important drivers of K transience during low discharge conditions. Specifically, pore clogging can be an important mechanism of transience over short durations (weeks to months) in immobile sediments subject to infrequent flows and minor reworking.  相似文献   

5.
Abstract

The term debris torrent was introduced in North America's Pacific Northwest region to describe coarse-grained, rapid channelized flows rich in organic debris. Environmental conditions favouring their occurrence and distinctive criteria for recognition of debris torrents are discussed. A classification of gravitational and fluvial sediment transporting processes is presented as a framework for papers presented at the Debris Torrent Workshop 6a held in Vancouver in August, 1987.  相似文献   

6.
Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ0) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time tx can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.  相似文献   

7.
Dividing rivers into homogeneous reaches is key for river processes and watershed management. In contrast to downstream fluvially dominated rivers, upstream debris-flow dominated torrents have steeper channel slopes and smaller valley width/depth ratios. Investigating transition reaches between torrents and fluvially dominated rivers, not only explores the structure of the landscape, but also contributes to hazard management. This study proposed a valley morphology index combining two variables, channel slope and valley width/depth ratio, to determine transition reaches between torrents and rivers. The methodology was applied to 41 mountain streams in Taiwan using a Geographic Information System (GIS)-based topographic analysis. Plots of valley width/depth ratio versus channel slope were used to determine boundary values of the valley morphology index (Iv) separating torrents from rivers. The plots showed that about 80% of the river basins present “L-shaped” curves, which indicate sharp decreases in slope for upstream sections and dramatic increases of valley width/depth ratio for downstream sections. Results further demonstrated an average value of Iv 0.0047 across the study sites. Spatial comparison between geographic regions indicated that transition reaches in eastern rivers tend to occur lower in the drainage basin due, in part, to higher terrain. Local factors, such as tributary confluences and landslides promote the transition from torrents to fluvially dominated rivers. Satellite images verified that the approach correctly identified transition reaches, suggesting that it may provide a useful reference for river management.  相似文献   

8.
The concept of stream channel grade – according to which a stream channel reach will adjust its gradient, S, in order to transport the imposed sediment load having magnitude Qb and characteristic grain size Db, with the available discharge Q (Mackin, 1948 , Geological Society of America Bulletin 59 : 463–512; Lane, 1955 , American Society of Civil Engineers, Proceedings 81 : 1–17) is one of the most influential ideas in fluvial geomorphology. Herein, we derive a scaling relation that describes how externally imposed changes in either Qb or Q can be accommodated by changes in the channel configuration, described by the energy gradient, mean flow depth, characteristic grain size and a parameter describing the effect of bed surface structures on grain entrainment. One version of this scaling relation is based on the dimensionless bed material transport parameter (W*) presented by Parker and Klingeman ( 1982 , Water Resources Research 18 : 1409–1423). An equivalent version is based on a new dimensionless transport parameter (E*) using dimensionless unit stream power. This version is nearly identical to the relation based on W*, except that it is independent of flow resistance. Both versions of the scaling relation are directly comparable to Lane's original relation. In order to generate this stream power‐based scaling relation, we derived an empirical transport function relation relating E* to dimensionless stream power using data from a wide range of stable, bed load‐dominated channels: the form of that transport function is based on the understanding that, while grain entrainment is related to the forces acting on the bed (described by dimensionless shear stress), sediment transport rate is related to the transfer of momentum from the fluid to the bed material (described by dimensionless stream power). Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
Two landsliding episodes between late 1973 and early 1975 delivered about 60000 m3 of sediment to six small deeply incised streams draining a 2·7 km2 area. About 4700 m3 of logs in the landslide debris formed major log jams in five streams, which impounded large volumes of landslide-derived sediment. Five years after the landsliding, 42 per cent (25000 m3) of sediment was still in storage behind 35 log jams ranging from 1·4–8·2 m high. The landsliding episodes have produced multi-stepped stream profiles, aggradation of channel reaches up to 150 m long to mean depths between 1·2 and 4·1 m, reductions in gradient, fining of bed material size, and related changes in bedforms and channel width:depth ratios that seem likely to persist for at least several decades. Sediment presently stored behind log jams is equivalent to between 50 and 220 years normal supply of sediment from hillslopes to stream channels. Long-delayed, large magnitude impacts on higher-order channels may occur if sudden failure of log jams is induced by a large storm at some future date.  相似文献   

10.
贺拿  杨建元  陈宁生  朱云华 《地震学刊》2013,(6):671-678,697
通过对四川凉山矮子沟流域的勘查,结合遥感解译及矮子沟流域的地质活动历史,分析“6·28”矮子沟泥石流的影响因素,在此基础上分析了矮子沟泥石流的成因。“6·28”矮子沟泥石流是在地震、人类活动及干旱影响下,受强降雨激发,坡面及支沟先发生泥石流,支沟泥石流堵塞主沟泄洪通道,强降雨使堵塞体瞬间溃决形成较大山洪,山洪的强大揭底作用使沟床松散固体物质启动,沿途不断有岸坡坍塌补给,中下游段有大块石加入,规模不断增大而形成。采用形态调查法计算的泥石流流量为925.91m3/s,泥石流总量为57.4×104m3,一次固体物质总量为27×104m3。通过对矮子沟流域的物源、地形及流域的气候特征研究,得出初步结论:近期泥石流的规模和频率都会降低,但在极端干旱、地震和强降雨的偶然联合作用下,依然可能发生大规模的泥石流灾害。  相似文献   

11.
Channel adjustments in the North Fork Toutle River and the Toutle River main stem were initiated by deposition of a 2.5 km3 debris avalanche and associated lahars that accompanied the catastrophic eruption of Mount St. Helens, Washington on 18 May 1980. Channel widening was the dominant process. In combination, adjustments caused average boundary shear stress to decrease non-linearly with time and critical shear stress to increase non-linearly with time. At the discharge that is equalled or exceeded 1 per cent of the time, these trends converged by 1991–1992 so that excess shear stress approached minimum values. Extremal hypotheses, such as minimization of unit stream power and minimization of the rate of energy dissipation (minimum stream power), are shown to be applicable to dynamic adjustments of the Toutle River system. Maximization of the Darcy–Weisbach friction factor did not occur, but increases in relative bed roughness, caused by the concomitant reduction in hydraulic depths and bed-material coarsening, were documented. Predictions of stable channel geometries using the minimum stream power approach were unsuccessful when compared to the 1991–1992 geometries and bed-material characteristics measured in the field. It is concluded that the predictions are not applicable because the study reaches are not truly stable and cannot become so until a new floodplain has been formed by renewed channel incision, retreat of stream-side hummocks, and establishment of riparian vegetation to limit the destabilizing effects of large floods. Further, prediction of energy slope (and consequently stream power) by the sediment transport equations is inaccurate because of the inability of the equations to account for significant contributions of finer grained (sand and gravel) bank materials (relative to the coarsened channel bed) from bank retreat and from upstream terrace erosion.  相似文献   

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

13.
Morphometric variables associated with 41 debris‐?ow and 18 ?uvial fans and their basins in the Southern Alps of New Zealand are examined. The results show statistically signi?cant differences in the area, maximum elevation, relief and ruggedness (Melton's R) of the basin and the area, gradient, and apex and toe elevations of the fan between debris‐?ow and ?uvial sites. Concavity of the fan longitudinal pro?le also differs between the two fan types, although this could not be tested statistically. Most of these morphometric differences re?ect differences in processes and environmental controls on them. Discriminant analysis indicates that basin area and fan gradient best differentiate the two fan types by process. Moderately strong correlations exist, on both debris‐?ow and ?uvial fans, between basin area or Melton's R and fan area. Correlations between basin area or Melton's R and fan gradient are generally weaker. The results of this study also indicate that on debris‐?ow‐prone fans the fan gradient and basin Melton's R have lower thresholds which overlap little with upper thresholds associated with basins where only stream?ow reaches the fan. These thresholds may therefore have value in preliminary morphometric assessments of debris‐?ow hazard on fans in the Southern Alps. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

14.
Spatial analysis of stream power using GIS: SLk anomaly maps   总被引:2,自引:0,他引:2  
The stream length‐gradient index (SL) shows the variation in stream power along river reaches. This index is very sensitive to changes in channel slope, thus allowing the evaluation of recent tectonic activity and/or rock resistance. Nevertheless, the comparison of SL values from rivers of different length is biased due to the manner in which the index is formulated, thus making correlations of SL anomalies along different rivers difficult. Therefore, when undertaking a comparison of SL values of rivers of different lengths, a normalization factor must be used. The graded river gradient (K) has already been used in some studies to normalize the SL index. In this work, we explore the relationships between the graded river gradient (K), the SL index and the stream power, proposing the use of a re‐named SLk index, which enables the comparison of variable‐length rivers, as well as the drawing of SLk anomaly maps. We present here a GIS‐based procedure to generate SLk maps and to identify SLk anomalies. In order to verify the advantages of this methodology, we compared an SLk map of the NE border of the Granada basin with both simple river profile–knickpoint identification and with an SL map. The results show that the SLk map supplies good results with defined anomalies and suitably reflects the main tectonic and lithological features of the study area. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

15.
16.
The impact of afforestation on stream bank erosion and channel form   总被引:1,自引:0,他引:1  
Modification of the land use of a small catchment through coniferous afforestation is shown to have influenced stream bank erosion and channel form. Field mapping and erosion pin measurements over a 19-month period provides evidence of more active bank erosion along forested channel reaches than along non-forested. Extrapolation of downstream increases in bankfull width, bankfull depth, and channel capacity with increasing basin area for the non-forested catchment has demonstrated that afforestation of the lower part of the catchment has had a marked effect on channel form. Channel widths within the forest are up to three times greater than that predicted from the regression. These changes in bankfull width have led to stream bed aggradation and the development of wide shallow channels within the forest, and channel capacities within the forest are over two times that predicted from the basin area. The relationship between channel sinuosity and valley gradient for non-forested reaches of the river also indicated decreased sinuosity resulting from afforestation. These changes in channel form result from active bank erosion within the forest with coarse material being deposited within the channel as point-bars and mid-channel bars. Active bank erosion is largely attributed to the suppression by the forest of a thick grass turf and its associated dense network of fine roots, and secondly to the river attempting to bypass log jams and debris dams in the stream channel.  相似文献   

17.
Rainfall was simulated on unconfined plots on regolith in debris flow source areas using a portable simulator. In total, 351 simulations were carried out on steep slopes (27–54°) with rainfall intensities of 28–291 mm/h. From these rainfall simulations the infiltration parameters sorptivity (S) and steady-state infiltration capacity (K) of the regolith, and a threshold for the occurrence of micro-scale mass movements, were obtained. Two evaluation methods were used to obtain the infiltration parameters K and S. The ‘infiltration envelope’ method uses rainfall intensity and time to ponding from multiple tests and fits an infiltration envelope through the data from which K and S can be obtained. The ‘constant runoff’ method uses rainfall intensity and overland flow intensity to calculate K, after which S can be calculated in several ways by using time to ponding. The constant runoff method produced K values of 16.6–128 mm/h, which usually show a log-normal distribution. K values depend on the regolith parent material and rainfall intensity. Using this method, S values are 0.088–0.381 cm/min1/2. The infiltration envelope method produced K values of 9.8–131 mm/h and S values of 0.14–0.32 cm/min1/2. It can be argued that both methods overestimate K as well as S, but quantitative relations between measured/calculated and actual values of K and S have not yet been obtained. At high rainfall intensities, typically 100 mm/h or more, micro-scale mass movements sometimes occur. A lower threshold curve for the occurrence of these micro-scale mass movements has been constructed. It is a function of both slope angle and rainfall intensity. The micro-scale mass movements could play an important part in the initiation of debris flows in the study area, possibly by delivering sediment to overland flow. On the very steep slopes, the sediment-rich overland flow can easily mobilize coarse material.  相似文献   

18.
Abstract

Piezometers and wells installed for water quality monitoring are frequently used to assess the saturated hydraulic conductivity (K) in the surrounding formation. A series of recovery tests was conducted to evaluate how purging, required to obtain representative water quality samples, affected measured values of hydraulic conductivity in 15 newly installed and undeveloped piezometers placed to between 2 and 15 m depth (in oxidized and unoxidized material) in a loamy glacial till (K range from 10?6 to 10?9 m s?1). Piezometers were purged between 9 and 11 times for sampling over a period of five months. The effect of the purgings on piezometer development was evaluated by changes in slope of the water level recovery curves which were used to calculate hydraulic conductivity. The first five purgings following piezometer installation increased K in the 15 piezometers by an average of 34%. The average increase in a value of K after 10 purgings was 44%. Values measured for hydraulic conductivity in a 75 mm diameter auger hole appeared stable after four purgings but piezometers installed in larger diameter boreholes (100 mm to 280 mm) snowed increases in K with up to 10 purgings. The hydraulic conductivity determined for piezometers installed at a 30° angle to the vertical showed greater variability than was observed in the adjacent vertically installed piezometers at the same depth.  相似文献   

19.
Abstract

Dissolved oxygen (DO) is one of the most useful indices of river's health and the stream re-aeration coefficient is an important input to computations related to DO. Normally, this coefficient is expressed as a function of several variables, such as mean stream velocity, shear stress velocity, bed slope, flow depth, and Froude number. However, in free surface flows, some of these variables are interrelated, and it is possible to obtain simplified stream re-aeration equations. In recent years, different functional forms have been advanced to represent the re-aeration coefficient for different data sets. In the present study, the artificial neural network (ANN) technique has been applied to estimate the re-aeration coefficient (K 2) using data sets measured at different reaches of the Kali River in India and values obtained from the literature. Observed stream/channel velocity, bed slope, flow depth, cross-sectional area and re-aeration coefficient data were used for the analysis. Different combinations of variables were tested to obtain the re-aeration coefficient using an ANN. The performance of the ANN was compared with other estimation methods. It was found that the re-aeration coefficient estimated by using an ANN was much closer to the observed values as compared with the other techniques.  相似文献   

20.
Xunhong Chen 《水文研究》2011,25(2):278-287
Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (Kv) of the channel sediments and the aquifer materials. This method was applied to produce Kv profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the Kv profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of Kv in the channel sediments. This depth‐dependent pattern of Kv distribution for the channel sediments is considered to be produced by hyporheic processes. This Kv‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

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

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