The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

Comparison of terrain following and cut cell grids using a non-hydrostatic model

Terrain following coordinates are widely used in operational models but the cut cell method has been proposed as an alternative that can more accurately represent atmospheric dynamics over steep orography. Because the type of grid is usually chosen during model implementation, it becomes necessary to use different models to compare the accuracy of different grids. In contrast, here a C-grid finite volume model enables a like-for-like comparison of terrain following and cut cell grids. A series of standard two-dimensional tests using idealised terrain are performed: tracer advection in a prescribed horizontal velocity field, a test starting from resting initial conditions, and orographically induced gravity waves described by nonhydrostatic dynamics. In addition, three new tests are formulated: a more challenging resting atmosphere case, and two new advection tests having a velocity field that is everywhere tangential to the terrain following coordinate surfaces. These new tests present a challenge on cut cell grids. The results of the advection tests demonstrate that accuracy depends primarily upon alignment of the flow with the grid rather than grid orthogonality. A resting atmosphere is well-maintained on all grids. In the gravity waves test, results on all grids are in good agreement with existing results from the literature, although terrain following velocity fields lead to errors on cut cell grids. Due to semi-implicit timestepping and an upwind-biased, explicit advection scheme, there are no timestep restrictions associated with small cut cells. We do not find the significant advantages of cut cells or smoothed coordinates that other authors find.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

Electroweak precision constraints on the Lee-Wick standard model

We perform an analysis of the electroweak precision observables in the Lee-Wick Standard Model. The most stringent restrictions come from the S and T parameters that receive important tree level and one loop contributions. In general the model predicts a large positive S and a negative T. To reproduce the electroweak data, if all the Lee-Wick masses are of the same order, the Lee-Wick scale is of order 5 TeV. We show that it is possible to find some regions in the parameter space with a fermionic state as light as 2.4-3.5 TeV, at the price of rising all the other masses to be larger than 5-8 TeV. To obtain a light Higgs with such heavy resonances a fine-tuning of order a few per cent, at least, is needed. We also propose a simple extension of the model including a fourth generation of Standard Model fermions with their Lee-Wick partners. We show that in this case it is possible to pass the electroweak constraints with Lee-Wick fermionic masses of order 0.4-1.5 TeV and Lee-Wick gauge masses of order 3 TeV.

A new spin-2 self-dual model in D=2+1

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

SU(5)circle times Z(13) grand unification model

We propose an SU(5) grand unified model with an invisible axion and the unification of the three coupling constants which is in agreement with the values, at M(Z), of alpha, alpha(s), and sin(2)theta(W). A discrete, anomalous, Z(13) symmetry implies that the Peccei-Quinn symmetry is an automatic symmetry of the classical Lagrangian protecting, at the same time, the invisible axion against possible semiclassical gravity effects. Although the unification scale is of the order of the Peccei-Quinn scale the proton is stabilized by the fact that in this model the standard model fields form the SU(5) multiplets completed by new exotic fields and, also, because it is protected by the Z(13) symmetry.

Godel type solution in teleparallel gravity

The stationary cosmological model without closed timelike curves of Godel type is obtained for the ideal dust matter source within the framework of the teleparallel gravity. For a specific choice of the teleparallel gravity parameters, this solution reproduces the causality violating stationary Godel solution in general relativity, in accordance with the teleparallel equivalent of general relativity. The relation between the axial-vector torsion and the cosmic vorticity is clarified. (C) 2004 Elsevier B.V. All rights reserved.

Metric-affine approach to teleparallel gravity

The teleparallel gravity theory, treated physically as a gauge theory of translations, naturally represents a particular case of the most general gauge-theoretic model based on the general affine group of spacetime. on the other hand, geometrically, the Weitzenbock spacetime of distant parallelism is a particular case of the general metric-affine spacetime manifold. These physical and geometrical facts offer a new approach to teleparallelism. We present a systematic treatment of teleparallel gravity within the framework of the metric-affine theory. The symmetries, conservation laws and the field equations are consistently derived, and the physical consequences are discussed in detail. We demonstrate that the so-called teleparallel GR-equivalent model has a number of attractive features which distinguishes it among the general teleparallel theories, although it has a consistency problem when dealing with spinning matter sources.

Fierz-Pauli higher derivative gravity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Cosmic acceleration from second order gauge gravity

We construct a phenomenological theory of gravitation based on a second order gauge formulation for the Lorentz group. The model presents a long-range modification for the gravitational field leading to a cosmological model provided with an accelerated expansion at recent times. We estimate the model parameters using observational data and verify that our estimative for the age of the Universe is of the same magnitude than the one predicted by the standard model. The transition from the decelerated expansion regime to the accelerated one occurs recently (at similar to 9.3 Gyr).

Higgsing M2 to D2 with gravity: N=6 chiral supergravity from topologically gauged ABJM theory

We present the higgsing of three-dimensional N = 6 superconformal ABJM type theories coupled to conformal supergravity, so called topologically gauged ABJM theory, thus providing a gravitational extension of previous work on the relation between N M2 and N D2-branes. The resulting N = 6 supergravity theory appears at a chiral point similar to that of three-dimensional chiral gravity introduced recently by Li, Song and Strominger, but with the opposite sign for the Ricci scalar term in the lagrangian. We identify the supersymmetry in the broken phase as a particular linear combination of the supersymmetry and special conformal supersymmetry in the original topologically gauged ABJM theory. We also discuss the higgsing procedure in detail paying special attention to the role played by the U(1) factors in the original ABJM model and the U(1) introduced in the topological gauging.

ELKO SPINOR FIELDS: LAGRANGIANS FOR GRAVITY DERIVED FROM SUPERGRAVITY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Gravity, antigravity and gravitational shielding in (2+1) dimensions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Investigations in massive 3D gravity

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Two-dimensional background field gravity: A Hamilton-Jacobi analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)