Mixing, sediment accumulation and focusing using 210Pb and 137Cs

The Portil lagoon is a natural freshwater reservoir located at the southwest of Spain, near the coast. In its surroundings an important tourist complex has been developed since the earlier 1970s. This has resulted in increased loads of nutrients, specially during summer months. In order to evaluate the impact from humans on the lagoon caused by anthropogenic activities, we have determined vertical profiles of excess ²¹⁰Pb and ¹³⁷Cs in a sediment core taken from its bottom. Vertical profile of excess ²¹⁰Pb seems to indicate that the sediment core was affected by large-scale mixing processes. Nevertheless this possibility was ruled out based on the shape of the ¹³⁷Cs vertical distribution after modelling the ¹³⁷Cs profile considering large-scale mixing. The chronology developed from ²¹⁰Pb and ¹³⁷Cs has allowed us to calculate the temporal evolution of the mass sedimentation rate during the last century. This mass sedimentation rate increased in 1973: from 0.08(2) g cm^-2 a^-1 in the period 1900-1973 to 0.17(5) g cm^-2 a^-1 in the period 1973-1995. This increase could be related to land-movements in the drainage area and to permitted sewage inputs into the lagoon, both with its origin in the growing/running of the surrounding tourist complex.     

An experimental study of nominally symmetric and asymmetric structures isolated with the FPS

This experimental investigation deals with the earthquake behaviour of a nominally symmetric and a mass‐asymmetric three‐storey structural model isolated with the frictional pendulum system (FPS). Both accidental and natural torsion are evaluated in the structure by using recorded accelerations in all building floors and measured deformations at the isolation level. A 3D‐shaking table was used to subject the model to five different ground motions, including impulsive as well as far‐field subduction‐zone type earthquakes. Results show that the analytical predictions of the earthquake behaviour of the isolated structure, as obtained from a physical model of the FPS, are in close agreement with the true complex inelastic measured behaviour of the FPS. Besides, experimental results also validate previous observations about the importance of accounting for the variability of the normal loads in modelling the earthquake behaviour of FPS isolators. Measured torsional deformation amplifications at the base of the building vary, in the mean, from 2.5% to 6% for the symmetric and asymmetric structural configurations, respectively. In relation to the fixed base structure, the reduction factors for the base shear of the isolated structure are, in the mean, about 3.9 for both configurations. Finally, it is concluded that the FPS is capable of controlling the lateral–torsional motions of mass‐asymmetric structures quite effectively by aligning the centre of mass of the superstructure with the centre of pendular and frictional resistance of the isolation system. Copyright © 2003 John Wiley & Sons, Ltd.     

Physical model for dynamic analysis of structures with FPS isolators

This article presents a physical model for frictional pendulum isolators (FPS) that is ready to be implemented in most commercial software. The model is capable of accounting for effects such as large deformations, sticking, and uplift and impact by sensing the normal loads in the isolators through a gap element. Sticking has been incorporated into the model by extending the Park–Wen hysteretic model to the case of large deformations. The proposed model has been tested against a theoretically ‘exact’ formulation leading to essentially identical results. To facilitate its use, the physical FPS model has been cast into a typical non‐linear structural element format, i.e. with deformation as input and restoring force as output. Examples of a building and a bridge have been chosen to show the potential of the element and to provide further insight into the earthquake response of structures with FPS isolators; in particular, in aspects such as the orientation in placement of the isolator, sticking, P? Δ, and other large deformation effects. Copyright © 2003 John Wiley & Sons, Ltd.     

Pre-Caving, Drilling and Blasting in the Esmeralda Sector of the El Teniente Mine

El Teniente Mine has developed the method of panel caving as an innovative pre-conditioning caving alternative for traditional mining of sectors located in primary rock. In applying this method, caving speed and area incorporation are the most relevant parameters to establish a desired production level in a specified sector. This is the context in which continuous improvements to several design alternatives have been developed and applied to optimize mine caving. Mine caving, by drilling and blasting, results in a narrow base cut of ore columns not larger than the height of the caving galleries, which usually measure 3.6 m. These confined blasting operations require drilling parallel blastholes that extend from the floor to the roof of the caving gallery. Panel caving aims at mitigating subvertical field stress transmission to the lower levels of the caving gallery. This improves the stability and helps maintain the advancement of the caved zone at the caving level, with respect to the development, preparation and ore mining activities in the production level. This paper analyzes and evaluates different panel caving alternatives. The result is the identification of the most adequate panel caving alternative for conditions in the Esmeralda Mine given safety, technical, operational and economic considerations.     

Preliminary Rock-Slope-Susceptibility Assessment Using GIS and the SMR Classification

This paper aims to examine the impact of large-scale structuraladjustments (like the Greater Dhaka Flood Protection Project, GDFPP) on local living environment.It focuses the importance of environmental factors in flood hazard mitigation, and examines theenvironmental attitudes of the floodplain residents arising from the large-scale structural adjustments.Based on `perceived natural hazard research perspectives', this paper examines: (i) the reasons for persistentfloodplain occupation, and (ii) the importance of environmental factors in the choice, motivations and decision-makingof floodplain residents.This research used data collected from 300 households situated inthe eastern part of Dhaka. The face-to-face household survey data provided individuals' responses to a structuredquestionnaire on hazards and environment. Survey concerned urban floodplains, and looked fordata on housing, household characteristics, and residents' attitudes. Results of interview surveys wereused to: (i) explore the reasons of floodplain occupation, and (ii) residents' attitudes to tolerable levelof flood risk and willingness to accept environmental change resulting from the proposed structural embankments inthe eastern perimeter of Dhaka City, Bangladesh.Findings revealed that floodplain occupation (by theindividuals' decision-making) was a result of overall reaction to the Government's structuraladjustment policies that resulted from institutional, locational and socio-economic factors. The attitude survey results provided residents' perception to hazards and environment to be dependenton the socio-economic factors – but in a complex manner, many factors are interrelated.In addition to support for structural embankments, the study sample displayed a common concernand widespread environmental awareness. In terms of any `trade-off' between thebenefits (resources) from the embankments and costs (hazards) due to the detrimental impact on environment, the residents of Dhaka, despite some concern forsacrificing embankments for environment, tended to show a generalconsensus for embankments.     

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.     

Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

Architecture of fault zones determined from outcrop, cores, 3-D seismic tomography and geostatistical modeling: example from the Albalá Granitic Pluton, SW Iberian Variscan Massif

The 3-D seismic tomographic data are used together with field, core and well log structural information to determine the detailed 3-D architecture of fault zones in a granitic massif of volume 500×575×168 m at Mina Ratones area in the Albalá Granitic Pluton. To facilitate the integration of the different data, geostatistical simulation algorithms are applied to interpolate the relatively sparse structural (hard) control data conditioned to abundant but indirect 3-D (soft) seismic tomographic data. To effectively integrate geologic and tomographic data, 3-D migration of the velocity model from the time domain into the depth domain was essential. The resulting 3-D model constitutes an image of the fault zone architecture within the granitic massif that honours hard and soft data and provides an evaluation of the spatial variability of structural heterogeneities based on the computation of 3-D experimental variograms of Fracture Index (fault intensity) data. This probabilistic quantitative 3-D model of spatially heterogeneous fault zones is suitable for subsequent fluid flow simulations. The modeled image of the 3-D fault distribution is consistent with the fault architecture in the Mina Ratones area, which basically consists of two families of subvertical structures with NNE–SSW and ENE–WSW trends that displaces the surfaces of low-angle faults (North Fault) and follows their seismically detected staircase geometry. These brittle structures cut two subvertical dykes (27 and 27′ Dykes) with a NNE–SSW to N–S trend. The faults present high FI (FI>12) adjacent bands of irregular geometry in detail that intersect in space delimiting rhombohedral blocks of relatively less fractured granite (FI<6). Both structural domains likely correspond with the protolith and the damaged zone/fault core in the widely accepted model for fault zone architecture. Therefore, the construction of 3-D grids of the FI in granitic areas affected by brittle tectonics permits the quantitative structural characterization of the rock massif.     

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.