Modified Markov Chain Monte Carlo Method for Dynamic Data Integration Using Streamline Approach

In this paper, the Markov Chain Monte Carlo (MCMC) approach is used for sampling of the permeability field conditioned on the dynamic data. The novelty of the approach consists of using an approximation of the dynamic data based on streamline computations. The simulations using the streamline approach allows us to obtain analytical approximations in the small neighborhood of the previously computed dynamic data. Using this approximation, we employ a two-stage MCMC approach. In the first stage, the approximation of the dynamic data is used to modify the instrumental proposal distribution. The obtained chain correctly samples from the posterior distribution; the modified Markov chain converges to a steady state corresponding to the posterior distribution. Moreover, this approximation increases the acceptance rate, and reduces the computational time required for MCMC sampling. Numerical results are presented.     

The circumstellar shells of carbon miras

The relation between mass loss rate and pulsation period in carbon Miras is discussed. The dust mass loss rate is very low (about 2 × 10^–10 M/yr) up to aboutP = 380 days, where there is a sudden increase. ForP > 400 days there is a linear relation between logM andP. The change in the mass loss rate near 380 days may be related to radiation pressure on dust becoming effective in driving the outflow.     

Retailing local food in the Scottish-English borders: A supply chain perspective

Local food is championed as one alternative response to industrial systems of food production and supply. While advocacy for local food is high, there is a lack of empirical evidence about the actual shape and scale of such food supply chains, especially from a retail perspective. Using supply chain diagrams, this paper presents a summary of ‘new’ agro-food geographies for five different retail types—farm shops, butchers, caterers, specialist shops, supermarkets/department stores—that all source local food from suppliers in the Scottish-English borders. Presented as five separate ‘shopping trips’, the paper examines where, how and why retailers source local food. Results reveal the complex nature of local food systems, especially in terms of intra-sector competitive dynamics (with a notable tension between direct forms of retail and established (independent) retailers), links and overlaps with ‘normal’ food retail systems and elastic notions of the ‘local’. The paper also draws a key distinction between locally produced and locally supplied food products.     

The Tornquist Zone, a north east inclining lithospheric transition at the south western margin of the Baltic Shield: Revealed through a nonlinear teleseismic tomographic inversion

From July 1996 to August 1997 the TOR project operated 130 seismographs in North Germany, Denmark and South Sweden, with the aim of collecting signals from local, regional and teleseismic earthquakes. This data set is particularly interesting since the seismic antenna crosses the most significant geological boundary in Europe, the Tornquist Zone, which in the northern part is the border between the Baltic Shield and the younger European lithosphere. Previous studies have shown significant physical changes in the crust and upper mantle across this transition zone, including two independent teleseismic tomographic studies of the TOR data set. But these two studies disagree on the orientation of the slope of the transition. Both studies used an iterative linearized inversion method. We will in this work Preprint submitted to Elsevier Science 27 July 2005 present an inversion based on Bayesian statistics, where the solution space is examined in order to study a very large number of tomographic solutions and to examine the solution uniqueness and uncertainty. The method is applied to measurements of 3345 relative teleseismic P-phase travel times from 48 teleseismic earthquakes with good azimuthal coverage with respect to the great circle arc of the TOR array. We find the lithospheric transition to be a north east inclination of around 30° to 45° off vertical.     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.     

GIS工作流系统的长事务处理研究

分析了Sagas事务模型及空间数据版本树技术,设计了一种灵活的长事务补偿控制方案,为集成GIS,MIS及OA的工作流系统提供了有力的事务处理技术支持。     

迭代折射静校正技术在鄂北地震资料处理中的应用

在低速层速度(地表直达波)已知的条件下,利用生产炮初至信息,用加权最小平方法和共轭梯度法建立的折射静校正数学模型,迭代求取长波长、短波长静校正量,在解决鄂尔多斯盆地北部的虚假小幅度构造问题及改善剖面质量方面取得了明显的效果。     

基于工作流技术的动态GIS服务链研究

提出了一种基于工作流技术的动态GIS服务链的实现方法,该方法使用抽象GIS服务组合服务链,由GIS服务动态选择机制将抽象服务映射到具体服务,并结合工作流技术,将工作流描述语言描述的GIS服务链流程送入工作流引擎,由工作流引擎管理和监控服务链的运行。     

GPS进行长边控制测量分析

采用基准和基准转换模型的不同会导致GPS观测边长与坐标计算边长不一致。因此由GPS测量的WGS－84坐标系向我国的北京54坐标系转换计算时，应尽可能选择参考点作对比分析。若满足精度要求，可采用全国精密参数作转换计算，若不能满足精度要求或达不到实测点的定位精度，可先将边长作高程异常改正，再对起算点作方位改正。