Geochemical and structural evolution of micas in the Rožná and Dobrá Voda pegmatites,Czech Republic

Summary The chemistry, structural parameters, polytypism, optical properties and Rb-Sr isotopes were examined in 11 to 60 samples of biotite, muscovite and lepidolite from the pegmatites at Roná (the type locality of lepidolite; 323 ± 4Ma) and Dobrá Voda (306 ± 9Ma) in western Moravia. At both localities, early endocontact biotite is followed inwards by muscovite and lepidolite, which is concentrated in and around the core. At Roná, a 1M lepidolite follows after 2M₁ muscovite but all later generations of lepidolite are 2M₂, close to Tri₅₀ Ply₅₀ and in part associated with muscovite 2M₁. At Dobrá Voda, all lepidolite types are 1M and free of muscovite, and the late varieties approximate Tri₃₀ Ply₇₀. At both localities, a trend of increasing HF is indicated during the progress of mica crystallization, culminating in precipitation of topaz. Polytypism of lepidolite is not correlatable with any compositional or growth feature, or their combination. Throughout the mica crystallization, Rb/Cs decreases but K/Rb becomes reversed after an initial decrease. Boron is partitioned preferentially into muscovite (up to 1.10 wt.% B₂O₃) but Be, Zn, Mn and Sc are enhanced in lepidolite. A slight increase in Fe, Ba and Cl in the last generation of lepidolite might be possibly due to mixing of residual pegmatite fluids with metamorphic pore solutions.

---

Geochemische und strukturelle Entwicklung der Glimmer in den Pegmatiten von Roná und Dobrá Voda, Tschechische Republik

Zusammenfassung In 11 bis 60 Proben von Biotit, Muskovit und Lepidolith aus den Pegmatiten von Roná (Typlokalität des Lepidoliths; 323 ± 9 Ma) in Westmähren wurden Chemie, Struktur-parameter, Polytypie, optische Eigenschaften und Rb-Sr-Isotopie untersucht. An beiden Lokalitäten wird früher Biotit an Endokontakten nach Innen von Muskovit und Lepidolith gefolgt, letzterer ist in und um den Kern konzentriert. In Roná folgt 1M-Lepidolith auf 2M₁-Muskovit, aber alle späteren Lepidolithgenerationen sind 2M₂, nahe Tri₅₀Ply₅₀ und zum Teil mit 2M₁-Muskovit vergesellschaftet. In Dobrá Voda sind alle Lepidolithe vom Typ 1M und frei von Muskovit, die späten Varietäten kommen Tri₅₀Ply₅₀ nahe. An beiden Lokalitäten ist während des Fortschreitens der Glimmerkristallisation eine Tendenz von steigendem HF angezeigt, die in der Ausfällung von Topas ihren Höhepunkt findet. Die Polytypie des Lepidoliths kann nicht mit irgendeiner Eigenheit der Zusammensetzung oder des Wachstums korreliert worden, auch nicht mit einer Kombination von diesen. Während der ganzen Glimmerkristallisation nimmt Rb/Cs ab, aber die Tendenz von K/Rb ändert sich nach anfänglichem Abfall. Das Bor verteilt sich bevorzugt auf den Muskovit (bis zu 1.10 Gew. -% B₂O₃), aber die Be-, Zn-, Mn- und Sc-Gehalte sind im Lepidolith erhöht. Ein leichter Ansteig von Fe, Ba und Cl in der letzten Lepidolithgeneration könnte vielleicht durch eine Mischung von pegmatitischen Restlösungen mit metamorphen Porenlösungen verursacht sein.

---

---

With 9 Figures     

Incorporating hydro-mechanical coupling in an analysis of the effects of rainfall patterns on unsaturated soil slope stability

Rainfall-induced landslides can cause loss of life and damage to property, infrastructure, and the environment. Rainfall patterns affect the pore-water pressure of unsaturated soil slopes, and are related to the slopes’ stability. Four rainfall patterns were chosen to represent natural rainfall patterns for an examination of rainfall infiltration into soil slopes using numerical models incorporating coupled water infiltration and deformation in unsaturated soils. Our analysis showed that rainfall patterns play a significant role in the distribution of the pore-water pressure in soil slopes, and influence the slope stability. The pore-water pressure profile of soil slopes and the factor of safety are affected by the ratio of rainfall intensity and the coefficient of permeability. The depth and shape of the shallow sliding plane of the landslide is closely related to the rainfall pattern; moreover, the results showed a correlation between the factor of safety of the slope and the rainfall intensity. This relationship can be described by a dimensionless rainfall intensity. The nonlinear relationship can be used to estimate the slope stability resulting from rainfall infiltration when the hydro-mechanical coupling in unsaturated soil slopes is considered.     

Suitability analysis for sand and gravel extraction site location in the context of a sustainable development in the surroundings of Zaragoza (Spain)

Zaragoza city is located in the central Ebro Basin, in the Iberian Peninsula. The fluvial terraces formed by the Ebro River present a valuable resource of sand and gravel deposits. However, taking advantage of these available resources implies conflicts with other land use interests like urban and industrial development as well as agricultural use, which has also traditionally occupied the alluvial terraces. These deposits represent a substantial groundwater resource that should be preserved for future generations. The development of spatial decision support systems (SDSS) has greatly assisted efforts for solving land-use conflicts. These systems combine the benefits of geographic information systems (GIS) and decision support methodologies and are therefore suitable to manage sustainable development of urban areas. In this contribution, an extraction suitability map taking into consideration a variety of environmental criteria is created with the help of a SDSS. The method used is the analytical hierarchy process which is integrated in ArcGIS. Areas most suitable to sand and gravel extraction are located in the high terraces, and in those terraces covered by pediments where the thickness of resource is relatively high. These areas are far from valuable natural areas, outside areas most vulnerable to groundwater contamination, and beneath soils with poor irrigation characteristics.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

GPR studies over the tsunami affected Karaikal beach,Tamil Nadu,south India

In this study, results of GPR profiling related to mapping of subsurface sedimentary layers at tsunami affected Karaikal beach are presented . A 400 MHz antenna was used for profiling along 262 m stretch of transect from beach to backshore areas with penetration of about 2.0 m depth (50 ns two-way travel time). The velocity analysis was carried out to estimate the depth information along the GPR profile. Based on the significant changes in the reflection amplitude, three different zones are marked and the upper zone is noticed with less moisture compared to other two (saturated) zones. The water table is noticed to vary from 0.5 to 0.75 m depth (12–15 ns) as moving away from the coastline. Buried erosional surface is observed at 1.5 m depth (40–42 ns), which represents the limit up to which the extreme event acted upon. In other words, it is the depth to which the tsunami sediments have been piled up to about 1.5 m thickness. Three field test pits were made along the transect and sedimentary sequences were recorded. The sand layers, especially, heavy mineral layers, recorded in the test pits indicate a positive correlation with the amplitude and velocity changes in the GPR profile. Such interpretation seems to be difficult in the middle zone due to its water saturation condition. But it is fairly clear in the lower zone located just below the erosional surface where the strata is comparatively more compact. The inferences from the GPR profile thus provide a lucid insight to the subsurface sediment sequences of the tsunami sediments in the Karaikal beach.     

Estimation of Regional Evapotranspiration by Combining Aircraftand Ground-Based Measurements

Two methods are examined for combining measurements from instrumented aircraftand towers to estimate regional scale evapotranspiration. Aircraft data provided spatially averaged values of properties of the surface, the evaporative fraction and maximum stomatal conductance. These quantities are less sensitive to meteorological conditions than the turbulent fluxes of heat and water vapour themselves. The methods allowed aircraft data collected over several days to be averaged and thus to reduce the random error associated with the temporal under-sampling inherent in aircraft measurements. Evaporative fraction is estimated directly from the aircraft data, while maximum stomatal conductance is estimated by coupling the Penman–Monteith equation to a simple model relating surface conductance to the incoming shortwave radiation and specific humidity saturation deficit. The spatial averages of evaporative fraction and maximum stomatal conductance can then be used with routine tower data to estimate the regional scale evapotranspiration. Data from aircraft flights and six ground based sites during the OASIS field campaign in south–east New South Wales in 1995 have been used to check the methods. Both the evaporative fraction and the maximum stomatal conductance derived from the aircraft data give information on the spatial variability of the surface energy budget at scales from 10 to 100 km. Daily averaged latent heat fluxes estimated using these methods for the OASIS study region agree with the available observations in quasi-stationary conditions or in weakly non-stationary conditions when the data from several aircraft flights are averaged to reduce the impact of short term imbalances in the surface energy budget.     

Comparison between geodetic and oceanographic approaches to estimate mean dynamic topography for vertical datum unification: evaluation at Australian tide gauges

The direct method of vertical datum unification requires estimates of the ocean’s mean dynamic topography (MDT) at tide gauges, which can be sourced from either geodetic or oceanographic approaches. To assess the suitability of different types of MDT for this purpose, we evaluate 13 physics-based numerical ocean models and six MDTs computed from observed geodetic and/or ocean data at 32 tide gauges around the Australian coast. We focus on the viability of numerical ocean models for vertical datum unification, classifying the 13 ocean models used as either independent (do not contain assimilated geodetic data) or non-independent (do contain assimilated geodetic data). We find that the independent and non-independent ocean models deliver similar results. Maximum differences among ocean models and geodetic MDTs reach >150 mm at several Australian tide gauges and are considered anomalous at the 99% confidence level. These differences appear to be of geodetic origin, but without additional independent information, or formal error estimates for each model, some of these errors remain inseparable. Our results imply that some ocean models have standard deviations of differences with other MDTs (using geodetic and/or ocean observations) at Australian tide gauges, and with levelling between some Australian tide gauges, of \({\sim }\pm 50\,\hbox {mm}\). This indicates that they should be considered as an alternative to geodetic MDTs for the direct unification of vertical datums. They can also be used as diagnostics for errors in geodetic MDT in coastal zones, but the inseparability problem remains, where the error cannot be discriminated between the geoid model or altimeter-derived mean sea surface.     

Geochemistry of the Adamello massif (northern Italy)

The Tertiary Adamello massif, outcropping over an area of more than 550 km² in the southern Alps (northern Italy) is composed mainly of granitoid rocks (granodiorite, tonalite, quartz diorite) with minor amounts of diorite and gabbro. The major and trace element composition of these rocks is comparable to calc-alkaline volcanic rocks of continental margins. The granitoid rocks display spatial and temporal variations in their composition, particularly in Na, P, Sr, La, Nb and Y contents and ⁸⁷Sr/⁸⁶Sr ratios. The variations were probably produced by concurrent contamination/wall-rock assimilation and fractional crystallization of high-alumina basaltic magma.     

Metal sources of the Navan carbonate-hosted base metal deposit, Ireland: Nd and Sr isotope evidence for deep hydrothermal convection

Nd and Sr isotope analyses are presented for gangue mineral samples from the giant carbonate-hosted Navan Zn–Pb deposit, Ireland, and for rocks from which Navan metals may have been derived. Analysis of gangue minerals spanning the Navan paragenetic sequence reveals systematic evolution in the composition of the mineralising fluid. Early fluid represented by replacive dolomite exhibits the lowest initial ⁸⁷Sr/⁸⁶Sr ratio (0.7083–0.7086), closest to that of the host limestone and to Lower Carboniferous seawater, and the highest ¹⁴³Nd/¹⁴⁴Nd ratio (0.51161–0.51176). Later generations of dolomite, barite and calcite, which encompass sulphide precipitation, have higher initial ⁸⁷Sr/⁸⁶Sr ratios (maximum 0.7105) and lower initial ¹⁴³Nd/¹⁴⁴Nd ratios (minimum 0.51157). All samples have initial Nd isotope ratios that are too low to have been acquired only from the host limestone. Drill core samples of presumed Ordovician volcanic and sedimentary rocks from beneath the Navan orebody have ¹⁴³Nd/¹⁴⁴Nd and ⁸⁷Sr/⁸⁶Sr ratios at the time of mineralisation of 0.51184–0.51217 and 0.7086–0.7138, respectively. The data are interpreted to indicate mixing of sulphide-rich, limestone-buffered brine, with a metal-bearing hydrothermal fluid, which had passed through sub-Carboniferous rocks, consistent with published fluid inclusion and S isotope data. The ¹⁴³Nd/¹⁴⁴Nd ratio of this basement-derived fluid is too low to have been imparted by flow through the Devonian Old Red Sandstone, as required in models of regional fluid flow in response to Hercynian uplift. Irrespective of whether such regional fluid flow occurred, the hydrothermal Nd must have been derived from sub-Devonian rocks. These conclusions broadly support the hydrothermal convection cell model in which brines, ultimately of surface origin, penetrated to a depth of several kilometres, leaching metals from the rocks through which they passed. The data also support increasing depth of penetration of convection cells with time. Metals were subsequently precipitated in carbonate rocks at sites of mixing with cooler, sulphide-rich fluids. However, comparison of the Navan hydrothermal gangue Nd–Sr isotope data with data from Lower Palaeozoic rocks strongly suggests that the latter cannot alone account for the “basement” signature. As the Navan deposit lies immediately north of the Iapetus Suture, this suggests that the Laurentian margin includes Precambrian basement.     

On the assessment of maximum earthquakes in the alps and adjacent areas

Strong tectonic earthquakes within the crust always occur on already existing faults, and they have the property of a shear rupture. Such earthquakes with surface-wave magnitudes M < 7 obviously have a geometric similarity. Because of this similarity and the validity of the Gutenberg and Richter's energy—magnitude relation, the expression M = 2 log₁₀ L + const., with L = focal length, is valid.The expression L_maxL^* for the maximum focal length, is also valid if L^* is the length of the rectilinear extent of the seismic line on which the maximum earthquake occurs. The bounds of L^* may be given by sharp bends and/or by traversing deep faults. Thus the maximum imaginable earthquake on a seismic line with the length L^* has the magnitude M_max = 2 log₁₀ L^* + const.For the investigated region — the Alps and adjacent areas — from the data of recent and historical strong earthquakes, it follows that M_max = 2 log₁₀ L^* + 1.7, if L^* is measured in kilometres. These limiting values lie in the centre-field of the magnitude range for maximum earthquakes, published by Shebalin in 1970. By the aid of this equation it is also possible to assess the upper limiting value of the accompanying maximum scale intensity.