An adaptive finite element water column model using the Mellor-Yamada level 2.5 turbulence closure scheme

A one-dimensional water column model using the Mellor and Yamada level 2.5 parameterization of vertical turbulent fluxes is presented. The model equations are discretized with a mixed finite element scheme. Details of the finite element discrete equations are given and adaptive mesh refinement strategies are presented. The refinement criterion is an "a posteriori" error estimator based on stratification, shear and distance to surface. The model performances are assessed by studying the stress driven penetration of a turbulent layer into a stratified fluid. This example illustrates the ability of the presented model to follow some internal structures of the flow and paves the way for truly generalized vertical coordinates. (c) 2005 Elsevier Ltd. All rights reserved.

Dynamical evolution of the 2003 southern hemisphere stratospheric winter using Envisat trace-gas observations

Data from the MIPAS instrument on Envisat, supplemented by meteorological analyses from ECMWF and the Met Office, are used to study the meteorological and trace-gas evolution of the stratosphere in the southern hemisphere during winter and spring 2003. A pole-centred approach is used to interpret the data in the physically meaningful context of the evolving stratospheric polar vortex. The following salient dynamical and transport features are documented and analysed: the merger of anticyclones in the stratosphere; the development of an intense, quasi-stationary anticyclone in spring; the associated top-down breakdown of the polar vortex; the systematic descent of air into the polar vortex; and the formation of a three-dimensional structure of a tracer filament on a planetary scale. The paper confirms and extends existing paradigms of the southern hemisphere vortex evolution. The quality of the MIPAS observations is seen to be generally good. though the water vapour retrievals are unrealistic above 10 hPa in the high-latitude winter.

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.

Atmospheric factors governing banded orographic convection

Atmospheric factors Governing Banded Orographic Convection The three-dimensional structure of shallow orographic convection is investigated through simulations performed with a cloud-resolving numerical model. In moist flows that overcome a given topographic barrier to form statically unstable cap clouds, the organization of the convection depends on both the atmospheric structure and the mechanism by which the convection is initiated. Convection initiated by background thermal fluctuations embedded in the flow over a smooth mountain (without any small-scale topographic features) tends to be cellular and disorganized except that shear-parallel bands may form in flows with strong unidirectional vertical shear. The development of well-organized bands is favored when there is weak static instability inside the cloud and when the dry air surrounding the cloud is strongly stable. These bands move with the flow and distribute their cumulative precipitation evenly over the mountain upslope. Similar shear-parallel bands also develop in flows where convection is initiated by small-scale topographic noise superimposed onto the main mountain profile, but in this case stronger circulations are also triggered that create stationary rainbands parallel to the low-level flow. This second dominant mode, which is less sensitive to the atmospheric structure and the strength of forcing, is triggered by lee waves that form over small-scale topographic bumps near the upstream edge of the main orographic cloud. Due to their stationarity, these flow-parallel bands can produce locally heavy precipitation amounts.

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.

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

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.

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.

Solvent-induced dynamic single-crystal-to-single-crystal transformation of a synthetic peptide-based cyclic compound

Solvent induced single-crystal-to-single-crystal transformation has been demonstrated indicating the dynamic behavior of one dimensional arrays obtained from a self-assembled new synthetic cyclic peptide.

The synthesis and characterisation of four new antimony sulphides incorporating transition-metal complexes

Four new antimony sulphides, [T(dien)(2)]Sb6S10 center dot xH(2)O [T = Ni (1), Co (2) x approximate to 0.45], [Co(en)(3)]SbsSI(3) (3) and [Ni(en)(3)]Sb12S19 (4), have been synthesised under solvothermal conditions. In compounds (1) - (3), Sb12S228- secondary building units are connected to form layered structures. In (1) and (2), Sb-6 S-2- layers containing Sb16S16 heterorings are separated by [T(dien]2](2+) cations, whilst in (3), Sb8 S2- layers 10 13 contain [Co(en)3]2+ cations within large Sb22S22 pores. Compound (4) adopts a three-dimensional structure in which [Ni(en)3 12 cations lie within ca. 5 A wide channels. (c) 2007 Elsevier Ltd. All rights reserved.

Structural studies on acridine derivatives binding to telomeric DNA

Acridine derivatives can inhibit a variety of nuclear enzymes by binding or intercalating to DNA. This class of compounds is of great interest in the development of novel anticancer agents. Despite the availability of crystallographic data for some of the compounds complexed with DNA, uncertainties remain about the mechanisms of action, binding preferences and biological targets. To investigate the intercalation of several acridine derivatives, a variety of techniques are being employed. Single-crystal X-ray diffraction is being used to determine the high resolution three-dimensional structure of short sequences of quadruplex telomeric DNA with bound drug. This will be compared to the effect of drug binding to long segments of double-stranded DNA using fibre diffraction, with neutron diffraction studies planned to analyse the hydrogen bonding patterns of the DNA-drug complexes. Small-angle neutron scattering (SANS) will also be applied to study drug binding to both short and long sequences of quadruplex and double-stranded DNA in solution. Initial SANS measurements of the telomeric repeat d(TGGGGT) imply that this hexamer is present as a quadruplex. (c) 2006 Elsevier B.V. All rights reserved.