Resonant Wave Interactions in the Presence of a Diurnally Varying Heat Source

Resonant interactions among equatorial waves in the presence of a diurnally varying heat source are studied in the context of the diabatic version of the equatorial beta-plane primitive equations for a motionless, hydrostatic, horizontally homogeneous and stably stratified background atmosphere. The heat source is assumed to be periodic in time and of small amplitude [i.e., O(epsilon)] and is prescribed to roughly represent the typical heating associated with deep convection in the tropical atmosphere. In this context, using the asymptotic method of multiple time scales, the free linear Rossby, Kelvin, mixed Rossby-gravity, and inertio-gravity waves, as well as their vertical structures, are obtained as leading-order solutions. These waves are shown to interact resonantly in a triad configuration at the O(e) approximation, and the dynamics of these interactions have been studied in the presence of the forcing. It is shown that for the planetary-scale wave resonant triads composed of two first baroclinic equatorially trapped waves and one barotropic Rossby mode, the spectrum of the thermal forcing is such that only one of the triad components is resonant with the heat source. As a result, to illustrate the role of the diurnal forcing in these interactions in a simplified fashion, two kinds of triads have been analyzed. The first one refers to triads composed of a k = 0 first baroclinic geostrophic mode, which is resonant with the stationary component of the diurnal heat source, and two dispersive modes, namely, a mixed Rossby-gravity wave and a barotropic Rossby mode. The other class corresponds to triads composed of two first baroclinic inertio-gravity waves in which the highest-frequency wave resonates with a transient harmonic of the forcing. The integration of the asymptotic reduced equations for these selected resonant triads shows that the stationary component of the diurnal heat source acts as an ""accelerator"" for the energy exchanges between the two dispersive waves through the excitation of the catalyst geostrophic mode. On the other hand, since in the second class of triads the mode that resonates with the forcing is the most energetically active member because of the energy constraints imposed by the triad dynamics, the results show that the convective forcing in this case is responsible for a longer time scale modulation in the resonant interactions, generating a period doubling in the energy exchanges. The results suggest that the diurnal variation of tropical convection might play an important role in generating low-frequency fluctuations in the atmospheric circulation through resonant nonlinear interactions.

Resonant Wave Interactions in the Equatorial Waveguide

Weakly nonlinear interactions among equatorial waves have been explored in this paper using the adiabatic version of the equatorial beta-plane primitive equations in isobaric coordinates. Assuming rigid lid vertical boundary conditions, the conditions imposed at the surface and at the top of the troposphere were expanded in a Taylor series around two isobaric surfaces in an approach similar to that used in the theory of surface-gravity waves in deep water and capillary-gravity waves. By adopting the asymptotic method of multiple time scales, the equatorial Rossby, mixed Rossby-gravity, inertio-gravity, and Kelvin waves, as well as their vertical structures, were obtained as leading-order solutions. These waves were shown to interact resonantly in a triad configuration at the O(epsilon) approximation. The resonant triads whose wave components satisfy a resonance condition for their vertical structures were found to have the most significant interactions, although this condition is not excluding, unlike the resonant conditions for the zonal wavenumbers and meridional modes. Thus, the analysis has focused on such resonant triads. In general, it was found that for these resonant triads satisfying the resonance condition in the vertical direction, the wave with the highest absolute frequency always acts as an energy source (or sink) for the remaining triad components, as usually occurs in several other physical problems in fluid dynamics. In addition, the zonally symmetric geostrophic modes act as catalyst modes for the energy exchanges between two dispersive waves in a resonant triad. The integration of the reduced asymptotic equations for a single resonant triad shows that, for the initial mode amplitudes characterizing realistic magnitudes of atmospheric flow perturbations, the modes in general exchange energy on low-frequency (intraseasonal and/or even longer) time scales, with the interaction period being dependent upon the initial mode amplitudes. Potential future applications of the present theory to the real atmosphere with the inclusion of diabatic forcing, dissipation, and a more realistic background state are also discussed.

Predictions from an anisotropic inflationary era

This paper investigates the predictions of an inflationary phase starting from a homogeneous and anisotropic universe of the Bianchi I type. After discussing the evolution of the background spacetime, focusing on the number of e-folds and the isotropization, we solve the perturbation equations and predict the power spectra of the curvature perturbations and gravity waves at the end of inflation. The main features of the early anisotropic phase is (1) a dependence of the spectra on the direction of the modes, (2) a coupling between curvature perturbations and gravity waves and (3) the fact that the two gravity wave polarizations do not share the same spectrum on large scales. All these effects are significant only on large scales and die out on small scales where isotropy is recovered. They depend on a characteristic scale that can, but a priori must not, be tuned to some observable scale. To fix the initial conditions, we propose a procedure that generalizes the one standardly used in inflation but that takes into account the fact that the WKB regime is violated at early times when the shear dominates. We stress that there exist modes that do not satisfy the WKB condition during the shear-dominated regime and for which the amplitude at the end of inflation depends on unknown initial conditions. On such scales, inflation loses its predictability. This study paves the way for the determination of the cosmological signature of a primordial shear, whatever the Bianchi I spacetime. It thus stresses the importance of the WKB regime to draw inflationary predictions and demonstrates that, when the number of e-folds is large enough, the predictions converge toward those of inflation in a Friedmann-Lemaitre spacetime but that they are less robust in the case of an inflationary era with a small number of e-folds.

Interaction of equatorial waves through resonance with the diurnal cycle of tropical heating

In this work, a new theoretical mechanism is presented in which equatorial Rossby and inertio-gravity wave modes may interact with each other through resonance with the diurnal cycle of tropical deep convection. We have adopted the two-layer incompressible equatorial primitive equations forced by a parametric heating that roughly represents deep convection activity in the tropical atmosphere. The heat source was parametrized in the simplest way according to the hypothesis that it is proportional to the lower-troposphere moisture convergence, with the background moisture state function mimicking the structure of the ITCZ. In this context, we have investigated the possibility of resonant interaction between equatorially trapped Rossby and inertio-gravity modes through the diurnal cycle of the background moisture state function. The reduced dynamics of a single resonant duo shows that when this diurnal variation is considered, a Rossby wave mode can undergo significant amplitude modulations when interacting with an inertio-gravity wave mode, which is not possible in the context of the resonant triad non-linear interaction. Therefore, the results suggest that the diurnal variation of the ITCZ can be a possible dynamical mechanism that leads the Rossby waves to be significantly affected by high frequency modes.

Alfven waves as a driving mechanism in stellar winds

Alfven waves have been invoked as an important mechanism of particle acceleration in stellar winds of cool stars. After their identification in the solar wind they started to be studied in winds of stars located in different regions of the FIR diagram. We discuss here some characteristics of these waves and we present a direct application in the acceleration of late-type stellar winds. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

Origin of Convectively Coupled Kelvin Waves over South America

Convectively coupled Kelvin waves over the South American continent are examined through the use of temporal and spatial filtering of reanalysis, satellite, and gridded rainfall data. They are most prominent from November to April, the season analyzed herein. The following two types of events are isolated: those that result from preexisting Kelvin waves over the eastern Pacific Ocean propagating into the continent, and those that apparently originate over Amazonia, forced by disturbances propagating equatorward from central and southern South America. The events with precursors in the Pacific are mainly upper-level disturbances, with almost no signal at the surface. Those events with precursors over South America, on the other hand, originate as upper-level synoptic wave trains that pass over the continent and resemble the ""cold surges`` documented by Garreaud and Wallace. As the wave train propagates over the Andes, it induces a southerly low-level wind that advects cold air to the north. Precipitation associated with a cold front reaches the equator a few days later and subsequently propagates eastward with the characteristics of a Kelvin wave. The structures of those waves originating over the Pacific are quite similar to those originating over South America as they propagate to eastern South America and into the Atlantic. South America Kelvin waves that originate over neither the Pacific nor the midlatitudes of South America can also be identified. In a composite sense, these form over the eastern slope of the Andes Mountains, close to the equator. There are also cases of cold surges that reach the equator yet do not form Kelvin waves. The interannual variability of the Pacific-originating events is related to sea surface temperatures in the central-eastern Pacific Ocean. When equatorial oceanic conditions are warm, there tends to be an increase in the number of disturbances that reach South America from the Pacific.

Wolf-Rayet optically thick winds with Alfven waves

The Wolf-Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfven waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.

Integrated marine gravity field in the Brazilian coast from altimeter-derived sea surface gradient and shipborne gravity

Sea surface gradients derived from the Geosat and ERS-1 satellite altimetry geodetic missions were integrated with marine gravity data from the National Geophysical Data Center and Brazilian national surveys. Using the least squares collocation method, models of free-air gravity anomaly and geoid height were calculated for the coast of Brazil with a resolution of 2` x 2`. The integration of satellite and shipborne data showed better statistical results in regions near the coast than using satellite data only, suggesting an improvement when compared to the state-of-the-art global gravity models. Furthermore, these results were obtained with considerably less input information than was used by those reference models. The least squares collocation presented a very low content of high-frequency noise in the predicted gravity anomalies. This may be considered essential to improve the high resolution representation of the gravity field in regions of ocean-continent transition. (C) 2010 Elsevier Ltd. All rights reserved.

Moho map of South America from receiver functions and surface waves

We estimate crustal structure and thickness of South America north of roughly 40 degrees S. To this end, we analyzed receiver functions from 20 relatively new temporary broadband seismic stations deployed across eastern Brazil. In the analysis we include teleseismic and some regional events, particularly for stations that recorded few suitable earthquakes. We first estimate crustal thickness and average Poisson`s ratio using two different stacking methods. We then combine the new crustal constraints with results from previous receiver function studies. To interpolate the crustal thickness between the station locations, we jointly invert these Moho point constraints, Rayleigh wave group velocities, and regional S and Rayleigh waveforms for a continuous map of Moho depth. The new tomographic Moho map suggests that Moho depth and Moho relief vary slightly with age within the Precambrian crust. Whether or not a positive correlation between crustal thickness and geologic age is derived from the pre-interpolation point constraints depends strongly on the selected subset of receiver functions. This implies that using only pre-interpolation point constraints (receiver functions) inadequately samples the spatial variation in geologic age. The new Moho map also reveals an anomalously deep Moho beneath the oldest core of the Amazonian Craton.

Magnetotelluric deep soundings, gravity and geoid in the south Sao Francisco craton: Geophysical indicators of cratonic lithosphere rejuvenation and crustal underplating

In the south Sao Francisco craton a circular and 8-m amplitude geoid anomaly coincides with the outcropping terrain of an Archean-Paleoproterozoic basement. Broadband magnetotelluric (MT) data inversions of two radial profiles within the positive geoid and Bouguer gravity anomaly yield geo-electrical crustal sections, whereby the lower crust is locally more conductive (10 to 100 Omega m) in spatial coincidence with a denser lower crust modeled by the gravity data. This anomalous lower crust may have resulted from magmatic underplating, associated with Mesoarchean and Proterozoic episodes of tholeiitic dike intrusion. Long-period MT soundings reveal a low electrical resistivity mantle (20 to 200 Omega m) from depths beyond 120 km. Forward geoid modeling, using the scope of the low electrical resistivity region within the mantle as a constraint, entails a density increase (40 to 50 kg/m(3)) possibly due to Fe enrichment of mantle minerals. However, this factor alone does not explain the observed resistivity. A supplemented presence of small amounts of percolated carbonatite melting (similar to 0.005 vol.%), dissolved water and enhanced oxygen fugacity within the peridotitic mantle are viable agents that could explain the less resistive upper mantle. We propose that metasomatic processes confined in the sub-continental lithospheric mantle foster the conditions for a low degree melting with variable CO(2), H(2)O and Fe content. Even though the precise age of this metasomatism is unknown it might be older than the Early Cretaceous based on the evidence that a high-degree of melting in a lithospheric mantle impregnated with carbonatites originated the tholeiitic dike intrusions dispersed from the southeastern border of the Sao Francisco craton, during the onset of the lithosphere extension and break-up of the western Gondwana. The proxies are the NE Parana and Espinhaco (130 Ma, Ar/Ar ages) tholeiitic dikes, which contain (similar to 3%) carbonatites in their composition. The occurrence of a positive geoid anomaly (+ 10 m) and pre-tholeiites (age > 138 Ma), carbonatites and kimberlites along the west African continental margin (Angola and Namibia) reinforces the presumed age of the Sao Francisco-Congo craton rejuvenation to be prior to its fragmentation in the Lower Cretaceous. (C) 2010 Elsevier B.V. All rights reserved.

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.

Programs to compute magnetization to density ratio and the magnetization inclination from 3-D gravity and magnetic anomalies

Vector field formulation based on the Poisson theorem allows an automatic determination of rock physical properties (magnetization to density ratio-MDR-and the magnetization inclination-MI) from combined processing of gravity and magnetic geophysical data. The basic assumptions (i.e., Poisson conditions) are: that gravity and magnetic fields share common sources, and that these sources have a uniform magnetization direction and MDR. In addition, the previously existing formulation was restricted to profile data, and assumed sufficiently elongated (2-D) sources. For sources that violate Poisson conditions or have a 3-D geometry, the apparent values of MDR and MI that are generated in this way have an unclear relationship to the actual properties in the subsurface. We present Fortran programs that estimate MDR and MI values for 3-D sources through processing of gridded gravity and magnetic data. Tests with simple geophysical models indicate that magnetization polarity can be successfully recovered by MDR-MI processing, even in cases where juxtaposed bodies cannot be clearly distinguished on the basis of anomaly data. These results may be useful in crustal studies, especially in mapping magnetization polarity from marine-based gravity and magnetic data. (c) 2007 Elsevier Ltd. All rights reserved.

Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.

BREAKING WAVES IN AN IDEAL QUARK GLUON PLASMA

We study the propagation of perturbations in the energy density in a quark gluon plasma. Expanding the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations we obtain a nonlinear differential equation called the breaking wave equation. We solve it numerically and follow the time-evolution of initially localized pulses. We find that, quite unexpectedly, these pulses live for a very long time (compared to the reaction time-scales) before breaking. In practice, they mimick the Korteweg-de Vries solitons. Their existence may have some observable consequences.

Sound waves and solitons in hot and dense nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density. The equation of state is derived from a relativistic mean field model, which is a variant of the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations leads to differential equations for the density perturbation. We solve them numerically for linear and spherical perturbations and follow the propagation of the initial pulses. For linear perturbations we find single soliton solutions and solutions with one or more solitons followed by ""radiation"". Depending on the equation of state a strong damping may occur. We consider also the evolution of perturbations in a medium without dispersive effects. In this case we observe the formation and breaking of shock waves. We study all these equations also for matter at finite temperature. Our results may be relevant for the analysis of RHIC data. They suggest that the shock waves formed in the quark gluon plasma phase may survive and propagate in the hadronic phase. (C) 2009 Elseiver. B.V. All rights reserved.