Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

Statistical evaluation of GLONASS amplitude scintillation over low latitudes in the Brazilian territory

The ionospheric scintillation, generated by the ionospheric plasma irregularities, affects the radio signals that pass through it. Their effects are widely studied in the literature with two different approaches. The first one deals with the use of radio signals to study and understand the morphology of this phenomenon, while the second one seeks to understand and model how much this phenomenon interferes in the radio signals and consequently in the services to which these systems work. The interest of several areas, particularly to those that are life critical, has increased using the concept of satellite multi-constellation, which consists of receiving, processing and using data from different navigation and positioning systems. Although there is a vast literature analyzing the effects of ionospheric scintillation on satellite navigation systems, the number of studies using signals received from the Russian satellite positioning system (named GLONASS) is still very rare. This work presents for the first time in the Brazilian low-latitude sector a statistical analysis of ionospheric scintillation data for all levels of magnetic activities obtained by a set of scintillation monitors that receive signals from the GLONASS system. In this study, data collected from four stations were used in the analysis; Fortaleza, Presidente Prudente, São José dos Campos and Porto Alegre. The GLONASS L-band signals were analyzed for the period from December 21, 2012 to June 20, 2016, which includes the peak of the solar cycle 24 that occurred in 2014. The main characteristics of scintillation presented in this study include: (1) the statistical evaluation of seasonal and solar activity, showing the chances that an user on similar geophysical conditions may be susceptible to the effects of ionospheric scintillation; (2) a temporal analysis based on the local time distribution of scintillation at different seasons and intensity levels; and (3) the evaluation of number of simultaneously affected channels and its effects on the dilution of precision (DOP) for GNSS users are also presented in order to alert the timetables in which navigation will be most susceptible to such effects, as well as statistics on simultaneously affected channels. Relevant results about these statistical characteristics of scintillation are presented and analyzed providing relevant information about availability of a navigation system.     

OSIRIS-REx Flight Dynamics and Navigation Design

OSIRIS-REx is the first NASA mission to return a sample of an asteroid to Earth. Navigation and flight dynamics for the mission to acquire and return a sample of asteroid 101955 Bennu establish many firsts for space exploration. These include relatively small orbital maneuvers that are precise to ～1 mm/s, close-up operations in a captured orbit about an asteroid that is small in size and mass, and planning and orbit phasing to revisit the same spot on Bennu in similar lighting conditions. After preliminary surveys and close approach flyovers of Bennu, the sample site will be scientifically characterized and selected. A robotic shock-absorbing arm with an attached sample collection head mounted on the main spacecraft bus acquires the sample, requiring navigation to Bennu’s surface. A touch-and-go sample acquisition maneuver will result in the retrieval of at least 60 grams of regolith, and up to several kilograms. The flight activity concludes with a return cruise to Earth and delivery of the sample return capsule (SRC) for landing and sample recovery at the Utah Test and Training Range (UTTR).     

On the detection of anomalous seismo-ionospheric behavior in the presence of space weather stimuli for large earthquakes

Anomalous behavior of ionospheric total electron content (TEC) prior to earthquake has been observed in many studies. Evidence of such seismo-ionospheric coupling effects suggests that it is plausible to rely on TEC signatures for early earthquake warning. However, the detection of pre-earthquake TEC anomalies (PETA) has not been adopted in practice due to two pertinent issues. Firstly, the effects of space weather activity can affect TEC levels and cause anomalous behavior in the TEC. Usually arbitrary thresholds are set for space weather indices to eliminate TEC anomaly due to space weather effects. Secondly, the choice regarding moving time-window length used to characterise background variation of TEC within the statistical envelope approach has an effect on detection of PETA. While the rule-of-thumb in selecting the moving window length is to have a time window capable of capturing background variability and short-term fluctuations, the length of the time window used in the literature varies with little justification. In this study, a critical examination is conducted on the statistical envelope approach and in particular, to eliminate the effect of space weather activity without the use of arbitrary space indices to detect PETA. A two-part PETA identification procedure is proposed, consisting of wavelet analyses isolating non-earthquake TEC contributions, followed by the statistical envelope method using a moving window length standardized based on observed periodicities and statistical implications is suggested. The approach is tested on a database of 30 large earthquakes (M?≥?7.0). The proposed procedure shows that PETA can be detected prior to earthquakes at higher confidence levels without the need to separately check for space weather activity. More importantly, the procedure was able to detect PETA for studies where it was previously reported that PETA could not be detected or detected convincingly.     

Characterization of the low latitude plasma density irregularities observed using C/NOFS and SCINDA data

Complex electrodynamic processes over the low latitude region often result in post sunset plasma density irregularities which degrade satellite communication and navigation. In order to forecast the density irregularities, their occurrence time, duration and location need to be quantified. Data from the Communication/Navigation Outage Forecasting System (C/NOFS) satellite was used to characterize the low latitude ion density irregularities from 2011 to 2013. This was supported by ground based data from the SCIntillation Network Decision Aid (SCINDA) receivers at Makerere (Geographic coordinate 32.6°E, 0.3°N, and dip latitude ?9.3°N) and Nairobi (Geographic coordinate 36.8°E, ?1.3°N, and dip latitude ?10.8°N). The results show that irregularities in ion density have a daily pattern with peaks from 20:00 to 24:00 Local Time (LT). Scintillation activity at L band and VHF over East Africa peaked in 2011 and 2012 from 20:00 to 24:00 LT, though in many cases scintillation at VHF persisted longer than that at L band. A longitudinal pattern in ion density irregularity occurrence was observed with peaks over 135–180°E and 270–300°E. The likelihood of ion density irregularity occurrence decreased with increasing altitude. Analysis of C/NOFS zonal ion drift velocities showed that the largest nighttime and daytime drifts were in 270–300°E and 300–330°E longitude regions respectively. Zonal irregularity drift velocities over East Africa were for the first time estimated from L-band scintillation indices. The results show that the velocity of plasma density irregularities in 2011 and 2012 varied daily, and hourly in the range of 50–150 m s^?1. The zonal drift velocity estimates from the L-band scintillation indices had good positive correlation with the zonal drift velocities derived from VHF receivers by the spaced receiver technique.     

Longitudinal variations of geomagnetic and ionospheric parameters in the Northern Hemisphere during magnetic storms according to multi-instrument observations

We present a joint analysis of longitude-temporal variations of ionospheric and geomagnetic parameters at middle and high latitudes in the Northern Hemisphere during the two severe magnetic storms in March and June 2015 by using data from the chains of magnetometers, ionosondes and GPS/GLONASS receivers. We identify the fixed longitudinal zones where the variability of the magnetic field is consistently high or low under quiet and disturbed geomagnetic conditions. The revealed longitudinal structure of the geomagnetic field variability in quiet geomagnetic conditions is caused by the discrepancy of the geographic and magnetic poles and by the spatial anomalies of different scales in the main magnetic field of the Earth. Variations of ionospheric parameters are shown to exhibit a pronounced longitudinal inhomogeneity with changing geomagnetic conditions. This inhomogeneity is associated with the longitudinal features of background and disturbed structure of the geomagnetic field. During the recovery phase of a storm, important role in dynamics of the mid-latitude ionosphere may belong to wave-like thermospheric disturbances of molecular gas, propagating westward for several days. Therefore, it is necessary to extend the time interval for studying the ionospheric effects of strong magnetic storms by a few days after the end of the magnetospheric source influence, while the disturbed regions in the thermosphere continues moving westward and causes the electron density decrease along the trajectories of propagation.