An inelastic stress integration algorithm for a rock mass containing sets of discontinuities

The mechanical behaviour of a rock mass is described in this paper using a method that combines experimental and numerical techniques. A new discontinuity (joint) constitutive model is proposed based on the non-associative yield law. Constitutive relationships for an anisotropic rock mass are determined by introducing an interaction between sets of discontinuities and the monolith in which they occur. The model is based on the concept that deformation of the monolith and its associated sets of discontinuities are coupled. An experimental investigation of the concrete foundations of the Bogovina Dam on the Crni Timok River in Serbia has been used to verify the presented numerical model. The tests included discontinuity and monolith sample investigations and large-scale field tests.     

Emission sources and health risk assessment of polycyclic aromatic hydrocarbons in ambient air during heating and non-heating periods in the city of Novi Sad,Serbia

The paper presents the assessment of emission sources and health risk of 16 PAHs in the city of Novi Sad, Serbia, with developed oil, food, chemical and cement industry. Although the observed region has been previously studied with regard to the occurrence of PAHs in ambient air, the sources of emission and health risk using multivariate statistical methods have never been studied before. The research was conducted during the heating and non-heating periods at three sampling sites located in industrial, residential and high-traffic area. The samples were collected by high-volume air sampler and analysed using the GC/MS system. The data were processed using the positive matrix factorization (PMF) and cluster analysis. To assess the cancer-risk of PAHs, potency equivalency factor approach was used. The results of PMF analysis have indicated that significant sources of PAH emissions are industrial processes and diesel combustion during the non-heating, and coal combustion, heating plants and diesel combustion during the heating season. Cluster analysis of experimental data has confirmed the results obtained using the PMF method. Estimated values of lifetime lung cancer risk were higher than those set by the World Health Organization and US Environmental Protection Agency which indicates that cancer risk due to exposure to PAHs in ambient air exists in both time periods on the territory of Novi Sad.     

The Western Carpathians—interaction of Hercynian and Alpine processes

The complicated structural and rheologic properties of Western Carpathian lithosphere reflect the complex geodynamic history of the Carpathian orogen. Based on critical analysis of earlier models, new interpolation of existing geophysical data and results of integrated modelling, a new map of the lithosphere thickness for the Carpathian–Pannonian region has been constructed. The map allows for the distinction of a frontal orogen collision zone in the NE (from increased lithosphere thickness) as well as a zone of oblique collision with the Bohemian Massif in the West, where lithosphere is not significantly thickened. The MOHO discontinuity beneath the Western Carpathian hinterland (Danube and East Slovak Basins), as defined by deep reflection seismic profiling, is relatively shallow. This probably reflects recent crustal extension related to oblique collision between the European plate and the ALCAPA block and an increase of the asthenospheric updoming from the Middle Miocene onward.Crustal thickness reflects the combined effects of deep-seated orogenic processes and mantle thermal evolution beneath the Pannonian Basin system. In this study, we focus particularly the structures of: (1) the Late Alpine collision and Neogene back arc basin development, including deep-seated contacts between colliding plates, a zone of slab detachment, the compressional accretionary wedge of the Outer Western Carpathian Flysch Belt, and extensional structures produced by subduction rollback and asthenosphere upwelling; (2) Early Alpine structures related to Cretaceous thrust-stacking, including subhorizontal reflection packages (interpreted as multi-generational extensional structures), the underplated intra-Penninic (Oravic) continental ribbon, and ophiolite traces of the Meliatic oceanic suture; and (3) north-dipping reflectors interpreted as remnant Hercynian lithotectonic fragments with opposed vergency to the subducted Alpine units.     

Petrology, geochemistry and zircon age for redwitzite at Abertamy, NW Bohemian Massif (Czech Republic): tracing the mantle component in Late Variscan intrusions

A small body of mafic texturally and compositionally varied igneous intrusive rocks corresponding to redwitzites occurs at Abertamy in the Western pluton of the Krušné hory/Erzgebirge granite batholith (Czech Republic). It is enclosed by porphyritic biotite granite of the older intrusive suite in the southern contact zone of the Nejdek-Eibenstock granite massif. We examined the petrology and geochemistry of the rocks and compared the data with those on redwitzites described from NE Bavaria and Western Bohemia.The redwitzites from Abertamy are coarse- to medium-grained rocks with massive textures and abundant up to 2 cm large randomly oriented biotite phenocrysts overgrowing the groundmass. They are high in MgO, Cr and Ni but have lower Rb and Li contents than the redwitzites in NE Bavaria. Compositional linear trends from redwitzites to granites at Abertamy indicate crystal fractionation and magma mixing in a magma chamber as possible mechanisms of magma differentiation. Plots of MgO versus SiO₂, TiO₂, Al₂O₃, FeO, CaO, Na₂O, and K₂O indicate mainly plagioclase and orthopyroxene fractionation as viable mechanisms for in situ differentiation of the redwitzites.The porphyritic biotite monzogranite enclosing the redwitzite is the typical member of the early granitic suite (Older Intrusive Complex, OIC ) with strongly developed transitional I/S-type features. The ages of zircons obtained by the single zircon Pb-evaporation method suggest that the redwitzites and granites at Abertamy originated during the same magmatic period of the Variscan plutonism at about 322 Ma.The granitic melts have been so far mainly interpreted to be formed by heat supply from a thickened crust or decompression melting accompanying exhumation and uplift of overthickened crust in the Krušné hory/Erzgebirge due to a previous collisional event at ca. 340 Ma. The presence of mafic bodies in the Western pluton of the Krušné hory/Erzgebirge batholith confirms a more significant role of mantle-derived mafic magmas in heating of the sources of granitic melts than previously considered.     

Changes of parameters of magnetic anisotropy of foliated rocks under pressure

Summary The changes of magnetic anisotropy under pressure were studied on a set of rocks with marked macroscopic stratification, collected from the Rjvíz borehole. The anisotropy parameters are relatively very stable under directional pressure, acting parallel with the rocks' bedding, and, on the contrary, significantly unstable under pressure acting perpendicular to the bedding. Systematically different magnitudes of stress sensitivity coefficients, , for parallel directional susceptibility were observed in both cases in the same types of rocks. The main cause is probably the different capability of transferring external stress to ferrimagnetics via the stratified non-magnetic matrix. With regard to magnetomechanical phenomena, therefore, a singel value of the stress sensitivity coefficient is insufficient to characterize rocks with a markedly anisotropic matrix.     

Geological interpretation of gravity profiles through the Karlovy Vary granite massif (Czech Republic)

We examined the shape of the Late Variscan Karlovy Vary granite massif located south of the Ohre/Eger graben in Northern Bohemia by reinterpretation of existing gravity data on two perpendicular profiles. The granite body of about 360 km² total outcrop size has the elongation ratio 0.35 with the major axis trending NE-SW. The SW part of the body was crossed in the nineties by the seismic profile 9HR which localized the bottom of granites in a depth of about 10 km. We used this value as a reference datum in our gravity profiles. We positioned one of our profiles along the seismic profile 9HR and the other one perpendicularly, i.e. parallel with the elongation of the outcrop surface. We interpret the shape of the main granite body in the vicinity of Karlovy Vary as a continuous desk whose floor is horizontal (or subhorizontal) and varies along its whole extension about a depth of 10 km. This thickness is approximately identical with that of the Saxothuringian nappes imaged by seismic reflection. The near surface upper contact of the granite body is mildly inclined, and outward dipping. It changes to steep sides or inward inclined contacts in deeper levels. The Lesny-Lysina (Kynžvart) massif is a separate granite body about 324 km thick, not continuously connected with the main Karlovy Vary massif. The gravity curve suggests that granites often enclose in their endocontact large blocks of country metasediments or metabasites the existence of which is partly evidenced by their outcrops outside the line of the profile. The granite body is found density-homogenous. Minor density differences between granite varieties are caused mainly by more intense hydrothermal alterations in younger suite granites. We interpret vertical conduits for the ascent of granitic magmas to be parallel to the Jáchymov-Gera and Ohře (Eger) lineaments or the Mariánské Lázně fault zone as indicated by the elongation of some outcrops. However, they are not clearly imaged from the gravity data. The effect of the depression of the Sokolov basin along the faults parallel with the Ohře (Eger) lineament is shallow and it is not indicated by any change in the floor depth of the granite body. Comparison of the seismicity distribution suggests that the hypocenters occur mostly outside of the granite bodies or near their contact with the country rock.     

Depth-Recursive Tomography Along the Eger Rift Using the S01 Profile Refraction Data: Tested at the KTB Super Drilling Hole,Structural Interpretation Supported by Magnetic,Gravity and Petrophysical Data

The refraction data from the SUDETES 2003 experiment were used for high-resolution tomography along the profile S01. The S01 profile crosses the zone Erbendorf-Vohenstrauss (ZEV) near the KTB site, then follows the SW–NE oriented Eger Rift in the middle part and continues toward the NE across the Elbe zone and the Sudetic structures as far as the Trans-European Suture Zone. To get the best resolution in the velocity image only the first arrivals of Pg waves with minimum picking errors were used. The previous depth-recursive tomographic method, based on Claerbout’s imaging principle, has been adapted to perform the linearized inversions in iterative mode. This innovative DRTG method (Depth-Recursive Tomography on Grid) uses a regular system of refraction rays covering uniformly the mapped domain. The DRTG iterations yielded a fine-grid velocity model with a required level of RMS travel-time fit and the model roughness. The travel-time residuals, assessed at single depth levels, were used to derive the statistical lateral resolution of “lens-shaped” velocity anomalies. Thus, for the 95% confidence level and 5% anomalies, one can resolve their lateral sizes from 15 to 40 km at the depths from 0 to 20 km. The DRTG tomography succeeded in resolving a significant low-velocity zone (LVZ) bound to the Franconian lineament nearby the KTB site. It is shown that the next optimization of the model best updated during the DRTG iterations tends to a minimum-feature model with sweeping out any LVZs. The velocities derived by the depth-recursive tomography relate to the horizontal directions of wave propagation rather than to the vertical. This was proved at the KTB site where pronounced anisotropic behavior of a steeply tilted metamorphic rock complex of the ZEV unit has been previously determined. Involving a ~7% anisotropy observed for the “slow” axis of symmetry oriented coincidentally in the horizontal SW–NE direction of the S01 profile, the DRTG velocity model agrees fairly well with the log velocities at the KTB site. Comparison with the reflectivity map obtained on the reflection seismic profile KTB8502 confirmed the validity of DRTG velocity model at maximum depths of ~16 km. The DRTG tomography enabled us to follow the relationship of major geological units of Bohemian Massif as they manifested in the obtained P-wave velocity image down to 15 km. Although the contact of Saxothuringian and the Teplá-Barrandian Unit (TBU) is collateral with the S01 profile direction, several major tectonic zones are rather perpendicular to the Variscan strike and so fairly imaged in the S01 cross-section. They exhibit a weak velocity gradient of sub-horizontal directions within the middle crust. In particular, the Moldanubian and TBU contact beneath the Western Krušné hory/Erzgebirge Pluton, the buried contact of the Lusatia unit and the TBU within the Elbe fault zone were identified. The maxima on the 6,100 ms⁻¹ isovelocity in the middle crust delimitated the known ultrabasic Erbendorf complex and implied also two next ultrabasic massifs beneath the Doupovské hory and the České středohoří volcanic complexes. The intermediate mid-crustal P-wave velocity lows are interpreted as granitic bodies. The presented geological model is suggested in agreement with available gravity, aeromagnetic and petrophysical data.     

The unique environment of the most acidified permanently meromictic lake in the Czech Republic

Changes in the water properties and biological characteristics of the highly acidic Hromnice Lake (Western Bohemia) were investigated. This 110-year-old lake, formed as a consequence of the mining of pyritic shales, is permanently meromictic. Two chemoclines separate an extremely acidic (pH ～ 2.6) mixolimnion from a metal-rich anoxic monimolimnion. The absence of spring mixolimnetic turnover due to ice melting and very slow heat propagation through the chemocline with a 6-month delay were observed. Extreme mixolimnetic oxygen maxima (up to 31 mg l^?1) in phosphorus-rich lake (PO₄^3? up to 1.6 mg l^?1) well correlated with outbursts of phytoplankton. Phytoplankton consist of several acido-tolerant species of the genera Coccomyxa, Lepocinclis, Chlamydomonas and Chromulina. Surface phytoplankton biomass expressed as chlorophyll-a varies from 2 to 140 μg l^?1. Multicellular zooplankton are almost absent with the exception of Cephalodella acidophila, a small rotifer occurring in low numbers. Large red larvae of the midge Chironomus gr. plumosus were found at the bottom close to the shore, with larvulae in the open water. Developmental stages (protonemata) of a moss, resembling filamentous algae, dwell in the otherwise plant-free littoral zone.     

Climate Change Impacts and Adaptation Strategies in Spring Barley Production in the Czech Republic

The crop model CERES-Barley was used to assess the impacts of increased concentration of atmospheric CO₂ on growth and development of the most important spring cereal in Central and Western Europe, i.e., spring barley, and to examine possible adaptation strategies. Three experimental regions were selected to compare the climate change impacts in various climatic and pedological conditions. The analysis was based on multi-year crop model simulations run with daily weather series obtained by stochastic weather generator and included two yield levels: stressed yields and potential yields. Four climate change scenarios based on global climate models and representing 2 × CO₂ climate were applied. Results: (i) The crop model is suitable for use in the given environment, e.g., the coefficient of determination between the simulated and experimental yields equals 0.88. (ii) The indirect effect related to changed weather conditions is mostly negative. Its magnitude ranges from ?19% to +5% for the four scenarios applied at the three regions. (iii) The magnitude of the direct effect of doubled CO₂ on the stressed yields for the three test sites is 35–55% in the present climate and 25–65% in the 2 × CO₂ climates. (iv) The stressed yields would increase in 2 × CO₂ conditions by 13–52% when both direct and indirect effects were considered. (v) The impacts of doubled CO₂ on potential yields are more uniform throughout the localities in comparison with the stressed yields. The magnitude of the indirect and direct effects ranges from ?1 to ?9% and from +31 to +33%, respectively. Superposition of both effects results in 19–30% increase of the potential yields. (vi) Application of the earlier planting date (up to 60 days) would result in 15–22% increase of the yields in 2 × CO₂ conditions. (vii) Use of a cultivar with longer vegetation duration would bring 1.5% yield increase per one extra day of the vegetation season. (viii) The initial water content in the soil water profile proved to be one of the key elements determining the spring barley yield. It causes the yields to increase by 54–101 kg.ha^?1 per 1% increase of the available soil water content on the sowing day.