Study of the radiation fields in LEO with the Space Application of Timepix Radiation Monitor (SATRAM)

We present the analysis of data taken by the Space Application of Timepix Radiation Monitor (SATRAM). It is centred on a Timepix detector (300?$\mathrm{\mu }$m thick silicon sensor, pixel pitch 55?$\mathrm{\mu }$m, 256?×?256 pixels). It was flown on Proba-V, an Earth observing satellite of the European Space Agency (ESA) from an altitude of 820?km on a sun-synchronous orbit, launched on May 7, 2013. A Monte Carlo simulation was conducted to determine the detector response to electrons (0.5–7?MeV) and protons (10–400?MeV) in an omnidirectional field taking into account the shielding of the detector housing and the satellite. With the help of the simulation, a strategy was developed to separate electrons, protons and ions in the data. The measured dose rate and stopping power distribution are presented as well as SATRAM’s capability to measure some of the stronger events in Earth’s magnetosphere. The stopping power, the cluster height and the shape of the particle tracks in the sensor were used to separate electrons, protons and ions. The results are presented as well. Finally, the pitch angles for a short period of time were extracted from the data and corrected with the angular response determined by the simulation.     

On the practical exploitation of perturbative effects in low Earth orbit for space debris mitigation

This paper presents the results of a numerical evaluation of the natural lifetime reduction in low Earth orbit, due to dynamical perturbations. The study considers two values for the area-to-mass ratio, a nominal ratio which resembles a typical value of spacecraft in orbit today, and an enhanced ratio which covers the surface augmentation. The results were obtained with two orbit propagators, one of a semi-analytical nature and the second one using non-averaged equations of motion. The simulations for both propagators were set up similarly to allow comparison. They both use the solar radiation pressure and the secular terms of the geopotential ($J_2,J_4$ and $J_6$). The atmospheric drag was turned on and off in both propagators to alternatively study the eccentricity build up and the residual lifetime. The non-averaging case also covers a validation with the full 6?×?6 geopotential. The results confirm the findings in previous publications, that is, the possibility for de-orbiting from altitudes above the residual atmosphere if a solar sail is deployed at the end-of-life, due to the combined effect of solar radiation pressure and the oblateness of the Earth. At near polar inclinations, shadowing effects can be exploited to the same end. The results obtained with the full, non-averaging propagator revealed additional de-orbiting corridors associated with solar radiation pressure which were not found by previous work on space debris mitigation. The results of both tools are compared for specific initial conditions. For nominal values of area-to-mass ratio, instead, it is confirmed that this resonance effect is negligible.The paper then puts the findings in the perspective of the current satellite catalogue. It identifies space missions which are currently close to a resonance corridor and shows the orbit evolution within the resonances with a significantly shorter residual orbital lifetime. The paper finishes with a discussion on the exploitation of these effects with regards to the long-term simulation of the space debris environment and a flux and collision probability comparison.     

Quantum dosimetry and online visualization of X-ray and charged particle radiation in commercial aircraft at operational flight altitudes with the pixel detector Timepix

We investigate the application of the hybrid semiconductor pixel detector Timepix for precise characterization, quantum sensitivity dosimetry and visualization of the charged particle radiation and X-ray field inside commercial aircraft at operational flight altitudes. The quantum counting capability and granularity of Timepix provides the composition and spectral-characteristics of the X-ray and charged-particle field with high sensitivity, wide dynamic range, high spatial resolution and particle type resolving power. For energetic charged particles the direction of trajectory and linear energy transfer can be measured. The detector is operated by the integrated readout interface FITPix for power, control and data acquisition together with the software package Pixelman for online visualization and real-time data processing. The compact and portable radiation camera can be deployed remotely being controlled simply by a laptop computer. The device performs continuous monitoring and accurate time-dependent measurements in wide dynamic range of particle fluxes, deposited energy, absorbed dose and equivalent dose rates. Results are presented for in-flight measurements at altitudes up to 12 km in various flights selected in the period 2006–2013.     

MONITORING OF ENERGETIC PARTICLE ENVIRONMENT INSIDE THE CHINA-BRAZIL EARTH RESOURCE SATELLITE

On 14 October 1999, the Chinese-Brazil earth resource satellite (CBERS-1) was launched in China. On board of the satellite there was an instrument designed at Peking University to detect the energetic particle radiation inside the satellite so the radiation fluxes of energetic particles in the cabin can be monitored continuously. Inside a satellite cabin, radiation environment consists of ether penetrated energetic particles or secondary radiation from satellite materials due to the interactions with primary cosmic rays.Purpose of the detectors are twofold, to monitor the particle radiation in the cabin and also to study the space radiation environment The data can be used to study the radiation environment and their effects on the electronics inside the satelhte cabin. On the other hand, the data are useful in study of geo-space energetic particle events such as solar proton events, particle precipitation and variations of the radiation belt since there should be some correlation between the radiation situation inside and outside the satellite.The instrument consists of two semi-conductor detectors for protons and electrons respectively. Each detector has two channels of energy ranges. They are 0.5-2MeV and ≥2MeV for electrons and 5-30MeV and 30-60MeV for protons. Counting rate for all channels are up to 104/(cm2@s)and power consumption is about 2.5 W. There are also the additional functions of CMOS TID (total integrated dose) effect and direct SEU monitoring. The data of CBMC was first sent back on Oct. 17 1999 and it's almost three years from then on. The detector has been working normally and the quality of data is good.The preliminary results of data analysis of CBMC not only reveal the effects of polar particle precipitation and radiation belt on radiation environment inside a satellite, but also show some important features of the geo-space energetic particle radiation.As one of the most important parameters of space weather, the energetic charged particles have great influences on space activities and ground tech nology. CBMC is perhaps the first long-term on-board special equipment to monitor the energetic particle radiation environment inside the satellite and the data it accnmulated are very useful in both satellite designing and space research.     

Depth dose measurements with the Liulin-5 experiment inside the spherical phantom of the MATROSHKA-R project onboard the International Space Station

The Liulin-5 experiment is a part of the international project MATROSHKA-R on the Russian segment of the ISS, which uses a tissue-equivalent spherical phantom equipped with a set of radiation detectors. The objective of the MATROSHKA-R project is to provide depth dose distribution of the radiation field inside the sphere in order to get more information on the distribution of dose in a human body. Liulin-5 is a charged particle telescope using three silicon detectors. It measures time resolved energy deposition spectra, linear energy transfer (LET) spectra, particle flux, and absorbed doses of electrons, protons and heavy ions, simultaneously at three depths along the radius of the phantom. Measurements during the minimum of the solar activity in cycle 23 show that the average absorbed daily doses at 40 mm depth in the phantom are between 180 μGy/day and 220 μGy/day. The absorbed doses at 165 mm depth in the phantom decrease by a factor of 1.6–1.8 compared to the doses at 40 mm depth due to the self-shielding of the phantom from trapped protons. The average dose equivalent at 40 mm depth is 590 ± 32 μSV/day and the galactic cosmic rays (GCR) contribute at least 70% of the total dose equivalent at that depth. Shown is that due to the South Atlantic Anomaly (SAA) trapped protons asymmetry and the direction of Liulin-5 lowest shielding zone the dose rates on ascending and descending nodes in SAA are different. The data obtained are compared to data from other radiation detectors on ISS.     

Simulations of absorbed dose on the phantom surface of MATROSHKA-R experiment at the ISS

The health risks associated with exposure to various components of space radiation are of great concern when planning manned long-term interplanetary missions, such as future missions to Mars. Since it is not possible to measure the radiation environment inside of human organs in deep space, simulations based on radiation transport/interaction codes coupled to phantoms of tissue equivalent materials are used. However, the calculated results depend on the models used in the codes, and it is therefore necessary to verify their validity by comparison with measured data. The goal of this paper is to compare absorbed doses obtained in the MATROSHKA-R experiment performed at the International Space Station (ISS) with simulations performed with the three-dimensional Monte Carlo Particle and Heavy-Ion Transport code System (PHITS). The absorbed dose was measured using passive detectors (packages of thermoluminescent and plastic nuclear track detectors) placed on the surface of the spherical tissue equivalent phantom MATROSHKA-R, which was exposed aboard the ISS in the Service Zvezda Module from December 2005 to September 2006. The data calculated by PHITS assuming an ISS shielding of 3 g/cm² and 5 g/cm² aluminum mass thickness were in good agreement with the measurements. Using a simplified geometrical model of the ISS, the influence of variations in altitude and wall mass thickness of the ISS on the calculated absorbed dose was estimated. The uncertainties of the calculated data are also discussed; the relative expanded uncertainty of absorbed dose in phantom was estimated to be 44% at a 95% confidence level.     

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.     

Time-varying solar radiation pressure on Ajisai in comparison with LAGEOS satellites

This study aims to investigate solar radiation pressure acting on the spherical geodetic satellites, Ajisai, LAGEOS-1, and LAGEOS-2. The solar radiation pressure coefficients ($C_R$) are derived in two independent ways: (a) through precise orbit determination (POD) using satellite laser ranging (SLR) data, and (b) through modeling using the optical properties of the satellite surface material. The average $C_R$ value of Ajisai (1.039), as calculated from the time series of $C_R$ POD estimates every 15?days, is consistently smaller than those of LAGEOS-1 (1.140) and LAGEOS-2 (1.103). This difference can be explained by the fact that the surface of Ajisai is mostly covered by mirrors. The Ajisai $C_R$ values estimated by POD show remarkable semi-annual variation, which disagrees with the results of a previous study (Sengoku et al., 1995) predicting that the $C_R$ of Ajisai varies almost annually. We attribute this semi-annual variation to two physical reasons: the non-spherical additional cross-sectional area due to the “attached fitting ring” and the low reflectivity of the surface material in the polar regions. Furthermore, the solar radiation pressure acting on Ajisai varies annually in a direction perpendicular to the sun-satellite vector. Finally, the two independent $C_R$ values of Ajisai agree more when we assume a total solar irradiance of 1361?W/m², whereas the value 1367?W/m² has been commonly used in POD.     

Fast Probabilistic Particle Identification algorithm using silicon strip detectors

Active detectors used as radiation monitors in space are not usually able to perform Particle Identification (PID). Common techniques need energy loss spectra with high statistics to estimate ion abundances. The ALTEA-space detector system is a set of silicon strip particle telescopes monitoring the radiation environment on board the International Space Station since July 2006 with real-time telemetry capabilities. Its large geometrical factor due to the concurrent use of six detectors permits the acquisition of good energy loss spectra even in a short period of observation. In this paper we present a novel Fast Probabilistic Particle Identification (FPPI) algorithm developed for the ALTEA data analysis in order to perform nuclear identification with low statistics and, with some limitations, also in real time.     

Optimal Earth’s reentry disposal of the Galileo constellation

Nowadays there is international consensus that space activities must be managed to minimize debris generation and risk. The paper presents a method for the end-of-life (EoL) disposal of spacecraft in Medium Earth Orbit (MEO). The problem is formulated as a multiobjective optimisation one, which is solved with an evolutionary algorithm. An impulsive manoeuvre is optimised to reenter the spacecraft in Earth’s atmosphere within 100?years. Pareto optimal solutions are obtained using the manoeuvre $\mathrm{\Delta }v$ and the time-to-reentry as objective functions to be minimised. To explore at the best the search space a semi-analytical orbit propagator, which can propagate an orbit for 100?years in few seconds, is adopted. An in-depth analysis of the results is carried out to understand the conditions leading to a fast reentry with minimum propellant. For this aim a new way of representing the disposal solutions is introduced. With a single 2D plot we are able to fully describe the time evolution of all the relevant orbital parameters as well as identify the conditions that enables the eccentricity build-up. The EoL disposal of the Galileo constellation is used as test case.     

Upgrading CCIR's foF2 maps using available ionosondes and genetic algorithms

We have developed a new approach towards a new database of the ionospheric parameter $f_{o}F2$. This parameter, being the frequency of the maximum of the ionospheric electronic density profile and its main modeller, is of great interest not only in atmospheric studies but also in the realm of radio propagation. The current databases, generated by CCIR (Committee Consultative for Ionospheric Radiowave propagation) and URSI (International Union of Radio Science), and used by the IRI (International Reference Ionosphere) model, are based on Fourier expansions and have been built in the 60s from the available ionosondes at that time. The main goal of this work is to upgrade the databases by using new available ionosonde data. To this end we used the IRI diurnal/spherical expansions to represent the $f_{o}F2$ variability, and computed its coefficients by means of a genetic algorithm (GA). In order to test the performance of the proposed methodology, we applied it to the South American region with data obtained by RAPEAS (Red Argentina para el Estudio de la Atmósfera Superior, i.e. Argentine Network for the Study of the Upper Atmosphere) during the years 1958–2009. The new GA coefficients provide a global better fit of the IRI model to the observed $f_{o}F2$ than the CCIR coefficients. Since the same formulae and the same number of coefficients were used, the overall integrity of IRI’s typical ionospheric feature representation was preserved. The best improvements with respect to CCIR are obtained at low solar activities, at large (in absolute value) modip latitudes, and at night-time. The new method is flexible in the sense that can be applied either globally or regionally. It is also very easy to recompute the coefficients when new data is available. The computation of a third set of coefficients corresponding to days of medium solar activity in order to avoid the interpolation between low and high activities is suggested. The same procedure as for $f_{o}F2$ can be perfomed to obtain the ionospheric parameter M(3000)F2.     

Direct comparison of transient radiation belt topology and dynamics in 1991 based on measurements onboard Mir space station and NOAA satellite.

In March 1991 the CRRES spacecraft measured a new transient radiation belt resulting from a solar proton event and subsequent geomagnetic disturbance. The presence of this belt was also noted by dosimeter-radiometers aboard the Mir space station (approx. 400 km, 51 degrees orbit) and by particle telescopes on the NOAA-10 spacecraft (850 km, 98 degrees). This event provides a unique opportunity to compare particle flux and dose measurements made by different instruments in different orbits under changing conditions. We present here a comparison of the measurements made by the different detectors. We discuss the topology and dynamics of the transient radiation belt over a period of more than one year.     

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 }$.     

GRACE accelerometer calibration by high precision non-gravitational force modeling

The precise modeling and knowledge of non-gravitational forces acting on satellites is of big interest to many scientific tasks and missions. Since 2002, the twin GRACE satellites have measured these forces in a low Earth orbit with highly precise accelerometers, for about 15?years. Besides the significance for the GRACE mission, these measurement data allow the evaluation of modeling approaches and the improvement of force models. Unfortunately, before any scientific usage, the accelerometer measurements need to be calibrated, namely scale factor and bias have to be regularly estimated.In this study we demonstrate an accelerometer calibration approach, solely based on high precision non-gravitational force modeling without any use of empirically or stochastically estimated parameters, using our in-house developed satellite simulation tool XHPS. The aim of this work is twofold, first we use the accelerometer data and the residuals resulting from the calibration to quantitatively analyze and validate different non-gravitational force model approaches. In a second step, we compare the calibration results to three different calibration methods from different authors, based on gravity field recovery, GPS-based precise orbit determination, and based on modeled accelerations.We consider atmospheric drag forces and winds, as well as radiation forces due to solar radiation pressure, albedo, Earth infrared and thermal radiation (TRP) of the satellite itself. For TRP, we investigate different transient temperature calculation approaches for the satellite surfaces with absorbed power from the aforementioned radiation sources. A detailed finite element model of the satellite is utilized for every force, considering orientation, material properties and shadowing conditions for each element.For cross-track and radial direction, which are mainly affected by the radiative forces, our calibration residuals are quite small when drag is not super dominant (1–3?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$50?$\text{nm}/{\text{s}}^2$). For these directions the calibration seems to perform better than the other compared methods, where some bigger differences were found. For the drag dominated along-track direction it is vice versa, here our method is not sensitive enough because the difference between modeled and measured drag is bigger (e.g. residuals around 10?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$70?$\text{nm}/{\text{s}}^2$ for low solar activity). In along-track direction the orbit determination based methods are more sensitive and produce more reliable results. Results for the complete GRACE mission time span from 2003 to 2017 are shown, covering different seasonal environmental conditions.     

Effects of dust polarity and nonextensive electrons on the dust-ion acoustic solitons and double layers in earth atmosphere

Propagation of dustion acoustic solitary waves (DIASWs) and double layers is discussed in earth atmosphere, using the Sagdeev potential method. The best model for distribution function of electrons in earth atmosphere is found by fitting available data on different distribution functions. The nonextensive function with parameter $q=0.58$ provides the best fit on observations. Thus we analyze the propagation of localized waves in an unmagnetized plasma containing nonextensive electrons, inertial ions, and negatively/positively charged stationary dust. It is found that both compressive and rarefactive solitons as well as double layers exist depending on the sign (and the value) of dust polarity. Characters of propagated waves are described using the presented model.     

Dynamical cartography of Earth satellite orbits

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term ($～$120?years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes ($～$400?km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded “ReDSHIFT” project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an “atlas” of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth’s oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device. We find that this would generally promote the deorbiting process, particularly at the transition region between LEO and MEO. Although direct reentry from very low eccentricities is very unlikely in most cases of interest, we find that a modest “delta-v” ($\mathrm{\Delta }V$) budget would be enough for satellites to be steered into a relatively short-lived resonance and achieve reentry into the Earth’s atmosphere within reasonable timescales ($～$50?years).     

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.     

Radiation measured during ISS-Expedition 13 with different dosimeters

Radiation in low Earth orbit (LEO) is mainly composed of galactic cosmic rays (GCR), solar energetic particles and particles in SAA (South Atlantic Anomaly). The biological impact of space radiation to astronauts depends strongly on the particles’ linear energy transfer (LET) and is dominated by high LET radiation. It is important to measure the LET spectrum for the space radiation field and to investigate the influence of radiation on astronauts. At present, the preferred active dosimeters sensitive to all LET are the tissue equivalent proportional counter (TEPC) and the silicon detectors in various configurations; the preferred passive dosimeters are CR-39 plastic nuclear track detectors (PNTDs) sensitive to high LET and thermoluminescence dosimeters (TLDs) as well as optically stimulated luminescence dosimeters (OSLDs) sensitive to low LET. The TEPC, CR-39 PNTDs, TLDs and OSLDs were used to investigate the radiation field for the ISS mission Expedition 13 (ISS-12S) in LEO. LET spectra and radiation quantities (fluence, absorbed dose, dose equivalent and quality factor) were measured for the space mission with different dosimeters. This paper introduces the role of high LET radiation in radiobiology, the operational principles for the different dosimeters, the LET spectrum method using CR-39 detectors, the method to combine the results measured with TLDs/OSLDs and CR-39 PNTDs, and presents the LET spectra and the radiation quantities measured and combined.     

Real-time spacecraft intercept strategy on J2-perturbed orbits

A real-time intercept strategy for spacecraft under the non-uniform gravitational perturbation of Earth is addressed in this paper. To intercept a target spacecraft on general conic sections, an interceptor considered in this work makes use of a thruster propelling the constant thrust which is comparable to unrealistic impulse-type thrust. The $J_2$ perturbation introduces critical dynamic variations of spacecraft orbiting the Earth, which results in a considerable amount of position error of the interceptor at the final intercept point. In order to release the burden of $J_2$ disturbance and make the miss distance between the target and interceptor small, a real-time intercept technique with an optimal intercept algorithm is suggested. The strategy proposed is to obtain an optimized output iteratively for a given time interval with previously obtained optimal values. These parameters are evaluated by the optimal intercept algorithm suggested. Once the optimal velocity change is obtained to satisfy intercept requirements, although the orbital system is perturbed, it is easy to regenerate a new solution by setting the previous solution as new initial guesses. This strategy is employed iteratively until the interceptor meets the target. Several numerical simulations are performed to highlight the proposed real-time strategy for spacecraft intercept missions.