Dynamics of baroclinic jets on extrasolar planets is studied using three-dimensional general circulation models (GCMs) which solve the traditional hydrostatic primitive equations. The focus is on: i ) baroclinic ow and instability on hot-Jupiters; ii ) detailed GCM intercomparison in a commonly used extrasolar planet setup; and, iii ) equatorial superrotation on Earth-like planets. Stability, non-linear evolution and equilibration of high-speed ageostrophic jets are studied under adiabatic condition relevant to hot-Jupiters. It is found that zonal jets can be baroclinically unstable, despite the planetary size of the Rossby deformation scale, and that high resolution is necessary to capture the process. Non-linear jet evolution is then used as a test case to assess model convergence in ve GCMs used in current hot-Jupiter simulations. The GCMs are also tested under a diabatic condition (thermal relaxation on a short timescale) similar to that used in many hot-Jupiter studies. In the latter case, in particular, the models show signi cant inter- and intra-model variability, limiting their quantitative prediction capability. Some models severely violate global angular momentum conservation. The generation of equatorial superrotation in Earth-like atmospheres, subject to \Held & Suarez-like" zonally-symmetric thermal forcing is also studied. It is shown that transition to superrotation occurs when the meridional gradient of the equilibrium surface entropy is weak in this setup. Two factors contribute to the onset of superrotation | suppression of breaking Rossby waves (generated by midlatitude baroclinic instability) that decelerate the equatorial ow, and, generation of inertial and barotropic instabilities in the equatorial region that provide the stirring to accelerate the equatorial ow. In summary, forcing condition and physical setup used in current hot-Jupiter simulations severely stretch model performance and predictive capability. Superrotation in Earth-like conditions may be common, but its strength decreases with resolution. Broadly, numerical convergence must be assessed in GCM experiments for each problem or setup considered.

In this work, we investigate the presence of stellar companions around hot Jupiter systems using data sets from the Clio and VISAO instruments on the Magellan Telescope. We observed eighteen targets of which eleven have known spin-orbit obliquity measurements. We detected eleven candidate companions of which five are new discoveries, five involved the validation and confirmation of previous studies, and one candidate proved to be a background star not bound to the transiting planet system. Out of eleven systems with known spin-orbit obliquity, seven systems have candidate companions. Due to the size of the sample, we could not find any correlation between the spin-orbit obliquity and the presence of a stellar companion. As future work, we will do follow up observations on the targets with candidate companions. We will increase our sample to one hundred systems to investigate if there is a correlation between spin-orbit obliquity and the presence of a distant stellar companion.

To date more than 1500 exoplanets have been discovered. A large number of these are hot Jupiters, Jupiter-sized planets orbiting < 0.1 au from their parent stars, due to limitations in observational techniques making them easier to detect than smaller planets in wider orbits. This is also, for the same reasons, the class of exoplanets with the most observational constraints. Due to the very large interaction between these planets and their parent stars they are believed to be tidally locked, causing a large temperature contrast between the permanently hot day side and colder night side. There are still many open questions about these planets. Many are observed to have inflated radii, i.e. the observed radius is larger for a given mass than evolutionary models predict. A mechanism that can transport some of the stellar heating into the interior of the planet may be able to explain this. The presence of hazes or clouds has been inferred on some planets, but their composition and distribution remain unknown. According to chemical equilibrium models TiO and VO should be present on the day side of the hottest of these planets, but these molecules have not yet been detected. Cold traps, where these molecules condense out on the night side, have been suggested to explain this. The efficiency of the heat redistribution from the day side to the night side has been found to vary significantly between different planets; the mechanism behind this is still unknown. To begin to answer many of these questions we need models capturing the three-dimensional nature of the atmospheres of these planets. General circulation models (GCMs) do this by solving the equations of fluid dynamics for the atmosphere coupled to a radiative transfer scheme. GCMs have previously been applied to several exoplanets, but many solve simplified fluid equations (shallow water or primitive equations) or highly parametrised radiation schemes (temperature-forcing, gray or band-averaged opacities). We here present an adaptation of the Met Office Unified Model (UM), a GCM used for weather predictions and climate studies for the Earth, to hot Jupiters. The UM solves the full 3D Euler equations for the fluid, and the radiation scheme uses the two-stream approximation and correlated-k method, which are state of the art for both Earth and exoplanet GCMs. This makes it ideally suited for the study of hot Jupiters. An important part of this work is devoted to the adaptation of the radiation scheme of the UM to hot Jupiters. This includes calculation of opacities for the main absorbers in these atmospheres from state-of-the-art high temperature line lists, the calculation of k-coefficients from these opacities, and making sure all aspects of the scheme perform satisfactorily at high temperatures and pressures. We have tested approximations made in previous works such as the two-stream approximation, use of band-averaged opacities and different treatments of gaseous overlap. Uncertainties in current models, such as the lack of high temperature line broadening parameters for these atmospheres, are discussed. We couple the adapted radiation scheme to the UM dynamical core, which has been tested independently. Our first application is devoted to one of the most well-observed hot Jupiters, HD 209458b. Differences between previous modelling works and our model are discussed, and we compare results from the full coupled model with results obtained using a temperature-forcing scheme. We have also developed a tool to calculate synthetic phase curves, and emission and transmission spectra from the output of our 3D model. This enables us to directly compare our model results to observations and test the effect of various parameters and model choices on observable quantities.

This dissertation explores the atmospheric circulation of extrasolar planets ranging from hot Jupiters to super Earths. For each of these studies, I utilize a three-dimensional circulation model coupled to a state-of-the-art, plane-parallel, two-stream, non-grey radiative transfer model dubbed the SPARC/MITgcm. First, I present models of the atmospheric circulation of eccentric hot Jupiters, a population which undergoes large variations in flux throughout their orbits. I demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits. For a select number of model integrations, I generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what is seen in infrared (IR) lightcurves. Next, I present circulation models of WASP-43b, a transiting hot Jupiter that is joining the ranks of HD 189733b and HD 209458b as a 'benchmark' hot Jupiter, with a wide array of observational constraints from the ground and space. Here I utilize the robust dataset of spectrophotometric observations taken with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST) to interpret my model results. I find that an atmospheric composition of 5x solar provides the best match to the data, particularly in emission. Lastly, I present atmospheric simulations of the super Earth GJ 1214b, exploring the planet's circulation as a function of atmospheric metallicity and composition. I find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity leads to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, I show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. In sum, these three studies explore exoplanet atmospheric circulation as a function of mass, radius, gravity, rotation rate, eccentricity and orbital distance.

Recent spectroscopic observations of transiting exoplanets have permitted the derivation of the thermal structure and molecular abundances of H₂O, CO, CO₂, CH₄, metallic oxides and alkali metals in these extreme atmospheres. Here, for the first time, a fully-fledged retrieval algorithm has been applied to exoplanet spectra to determine the thermal structure and composition. The development of a suite of radiative transfer and retrieval tools for exoplanet atmospheres is described, building upon an optimal estimation retrieval algorithm extensively used in solar system studies. Firstly, the collection of molecular line lists and the pre-tabulation of the absorption coefficients (k-distribution tables) for high temperature application are discussed. Secondly, the best-fit spectra for hot Jupiters are demonstrated and discussed case by case. Available sets of primary and secondary transit observations of exoplanets are used to retrieve atmospheric properties from these spectra, quantifying the limits of our knowledge of exoplanetary atmospheres based on the current quality of the data. The contribution functions and the vertical sensitivity to the molecules are fully utilised to interpret these spectra, probing the structure and composition of the atmosphere. Finally, the retrievals provide our best estimates of the thermal and compositional structure to date, using the covariance matrices to properly assess the degeneracy between different parameters and the uncertainties on derived quantities for the first time. This sheds light on the range of diverse interpretations offered by other authors so far, and allows us to scrutinise further atmospheric features by maximising the capability of the current retrieval algorithm and to demonstrate the need for broadband spectroscopy from future missions.

This report is exploring the possibility of characterizing hot Jupiter atmospheres using ground-based high-resolution spectroscopy. The ESO CRIRES infrared spectrometer is selected as the observing tool. Simulated observations are computed for known transiting systems. The properties of observations (noise, spectral coverage, resolution) are estimated with the CRIRES Exposure Time Calculator. An inverse method is used for reconstructing the transmission spectra of exoplanetary atmospheres and identifying spectral features. The possibility of using this method for non-transiting systems is examined. Three exoplanets are deemed possible to reconstruct the spectrum of. / Den här rapporten studerar möjligheten att karaktärisera heta Jupiter atmosfärer genom att använda markbaserad högupplöst spektroskopi. ESO:s CRIRES infraröda spektrometer används som observeringsvektyg. Simulerade observationer beräknas för kända system där planeten genomgår en transit. Flera aspekter kring observationer (brus, våglängdstäckning, upplösning) uppskattas med CRIRES Exposure Time Calculator. En invers metod används för att rekonstruera exoplanetatmosfärers transmissionsspektrum och för att identifiera spektraldrag. Möjligheten att använda den här metoden för system utan transit utforskas. Tre exoplaneter antas vara möjliga att rekonstruera spektrat från.

No description available.

Many known extra-solar giant planets lie close to their host stars. Around 60 have their semi-major axes smaller than 0.05 AU. In contrast to planets further out, the vast majority of these close-in planets have low eccentricity orbits. This suggests that their orbits have been circularized likely due to tidal dissipation inside the planets. These exoplanets share with our own Jupiter at least one trait in common: when they are subject to periodic tidal forcing, they behave like a lossy spring, with a tidal ``quality factor'', Q, of order 10^5. This parameter is the ratio between the energy in the tide and the energy dissipated per period. To explain this, a possible solution is resonantly forced internal oscillation. If the frequency of the tidal forcing happens to land on that of an internal eigenmode, this mode can be resonantly excited to a very large amplitude. The damping of such a mode inside the planet may explain the observed Q value. The only normal modes that fall in the frequency range of the tidal forcing (~ few days) are inertial modes, modes restored by the Coriolis force. We present a new numerical technique to solve for inertial modes in a convective, rotating sphere. This technique combines the use of an ellipsoidal coordinate system with a pseudo-spectral method to solve the partial differential equation that governs the inertial oscillations. We show that, this technique produces highly accurate solutions when the density profile is smooth. In particular, the lines of nodes are roughly parallel to the ellipsoidal coordinate axes. In particular, using these accurate solutions, we estimate the resultant tidal dissipation for giant planets, and find that turbulent dissipation of inertial modes in planets with smooth density profiles do not give rise to dissipation as strong as the one observed. We also study inertial modes in density profiles that exhibit discontinuities, as some recent models of Jupiter show. We found that, in this case, our method could not produce convergent solutions for the inertial modes. Additionally, we propose a way to observe inertial modes inside Saturn indirectly, by observing waves in its rings that may be excited by inertial modes inside Saturn.

Many known extra-solar giant planets lie close to their host stars. Around 60 have their semi-major axes smaller than 0.05 AU. In contrast to planets further out, the vast majority of these close-in planets have low eccentricity orbits. This suggests that their orbits have been circularized likely due to tidal dissipation inside the planets. These exoplanets share with our own Jupiter at least one trait in common: when they are subject to periodic tidal forcing, they behave like a lossy spring, with a tidal ``quality factor'', Q, of order 10^5. This parameter is the ratio between the energy in the tide and the energy dissipated per period. To explain this, a possible solution is resonantly forced internal oscillation. If the frequency of the tidal forcing happens to land on that of an internal eigenmode, this mode can be resonantly excited to a very large amplitude. The damping of such a mode inside the planet may explain the observed Q value. The only normal modes that fall in the frequency range of the tidal forcing (~ few days) are inertial modes, modes restored by the Coriolis force. We present a new numerical technique to solve for inertial modes in a convective, rotating sphere. This technique combines the use of an ellipsoidal coordinate system with a pseudo-spectral method to solve the partial differential equation that governs the inertial oscillations. We show that, this technique produces highly accurate solutions when the density profile is smooth. In particular, the lines of nodes are roughly parallel to the ellipsoidal coordinate axes. In particular, using these accurate solutions, we estimate the resultant tidal dissipation for giant planets, and find that turbulent dissipation of inertial modes in planets with smooth density profiles do not give rise to dissipation as strong as the one observed. We also study inertial modes in density profiles that exhibit discontinuities, as some recent models of Jupiter show. We found that, in this case, our method could not produce convergent solutions for the inertial modes. Additionally, we propose a way to observe inertial modes inside Saturn indirectly, by observing waves in its rings that may be excited by inertial modes inside Saturn.

Les hétérogénéités de surface des étoiles actives, telles que les taches et les facules, peuvent compliquer l'interprétation des spectres de transmission en introduisant des caractéristiques spectrales qui chevauchent celles d'atmosphère d'exoplanètes. Les courbes de lumière de transit d'HAT-P-18\(\,\)b et de WASP-52\(\,\)b, observées avec le mode SOSS de l'instrument NIRISS à bord du JWST, sont déformées par des occultations de taches. Avant le déploiement du JWST, ces régions actives étaient souvent simplement masquées, toutefois ceci peut mener à des mesures incorrectes des paramètres du transit. J'ai adapté et implémenté \(\texttt{spotrod}\), un modèle de transit avec occultation de taches, dans l'outil \(\texttt{Juliet}\) pour inférer conjointement les paramètres du transit et des taches occultées. J'ai ainsi ajusté les courbes de lumière de transit de ces deux Jupiters chaudes et récupéré la position de chaque tache, leur rayon et leur spectre de contraste, c'est-à-dire le rapport du flux de la tache sur le flux stellaire. J'ai contraint la température des taches et leur gravité de surface (pour prendre en compte les effets du champ magnétique local) en ajustant chaque spectre de contraste avec des spectres de modèles stellaires PHOENIX. Cependant, un certain degré de dégénérescence est présent, conduisant à une solution plus probable pour chaque tache, mais aussi à d'autres solutions qui ne peuvent être exclues. Le spectre de transmission d'HAT-P-18\(\,\)b nous a permis de détecter de l'H\(_2\)O (12,5\(\,\sigma\)) avec une abondance sub-solaire de \(\log\) H\(_2\)O \(\approx\) -4,4 \(\pm\) 0,3, des nuages (7,4\(\,\sigma\)) et du CO\(_2\) (7,3\(\,\sigma\)) dans l'atmosphère planétaire ainsi que des régions actives non occultées (5,8\(\,\sigma\)) qui imitent une pente de diffusion Rayleigh. / Surface heterogeneities on active stars, such as starspots and faculae, can complicate the interpretation of transmission spectra and introduce spectral features that overlap those of exoplanetary atmospheres. The transit light curves of HAT-P-18\(\,\)b and WASP-52\(\,\)b, observed in the SOSS mode of the NIRISS instrument aboard the JWST, are deformed by spot-crossings. These active regions were often simply masked before the launch of the JWST; however, this can prevent the correct measure of transit parameters. I adapted and implemented \(\texttt{spotrod}\), a model for transits of spotted stars, into the \(\texttt{Juliet}\) tool to simultaneously infer the transit and occulted starspots parameters. I fitted the transit light curves of these two hot Jupiters and retrieved for each spot its position, radius and spot-to-stellar flux contrast spectrum. I constrained the spots' temperature and surface gravity \(-\) attempting to capture the effects of the local magnetic pressure \(-\) by fitting each contrast spectrum with PHOENIX stellar model spectra. However, some degree of degeneracy is present, leading to a most likely solution for each starspot and other solutions that cannot be excluded. The transmission spectrum of HAT-P-18\(\,\)b enabled us to detect H\(_2\)O (12.5\(\,\sigma\)) with a sub-solar abundance of \(\log\) H\(_2\)O \(\approx\) -4.4 \(\pm\) 0.3, a cloud deck (7.4\(\,\sigma\)) and CO\(_2\) (7.3\(\,\sigma\)) in the planetary atmosphere as well as unocculted active regions (5.8\(\,\sigma\)) which mimic a Rayleigh scattering slope.