Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequent in practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of the sample, thus producing a transfer of mass from some components to other ones, but preserving the total mass present in the system. This evolution is traditionally modelled as a system of ordinary di erential equations of the mass of each component. However, this kind of evolution can be decomposed into a compositional change, expressed in terms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despite of some subcompositions behaving linearly. The goal is to study the characteristics of such simplicial systems of di erential equa- tions such as linearity and stability. This is performed extracting the compositional dif ferential equations from the mass equations. Then, simplicial derivatives are expressed in coordinates of the simplex, thus reducing the problem to the standard theory of systems of di erential equations, including stability. The characterisation of stability of these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and the associated behaviour of the orbits are the main tools. For a three component system, these orbits can be plotted both in coordinates of the simplex or in a ternary diagram. A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is a radioactive decay

Validation of ocean model syntheses against hydrography using a new web application

Results are presented from a new web application called OceanDIVA - Ocean Data Intercomparison and Visualization Application. This tool reads hydrographic profiles and ocean model output and presents the data on either depth levels or isotherms for viewing in Google Earth, or as probability density functions (PDFs) of regional model-data misfits. As part of the CLIVAR Global Synthesis and Observations Panel, an intercomparison of water mass properties of various ocean syntheses has been undertaken using OceanDIVA. Analysis of model-data misfits reveals significant differences between the water mass properties of the syntheses, such as the ability to capture mode water properties.

The release of silicon from amorphous silica and rice straw in Sri Lankan soils

Silicon release from rice straw and amorphous silica when shaken in solution with five Sri Lankan soils was studied indirectly using sorption isotherms and changes in concentration and directly using straw in dialysis bags examined using electron microscopy. The aim was to further our understanding of the processes and factors affecting the release of straw-Si in soils and its availability to rice. The soils (alfisols and ultisols) shaken with 0.1 M NaCl (5 g per 125 mL for 250 days) produced concentrations of 1 - 4 mg L-1 of monosilicic acid-Si. Amorphous silica added to these suspensions (36.5 mg, containing 17 mg Si) raised the concentrations to 20 - 40 mg L-1, and added rice straw (0.5 g, containing 17 mg Si) gave 10 - 25 mg L-1. Sorption isotherms (7 days equilibrations) were used to calculate from the concentrations the amounts of Si released ( 24 - 38% and 8 - 21%, respectively). Both materials gave about 40 mg L-1 of monosilicic acid-Si plus 30 mg L-1 of disilicic acid-Si when shaken in solution alone (5 g per 125 mL). Straw in dialysis bags ( 0.5 g per 25 mL in 0.1 M NaCl) was shaken in soil suspension ( 5 g per 100 mL) for 60 days. Similar concentrations and releases were measured to those obtained above. About one fifth of the mass of straw was lost by decomposition in the first 15 days. A chloroform treatment prevented decomposition, but Si release was unaffected. Disintegration continued throughout the experiments, with phytoliths being exposed and dissolved. Compared to the rate of release from straw into solution without soil, the release of Si into soil suspensions was increased during the first 20 days by adsorption on the soil, but was then reduced probably through the effect of Fe and Al on the phytolith surfaces. The extent of this blocking effect varied between soils and was not simply related to soil pH.

A review of dissolved oxygen modelling techniques for lowland rivers

This review introduces the methods used to simulate the processes affecting dissolved oxygen (DO) in lowland rivers. The important processes are described and this provides a modelling framework to describe those processes in the context of a mass-balance model. The process equations that are introduced all require (reaction) rate parameters and a variety of common procedures for identifying those parameters are reviewed. This is important because there is a wide range of estimation techniques for many of the parameters. These different techniques elicit different estimates of the parameter value and so there is the potential for a significant uncertainty in the model's inputs and therefore in the output too. Finally, the data requirements for modelling DO in lowland rivers are summarised on the basis of modelling the processes described in this review using a mass-balance model. This is reviewed with regard to what data are available and from where they might be obtained. (C) 2003 Elsevier Science B.V. All rights reserved.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Causes of the apparent scale independence of fractal indices associated with forest fragmentation in Bolivia

In a recent investigation, Landsat TM and ETM+ data were used to simulate different resolutions of remotely-sensed images (from 30 to 1100 m) and to analyze the effect of resolution on a range of landscape metrics associated with spatial patterns of forest fragmentation in Chapare, Bolivia since the mid-1980s. Whereas most metrics were found to be highly dependent on pixel size, several fractal metrics (DLFD, MPFD, and AWMPFD) were apparently independent of image resolution, in contradiction with a sizeable body of literature indicating that fractal dimensions of natural objects depend strongly on image characteristics. The present re-analysis of the Chapare images, using two alternative algorithms routinely used for the evaluation of fractal dimensions, shows that the values of the box-counting and information fractal dimensions are systematically larger, sometimes by as much as 85%, than the "fractal" indices DLFD, MPFD, and AWMFD for the same images. In addition, the geometrical fractal features of the forest and non-forest patches in the Chapare region strongly depend on the resolution of images used in the analysis. The largest dependency on resolution occurs for the box-counting fractal dimension in the case of the non-forest patches in 1993, where the difference between the 30 and I 100 m-resolution images corresponds to 24% of the full theoretical range (1.0 to 2.0) of the mass fractal dimension. The observation that the indices DLFD, MPFD, and AWMPFD, unlike the classical fractal dimensions, appear relatively unaffected by resolution in the case of the Chapare images seems due essentially to the fact that these indices are based on a heuristic, "non-geometric" approach to fractals. Because of their lack of a foundation in fractal geometry, nothing guarantees that these indices will be resolution-independent in general. (C) 2006 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.

DISTRIBUTION OF Cr, Pb, Cu, Cd AND Zn IN SEDIMENTS OF THE DOBCZYCE DAM RESERVOIR (SOUTHERN POLAND)

Sediments play a fundamental role in the behaviour of contaminants in aquatic systems. Various processes in sediments, eg adsorption-desorption, oxidation-reduction, ion exchange or biological activities, can cause accumulation or release of metals and anions from the bottom of reservoirs, and have been recently studied in Polish waters [1-3]. Sediment samples from layer A: (1 divided by 6 cm depth in direct contact with bottom water); layer B: (7 divided by 12 cm depth moderate contact); and layer C: (12+ cm depth, in theory an inactive layer) were collected in September 2007 from six sites representing different types of hydrological conditions along the Dobczyce Reservoir (Fig. l). Water depths at the sampling points varied from 3.5 to 21 m. We have focused on studying the distribution and accumulation of several heavy metals (Cr, Pb, Cd, Cu and Zn) in the sediments. The surface, bottom and pore water (extracted from sediments by centrifugation) samples were also collected. Possible relationships between the heavy-metal distribution in sediments and the sediment characteristics (mineralogy, organic matter) as well as the Fe, Mn and Ca content of sediments, have been studied. The 02 concentrations in water samples were also measured. The heavy metals in sediments ranged from 19.0 to 226.3 mg/kg of dry mass (ppm). The results show considerable variations in heavy-metal concentrations between the 6 stations, but not in the individual layers (A, B, C). These variations are related to the mineralogy and chemical composition of the sediments and their pore waters.

Long-term change, interannual and intra-seasonal variability in climate and glacier mass balance in the Central Greater Caucasus, Russia

Long-term trends, interannual and intra-seasonal variability in the mass-balance record from Djankuat glacier, central Greater Caucasus, Russia, are related to local climate change, synoptic and large-scale anomalies in atmospheric circulation. A clear warming signal emerged in the central Greater Caucasus in the early 1990s, leading to a strong increase in ablation. In the absence of a compensating change in winter accumulation, the net mass balance of Djankuat has declined. The highest value of seasonal ablation on record was registered in the summer of 2000. At the beginning of the 21st century these trends reversed. Ablation was below average even in the summer of 2003, which was unusually warm in western Europe. Precipitation and winter accumulation were high, allowing for a partial recovery of net mass balance. The interannual variability in the components of mass balance is weakly related to the North Atlantic Oscillation (NAO) and the Scandinavian teleconnection patterns, but there is a clear link with the large-scale circulation anomalies represented by the Rossby pattern. Five synoptic categories have been identified for the ablation season of 2005, revealing a strong separation between components of radiation budget, air temperature and daily melt. Air temperature is the main control over melt. The highest values of daily ablation are related to the strongly positive NAO which forces high net radiation, and to the warm and moist advection from the Black Sea.

The release of silicon from amorphous silica and rice straw in Sri Lankan soils

Silicon release from rice straw and amorphous silica when shaken in solution with five Sri Lankan soils was studied indirectly using sorption isotherms and changes in concentration and directly using straw in dialysis bags examined using electron microscopy. The aim was to further our understanding of the processes and factors affecting the release of straw-Si in soils and its availability to rice. The soils (alfisols and ultisols) shaken with 0.1 M NaCl (5 g per 125 mL for 250 days) produced concentrations of 1 - 4 mg L-1 of monosilicic acid-Si. Amorphous silica added to these suspensions (36.5 mg, containing 17 mg Si) raised the concentrations to 20 - 40 mg L-1, and added rice straw (0.5 g, containing 17 mg Si) gave 10 - 25 mg L-1. Sorption isotherms (7 days equilibrations) were used to calculate from the concentrations the amounts of Si released ( 24 - 38% and 8 - 21%, respectively). Both materials gave about 40 mg L-1 of monosilicic acid-Si plus 30 mg L-1 of disilicic acid-Si when shaken in solution alone (5 g per 125 mL). Straw in dialysis bags ( 0.5 g per 25 mL in 0.1 M NaCl) was shaken in soil suspension ( 5 g per 100 mL) for 60 days. Similar concentrations and releases were measured to those obtained above. About one fifth of the mass of straw was lost by decomposition in the first 15 days. A chloroform treatment prevented decomposition, but Si release was unaffected. Disintegration continued throughout the experiments, with phytoliths being exposed and dissolved. Compared to the rate of release from straw into solution without soil, the release of Si into soil suspensions was increased during the first 20 days by adsorption on the soil, but was then reduced probably through the effect of Fe and Al on the phytolith surfaces. The extent of this blocking effect varied between soils and was not simply related to soil pH.

Correlation of magnetic field intensities and solar wind speeds of events observed by ACE

The relationship between the magnetic field intensity and speed of solar wind events is examined using ∼3 years of data from the ACE spacecraft. No preselection of coronal mass ejections (CMEs) or magnetic clouds is carried out. The correlation between the field intensity and maximum speed is shown to increase significantly when |B| > 18 nT for 3 hours or more. Of the 24 events satisfying this criterion, 50% are magnetic clouds, the remaining half having no ordered field structure. A weaker correlation also exists between southward magnetic field and speed. Sixteen of the events are associated with halo CMEs leaving the Sun 2 to 4 days prior to the leading edge of the events arriving at ACE. Events selected by speed thresholds show no significant correlation, suggesting different relations between field intensity and speed for fast solar wind streams and ICMEs.

Airborne measurements of the spatial distribution of aerosol chemical composition across Europe and evolution of the organic fraction

The spatial distribution of aerosol chemical composition and the evolution of the Organic Aerosol (OA) fraction is investigated based upon airborne measurements of aerosol chemical composition in the planetary boundary layer across Europe. Sub-micron aerosol chemical composition was measured using a compact Time-of-Flight Aerosol Mass Spectrometer (cToF-AMS). A range of sampling conditions were evaluated, including relatively clean background conditions, polluted conditions in North-Western Europe and the near-field to far-field outflow from such conditions. Ammonium nitrate and OA were found to be the dominant chemical components of the sub-micron aerosol burden, with mass fractions ranging from 20--50% each. Ammonium nitrate was found to dominate in North-Western Europe during episodes of high pollution, reflecting the enhanced NO_x and ammonia sources in this region. OA was ubiquitous across Europe and concentrations generally exceeded sulphate by 30--160%. A factor analysis of the OA burden was performed in order to probe the evolution across this large range of spatial and temporal scales. Two separate Oxygenated Organic Aerosol (OOA) components were identified; one representing an aged-OOA, termed Low Volatility-OOA and another representing fresher-OOA, termed Semi Volatile-OOA on the basis of their mass spectral similarity to previous studies. The factors derived from different flights were not chemically the same but rather reflect the range of OA composition sampled during a particular flight. Significant chemical processing of the OA was observed downwind of major sources in North-Western Europe, with the LV-OOA component becoming increasingly dominant as the distance from source and photochemical processing increased. The measurements suggest that the aging of OA can be viewed as a continuum, with a progression from a less oxidised, semi-volatile component to a highly oxidised, less-volatile component. Substantial amounts of pollution were observed far downwind of continental Europe, with OA and ammonium nitrate being the major constituents of the sub-micron aerosol burden. Such anthropogenically perturbed air masses can significantly perturb regional climate far downwind of major source regions.

Redefinition of the kilogram, ampere, kelvin and mole: a proposed approach to implementing CIPM recommendation 1 (CI-2005)

The International System of Units (SI) is founded on seven base units, the metre, kilogram, second, ampere, kelvin, mole and candela corresponding to the seven base quantities of length, mass, time, electric current, thermodynamic temperature, amount of substance and luminous intensity. At its 94th meeting in October 2005, the International Committee for Weights and Measures (CIPM) adopted a recommendation on preparative steps towards redefining the kilogram, ampere, kelvin and mole so that these units are linked to exactly known values of fundamental constants. We propose here that these four base units should be given new definitions linking them to exactly defined values of the Planck constant h, elementary charge e, Boltzmann constant k and Avogadro constant NA, respectively. This would mean that six of the seven base units of the SI would be defined in terms of true invariants of nature. In addition, not only would these four fundamental constants have exactly defined values but also the uncertainties of many of the other fundamental constants of physics would be either eliminated or appreciably reduced. In this paper we present the background and discuss the merits of these proposed changes, and we also present possible wordings for the four new definitions. We also suggest a novel way to define the entire SI explicitly using such definitions without making any distinction between base units and derived units. We list a number of key points that should be addressed when the new definitions are adopted by the General Conference on Weights and Measures (CGPM), possibly by the 24th CGPM in 2011, and we discuss the implications of these changes for other aspects of metrology.

The propagation mechanism of a vortex on the β plane

The propagation velocity and propagation mechanism for vortices on a β plane are determined for a reduced-gravity model by integrating the momentum equations over the β plane. Isolated vortices, vortices in a background current, and initial vortex propagation from rest are studied. The propagation mechanism for isolated anticyclones as well as cyclones, which has been lacking up to now, is presented. It is shown that, to first order, the vortex moves to generate a Coriolis force on the mass anomaly of the vortex to compensate for the force on the vortex due to the variation of the Coriolis parameter. Only the mass anomaly of the vortex is of importance, because the Coriolis force due to the motion of the bulk of the layer moving with the vortex is almost fully compensated by the Coriolis force on the motion of the exterior flow. Because the mass anomaly of a cyclone is negative the force and acceleration have opposite sign. The role of dipolar structures in steadily moving vortices is discussed, and it is shown that their overall structure is fixed by the steady westward motion of the mass anomaly. Furthermore, it is shown that reduced-gravity vortices are not advected with a background flow. The reason for this behavior is that the background flow changes the ambient vorticity gradient such that the vortex obtains an extra self-propagation term that exactly cancels the advection by the background flow. Last, it is shown that a vortex initially at rest will accelerate equatorward first, after which a westward motion is generated. This result is independent of the sign of the vortex.

Global estimates of snow water equivalent from passive microwave instruments: history, challenges and future developments

Snow properties have been retrieved from satellite data for many decades. While snow extent is generally felt to be obtained reliably from visible-band data, there is less confidence in the measurements of snow mass or water equivalent derived from passive microwave instruments. This paper briefly reviews historical passive microwave instruments and products, and compares the large-scale patterns from these sources to those of general circulation models and leading reanalysis products. Differences are seen to be large between the datasets, particularly over Siberia. A better understanding of the errors in both the model-based and measurement-based datasets is required to exploit both fully. Techniques to apply to the satellite measurements for improved large-scale snow data are suggested.

Forward modelling to determine the observational signatures of white-light imaging and interplanetary scintillation for the propagation of an interplanetary shock in the ecliptic plane

Recent coordinated observations of interplanetary scintillation (IPS) from the EISCAT, MERLIN, and STELab, and stereoscopic white-light imaging from the two heliospheric imagers (HIs) onboard the twin STEREO spacecraft are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun–Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach.