首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 26 毫秒
1.
ABSTRACT

The non-parametric mathematical framework of bilinear surface smoothing (BSS) methodology provides flexible means for spatial (two dimensional) interpolation of variables. As presented in a companion paper, interpolation is accomplished by means of fitting consecutive bilinear surface into a regression model with known break points and adjustable smoothing terms defined by means of angles formed by those bilinear surface. Additionally, the second version of the methodology (BSSE) incorporates, in an objective manner, the influence of an explanatory variable available at a considerably denser dataset. In the present study, both versions are explored and illustrated using both synthesized and real world (hydrological) data, and practical aspects of their application are discussed. Also, comparison and validation against the results of commonly used spatial interpolation methods (inverse distance weighted, spline, ordinary kriging and ordinary cokriging) are performed in the context of the real world application. In every case, the method’s efficiency to perform interpolation between data points that are interrelated in a complicated manner was confirmed. Especially during the validation procedure presented in the real world case study, BSSE yielded very good results, outperforming those of the other interpolation methods. Given the simplicity of the approach, the proposed mathematical framework’s overall performance is quite satisfactory, indicating its applicability for diverse tasks of scientific and engineering hydrology and beyond.
Editor Z. W. Kundzewicz; Associate editor A. Carsteanu  相似文献   

2.
Abstract

Gridded meteorological data are available for all of Norway as time series dating from 1961. A new way of interpolating precipitation in space from observed values is proposed. Based on the criteria that interpolated precipitation fields in space should be consistent with observed spatial statistics, such as spatial mean, variance and intermittency, spatial fields of precipitation are simulated from a gamma distribution with parameters determined from observed data, adjusted for intermittency. The simulated data are distributed in space, using the spatial pattern derived from kriging. The proposed method is compared to indicator kriging and to the current methodology used for producing gridded precipitation data. Cross-validation gave similar results for the three methods with respect to RMSE, temporal mean and standard deviation, whereas a comparison on estimated spatial variance showed that the new method has a near perfect agreement with observations. Indicator kriging underestimated the spatial variance by 60–80% and the current method produced a significant scatter in its estimates.

Citation Skaugen, T. & Andersen, J. (2010) Simulated precipitation fields with variance-consistent interpolation. Hydrol. Sci. J. 55(5), 676–686.  相似文献   

3.
Abstract

New mathematical programming models are proposed, developed and evaluated in this study for estimating missing precipitation data. These models use nonlinear and mixed integer nonlinear mathematical programming (MINLP) formulations with binary variables. They overcome the limitations associated with spatial interpolation methods relevant to the arbitrary selection of weighting parameters, the number of control points within a neighbourhood, and the size of the neighbourhood itself. The formulations are solved using genetic algorithms. Daily precipitation data obtained from 15 rain gauging stations in a temperate climatic region are used to test and derive conclusions about the efficacy of these methods. The developed methods are compared with some naïve approaches, multiple linear regression, nonlinear least-square optimization, kriging, and global and local trend surface and thin-plate spline models. The results suggest that the proposed new mathematical programming formulations are superior to those obtained from all the other spatial interpolation methods tested in this study.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi

Citation Teegavarapu, R.S.V., 2012. Spatial interpolation using nonlinear mathematical programming models for estimation of missing precipitation records. Hydrological Sciences Journal, 57 (3), 383–406.  相似文献   

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

5.
Rainfall data are a fundamental input for effective planning, designing and operating of water resources projects. A well‐designed rain gauge network is capable of providing accurate estimates of necessary areal average and/or point rainfall estimates at any desired ungauged location in a catchment. Increasing network density with additional rain gauge stations has been the main underlying criterion in the past to reduce error and uncertainty in rainfall estimates. However, installing and operation of additional stations in a network involves large cost and manpower. Hence, the objective of this study is to design an optimal rain gauge network in the Middle Yarra River catchment in Victoria, Australia. The optimal positioning of additional stations as well as optimally relocating of existing redundant stations using the kriging‐based geostatistical approach was undertaken in this study. Reduction of kriging error was considered as an indicator for optimal spatial positioning of the stations. Daily rainfall records of 1997 (an El Niño year) and 2010 (a La Niña year) were used for the analysis. Ordinary kriging was applied for rainfall data interpolation to estimate the kriging error for the network. The results indicate that significant reduction in the kriging error can be achieved by the optimal spatial positioning of the additional as well as redundant stations. Thus, the obtained optimal rain gauge network is expected to be appropriate for providing high quality rainfall estimates over the catchment. The concept proposed in this study for optimal rain gauge network design through combined use of additional and redundant stations together is equally applicable to any other catchment. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.  相似文献   

6.
Abstract

This paper compares the performance of three geostatistical algorithms, which integrate elevation as an auxiliary variable: kriging with external drift (KED); kriging combined with regression, called regression kriging (RK) or kriging after detrending; and co-kriging (CK). These three methods differ by the way by in which the secondary information is introduced into the prediction procedure. They are applied to improve the prediction of the monthly average rainfall observations measured at 106 climatic stations in Tunisia over an area of 164 150 km2 using the elevation as the auxiliary variable. The experimental sample semivariograms, residual semivariograms and cross-variograms are constructed and fitted to estimate the rainfall levels and the estimation variance at the nodes of a square grid of 20 km?×?20 km resolution and to develop corresponding contour maps. Contour diagrams for KED and RK were similar and exhibited a pattern corresponding more closely to local topographic features when (a) the network is sparse and (b) the rainfall–elevation correlation is poor, while CK showed a smooth zonal pattern. Smaller prediction variances are obtained for the RK algorithm. The cross-validation showed that the RMSE obtained for CK gave better results than for KED or RK.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Feki, H., Slimani, M., and Cudennec, C., 2012. Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods. Hydrological Sciences Journal, 57 (7), 1294–1314.  相似文献   

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

8.
Abstract

Because the properties of eroded soil affect the deposition phenomena and transport capacity of chemical materials by eroded particles, recent research is trying to link the grain-size distribution of the eroded sediment to that of the original soil in order to explain the enrichment of chemical content of the sediment with the respect to the parent soil. In this study, the spatial distribution of nitrogen, phosphorus and total organic carbon was firstly deduced using the measurements carried out in 47 soil samples distributed over a forested basin together with a kriging interpolation method. Then the load of each chemical was calculated at morphological unit and basin scales using the above-mentioned spatial distributions and sediment yield values calculated by the SEDD (SEdiment Delivery Distributed) model, which couples the universal soil loss equation with a spatial disaggregation criterion of sediment delivery processes. Finally, at basin scale, a new expression of the enrichment ratio of a given chemical was applied.  相似文献   

9.
In this technical note, we investigate the hypothesis that ‘non-linearity matters in the spatial mapping of complex patterns of groundwater arsenic contamination’. The spatial mapping pertained to data-driven techniques of spatial interpolation based on sampling data at finite locations. Using the well known example of extensive groundwater contamination by arsenic in Bangladesh, we find that the use of a highly non-linear pattern learning technique in the form of an artificial neural network (ANN) can yield more accurate results under the same set of constraints when compared to the ordinary kriging method. One ANN and a variogram model were used to represent the spatial structure of arsenic contamination for the whole country. The probability for successful detection of a well as safe or unsafe was found to be atleast 15% larger than that by kriging under the country-wide scenario. The probability of false hopes, which is a serious issue in public health monitoring was found to be significantly lower (by more than 10%) than that by kriging.  相似文献   

10.
ABSTRACT

Flow–duration curves (FDCs) are essential to support decisions on water resources management, and their regionalization is fundamental for the assessment of ungauged basins. In comparison with calibrated rainfall–runoff models, statistical methods provide data-driven estimates representing a useful benchmark. The objective of this work is the interpolation of FDCs from ~500 discharge gauging stations in the Danube. To this aim we use total negative deviation top-kriging (TNDTK), as multi-regression models are shown to be unsuitable for representing FDCs across all durations and sites. TNDTK shows a high accuracy for the entire Danube region, with overall Nash-Sutcliffe efficiency values computed in a leave-p-out cross-validation scheme (p equal to one site, one-third and half of the sites), all above 0.88. A reliability measure based on kriging variance is attached to each interpolated FDC at ~4000 prediction nodes. The GIS layer of regionalized FDCs is made available for broader use in the region.  相似文献   

11.
Abstract

A case study is presented for the application of statistical and geostatistical methods to the problem of estimating groundwater quality variables. This methodology has been applied to the investigation of the detrital aquifer of the Bajo Andarax (Almería, Spain). The use of principal components analysis is proposed, as a first step, for identifying relevant types of groundwater and the processes that bring about a change in their quality. As a result of this application, three factors were obtained, which were used as three new variables (VI: sulphate influence; V2: thermal influence; and V3: marine influence). Analysis of their spatial distribution was performed through the calculation of experimental and theoretical variograms, which served as input for geostatistical modelling using ordinary block kriging. This analysis has allowed a probabilistic representation of the data to be obtained by mapping the three variables throughout the aquifer for each sampling point. In this way, one can evaluate the spatial and temporal variation of the principal physico-chemical processes associated with the three variables VI, V2 and V3 implicated in the groundwater quality of the detrital aquifer.  相似文献   

12.
Abstract

Methods were evaluated for interpolating precipitation (P), evapotranspiration (ET), and runoff (RO) at ungauged points on Shikoku Island, Japan, using data gathered from gauged stations on the same island. Two methods were examined: a “local” cubic spline interpolator, which, for a given point, fitted the function exactly to nearby gauged data points; and a “global” multivariate regression interpolator, which fitted the function to all gauged data points based on their topographic positions (i.e. latitude, longitude, altitude). Local and global interpolators did not generate similar results for P and temperature (T). The spatial density of gauged data points used in the interpolation affected the performance of the interpolators. With any given density of gauged data points included in the interpolation, the local interpolator outperformed the global interpolator. The findings indicate that local interpolators are more accurate predictors of the spatial distribution of water balance components in mountainous regions such as Shikoku Island.  相似文献   

13.
ABSTRACT

A new deep extreme learning machine (ELM) model is developed to predict water temperature and conductivity at a virtual monitoring station. Based on previous research, a modified ELM auto-encoder is developed to extract more robust invariance among the water quality data. A weighted ELM that takes seasonal variation as the basis of weighting is used to predict the actual value of water quality parameters at sites which only have historical data and no longer generate new data. The performance of the proposed model is validated against the monthly data from eight monitoring stations on the Zengwen River, Taiwan (2002–2017). Based on root mean square error, mean absolute error, mean absolute percentage error and correlation coefficient, the experimental results show that the new model is better than the other classical spatial interpolation methods.  相似文献   

14.
Spatial interpolation methods for nonstationary plume data   总被引:1,自引:0,他引:1  
Plume interpolation consists of estimating contaminant concentrations at unsampled locations using the available contaminant data surrounding those locations. The goal of ground water plume interpolation is to maximize the accuracy in estimating the spatial distribution of the contaminant plume given the data limitations associated with sparse monitoring networks with irregular geometries. Beyond data limitations, contaminant plume interpolation is a difficult task because contaminant concentration fields are highly heterogeneous, anisotropic, and nonstationary phenomena. This study provides a comprehensive performance analysis of six interpolation methods for scatter-point concentration data, ranging in complexity from intrinsic kriging based on intrinsic random function theory to a traditional implementation of inverse-distance weighting. High resolution simulation data of perchloroethylene (PCE) contamination in a highly heterogeneous alluvial aquifer were used to generate three test cases, which vary in the size and complexity of their contaminant plumes as well as the number of data available to support interpolation. Overall, the variability of PCE samples and preferential sampling controlled how well each of the interpolation schemes performed. Quantile kriging was the most robust of the interpolation methods, showing the least bias from both of these factors. This study provides guidance to practitioners balancing opposing theoretical perspectives, ease-of-implementation, and effectiveness when choosing a plume interpolation method.  相似文献   

15.
ABSTRACT

Several satellite-based precipitation estimates are becoming available at a global scale, providing new possibilities for water resources modelling, particularly in data-sparse regions and developing countries. This work provides a first validation of five different satellite-based precipitation products (TRMM-3B42 v6 and v7, RFE 2.0, PERSIANN-CDR, CMORPH1.0 version 0.x) in the 1785 km2 Makhazine catchment (Morocco). Precipitation products are first compared against ground observations. Ten raingauges and four different interpolation methods (inverse distance, nearest neighbour, ordinary kriging and residual kriging with altitude) were used to compute a set of interpolated precipitation reference fields. Second, a parsimonious conceptual hydrological model is considered, with a simulation approach based on the random generation of model parameters drawn from existing parameter set libraries, to compare the different precipitation inputs. The results indicate that (1) all four interpolation methods, except the nearest neighbour approach, give similar and valid precipitation estimates at the catchment scale; (2) among the different satellite-based precipitation estimates verified, the TRMM-3B42 v7 product is the closest to observed precipitation, and (3) despite poor performance at the daily time step when used in the hydrological model, TRMM-3B42 v7 estimates are found adequate to reproduce monthly dynamics of discharge in the catchment. The results provide valuable perspectives for water resources modelling of data-scarce catchments with satellite-based rainfall data in this region.
Editor M.C. Acreman; Associate editor N. Verhoest  相似文献   

16.
Abstract

Since droughts are natural phenomena, their occurrence cannot be predicted with certainty and thus it must be treated as a random variable. Once drought duration and magnitude have been found objectively, it is possible to plan for the transport of water in known quantities to drought-stricken areas either from alternative water resources or from water stored during wet periods. The summation of deficits over a particular period is referred to as the drought magnitude. Drought intensity is the ratio of drought magnitude to its duration. These drought properties at different truncation levels provide significant hydrological and hydrometeorological design quantities. In this study, the run analysis and z-score are used for determining drought properties of given hydrological series. In addition, kriging is used as a spatial drought analysis for mapping. This study is applied to precipitation records for Istanbul, Edirne, Tekirdag and Kirklareli in the Trakya region, Turkey and then the drought period, magnitude and standardized precipitation index (SPI) values are presented to depict the relationships between drought duration and magnitude.  相似文献   

17.
ABSTRACT

The applicability of multivariate interpolation and information entropy to optimize the raingauge network in the Mekong River Basin (MRB) is investigated. Three different spatial interpolation methods are tested: inverse distance squared (IDS), ordinary kriging (OK) and gradient plus inverse distance squared (GIDS). The validated results confirm that the GIDS method outperformed IDS and OK. The application of information entropy together with GIDS on a network of 57 gauges provided the same information content (7.34 nat) as could be obtained using all 6788 gauges in the MRB. Combining this result with meteorological and hydrological indicators revealed that the number of gauges for the optimum raingauge network could be reduced to 40. The results imply good applicability of the proposed method, which may be used to help prioritize efforts and funds to maintain the raingauge network in a given river basin.  相似文献   

18.
ABSTRACT

The current state of kriging in subsurface hydrology is critically reviewed. In an application to a region where boreholes already exist, methods of optimal location of additional observation wells for geophysical parameter investigation and optimal interpolation for the purpose of solving the inverse problem are investigated. The particular case of the location of wells for the measurements of transmissivity and hydraulic head in the Kennet Valley Chalk aquifer, UK, is examined. Results of interpolation of measured hydraulic conductivity values by kriging are compared with results from a standard graphical package for interpolation. Reference is also made to the distribution obtained by the inverse method (in which the conductivity distribution is obtained from the head distribution). On the basis of the application, the conditional simulation (in which the generated data are both consistent with field values and the field statistical structure) is deemed to be the best. It is also found that different methods of interpolation give widely different distributions in the case of hydraulic conductivity. It is suggested that the kriged map or conditional map of the transmissivity should serve as the basis for regional discretization to which corrections via the inverse model may be made.  相似文献   

19.
Interpolations of groundwater table elevation in dissected uplands   总被引:3,自引:0,他引:3  
Chung JW  Rogers JD 《Ground water》2012,50(4):598-607
The variable elevation of the groundwater table in the St. Louis area was estimated using multiple linear regression (MLR), ordinary kriging, and cokriging as part of a regional program seeking to assess liquefaction potential. Surface water features were used to determine the minimum water table for MLR and supplement the principal variables for ordinary kriging and cokriging. By evaluating the known depth to the water and the minimum water table elevation, the MLR analysis approximates the groundwater elevation for a contiguous hydrologic system. Ordinary kriging and cokriging estimate values in unsampled areas by calculating the spatial relationships between the unsampled and sampled locations. In this study, ordinary kriging did not incorporate topographic variations as an independent variable, while cokriging included topography as a supporting covariable. Cross validation suggests that cokriging provides a more reliable estimate at known data points with less uncertainty than the other methods. Profiles extending through the dissected uplands terrain suggest that: (1) the groundwater table generated by MLR mimics the ground surface and elicits a exaggerated interpolation of groundwater elevation; (2) the groundwater table estimated by ordinary kriging tends to ignore local topography and exhibits oversmoothing of the actual undulations in the water table; and (3) cokriging appears to give the realistic water surface, which rises and falls in proportion to the overlying topography. The authors concluded that cokriging provided the most realistic estimate of the groundwater surface, which is the key variable in assessing soil liquefaction potential in unconsolidated sediments.  相似文献   

20.
Kriging with external drift for functional data for air quality monitoring   总被引:3,自引:2,他引:1  
Functional data featured by a spatial dependence structure occur in many environmental sciences when curves are observed, for example, along time or along depth. Recently, some methods allowing for the prediction of a curve at an unmonitored site have been developed. However, the existing methods do not allow to include in a model exogenous variables that, for example, bring meteorology information in modeling air pollutant concentrations. In order to introduce exogenous variables, potentially observed as curves as well, we propose to extend the so-called kriging with external drift—or regression kriging—to the case of functional data by means of a three-step procedure involving functional modeling for the trend and spatial interpolation of functional residuals. A cross-validation analysis allows to choose smoothing parameters and a preferable kriging predictor for the functional residuals. Our case study considers daily PM10 concentrations measured from October 2005 to March 2006 by the monitoring network of Piemonte region (Italy), with the trend defined by meteorological time-varying covariates and orographical constant-in-time variables. The performance of the proposed methodology is evaluated by predicting PM10 concentration curves on 10 validation sites, even with simulated realistic datasets on a larger number of spatial sites. In this application the proposed methodology represents an alternative to spatio-temporal modeling but it can be applied more generally to spatially dependent functional data whose domain is not a time interval.  相似文献   

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

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