Properties of singular vectors using convective available potential energy as final time norm

We study the feasibility of using the singular vector technique to create initial condition perturbations for short-range ensemble prediction systems (SREPS) focussing on predictability of severe local storms and in particular deep convection. For this a new final time semi-norm based on the convective available potential energy (CAPE) is introduced. We compare singular vectors using the CAPE-norm with SVs using the more common total energy (TE) norm for a 2-week summer period in 2007, which includes a case of mesoscale extreme rainfall in the south west of Finland. The CAPE singular vectors perturb the CAPE field by increasing the specific humidity and temperature of the parcel and increase the lapse rate above the parcel in the lower troposphere consistent with physical considerations. The CAPE-SVs are situated in the lower troposphere. This in contrast to TE-SVs with short optimization times which predominantly remain in the high troposphere. By examining the time evolution of the CAPE singular values we observe that the convective event in the south west of Finland is clearly associated with high CAPE singular values.

The importance of friction in mountain wave drag amplification by Scorer parameter resonance

A mechanism for amplification of mountain waves, and their associated drag, by parametric resonance is investigated using linear theory and numerical simulations. This mechanism, which is active when the Scorer parameter oscillates with height, was recently classified by previous authors as intrinsically nonlinear. Here it is shown that, if friction is included in the simplest possible form as a Rayleigh damping, and the solution to the Taylor-Goldstein equation is expanded in a power series of the amplitude of the Scorer parameter oscillation, linear theory can replicate the resonant amplification produced by numerical simulations with some accuracy. The drag is significantly altered by resonance in the vicinity of n/l_0 = 2, where l_0 is the unperturbed value of the Scorer parameter and n is the wave number of its oscillation. Depending on the phase of this oscillation, the drag may be substantially amplified or attenuated relative to its non-resonant value, displaying either single maxima or minima, or double extrema near n/l_0 = 2. Both non-hydrostatic effects and friction tend to reduce the magnitude of the drag extrema. However, in exactly inviscid conditions, the single drag maximum and minimum are suppressed. As in the atmosphere friction is often small but non-zero outside the boundary layer, modelling of the drag amplification mechanism addressed here should be quite sensitive to the type of turbulence closure employed in numerical models, or to computational dissipation in nominally inviscid simulations.

Assessment of wheat cultivars for drought tolerance via osmotic stress imposed at early seedling growth stages

A study was conducted in the Department of Plant Breeding and Genetics,Sindh Agriculture University, Tandojam, Pakistan during the year 2009. Sixteen spring wheat cultivars (Triticum aestivum L.) were screened under osmotic stress with three treatments i.e. control-no PEG (polyethylene glycol), 15 percent and 25 percent PEG-6000 solution. The analysis of variance indicated significant differences among treatments for all seedling traits except seed germination percentage. Varieties also differed significantly in germination percentage, coleoptile length, shoot root length, shoot weight, root/shoot ratio and seed vigour index. However, shoot and root weights were non-significant. Significant interactions revealed that cultivars responded variably to osmotic stress treatments; hence provided better opportunity to select drought tolerant cultivars at seedling growth stages. The relative decrease over averages due to osmotic stress was 0.8 percent in seed germination, 53 percent in coleoptile length 62.9 percent in shoot length, 74.4 percent in root length, 50.6 percent in shoot weight, 45.1 percent in root weight, 30.2 percent in root/shoot ratio and 68.5 percent in seed vigour index. However, relative decrease of individual variety for various seedling traits could be more meaningful which indicated that cultivar TD-1 showed no reduction in coleoptile length, while minimum decline was noted in Anmol. For shoot length, cultivar Sarsabz expressed minimum reduction followed by Anmol. However, cultivars Anmol, Moomal, Inqalab-91, and Pavan gave almost equally lower reductions for root length suggesting their higher stress tolerance. In other words, cultivars Anmol, Moomal, Inqalab-91, Sarsabz, TD-1, ZA-77 and Pavan had relatively longer coleoptiles, shoots and roots, and were regarded as drought tolerant. Correlation coefficients among seedlings traits were significant and positive for all traits except germination percentage which had no significant correlation with any of other trait. The results indicated that increase in one trait may cause simultaneous increase in other traits; hence selection for any of these seedling attributes will lead to develop drought tolerant wheat cultivars.

A rare case of solution and solid state inter-conversion of two copper(II) dimers and a copper(II) chain

Three Cu(II)-azido complexes of formula [Cu2L2(N-3)(2)] (1), [Cu2L2(N-3)(2)]center dot H2O (2) and [CuL(N-3)](n) (3) have been synthesized using the same tridentate Schiff base ligand HL (2-[(3-methylaminopropylimino)-methyl]-phenol), the condensation product of N-methyl-1,3-propanediamine and salicyldehyde). Compounds 1 and 2 are basal-apical mu-1,1 double azido bridged dimers. The dimeric structure of 1 is centro-symmetric but that of 2 is non-centrommetric. Compound 3 is a mu-1,1 single azido bridged 1D chain. The three complexes interconvert in solution and can be obtained in pure form by carefully controlling the synthetic conditions. Compound 2 undergoes an irreversible transformation to 1 upon dehydration in the solid state. The magnetic properties of compounds 1 and 2 show the presence of weak antiferromagnetic exchange interactions mediated by the double 1,1-N-3 azido bridges (J = -2.59(4) and -0.10(1) cm-(1), respectively). The single 1,1-N-3 bridge in compound 3 mediates a negligible exchange interaction.

Homoleptic copper(II) and copper(I) complexes of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline. Six-coordinate copper(I) in solution

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Cu(ClO(4))(2)center dot 6H(2)O in methanol in 3:1 M ratio at room temperature yields light green [CuL(3)](ClO(4))(2)center dot H(2)O (1). The X-ray crystal structure of the hemi acetonitrile solvate [CuL(3)](ClO(4))(2)center dot 0.5CH(3)CN has been determined which shows Jahn-Teller distortion in the CuN(6) core present in the cation [CuL(3)](2+). Complex 1 gives an axial EPR spectrum in acetonitrile-toluene glass with g(parallel to) = 2.262 (A(parallel to) = 169 x 10 (4) cm (1)) and g(perpendicular to) = 2.069. The Cu(II/I) potential in 1 in CH(2)Cl(2) at a glassy carbon electrode is 0.32 V versus NHE. This potential does not change with the addition of extra L in the medium implicating generation of a six-coordinate copper(I) species [CuL(3)](+) in solution. B3LYP/LanL2DZ calculations show that the six Cu-N bond distances in [CuL(3)](+) are 2.33, 2.25, 2.32, 2.25, 2.28 and 2.25 angstrom while the ideal Cu(I)-N bond length in a symmetric Cu(I)N(6) moiety is estimated as 2.25 angstrom. Reaction of L with Cu(CH(3)CN)(4)ClO(4) in dehydrated methanol at room temperature even in 4:1 M proportion yields [CuL(2)]ClO(4) (2). Its (1)H NMR spectrum indicates that the metal in [CuL(2)](+) is tetrahedral. The Cu(II/I) potential in 2 is found to be 0.68 V versus NHE in CH(2)Cl(2) at a glassy carbon electrode. In presence of excess L, 2 yields the cyclic voltammogram of 1. From (1)H NMR titration, the free energy of binding of L to [CuL(2)](+) to produce [CuL(3)](+) in CD(2)Cl(2) at 298 K is estimated as -11.7 (+/-0.2) kJ mol (1).

The (de)merits of minimum-variance hedging: application to the crack spread

We study the empirical performance of the classical minimum-variance hedging strategy, comparing several econometric models for estimating hedge ratios of crude oil, gasoline and heating oil crack spreads. Given the great variability and large jumps in both spot and futures prices, considerable care is required when processing the relevant data and accounting for the costs of maintaining and re-balancing the hedge position. We find that the variance reduction produced by all models is statistically and economically indistinguishable from the one-for-one “naïve” hedge. However, minimum-variance hedging models, especially those based on GARCH, generate much greater margin and transaction costs than the naïve hedge. Therefore we encourage hedgers to use a naïve hedging strategy on the crack spread bundles now offered by the exchange; this strategy is the cheapest and easiest to implement. Our conclusion contradicts the majority of the existing literature, which favours the implementation of GARCH-based hedging strategies.

'Quasi'-norm of an arithmetical convolution operator and the order of the Riemann zeta function

In this paper we study Dirichlet convolution with a given arithmetical function f as a linear mapping 'f that sends a sequence (an) to (bn) where bn = Pdjn f(d)an=d. We investigate when this is a bounded operator on l2 and ¯nd the operator norm. Of particular interest is the case f(n) = n¡® for its connection to the Riemann zeta function on the line 1, 'f is bounded with k'f k = ³(®). For the unbounded case, we show that 'f : M2 ! M2 where M2 is the subset of l2 of multiplicative sequences, for many f 2 M2. Consequently, we study the `quasi'-norm sup kak = T a 2M2 k'fak kak for large T, which measures the `size' of 'f on M2. For the f(n) = n¡® case, we show this quasi-norm has a striking resemblance to the conjectured maximal order of j³(® + iT )j for ® > 12 .

Sparse probability density function estimation using the minimum integrated square error

We develop a new sparse kernel density estimator using a forward constrained regression framework, within which the nonnegative and summing-to-unity constraints of the mixing weights can easily be satisfied. Our main contribution is to derive a recursive algorithm to select significant kernels one at time based on the minimum integrated square error (MISE) criterion for both the selection of kernels and the estimation of mixing weights. The proposed approach is simple to implement and the associated computational cost is very low. Specifically, the complexity of our algorithm is in the order of the number of training data N, which is much lower than the order of N2 offered by the best existing sparse kernel density estimators. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with comparable accuracy to those of the classical Parzen window estimate and other existing sparse kernel density estimators.

Scattering of electromagnetic waves by rough interfaces and inhomogeneous layers

We consider a two-dimensional problem of scattering of a time-harmonic electromagnetic plane wave by an infinite inhomogeneous conducting or dielectric layer at the interface between semi-infinite homogeneous dielectric half-spaces. The magnetic permeability is assumed to be a fixed positive constant. The material properties of the media are characterized completely by an index of refraction, which is a bounded measurable function in the layer and takes positive constant values above and below the layer, corresponding to the homogeneous dielectric media. In this paper, we examine only the transverse magnetic (TM) polarization case. A radiation condition appropriate for scattering by infinite rough surfaces is introduced, a generalization of the Rayleigh expansion condition for diffraction gratings. With the help of the radiation condition the problem is reformulated as an equivalent mixed system of boundary and domain integral equations, consisting of second-kind integral equations over the layer and interfaces within the layer. Assumptions on the variation of the index of refraction in the layer are then imposed which prove to be sufficient, together with the radiation condition, to prove uniqueness of solution and nonexistence of guided wave modes. Recent, general results on the solvability of systems of second kind integral equations on unbounded domains establish existence of solution and continuous dependence in a weighted norm of the solution on the given data. The results obtained apply to the case of scattering by a rough interface between two dielectric media and to many other practical configurations.

Acoustic scattering by an inhomogeneous layer on a rigid plate

The problem of scattering of time-harmonic acoustic waves by an inhomogeneous fluid layer on a rigid plate in R2 is considered. The density is assumed to be unity in the media: within the layer the sound speed is assumed to be an arbitrary bounded measurable function. The problem is modelled by the reduced wave equation with variable wavenumber in the layer and a Neumann condition on the plate. To formulate the problem and prove uniqueness of solution a radiation condition appropriate for scattering by infinite rough surfaces is introduced, a generalization of the Rayleigh expansion condition for diffraction gratings. With the help of the radiation condition the problem is reformulated as a system of two second kind integral equations over the layer and the plate. Under additional assumptions on the wavenumber in the layer, uniqueness of solution is proved and the nonexistence of guided wave solutions of the homogeneous problem established. General results on the solvability of systems of integral equations on unbounded domains are used to establish existence and continuous dependence in a weighted norm of the solution on the given data.

Approximate solution of second kind integral equations on infinite cylindrical surfaces

The paper considers second kind integral equations of the form $\phi (x) = g(x) + \int_S {k(x,y)} \phi (y)ds(y)$ (abbreviated $\phi = g + K\phi $), in which S is an infinite cylindrical surface of arbitrary smooth cross section. The “truncated equation” (abbreviated $\phi _a = E_a g + K_a \phi _a $), obtained by replacing S by $S_a $, a closed bounded surface of class $C^2 $, the boundary of a section of the interior of S of length $2a$, is also discussed. Conditions on k are obtained (in particular, implying that K commutes with the operation of translation in the direction of the cylinder axis) which ensure that $I - K$ is invertible, that $I - K_a $ is invertible and $(I - K_a )^{ - 1} $ is uniformly bounded for all sufficiently large a, and that $\phi _a $ converges to $\phi $ in an appropriate sense as $a \to \infty $. Uniform stability and convergence results for a piecewise constant boundary element collocation method for the truncated equations are also obtained. A boundary integral equation, which models three-dimensional acoustic scattering from an infinite rigid cylinder, illustrates the application of the above results to prove existence of solution (of the integral equation and the corresponding boundary value problem) and convergence of a particular collocation method.

Some uniform stability and convergence results for integral equations on the real line and projection methods for their solution

The paper considers second kind equations of the form (abbreviated x=y + K2x) in which and the factor z is bounded but otherwise arbitrary so that equations of Wiener-Hopf type are included as a special case. Conditions on a set are obtained such that a generalized Fredholm alternative is valid: if W satisfies these conditions and I − Kz, is injective for each z ε W then I − Kz is invertible for each z ε W and the operators (I − Kz)−1 are uniformly bounded. As a special case some classical results relating to Wiener-Hopf operators are reproduced. A finite section version of the above equation (with the range of integration reduced to [−a, a]) is considered, as are projection and iterated projection methods for its solution. The operators (where denotes the finite section version of Kz) are shown uniformly bounded (in z and a) for all a sufficiently large. Uniform stability and convergence results, for the projection and iterated projection methods, are obtained. The argument generalizes an idea in collectively compact operator theory. Some new results in this theory are obtained and applied to the analysis of projection methods for the above equation when z is compactly supported and k(s − t) replaced by the general kernel k(s,t). A boundary integral equation of the above type, which models outdoor sound propagation over inhomogeneous level terrain, illustrates the application of the theoretical results developed.

Runge-Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations

Many operational weather forecasting centres use semi-implicit time-stepping schemes because of their good efficiency. However, as computers become ever more parallel, horizontally explicit solutions of the equations of atmospheric motion might become an attractive alternative due to the additional inter-processor communication of implicit methods. Implicit and explicit (IMEX) time-stepping schemes have long been combined in models of the atmosphere using semi-implicit, split-explicit or HEVI splitting. However, most studies of the accuracy and stability of IMEX schemes have been limited to the parabolic case of advection–diffusion equations. We demonstrate how a number of Runge–Kutta IMEX schemes can be used to solve hyperbolic wave equations either semi-implicitly or HEVI. A new form of HEVI splitting is proposed, UfPreb, which dramatically improves accuracy and stability of simulations of gravity waves in stratified flow. As a consequence it is found that there are HEVI schemes that do not lose accuracy in comparison to semi-implicit ones. The stability limits of a number of variations of trapezoidal implicit and some Runge–Kutta IMEX schemes are found and the schemes are tested on two vertical slice cases using the compressible Boussinesq equations split into various combinations of implicit and explicit terms. Some of the Runge–Kutta schemes are found to be beneficial over trapezoidal, especially since they damp high frequencies without dropping to first-order accuracy. We test schemes that are not formally accurate for stiff systems but in stiff limits (nearly incompressible) and find that they can perform well. The scheme ARK2(2,3,2) performs the best in the tests.