首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 336 毫秒
1.
Radar‐based estimates of rainfall are affected by many sources of uncertainties, which would propagate through the hydrological model when radar rainfall estimates are used as input or initial conditions. An elegant solution to quantify these uncertainties is to model the empirical relationship between radar measurements and rain gauge observations (as the ‘ground reference’). However, most current studies only use a fixed and uniform model to represent the uncertainty of radar rainfall, without consideration of its variation under different synoptic regimes. Wind is such a typical weather factor, as it not only induces error in rain gauge measurements but also causes the raindrops observed by weather radar to drift when they reach the ground. For this reason, as a first attempt, this study introduces the wind field into the uncertainty model and designs the radar rainfall uncertainty model under different wind conditions. We separate the original dataset into three subsamples according to wind speed, which are named as WDI (0–2 m/s), WDII (2–4 m/s) and WDIII (>4 m/s). The multivariate distributed ensemble generator is introduced and established for each subsample. Thirty typical events (10 at each wind range) are selected to explore the behaviours of uncertainty under different wind ranges. In each time step, 500 ensemble members are generated, and the values of 5th to 95th percentile values are used to produce the uncertainty bands. Two basic features of uncertainty bands, namely dispersion and ensemble bias, increase significantly with the growth of wind speed, demonstrating that wind speed plays a considerable role in influencing the behaviour of the uncertainty band. On the basis of these pieces of evidence, we conclude that the radar rainfall uncertainty model established under different wind conditions should be more realistic in representing the radar rainfall uncertainty. This study is only a start in incorporating synoptic regimes into rainfall uncertainty analysis, and a great deal of more effort is still needed to build a realistic and comprehensive uncertainty model for radar rainfall data. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.  相似文献   

3.
Automated rainfall simulator for variable rainfall on urban green areas   总被引：1，自引：0，他引：1  
Rainfall simulators can enhance our understanding of the hydrologic processes affecting the total runoff to urban drainage systems. This knowledge can be used to improve urban drainage designs. In this study, a rainfall simulator is developed to simulate rainfall on urban green surfaces. The rainfall simulator is controlled by a microcomputer programmed to replicate the temporal variations in rainfall intensity of both historical and synthetic rainfall events with constant rainfall intensity on an area of 1 m2. The performance of the rainfall simulator is tested under laboratory conditions with regard to spatial uniformity of the rainfall, the kinetic energy of the raindrops, and the ability to replicate historical and synthetic rainfall events with temporally varying intensity. The rainfall simulator is applied in the field to evaluate its functionality under field conditions and the influence of wind on simulated rainfall. Finally, a field study is carried out on the relationship between runoff, soil volumetric water content, and surface slope. Performance and field tests show that the simulated rainfall has a uniform spatial distribution, whereas the kinetic energy of the raindrops is slightly higher than that of other comparable rainfall simulators. The rainfall simulator performs best in low wind speed conditions. The simulator performs well in replicating historical and synthetic rainfall events by matching both intensity variations and accumulated rainfall depth. The field study shows good correlation between rainfall, runoff, infiltration, soil water content, and surface slope.  相似文献   

4.
Rain‐gauge catch efficiencies are affected by wind. Wind makes raindrops fall at an angle of inclination and the effective diameter of the rain gauge orifice smaller than if raindrops fall into the gauge vertically. Two spherical and two semi‐spherical orifices were designed to modify standard gauges and others in use today. The two spherical orifices catch rain with an effective diameter always equal to the actual diameter regardless of wind speed and direction. The semi‐spherical orifices, used side‐by‐side with a standard gauge, correct 50% of catch deficiencies made by the standard gauge. Tests based on 115 storms show that the four new gauges caught more rainfall than the standard gauge, with an average catch increase ranging from 8% to 16%. Compared with the pit gauge, average deficiency in catch ranged from ?1% (spherical rain gauge orifice with cylinders) to 4%, whereas the deficiency for the standard gauge was ?10%. Percentage deficiencies of the new gauges were positively affected by wind speed, raindrop inclination and rainfall intensity. Although the new gauges tended to underestimate the standard gauge in small storms (<0·25 cm) and overestimated the pit gauge under strong winds, their deviations are small. Underestimates for small storms could be improved by using gauge materials that reduce surface temperature, evaporation and water retention. The gauges are simple in design, easy to operate and inexpensive. In order to maintain a historically consistent set of rainfall data, a dual‐gauge (standard gauge + spherical gauge) is recommended for existing rainfall stations. The new rain gauge orifices are suitable for large‐scale applications. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
This paper reports the results of an investigation into flood simulation by areal rainfall estimated from the combination of gauged and radar rainfalls and a rainfall–runoff model on the Anseong‐cheon basin in the southern part of Korea. The spatial and temporal characteristics and behaviour of rainfall are analysed using various approaches combining radar and rain gauges: (1) using kriging of the rain gauge alone; (2) using radar data alone; (3) using mean field bias (MFB) of both radar and rain gauges; and (4) using conditional merging technique (CM) of both radar and rain gauges. To evaluate these methods, statistics and hyetograph for rain gauges and radar rainfalls were compared using hourly radar rainfall data from the Imjin‐river, Gangwha, rainfall radar site, Korea. Then, in order to evaluate the performance of flood estimates using different rainfall estimation methods, rainfall–runoff simulation was conducted using the physics‐based distributed hydrologic model, Vflo?. The flood runoff hydrograph was used to compare the calculated hydrographs with the observed one. Results show that the rainfall field estimated by CM methods improved flood estimates, because it optimally combines rainfall fields representing actual spatial and temporal characteristics of rainfall. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

6.
This paper provides a comparison of gauge and radar precipitation data sources during an extreme hydrological event. November–December 2006 was selected as a time period of intense rainfall and large river flows for the Severn Uplands, an upland catchment in the United Kingdom. A comparison between gauge and radar precipitation time‐series records for the event indicated discrepancies between data sources, particularly in areas of higher elevation. The HEC‐HMS rainfall‐runoff model was selected to assess the accuracy of the precipitation to simulate river flows for the extreme event. Gauge, radar and gauge‐corrected radar rainfall were used as model inputs. Universal cokriging was used to geostatistically interpolate gauge data with radar and elevation data as covariates. This interpolated layer was used to calculate the mean‐field bias and correct the radar composites. Results indicated that gauge‐ and gauge‐corrected radar‐driven models replicated flows adequately for the extreme event. Gauge‐corrected flow predictions produced an increase in flow prediction accuracy when compared with the raw radar, yet predictions were comparative in accuracy to those using the interpolated gauge network. Subsequent investigations suggested this was due to an adequate spatial and temporal resolution of the precipitation gauge network within the Severn Uplands. Results suggested that the six rain gauges could adequately represent precipitation variability of the Severn Uplands to predict flows at an approximately equal accuracy to that obtained by radar. Temporally, radar produced an increase in flow prediction accuracy in mountainous reaches once the gauge time step was in excessive of an hourly interval. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

7.
The feasibility of linear and nonlinear geostatistical estimation techniques for optimal merging of rainfall data from raingage and radar observations is investigated in this study by use of controlled numerical experiments. Synthetic radar and raingage data are generated with their hypothetical error structures that explicitly account for sampling characteristics of the two sensors. Numerically simulated rainfall fields considered to be ground-truth fields on 4×4 km grids are used in the generation of radar and raingage observations. Ground-truth rainfall fields consist of generated rainfall fields with various climatic characteristics that preserve the space-time covariance function of rainfall events in extratropical cyclonic storms. Optimal mean areal precipitation estimates are obtained based on the minimum variance, unbiased property of kriging techniques under the second order homogeneity assumption of rainfall fields. The evaluation of estimated rainfall fields is done based on the refinement of spatial predictability over what would be provided from each sensor individually. Attention is mainly given to removal of measurement error and bias that are synthetically introduced to radar measurements. The influence of raingage network density on estimated rainfall fields is also examined.  相似文献   

8.
The feasibility of linear and nonlinear geostatistical estimation techniques for optimal merging of rainfall data from raingage and radar observations is investigated in this study by use of controlled numerical experiments. Synthetic radar and raingage data are generated with their hypothetical error structures that explicitly account for sampling characteristics of the two sensors. Numerically simulated rainfall fields considered to be ground-truth fields on 4×4 km grids are used in the generation of radar and raingage observations. Ground-truth rainfall fields consist of generated rainfall fields with various climatic characteristics that preserve the space-time covariance function of rainfall events in extratropical cyclonic storms. Optimal mean areal precipitation estimates are obtained based on the minimum variance, unbiased property of kriging techniques under the second order homogeneity assumption of rainfall fields. The evaluation of estimated rainfall fields is done based on the refinement of spatial predictability over what would be provided from each sensor individually. Attention is mainly given to removal of measurement error and bias that are synthetically introduced to radar measurements. The influence of raingage network density on estimated rainfall fields is also examined.  相似文献   

9.
In this study, the correction problem of mean‐field bias of radar rain rate was investigated using the concept of linear regression. Three different relationships were reviewed for their slopes to be used as the bias correction factor: Relationship 1 (R1) is based on the conventional linear regression, relationship 2 (R2) is forced to pass the origin and relationship 3 (R3) is the line whose slope is the G/R ratio. In other words, R1 is the regression line connecting the intercept and the mass centre of measurement pairs, R2 is the regression line forced to pass the origin, and R3 is the line connecting the origin and the mass centre. The slopes of all three relationships were reviewed analytically to compare them, and thereby, the effect of zero measurements could be evaluated. Additionally, the effect of using switched independent and dependent variables on the derived slopes was also evaluated. The theoretically derived results were then verified by analysing the rainfall event on 10–11 August 2010 in Korea. Finally, the difference between the bias‐corrected radar rain rate and the rain gauge rain rate was quantified by root mean square error and mean error so that it could be used as a measure for the evaluation of bias correction factors. In conclusion, the slope of R2 was found to be the best for the bias correction factor. However, when deciding the slope of this R2, the radar rain rate should be used as the independent variable in the low rain rate region, and the rain gauge rain rate in the high rain rate region above a certain threshold. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

10.
Flash-flood warning models can save lives and protect various kinds of infrastructure. In dry climate regions, rainfall is highly variable and can be of high-intensity. Since rain gauge networks in such areas are sparse, rainfall information derived from weather radar systems can provide useful input for flash-flood models. This paper presents a flash-flood warning model which utilizes radar rainfall data and applies it to two catchments that drain into the dry Dead Sea region. Radar-based quantitative precipitation estimates (QPEs) were derived using a rain gauge adjustment approach, either on a daily basis (allowing the adjustment factor to change over time, assuming available real-time gauge data) or using a constant factor value (derived from rain gauge data) over the entire period of the analysis. The QPEs served as input for a continuous hydrological model that represents the main hydrological processes in the region, namely infiltration, flow routing and transmission losses. The infiltration function is applied in a distributed mode while the routing and transmission loss functions are applied in a lumped mode. Model parameters were found by calibration based on the 5 years of data for one of the catchments. Validation was performed for a subsequent 5-year period for the same catchment and then for an entire 10-year record for the second catchment. The probability of detection and false alarm rates for the validation cases were reasonable. Probabilistic flash-flood prediction is presented applying Monte Carlo simulations with an uncertainty range for the QPEs and model parameters. With low probability thresholds, one can maintain more than 70% detection with no more than 30% false alarms. The study demonstrates that a flash-flood warning model is feasible for catchments in the area studied.  相似文献   

11.
This study evaluated four possible cases of comparing radar and rain gauge rain rate for the detection of mean‐field bias. These four cases, or detection designs, consider in this study are: (1) design 1‐uses all the data sets available, including zero radar rain rate and zero rain gauge rain rate, (2) design 2—uses the data sets of positive radar rain rate and zero or positive rain gauge rain rate, (3) design 3—uses the data sets of zero or positive radar rain rate and positive rain gauge rain rate and (4) design 4—uses the data sets of positive radar rain rate and positive rain gauge rain rate. A theoretical review of these four detection designs showed that only the design 1 causes no design bias, but designs 2, 3 and 4 can cause positive, negative and negative design biases, respectively. This theoretical result was also verified by applying these four designs to the rain rate field generated by a multi‐dimensional rain rate model, as well as to that of the Mt Gwanak radar in Korea. The results from both applications showed that especially the design 4, which is generally used for the detection of mean‐field bias of radar rain rate, causes a serious design bias; therefore, is inappropriate as a design for detecting the mean‐field bias of radar rain rate. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

12.
The final stage in processing radar data so as to arrive at an estimated rain field typically involves a comparison of the preliminary radar-derived estimates of hourly rainfall with those observed by ground-based gauges. Often a mean field bias adjustment will then be applied using an age-weighted average of the individual gauge–radar comparisons. In this paper, a mean field bias adjustment is presented that uses the path-integrated rainfall estimates provided by microwave links together with information from gauges. It is shown to be at least as efficient as the current gauge-based procedure used by the UK Met Office to improve the accuracy of radar-based estimates of rainfall at the ground.  相似文献   

13.
Uncertainty analysis of radar rainfall enables stakeholders and users have a clear knowledge of the possible uncertainty associated with the rainfall products. Long-term empirical modeling of the relationship between radar and gauge measurements is an efficient and practical method to describe the radar rainfall uncertainty. However, complicated variation of synoptic conditions makes the radar-rainfall uncertainty model based on historical data hard to extend in the future state. A promising solution is to integrate synoptic regimes with the empirical model and explore the impact of individual synoptic regimes on radar rainfall uncertainty. This study is an attempt to introduce season, one of the most important synoptic factor, into the radar rainfall uncertainty model and proposes a seasonal ensemble generator for radar rainfall using copula and autoregressive model. We firstly analyze the histograms of rainfall-weighted temperature, the radar-gauge relationships, and Box and Whisker plots in different seasons and conclude that the radar rainfall uncertainty has strong seasonal dependence. Then a seasonal ensemble generator is designed and implemented in a UK catchment under a temperate maritime climate, which can fully model marginal distribution, spatial dependence, temporal dependence and seasonal dependence of radar rainfall uncertainty. To test its performance, 12 typical rainfall events (4 for each season) are chosen to generate ensemble rainfall values. In each time step, 500 ensemble members are produced and the values of 5th to 95th percentiles are used to derive the uncertainty bands. Except several outliers, the uncertainty bands encompass the observed gauge rainfall quite well. The parameters of the ensemble generator vary considerably for each season, indicating the seasonal ensemble generator reflects the impact of seasons on radar rainfall uncertainty. This study is an attempt to simultaneously consider four key features of radar rainfall uncertainty and future study will investigate their impacts on the outputs of hydrological models with radar rainfall as input or initial conditions.  相似文献   

14.
司伟  包为民  瞿思敏  石朋 《湖泊科学》2018,30(2):533-541
空间集总式水文模型的洪水预报精度会受到面平均雨量估计误差的严重影响.点雨量监测值的误差类型、误差大小以及流域的雨量站点密度和站点的空间分布都会影响到面平均雨量的计算.为提高实时洪水预报精度,本文提出了一种基于降雨系统响应曲线洪水预报误差修正方法.通过此方法估计降雨输入项的误差,从而提高洪水预报精度.此方法将水文模型做为输入和输出之间的响应系统,用实测流量和计算流量之间的差值做为信息,通过降雨系统响应曲线,使用最小二乘估计原理,对面平均雨量进行修正,再用修正后的面平均雨量重新计算出流过程.将此修正方法结合新安江模型使用理想案例进行检验,并应用于王家坝流域的16场历史洪水以及此流域不同雨量站密度的情况下,结果证明均有明显修正效果,且在雨量站密度较低时修正效果更加明显.该方法是一种结构简单且不增加模型参数和复杂度的实时洪水修正的新方法.  相似文献   

15.
A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free‐falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of 1 m × 1 m, the expected number of received free‐falling raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re‐detachment amount. The re‐detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free‐falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re‐detachment amount were small parts of the total splash amount. Their proportions were 0·15% and 2·6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil‐splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

16.
Rainfall is a phenomenon difficult to model and predict, for the strong spatial and temporal heterogeneity and the presence of many zero values. We deal with hourly rainfall data provided by rain gauges, sparsely distributed on the ground, and radar data available on a fine grid of pixels. Radar data overcome the problem of sparseness of the rain gauge network, but are not reliable for the assessment of rain amounts. In this work we investigate how to calibrate radar measurements via rain gauge data and make spatial predictions for hourly rainfall, by means of Monte Carlo Markov Chain algorithms in a Bayesian hierarchical framework. We use zero-inflated distributions for taking zero-measurements into account. Several models are compared both in terms of data fitting and predictive performances on a set of validation sites. Finally, rainfall fields are reconstructed and standard error estimates at each prediction site are shown via easy-to-read spatial maps.  相似文献   

17.
ABSTRACT

Multisource rainfall products can be used to overcome the absence of gauged precipitation data for hydrological applications. This study aims to evaluate rainfall estimates from the Chinese S-band weather radar (CINRAD-SA), operational raingauges, multiple satellites (CMORPH, ERA-Interim, GPM, TRMM-3B42RT) and the merged satellite–gauge rainfall products, CMORPH-GC, as inputs to a calibrated probability distribution model (PDM) on the Qinhuai River Basin in Nanjing, China. The Qinhuai is a middle-sized catchment with an area of 799 km2. All sources used in this study are capable of recording rainfall at high spatial and temporal resolution (3 h). The discrepancies between satellite and radar data are analysed by statistical comparison with raingauge data. The streamflow simulation results from three flood events suggest that rainfall estimates using CMORPH-GC, TRMM-3B42RT and S-band radar are more accurate than those using the other rainfall sources. These findings indicate the potential to use satellite and radar data as alternatives to raingauge data in hydrological applications for ungauged or poorly gauged basins.  相似文献   

18.
ABSTRACT

Assessment of forecast precipitation is required before it can be used as input to hydrological models. Using radar observations in southeastern Australia, forecast rainfall from the Australian Community Climate Earth-System Simulator (ACCESS) was evaluated for 2010 and 2011. Radar rain intensities were first calibrated to gauge rainfall data from four research rainfall stations at hourly time steps. It is shown that the Australian ACCESS model (ACCESS-A) overestimated rainfall in low precipitation areas and underestimated elevated accumulations in high rainfall areas. The forecast errors were found to be dependent on the rainfall magnitude. Since the cumulative rainfall observations varied across the area and through the year, the relative error (RE) in the forecasts varied considerably with space and time, such that there was no consistent bias across the study area. Moreover, further analysis indicated that both location and magnitude errors were the main sources of forecast uncertainties on hourly accumulations, while magnitude was the dominant error on the daily time scale. Consequently, the precipitation output from ACCESS-A may not be useful for direct application in hydrological modelling, and pre-processing approaches such as bias correction or exceedance probability correction will likely be necessary for application of the numerical weather prediction (NWP) outputs.
EDITOR M.C. Acreman ASSOCIATE EDITOR A. Viglione  相似文献   

19.
The radar reflectivity (Z)–rain intensity (R) relationship fluctuates in both temporal and spatial scales. The dynamic factor analysis (DFA) and min/max autocorrelation factor analysis (MAFA) was specifically designed for considering various space–time integrations of gauge rainfall and radar reflectivity. We detect representative radar reflectivity observed around rainfall stations that were most responsible for rainfall intensity and identify the crucial patterns of the radar reflectivity in the Kaoping River watershed during Typhoon Morakot. Result shows that the MAFA and DFA can reduce the uncertainty of the dynamic Z‐R relationship effectively. The MAFA separates an entire area into two subareas (southern and northern areas) according to the relationships between the radar reflectivity and min/max autocorrelation factor (MAF) axes. For both areas, the different extents of temporal rainfall correlated with the radar reflectivity were determined using DFA. Especially in the northern area, the radar reflectivity was significantly related to the rainfall intensity for most stations without mountain blockage. Mountain blockages associated with the presence of terrain and wind direction were inferred the major factors that affected the relationship between radar reflectivity and rainfall intensity in the mountainous watershed. Further study can consider the terrain effect and meteorological information, such as wind speed and direction in the DFA model, with the dominant radar reflectivity to estimate the temporal rainfall patterns. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
Precipitation is a key control on watershed hydrologic modelling output, with errors in rainfall propagating through subsequent stages of water quantity and quality analysis. Most watershed models incorporate precipitation data from rain gauges; higher‐resolution data sources are available, but they are associated with greater computational requirements and expertise. Here, we investigate whether the Multisensor Precipitation Estimator (MPE or Stage IV Next‐Generation Radar) data improve the accuracy of streamflow simulations using the Soil and Water Assessment Tool (SWAT), compared with rain gauge data. Simulated flows from 2002 to 2010 at five timesteps were compared with observed flows for four nested subwatersheds of the Neuse River basin in North Carolina (21‐, 203‐, 2979‐, and 10 100‐km2 watershed area), using a multi‐objective function, informal likelihood‐weighted calibration approach. Across watersheds and timesteps, total gauge precipitation was greater than radar precipitation, but radar data showed a conditional bias of higher rainfall estimates during large events (>25–50 mm/day). Model parameterization differed between calibrations with the two datasets, despite the fact that all watershed characteristics were the same across simulation scenarios. This underscores the importance of linking calibration parameters to realistic processes. SWAT simulations with both datasets underestimated median and low flows, whereas radar‐based simulations were more accurate than gauge‐based simulations for high flows. At coarser timesteps, differences were less pronounced. Our results suggest that modelling efforts in watersheds with poor rain gauge coverage can be improved with MPE radar data, especially at short timesteps. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.  相似文献   

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

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