首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A rising exposure to flood risk is a predicted consequence of increased development in vulnerable areas and an increase in the frequency of extreme weather events due to climate change. In the face of this challenge, a continued reliance on engineered at‐a‐point flood defences is seen as both unrealistic and undesirable. The contribution of ‘soft engineering’ solutions (e.g. riparian forests, wood in rivers) to integrated, catchment scale flood risk management has been demonstrated at small scales but not larger ones. In this study we use reduced complexity hydrological modelling to analyse the effects of land use and channel changes resulting from river restoration upon flood flows at the catchment scale. Results show short sections of river‐floodplain restoration using engineered logjams, typical of many current restoration schemes, have highly variable impacts on catchment‐scale flood peak magnitude and so need to be used with caution as a flood management solution. Forested floodplains have a more general impact upon flood hydrology, with areas in the middle and upper catchment tending to show reductions in peak magnitude at the catchment outflow. The most promising restoration scenarios for flood risk management are for riparian forest restoration at the sub‐catchment scale, representing 20–40% of the total catchment area, where reductions in peak magnitude of up to 19% are observed through de‐synchronization of the timings of sub‐catchment flood waves. Sub‐catchment floodplain forest restoration over 10–15% of total catchment area can lead to reductions in peak magnitude of 6% at 25 years post‐restoration. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
Regulated rivers generally incise below dams that cut off sediment supply, but how that happens and what the consequences are at different spatial scales is poorly understood. Modern topographic mapping at meter‐scale resolution now enables investigation of the details of spatial processes. In this study, spatial segregation was applied to a meter‐scale raster map of topographic change from 1999 to 2008 on the gravel‐cobble, regulated lower Yuba River in California to answer specific scientific questions about how a decadal hydrograph that included a flood peak of 22 times bankfull discharge affected the river at segment, reach, and morphological unit scales. The results show that the river preferentially eroded sediment from floodplains compared to the channel, and this not only promoted valley‐wide sediment evacuation, but also facilitated the renewal and differentiation of morphological units, especially in the channel. At the reach scale, area of fill and mean net rate of elevational change were directly correlated with better connectivity between the channel and floodplain, while the mean rate of scour in scour areas was influenced by the ratio of slope to bankfull Froude number, a ratio indicative of lateral migration versus vertical downcutting. Hierarchical segregation of topographic change rasters proved useful for understanding multi‐scalar geomorphic dynamics. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

3.
  相似文献   

4.
Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

5.
Most analyses of river adjustment have focused on parts of catchments where metamorphosis has occurred. This provides a non‐representative view of river responses to human‐disturbance. Although many rivers have been subjected to systematic land‐use change and disturbance, significant variability is evident in the form, extent and consequences of adjustment. This study documents the catchment‐wide distribution of river sensitivity and adjustment in the upper Hunter catchment, New South Wales, Australia in the period since European settlement. The spatial distribution and timing of lateral, vertical and wholesale river adjustments are used to assess river sensitivity to change. The type and pattern of rivers, influenced largely by valley setting, have induced a fragmented pattern of river adjustment in the upper Hunter catchment. Adjustments have been largely non‐uniform and localized, reflecting the predominance of bedrock‐controlled rivers which have limited capacity to adjust and are resilient to change. Less than 20% of river courses have experienced metamorphosis. Phases of reach‐scale geomorphic adjustment to human disturbance are characterized as a gradient of primary, secondary and tertiary responses. In general terms, primary responses such as cutoffs or straightening were followed by secondary responses such as channel expansion. These secondary responses occurred between 50–70 years after initial disturbance. A subsequent tertiary phase of river recovery, denoted as a transition from predominantly erosional to predominantly depositional geomorphic processes such as channel contraction, occurred around 70–120 years after initial disturbance. Such responses are ongoing across much of the upper Hunter catchment. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

6.
The current benchmark approach for mathematical modelling of floodplain hydrologic regime consists of dynamically coupling one‐dimensional (1D) and two‐dimensional (2D) models for flow routing along the main channel and the floodplain, respectively. For large‐scale sites, floodplain inundation may spread over hundreds of square kilometres and may last for many months and even influence seasonal floods in following years. This paper aims at evaluating the effect of vertical water balance representation on modelling a large‐scale floodplain. The Pantanal wetland (140 000 km2; Brazil) is simulated using a 1D/2D coupled model approach, which also considers the representation of vertical water processes over the floodplain. Four scenarios are simulated: Baseline (the reference scenario), NoVertBal (in which the vertical water balance over floodplain is turned off) and ETp+1 and ETp?1 scenarios, characterized by artificially increasing or decreasing daily potential evapotranspiration (ETp) by 1 mm, respectively. The results showed that the effect of the vertical water processes scenarios on channel flow is directly dependent on the lateral exchange of water between the channel and floodplain in the upstream river reach. This influence is stronger when there is a gain of water from the floodplain to the channel. The inclusion of these vertical water processes into floodplain modelling was essential to represent the process of wetting and drying, this effect being more relevant for areas not directly connected to main channels, which is a characteristic of the Pantanal region. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
Bankfull discharge is a key parameter in the context of river engineering and geomorphology, as an indicator of flood discharge capacity in alluvial rivers, and varying in response to the incoming flow and sediment regimes. Bankfull channel dimensions have significantly adjusted along the Lower Yellow River (LYR) due to recent channel degradation, caused by the operation of the Xiaolangdi Reservoir, which has led to longitudinal variability in cross‐sectional bankfull discharges. Therefore, it is more representative to describe the flood discharge capacity of the LYR, using the concept of reach‐averaged bankfull discharge. Previous simple mean methods to estimate reach‐scale bankfull discharge cannot meet the condition of flow continuity or account for the effect of different spacing between two sections. In this study, a general method to calculate cross‐sectional bankfull discharge using the simulated stage‐discharge relation is outlined briefly, and an integrated method is then proposed for estimating reach‐scale bankfull discharge. The proposed method integrates a geometric mean based on the log‐transformation with a weighted average based on the spacing between two consecutive sections, which avoids the shortcomings of previous methods. The post‐flood reach‐scale bankfull discharges in three different channel‐pattern reaches of the LYR were estimated annually during the period from 1999 to 2010 using the proposed method, based on surveyed post‐flood profiles at 91 sedimentation sections and the measured hydrological data at seven hydrometric sections. The calculated results indicate that: (i) the estimated reach‐scale bankfull discharges can effectively represent the flood discharge capacity of different reaches, with their ranges of variation being less than those of typical cross‐sectional bankfull discharges; and (ii) the magnitude of the reach‐scale bankfull discharge in each reach can respond well to the accumulative effect of incoming flow and sediment conditions. Finally, empirical relationships for different reaches in the LYR were developed between the reach‐scale bankfull discharge and the previous four‐year average discharge and incoming sediment coefficient during flood seasons, with relatively high correlation coefficients between them being obtained, and the reach‐scale bankfull discharges in different reaches predicted by the delayed response model were also presented for a comparison. These relations for the prediction of reach‐scale bankfull discharges were validated using the cross‐sectional profiles and hydrological data measured in 2011. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
Digital elevation models (DEMs) of river channel bathymetries are developed by interpolating elevations between data collected at discrete points or along transects. The accuracy of interpolated bathymetries depends on measurement error, the density and distribution of point data, and the interpolation method. Whereas point measurement errors can be minimized by selecting the most efficient equipment, the effect of data density and interpolation method on river bathymetry is relatively unknown. Thus, this study focuses on transect‐based collection methods and investigates the effects of transect location, the spacing between transects, and interpolation methods on the accuracy of interpolated bathymetry. This is accomplished by comparing four control bathymetries generated from accurate and high resolution, sub‐meter scale data to bathymetries interpolated from transect data extracted from the control bathymetries using two transect locating methods and four interpolation methods. The transect locating methods are a morphologically‐spaced and an equally‐spaced model. The four interpolation methods are Ordinary Kriging, Delaunay Triangulation, and Simple Linear, which are applied in curvilinear coordinates (Delaunay Triangulation is also applied in Cartesian coordinates), and Natural Neighbor only in Cartesian Coordinates. The bathymetric data were obtained from morphologically simple and complex reaches of a large (average bankfull width = 90 m) and a small (average bankfull width = 17 m) river. The accuracy of the developed DEMs is assessed using statistical analysis of the differences between the control and interpolated bathymetries and hydraulic parameters assessed from bankfull water surface elevations. Results indicate that DEM accuracy is not influenced by the choice of transect location method (with same averaged cross‐section spacing) or a specific interpolation method, but rather by the coordinate system for which the interpolation method is applied and the spacing between transects. They also show negligible differences between the mean depths and surface areas calculated from bathymetries with dense or coarse spacing. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

9.
Floodplains and terraces in river valleys play important roles in the transport dynamics of water and sediment. While flat areas in river valleys can be identified from LiDAR data, directly characterizing them as either floodplain or terraces is not yet possible. To address this challenge, we hypothesize that, since geomorphic features are strongly coupled to hydrological and hydraulic dynamics and their associated variability, there exists a return frequency, or possibly a narrow band of return frequencies, of flow that is associated with floodplain formation; and this association can provide a distinctive signature for distinguishing them from terraces. Based on this hypothesis we develop a novel approach for distinguishing between floodplains and terraces that involves transforming the transverse cross‐sectional geometry of a river valley into a curve, named a river valley hypsometric (RVH) curve, and linking hydraulic inundation frequency with the features of this curve. Our approach establishes that the demarcation between floodplains and terraces can be established from the structure of steps and risers in the RVH curves which can be obtained from the DEM data. Further, it shows that these transitions may themselves be shaped by floods with 10‐ to 100‐year recurrence. We additionally show that, when floodplain width and height (above channel bottom) are normalized by bankfull width and depth, the ratio lies in a narrow range independent of the scale of the river valley. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

10.
  相似文献   

11.
This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km2) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km2 reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

12.
There is increasing demand for models that can accurately predict river temperature at the large spatial scales appropriate to river management. This paper combined summer water temperature data from a strategically designed, quality controlled network of 25 sites, with recently developed flexible spatial regression models, to understand and predict river temperature across a 3,000 km2 river catchment. Minimum, mean and maximum temperatures were modelled as a function of nine potential landscape covariates that represented proxies for heat and water exchange processes. Generalised additive models were used to allow for flexible responses. Spatial structure in the river network data (local spatial variation) was accounted for by including river network smoothers. Minimum and mean temperatures decreased with increasing elevation, riparian woodland and channel gradient. Maximum temperatures increased with channel width. There was greater between‐river and between‐reach variability in all temperature metrics in lower‐order rivers indicating that increased monitoring effort should be focussed at these smaller scales. The combination of strategic network design and recently developed spatial statistical approaches employed in this study have not been used in previous studies of river temperature. The resulting catchment scale temperature models provide a valuable quantitative tool for understanding and predicting river temperature variability at the catchment scales relevant to land use planning and fisheries management and provide a template for future studies.  相似文献   

13.
Field, laboratory, and numerical modelling research are increasingly demonstrating the potential of riparian tree colonization and growth to influence fluvial dynamics and the evolution of fluvial landforms. This paper jointly analyses multi‐temporal, multispectral ASTER data, continuous river stage and discharge data, and field observations of the growth rates of the dominant riparian tree species (Populus nigra) along a 21 km reach of the Tagliamento River, Italy. Research focuses on the period 2004–2009, during which there was a bankfull flood on 24 October 2004, followed by 2 years with low water levels, nearly 2 years with only modest flow pulses, and then a final period from 15 August 2008 that included several intermediate to bankfull flow events. This study period of increasing flow disturbance allows the exploration of vegetation dynamics within the river's active corridor under changing flow conditions. The analysis demonstrates the utility of ASTER data for investigating vegetation dynamics along large fluvial corridors and reveals both spatial and temporal variations in the expansion, coalescence, and erosion of vegetated patches within the study reach. Changes in the extent of the vegetated area and its dynamics vary along the study reach. In sub‐reaches where riparian tree growth is vigorous, the vegetated area expands rapidly during time periods without channel‐shaping flows, and is subsequently able to resist erosion by bankfull floods. In contrast, in sub‐reaches where tree growth is less vigorous, the vegetated area expands at a slower rate and is more readily re‐set by bankfull flood events. This illustrates that the rate of growth of riparian trees is crucial to their ability to contribute actively to river corridor dynamics. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

14.
Bankfull discharge is identified as an important parameter for studying river morphology, sediment motion, flood dynamics and their ecological impacts. In practice, the determination of this discharge and its hydrological characteristics is not easy, and a choice has to be made between several existing methods. To evaluate the impact of the choice of methods, five bankfull elevation definitions and four hydrological characterizations (determination of duration and frequency of exceedance applied to instantaneous or mean daily data) were compared on 16 gravel‐bed river reaches located in France (the catchment sizes vary from 10 km2 to 1700 km2). The consistency of bankfull discharge estimated at reach scale and the hydraulic significance of the five elevation definitions were examined. The morphological definitions (Bank Inflection, Top of Bank) were found more relevant than the definitions based on a geometric criterion. The duration of exceedance was preferred to recurrence intervals (partial duration series approach) because it is not limited by the independency of flood events, especially for low discharges like those associated with the Bank Inflection definition. On average, the impacts of the choice of methods were very important for the bankfull discharge magnitude (factor of 1·6 between Bank Inflection and Top of Bank) and duration of exceedance or frequency (respectively a factor 1·8 and 1·9 between mean daily and instantaneous discharge data). The choice of one combination of methods rather than another can significantly modify the conclusions of a comparative analysis in terms of bankfull discharge magnitude and its hydrological characteristics, so that one must be cautious when comparing results from different studies that use different methods. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
Lu Zhuo  Qiang Dai  Dawei Han 《水文研究》2015,29(11):2463-2477
Hydrological models play a significant role in modelling river flow for decision making support in water resource management. In the past decades, many researchers have made a great deal of efforts in calibrating and validating various models, with each study being focused on one or two models. As a result, there is a lack of comparative analysis on the performance of those models to guide hydrologists to choose appropriate models for the individual climate and physical conditions. This paper describes a two‐level meta‐analysis to develop a matching system between catchment complexity (based on catchment significant features (CSFs)) and model types. The intention is to use the available CSF information for choosing the most suitable model type for a given catchment. In this study, the CSFs include the elements of climate, soil type, land cover and catchment scale. Specific choices of model types in small and medium catchments are further explored with all CSF information obtained. In particular, it is interesting to find that semi‐distributed models are the most suitable model type for catchments with the area over 3000 km2, regardless of other CSFs. The potential methodology for expanding the matching system between catchment complexity and model complexity is discussed. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

16.
The channel boundary conditions along the Lower Yellow River (LYR) have been altered significantly since the 1950s with the continual reinforcement and construction of both main and secondary dykes and river training works. To evaluate how the confined complex channel–floodplain system of the LYR responds to floods, this study presents a detailed investigation of the relationship between the tempo‐spatial distribution of sedimentation/erosion and overbank floods occurred in the LYR. For large overbank floods, we found that when the sediment transport coefficient (ratio of sediment concentration of flow to flow discharge) is less than 0.034, the bankfull channel is subject to significant erosion, whereas the main and secondary floodplains both accumulate sediment. The amount of sediment deposited on the main and secondary floodplains is closely related to the ratio of peak discharge to bankfull discharge, volume of water flowing over the floodplains, and sediment concentration of overbank flow, whereas the degree of erosion in the bankfull channel is related to the amount of sediment deposited on the main and secondary floodplains, water volume, and sediment load in flood season. The significant increase in erosion in the bankfull channel is due to the construction of the main and secondary dykes and river training works, which are largely in a wide and narrow alternated pattern along the LYR such that the water flowing over wider floodplains returns to the channel downstream after it drops sediment. For small overbank floods, the bankfull channel is subject to erosion when the sediment transport coefficient is less than 0.028, whereas the amount of sediment deposited on the secondary floodplain is associated closely with the sediment concentration of flow. Over the entire length of the LYR, the situation of erosion in the bankfull channel and sediment deposition on the main and secondary floodplains occurred mainly in the upper reach of the LYR, in which a channel wandering in planform has been well developed.  相似文献   

17.
Past research investigated the surpassing of mean velocity at riffle cross sections by that at pool cross sections for flows up to bankfull, termed ‘velocity reversals’, to understand one mechanism by which riffle–pool relief is maintained. This study reenvisioned the classic velocity reversal by documenting stage‐dependent changes to the locations of peak velocity without cross sections. Instead, the dynamics of peak velocity patches were considered for flows spanning 0.2 to 22 times bankfull discharge through the use of a high‐resolution DEM and two‐dimensional hydrodynamic modelling. A remarkable diversity in peak velocity patch behaviour was found across discharges, including gradual expansion and shifting as well as abrupt disappearance and emergence relative to the low‐flow patch locations. These behaviours blended together to varying degrees to produce many reversals in peak velocity across morphological units, but it took substantially higher than bankfull discharge for peak velocities to move from riffles and chutes to fast glides and pools. The discharges at which reversals occurred among morphological units were significantly higher for the valley‐confined reach than for the anastomosing reach. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

18.
Stream-water temperature is a key variable controlling chemical, biological, and ecological processes in freshwater environments. Most models focus on a single river cross-section; however, temperature gradients along stretches and tributaries of a river network are crucial to assess ecohydrological features such as aquatic species suitability, growth and feeding rates, or disease transmission. We propose SESTET, a deterministic, spatially explicit stream temperature model for a whole river network, based on water and energy budgets at a reach scale and requiring only commonly available spatially distributed datasets, such as morphology and air temperature, as input. Heat exchange processes at the air–water interface are modelled via the widely used equilibrium temperature concept, whereas the effects of network structure are accounted for through advective heat fluxes. A case study was conducted on the prealpine Wigger river (Switzerland), where water temperatures have been measured in the period 2014–2018 at 11 spatially distributed locations. The results show the advantages of accounting for water and energy budgets at the reach scale for the entire river network, compared with simpler, lumped formulations. Because our approach fundamentally relies on spatially distributed air temperature fields, adequate spatial interpolation techniques that account for the effects of both elevation and thermal inversion in air temperature are key to a successful application of the model. SESTET allows the assessment of the magnitude of the various components of the heat budget at the reach scale and the derivation of reliable estimates of spatial gradients of mean daily stream temperatures for the whole catchment based on a limited number of conveniently located (viz., spanning the largest possible elevation range) measuring stations. Moreover, accounting for mixing processes and advective fluxes through the river network allows one to trust regionalized values of the parameters controlling the relationship between equilibrium and air temperature, a key feature to generalize the model to data-scarce catchments.  相似文献   

19.
Stream invertebrate distribution patterns reflect local sedimentary and hydraulic conditions, which in turn are influenced by a range of factors operating at larger scales. We assessed whether spatial variation in invertebrate assemblages across a meso-scale catchment is best understood in terms of the characteristics of the study reaches themselves or the characteristics of respective upstream catchment areas. The study river experiences naturally high fine sediment loads as a result of the extraordinary supply of sediment from high erodible marls in its catchment. We hypothesized that between-reach variation in the volume of fine sediment stored within the channel results from a combination of reach and upstream catchment characteristics, and that these characteristics help explain variation in invertebrate assemblages. The storage of fine sediment in study reaches correlated with a number of upstream catchment characteristics, as well as reach-scale hydraulic conditions. Variability in invertebrate assemblages correlated most strongly (62% of variance explained) with the characteristics of the catchment upstream from each reach (area of contribution), with the characteristics of the reaches accounting for only 35% of the variability. The explanatory power of the reach-scale habitat variables was reduced when the effect of upstream catchment conditions was removed. This suggests inbuilt effects of larger scale conditions on reach habitat and invertebrate assemblages. Results lend support to theories of scale hierarchy within river systems and help emphasize the need to target management at upstream catchment areas.  相似文献   

20.
This article presents results from an investigation of the hydraulic characteristics of overbank flows on topographically‐complex natural river floodplains. A two‐dimensional hydraulic model that solves the depth‐averaged shallow water form of the Navier–Stokes equations is used to simulate an overbank flow event within a multiple channel reach of the River Culm, Devon, UK. Parameterization of channel and floodplain roughness by the model is evaluated using monitored records of main channel water level and point measurements of floodplain flow depth and unit discharge. Modelled inundation extents and sequences are assessed using maps of actual inundation patterns obtained using a Global Positioning System, observational evidence and ground photographs. Simulation results suggest a two‐phase model of flooding at the site, which seems likely to be representative of natural floodplains in general. Comparison of these results with previous research demonstrates the complexity of overbank flows on natural river floodplains and highlights the limitations of laboratory flumes as an analogue for these environments. Despite this complexity, frequency distributions of simulated depth, velocity and unit discharge data closely follow a simple gamma distribution model, and are described by a shape parameter (α) that exhibits clear systematic trends with changing discharge and floodplain roughness. Such statistical approaches have the potential to provide the basis for computationally efficient flood routing and overbank sedimentation models. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

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

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