Spacecraft survivability in the natural debris environment near the stable Earth-Moon Lagrange points

The theoretical analysis of the motion of natural space debris near the stable Earth-Moon Lagrange Points, $L_4$ and $L_5$, is presented with a focus on the potential debris risks to spacecraft operating near these points. Specifically, the research formulates a debris propagation model using four-body dynamics, then applies candidate probabilistic survivability models to a notional spacecraft operating at the $L_4$ and $L_5$ Lagrange points to quantify the collision risks to the spacecraft from natural debris particles. Of the survivability models implemented, the natural debris collision risks to spacecraft survivability are found to be incredibly low, but mitigation strategies to reduce the risk further are identified in this study. Overall, research into stable Lagrange point natural debris propagation improves understanding of the collision risks posed by the naturally occurring Kordylewski clouds and enhances operational planning for Lagrange point space missions.     

Where is the base of the Transition Region? Evidence from TRACE,SDO, IRIS and ALMA observations

Classic solar atmospheric models put the Chromosphere-Corona Transition Region (CCTR) at $～2$ Mm above the ${\tau }_{5000}=1$ level, whereas radiative MHD (rMHD) models place the CCTR in a wider range of heights. However, observational verification is scarce. In this work we review and discuss recent results from various instruments and spectral domains. In SDO and TRACE images spicules appear in emission in the 1600, 1700 and 304 Å bands and in absorption in the EUV bands; the latter is due to photo-ionization of H i and He i, which increases with wavelength. At the shortest available AIA wavelength and taking into account that the photospheric limb is $～0.34$ Mm above the ${\tau }_{5000}=1$ level, we found that CCTR emission starts at $～3.7$ Mm; extrapolating to $\lambda =0$, where there is no chromospheric absorption, we deduced a height of $3.0\pm 0.5$ Mm, which is above the value of 2.14 Mm of the Avrett and Loeser model. Another indicator of the extent of the chromosphere is the height of the network structures. Height differences produce a limbward shift of features with respect to the position of their counterparts in magnetograms. Using this approach, we measured heights of $0.14\pm 0.04$ Mm (at 1700 Å), $0.31\pm 0.09$ Mm (at 1600 Å) and $3.31\pm 0.18$ Mm (at 304 Å) for the center of the solar disk. A previously reported possible solar cycle variation is not confirmed. A third indicator is the position of the limb in the UV, where IRIS observations of the Mg ii triplet lines show that they extend up to $～2.1$ Mm above the 2832 Å limb, while AIA/SDO images give a limb height of $1.4\pm 0.2$ Mm (1600 Å) and $5.7\pm 0.2$ Mm (304 Å). Finally, ALMA mm-$\lambda$ full-disk images provide useful diagnostics, though not very accurate, due to their relatively low resolution; values of $2.4\pm 0.7$ Mm at 1.26 mm and $4.2\pm 2.5$ Mm at 3 mm were obtained. Putting everything together, we conclude that the average chromosphere extends higher than homogeneous models predict, but within the range of rMHD models..     

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.     

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.     

A theoretical assessment of the feasibility of potential Lunar Reconnaissance Orbiter radio occultation observations of the lunar ionosphere

The existence of a “dense” lunar ionosphere has been controversial for decades. Positive ions produced from the lunar surface and exosphere are inferred to have densities that are $?$${10}^6$ – ${10}^7$ m^?3 near the surface and smaller at higher altitudes, yet electron densities derived from radio occultation measurements occasionally exceed these values by orders of magnitude. For example, about 4% of the single-spacecraft radio occultation measurements from Kaguya/SELENE were consistent with peak electron densities of $~$$3\times {10}^8$ m^?3. Space plasmas should be neutral on macroscopic scales, so this represents a substantial discrepancy. Aditional observations of electron densities in the lunar ionosphere are critical to resolving this longstanding paradox. Here we theoretically assess whether radio occultation observations using two-way coherent S-band radio signals from the Lunar Reconnaissance Orbiter (LRO) spacecraft could provide useful measurements of electron densities in the lunar ionosphere. We predict the uncertainty in a single LRO radio occultation measurement of electron density to be $~$$3\times {10}^8$ m^?3, comparable to occasional observations by Kaguya/SELENE of a dense lunar ionosphere. Thus an individual profile from LRO is unlikely to reliably detect the lunar ionosphere; however, averages of multiple ($~$10) LRO profiles acquired under similar geophysical and viewing conditions should be able to make reliable detections. An observing rate of six ingress occultations per day ($~$2000 per year) could be achieved with minimal impact on current LRO operations. This rate compares favorably with the 378 observations reported from the single-spacecraft experiment on Kaguya/SELENE between November 2007 and June 2009. The large number of observations possible for LRO would be sufficient to permit wide-ranging investigations of spatial and temporal variations in the poorly understood lunar ionosphere. These findings strengthen efforts to conduct such observations with LRO.     

Study of aerosol anomaly associated with large earthquakes (M>6)

We present the variation of unusual atmospheric phenomena, aerosols, to understand the preseismic irregularities for two major earthquakes in Japan. We consider aerosol optical depth and Angstrom exponent data retrieved from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument onboard the Terra satellite to establish possible connections between earthquakes and the generation of aerosols. Variation of the aerosol parameters shows significant changes before the April 15, 2016, Kumamoto earthquake ($M=7.0,h=10$ km) and the November 21, 2016, Fukushima earthquake ($M=6.9$ and $h=9$ km), where M indicates the Richter magnitude and h indicates the focal depth. To identify the source of the aerosol particles, we use the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT-4). This model uses both Lagrangian and Eulerian approaches to compute trajectories and establish a source-receptor relationship. We compute backward trajectories to check whether the aerosol generated near the epicenter is due to the preseismic processes or is transported from other areas. From our results, we conclude the fine-mode aerosols are generated in the vicinity of the epicenter, 3–7 days before the earthquakes.     

Tumbling object deorbiting using spaceborne laser engagement – A CubeSat case study

This paper demonstrates active space debris removal using spaceborne laser systems. The laser beam and the surface of the target are discretised into multiple rays and finite elements, respectively, for laser-target interaction modelling, in which the laser ablation process is investigated. A high-fidelity attitude/orbit propagator tool is developed to account for both the linear impulse and angular impulse induced by the laser engagement and other perturbations. The laser system is activated only when three switch criteria are satisfied. In numerical simulations, laser pulses from international space station are generated to deorbit a 3U CubeSat with initially tumbling modes. The results validate the effectiveness of deorbiting tumbling CubeSats using spaceborne laser engagement, with the perigee height lowered by approximately $2.4\mathrm{km}$ in around $30\min$ after $2\mathrm{h}$ propagation. It is also found that the laser engagement becomes more effective for an initially faster rotating object.     

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.     

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.