首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Rainfall data in continuous space provide an essential input for most hydrological and water resources planning studies. Spatial distribution of rainfall is usually estimated using ground‐based point rainfall data from sparsely positioned rain‐gauge stations in a rain‐gauge network. Kriging has become a widely used interpolation method to estimate the spatial distribution of climate variables including rainfall. The objective of this study is to evaluate three geostatistical (ordinary kriging [OK], ordinary cokriging [OCK], kriging with an external drift [KED]), and two deterministic (inverse distance weighting, radial basis function) interpolation methods for enhanced spatial interpolation of monthly rainfall in the Middle Yarra River catchment and the Ovens River catchment in Victoria, Australia. Historical rainfall records from existing rain‐gauge stations of the catchments during 1980–2012 period are used for the analysis. A digital elevation model of each catchment is used as the supplementary information in addition to rainfall for the OCK and kriging with an external drift methods. The prediction performance of the adopted interpolation methods is assessed through cross‐validation. Results indicate that the geostatistical methods outperform the deterministic methods for spatial interpolation of rainfall. Results also indicate that among the geostatistical methods, the OCK method is found to be the best interpolator for estimating spatial rainfall distribution in both the catchments with the lowest prediction error between the observed and estimated monthly rainfall. Thus, this study demonstrates that the use of elevation as an auxiliary variable in addition to rainfall data in the geostatistical framework can significantly enhance the estimation of rainfall over a catchment.  相似文献   

2.
Compared to other estimation techniques, one advantage of geostatistical techniques is that they provide an index of the estimation accuracy of the variable of interest with the kriging estimation standard deviation (ESD). In the context of radar–raingauge quantitative precipitation estimation (QPE), we address in this article the question of how the kriging ESD can be transformed into a local spread of error by using the dependency of radar errors to the rain amount analyzed in previous work. The proposed approach is implemented for the most significant rain events observed in 2008 in the Cévennes-Vivarais region, France, by considering both the kriging with external drift (KED) and the ordinary kriging (OK) methods. A two-step procedure is implemented for estimating the rain estimation accuracy: (i) first kriging normalized ESDs are computed by using normalized variograms (sill equal to 1) to account for the observation system configuration and the spatial structure of the variable of interest (rainfall amount, residuals to the drift); (ii) based on the assumption of a linear relationship between the standard deviation and the mean of the variable of interest, a denormalization of the kriging ESDs is performed globally for a given rain event by using a cross-validation procedure. Despite the fact that the KED normalized ESDs are usually greater than the OK ones (due to an additional constraint in the kriging system and a weaker spatial structure of the residuals to the drift), the KED denormalized ESDs are generally smaller the OK ones, a result consistent with the better performance observed for the KED technique. The evolution of the mean and the standard deviation of the rainfall-scaled ESDs over a range of spatial (5–300 km2) and temporal (1–6 h) scales demonstrates that there is clear added value of the radar with respect to the raingauge network for the shortest scales, which are those of interest for flash-flood prediction in the considered region.  相似文献   

3.
Using auxiliary information to improve the prediction accuracy of soil properties in a physically meaningful and technically efficient manner has been widely recognized in pedometrics. In this paper, we explored a novel technique to effectively integrate sampling data and auxiliary environmental information, including continuous and categorical variables, within the framework of the Bayesian maximum entropy (BME) theory. Soil samples and observed auxiliary variables were combined to generate probability distributions of the predicted soil variable at unsampled points. These probability distributions served as soft data of the BME theory at the unsampled locations, and, together with the hard data (sample points) were used in spatial BME prediction. To gain practical insight, the proposed approach was implemented in a real-world case study involving a dataset of soil total nitrogen (TN) contents in the Shayang County of the Hubei Province (China). Five terrain indices, soil types, and soil texture were used as auxiliary variables to generate soft data. Spatial distribution of soil total nitrogen was predicted by BME, regression kriging (RK) with auxiliary variables, and ordinary kriging (OK). The results of the prediction techniques were compared in terms of the Pearson correlation coefficient (r), mean error (ME), and root mean squared error (RMSE). These results showed that the BME predictions were less biased and more accurate than those of the kriging techniques. In sum, the present work extended the BME approach to implement certain kinds of auxiliary information in a rigorous and efficient manner. Our findings showed that the BME prediction technique involving the transformation of variables into soft data can improve prediction accuracy considerably, compared to other techniques currently in use, like RK and OK.  相似文献   

4.
Abstract

Radar quantitative precipitation estimates (QPEs) were assessed using reference values established by means of a geostatistical approach. The reference values were estimated from raingauge data using the block kriging technique, and the reference meshes were selected on the basis of the kriging estimation variance. Agreement between radar QPEs and reference rain amounts was shown to increase slightly with the space–time scales. The statistical distributions of the errors were modelled conditionally with respect to several factors using the GAMLSS approach. The conditional bias of the errors presents a complex structure that depends on the space–time scales and the considered geographical sub-domains, while the standard deviation of the errors has a more homogeneous behaviour. The estimation standard deviation of the reference rainfall and the standard deviation of the errors between radar and reference rainfall were found to have the same magnitude, indicating the limitations of the available network in terms of providing accurate reference values for the spatial scales considered (5–100 km2).
Editor D. Koutsoyiannis; Guest editor R.J. Moore

Citation Delrieu, G., Bonnifait, L., Kirstetter, P.-E., and Boudevillain, B., 2013. Dependence of radar quantitative precipitation estimation error on the rain intensity in the Cévennes region, France. Hydrological Sciences Journal, 59 (7), 1300–1311. http://dx.doi.org/10.1080/02626667.2013.827337  相似文献   

5.
In this study, we link and compare the geographically weighted regression (GWR) model with the kriging with an external drift (KED) model of geostatistics. This includes empirical work where models are performance tested with respect to prediction and prediction uncertainty accuracy. In basic forms, GWR and KED (specified with local neighbourhoods) both cater for nonstationary correlations (i.e. the process is heteroskedastic with respect to relationships between the variable of interest and its covariates) and as such, can predict more accurately than models that do not. Furthermore, on specification of an additional heteroskedastic term to the same models (now with respect to a process variance), locally-accurate measures of prediction uncertainty can result. These heteroskedastic extensions of GWR and KED can be preferred to basic constructions, whose measures of prediction uncertainty are only ever likely to be globally-accurate. We evaluate both basic and heteroskedastic GWR and KED models using a case study data set, where data relationships are known to vary across space. Here GWR performs well with respect to the more involved KED model and as such, GWR is considered a viable alternative to the more established model in this particular comparison. Our study adds to a growing body of empirical evidence that GWR can be a worthy predictor; complementing its more usual guise as an exploratory technique for investigating relationships in multivariate spatial data sets.  相似文献   

6.
The spatial variability of precipitation has often been a topic of research, since accurate modelling of precipitation is a crucial condition for obtaining reliable results in hydrology and geomorphology. In mountainous areas, the sparsity of the measurement networks makes an accurate and reliable spatialization of rainfall amounts at the local scale difficult. The purpose of this paper is to show how the use of a digital elevation model can improve interpolation processes at the subregional scale for mapping the mean annual and monthly precipitation from rainfall observations (40 years) recorded in a region of 1400 km2 in southern Italy. Besides linear regression of precipitation against elevation, two methods of interpolation are applied: inverse squared distance and ordinary cokriging. Cross‐validation indicates that the inverse distance interpolation, which ignores the information on elevation, yields the largest prediction errors. Smaller prediction errors are produced by linear regression and ordinary cokriging. However, the results seem to favour the multivariate geostatistical method including auxiliary information (related to elevation). We conclude that ordinary cokriging is a very flexible and robust interpolation method because it can take into account several properties of the landscape; it should therefore be applicable in other mountainous regions, especially where precipitation is an important geomorphological factor. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

7.
《水文科学杂志》2013,58(5):1038-1050
Abstract

Universal kriging is applied to the command area of a set of canal irrigation projects in north-western India to show its applicability for optimal contour mapping of groundwater levels. This command area faces the problem of a rising groundwater table and manifestation of waterlogging over the years at many places. With the use of measured elevations of the water table in September 1990 at 143 observation sites in an area of 4500 km2, an omni-directional experimental semivariogram was constructed. The concave upward shape of omni-directional experimental semivariogram indicated the non-stationarity of the data and so the need for universal kriging. Directional semivariograms were calculated to find out the direction along which there is least drift. This directional semivariogram was considered as the underlying semivariogram and various theoretical semivariogram models were fitted to it. The drift order was estimated from the cross-validation procedure. The model and drift order finally selected were used to estimate the groundwater levels and corresponding estimation variances at the nodes of a square grid of 2 km × 2 km, and to develop the corresponding contour map.  相似文献   

8.
Abstract

The present research study investigates the application of nonlinear normalizing data transformations in conjunction with ordinary kriging (OK) for the accurate prediction of groundwater level spatial variability in a sparsely-gauged basin. We investigate three established normalizing methods, Gaussian anamorphosis, trans-Gaussian kriging and the Box-Cox method to improve the estimation accuracy. The first two are applied for the first time to groundwater level data. All three methods improve the mean absolute prediction error compared to the application of OK to the non-transformed data. In addition, a modified Box-Cox transformation is proposed and applied to normalize the hydraulic heads. The modified Box-Cox transformation in conjunction with OK is found to be the optimal spatial model based on leave-one-out cross-validation. The recently established Spartan semivariogram family provides the optimal model fit to the transformed data. Finally, we present maps of the groundwater level and the kriging variance based on the optimal spatial model.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Varouchakis, E.A., Hristopoulos, D.T., and Karatzas, G.P., 2012. Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application. Hydrological Sciences Journal, 57 (7), 1404–1419.  相似文献   

9.
In humid, well-vegetated areas, such as in the northeastern US, runoff is most commonly generated from relatively small portions of the landscape becoming completely saturated, however, little is known about the spatial and temporal behavior of these saturated regions. Indicator kriging provides a way to use traditional water table data to quantify probability of saturation to evaluate predicted spatial distributions of runoff generation risk, especially for the new generation of water quality models incorporating saturation excess runoff theory. When spatial measurements of a variable are transformed to binary indicators (i.e., 1 if above a given threshold value and 0 if below) and the resulting indicator semivariogram is modeled, indicator kriging produces the probability of the measured variable to exceed the threshold value. Indicator kriging gives quantified probability of saturation or, consistent with saturation excess runoff theory, runoff generation risk with depth to water table as the variable and the threshold set near the soil surface. The probability of saturation for a 120 m × 180 m hillslope based upon 43 measurements of depth to water table is investigated with indicator semivariograms for six storm events. The indicator semivariograms show high spatial structure in saturated regions with large antecedent rainfall conditions. The temporal structure of the data is used to generate interpolated (soft) data to supplement measured (hard) data. This improved the spatial structure of the indicator semivariograms for lower antecedent rainfall conditions. Probability of saturation was evaluated through indicator kriging incorporating soft data showing, based on this preliminary study, highly connected regions of saturation as expected for the wet season (April through May) in the Catskill Mountain region of New York State. Supplementation of hard data with soft data incorporates physical hydrology of the hillslope to capture significant patterns not available when using hard data alone for indicator kriging. With the need for water quality models incorporating appropriate runoff generation risk estimates on the rise, this manner of data will lay the groundwork for future model evaluation and development.  相似文献   

10.
Abstract

This paper refers to the quantification and prediction of the sedimentation rate of 26 hillside-dam reservoirs in Central Tunisia. The objectives of the study are to develop a simple and practical methodology to identify controlling factors of sedimentation, and to propose a regionalization from the study sites. Principal component analysis (PCA) and complementary multi-dimensional statistical methods are used to relate highly variable area-specific sediment yield to hydro-morphometric, lithological, geomorphological and anthropogenic characteristics of catchments. It appears that catchment area is not the main controlling factor of sedimentation in the studied area. The overall slope index, drainage network characteristics and runoff parameters are also important in characterizing sediment yield. Applied to the annual sedimentation rate series, PCA resulted in retaining the first three principal axes, explaining 65% of the total variance. Statistical methods showed that the overall slope index, the total drainage length, the compacity index and the runoff parameters are as important for the sedimentation quantification. This allowed a graphical clustering of the study zone into three distinct groups having similar behaviours: (i) watersheds characterized by high sediment transport rates and high runoff coefficients, (ii) basins distinguished by relatively low values of both flow discharge and sediment transport rates, and (iii) watersheds with an intermediate sediment yield, especially characterized by relatively high relief. In a second step, a multiple regression model including the four characteristic catchment properties was developed, presenting a valuable tool to predict area-specific sediment yield from catchments in central Tunisia. This model shows reasonable efficiency with an absolute prediction error of 81%.

Citation Ayadi, I., Abida, H., Djebbar, Y. & Mahjoub, M. R. (2010) Sediment yield variability in central Tunisia: a quantitative analysis of its controlling factors. Hydrol. Sci. J. 55(3), 446–458.  相似文献   

11.
Abstract

A method is introduced for probabilistic forecasting of hydrological events based on geostatistical analysis. In this method, the predictors of a hydrological variable define a virtual field such that, in this field, the observed dependent variables are considered as measurement points. Variography of the measurement points enables the use of the system of kriging equations to estimate the value of the variable at non-measured locations of the field. Non-measured points are the forecasts associated with specific predictors. Calculation of the estimation variance facilitates probabilistic analysis of the forecast variables. The method is applied to case studies of the Red River in Manitoba, Canada and Karoon River in Khoozestan, Iran. The study analyses the advantages and limitations of the proposed method in comparison with a K-nearest neighbour approach and linear and nonlinear multiple regression. The utility of the proposed method for forecasting hydrological variables with a conditional probability distribution is demonstrated.  相似文献   

12.
ABSTRACT

Predicting the impacts of climate change on water resources remains a challenging task and requires a good understanding of the dynamics of the forcing terms in the past. In this study, the variability of precipitation and drought patterns is studied over the Mediterranean catchment of the Medjerda in Tunisia based on an observed rainfall dataset collected at 41 raingauges during the period 1973–2012. The standardized precipitation index and the aridity index were used to characterize drought variability. Multivariate and geostatistical techniques were further employed to identify the spatial variability of annual rainfall. The results show that the Medjerda is marked by a significant spatio-temporal variability of drought, with varying extreme wet and dry events. Four regions with distinct rainfall regimes are identified by utilizing the K-means cluster analysis. A principal component analysis identifies the variables that are responsible for the relationships between precipitation and drought variability.  相似文献   

13.
Abstract

A novel approach is presented for combining spatial and temporal detail from newly available TRMM-based data sets to derive hourly rainfall intensities at 1-km spatial resolution for hydrological modelling applications. Time series of rainfall intensities derived from 3-hourly 0.25° TRMM 3B42 data are merged with a 1-km gridded rainfall climatology based on TRMM 2B31 data to account for the sub-grid spatial distribution of rainfall intensities within coarse-scale 0.25° grid cells. The method is implemented for two dryland catchments in Tunisia and Senegal, and validated against gauge data. The outcomes of the validation show that the spatially disaggregated and intensity corrected TRMM time series more closely approximate ground-based measurements than non-corrected data. The method introduced here enables the generation of rainfall intensity time series with realistic temporal and spatial detail for dynamic modelling of runoff and infiltration processes that are especially important to water resource management in arid regions.

Editor D. Koutsoyiannis

Citation Tarnavsky, E., Mulligan, M. and Husak, G., 2012. Spatial disaggregation and intensity correction of TRMM-based rainfall time series for hydrological applications in dryland catchments. Hydrological Sciences Journal, 57 (2), 248–264.  相似文献   

14.
Abstract

The scale invariance of rainfall series in the Tunis area, Tunisia (semi-arid Mediterranean climate) is studied in a mono-fractal framework by applying the box counting method to four series of observations, each about 2.5 years in length, based on a time resolution of 5 min. In addition, a single series of daily rainfall records for the period 1873–2009 was analysed. Three self-similar structures were identified: micro-scale (5 min to 2 d) with fractal dimension 0.44, meso-scale (2 d to one week) and synoptic-scale (one week to eight months) with fractal dimension 0.9. Interpretation of these findings suggests that only the micro-scale and transition to saturation are consistent, while the high fractal dimension relating to the synoptic scale might be affected by the tendency to saturation. A sensitivity analysis of the estimated fractal dimension was performed using daily rainfall data by varying the series length, as well as the intensity threshold for the detection of rain.

Editor Z.W. Kundzewicz; Associate editor S. Grimaldi

Citation Ghanmi, H., Bargaoui, Z., and Mallet, C., 2013. Investigation of the fractal dimension of rainfall occurrence in a semi-arid Mediterranean climate. Hydrological Sciences Journal, 58 (3), 483–497.  相似文献   

15.
 This paper deals with the problem of spatial data mapping. A new method based on wavelet interpolation and geostatistical prediction (kriging) is proposed. The method – wavelet analysis residual kriging (WARK) – is developed in order to assess the problems rising for highly variable data in presence of spatial trends. In these cases stationary prediction models have very limited application. Wavelet analysis is used to model large-scale structures and kriging of the remaining residuals focuses on small-scale peculiarities. WARK is able to model spatial pattern which features multiscale structure. In the present work WARK is applied to the rainfall data and the results of validation are compared with the ones obtained from neural network residual kriging (NNRK). NNRK is also a residual-based method, which uses artificial neural network to model large-scale non-linear trends. The comparison of the results demonstrates the high quality performance of WARK in predicting hot spots, reproducing global statistical characteristics of the distribution and spatial correlation structure.  相似文献   

16.
Abstract

Abstract The construction of the Gabcikovo hydropower plant and the diversion of the Danube River over 25 km into an artificial channel in 1992 influenced the groundwater regime of the region considerably. Statistical and geostatistical methods are used to quantify changes of different groundwater characteristics on the Hungarian side of the river based on observations in the time period 1960–2000. External drift kriging was used to interpolate groundwater levels and the other related variables. While mean groundwater levels did not change appreciably, there are significant changes in the variability. Standard deviations of the groundwater levels and the amplitude of the annual cycle decreased near the old river bed of the Danube. The water-level fluctuations of the Danube have a decreased influence on the groundwater dynamics. Interrelationships between water levels in wells have also changed.  相似文献   

17.
Abstract

The generation of reliable quantitative precipitation estimations (QPEs) through use of raingauge and radar data is an important issue. This study investigates the impacts of radar QPEs with different densities of raingauge networks on rainfall–runoff processes through a semi-distributed parallel-type linear reservoir rainfall–runoff model. The spatial variation structures of the radar QPE, raingauge QPE and radar-gauge residuals are examined to review the current raingauge network, and a compact raingauge network is identified via the kriging method. An analysis of the large-scale spatial characteristics for use with a hydrological model is applied to investigate the impacts of a raingauge network coupled with radar QPEs on the modelled rainfall–runoff processes. Since the precision in locating the storm centre generally represents how well the large-scale variability is reproduced; the results show not only the contribution of kriging to identify a compact network coupled with radar QPE, but also that spatial characteristics of rainfalls do affect the hydrographs.
Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Pan, T.-Y., Li, M.-Y., Lin, Y.-J., Chang, T.-J., Lai, J.-S., and Tan, Y.-C., 2014. Sensitivity analysis of the hydrological response of the Gaping River basin to radar-raingauge quantitative precipitation estimates. Hydrological Sciences Journal, 59 (7), 1335–1352. http://dx.doi.org/10.1080/02626667.2014.923969  相似文献   

18.
Abstract

Modelling and prediction of hydrological processes (e.g. rainfall–runoff) can be influenced by discontinuities in observed data, and one particular case may arise when the time scale (i.e. resolution) is coarse (e.g. monthly). This study investigates the application of catastrophe theory to examine its suitability to identify possible discontinuities in the rainfall–runoff process. A stochastic cusp catastrophe model is used to study possible discontinuities in the monthly rainfall–runoff process at the Aji River basin in Azerbaijan, Iran. Monthly-averaged rainfall and flow data observed over a period of 20 years (1981–2000) are analysed using the Cuspfit program. In this model, rainfall serves as a control variable and runoff as a behavioural variable. The performance of this model is evaluated using four measures: correlation coefficient, log-likelihood, Akaike information criterion (AIC) and Bayesian information criterion (BIC). The results indicate the presence of discontinuities in the rainfall–runoff process, with a significant sudden jump in flow (cusp signal) when rainfall reaches a threshold value. The performance of the model is also found to be better than that of linear and logistic models. The present results, though preliminary, are promising in the sense that catastrophe theory can play a possible role in the study of hydrological systems and processes, especially when the data are noisy.

Citation Ghorbani, M. A., Khatibi, R., Sivakumar, B. & Cobb, L. (2010) Study of discontinuities in hydrological data using catastrophe theory. Hydrol. Sci. J. 55(7), 1137–1151.  相似文献   

19.
Abstract

This study examines relationships between model parameters and urbanization variables for evaluating urbanization effects in a watershed. Rainfall–runoff simulation using the Nash model is the main basis of the study. Mean rainfall and excesses resulting from time-variant losses were completed using the kriging and nonlinear programming methods, respectively. Calibrated parameters of 47 events were related to urbanized variables, change of shape parameter responds more sensitively than that of scale parameter based on comparisons between annual average and optimal interval methods. Regression equations were used to obtain four continuous correlations for linking shape parameter with urbanization variables. Verification of 10 events demonstrates that shape parameter responds more strongly to imperviousness than to population, and a power relationship is suitable. Therefore, an imperviousness variable is a major reference for analysing urbanization changes to a watershed. This study found that time to peak of IUH was reduced from 11.76 to 3.97 h, whereas peak discharge increased from 44.79 to 74.92 m3/s.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi

Citation Huang, S.-Y., Cheng, S.-J., Wen, J.-C. and Lee, J.-H., 2012. Identifying hydrograph parameters and their relationships to urbanization variables. Hydrological Sciences Journal, 57 (1), 144–161.  相似文献   

20.
Scale properties of daily areal rainfall in south-eastern Norway are investigated, and the events are classified into small- and large-scale events by statistical pattern recognition. The coefficient of variation and the standard deviation, computed from all precipitation gauges for one event, are considered to be mutually independent and are assigned gamma densities. The two statistical parameters form a feature vector which is used as the discriminatory variable in the classification algorithm. Empirical and theoretical semivariograms were fitted to the classified sets of events, and the parameters of the semivariogram (nugget, sill and range) were investigated for each set. Different values for all parameters were found, and the differences in range, describing the scale of the phenomena, were further investigated by the use of cross validation with the MSIE (mean squared interpolation error) criterion. Optimal range of daily small-scale precipitation was found to be 16 km (19 km without nugget effect), and 33 km (70 km without nugget effect) for daily large-scale precipitation. These optimal values for range, however, proved to be of small consequence for the estimated areal precipitation using kriging interpolation. Areal precipitation was computed for catchments of size 50 km2 and 1200 km2, and insignificant differences in the computed areal precipitation were found when applying the class-specific semivariograms instead of the mean semivariogram calculated from all events regardless of class. However, when areal reduction factors (ARFs) were derived from the analytical expressions of the class-specific semivariograms, considerable differences were found for the ARF curves of the two precipitation types. The difference became more marked for precipitation events of low probability of occurrence.  相似文献   

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

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