A comparison of radiation shielding effectiveness of materials for highly elliptical orbits

The Canadian Space Agency (CSA) has proposed a Polar Communications and Weather (PCW) satellite mission, in conjunction with other partners. The PCW will provide essential communications and meteorological services to the Canadian Arctic, as well as space weather observations of in situ ionizing radiation along the orbit. The CSA has identified three potential Highly Elliptical Orbits (HEOs) for a PCW satellite constellation, Molniya, Modified Tundra and Triple Apogee (TAP), each having specific merits, which would directly benefit the performance/longevity of a PCW spacecraft. Radiation shielding effectiveness of various materials was studied for the three PCW orbit options to determine the feasibility of employing materials other than conventional aluminium to achieve a specified spacecraft shielding level with weight savings over aluminium. It was found that, depending on the orbit-specific radiation environment characteristics, the benefits of using polyethylene based materials is significant enough (e.g., 22% in Molniya for PE at 50 krad TID) to merit further investigation.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

太阳同步轨道真空紫外光电倍增管的高能粒子屏蔽方法研究

空间辐射环境是导致光电倍增管性能衰变甚至失效的主要原因, 由于空间辐射环境的复杂性以及飞行器在空间活动的轨道、设计寿命和光电倍增管的类型不同, 光电倍增管需要的抗辐射屏蔽要求不同. 以太阳同步轨道高能粒子辐射环境为例, 对真空紫外光电倍增管进行了辐射屏蔽的Monte-Carlo计算, 在此基础上提出了运行于该轨道上的真空紫外光电倍增管的高能粒子屏蔽优化组合. 模拟计算表明, 采用Al-W-Al三层结构的屏蔽组合形式, 其中Al层厚2mm, W厚 0.5mm,可将辐射剂量从原本的9.98× 106rad (Si)降低到 8.41× 102rad (Si), 能较好地满足实际使用要求.      

星载辐照剂量深度分布探测仪设计

空间带电粒子对在轨航天器会产生辐照剂量效应,严重时可导致星载设备及航天器等的性能衰减及寿命下降,因此采用了星内多点多方位的辐照总剂量探测技术.在中国首次采用深度剂量监测方案.每个探头设置5个剂量监测点,对应屏蔽厚度分别为0（开窗）,1mm,2.5mm,3.5mm,6mm的铝结构,探测剂量总量程为2×106rad （Si）,最高灵敏度为1rad （Si）.所设计的星载辐照剂量深度分布探测仪可以对卫星在轨遭受的辐照剂量进行实时监测和预警,其探测结果用于研究星内剂量深度分布,对卫星载荷在轨工作状态进行分析评估,并为后续卫星的抗辐照加固设计提供依据.      

Accurate orbit predictions for debris orbit manoeuvre using ground-based lasers

Orbit manoeuvre of low Earth orbiting (LEO) debris using ground-based lasers has been proposed as a cost-effective means to avoid debris collisions. This requires the orbit of the debris object to be determined and predicted accurately so that the laser beam can be locked on the debris without the loss of valuable laser operation time. This paper presents the method and results of a short-term accurate LEO (<900 km in altitude) debris orbit prediction study using sparse laser ranging data collected by the EOS Space Debris Tracking System (SDTS). A main development is the estimation of the ballistic coefficients of the LEO objects from their archived long-term two line elements (TLE). When an object is laser tracked for two passes over about 24 h, orbit prediction (OP) accuracy of 10–20 arc seconds for the next 24–48 h can be achieved – the accuracy required for laser debris manoeuvre. The improvements in debris OP accuracy are significant in other applications such as debris conjunction analyses and the realisation of daytime debris laser tracking.     

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.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Comparison of the transport codes HZETRN,HETC and FLUKA for a solar particle event

The protection of astronauts and instrumentation from galactic cosmic rays and solar particle events is one of the primary constraints associated with mission planning in low earth orbit or deep space. To help satisfy this constraint, several computational tools have been developed to analyze the effectiveness of various shielding materials and structures exposed to space radiation. These tools are now being carefully scrutinized through a systematic effort of verification, validation, and uncertainty quantification. In this benchmark study, the deterministic transport code HZETRN is compared to the Monte Carlo transport codes HETC-HEDS and FLUKA for a 30 g/cm² water target protected by a 20 g/cm² aluminum shield exposed to a parameterization of the February 1956 solar particle event. Neutron and proton fluences as well as dose and dose equivalent are compared at various depths in the water target. The regions of agreement and disagreement between the three codes are quantified and discussed, and recommendations for future work are given.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

Method for the prediction of the effective dose equivalent to the crew of the International Space Station

This paper describes a methodology for assessing the pre-mission exposure of space crew aboard the International Space Station (ISS) in terms of an effective dose equivalent. In this approach, the PHITS Monte Carlo code was used to assess the particle transport of galactic cosmic radiation (GCR) and trapped radiation for solar maximum and minimum conditions through an aluminum shield thickness. From these predicted spectra, and using fluence-to-dose conversion factors, a scaling ratio of the effective dose equivalent rate to the ICRU ambient dose equivalent rate at a 10 mm depth was determined. Only contributions from secondary neutrons, protons, and alpha particles were considered in this analysis.     

Heliotropic dust rings for Earth climate engineering

This paper examines the concept of a Sun-pointing elliptical Earth ring comprised of dust grains to offset global warming. A new family of non-Keplerian periodic orbits, under the effects of solar radiation pressure and the Earth’s J₂ oblateness perturbation, is used to increase the lifetime of the passive cloud of particles and, thus, increase the efficiency of this geoengineering strategy. An analytical model is used to predict the orbit evolution of the dust ring due to solar-radiation pressure and the J₂ effect. The attenuation of the solar radiation can then be calculated from the ring model. In comparison to circular orbits, eccentric orbits yield a more stable environment for small grain sizes and therefore achieve higher efficiencies when the orbit decay of the material is considered. Moreover, the novel orbital dynamics experienced by high area-to-mass ratio objects, influenced by solar radiation pressure and the J₂ effect, ensure the ring will maintain a permanent heliotropic shape, with dust spending the largest portion of time on the Sun facing side of the orbit. It is envisaged that small dust grains can be released from a circular generator orbit with an initial impulse to enter an eccentric orbit with Sun-facing apogee. Finally, a lowest estimate of 1 × 10¹² kg of material is computed as the total mass required to offset the effects of global warming.     

Trapped proton fluxes measured on board LEO satellites in comparison with models

The proton fluxes from the low-Earth orbital satellites databases (NPOES-17 and CORONAS-F) were analyzed for the quiet geomagnetic period in April 2005. The satisfactory consent was found between the experimental and the AP8 model fluxes of the trapped protons with energy more than ∼10 MeV. At the same time, trapped proton fluxes with energy less than 10 MeV measured by LEO satellites were higher than the ones predicted by the AP8 model in the region of the SAA (drift shell L < 1.5).     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

Analysis of the space radiation doses obtained simultaneously at two different locations outside the ISS

Space weather and related ionizing radiation has been recognized as one of the main health concerns for the International Space Station (ISS) crew. The estimation of the radiation effect on humans outside the ISS requires at first order accurate knowledge of their accumulated absorbed dose rates, which depend on the global space radiation distribution, solar cycle and local variations generated by the 3D mass distribution surrounding the ISS. The R3DE (Radiation Risks Radiometer-Dosimeter for the EXPOSE-E platform) on the European Technological Exposure Facility (EuTEF) worked successfully outside of the European Columbus module between February 2008 and September 2009. A very similar instrument named R3DR for the EXPOSE-R platform worked outside the Russian Zvezda module of the ISS between March 2009 and August 2010. Both are Liulin-type detectors, Bulgarian-built miniature spectrometer-dosimeters. The acquired approximately 5 million deposited energy spectra from which the flux and absorbed dose rate were calculated with 10 s resolution behind less than 0.41 g cm⁻² shielding. This paper analyses the spectra collected in 2009 by the R3DE/R instruments and the long-term variations in the different radiation environments of Galactic Cosmic Rays (GCR), inner radiation belt trapped protons in the region of the South Atlantic Anomaly (SAA) and relativistic electrons from the Outer Radiation Belt (ORB). The R3DE instrument, heavily shielded by the surrounding structures, measured smaller primary fluxes and dose rates from energetic protons from the SAA and relativistic electrons from the ORB but higher values from GCRs because of the contribution from secondary particles. The main conclusion from this investigation is that the dose rates from different radiation sources around the International Space Station (ISS) have a large special and temporal dynamic range. The collected data can be interpreted as possible doses obtained by the cosmonauts and astronauts during Extra Vehicular Activities (EVA) because the R3DE/R instruments shielding is very similar to the Russian and American space suits average shielding (,  and ). Fast, active measurements are required to assess accurately the dose accumulated by astronauts during EVA.     

Debris flux comparisons from the Goldstone Radar,Haystack Radar,and Hax Radar prior,during, and after the last solar maximum

The continual monitoring of the low Earth orbit (LEO) debris environment using highly sensitive radars is essential for an accurate characterization of these dynamic populations. Debris populations are continually evolving since there are new debris sources, previously unrecognized debris sources, and debris loss mechanisms that are dependent on the dynamic space environment. Such radar data are used to supplement, update, and validate existing orbital debris models. NASA has been utilizing radar observations of the debris environment for over a decade from three complementary radars: the NASA JPL Goldstone radar, the MIT Lincoln Laboratory (MIT/LL) Long Range Imaging Radar (known as the Haystack radar), and the MIT/LL Haystack Auxiliary radar (HAX). All of these systems are highly sensitive radars that operate in a fixed staring mode to statistically sample orbital debris in the LEO environment. Each of these radars is ideally suited to measure debris within a specific size region. The Goldstone radar generally observes objects with sizes from 2 mm to 1 cm. The Haystack radar generally measures from 5 mm to several meters. The HAX radar generally measures from 2 cm to several meters. These overlapping size regions allow a continuous measurement of cumulative debris flux versus diameter from 2 mm to several meters for a given altitude window. This is demonstrated for all three radars by comparing the debris flux versus diameter over 200 km altitude windows for 3 nonconsecutive years from 1998 to 2003. These years correspond to periods before, during, and after the peak of the last solar cycle. Comparing the year to year flux from Haystack for each of these altitude regions indicate statistically significant changes in subsets of the debris populations. Potential causes of these changes are discussed. These analysis results include error bars that represent statistical sampling errors.     

BeiDou-3 orbit and clock quality of the IGS Multi-GNSS Pilot Project

Within the Multi-GNSS Pilot Project (MGEX) of the International GNSS Service (IGS), precise orbit and clock products for the BeiDou-3 global navigation satellite system (BDS-3) are routinely generated by a total of five analysis centers. The processing standards and specific properties of the individual products are reviewed and the BDS-3 orbit and clock product performance is assessed through direct inter-comparison, satellite laser ranging (SLR) residuals, clock stability analysis, and precise point positioning solutions. The orbit consistency evaluated by the signal-in-space range error is on the level of 4–8 cm for the medium Earth orbit satellites whereas SLR residuals have RMS values between 3 and 9 cm. The clock analysis reveals sytematic effects related to the elevation of the Sun above the orbital plane for all ACs pointing to deficiencies in solar radiation pressure modeling. Nevertheless, precise point positioning with the BDS-3 MGEX orbit and clock products results in 3D RMS values between 7 and 8 mm.     

On the validity of the aluminum equivalent approximation in space radiation shielding applications

The origin of the aluminum equivalent shield approximation in space radiation analysis can be traced back to its roots in the early years of the NASA space programs (Mercury, Gemini and Apollo) wherein the primary radiobiological concern was the intense sources of ionizing radiation causing short term effects which was thought to jeopardize the safety of the crew and hence the mission. Herein, it is shown that the aluminum equivalent shield approximation, although reasonably well suited for that time period and to the application for which it was developed, is of questionable usefulness to the radiobiological concerns of routine space operations of the 21st century which will include long stays onboard the International Space Station (ISS) and perhaps the moon. This is especially true for a risk based protection system, as appears imminent for deep space exploration where the long-term effects of Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis demonstrates that sufficiently large errors in the interior particle environment of a spacecraft result from the use of the aluminum equivalent approximation, and such approximations should be avoided in future astronaut risk estimates. In this study, the aluminum equivalent approximation is evaluated as a means for estimating the particle environment within a spacecraft structure induced by the GCR radiation field. For comparison, the two extremes of the GCR environment, the 1977 solar minimum and the 2001 solar maximum, are considered. These environments are coupled to the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport (HZETRN), which propagates the GCR spectra for elements with charges (Z) in the range 1 ? Z ? 28 (H–Ni) and secondary neutrons through selected target materials. The coupling of the GCR extremes to HZETRN allows for the examination of the induced environment within the interior of an idealized spacecraft as approximated by a spherical shell shield, and the effects of the aluminum equivalent approximation for a good polymeric shield material such as generic polyethylene (PE). The shield thickness is represented by a 25 g/cm² spherical shell. Although, one could imagine the progression to greater thickness, the current range will be sufficient to evaluate the qualitative usefulness of the aluminum equivalent approximation. Upon establishing the inaccuracies of the aluminum equivalent approximation through numerical simulations of the GCR radiation field attenuation for PE and aluminum equivalent PE spherical shells, we further present results for a limited set of commercially available, hydrogen rich, multifunctional polymeric constituents to assess the effect of the aluminum equivalent approximation on their radiation attenuation response as compared to the generic PE.     

Near-Earth radiation model deficiencies as seen on CRRES.

The Space Radiation (SPACERAD) experiments on the Combined Release and Radiation Effects Satellite (CRRES) gathered 14 months of radiation particle data in an 18 degrees inclination orbit between 350 km and 36000 km from July 1990 to October 1991. When compared to the NASA radiation belt models AP8 and AE8, the data show the proton model (AP8) does not take into account a second belt formed after major solar flare/shock injection events, and the electron model (AE8) is misleading, at best, in calculating dose in near-Earth orbits. The second proton belt, although softer in energy than the main proton belt, can produce upsets in proton sensitive chips and would produce significant dose in satellites orbiting in it. The MeV electrons observed on CRRES show a significant particle population above 5 MeV (not in the AE8 model) which must be included in any meaningful dose predictions for satellites operating between L-shells of 1.7 and 3.0 RE.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.