Solar Cycle Effects on the Dynamics of Jupiter’s and Saturn’s Magnetospheres

The giant planetary magnetospheres surrounding Jupiter and Saturn respond in quite different ways, compared to Earth, to changes in upstream solar wind conditions. Spacecraft have visited Jupiter and Saturn during both solar cycle minima and maxima. In this paper we explore the large-scale structure of the interplanetary magnetic field (IMF) upstream of Saturn and Jupiter as a function of solar cycle, deduced from solar wind observations by spacecraft and from models. We show the distributions of solar wind dynamic pressure and IMF azimuthal and meridional angles over the changing solar cycle conditions, detailing how they compare to Parker predictions and to our general understanding of expected heliospheric structure at 5 and 9 AU. We explore how Jupiter’s and Saturn’s magnetospheric dynamics respond to varying solar wind driving over a solar cycle under varying Mach number regimes, and consider how changing dayside coupling can have a direct effect on the nightside magnetospheric response. We also address how solar UV flux variability over a solar cycle influences the plasma and neutral tori in the inner magnetospheres of Jupiter and Saturn, and estimate the solar cycle effects on internally driven magnetospheric dynamics. We conclude by commenting on the effects of the solar cycle in the release of heavy ion plasma into the heliosphere, ultimately derived from the moons of Jupiter and Saturn.     

On the thermal history of Saturn’s satellites Titan and Enceladus

The thermal histories of two geologically active satellites of Saturn—Titan and Enceladus—are discussed. During the Cassini mission, it was found that there are both nitrogen-containing compounds—NH₃ and N₂-and CO₂ and CH₄ in the water plumes of Enceladus; at that, ammonia is the prevailing form. This may testify that during evolution, the material of the satellite was warmed up to T ∼ 500–600 K, when NH₃ (the form of nitrogen capable of being accreted) could only be partly converted into N₂. Contrary to Enceladus, the temperature inside Titan probably reached values higher than 800 K or even higher than 1000 K, since the process of the chemical dissociation of ammonia was completely finished on this satellite and its atmosphere contains only molecular nitrogen. While the internal heating of Titan up to high temperatures can be explained by its large mass, the heating source for Enceladus’ interior is far from evident. Such traditional heating sources as the energy of gravitational differentiation and the radiogenic heating due to shortliving ²⁶Al and ⁶⁰Fe could not be effective. The first one is because of the small size of Enceladus (RE ≈ 250 km), and the inefficiency of the second one is caused by the fact that the satellite was formed not earlier than 8–10 Myr after the formation of calcium and aluminum-enriched inclusions in carbonaceous chondrites (CAI), i.e., after ²⁶Al had completely decayed. In the present paper, we propose other heating mechanisms-the heat of long-living radioactive elements and tidal heat, which could provide the observed chemical composition of the water plumes of Enceladus rather than only the differentiation of its protomatter into the ironstone core and the ice mantle.     

A Comparison of Solar Cycle Variations in the Equatorial Rotation Rates of the Sun’s Subsurface, Surface, Corona, and Sunspot Groups

Using the Solar Optical Observing Network (SOON) sunspot-group data for the period 1985?–?2010, the variations in the annual mean equatorial-rotation rates of the sunspot groups are determined and compared with the known variations in the solar equatorial-rotation rates determined from the following data: i) the plasma rotation rates at 0.94R_⊙,0.95R_⊙,…,1.0R_⊙ measured by the Global Oscillation Network Group (GONG) during the period 1995?–?2010, ii) the data on the soft-X-ray corona determined from Yohkoh/SXT full-disk images for the years 1992?–?2001, iii) the data on small bright coronal structures (SBCS) that were traced in Solar and Heliospheric Observatory (SOHO)/EIT images during the period 1998?–?2006, and iv) the Mount Wilson Doppler-velocity measurements during the period 1986?–?2007. A large portion (up to ≈?30^° latitude) of the mean differential-rotation profile of the sunspot groups lies between those of the internal differential-rotation rates at 0.94R_⊙ and 0.98R_⊙. The variation in the yearly mean equatorial-rotation rate of the sunspot groups seems to be lagging behind that of the equatorial-rotation rate determined from the GONG measurements by one to two years. The amplitude of the GONG measurements is very small. The solar-cycle variation in the equatorial-rotation rate of the solar corona closely matches that determined from the sunspot-group data. The variation in the equatorial-rotation rate determined from the Mount Wilson Doppler-velocity data closely resembles the corresponding variation in the equatorial-rotation rate determined from the sunspot-group data that included the values of the abnormal angular motions (>?|3^°|?day^?1) of the sunspot groups. Implications of these results are pointed out.     

On the Possibility of Gly and Ala Amino Acids on Titan’s Surface

In view of the possible production of formaldehyde and acetaldehyde in Titan’s atmosphere, the production of α-amino nitriles, the precursors of Glycine and Alanine amino acids, is explored in the upper atmosphere of Titan. The presence of Glycine and Alanine amino acids or their precursor amino nitriles can be used as a diagnostic, respectively for the presence or absence of water locally on the surface of Titan.     

A Simulation Study of the Plasma Wave Enhancements in the Earth’s Equatorial Plasmasphere

In the equatorial plasmasphere, plasma waves are frequently observed. To improve our understanding of the mechanism generating plasma waves from instabilities, a comparison of observations, linear growth-rate calculations, and simulation results is presented. To start the numerical experiments from realistic initial plasma conditions, we use the initial parameters inferred from observational data obtained around the plasma-wave generation region by the Akebono satellite. The linear growth rates of waves of different modes are calculated under resonance conditions, and compared with simulation results and observations. By employing numerical experiments by a particle code, we first show that upper hybrid-, Z-, and whistler-mode waves are excited through instabilities driven by a ring-type velocity distribution. The simulation results suggest a possibility that energetic electrons with energies of some tens of keV confined around the geomagnetic equator are responsible for the observed enhancements of Z- and whistler-mode waves. While the comparison between linear growth-rate calculations and observations shows the different tendency of wave amplitude of Z-mode and whistler-mode waves, the wave amplitude of these wave modes in the simulation results is consistent with the observation.     

A Comparison of Values of the Imaginary Part of Aerosol Refractive Index in the Latitudinal Belts of Saturn’s Northern Hemisphere

Using the spectrophotometric measurements data of 2015, the relation of values of the imaginary part n _i of aerosol refractive index was determined for latitudinal belts 17° N, 33° N, 49° N, and 66° N of Saturn’s disc. A steadily decreasing tendency in the relative n _i values when moving northward from the equatorial region of the disk to the latitude 49° N, inclusive, was revealed. The n _i values in the 17° N and 49° N belts were found to differ significantly from other latitudinal regions of the giant planet’s disk.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.     

Processes forming and sustaining Saturn’s proton radiation belts

Saturn’s proton radiation belts extend over the orbits of several moons that split this region of intense radiation into several distinct belts. Understanding their distribution requires to understand how their particles are created and evolve. High-energy protons are thought to be dominantly produced by cosmic ray albedo neutron decay (CRAND). The source of the lower energies and the role of other effects such as charge exchange with the gas originating from Enceladus is still an open question. There is also no certainty so far if the belts exist independently from each other and the rest of the magnetosphere or if and how particles are exchanged between these regions. We approach these problems by using measurements acquired by the MIMI/LEMMS instrument onboard the Cassini spacecraft. Protons in the range from 500 keV to 40 MeV are considered. Their intensities are averaged over 7 years of the mission and converted to phase space densities at constant first and second adiabatic invariant. We reproduce the resulting radial profiles with a numerical model that includes radial diffusion, losses from moons and interactions with gas, and a phenomenological source. Our results show that the dominating effects away from the moon sweeping corridors are diffusion and the source, while interactions with gas are secondary. Based on a GEANT4 simulation of the interaction of cosmic rays with Saturn’s rings, we conclude that secondary particles produced within the rings can only account for the high-energy part of the source. A comparison with the equivalent processes within Earth’s atmosphere shows that Saturn’s atmosphere can contribute to the production of the lower energies and might be even dominating at the higher energies. Other possibilities to supply the belts and exchange particles between them, as diffusion and injections from outside the belts, or stripping of ENAs, can be excluded.     

The origin and evolution of Saturn’s 2011–2012 stratospheric vortex

The planet-encircling springtime storm in Saturn’s troposphere (December 2010–July 2011, Fletcher, L.N. et al. [2011]. Science 332, 1413–1414; Sánchez-Lavega, A. et al. [2011]. Nature 475, 71–74; Fischer, G. et al. [2011]. Nature 475, 75–77) produced dramatic perturbations to stratospheric temperatures, winds and composition at mbar pressures that persisted long after the tropospheric disturbance had abated. Thermal infrared (IR) spectroscopy from the Cassini Composite Infrared Spectrometer (CIRS), supported by ground-based IR imaging from the VISIR instrument on the Very Large Telescope and the MIRSI instrument on NASA’s IRTF, is used to track the evolution of a large, hot stratospheric anticyclone between January 2011 and March 2012. The evolutionary sequence can be divided into three phases: (I) the formation and intensification of two distinct warm airmasses near 0.5 mbar between 25 and 35°N (B1 and B2) between January–April 2011, moving westward with different zonal velocities, B1 residing directly above the convective tropospheric storm head; (II) the merging of the warm airmasses to form the large single ‘stratospheric beacon’ near 40°N (B0) between April and June 2011, disassociated from the storm head and at a higher pressure (2 mbar) than the original beacons, a downward shift of 1.4 scale heights (approximately 85 km) post-merger; and (III) the mature phase characterised by slow cooling (0.11 ± 0.01 K/day) and longitudinal shrinkage of the anticyclone since July 2011. Peak temperatures of 221.6 ± 1.4 K at 2 mbar were measured on May 5th 2011 immediately after the merger, some 80 K warmer than the quiescent surroundings. From July 2011 to the time of writing, B0 remained as a long-lived stable stratospheric phenomenon at 2 mbar, moving west with a near-constant velocity of 2.70 ± 0.04 deg/day (?24.5 ± 0.4 m/s at 40°N relative to System III longitudes). No perturbations to visible clouds and hazes were detected during this period.With no direct tracers of motion in the stratosphere, we use thermal windshear calculations to estimate clockwise peripheral velocities of 200–400 m/s at 2 mbar around B0. The peripheral velocities of the two original airmasses were smaller (70–140 m/s). In August 2011, the size of the vortex as defined by the peripheral collar was 65° longitude (50,000 km in diameter) and 25° latitude. Stratospheric acetylene (C₂H₂) was uniformly enhanced by a factor of three within the vortex, whereas ethane (C₂H₆) remained unaffected. The passage of B0 generated a new band of warm stratospheric emission at 0.5 mbar at its northern edge, and there are hints of warm stratospheric structures associated with the beacons at higher altitudes (p < 0.1 mbar) than can be reliably observed by CIRS nadir spectroscopy. Analysis of the zonal windshear suggests that Rossby wave perturbations from the convective storm could have propagated vertically into the stratosphere at this point in Saturn’s seasonal cycle, one possible source of energy for the formation of these stratospheric anticyclones.     

ABCAnalytical and Semi-Analytical Analysis of an Artificial Satellite’s Rotational Motion

The rotational motion of an artificial satellite is studied by considering torques produced by gravity gradient and direct solar radiation pressure. A satellite of circular cylinder shape is considered here, and Andoyers variables are used to describe the rotational motion. Expressions for direct solar radiation torque are derived. When the earths shadow is not considered, an analytical solution is obtained using Lagranges method of variation of parameters. A semi-analytical procedure is proposed to predict the satellites attitude under the influence of the earths shadow. The analytical solution shows that angular variables are linear and periodic functions of time while their conjugates suffer only periodic variations. When compared, numerical and analytical solutions have a good agreement during the time range considered.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation when the two planets are in migration

In a previous paper (Hou et al. in Celest Mech Dyn Astron 119:119–142, 2014a), the problem of dynamical symmetry between two Jupiter triangular libration points (TLPs) with Saturn’s perturbation in the present configuration of the two planets was studied. A small short-time scale spatial asymmetry exists but gradually disappears with the time going, so the planar stable regions around the two Jupiter TLPs should be dynamically symmetric from a longtime perspective. In this paper, the symmetry problem is studied when the two planets are in migration. Several mechanisms that can cause asymmetries are discussed. Studies show that three important ones are the large short-time scale spatial asymmetry when Jupiter and Saturn are in resonance, the changing orbits of Jupiter and Saturn in the planet migration process, and the chaotic nature of Trojan orbits during the planet migration process. Their joint effects can cause an observable difference to the two Jupiter Trojan swarms. The thermal Yarkovsky effect is also found to be able to cause dynamical differences to the two TLPs, but generally they are too small to be practically observed.     

Photographic observations of Saturn’s moons at the MAO NAS of Ukraine in 1961–1990

A catalog of 1385 astrometric positions of Saturn’s moons S2–S9 has been compiled with Tycho-2 as a reference frame from photographic observations obtained at the Main Astronomical Observatory, National Academy of Sciences of Ukraine, in 1961–1990. Astronegatives have been digitized with an Epson Expression 10000XL commercial scanner in 16-bit grayscale with a resolution of 1200 dpi. Reduction has been performed in the LINUX-MIDAS-ROMAFOT software supplemented with additional modules. The internal positional accuracy of the reduction is 0.09…0.23′′ for both coordinates and 0.27…0.37^m for the photographic magnitudes of the Tycho-2 catalog. The calculated topocentric positions of the moons are compared online with the IMCCE ephemeris data (DE405 + TASS1,7). Moon-minus-moon differential coordinates are found for most of the moons and compared with theoretical data (http://lnfm1.sai.msu.ru/neb/nss/nssephmr.htm).     

Optimization of the Return Trajectories from the Moon for Delivering Soil to a Specified Region on the Earth’s Surface

Solar System Research - The construction of the trajectories for spacecraft return from the Moon to the Earth for delivering lunar soil to the polygon specified by the corresponding geographical...     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.     

The effect of the Earth’s oblateness on the Moon’s physical libration in latitude

The Moon’s physical libration in latitude generated by gravitational forces caused by the Earth’s oblateness has been examined by a vector analytical method. Libration oscillations are described by a close set of five linear inhomogeneous differential equations, the dispersion equation has five roots, one of which is zero. A complete solution is obtained. It is revealed that the Earth’s oblateness: a) has little effect on the instantaneous axis of Moon’s rotation, but causes an oscillatory rotation of the body of the Moon with an amplitude of 0.072″ and pulsation period of 16.88 Julian years; b) causes small nutations of poles of the orbit and of the ecliptic along tight spirals, which occupy a disk with a cut in a center and with radius of 0.072″. Perturbations caused by the spherical Earth generate: a) physical librations in latitude with an amplitude of 34.275″; b) nutational motion for centers of small spiral nutations of orbit (ecliptic) pole over ellipses with semi-major axes of 113.850″ (85.158″) and the first pole rotates round the second one along a circle with radius of 28.691″; c) nutation of the Moon’s celestial pole over an ellipse with a semi-major axis of 45.04″ and with an axes ratio of about 0.004 with a period of T = 27.212 days. The principal ellipse’s axis is directed tangentially with respect to the precession circumference, along which the celestial pole moves nonuniformly nearly in one dimension. In contrast to the accepted concept, the latitude does not change while the Moon’s poles of rotation move. The dynamical reason for the inclination of the Moon’s mean equator with respect to the ecliptic is oblateness of the body of the Moon.     

A predator–prey model for moon-triggered clumping in Saturn’s rings

UVIS occultation data show clumping in Saturn’s F ring and at the B ring outer edge, indicating aggregation and disaggregation at these locations that are perturbed by Prometheus and by Mimas. The inferred timescales range from hours to months. Occultation profiles of the edge show wide variability, indicating perturbations by local mass aggregations. Structure near the B ring edge is seen in power spectral analysis at scales 200–2000 m. Similar structure is also seen at the strongest density waves, with significance increasing with resonance strength. For the B ring outer edge, the strongest structure is seen at longitudes 90° and 270° relative to Mimas. This indicates a direct relation between the moon and the ring clumping. We propose that the collective behavior of the ring particles resembles a predator–prey system: the mean aggregate size is the prey, which feeds the velocity dispersion; conversely, increasing dispersion breaks up the aggregates. Moons may trigger clumping by streamline crowding, which reduces the relative velocity, leading to more aggregation and more clumping. Disaggregation may follow from disruptive collisions or tidal shedding as the clumps stir the relative velocity. For realistic values of the parameters this yields a limit cycle behavior, as for the ecology of foxes and hares or the “boom-bust” economic cycle. Solving for the long-term behavior of this forced system gives a periodic response at the perturbing frequency, with a phase lag roughly consistent with the UVIS occultation measurements. We conclude that the agitation by the moons in the F ring and at the B ring outer edge drives aggregation and disaggregation in the forcing frame. This agitation of the ring material may also allow fortuitous formation of solid objects from the temporary clumps, via stochastic processes like compaction, adhesion, sintering or reorganization that drives the denser parts of the aggregate to the center or ejects the lighter elements. Any of these more persistent objects would then orbit at the Kepler rate. We would also expect the formation of clumps and some more permanent objects at the other perturbed regions in the rings… including satellite resonances, shepherded ring edges, and near embedded objects like Pan and Daphnis (where the aggregation/disaggregation cycles are forced similar to Prometheus forcing of the F ring).     

Atronomical CCD observations of the main Saturn’s satellites at Pulkovo Observatory in 2004–2007

The results of astrometric observations of Saturn’s satellites (S1–S8) obtained using a 26-inch refractor and a normal astrograph at Pulkovo Observatory in 2004–2007 are given. High-accuracy equatorial coordinates of Saturn’s satellites in the system of the UCAC2 reference catalog and the relative “satellite-satellite” positions suitable for specifying their motion theories are obtained. The observations are compared with the DE405 + TASS1.7 and INPOP06 + TASS1.7 theories of motion. The root-mean-square errors of the obtained satellite positions lie within the range of 10–50 mas, as far as the intrinsic convergence is concerned, and 20–70 mas, as far as the extrinsic one is concerned. The observation results are included into the astrometrical database of the Pulkovo Observatory (www.puldb.ru).     

Statistical interrelation of 22-yr-long variations in the data on parameters of the Earth’s rotation and magnetic fields of the Earth and the Sun

By means of spectral analysis, oscillations have been detected in many-year time series of deviations in the duration of days from the standard that cannot be explained within the framework of existing gravitational theory. The solution of the problem of the origin and structure of these oscillations is associated with the essence of the phenomena taking place inside the Earth and the mechanisms of energy transfer by external fields and the Sun. The effect of the Sun on the Earth’s rotation also leads to the formation of a unified cyclic background, resulting in correlated oscillations in all shells of the Earth, as well as in its atmosphere and its nucleus. Ground magnetic fields of the Earth and the Sun play the leading role in the abovementioned unified cyclic process. The results of spectral autoregressive and wavelet analyses of experimental data concerning deviations in the duration of days from the standard in the years 1832–2006, ground geomagnetic field intensities in 1832–2006, and the mean magnetic field of the Sun during the period from 1975 to 2005 have been discussed in this paper to reveal and compare correlating oscillations. To analyze a short-period (31-yr-long) series of daily data on the mean magnetic field of the Sun, the results of wavelet transformations (the Morlet wavelet) of the detected amplitude of a burst envelope with a carrier frequency of 13 cycles/yr have been obtained.     

Interpretation of the 6818.9-Å methane feature observed on Jupiter,Saturn, and Uranus

High-resolution (0.1-Å) spectra of the 6818.9-Å methane feature obtained for Jupiter, Saturn, and Uranus by K. H. Baines, W. V. Schempp, and W. H. Smith ((1983). Icarus56, 534–542) are modeled using a doubling and adding code after J. H. Hansen ((1969). Astrophys. J.155, 565–573). The feature's rotational quantum number is estimated using the relatively homogeneous atmosphere of Saturn, with only J = 0 and J = 1 fitting the observational constraints. The aerosol content within Saturn's northern temperate region is shown to be substantially less than at the equator, indicating a haze only half as optically thick. Models of Jupiter's atmosphere are consistent with the rotational quantum-number assignment. Synthetic line profiles of the 6818.9-Å feature observed on Uranus reveal that a substantial haze exists at or above the methane condensation region with an optical depth eight times greater than previously reported. Seasonal effects are indicated. The methane column abundance is 5 ± 1 km-am. The mixing ratio of methane to hydrogen within the deep unsaturated region of the planet is 0.045 ± 0.025, based on an H₂ column abundance of 240 ± 60 km-am (W. H. Smith, W. Macy, and C. B. Pilcher (1980). Icarus43, 153–160), thus indicating that the methane comprises between one-sixth and one-half of the planet's mass. However, proper reevaluation of H₂ quadrupole features accounting for the haze reported here may significantly reduce the relative methane abundance.