This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Observations of active M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. We detail the analogy between some anelastic dynamo simulations and spectropolarimetric observations of 23 M stars. In numerical models, the relative contribution of inertia and Coriolis force in the global force balance -estimated by the so-called local Rossby number- is known to have a strong impact on the magnetic field geometry. We discuss the relevance of this parameter in setting the large-scale magnetic field of M dwarfs.

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.

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.

Context: observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance -- quantified by the local Rossby number -- is known to have a strong impact on the magnetic field geometry. Aims: we aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. Methods: we explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to reflect the diagnostic parameters employed in numerical models. Results: the transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby number. By analogy with numerical models, we expect different amplitudes of differential rotation on the two dynamo branches.

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.

A didactic introduction to current thinking on some aspects of the solar dynamo is given for geophysicists and planetary scientists.

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.

A comparison is made between mean-field models and direct numerical simulations of rotating magnetoconvection and the geodynamo. The mean-field coefficients are calculated with the fluid velocity taken from the direct numerical simulations. The magnetic fields resulting from mean-field models are then compared with the mean magnetic field from the direct numerical simulations.