Search

The outer areas of Jupiter and Saturn have multiple zonal winds, reaching the high latitudes, that penetrate deep into the planets' interiors, as suggested by gravity measurements. These characteristics are replicable in numerical simulations by including both a shallow stably stratified layer, below a convecting envelope, and increasing electrical conductivity. A dipolar magnetic field, assumed to be generated by a dynamo below our model, is imposed. We find that the winds' depth into the stratified layer depends on the local product of the squared magnetic field strength and electrical conductivity. The key for the drop-off of the zonal winds is a meridional circulation which perturbs the density structure in the stable layer. In the stable region its dynamics is governed by a balance between Coriolis and electromagnetic forces. Our models suggest that a stable layer extending into weakly conducting regions could account for the observed deep zonal wind structures., Comment: Accepted by JGR - Planets

The last 22.5 orbits of the Cassini mission brought the spacecraft to less than 3000 km from Saturn's 1-bar surface. These close encounters offered an unprecedented view of Saturn's magnetic field, including contributions from the internal dynamo, the ionosphere, and the magnetosphere. In this chapter, we highlight the new picture of Saturn's magnetic field from the Cassini mission including the persistent yet time-varying low-latitude field-aligned currents, Alfv\'en waves planet-ward of the D-ring, extreme axisymmetry, and high-degree magnetic moments. We then discuss the implications and new questions raised for Saturn's innermost magnetosphere, equatorial ionosphere, and interior. We conclude this chapter with an outlook for the future exploration of Saturn and other giant planets., Comment: Peer reviewed and accepted for publication as Chapter 5 of a multi-volume work edited by Kevin Baines, Michael Flasar, Norbert Krupp, and Thomas Stallard, entitled Cassini at Saturn: The Grand Finale, to be published by Cambridge University Press. 15 pages, 10 figures

The precise mechanism that forms jets and large-scale vortices on the giant planets is unknown. An inverse cascade has been suggested. Alternatively, energy may be directly injected by small-scale convection. Our aim is to clarify whether an inverse cascade feeds zonal jets and large-scale eddies in a system of rapidly rotating, deep, geostrophic spherical-shell convection. We analyze the nonlinear scale-to-scale transfer of kinetic energy in such simulations as a function of the azimuthal wave number, m. We find that the main driving of the jets is associated with upscale transfer directly from the small convective scales to the jets. This transfer is very nonlocal in spectral space, bypassing large-scale structures. The jet formation is thus not driven by an inverse cascade. Instead, it is due to a direct driving by Reynolds stresses from small-scale convective flows. Initial correlations are caused by the effect of uniform background rotation and shell geometry on the flows. While the jet growth suppresses convection, it increases the correlation of the convective flows, which further amplifies the jet growth until it is balanced by viscous dissipation. To a much smaller extent, energy is transferred upscale to large-scale vortices directly from the convective scales, mostly outside the tangent cylinder. There, large-scale vortices are not driven by an inverse cascade either. Inside the tangent cylinder, the transfer to large-scale vortices is weaker, but more local in spectral space, leaving open the possibility of an inverse cascade as a driver of large-scale vortices. In addition, large-scale vortices receive kinetic energy from the jets via forward transfer. We therefore suggest a jet instability as an alternative formation mechanism of largescale vortices. Finally, we find that the jet kinetic energy scales as $\ell^{-5}$, the same as for the zonostrophic regime., Comment: 15 pages, 14 figures, accepted for publication in A&A

The recent gravity field measurements of Jupiter (Juno) and Saturn (Cassini) confirm the existence of deep zonal flows reaching to a depth of 5\% and 15\% of the respective radius. Relating the zonal wind induced density perturbations to the gravity moments has become a major tool to characterise the interior dynamics of gas giants. Previous studies differ with respect to the assumptions made on how the wind velocity relates to density anomalies, on the functional form of its decay with depth, and on the continuity of antisymmetric winds across the equatorial plane. Most of the suggested vertical structures exhibit a rather smooth radial decay of the zonal wind, which seems at odds with the observed secular variation of the magnetic field and the prevailing geostrophy of the zonal winds. Moreover, the results relied on an artificial equatorial regularisation or ignored the equatorial discontinuity altogether. We favour an alternative structure, where the equatorially antisymmetric zonal wind in an equatorial latitude belt between $\pm 21^\circ$ remains so shallow that it does not contribute to the gravity signal. The winds at higher latitudes suffice to convincingly explain the measured gravity moments. Our results indicate that the winds are geostrophic, i.e. constant along cylinders, in the outer $3000\,$ km and decay rapidly below. The preferred wind structure is 50\% deeper than previously thought, agrees with the measured gravity moment, is compliant with the magnetic constraints and the requirement of an adiabatic atmosphere and unbiased by the treatment of the equatorial discontinuity.

Recent precise measurements at Jupiter's and Saturn's gravity fields constrain the properties of the zonal flows in the outer envelopes of these planets. A simplified dynamic equation, sometimes called the thermal wind or thermo-gravitational wind equation, establishes a link between zonal flows and the related buoyancy perturbation, which in turn can be exploited to yield the dynamic gravity perturbation. Whether or not the action of the dynamic gravity perturbation needs to be explicitly included in this equation, an effect we call the Dynamic Self Gravity (DSG), has been a matter of intense debate. We show that, under reasonable assumptions, the equation can be solved (semi) analytically. This allows us to quantify the impact of the DSG on each gravity harmonic, practically independent of the zonal flow or the details of the planetary interior model. The impact decreases with growing spherical harmonic degree l. For degrees l=2 to about l=4, the DSG is a first order effect and should be taken into account in any attempt of inverting gravity measurements for zonal flow properties. For degrees of about l=5 to roughly l=10, the relative impact of DSG is about 10% and thus seems worthwhile to include, in particular since this comes at little extra costs with the method presented here. For yet higher degrees, is seems questionable whether gravity measurements or interior models will ever reach the required precision equivalent of the DSG impact of only a few percent of less., Comment: 28 pages, 6 figures

The Sun's complex dynamics is controlled by buoyancy and rotation in the convection zone and by magnetic forces in the atmosphere and corona. While small-scale solar convection is well understood, the dynamics of large-scale flows in the solar convection zone is not explained by theory or simulations. Waves of vorticity due to the Coriolis force, known as Rossby waves, are expected to remove energy out of convection at the largest scales. Here we unambiguously detect and characterize retrograde-propagating vorticity waves in the shallow subsurface layers of the Sun at angular wavenumbers below fifteen, with the dispersion relation of textbook sectoral Rossby waves. The waves have lifetimes of several months, well-defined mode frequencies below 200 nHz in a co-rotating frame, and eigenfunctions of vorticity that peak at the equator. Rossby waves have nearly as much vorticity as the convection at the same scales, thus they are an essential component of solar dynamics. We find a transition from turbulence-like to wave-like dynamics around the Rhines scale of angular wavenumber of twenty; this might provide an explanation for the puzzling deficit of kinetic energy at the largest spatial scales., Comment: This is the submitted version of the paper published in Nature Astronomy. 23 pages, 8 figures, 1 table

Earth sustains its magnetic field by a dynamo process driven by convection in the liquid outer core. Geodynamo simulations have been successful in reproducing many observed properties of the geomagnetic field. However, while theoretical considerations suggest that flow in the core is governed by a balance between Lorentz force, rotational force and buoyancy (called MAC balance for Magnetic, Archimedean, Coriolis) with only minute roles for viscous and inertial forces, dynamo simulations must employ viscosity values that are many orders of magnitude larger than in the core due to computational constraints. In typical geodynamo models viscous and inertial forces are not much smaller than the Coriolis force and the Lorentz force plays a sub-dominant role. This has led to conclusions that these simulations are viscously controlled and do not represent the physics of the geodynamo. Here we show by a direct analysis of the relevant forces that a MAC balance can be achieved when the viscosity is reduced to values close to the current practical limit. Lorentz force, buoyancy and the uncompensated (by pressure) part of the Coriolis force are of very similar strength, whereas viscous and inertia are smaller by a factor of at least 20 in the bulk of the fluid volume. Compared to non-magnetic convection at otherwise identical parameters, the dynamo flow is of larger scale, less invariant parallel to the rotation axis (less geostrophic) and convection transports twice as much heat, all of which is expected when the Lorentz force strongly influences the convection properties., Comment: Accepted for publication in the Proceedings of the National Academy of Sciences (PNAS)

The recent discovery of an Earth-like exoplanet around Proxima Centauri has shined a spot light on slowly rotating fully convective M-stars. When such stars rotate rapidly (period $\lesssim 20$ days), they are known to generate very high levels of activity that is powered by a magnetic field much stronger than the solar magnetic field. Recent theoretical efforts are beginning to understand the dynamo process that generates such strong magnetic fields. However, the observational and theoretical landscape remains relatively uncharted for fully convective M-stars that rotate slowly. Here we present an anelastic dynamo simulation designed to mimic some of the physical characteristics of Proxima Centauri, a representative case for slowly rotating fully convective M-stars. The rotating convection spontaneously generates differential rotation in the convection zone which drives coherent magnetic cycles where the axisymmetric magnetic field repeatedly changes polarity at all latitudes as time progress. The typical length of the `activity' cycle in the simulation is about nine years, in good agreement with the recently proposed activity cycle length of about seven years for Proxima Centauri. Comparing our results with earlier work, we hypothesis that the dynamo mechanism undergoes a fundamental change in nature as fully convective stars spin down with age., Comment: 8 pages, 4 figures, double column; Accepted in ApJ Letters

Despite the lack of a shear-rich tachocline region low-mass fully convective stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three dimensional model of magnetic field generation in low-mass fully convective stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman-Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass fully convective stars., Comment: Double column format, 6 pages, 5 figures. To appear in ApJL

We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime., Comment: 15 pages, double column format, 10 figures. Substantial modifications in version 2. Data for shells with aspect ratio 0.35 ("Supple_data") can be found in the source. To appear in "Geophysical Journal International"

We present an analysis of the olivine-rich exposures at Bellicia and Arruntia craters using Dawn Framing Camera (FC) color data. Our results confirm the existence of olivine-rich materials at these localities as described by Ammannito et al. (2013a) using Visual Infrared Spectrometer (VIR) data. Analyzing laboratory spectra of various Howardite-Eucrite-Diogenite meteorites, high-Ca pyroxenes, olivines and olivine-orthopyroxene mixtures, we derive three FC spectral band parameters that are indicators of olivine-rich materials. Combining the three band parameters allows us, for the first time, to reliably identify sites showing modal olivine contents >40%. The olivine-rich exposures at Bellicia and Arruntia are mapped using higher spatial resolution FC data. The exposures are located on the slopes of outer/inner crater walls, on the floor of Arruntia, in the ejecta, as well as in nearby fresh small impact craters. The spatial extent of the exposures ranges from a few hundred meters to few kilometers. The olivine-rich exposures are in accordance with both the magma ocean and the serial magmatism model (e.g., Righter and Drake 1997; Yamaguchi et al. 1997). However, it remains unsolved why the olivine-rich materials are mainly concentrated in the northern hemisphere (~36-42{\deg} N, 46-74{\deg} E) and are almost absent in the Rheasilvia basin., Comment: Accepted for publication in Meteoritics and Planetary Science

Observations of cool stars reveal dark spot-like features on their surfaces. Compared to sunspots, starspots can be bigger or cover a larger fraction of the stellar surface. While sunspots appear only at low latitudes, starspots are also found in polar regions, in particular on rapidly rotating stars. Sunspots are believed to result from the eruption of magnetic flux-tubes rising from the deep interior of the Sun. The strong magnetic field locally reduces convective heat transport to the solar surface. Such flux-tube models have also been invoked to explain starspot properties. However, these models use several simplifications and so far the generation of either sunspots or starspots has not been demonstrated in a self-consistent simulation of stellar magnetic convection. Here we show that direct numerical simulations of a distributed dynamo operating in a density-stratified rotating spherical shell can spontaneously generate cool spots. Convection in the interior of the model produces a large scale magnetic field which interacts with near surface granular convection leading to strong concentrations of magnetic flux and formation of starspots. Prerequisites for the formation of sizeable high-latitude spots in the model are sufficiently strong density stratification and rapid rotation. Our model presents an alternate mechanism for starspot formation by distributed dynamo action., Comment: 14 pages; Important additions in version 2; To appear in A&A

Dynamos operating in the interiors of rapidly rotating planets and low-mass stars might belong to a similar category where rotation plays a vital role. We quantify this similarity using scaling laws. We analyse direct numerical simulations of Boussinesq and anelastic spherical shell dynamos. These dynamos represent simplified models which span from Earth-like planets to rapidly rotating low-mass stars. We find that magnetic field and velocity in these dynamos are related to the available buoyancy power via a simple power law which holds over wide variety of control parameters., Comment: 2 pages; Proceedings of IAUS 302: Magnetic fields throughout stellar evolution (August 2013, Biarritz, France)

Numerical dynamo models always employ parameter values that differ by orders of magnitude from the values expected in natural objects. However, such models have been successful in qualitatively reproducing properties of planetary and stellar dynamos. This qualitative agreement fuels the idea that both numerical models and astrophysical objects may operate in the same asymptotic regime of dynamics. This can be tested by exploring the scaling behavior of the models. For convection-driven incompressible spherical shell dynamos with constant material properties, scaling laws had been established previously that relate flow velocity and magnetic field strength to the available power. Here we analyze 273 direct numerical simulations using the anelastic approximation, involving also cases with radius-dependent magnetic, thermal and viscous diffusivities. These better represent conditions in gas giant planets and low-mass stars compared to Boussinesq models. Our study provides strong support for the hypothesis that both mean velocity and mean magnetic field strength scale as a function of power generated by buoyancy forces in the same way for a wide range of conditions., Comment: 9 pages, 4 figures, 1 table; data used in the paper can be found in "Dataset.txt" file available in the source; to appear in ApJ

Numerical simulations of convection driven rotating spherical shell dynamos have often been performed with rigid boundary conditions, as is appropriate for the metallic cores of terrestrial planets. Free-slip boundaries are more appropriate for dynamos in other astrophysical objects, such as gas-giants or stars. Using a set of 57 direct numerical simulations, we investigate the effect of free-slip boundary conditions on the scaling properties of heat flow, flow velocity and magnetic field strength and compare it with earlier results for rigid boundaries. We find that the nature of the mechanical boundary condition has only a minor influence on the scaling laws. We also find that although dipolar and multipolar dynamos exhibit approximately the same scaling exponents, there is an offset in the scaling pre-factors for velocity and magnetic field strength. We argue that the offset can be attributed to the differences in the zonal flow contribution between dipolar and multipolar dynamos., Comment: 10 pages, 9 figures, 1 table. To appear in ICARUS

Very little is known about magnetic fields of extrasolar planets and brown dwarfs. We use the energy flux scaling law presented by Christensen et al. (2009) to calculate the evolution of average magnetic fields in extrasolar planets and brown dwarfs under the assumption of fast rotation, which is probably the case for most of them. We find that massive brown dwarfs of about 70 M_Jup can have fields of a few kilo-Gauss during the first few hundred Million years. These fields can grow by a factor of two before they weaken after deuterium burning has stopped. Brown dwarfs with weak deuterium burning and extrasolar giant planets start with magnetic fields between ~100G and ~1kG at the age of a few Myr, depending on their mass. Their magnetic field weakens steadily until after 10Gyr it has shrunk by about a factor of 10. We use observed X-ray luminosities to estimate the age of the known extrasolar giant planets that are more massive than 0.3M_Jup and closer than 20pc. Taking into account the age estimate, and assuming sun-like wind-properties and radio emission processes similar to those at Jupiter, we calculate their radio flux and its frequency. The highest radio flux we predict comes out as 700mJy at a frequency around 150MHz for $\tau$Boob, but the flux is below 60mJy for the rest. Most planets are expected to emit radiation between a few Mhz and up to 100MHz, well above the ionospheric cutoff frequency., Comment: 7 pages, accepted by A&A

We have measured the surface magnetic flux on four accreting young brown dwarfs and one non-accreting young very low-mass star utilizing high resolution spectra of absorption lines of the FeH molecule. A magnetic field of 1-2 kG had been proposed for one of the brown dwarfs, 2MASS J1207334$-$393254, because of its similarities to higher mass T Tauri stars as manifested in accretion and the presence of a jet. We do not find clear evidence for a kilo-Gauss field in any of our young brown dwarfs but do find a 2 kG field on the young VLM star. Our 3-$\sigma$ upper limit for the magnetic flux in 2MASS J1207334$-$393254 just reaches 1 kG. We estimate the magnetic field required for accretion in young brown dwarfs given the observed rotations, and find that fields of only a few hundred Gauss are sufficient for magnetospheric accretion. This predicted value is less than our observed upper limit. We conclude that magnetic fields in young brown dwarfs are a factor of five or more lower than in young stars of about one solar mass, and in older stars with spectral types similar to our young brown dwarfs. It is interesting that, during the first few million years, the fields scale down with mass in line with what is needed for magnetospheric accretion, yet no such scaling is observed at later ages within the same effective temperature range. This scaling is opposite to the trend in rotation, with shorter rotation periods for very young accreting brown dwarfs compared with accreting solar-mass objects (and very low Rossby numbers in all cases). We speculate that in young objects a deeper intrinsic connection may exist between magnetospheric accretion and magnetic field strength, or that magnetic field generation in brown dwarfs may be less efficient than in stars. Neither of these currently have an easy physical explanation., Comment: accepted by ApJ

The six sensors of SEIS (The Seismic Experiment for Interior Structure) [- one of three primary instruments on NASA's Mars Lander Insight] cover a broad range of the seismic bandwidth, from 0.01 hertz to 50 hertz, with possible extension to longer periods. Data are transmitted in the form of three continuous VBB (Very Broad-Band) components at 2 samples per second (sps), an estimation of the short period (SP) energy content from the SP at 1 sps, and a continuous compound VBB/SP vertical axis at 10 sps. The continuous streams are augmented by requested event data with sample rates from 20 to 100 sps. SEIS data products are downlinked from the spacecraft in raw CCSDS (Consultative Committee for Space Data Systems) packets and converted to both the Standard for the Exchange of Earthquake Data (SEED) format files and ASCII tables (GeoCSV) for analysis and archiving. Metadata are available in dataless SEED and StionXML. Time series data (waveforms) are available in miniseed and GeoCSV. Data are distributed according to FDSN (Federation of Digital Seismograph Networks - http://www.fdsn.org) formats and interfaces. Wind, pressure and temperature data from the Auxiliary Payload Sensor Suite (APSS) will also be available in SEED format, and can be used for decorrelation and diagnostic purposes on SEIS.