An estimation of Envisat’s rotational state accounting for the precession of its rotational axis caused by gravity-gradient torque

The rotational state of Envisat is re-estimated using the specular glint times in optical observation data obtained from 2013 to 2015. The model is simplified to a uniaxial symmetric model with the first order variation of its angular momentum subject to a gravity-gradient torque causing precession around the normal of the orbital plane. The sense of Envisat’s rotation can be derived from observational data, and is found to be opposite to the sense of its orbital motion. The rotational period is estimated to be $(120.674\pm 0.068)·\exp \left((4.5095\pm 0.0096)\times {10}^{-4}·t\right)\text{s}$, where t is measured in days from the beginning of 2013. The standard deviation is 0.760?s, making this the best fit obtained for Envisat in the literature to date. The results demonstrate that the angle between the angular momentum vector and the negative normal of the orbital plane librates around a mean value of $8.53°\pm 0.42°$ with an amplitude from about $0.7°$ (in 2013) to $0.5°$ (in 2015), with the libration period equal to the precession period of the angular momentum, from about 4.8?days (in 2013) to 3.4?days (in 2015). The ratio of the minimum to maximum principal moments of inertia is estimated to be $0.0818\pm 0.0011$, and the initial longitude of the angular momentum in the orbital coordinate system is $40.5°\pm 9.3°$. The direction of the rotation axis derived from our results at September 23, 2013, UTC 20:57 is similar to the results obtained from satellite laser ranging data but about $20°$ closer to the negative normal of the orbital plane.     

Vertical and radial metallicity gradients in high latitude galactic fields with SDSS

We used the ugr magnitudes of 1437467 F-G type main-sequence stars with metal abundance $-2?[\mathrm{Fe}/\mathrm{H}]?+0.2$ dex and estimated radial and vertical metallicity gradients for high Galactic-latitude fields, $50°<b?90°$ and $0°<l?360°$, of the Milky Way Galaxy. The radial metallicity gradient $d[\mathrm{Fe}/\mathrm{H}]/\mathit{dR}=-0.042\pm 0.011$ dex kpc^?1 estimated for the stars with $1.31<z\le 1.74$ kpc is attributed to the thin-disc population. While, the radial gradients evaluated for stars at higher vertical distances are close to zero indicating that the thick disc and halo have not undergone a radial collapse phase at least at high Galactic latitudes. The vertical metallicity gradients estimated for stars with three different Galactic latitudes, $50°<b?65°,65°<b?80°$ and $80°<b?90°$ do not show a strong indication for Galactic latitude dependence of our gradients. The thin disc, $0.5<z?2$ kpc, with a vertical metallicity gradient $d\left[\mathrm{Fe}/\mathrm{H}\right]/\mathit{dz}=-0.308\pm 0.018$ dex kpc^?1, is dominant only in galactocentric distance interval $6<R?10$ kpc, while the thick disc ($2<z?5$ kpc) could be observed in the intervals $6<R?10$ and $10<R?15$ kpc with compatible vertical metallicity gradients, i.e. $d\left[\mathrm{Fe}/\mathrm{H}\right]/\mathit{dz}=-0.164\pm 0.014$ dex kpc^?1 and $d\left[\mathrm{Fe}/\mathrm{H}\right]/\mathit{dz}=-0.172\pm 0.016$ dex kpc^?1. Five vertical metallicity gradients are estimated for the halo ($z>5$ kpc) in three galactocentric distance intervals, $6<R?10,10<R?15$ and $15<R?20$ kpc. The first one corresponding to the interval $6<R?10$ kpc is equal to $d\left[\mathrm{Fe}/\mathrm{H}\right]/\mathit{dz}=-0.023\pm 0.006$ dex kpc^?1, while the others at larger galactocentric distances are close to zero. We derived synthetic vertical metallicity gradients for 2,230,167 stars and compared them with the observed ones. There is a good agreement between the two sets of vertical metallicity gradients for the thin disc, while they are different for the thick disc. For the halo, the conspicuous difference corresponds to the galactocentric distance interval $6<R?10$ kpc, while they are compatible at higher galactocentric distance intervals.     

Near-infrared study of open clusters Teutsch 10 and Teutsch 25

The astrophysical parameters have been estimated for two unstudied open star clusters Teutsch 10 and Teutsch 25 using the Two Micron All Sky Survey ($2\mathit{MASS}$) database. Radius is estimated as 4.5 arcmin for both clusters using radial density profiles. We have estimated proper motion values in both RA and DEC directions as $2.28\pm 0.3$ and $-0.38\pm 0.11$?mas?yr^?1 for Teutsch 10 and $0.48\pm 0.3$ and $3.35\pm 0.16$?mas?yr^?1 for Teutsch 25 using $\mathit{PPMXL}$¹ catalog. By estimating the stellar membership probabilities, we have identified 30 and 28 most likely members for Teutsch 10 and Teutsch 25 respectively. We have estimated the reddening as $E(B-V)=0.96\pm 0.3$?mag for Teutsch 10 and $0.58\pm 0.2$?mag for Teutsch 25, while the corresponding distances are $2.4\pm 0.2$ and $1.9\pm 0.1$?kpc. Ages of $70\pm 10$?Myr for Teutsch 10 and $900\pm 100$?Myr for Teutsch 25 are estimated using the theoretical isochrones of metallicity Z?=?0.019. The mass function slopes are derived as $1.23\pm 0.30$ and $1.09\pm 0.35$ for Teutsch 10 and Teutsch 25 respectively. Estimated mass function slope for both the clusters are close to the Salpeter value ($x=1.35$) within the errors. Estimated values of dynamical relaxation time are found to be less than cluster’s age for these objects. This concludes that both objects are dynamically relaxed. The possible reason for relaxation may be due to dynamical evolution or imprint of star formation or both.     

Physical characterization of the deep-space debris WT1190F: A testbed for advanced SSA techniques

We report on extensive ${\mathit{BVR}}_c{I}_c$ photometry and low-resolution ($\lambda /\mathrm{\Delta }\lambda ～250$) spectroscopy of the deep-space debris WT1190F, which impacted Earth offshore from Sri Lanka, on 2015 November 13. In spite of its likely artificial origin (as a relic of some past lunar mission), the case offered important points of discussion for its suggestive connection with the envisaged scenario for a (potentially far more dangerous) natural impactor, like an asteroid or a comet.Our observations indicate for WT1190F an absolute magnitude $R_c={32.45}_{\pm 0.31}$, with a flat dependence of reflectance on the phase angle, such as ${\mathit{dR}}_c/d?～{0.007}_{\pm 2}$?mag?deg^?1. The detected short-timescale variability suggests that the body was likely spinning with a period twice the nominal figure of $P_{\mathrm{flash}}={1.4547}_{\pm 0.0005}\mathrm{s}$, as from the observed lightcurve. In the ${\mathit{BVR}}_c{I}_c$ color domain, WT1190F closely resembled the Planck deep-space probe. This match, together with a depressed reflectance around 4000 and 8500 Å may be suggestive of a “grey” (aluminized) surface texture.The spinning pattern remained in place also along the object fiery entry in the atmosphere, a feature that may have partly shielded the body along its fireball phase perhaps leading a large fraction of its mass to survive intact, now lying underwater along a tight ($～1\times 80$?km) strip of sea, at a depth of 1500?m or less.Under the assumption of Lambertian scatter, an inferred size of ${216}_{\pm 30}/\sqrt{\alpha /0.1}$?cm is obtained for WT1190F. By accounting for non-gravitational dynamical perturbations, the Area-to-Mass ratio of the body was in the range $(0.006?\mathit{AMR}?0.011)$?m²?kg^?1.Both these figures resulted compatible with the two prevailing candidates to WT1190F’s identity, namely the Athena II Trans-Lunar Injection Stage of the Lunar Prospector mission, and the ascent stage of the Apollo 10 lunar module, callsign “Snoopy”. Both candidates have been analyzed in some detail here through accurate 3D CAD design mockup modelling and BRDF reflectance rendering to derive the inherent photometric properties to be compared with the observations.     

Medium Earth Orbit dynamical survey and its use in passive debris removal

The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable $\mathrm{\Delta }V$ cost. We simulated the dynamics over 120–200?years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90°, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed $～40$ years, mainly around three specific inclination values: 46°, 56°, and 68°; a result compatible with our understanding of MEO secular dynamics. For $\mathrm{\Delta }V?300$ m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with $e<0.3$), reentry times from GNSS altitudes exceed $～70$ years, while low-cost ($\mathrm{\Delta }V?5–35$ m/s) graveyard orbits, stable for at lest 200?years, are found for eccentricities up to $e\approx 0.018$. This investigation was carried out in the framework of the EC-funded “ReDSHIFT” project.     

Modelling ionospheric vertical drifts over Africa low latitudes using Empirical Orthogonal functions and comparison with climatological model

For the first time, empirical model of daytime vertical $\mathbf{E}\times \mathbf{B}$ drift based on Empirical Orthogonal functions (EOF) decomposition technique is presented. Day-to-day variability of $\mathbf{E}\times \mathbf{B}$ drift inferred from horizontal (H) geomagnetic field data around dip latitude for the period of 2008–2013 is used to both develop and validate the model. Results show that the EOF technique is promising with modelled values and data giving correlation coefficient values of at least 0.90 for geomagnetic conditions of both $K_p?3$ and $K_p>3$ within 2008–2013. Independent model validation shows that in situ $\mathbf{E}\times \mathbf{B}$ values from ion velocity meter (IVM) instrument on-board C/NOFS satellite are closer to model $\mathbf{E}\times \mathbf{B}$ estimates than the climatological Scherliess-Fejer (SF) model incorporated within the International Reference Ionosphere (IRI).     

Velocity shear Kelvin-Helmholtz instability with inhomogeneous DC electric field in the magnetosphere of Saturn

In this paper parallel flow velocity shear Kelvin-Helmholtz instability has been studied in two different extended regions of the inner magnetosphere of Saturn. The method of the characteristic solution and kinetic approach has been used in the mathematical calculation of dispersion relation and growth rate of K-H waves. Effect of magnetic field (B), inhomogeneity (P/a), velocity shear scale length ($A_i$), temperature anisotropy ($T_{\perp }/T_{||}$), electric field (E), ratio of electron to ion temperature ($T_e/T_i$), density gradient (${\epsilon }_n{\rho }_i$) and angle of propagation ($\theta$) on the dimensionless growth rate of K-H waves in the inner magnetosphere of Saturn has been observed with respect to $k_{\perp }{\rho }_i$. Calculations of this theoretical analysis have been done taking the data from the Cassini in the inner magnetosphere of Saturn in the two extended regions of Rs ～4.60–4.01 and Rs ～4.82–5.0. In our study velocity shear, temperature anisotropy and magnitude of the electric field are observed to be the major sources of free energy for the K-H instability in both the regions considered. The inhomogeneity of electric field, electron-ion temperature ratio, and density gradient have been observed playing stabilizing effect on K-H instability. This study also indicates the effect of the vicinity of icy moon Enceladus on the growth of K-H instability.     

Impact of rogue active regions on hemispheric asymmetry

The solar dipole moment at activity minimum is a good predictor of the strength of the subsequent solar cycle. Through a systematic analysis using a state-of-the-art $2\times 2$D solar dynamo model, we found that bipolar magnetic regions (BMR) with atypical characteristics can modify the strength of the next cycle via their impact on the buildup of the dipole moment as a sunspot cycle unfolds. In addition to summarizing these results, we present further effects of such “rogue” BMRs. These have the ability to generate hemispheric asymmetry in the subsequent sunspot cycle, since they modify the polar cap flux asymmetry of the ongoing cycle. We found strong correlation between the polar cap flux asymmetry of cycle i and the total pseudo sunspot number asymmetry of cycle $i+1$. Good correlation also appears in the case of the time lag of the hemispheres of cycle $i+1$.     

Kelvin–Helmholtz instability in a twisting solar polar coronal hole jet observed by SDO/AIA

We investigate the conditions under which the fluting ($m=2$), $m=3$, and $m=12$ magnetohydrodynamic (MHD) modes in a uniformly twisted flux tube moving along its axis become unstable in order to model the Kelvin–Helmholtz (KH) instability in a twisting solar coronal hole jet near the northern pole of the Sun. We employed the dispersion relations of MHD modes derived from the linearized MHD equations. We assumed real wavenumbers and complex angular wave frequencies, namely complex wave phse velocities. The dispersion relations were solved numerically at fixed input parameters (taken from observational data) and varying degrees of torsion of the internal magnetic field. It is shown that the stability of the modes depends upon five parameters: the density contrast between the flux tube and its environment, the ratio of the external and internal axial magnetic fields, the twist of the magnetic field lines inside the tube, the ratio of transverse and axial jet’s velocities, and the value of the Alfvén Mach number (the ratio of the tube axial velocity to Alfvén speed inside the flux tube). Using a twisting jet of 2010 August 21 by SDO/AIA and other observations of coronal jets we set the parameters of our theoretical model and have obtained that in a twisted magnetic flux tube of radius of $9.8$ Mm, at a density contrast of $0.474$ and fixed Alfvén Mach number of $?0.76$, for the three MHD modes there exist instability windows whose width crucially depends upon the internal magnetic field twist. It is found that for the considered modes an azimuthal magnetic field of $1.3$–$1.4$ G (computed at the tube boundary) makes the width of the instability windows equal to zero, that is, it suppress the KH instability onset. On the other hand, the times for developing KH instability of the $m=12$ MHD mode at instability wavelengths between 15 and 12 Mm turn out to be in the range of $1.9$–$4.7$ min that is in agreement with the growth rates estimated from the temporal evolution of the observed unstable jet’s blobs in their initial stage.     

A combination method using evolutionary algorithms in initial orbit determination for too short arc

With the combination of two evolutionary algorithms EDA and DE, a new method of initial orbit determination for satellites based on ground-based too-short-arc is established. Compared with other algorithms, the proposed method focuses on the most densely populated region in the solution space rather than the individual with best fitness value. Both the global information and local information are well fused in the search of optimum. In the method $(a,e,M)$ are treated as variables of the optimization, and the optimization procedure is carried out as a two-stage hierarchical optimization problem which has three variables for each stage. Kernel density estimation is applied to build the probability distribution model without any assumptions of the specified distribution, accompanied by handling semi-major axis and eccentricity as a pair of dependent variables in the construction of the probability for the correlation between them in the practice. Numerical experiments with real ground-based observations show that the proposed method is applicable to too-short-arc with even 3?s, and the result of bias in several kilometers can be achieved with $5^{″}$ error added to angular measurements.     

First solar observations with ALMA

The Atacama Large Millimeter-Submillimeter Array (ALMA) has opened a new window for studying the Sun via high-resolution high-sensitivity imaging at millimeter wavelengths. In this contribution I review the capabilities of the instrument for solar observing and describe the extensive effort taken to bring the possibility of solar observing with ALMA to the scientific community. The first solar ALMA observations were carried out during 2014 and 2015 in two ALMA bands, Band 3 ($\lambda =3$?mm) and Band 6 ($\lambda =1.3$?mm), in single-dish and interferometric modes, using single pointing and mosaicing observing techniques, with spatial resolution up to $～$2″ and $～$1″ in the two bands, respectively. I overview several recently published studies which made use of the first solar ALMA observations, describe current status of solar observing with ALMA and briefly discuss the future capabilities of the instrument.     

Effects of space weather on the ionosphere and LEO satellites’ orbital trajectory in equatorial,low and middle latitude

We study the effects of space weather on the ionosphere and low Earth orbit (LEO) satellites’ orbital trajectory in equatorial, low- and mid-latitude (EQL, LLT and MLT) regions during (and around) the notable storms of October/November, 2003. We briefly review space weather effects on the thermosphere and ionosphere to demonstrate that such effects are also latitude-dependent and well established. Following the review we simulate the trend in variation of satellite’s orbital radius (r), mean height (h) and orbit decay rate (ODR) during 15 October–14 November 2003 in EQL, LLT and MLT. Nominal atmospheric drag on LEO satellite is usually enhanced by space weather or solar-induced variations in thermospheric temperature and density profile. To separate nominal orbit decay from solar-induced accelerated orbit decay, we compute $r\text{,}h$ and ODR in three regimes viz. (i) excluding solar indices (or effect), where $r=r_0\text{,}h=h_0$ and $\mathit{ODR}={\mathit{ODR}}_0$ (ii) with mean value of solar indices for the interval, where $r=r_m\text{,}h=h_m$ and $\mathit{ODR}={\mathit{ODR}}_m$ and (iii) with actual daily values of solar indices for the interval ($r\text{,}h$ and ODR). For a typical LEO satellite at h?=?450?km, we show that the total decay in r during the period is about 4.20?km, 3.90?km and 3.20?km in EQL, LLT and MLT respectively; the respective nominal decay ($r_0$) is 0.40?km, 0.34?km and 0.22?km, while solar-induced orbital decay ($r_m$) is about 3.80?km, 3.55?km and 2.95?km. h also varied in like manner. The respective nominal ${\mathit{ODR}}_0$ is about 13.5?m/day, 11.2?m/day and 7.2?m/day, while solar-induced ${\mathit{ODR}}_m$ is about 124.3?m/day, 116.9?m/day and 97.3?m/day. We also show that severe geomagnetic storms can increase ODR by up to 117% (from daily mean value). However, the extent of space weather effects on LEO Satellite’s trajectory significantly depends on the ballistic co-efficient and orbit of the satellite, and phase of solar cycles, intensity and duration of driving (or influencing) solar event.     

Simulation-based evaluation of a cold atom interferometry gradiometer concept for gravity field recovery

The prospects of future satellite gravimetry missions to sustain a continuous and improved observation of the gravitational field have stimulated studies of new concepts of space inertial sensors with potentially improved precision and stability. This is in particular the case for cold-atom interferometry (CAI) gradiometry which is the object of this paper. The performance of a specific CAI gradiometer design is studied here in terms of quality of the recovered gravity field through a closed-loop numerical simulation of the measurement and processing workflow. First we show that mapping the time-variable field on a monthly basis would require a noise level below $5\text{mE}/\sqrt{\text{Hz}}$. The mission scenarios are therefore focused on the static field, like GOCE. Second, the stringent requirement on the angular velocity of a one-arm gradiometer, which must not exceed ${10}^{-6}$?rad/s, leads to two possible modes of operation of the CAI gradiometer: the nadir and the quasi-inertial mode. In the nadir mode, which corresponds to the usual Earth-pointing satellite attitude, only the gradient $V_{\mathit{yy}}$, along the cross-track direction, is measured. In the quasi-inertial mode, the satellite attitude is approximately constant in the inertial reference frame and the 3 diagonal gradients $V_{\mathit{xx}}\text{,}{V}_{\mathit{yy}}$ and $V_{\mathit{zz}}$ are measured. Both modes are successively simulated for a 239?km altitude orbit and the error on the recovered gravity models eventually compared to GOCE solutions. We conclude that for the specific CAI gradiometer design assumed in this paper, only the quasi-inertial mode scenario would be able to significantly outperform GOCE results at the cost of technically challenging requirements on the orbit and attitude control.     

Constraining the mass of the black hole GX 339-4 using spectro-temporal analysis of multiple outbursts

We carried out spectro-temporal analysis of the archived data from multiple outbursts spanning over the last two decades from the black hole X-ray binary GX 339-4. In this paper, the mass of the compact object in the X-ray binary system GX 339-4 is constrained based on three indirect methods. The first method uses broadband spectral modelling with a two component flow structure of the accretion around the black hole. The broadband data are obtained from RXTE (Rossi X-ray Timing Explorer) in the range 3.0 to 150.0?keV and from Swift and NuSTAR (Nuclear Spectroscopic Telescope Array) simultaneously in the range 0.5 to 79.0?keV. In the second method, we model the time evolution of Quasi-periodic Oscillation (QPO) frequencies, considering it to be the result of an oscillating shock that radially propagates towards or away from the compact object. The third method is based on scaling a mass dependent parameter from an empirical model of the photon index ($\mathrm{\Gamma }$) – QPO ($\nu$) correlation. We compare the results at 90 percent confidence from the three methods and summarize the mass estimate of the central object to be in the range $8.28–11.89M_{\odot }$.     

Highlights about the performances of storm-time TEC modelling techniques for low/equatorial and mid-latitude locations

A statistical evaluation of storm-time total electron content (TEC) modelling techniques over various latitudes of the African sector and surrounding areas is presented. The source of observational TEC data used in this study is the Global Navigation Satellite Systems (GNSS), specifically the Global Positioning Systems (GPS) receiver networks. For each selected receiver station, three different storm-time models based on empirical orthogonal functions (EOF) analysis, non-linear regression analysis (NLRA) and Artificial neural networks (ANN), were implemented. Storm-time GPS TEC data used for both development and validation of the models was selected based on the storm criterion of $\mathit{Dst}?-50$ nT or $K_p?4$ to take into account both coronal mass ejections (CMEs) and co-rotating interaction regions (CIRs) driven storms, respectively. To make an independent test of the models, storm periods considered for validation were excluded from datasets used during the implementation of the models and results are compared with observations, monthly median values, and International Reference Ionosphere (IRI-2016) predictions. Considering GPS TEC as reference, a statistical analysis performed over six storm periods reserved for validation revealed that ANN model is about 10%, 26%, and 58% more accurate than EOF, NLRA, and IRI models, respectively. It was further found that, EOF model performs 15%, and 44% better than NLRA, and IRI models, respectively, while NLRA is 25% better than IRI. On the other hand, results are also discussed referring to the background ionosphere represented by monthly median TEC (MM TEC) and statistics are provided. Moreover, strengths and weaknesses of each model are highlighted.