排序方式： 共有24条查询结果，搜索用时 31 毫秒
1.
Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data 
下载免费PDF全文
下载免费PDF全文 Antonio‐Juan Collados‐Lara Eulogio Pardo‐Igúzquiza David Pulido‐Velazquez 《水文研究》2017,31(10):1966-1982
Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources. 相似文献
2.
Compositional Bayesian indicator estimation 总被引：1，自引：1，他引：0
Carolina Guardiola-Albert Eulogio Pardo-Ig��zquiza 《Stochastic Environmental Research and Risk Assessment (SERRA)》2011,25(6):835-849
Indicator kriging is widely used for mapping spatial binary variables and for estimating the global and local spatial distributions
of variables in geosciences. For continuous random variables, indicator kriging gives an estimate of the cumulative distribution
function, for a given threshold, which is then the estimate of a probability. Like any other kriging procedure, indicator
kriging provides an estimation variance that, although not often used in applications, should be taken into account as it
assesses the uncertainty of the estimate. An alternative approach to indicator estimation is proposed in this paper. In this
alternative approach the complete probability density function of the indicator estimate is evaluated. The procedure is described
in a Bayesian framework, using a multivariate Gaussian likelihood and an a priori distribution which are both combined according
to Bayes theorem in order to obtain a posterior distribution for the indicator estimate. From this posterior distribution,
point estimates, interval estimates and uncertainty measures can be obtained. Among the point estimates, the median of the
posterior distribution is the maximum entropy estimate because there is a fifty-fifty chance of the unknown value of the estimate
being larger or smaller than the median; that is, there is maximum uncertainty in the choice between two alternatives. Thus
in some sense, the latter is an indicator estimator, alternative to the kriging estimator, that includes its own uncertainty.
On the other hand, the mode of the posterior distribution estimator, assuming a uniform prior, is coincidental with the simple
kriging estimator. Additionally, because the indicator estimate can be considered as a two-part composition which domain of
definition is the simplex, the method is extended to compositional Bayesian indicator estimation. Bayesian indicator estimation
and compositional Bayesian indicator estimation are illustrated with an environmental case study in which the probability
of the content of a geochemical element in soil being over a particular threshold is of interest. The computer codes and its
user guides are public domain and freely available. 相似文献
3.
Geostatistical estimation (kriging) and geostatistical simulation are routinely used in ground water hydrology for optimal spatial interpolation and Monte Carlo risk assessment, respectively. Both techniques are based on a model of spatial variability (semivariogram or covariance) that generally is not known but must be inferred from the experimental data. Where the number of experimental data is small (say, several tens), as is not unusual in ground water hydrology, the model fitted to the empirical semivariogram entails considerable uncertainty. If all the practical results are based on this unique fitted model, the final results will be biased. We propose that, instead of using a unique semivariogram model, the full range of models that are inside a given confidence region should be used, and the weight that each semivariogram model has on the final result should depend on its plausibility. The first task, then, is to evaluate the uncertainty of the model, which can be efficiently done by using maximum likelihood inference. The second task is to use the range of plausible models in applications and to show the effect observed on the final results. This procedure is put forth here with kriging and simulation applications, where the uncertainty in semivariogram parameters is propagated into the final results (e.g., the prediction of ground water head). A case study using log-transmissivity data from the Vega de Granada aquifer, in southern Spain, is given to illustrate the methodology. 相似文献
4.
Eulogio?Pardo-Igúzquiza Pedro?Martínez-SantosEmail author Miguel?Martín-Loeches 《Stochastic Environmental Research and Risk Assessment (SERRA)》2018,32(8):2433-2444
This paper deals with the design of optimal spatial sampling of water quality variables in remote regions, where logistics are complicated and the optimization of monitoring networks may be critical to maximize the effectiveness of human and material resources. A methodology that combines the probability of exceeding some particular thresholds with a measurement of the information provided by each pair of experimental points has been introduced. This network optimization concept, where the basic unit of information is not a single spatial location but a pair of spatial locations, is used to emphasize the locations with the greatest information, which are those at the border of the phenomenon (for example contamination or a quality variable exceeding a given threshold), that is, where the variable at one of the locations in the pair is above the threshold value and the other is below the threshold. The methodology is illustrated with a case of optimizing the monitoring network by optimal selection of the subset that best describes the information provided by an exhaustive survey done at a given moment in time but which cannot be repeated systematically due to time or economic constrains. 相似文献
5.
Pardo-Igúzquiza Eulogio Dowd Peter A. Rizo-Decelis Luis David 《Mathematical Geosciences》2020,52(5):639-650
Mathematical Geosciences - The universality of fractals implies that very different physical processes can give rise to similar complex spatial patterns. Sinkholes (dolines) and galaxies provide a... 相似文献
6.
Empirical Maximum Likelihood Kriging: The General Case 总被引：4，自引：0，他引：4
Although linear kriging is a distribution-free spatial interpolator, its efficiency is maximal only when the experimental data follow a Gaussian distribution. Transformation of the data to normality has thus always been appealing. The idea is to transform the experimental data to normal scores, krige values in the “Gaussian domain” and then back-transform the estimates and uncertainty measures to the “original domain.” An additional advantage of the Gaussian transform is that spatial variability is easier to model from the normal scores because the transformation reduces effects of extreme values. There are, however, difficulties with this methodology, particularly, choosing the transformation to be used and back-transforming the estimates in such a way as to ensure that the estimation is conditionally unbiased. The problem has been solved for cases in which the experimental data follow some particular type of distribution. In general, however, it is not possible to verify distributional assumptions on the basis of experimental histograms calculated from relatively few data and where the uncertainty is such that several distributional models could fit equally well. For the general case, we propose an empirical maximum likelihood method in which transformation to normality is via the empirical probability distribution function. Although the Gaussian domain simple kriging estimate is identical to the maximum likelihood estimate, we propose use of the latter, in the form of a likelihood profile, to solve the problem of conditional unbiasedness in the back-transformed estimates. Conditional unbiasedness is achieved by adopting a Bayesian procedure in which the likelihood profile is the posterior distribution of the unknown value to be estimated and the mean of the posterior distribution is the conditionally unbiased estimate. The likelihood profile also provides several ways of assessing the uncertainty of the estimation. Point estimates, interval estimates, and uncertainty measures can be calculated from the posterior distribution. 相似文献
7.
Mathematical Geosciences - Fractal analysis has, for some time, been used to evaluate the roughness of both natural and engineered surfaces. Since the introduction of fractal geometry, evidence of... 相似文献
8.
Eulogio Pardo-Igúzquiza 《Mathematical Geology》1998,30(1):95-108
In this paper, the maximum likelihood method for inferring the parameters of spatial covariances is examined. The advantages of the maximum likelihood estimation are discussed and it is shown that this method, derived assuming a multivariate Gaussian distribution for the data, gives a sound criterion of fitting covariance models irrespective of the multivariate distribution of the data. However, this distribution is impossible to verify in practice when only one realization of the random function is available. Then, the maximum entropy method is the only sound criterion of assigning probabilities in absence of information. Because the multivariate Gaussian distribution has the maximum entropy property for a fixed vector of means and covariance matrix, the multinormal distribution is the most logical choice as a default distribution for the experimental data. Nevertheless, it should be clear that the assumption of a multivariate Gaussian distribution is maintained only for the inference of spatial covariance parameters and not necessarily for other operations such as spatial interpolation, simulation or estimation of spatial distributions. Various results from simulations are presented to support the claim that the simultaneous use of maximum likelihood method and the classical nonparametric method of moments can considerably improve results in the estimation of geostatistical parameters. 相似文献
9.
Marta Rodrigo-Gámiz Francisca Martínez-Ruiz Francisco J. Rodríguez-Tovar Francisco J. Jiménez-Espejo Eulogio Pardo-Igúzquiza 《Quaternary Research》2014
Cyclostratigraphic analysis conducted on a continuous high-resolution marine record from the western most Mediterranean reveals well-identified paleoclimate cycles for the last 20,000 yr. The detrital proxies used (Si/Al, Ti/Al, Zr/Al, Mg/Al, K/Al, Rb/Al) are related to different sediment-transport mechanisms, including eolian dust and fluvial runoff, which involve fluctuations in the atmosphere–hydrosphere systems. These fluctuations are accompanied by changes in marine productivity (supported by Ba/Al) and bottom-water redox conditions (Cu/Al, V/Al, Zn/Al, Fe/Al, Mn/Al, U/Th). Spectral analysis conducted using the Lomb–Scargle periodogram and the achieved significance level implemented with the permutation test allowed us to establish major periodicities at 1300, 1515, 2000, and 5000 yr, and secondary peaks at 650, 1087, and 3000 yr. Some of these cycles also agree with those previously described in the North Atlantic Ocean and circum-Mediterranean records. The periodicities obtained at 2000 and 5000 yr support a global connection with records distributed at high, mid, and low latitudes associated with solar activity, monsoonal regime and orbital forcing. The 1300- and 1515-yr cycles appear to be linked with North Atlantic climate variability and the African monsoon system. Thus, the analyzed record provides evidence of climate cycles and plausible forcing mechanisms coupled with ocean–atmosphere fluctuations. 相似文献
10.
Generalized Bootstrap Method for Assessment of?Uncertainty in Semivariogram Inference 总被引：1，自引：0，他引：1
The semivariogram and its related function, the covariance, play a central role in classical geostatistics for modeling the average continuity of spatially correlated attributes. Whereas all methods are formulated in terms of the true semivariogram, in practice what can be used are estimated semivariograms and models based on samples. A generalized form of the bootstrap method to properly model spatially correlated data is used to advance knowledge about the reliability of empirical semivariograms and semivariogram models based on a single sample. Among several methods available to generate spatially correlated resamples, we selected a method based on the LU decomposition and used several examples to illustrate the approach. The first one is a synthetic, isotropic, exhaustive sample following a normal distribution, the second example is also a synthetic but following a non-Gaussian random field, and a third empirical sample consists of actual raingauge measurements. Results show wider confidence intervals than those found previously by others with inadequate application of the bootstrap. Also, even for the Gaussian example, distributions for estimated semivariogram values and model parameters are positively skewed. In this sense, bootstrap percentile confidence intervals, which are not centered around the empirical semivariogram and do not require distributional assumptions for its construction, provide an achieved coverage similar to the nominal coverage. The latter cannot be achieved by symmetrical confidence intervals based on the standard error, regardless if the standard error is estimated from a parametric equation or from bootstrap. 相似文献