The Rožná uranium deposit (Bohemian Massif, Czech Republic): shear zone-hosted, late Variscan and post-Variscan hydrothermal mineralization

Three major mineralization events are recorded at the Rožná uranium deposit (total mine production of 23,000 t U, average grade of 0.24% U): (1) pre-uranium quartz-sulfide and carbonate-sulfide mineralization, (2) uranium, and (3) post-uranium quartz-carbonate-sulfide mineralization. (1) K–Ar ages for white mica from wall rock alteration of the pre-uranium mineralization style range from 304.5 ± 5.8 to 307.6 ± 6.0 Ma coinciding with the post-orogenic exhumation of the Moldanubian orogenic root and retrograde-metamorphic equilibration of the high-grade metamorphic host rocks. The fluid inclusion record consists of low-salinity aqueous inclusions, together with H₂O-CO₂-CH₄, CO₂-CH₄, and pure CH₄ inclusions. The fluid inclusion, paragenetic, and isotope data suggest that the pre-uranium mineralization formed from a reduced low-salinity aqueous fluid at temperatures close to 300°C. (2) The uraniferous hydrothermal event is subdivided into the pre-ore, ore, and post-ore substages. K–Ar ages of pre-ore authigenic K-feldspar range from 296.3 ± 7.5 to 281.0 ± 5.4 Ma and coincide with the transcurrent reorganization of crustal blocks of the Bohemian Massif and with Late Stephanian to Early Permian rifting. Massive hematitization, albitization, and desilicification of the pre-ore altered rocks indicate an influx of oxidized basinal fluids to the crystalline rocks of the Moldanubian domain. The wide range of salinities of fluid inclusions is interpreted as a result of the large-scale mixing of basinal brines with meteoric water. The cationic composition of these fluids indicates extensive interaction with crystalline rocks. Chlorite thermometry yielded temperatures of 260°C to 310°C. During this substage, uranium was probably leached from the Moldanubian crystalline rocks. The hydrothermal alteration of the ore substage followed, or partly overlapped in time, the pre-ore substage alteration. K–Ar ages of illite from ore substage alteration range from 277.2 ± 5.5 to 264.0 ± 4.3 Ma and roughly correspond with the results of chemical U–Pb dating of authigenic monazite (268 ± 50 Ma). The uranium ore deposition was accompanied by large-scale decomposition of biotite and pre-ore chlorite to Fe-rich illite and iron hydrooxides. Therefore, it is proposed that the deposition of uranium ore was mostly in response to the reduction of the ore-bearing fluid by interaction with ferrous iron-bearing silicates (biotite and pre-ore chlorite). The Th data on primary, mostly aqueous, inclusions trapped in carbonates of the ore substage range between 152°C and 174°C and total salinity ranges over a relatively wide interval of 3.1 to 23.1 wt% NaCl eq. Gradual reduction of the fluid system during the post-ore substage is manifested by the appearance of a new generation of authigenic chlorite and pyrite. Chlorite thermometry yielded temperatures of 150°C to 170°C. Solid bitumens that post-date uranium mineralization indicate radiolytic polymerization of gaseous and liquid hydrocarbons and their derivatives. The origin of the organic compounds can be related to the diagenetic and catagenetic transformation of organic matter in Upper Stephanian and Permian sediments. (3) K–Ar ages on illite from post-uranium quartz-carbonate-sulfide mineralization range from 233.7 ± 4.7 to 227.5 ± 4.6 Ma and are consistent with the early Tethys-Central Atlantic rifting and tectonic reactivation of the Variscan structures of the Bohemian Massif. A minor part of the late Variscan uranium mineralization was remobilized during this hydrothermal event.     

The Western Bohemia uppermost crust shear wave velocities from Love wave dispersion

We present a detailed study of Western Bohemia Love waves generated by blasts with an intention to estimate the uppermost crust structure for a more detailed layer distribution than previous studies have used. The use of short-period (4 s and shorter) Love waves represents a new approach in the studied region. Properties of multiple filtering as a tool of frequency–time analysis are discussed. A new method of selecting the dispersion ridges is introduced. Tests of filtering are provided by analyzing signals with analytically known dispersion. The isometric algorithm for the inversion problem is applied, the problem of non-uniqueness is discussed, and tests of reliability of the inversion are presented. During the inversion, the forward problem is solved by use of the matrix method. Six records of blasts from the Western Bohemia region are analyzed to separate the fundamental modes of Love wavegroups, and shear wave velocity distributions down to a depth of 3.0 km are inferred. Modal summation is used to compute synthetic velocigrams, which are compared to measured ones. The lateral heterogeneity of the region is discussed and the presented models are compared to those of previous studies and to the geological setting of the region.     

Dynamics of the Subsurface Temperature Field Caused by Mining-Related Groundwater Management

The subsurface temperature field was studied on a set of 46 borehole logs measured in the vicinity of uranium deposits in the Bohemian Cretaceous Basin. Vertical variations of the steady state temperature and the temperature gradient are governed by thermal conductivity which strongly varies in dependence on lithology. Large departures from undisturbed temperature detected in many holes are associated with uranium mining. A positive anomaly is observed in leaching fields where large amounts of acid are injected into the uranium-bearing Cenomanian. A negative anomaly is linked to the operation of hydraulic barriers which enclose the mining area and helps to contain pollution by pumping clean water into the Cenomanian aquifer. The spatial distribution of the observed temperature anomalies helps to map the migration of the fluids used in both processes. The temperature disturbance is propagated from the Cenomanian aquifer up through overlying impermeable sediments. The good fit of transient conductive models to the measured temperatures rules out heat advection and hence upward water flow towards and contamination of the upper Middle Turonian aquifer in the vicinity of the holes studied.     

The telluric field in a halfspace with a spherical inhomogeneity

Summary An exact solution of the distribution of the intensity of the telluric field in a halfspace with a spherical inhomogeneity has been obtained by solving the Laplace equation in bipolar coordinates. This exact solution is compared with the presently known approximate solution, obtained by the method of images, and the region of their coincidence has been determined. For the image solution theoretical relations have been derived for computing the anomaly of the magnetic field. Craphs depict many of the properties of the telluric and magnetic fields on the surface of the considered halfspace for various parameters of the given problem.Dedicated to 90th Birthday of Professor Frantiek Fiala     

Imprint of the Atlantic multi-decadal oscillation and Pacific decadal oscillation on southwestern US climate: past,present, and future

The surface air temperature increase in the southwestern United States was much larger during the last few decades than the increase in the global mean. While the global temperature increased by about 0.5 °C from 1975 to 2000, the southwestern US temperature increased by about 2 °C. If such an enhanced warming persisted for the next few decades, the southwestern US would suffer devastating consequences. To identify major drivers of southwestern climate change we perform a multiple-linear regression of the past 100 years of the southwestern US temperature and precipitation. We find that in the early twentieth century the warming was dominated by a positive phase of the Atlantic multi-decadal oscillation (AMO) with minor contributions from increasing solar irradiance and concentration of greenhouse gases. The late twentieth century warming was about equally influenced by increasing concentration of atmospheric greenhouse gases (GHGs) and a positive phase of the AMO. The current southwestern US drought is associated with a near maximum AMO index occurring nearly simultaneously with a minimum in the Pacific decadal oscillation (PDO) index. A similar situation occurred in mid-1950s when precipitation reached its minimum within the instrumental records. If future atmospheric concentrations of GHGs increase according to the IPCC scenarios (Solomon et al. in Climate change 2007: working group I. The Physical Science Basis, Cambridge, 996 pp, 2007), climate models project a fast rate of southwestern warming accompanied by devastating droughts (Seager et al. in Science 316:1181–1184, 2007; Williams et al. in Nat Clim Chang, 2012). However, the current climate models have not been able to predict the behavior of the AMO and PDO indices. The regression model does support the climate models (CMIP3 and CMIP5 AOGCMs) projections of a much warmer and drier southwestern US only if the AMO changes its 1,000 years cyclic behavior and instead continues to rise close to its 1975–2000 rate. If the AMO continues its quasi-cyclic behavior the US SW temperature should remain stable and the precipitation should significantly increase during the next few decades.     

New Techniques and Limitations of Light Curve Analysis

The solution of light curves of eclipsing binaries provides the simplest way of measuring the basic physical parameters of stars. The development of computers and numerical methods enables us to perform an analysis of large amount of photometric data by comparing it to more and more sophisticated physical models which are also suitable for stellar systems with different peculiarities. Nevertheless, some parameters of stellar systems are poorly determined by light curves alone, and simultaneous fitting together with other data is preferable. The progress of instrumentation not only improves the precision and resolution of data, but also overcomes the classical gaps between stellar systems observable as visual, photometric and spectroscopic binaries. This opens new problems in the development of techniques of data analysis, which can be treated as a generalization of light curve modeling. Some of these directions are outlined here. In particular, the disentangling of spectra with line-profile variability due to proximity effects is discussed.     

Structure and stationarity of quasi-perpendicular shocks: Numerical simulations

This paper presents an overview of numerical simulation studies of fast collisionless shocks and compares these simulation results with observations of the Earth's bow shock and theoretical works. Especially, we review the structure and stationarity of the supercritical quasi-perpendicular shocks. In situ observations indicate that these shocks are generally quasi-stationary whereas full particle simulations as well as hybrid simulations often present a strong nonstationary behavior, a shock self-reformation. The simulation results, along with theoretical and observational works, suggest that the classical models of the quasi-stationary structure generated by reflected protons or by dispersive whistlers are not generally applicable for the supercritical quasi-perpendicular shocks and other phenomena are to be included into the model to ensure the observed quasi-stationarity: The role of a small scale turbulence and shock ripples is investigated. The downstream turbulence and the electron dynamics in the quasi-perpendicular shocks are also discussed.     

A combined statistical approach for evaluation of the effects of land use,agricultural and urban activities on stream water chemistry in small tile-drained catchments of south Bohemia,Czech Republic

This work evaluates the changes of nitrate-nitrogen (NO₃-N), ammonium-nitrogen (NH₄-N), total phosphorus (P) and chemical oxygen demand (COD) concentrations in stream waters as related to the land use/land cover (LULC) alterations within eight small (5–39 km²) tile-drained catchments in the southern part of The Czech Republic in the period 1993–2010, when massive grassing of arable land took place. The robust tools of seasonal Mann–Kendall trend test and LOcally WEighted Scatterplot Smoothing methods were employed to reveal trends of the monitored parameters with adjustment to hydrology. Using principal component analysis and multiple regressions, statistically significant factors with highest impacts on the assessed water quality parameters were identified. Besides indicators of LULC changes in the catchments and their various zones, information of built tile drainage systems were used along with factors reflecting point pollution sources such as the population number, sewerage type and proximity to a watercourse, effectiveness of wastewater treatment, and number of livestock units. The change in LULC was essential only for NO₃-N concentrations, when grassing of arable land, presence of water ponds, areas of permanent cultures and also areas of drained land explained up to 90.6 % NO₃-N variability and nitrate-nitrogen concentrations showed a significantly decreasing trend in all monitored catchments during the evaluated period. LULC changes within infiltration vulnerable zones were discovered as less important for the assessed water quality parameters compared to LULC changes in the whole catchment area. However, for NH₄-N, P and COD, the results did not enable a definite quantification of the effects of LULC changes. The influence of non-point pollution sources on these parameters was revealed as uncertain and was heavily overshadowed by point sources, in particular by wastewater management, and livestock numbers, although the proportion of arable land in tile drainage subcatchments was discovered fundamental in case of the COD. The increasing numbers of livestock, population, and changes in sewage treatment led in some catchments to significant worsening of water quality. Achieved findings may be critical for supporting water quality policy and management decisions.     

Electron transport and precipitation at Mercury during the MESSENGER flybys: Implications for electron-stimulated desorption

To examine electron transport, energization, and precipitation in Mercury's magnetosphere, a hybrid simulation study has been carried out that follows electron trajectories within the global magnetospheric electric and magnetic field configuration of Mercury. We report analysis for two solar-wind parameter conditions corresponding to the first two MESSENGER Mercury flybys on January 14, 2008, and October 6, 2008, which occurred for similar solar wind speed and density but contrasting interplanetary magnetic field (IMF) directions. During the first flyby the IMF had a northward component, while during the second flyby the IMF was southward. Electron trajectories are traced in the fields of global hybrid simulations for the two flybys. Some solar wind electrons follow complex trajectories at or near where dayside reconnection occurs and enter the magnetosphere at these locations. The entry locations depend on the IMF orientation (north or south). As the electrons move through the entry regions they can be energized as they execute non-adiabatic (demagnetized) motion. Some electrons become magnetically trapped and drift around the planet with energies on the order of 1–10 keV. The highest energy of electrons anywhere in the magnetosphere is about 25 keV, consistent with the absence of high-energy (>35 keV) electrons observed during either MESSENGER flyby. Once within the magnetosphere, a fraction of the electrons precipitates at the planetary surface with fluxes on the order of 10⁹ cm⁻² s⁻¹ and with energies of hundreds of eV. This finding has important implications for the viability of electron-stimulated desorption (ESD) as a mechanism for contributing to the formation of the exosphere and heavy ion cloud around Mercury. From laboratory estimates of ESD ion yields, a calculated ion production rate due to ESD at Mercury is found to be on par with ion sputtering yields.