Seismic Moment Tensor Resolution on a Local Scale: Simulated Rockburst and Mine-induced Seismic Events in the Kopanang Gold Mine, South Africa

Seismic records contain information about the effect of the source as well as the effect of wave propagation through the rock mass. The effect of wave propagation is usually not well known as only simplified models of geological structures are available. Therefore, the information about the source retrieved by inverting seismograms may include errors due to incomplete knowledge of the rock mass along the propagation path, which in turn cause a distortion in the calculated moment tensor (MT). The distortion of the MT on a local scale was observed by inverting records of a simulated rockburst conducted at the Kopanang gold mine in South Africa. A dominant isotropic component of the explosive characteristics was found from the inversion. The deviatoric components retrieved from the blast are spurious. A test of their stability indicated that they are not significant, assuming an uncertainty above 5% for velocities and 10% for attenuation within the homogeneous model available for the mine. Thus, the retrieval of the MT from records of local networks in mines using a homogeneous model of the rock mass seems to be feasible. However, the homogeneous model of the rock mass can only be applied to close stations, within a few kilometers of the source. The seismic records from distant stations were too complex to be modelled by a homogeneous rock mass. Records of six mine-induced seismic events recorded at the Kopanang gold mine were also inverted. A vertical linear dipole along the pressure (P) axis was found for three of the events, suggesting a pillar burst. The mechanism of two events contains an isotropic implosion together with a nearly vertical dip-slip, and seems to indicate a combination of a cavity collapse with a down dip-slip along a nearly vertical fault. One event corresponds to a dipole along the tensional (T) axis. However, it is vertical, thus its association with tensile faulting of the hangingwall is uncertain.     

Height transformation models from ellipsoidal into the normal orthometric height system for the territory of the City of Zagreb

The paper presents the testing of the possibility of determining the heights of GPS points in the homogeneous field in the new Croatian Height Reference System (HVRS71) by using the method of height transformation. The testing was made in the area of Zagreb. As part of the field works, normal orthometric heights of 27 GPS points were determined according to the new height system, by transferring the benchmark heights using the geometric levelling method, thus obtaining GPS/levelling points of known ellipsoidal and normal orthometric heights. The GPS/levelling points served as the basis for determining the transformation models that enabled the computation of normal orthometric heights from ellipsoidal heights of any GPS point in the observed area. The empirical data used for modelling were reduced undulation dN values of GPS/levelling points. As part of the dN modelling with parametric functions, the approximation surfaces were obtained on the basis of three polynomials: FN310, FN312 and FN318. The transformation models were also tested using non-parametric Watson and Loess algorithms. The FN318 and Loess models yielded the best results.     

A parallel mesh-free contaminant transport model based on the Analytic Element and Streamline Methods

This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources.     

Discontinuous Nature of the Lower Part of the South American Wadati-Benioff Zone in The Arica Elbow Region

The Arica Elbow region represents that part of Andean South America where the azimuth of the strike of the Peru-Chile trench changes from 150° to 190°. The area under study is roughly bounded by latitudes 17 °S and 23 °S, The shape of the Wadati-Benioff zone was studied in terms of the distribution of ISC hypocentres dated between 1964 and 1993. A system of 22 vertical cross-sections, perpendicular to the trench axis, and a map of epicentres was used to derive the detailed shape of the Wadati-Benioff zone of the presently descending slab. The distribution of earthquake foci indicates a fingerlike shape of the lower part of the Wadati-Benioff zone beneath the aseismic gap. The slab length shows small changes around 350 km in the northern sections, pronounced length oscillations between 350 and 750 km in the neighbouring central sections and a constant value of 650 km in the southern sections. The dip and thickness of the Wadati-Benioff zone are practically constant in all sections. Fault plane solutions, separated spatially into three zones, were used to estimate the state of stress in the slab.     

Lysosomal membrane stability and metallothioneins in digestive gland of Mussels (Mytilus galloprovincialis Lam.) as biomarkers in a field study

The lysosomal membrane destabilization and the metallothionein content in the digestive gland cells of mussels (Mytilus galloprovincialis Lam.), collected along the east coast of the North Adriatic (Istrian and Kvarnerine coast, Croatia), were examined over a period of four years (1996–1999). The lysosomal membrane stability, as a biomarker of general stress, showed that the membrane labilization period in mussels from polluted, urban- and industrial-related areas was significantly decreased (p<0.05) when compared to mussels from control, clean sea water sites. In the harbour of Rijeka, the most contaminated site, the lysosomal membrane stability was reduced by more than 70% compared to the control. This method also proved to be a useful biomarker for detection of additional stress caused by short-term hypoxia that occurred once during this study inside the polluted and periodically quite eutrophic Pula Harbour. The concentration of metallothioneins in the mussel digestive gland, as a specific biomarker of exposure to heavy metals, did not reveal significant differences (p<0.05) between sites covered by this study.     

Comparison of surface and borehole locations of induced seismicity

Monitoring of induced microseismic events has become an important tool in hydraulic fracture diagnostics and understanding fractured reservoirs in general. We compare microseismic event and their uncertainties using data sets obtained with surface and downhole arrays of receivers. We first model the uncertainties to understand the effect of different acquisition geometries on location accuracy. For a vertical array of receivers in a single monitoring borehole, we find that the largest part of the final location uncertainty is related to estimation of the backazimuth. This is followed by uncertainty in the vertical position and radial distance from the receivers. For surface monitoring, the largest uncertainty lies in the vertical position due to the use of only a single phase (usually P‐wave) in the estimation of the event location. In surface monitoring results, lateral positions are estimated robustly and are not sensitive to the velocity model. In this case study, we compare event location solutions from two catalogues of microseismic events; one from a downhole array and the second from a surface array of 1C geophone. Our results show that origin time can be reliably used to find matching events between the downhole and surface catalogues. The locations of the corresponding events display a systematic shift consistent with a poorly calibrated velocity model for downhole dataset. For this case study, locations derived from surface monitoring have less scatter in both vertical and horizontal directions.     

Changes in Components of the Water Balance in the Croatian Lowlands

Summary ?This paper deals with variations and trends in some components of the water balance: the soil water content; evaporation loss from the soil (from the surface and underlying layers); transpiration; ground water recharge and runoff. These components are calculated by means of the Palmer procedure. This analysis is based on data from Osijek, Croatia from this century (1900–1995). Besides the meteorological input parameters necessary for the water balance calculations, i.e. precipitation, temperature and relative humidity, the pedological characteristics of this area have also been taken into account. Fluctuations have been considered by means of the 11-year binomial filtered series and linear trends were tested by means of the Mann-Kendall rank test. For a closer look on the trends of water balance components, a progressive analysis of the time series was performed, too. The results show a significant increase in potential evapotranspiration and evapotranspiration and decrease in runoff and soil water content has occurred during the century. Received February 22, 1999/Revised August 3, 1999     

Competitive Adjustment and Advancement in Global Commodity Chains II: The Case of the Singapore Garment Industry

This is the second part of a contribution to the debate on the possibilities of leveraging economic globalisation — through incorporation in international production networks and global commodity chains (GCCs) — as a strategy for developing countries to industrialise and advance their position in the world economy. In the first part, we analysed the experience of the East Asian apparel industry and addressed the issues of sustaining positions, upgrading to more rewarding roles, and advancing to less dependent positions within production networks and GCCs. We developed an analytical framework at both the firm and industry levels, and subsequently identified alternative firm- and industry-level strategies and trajectories. The present paper deals with these issues in the context of the Singapore apparel industry. Based on extensive empirical research, we demonstrate that although the East Asian experience of upgrading and repositioning within the GCC is to some extent emulated in the case of the Singapore apparel industry, the outcomes have been less favourable in terms of the depth, extent and strength of these trajectories. The differing outcomes can be explained in terms of different (systemic) conditions in the Singapore business environment, including the agency of local players, the geography of sourcing networks, and the role of the state and prevailing business attitudes. Our conclusions merit continued attention in both research and policy circles on the development of capabilities at the firm level, and the role of local business and institutional environments in local industry development processes under globalisation.     

On forecasting of cyclogenesis in the West Mediterranean and other areas bounded by mountain ranges by baroclinic model

Summary Changes occurring in the thermal structure of air masses in the Alps area and the West Mediterranean region, during cyclonic developments in the Gulf of Genoa, have been considered. The effect of the height and the form of the mountain ranges as well as that of the nonadiabatic processes, on the deformation of the thickness pattern, have been discussed. These factors have been further taken into consideration for a modification of the equations used inEstoque's method for the baroclinic model, and a new method has been developed for forecasting cyclogenesis in regions bounded by mountain ranges.

---

Zusasmmenfassung Es werden die Änderungen der thermischen Struktur der Luftmassen im Alpengebiet und im westlichen Mittelmeerraum bei der Entwicklung von Zyklonen im Golf von Genua aufgezeigt und die Einflüsse von Höhe und Form der Gebirgszüge sowie von nichtadiabatischen Vorgängen auf die Deformation der relativen Topographie besprochen. Diese Faktoren werden auch zur Modifikation der in der Methode vonEstoque für ein baroklines Modell verwendeten Gleichungen in Betracht gezogen und damit wird eine neue Methode zur Vorhersage der Zyklogenese unter dem Einfluß von Gebirgszügen entwickelt.

---

Résumé On expose dans ce mémoire à quelles modifications la structure thermique des masses d'air est soumise dans la région des Alpes et dans le bassin occidental de la Méditerranée lors du développement de dépressions sur le Golfe de Gênes. On y discute en outre l'influence de la hauteur et de la forme des chaînes de montagnes ainsi que celle des processus non-adiabatiques sur des déformations de la topographie relative. Tous ces éléments sont pris en considération pour leur adapter les équations utilisées dans le modèle baroclined'Estoque. Par là, on accède à une nouvelle méthode de prévision de la cyclogénèse, méthode qui tient compte de l'influence des chaînes de montagnes.

---

---

With 5 Figures     

A note on the performance of the multiply-upstream semi-Lagrangian advection schemes for one-dimensional nonlinear momentum conservation equation

Summary Two-time-level multiply-upstream semi-Lagrangian schemes were examined in the case of the self-advecting, one-dimensional nonlinear momentum conservation equation. The shock formation process was analyzed. It is pointed out that the shocks cannot be created in the truncated systems satisfying the Pudykiewicz, Benoit and Staniforth criterion.The numerical integrations were restricted to 12 h. It was shown that, at least in the sub-CFL range, increased complexity of the scheme can compensate reduced horizontal resolution. A considerable sensitivity of the schemes with respect to the time step was detected. In the super-CFL mode, several windows on various time scales were found within which the Pudykiewicz, Benoit and Staniforth criterion was satisfied. The time step of 1.44 times the maximum time step allowed by the CFL criterion was used in the semi-Lagrangian runs.The super-CFL, semi-Lagrangian solutions were diverging progressively from the sub-CFL ones as the forecasts advanced. This was also reflected in the energy spectra.Unacceptably large energy losses were encountered in the super-CFL, semi-Lagrangian runs. Most of these losses could be explained by the reduced mean wind speed, i.e., the amplitude of the zero wavenumber wave. At the same time, the energy content in the shorter waves increased. In a more complex model, such a situation would resemble a loss of zonal, and an increase of transient eddy kinetic energy.A trajectory error measure was defined as the maximum absolute value of the distance between the actual arriving point of the particle originating at the estimated departure point, and the grid point assumed to be the arrival point in the semi-Lagrangian procedure. In contrast to the sub-CFL regime, this measure could reach a considerable fraction of the grid distance in the computations with the super-CFL time steps.In the physical system considered, the trajectories are determined only by the velocities at the departure points. With the semi-Lagrangian schemes the distances traveled by the particles are estimated on the basis of the velocities at the points downstream with respect to the departure points. Thus, unless the solution is smooth (in space and time) on the scales of the extrapolation distances/times, the upstream extrapolation does not promise the convergence of the solution.With 16 Figures