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.     

Driving forces of tropical deforestation: The role of remote sensing and spatial models

Remote sensing technologies are increasingly used to monitor landscape change in many parts of the world. While the availability of extensive and timely imagery from various satellite sensors can aid in identifying the rates and patterns of deforestation, modelling techniques can evaluate the socioeconomic and biophysical forces driving deforestation processes. This paper briefly reviews some emerging spatial methodologies aimed at identifying driving forces of land use change and applies one such methodology to understand deforestation in Mexico. Satellite image classification, change analysis and econometric modelling are used to identify the rates, hotspots and drivers of deforestation in a case study of the southern Yucatán peninsular region, an enumerated global hotspot of biodiversity and tropical deforestation. In particular, the relative roles of biophysical and socioeconomic factors in driving regional deforestation rates are evaluated. Such methodological approaches can be applied to other regions of the forested tropics and contribute insights to conservation planning and policy.     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

The Physico-Chemical Properties of a Subducted Slab from Garnet Zonation Patterns (Sesia Zone, Western Alps)

Garnets in continentally derived high-pressure (HP) rocks ofthe Sesia Zone (Western Alps) exhibit three different chemicalzonation patterns, depending on sample locality. Comparisonof observed garnet zonation patterns with thermodynamicallymodelled patterns shows that the different patterns are causedby differences in the water content of the subducted protolithsduring prograde metamorphism. Zonation patterns of garnets inwater-saturated host rocks show typical prograde chemical zonationswith steadily increasing pyrope content and increasing XMg,together with bell-shaped spessartine patterns. In contrast,garnets in water-undersaturated rocks have more complex zonationpatterns with a characteristic decrease in pyrope and XMg betweencore and inner rim. In some cases, garnets show an abrupt compositionalchange in core-to-rim profiles, possibly due to water-undersaturationprior to HP metamorphism. Garnets from both water-saturatedand water-undersaturated rocks show signs of intervening growthinterruptions and core resorption. This growth interruptionresults from bulk-rock depletion caused by fractional garnetcrystallization. The water content during burial influences significantly thephysical properties of the subducted rocks. Due to enhancedgarnet crystallization, water-undersaturated rocks, i.e. thoselacking a free fluid phase, become denser than their water-saturatedequivalents, facilitating the subduction of continental material.Although water-bearing phases such as phengite and epidote arestable up to eclogite-facies conditions in these rocks, dehydrationreactions during subduction are lacking in water-undersaturatedrocks up to the transition to the eclogite facies, due to thethermodynamic stability of such hydrous phases at high P–Tconditions. Our calculations show that garnet zonation patternsstrongly depend on the mineral parageneses stable during garnetgrowth and that certain co-genetic mineral assemblages causedistinct garnet zonation patterns. This observation enablesinterpretation of complex garnet growth zonation patterns interms of garnet-forming reactions and water content during HPmetamorphism, as well determination of detailed P–T paths. KEY WORDS: dehydration; high-pressure metamorphism; Sesia Zone; subduction; thermodynamic modelling     

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

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