Constraining the Atmospheric Composition of the Day-Night Terminators of HD 189733b : Atmospheric Retrieval with Aerosols

A number of observations have shown that Rayleigh scattering by aerosols dominates the transmission spectrum of HD 189733b at wavelengths shortward of 1 $\mu$m. In this study, we retrieve a range of aerosol distributions consistent with transmission spectroscopy between 0.3-24 $\mu$m that were recently re-analyzed by Pont et al. (2013). To constrain the particle size and the optical depth of the aerosol layer, we investigate the degeneracies between aerosol composition, temperature, planetary radius, and molecular abundances that prevent unique solutions for transit spectroscopy. Assuming that the aerosol is composed of MgSiO$_3$, we suggest that a vertically uniform aerosol layer over all pressures with a monodisperse particle size smaller than about 0.1 $\mu$m and an optical depth in the range 0.002-0.02 at 1 $\mu$m provides statistically meaningful solutions for the day/night terminator regions of HD 189733b. Generally, we find that a uniform aerosol layer provide adequate fits to the data if the optical depth is less than 0.1 and the particle size is smaller than 0.1 $\mu$m, irrespective of the atmospheric temperature, planetary radius, aerosol composition, and gaseous molecules. Strong constraints on the aerosol properties are provided by spectra at wavelengths shortward of 1 $\mu$m as well as longward of 8 $\mu$m, if the aerosol material has absorption features in this region. We show that these are the optimal wavelengths for quantifying the effects of aerosols, which may guide the design of future space observations. The present investigation indicates that the current data offer sufficient information to constrain some of the aerosol properties of HD189733b, but the chemistry in the terminator regions remains uncertain.Comment: Transferred to ApJ and accepted. 11 pages, 10 figures, 1 tabl

Seasonal Variability of Saturn's Tropospheric Temperatures, Winds and Para-H2_2 from Cassini Far-IR Spectroscopy

Far-IR 16-1000 $\mu$m spectra of Saturn's hydrogen-helium continuum measured by Cassini's Composite Infrared Spectrometer (CIRS) are inverted to construct a near-continuous record of upper tropospheric (70-700 mbar) temperatures and para-H$_2$ fraction as a function of latitude, pressure and time for a third of a Saturnian year (2004-2014, from northern winter to northern spring). The thermal field reveals evidence of reversing summertime asymmetries superimposed onto the belt/zone structure. The temperature structure that is almost symmetric about the equator by 2014, with seasonal lag times that increase with depth and are qualitatively consistent with radiative climate models. Localised heating of the tropospheric hazes (100-250 mbar) create a distinct perturbation to the temperature profile that shifts in magnitude and location, declining in the autumn hemisphere and growing in the spring. Changes in the para-H$_2$ ($f_p$) distribution are subtle, with a 0.02-0.03 rise over the spring hemisphere (200-500 mbar) perturbed by (i) low-$f_p$ air advected by both the springtime storm of 2010 and equatorial upwelling; and (ii) subsidence of high-$f_p$ air at northern high latitudes, responsible for a developing north-south asymmetry in $f_p$. Conversely, the shifting asymmetry in the para-H$_2$ disequilibrium primarily reflects the changing temperature structure (and the equilibrium distribution of $f_p$), rather than actual changes in $f_p$ induced by chemical conversion or transport. CIRS results interpolated to the same point in the seasonal cycle as re-analysed Voyager-1 observations show qualitative consistency, with the exception of the tropical tropopause near the equatorial zones and belts, where downward propagation of a cool temperature anomaly associated with Saturn's stratospheric oscillation could potentially perturb tropopause temperatures, para-H$_2$ and winds. [ABRIDGED]Comment: Preprint accepted for publication in Icarus, 29 pages, 18 figure

On the potential of the EChO mission to characterise gas giant atmospheres

Space telescopes such as EChO (Exoplanet Characterisation Observatory) and JWST (James Webb Space Telescope) will be important for the future study of extrasolar planet atmospheres. Both of these missions are capable of performing high sensitivity spectroscopic measurements at moderate resolutions in the visible and infrared, which will allow the characterisation of atmospheric properties using primary and secondary transit spectroscopy. We use the NEMESIS radiative transfer and retrieval tool (Irwin et al. 2008, Lee et al. 2012) to explore the potential of the proposed EChO mission to solve the retrieval problem for a range of H2-He planets orbiting different stars. We find that EChO should be capable of retrieving temperature structure to ~200 K precision and detecting H2O, CO2 and CH4 from a single eclipse measurement for a hot Jupiter orbiting a Sun-like star and a hot Neptune orbiting an M star, also providing upper limits on CO and NH3. We provide a table of retrieval precisions for these quantities in each test case. We expect around 30 Jupiter-sized planets to be observable by EChO; hot Neptunes orbiting M dwarfs are rarer, but we anticipate observations of at least one similar planet.Comment: 22 pages, 30 figures, 4 tables. Accepted for publication in MNRA

New insights on Saturn's formation from its nitrogen isotopic composition

The recent derivation of a lower limit for the $^{14}$N/$^{15}$N ratio in Saturn's ammonia, which is found to be consistent with the Jovian value, prompted us to revise models of Saturn's formation using as constraints the supersolar abundances of heavy elements measured in its atmosphere. Here we find that it is possible to account for both Saturn's chemical and isotopic compositions if one assumes the formation of its building blocks at $\sim$45 K in the protosolar nebula, provided that the O abundance was $\sim$2.6 times protosolar in its feeding zone. To do so, we used a statistical thermodynamic model to investigate the composition of the clathrate phase that formed during the cooling of the protosolar nebula and from which the building blocks of Saturn were agglomerated. We find that Saturn's O/H is at least $\sim$34.9 times protosolar and that the corresponding mass of heavy elements ($\sim$43.1 \Mearth) is within the range predicted by semi-convective interior models.Comment: Accepted for publication in Astrophysical Journal Letter

Optimal Estimation Retrievals of the Atmospheric Structure and Composition of HD 189733b from Secondary Eclipse Spectroscopy

Recent spectroscopic observations of transiting hot Jupiters have permitted the derivation of the thermal structure and molecular abundances of H2O, CO2, CO, and CH4 in these extreme atmospheres. Here, for the first time, we apply the technique of optimal estimation to determine the thermal structure and composition of an exoplanet by solving the inverse problem. The development of a suite of radiative transfer and retrieval tools for exoplanet atmospheres is described, building upon a retrieval algorithm which is extensively used in the study of our own solar system. First, we discuss the plausibility of detection of different molecules in the dayside atmosphere of HD 189733b and the best-fitting spectrum retrieved from all publicly available sets of secondary eclipse observations between 1.45 and 24 {\mu}m. Additionally, we use contribution functions to assess the vertical sensitivity of the emission spectrum to temperatures and molecular composition. Over the altitudes probed by the contribution functions, the retrieved thermal structure shows an isothermal upper atmosphere overlying a deeper adiabatic layer (temperature decreasing with altitude), which is consistent with previously-reported dynamical and observational results. The formal uncertainties on retrieved parameters are estimated conservatively using an analysis of the cross-correlation functions and the degeneracy between different atmospheric properties. The formal solution of the inverse problem suggests that the uncertainties on retrieved parameters are larger than suggested in previous studies, and that the presence of CO and CH4 is only marginally supported by the available data. Nevertheless, by including as broad a wavelength range as possible in the retrieval, we demonstrate that available spectra of HD 189733b can constrain a family of potential solutions for the atmospheric structure.Comment: 13 pages, 10 figures, 1 table, Accepted for publication in MNRA

Exoplanet atmospheres with EChO: spectral retrievals using EChOSim

We demonstrate the effectiveness of the Exoplanet Characterisation Observatory mission concept for constraining the atmospheric properties of hot and warm gas giants and super Earths. Synthetic primary and secondary transit spectra for a range of planets are passed through EChOSim (Waldmann & Pascale 2014) to obtain the expected level of noise for different observational scenarios; these are then used as inputs for the NEMESIS atmospheric retrieval code and the retrieved atmospheric properties (temperature structure, composition and cloud properties) compared with the known input values, following the method of Barstow et al. (2013a). To correctly retrieve the temperature structure and composition of the atmosphere to within 2 {\sigma}, we find that we require: a single transit or eclipse of a hot Jupiter orbiting a sun-like (G2) star at 35 pc to constrain the terminator and dayside atmospheres; 20 transits or eclipses of a warm Jupiter orbiting a similar star; 10 transits/eclipses of a hot Neptune orbiting an M dwarf at 6 pc; and 30 transits or eclipses of a GJ1214b-like planet.Comment: 13 pages, 15 figures, 1 table. Accepted by Experimental Astronomy. The final publication will shortly be available at Springer via http://dx.doi.org/10.1007/s10686-014-9397-

Time variability of Neptune's horizontal and vertical cloud structure revealed by VLT/SINFONI and Gemini/NIFS from 2009 to 2013

New observations of Neptune's clouds in the near infrared were acquired in October 2013 with SINFONI on ESO's Very Large Telescope (VLT) in Chile. SINFONI is an Integral Field Unit spectrometer returning a 64 × 64 pixel image with 2048 wavelengths. Image cubes in the J-band (1.09-1.41 μm) and H-band (1.43-1.87 μm) were obtained at spatial resolutions of 0.1″and 0.025″per pixel, while SINFONI's adaptive optics provided an effective resolution of approximately 0.1″. Image cubes were obtained at the start and end of three successive nights to monitor the temporal development of discrete clouds both at short timescales (i.e. during a single night) as well as over the longer period of the three-day observing run. These observations were compared with similar H-band observations obtained in September 2009 with the NIFS Integral Field Unit spectrometer on the Gemini-North telescope in Hawaii, previously reported by Irwin et al. (2011) [Icarus, 216, 141-158], and previously unreported Gemini/NIFS observations at lower spatial resolution made in 2011.We find both similarities and differences between these observations, spaced over four years. The same overall cloud structure is seen with high, bright clouds visible at mid-latitudes (30-40°N,S), with slightly lower clouds observed at lower latitudes, together with small discrete clouds seen circling the pole at a latitude of approximately 60°S. However, while discrete clouds were visible at this latitude at both the main cloud deck level (at 2-3 bar) and in the upper troposphere (100-500 mb) in 2009, no distinct deep (2-3 bar), discrete circumpolar clouds were visible in 2013, although some deep clouds were seen at the southern edge of the main cloud belt at 30-40°S, which have not been observed before. The nature of the deep sub-polar discrete clouds observed in 2009 is intriguing. While it is possible that in 2013 these deeper clouds were masked by faster moving, overlying features, we consider that it is unlikely that this should have happened in 2013, but not in 2009 when the upper-cloud activity was generally similar. Meanwhile, the deep clouds seen at the southern edge of the main cloud belt at 30-40°S in 2013, should also have been detectable in 2009, but were not seen. Hence, these observations may have detected a real temporal variation in the occurrence of Neptune's deep clouds, pointing to underlying variability in the convective activity at the pressure of the main cloud deck at 2-3 bar near Neptune's south pole and also in the main observable cloud belt at 30-40°S.</p