Results of astrometric observations of Jupiter’s Galilean satellites at the Pulkovo Observatory from 1986 to 2005

The results of photographic observations of Jupiter’s Galilean satellites made with the 26-inch refractor at the Pulkovo Observatory from 1986 to 2005 are given. Satellite coordinates with respect to Jupiter and the mutual distances between the satellites have been determined. A scale-trale technique that does not require reference stars for the astrometric reduction of measurements has been used. The effect of the Jupiter phase has been taken into account in the jovicentric coordinates. The observation results have been compared with a modern theory of the Galilean satellites’ motions. Systematic observation errors depending on the observation technique have been studied. The intrinsic observation accuracy in the random quotient is characterized by the values 0.041″ over X and Y. The external accuracy of the relative Galilean satellite coordinates determined by comparing the observations with modern ephemerides turned out to be equal to 0.165″, 0.213″ for the Jovicentric coordinates and 0.134″, 0.170″ for the “satellite-satellite” coordinates. The highest accuracy of the relative satellite coordinates is reached at small distances between the satellites which are less than 100″: the corresponding mean-square errors of one observation are equal in to the external convergence to 0.050″, 0.070″. The results of photographic observations have been compared with the first CCD observations of the Jupiter satellites made in 2004 with the 26-inch refractor.     

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).     

Astrometric results of observations at Russian observatories of mutual occultations and eclipses of Jupiter’s Galilean satellites in 2009

In 2009, in five Russian observatories photometric observations of Jupiter’s Galilean satellites during their mutual occultations and eclipses were carried out. Based on these observations, an original method was used to ascertain astrometric results such as the difference between the coordinates of pairs of satellites. Fifty-three phenomena were successfully observed. A total of 94 light curves of satellites were measured. The error in the coordinates of satellites due to random errors in photometry, calculated on all data obtained, was 0.041″ in right ascension and 0.046″ in declination. The discrepancies between the theory and observations in these coordinates was found to be 0.060″ and 0.057″, respectively. The results were uploaded to the common database for all observations of natural satellites of planets at the Natural Satellites Data Center (NSDC), which is available online at http://www.sai.msu.ru/neb/nss/index.htm. For the first time in the practice of photometric observations of satellites in epochs of mutual occultations and eclipses a new method of observation was tested, which eliminates from astrometric results the major systematic errors caused by an inaccurate account of the background level. The tests were conducted in the Terskol Observatory and the observatory of the Crimean laboratory of the Sternberg State Astronomical Institute of the Moscow State University. The application of the new method showed that the elimination of the background level at these observatories was carried out correctly.     

Hard X-ray and γ-ray observations of solar flares with the Phebus experiment

The Phebus experiment on board the GRANAT satellite provides temporal and spectral observations of solar and cosmic -ray bursts in the 0.1 100 MeV nominal energy range. The experiment was turned on January 8, 1990 and is still in operation. In this paper we present the main characteristics of the Phebus experiment and we describe and discuss some of the observational properties of the 18 solar hard X-ray/-ray events detected during the first semester of the Phebus operation. It is found that: (i) events of a few minutes duration, detected above 100 keV, systematically show subsecond time variations; (ii) longer duration events (>5 min) do not exhibit fast time variations and generally consist of 1-min peaks superimposed on a less intense, sometimes harder, slowly varying component. In addition to these general trends we discuss in more detail three events for which significant count-rates have been detected above 10 MeV.     

Robe’s restricted problem of 2+2 bodies with one of the primaries an oblate body

In this problem, one of the primaries of mass \(m^{*}_{1}\) is a rigid spherical shell filled with a homogeneous incompressible fluid of density ρ ₁. The smaller primary of mass m ₂ is an oblate body outside the shell. The third and the fourth bodies (of mass m ₃ and m ₄ respectively) are small solid spheres of density ρ ₃ and ρ ₄ respectively inside the shell, with the assumption that the mass and the radius of the third and the fourth body are infinitesimal. We assume that m ₂ is describing a circle around \(m^{*}_{1}\) . The masses m ₃ and m ₄ mutually attract each other, do not influence the motions of \(m^{*}_{1}\) and m ₂ but are influenced by them. We also assume that masses m ₃ and m ₄ are moving in the plane of motion of mass m ₂. In the paper, equilibrium solutions of m ₃ and m ₄ and their linear stability are analyzed. There are two collinear equilibrium solutions for the system. The non collinear equilibrium solutions exist only when ρ ₃=ρ ₄. There exist an infinite number of non collinear equilibrium solutions of the system, provided they lie inside the spherical shell. In a system where the primaries are considered as earth-moon and m ₃,m ₄ as submarines, the collinear equilibrium solutions thus obtained are unstable for the mass parameters μ,μ ₃,μ ₄ and oblateness factor A. In this particular case there are no non-collinear equilibrium solutions of the system.     

Astrometric Results for Observations of Jupiter’s Galilean Satellites During Mutual Occultations and Eclipses in 2009 and 2014–2015

Solar System Research - Photometric observations of satellites during their mutual occultations and eclipses are a valuable source of astrometric data for studying the motion of natural planetary...     

A photometric and spectroscopic study of WW And – an Algol-type,long period binary system with an accretion disc

We have analysed the available spectra of WW And and for the first time obtained a reasonably well defined radial velocity curve of the primary star. Combined with the available radial velocity curve of the secondary component, these data led to the first determination of the spectroscopic mass ratio of the system at q_spec = 0.16 ± 0.03. We also determined the radius of the accretion disc from analysis of the double-peaked H_α emission lines. Our new, high-precision, Johnson VRI and the previously-available Strömgren vby light curves were modelled with stellar and accretion disc models. A consistent model for WW And – a semidetached system harbouring an accretion disc which is optically thick in its inner region, but optically thin in the outer parts – agrees well with both spectroscopic and photometric data.     

Highly accurate determination of the coordinates and the Earth’s rotation parameters involving the Svetloe VLBI Observatory

We present the results of our processing of the first observations of extragalactic radio sources obtained with the eight-element International VLBI Network, which includes the Svetloe Russian Radio Astronomy Observatory equipped with a Mark 3A recording terminal. Our observations and their processing yielded highly accurate coordinates (in meters) of the Svetloe Observatory in the ITRF 2000 system: X = 2730173.854 ± 0.002, Y = 1562442.668 ± 0.004, Z = 5529969.069 ± 0.007. We also show that including the Svetloe Observatory in the International Network led to an appreciable improvement in the accuracy of determining the Earth’s rotation parameters (microarcseconds for the coordinates of the pole and nutation angles, microseconds for Universal Time): X_p = ?154683 ± 77, Y_p = 361809 ± 59, UT1-UTC = ?325162.9 ± 2.5, Δψ = ?53147 ± 114, Δε = ?2286 ± 47.     

Changes in the Sun’s mass and gravitational constant estimated using modern observations of planets and spacecraft

More than 635 thousand positional observations of planets and spacecraft of various types (mostly radiotechnical ones, 1961–2010) were used to estimate possible changes in the gravitational constant, Sun’s mass, and semi-major axes of planetary orbits, as well as the associated value of the astronomical unit. The observations were analyzed based on the EPM2010 ephemerides constructed at the Institute of Applied Astronomy (Russian Academy of Sciences) in a post-Newtonian approximation as a result of simultanious numerical integration of the equations of motion of nine major planets, the Sun, the Moon, asteroids, and trans-Neptunian objects. The heliocentric gravitational constant GM _⊙ was found to vary with a rate of (GṀ _⊙/GM _⊙ = (−5.0 ± 4.1)) × 10⁻¹⁴ per year (at the 3σ level). The positive secular changes in the semimajor axes ȧ _i /a _i were found for Mercury, Venus, Mars, Jupiter, and Saturn provided by high-precision observations. These changes also correspond to the decrease in the heliocentric gravitational constant. The changing of GM _⊙, itself is probably caused by the loss of the mass M _⊙ of the Sun due to its radiation and solar wind; these effects are partly compensated by the material falling onto the Sun. Allowing for the maximum bounds on the possible change in the Sun’s mass M _⊙, it has been found from the change obtained in GM _⊙ that the annual change Ġ/G of the gravitational constant G falls within the interval −4.2 × 10⁻¹⁴ < ȧ/G < +7.5 × 10⁻¹⁴ with a 95% probability. The astronomical unit (AU) is connected by its definition only with the heliocentric gravitational constant. The decrease of GM _⊙ obtained in this paper should correspond to a secular decrease in the AU. It is shown, however, that the modern level of accuracy does not allow us to determine a change in the AU. The attained posibility of determining changes in GM _⊙ using high-accuracy observations encourages us to have a relation between GM _⊙ and the AU fixed for a certain moment in time, since it is inconvenient to have a time-dependent length for the AU.     

Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models

The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.     

Temporal Pointing Variations of the Solar Dynamics Observatory’s HMI and AIA Instruments on Subweekly Time Scales

Achieving subarcsecond co-registration across varying time-lines of multi-wavelength and instrument images is difficult and requires an accurate characterization of the instrument pointing jitter. We investigated the internal pointing errors on daily and yearly time-scales that occur across the Solar Dynamics Observatory’s (SDO) Atmospheric Imaging Assembly (AIA) and Helioseismic Magnetic Imager (HMI). Using cross-correlation techniques on the AIA 1700 Å passband and the HMI line-of-sight magnetograms from three years of observational image pairs at approximately three-day intervals, internal pointing errors were quantified. Pointing variations of ±?0.26″ (jitter-limited) and ±?0.50″ in the solar East–West (x) and North–South (y) directions, respectively, were measured. AIA observations of the Venus transit in June 2012 were used to measure existing coalignment offsets in all passbands. We found that the AIA passband pointing variations are 〈ΔX _CO〉=1.10″±1.41″ and 〈ΔY _CO〉=1.25″±1.24″ when aligned to the HMI nominal image center, referred to here as the CutOut technique. Minimal long-term pointing variations found between limb and correlation derived pointings provide evidence that the image-center positions provided by the instrument teams achieve single-pixel accuracy on time scales shorter than their characterization. However, daily AIA passband pointing variations of ??1.18″ indicate that autonomous subarcsecond co-registration is not fully achieved yet.     

An analysis of apparent r-mode oscillations in solar activity,the solar diameter,the solar neutrino flux,and nuclear decay rates,with implications concerning the Sun’s internal structure and rotation,and neutrino processes

This article presents a comparative analysis of solar activity data, Mt Wilson diameter data, Super-Kamiokande solar neutrino data, and nuclear decay data acquired at the Lomonosov Moscow State University (LMSU). We propose that salient periodicities in all of these datasets may be attributed to r-mode oscillations. Periodicities in the solar activity data and in Super-Kamiokande solar neutrino data may be attributed to r-mode oscillations in the known tachocline, with normalized radius in the range 0.66–0.74, where the sidereal rotation rate is in the range 13.7–14.6 year⁻¹. We propose that periodicities in the Mt Wilson and LMSU data may be attributed to similar r-mode oscillations where the sidereal rotation rate is approximately 12.0 year⁻¹, which we attribute to a hypothetical “inner” tachocline separating a slowly rotating core from the radiative zone. We also discuss the possible role of the Resonant Spin Flavor Precession (RSFP) process, which leads to estimates of the neutrino magnetic moment and of the magnetic field strength in or near the solar core.     

Families of periodic orbits in Hill’s problem with solar radiation pressure: application to Hayabusa 2

This work studies periodic solutions applicable, as an extended phase, to the JAXA asteroid rendezvous mission Hayabusa 2 when it is close to target asteroid 1999 JU3. The motion of a spacecraft close to a small asteroid can be approximated with the equations of Hill’s problem modified to account for the strong solar radiation pressure. The identification of families of periodic solutions in such systems is just starting and the field is largely unexplored. We find several periodic orbits using a grid search, then apply numerical continuation and bifurcation theory to a subset of these to explore the changes in the orbit families when the orbital energy is varied. This analysis gives information on their stability and bifurcations. We then compare the various families on the basis of the restrictions and requirements of the specific mission considered, such as the pointing of the solar panels and instruments. We also use information about their resilience against parameter errors and their ground tracks to identify one particularly promising type of solution.     

Collinear equilibrium points of Hill’s problem with radiation and oblateness and their fractal basins of attraction

We discuss the existence, location, and stability of the collinear equilibrium points of a generalized Hill problem with radiation of the primary (the Sun) and oblateness of the secondary (the planet), and present some remarkable fractals created as basins of attraction of Newton’s method applied for their computation in several cases of the parameters.     

Bicentennial decrease of the solar constant leads to the Earth’s unbalanced heat budget and deep climate cooling

Long-wave energy emitted by the Earth-atmosphere into space is characterized by changes in power over time that always lag behind the changes in power of the absorbed solar radiation due to slow variation in enthalpy of the Earth-atmosphere system. Long-term variation of the solar energy radiation absorbed by the Earth remains uncompensated by the energy radiated into space over the interval of time that is determined by the thermal inertia. The basic state of the climate system is when the debit and credit sides in the Earth’s global annual mean energy budget (including the air and water envelopes) are almost always unbalanced. The annual mean balance of the heat budget of the Earth-atmosphere over a long time period will reliably define the behavior and magnitude of the energy excess accumulated by the Earth or energy deficit to allow us to determine adequately and to predict beforehand the trend and amplitude of the forthcoming climate change using the prognosis of variations in the total solar irradiance (solar constant). The decrease in solar constant has been observed since the early 1990s. The Earth as a planet will have a negative balance in the energy budget in the future as well, because the Sun is entering the decline phase of the bicentennial luminosity changes. This will lead to a drop in temperature in approximately 2014. The increase in albedo and decrease in greenhouse gas concentration in the atmosphere will result in the additional decrease in absorbed portion of the solar energy and reduced greenhouse effect. The additional drop in temperature exceeding the effect of decreased solar constant can occur as a result of successive feedback effects. A deep bicentennial minimum in solar constant is to be anticipated in 2042 ± 11 and the 19th Little Ice Age (for the last 7500 years) may occur in 2055 ± 11.     

Diurnal variation of electron density in Saturn’s ionosphere: Model comparisons with Saturn Electrostatic Discharge (SED) observations

Using the Saturn Thermosphere Ionosphere Model (STIM), we present a study of the diurnal variation of electron density, with a focus on comparisons with peak electron densities (N_MAX) inferred from the low-frequency cutoff of radio emission due to lightning in the lower atmosphere, called Saturn Electrostatic Discharges (SEDs). It is demonstrated that photochemistry in Saturn’s ionosphere cannot reproduce the SED-inferred diurnal variation in N_MAX unless additional production and loss sources outside of the current best estimates are considered. Additional explanations of the SED-inferred diurnal variation of N_MAX are presented and analyzed, such as the possibility that the low-frequency cutoff seen in SEDs is due to the presence of sharp low-altitude layers of plasma, as frequently seen in radio occultation measurements. Finally, we outline the observational constraints that must be fulfilled by any candidate explanations of the SED-inferred diurnal variation of N_MAX.     

Cross-correlation of solar wind parameters with sunspots (‘Long-term variations’) at 1 AU during cycles 21 and 22

Long-term variations of solar wind parameters at 1 AU are correlated with sunspots for the time interval 1973 to 1993 (solar cycles 21, 22). Using theNear-Earth Heliosphere Data OMNI the plasma density, the magnitude of the interplanetary magnetic field, the solar wind velocity and the solar wind temperature show consistent long-term variations in each cycle (21 and 22) — pointing to specifictime-lags in the coupling between sunspots (and the underlying convection zone), the solar corona and the solar wind parameters at 1 AU (ecliptic).     

Comparison of irregular variations in the Earth’s rotation rate and parameters of geophysical and atmospheric processes with events on the sun

A study of the Earth’s rotation in space reveals a complex pattern of variations in its orientation, the excitation mechanisms of these variations, and their manifestations in various natural processes. The Earth’s rotation rate is not constant but exhibits complex fluctuations that account for some fraction of 10⁸ (corresponding to variations of several milliseconds (ms) in the length of the day). These variations span a wide spectrum of time scales, from hours to centuries or longer, reflecting the fact that they are produced by a wide variety of geophysical and astronomical processes. We discuss the results of our statistical comparison of long series of observations to reveal the most coherent variations. The spectral composition of the experimental time series has been determined using modified periodogram and single-channel autoregression methods. A comparative analysis has been performed by a two-channel autoregression spectral estimation method. The results of our comparison of the time series suggest that the fluctuations with periods of about 73 years are highly coherent.     

Robe’s restricted problem of 2+2 bodies when the bigger primary is a Roche ellipsoid and the smaller primary is an oblate body

In this problem, one of the primaries of mass m ₁ is a Roche ellipsoid filled with a homogeneous incompressible fluid of density ρ ₁. The smaller primary of mass m ₂ is an oblate body outside the Ellipsoid. The third and the fourth bodies (of mass m ₃ and m ₄ respectively) are small solid spheres of density ρ ₃ and ρ ₄ respectively inside the Ellipsoid, with the assumption that the mass and the radius of the third and the fourth body are infinitesimal. We assume that m ₂ is describing a circle around m ₁. The masses m ₃ and m ₄ mutually attract each other, do not influence the motions of m ₁ and m ₂ but are influenced by them. We have extended the Robe’s restricted three-body problem to 2+2 body problem under the assumption that the fluid body assumes the shape of the Roche ellipsoid (Chandrashekhar in Ellipsoidal figures of equilibrium, Chap. 8, Dover, New York, 1987). We have taken into consideration all the three components of the pressure field in deriving the expression for the buoyancy force viz (i) due to the own gravitational field of the fluid (ii) that originating in the attraction of m ₂ (iii) that arising from the centrifugal force. In this paper, equilibrium solutions of m ₃ and m ₄ and their linear stability are analyzed. We have proved that there exist only six equilibrium solutions of the system, provided they lie within the Roche ellipsoid. In a system where the primaries are considered as Earth-Moon and m ₃,m ₄ as submarines, the equilibrium solutions of m ₃ and m ₄ respectively when the displacement is given in the direction of x ₁-axis or x ₂-axis are unstable.