Separate lead isotope analyses of leachate and residue rock fractions: implications for metal source tracing in ore deposit studies

Separate lead isotope analyses of leachate and residue fractions are applied to a broad spectrum of rocks commonly investigated in metallogenic studies. Resulting data highlight a systematic behavior of leachate and residue fractions with respect to lead isotope compositions, which essentially depends on the mineralogical composition of the rock. Granitoid and high-grade metamorphic rocks have residue compositions virtually identical to common lead. In contrast, low-grade metasedimentary rocks may have residue compositions swamped by radiogenic lead of leach-resistant zircons. Mafic magmatic rocks have residues that are often more radiogenic than leachates, depending on the ratio of leach-refractory zircons to common lead in the residual fraction of these rocks. Separate leachate and residue analyses of source rocks provide two lead isotope end members whose mixture may represent lead with the appropriate ore fluid composition. Our leaching experiments indicate that hot acid solutions (and by inference hydrothermal fluids) preferentially leach radiogenic lead from medium- to high-grade metamorphic and granitoid rocks, whereas they preferentially leach common lead from low-grade metasedimentary and mafic magmatic rocks. The method presented in this study provides a reliable alternative to other methods (i.e., age-correction of bulk-rock compositions) for the determination of the common lead signature of felsic to intermediate magmatic rocks. This may be preferable to age-corrected bulk-rock analyses, where ages to apply for corrections of bulk-rock data are not known or where moderately to highly altered rocks must be used. Case studies of orogenic gold and MVT districts of Peru (Pataz and San Vicente, respectively) show that separate leachate and residue lead isotope analyses carried out systematically on whole rocks allow a more thorough evaluation of metal source reservoirs than does the standard method of age-corrected or uncorrected bulk-rock analyses.     

Fault plane solutions of the january 26th, 2001 bhuj earthquake sequence

A 12-station temporary microearthquake network was established by the Geological Survey of India for aftershock monitoring of the January 26th, 2001 Bhuj earthquake (M _w 7.6) in the Kutch district of Gujarat state, western India. The epicentres of the aftershocks show two major trends: one in the NE direction and the other in the NW direction. Fault-plane solutions of the best-located and selected cluster of events that occurred along the NE trend, at a depth of 15–38 km, show reverse faulting with a large left-lateral strike-slip motion, which are comparable with the main-shock solution. The NW trending upper crustal aftershocks at depth <10 km, on the other hand, show reverse faulting with right-lateral strike-slip motion, and the mid crustal and lower crustal aftershocks, at a depth of 15–38 km, show pure reverse faulting as well as reverse faulting with right-lateral and left-lateral strike-slip motions; these solutions are not comparable with the main-shock solution. It is inferred that the intersection of two faults has been the source area for stress concentration to generate the main shock and the aftershocks.     

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

Nest microhabitat selection in larks from a European semi-arid shrub-steppe: the role of sunlight and predation

We analysed nest site selection in two shrub-steppe larks, the lesser short-toed larkCalandrella rufescensand the Thekla larkGalerida theklae. Both species nest on the ground next to a bush preferably facing north or north-east. Nestling growth and fledging weight are greater when the nest faces these directions. Both species prefer to nest next to camephytes with little foliage, and the lesser short-toed lark chooses sites with high visibility. We conclude that sunlight and nest predation are fundamental selective pressures determining nest microhabitat selection in larks in semi-arid areas.     

Synthetic Seismograms for a Layered Earth Geological Model Using the Absorption and Dispersion Phenomena

—An algorithm has been developed to compute the dispersive and dissipative seismic response using FUTTERMAN’S (1962) third attenuation-dispersion relationship. In the computation, frequency-dependent velocity and quality factor Q have been used but in the case of the nondispersive synthetic seismogram, frequency-independent velocity has been used. The model’s parameters are density, phase velocity, quality factors and thicknesses of the layers. Dispersive and nondispersive synthetic seismograms have been computed with and without absorption for a layered earth geological model. Fast Fourier transform (FFT) technique has been adopted for converting the frequency domain response into the time domain. The frequency spacing, Δf = 0.976?Hz, has been considered to avoid the aliasing effect. The results have revealed changes in the reflected waveforms in the frequency domain as well as in the time domain for absorption and dispersion cases. It is also concluded that dispersion reduces the arrival time and this effect is increasing with the travel time. The effect of constant Q on the seismic response has also been studied.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.     

Remote sensing detection of atmospheric pollutants by differential absorption LIDAR 510M/SODAR PA2 mobile system

Summary Air pollution problem in the Czech Republic is very complex due to large number of factors as turbulent dispersion and chemical reactions concurring in the status of the PBL. For this 3D distribution of pollutants is a key information that cannot be given by conventional monitoring stations. Combined use of DIAL lidar and sodar can give access to such information. Case studies of air pollutants transport are presented using 3D concentration fields of NO₂ and ozone measured by lidar and meteorological conditions monitored by Doppler Sodar.     

Optimal design of viscoelastic dampers using eigenvalue assignment

In this study a procedure for determining the optimum size and location of viscoelastic dampers is proposed using the eigenvalue assignment technique. Natural frequencies and modal damping ratios, required to realize a given target response, are determined first by the convex model. Then the desired dynamic structural properties are realized by optimally distributing the damping and stiffness coefficients of viscoelastic dampers using non‐linear programming based on the gradient of eigenvalues. This optimization method provides information on the optimal location as well as the magnitude of the damper parameters. The proposed procedure is applied to the retrofit of a 10‐story shear frame, and to a three‐dimensional structure with an asymmetric plan. The analysis results confirm that the responses of model structures retrofitted by the proposed method correspond well with the given target response. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic structural control using an electric servomotor active mass driver system

This study investigates an electric‐type active mass driver (AMD) system for structural vibration control. Composed primarily of an electric servomotor and a ball screw, the electrical AMD system is free from noise problems, oil leakage, and labor‐intensive maintenance that commonly are associated with hydraulic AMD systems. The desired stroke amplification of the mass and the power demand of the servomotor can be adjusted via the ball screw pitch, which in turn affects the effectiveness and efficiency of the system. Meanwhile, an instantaneous optimal direct output feedback control algorithm is adopted. Numerical simulation is performed using a five‐story steel frame as the object structure under the conditions of the 1940 El Centro earthquake. The AMD system proves to be effective and efficient within a certain range of the ball screw pitch. The reductions of the peak responses can reach as high as 70% if properly designed. Requiring only the velocity measurement of the top floor for on‐line feedback control, the proposed control algorithm is recommended for practical implementation. Copyright © 2004 John Wiley & Sons, Ltd.