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.

Super Earths are a class of exoplanets with masses between 1-10 M⊕. Comprising nearly 70 % of the discovered planet population, they are largest class of exoplanets known. Super Earths exhibit an interesting variety of compositions, as their densities imply that they range from dense, rocky planets to those with substantial amounts of water. This thesis aims to understand why super Earths form so frequently, and to connect the final compositions of super Earths to the regions where they form in protoplanetary disks. To do this, we combine a model that calculates the physical and chemical conditions within a protoplanetary disk with a core accretion model of planet formation. A key feature of our planet formation model is planet traps that act as barriers to rapid type-I migration. The traps we include in our model are the dead zone, which can be caused by cosmic ray or X-ray ionization, the ice line, and the heat transition. In disks with lifetimes >􏰐 4 Myr we find that planet formation in all traps results in Jovian planets. Typically, the X-ray dead zone and heat transition traps produce hot Jupiters orbiting near 0.05 AU while the cosmic ray dead zone and ice line traps produce Jupiters near 1 AU. Super Earths are found to form in disks with short lifetimes 􏰑< 2 Myr that photoevaporate prior to planets undergoing runaway gas accretion. Additionally, we find that super Earth formation takes place in low-mass disks (<􏰑 0.05 M⊙), where planet formation timescales exceed disk lifetimes inferred through observations. The location of various traps throughout the disk play a key role in allowing super Earths to achieve a range of compositions. Super Earths forming in the ice line or heat transition accrete solids from cold regions of the disk, resulting in planets with large ice contents (up to 50 % by mass). Conversely, super Earths formed in the dead zone trap accrete solids from warm regions of the disk and are therefore composed of mostly rocky materials (less than 5 % ice by mass). / Thesis / Master of Science (MSc)

Atmospheres play an important role in the habitability of a planet, so understanding and modelling them is an important step in the search for life on other planets. This thesis presents a 1D frequency-dependent radiative-convective code that was written to help determine the temperature-pressure structure of potentially-habitable exoplanets. This code pairs with a chemistry model to determine the chemical composition of these planets' atmospheres. This code is applied to the planets in the TRAPPIST-1 system. The initial atmospheric compositions of the TRAPPIST-1 planets are determined through planet formation history and considered for both outgassed and accreted atmospheres. An interesting result is found when running these initial atmospheric compositions through the chemistry model: when the atmosphere equilibrates, it can change its C/O ratio from equal to that found in the accreted or outgassed volatiles to something lower, because, in temperate conditions, CO$_2$ is favoured over CO. This has the consequence that observed C/O ratios in terrestrial atmospheres cannot be relied on to infer the C/O ratio of the protoplanetary disc in which the planet formed. The initial results of atmospheric modelling for TRAPPIST-1 planets indicate that these planets are likely to have relatively warmer upper atmospheres due to the fact that their host star emits primarily in the infrared, and a portion of this radiation is then absorbed as it enters the top of the atmosphere. These initial results have not been seen in previous work. These initial results are the beginning of a database of potential atmospheres on the TRAPPIST-1 planets. It is hoped that these atmospheres can be compared with observations from future observing missions like the James Webb Space Telescope to help constrain the surface conditions of these potentially-habitable planets and ultimately, to help in the search for life. / Thesis / Master of Science (MSc)

Plane-layer geometry convection models remain useful for modelling planetary mantle dynamics however they yield significantly warmer mean temperatures than spherical shell models. For example, in a uniform property spherical shell with the same radius ratio, f, as the Earth's mantle; a bottom heating Rayleigh number, Ra, of 10^7 and a nondimensional internal heating rate, H, of 23 (arguably Earth-like values) are insufficient to heat the mean temperature, θ, above the mean of the non-dimensional boundary value temperatures (0.5), the temperature in a plane-layer model with no internal heating. This study investigates the impact of this geometrical effect in convection models featuring uniform and stratified viscosity. To address the effect of geometry, heat sinks are implemented to lower the mean temperature in 3D plane-layer isoviscous convection models. Over 100 models are analyzed, and their mean temperatures are used to derive a single equation for predicting θ, as a function of Ra, H and f in spherical and plane-layer systems featuring free-slip surfaces. The inclusion of first-order terrestrial characteristics is introduced to quantitatively assess the influence of system geometry on planetary scale simulations. Again, over 100 models are analyzed featuring a uniform upper mantle viscosity and a lower mantle viscosity that increases by a factor of 30 or 100. An effective Rayleigh number, Raη, is defined based on the average viscosity of the mantle. Equations for the relationship between θ, Raη, and H are derived for convection in a spherical shell with f = 0.547 and plane-layer geometries. These equations can be used to determine the appropriate heating rate for a plane-layer convection model to emulate spherical shell convection mean temperatures for effective Rayleigh numbers comparable to the Earth’s value and greater. Comparing cases with the same H and Raη, the increased lower mantle viscosity amplifies the mismatch in mean temperatures between spherical shell and plane-layer models. These findings emphasize the importance of adjusting heating rates in plane-layer geometry models and have important implications for studying convection with temperature-dependent parameters in plane-layer systems. The findings are particularly relevant to the study of convection in super-Earths where full spherical shell calculations remain intractable.

Plane-layer geometry convection models remain useful for modelling planetary mantle dynamics however they yield significantly warmer mean temperatures than spherical shell models. For example, in a uniform property spherical shell with the same radius ratio, f, as the Earth's mantle; a bottom heating Rayleigh number, Ra, of 10^7 and a nondimensional internal heating rate, H, of 23 (arguably Earth-like values) are insufficient to heat the mean temperature, θ, above the mean of the non-dimensional boundary value temperatures (0.5), the temperature in a plane-layer model with no internal heating. This study investigates the impact of this geometrical effect in convection models featuring uniform and stratified viscosity. To address the effect of geometry, heat sinks are implemented to lower the mean temperature in 3D plane-layer isoviscous convection models. Over 100 models are analyzed, and their mean temperatures are used to derive a single equation for predicting θ, as a function of Ra, H and f in spherical and plane-layer systems featuring free-slip surfaces. The inclusion of first-order terrestrial characteristics is introduced to quantitatively assess the influence of system geometry on planetary scale simulations. Again, over 100 models are analyzed featuring a uniform upper mantle viscosity and a lower mantle viscosity that increases by a factor of 30 or 100. An effective Rayleigh number, Raη, is defined based on the average viscosity of the mantle. Equations for the relationship between θ, Raη, and H are derived for convection in a spherical shell with f = 0.547 and plane-layer geometries. These equations can be used to determine the appropriate heating rate for a plane-layer convection model to emulate spherical shell convection mean temperatures for effective Rayleigh numbers comparable to the Earth’s value and greater. Comparing cases with the same H and Raη, the increased lower mantle viscosity amplifies the mismatch in mean temperatures between spherical shell and plane-layer models. These findings emphasize the importance of adjusting heating rates in plane-layer geometry models and have important implications for studying convection with temperature-dependent parameters in plane-layer systems. The findings are particularly relevant to the study of convection in super-Earths where full spherical shell calculations remain intractable.

The detection of low-mass extra-solar planets through radial-velocity searches is currently limited by the intrinsic magnetic activity of the host stars. The correlated noise that arises from their natural radial-velocity variability can easily mimic or conceal the orbital signals of super-Earth and Earth-mass extra-solar planets. I developed an intuitive and robust data analysis framework in which the activity-induced variations are modelled with a Gaussian process that has the frequency structure of the photometric variations of the star, thus allowing me to determine precise and reliable planetary masses. I applied this technique to three recently discovered planetary systems: CoRoT-7, Kepler-78 and Kepler-10. I determined the masses of the transiting super-Earth CoRoT-7b and the small Neptune CoRoT-7c to be 4.73 ± 0.95 M⊕ and 13.56 ± 1.08 M⊕, respectively. The density of CoRoT-7b is 6.61 ± 1.72 g.cm⁻³, which is compatible with a rocky composition. I carried out Bayesian model selection to assess the nature of a previously identified signal at 9 days, and found that it is best interpreted as stellar activity. Despite the high levels of activity of its host star, I determined the mass of the Earth-sized planet Kepler-78b to be 1.76 ± 0.18 M⊕. With a density of 6.2(+1.8:-1.4) g.cm⁻³, it is also a rocky planet. I found the masses of Kepler-10b and Kepler-10c to be 3.31 ± 0.32 M⊕ and 16.25 ± 3.66 M⊕, respectively. Their densities, of 6.4(+1.1:-0.7) g.cm⁻³ and 8.1 ± 1.8 g.cm⁻³, imply that they are both of rocky composition – even the 2 Earth-radius planet Kepler-10c! In parallel, I deepened our understanding of the physical origin of stellar radial-velocity variability through the study of the Sun, which is the only star whose surface can be imaged at high resolution. I found that the full-disc magnetic flux is an excellent proxy for activity-induced radial-velocity variations; this result may become key to breaking the activity barrier in coming years. I also found that in the case of CoRoT-7, the suppression of convective blueshift leads to radial-velocity variations with an rms of 1.82 m.s⁻¹, while the modulation induced by the presence of dark spots on the rotating stellar disc has an rms of 0.46 m.s⁻¹. For the Sun, I found these contributions to be 2.22 m.s⁻¹ and 0.14 m.s⁻¹, respectively. These results suggest that for slowly rotating stars, the suppression of convective blueshift is the dominant contributor to the activity-modulated radial-velocity signal, rather than the rotational Doppler shift of the flux blocked by starspots.