SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

We present the design of a point-and-shoot non-imaging full-Stokes spectropolarimeter dedicated to detecting life on Earth from an orbiting platform like the ISS. We specifically aim to map circular polarization in the spectral features of chlorophyll and other biopigments for our planet as a whole. These non-zero circular polarization signatures are caused by homochirality of the molecular and supramolecular configurations of organic matter, and are considered the most unambiguous biomarker. To achieve a fully solid-state snapshot design, we implement a novel spatial modulation that completely separates the circular and linear polarization channels. The polarization modulator consists of a patterned liquid-crystal quarter-wave plate inside the spectrograph slit, which also constitutes the first optical element of the instrument. This configuration eliminates cross-talk between linear and circular polarization, which is crucial because linear polarization signals are generally much stronger than the circular polarization signals. This leads to a quite unorthodox optical concept for the spectrograph, in which the object and the pupil are switched. We discuss the general design requirements and trade-offs of LSDpol (Life Signature Detection polarimeter), a prototype instrument that is currently under development.

One of the most exciting scientific challenges is to detect Earth-like planets in the habitable zones of other stars in the galaxy and search for evidence of life. The ability to observe and characterise dozens of potentially Earth-like planets now lies within the realm of possibility due to rapid advances in key space and imaging technologies. The associated challenge of directly imaging very faint planets in orbit around nearby very bright stars is now well understood, with the key instrumentation also being perfected and developed. Such advances will allow us to develop large transformative telescopes, covering a broad UV-optical-IR spectral range, which can carry out the detailed research programmes designed to answer the questions we wish to answer: Carry out high contrast imaging surveys of nearby stars to search for planets within their habitable zones. Characterise the planets detected to determine masses and radii from photometric measurements. Through spectroscopic studies of their atmospheres and surfaces, search for habitability indicators and for signs of an environment that has been modified by the presence of life. Active studies of potential missions have been underway for a number of years. The latest of these is the Large UV Optical IR space telescope (LUVOIR), one of four flagship mission studies commissioned by NASA in support of the 2020 US Decadal Survey. LUVOIR, if selected, will be of interest to a wide scientific community and will be the only telescope capable of searching for and characterizing a sufficient number of exoEarths to provide a meaningful answer to the question - Are we alone?. This paper is a submission to the European Space Agency Voyage 2050 call for white papers outlining the case for an ESA contribution to a Large UVOIR telescope.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

An exoplanet's habitability will depend strongly on the presence of liquid water. Flux and/or polarization measurements of starlight that is reflected by exoplanets could help to identify exo-oceans. We investigate which broadband spectral features in flux and polarization phase functions of reflected starlight uniquely identify exo-oceans. We compute total fluxes F and polarized fluxes Q of starlight reflected by cloud-free and (partly) cloudy exoplanets, for wavelengths from 350 to 865 nm. The ocean surface has waves composed of Fresnel reflecting wave facets and whitecaps, and scattering within the water body is included. Total flux F, polarized flux Q, and degree of polarization P of ocean planets change color from blue, through white, to red at phase angles alpha ranging from 134-108 deg for F, and from 123-157 deg for Q, with cloud coverage fraction fc increasing from 0.0 to 1.0 for F, and to 0.98 for Q. The color change in P only occurs for fc ranging from 0.03-0.98, with the color crossing angle alpha ranging from 88-161 deg. The total flux F of a cloudy, zero surface albedo planet can also change color, and for fc=0.0, an ocean planet's F will not change color for surface pressures ps > 8 bars. Polarized flux Q of a zero surface albedo planet does not change color for any fc. The color change of P of starlight reflected by an exoplanet, from blue, through white, to red with increasing alpha above 88 deg, appears to identify a (partly) cloudy exo-ocean. The color change of polarized flux Q with increasing alpha above 123 deg appears to uniquely identify an exo-ocean, independent of surface pressure or cloud fraction. At the color changing phase angle, the angular distance between a star and its planet is much larger than at the phase angle where the glint appears in reflected light. The color change in polarization thus offers better prospects for detecting exo-oceans.

The polarization state of starlight reflected by a planetary atmosphere uniquely reveals coverage, particle size, and composition of aerosols as well as changing cloud patterns. It is not possible to obtain a comparable level of detailed from flux-only observations. Furthermore, polarization observations can probe the atmosphere of planets independently of the orbital geometry (i.e., transiting and non-transiting planets). We show that a high-resolution spectropolarimeter with a broad wavelength coverage, particularly if attached to a large space telescope, would enable simultaneous study of the polarimetric exoplanet properties of the continuum and to look for and characterize the polarimetric signal due to scattering from single molecules.

We present a comprehensive set of spectropolarimetric observations of Earthshine as obtained by FORS2 at the VLT for phase angles from 50degree to 135degree (Sun-Earth-Moon angle), covering a spectral range from 430nm to 920nm. The degree of polarization in BVRI passbands, the differential polarization vegetation index, and the equivalent width of the O2A polarization band around 760nm are determined with absolute errors around 0.1 percent in the degree of polarization. Earthshine polarization spectra are corrected for the effect of depolarization introduced by backscattering on the lunar surface, introducing systematic errors of the order of 1 percent in the degree of polarization. Distinct viewing sceneries such as observing the Atlantic or Pacific side in Earthshine yield statistically different phase curves. The equivalent width defined for the O2A band polarization is found to vary from -5nm to +2nm. A differential polarized vegetation index is introduced and reveals a larger vegetation signal for those viewing sceneries that contain larger fractions of vegetated surface areas. We corroborate the observed correlations with theoretical models from the literature, and conclude that the Vegetation Red Edge(VRE) is a robust and sensitive signature in polarization spectra of planet Earth. The overall behaviour of polarization of planet Earth in the continuum and in the O2A band can be explained by existing models. Biosignatures such as the O2A band and the VRE are detectable in Earthshine polarization with a high degree of significance and sensitivity. An in-depth understanding of Earthshines temporal and spectral variability requires improved models of Earths biosphere, as a prerequisite to interpret possible detections of polarised biosignatures in earthlike exoplanets in the future.

Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere.

Context: Detecting moons around exoplanets is a major goal of current and future observatories. Moons are suspected to influence rocky exoplanet habitability, and gaseous exoplanets in stellar habitable zones could harbour abundant and diverse moons to target in the search for extraterrestrial habitats. Exomoons will contribute to exoplanetary signals but are virtually undetectable with current methods. Aims: We identify and analyse traces of exomoons in the temporal variation of total and polarised fluxes of starlight reflected by an Earth-like exoplanet and its spatially unresolved moon across all phase angles, with both orbits viewed in an edge-on geometry. Methods: We compute the total and linearly polarised fluxes, and the degree of linear polarization P of starlight that is reflected by the exoplanet with its moon along their orbits, accounting for the temporal variation of the visibility of the planetary and lunar disks, and including effects of mutual transits and mutual eclipses. Our computations pertain to a wavelength of 450 nm. Results: Total flux F shows regular dips due to planetary and lunar transits and eclipses. Polarization P shows regular peaks due to planetary transits and lunar eclipses, and P can increase and/or slightly decrease during lunar transits and planetary eclipses. Changes in F and P will depend on the radii of the planet and moon, on their reflective properties, and their orbits, and are about one magnitude smaller than the smooth background signals. The typical duration of a transit or an eclipse is a few hours. Conclusions: Traces of an exomoon due to planetary and lunar transits and eclipses show up in F and P of sunlight reflected by planet-moon systems and could be searched for in exoplanet flux and/or polarisation phase functions.

The circular polarization of light that planets reflect is often neglected because it is very small compared to the linear polarization. It could, however, provide information on a planet's atmosphere and surface, and on the presence of life, because homochiral molecules that are the building blocks of life on Earth are known to reflect circularly polarized light. We compute $P_c$, the degree of circular polarization, for light that is reflected by rocky (exo)planets with liquid water or sulfuric acid solution clouds, both spatially resolved across the planetary disk and, for planets with patchy clouds, integrated across the planetary disk, for various planetary phase angles $\alpha$. The optical thickness and vertical distribution of the atmospheric gas and clouds, the size parameter and refractive index of the cloud particles, and $\alpha$ all influence $P_c$. Spatially resolved, $P_c$ varies between $\pm 0.20\%$ (the sign indicates the polarization direction). Only for small gas optical thicknesses above the clouds do significant sign changes (related to cloud particle properties) across the planets' hemispheres occur. For patchy clouds, the disk--integrated $P_c$ is typically smaller than $\pm 0.025\%$, with maximums for $\alpha$ between $40^\circ$ and $70^\circ$, and $120^\circ$ to $140^\circ$. As expected, the disk--integrated $P_c$ is virtually zero at $\alpha=0^\circ$ and 180$^\circ$. The disk--integrated $P_c$ is also very small at $\alpha \approx 100^\circ$. Measuring circular polarization signals appears to be challenging with current technology. The small atmospheric circular polarization signal could, however, allow the detection of circular polarization due to homochiral molecules. Confirmation of the detectability of such signals requires better knowledge of the strength of circular polarization signals of biological sources.

PyMieDAP (the Python Mie Doubling-Adding Programme) is a Python--based tool for computing the total, linearly, and circularly polarized fluxes of incident unpolarized sun- or starlight that is reflected by, respectively, Solar System planets or moons, or exoplanets at a range of wavelengths. The radiative transfer computations are based on an adding--doubling Fortran algorithm and fully include polarization for all orders of scattering. The model (exo)planets are described by a model atmosphere composed of a stack of homogeneous layers containing gas and/or aerosol and/or cloud particles bounded below by an isotropically, depolarizing surface (that is optionally black). The reflected light can be computed spatially--resolved and/or disk--integrated. Spatially--resolved signals are mostly representative for observations of Solar System planets (or moons), while disk--integrated signals are mostly representative for exoplanet observations. PyMieDAP is modular and flexible, and allows users to adapt and optimize the code according to their needs. PyMieDAP keeps options open for connections with external programs and for future additions and extensions. In this paper, we describe the radiative transfer algorithm that PyMieDAP is based on and the code's principal functionalities. And we provide benchmark results of PyMieDAP that can be used for testing its installation and for comparison with other codes. PyMieDAP is available online under the GNU GPL license at http://gitlab.com/loic.cg.rossi/pymiedap

Young giant exoplanets emit infrared radiation that can be linearly polarized up to several percent. This linear polarization can trace: 1) the presence of atmospheric cloud and haze layers, 2) spatial structure, e.g. cloud bands and rotational flattening, 3) the spin axis orientation and 4) particle sizes and cloud top pressure. We introduce a novel high-contrast imaging scheme that combines angular differential imaging (ADI) and accurate near-infrared polarimetry to characterize self-luminous giant exoplanets. We implemented this technique at VLT/SPHERE-IRDIS and developed the corresponding observing strategies, the polarization calibration and the data-reduction approaches. By combining ADI and polarimetry we can characterize planets that can be directly imaged with a very high signal-to-noise ratio. We use the IRDIS pupil-tracking mode and combine ADI and principal component analysis to reduce speckle noise. We take advantage of IRDIS' dual-beam polarimetric mode to eliminate differential effects that severely limit the polarimetric sensitivity (flat-fielding errors, differential aberrations and seeing), and thus further suppress speckle noise. To correct for instrumental polarization effects, we apply a detailed Mueller matrix model that describes the telescope and instrument and that has an absolute polarimetric accuracy $\leq0.1\%$. Using this technique we have observed the planets of HR 8799 and the (sub-stellar) companion PZ Tel B. Unfortunately, we do not detect a polarization signal in a first analysis. We estimate preliminary $1\sigma$ upper limits on the degree of linear polarization of $\sim1\%$ and $\sim0.1\%$ for the planets of HR 8799 and PZ Tel B, respectively. The achieved sub-percent sensitivity and accuracy show that our technique has great promise for characterizing exoplanets through direct-imaging polarimetry.

Context. Clouds have already been detected in exoplanetary atmospheres. They play crucial roles in a planet's atmosphere and climate and can also create ambiguities in the determination of atmospheric parameters such as trace gas mixing ratios. Knowledge of cloud properties is required when assessing the habitability of a planet. Aims. We aim to show that various types of cloud cover such as polar cusps, subsolar clouds, and patchy clouds on Earth-like exoplanets can be distinguished from each other using the polarization and flux of light that is reflected by the planet. Methods. We have computed the flux and polarization of reflected starlight for different types of (liquid water) cloud covers on Earth-like model planets using the adding-doubling method, that fully includes multiple scattering and polarization. Variations in cloud-top altitudes and planet-wide cloud cover percentages were taken into account. Results. We find that the different types of cloud cover (polar cusps, subsolar clouds, and patchy clouds) can be distinguished from each other and that the percentage of cloud cover can be estimated within 10%. Conclusions. Using our proposed observational strategy, one should be able to determine basic orbital parameters of a planet such as orbital inclination and estimate cloud coverage with reduced ambiguities from the planet's polarization signals along its orbit.

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmosphere and surface can be derived by analyzing light of the parent star reflected by the planet. We investigate the influence of the surface albedo $A_{\rm s}$, the optical thickness $b_{\rm cloud}$ and altitude of water clouds, and the mixing ratio $\eta$ of biosignature O$_2$ on the strength of the O$_2$ A-band (around 760 nm) in flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ($\eta$ < 0.4) will yield moderately deep bands in flux and moderate to small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on their properties such as optical thickness, top altitude, particle phase, coverage fraction, horizontal distribution. Depending on the surface albedo, and cloud properties, different O$_2$ mixing ratios $\eta$ can give similar absorption band depths in flux and band strengths in polarization, in particular if the clouds have moderate to high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, in particular not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O$_2$ or any other absorbing gas.

The Hera Saturn entry probe mission is proposed as an M--class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-Relay Spacecraft is powered by solar panels and batteries. We anticipate two major subsystems to be supplied by the United States, either by direct procurement by ESA or by contribution from NASA: the solar electric power system (including solar arrays and the power management and distribution system), and the probe entry system (including the thermal protection shield and aeroshell). Hera is designed to perform in situ measurements of the chemical and isotopic compositions as well as the dynamics of Saturn's atmosphere using a single probe, with the goal of improving our understanding of the origin, formation, and evolution of Saturn, the giant planets and their satellite systems, with extrapolation to extrasolar planets. Hera's aim is to probe well into the cloud-forming region of the troposphere, below the region accessible to remote sensing, to the locations where certain cosmogenically abundant species are expected to be well mixed. By leading to an improved understanding of the processes by which giant planets formed, including the composition and properties of the local solar nebula at the time and location of giant planet formation, Hera will extend the legacy of the Galileo and Cassini missions by further addressing the creation, formation, and chemical, dynamical, and thermal evolution of the giant planets, the entire solar system including Earth and the other terrestrial planets, and formation of other planetary systems.

We present numerically calculated, disk--integrated, spectropolarimetric signals of starlight that is reflected by vertically and horizontally inhomogeneous gaseous exoplanets. We include various spatial features that are present on Solar System's gaseous planets: belts and zones, cyclonic spots, and polar hazes, to test whether such features leave traces in the disk--integrated flux and polarization signals. Broadband flux and polarization signals of starlight that is reflected by gaseous exoplanets are calculated using an efficient, adding--doubling radiative transfer code, that fully includes single and multiple scattering and polarization. The planetary model atmospheres are vertically inhomogeneous and can be horizontally inhomogeneous, and contain gas molecules and/or cloud and/or aerosol particles.

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

Looking for the primary rainbow in starlight that is reflected by exoplanets appears to be a promising method to search for liquid water clouds in exoplanetary atmospheres. Ice water clouds, that consist of water crystals instead of water droplets, could potentially mask the rainbow feature in the planetary signal by covering liquid water clouds. Here, we investigate the strength of the rainbow feature for exoplanets that have liquid and icy water clouds in their atmosphere, and calculate the rainbow feature for a realistic cloud coverage of Earth. We calculate flux and polarization signals of starlight that is reflected by horizontally and vertically inhomogeneous Earth--like exoplanets, covered by patchy clouds consisting of liquid water droplets or water ice crystals. The planetary surfaces are black. On a planet with a significant coverage of liquid water clouds only, the total flux signal shows a weak rainbow feature. Any coverage of the liquid water clouds by ice clouds, however, dampens the rainbow feature in the total flux, and thus the discovery of liquid water in the atmosphere. On the other hand, detecting the primary rainbow in the polarization signal of exoplanets appears to be a powerful tool for detecting liquid water in exoplanetary atmospheres, even when these clouds are partially covered by ice clouds. In particular, liquid water clouds covering as little as 10%-20% of the planetary surface, with more than half of these covered by ice clouds, still create a polarized rainbow feature in the planetary signal. Indeed, calculations of flux and polarization signals of an exoplanet with a realistic Earth--like cloud coverage, show a strong polarized rainbow feature.

We present modeled flux and linear polarisation signals of starlight that is reflected by spatially unresolved, horizontally inhomogeneous planets and discuss the effects of including horizontal inhomogeneities on the flux and polarisation signals of Earth-like exoplanets. Our code is based on an efficient adding--doubling algorithm, which fully includes multiple scattering by gases and aerosol/cloud particles. We divide a model planet into pixels that are small enough for the local properties of the atmosphere and surface (if present) to be horizontally homogeneous. Given a planetary phase angle, we sum up the reflected total and linearly polarised fluxes across the illuminated and visible part of the planetary disk, taking care to properly rotate the polarized flux vectors towards the same reference plane. We compared flux and polarisation signals of simple horizontally inhomogeneous model planets against results of the weighted sum approximation, in which signals of horizontally homogeneous planets are combined. Apart from cases in which the planet has only a minor inhomogeneity, the signals differ significantly. In particular, the shape of the polarisation phase function appears to be sensitive to the horizontal inhomogeneities. The same holds true for Earth-like model planets with patchy clouds above an ocean and a sandy continent. Our simulations clearly show that horizontal inhomogeneities leave different traces in flux and polarisation signals. Combining flux with polarisation measurements would help retrieving the atmospheric and surface patterns on a planet.

[Abridged] The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the planetary atmosphere, the amount of transmitted flux not only depends on the total optical thickness of the slant light path that is probed, but also on the amount of forward-scattering by the scattering particles. Here, we present results of calculations with a three-dimensional Monte Carlo code that simulates the transmitted flux during occultations or transits. For isotropically scattering particles, like gas molecules, the transmitted flux appears to be well-described by the total atmospheric optical thickness. Strongly forward-scattering particles, however, such as commonly found in atmospheres of Solar System planets, can increase the transmitted flux significantly. For exoplanets, such added flux can decrease the apparent radius of the planet by several scale heights, which is comparable to predicted and measured features in exoplanet transit spectra. We performed detailed calculations for Titan's atmosphere between 2.0 and 2.8 micron and show that haze and gas abundances will be underestimated by about 8% if forward-scattering is ignored in the retrievals. At shorter wavelengths, errors in the gas and haze abundances and in the spectral slope of the haze particles can be several tens of percent, also for other Solar System planetary atmospheres. We also find that the contribution of forward-scattering can be fairly well described by modelling the atmosphere as a plane-parallel slab.

Context. High-contrast imaging is currently the only available technique for the study of the thermodynamical and compositional properties of exoplanets in long-period orbits. The SPICES project is a coronagraphic space telescope dedicated to the spectro-polarimetric analysis of gaseous and icy giant planets as well as super-Earths at visible wavelengths. So far, studies for high-contrast imaging instruments have mainly focused on technical feasibility because of the challenging planet/star flux ratio of 10-8-10-10 required at short separations (200 mas or so) to image cold exoplanets. However, the analysis of planet atmospheric/surface properties has remained largely unexplored. Aims. The aim of this paper is to determine which planetary properties SPICES or an equivalent direct imaging mission can measure, considering realistic reflected planet spectra and instrument limitation. Methods. We use numerical simulations of the SPICES instrument concept and theoretical planet spectra to carry out this performance study. Results. We find that the characterization of the main planetary properties (identification of molecules, effect of metallicity, presence of clouds and type of surfaces) would require a median signal-to-noise ratio of at least 30. In the case of a solar-type star \leq 10 pc, SPICES will be able to study Jupiters and Neptunes up to ~5 and ~2 AU respectively. It would also analyze cloud and surface coverage of super-Earths of radius 2.5 RE at 1 AU. Finally, we determine the potential targets in terms of planet separation, radius and distance for several stellar types. For a Sun analog, we show that SPICES could characterize Jupiters (M \geq 30 ME) as small as 0.5 Jupiter radii at ~2 AU up to 10 pc, and super-Earths at 1-2 AU for the handful of stars that exist within 4-5 pc. Potentially, SPICES could perform analysis of a hypothetical Earth-size planet around alpha Cen A and B.

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System.

The detections of small, rocky exoplanets have surged in recent years and will likely continue to do so. To know whether a rocky exoplanet is habitable, we have to characterise its atmosphere and surface. A promising characterisation method for rocky exoplanets is direct detection using spectropolarimetry. This method will be based on single pixel signals, because spatially resolving exoplanets is impossible with current and near-future instruments. Well-tested retrieval algorithms are essential to interpret these single pixel signals in terms of atmospheric composition, cloud and surface coverage. Observations of Earth itself provide the obvious benchmark data for testing such algorithms. The observations should provide signals that are integrated over the Earth's disk, that capture day and night variations, and all phase angles. The Moon is a unique platform from where the Earth can be observed as an exoplanet, undisturbed, all of the time. Here, we present LOUPE, the Lunar Observatory for Unresolved Polarimetry of Earth, a small and robust spectropolarimeter to observe our Earth as an exoplanet.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g. due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1 %, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line-of-sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes.

(abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scattering in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering particles will show deeper absorption features than absorbing (e.g. iron) particles and particle size will also affect the calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can yield retrieved properties that differ significantly from the real atmospheric properties.

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 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

After Earth's origin, our host star, the Sun, was shining 20 to 25 percent less brightly than today. Without greenhouse-like conditions to warm the atmosphere, our early planet would have been an ice ball and life may never have evolved. But life did evolve, which indicates that greenhouse gases must have been present on early Earth to warm the planet. Evidence from the geologic record indicates an abundance of the greenhouse gas CO2. CH4 was probably present as well, and in this regard methanogenic bacteria, which belong to a diverse group of anaerobic procaryotes that ferment CO 2 plus H2 to CH4, may have contributed to modification of the early atmosphere. Molecular oxygen was not present, as is indicated by the study of rocks from that era, which contain iron carbonate rather than iron oxide. Multicellular organisms originated as cells within colonies that became increasingly specialized. The development of photosynthesis allowed the Sun's energy to be harvested directly by life forms. The resultant oxygen accumulated in the atmosphere and formed the ozone layer in the upper atmosphere. Aided by the absorption of harmful UV radiation in the ozone layer, life colonized Earth's surface. Our own planet is a very good example of how life forms modified the atmosphere over the planets' life time. We show that these facts have to be taken into account when we discover and characterize atmospheres of Earth-like exoplanets. If life has originated and evolved on a planet, then it should be expected that a strong co-evolution occurred between life and the atmosphere, the result of which is the planets' climate.

In this outlook we describe what could be the next steps of the direct characterization of habitable exoplanets after first the medium and large mission projects and investigate the benefits of the spectroscopic and direct imaging approaches. We show that after third and fourth generation missions foreseeable for the next 100 years, we will face a very long era before being able to see directly the morphology of extrasolar organisms.

The physical characterization of exoplanets will require to take spectra at several orbital positions. For that purpose, a direct imaging capability is necessary. Direct imaging requires an efficient stellar suppression mechanism, associated with an ultrasmooth telescope. We show that before future large space missions (interferometer, 4-8 m class coronograph, external occulter or Fresnel imager), direct imaging of giant planets and close-by super-Earth are at the cross-road of a high scientific interest and a reasonable feasibility. The scientific interest lies in the fact that super-Earths share common geophysical attributes with Earths. They already begin to be detected by radial velocity (RV) and, together with giant planets, they have a larger area than Earths, making them detectable with a 1.5-2 m class telescope in reflected light. We propose such a (space) telescope be a first step before large direct imaging missions.

The Super-Earth Explorer is an Off-Axis Space Telescope (SEE-COAST) designed for high contrast imaging. Its scientific objective is to make the physico-chemical characterization of exoplanets possibly down to 2 Earth radii >. For that purpose it will analyze the spectral and polarimetric properties of the parent starlight reflected by the planets, in the wavelength range 400-1250 nm

We present measurements of ratios of elements of the scattering matrix of Martian analogue palagonite particles for scattering angles ranging from 3 to 174 degrees and a wavelength of 632.8 nm. To facilitate the use of these measurements in radiative transfer calculations we have devised a method that enables us to obtain, from these measurements, a normalized synthetic scattering matrix covering the complete scattering angle range from 0 to 180 degrees. Our method is based on employing the coefficients of the expansions of scattering matrix elements into generalized spherical functions. The synthetic scattering matrix elements and/or the expansion coefficients obtained in this way, can be used to include multiple scattering by these irregularly shaped particles in (polarized) radiative transfer calculations, such as calculations of sunlight that is scattered in the dusty Martian atmosphere.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

We present results of numerical simulations of the flux (irradiance), F, and the degree of polarization (i.e. the ratio of polarized to total flux), P, of light that is reflected by Earth-like extrasolar planets orbiting solar-type stars, as functions of the wavelength (from 0.3 to 1.0 micron, with 0.001 micron spectral resolution) and as functions of the planetary phase angle. We use different surface coverages for our model planets, including vegetation and a Fresnel reflecting ocean, and clear and cloudy atmospheres. Our adding-doubling radiative transfer algorithm, which fully includes multiple scattering and polarization, handles horizontally homogeneous planets only; we simulate fluxes and polarization of horizontally inhomogeneous planets by weighting results for homogeneous planets. Like the flux, F, the degree of polarization, P, of the reflected starlight is shown to depend strongly on the phase angle, on the composition and structure of the planetary atmosphere, on the reflective properties of the underlying surface, and on the wavelength, in particular in wavelength regions with gaseous absorption bands. The sensitivity of P to a planet's physical properties appears to be different than that of F. Combining flux with polarization observations thus makes for a strong tool for characterizing extrasolar planets. The calculated total and polarized fluxes will be made available through the CDS.