Some Results on the Spatial Breakdown Point of Robust Point Estimates of the Variogram

The effect of outliers on estimates of the variogram depends on how they are distributed in space. The ‘spatial breakdown point’ is the largest proportion of observations which can be drawn from some arbitrary contaminating process without destroying a robust variogram estimator, when they are arranged in the most damaging spatial pattern. A numerical method is presented to find the spatial breakdown point for any sample array in two dimensions or more. It is shown by means of some examples that such a numerical approach is needed to determine the spatial breakdown point for two or more dimensions, even on a regular square sample grid, since previous conjectures about the spatial breakdown point in two dimensions do not hold. The ‘average spatial breakdown point’ has been used as a basis for practical guidelines on the intensity of contaminating processes that can be tolerated by robust variogram estimators. It is the largest proportion of contaminating observations in a data set such that the breakdown point of the variance estimator used to obtain point estimates of the variogram is not exceeded by the expected proportion of contaminated pairs of observations over any lag. In this paper the behaviour of the average spatial breakdown point is investigated for cases where the contaminating process is spatially dependent. It is shown that in two dimensions the average spatial breakdown point is 0.25. Finally, the ‘empirical spatial breakdown point’, a tool for the exploratory analysis of spatial data thought to contain outliers, is introduced and demonstrated using data on metal content in the soils of Sheffield, England. The empirical spatial breakdown point of a particular data set can be used to indicate whether the distribution of possible contaminants is likely to undermine a robust variogram estimator.     

Influence of initial plastic anisotropy of overconsolidated clays on ground behaviour during tunneling

In order to evaluate the influence of the initial plastic anisotropy in the excavation of a tunnel, a “bubble” bounding surface model for structured soils, formulated by Kavvadas and Belokas (Proc. 10th IACMAG Conf., 2001), was implemented in the explicit finite difference code FLAC. Two different initial stress (K ₀) conditions were considered. The size and shape of the initial bounding surfaces were specified to be consistent with the initial stress field. The distorted and rotated shape of the bounding surface, supported by experimental results, defines the anisotropy of shear strength, which is shown to have a significant influence on the displacements. There is also considerable sensitivity of the soil model to the initial stress field.     

Upper Cretaceous carbonate hosted zinc–lead–barite deposits in Northern Thrust Zone, northern Iraq: petrography and geochemistry

Zinc–lead–barite deposits located in Lefan and Lower Banik localities of about 25 km northeast of Zakho City, Northern Iraq consist of a group of strata-bound sulfides hosted in Upper Cretaceous (Upper Campanian–Maastrichtian) dolomitic limestone. Carbonate-hosted ores contain 3.77% Zn, 2% Pb, and 5% Fe, while in lower Banik, they contain 1.5% Zn, 0.37% Pb, and 1.4% Fe. Diagenetic processes, such as dolomitization and recrystalization in addition to the type of microfacies, provided appropriate physical and chemical conditions that permitted the passage of ore-bearing fluids and participated in precipitation and ore localization. These deposits are precipitated in a platform and developed within the Foreland Thrust Belt. Ore precipitated as infill of intergranular dolomite porosity with replaced dolomite and rudist shells forming disseminated crystals that occupy intergranular pore spaces around dolomite and calcite and as infill of dissolution spaces and fractures.     

Strain-softening behaviour of sand in strain path testing under plane-strain conditions

Experimental data are presented in this paper to study the strain-softening behaviour of sand under plane-strain conditions. K ₀ consolidated strain path tests were conducted using a new plane-strain apparatus. The stress–strain behaviour of medium dense sand under plane-strain conditions was characterized. The test results show that the occurrence of pre-failure strain softening under plane-strain conditions is affected by the void ratio, the strain increment ratio and the initial effective confining stress. This is consistent with previous findings established under axisymmetric conditions. However, a pre-failure strain-softening behaviour in plane-strain tests conducted under high-confining stresses may consist of three stages, namely, material softening, banding softening, and ultimate state. This observation is different from that in triaxial tests where banding softening does not normally occur.     

Attenuation relations for strong seismic ground motion in the Himalayan region

Strong motion data from various regions of India have been used to study attenuation characteristics of horizontal peak acceleration and velocity. The strong ground motion data base considered in the present work consists of various earthquakes recorded in the northern part of India since 1986 with magnitudes 5.7 to 7.2. Using these data, relations for horizontal peak acceleration and velocity, which are $$\begin{gathered} log_{10} a = 1.14 + 0.31M + 0.65log_{10} R \hfill \\ log_{10} v = 0.571 + 0.41M + 0.768log_{10} R \hfill \\ \end{gathered} $$ have been proposed wherea is the peak horizontal acceleration in cm/sec²,v is the peak horizontal velocity in mm/sec,M is body wave magnitude, andR is the hypocentral distance in km. The proposed relations are in reasonable agreement with the small amount of strong ground motion data available for the northern part of India. The present results will be useful in estimating strong ground motion parameters and in the earthquake resistant design in the Himalayan region.     

The qualitative zoning record of minerals. A method for determining the duration of metamorphic events?

Summary The efficacy of ion exchange in petrological systems, for example the Fe-Mg exchange between garnet and biotite, is a function of grainsize, temperature, rates of temperature change and diffusion parameters. The combination of these variables determines the final zoning profile of minerals. Therefore, zoning profiles may be used to derive one of these variables if the others are known. For example, a mineral grain that experienced a short thermal event may still preserve a zoning profile characteristic of the heating path. In contrast, if that grain experienced a long thermal event it may develop a zoning profile typical of the cooling path. Conversely, for the same temperature-time cycle, large grains may not completely equilibrate at the metamorphic peak, and in smaller grains any previous record may have been erased. This is commonly observed in natural rocks where different grains within one thin section preserve often qualitatively different zoning profiles. Thus, a critical grainsize,l _crit, may exist that separates grains with qualitatively different zoning profiles so that grains of sizel > l _crit still retain zoning information about the heating path and grains of sizel <l _crit contain only information about the cooling path. If the critical grainsize can be measured and an independent estimate for the peak metamorphic temperature exists, the duration of the thermal event may, in principle, be estimated.The applicability of this method to natural garnets is hampered by a range of uncertainties. However, a parameterisation of the critical grain size may be used to illustrate the relative importance of grain size, temperature and event duration to the equilibration of minerals. In this paper, the critical grainsize is parameterized for various temperature-time cycles and its dependence on diffusion parameters is discussed. It is shown that, for Barrovian conditions,l _crit between 0.1 and 1 mm separates garnets retaining prograde information from garnets retaining cooling path information. For illustration, we compare the results with the critical grainsize of garnets from two metamorphic terranes, the Prydz Bay region (Antarctica) and the Koralm complex (Eastern Alps). Despite the large range of uncertainties attached to the method, it is shown that the critical grainsize of garnets in both terranes is consistent with a very short duration of the last thermal event that affected the two regions.

---

Die qualitative Form von Zonierungsprofilen in Mineralen. Information über die Dauer metamorpher Ereignisse?

Zusammenfassung Zonierungsprofile von Paragenesen im Ionenaustausch, zum Beispiel das Fe-Mg Austauschsgleichgewicht zwischen Granat und Biotit, sind eine Funktion von (i) Korngröße, (ii) Temperatur, (iii) Kühl- und Heizgeschwindigkeit sowie (iv) den Diffusionskonstanen. Es sollte daher möglich sein, Zonierungsprofile dazu zu verwenden, einen dieser Parameter zu bestimmen, wenn die anderen bekannt sind. Zonierungsprofile von Kristallen die ein kurzes thermisches Ereignis erfahren haben, mögen daher prograde Information beinhalten, wogegen derselbe Kristall ein retrogrades Profil aufweisen mag, wenn er einem langlebigen thermischen Ereignis unterlag (prograd und retrograd ist hier als: bei ansteigender Temperatur und bei abfallender Temperatur definiert). In Kristallen die dem gleichen Temperatur-Zeitpfad unterworfen waren, könnten kleine Körner nur die Kühlgeschichte dokumentieren, wogegen große noch Zonierungsinformation vom Heizpfad aufweisen. In natürlichen Gesteinen wird das oft dadurch beobachtet, daß verschiedene Körner auseinem Dünnschliff qualitativ verschiedene Zonierungsprofile aufweisen. Es ist daher möglich eine kritische Korngöße,l _crit, zu definieren, die Korngrößen mit qualitativ verschiedenen Zonierungsprofilen voneinander trennt. Körner mit einem Durchmesserl > l _crit haben, zumindest teilweise, noch prograde Profile, wogegen Körner mitl <l _crit nur retrograde Information dokumentieren. Wenn man diese kritische Korngröße messen kann, sollte sie dazu benutzt werden können, etwas über die Dauer des thermischen Ereignisses auszusagen.Die Anwendbarkeit dieser Methode ist durch eine Reihe von Fehlern limitiert. Nichtsdestotrotz, ist eine Parameterisierung der kritischen Korngröße nützlich um die relative Wichtigkeit von Korngröße, Temperatur und Zeit, für die Entwicklung von Zonierungsprofilen, zu illustrieren. Unser Modell kann dazu benutzt werden, um die Größenordnung von Granaten abzuschätzen, die noch prograde Information dokumentieren können. Die Abhängigkeit der kritischen Korngröße von Diffusionsparametern und der Form des Temperatur-Zeit Pfades wird ebenfalls diskutiert. Es wird gezeigt, daß für Metamorphose in der mittleren Grünschiefer und Amphibolit fazies,l _crit zwischen 0.1 and 1 mm prograd zonierte von retrograd zonierten Granaten trennt. Um das Modell zu illustrieren, werden die Ergebnisse auf die Prydz Bay Region (Antarktis) und den Koralm Komplex (Ostalpen) angewendet. Trotz der großen Fehlergrenzen der Methode glauben wir zeigen zu können, daß die kritische Korngröße beider Terrains andeutet, daß das letzte thermische Ereignis in beiden Gebieten kurzlebig war.

---

---

With 5 Figures     

Transition from eclogite to amphibolite-facies metamorphism in the Austroalpine Ulten Zone

Summary The Ulten Zone of the Austroalpine crystalline basement south-west of Meran (Italy) contains metapelitic schists and granoblastic paragneisses, leucocratic orthogneisses, migmatites (in both gneiss-lithologies), metabasites and ultramafic lenses. Metamorphic textures of the metapelitic schists and granoblastic paragneisses indicate two different metamorphic events, characterized by two mineral assemblages, which differ in mineral chemistry: (1) an eclogite facies mineral assemblage (M1) comprising Grt-Ky I-Bt. Ms-Kfs-PI-Qtz-Rt, and (2) an amphibolite facies mineral assemblage (M2) comprising Grt-KyII-Bt-Ms-PI-Qtz-Ilm±St. For the M1 event, pressures of at least 15kbar and temperatures of about 700°±50°C can be estimated. The later amphibolite facies overprint occurred at pressures of 6 to 8kbar and about 600°±50°C. The M1 and M2 assemblages belong to a continuous clockwise metamorphic evolution during the Variscan orogeny. Evidence for Alpine metamorphism can only be detected by sericite rims around kyanite and reset biotite ages. The migmatites, which contribute about 15–30vol.% of all rocks in the investigated area, were formed on the prograde path during the M1 event. Dissolution of H₂O in the melted part of the migmatites resulted in a CO₂dominated fluid, which was trapped in primary kyanite (M1) fluid inclusions. Secondary H₂O-rich fluid inclusions are found in quartz grains and may represent the fluid which enabled a pervasive equilibration during M2.

---

Übergang von eklogit-zu amphibolitfazieller Matamorphose in der austroalpinen Ultenzone

Zusammenfassung Die Ulten Zone, ein Teil des ostalpinen kristallinen basements, südwestlich von Meran, wird aus Metapeliten and granoblastischen Paragneisen, leukokraten Orthogneisen, Migmatiten (in beiden Lithologien), Metabasiten and ultramafischen Linsen aufgebaut. Metamorphe Texturen der Metapelite und granoblastischen Paragneise lassen auf zwei verschiedene metamorphe Ereignisse schließen, die durch unterschiedliche Mineral-chemismen und Paragenesen charakterisiert sind: (1) eine eklogitfazielle Paragenese (M1), bestehend aus Grt-KyI-Bt-Ms-Kfs-P1-Qtz-Rt und (2) eine amphibolitfazielle Paragenese (M2), bestehend aus Grt-KyII-Bt-Ms-P1-Qtz-Ilm±St. Für M1 konnten Minimaldrucke von 15kbar und Temperaturen von 700°±50°C abgeleitet werden. Die spätere amphibolitfazielle Überprägung fand bei 6 bis 8kbar und 600°±50°C statt. M1 und M2 gehören einer kontinuierlichen Metamorphoseentwicklung während der variszischen Orogenese an.Die Migmatite, ungefähr 15–30vol.% der Gesteine im untersuchten Gebiet, wurden am prograden Pfad während des M1 Ereignisses gebildet. Aufgrund der höheren Löslichkeit von H₂0 in der Schmelze, blieb ein CO₂, reiches Fluid zurück, das im primären Kyanit (M1) eingeschlossen wurde. Wässrige Flüssigkeitseinschlüsse können in Quarzkörnern gefunden werden. Dieses Fluid ist wahrscheinlich für die Reequilibrierung zu amphibolitfaziellen Bedingungen verantwortlich.

---

---

With 5 Figures     

Spatial orientation of the spin vectors of the disk galaxies in the Virgo cluster I. A catalogue of the disk galaxies in the Virgo cluster

Based on four UK Schmidt Telescope (UKST) IIIa-j direct plates which had been digitized by the Automated Plate Measuring System (APM), a detailed identification of the disk galaxies in the Virgo cluster has been carried out. The certain and possible member galaxies are mainly selected from the updated Catalogue of 2096 Galaxies in the Virgo Cluster Area (VCC). The area of this identification covers a 85 square degree region centered on this cluster. As a result, we compiled a new catalogue of 340 certain and possible member disk galaxies, including 85 lenticulars (SO), 136 spirals (S) and 119 irregulars (Irr), where the APM-measured values of the position angle (P.A.) and diameters of 299 disk galaxies are provided. For at least 149 galaxies, the values of the P.A. and diameters are published for the first time.Comparison between the APM-measured P.A. and diameters of galaxies and those given by the Uppsala General Catalog of Galaxies (UGC) demonstrates the reliability of our data. The statistical properties and completeness of this database are also discussed. Compared with the UGC-given P.A. and diameters which were determined by eyes, our data has no selection effect owing to the instrument measurements. Within our searching region, this might be the largest and relatively complete photometric database published so far for investigating the spatial orientation of the spin vectors of galaxies in the Virgo cluster.     

Gravitational mass of magnetostatic field

Maxwell's displacement current equation is interpreted in the light of recent work to show that static magnetic field in free space should have a colocated and contemporaneous mass that is neither embodied in, nor can be anticipated from, the mass-energy relation. Thus magnetostatic field in the universe represents an invisible mass. Some consequences are discussed.     

Towards a sustainable economy for the Shetland Islands: Development and management issues in fishing and fish farming

Shetland, as an isolated group of islands on the north-west margin of Europe, has always depended in its economy on marine-based resources, and current trends suggest that this will continue for the future. Up to the oil era Shetland depended mainly on the fisheries around it, although for centuries externally-based fishermen have profited more from these fisheries than have the Shetland people. As well as the development of oil from the 1970s there has been the other major economic innovation of fish farming. The oil industry is on massive scale and has greatly changed the economy of the islands and brought notable social benefits; but the direct employment it provides is considerably exceeded by the fish-related industries, which include employment in fish processing and packing plants as well as in fish catching and fish farming. With the oil industry in Shetland now on a settled downward trajectory, the islands are inevitably now having to look more to fish-related employment for their future. In the internal structure of Shetland, fish farming has a notable effect in distributing employment throughout the islands. Activity in fishing and fish processing is more selectively located: fishing is mainly concentrated in certain islands in the archipelago, particularly Whalsay and Burra; and fish processing is mainly in central locations, especially Lerwick and Scalloway. A prominent problem for fish processors is getting reliable supplies in face of higher fish prices at markets on the Scottish mainland. The processors also now face the challenge of putting more value added on their products in an outlying location. For the future there is a major issue in getting more effective measures for conservation of fish stocks in the Shetland area; there is also the related issue of whether Shetland fishermen should get preference in the allocation of fish quotas in this area. For fish farming there is the basic issue of maintaining a pollution-free environment: this was emphasised by the oil spill from the Braer oil tanker, and there is also concern about potential polluution from the prospective development of oil fields in the more difficult deep-water environment to the west of the islands.