首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Upland agricultural land management activities such as grazing, vegetation burning, and bare ground restoration impact hydrological elements of headwater catchments, many of which may be important for downstream flood peaks (e.g., overland flow and soil water storage). However, there is poor understanding of how these management practices affect river flow peaks during high magnitude rainfall events. Using the distributed TOPMODEL, spatial configurations of land management were modelled to predict flood response in an upland catchment, which contains different regions operating subsidized agricultural stewardship schemes. Heavy grazing leading to soil compaction and loss of vegetation cover in stewardship regions covering 79.8% of the catchment gave a 42‐min earlier flow peak, which was 82.2% higher (under a 1‐hr 15‐mm storm) than the current simulated hydrograph. Light grazing over the same regions of the catchment had much less influence on river flow peaks (18 min earlier and 32.9% increase). Rotational burning (covering 8.8% of the catchment), most of which is located in the headwater areas, increased the peak by 3.2% in the same rainfall event. Vegetation restoration with either Eriophorum or Sphagnum (higher density) in bare areas (5.8%) of the catchment provided a reduction of flood peak (3.9% and 5.2% in the 15‐mm storm event), whereas the same total area revegetated with Sphagnum in riparian regions delivered a much larger decrease (15.0%) in river flow peaks. We show that changes of vegetation cover in highly sensitive areas (e.g., near‐stream zones) generate large impacts on flood peaks. Thus, it is possible to design spatially distributed management systems for upland catchments, which reduce flood peaks while at the same time ensuring economic viability for upland farmers.  相似文献   

2.
Floods in small mountainous watersheds cover a wide spectrum of flow. They can range from clear water flows and hyperconcentrated flows to debris floods and debris flows, and calculation of the peak discharge is crucial for predicting and mitigating such hazards. To determine the optimal approach for discharge estimation, this study compared water flow monitoring hydrographs to investigate the performance of five hydrological models that incorporate different runoff yields and influx calculation methods. Two of the models performed well in simulating the peak discharge, peak time, and total flow volume of the water flood. The ratio (γ) of the monitored debris flood discharge (Qd) to the simulated water flow discharge (Qw) was investigated. Qualitatively, γ initially increased with Qw but then decreased when Qw exceeded a certain threshold, which corresponded to rainfall of 95 and 120 mm in a 6- and 24-h event with a normal distribution of precipitation, respectively. The decrease might be attributable to a threshold of sediment availability being reached, beyond which increased flow rate is not matched by increased sediment input in the large watershed. Uncertainty of hydrological calculation was evaluated by dividing the catchment into sub-basins and adopting different rainfall time steps as input. The efficiency of using a distributed simulation exhibited marginal improvement potential compared with a lumped simulation. Conversely, the rainfall time step input significantly affected the simulation results by delaying the peak time and decreasing the peak discharge. This research demonstrates the applicability of a discharge estimation method that combines a hydrological water flow simulation and an estimation of γ. The results were verified on the basis of monitored flow densities and videos obtained in two watersheds with areas of 2.34 and 32.4 km2.  相似文献   

3.
Abstract

Flood frequency estimation is crucial in both engineering practice and hydrological research. Regional analysis of flood peak discharges is used for more accurate estimates of flood quantiles in ungauged or poorly gauged catchments. This is based on the identification of homogeneous zones, where the probability distribution of annual maximum peak flows is invariant, except for a scale factor represented by an index flood. The numerous applications of this method have highlighted obtaining accurate estimates of index flood as a critical step, especially in ungauged or poorly gauged sections, where direct estimation by sample mean of annual flood series (AFS) is not possible, or inaccurate. Therein indirect methods have to be used. Most indirect methods are based upon empirical relationships that link index flood to hydrological, climatological and morphological catchment characteristics, developed by means of multi-regression analysis, or simplified lumped representation of rainfall–runoff processes. The limits of these approaches are increasingly evident as the size and spatial variability of the catchment increases. In these cases, the use of a spatially-distributed, physically-based hydrological model, and time continuous simulation of discharge can improve estimation of the index flood. This work presents an application of the FEST-WB model for the reconstruction of 29 years of hourly streamflows for an Alpine snow-fed catchment in northern Italy, to be used for index flood estimation. To extend the length of the simulated discharge time series, meteorological forcings given by daily precipitation and temperature at ground automatic weather stations are disaggregated hourly, and then fed to FEST-WB. The accuracy of the method in estimating index flood depending upon length of the simulated series is discussed, and suggestions for use of the methodology provided.
Editor D. Koutsoyiannis  相似文献   

4.
Abstract

A modelling experiment is used to examine different land-use scenarios ranging from extreme deforestation (31% forest cover) to pristine (95% forest cover) conditions and related Payment for Ecosystem Services (PES) schemes to assess whether a change in streamflow dynamics, discharge extremes and mean annual water balance of a 73.4-km2 tropical headwater catchment in Costa Rica could be detected. A semi-distributed, conceptual rainfall–runoff model was adapted to conceptualize the empirically-based, dominant hydrological processes of the study area and was multi-criteria calibrated using different objective functions and empirical constraints on model simulations in a Monte Carlo framework to account for parameter uncertainty. The results suggest that land-use change had relatively little effect on the overall mean annual water yield (<3%). However, streamflow dynamics proved to be sensitive in terms of frequency, timing and magnitude of discharge extremes. For low flows and peak discharges of return periods greater than one year, land use had a minor influence on the runoff response. Below these thresholds (<1-year return period), forest cover potentially decreased runoff peaks and low flows by as much as 10%, and non-forest cover increased runoff peaks and low flows by up to 15%. The study demonstrated the potential for using hydrological modelling to help identify the impact of protection and reforestation efforts on ecosystem services.

Editor Z.W. Kundzewicz

Citation Birkel, C., Soulsby, C., and Tetzlaff, D., 2012. Modelling the impacts of land-cover change on streamflow dynamics of a tropical rainforest headwater catchment. Hydrological Sciences Journal, 57 (8), 1543–1561.  相似文献   

5.
Abstract

The design and construction of a special-purpose laboratory catchment and rainfall simulator is described. The equipment consists of a soil catchment area that can be inclined at various angles. Additional instrumentation then measures the flow of water across the surface of, and through, the soil bed. Precipitation is provided by a unit that simulates rainfall at particular rates with uniform distribution.

The equipment was used to examine infiltration, runoff and other hydrological properties of a number of soils under different rainfall intensities and with different catchment slopes. Correlations were obtained for these variables.  相似文献   

6.
Abstract

Agricultural use and related water erosion may lead to significant changes in the sedimentological and hydrological characteristics of watersheds, and therefore negative consequences for rural development. This research aimed to put present-day soil erosion of the important Mejerda catchment into a historical context. The catchment of Wadi Mejerda in northern Tunisia has experienced soil erosion due to weather and human impacts for thousands of years. We used historical texts and results from archaeological research that go back to 1000 BC, as well as data collected during the last century. Soil erosion from different types of agricultural landscape management was analysed together with information on the soils' production potential, the hydrographic network and flood frequency. The results showed that water erosion has increased the hydrographic network by 65 km and increased the deltaic plain by as much as 15 km2/century. However, soil productivity has decreased significantly. Moreover, due to in channel sedimentation and river choking, the number of flooding occurrences has multiplied over the last century. Finally, it is shown that water erosion follows a specific cycle of degradation throughout the watershed. These findings should be considered for better water and soil management in the context of semi-arid areas.

Editor Z.W. Kundzewicz

Citation Jebari, S., Berndtsson, R., Lebdi, F., and Bahri, A., 2012. Historical aspects of soil erosion in the Mejerda catchment. Hydrological Sciences Journal, 57 (5), 901–912.  相似文献   

7.
High resolution radar rainfall fields and a distributed hydrologic model are used to evaluate the sensitivity of flood and flash flood simulations to spatial aggregation of rainfall and soil properties at catchment scales ranging from 75 to 983 km2. Hydrologic modeling is based on a Hortonian infiltration model and a network-based representation of hillslope and channel flow. The investigation focuses on three extreme flood and flash flood events occurred on the Sesia river basin, North Western Italy, which are analysed by using four aggregation lengths ranging from 1 to 16 km. The influence of rainfall spatial aggregation is examined by using the flow distance as a spatial coordinate, hence emphasising the role of river network in the averaging of space–time rainfall. The effects of reduced and distorted rainfall spatial variability on peak discharge have been found particularly severe for the flash flood events, with peak errors up to 35% for rainfall aggregation of 16 km and at 983 km2 catchment size. Effects are particularly remarkable when significant structured rainfall variability combines with relatively important infiltration volumes due to dry initial conditions, as this emphasises the non-linear character of the rainfall–runoff relationship. In general, these results confirm that the correct estimate of rainfall volume is not enough for the accurate reproduction of flash flood events characterised by large and structured rainfall spatial variability, even at catchment scales around 250 km2. However, accurate rainfall volume estimation may suffice for less spatially variable flood events. Increasing the soil properties aggregation length exerts similar effects on peak discharge errors as increasing the rainfall aggregation length, for the cases considered here and after rescaling to preserve the rainfall volume. Moreover, peak discharge errors are roughly proportional to runoff volume errors, which indicates that the shape of the flood wave is influenced in a limited way by modifying the detail of the soil property spatial representation. Conversely, rainfall aggregation may exert a pronounced influence on the discharge peak by reshaping the spatial organisation of the runoff volumes and without a comparable impact on the runoff volumes.  相似文献   

8.
In this paper, the controls of different indicators on the statistical moments (i.e. mean annual flood (MAF), coefficient of variation (CV) and skewness (CS)) of the maximum annual flood records of 459 Austrian catchments are analysed. The process controls are analysed in terms of the correlation of the flood moments within five hydrologically homogeneous regions to two different types of indicators. Indicators of the first type are static catchment attributes, which are associated with long‐term observations such as mean annual precipitation, the base flow index, and the percentage of catchment area covered by a geological unit or soil type. Indicators of the second type are dynamic catchment attributes that are associated with the event scale. Indicators of this type used in the study are event runoff coefficients and antecedent rainfall. The results indicate that MAF and CV are strongly correlated with indicators characterising the hydro‐climatic conditions of the catchments, such as mean annual precipitation, long‐term evaporation and the base flow index. For the catchments analysed, the flood moments are not significantly correlated with static catchment attributes representing runoff generation, such as geology, soil types, land use and the SCS curve number. Indicators of runoff generation that do have significant predictive power for flood moments are dynamic catchment attributes such as the mean event runoff coefficients and mean antecedent rainfall. The correlation analysis indicates that flood runoff is, on average, more strongly controlled by the catchment moisture state than by event rainfall. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

9.
Book reviews     
Abstract

Statistical and deterministic modelling estimates of flood magnitudes and frequencies that can affect flood-plain ecology in the upper Ahuriri River catchment, a mountainous high country catchment in the New Zealand Southern Alps, were evaluated. Statistical analysis of 46 years of historical data showed that floods are best modelled by the generalized extreme value and lognormal distributions. We evaluated application of the HEC-HMS model to this environment by modelling flood events of various frequencies. Model results were validated and compared with the statistical estimates. The SCS curve number method was used for losses and runoff generation, and the model was very sensitive to curve number. The HEC-HMS flood estimates matched the statistical estimates reasonably well, and, over all return periods, were on average approximately 1% greater. However, the model generally underestimated flood peaks up to the 25-year event and overestimated magnitudes above this. The results compared well with other regional estimates, including studies based on L-moments, and showed that this catchment has smaller floods than other similarly-sized catchments in the Southern Alps.

Editor D. Koutsoyiannis; Associate editor H. Aksoy

Citation Caruso, B.S., Rademaker, M., Balme, A., and Cochrane, T.A., 2013. Flood modelling in a high country mountain catchment, New Zealand: comparing statistical and deterministic model estimates for ecological flows. Hydrological Sciences Journal, 58 (2), 328–341.  相似文献   

10.
E. Morin  H. Yakir 《水文科学杂志》2014,59(7):1353-1362
Abstract

t Spatio-temporal storm properties have a large impact on catchment hydrological response. The sensitivity of simulated flash floods to convective rain-cell characteristics is examined for an extreme storm event over a 94 km2 semi-arid catchment in southern Israel. High space–time resolution weather radar data were used to derive and model convective rain cells that then served as input into a hydrological model. Based on alterations of location, direction and speed of a major rain cell, identified as the flooding cell for this case, the impacts on catchment rainfall and generated flood were examined. Global sensitivity analysis was applied to identify the most important factors affecting the flash flood peak discharge at the catchment outlet. We found that the flood peak discharge could be increased three-fold by relatively small changes in rain-cell characteristics. We assessed that the maximum flash flood magnitude that this single rain cell can produce is 175 m3/s, and, taking into account the rest of the rain cells, the flash flood peak discharge can reach 260 m3/s.
Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Morin, E. and Yakir, H., 2013. Hydrological impact and potential flooding of convective rain cells in a semi-arid environment. Hydrological Sciences Journal, 59 (7), 1275–1284. http://dx.doi.org/10.1080/02626667.2013.841315  相似文献   

11.
An essential part of hydrological research focuses on hydrological extremes, such as river peak flows and associated floods, because of their large impact on economy, environment, and human life. These extremes can be affected by potential future environmental change, including global climate change and land cover change. In this paper, the relative impact of both climate change and urban expansion on the peak flows and flood extent is investigated for a small‐scale suburban catchment in Belgium. A rainfall‐runoff model was coupled to a hydrodynamic model in order to simulate the present‐day and future river streamflow. The coupled model was calibrated based on a series of measured water depths and, after model validation, fed with different climate change and urban expansion scenarios in order to evaluate the relative impact of both driving factors on the peak flows and flood extent. The three climate change scenarios that were used (dry, wet winter, wet summer) were based on a statistical downscaling of 58 different RCM and GCM scenario runs. The urban expansion scenarios were based on three different urban growth rates (low, medium, high urban expansion) that were set up by means of an extrapolation of the observed trend of urban expansion. The results suggest that possible future climate change is the main source of uncertainty affecting changes in peak flow and flood extent. The urban expansion scenarios show a more consistent trend. The potential damage related to a flood is, however, mainly influenced by land cover changes that occur in the floodplain. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

12.
Abstract

This article presents a comparison between real-time discharges calculated by a flash-flood warning system and post-event flood peak estimates. The studied event occurred on 15 and 16 June 2010 at the Argens catchment located in the south of France. Real-time flood warnings were provided by the AIGA (Adaptation d’Information Géographique pour l’Alerte en Crue) warning system, which is based on a simple distributed hydrological model run at a 1-km2 resolution using radar rainfall information. The timing of the warnings (updated every 15 min) was compared to the observed flood impacts. Furthermore, “consolidated” flood peaks estimated by an intensive post-event survey were used to evaluate the AIGA-estimated peak discharges. The results indicated that the AIGA warnings clearly identified the most affected areas. However, the effective lead-time of the event detection was short, especially for fast-response catchments, because the current method does not take into account any rainfall forecast. The flood peak analysis showed a relatively good correspondence between AIGA- and field-estimated peak values, although some differences were due to the rainfall underestimation by the radar and rainfall–runoff model limitations.
Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Javelle, P., Demargne, J., Defrance, D., Pansu, J. and Arnaud, P., 2014. Evaluating flash-flood warnings at ungauged locations using post-event surveys: a case study with the AIGA warning system. Hydrological Sciences Journal, 59 (7), 1390–1402. http://dx.doi.org/10.1080/02626667.2014.923970  相似文献   

13.
ABSTRACT

The spread of impervious surfaces in urban areas combined with the rise in the intensity of rainfall events as a result of climate change has led to dangerous increases in storm water flows. This paper discusses a new implementation of the fully distributed hydrological model Multi-Hydro (developed at École des Ponts ParisTech), when operating storage basins, and its ability to deal with high-resolution radar rainfall data. The peri-urban area of Massy (south of Paris, France) was selected as a case study for having six of these drainage facilities, used extensively in flood control. Two radar rainfall datasets with different spatiotemporal resolutions were used: Météo-France’s PANTHER rainfall product (C-band) and ENPC’s X-band DPSRI. The rainfall spatiotemporal variability was analysed statistically using Universal Multifractals (UM). Finally, to validate the application, the water level simulations were compared with local measurements in the Cora storage basin located next to the catchment’s single outlet.  相似文献   

14.
Abstract

The physically-based flood frequency models use readily available rainfall data and catchment characteristics to derive the flood frequency distribution. In the present study, a new physically-based flood frequency distribution has been developed. This model uses bivariate exponential distribution for rainfall intensity and duration, and the Soil Conservation Service-Curve Number (SCS-CN) method for deriving the probability density function (pdf) of effective rainfall. The effective rainfall-runoff model is based on kinematic-wave theory. The results of application of this derived model to three Indian basins indicate that the model is a useful alternative for estimating flood flow quantiles at ungauged sites.  相似文献   

15.
ABSTRACT

When discharge measurements are not available, design of water structures relies on using frequency analysis of rainfall data and applying a rainfall–runoff model to estimate a hydrograph. The Soil Conservation Service (SCS) method estimates the design hydrograph first through a rainfall–runoff transformation and next by propagating runoff to the basin outlet via the SCS unit hydrograph (UH) method. The method uses two parameters, the Curve Number (CN) and the time of concentration (Tc). However, in data-scarce areas, the calibration of CN and Tc from nearby gauged watersheds is limited and subject to high uncertainties. Therefore, the inherent uncertainty/variability of the SCS parameters may have considerable ramifications on the safety of design. In this research, a reliability approach is used to evaluate the impact of incorporating the uncertainty of CN and Tc in flood design. The sensitivity of the probabilistic outcome against the uncertainty of input parameters is calculated using the First Order Reliability Method (FORM). The results of FORM are compared with the conventional SCS results, taking solely the uncertainty of the rainfall event. The relative importance of the uncertainty of the SCS parameters is also estimated. It is found that the conventional approach, used by many practitioners, might grossly underestimate the risk of failure of water structures, due to neglecting the probabilistic nature of the SCS parameters and especially the Curve Number. The most predominant factors against which the SCS-CN method is highly uncertain are when the average rainfall value is low (less than 20 mm) or its coefficient of variation is not significant (less than 0.5), i.e. when the resulting rainfall at the design return period is low. A case study is presented for Egypt using rainfall data and CN values driven from satellite information, to determine the regions of acceptance of the SCS-CN method.
EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR A. Efstratiadis  相似文献   

16.
Yi-Ru Chen  Bofu Yu 《水文科学杂志》2013,58(10):1759-1769
Abstract

Over the past century, land-use has changed in southeast Queensland, and when coupled with climatic change, the risk of flooding has increased. This research aims to examine impacts of climate and land-use changes on flood runoff in southeast Queensland, Australia. A rainfall–runoff model, RORB, was calibrated and validated using observed flood hydrographs for one rural and one urbanized catchment, for 1961–1990. The validated model was then used to generate flood hydrographs using projected rainfall based on two climate models: the Geophysical Fluid Dynamics Laboratory Climate Model 2.1 (GFDL CM2.1) and the Conformal-Cubic Atmospheric Model (CCAM), for 2016–2045. Projected daily rainfall for the two contrasting periods was used to derive adjustment factors for a given frequency of occurrence. Two land-use change scenarios were used to evaluate likely impacts. Based on the projected rainfall, the results showed that, in both catchments, future flood magnitudes are unlikely to increase for large flood events. Extreme land-use change would significantly impact flooding in the rural catchment, but not the urbanized catchment.
Editor Z.W. Kundzewicz; Associate editor Y. Gyasi-Agyei  相似文献   

17.
《水文科学杂志》2013,58(6):1208-1220
Abstract

Betel nut or betel palm (Areca catechu Linn.) has become a major cash crop in southern and central Taiwan since the mid-1980s. Many slopeland fruit orchards and forests have been converted to betel nut plantations. The total area of betel nut plantations has increased over the past 30 years, reaching a peak of 56 542 ha in 1997. The public and conservation groups frequently express great concern over the potential negative hydrological impacts of betel nut plantations on steep slopes. This investigation in central Taiwan examines the effects of hillslope betel nut plantations on hydrological processes. Differences in hydrological characteristics between plots planted with betel nut trees and those with other ground cover types were evaluated at two study sites. The tall, single-layer canopy and wide spacing between planted betel nut trees led to low interception losses, high throughfall and high net rainfall. Plots planted with betel nut trees had lower infiltration, higher surface runoff and higher erosion than forested sites. These hydrological characteristics can be related to factors such as crown cover, soil organic content and soil porosity in betel nut plantations. Streamflow data from three gauged watersheds with different proportions of total area planted with betel nut showed that where greater proportions of total area were planted with betel nut trees, there were higher annual streamflow/rainfall ratios, higher specific peak flows, steeper recessions and higher peak flow/baseflow ratios. The results from this study suggest that, in general, betel nut trees are less desirable from the soil and water conservation viewpoints than natural forests.  相似文献   

18.
A statistically based runoff‐yield model is proposed in this paper. The model considers spatial heterogeneities of rainfall, soil infiltration capacity and soil water storage capacity that are main factors controlling runoff‐yield process. It assumes that the spatial variation of rainfall intensity at each time step can be characterized by a probability density function, which is estimated by matching the hyetograph through goodness‐of‐fit measure, whereas the spatial heterogeneities of soil infiltration capacity and soil water storage capacity are described by parabola‐type functions. Surface runoff is calculated according to infiltration excess mechanism; the statistical distribution of surface runoff rate can be deduced with the joint distribution of rainfall intensity and soil infiltration rate, thus obtaining a quasi‐analytical solution for surface runoff. Based on saturation excess mechanism, the groundwater flow (flows below the ground are collectively referred to as groundwater flow) is calculated by infiltration and the probability distribution of soil water storage capacity. Consequently, the total runoff is composed of infiltration excess and saturation excess runoff components. As an example, this model is applied to flood event simulation in Dongwan catchment, a semi‐humid region and a tributary of Yellow River in China. It indicates that the proposed runoff‐yield model could achieve acceptable accuracy. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
Abstract

The hydrological response of a small agroforestry catchment in northwest Spain (Corbeira catchment, 16 km2) is analysed, with particular focus on rainfall events. Fifty-four rainfall–runoff events, from December 2004 to September 2007, were used to analyse the principal hydrological patterns and show which factors best explain the hydrological response. The nonlinearity between rainfall and runoff showed that the variability in the hydrological response of the catchment was linked to the seasonal dynamics of the rainfall and, to a lesser extent, to evapotranspiration. The runoff coefficient, estimated as the ratio between direct runoff and rainfall volume, on an event basis, was analysed as a function of rainfall characteristics (amount and intensity) and the initial catchment state conditions prior to an event, such as pre-event baseflow and antecedent rainfall index. The results revealed that the hydrological response depends both on the soil humidity conditions at the start of the event and on rainfall amount, whereas rainfall intensity presented only a significant correlation with discharge increment. The antecedent conditions seem to be a key point in runoff production, and they explain much of the response. The hydrographs are characterized by a steep rising limb, a relatively narrow peak discharge and slow recession limb. These data and the observations suggest that the subsurface flow is the dominant runoff process.

Editor Z.W. Kundzewicz; Associate editor T. Wagener

Citation Rodríguez-Blanco, M.L., Taboada-Castro, M.M. and Taboada-Castro, M.T., 2012. Rainfall–runoff response and event-based runoff coefficients in a humid area (northwest Spain). Hydrological Sciences Journal, 57 (3), 445–459.  相似文献   

20.
Runoff and peak flows in three experimental catchments with different forest conditions were analysed in a rainy temperate climate in southern Chile. The hydrological effects of clearcutting a Pinus radiata plantation covering 79·4% of the La Reina catchment were studied by analysing runoff and peak flows in the pre‐ and post‐harvesting periods. Mean annual runoff increased 110% after timber harvesting. Clearcutting generated a 32% mean increase in peak flows. Peak flow and runoff were examined in two adjacent catchments with different forest conditions. The older plantation in Los Ulmos 1 increasingly consumed more water than the younger plantation established at Los Ulmos 2, whereas differences in peak flows between these two catchments were not significant. Runoff and peak flows comparisons not only reflected changes in forest cover, but also the effect of rainfall characteristics during the study periods and the basins' morphologies. Comparisons between pre‐ and post‐harvesting peak discharges must be undertaken with caution, because a previous analysis at La Reina using a partial set of data overestimated changes in peak flows after timber harvesting. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

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

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