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.     

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.