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 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.

Comparison of EDTA and EDDS as potential soil amendments for enhanced phytoextraction of heavy metals

Phytoextraction has been proposed as an alternative remediation technology for soils polluted with heavy metals or radionuclides, but is generally conceived as too slow working. Enhancing the accumulation of trace pollutants in harvestable plant tissues is a prerequisite for the technology to be practically applicable. The chelating aminopolycarboxylic acid, ethylene diamine tetraacetate (EDTA), has been found to enhance shoot accumulation of heavy metals. However, the use of EDTA in phytoextraction may not be suitable due to its high environmental persistence, which may lead to groundwater contamination. This paper aims to assess whether ethylene diamine disuccinate (EDDS), a biodegradable chelator, can be used for enhanced phytoextraction purposes. A laboratory experiment was conducted to examine mobilisation of Cd, Cu, Cr, Ni, Pb and Zn into the soil solution upon application of EDTA or EDDS. The longevity of the induced mobilisation was monitored for a period of 40 days after application. Estimated effect half lives ranged between 3.8 and 7.5 days for EDDS, depending on the applied dose. The minimum observed effect half life of EDTA was 36 days, while for the highest applied dose no decrease was observed throughout the 40 day period of the mobilisation experiment. Performance of EDTA and EDDS for phytoextraction was evaluated by application to Helianthus annuus. Two other potential chelators, known for their biodegradability in comparison to EDTA, were tested in the plant experiment: nitrilo acetic acid (NTA) and citric acid. Uptake of heavy metals was higher in EDDS-treated pots than in EDTA-treated pots. The effects were still considered insufficiently high to consider efficient remediation. This may be partly due to the choice of timing for application of the soil amendment. Fixing the time of application at an earlier point before harvest may yield better results. NTA and citric acid induced no significant effects on heavy metal uptake. (C) 2004 Elsevier Ltd. All rights reserved.

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.

A well-posed integral equation formulation for three-dimensional rough surface scattering

We consider the problem of scattering of time-harmonic acoustic waves by an unbounded sound-soft rough surface. Recently, a Brakhage Werner type integral equation formulation of this problem has been proposed, based on an ansatz as a combined single- and double-layer potential, but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half-space Green's function. Moreover, it has been shown in the three-dimensional case that this integral equation is uniquely solvable in the space L-2 (Gamma) when the scattering surface G does not differ too much from a plane. In this paper, we show that this integral equation is uniquely solvable with no restriction on the surface elevation or slope. Moreover, we construct explicit bounds on the inverse of the associated boundary integral operator, as a function of the wave number, the parameter coupling the single- and double-layer potentials, and the maximum surface slope. These bounds show that the norm of the inverse operator is bounded uniformly in the wave number, kappa, for kappa > 0, if the coupling parameter h is chosen proportional to the wave number. In the case when G is a plane, we show that the choice eta = kappa/2 is nearly optimal in terms of minimizing the condition number.

Approximate factorization constraint preconditioners for saddle-point matrices

We consider the application of the conjugate gradient method to the solution of large, symmetric indefinite linear systems. Special emphasis is put on the use of constraint preconditioners and a new factorization that can reduce the number of flops required by the preconditioning step. Results concerning the eigenvalues of the preconditioned matrix and its minimum polynomial are given. Numerical experiments validate these conclusions.

Acoustic scattering by mildly rough unbounded surfaces in three dimensions

For a nonlocally perturbed half- space we consider the scattering of time-harmonic acoustic waves. A second kind boundary integral equation formulation is proposed for the sound-soft case, based on a standard ansatz as a combined single-and double-layer potential but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half- space Green's function. Due to the unboundedness of the surface, the integral operators are noncompact. In contrast to the two-dimensional case, the integral operators are also strongly singular, due to the slow decay at infinity of the fundamental solution of the three-dimensional Helmholtz equation. In the case when the surface is sufficiently smooth ( Lyapunov) we show that the integral operators are nevertheless bounded as operators on L-2(Gamma) and on L-2(Gamma G) boolean AND BC(Gamma) and that the operators depend continuously in norm on the wave number and on G. We further show that for mild roughness, i.e., a surface G which does not differ too much from a plane, the boundary integral equation is uniquely solvable in the space L-2(Gamma) boolean AND BC(Gamma) and the scattering problem has a unique solution which satisfies a limiting absorption principle in the case of real wave number.

Upwind solution of singular differential equations arising from steady channel flows

We study ordinary nonlinear singular differential equations which arise from steady conservation laws with source terms. An example of steady conservation laws which leads to those scalar equations is the Saint–Venant equations. The numerical solution of these scalar equations is sought by using the ideas of upwinding and discretisation of source terms. Both the Engquist–Osher scheme and the Roe scheme are used with different strategies for discretising the source terms.