The Cassini mission entered the Grand Finale phase in April 2017 and executed 22.5 highly inclined, close-in orbits around Saturn before diving into the planet on September 15th 2017. Here we present our analysis of the Cassini Grand Finale magnetometer (MAG) dataset, focusing on Saturn's internal magnetic field. These measurements demonstrate that Saturn's internal magnetic field is exceptionally axisymmetric, with a dipole tilt less than 0.007 degrees (25.2 arcsecs). Saturn's magnetic equator was directly measured to be shifted northward by ~ 0.0468 +/- 0.00043 (1-sigma) $R_S$, 2820 +/- 26 km, at cylindrical radial distances between 1.034 and 1.069 $R_S$ from the spin-axis. Although almost perfectly axisymmetric, Saturn's internal magnetic field exhibits features on many characteristic length scales in the latitudinal direction. Examining Br at the a=0.75 $R_S$, c=0.6993 $R_S$ isobaric surface, the degrees 4 to 11 contributions correspond to latitudinally banded magnetic perturbations with characteristic width similar to that of the off-equatorial zonal jets observed in the atmosphere of Saturn. Saturn's internal magnetic field beyond 60 degrees latitude, in particular the small-scale features, are less well constrained by the available measurements, mainly due to incomplete spatial coverage in the polar region. A stably stratified layer thicker than 2500 km likely exists above Saturn's deep dynamo to filter out the non-axisymmetric internal magnetic field. A heat transport mechanism other than pure conduction, e.g. double diffusive convection, must be operating within this layer to be compatible with Saturn's observed luminosity. The latitudinally banded magnetic perturbations likely arise from a shallow secondary dynamo action with latitudinally banded differential rotation in the semi-conducting layer.

The nature and structure of the observed east-west flows on Jupiter and Saturn has been one of the longest-lasting mysteries in planetary science. This mystery has been recently unraveled due to the accurate gravity measurements provided by the Juno mission to Jupiter and the Grand Finale of the Cassini mission to Saturn. These two experiments, which coincidentally happened around the same time, allowed determination of the vertical and meridional profiles of the zonal flows on both planets. This paper reviews the topic of zonal jets on the gas giants in light of the new data from these two experiments. The gravity measurements not only allow the depth of the jets to be constrained, yielding the inference that the jets extend roughly 3000 and 9000 km below the observed clouds on Jupiter and Saturn, respectively, but also provide insights into the mechanisms controlling these zonal flows. Specifically, for both planets this depth corresponds to the depth where electrical conductivity is within an order of magnitude of 1 S/m, implying that the magnetic field likely plays a key role in damping the zonal flows.

A suite of discoveries in the last two decades demonstrate that we are now at a point where incorporating magnetic behavior is key for advancing our ability to characterize substellar and planetary systems. The next decade heralds the exciting maturation of the now-burgeoning field of brown dwarf magnetism, and investing now in brown dwarf magnetism will provide a key platform for exploring exoplanetary magnetism and habitability beyond the solar system. We anticipate significant discoveries including: the nature of substellar and planetary magnetic dynamos, the characterization of exo-aurora physics and brown dwarf magnetospheric environments, and the role of satellites in manifestations of substellar magnetic activity. These efforts will require significant new observational capabilities at radio and near infrared wavelengths, dedicated long-term monitoring programs, and committed support for the theoretical modeling efforts underpinning the physical processes of the magnetic phenomena

The higher order gravitational moments of a differentially rotating planet are greatly affected or even dominated by the near surface differential flow. Unlike the contribution from rigid body rotation, this part of the gravity field can be well estimated by linear theory when the differential rotation is on cylinders and the corresponding gravity field arises from the higher order moments directly introduced by that flow. In the context of an n = 1 polytrope, we derive approximate analytical formulas for the gravity moments. We find that $\Delta J_{2\mathrm{n}}$ is typically at most a few times $(-1)^{n+1} \; a \; q \; d^{5/2}$, where $a$ is the amplitude of the differential rotation (as a fraction of the background rigid body rotation), $q = \Omega^2 R^3 / G M$ is the usual dimensionless measure of rotation for the planet (mass $M$, radius $R$), and $d << 1$ is the characteristic depth of the flow as a fraction of the planetary radius. Applied to Saturn, with $a$ set by the observed surface wind amplitude, we find first that the observed signs of the $\Delta J_{2\mathrm{n}}$ are a trivial consequence of the definition of the corresponding Legendre polynomials, but the {\it Cassini} observations can not be explained by a simple exponentially decaying flow and instead require a substantial retrograde flow at depth and a larger $d$ than the simple scaling suggests. This is consistent with the results reported by \citet{Iess2019}. However, there is no fluid dynamical requirement that the flows observed in the atmosphere are a guide to the flows thousands of km deeper. We explore a wide range of flow depths and amplitudes which yield values for $\Delta J_{\mathrm{n}}$ that are acceptable within the error estimates and thus highlight the inherent non-uniqueness of inferences made from the higher order gravity moments.

The Earth's Moon is thought to have formed from a circumterrestrial disk generated by a giant impact between the proto-Earth and an impactor approximately 4.5 billion years ago. Since the impact was energetic, the disk would have been hot and partially vaporized. This formation process is thought to be responsible for the geochemical observation that the Moon is depleted in volatiles. This model predicts that the Moon should be significantly depleted in water as well, but this appears to contradict some of the recently measured lunar water abundances and D/H ratios that suggest that the Moon is more water-rich than previously thought. Alternatively, the Moon could have retained its water if the upper of the disk were dominated by heavier species because hydrogen would have had to diffuse out from the heavy-element rich disk, and therefore the escape rate would have been limited by this slow diffusion process (diffusion-limited escape). To identify which escape the disk would have experienced and to quantify volatiles loss from the disk, we compute the thermal structure of the Moon-forming disk considering various bulk water abundances and mid-plane disk temperatures. Our calculations show that the upper parts of the Moon-forming disk are dominated by heavy atoms or molecules and hydrogen is a minor species. This indicates that hydrogen escape would have been diffusion-limited, and therefore the amount of lost water and hydrogen would have been small compared to the initial abundance assumed. This result indicates that the giant impact hypothesis can be consistent with the water-rich Moon. Furthermore, since the hydrogen wind would have been weak, the other volatiles would not have escaped either. Thus, the observed volatile depletion of the Moon requires another mechanism.

Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M$_\oplus$, and radius, 3.68 R$_\oplus$, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.

We have used NSF's Karl G. Jansky Very Large Array (VLA) to observe a sample of five known radio-emitting late L and T dwarfs ranging in age from ~0.2 - 3.4Gyr. We observed each target for seven hours, extending to higher frequencies than previously attempted and establishing proportionally higher limits on maximum surface magnetic field strengths. Detections of circularly polarized pulses at 8 - 12~GHz yield measurements of 3.2 - 4.1 kG localized magnetic fields on four of our targets, including the archetypal cloud variable and likely planetary-mass object T2.5 dwarf SIMP J01365663+0933473. We additionally detect a pulse at 15 - 16.5 GHz for the T6.5 dwarf 2MASS 10475385+2124234, corresponding to a localized 5.6 kG field strength. For the same object, we tentatively detect a 16.5 - 18 GHz pulse, corresponding to a localized 6.2 kG field strength. We measure rotation periods between ~1.47 - 2.28 hr for 2MASS J10430758+2225236, 2MASS J12373919+6526148, and SDSS J04234858-0414035, supporting (i) an emerging consensus that rapid rotation may be important for producing strong dipole fields in convective dynamos and/or (ii) rapid rotation is a key ingredient for driving the current systems powering auroral radio emission. We observe evidence of variable structure in the frequency-dependent time series of our targets on timescales shorter than a rotation period, suggesting a higher degree of variability in the current systems near the surfaces of brown dwarfs. Finally, we find that age, mass, and temperature together cannot account for the strong magnetic fields produced by our targets.

Recent (2011) results from the Nice Observatory indicate the existence of global seismic modes on Jupiter in the frequency range between 0.7 and 1.5mHz with amplitudes of tens of cm/s. Currently, the driving force behind these modes is a mystery; the measured amplitudes are many orders of magnitude larger than anticipated based on theory analogous to heliosiesmology (that is, turbulent convection as a source of stochastic excitation). One of the most promising hypotheses is that these modes are driven by Jovian storms. This work constructs a framework to analytically model the expected equilibrium normal mode amplitudes arising from convective columns in storms. We also place rough constraints on Jupiter's seismic modal quality factor. Using this model, neither meteor strikes, turbulent convection, nor water storms can feasibly excite the order of magnitude of observed amplitudes. Next we speculate about the potential role of rock storms deeper in Jupiter's atmosphere, because the rock storms' expected energy scales make them promising candidates to be the chief source of excitation for Jovian seismic modes, based on simple scaling arguments. We also suggest some general trends in the expected partition of energy between different frequency modes. Finally we supply some commentary on potential applications to gravity, Juno, Cassini and Saturn, and future missions to Uranus and Neptune.

Low-mass low-density planets discovered by Kepler in the super-Earth mass regime typically have large radii for their inferred masses, implying the presence of H$_2$-He atmospheres. These planets are vulnerable to atmospheric mass loss due to heating by the parent star's XUV flux. Models coupling atmospheric mass loss with thermal evolution predicted a bimodal distribution of planetary radii, which has gained observational support. However, a key component that has been ignored in previous studies is the dissolution of these gases into the molten core of rock and iron that constitute most of their mass. Such planets have high temperatures ($>$2000 K) and pressures ($\sim$kbars) at the core-envelope boundary, ensuring a molten surface and a subsurface reservoir of hydrogen that can be 5-10 times larger than the atmosphere. This study bridges this gap by coupling the thermal evolution of the planet and the mass loss of the atmosphere with the thermodynamic equilibrium between the dissolved H$_2$ and the atmospheric H$_2$ (Henry's law). Dissolution in the interior allows a planet to build a larger hydrogen repository during the planet formation stage. We show that the dissolved hydrogen outgasses to buffer atmospheric mass loss. The slow cooling of the planet also leads to outgassing because solubility decreases with decreasing temperature. Dissolution of hydrogen in the interior therefore increases the atmosphere retention ability of super-Earths. The study highlights the importance of including the temperature- and pressure-dependent solubility of gases in magma oceans and coupling outgassing to planetary evolution models.

Giant planets are thought to have cores in their deep interiors, and the division into a heavy-element core and hydrogen-helium envelope is applied in both formation and structure models. We show that the primordial internal structure depends on the planetary growth rate, in particular, the ratio of heavy elements accretion to gas accretion. For a wide range of likely conditions, this ratio is in one-to-one correspondence with the resulting post-accretion profile of heavy elements within the planet. This flux ratio depends sensitively on the assumed solid surface density in the surrounding nebula. We suggest that giant planets' cores might not be distinct from the envelope and includes some hydrogen and helium, and the deep interior can have a gradual heavy-element structure. Accordingly, Jupiter's core may not be well-defined. Accurate measurements of Jupiter's gravitational field by Juno could put constraints on Jupiter's core mass. However, as we suggest here, the definition of Jupiter's core is complex, and the core's physical properties (mass, density) depend on the actual definition of the core and on its growth history.

Any nonspherical distribution of density inside planets and stars gives rise to a non-spherical external gravity and change of shape. If part or all of the observed zonal flows at the cloud deck of Jupiter and Saturn represent deep interior dynamics, then the density perturbations associated with the deep zonal flows could generate gravitational signals detectable by the Juno mission and the Cassini Grand Finale. Here we present a critical examination of the applicability of the thermal wind equation to calculate the wind induced gravity moments. Our analysis shows that wind induced gravity moments calculated from TWE are in overall agreement with the full solution to the Euler equation. However, the accuracy of individual high-degree moments calculated from TWE depends crucially on retaining the nonsphericity of the background density and gravity. Only when the background nonsphericity of the planet is taken into account, does the TWE make accurate enough prediction (with a few tens of percent errors) for individual high-degree gravity moments associated with deep zonal flows. Since the TWE is derived from the curl of the Euler equation and is a local relation, it necessarily says nothing about any density perturbations that contribute irrotational terms to the Euler equation and that have a non-local origin. However, the predicted corrections from these density contributions to the low harmonic degree gravity moments are not discernible from insignificant changes in interior models while the corrections at high harmonic degree are very small, a few percent or less.

All four giant planets in the Solar System feature zonal flows on the order of 100 m/s in the cloud deck, and large-scale intrinsic magnetic fields on the order of 1 Gauss near the surface. The vertical structure of the zonal flows remains obscure. The end-member scenarios are shallow flows confined in the radiative atmosphere and deep flows throughout the entire planet. The electrical conductivity increases rapidly yet smoothly as a function of depth inside Jupiter and Saturn. Deep zonal flows will inevitably interact with the magnetic field, at depth with even modest electrical conductivity. Here we investigate the interaction between zonal flows and magnetic fields in the semi-conducting region of giant planets. Employing mean-field electrodynamics, we show that the interaction will generate detectable poloidal magnetic field perturbations spatially correlated with the deep zonal flows. Assuming the peak amplitude of the dynamo alpha-effect to be 0.1 mm/s, deep zonal flows on the order of 0.1 - 1 m/s in the semi-conducting region of Jupiter and Saturn would generate poloidal magnetic perturbations on the order of 0.01% - 1% of the background dipole field. These poloidal perturbations should be detectable with the in-situ magnetic field measurements from the Juno mission and the Cassini Grand Finale. This implies that magnetic field measurements can be employed to constrain the properties of deep zonal flows in the semi-conducting region of giant planets.

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

We have observed 6 late-L and T dwarfs with the Karl G. Jansky Very Large Array (VLA) to investigate the presence of highly circularly polarized radio emission, associated with large-scale auroral currents. Previous surveys encompassing ~60 L6 or later targets in this spectral range have yielded only one detection. Our sample includes the previously detected T6.5 dwarf 2MASS 10475385+2124234 as well as 5 new targets selected for the presence of H-alpha emission or optical/infrared photometric variability, which are possible manifestations of auroral activity. We detect 2MASS 10475385+2124234, as well as 4 of the 5 targets in our biased sample, including the strong IR variable SIMP J01365662+0933473 and bright H-alpha emitter 2MASS 12373919+6526148, reinforcing the possibility that activity at these disparate wavelengths is related. The radio emission frequency corresponds to a precise determination of the lower-bound magnetic field strength near the surface of each dwarf and this new sample provides robust constraints on dynamo theory in the low mass brown dwarf regime. Magnetic fields >2.5 kG are confirmed for 5/6 targets. Our results provide tentative evidence that the dynamo operating in this mass regime may be inconsistent with predicted values from a recently proposed model. Further observations at higher radio frequencies are essential for verifying this assertion.

The Earth's Moon is thought to have formed by an impact between the Earth and an impactor around 4.5 billion years ago. This impact could have been so energetic that it could have mixed and homogenized the Earth's mantle. However, this view appears to be inconsistent with geochemical studies that suggest that the Earth's mantle was not mixed by the impact. Another plausible outcome is that this energetic impact melted the whole mantle, but the extent of mantle melting is not well understood even though it must have had a significant effect on the subsequent evolution of the Earth's interior and atmosphere. To understand the initial state of the Earth's mantle, we perform giant impact simulations using smoothed particle hydrodynamics (SPH) for three different models: (a) standard: a Mars-sized impactor hits the proto-Earth, (b) fast-spinning Earth: a small impactor hits a rapidly rotating proto-Earth, and (c) sub-Earths: two half Earth-sized planets collide. We use two types of equations of state (MgSiO3 liquid and forsterite) to describe the Earth's mantle. We find that the mantle remains unmixed in (a), but it may be mixed in (b) and (c). The extent of mixing is most extensive in (c). Therefore, (a) is most consistent and (c) may be least consistent with the preservation of the mantle heterogeneity, while (b) may fall between. We determine that the Earth's mantle becomes mostly molten by the impact in all of the models. The choice of the equations of state does not affect these outcomes. Additionally, our results indicate that entropy gains of the mantle materials by a giant impact cannot be predicted well by the Rankine-Hugoniot equations. Moreover, we show that the mantle can remain unmixed on a Moon-forming timescale if it does not become mixed by the impact.

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 {\mu}m, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

According to the standard giant impact hypothesis, the Moon formed from a partially vaporized disk generated by a collision between the proto Earth and a Mars sized impactor. The initial structure of the disk significantly affects the Moon forming process, including the Moons mass, its accretion time scale, and its isotopic similarity to Earth. The dynamics of the impact event determines the initial structure of a nearly hydrostatic Moon forming disk. However, the hydrostatic and hydrodynamic models have been studied separately and their connection has not previously been well quantified. Here, we show the extent to which the properties of the disk can be inferred from Smoothed Particle Hydrodynamic (SPH) simulations. By using entropy, angular momentum and mass distributions of the SPH outputs as approximately conserved quantities, we compute the two dimensional disk structure. We investigate four different models: (a) standard, the canonical giant impact model, (b) fast spinning Earth, a collision between a fast spinning Earth and a small impactor, (c) sub Earths, a collision between two objects with half Earths mass, and (d) intermediate, a collision of two bodies whose mass ratio is 7:3. Our SPH calculations show that the initial disk has approximately uniform entropy. The disks of the fast spinning Earth and sub Earths cases are hotter and more vaporized (80-90% vapor) than the standard case (20%). The intermediate case falls between these values. In the highly vaporized cases, our procedure fails to establish a unique surface density profile of the disk because the disk is unstable according to the Rayleigh criterion. In these cases, we estimate non-unique disk models by conserving global quantities. We also develop a semi analytic model for the thermal structure of the disk, which requires only two inputs: the average entropy and the surface density of the disk.

Titan's moment of inertia, calculated assuming hydrostatic equilibrium from gravity field data obtained during the Cassini-Huygens mission, implies an internal mass distribution that may be incompatible with complete differentiation. This suggests that Titan may have a mixed ice/rock core, possibly consistent with slow accretion in a gas-starved disk, which may initially spare Titan from widespread ice melting and subsequent differentiation. A partially differentiated Titan, however, must still efficiently remove radiogenic heat over geologic time. We argue that compositional heterogeneity in the major Saturnian satellites indicates that Titan formed from planetesimals with disparate densities. The resulting compositional anomalies would quickly redistribute to form a vertical density gradient that would oppose thermal convection. We use elements of the theory of double-diffusive convection to create a parameterized model for the thermal evolution of ice/rock mixtures with a stabilizing compositional gradient. Simulations are performed for a wide range of initial conditions to account for large uncertainties in material properties and accretionary processes. Ultimately, for realistic density gradients, double-diffusive convection in the ice/rock interior can delay, but not prevent, ice melting and differentiation, even if a substantial fraction of potassium is leached from the rock component. Consequently, Titan is not partially differentiated.

We compare the moment of inertia (MOI) of a simple hydrostatic, two layer body as determined by the Radau-Darwin Approximation (RDA) to its exact hydrostatic MOI calculated to first order in the parameter q = w^2R^3/GM, where w, R, and M are the spin angular velocity, radius, and mass of the body, and G is the gravitational constant. RDA is in error by less than 1% for many configurations of core sizes and layer densities congruent with those of solid bodies in the Solar System. We determine the error in the MOI of icy satellites calculated with the RDA due to nonhydrostatic effects by using a simple model in which the core and outer shell have slight degree 2 distortions away from their expected hydrostatic shapes. Since the hydrostatic shape has an associated stress of order pw^2R^2 (where p is density) it follows that the importance of nonhydrostatic effects scales with the dimensionless number s/pw^2R^2, where s is the nonhydrostatic stress. This highlights the likely importance of this error for slowly rotating bodies (e.g., Titan and Callisto) and small bodies (e.g., Saturn moons other than Titan). We apply this model to Titan, Callisto, and Enceladus and find that the RDA-derived MOI can be 10% greater than the actual MOI for nonhydrostatic stresses as small as ~0.1 bars at the surface or ~1 bar at the core-mantle boundary, while for Ganymede the stresses necessary to produce the same MOI errors are an order of magnitude greater due to its faster rotation. If satellites can reorient to the lowest energy state then RDA will always give an overestimate of the true MOI. Observations have shown that small nonhydrostatic gravity anomalies exist on Ganymede and Titan. We conclude that nonhydrostatic effects could be present to an extent that allows Callisto and Titan to be fully differentiated.

Through the process of thermal ionization, intense stellar irradiation renders Hot Jupiter atmospheres electrically conductive. Simultaneously, lateral variability in the irradiation drives the global circulation with peak wind speeds of order ~ km/s. In turn, the interactions between the atmospheric flows and the background magnetic field give rise to Lorentz forces that can act to perturb the flow away from its purely hydrodynamical counterpart. Using analytical theory and numerical simulations, here we show that significant deviations away from axisymmetric circulation are unstable in presence of a non-negligible axisymmetric magnetic field. Specifically, our results suggest that dayside-to-nightside flows, often obtained within the context of three-dimensional circulation models, only exist on objects with anomalously low magnetic fields, while the majority of highly irradiated exoplanetary atmospheres are entirely dominated by zonal jets.

Recently, the Kepler spacecraft has detected a sizable aggregate of objects, characterized by giant-planet-like radii and modest levels of stellar irradiation. With the exception of a handful of objects, the physical nature, and specifically the average densities, of these bodies remain unknown. Here, we propose that the detected giant planet radii may partially belong to planets somewhat less massive than Uranus and Neptune. Accordingly, in this work, we seek to identify a physically sound upper limit to planetary radii at low masses and moderate equilibrium temperatures. As a guiding example, we analyze the interior structure of the Neptune-mass planet Kepler-30d and show that it is acutely deficient in heavy elements, especially compared with its solar system counterparts. Subsequently, we perform numerical simulations of planetary thermal evolution and in agreement with previous studies, show that generally, 10 - 20 Earth mass, multi-billion year old planets, composed of high density cores and extended H/He envelopes can have radii that firmly reside in the giant planet range. We subject our results to stability criteria based on extreme ultraviolet radiation, as well as Roche-lobe overflow driven mass-loss and construct mass-radius relationships for the considered objects. We conclude by discussing observational avenues that may be used to confirm or repudiate the existence of putative low mass, gas-dominated planets.

We present calculations of thermal evolution of Hot Jupiters with various masses and effective temperatures under Ohmic dissipation. The resulting evolutionary sequences show a clear tendency towards inflated radii for effective temperatures that give rise to significant ionization of alkali metals in the atmosphere, compatible with the trend of the data. The degree of inflation shows that Ohmic dissipation, along with the likely variability in heavy element content can account for all of the currently detected radius anomalies. Furthermore, we find that in absence of a massive core, low-mass hot Jupiters can over-flow their Roche-lobes and evaporate on Gyr time-scales, possibly leaving behind small rocky cores.

The moment of inertia of a giant planet reveals important information about the planet's internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter's normalized moment of inertia NMoI=C/MR^2 by measuring Jupiter's pole precession and the Lense-Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter's mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J_2 and J_4 are perfectly known and are equal to their measured values. The model suggests that for fixed values of J_2 and J_4 a range of NMOI values between 0.2629 and 0.2645 can be found. The Radau-Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ~ 0.236, inferred from a dynamical model (Ward & Canup, 2006, ApJ, 640, L91) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter's NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter's internal structure. We carry out a simplified assessment of the error involved in Juno's possible determination of Jupiter's NMoI.

In order to accurately model giant planets, a whole set of observational constraints is needed. As the conventional constraints for extrasolar planets like mass, radius, and temperature allow for a large number of acceptable models, a new planetary parameter is desirable in order to further constrain planetary models. Such a parameter may be the tidal Love number k2. In this paper we aim to study the capability of k2 to reveal further information about the interior structure of a planet. With theoretical planetary models we investigate how the tidal Love number k2 responds to the internal density distribution of a planet. In particular, we demonstrate the effect of the degeneracy of k2 due to a density discontinuity in the envelope of a three-layer planetary model. The effect of a possible outer density discontinuity masks the effect of the core mass on the Love number k2. Hence, there is no unique relationship between the Love number k2 and the core mass of a planet. We show that the degeneracy of k2 with respect to a layer boundary in the envelope also occurs in existing planets, e.g. Saturn and the Hot Neptune GJ436b. As a result of the degeneracy, the planetary parameter k2 cannot be used to further constrain Saturnian models and for GJ436b only a maximum possible core mass can be derived from a given k2. To significantly narrow the uncertainty about the core mass of GJ436b the combined knowledge of k2 and atmospheric metallicity and temperature profile is necessary.

Despite its importance to questions of lunar origin, the chemical composition of the Moon is not precisely known. In recent years, however, the isotopic composition of lunar samples has been determined to high precision and found to be indistinguishable from the terrestrial mantle despite widespread isotopic heterogeneity in the Solar System. In the context of the giant-impact hypothesis, this level of isotopic homogeneity can evolve if the proto-lunar disk and post-impact Earth undergo turbulent mixing into a single uniform reservoir while the system is extensively molten and partially vaporized. In the absence of liquid-vapor separation, such a model leads to the lunar inheritance of the chemical composition of the terrestrial magma ocean. Hence, the turbulent mixing model raises the question of how chemical differences arose between the silicate Earth and Moon. Here we explore the consequences of liquid-vapor separation in one of the settings relevant to the lunar composition: the silicate vapor atmosphere of the post-giant-impact Earth. We use a model atmosphere to quantify the extent to which rainout can generate chemical differences by enriching the upper atmosphere in the vapor, and show that plausible parameters can generate the postulated enhancement in the FeO/MgO ratio of the silicate Moon relative to the Earth's mantle. Moreover, we show that liquid-vapor separation also generates measurable mass-dependent isotopic offsets between the silicate Earth and Moon and that precise silicon isotope measurements can be used to constrain the degree of chemical fractionation during this earliest period of lunar history. An approach of this kind has the potential to resolve long-standing questions on the lunar chemical composition.

Simulations of the moon-forming impact suggest that most of the lunar material derives from the impactor rather than the Earth. Measurements of lunar samples, however, reveal an oxygen isotope composition that is indistinguishable from terrestrial samples, and clearly distinct from meteorites coming from Mars and Vesta. Here we explore the possibility that the silicate Earth and impactor were compositionally distinct with respect to oxygen isotopes, and that the terrestrial magma ocean and lunar-forming material underwent turbulent mixing and equilibration in the energetic aftermath of the giant impact. This mixing may arise in the molten disk epoch between the impact and lunar accretion, lasting perhaps 10^2-10^3 yr. The implications of this idea for the geochemistry of the Moon, the origin of water on Earth, and constraints on the giant impact are discussed.

We present a new, magnetohydrodynamic mechanism for inflation of close-in giant extrasolar planets. The idea behind the mechanism is that current, which is induced through interaction of atmospheric winds and the planetary magnetic field, results in significant Ohmic dissipation of energy in the interior. We develop an analytical model for computation of interior Ohmic dissipation, with a simplified treatment of the atmosphere. We apply our model to HD209458b, Tres-4b and HD189733b. With conservative assumptions for wind speed and field strength, our model predicts a generated power that appears to be large enough to maintain the transit radii, opening an unexplored avenue towards solving a decade-old puzzle of extrasolar gas giant radius anomalies.

The atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists of two parts. The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) shows that the observed zonal winds cannot penetrate below a depth at which the electrical conductivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum penetration of 0.96R_J and 0.86R_S, with uncertainties of a few percent of R. At these radii, the magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent convection. They permit complete penetration along cylinders of the equatorial jets observed in the atmospheres of Jupiter and Saturn. The second part investigates how deep the observed zonal winds actually do penetrate. Truncation of the winds in the planet's convective envelope would involve breaking the Taylor-Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration which none of the obvious candidates appears able to provide. Accelerations arising from entropy gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These considerations suggest that strong zonal winds are confined to shallow, stably stratified layers, with equatorial jets being the possible exception.