Calculating the ocean's mean dynamic topography from a mean sea surface and a geoid

In principle the global mean geostrophic surface circulation of the ocean can be diagnosed by subtracting a geoid from a mean sea surface (MSS). However, because the resulting mean dynamic topography (MDT) is approximately two orders of magnitude smaller than either of the constituent surfaces, and because the geoid is most naturally expressed as a spectral model while the MSS is a gridded product, in practice complications arise. Two algorithms for combining MSS and satellite-derived geoid data to determine the ocean’s mean dynamic topography (MDT) are considered in this paper: a pointwise approach, whereby the gridded geoid height field is subtracted from the gridded MSS; and a spectral approach, whereby the spherical harmonic coefficients of the geoid are subtracted from an equivalent set of coefficients representing the MSS, from which the gridded MDT is then obtained. The essential difference is that with the latter approach the MSS is truncated, a form of filtering, just as with the geoid. This ensures that errors of omission resulting from the truncation of the geoid, which are small in comparison to the geoid but large in comparison to the MDT, are matched, and therefore negated, by similar errors of omission in the MSS. The MDTs produced by both methods require additional filtering. However, the spectral MDT requires less filtering to remove noise, and therefore it retains more oceanographic information than its pointwise equivalent. The spectral method also results in a more realistic MDT at coastlines. 1. Introduction An important challenge in oceanography is the accurate determination of the ocean’s time-mean dynamic topography (MDT). If this can be achieved with sufficient accuracy for combination with the timedependent component of the dynamic topography, obtainable from altimetric data, then the resulting sum (i.e., the absolute dynamic topography) will give an accurate picture of surface geostrophic currents and ocean transports.

Performance of the OPA/ARPEGE-T21 global ocean-atmosphere coupled model

The climatology of the OPA/ARPEGE-T21 coupled general circulation model (GCM) is presented. The atmosphere GCM has a T21 spectral truncation and the ocean GCM has a 2°×1.5° average resolution. A 50-year climatic simulation is performed using the OASIS coupler, without flux correction techniques. The mean state and seasonal cycle for the last 10 years of the experiment are described and compared to the corresponding uncoupled experiments and to climatology when available. The model reasonably simulates most of the basic features of the observed climate. Energy budgets and transports in the coupled system, of importance for climate studies, are assessed and prove to be within available estimates. After an adjustment phase of a few years, the model stabilizes around a mean state where the tropics are warm and resemble a permanent ENSO, the Southern Ocean warms and almost no sea-ice is left in the Southern Hemisphere. The atmospheric circulation becomes more zonal and symmetric with respect to the equator. Once those systematic errors are established, the model shows little secular drift, the small remaining trends being mainly associated to horizontal physics in the ocean GCM. The stability of the model is shown to be related to qualities already present in the uncoupled GCMs used, namely a balanced radiation budget at the top-of-the-atmosphere and a tight ocean thermocline.

Equivalence of truncated count mixture distributions and mixtures of truncated count distributions

This article is about modeling count data with zero truncation. A parametric count density family is considered. The truncated mixture of densities from this family is different from the mixture of truncated densities from the same family. Whereas the former model is more natural to formulate and to interpret, the latter model is theoretically easier to treat. It is shown that for any mixing distribution leading to a truncated mixture, a (usually different) mixing distribution can be found so. that the associated mixture of truncated densities equals the truncated mixture, and vice versa. This implies that the likelihood surfaces for both situations agree, and in this sense both models are equivalent. Zero-truncated count data models are used frequently in the capture-recapture setting to estimate population size, and it can be shown that the two Horvitz-Thompson estimators, associated with the two models, agree. In particular, it is possible to achieve strong results for mixtures of truncated Poisson densities, including reliable, global construction of the unique NPMLE (nonparametric maximum likelihood estimator) of the mixing distribution, implying a unique estimator for the population size. The benefit of these results lies in the fact that it is valid to work with the mixture of truncated count densities, which is less appealing for the practitioner but theoretically easier. Mixtures of truncated count densities form a convex linear model, for which a developed theory exists, including global maximum likelihood theory as well as algorithmic approaches. Once the problem has been solved in this class, it might readily be transformed back to the original problem by means of an explicitly given mapping. Applications of these ideas are given, particularly in the case of the truncated Poisson family.

Peptide inhibitors of G protein-coupled receptor kinases

G protein-coupled receptor kinases (GRKs) are regulatory enzymes involved in the modulation of seven-transmembrane-helix receptors. In order to develop specific inhibitors for these kinases, we synthesized and investigated peptide inhibitors derived from the sequence of the first intracellular loop of the beta(2)-adrenergic receptor. Introduction of changes in the sequence and truncation of N- and C-terminal amino acids increased the inhibitory potency by a factor of 40. These inhibitors not only inhibited the prototypical GRK2 but also GRK3 and GRK5. In contrast there was no inhibition of protein kinase C and protein kinase A even at the highest concentration tested. The peptide with the sequence AKFERLQTVTNYFITSE inhibited GRK2 with an IC50 of 0.6 mu M, GRK3 with 2.6 mu M and GRK5 with 1.6 mu M. The peptide inhibitors were non-competitive for receptor and ATP. These findings demonstrate that specific peptides can inhibit GRKs in the submicromolar range and suggest that a further decrease in size is possible without losing the inhibitory potency. (c) 2005 Published by Elsevier Inc.

Parallel Monte Carlo approach for integration of the rendering equation

This paper is addressed to the numerical solving of the rendering equation in realistic image creation. The rendering equation is integral equation describing the light propagation in a scene accordingly to a given illumination model. The used illumination model determines the kernel of the equation under consideration. Nowadays, widely used are the Monte Carlo methods for solving the rendering equation in order to create photorealistic images. In this work we consider the Monte Carlo solving of the rendering equation in the context of the parallel sampling scheme for hemisphere. Our aim is to apply this sampling scheme to stratified Monte Carlo integration method for parallel solving of the rendering equation. The domain for integration of the rendering equation is a hemisphere. We divide the hemispherical domain into a number of equal sub-domains of orthogonal spherical triangles. This domain partitioning allows to solve the rendering equation in parallel. It is known that the Neumann series represent the solution of the integral equation as a infinity sum of integrals. We approximate this sum with a desired truncation error (systematic error) receiving the fixed number of iteration. Then the rendering equation is solved iteratively using Monte Carlo approach. At each iteration we solve multi-dimensional integrals using uniform hemisphere partitioning scheme. An estimate of the rate of convergence is obtained using the stratified Monte Carlo method. This domain partitioning allows easy parallel realization and leads to convergence improvement of the Monte Carlo method. The high performance and Grid computing of the corresponding Monte Carlo scheme are discussed.

The advection of high-resolution tracers by low-resolution winds

The usefulness of any simulation of atmospheric tracers using low-resolution winds relies on both the dominance of large spatial scales in the strain and time dependence that results in a cascade in tracer scales. Here, a quantitative study on the accuracy of such tracer studies is made using the contour advection technique. It is shown that, although contour stretching rates are very insensitive to the spatial truncation of the wind field, the displacement errors in filament position are sensitive. A knowledge of displacement characteristics is essential if Lagrangian simulations are to be used for the inference of airmass origin. A quantitative lower estimate is obtained for the tracer scale factor (TSF): the ratio of the smallest resolved scale in the advecting wind field to the smallest “trustworthy” scale in the tracer field. For a baroclinic wave life cycle the TSF = 6.1 ± 0.3 while for the Northern Hemisphere wintertime lower stratosphere the TSF = 5.5 ± 0.5, when using the most stringent definition of the trustworthy scale. The similarity in the TSF for the two flows is striking and an explanation is discussed in terms of the activity of potential vorticity (PV) filaments. Uncertainty in contour initialization is investigated for the stratospheric case. The effect of smoothing initial contours is to introduce a spinup time, after which wind field truncation errors take over from initialization errors (2–3 days). It is also shown that false detail from the proliferation of finescale filaments limits the useful lifetime of such contour advection simulations to 3σ−1 days, where σ is the filament thinning rate, unless filaments narrower than the trustworthy scale are removed by contour surgery. In addition, PV analysis error and diabatic effects are so strong that only PV filaments wider than 50 km are at all believable, even for very high-resolution winds. The minimum wind field resolution required to accurately simulate filaments down to the erosion scale in the stratosphere (given an initial contour) is estimated and the implications for the modeling of atmospheric chemistry are briefly discussed.

Conserved hydrophobic clusters on the surface of the Caf1A usher C-terminal domain are important for F1 antigen assembly

The outer membrane usher protein Caf1A of the plague pathogen Yersinia pestis is responsible for the assembly of a major surface antigen, the F1 capsule. The F1 capsule is mainly formed by thin linear polymers of Caf1 (capsular antigen fraction 1) protein subunits. The Caf1A usher promotes polymerization of subunits and secretion of growing polymers to the cell surface. The usher monomer (811 aa, 90.5 kDa) consists of a large transmembrane β-barrel that forms a secretion channel and three soluble domains. The periplasmic N-terminal domain binds chaperone-subunit complexes supplying new subunits for the growing fiber. The middle domain, which is structurally similar to Caf1 and other fimbrial subunits, serves as a plug that regulates the permeability of the usher. Here we describe the identification, characterization, and crystal structure of the Caf1A usher C-terminal domain (Caf1A(C)). Caf1A(C) is shown to be a periplasmic domain with a seven-stranded β-barrel fold. Analysis of C-terminal truncation mutants of Caf1A demonstrated that the presence of Caf1A(C) is crucial for the function of the usher in vivo, but that it is not required for the initial binding of chaperone-subunit complexes to the usher. Two clusters of conserved hydrophobic residues on the surface of Caf1A(C) were found to be essential for the efficient assembly of surface polymers. These clusters are conserved between the FGL family and the FGS family of chaperone-usher systems.

Blind channel estimation for multi-rate multicarrier DS/CDMA communications

Multi-rate multicarrier DS/CDMA is a potentially attractive multiple access method for future wireless communications networks that must support multimedia, and thus multi-rate, traffic. Several receiver structures exist for single-rate multicarrier systems, but little has been reported on multi-rate multicarrier systems. Considering that high-performance detection such as coherent demodulation needs the explicit knowledge of the channel, based on the finite-length chip waveform truncation, this paper proposes a subspace-based scheme for timing and channel estimation in multi-rate multicarrier DS/CDMA systems, which is applicable to both multicode and variable spreading factor systems. The performance of the proposed scheme for these two multi-rate systems is validated via numerical simulations. The effects of the finite-length chip waveform truncation on the performance of the proposed scheme is also analyzed theoretically.

Symmetry breaking, mixing, instability, and low frequency variability in a minimal Lorenz-like system

Starting from the classical Saltzman two-dimensional convection equations, we derive via a severe spectral truncation a minimal 10 ODE system which includes the thermal effect of viscous dissipation. Neglecting this process leads to a dynamical system which includes a decoupled generalized Lorenz system. The consideration of this process breaks an important symmetry and couples the dynamics of fast and slow variables, with the ensuing modifications to the structural properties of the attractor and of the spectral features. When the relevant nondimensional number (Eckert number Ec) is different from zero, an additional time scale of O(Ec−1) is introduced in the system, as shown with standard multiscale analysis and made clear by several numerical evidences. Moreover, the system is ergodic and hyperbolic, the slow variables feature long-term memory with 1/f3/2 power spectra, and the fast variables feature amplitude modulation. Increasing the strength of the thermal-viscous feedback has a stabilizing effect, as both the metric entropy and the Kaplan-Yorke attractor dimension decrease monotonically with Ec. The analyzed system features very rich dynamics: it overcomes some of the limitations of the Lorenz system and might have prototypical value in relevant processes in complex systems dynamics, such as the interaction between slow and fast variables, the presence of long-term memory, and the associated extreme value statistics. This analysis shows how neglecting the coupling of slow and fast variables only on the basis of scale analysis can be catastrophic. In fact, this leads to spurious invariances that affect essential dynamical properties (ergodicity, hyperbolicity) and that cause the model losing ability in describing intrinsically multiscale processes.

Evolutionary dynamics of Staphylococcus aureus during progression from carriage to disease

Whole-genome sequencing offers new insights into the evolution of bacterial pathogens and the etiology of bacterial disease. Staph- ylococcus aureus is a major cause of bacteria-associated mortality and invasive disease and is carried asymptomatically by 27% of adults. Eighty percent of bacteremias match the carried strain. How- ever, the role of evolutionary change in the pathogen during the progression from carriage to disease is incompletely understood. Here we use high-throughput genome sequencing to discover the genetic changes that accompany the transition from nasal carriage to fatal bloodstream infection in an individual colonized with meth- icillin-sensitive S. aureus. We found a single, cohesive population exhibiting a repertoire of 30 single-nucleotide polymorphisms and four insertion/deletion variants. Mutations accumulated at a steady rate over a 13-mo period, except for a cluster of mutations preceding the transition to disease. Although bloodstream bacteria differed by just eight mutations from the original nasally carried bacteria, half of those mutations caused truncation of proteins, including a prema- ture stop codon in an AraC-family transcriptional regulator that has been implicated in pathogenicity. Comparison with evolution in two asymptomatic carriers supported the conclusion that clusters of pro- tein-truncating mutations are highly unusual. Our results demon- strate that bacterial diversity in vivo is limited but nonetheless detectable by whole-genome sequencing, enabling the study of evolutionary dynamics within the host. Regulatory or structural changes that occur during carriage may be functionally important for pathogenesis; therefore identifying those changes is a crucial step in understanding the biological causes of invasive bacterial disease.

Efficient calculation of the green function for acoustic propagation above a homogeneous impedance plane

This paper is concerned with the problem of propagation from a monofrequency coherent line source above a plane of homogeneous surface impedance. The solution of this problem occurs in the kernel of certain boundary integral equation formulations of acoustic propagation above an impedance boundary, and the discussion of the paper is motivated by this application. The paper starts by deriving representations, as Laplace-type integrals, of the solution and its first partial derivatives. The evaluation of these integral representations by Gauss-Laguerre quadrature is discussed, and theoretical bounds on the truncation error are obtained. Specific approximations are proposed which are shown to be accurate except in the very near field, for all angles of incidence and a wide range of values of surface impedance. The paper finishes with derivations of partial results and analogous Laplace-type integral representations for the case of a point source.

Crooks' fluctuation theorem for a process on a two dimensional fluid field

We investigate the behavior of a two-dimensional inviscid and incompressible flow when pushed out of dynamical equilibrium. We use the two-dimensional vorticity equation with spectral truncation on a rectangular domain. For a sufficiently large number of degrees of freedom, the equilibrium statistics of the flow can be described through a canonical ensemble with two conserved quantities, energy and enstrophy. To perturb the system out of equilibrium, we change the shape of the domain according to a protocol, which changes the kinetic energy but leaves the enstrophy constant. We interpret this as doing work to the system. Evolving along a forward and its corresponding backward process, we find numerical evidence that the distributions of the work performed satisfy the Crooks relation. We confirm our results by proving the Crooks relation for this system rigorously.

Severe acute respiratory syndrome coronavirus nonstructural proteins 3, 4, and 6 induce double-membrane vesicles

Coronaviruses (CoV), like other positive-stranded RNA viruses, redirect and rearrange host cell membranes for use as part of the viral genome replication and transcription machinery. Specifically, coronaviruses induce the formation of double-membrane vesicles in infected cells. Although these double-membrane vesicles have been well characterized, the mechanism behind their formation remains unclear, including which viral proteins are responsible. Here, we use transfection of plasmid constructs encoding full-length versions of the three transmembrane-containing nonstructural proteins (nsps) of the severe acute respiratory syndrome (SARS) coronavirus to examine the ability of each to induce double-membrane vesicles in tissue culture. nsp3 has membrane disordering and proliferation ability, both in its full-length form and in a C-terminal-truncated form. nsp3 and nsp4 working together have the ability to pair membranes. nsp6 has membrane proliferation ability as well, inducing perinuclear vesicles localized around the microtubule organizing center. Together, nsp3, nsp4, and nsp6 have the ability to induce double-membrane vesicles that are similar to those observed in SARS coronavirus-infected cells. This activity appears to require the full-length form of nsp3 for action, as double-membrane vesicles were not seen in cells coexpressing the C-terminal truncation nsp3 with nsp4 and nsp6. IMPORTANCE Although the majority of infections caused by coronaviruses in humans are relatively mild, the SARS outbreak of 2002 to 2003 and the emergence of the human coronavirus Middle Eastern respiratory syndrome (MERS-CoV) in 2012 highlight the ability of these viruses to cause severe pathology and fatality. Insight into the molecular biology of how coronaviruses take over the host cell is critical for a full understanding of any known and possible future outbreaks caused by these viruses. Additionally, since membrane rearrangement is a tactic used by all known positive-sense single-stranded RNA viruses, this work adds to that body of knowledge and may prove beneficial in the development of future therapies not only for human coronavirus infections but for other pathogens as well.

The effect of the magnetic field of the secondary star in dwarf novae

We investigate the effect of a secondary star magnetic field on the accretion disc dynamics of dwarf novae. Simulations have been carried out with a particle code and a dipolar magnetic field structure. The magnetic field acts to remove angular momentum from the disc material, increasing the inward mass flow. This makes the accretion disc more centrally condensed, causing a reduction in the recurrence time for dwarf nova outbursts. We have produced Doppler tomograms and light curves which may be compared with observations. These tomograms are significantly different from those produced in the absence of a magnetic field on the secondary. We derive an upper limit to the magnetic moment of the secondary star in UGem of mu_2<2x10^32 A m^2. The magnetic truncation of the accretion disc produces resonance phenomena similar to those seen in the superoutbursts of SUUMa systems. While these have not been observed for systems like UGem, observations of the SUUMa systems provide us with a useful diagnostic of the disc-field interaction. We are able to place an upper limit on the magnetic moment of the secondary in ZCha of mu_2<1x10^30 A m^2.

The composite-tendency RAW filter in semi-implicit integrations

The time discretization in weather and climate models introduces truncation errors that limit the accuracy of the simulations. Recent work has yielded a method for reducing the amplitude errors in leapfrog integrations from first-order to fifth-order. This improvement is achieved by replacing the Robert--Asselin filter with the RAW filter and using a linear combination of the unfiltered and filtered states to compute the tendency term. The purpose of the present paper is to apply the composite-tendency RAW-filtered leapfrog scheme to semi-implicit integrations. A theoretical analysis shows that the stability and accuracy are unaffected by the introduction of the implicitly treated mode. The scheme is tested in semi-implicit numerical integrations in both a simple nonlinear stiff system and a medium-complexity atmospheric general circulation model, and yields substantial improvements in both cases. We conclude that the composite-tendency RAW-filtered leapfrog scheme is suitable for use in semi-implicit integrations.