The BanskáŠtiavnica ore district: relationship between metallogenetic processes and the geological evolution of a stratovolcano

The Banská?tiavnica ore district is in the central zone of the largest stratovolcano in the Central Slovakia Neogene Volcanic Field, which is situated at the inner side of the Carpathian arc over the Hercynian basement with the Late Paleozoic and Mesozoic sedimentary cover. Volcanic rocks of the High-K orogenic suite are of the Badenian through Pannonian age (16.5–8.5?Ma). Their petrogenesis is closely related to subduction of flysch belt oceanic basement underneath the advancing Carpathian arc and to back-arc extension processes. The stratovolcano includes a large caldera 20?km in diameter and a late-stage resurgent horst in its centre, exposing a basement and extensive subvolcanic intrusive complex. The following stages have been recognized in the evolution of the stratovolcano: (1)?formation of a large pyroxene/hornblende-pyroxene andesite stratovolcano; (2)?denudation, emplacement of a diorite intrusion; (3) emplacement of a large granodiorite bell-jar pluton within the basement; (4) emplacement of granodiorite/quartz-diorite porphyry stocks and dyke clusters around the pluton; (5) caldera subsidence and its filling by biotite-hornblende andesite volcanics, emplacement of quartz-diorite porphyry sills and dykes at the subvolcanic level; (6)?renewed activity of andesites from dispersed centres on slopes of the volcano; (7) uplift of a resurgent horst accompanied by rhyolite volcanics and granite porphyry dykes. The following types of ore deposits (mineralizations) have been identified in the Banská?tiavnica ore district: 1. Quartz-pyrophyllite-pyrite high-sulphidation system at ?obov, related to the diorite intrusion. 2. Magnetite skarn deposits and occurrences?at contacts of the granodiorite pluton with Mesozoic carbonate rocks. Magnetite ores occur as lenses in the calcic skarns. 3.?Stockwork/disseminated base metal deposit along an irregular network of fractures in apical parts of the granodiorite pluton and in remnants of basement rocks. Mineral paragenesis is simple, with leading sphalerite and galena and minor chalcopyrite and pyrite. In overlying andesites the mineralization is accompanied by metasomatic quartzites and argillites with pyrophyllite, kaolinite, illite and pyrite. 4. Porphyry/skarn copper deposits and occurrences related to granodiorite/quartz-diorite porphyry dyke clusters and stocks around the granodiorite intrusion. The mineralized zone is represented by accumulations of chalcopyrite in exo- and endo-skarns, usually of the magnesian type affected by serpentinization. Besides chalcopyrite, pyrhotite, minor bornite, chalcosite, tennantite and magnetite, rare molybdenite and gold are present. The alteration pattern around productive intrusions includes an external zone of propylitization, a zone of argillitic alteration (kaolinite – illite – pyrite) and an internal zone of phyllic alteration (quartz – sericite – pyrite). Biotitization is rare and limited to porphyry intrusions. 5. Intrusion related “mesothermal” gold deposit in an andesitic environment just above the granodiorite intrusion. Gold of high fineness with base metal mineralization is contained in brecciated and/or banded quartz veins of subhorizontal orientation, parallel to the surface of granodiorite pluton. At least the first phase of mineralization is older than quartz-diorite porphyry sills, which separate granodiorite and blocks of mineralized andesite. 6. Hot spring type advanced argillic systems in the caldera filling. Silicites and opalites accompanied by kaolinite, alunite and pyrite grade downward into smectite dominated argillites. 7. Vein type epithermal precious/base metal deposits and occurrences as a result of the long lasting interaction among structural evolution of the resurgent horst and evolving hydrothermal system, extensive intrusive complex and deep seated siliceous magma chamber serving as heat and magmatic fluid source. Three types of epithermal veins occur in a zonal arrangement: (a) base metal veins ± Au with transition to Cu?±?Bi mineralization at depth in the east/central part of the horst, (b)?Ag – Au veins with minor base metal mineralization and (c) Au – Ag veins located at marginal faults of the horst. Isotopic composition of oxygen and hydrogen in hydrothermal fluids indicate mixing of magmatic and meteoric component (with generally increasing proportion of meteoric component towards younger mineralization periods?). Veins are accompanied by zones of silicification, adularization and sericitization, indicating a low sulphidation environment. 8.?Replacement base metal mineralization of a limited extent in the Mesozoic carbonate rocks next to sulphide rich epithermal base metal veins.     

Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW”

The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.     

Prediction of magnitude of the largest potentially induced seismic event

We propose a method for determining the possible magnitude of a potentially largest induced seismic event derived from the Gutenberg–Richter law and an estimate of total released seismic moment. We emphasize that the presented relationship is valid for induced (not triggered) seismicity, as the total seismic moment of triggered seismicity is not bound by the injection. The ratio of the moment released by the largest event and weaker events is determined by the constants a and b of the Gutenberg–Richter law. We show that for a total released seismic moment, it is possible to estimate number of events greater than a given magnitude. We determine the formula for the moment magnitude of a probable largest seismic event with one occurrence within the recurrence interval (given by one volumetric change caused by mining or injecting). Finally, we compare theoretical and measured values of the moment magnitudes of the largest induced seismic events for selected geothermal and hydraulic fracturing projects.     

Dust ejection from planetary bodies by temperature gradients: Laboratory experiments

Laboratory experiments show that dusty bodies in a gaseous environment eject dust particles if they are illuminated. We find that even more intense dust eruptions occur when the light source is turned off. We attribute this to a compression of gas by thermal creep in response to the changing temperature gradients in the top dust layers. The effect is studied at a light flux of 13 kW/m² and 1 mbar ambient pressure. The effect is applicable to protoplanetary disks and Mars. In the inner part of protoplanetary disks, planetesimals can be eroded especially at the terminator of a rotating body. This leads to the production of dust which can then be transported towards the disk edge or the outer disk regions. The generated dust might constitute a significant fraction of the warm dust observed in extrasolar protoplanetary disks. We estimate erosion rates of about 1 kg s^?1 for 100 m parent bodies. The dust might also contribute to subsequent planetary growth in different locations or on existing protoplanets which are large enough not to be susceptible to particle loss by light induced ejection. Due to the ejections, planetesimals and smaller bodies will be accelerated or decelerated and drift outward or inward, respectively. The effect might also explain the entrainment of dust in dust devils on Mars, especially at high altitudes where gas drag alone might not be sufficient.     

Thomsen's parameters from p‐wave measurements in a spherical sample‡

The aim of this paper is to understand the seismic anisotropy of the overburden shale in an oilfield in the North West Shelf of Western Australia. To this end, we first find the orientation of the symmetry axis of a spherical shale sample from measurements of ultrasonic P‐wave velocities in 132 directions at the reservoir pressure. After transforming the data to the symmetry axis coordinates, we find Thomsen's anisotropy parameters δ and ? using these measurements and measurements of the shear‐wave velocity along the symmetry axis from a well log. To find these anisotropy parameters, we use a very fast simulated re‐annealing algorithm with an objective function that contains only the measured ray velocities, their numerical derivatives and the unknown elasticity parameters. The results show strong elliptical anisotropy in the overburden shale. This approach produces smaller uncertainty of Thomsen parameter δ than more direct approaches.     

Validation of 3D synthetic seismograms based on the ray-Born approximation

The first-order Born approximation is a weak scattering perturbation method which is a powerful tool. The combination of the Born approximation and the ray theory enables to extend the applicability of the ray theory in terms of the required smoothness of the model and ensures faster computations than with, e.g., the finite difference method. We are motivated to describe and explain the effects of the numerical discretization of the Born integral on the resulting seismograms. We focus on forward modelling and study the cases in which perturbation from the background model contains the interface. We restrict ourselves to isotropic models that contain two homogeneous layers. We compare the 2D and 3D ray-based Bornapproximation seismograms with the ray-theory seismograms. The Born seismograms are computed using a grid of finite extent. We anticipate that the computational grid should contain an appropriate number of gridpoints, otherwise the seismogram would be inaccurate. We also anticipate that the limited size of the computational grid can cause problems. We demonstrate numerically that an incorrect grid can produce significant errors in the amplitude of the wave, or it can shift the seismogram in time. Moreover, the grid boundaries work as interfaces, where spurious waves can be generated. We also attempt to explain these phenomena theoretically. We give and test the options of removing the spurious waves. We show that it is possible to compute the Born approximation in a sparser grid, if we use elastic parameters averaged from some dense grid.     

Small-scale patterns of meiofauna in a bryophyte covered tufa barrier (Plitvice Lakes,Croatia)

Meiobenthic fauna is complex and essential part of the stream benthos and it is known as intermediaries from bacteria to higher trophic levels. However, the dynamics and mechanisms regulating this community at small-scale largely have been neglected. This study was carried out to determine meiofaunal dynamics in a small scale-patterns influenced by flow velocity and other abiotic and biotic factors. We examined meiofauna within bryophytes on a tufa barrier in Plitvice Lakes National Park using a detailed taxonomic approach of various meiofaunal groups and their feeding guilds. We choose three microhabitats (slow, medium, fast) differing significantly in flow velocity above bryophytes. Bdelloid rotifers were the most abundant group in microhabitat with highest flow velocity, while in other two microhabitats nematodes and monogonont rotifers prevailed in abundance. Data on environmental variables and main meiofaunal taxa and feeding guilds were analyzed using redundancy analysis. This analysis indicated that microfilter feeding guild (e.g. bdelloid rotifers) was strongly affected by interaction of flow velocity and POM fractions. Other feeding guilds were influenced by temperature, oxygen and/or pH and did not prefer high flow velocity. Suction-feeder nematodes and microfilter-feeder rotifers were dominant on temporal and spatial scale, indicating their good adaptations on frequently disturbed conditions that prevailed on bryophyte covered tufa barrier. Our results provide comprehensive survey of diversity, density as well as trophic structure of meiofauna in aquatic bryophytes. Differences in meiofaunal composition and density between three microhabitats suggest that the meiofauna is relevant indicator of environmental changes even at small-scale pattern.