Selenium cycling in volcanic environments: the role of soils as reactive interfaces

Selenium (Se) is an element with important health implications that is emitted in significant amounts from volcanoes. Attracted by the fertility of volcanic soils, around 10% of the world population lives within 100 km of an active volcano. Nevertheless, the behaviour of Se in volcanic environments is poorly understood. Therefore, the main aim of this thesis is to investigate the role of soils in the Se cycling in volcanic environments. Prior to the geochemical studies, precise and accurate methods for the determination of Se contents, speciation and isotopic signatures were developed. Afterwards, a combination of field studies and lab controlled experiments were performed with soils from two contrasting European volcanic settings: Mount Etna in Sicily (Italy) and Mount Teide in Tenerife (Spain). The results showed a strong link between Se behaviour and soil development, indicating that Se mobility in volcanic soils is controlled by sorption processes and soil mineralogy.

A review of the theoretical basis for bulk mass flux convective parameterization

Most parameterizations for precipitating convection in use today are bulk schemes, in which an ensemble of cumulus elements with different properties is modelled as a single, representative entraining-detraining plume. We review the underpinning mathematical model for such parameterizations, in particular by comparing it with spectral models in which elements are not combined into the representative plume. The chief merit of a bulk model is that the representative plume can be described by an equation set with the same structure as that which describes each element in a spectral model. The equivalence relies on an ansatz for detrained condensate introduced by Yanai et al. (1973) and on a simplified microphysics. There are also conceptual differences in the closure of bulk and spectral parameterizations. In particular, we show that the convective quasi-equilibrium closure of Arakawa and Schubert (1974) for spectral parameterizations cannot be carried over to a bulk parameterization in a straightforward way. Quasi-equilibrium of the cloud work function assumes a timescale separation between a slow forcing process and a rapid convective response. But, for the natural bulk analogue to the cloud-work function (the dilute CAPE), the relevant forcing is characterised by a different timescale, and so its quasi-equilibrium entails a different physical constraint. Closures of bulk parameterization that use the non-entraining parcel value of CAPE do not suffer from this timescale issue. However, the Yanai et al. (1973) ansatz must be invoked as a necessary ingredient of those closures.

Physico-chemical controls on phosphorus cycling in two lowland streams. Part 1 - the water column

This paper investigates phosphorus (P) transport and transformation dynamics in two contrasting sub-catchments of the River Kennel, England. Samples were collected daily under baseflow and hourly under stormflow conditions using autosamplers for 2 years and analysed for a range of determinands (full P fractionation, suspended sediment (SS), cations, pH, alkalinity, temperature and oxygen). Concentrations of SRP, SUP, PP and SS were higher in the flashy River Enborne (means of 0.186, 0.071, 0.101 and 34 mg l(-1), respectively) than the groundwater-fed River Lambourn (0.079, 0.057, 0.028 and 9 mg l(-1), respectively). A seasonal trend in the daily P dataset was evident, with lower concentrations during intermediate flows and the spring (caused by a dilution effect and macrophyte uptake) than during baseflow conditions. However, in the hourly P dataset, highest concentrations were observed during storm events in the autumn and winter (reflecting higher scour with increased capacity to entrain particles). Storm events were more significant in contributing to the total P load in the River Enborne than the River Lambourn, especially during August to October, when dry antecedent conditions were observed in the catchment. Re-suspension of P-rich sediment that accumulated within the channel during summer low flows might account for these observations. It is suggested that a P-calcite co-precipitation mechanism was operating during summer in the River Lambourn, while adsorption by metal oxyhydroxide groups was an important mechanism controlling P fractionation in the River Enborne. The influence of flow conditions and channel storage/release mechanisms on P dynamics in these two lowland rivers is assessed. (C) 2004 Elsevier B.V. All rights reserved.

Physico-chemical controls on phosphorus cycling in two lowland streams. Part 2 - The sediment phase

This article investigates the temporal and spatial controls on sediment-phosphorus (P) dynamics in two contrasting sub-catchments of the River Kennet, England. Suspended sediment (collected under representative flow conditions) and size-fractionated bedload (collected weekly for one year) from the Rivers Lambourn and Enborne was analysed for a range of physico-chemical determinands. Total P concentrations were highest in the most mobile fractions of sediment: suspended sediment, fine silt and clay and organic matter (mean concentrations of 1758, 1548 and 1440 mug P g(-1) dry sediment, respectively). Correlation analysis showed significant relationships between total P and total iron (n = 110), total manganese (n = 110), organic matter (n = 110) and specific surface area (n = 28) in the Lambourn (r(2) 0.71, 0.68, 0.62 and 0.52, respectively) and between total P and total iron (n = 110), total manganese (n = 110) and organic matter (n = 110) in the Enborne (r(2) 0.74, 0.85 and 0.68, respectively). These data highlight the importance of metal oxyhydroxide adsorption of P on fine particulates and organic matter. However, high total P concentrations in the granule gravel and coarse sand size fraction during the summer period (mean concentration 228 mug P g(-1) dry sediment) also highlight the role of calcite co-precipitation on P dynamics in the Lambourn. P to cation ratios in Lambourn sediment indicated that fine silt and clay and granule gravel and coarse sand size fractions were potential sources of P release to the water column during specific periods of the summer and autumn. In the Enborne, however, only the granule gravel and coarse sand size fraction had high ratios and a slow, constant release of P was observed. In addition, scanning electron microscopy work confirmed the association of P with calcite in the Lambourn and P with iron on clay particles in the Enborne. The study highlighted the importance of the chemical and physical properties of the sediment in influencing the mechanisms controlling P storage and release within river channels. (C) 2004 Elsevier B.V. All rights reserved.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

Remote estimation of thermal inertia and soil heat flux for bare soil

A method is presented which allows thermal inertia (the soil heat capacity times the square root of the soil thermal diffusivity, C(h)rootD(h)), to be estimated remotely from micrometeorological observations. The method uses the drop in surface temperature, T-s, between sunset and sunrise, and the average night-time net radiation during that period, for clear, still nights. A Fourier series analysis was applied to analyse the time series of T-s . The Fourier series constants, together with the remote estimate of thermal inertia, were used in an analytical expression to calculate diurnal estimates of the soil heat flux, G. These remote estimates of C(h)rootD(h) and G compared well with values derived from in situ sensors. The remote and in situ estimates of C(h)rootD(h) both correlated well with topsoil moisture content. This method potentially allows area-average estimates of thermal inertia and soil heat flux to be derived from remote sensing, e.g. METEOSAT Second Generation, where the area is determined by the sensor's height and viewing angle. (C) 2003 Elsevier B.V. All rights reserved.