The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.     

Metamorphic kyanite eclogites in the lufilian arc of Zambia

Eclogites formed by the Lufilian (post-Katanga) metamorphism of gabbros intruded into rocks of the Katanga System (Upper Proterozoic) occur in the Central Province of Zambia. Typical constituents of these rocks are omphacite, almandine garnet with significant contents of the pyrope and grossular components, kyanite and rutile. Eclogites from some localities display well preserved, relict ophitic texture inherited from the original gabbro. Hornblende-scapolite metagabbro and coronite metagabbro often accompany the eclogites, and metalherzolite and serpentinite occur at a few localities. The Zambian eclogites are broadly associated with a variety of relatively high-pressure, kyanite-bearing mineral assemblages some of which are possibly in the eclogite facies. It is suggested that, in addition to temperature and load pressure, fluid pressure, fluid composition (H₂O, CO₂, Cl, SO₃) and permeability of the gabbro due to deformation were important factors for the transformation of gabbro to eclogite in the environment of the Lufilian Arc.     

Response of single benthic metrics and multi-metric methods to anthropogenic pressure gradients, in five distinct European coastal and transitional ecosystems

In recent times many benthic indices have been proposed to assess the ecological quality of marine waters worldwide. In this study we compared single metrics and multi-metric methods to assess coastal and transitional benthic status along human pressure gradients in five distinct environments across Europe: Varna bay and lake (Bulgaria), Lesina lagoon (Italy), Mondego estuary (Portugal), Basque coast (Spain) and Oslofjord (Norway). Hence, 13 single metrics (abundance, number of taxa, and several diversity and sensitivity indices) and eight of the most common indices used within the European Water Framework Directive (WFD) for benthic assessment were selected: index of size spectra (ISS), Benthic assessment tool (BAT), Norwegian quality index (NQI), Multivariate AMBI (M-AMBI), Benthic quality index (BQI), (Benthic ecosystem quality index (BEQI), Benthic index based on taxonomic sufficiency (BITS), and infaunal quality index (IQI). Within each system, sampling sites were ordered in an increasing pressure gradient according to a preliminary classification based on professional judgement. The different indices are largely consistent in their response to pressure gradient, except in some particular cases (i.e. BITS, in all cases, or ISS when a low number of individuals is present). Inconsistencies between indicator responses were most pronounced in transitional waters (i.e. IQI, BEQI), highlighting the difficulties of the generic application of indicators to all marine, estuarine and lagoonal environments. However, some of the single (i.e. ecological groups approach, diversity, richness) and multi-metric methods (i.e. BAT, M-AMBI, NQI) were able to detect such gradients both in transitional and coastal environments, being these multi-metric methods more consistent in the detection than single indices. This study highlights the importance of survey design and good reference conditions for some indicators. The agreement observed between different methodologies and their ability to detect quality trends across distinct environments constitutes a promising result for the implementation of the WFD’s monitoring plans. Moreover, these results have management implications, regarding the dangers of misclassification, uncertainty in the assessment, use of conflicting indices, and testing and validation of indices.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

North Atlantic atmospheric circulation and surface wind in the Northeast of the Iberian Peninsula: uncertainty and long term downscaled variability

The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.     

Objective method for classifying air masses: an application to the analysis of Buenos Aires’ (Argentina) urban heat island intensity

Summary ?During recent years, numerous studies have examined the Buenos Aires urban climate, but the relationship between large-scale weather conditions and the Buenos Aires urban heat island (UHI) intensity has not been studied. The goal of this paper is to apply an objective synoptic climatological method to identify homogeneous air masses or weather types affecting Buenos Aires during winter, and to relate the results to the UHI intensity. A K-means clustering method was used to define six different air masses considering the 03:00, 09:00, 15:00 and 21:00 LT surface observations of dry bulb temperature, dew point, cloud cover, atmospheric pressure and wind direction and velocity at Ezeiza, the most rural meteorological station of the Buenos Aires metropolitan area (Fig. 1). Results show that the mean UHI intensity is at its maximum (2.8 °C) a few hours before sunrise when conditions are dominated by cold air masses associated with cold-core anticyclones, weak winds and low cloud cover. Inverse heat islands are found during the afternoon for all air masses indicating that surface processes are not dominant at that time. The relatively infrequent and warmest air mass is the only one that presents a mean negative urban-rural temperature difference (−0.1 °C) during the afternoon with the smallest diurnal cycle of the UHI intensity probably due to the prevailing high humidity and cloudy sky conditions. The paper provides an insight into the Buenos Aires urban–rural temperature difference under a variety of winter weather types and results could be useful to improve local daily temperature forecasts for the metropolitan area of Buenos Aires on the basis of the routine forecasts of weather types. Received October 24, 2001; revised June 12, 2002; accepted October 10, 2002     

Surface Synoptic Circulation and Daily Precipitation in Catalonia

Summary Catalonian daily precipitation patterns were obtained by means of S-mode Principal Component Analysis (PCA) applied to subsets built up according to synoptic surface airflow. Daily data from November to April 1990–1994 from 45 raingauge sites in Catalonia without missing data were used. Eight circulation classes for days in this period were considered. A Precipitation Activity Index (PAI) was calculated for each station and class, leading to the determination of the highest activity regions for each subset. In addition, we determined which circulation type produced most precipitation at each station of the network. Separate PCAs were performed for each of the eight classes considered and the first four Varimax rotated solutions were interpreted for each circulation type. The physical meaning of the rest of the factors were omitted as they were associated with local behaviour. Finally, some concluding remarks on the cyclogenic properties of the Western Mediterranean Basin and the diversity of the extracted patterns are presented and supported by the authors experience in forecasting in Catalonia. Received February, 11, 1997 Revised May 23, 1997     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

Estimating Buoyancy Heat Flux Using the Surface Renewal Technique over Four Amazonian Forest Sites in Brazil

The surface renewal (SR) method was applied for the first time to measurements of air temperature over four Amazonian forest sites and different seasons in order to obtain estimates of buoyancy heat flux. The required calibration of this method against eddy covariance resulted in a value for a specific parameter that is close to the range reported in other studies, contributing to the generalization of the SR method to different kinds of canopies. The comparison with fluxes obtained using the eddy covariance technique revealed a good match between the two methods for different sites, heights and seasons. Sites with high levels of non-stationarity in the signals of temperature and wind speed presented higher scatter in the regression with fluxes from eddy covariance. For a particular site with previously reported influence of low-frequency motions, the regression was only satisfactory, i.e., slope parameter close to unity and small offset, when oscillations with periods longer than $\approx $ ≈ 13 min were filtered out. The SR method has a great potential due to the simplicity of the instrumentation required. However, care should be taken when measuring under the influence of mesoscale motions, which can lead to high levels of non-stationarity, compromising the fundamental concepts of the SR theory.