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.     

Use of Fuzzy Membership Input Layers to Combine Subjective Geological Knowledge and Empirical Data in a Neural Network Method for Mineral-Potential Mapping

Use of GIS layers, in which the cell values represent fuzzy membership variables, is an effective method of combining subjective geological knowledge with empirical data in a neural network approach to mineral-prospectivity mapping. In this study, multilayer perceptron (MLP), neural networks are used to combine up to 17 regional exploration variables to predict the potential for orogenic gold deposits in the form of prospectivity maps in the Archean Kalgoorlie Terrane of Western Australia. Two types of fuzzy membership layers are used. In the first type of layer, the statistical relationships between known gold deposits and variables in the GIS thematic layer are used to determine fuzzy membership values. For example, GIS layers depicting solid geology and rock-type combinations of categorical data at the nearest lithological boundary for each cell are converted to fuzzy membership layers representing favorable lithologies and favorable lithological boundaries, respectively. This type of fuzzy-membership input is a useful alternative to the 1-of-N coding used for categorical inputs, particularly if there are a large number of classes. Rheological contrast at lithological boundaries is modeled using a second type of fuzzy membership layer, in which the assignment of fuzzy membership value, although based on geological field data, is subjective. The methods used here could be applied to a large range of subjective data (e.g., favorability of tectonic environment, host stratigraphy, or reactivation along major faults) currently used in regional exploration programs, but which normally would not be included as inputs in an empirical neural network approach.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.     

Extracting flow structures from tracer data

A method of visualizing structures in closed chaotic flows out of homogenous particle distributions is presented in the example of models of a meandering jet. To this end, the system will be leaked or opened up by defining a region of the flow, so that a particle is considered to be escaped if it leaves this region. By applying this method to an ensemble of nonescaped tracers, we are able to characterize mixing processes by visualizing the converging and stretching filamentations (stable and unstable manifolds) in the flow without using additional mathematical tools. The possibility of applying the algorithm to analyze buoy data, and a comparison with the finite time manifolds are discussed.     

Lower crustal granulite xenoliths from the Pannonian Basin,Hungary, Part 2: Sr–Nd–Pb–Hf and O isotope evidence for formation of continental lower crust by tectonic emplacement of oceanic crust

Mafic granulite xenoliths from the lower crust of the Pannonian Basin are dominated by LREE-depleted bulk-rock compositions. Many of these have MORB-like ¹⁴³Nd/¹⁴⁴Nd but ⁸⁷Sr/⁸⁶Sr is elevated relative to most MORBs. Their '¹⁸O values cover a wide range from +3.8 to +9.5‰. A group of LREE-enriched mafic granulites have higher ⁸⁷Sr/⁸⁶Sr (0.704-0.708) and lower ¹⁴³Nd/¹⁴⁴Nd (0.5128-0.5124), with higher '¹⁸O values on average (+7.8 to +10.6‰) than the LREE-depleted granulites. The LREE-enriched granulites are, however, isotopically similar to newly discovered metasedimentary granulite xenoliths. A sublinear correlation in )Hf-)Nd isotope space has a shallower slope than the crust-mantle array, with the metasedimentary rocks forming the low )Hf end member; the radiogenic end is restricted to the LREE-depleted granulites and these overlap the field of MORB. Pb isotopes for the LREE-depleted samples are less radiogenic on average than those of the LREE-enriched and metasedimentary xenoliths, and metasedimentary granulites have consistently higher ²⁰⁸Pb/²⁰⁴Pb. The wide range in '¹⁸O over a restricted range in Nd and Sr isotope values, in combination with the predominance of LREE-depleted trace-element compositions, is consistent with an origin as a package of hydrothermally altered oceanic crust. The existence of '¹⁸O values lower than average MORB and/or mantle peridotite requires that at least some of these rocks were hydrothermally altered at high temperature, presumably in the oceanic lower crust. The low ¹⁴³Nd/¹⁴⁴Nd of the LREE-enriched mafic granulites cannot be explained by simple mixing between a LREE-depleted melt and an enriched component, represented by the recovered metasediments. Instead, we interpret these rocks as the metamorphic equivalent of the shallowest levels of the ocean crust where pillow basalts are intimately intercalated with oceanic sediments. A possible model is accretion of oceanic crustal slices during subduction and convergence, followed by high-grade metamorphism during the Alpine orogeny.     

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.     

Surface-water quality for the Salí River watershed in NW Argentina

Spatial and temporal variations of Na, Mg, Si, K, Ca, SO₄^2-, Cl^-, HCO₃^-, Ag, Al, As, Au, B, Ba, Be, Bi, Br, Ce, Cd, Co, Cr, Cs, Cu, Dy, Er, Eu, F, Fe, Ga, Ge, Gd, Hf, Hg, Ho, I, La, Li, Mn, Mo, Nb, Nd, Ni, P, Pb, Pr, Rb, Sb, Sm, Sn, Sr, Ta, Tb, Te, Th, Ti, Tl, U, V, W, Y, Yb, Zn, and Zr were monitored through 37 sampling stations to determine the main aspects influencing the surface-water quality in the Salí River watershed (Tucumán Province, NW Argentina). The influence of the regional geological setting on water chemistry allows to distinguish three sub-basins. The interaction with sedimentary rocks was found to be dominant in the northern and central-eastern sub-basin as well as in the southern sub-basin, whereas the metamorphic-granitic basement of the Pampean Ranges was noted in the central sub-basin. In addition, anthropogenic activities affect the spatial variation of K, P, Mn, Rb, and Pb as well as dissolved oxygen concentrations and Eh. Temporal water-quality variation is related to the spatial distribution of precipitation and to the seasonal character of the main local industries (sugar cane, alcohol, citrus), increasing P and K concentrations and decreasing dissolved oxygen concentration and Eh in winter. Cl^-, Na, SO₄^2-, Al, As, B, Fe, Mn, Se and U concentrations exceed the regulated drinking-water thresholds at several sampling stations.     

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.     

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.