Hurricanes are one of the most extreme storm systems that occur on Earth, characterized by strong rainfall and fast winds. The terrestrial exoplanets that will be characterized with future infrared space telescopes orbit M dwarf stars. As a result, the best observable terrestrial exoplanets have vastly different climates than Earth, with a large dayside-to-nightside irradiation contrast and relatively slow rotation. Hurricanes may affect future observations of terrestrial exoplanets because they enhance the vertical transport of water vapor and could influence ocean heat transport. In this work, we explore how the environment of terrestrial exoplanets orbiting M dwarf stars affects the favorability of hurricane genesis (formation). To do so, we apply metrics developed to understand hurricane genesis on Earth to three-dimensional climate models of ocean-covered exoplanets orbiting M dwarf stars. We find that hurricane genesis is most favorable on intermediate-rotating tidally locked terrestrial exoplanets with rotation periods of $\sim 8-10~\mathrm{days}$. As a result, hurricane genesis is most favorable for terrestrial exoplanets in the habitable zones of late-type M dwarf stars. The peak in the favorability of hurricane genesis at intermediate rotation occurs because sufficient spin is required for hurricane genesis, but the vertical wind shear on fast-rotating terrestrial exoplanets disrupts hurricane genesis. We find that hurricane genesis is less favorable on slowly rotating terrestrial exoplanets, which agrees with previous work. Future work using simulations that resolve hurricane genesis and evolution can test our expectations for how the environment affects the favorability of hurricane genesis on tidally locked terrestrial exoplanets.

The potential habitability of tidally locked planets orbiting M-dwarf stars has been widely investigated in recent work, typically with a non-dynamic ocean and without continents. On Earth, ocean dynamics are a primary means of heat and nutrient distribution. Continents are a critical source of nutrients, strongly influence ocean dynamics, and participate in climate regulation. In this work, we investigate how the size of a substellar land mass affects the oceans ability to transport heat and upwell nutrients on the tidally locked planet Proxima Centauri b using the ROCKE-3D coupled ocean-atmosphere General Circulation Model (GCM). We find that dayside ice-free ocean and nutrient delivery to the mixed layer via upwelling are maintained across all continent sizes. We also find that Proxima Centauri bs climate is more sensitive to differences among atmospheric GCMs than to the inclusion of ocean dynamics in ROCKE-3D. Finally, we find that Proxima Centauri b transitions from a lobster state where ocean heat transport distributes heat away from the substellar point to an eyeball state where heat transport is restricted and surface temperature decreases symmetrically from the substellar point when the continent size exceeds about 20 percent of the surface area. Our work suggests that both a dynamic ocean and continents are unlikely to decrease the habitability prospects of nearby tidally locked targets like Proxima Centauri b that could be investigated with future observations by the James Webb Space Telescope (JWST).

Understanding the anomalous radii of many transiting hot gas giant planets is a fundamental problem of planetary science. Recent detections of re-inflated warm Jupiters orbiting post-main-sequence stars and the re-inflation of hot Jupiters while their host stars evolve on the main-sequence may help constrain models for the anomalous radii of hot Jupiters. In this work, we present evolution models studying the re-inflation of gas giants to determine how varying the depth and intensity of deposited heating affects both main-sequence re-inflation of hot Jupiters and post-main-sequence re-inflation of warm Jupiters. We find that deeper heating is required to re-inflate hot Jupiters than is needed to suppress their cooling, and that the timescale of re-inflation decreases with increasing heating rate and depth. We find a strong degeneracy between heating rate and depth, with either strong shallow heating or weak deep heating providing an explanation for main-sequence re-inflation of hot Jupiters. This degeneracy between heating rate and depth can be broken in the case of post-main-sequence re-inflation of warm Jupiters, as the inflation must be rapid to occur within post-main-sequence evolution timescales. We also show that the dependence of heating rate on incident stellar flux inferred from the sample of hot Jupiters can explain re-inflation of both warm and hot Jupiters. TESS will obtain a large sample of warm Jupiters orbiting post-main-sequence stars, which will help to constrain the mechanism(s) causing the anomalous radii of gas giant planets.

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

Hot Jupiters receive intense incident stellar light on their daysides, which drives vigorous atmospheric circulation that attempts to erase their large dayside-to-nightside flux contrasts. Propagating waves and instabilities in hot Jupiter atmospheres can cause emergent properties of the atmosphere to be time-variable. In this work, we study such weather in hot Jupiter atmospheres using idealized cloud-free general circulation models with double-grey radiative transfer. We find that hot Jupiter atmospheres can be time-variable at the $\sim 0.1-1\%$ level in globally averaged temperature and at the $\sim 1-10\%$ level in globally averaged wind speeds. As a result, we find that observable quantities are also time variable: the secondary eclipse depth can be variable at the $\lesssim 2\%$ level, the phase curve amplitude can change by $\lesssim 1\%$, the phase curve offset can shift by $\lesssim 5^{\circ}$, and terminator-averaged wind speeds can vary by $\lesssim 2~ \mathrm{km}~\mathrm{s}^{-1}$. Additionally, we calculate how the eastern and western limb-averaged wind speeds vary with incident stellar flux and the strength of an imposed drag that parameterizes Lorentz forces in partially ionized atmospheres. We find that the eastern limb is blueshifted in models over a wide range of equilibrium temperature and drag strength, while the western limb is only redshifted if equilibrium temperatures are $\lesssim1500~\mathrm{K}$ and drag is weak. Lastly, we show that temporal variability may be observationally detectable in the infrared through secondary eclipse observations with JWST, phase curve observations with future space telescopes (e.g., ARIEL), and/or Doppler wind speed measurements with high-resolution spectrographs.

We are on the verge of characterizing the atmospheres of terrestrial exoplanets in the habitable zones of M dwarf stars. Due to their large planet-to-star radius ratios and higher frequency of transits, terrestrial exoplanets orbiting M dwarf stars are favorable for transmission spectroscopy. In this work, we quantify the effect that water clouds have on the amplitude of water vapor transmission spectral features of terrestrial exoplanets orbiting M dwarf stars. To do so, we make synthetic transmission spectra from general circulation model (GCM) experiments of tidally locked planets. We improve upon previous work by considering how varying a broad range of planetary parameters affects transmission spectra. We find that clouds lead to a 10-100 times increase in the number of transits required to detect water features with the James Webb Space Telescope (JWST) with varying rotation period, incident stellar flux, surface pressure, planetary radius, and surface gravity. We also find that there is a strong increase in the dayside cloud coverage in our GCM simulations with rotation periods $\gtrsim 12 \ \mathrm{days}$ for planets with Earth's radius. This increase in cloud coverage leads to even stronger muting of spectral features for slowly rotating exoplanets orbiting M dwarf stars. We predict that it will be extremely challenging to detect water transmission features in the atmospheres of terrestrial exoplanets in the habitable zone of M dwarf stars with JWST. However, species that are well-mixed above the cloud deck (e.g., CO$_2$ and CH$_4$) may still be detectable on these planets with JWST.

Recent observations of ultra-hot Jupiters with dayside temperatures in excess of $2500~\mathrm{K}$ have found evidence for new physical processes at play in their atmospheres. In this work, we investigate the effects of the dissociation of molecular hydrogen and recombination of atomic hydrogen on the atmospheric circulation of ultra-hot Jupiters. To do so, we incorporate these effects into a general circulation model (GCM) for hot Jupiter atmospheres, and run a large suite of models varying the incident stellar flux, rotation period, and strength of frictional drag. We find that including hydrogen dissociation and recombination reduces the fractional day-to-night temperature contrast of ultra-hot Jupiter atmospheres and causes the speed of the equatorial jet to decrease in simulations with fixed rotation. This is because the large energy input required for hydrogen dissociation cools the dayside of the planet, and the energy released due to hydrogen recombination warms the nightside. The resulting decrease in the day-to-night temperature contrast reduces the day-to-night pressure gradient that drives the circulation, resulting in weaker wind speeds. The results from our GCM experiments qualitatively agree with previous theory which found that the fractional day-night temperature contrast of ultra-hot Jupiters should decrease with increasing equilibrium temperature due to hydrogen dissociation and recombination. Lastly, we compute full-phase light curves from our suite of GCM experiments, finding that the reduced day-to-night temperature contrast in ultra-hot Jupiter atmospheres causes a smaller phase curve amplitude. The reduction in phase curve amplitude due to hydrogen dissociation and recombination could explain the relatively small phase curve amplitudes of observed ultra-hot Jupiters.

The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari (2013), and determines how unstable the atmosphere is to baroclinic instability (the baroclinic "criticality") through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes.

Aerosols appear to be ubiquitous in close-in gas giant atmospheres, and disequilibrium chemistry likely impacts the emergent spectra of these planets. Lofted aerosols and disequilibrium chemistry are caused by vigorous vertical transport in these heavily irradiated atmospheres. Here we numerically and analytically investigate how vertical transport should change over the parameter space of spin-synchronized gas giants. In order to understand how tracer transport depends on planetary parameters, we develop an analytic theory to predict vertical velocities and mixing rates ($K_\mathrm{zz}$) and compare the results to our numerical experiments. We find that both our theory and numerical simulations predict that, if the vertical mixing rate is described by an eddy diffusivity, then this eddy diffusivity $K_\mathrm{zz}$ should increase with increasing equilibrium temperature, decreasing frictional drag strength, and increasing chemical loss timescales. We find that the transition in our numerical simulations between circulation dominated by a superrotating jet and that with solely day-to-night flow causes a marked change in the vertical velocity structure and tracer distribution. The mixing ratio of passive tracers is greatest for intermediate drag strengths that corresponds to this transition between a superrotating jet with columnar vertical velocity structure and day-to-night flow with upwelling on the dayside and downwelling on the nightside. Lastly, we present analytic solutions for $K_\mathrm{zz}$ as a function of planetary effective temperature, chemical loss timescales, and other parameters, for use as input to one-dimensional chemistry models of spin-synchronized gas giant atmospheres.

Launching in 2028, ESA's Atmospheric Remote-sensing Exoplanet Large-survey (ARIEL) survey of $\sim$1000 transiting exoplanets will build on the legacies of Kepler and TESS and complement JWST by placing its high precision exoplanet observations into a large, statistically-significant planetary population context. With continuous 0.5--7.8~$\mu$m coverage from both FGS (0.50--0.55, 0.8--1.0, and 1.0--1.2~$\mu$m photometry; 1.25--1.95~$\mu$m spectroscopy) and AIRS (1.95--7.80~$\mu$m spectroscopy), ARIEL will determine atmospheric compositions and probe planetary formation histories during its 3.5-year mission. NASA's proposed Contribution to ARIEL Spectroscopy of Exoplanets (CASE) would be a subsystem of ARIEL's FGS instrument consisting of two visible-to-infrared detectors, associated readout electronics, and thermal control hardware. FGS, to be built by the Polish Academy of Sciences' Space Research Centre, will provide both fine guiding and visible to near-infrared photometry and spectroscopy, providing powerful diagnostics of atmospheric aerosol contribution and planetary albedo, which play a crucial role in establishing planetary energy balance. The CASE team presents here an independent study of the capabilities of ARIEL to measure exoplanetary metallicities, which probe the conditions of planet formation, and FGS to measure scattering spectral slopes, which indicate if an exoplanet has atmospheric aerosols (clouds and hazes), and geometric albedos, which help establish planetary climate. Our design reference mission simulations show that ARIEL could measure the mass-metallicity relationship of its 1000-planet single-visit sample to $>7.5\sigma$ and that FGS could distinguish between clear, cloudy, and hazy skies and constrain an exoplanet's atmospheric aerosol composition to $>5\sigma$ for hundreds of targets, providing statistically-transformative science for exoplanet atmospheres.

A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the habitable zone for more massive planets occurs at higher stellar irradiation, as found in previous one-dimensional simulations. We also determine the relative effects of varying planetary radius and surface gravity. Increasing the planetary radius leads to a lower planetary albedo and warmer climate, pushing the inner edge of the habitable zone to lower stellar irradiation. This results from a change in circulation regime that leads to the disruption of the thick, reflective cloud deck around the substellar point. Increasing gravity increases the outgoing longwave radiation, which moves the inner edge of the habitable zone to higher stellar irradiation. This is because the column mass of water vapor decreases with increasing gravity, leading to a reduction in the greenhouse effect. The effect of gravity on the outgoing longwave radiation is stronger than the effect of radius on the planetary albedo, so that increasing gravity and radius together causes the inner edge of the habitable zone to move to higher stellar irradiation. Our results show that the inner edge of the habitable zone for more massive terrestrial planets occurs at a larger stellar irradiation.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

The recent detections of temperate terrestrial planets orbiting nearby stars and the promise of characterizing their atmospheres motivates a need to understand how the diversity of possible planetary parameters affects the climate of terrestrial planets. In this work, we investigate the atmospheric circulation and climate of terrestrial exoplanets orbiting both Sun-like and M-dwarf stars over a wide swath of possible planetary parameters, including the planetary rotation period, surface pressure, incident stellar flux, surface gravity, planetary radius, and cloud particle size. We do so using a general circulation model (GCM) that includes non-grey radiative transfer and the effects of clouds. The results from this suite of simulations generally show qualitatively similar dependencies of circulation and climate on planetary parameters as idealized GCMs, with quantitative differences due to the inclusion of additional model physics. Notably, we find that the effective cloud particle size is a key unknown parameter that can greatly affect the climate of terrestrial exoplanets. We confirm a transition between low and high dayside cloud coverage of synchronously rotating terrestrial planets with increasing rotation period. We determine that this cloud transition is due to eddy-driven convergence near the substellar point and should not be parameterization-dependent. Finally, we compute full-phase light curves from our simulations of planets orbiting M-dwarf stars, finding that changing incident stellar flux and rotation period affect observable properties of terrestrial exoplanets. Our GCM results can guide expectations for planetary climate over the broad range of possible terrestrial exoplanets that will be observed with future space telescopes.

Secondary eclipse observations of ultra-hot Jupiters have found evidence that hydrogen is dissociated on their daysides. Additionally, full-phase light curve observations of ultra-hot Jupiters show a smaller day-night emitted flux contrast than that expected from previous theory. Recently, it was proposed by Bell & Cowan (2018) that the heat intake to dissociate hydrogen and heat release due to recombination of dissociated hydrogen can affect the atmospheric circulation of ultra-hot Jupiters. In this work, we add cooling/heating due to dissociation/recombination into the analytic theory of Komacek & Showman (2016) and Zhang & Showman (2017) for the dayside-nightside temperature contrasts of hot Jupiters. We find that at high values of incident stellar flux, the day-night temperature contrast of ultra-hot Jupiters may decrease with increasing incident stellar flux due to dissociation/recombination, the opposite of that expected without including the effects of dissociation/recombination. We propose that a combination of a greater number of full-phase light curve observations of ultra-hot Jupiters and future General Circulation Models that include the effects of dissociation/recombination could determine in detail how the atmospheric circulation of ultra-hot Jupiters differs from that of cooler planets.

Because of their intense incident stellar irradiation and likely tidally locked spin states, hot Jupiters are expected to have wind speeds that approach or exceed the speed of sound. In this work we develop a theory to explain the magnitude of these winds. We model hot Jupiters as planetary heat engines and show that hot Jupiters are always less efficient than an ideal Carnot engine. Next, we demonstrate that our predicted wind speeds match those from three-dimensional numerical simulations over a broad range of parameters. Finally, we use our theory to evaluate how well different drag mechanisms can match the wind speeds observed with Doppler spectroscopy for HD 189733b and HD 209458b. We find that magnetic drag is potentially too weak to match the observations for HD 189733b, but is compatible with the observations for HD 209458b. In contrast, shear instabilities and/or shocks are compatible with both observations. Furthermore, the two mechanisms predict different wind speed trends for hotter and colder planets than currently observed. As a result, we propose that a wider range of Doppler observations could reveal multiple drag mechanisms at play across different hot Jupiters.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

Hot Jupiters receive strong stellar irradiation, producing equilibrium temperatures of $1000 - 2500 \ \mathrm{Kelvin}$. Incoming irradiation directly heats just their thin outer layer, down to pressures of $\sim 0.1 \ \mathrm{bars}$. In standard irradiated evolution models of hot Jupiters, predicted transit radii are too small. Previous studies have shown that deeper heating -- at a small fraction of the heating rate from irradiation -- can explain observed radii. Here we present a suite of evolution models for HD 209458b where we systematically vary both the depth and intensity of internal heating, without specifying the uncertain heating mechanism(s). Our models start with a hot, high entropy planet whose radius decreases as the convective interior cools. The applied heating suppresses this cooling. We find that very shallow heating -- at pressures of $1 - 10 \ \mathrm{bars}$ -- does not significantly suppress cooling, unless the total heating rate is $\gtrsim 10\%$ of the incident stellar power. Deeper heating, at $100 \ \mathrm{bars}$, requires heating at only $1\%$ of the stellar irradiation to explain the observed transit radius of $1.4 R_{\rm Jup}$ after 5 Gyr of cooling. In general, more intense and deeper heating results in larger hot Jupiter radii. Surprisingly, we find that heat deposited at $10^4 \ \mathrm{bars}$ -- which is exterior to $\approx 99\%$ of the planet's mass -- suppresses planetary cooling as effectively as heating at the center. In summary, we find that relatively shallow heating is required to explain the radii of most hot Jupiters, provided that this heat is applied early and persists throughout their evolution.

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing fractional dayside-nightside brightness temperature difference with increasing incident stellar flux, both averaged across the infrared and in each individual wavelength band. The analytic theory of Komacek & Showman (2016) shows that this trend is due to the decreasing ability with increasing incident stellar flux of waves to propagate from day to night and erase temperature differences. Here, we compare the predictions of this theory to observations, showing that it explains well the shape of the trend of increasing dayside-nightside temperature difference with increasing equilibrium temperature. Applied to individual planets, the theory matches well with observations at high equilibrium temperatures but, for a fixed photosphere pressure of $100 \ \mathrm{mbar}$, systematically under-predicts the dayside-nightside brightness temperature differences at equilibrium temperatures less than $2000 \ \mathrm{K}$. We interpret this as due to as the effects of a process that moves the infrared photospheres of these cooler hot Jupiters to lower pressures. We also utilize general circulation modeling with double-grey radiative transfer to explore how the circulation changes with equilibrium temperature and drag strengths. As expected from our theory, the dayside-nightside temperature differences from our numerical simulations increase with increasing incident stellar flux and drag strengths. We calculate model phase curves using our general circulation models, from which we compare the broadband infrared offset from the substellar point and dayside-nightside brightness temperature differences against observations, finding that strong drag or additional effects (e.g. clouds and/or supersolar metallicities) are necessary to explain many observed phase curves.

Terrestrial exoplanets in the canonical habitable zone may have a variety of initial water fractions due to random volatile delivery by planetesimals. If the total planetary water complement is high, the entire surface may be covered in water, forming a "waterworld." On a planet with active tectonics, competing mechanisms act to regulate the abundance of water on the surface by determining the partitioning of water between interior and surface. Here we explore how the incorporation of different mechanisms for the degassing and regassing of water changes the volatile evolution of a planet. For all of the models considered, volatile cycling reaches an approximate steady-state after $\sim 2 \ \mathrm{Gyr}$. Using these steady-states, we find that if volatile cycling is either solely dependent on temperature or seafloor pressure, exoplanets require a high abundance ($\gtrsim 0.3\%$ of total mass) of water to have fully inundated surfaces. However, if degassing is more dependent on seafloor pressure and regassing mainly dependent on mantle temperature, the degassing rate is relatively large at late times and a steady-state between degassing and regassing is reached with a substantial surface water fraction. If this hybrid model is physical, super-Earths with a total water fraction similar to that of the Earth can become waterworlds. As a result, further understanding of the processes that drive volatile cycling on terrestrial planets is needed to determine the water fraction at which they are likely to become waterworlds.

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing dayside-to-nightside brightness temperature difference with increasing equilibrium temperature. Here we present a three-dimensional model that explains this relationship, in order to shed insight on the processes that control heat redistribution in tidally-locked planetary atmospheres. This three-dimensional model combines predictive analytic theory for the atmospheric circulation and dayside-nightside temperature differences over a range of equilibrium temperature, atmospheric composition, and potential frictional drag strengths with numerical solutions of the circulation that verify this analytic theory. This analytic theory shows that the longitudinal propagation of waves mediates dayside-nightside temperature differences in hot Jupiter atmospheres, analogous to the wave adjustment mechanism that regulates the thermal structure in Earth's tropics. These waves can be damped in hot Jupiter atmospheres by either radiative cooling or potential frictional drag. This frictional drag would likely be caused by Lorentz forces in a partially ionized atmosphere threaded by a background magnetic field, and would increase in strength with increasing temperature. Additionally, the amplitude of radiative heating and cooling increases with increasing temperature, and hence both radiative heating/cooling and frictional drag damp waves more efficiently with increasing equilibrium temperature. Radiative heating and cooling play the largest role in controlling dayside-nightside temperature temperature differences in both our analytic theory and numerical simulations, with frictional drag only important if it is stronger than the Coriolis force. As a result, dayside-nightside temperature differences in hot Jupiter atmospheres increase with increasing stellar irradiation and decrease with increasing pressure.

We present magnetohydrodynamic (MHD) simulations of the atmospheres of hot Jupiters ranging in temperature from 1100-1800K. Magnetic effects are negligible in atmospheres with temperatures $\lesssim$ 1400K. At higher temperatures winds are variable and in many cases, mean equatorial flows can become westward, opposite to their hydrodynamic counterparts. Ohmic dissipation peaks at temperatures $\sim$1500-1600K, depending on field strength, with maximum values $\sim 10^{18}$W at 10bar, substantially lower than previous estimates. Based on the limited parameter study done, this value can not be increased substantially with increasing winds, higher temperatures, higher field strengths, different boundary conditions or lower diffusivities. Although not resolved in these simulations there is modest evidence that a magnetic buoyancy instability may proceed in hot atmospheres.