A cell by cell anisotropic adaptive mesh ALE scheme for the numerical solution of the Euler equations

In this paper a cell by cell anisotropic adaptive mesh technique is added to an existing staggered mesh Lagrange plus remap finite element ALE code for the solution of the Euler equations. The quadrilateral finite elements may be subdivided isotropically or anisotropically and a hierarchical data structure is employed. An efficient computational method is proposed, which only solves on the finest level of resolution that exists for each part of the domain with disjoint or hanging nodes being used at resolution transitions. The Lagrangian, equipotential mesh relaxation and advection (solution remapping) steps are generalised so that they may be applied on the dynamic mesh. It is shown that for a radial Sod problem and a two-dimensional Riemann problem the anisotropic adaptive mesh method runs over eight times faster.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

GDASH: a grid-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Very large scale increases in speed of execution can therefore be achieved by distributing individual DASH runs over a network of computers. The GDASH program achieves this by packaging DASH in a form that enables it to run under the Univa UD Grid MP system, which harnesses networks of existing computing resources to perform calculations.

MDASH: a multi-core-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Modest increases in speed of execution can therefore be achieved by executing individual DASH runs on the individual cores of CPUs.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

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.

The influence of mineral solubility and soil solution concentration on the toxicity of copper to Eisenia fetida Savigny

The OECD 14 d earthworm acute toxicity test was used to determine the toxicity of copper added as copper nitrate (Cu(NO3)(2)), copper sulphate (CuSO4) and malachite (Cu-2(OH)(2)(CO3)) to Eisenia fetida Savigny. Cu(NO3)(2), and CuSO4 were applied in both an aqueous (aq) and solid (s) form, Cu-2(OH)(2)(CO3) was added as a solid. Soil solution was extracted by centrifugation, and analysed for copper. Two extractants [0.01 M CaCl2 and 0.005 M diethylenetriminpentaacetic acid (DTPA)] were used as a proxy of the bioavailable copper fraction in the soil. For bulk soil copper content the calculated copper toxicity decreased in the order nitrate > sulphide > carbonate, the same order as decreasing solubility of the metal compounds. For Cu(NO3)(2) and CuSO4, the LC50s obtained were not significantly different when the compound was added in solution or solid form. There was a significant correlation between the soil solution copper concentration and the percentage earthworm mortality for all 3 copper compounds (P less than or equal to 0.05) indicating that the soil pore water copper concentration is important for determining copper availability and toxicity to E. fetida. In soil avoidance tests the earthworms avoided the soils treated with Cu(NO3)(2) (aq and s) and CuSO4 (aq and s), at all concentrations used (110-8750 mug Cu g(-1), and 600-8750 mug Cu g(-1) respectively). In soils treated with Cu-2(OH2)CO3, avoidance behaviour was exhibited at all concentrations greater than or equal to3500 mug Cu g(-1). There was no significant correlation between the copper extracted by either CaCl2 or DTPA and percentage mortality. These two extractants are therefore not useful indicators of copper availability and toxicity to E. fetida.

Zinc deficiency in selected cultivars of wheat and barley as tested in solution culture

The relative zinc (Zn) efficiencies of 33 wheat and 3 barley cultivars were determined by growing them in chelate-buffered culture solutions. Zn efficiency, determined by growth in a Zn-deficient solution relative to that in a medium containing an adequate concentration of Zn, was found to vary between 10% and 63% among the cultivars tested. Out of the 36 cultivars tested, 12 proved to be Zn efficient, 10 were Zn inefficient, and the remaining 14 varieties were classed as intermediate. The most Zn-efficient cultivars included Bakhtawar, Gatcher S61, Wilgoyne, and Madrigal, and the most Zn inefficient included Durati, Songlen, Excalibur, and Chakwal-86. Zn-efficient cultivars accumulated greater amounts of Zn in their shoots than inefficient cultivars, but the correlation between shoot Zn and shoot dry matter production was poor. All the cultivars accumulated higher concentrations of iron (Fe), copper (Cu), manganese (Mn), and phosphorus (P) at deficient levels of Zn, compared with adequate Zn concentrations. The Zn-inefficient cultivars accumulated higher concentrations of these other elements compared to efficient cultivars.

The role of calcite and gypsum in the leaching of potassium in a sandy soil

Pulses of potassium (K+) applied to columns of repacked calcium (Ca2+) saturated soil were leached with distilled water or calcium chloride (CaCl2) solutions of various concentrations at a rate of 12 mm h(-1). With increased Ca2+ concentration, the rate of movement of K+ increased, as did the concentration of K+ in the displaced pulse, which was less dispersed. The movement of K+ in calcite-amended soil leached with water was at a similar rate to that of the untreated soil leached with 1 mM CaCl2, and in soil containing gypsum, movement was similar to that leached with 15 mM CaCl2. The Ca2+ concentrations in the leachates were about 0.4 and 15 mM respectively the expected values for the dissolution of the two amendments. Soil containing native K+ was leached with distilled water or CaCl2 solutions. The amount of K+ leached increased as Ca2+ concentration increased, with up to 34% of the exchangeable K+ being removed in five pore volumes of 15 mM CaCl2. Soil amended with calcite and leached with water lost K+ at a rate between that for leaching the unamended soil with 1 mM CaCl2 and that with water. Soil containing gypsum and leached with water lost K+ at a similar rate to unamended soil leached with 15 mM CaCl2. The presence of Ca2+ in irrigation water and of soil minerals able to release Ca2+ are of importance in determining the amounts of K+ leached from soils. The LEACHM model predicted approximately the displacement of K+, and was more accurate with higher concentrations of displacing solution. The shortcomings of this model are its inability to account for rate-controlled processes and the assumption that K+:Ca2+ exchange during leaching can be described using a constant adsorption coefficient. As a result, the pulse is predicted to appear a little earlier and the following edge has less of a tail than chat measured. In practical agriculture, the model will be more useful in soils containing gypsum or leached with saline water than in either calcareous or non-calcareous soils leached with rainwater.

Wetland restoration within agricultural watersheds: Balancing water quality protection with habitat conservation

Peat wetlands that have been restored from agricultural Land have the potential to act as Long term sources of phosphorus (P) and, therefore have to potenital to accelerate freshwater eutrophication. During a two-year study the water table in a eutrophic fen peat that was managed by pump drainage fluctuated annually between +20 cm and -60 cm relative to ground Level. This precise management was facilitated by the high hydraulic conductivity (K) of the humified peat (1.1 x 10(-5) m s(-1)) below around 60 cm depth. However, during one week of intermittent pumping, as much as 50 g ha(-1) dissolved P entered the pumped ditch. Summer. rainfall events and autumn reflooding also triggered P losses. The P Losses were attributed to the low P sorption capacity (217 mg kg(-1)) of the saturated peat below 60 cm, combined with its high K and the reductive dissolution of Fe bound P.

Plant roots release phospholipid surfactants that modify the physical and chemical properties of soil

Plant root mucilages contain powerful surfactants that will alter the interaction of soil solids with water and ions, and the rates of microbial processes. The lipid composition of maize, lupin and wheat root mucilages was analysed by thin layer chromatography and gas chromatography-mass spectrometry. A commercially available phosphatidylcholine (lecithin), chemically similar to the phospholipid surfactants identified in the mucilages, was then used to evaluate its effects on selected soil properties. The lipids found in the mucilages were principally phosphatidylcholines, composed mainly of saturated fatty acids, in contrast to the lipids extracted from root tissues. In soil at low tension, lecithin reduced the water content at any particular tension by as much as 10 and 50% in soil and acid-washed sand, respectively. Lecithin decreased the amount of phosphate adsorption in soil and increased the phosphate concentration in solution by 10%. The surfactant also reduced net rates of ammonium consumption and nitrate production in soil. These experiments provide the first evidence we are aware of that plant-released surfactants will significantly modify the biophysical environment of the rhizosphere.

Morphology, biomass and nutrient status of fine roots of Scots pine (Pinus sylvestris) as influenced by seasonal fluctuations in soil moisture and soil solution chemistry

A field monitoring study was carried out to follow the changes of fine root morphology, biomass and nutrient status in relation to seasonal changes in soil solution chemistry and moisture regime in a mature Scots pine stand on acid soil. Seasonal and yearly fluctuations in soil moisture and soil solution chemistry have been observed. Changes in soil moisture accounted for some of the changes in the soil solution chemistry. The results showed that when natural acidification in the soil occurs with low pH (3.5-4.2) and high aluminium concentration in the soil solution (> 3-10 mg l(-1)), fine root longevity and distribution could be affected. However, fine root growth of Scots pine may not be negatively influenced by adverse soil chemical conditions if soil moisture is not a limiting factor for root growth. In contrast, dry soil conditions increase Scots pine susceptibility to soil acidification and this could significantly reduce fine root growth and increase root mortality. It is therefore important to study seasonal fluctuations of the environmental variables when investigating and modelling cause-effect relationships.

Short-term effects of manipulated increase in acid deposition on soil, soil solution chemistry and fine roots in Scots pine (Pinus sylvestris) stand on a podzol

A manipulated increase in acid deposition (15 kg S ha(-1)), carried out for three months in a mature Scots pine (Pinus sylvestris) stand on a podzol, acidified the soil and raised dissolved Al at concentrations above the critical level of 5 mg l(-1) previously determined in a controlled experiment with Scots pine seedlings. The induced soil acidification reduced tree fine root density and biomass significantly in the top 15 cm of soil in the field. The results suggested that the reduction in fine root growth was a response not simply to high Al in solution but to the depletion of exchangeable Ca and Mg in the organic layer, K deficiency, the increase in NH4:NO3 ratio in solution and the high proton input to the soil by the acid manipulation. The results from this study could not justify the hypothesis of Al-induced root damage under field conditions, at least not in the short term. However, the study suggests that a short exposure to soil acidity may affect the fine root growth of mature Scots pine.