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.

Consistent approximate models of the global atmosphere: shallow, deep, hydrostatic, quasi-hydrostatic and non-hydrostatic

We study global atmosphere models that are at least as accurate as the hydrostatic primitive equations (HPEs), reviewing known results and reporting some new ones. The HPEs make spherical geopotential and shallow atmosphere approximations in addition to the hydrostatic approximation. As is well known, a consistent application of the shallow atmosphere approximation requires omission of those Coriolis terms that vary as the cosine of latitude and of certain other terms in the components of the momentum equation. An approximate model is here regarded as consistent if it formally preserves conservation principles for axial angular momentum, energy and potential vorticity, and (following R. Müller) if its momentum component equations have Lagrange's form. Within these criteria, four consistent approximate global models, including the HPEs themselves, are identified in a height-coordinate framework. The four models, each of which includes the spherical geopotential approximation, correspond to whether the shallow atmosphere and hydrostatic (or quasi-hydrostatic) approximations are individually made or not made. Restrictions on representing the spatial variation of apparent gravity occur. Solution methods and the situation in a pressure-coordinate framework are discussed. © Crown copyright 2005.

The influence of physical processes on extratropical singular vectors

An investigation is made of the impact of a full linearized physical (moist) parameterization package on extratropical singular vectors (SVs) using the ECMWF integrated forecasting system (IFS). Comparison is made for one particular period with a dry physical package including only vertical diffusion and surface drag. The crucial extra ingredient in the full package is found to be the large-scale latent heat release. Consistent with basic theory, its inclusion results in a shift to smaller horizontal scales and enhanced growth for the SVs. Whereas, for the dry SVs, T42 resolution is sufficient, the moist SVs require T63 to resolve their structure and growth. A 24-h optimization time appears to be appropriate for the moist SVs because of the larger growth of moist SVs compared with dry SVs. Like dry SVs, moist SVs tend to occur in regions of high baroclinicity, but their location is also influenced by the availability of moisture. The most rapidly growing SVs appear to enhance or reduce large-scale rain in regions ahead of major cold fronts. The enhancement occurs in and ahead of a cyclonic perturbation and the reduction in and ahead of an anticyclonic perturbation. Most of the moist SVs for this situation are slightly modified versions of the dry SVs. However, some occur in new locations and have particularly confined structures. The most rapidly growing SV is shown to exhibit quite linear behavior in the nonlinear model as it grows from 0.5 to 12 hPa in 1 day. For 5 times this amplitude the structure is similar but the growth is about half as the perturbation damps a potential vorticity (PV) trough or produces a cutoff, depending on its sign.

The structure of crystallisable copolymers of L-lactide, epsilon-caprolactone and glycolide

We apply a new X-ray scattering approach to the study of melt-spun filaments of tri-block and random terpolymers prepared from lactide, caprolactone and glycolide. Both terpolymers contain random sequences, in both cases the overall fraction of lactide units is similar to 0.7 and C-13 and H-1 NMR shows the lactide sequence length to be similar to 9-10. A novel representation of the X-ray fibre pattern as series of spherical harmonic functions considerably facilitates the comparison of the scattering from the minority crystalline phase with hot drawn fibres prepared from the poly(L-lactide) homopolymer. Although the fibres exhibit rather disordered structures we show that the crystal structure is equivalent to that displayed by poly(L-lactide) for both the block and random terpolymers. There are variations in the development of a two-phase structure which reflect the differences in the chain architectures. There is evidence that the random terpolymer includes non-lactide units in to the crystal interfaces to achieve a well defined two-phase structure. (c) 2005 Published by Elsevier Ltd.

Rotational spectra induced by vibrations

Mizushima and Venkateswarlu showed in 1953 that certain molecules have the property that excited vibrational states may possess rotational spectra even when the rotational spectrum of the ground vibrational state is forbidden by symmetry. We call such a spectrum a vibrationally induced rotational spectrum, and have made a systematic examination of the point groups which permit such behaviour. We also give formulae for the approximate line frequencies and intensities in these spectra, and discuss some of the problems involved in observing them. The spectra can only arise from degenerate vibrational states, and are of three possible types: i) symmetric top perpendicular spectra, shown by molecules belonging to the point groups Dnh, Dn and Cnh, where n is odd; (ii) symmetric top parallel spectra, shown by molecules belonging to Dnd and S2n, where n is even; and (iii) spherical top spectra, shown by molecules belonging to T or Td. Excited vibrational states of polar molecules of point groups Cnv or Cn, where n is odd, may also possess vibrationally induced perpendicular components of type (i), in addition to their ordinary parallel spectra. In addition to the above limitations on the point groups there are, in general, limitations on the symmetry species of the degenerate vibrational states.

Wave functions for the methane molecule

If the potential field due to the nuclei in the methane molecule is expanded in terms of a set of spherical harmonics about the carbon nucleus, only the terms involving s, f, and higher harmonic functions differ from zero in the equilibrium configuration. Wave functions have been calculated for the equilibrium configuration, first including only the spherically symmetric s term in the potential, and secondly including both the s and the f terms. In the first calculation the complete Hartree-Fock S.C.F. wave functions were determined; in the second calculation a variation method was used to determine the best form of the wave function involving f harmonics. The resulting wave functions and electron density functions are presented and discussed

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.

Monthly average temperature modelling

This research is associated with the goal of the horticultural sector of the Colombian southwest, which is to obtain climatic information, specifically, to predict the monthly average temperature in sites where it has not been measured. The data correspond to monthly average temperature, and were recorded in meteorological stations at Valle del Cauca, Colombia, South America. Two components are identified in the data of this research: (1) a component due to the temporal aspects, determined by characteristics of the time series, distribution of the monthly average temperature through the months and the temporal phenomena, which increased (El Nino) and decreased (La Nina) the temperature values, and (2) a component due to the sites, which is determined for the clear differentiation of two populations, the valley and the mountains, which are associated with the pattern of monthly average temperature and with the altitude. Finally, due to the closeness between meteorological stations it is possible to find spatial correlation between data from nearby sites. In the first instance a random coefficient model without spatial covariance structure in the errors is obtained by month and geographical location (mountains and valley, respectively). Models for wet periods in mountains show a normal distribution in the errors; models for the valley and dry periods in mountains do not exhibit a normal pattern in the errors. In models of mountains and wet periods, omni-directional weighted variograms for residuals show spatial continuity. The random coefficient model without spatial covariance structure in the errors and the random coefficient model with spatial covariance structure in the errors are capturing the influence of the El Nino and La Nina phenomena, which indicates that the inclusion of the random part in the model is appropriate. The altitude variable contributes significantly in the models for mountains. In general, the cross-validation process indicates that the random coefficient model with spatial spherical and the random coefficient model with spatial Gaussian are the best models for the wet periods in mountains, and the worst model is the model used by the Colombian Institute for Meteorology, Hydrology and Environmental Studies (IDEAM) to predict temperature.

The extrapolation problem and how population modeling can help

We argue that population modeling can add value to ecological risk assessment by reducing uncertainty when extrapolating from ecotoxicological observations to relevant ecological effects. We review other methods of extrapolation, ranging from application factors to species sensitivity distributions to suborganismal (biomarker and "-omics'') responses to quantitative structure activity relationships and model ecosystems, drawing attention to the limitations of each. We suggest a simple classification of population models and critically examine each model in an extrapolation context. We conclude that population models have the potential for adding value to ecological risk assessment by incorporating better understanding of the links between individual responses and population size and structure and by incorporating greater levels of ecological complexity. A number of issues, however, need to be addressed before such models are likely to become more widely used. In a science context, these involve challenges in parameterization, questions about appropriate levels of complexity, issues concerning how specific or general the models need to be, and the extent to which interactions through competition and trophic relationships can be easily incorporated.

Electrochemical and sonoelectrochemical monitoring of indigo reduction by glucose

The reduction of indigo (dispersed in water) to leuco-indigo (dissolved in water) is an important industrial process and investigated here for the case of glucose as an environmentally benign reducing agent. In order to quantitatively follow the formation of leuco-indigo two approaches based on (i) rotating disk voltammetry and (ii) sonovoltammetry are developed. Leuco-indigo, once formed in alkaline solution, is readily monitored at a glassy carbon electrode in the mass transport limit employing hydrodynamic voltammetry. The presence of power ultrasound further improves the leuco-indigo determination due to additional agitation and homogenization effects. While inactive at room temperature, glucose readily reduces indigo in alkaline media at 65 degrees C. In the presence of excess glucose, a surface dissolution kinetics limited process is proposed following the rate law d eta(leuco-indigo)/dt = k x c(OH-) x S-indigo where eta(leuco-indigo) is the amount of leuco-indigo formed, k = 4.1 x 10(-9) m s(-1) (at 65 degrees C, assuming spherical particles of I gm diameter) is the heterogeneous dissolution rate constant,c(OH-) is the concentration of hydroxide, and Sindigo is the reactive surface area. The activation energy for this process in aqueous 0.2 M NaOH is E-A = 64 U mol(-1) consistent with a considerable temperature effects. The redox mediator 1,8-dihydroxyanthraquinone is shown to significantly enhance the reaction rate by catalysing the electron transfer between glucose and solid indigo particles. (c) 2006 Elsevier Ltd. All fights reserved.

Malacosporean-like spores in urine of rainbow trout react with antibody and DNA probes to Tetracapsuloides bryosalmonae

Tetracapsuloides bryosalmonae is the myxozoan parasite causing proliferative kidney disease (PKD) of salmonid fishes in Europe and North America. The complete life cycle of the parasite remains unknown despite recent discoveries that the stages infectious for fish develop in freshwater bryozoans. During the course of examinations of the urine of rainbow trout (Oncorhynchus mykiss) with or recovering from PKD we identified spores with features similar to those of T. bryosalmonae found in the bryozoan host. Spores found in the urine were subspherical, with a width of 16 mum and height of 14 mum, and possessed two soft valves surrounding two spherical polar capsules (2 mum in diameter) and a single sporoplasm. The absence of hardened valves is a distinguishing characteristic of the newly established class Malacosporea that includes T. bryosalmonae as found in the bryozoan host. The parasite in the urine of rainbow trout possessed only two polar capsules and two valve cells compared to the four polar capsules and four valves observed in the spherical spores of 19 mum in diameter from T. bryosalmonae from the bryozoan host. Despite morphological differences, a relationship between the spores in the urine of rainbow trout and T. bryosalmonae was demonstrated by binding of monoclonal and polyclonal antibodies and DNA probes specific to T. bryosalmonae.

Insights into the effects on metal binding of the systematic substitution of five key glutamate ligands in the ferritin of Escherichia coli

Ferritins are nearly ubiquitous iron storage proteins playing a fundamental role in iron metabolism. They are composed of 24 subunits forming a spherical protein shell encompassing a central iron storage cavity. The iron storage mechanism involves the initial binding and subsequent O-2-dependent oxidation of two Fe2+ ions located at sites A and B within the highly conserved dinuclear "ferroxidase center" in individual subunits. Unlike animal ferritins and the heme-containing bacterioferritins, the Escherichia coli ferritin possesses an additional iron-binding site (site C) located on the inner surface of the protein shell close to the ferroxidase center. We report the structures of five E. coli ferritin variants and their Fe3+ and Zn2+ (a redox-stable alternative for Fe2+) derivatives. Single carboxyl ligand replacements in sites A, B, and C gave unique effects on metal binding, which explain the observed changes in Fe2+ oxidation rates. Binding of Fe2+ at both A and B sites is clearly essential for rapid Fe2+ oxidation, and the linking of Fe-B(2+) to Fe-C(2+) enables the oxidation of three Fe2+ ions. The transient binding of Fe2+ at one of three newly observed Zn2+ sites may allow the oxidation of four Fe2+ by one dioxygen molecule.

LDA or GGA? A combined experimental inelastic neutron scattering and ab initio lattice dynamics study of alkali metal hydrides

In a previous work, we carried out inelastic neutron scattering (INS) spectroscopy experiments and preliminary first principles calculations on alkali metal hydrides. The complete series of alkali metal hydrides, LiH, NaH, KH, RbH and CsH was measured in the high-resolution TOSCA INS spectrometer at ISIS. Here, we present the results of ab initio electronic structure calculations of the properties of the alkali metal hydrides using both the local density approximation (LDA) and the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) parameterization. Properties calculated were lattice parameters, bulk moduli, dielectric constants, effective charges, electronic densities and inelastic neutron scattering (INS) spectra. We took advantage of the currently available computer power to use full lattice dynamics theory to calculate thermodynamic properties for these materials. For the alkali metal hydrides (LiH, NaH, KH, RbH and CsH) using lattice dynamics, we found that the INS spectra calculated using LDA agreed better with the experimental data than the spectra calculated using GGA. Both zero-point effects and thermal contributions to free energies had an important effect on INS and several thermodynamic properties.

Fractionated crystallization and fractionated melting of confined PEO microdomains in PB-b-PEO and PE-b-PEO diblock copolymers

The confined crystallization of poly(ethylene oxide) (PEO) in predominantly spherical microdomains formed by several diblock copolymers was studied and compared. Two polybutadiene-b-poly(ethylene oxide) diblock copolymers were prepared by sequential anionic polymerization (with approximately 90 and 80 wt % polybutadiene (PB)). These were compared to equivalent samples after catalytic hydrogenation that produced double crystalline polyethylene-b-poly(ethylene oxide) diblock copolymers. Both systems are segregated into microdomains as indicated by small-angle X-ray scattering (SAXS) experiments performed in the melt and at lower temperatures. However, the PB-b-PEO systems exhibited a higher degree of order in the melt. A predominantly spherical morphology of PEO in a PB or a PE matrix was observed by both SAXS and transmission electron microscopy, although a possibly mixed morphology (spheres and cylinders) was formed when the PEO composition was close to the cylinder-sphere domain transitional composition as indicated by SAXS. Differential scanning calorimetry experiments showed that a fractionated crystallization process for the PEO occurred in all samples, indicating that the PE cannot nucleate PEO in these diblock copolymers. A novel result was the observation of a subsequent fractionated melting that reflected the crystallization process. Sequential isothermal crystallization experiments allowed us to thermally separate at least three different crystallization and melting peaks for the PEO microdomains. The lowest melting point fraction was the most important in terms of quantity and corresponded to the crystallization of isolated PEO spheres (or cylinders) that were either superficially or homogeneously nucleated. This was confirmed by Avrami index values of approximately 1. The isothermal crystallization results indicate that the PE matrix restricts the crystallization of the covalently bonded PEO to a higher degree compared to PB.

Effect of ethanol on the micellization and gelation of Pluronic P123

In certain applications copolymer P123 (E21P67E21) is dissolved in water-ethanol mixtures, initially to form micellar solutions and eventually to gel. For P123 in 10, 20, and 30 wt % aqueous ethanol we used dynamic light scattering from dilute solutions to confirm micellization, oscillatory rheometry, and visual observation of mobility (tube inversion) to determine gel formation in concentrated solutions and small-angle X-ray scattering (SAXS) to determine gel structure. Except for solutions in 30 wt % aqueous ethanol, a clear-turbid transition was encountered on heating dilute and concentrated micellar solutions alike, and as for solutions in water alone (Chaibundit et al. Langmuir 2007, 23, 9229) this could be ascribed to formation of wormlike micelles. Dense clouding, typical of phase separation, was observed at higher temperatures. Regions of isotropic and birefringent gel were defined for concentrated solutions and shown (by SAXS) to have Cubic (fcc and hcp) and hexagonal structures, consistent with packed spherical and elongated micelles, respectively. The cubic gels (0, 10, and 20 wt % ethanol) were clear, while the hex gels were either turbid (0 and 10 wt % ethanol), turbid enclosing a clear region (20 wt % ethanol), or entirely clear (30 wt % ethanol). The SAXS profile was unchanged between turbid and clear regions of the 20 wt % ethanol gel. Temperature scans of dynamic moduli showed (as expected) a clear distinction between high-modulus cubic gels (G'(max) approximate to 20-30 kPa) and lower modulus hex gels (G'(max) < 10 kPa).