A comparison of two methods for developing the linearization of a shallow-water model

We develop the linearization of a semi-implicit semi-Lagrangian model of the one-dimensional shallow-water equations using two different methods. The usual tangent linear model, formed by linearizing the discrete nonlinear model, is compared with a model formed by first linearizing the continuous nonlinear equations and then discretizing. Both models are shown to perform equally well for finite perturbations. However, the asymptotic behaviour of the two models differs as the perturbation size is reduced. This leads to difficulties in showing that the models are correctly coded using the standard tests. To overcome this difficulty we propose a new method for testing linear models, which we demonstrate both theoretically and numerically. © Crown copyright, 2003. Royal Meteorological Society

Case study of the use of remotely sensed data for modeling flood inundation on the river Severn, UK.

A methodology for using remotely sensed data to both generate and evaluate a hydraulic model of floodplain inundation is presented for a rural case study in the United Kingdom: Upton-upon-Severn. Remotely sensed data have been processed and assembled to provide an excellent test data set for both model construction and validation. In order to assess the usefulness of the data and the issues encountered in their use, two models for floodplain inundation were constructed: one based on an industry standard one-dimensional approach and the other based on a simple two-dimensional approach. The results and their implications for the future use of remotely sensed data for predicting flood inundation are discussed. Key conclusions for the use of remotely sensed data are that care must be taken to integrate different data sources for both model construction and validation and that improvements in ground height data shift the focus in terms of model uncertainties to other sources such as boundary conditions. The differences between the two models are found to be of minor significance.

On the evolution of the solar wind between 1 and 5 AU at the time of the Cassini Jupiter flyby: Multispacecraft observations of interplanetary coronal mass ejections including the formation of a merged interaction region

The Cassini flyby of Jupiter occurred at a time near solar maximum. Consequently, the pre-Jupiter data set reveals clear and numerous transient perturbations to the Parker Spiral solar wind structure. Limited plasma data are available at Cassini for this period due to pointing restrictions imposed on the instrument. This renders the identification of the nature of such structures ambiguous, as determinations based on the magnetic field data alone are unreliable. However, a fortuitous alignment of the planets during this encounter allowed us to trace these structures back to those observed previously by the Wind spacecraft near the Earth. Of the phenomena that we are satisfactorily able to trace back to their manifestation at 1 AU, two are identified as being due to interplanetary coronal mass ejections. One event at Cassini is shown to be a merged interaction region, which is formed from the compression of a magnetic cloud by two anomalously fast solar wind streams. The flux-rope structure associated with this magnetic cloud is not as apparent at Cassini and has most likely been compressed and deformed. Confirmation of the validity of the ballistic projections used here is provided by results obtained from a one-dimensional magnetohydrodynamic projection of solar wind parameters measured upstream near the Earth. It is found that when the Earth and Cassini are within a few tens of degrees in heliospheric longitude, the results of this one-dimensional model predict the actual conditions measured at 5 AU to an impressive degree. Finally, the validity of the use of such one-dimensional projections in obtaining quasi-solar wind parameters at the outer planets is discussed.

Moving boundary value problems for the wave equation

We study certain boundary value problems for the one-dimensional wave equation posed in a time-dependent domain. The approach we propose is based on a general transform method for solving boundary value problems for integrable nonlinear PDE in two variables, that has been applied extensively to the study of linear parabolic and elliptic equations. Here we analyse the wave equation as a simple illustrative example to discuss the particular features of this method in the context of linear hyperbolic PDEs, which have not been studied before in this framework.

Boundary value problems for the N-wave interaction equations

We consider boundary value problems for the N-wave interaction equations in one and two space dimensions, posed for x [greater-or-equal, slanted] 0 and x,y [greater-or-equal, slanted] 0, respectively. Following the recent work of Fokas, we develop an inverse scattering formalism to solve these problems by considering the simultaneous spectral analysis of the two ordinary differential equations in the associated Lax pair. The solution of the boundary value problems is obtained through the solution of a local Riemann–Hilbert problem in the one-dimensional case, and a nonlocal Riemann–Hilbert problem in the two-dimensional case.

Characterising the background errors for the boundary-layer capping inversion

The one-dimensional variational assimilation of vertical temperature information in the presence of a boundary-layer capping inversion is studied. For an optimal analysis of the vertical temperature profile, an accurate representation of the background error covariances is essential. The background error covariances are highly flow-dependent due to the variability in the presence, structure and height of the boundary-layer capping inversion. Flow-dependent estimates of the background error covariances are shown by studying the spread in an ensemble of forecasts. A forecast of the temperature profile (used as a background state) may have a significant error in the position of the capping inversion with respect to observations. It is shown that the assimilation of observations may weaken the inversion structure in the analysis if only magnitude errors are accounted for as is the case for traditional data assimilation methods used for operational weather prediction. The positional error is treated explicitly here in a new data assimilation scheme to reduce positional error, in addition to the traditional framework to reduce magnitude error. The distribution of the positional error of the background inversion is estimated for use with the new scheme.

Lagrangian validation of numerical drifter trajectories using drifting buoys: Application to the Agulhas system

The skill of numerical Lagrangian drifter trajectories in three numerical models is assessed by comparing these numerically obtained paths to the trajectories of drifting buoys in the real ocean. The skill assessment is performed using the two-sample Kolmogorov–Smirnov statistical test. To demonstrate the assessment procedure, it is applied to three different models of the Agulhas region. The test can either be performed using crossing positions of one-dimensional sections in order to test model performance in specific locations, or using the total two-dimensional data set of trajectories. The test yields four quantities: a binary decision of model skill, a confidence level which can be used as a measure of goodness-of-fit of the model, a test statistic which can be used to determine the sensitivity of the confidence level, and cumulative distribution functions that aid in the qualitative analysis. The ordering of models by their confidence levels is the same as the ordering based on the qualitative analysis, which suggests that the method is suited for model validation. Only one of the three models, a 1/10° two-way nested regional ocean model, might have skill in the Agulhas region. The other two models, a 1/2° global model and a 1/8° assimilative model, might have skill only on some sections in the region

Orography in a contour dynamics model of large-scale atmospheric flow

The influence of orography on the structure of stationary planetary Rossby waves is studied in the context of a contour dynamics model of the large-scale atmospheric flow. Orography of infinitesimal and finite amplitude is studied using analytical and numerical techniques. Three different types of orography are considered: idealized orography in the form of a global wave, idealized orography in the form of a local table mountain, and the earth's orography. The study confirms the importance of resonances, both in the infinitesimal orography and in the finite orography cases. With finite orography the stationary waves organize themselves into a one-dimensional set of solutions, which due to the resonances, is piecewise connected. It is pointed out that these stationary waves could be relevant for atmospheric regimes.

Stabilization of a helical water chain in a metal-organic host of a trinuclear Schiff base complex

Three heterometallic trinuclear Schiff base complexes, [{GuL(1)(H2O)}(2)Ni(CN)(4)]center dot 4H(2)O (1), [{CuL2(H2O)}(2)Ni(CN)(4)] (2), and [{CuL3(H2O)}(2)Ni(CN)(4)] (3) (HL1 = 7-amino-4-methyl-5-azahept-3-en-2-one, HL2 = 7-methylamino-4-methyl-5-azahept-3-en-2-one, and HL3 = 7-dimethylamino-4-methyl-5-azahept-3-en-2-one), were synthesized. All three complexes were characterized by elemental analysis, IR and UV spectroscopies, and thermal analysis. Two of them (1 and 3) were also characterized by single crystal X-ray crystallography. Complex 1 forms a hydrogen-bonded one-dimensional metal-organic framework that stabilizes a helical water chain into its cavity, but when any of the amine hydrogen atoms of the Schiff base are replaced by methyl groups, as in L 2 and L 3, the water chain, vanishes, showing explicitly the importance of the host-guest H-bonding interactions for the stabilization of a water cluster.

In vivo transformations of dihydroartemisinic acid in Artemisia annua plants

[15-(CH3)-C-13-H-2]-dihydroartemisinic acid (2a) and [15-(CH3)-H-2]-dihydroartemisinic acid (2b) have been fed via the root to intact Artemisia annua plants and their transformations studied in vivo by one-dimensional H-2 NMR spectroscopy and two-dimensional, C-13-H-2 correlation NMR spectroscopy (C-13-(2) H COSY). Labelled dihydroartemisinic acid was transformed into 16 12-carboxy-amorphane and cadinane sesquiterpenes within a few days in the aerial parts of A. annua, although transformations in the root were much slower and more limited. Fifteen of these 16 metabolites have been reported previously as natural products from A. annua. Evidence is presented that the first step in the transformation of dihydroartemisinic acid in vivo is the formation of allylic hydroperoxides by the reaction of molecular oxygen with the Delta(4,5)-double bond in this compound. The origin of all 16 secondary metabolites might then be explained by the known further reactions of such hydroperoxides. The qualitative pattern for the transformations of dihydroartemisinic acid in vivo was essentially unaltered when a comparison was made between plants, which had been kept alive and plants which were allowed to die after feeding of the labelled precursor. This, coupled with the observation that the pattern of transformations of 2 in vivo demonstrated very close parallels with the spontaneous autoxidation chemistry for 2, which we have recently demonstrated in vitro, has lead us to conclude that the main 'metabolic route' for dihydroartemisinic acid in A. annua involves its spontaneous autoxidation and the subsequent spontaneous reactions of allylic hydroperoxides which are derived from 2. There may be no need to invoke the participation of enzymes in any of the later biogenetic steps leading to all 16 of the labelled 11,13-dihydro-amorphane sesquiterpenes which are found in A. annua as natural products. (C) 2003 Elsevier Ltd. All rights reserved.

Facile synthesis of Cu(II) complexes of monocondensed N,N,N donor Schiff base ligands: crystal structure, spectroscopic and magnetic properties

Four new copper(II) complexes, [((CuLN3)-N-1)(2)](ClO4)(2) (1), [(CuL2 N-3)(2)](ClO4)(2) (2), [CuL3(N-3)ClO4)](n) (3) and [CuL4(mu-1,1-N-3)(mu-1,3-N-3)(ClO4)](n) (4) where L-1 = N-1-pyridin-2-yl-methylene-propane-1,3-diamine, L-2 = N-1-(1-pyridin-2-yl-ethylidene)propane-1,3-diamine, L-3 =N-1-(1-pyridin-2-yl-ethylidene)ethane-1,2-diamine and L-4=N-1-(1-pyridin-2-yl-ethylidene)propane-1,2-diamine are four tridentate N,N,N donor Schiff base ligands, have been derived and structurally characterized by X-ray crystallography. Compounds 1 and 2 consist of double basal-apical end-on (EO) azide bridged dinuclear Cu-II complexes with square-pyramidal geometry. In complex 3 the square planar mononuclear [CuL3 (N-3)] units are linked by weakly coordinated perchlorate ions in the axial positions of Cu-II to form a one-dimensional chain. Two such chains are connected by hydrogen bonds involving perchlorate ions and azide groups. Compound 4 consists of 1-D chains in which the Cu-II ions with a square-pyramidal geometry are alternately bridged by single EO and end-to-end (EE) azido ligands, both adopting a basal-apical disposition. Variable temperature (300-2 K) magnetic susceptibility measurements and magnetization measurements at 2 K have been performed. The results reveal that complexes 1 and 2 are antiferromagnetically coupled through azido bridges (J= -12.18 +/- 0.09 and -4.43 +/- 0.1 cm(-1) for 1 and 2, respectively). Complex 3 shows two different magnetic interactions through the two kinds of hydrogen bonds; one is antiferromagnetic (J(1) = - 9.69 +/- 0.03 cm(-1)) and the other is ferromagnetic (J(2) = 1.00 +/- 0.01 cm(-1)). From a magnetic point of view complex 4 is a ferromagnetic dinuclear complex (J= 1.91 +/- 0.01 cm(-1)) coupled through the EO bridge only. The coupling through the EE bridge is practically nil as the N(azido)-Cu-II (axial) distance (2.643 angstrom) is too long. (C) 2006 Elsevier Ltd. All rights reserved.

A case of four-coordinate silver(I) adopting a high energy planar geometry. Experiment and theory

The ligands PhL and MeL are obtained by condensing 2-formylpyridine with benzil dihydrazone and diacetyl dihydrazone, respectively, in 2: 1 molar proportion. With silver( I), PhL yields a double-stranded dinuclear cationic helicate 1 in which the metal is tetrahedral but MeL gives a cationic one-dimensional polymeric complex 2 where silver( I) is distorted square planar and the ligand backbone is nearly planar. In both complexes, metal: ligand ratio is 1: 1. Ab initio calculations on the ligands at the HF/6-31+G* level reveal that while PhL strongly prefers a helical conformation, MeL has a natural inclination to remain in a planar conformation. Density functional theory calculations on model silver( I) complexes show that formation of the linear polymer in the case of MeL is also an important factor in imposing the planar geometry of Ag(I) in 2.

Zinc(II) mediated synthesis of an N-substituted 2-(2-pyridyl)imidazole from a 1,2-diketone and 2-(aminomethyl)pyridine and its ligational behaviour

Reaction of anhydrous ZnCl2 with the 1:2 condensate (L) of benzil and 2-(aminomethyl)pyridine in methanol gives monomeric ZnL'Cl-2 (1) where L' is 2-[(4,5-diphenyl-2-pyridin-2-yl-1H-imidazol-1-yl)-methyl]pyridine. In the X-ray crystal structure, 1 is found to contain tetrahedral zinc with an N2Cl2 coordination sphere and the N-substituent methylpyridine fragment hanging as a free arm. A tentative mechanism is proposed for the zinc mediated conversion of L-->L'. Demetallation of 1 by the action of aqueous NaOH yields L' in the free state. When L' is reacted with Zn(ClO4)(2).6H(2)O in a 1:2 molar proportion, [Zn(L')(2)](n)(ClO4)(2n).(H2O)(n/2).(CH2Cl2)(n/2) (2) is obtained. The zinc atom in 2, as revealed by X-ray crystallography, has a trigonal bipyramidal N-5 coordination sphere. There are two independent ligands in the asymmetric unit of 2. One of them bonds only to one zinc atom in a bidentate mode with the N-substituent methylpyridine hanging free while the other ligand binds to two different zinc atoms in a tridentate fashion, employing the N-substituent methylpyridine nitrogen atom to form the polymeric one-dimensional chain cation.

A novel cation induced polymeric chain in Na-8[{Cu(gly)(2)}(2)]{H-2(H2W12O42)}] center dot 24H(2)O: hydrothermal synthesis, spectroscopic characterization and X-ray structure analysis

A novel bis(glycinato) copper(II) paradodecatungstate Na-8[{Cu(gly)(2)}(2)]-{H-2(H2W12O42)}] center dot 24H(2)O (1) has been synthesized under hydrothermal conditions. The crystal structure of 1 reveals an infinite one-dimensional chain along the [100] direction and is built from paradodecatungstate (H2W12O42)(10-) clusters joined through [Cu(gly)(2)] moieties. Parallel chains are interlinked by NaO6 octahedra to generate a two-dimensional network.