A robotic telescope network for space debris identification and tracking

The two TAROT (Télescopes à Action Rapide pour les Objets Transitoires; Rapid Action Telescopes for Transient Objects) installations are fully robotic optical observatories with optimized observation scheduling, data processing and archiving. Zadko is a 1 m telescope in Western Australia. The fully robotisation of the Zadko telescope has just been completed; it is now included in the TAROT network. In this paper we provide an overview of this international network of robotic optical telescopes. We discuss the advantages of using the network to participate in a satellite and space debris tracking program. This network will access almost all geostationary belt objects, and provide the first real-time satellite positioning capability. The inclusion of the 1 m Zadko telescope into the network significantly extends the efficiency and sensitivity of the existing two telescope configuration.     

Airglow OH emission height inferred from the OH temperature and meteor trail diffusion coefficient

Simultaneous observations of the airglow OH(6,2) band rotational temperature, TOH, and meteor trail ambipolar diffusion coefficient, D, were carried out at Shigaraki (35°N, 136°E), during PSMOS 2003 Campaign, January 28 to February 8, 2003. The OH emission height was estimated by cross correlation analysis of the TOH and D nocturnal variations. A good correlation between TOH and D was obtained at 85 km of altitude. From the nocturnal variations of TOH and D, it is found that the OH emission peak height varied from 88 km before the midnight to 84 km in the early morning. The height variation could be caused by an atmospheric tidal effect in the emission height.     

瑞利激光雷达反演中层大气温度算法

瑞利激光雷达是探测中层大气温度的重要设备,具有高时空分辨率且不存在探测盲区的优点.为充分利用激光雷达探测到的原始数据,改进了传统的Chanin-Haunchecorne方法,采用由均匀搜索生成温度初值的方式反演大气温度.位于北京延庆（40.3°N,116.2°E）的瑞利激光雷达具有589nm和532nm两个通道,将589nm通道用于计算,532nm通道作为参考.在150~250K区间等间隔选取多个温度初值,利用589nm通道反演60~70km高度范围内的大气温度廓线和大气密度廓线,利用参考密度廓线选取准确的温度初值,反演得到准确的大气温度廓线.将589nm通道和532nm通道反演的温度廓线相比较,发现二者具有较高的一致性.改进后的方法有效利用了信噪比较差的光信号,使589nm通道的温度探测上限从60km提高到70km.      

Ground-based measurement of strato–mesospheric CO by a FTIR spectrometer over Poker Flat,Alaska

Upper atmospheric CO above 24 km has been observed over Poker Flat (147°W, 65°N, altitude 0.61 km), Alaska using ground-based solar absorption infrared spectroscopy. This is the first reported detection of stratospheric–mesospheric CO using this method from the ground. The results clearly indicate that there is a seasonal variation of the CO profile with enhanced abundances in spring while remaining low from May onwards.The Poker Flat Research Range is one of the many measurement sites that constitute the Network for the Detection of Stratospheric Change (NDSC). The method used in this work, estimating the CO partial column abundances above the middle stratosphere, can be applied to spectra observed using FTIR spectrometry at many other NDSC sites. This suggests the availability of this established technique as a new method for CO measurements in the upper atmosphere.     

A future Japanese program in far-ultraviolet astronomy and solar physics

A design study of an ultraviolet-telescope satellite (UVSAT) has been carried out by the Japanese astronomical community. The main purpose of this satellite mission would be to investigate (i) the distribution and nature of ultraviolet sources in star clusters, galaxies, and clusters of galaxies, (ii) the physical structures of galactic nebulae, and (iii) the dynamics of stellar and galactic activity. A 60 cm, f/4 Cassegrain telescope will be launched into a semi-circular orbit of ～500 km altitude and of ～30° inclination, by an ISAS/Japan M3S-III rocket. An intensified CCD camera and/or a concave-grating spectrograph will be operated for the wavelength range λ?1200 Å. Various possible auxiliary instruments and sub-telescopes are considered.     

The Zadko telescope: A resource for science education enrichment

This paper describes plans for employing a new fully robotic optical telescope in Western Australia, the Zadko telescope, for evaluating education-outreach. A key feature of the telescope is the optimized observation-scheduling program, developed by our French collaborators who operate the TAROT robotic telescope network. It provides a simple interface for requesting observation time remotely, and has the potential for school students to participate in real astronomical research. The University of Western Australia and Curtin University are commencing a study for evaluating changes in student perceptions of science by participation in our astronomy research. Other areas of interest include broadening the program to include access of students from European countries, and exploring how remote access astronomy can be used to encourage awareness between different cultures.     

Structure of the Martian atmosphere from ε-Gem occultation observations

We have collected information about Martian atmospheric scale heights derived from observations of the occultation of ? Gem by Mars on April 8, 1976. The observations give data in the altitude range ～50 to 80 km. A rough, unweighted average of results so far available yields a temperature of ～165°K. Excursions of ～±40°K about this mean may be present as a function of both altitude and areographic coordinates.     

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.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.     

Long-period meteor radar temperature variations over Collm (51°N, 13°E) and Kazan (56°N, 49°E)

We have estimated temperatures from meteor radar measurements using the gradient method and the full width at half maximum method over Kazan (56°N 49°E) and Collm (51°N, 13°E). The time series cover the period 2016–2019. The temperature gradient model is constructed from SABER temperature observations. We demonstrate that annual mean, amplitudes and phases of the annual and semiannual oscillations of the radar temperatures are close to those of the MLS and SABER temperatures. The annual mean temperatures over Kazan and Collm differ non-significantly. The seasonal variability of the radar temperature is mostly due to the annual cycle which tends to grow with latitude. The gradient method produces temperatures which agree with the SABER temperatures better than with the MLS ones. The harmonics of the annual oscillations from periods of 73 days up to periods of about 40 days are the most significant day-to-day temperature oscillations and have zonal wavenumber zero. Their periods and phases are in good correspondence with those of the MLS and SABER ones. We also show some results which demonstrate that at 56°N the FWHM method is not as robust as the gradient method.     

Deorbiting spacecraft with passively stabilised attitude using a simplified quasi-rhombic-pyramid sail

This paper studies deorbiting using an analogue to the quasi-rhombic-pyramid concept for planar motion. The focus is on maintaining a stable (meaning oscillatory) attitude close to the direction of the velocity of the spacecraft relative to the atmosphere. The study consists of a massive computation of deorbit times chosen in a region of the phase space where atmospheric drag plays a leading role. Here, no damping effects are considered. Thus, any passive stabilisation observed is either due to solar radiation pressure or atmospheric drag. The results show that such stable deorbiting is feasible up to a threshold that depends upon the physical parameters of the sail. This threshold is around 500 km of altitude. Stable deorbiting is also shown to reduce the unpredictability that appears due to tumbling.     

Seasonal variations in vertical distribution of meteor decay time as observed from meteor radars at 8.5°N and 80°N

The decay times of meteor radar echoes have been used for decades to investigate characteristics of the mesosphere and lower thermosphere (MLT) region. As the meteor echo decay time depends on background atmospheric parameters, in the present communication, we examine the seasonal variation of the vertical distributions of underdense meteor echo decay times with respect to echo strength. Observations from two similar radars located at two distinct geographical locations, Thumba (8.5°N, 77°E) and Eureka (80°N, 85.8°W) were used for the present study. Here, the radar received signal power is categorized into strong and weak echoes and vertical profiles of their decay times are constructed. It has been noticed that the monthly mean decay time vertical profile turning altitude (i.e., inflection point) varies in the range of 80–87?km of altitude depending on latitude. The turning altitude is observed at relatively lower heights in the winter than in summer at both the latitudes. The present analysis shows that the meteor decay time below the mean turning altitude follows a decreasing trend with decreasing altitude, which is quite distinct to the behaviour of ambipolar diffusion. It is also observed that there is a difference in mean decay time of strong and weak echoes below 90?km of altitude, which is very prominently seen at lower altitudes. This difference shows a seasonal pattern at high latitude, but does not show any seasonal variation at low latitude. The present results are discussed in light of current understanding of the meteor decay time.     

Observation of aphelion cloud belt over Martian tropics,its evolution,and associated dust distribution from MCS data

The present study uses five Martian years of observations from Mars Climate Sounder onboard Mars Reconnaissance Orbiter for investigating the Aphelion Cloud Belt (ACB) over the tropics. Analysis of zonal mean water ice column opacity suggests that the spatial extension of the ACB is mainly confined over the tropics and mid-latitudes (-20 – 40°N) during L_S ~ 45 – 135° (L_S = 0° signifies northern spring equinox). The ACB is seen primarily in the nighttime only due to the truncation of the daytime profile observations at significantly higher altitudes (at ~30 km). Zonal mean ice extinction profiles show ACB’s altitudinal range within ~10 – 40 km, and the existence of a thin cloud band in the absence of a thick ACB during aphelion season. Three phases of the ACB could be identified as the formation phase during L_S = 45 – 75° (phase 1), the peak phase during L_S = 76 – 105° (phase 2), and the decaying phase during L_S = 106 – 135° (phase 3). Observation of the cloud latitude belt shows a northward movement starting from phase 2, prominent over regions nearby Lunae Planum and Xanthe Terra. During this phase, the top level of thick clouds within the ACB decreases to ~20 km in the southern hemisphere, while it increases a little over the northern hemisphere (NH). The decreasing tendency continues in phase 3 over the entire region ?10 – 10°N, and the thick cloud base moves higher over the NH, though the vertical depth of it becomes narrower than phase 2. Temperature profiles do not show any noticeable influence on the northward evolution of the ACB. However, the study at a regional level indicates a possible association of upper tropospheric dustiness with the ACB’s evolution. The mechanism is evident in the correlation analysis mostly at an altitude range of ~18 – 35 km. The migrating semidiurnal tide (SMD) as a proxy of dust or water ice forcing, and the calculated upper tropospheric dust radiative heating, shows an apparent northward movement of their peak amplitude within the three phases of the ACB. This match between the spatiotemporal variations of the SMD and the water ice was not observed previously. However, the correlating behavior seems to be prominent in the areas nearby Lunae Planum and Xanthe Terra and the upper-tropospheric region of the atmosphere.     

海口上空中高层大气强烈的突发钠层观测

海南激光雷达站（20°N,110°E）2011年11月2日观测到强烈的突发钠层（SSL）.该SSL峰值密度达37087cm-3,半峰全宽仅为0.9km.对该SSL发生前后的钠层峰值密度和高度特性进行分析.统计激光雷达2010年5月至2013年12月共计377天观测到的SSL事件222次,其中仅有1次SSL峰值密度超过30000cm-3.对比距离最近的海南儋州（19.5°N,109.1°E）测高仪和VHF雷达观测到的突发E层（Es）事件,分析事件的相关性.Es最低高度与SSL峰值高度差均在5km以内,约有75%的Es与SSL时间差在±30min之内.实验中使用continuum的Nd:YAG激光器和泵浦染料激光器产生589nm激光,能量为45mJ,使用直径为1000mm的望远镜接收钠层的光子回波.      

A balloon-borne sun-tracking multichannel photometer for atmospheric aerosol measurements

A balloon borne multichannel photometer for measurement of atmospheric scattering in the near ultraviolet and the visible wavelength regions has been developed at the Physical Research Laboratory, Ahmedabad for study of the size distribution and number density of aerosols at tropospheric and lower stratospheric altitudes. The instrumentation involves tracking the sun in elevation and scanning in azimuth. The payload was recently flown on a 100 kg. balloon from the Hyderabad Balloon Facility on 18 April 1984. The balloon reached a float altitude of 35 km and good quality data has been obtained from an altitude of 6 km upto float altitude. Data analysis is still in progress. The present paper details the instrument design and presents a few illustrations of the instrument performance from this flight.     

Scientific ballooning in Japan – An over view of recent activities

Current status of scientific ballooning in Japan is reviewed. First, I describe successful application of balloon technologies to construct a vessel of transparent plastic film, to contain about 1000 tons of liquid scintillator in Kamioka Liquid Scintillator Anti-Neutrino Detector (KamLAND). KamLAND is a project to study neutrino oscillation phenomena, by detecting anti-neutrinos produced in distant nuclear reactors. Next, I describe high altitude balloons developed by the ISAS balloon group. They developed balloons made from ultra-thin polyethylene film, producing a balloon of volume 60,000 m³ which successfully reached an altitude of 53 km in 2002. This is a world record, the greatest altitude that a balloon has ever achieved. ISAS is applying further effort to develop balloons with even thinner films, to achieve a higher altitude than 53 km. Other recent activities by the ISAS balloon group are briefly described.I also review scientific ballooning projects now operating in Japan, particularly focusing on the Balloon-Borne Experiment with a Superconducting Spectrometer (BESS) program. This is a US–Japan collaborative program that has carried out very precise measurements of antiprotons, protons and other components in primary cosmic rays, as well as measuring the fluxes of atmospheric muons and other components. The results of these observations give us important information to improve our understanding of the production mechanism of antiprotons observed in the primary cosmic rays. The data are also important for analysis of atmospheric neutrino events observed by Super-Kamiokande and other ground-based neutrino detectors. Future prospects of BESS and other balloon-borne cosmic-ray research programs are also presented.     

Atmosphere sounding by GPS radio occultation: First results from TanDEM-X and comparison with TerraSAR-X

On 21 June 2010 the TerraSAR-X satellite was joined by the TanDEM-X satellite. A Global Positioning System (GPS) radio occultation (RO) experiment using the twin satellites has been carried out to estimate the precision of GPS atmospheric soundings. For the Day Of Year (DOY) 330–336, 2011, we analyze phase and amplitude data recorded by GPS receivers separated by a few hundred meters in a low earth orbit and derive collocated atmospheric refractivity profiles. In the altitude range 10–20 km the standard deviation between TerraSAR-X and TanDEM-X refractivity does not exceed 0.15%. The standard deviation is rapidly increasing for lower and higher altitudes; close to the surface and at an altitude of 30 km the standard deviation reaches 0.8% and 0.5%, respectively. Systematic deviations between TerraSAR-X and TanDEM-X refractivity in the considered altitude range (0–30 km) are negligible. The results confirm the anticipated high precision of the GPS RO technique. However, the difference in the retrieved refractivity in the lower troposphere for different Open Loop (OL) signal tracking parameters, altered onboard TanDEM-X for DOY 49–55, 2012, calls for an in depth analysis. At the moment we can not exclude that a potential bias in the OL Doppler model introduces a bias in our retrieved refractivity at altitudes <8$<8$ km.     

Satellite availability and positioning performance of uncombined precise point positioning using BeiDou-2 and BeiDou-3 multi-frequency signals

To realize the smooth transition from regional BeiDou Navigation Satellite System (BDS-2) to the global one (BDS-3), the integration of BDS-2 and BDS-3 is important for providing continuous, stable and reliable positioning, navigation and timing (PNT) services for global users. This work used 154 globally distributed multi-GNSS (Global Navigation Satellite System) experiment stations spanning 30 days to analyze the satellite availability and positioning performance of uncombined precise point positioning (UC-PPP) under current BDS-2 and BDS-3 constellations. We focused on three issues: the influence of BDS-3 receiver tracking abilities, the positioning performance among different areas, and the benefit of multi-frequency observations. The results show that the elliptical zone caused by poor BDS-2 satellite visibility is disappeared when the evenly distributed BDS-3 medium earth orbit satellites are introduced. When BDS-3 are integrated with BDS-2, the area with the Position Dilution of Precision (PDOP) less than 2 can be expanded to 75° S-75° N and 30° E-150° W. The positioning performance of BDS-3 and BDS-2/BDS-3 UC-PPP are seriously affected by the receiver tracking abilities of BDS-3 signals. When the maximum pseudo-random noise sequences (PRNs) of BDS-3 satellites tracked by stations are within 30 or 37, the positioning accuracy of static UC-PPP can be improved by 22.94% or 8.27% due to the integration of BDS-2 and BDS-3. Besides, the most improvement of BDS-2 and BDS-3 integration is achieved in Asia-Pacific regions, especially for the kinematic UC-PPP or the poor receiver tracking abilities of BDS-3. Similar to the multi-frequency BDS-2 UC-PPP, the benefit of multi-frequency signals for BDS-3 or BDS-2/BDS-3 UC-PPP is also non-vital. The three-dimensional positioning accuracy of BDS-2/BDS-3 multi-frequency UC-PPP in static mode and kinematic mode are 2.24 cm and 5.39 cm, while the corresponding convergence time are 49.62 min and 73.80 min, respectively. Compared with BDS-2, both the positioning accuracy and the convergence time of BDS-2/BDS-3 joint UC-PPP are improved by approximately over 50%, which indicates that BDS-3 has a great potential to provide high-quality PNT services as other global navigation satellite systems.     

World grid of cosmic ray vertical cut-off rigidity for the last decade

Cosmic ray cut-off rigidity tables and maps over the world concerning the epochs 2010, 2015 and the current one 2020 have been constructed. These maps display the effective cut-off rigidity in every five degrees in latitude and in longitude at the altitude of 20 km above the surface of the international reference ellipsoid. The values of the geomagnetic cut-off rigidity were calculated in every 5° in latitude and in every 15° in longitude applying the well-known method of particle trajectory calculations resulted from the theory of the particle motion in the Earth's magnetic field. The applied software employed the 12th Generation of the International Geomagnetic Reference Field (IGRF 12) and trajectories were calculated at 0.01 GV intervals in order to determine the vertical cut-off rigidity for each location. Beyond the use of the calculated cut-off rigidity values as a basic reference of charged particle access to different geographical locations during quiet and/or more intense geomagnetic periods, these results can be used for a long- term forecasting of the geomagnetic conditions variations.