Characteristics of planetary-scale waves simulated by a new venusian mesosphere and thermosphere general circulation model

We have developed a new general circulation model (GCM) for the venusian mesosphere and thermosphere (80-about 180 km). Our GCM simulations show that winds in the subsolar-to-antisolar direction (SS–AS) are predominant above about 90 km. A weak return flow of the SS–AS is seen below about 90 km. We performed GCM simulations imposing the planetary-scale waves (thermal tides, Rossby wave, and Kelvin wave) at the lower boundary. Although the diurnal and semidiurnal tides are damped below 95 km, the Rossby wave propagates up to around 130 km. However, the amplitude of the Rossby wave is too small (<1 m/s) to affect the general circulation. On the other hand, the Kelvin wave propagates up to about 130 km with a maximum zonal wind fluctuation of approximately 5.9 m/s on average. The amplitude of the Kelvin wave sometimes exceeds 10 m/s around the terminator. The Kelvin wave causes a temporal variation in the wind velocity at the altitude of the O₂-1.27 μm nightglow emission (about 95 km). Using a newly developed 1-D nightglow model and the composition distribution calculated from our GCM, we investigated the impact of the Kelvin wave on the nightglow distribution. Our results suggest that the Kelvin wave would cause temporal variations in the nightglow emission in the 23:50–00:20 LT region with an intensity of 1.1–1.3 MR and a period of approximately 4 days.     

Fourier Analysis of Gapped Time Series: Improved Estimates of Solar and Stellar Oscillation Parameters

Quantitative helioseismology and asteroseismology require very precise measurements of the frequencies, amplitudes, and lifetimes of the global modes of stellar oscillation. The precision of these measurements depends on the total length (T), quality, and completeness of the observations. Except in a few simple cases, the effect of gaps in the data on measurement precision is poorly understood, in particular in Fourier space where the convolution of the observable with the observation window introduces correlations between different frequencies. Here we describe and implement a rather general method to retrieve maximum likelihood estimates of the oscillation parameters, taking into account the proper statistics of the observations. Our fitting method applies in complex Fourier space and exploits the phase information. We consider both solar-like stochastic oscillations and long-lived harmonic oscillations, plus random noise. Using numerical simulations, we demonstrate the existence of cases for which our improved fitting method is less biased and has a greater precision than when the frequency correlations are ignored. This is especially true of low signal-to-noise solar-like oscillations. For example, we discuss a case where the precision of the mode frequency estimate is increased by a factor of five, for a duty cycle of 15%. In the case of long-lived sinusoidal oscillations, a proper treatment of the frequency correlations does not provide any significant improvement; nevertheless, we confirm that the mode frequency can be measured from gapped data with a much better precision than the 1/T Rayleigh resolution.     

Inference of dust opacities for the 1977 Martian great dust storms from Viking Lander 1 pressure data

A δ-Eddington radiative transfer algorithm is used to compute the thermal tidal heating of a dusty Martian atmosphere for a given set of dust optical depth, effective single scattering albedo, and phase function asymmetry parameter. The resulting thermal tidal forcing is used in a classical atmospheric tidal model to compute the amplitudes of the surface pressure oscillations at the Viking Lander 1 site for the two 1977 Martian great dust storms. Parametric studies show that the dust opacities and optical parameters derived from the Viking Lander imaging data are roughly representative of the global dust haze needed to reproduce the tidal surface pressure amplitudes also observed at Lander 1, except that the model-inferred asymmetry parameter is smaller during the onset of a great storm. The observed preferential enhancement during dust-storm onset of the semidiurnal tide at Viking Lander 1 relative to its diurnal counterpart is shown to be due primarily to the elevation of the tidal heating source in a very dusty atmosphere, although resonant enhancement of the main semidiurnal tidal mode makes an important secondary contribution.     

Period investigation of two RS CVn-type binary stars: RU Cancri and AW Herculis

Orbital period variations of two RS CVn-type binaries, RU Cnc and AW Her, are presented based on the analysis of all available times of light minima. It is discovered that the orbital period of RU Cnc shows two possible period oscillations with periods of 13.38(±0.23) and 37.6(±3.4) years. The corresponding amplitudes of the oscillations are about 0.0098(±0.0023) and 0.0119(±0.0017) days. For AW Her, it is found that the period shows a cyclic variation with a period of 12.79(±0.34) years and an amplitude of about 0.0327(±0.0063) days. Since RU Cnc and AW Her are two RS CVn-type systems, the cyclic period oscillations are more likely to be caused by the magnetic activity cycles.     

Nonlinear Numerical Simulations of Magneto-Acoustic Wave Propagation in Small-Scale Flux Tubes

We present results of nonlinear, two-dimensional, numerical simulations of magneto-acoustic wave propagation in the photosphere and chromosphere of small-scale flux tubes with internal structure. Waves with realistic periods of three to five minutes are studied, after horizontal and vertical oscillatory perturbations are applied to the equilibrium model. Spurious reflections of shock waves from the upper boundary are minimized by a special boundary condition. This has allowed us to increase the duration of the simulations and to make it long enough to perform a statistical analysis of oscillations. The simulations show that deep horizontal motions of the flux tube generate a slow (magnetic) mode and a surface mode. These modes are efficiently transformed into a slow (acoustic) mode in the v _A<c _S atmosphere. The slow (acoustic) mode propagates vertically along the field lines, forms shocks, and remains always within the flux tube. It might effectively deposit the energy of the driver into the chromosphere. When the driver oscillates with a high frequency, above the cutoff, nonlinear wave propagation occurs with the same dominant driver period at all heights. At low frequencies, below the cutoff, the dominant period of oscillations changes with height from that of the driver in the photosphere to its first harmonic (half period) in the chromosphere. Depending on the period and on the type of the driver, different shock patterns are observed.     

On searching for observational manifestations of Alfvén waves in solar faculae

In an effort to detect torsional oscillations, we have studied the periodic half-width variations for several spectral lines in solar faculae. The duration of the series being analyzed was from 40 to 150 min. We have determined the dominant frequencies and amplitudes of the half-width oscillations and considered their phase relations to the intensity and line-of-sight velocity oscillations. Five-minute profile halfwidth oscillations with a peak-to-peak amplitude of ～10 m ?A are recorded with confidence in the upperphotospheric Si I 10 827 ?A line in faculae. The chromospheric He I 10 830 A? and Hα line profiles shows ～40–60 m ?A variations in two frequency bands, 2.5–4 and 1–1.9 mHz. No center-to-limb dependence that, according to the theory, must accompany the torsional oscillations has been revealed in the behavior of the oscillation amplitudes. According to present views, these variations cannot be caused by periodic temperature and magnetic field changes. Our observations do not allow us to explain these variations by the sausage mode action either, which should manifest itself at the double frequency.     

Radio occultation measurements and MGCM simulations of Kelvin waves on Mars

We have derived new results concerning thermal tides on Mars from a combination of radio occultation measurements and numerical simulations by a Mars General Circulation Model (MGCM). This investigation exploits a set of concurrent observations by Mars Express (MEX) and Mars Global Surveyor (MGS) in mid-2004, when the season on Mars was midspring in the northern hemisphere. The MEX occultations sampled the atmosphere near the evening terminator at latitudes ranging from 54° N to 15° S. The MGS occultations provided complementary coverage near the morning terminator at latitudes of 35° N and 71° S. The geopotential field derived from these measurements contains distinctive modulation caused by solar-asynchronous thermal tides. Through careful analysis of the combined observations, we characterized two prominent wave modes, obtaining direct solutions for some properties, such as the amplitude and phase, as well as constraints on others, such as the period, zonal wave number, and meridional structure. We supplemented these observations with MGCM simulations. After evaluating the performance of the MGCM against the measurements, we used the validated simulation to deduce the identity of the two tidal modes and to explore their behavior. One mode is a semidiurnal Kelvin wave with a zonal wave number of 2 (SK2), while the other is a diurnal Kelvin wave with a zonal wave number of 1 (DK1). Both modes are known to be close to resonance in the martian atmosphere. Our observations of the SK2 are more complete and less ambiguous than any previous measurement. The well-known DK1 is the dominant solar-asynchronous tide in the martian atmosphere, and our results confirm and extend previous observations by diverse instruments.     

Interannual variability of regional climate and its change due to the greenhouse effect

Interannual variability of regional climate was investigated on a seasonal basis. Observations and two global climate model (GCM) simulations were intercompared to identify model biases and climate change signals due to the enhanced greenhouse effect. Observed record length varies from 40 to 100 years, while the model output comes from two 100-year equilibrium climate simulations corresponding to atmospheric greenhouse gas concentrations at observed 1990 and projected 2050 levels. The GCM includes an atmosphere based on the NCAR CCM1 with the addition of the radiative effects of CH₄, N₂O and CFCs, a bulk layer land surface and a mixed-layer ocean with thermodynamic sea-ice and fixed meridional oceanic heat transport.Because comparisons of interannual variability are sensitive to the time period chosen, a climate ensemble technique has been developed. This technique provides comparisons between variance ratios of two time series for all possible contiguous sub-periods of a fixed length. The time autocorrelation is thus preserved within each sub-period. The optimal sub-period length was found to be 30 years, based on which robust statistics of the ensemble were obtained to identify substantial differences in interannual variability that are both physically important and statistically significant.Several aspects of observed interannual variability were reproduced by the GCM. These include: global surface air temperature; Arctic sea-ice extent; and regional variability of surface air temperature, sea level pressure and 500 mb height over about one quarter of the observed data domains. Substantial biases, however, exist over broad regions, where strong seasonality and systematic links between variables were identified. For instance, during summer substantially greater model variability was found for both surface air temperature and sea-level pressure over land areas between 20–50°N, while this tendency was confined to 20–30°N in other seasons. When greenhouse gas concentrations increase, atmospheric moisture variability is substantially larger over areas that experience the greatest surface warming. This corresponds to an intensified hydrologic cycle and, hence, regional increases in precipitation variability. Surface air temperature variability increases where hydrologic processes vary greatly or where mean soil moisture is much reduced. In contrast, temperature variability decreases substantially where sea-ice melts completely. These results indicate that regional changes in interannual variability due to the enhanced greenhouse effect are associated with mechanisms that depend on the variable and season.     

The vector approach to the problem of physical libration of the Moon: the linearized problem

A flexible and informative vector approach to the problem of physical libration of the rigid Moon has been developed in which three Euler differential equations are supplemented by 12 kinematic ones. A linearized system of equations can be split into an even and odd systems with respect to the reflection in the plane of the lunar equator, and rotational oscillations of the Moon are presented by superposition of librations in longitude and latitude. The former is described by three equations and consists of unrestricted oscillations with a period of T ₁ = 2.878 Julian years (amplitude of 1.855″) and forced oscillations with periods of T ₂ = 27.201 days (15.304″), one stellar year (0.008″), half a year (0.115″), and the third of a year (0.0003″) (five harmonics altogether). A zero frequency solution has also been obtained. The effect of the Sun on these oscillations is two orders of magnitude less than that of the Earth. The libration in latitude is presented by five equations and, at pertrubations from the Earth, is described by two harmonics of unrestricted oscillations (T ₅ ≈ 74.180 Julian years, T ₆ ≈ 27.347 days) and one harmonic of forced oscillations (T ₃ = 27.212 days). The motion of the true pole is presented by the same harmonics, with the maximum deviation from the Cassini pole being 45.3″. The fifth (zero) frequency yields a stationary solution with a conic precession of the rotation axis (previously unknown). The third Cassini law has been proved. The amplitudes of unrestricted oscillations have been determined from comparison with observations. For the ratio $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 , the theory gives 0.2319, which confirms the adequacy of the approach. Some statements of the previous theory are revised. Poinsot’s method is shown to be irrelevant in describing librations of the Moon. The Moon does not have free (Euler) oscillations; it has oscillations with a period of T ₅ ≈ 74.180 Julian years rather than T ≈ 148.167 Julian years.     

Multisite campaign on the open cluster M67 – II. Evidence for solar-like oscillations in red giant stars

Measuring solar-like oscillations in an ensemble of stars in a cluster, holds promise for testing stellar structure and evolution more stringently than just fitting parameters to single field stars. The most-ambitious attempt to pursue these prospects was by Gilliland et al. who targeted 11 turn-off stars in the open cluster M67 (NGC 2682), but the oscillation amplitudes were too small (<20 μmag) to obtain unambiguous detections. Like Gilliland et al. we also aim at detecting solar-like oscillations in M67, but we target red giant stars with expected amplitudes in the range 50–  500 μmag  and periods of 1 to 8 h. We analyse our recently published photometry measurements, obtained during a six-week multisite campaign using nine telescopes around the world. The observations are compared with simulations and with estimated properties of the stellar oscillations. Noise levels in the Fourier spectra as low as  27 μmag  are obtained for single sites, while the combined data reach  19 μmag  , making this the best photometric time series of an ensemble of red giant stars. These data enable us to make the first test of the scaling relations (used to estimate frequency and amplitude) with an homogeneous ensemble of stars. The detected excess power is consistent with the expected signal from stellar oscillations, both in terms of its frequency range and amplitude. However, our results are limited by apparent high levels of non-white noise, which cannot be clearly separated from the stellar signal.     

Flow around sunspots from it p-mode frequency shift measurements

We use five-day helioseimic data from the Taiwan Oscillation Network to study the flow around a sunspot, NOAA 7887. The p-mode oscillations in an annular region centered at the sunspot are decomposed into the modes propagating toward and away from the sunspot. We find that the frequency of an outgoing mode is greater than that of the corresponding incoming mode. This indicates that the plasma is flowing outward from the sunspot. The outflow velocity is estimated to be about 40–80 m s^-1.     

Titan global climate model: A new 3-dimensional version of the IPSL Titan GCM

We have developed a new 3-dimensional climate model for Titan’s atmosphere, using the physics of the IPSL Titan 2-dimensional climate model with the current version of the LMDZ General Circulation Model dynamical core. Microphysics and photochemistry are still computed as zonal averages. This GCM covers altitudes from surface to 500 km altitude, with barotropic waves now being resolved and the diurnal cycle included. The boundary layer scheme has been changed, yielding a strong improvement in the tropospheric zonal wind profile modeled at Huygens descent position and season. The potential temperature profile is fairly consistent with Huygens observations in the lowest 10 km. The latitudinal profile of the near-surface temperature is close to observed values. The minimum of zonal wind observed by the Huygens probe just above the tropopause is also present in these simulations, and its origin is discussed by comparing solar heating and dynamical transport of energy. The stratospheric temperature and wind fields are consistent with our previous works. Compared to observations, the zonal wind peak is too weak (around 120 m/s) and too low (around 200 km). The temperature structures appear to be compressed in altitude, and depart strongly from observations in the upper stratosphere. These discrepancies are correlated, and most probably related to the altitude of the haze production. The model produces a detached haze layer located more than 150 km lower than observed by the Cassini instruments. This low production altitude is due to the current position of the GCM upper boundary. However, the temporal behaviour of the detached haze layer in the model may explain the seasonal differences observed between Cassini and Voyager 1. The waves present in the GCM are analyzed, together with their respective roles in the angular momentum budget. Though the role of the mean meridional circulation in momentum transport is similar to previous work, and the transport by barotropic waves is clearly seen in the stratosphere, a significant part of the transport at high latitudes is done all year long through low-frequency tropospheric waves that may be baroclinic waves.     

Measurement of the solar diurnal anisotropy of the cosmic-ray albedo neutron flux

The solar diurnal anisotropy of the cosmic-ray albedo neutron flux has been measured by a neutron detector on board the OGO-6 satellite. On the average the diurnal amplitudes and phases of the cosmic ray albedo neutron flux (10 MeV) were respectively 0.18 (±0.02)% and 15(±1) hr LT though there were substantial fluctuations of a few days duration which did not depend on the solar sector structure polarity and a 27-day periodicity in the diurnal amplitudes which was associated with the Sun's rotation.     

Gaussian process modelling of asteroseismic data

The measured properties of stellar oscillations can provide powerful constraints on the internal structure and composition of stars. To begin this process, oscillation frequencies must be extracted from the observational data, typically time series of the star's brightness or radial velocity. In this paper, a probabilistic model is introduced for inferring the frequencies and amplitudes of stellar oscillation modes from data, assuming that there is some periodic character to the oscillations, but that they may not be exactly sinusoidal. Effectively, we fit damped oscillations to the time series, and hence the mode lifetime is also recovered. While this approach is computationally demanding for large time series (>1500 points), it should at least allow improved analysis of observations of solar-like oscillations in subgiant and red giant stars, as well as sparse observations of semiregular stars, where the number of points in the time series is often low. The method is demonstrated on simulated data and then applied to radial velocity measurements of the red giant star  ξ Hydrae  , yielding a mode lifetime between 0.41 and 2.65 d with 95 per cent posterior probability. The large frequency separation between modes is ambiguous, however we argue that the most plausible value is 6.3 μHz, based on the radial velocity data and the star's position in the Hertzsprung–Russell diagram.     

ALMA as the ideal probe of the solar chromosphere

The very nature of the solar chromosphere, its structuring and dynamics, remains far from being properly understood, in spite of intensive research. Here we point out the potential of chromospheric observations at millimeter wavelengths to resolve this long-standing problem. Computations carried out with a sophisticated dynamic model of the solar chromosphere due to Carlsson and Stein demonstrate that millimeter emission is extremely sensitive to dynamic processes in the chromosphere and the appropriate wavelengths to look for dynamic signatures are in the range 0.8–5.0 mm. The model also suggests that high resolution observations at mm wavelengths, as will be provided by ALMA, will have the unique property of reacting to both the hot and the cool gas, and thus will have the potential of distinguishing between rival models of the solar atmosphere. Thus, initial results obtained from the observations of the quiet Sun at 3.5 mm with the BIMA array (resolution of 12″) reveal significant oscillations with amplitudes of 50–150 K and frequencies of 1.5–8 mHz with a tendency toward short-period oscillations in internetwork and longer periods in network regions. However higher spatial resolution, such as that provided by ALMA, is required for a clean separation between the features within the solar atmosphere and for an adequate comparison with the output of the comprehensive dynamic simulations.     

Short-period intensity oscillations in the solar corona observed during the total solar eclipse of 26 February 1998

Encouraged by the detection of high-frequency, low-amplitude continuum intensity oscillations in the solar corona during the total solar eclipse of 1995, we designed and fabricated a six-channel photometer incorporating low-noise Hamamatsu R647 photomultipliers. Fast photometry at five different locations in the solar corona was performed at Don Bosco Mission, Venezuela during the total solar eclipse of 26 February 1998. Three interference filters with passbands of about 150 Å and centered around 4700, 4900, and 5000 Å were used. The photometric data were recorded at a rate of 20 Hz in three channels and 50 Hz in the remaining three channels. The power spectrum analysis of one of the channels that recorded appreciable counts indicates the existence of intensity oscillations in the frequency range 0.01–0.2 Hz. A least-squares analysis yields 90.1, 25.2, and 6.9 s periods for the three prominent components which have amplitudes in the range 0.5–3.5% of the coronal brightness. These periods and their amplitudes are similar to those detected in the coronal intensity oscillations during the 1995 eclipse.     

Fluctuations in tidal (24-, 12-h) characteristics and oscillations (8-h-5-d) in the mesosphere and lower thermosphere (70–110 km): Saskatoon (52°N, 107°W), 1979–1981

An M.F. radar (2.2 MHz) operating at Saskatoon, Canada (52°N, 107°W) has been used to produce continuous wind data ( 80–110km) from September 1978–April 1981. The 24-, 12-h tidal oscillations reveal regular summer-winter transitions; in particular the semi-diurnal tide demonstrates strikingly regular and rapid equinoctial changes over the three years. The vernal and autumnal equinox changes are clearly different in morphology. Shorter term tidal fluctuations (2d τ 10d) are compared with mean winds and gravity wave amplitudes, as well as with satellite-derived stratospheric temperatures.

Spectral analysis of monthly data sets for 1980, from 90–105 km, reveal oscillations of the expected 8-, 12-, 24-h periods, but also of 10-, 16- and 2-, 5/6d. A modulation of the “2-d” wave by the 12-h wave is suggested as a possible cause of these surprisingly regular oscillations.     

Semi-regular and alternating pulsations in hydrodynamical models for low-mass giants

On the basis of hydrodynamical simulations the semi-regular and alternating pulsations of the W Vir and RV Tau stars are studied including non-steady radiative transfer effects in gray sphericalsymmetric atmospheres. It is shown that the transition from regular to semi-regular oscillations occurs at periodsP13–15 days due to pulsation energy imbalance provoqued by strong radiating shocks. The differences between the periods and amplitudes of W Vir and RV Tau variables are explained as a result of this imbalance. It is suggested that the semi-regular alternating pulsations of the RV Tau stars appear following a period-doubling phenomenon, which can be also understood in terms of the pulsation energy variations between odd and even fundamental cycles. The pulsational characteristics of the period-doubled models are in general accordance with those observed in the RV Tau stars. On the basis of numberical results, a theoretical estimation of the upper limit for the luminosity of the RV Tau stars is derived asL10³ L .