Application of Regional Phase Amplitude Tomography to Seismic Verification

v--vWe have applied tomographic techniques to amplitude data to quantify regional phase path effects for use in source discrimination studies. Tomography complements interpolation (kriging) methods by extending our ability to apply path corrections into regions devoid of calibration events, as well as raising levels of confidence in the corrections because of their more physical basis. Our tomography technique solves for resolvable combinations of attenuation, source-generation, site and spreading terms. First difference regularization is used to remove singularities and reduce noise effects.¶In initial tests the technique was applied to a data set of 1488, 1.0 Hz, P_g/L_g amplitude ratios from 13 stations for paths inside a 30° by 40° box covering western China and surrounding regions. Tomography reduced variance 60%, relative to the power-law distance correction traditionally applied to amplitude ratios. Relative P_g/L_g attenuation varied with geologic region, with low values in Tibet, intermediate values in basins and high values for platforms and older crust. Spatial patterns were consistent with previous path effect studies in Asia, especially local earthquake coda-Q. Relative spreading was consistent with expected values for P_g and L_g. Relative site terms were similar to one another, yet some tradeoff with attenuation was evident.¶Tomography residuals followed systematic trends with distance, which may result from the evolution from direct to coda phases, focusing, model tradeoff or data windowing effects. Examination of the residuals using a kriging interpolator showed coherent geographical variations, indicating unmodeled path effects. The residual patterns often follow geological boundaries, which could result from attenuating zones or minor blockages that are too thin to be resolved, or that have anisotropic effect on regional phases. These results will guide efforts to reparameterize tomography models to more effectively represent regional wave attenuation and blockage. The interpolated residuals also can be combined with predictions of the tomographic model to account for path effects in discrimination studies on a station by station basis.     

Declustering experimental variograms by global estimation with fourth order moments

The variogram is a key parameter for geostatistical estimation and simulation. Preferential sampling may bias the spatial structure and often leads to noisy and unreliable variograms. A novel technique is proposed to weight variogram pairs in order to compensate for preferential or clustered sampling . Weighting the variogram pairs by global kriging of the quadratic differences between the tail and head values gives each pair the appropriate weight, removes noise and minimizes artifacts in the experimental variogram. Moreover, variogram uncertainty could be computed by this technique. The required covariance between the pairs going into variogram calculation, is a fourth order covariance that must be calculated by second order moments. This introduces some circularity in the calculation whereby an initial variogram must be assumed before calculating how the pairs should be weighted for the experimental variogram. The methodology is assessed by synthetic and realistic examples. For synthetic example, a comparison between the traditional and declustered variograms shows that the declustered variograms are better estimates of the true underlying variograms. The realistic example also shows that the declustered sample variogram is closer to the true variogram.     

A brief review of kriging and its application to optimal interpolation and observation well selection

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.     

Studying the effect of different approaches to estimating model parameters on the accuracy of river runoff simulation by SWAP model

Different approaches to estimating the parameters of SWAP physically based model, which describes heat and water transfer processes in the soil-vegetaion (snow) cover-atmosphere system are examined. In particular, two methods of a priori estimation of parameter values and two variants of their calibration are discussed. The parameter sets obtained by different methods were used to simulate the runoff from 12 experimental catchments in the eastern USA. The calculations were conducted for a 39-year period (1960–1998) with a 3-hour step. The results of calculations were compared with each other and with measured river runoff values in order to identify the parameter set that is optimal for runoff evaluation. A strategy is proposed for a priori parameter estimation in the case of basins where observational data are too poor to enable parameter calibration.     

Data time-step dependency of conceptual rainfall—streamflow model parameters: an empirical study with implications for regionalisation

Abstract

Employing 1-, 2-, 4-, 6-, 12- and 24-hourly data sets for two catchments (10.6 and 298 km²) in Wales, the calibrated parameters of a unit hydrograph-based model are shown to change substantially over that range of data time steps. For the smaller basin, each model parameter reaches, or approaches, a stable value as the data time step decreases, providing a straightforward method of estimating time-step independent model parameter values. For the larger basin, the model parameters also reach, or approach, stable values using hourly data, but, for reasons given in the paper, interpretation of the results is more difficult. Model parameter sensitivity analyses are presented that give insights into the relative precision on the parameters for both catchments. The paper discusses the importance of accounting for model parameter data time-step dependency in pursuit of a reduction in the uncertainty associated with estimates of flow in ungauged basins, and suggests that further work along these lines be undertaken using different catchments and models.     

Generation of runoff series at ungauged locations using empirical orthogonal functions in combination with kriging

The method of Empirical Orthogonal Functions (EOF method) is combined with an objective interpolation technique, kriging, to generate runoff series at ungauged locations. In a case study the results are compared to series interpolated by a combination of EOF analysis and regression using catchment characteristics as independent variables. The results are also compared to linear weighting of an existing runoff series, a commonly used method for spatial interpolation. The influence of altitude on the runoff is studied comparing kriging based on 2 and 3 coordinates. The study showed that the capacity of EOF analysis combined with kriging is as good as the traditionally used linear weighting. The results, when altitude is included in the kriging, are improved.     

Application of distributed hydrological models for predictions in ungauged basins: a method to quantify predictive uncertainty

Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparison of approaches for predicting solute transport: sandbox experiments

The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.     

Comparison of ordinary kriging and artificial neural network for spatial mapping of arsenic contamination of groundwater

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.     

Minimum/maximum autocorrelation factors applied to grade estimation

There are many situations in the mining industry where grade estimation of multiple correlated variables is required. The resulting model is expected to reproduce the data correlation, but there is no guarantee that the correlation observed among data will be reproduced by the model if the variables are independently estimated by kriging, and the correlation is not explicitly taken into account. The best geostatistical approach to address this estimation problem is to use co-kriging, which requires both cross and direct covariance modeling of all variables. However, the co-kriging method is labor-intensive when the problem involves more than three attributes. An alternative is to decorrelate the variables and estimate each one independently, using, for instance, the minimum/maximum autocorrelation factors (MAF) approach. This method involves the application of a linear transformation to the correlated variables, transforming the original data into a space where they are uncorrelated. The resulting transformed data can be individually estimated using kriging, avoiding the use of the linear model of coregionalization. Once the kriging has been performed, the MAF estimates are back-transformed to the original data space, re-establishing their correlation.The methodology is illustrated in a case study where there are two variables with correlation coefficient, ρ = ?0.98. The MAF transformation was applied in combination with ordinary kriging (herein denoted as KMAF). Co-kriging was performed to provide a benchmark for comparing the results obtained through KMAF. The results obtained by co-kriging and KMAF showed less than 1 % average deviation between the two block models.     

Simulated precipitation fields with variance-consistent interpolation

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.     

Spatial prediction of river channel topography by kriging

Topographic information is fundamental to geomorphic inquiry, and spatial prediction of bed elevation from irregular survey data is an important component of many reach‐scale studies. Kriging is a geostatistical technique for obtaining these predictions along with measures of their reliability, and this paper outlines a specialized framework intended for application to river channels. Our modular approach includes an algorithm for transforming the coordinates of data and prediction locations to a channel‐centered coordinate system, several different methods of representing the trend component of topographic variation and search strategies that incorporate geomorphic information to determine which survey data are used to make a prediction at a specific location. For example, a relationship between curvature and the lateral position of maximum depth can be used to include cross‐sectional asymmetry in a two‐dimensional trend surface model, and topographic breaklines can be used to restrict which data are retained in a local neighborhood around each prediction location. Using survey data from a restored gravel‐bed river, we demonstrate how transformation to the channel‐centered coordinate system facilitates interpretation of the variogram, a statistical model of reach‐scale spatial structure used in kriging, and how the choice of a trend model affects the variogram of the residuals from that trend. Similarly, we show how decomposing kriging predictions into their trend and residual components can yield useful information on channel morphology. Cross‐validation analyses involving different data configurations and kriging variants indicate that kriging is quite robust and that survey density is the primary control on the accuracy of bed elevation predictions. The root mean‐square error of these predictions is directly proportional to the spacing between surveyed cross‐sections, even in a reconfigured channel with a relatively simple morphology; sophisticated methods of spatial prediction are no substitute for field data. Copyright © 2007 John Wiley & Sons, Ltd.     

The role of debris supply conditions in predicting debris flow activity

Debris flow frequency and magnitude were determined for 33 basins in southwest British Columbia. Basins were first classified as either weathering-limited or transport-limited using a discriminant function based on debris-contributing area, an area-weighted terrain stability number, and drainage density. Multiple regression was used to predict magnitude, peak discharge, frequency and activity (frequency times magnitude) within each group of basins. Model performance was improved by stratifying the total sample of debris flow basins into weathering-and transport-limited groups. Explained variance increased by an average of 15 per cent in the transport-limited sample, indicating that sediment supply conditions in the more active basins are fundamental in predicting debris flow activity. An independent test of the regression models with 11 basins yielded generally good results for debris flow magnitude and peak discharge. Prediction of debris flow frequency proved problematical in weathering-limited basins. The methods developed here provide estimates of debris flow attributes in basins for which few data on past events are available. Copyright © 1999 John Wiley & Sons, Ltd.     

Multigaussian kriging for point-support estimation: incorporating constraints on the sum of the kriging weights

In the geostatistical analysis of regionalized data, the practitioner may not be interested in mapping the unsampled values of the variable that has been monitored, but in assessing the risk that these values exceed or fall short of a regulatory threshold. This kind of concern is part of the more general problem of estimating a transfer function of the variable under study. In this paper, we focus on the multigaussian model, for which the regionalized variable can be represented (up to a nonlinear transformation) by a Gaussian random field. Two cases are analyzed, depending on whether the mean of this Gaussian field is considered known or not, which lead to the simple and ordinary multigaussian kriging estimators respectively. Although both of these estimators are theoretically unbiased, the latter may be preferred to the former for practical applications since it is robust to a misspecification of the mean value over the domain of interest and also to local fluctuations around this mean value. An advantage of multigaussian kriging over other nonlinear geostatistical methods such as indicator and disjunctive kriging is that it makes use of the multivariate distribution of the available data and does not produce order relation violations. The use of expansions into Hermite polynomials provides three additional results: first, an expression of the multigaussian kriging estimators in terms of series that can be calculated without numerical integration; second, an expression of the associated estimation variances; third, the derivation of a disjunctive-type estimator that minimizes the variance of the error when the mean is unknown.     

南海北部大陆边缘盆地地热特征与油气富集

本文报道了莺歌海盆地、北部湾盆地共148个新测热导率数据,根据收集的钻井温度数据新增计算65个大地热流数据;结合前人研究成果绘制了南海北部大陆边缘沉积盆地的地温梯度图、大地热流分布图;系统归纳了南海北部大陆边缘油气勘探成果.结果表明,南海北部大陆边缘珠江口盆地、琼东南盆地、北部湾盆地、莺歌海盆地的平均热流值分别为68.7±11 mW/m²、71.1±13 mW/m²、65.7±8.9 mW/m²、74.7±10 mW/m²,属于典型的“热盆”.热流区域分布特征总体上受大地构造背景控制,随地壳厚度从北向南,由陆架到陆坡区逐渐减薄而增高,水热活动与岩浆活动等是引起局部高热流异常的原因.盆地地温场的差异控制和约束了油气分布富集规律,从研究区油气勘探成果中可以发现,该区域的气田多发育于高热流盆地（凹陷）,而中-低热流盆地（凹陷）则多孕育油田,油气田具有“北油南气”的分布特征.     

Transferability of hydrological model parameters between basins in data-sparse areas, subarctic Canada

Hydrological models used for the simulation of runoff are often calibrated only on the basis of data obtained at the catchment outlet but the parameters thus derived are then applied to the simulations for the subbasins. Such a practice is common for the data-sparse areas such as the subarctic. However, it may yield erroneous results when the calibrated model parameters are applied to basins of various sizes, or with divergent physical characteristics. This study assesses the feasibility of transferring parameter estimates derived for one basin of a particular size to other basins of different dimensions, using the SLURP model for simulation and the Liard and two of its subbasins as an example. Results indicate that other than the snowmelt factor, the parameter values obtained from the subbasins are similar, but values of several parameters (e.g. maximum capacity of the soil water and groundwater storage, and snowmelt factor) are different from those derived for the large basin. Compared with applying the Liard basin parameters, the subbasins parameter sets generate higher evapotranspiration, earlier termination of the snowmelt period, more soil water storage, a shorter period with significant soil water storage and a better overall agreement between the observed and simulated runoff. It is recommended that adequate attention be given to the transferability of the parameter values to improve the simulation of subbasins hydrology.     

Observed hydrographs: on their ability to infer a time‐invariant hydrological transfer function for flow prediction in ungauged basins

Streamflow measurements provide information about the flow generation characteristics of land surfaces as well as the flow transferring nature of the channel network. In this study, such flow transferring properties of the channel network that were obtained from downstream flow observations were used for predicting flow in ungauged basins. A temporally averaged transfer function (ATF) of the channel segments of Kentucky River Basin (KRB) in Kentucky, USA, was extracted from observed hydrographs in a time‐invariant system as a function of drainage area. The ATF was regionalized through multiple regression analysis for 194 combinations of drainage areas that differ in topography, terrain, and geology. The application of ATF for flow prediction in ungauged basins was performed for Goose Creek, a subbasin of KRB by integrating ATF with the TOPMODEL. In addition, the ATF was shown to be capable of providing calibration and validation data for ungauged basins in a backward technique from a measured stream gauge downstream, with minimal data requirement of drainage area. The applicability of ATF was illustrated across a range of streamflow conditions from watersheds that varied greatly in their terrain and geology. Nash–Sutcliffe efficiency of the proposed method, as a function of drainage areas of the corresponding basins, to predict daily streamflow from ungauged basins ranged from 0.83 to 0.92. The results of the study concluded that the ATF obtained from measured streamflow thus proved to be a quick and simple tool for assessment of streamflow in both operational and modeling hydrology. Copyright © 2012 John Wiley & Sons, Ltd.     

Predictions in ungauged basins: an approach for regionalization of hydrological models considering the probability distribution of model parameters

Regionalization of model parameters by developing appropriate functional relationship between the parameters and basin characteristics is one of the potential approaches to employ hydrological models in ungauged basins. While this is a widely accepted procedure, the uniqueness of the watersheds and the equifinality of parameters bring lot of uncertainty in the simulations in ungauged basins. This study proposes a method of regionalization based on the probability distribution function of model parameters, which accounts the variability in the catchment characteristics. It is envisaged that the probability distribution function represents the characteristics of the model parameter, and when regionalized the earlier concerns can be addressed appropriately. The method employs probability distribution of parameters, derived from gauged basins, to regionalize by regressing them against the catchment attributes. These regional functions are used to develop the parameter characteristics in ungauged basins based on the catchment attributes. The proposed method is illustrated using soil water assessment tool model for an ungauged basin prediction. For this numerical exercise, eight different watersheds spanning across different climatic settings in the USA are considered. While all the basins considered in this study were gauged, one of them was assumed to be ungauged (pseudo-ungauged) in order to evaluate the effectiveness of the proposed methodology in ungauged basin simulation. The process was repeated by considering representative basins from different climatic and landuse scenarios as pseudo-ungauged. The results of the study indicated that the ensemble simulations in the ungauged basins were closely matching with the observed streamflow. The simulation efficiency varied between 57 and 61 % in ungauged basins. The regional function was able to generate the parameter characteristics that were closely matching with the original probability distribution derived from observed streamflow data.     

Prediction of streamflow regimes over large geographical areas: interpolated flow–duration curves for the Danube region

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.     

Modeling locally varying anisotropy of CO2 emissions in the United States

In a spatial property modeling context, the variables of interest to be modeled often display complex nonlinear features. Techniques to incorporate these nonlinear features, such as multiple point statistics or cummulants, are often complex with input parameters that are difficult to infer. The methodology proposed in this paper uses a classical vector-based definition of locally varying anisotropy to characterize nonlinear features and incorporate locally varying anisotropy into numerical property models. The required input is an exhaustive field of anisotropy orientation and magnitude. The methodology consists of (1) using the shortest path distance between locations to define the covariance between points in space (2) multidimensional scaling of the domain to ensure positive definite kriging equations and (3) estimation or simulation with kriging or sequential Gaussian simulation. The only additional parameter required when kriging or simulating with locally varying anisotropy is the number of dimensions to retain in multidimensional scaling. The methodology is demonstrated on a CO₂ emissions data set for the United States in 2002 and shows an improvement in cross validation results as well as a visual reproduction of nonlinear features.