Geochronology and geochemistry of eclogites from the Mariánské Lázně Complex,Czech Republic: Implications for Variscan orogenesis

The Mariánské Lázn complex (MLC) is located in the Bohemian Massif along the north-western margin of the Teplá-Barrandian microplate and consists of metagabbro, amphibolite and eclogite, with subordinate amounts of serpentinite, felsic gneiss and calcsilicate rocks. The MLC is interpreted as a metaophiolite complex that marks the suture zone between the Saxothuringian rocks to the north-west and the Teplá-Barrandian microplate to the south-east. Sm-Nd geochronology of garnet-omphacite pairs from two eclogite samples yields ages of 377±7, and 367±4 Ma. Samples of eclogite and amphibolite do not define a whole rock Sm-Nd isochron, even though there is a large range in Sm/Nd ratio, implying that the suite of samples may not be cogenetic. Eclogites do not have correlated _Nd values and initial ⁸⁷Sr/⁸⁶Sr ratios. Five of the eight eclogite samples have high _Nd values (+10.2 to +7.1) consistent with derivation from a MORB-like source, but variable ⁸⁷Sr/⁸⁶Sr ratios (0.7033 to 0.7059) which probably reflect hydrothermal seawater alteration. Three other eclogite samples have lower _Nd values (+ 5.4 to –0.8) and widely variable ⁸⁷Sr/⁸⁶Sr ratios (0.7033 to 0.7096). Such low _Nd values are inconsistent with derivation from a MORB, source and may reflect a subduction or oceanic island basalt component in their source. The MLC is an important petrotectonic element in the Bohemian Massif, providing evidence for Cambro-Ordovician formation of oceanic crust and interaction with seawater, Late Devonian (Frasnian-Famennian) high- and medium-pressure metamorphism related to closure of a Saxothuringian ocean basin, Early Carboniferous (Viséan) thrusting of the Teplá terrane over Saxothuringian rocks and Late Viséan extension.     

Comparison between a radiative wind model for be stars and high angular resolution data from the GI2T interferometer

I present theoretical line profiles and intensity maps from an axi-symmetric radiative wind model from a rapidly rotating Be star. The introduction of a viscosity parameter in the latitude-dependent hydrodynamic code enables us to consider the effects of the viscous force in the azimuthal component of momentum equations (Araújo et al. 1994). Both velocity field and density law derived from the hydrodynamic equations have been used for solving the statistical equilibrium equations. By adopting the Sobolev approximation, we could easily obtain a good estimate of both electronic density and hydrogen level populations throughout the envelope. The numerical calculation was performed for parameters characterisic of the Be star Cassiopeiae.     

Structure of the satellitary accretion disk of Saturn

A detailed computation on the equilibrium structure of an accretion disk around Saturn from which the regular satellites presumably originated is reported. Such a disk is the predecessor of the self-dissipating disk that is formed when the mass infall stops (Cassen and Moosman, 1981, Icarus48, 353–376). When determining the disk structure local energy balance was assumed. Convention was taken into account by introducing local energy dissipation and, in an approximate manner, sonic convection. Changes in the disk structure were investigated by varying the free parameters, i.e., the external flux from both the protosun and the protoplanet, the abundance of dust and the strength of turbulence. It has been verified that the external energy flux does not play an important role in the evolution of the disk structure. Models characterized by either longer times (?3 10³ year) or a noticeable depletion of condensable elements (10^?2 times less than the solar value) have a total mass of the order of 0.34?0.1 times the mass of the regular satellites increased by the mass of the light elements. Low turbulence models (Reynolds critical number $\text{Re}{}^1\text{= 150}$) are characterized approximately by a total mass twice as large the mass of the regular satellites. All the studied models present a temperature distribution that allows the condensation of iron, silicate, and, in the outer regions, ice grains. All models but the one with 10^?2 of the solar value of condensable elements are characterized by a wide convective region that contains the formation zone of the regular satellites.     

Using numerical modelling to evaluate tsunami hazard near the Kuril Island

A sequence of computer experiments is used to study questions concerning the tsunami problem as a quantitative estimate of tsunami danger, detailed geographical tsunami classification, determination of the parameters of critical tsunami waves, and the conditions of their development. We call a wave critical, if its impact on the coast is most hazardous.Using the Middle Kuril Island as an example, we present the results of a computer experiment which includes determining the wavefields on the shelf and estimating the effects connected with the deep-water Bussol and Diana Straits.Numerical simulation of tsunami waves of different sources permits the assessment of the extent of tsunami danger in different areas of the coastal zone of Simushir Island, depending on the location of the focus zone and their geometry.The major singularities of the wavefield arise in the zones of the deep-water straits. The distribution of the amplification factors is determined by both the global parameters of the wavefields and the local properties of individual harbours. The results obtained for a particular harbour in the northern part of Simushir Island, formed the basis for the quantitative estimate of tsunami danger for this area.     

A chlorine isotope study of DSDP/ODP serpentinized ultramafic rocks: Insights into the serpentinization process

Eight DSDP/ODP cores were analyzed for major ion concentrations and δ³⁷Cl values of water-soluble chloride (δ³⁷Cl_WSC) and structurally bound chloride (δ³⁷Cl_SBC) in serpentinized ultramafic rocks. This diverse set of cores spans a wide range in age, temperature of serpentinization, tectonic setting, and geographic location of drilled serpentinized oceanic crust. Three of the cores were sampled at closely spaced intervals to investigate downhole variation in Cl concentration and chlorine isotope composition.The average total Cl content of all 86 samples is 0.26 ± 0.16 wt.% (0.19 ± 0.10 wt.% as water-soluble Cl (X_WSC) and 0.09 ± 0.09 wt.% as structurally bound Cl (X_SBC)). Structurally bound Cl concentration nearly doubles with depth in all cores; there is no consistent trend in water-soluble Cl content among the cores. Chlorine isotope fractionation between the structurally bound Cl⁻ site and the water-soluble Cl⁻ site varies from − 1.08‰ to + 1.16‰, averaging to + 0.21‰. Samples with negative fractionations may be related to reequilibration of the water-soluble chloride with seawater post-serpentinite formation. Six of the cores have positive bulk δ³⁷Cl values (+ 0.05‰ to + 0.36‰); the other two cores (173-1068A (Leg-Hole) and 84-570) have negative bulk δ³⁷Cl values (− 1.26‰ and − 0.54‰). The cores with negative δ³⁷Cl values also have variable Cl⁻ / SO₄² ratios, in contrast to all other cores. The isotopically positive cores (153-920D and 147-895E) show no isotopic variation with depth; the isotopically negative core (173-1068A) decreases by ∼1‰ with depth for both the water-soluble and structurally bound Cl fractions.Non-zero bulk δ³⁷Cl values indicate Cl in serpentinites was incorporated during original hydration and is not an artifact of seawater infiltration during drilling. Cores with positive δ³⁷Cl values are most likely explained by open system fractionation during hydrothermal alteration, with preferential incorporation of ³⁷Cl from seawater into the serpentinite and loss of residual light Cl back to the ocean. Fluid / rock ratios were probably low as evidenced by the presence of water-soluble salts. The two isotopically negative cores are characterized by a thick overlying sedimentary package that was in place prior to serpentinization. We believe the low δ³⁷Cl values of these cores are a result of hydration of ultramafic rock by infiltrating aqueous pore fluids from the overlying sediments. The resulting serpentinites inherit the characteristic negative δ³⁷Cl values of the pore waters. Chlorine stable isotopes can be used to identify the source of the serpentinizing fluid and ultimately discern chemical and tectonic processes involved in serpentinization.     

Fluid sources for the La Guitarra epithermal deposit (Temascaltepec district, Mexico): Volatile and helium isotope analyses in fluid inclusions

The La Guitarra deposit (Temascaltepec district, South-Central Mexico), belongs to the low/intermediate sulfidation epithermal type, has a polymetallic character although it is currently being mined for Ag and Au. The mineralization shows a polyphasic character and formed through several stages and sub-stages (named I, IIA, IIB, IIC, IID, and III). The previous structural, mineralogical, fluid inclusion and stable isotope studies were used to constrain the selection of samples for volatile and helium isotope analyses portrayed in this study. The N₂/Ar overall range obtained from analytical runs on fluid inclusion volatiles, by means of Quadrupole Mass Spectrometry (QMS), is 0 to 2526, and it ranges 0 to 2526 for stage I, 0 to 1264 for stage IIA, 0 to 1369 for stage IIB, 11 to 2401 for stage IIC, 19 to 324 for stage IID, and 0 to 2526 for stage III. These values, combined with the CO₂/CH₄ ratios, and N₂-He-Ar and N₂-CH₄-Ar relationships, suggest the occurrence of fluids from magmatic, crustal, and shallow meteoric sources in the forming epithermal vein deposit. The helium isotope analyses, obtained by means of Noble Gas Mass Spectrometry, display R/Ra average values between 0.5 and 2, pointing to the occurrence of mantle-derived helium that was relatively diluted or “contaminated” by crustal helium. These volatile analyses, when correlated with the stable isotope data from previous works and He isotope data, show the same distribution of data concerning sources for mineralizing fluids, especially those corresponding to magmatic and crustal sources. Thus, the overall geochemical data from mineralizing fluids are revealed as intrinsically consistent when compared to each other.The three main sources for mineralizing fluids (magmatic, and both deep and shallow meteoric fluids) are accountable at any scale, from stages of mineralization down to specific mineral associations. The volatile and helium isotope data obtained in this paper suggest that the precious metal-bearing mineral associations formed after hydrothermal pulses of predominantly oxidized magmatic fluids, and thus it is likely that precious metals were carried by fluids with such origin. Minerals from base-metal sulfide associations record both crustal and magmatic sources for mineralizing fluids, thus suggesting that base metals could be derived from deep leaching of crustal rocks. At the La Guitarra epithermal deposit there is no evidence for an evolution of mineralizing fluids towards any dominant source. Rather than that, volatile analyses in fluid inclusions suggest that this deposit formed as a pulsing hydrothermal system where each pulse or set of pulses accounts for different compositions of mineralizing fluids.The positive correlation between the relative content of magmatic fluids (high N₂/Ar ratios) and H₂S suggests that the necessary sulfur to carry mostly gold as bisulfide complexes came essentially from magmatic sources. Chlorine necessary to carry silver and base metals was found to be abundant in inclusion fluids and although there is no evidence about its source, it is plausible that it may come from magmatic sources as well.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.     

Temperature–composition–time (T–X–t) data from authigenic K-feldspar: An integrated methodology for dating fluid flow events

Integration of fluid inclusion analysis with high spatial resolution Ar–Ar dating of K-feldspar cements has been used to resolve and reconstruct palaeo-fluid flow. Fluid inclusion analysis allows discrimination of distinct cement phases, thereby identifying discrete episodes of fluid flow. Ar–Ar dating of the same cements via high spatial resolution laserprobe establishes absolute age constraints on the framework previously constructed. Integration of these two datasets yields temperature–composition–time data.     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).