A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

Moment Tensor Inversion of Regional Phases: Application to a Mine Collapse

—?The aim of our study consists of analyzing potentially non-double-couple seismic events recorded at regional distances. In order to define the nature of the seismic source, a moment tensor inversion is carried out as this method is general enough not to initially constrain the source mechanism. In this paper we present an application to a seismic event induced by a mine collapse which occurred near the town of Halle in Germany. Because of its induced nature, many parameters such as the location and geometry of this seismic source are known. This information allows us to test the influence of inadequate propagation modeling on the moment tensor obtained from the inversion. Green's functions have been computed with the reflectivity method in a flat layered medium, using the European model EurID (Du et? al., 1998; Dufumier et al., 1997). From the inversion of P-wave seismograms recorded by the German Regional Seismic Network will, we obtained a source time function which can be decomposed into two subevents. The first one has a large isotropic part and a deviatoric mechanism with near vertical nodal planes. No volume change is observed for the second subevent, but a deviatoric component opposite of the first one. The addition of S-waves does not change the results of the inversion which are stable. Surface waves were not used because of their poor dispersion curves. Based on the moment tensor obtained from these inversions, the physical process at the source is compatible with a large cavity collapse.     

Preliminary seismic hazard maps of Slovenia

The preparation of the preliminary seismic hazard maps of the territory of Slovenia has been based on an expansion of the basic approach laid out by Cornell in 1968. Three seismic source models were prepared. Two of them are based mainly on the earthquake catalogue using the Poissonian probability model. A map of seismic energy release and a map of earthquake epicenter density are used to delineate seismic sources in these models. The geometry of the third model which is based on a rough estimate of seismotectonic setting is taken from the probabilistic seismic hazard analysis of a nuclear power plant in Slovenia. Published ground motion attenuation models based on strong motion records of recent strong earthquakes in Italy are used. Test maps for variable and uniform b-values are presented. The computer program, Seisrisk III, developed by the U.S. Geological Survey is used.     

Repeated Geodetic Measurements in the West Bohemia Seismoactive Region: Period 1993 - 1996

The western part of the Bohemian Massif is characterized by repeated occurrences of intraplate earthquake swarms. To study surface deformations of this anomalous region, a network covering about 2000 square kilometres for repeated geodetic measurements was established in 1993 - 1994. The positions of the individual points of the network were carefully picked with respect to local tectonic structure and earthquake foci distribution. GPS and precise levelling measurements were performed 1-2 times a year. The GPS data were processed by Bernese GPS software. No tendency to any displacement - either horizontal or vertical - of geological blocks was derived from the geodetic data for the period 1993 - 96. Only displacements of less than 5 mm/year in average could occur in the whole region in that period; larger displacements would have been revealed by our measurements.     

Argon isotope fractionation induced by stepwise heating

Noble gas isotopes are widely used to elucidate the history of the rocks in which they have been trapped, either from distinct reservoirs or by accumulation following radioactive decay. To extract noble gases from their host rocks, stepwise heating is the most commonly used technique to deconvolve isotopically different components, e.g., atmospheric, in situ radiogenic, or excess radiogenic from mantle or crustal reservoirs. The accurate determination of the isotopic composition of these different components is of crucial importance, e.g., for ages obtained by ⁴⁰Ar-³⁹Ar stepheating plateaus. However, diffusion theory-based model calculations predict that the stepwise thermal extraction process from mineral phases induces isotope fractionation and, hence, adulterates the original composition. Such effects are largely unconsidered, as they are small and a compelling experimental observation is lacking. We report the first unequivocal evidence for significant mass fractionation of argon isotopes during thermal extraction, observed on shungite, a carbon-rich Precambrian sedimentary rock. The degree of fractionation, as monitored by ³⁸Ar/³⁶Ar and ⁴⁰Ar/³⁶Ar ratios, very well agrees with theoretical predictions assuming an inverse square root dependence of diffusion coefficient and atomic mass, resulting in easier extraction of lighter isotopes. Hence, subatmospheric ⁴⁰Ar/³⁶Ar ratios obtained for argon extracted at low temperatures may not represent paleoatmospheric argon. Shungite argon resembles modern atmospheric composition, but constraints on the timing of trapping appear difficult to obtain, as shungites are multicomponent systems.In ⁴⁰Ar-³⁹Ar stepwise heating, the isotope fractionation effect could cause systematic underestimations of plateau ages, between 0.15 and 0.4% depending on age, or considerably higher if samples contain appreciable atmospheric Ar. The magnitude of this effect is similar to the presently achieved uncertainties of this increasingly precise dating technique. Our results also indicate the importance of thermally activated diffusion as a possible fractionation mechanism, e.g., for hydrothermal gas exhalations, or for carbonaceous carrier phases such as “Q” in meteorites that have been suggested as carriers of highly fractionated noble gas residues from the early solar nebula.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

Forecasting ozone concentrations in the east of Croatia using nonparametric Neural Network Models

Ozone is one of the most significant secondary pollutants with numerous negative effects on human health and environment including plants and vegetation. Therefore, more effort is made recently by governments and associations to predict ozone concentrations which could help in establishing better plans and regulation for environment protection. In this study, we use two Artificial Neural Network based approaches (MPL and RBF) to develop, for the first time, accurate ozone prediction models, one for urban and another one for rural area in the eastern part of Croatia. The evaluation of actual against the predicted ozone concentrations revealed that MLP and RBF models are very competitive for the training and testing data in the case of Kopa?ki Rit area whereas in the case of Osijek city, MLP shows better evaluation results with 9% improvement in the correlation coefficient. Furthermore, subsequent feature selection process has improved the prediction power of RBF network.