The diabatic contour advective semi-Lagrangian model

This article describes a novel algorithmic development extending the contour advective semi-Lagrangian model to include nonconservative effects. The Lagrangian contour representation of finescale tracer fields, such as potential vorticity, allows for conservative, nondiffusive treatment of sharp gradients allowing very high numerical Reynolds numbers. It has been widely employed in accurate geostrophic turbulence and tracer advection simulations. In the present, diabatic version of the model the constraint of conservative dynamics is overcome by including a parallel Eulerian field that absorbs the nonconservative ( diabatic) tendencies. The diabatic buildup in this Eulerian field is limited through regular, controlled transfers of this field to the contour representation. This transfer is done with a fast newly developed contouring algorithm. This model has been implemented for several idealized geometries. In this paper a single-layer doubly periodic geometry is used to demonstrate the validity of the model. The present model converges faster than the analogous semi-Lagrangian models at increased resolutions. At the same nominal spatial resolution the new model is 40 times faster than the analogous semi-Lagrangian model. Results of an orographically forced idealized storm track show nontrivial dependency of storm-track statistics on resolution and on the numerical model employed. If this result is more generally applicable, this may have important consequences for future high-resolution climate modeling.

On the spectra and pseudospectra of a class of non-self-adjoint random matrices and operators

In this paper we develop and apply methods for the spectral analysis of non-selfadjoint tridiagonal infinite and finite random matrices, and for the spectral analysis of analogous deterministic matrices which are pseudo-ergodic in the sense of E. B. Davies (Commun. Math. Phys. 216 (2001), 687–704). As a major application to illustrate our methods we focus on the “hopping sign model” introduced by J. Feinberg and A. Zee (Phys. Rev. E 59 (1999), 6433–6443), in which the main objects of study are random tridiagonal matrices which have zeros on the main diagonal and random ±1’s as the other entries. We explore the relationship between spectral sets in the finite and infinite matrix cases, and between the semi-infinite and bi-infinite matrix cases, for example showing that the numerical range and p-norm ε - pseudospectra (ε > 0, p ∈ [1,∞] ) of the random finite matrices converge almost surely to their infinite matrix counterparts, and that the finite matrix spectra are contained in the infinite matrix spectrum Σ. We also propose a sequence of inclusion sets for Σ which we show is convergent to Σ, with the nth element of the sequence computable by calculating smallest singular values of (large numbers of) n×n matrices. We propose similar convergent approximations for the 2-norm ε -pseudospectra of the infinite random matrices, these approximations sandwiching the infinite matrix pseudospectra from above and below.

Spatial variability of flow statistics within regular building arrays

Turbulence statistics obtained by direct numerical simulations are analysed to investigate spatial heterogeneity within regular arrays of building-like cubical obstacles. Two different array layouts are studied, staggered and square, both at a packing density of λp=0.25 . The flow statistics analysed are mean streamwise velocity ( u− ), shear stress ( u′w′−−−− ), turbulent kinetic energy (k) and dispersive stress fraction ( u˜w˜ ). The spatial flow patterns and spatial distribution of these statistics in the two arrays are found to be very different. Local regions of high spatial variability are identified. The overall spatial variances of the statistics are shown to be generally very significant in comparison with their spatial averages within the arrays. Above the arrays the spatial variances as well as dispersive stresses decay rapidly to zero. The heterogeneity is explored further by separately considering six different flow regimes identified within the arrays, described here as: channelling region, constricted region, intersection region, building wake region, canyon region and front-recirculation region. It is found that the flow in the first three regions is relatively homogeneous, but that spatial variances in the latter three regions are large, especially in the building wake and canyon regions. The implication is that, in general, the flow immediately behind (and, to a lesser extent, in front of) a building is much more heterogeneous than elsewhere, even in the relatively dense arrays considered here. Most of the dispersive stress is concentrated in these regions. Considering the experimental difficulties of obtaining enough point measurements to form a representative spatial average, the error incurred by degrading the sampling resolution is investigated. It is found that a good estimate for both area and line averages can be obtained using a relatively small number of strategically located sampling points.

Sufficiency of Favard's condition for a class of band-dominated operators on the axis

The purpose of this paper is to show that, for a large class of band-dominated operators on $\ell^\infty(Z,U)$, with $U$ being a complex Banach space, the injectivity of all limit operators of $A$ already implies their invertibility and the uniform boundedness of their inverses. The latter property is known to be equivalent to the invertibility at infinity of $A$, which, on the other hand, is often equivalent to the Fredholmness of $A$. As a consequence, for operators $A$ in the Wiener algebra, we can characterize the essential spectrum of $A$ on $\ell^p(Z,U)$, regardless of $p\in[1,\infty]$, as the union of point spectra of its limit operators considered as acting on $\ell^p(Z,U)$.

Condition number estimates for combined potential boundary integral operators in acoustic scattering

We study the classical combined field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle, namely the indirect formulation due to Brakhage-Werner/Leis/Panic, and the direct formulation associated with the names of Burton and Miller. We obtain lower and upper bounds on the condition numbers for these formulations, emphasising dependence on the frequency, the geometry of the scatterer, and the coupling parameter. Of independent interest we also obtain upper and lower bounds on the norms of two oscillatory integral operators, namely the classical acoustic single- and double-layer potential operators.

The timing of Quaternary calcrete development in semi-arid southeast Spain: Investigating the role of climate on calcrete genesis

Quaternary-aged calcrete horizons are common weathering products in arid and semi-arid regions. It is, however, unclear how calcrete forming processes respond to the major oscillations in climate that occur over the Quaternary period. This paper presents a U-series-based calcrete age database from the Sorbas basin, southeast Spain. The study constructs an age frequency distribution of these ages which is consequently compared to a range of palaeoenvironmental records from the Mediterranean. The age distribution presented here suggests that the formation of pedogenic calcrete horizons in the Sorbas basin primarily occurs during 'warm' isotope stages (MIS 1 and 5), with very few calcrete ages occurring during cold glacial/stadial stages (MIS 2, 3 and 4). It is suggested that this is a function of the environments that existed during 'warm' isotope stages being more conducive to calcrete development than those that existed during cold climate episodes. In a semi-arid region such as the Sorbas basin it is likely that increased aridity during glacial stages, coupled with reduced vegetation and accelerated landscape instability, was crucial in reducing rates of calcrete formation. In a semi-arid region such as southeast Spain, calcrete formation during the Quaternary, therefore, oscillates with climate change but is primarily a "warm" episode phenomenon. It is suggested that further studies are required to see how calcrete genesis responds to environmental change in more humid parts of the Mediterranean. (C) 2009 Elsevier B.V. All rights reserved.

A semi-quantitative approach to deriving a model structure for the uptake of organic chemicals by vegetation

An expert elicitation exercise was undertaken to determine those components and processes that are most important for modeling plant uptake of organic chemicals. The state of our knowledge of these processes was also assessed. This semi-quantitative analysis allowed the construction of an idealized model with seven compartments; soil bulk, soil water, roots, stem, leaves, fruit, and air. Three main areas were identified further research: 1) the uptake of organic chemicals by fruit; 2) the internal transfer of organic chemicals between plant structures (e.g., stem and leaves); and 3) the transfer via the soil-air-plant pathway. Until new data becomes available to quantify these processes, it is proposed that an equilibrium partitioning approach is used between plant components other than fruit or that models consist of both an edible and inedible compartment.

Effects of manure and fertilizer on grain yield, soil carbon and phosphorus in a 13-year field trial in semi-arid Kenya

Long-term indicators of soil fertility were assessed by measuring grain yield, soil organic carbon (SOC) and soil Olsen phosphorous for a P-deficient soil. In one set of treatments, goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1), and intercrops of sorghum/cowpea, millet/green gram and maize/pigeonpea were grown. Yield depended on rainfall and trends with time were not identifiable. Manure caused an upward trend in SOC, but 10 t ha(-1) manure did not give significantly more SOC than 5 t ha(-1). Only 10 t ha(-1) manure increased Olsen P. Measurements of both SOC and Olsen P are recommended. In another set of treatments, manure was applied for four years; the residual effect lasted another seven to eight years when assessed by yield, SOC and Olsen P Treatment with mineral fertilizers provided the same rates of N and P as 5 t hat manure and yields from manure and fertilizer were similar. Fertilizer increased Olsen P but not SOC. Management systems with occasional manure application and intermediate fertilizer applications should be assessed. Inputs and offtakes of C, N and P were measured for three years. Approximately 16, 25 and 11% of C, N and P respectively were stabilized into soil organic matter from 5 t ha(-1) a(-1) manure. The majority of organic P was fixed as soil inorganic P.

Effect of manure application on crop yield and soil chemical properties in a long-term field trial of semi-arid Kenya

The sustainability of cereal/legume intercropping was assessed by monitoring trends in grain yield, soil organic C (SOC) and soil extractable P (Olsen method) measured over 13 years at a long-term field trial on a P-deficient soil in semi-arid Kenya. Goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1) and trends in grain yield were not identifiable because of season-to-season variations. SOC and Olsen P increased for the first seven years of manure application and then remained constant. The residual effect of manure applied for four years only lasted another seven to eight years when assessed by yield, SOC and Olsen P. Mineral fertilizers provided the same annual rates of N and P as in 5 t ha(-1) manure and initially ,gave the same yield as manure, declining after nine years to about 80%. Therefore, manure applications could be made intermittently and nutrient requirements topped-up with fertilizers. Grain yields for sorghum with continuous manure were described well by correlations with rainfall and manure input only, if data were excluded for seasons with over 500 mm rainfall. A comprehensive simulation model should correctly describe crop losses caused by excess water.

A contour-advective semi-lagrangian numerical algorithm for simulating fine-scale conservative dynamical fields

This paper describes a novel numerical algorithm for simulating the evolution of fine-scale conservative fields in layer-wise two-dimensional flows, the most important examples of which are the earth's atmosphere and oceans. the algorithm combines two radically different algorithms, one Lagrangian and the other Eulerian, to achieve an unexpected gain in computational efficiency. The algorithm is demonstrated for multi-layer quasi-geostrophic flow, and results are presented for a simulation of a tilted stratospheric polar vortex and of nearly-inviscid quasi-geostrophic turbulence. the turbulence results contradict previous arguments and simulation results that have suggested an ultimate two-dimensional, vertically-coherent character of the flow. Ongoing extensions of the algorithm to the generally ageostrophic flows characteristic of planetary fluid dynamics are outlined.