We present results of simulations with the Max Planck Institute general circulation model (MPI-MGCM) implementing a hydrological cycle scheme. The simulations reveal a seasonal water "pump" mechanism responsible for the upward transport of water vapor. This mechanism occurs in high latitudes above 60$^\circ$ of the southern hemisphere at perihelion, when the upward branch of the meridional circulation is particularly strong. A combination of the mean vertical flux with variations induced by solar tides facilitates penetration of water across the "bottleneck" at approximately 60 km. The meridional circulation then transports water across the globe to the northern hemisphere. Since the intensity of the meridional cell is tightly controlled by airborne dust, the water abundance in the thermosphere strongly increases during dust storms.

Carbon dioxide (CO$_2$) ice clouds have been routinely observed in the middle atmosphere of Mars. However, there are still uncertainties concerning physical mechanisms that control their altitude, geographical, and seasonal distributions. Using the Max Planck Institute Martian General Circulation Model (MPI-MGCM), incorporating a state-of-the-art whole atmosphere subgrid-scale gravity wave parameterization [Yi\u{g}it et al., 2008], we demonstrate that internal gravity waves generated by lower atmospheric weather processes have wide reaching impact on the Martian climate. Globally, GWs cool the upper atmosphere of Mars by $\sim$10 % and facilitate high-altitude CO$_2$ ice cloud formation. CO$_2$ ice cloud seasonal variations in the mesosphere and the mesopause region appreciably coincide with the spatio-temporal variations of GW effects, providing insight into the observed distribution of clouds. Our results suggest that GW propagation and dissipation constitute a necessary physical mechanism for CO$_2$ ice cloud formation in the Martian upper atmosphere during all seasons.

We review the measurements of magnetic fields of OBA stars. Based on these data we confirm that magnetic fields are distributed according to a lognormal law with a mean log(B)=-0.5 (B in kG) with a standard deviation sigma=0.5. The shape of the magnetic field distribution is similar to that for neutron stars. This finding is in favor of the hypothesis that the magnetic field of a neutron star is determined mainly by the magnetic field of its predecessor, the massive OB star. Further, we model the evolution of an ensemble of magnetic massive stars in the Galaxy. We use our own population synthesis code to obtain the distribution of stellar radii, ages, masses, temperatures, effective magnetic fields and magnetic fluxes from the pre-main sequence (PMS) via zero age main sequence (ZAMS) up to the terminal age main sequence (TAMS) stages. A comparison of the obtained in our model magnetic field distribution (MFD) with that obtained from the recent measurements of the stellar magnetic field allows us to conclude that the evolution of magnetic fields of massive stars is slow if not absent. The shape of the real MFD shows no indications of the magnetic desert proposed previously. Based on this finding we argue that the observed fraction of magnetic stars is determined by physical conditions at the PMS stage of stellar evolution.

Effects of subgrid-scale gravity waves (GWs) on the diurnal migrating tides are investigated from the mesosphere to the upper thermosphere for September equinox conditions, using a general circulation model coupled with the extended spectral nonlinear GW parameterization of Yi\u{g}it et al (2008). Simulations with GW effects cut-off above the turbopause and included in the entire thermosphere have been conducted. GWs appreciably impact the mean circulation and cool the thermosphere down by up to 12-18%. GWs significantly affect the winds modulated by the diurnal migrating tide, in particular in the low-latitude mesosphere and lower thermosphere and in the high-latitude thermosphere. These effects depend on the mutual correlation of the diurnal phases of the GW forcing and tides: GWs can either enhance or reduce the tidal amplitude. In the low-latitude MLT, the correlation between the direction of the deposited GW momentum and the tidal phase is positive due to propagation of a broad spectrum of GW harmonics through the alternating winds. In the Northern Hemisphere high-latitude thermosphere, GWs act against the tide due to an anti-correlation of tidal wind and GW momentum, while in the Southern high-latitudes they weakly enhance the tidal amplitude via a combination of a partial correlation of phases and GW-induced changes of the circulation. The variable nature of GW effects on the thermal tide can be captured in GCMs provided that a GW parameterization (1) considers a broad spectrum of harmonics, (2) properly describes their propagation, and (3) correctly accounts for the physics of wave breaking/saturation.

We review the measurements of magnetic fields of OB stars and compile a catalog of magnetic OB stars. Based on available data we confirm that magnetic field values are distributed according to a log--normal law with a mean log(B)=2.53 and a standard deviation $\sigma=0.54$. We also investigate the formation of the magnetic field of OBA stars before the Main Sequence (MS).

Gravity waves are primarily generated in the lower atmosphere, and can reach thermospheric heights in the course of their propagation. This paper reviews the recent progress in understanding the role of gravity waves in vertical coupling during sudden stratospheric warmings. Modeling of gravity wave effects is briefly reviewed, and the recent developments in the field are presented. Then, the impact of these waves on the general circulation of the upper atmosphere is outlined. Finally, the role of gravity waves in vertical coupling between the lower and the upper atmosphere is discussed in the context of sudden stratospheric warmings.

Global characteristics of the small-scale gravity wave (GW) field in the Martian atmosphere obtained from a high-resolution general circulation model (GCM) are presented for the first time. The simulated GW-induced temperature variances are in a good agreement with available radio occultation data in the lower atmosphere between 10 and 30 km. The model reveals a latitudinal asymmetry with stronger wave generation in the winter hemisphere, and two distinctive sources of GWs: mountainous regions and the meandering winter polar jet. Orographic GWs are filtered while propagating upward, and the mesosphere is primarily dominated by harmonics with faster horizontal phase velocities. Wave fluxes are directed mainly against the local wind. GW dissipation in the upper mesosphere generates body forces of tens of m~s$^{-1}$~sol$^{-1}$, which tend to close the simulated jets. The results represent a realistic surrogate for missing observations, which can be used for constraining GW parameterizations and validating GCM simulations.

First high-altitude observations of gravity wave (GW)-induced CO$_2$ density perturbations in the Martian thermosphere retrieved from NASA's NGIMS instrument on board the MAVEN satellite are presented and interpreted using the extended GW parameterization of Yi\u{g}it et al. [2008] and the Mars Climate Database as an input. Observed relative density perturbations between 180-220 km of 20-40 % demonstrate appreciable local time, latitude, and altitude variations. Modeling for the spatiotemporal conditions of the MAVEN observations suggests that GWs can directly propagate from the lower atmosphere to the thermosphere, produce appreciable dynamical effects, and likely contribute to the observed fluctuations. Modeled effects are somewhat smaller than the observed but their highly variable nature is in qualitative agreement with observations. Possible reasons for discrepancies between modeling and measurements are discussed.

We present the first general circulation model simulations that quantify and reproduce patches of extremely cold air required for CO$_2$ condensation and cloud formation in the Martian mesosphere. They are created by subgrid-scale gravity waves (GWs) accounted for in the model with the interactively implemented spectral parameterization. Distributions of GW-induced temperature fluctuations and occurrences of supersaturation conditions are in a good agreement with observations of high-altitude CO$_2$ ice clouds. Our study confirms the key role of GWs in facilitating CO$_2$ cloud formation, discusses their tidal modulation, and predicts clouds at altitudes higher than have been observed to date.

Observations show that the lower thermosphere of Mars ($\sim$100--140 km) is up to 40 K colder than the current general circulation models (GCMs) can reproduce. Possible candidates for physical processes missing in the models are larger abundances of atomic oxygen facilitating stronger CO$_2$ radiative cooling, and thermal effects of gravity waves. Using two state-of-the-art Martian GCMs, the Laboratoire de M\'et\'eorologie Dynamique and Max Planck Institute models that self-consistently cover the atmosphere from the surface to the thermosphere, these physical mechanisms are investigated. Simulations demonstrate that the CO$_2$ radiative cooling with a sufficiently large atomic oxygen abundance, and the gravity wave-induced cooling can alone result in up to 40 K colder temperature in the lower thermosphere. Accounting for both mechanisms produce stronger cooling at high latitudes. However, radiative cooling effects peak above the mesopause, while gravity wave cooling rates continuously increase with height. Although both mechanisms act simultaneously, these peculiarities could help to further quantify their relative contributions from future observations.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

This paper presents a contemporary review of vertical coupling in the atmosphere and ionosphere system induced by internal waves of lower atmospheric origin. Atmospheric waves are primarily generated by meteorological processes, possess a broad range of spatial and temporal scales, and can propagate to the upper atmosphere. A brief summary of internal wave theory is given, focusing on gravity waves, solar tides, planetary Rossby and Kelvin waves. Observations of wave signatures in the upper atmosphere, their relationship with the direct propagation of waves into the upper atmosphere, dynamical and thermal impacts as well as concepts, approaches, and numerical modeling techniques are outlined. Recent progress in studies of sudden stratospheric warming and upper atmospheric variability are discussed in the context of wave-induced vertical coupling between the lower and upper atmosphere.

We study a stochastic variability of SS433 in the $10^{-4} - 5\times 10^{-2}$ Hz frequency range based on RXTE data, and on simultaneous observations with RXTE and optical telescopes. We find that the cross-correlation functions and power spectra depend drastically on the precession phase of the supercritical accretion disc. When the wind funnel of the disc is maximally open to the observer, a flat part emerges in the power spectrum; a break is observed at the frequency $1.7\times10^{-3}$ Hz, with a power-law index $\beta \approx 1.67$ at higher frequencies. The soft emission forming mostly in the jets, lags behind the hard and optical emission. When the observer does not see the funnel and jets (the `edge-on' disc), the power spectrum is described by a single power-law with $\beta \approx 1.34$ and no correlations between X-ray ranges are detected. We investigated two mechanisms to explain the observed variability at the open disc phase, 1) reflection of radiation at the funnel wall (X-rays and optical) and 2) the gas cooling in the jets (X-rays only). The X-ray variability is determined by the contribution of both mechanisms, however the contribution of the jets is much higher. We found that the funnel size is $(2-2.5)\times10^{12}$ cm, and the opening angle is $\vartheta_f\sim 50^\circ$. X-ray jets may consist of three fractions with different densities: $8\times10^{13}$, $3\times10^{13}$ and $5\times10^{11}$ cm$^{-3}$, with most of the jet's mass falling within the latter fraction. We suppose that revealed flat part in the power spectrum may be related to an abrupt change in the disc structure and viscous time-scale at the spherization radius, because the accretion disc becomes thick at this radius, $h/r \sim 1$. The extent of the flat spectrum depends on the variation of viscosity at the spherization radius.

The aim of the present paper is to investigate a possible contribution of the rotation-powered pulsars and pulsar wind nebulae to the population of ultraluminous X-ray sources (ULXs). We first develop an analytical model for the evolution of the distribution function of pulsars over the spin period and find both the steady-state and the time-dependent solutions. Using the recent results on the X-ray efficiency dependence on pulsar characteristic age, we then compute the X-ray luminosity function (XLF) of rotation-powered pulsars. In a general case it has a broken power-law shape with a high luminosity cutoff, which depends on the distributions of the birth spin period and the magnetic field. Using the observed XLF of sources in the nearby galaxies and the condition that the pulsar XLF does not exceed that, we find the allowed region for the parameters describing the birth period distribution. We find that the mean pulsar period should be greater than 10-40 ms. These results are consistent with the constraints obtained from the X-ray luminosity of core-collapse supernovae. We estimate that the contribution of the rotation-powered pulsars to the ULX population is at a level exceeding 3 per cent. For a wide birth period distribution, this fraction grows with luminosity and above 10E40 erg/s pulsars can dominate the ULX population.

The observed X-ray luminosity of SS 433 is ~10^36 erg/s, it is known that all the radiation is formed in the famous SS 433 jets. The bolometric luminosity of SS 433 is ~10^40 erg/s, and originally the luminosity must be realized in X-rays. The original radiation is probably thermalized in the supercritical accretion disk wind, however the missing more than four orders of magnitude is surprising. We have analysed the XMM-Newton spectra of SS 433 using a model of adiabatically and radiatively cooling X-ray jets. The multi-temperature thermal jet model reproduces very well the strongest observed emission lines, but it can not reproduce the continuum radiation and some spectral features. We have found a notable contribution of ionized reflection to the spectrum in the energy range from 3 to 12 keV. The reflected spectrum is an evidence of the supercritical disk funnel, where the illuminating radiation comes from deeper funnel regions, to be further reflected in the outer visible funnel walls. The illuminating spectrum is similar to that observed in ULXs, its luminosity has to be no less than ~10^39 erg/s. A soft excess has been detected, that does not depend on the thermal jet model details. It may be represented as a BB with a temperature of ~0.1 keV and luminosity of ~3*10^37 erg/s. The soft spectral component has about the same parameters as those found in ULXs.

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

We present results of optical spectroscopic observations of the mass donor star in SS 433 with Subaru and Gemini, with an aim to best constrain the mass of the compact object. Subaru/FOCAS observations were performed on October 6-8 and 10, 2007, covering the orbital phase of phi=0.96-0.26. We first calculate cross correlation function of these spectra with that of the reference star HD 9233 in the wavelength range of 4740-4840 Angstrom. This region is selected to avoid 'strong' absorption lines accompanied with contaminating emission components. The same analysis is applied to archive data of Gemini/GMOS taken at phi=0.84-0.30 by Hillwig & Gies (2008). From the Subaru and Gemini CCF results, the amplitude of radial velocity curve of the donor star is determined to be 58.3+/-3.8 km s-1 with a systemic velocity of 59.2+/-2.5 km s-1. Together with the radial velocity curve of the compact object, we derive the mass of the donor star and compact object to be M_O=12.4+/-1.9 M_sun and M_X=4.3+/-0.6 M_sun, respectively. We conclude, however, that these values should be taken as upper limits. From the analysis of the averaged absorption line profiles of strong lines and weak lines observed with Subaru, we find evidence for heating effects from the compact object. Using a simple model, we find that the true radial velocity amplitude of the donor star could be as low as 40+/-5 km s-1 in order to produce the observed absorption-line profiles. Taking into account the heating of the donor star may lower the derived masses to M_O=10.4 +2.3/-1.9 M_sun and M_X=2.5 +0.7/-0.6 M_sun. Our final constraint, 1.9 M_sun< M_X <4.9 M_sun, indicates that the compact object in SS 433 is most likely a low mass black hole, although the possibility of a massive neutron star cannot be firmly excluded.

We have analysed the XMM-Newton spectra of SS 433 using a standard model of adiabatically and radiatively cooling X-ray jets. The multi-temperature thermal jet model reproduces well the strongest observed emission line fluxes. The thermal model alone can not reproduce the continuum radiation in the XMM spectral range, the fluorescent iron line and some broad spectral features. Using the thermal jet-plus-reflection model, we find a notable contribution of ionized reflection to the spectrum in the energy range from ~3 to 12 keV. The reflecting surface is highly ionized (xi ~300), the illuminating radiation spectrum is flat and probably strongly absorbed. We conclude that the reflected spectrum is an evidence of the supercritical disc funnel, where the illuminating radiation comes from deeper funnel regions, to be further reflected in the outer visible funnel walls (r >~2*10^11 cm). We have not found any evidences of reflection in the soft 0.8-2 keV energy range, instead, a soft excess is detected, that does not depend on the thermal jet model details. This soft component might prove to be the direct radiation of the visible funnel wall. It is represented here either as black body radiation with the temperature of ~0.1 keV and luminosity of L_bb ~3*10^37 erg/s, or with a multicolour funnel (MCF) model. The soft spectral component has about the same parameters as those found in ULXs.

We present the first results of time-resolved spectroscopic observations of little-studied dwarf nova BZ Ursae Majoris in quiescence.