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.     

Goddard Robotic Telescope – Optical follow-up of GRBs and coordinated observations of AGNs

Since it is not possible to predict when a Gamma-Ray Burst (GRB) will occur or when Active Galactic Nucleus (AGN) flaring activity starts, follow-up/monitoring ground telescopes must be located as uniformly as possible all over the world in order to collect data simultaneously with Fermi and Swift detections. However, there is a distinct gap in follow-up coverage of telescopes in the eastern U.S. region based on the operations of Swift. Motivated by this fact, we have constructed a 14″ fully automated optical robotic telescope, Goddard Robotic Telescope (GRT), at the Goddard Geophysical and Astronomical Observatory. The aims of our robotic telescope are (1) to follow-up Swift/Fermi GRBs and (2) to perform the coordinated optical observations of Fermi Large Area Telescope (LAT) AGN. Our telescope system consists of off-the-shelf hardware. With the focal reducer, we are able to match the field of view of Swift narrow instruments (20′ × 20′). We started scientific observations in mid-November 2008 and GRT has been fully remotely operated since August 2009. The 3σ upper limit in a 30 s exposure in the R filter is ∼15.4 mag; however, we can reach to ∼18 mag in a 600 s exposures. Due to the weather condition at the telescope site, our observing efficiency is 30–40% on average.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.     

A new algorithm for optical observations of space debris with the TAROT telescopes

Since 2004, we observe satellites in the geostationary orbit with a network of robotic ground based fully automated telescopes called TAROT. One of them is located in France and the second at ESO, La Silla, Chile. The system processes the data in real time. Its wide field of view is useful for the discovery, the systematic survey and for the tracking of both catalogued and un-catalogued objects. We present a new source extraction algorithm based on morphological mathematic, which has been tested and is currently under implementation in the standard pipeline. Using this method, the observation strategy will correlate the measurements of the same object on successive images and give better detection rate and false alarm rate than the previous one. The overall efficiency and quality of the survey of the geostationary orbit has drastically improved and we can now detect satellites and debris in different orbits like Geostationary Transfer Orbit (GTO). Results obtained in real conditions with TAROT are presented.     

First astrometric observations of space debris with the MéO telescope

The MéO (for Métrologie Optique) telescope is the Satellite and Lunar Laser Ranging (SLR) dedicated telescope of Observatoire de la Côte d’Azur (France) located at plateau de Calern. The telescope uses an altazimuth mount. The motorization of the mount has a capability of 6 deg/s allowing the follow up of Low Earth Orbits (LEO) satellites, as well as Medium Earth Orbits (MEO) and geostationary (GEO) satellites, and the Moon. The telescope has a primary mirror of 1.54 m. It uses a Nasmyth focus equipped with an EMCCD camera. The telescope field of view, defined by the equivalent focal length and the size of the camera, is currently 3.4 arcmin × 3.4 arcmin.     

Status of MAGIC and recent results

Ground-based γ-ray astronomy is part of a new field of fundamental research of Astroparticle Physics, that recently made spectacular discoveries mostly thanks to Imaging Air Cherenkov Telescopes (IACT). The MAGIC telescope is a IACT located at La Palma, Canary Islands, Spain. Composed of two telescopes with 17 m diameter each, MAGIC is equipped with the largest optical reflectors in the world, and it has the lowest threshold energy (25 GeV). MAGIC started operations in 2004 in the single-detector configuration, and in 2009 as a stereo detector. Since then, it has discovered many new sources and classes of sources, both galactic and extragalactic. Here some highlights from the most recent results are presented.     

In orbit performance of butane propulsion system

The present work focuses on the determination of the in orbit performance of the Alsat-1 microsatellite propulsion system. The satellite mass is 90 kg, of which 6.2 kg is the propulsion system dry mass. The system is a butane propulsion system using low power resistojet thruster with 2.3 kg of propellant. The liquefied butane gas was selected due to its higher storage density and safety compared to the other propellants used for microsatellites. The purpose of this paper is the analysis of the firings performed after the launch of the satellite and to evaluate the system specific impulse and thrust level during the system lifetime. A total of 273 firings were performed on the Alsat-1 propulsion system in the period between the end of 2002 and mid 2009, the cumulated firing time is more than 12 h 49 min. The analysis of all the propulsion telemetry data shows that the system provides a total mission delta V of 25.3 m/s which is more than the 10 m/s specified for this mission. Furthermore, the mission average specific impulse and thrust are respectively 99.9 s and 48.8 mN.     

Shape and motion estimate of LEO debris using light curves

We have derived a tri-axial ellipsoidal model of an LEO object, a Cosmos 2082 rocket body, including its rotational axis direction, rotation period, precession, and a compositional parameter, using only light curve data from an optical telescope. The brightness of the object was monitored for two days and least-squares fitting was used to determine these values. The derived axial ratios are 100:18:18, the coordinates of the rotational axis direction on the celestial sphere are R.A. = 305.8° and Dec. = 2.6°, and its observed average rotation period is 41 s. When precession is considered, its amplitude and precession period are 30.5° and 29.4 min. These results show that optical light curve data are sufficient to determine the shapes and the motions of LEO objects.     

Synergetic retrieval of terrestrial AOD from MODIS images of twin satellites Terra and Aqua

Aerosol optical depth (AOD) is one of the most important indicators of atmospheric pollution. It can be retrieved from satellite imagery using several established methods, such as the dark dense vegetation method and the deep blue algorithm. All of these methods require estimation of surface reflectance prior to retrieval, and are applicable to a certain pre-designated type of surface cover. Such limitations can be overcome by using a synergetic method of retrieval proposed in this study. This innovative method is based on the fact that the ratio K of surface reflectance at different angles/geometries is independent of wavelength as reported by Flowerdew and Haigh (1995). An atmospheric radiative transfer model was then established and resolved with the assistance of the ratio K obtained from two Moderate Resolution Imaging Spectroradiometer (MODIS) spectral bands acquired from the twin satellites of Terra and Aqua whose overpass is separated by three hours. This synergetic method of retrieval was tested with 20 pairs of MODIS images. The retrieved AOD was validated against the ground observed AOD at the Taihu station of the AErosol RObotic NETwork (AERONET). It is found that they are correlated with the observations at a coefficient of 0.828 at 0.47 μm and 0.921 at 0.66 μm wavelengths. The retrieved AOD has a mean relative error of 25.47% at 0.47 μm and 24.3% at 0.66 μm. Of the 20 samples, 15 and 17 fall within two standard error of the line based observed AOD data on the ground at the 0.47 μm and 0.66 μm, respectively. These results indicate that this synergetic method can be used to reliably retrieve AOD from the twin satellites MODIS images, namely Terra and Aqua. It is not necessary to determine surface reflectance first.     

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.     

Properties of the high area-to-mass ratio space debris population at high altitudes

In the framework of its space debris research activities ESA established an optical survey program to study the space debris environment at high altitudes, in particular in the geostationary ring and in the geostationary transfer orbit region. The Astronomical Institute of the University of Bern (AIUB) performs these surveys on behalf of ESA using ESA’s 1-m telescope in Tenerife. Regular observations were started in 1999 and are continued during about 120–140 nights per year. Results from these surveys revealed a substantial amount of space debris at high altitudes in the size range from 0.1 to 1 m. Several space debris populations with different dynamical properties were identified in the geostationary ring. During the searches for debris in the geostationary transfer orbit region a new population of objects in unexpected orbits, where no potential progenitors exist, was found. The orbital periods of these objects are clustered around one revolution per day; the eccentricities, however, are scattered between 0 and 0.6. By following-up some of these objects using the ESA telescope and AIUB’s 1-m telescope in Zimmerwald, Switzerland, it was possible to study the properties of this new population. One spectacular finding from monitoring the orbits over time spans of days to months is the fact that these objects must have extreme area-to-mass ratios, which are by several orders of magnitudes higher than for ‘normal-type’ debris. This in turn supports the hypothesis that the new population actually is debris generated in or near the geostationary ring and which is in orbits with periodically varying eccentricity and inclination due to perturbations by solar radiation pressure. In order to further study the nature of these debris, multi-color and temporal photometry (light curves) were acquired with the Zimmerwald telescope. The light curves show strong variations over short time intervals, including signals typical for specular reflections. Some objects exhibit distinct periodic variations with periods ranging from 10 to several 100 s. All this is indicative for objects with complicated shapes and some highly reflective surfaces.     

The study of enhanced earth observations on a satellite image chain

The Multi-Spectral Camera (MSC) on the KOrea Multi-Propose SATellite (KOMPSAT)-2 was developed and launched as a main payload to provide a One(1) m panchromatic image and four(4) band four(4) m multi-spectral images at an altitude of 685 km covering a swath width of 15 km. These images, archived around the world, are a useful resource for space applications in agriculture, cartography, geology, forestry, regional planning, surveillance, and national security. The image quality of KOMPSAT-2 depends upon its image chain, which is comprised of an on-board system in the satellite and a processing system at the ground station. Therefore, in this study we determine the factors that have a major impact on the image quality through an investigation of the entire image chain. Consequently, two methods, involving a compression algorithm and a deconvolution technique, were determined as having a significant influence on the KOMPSAT-2 image quality. The compression algorithm of KOMPSAT-2 is rate-controlled JPEG-like algorithm that controls the mismatch between the input and output data rate. The ability to control the input/output data rate may be useful during the operation of the satellite but can also lower the overall image quality. The deconvolution technique may increase the sharpness of images, but it can also amplify the image noise level. Therefore, we propose methods of wavelet-based compression and denoising as an alternative to currently existing algorithms. Satisfactory results were obtained through experimentation with these two algorithms, and they are expected to be successfully implemented into the future KOMPSAT series to yield high-quality images for enhanced earth observation.     

A genetic algorithm for Initial Orbit Determination from a too short arc optical observation

Optical observations constitute a source of angular measurements of a satellite pass. Commonly, these observations have short durations with respect to the satellite orbit period. As a consequence, the use of classical orbit determination algorithms, as Laplace, Gauss or Escobal methods, give very poor results. The present work faces with the problem of estimating the orbital parameters of an orbiting object using its optical streak acquired by a telescope or a high accuracy camera. In the paper a new technique is developed for the Initial Orbit Determination from optical data by exploiting the genetic algorithms. The algorithm works without restrictions on the observer location. A recent challenging problem is the Initial Orbit Determination with space-based observations. This work focuses on the problem of determinating the orbital parameters of a satellite from an orbiting observer in LEO, using short time observations. We present the results based on a simulation with the observer on a sun-synchronous orbit with a single observation of just 60 s. Monte Carlo simulations are presented with different levels of sensor accuracy to show the reliability of the algorithm. The algorithm is able to yield a candidate solution for each observation. The coplanar case is analyzed and discussed as well.     

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

Low-latitude ionospheric scintillations and total electron content obtained with the CITRIS instrument on STPSat1 using radio transmissions from DORIS ground beacons

The primary objective of the Scintillation and Tomography Receiver in Space (CITRIS) is to detect ionospheric irregularities from space at low latitude. For this purpose, the satellite receiver uses the UHF and S-Band transmissions of the ground network of Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) beacons. CITRIS, developed at the Naval Research Laboratory, differs from the normal DORIS receiver by being able to capture and store the complex amplitude of the 401.25 and 2036.25 MHz transmissions at 200 Hz sample rate. Ground processing of the CITRIS data yields total electron content (TEC) and both phase and amplitude scintillations. With CITRIS flying on the US Space Test Program (STP) satellite STPSat1, 2 years of data were collected and processed to determine the fluctuations in ionospheric TEC and radio scintillations associated with equatorial irregularities. CITRIS flights over DORIS transmitters yield direct measurements of the horizontal plasma density fluctuations associated with equatorial plasma bubbles. Future flights of CITRIS can provide valuable complements to other satellite instruments such as GPS occultation receivers used to estimate vertical electron density profiles in the ionosphere.     

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.     

Development of simultaneous imaging polarimeter for asteroids

We are currently developing a polarimeter to study surface physical properties of asteroids. To enhance polarimetric accuracy and observational efficiency, we newly devised the polarimeter whose measurements can provide the two Stokes parameters Q and U, simultaneously. The test-observations of the prototype polarimeter have been carried out in December 2003 and January 2004, mounted on the 101 cm telescope at Bisei Astronomical Observatory, Okayama, Japan. In the observations, unpolarized and polarized standard stars were observed to measure an instrumental polarization and its uncertainty. As a result, an instrumental polarization of 5.06 ± 0.18% has been measured.     

Feasibility of performing space surveillance tasks with a proposed space-based optical architecture

Under ESA contract an industrial consortium including Aboa Space Research Oy (ASRO), the Astronomical Institute of the University of Bern (AIUB), and the Dutch National Aerospace Laboratory (NLR), proposed the observation concept, developed a suitable sensor architecture, and assessed the performance of a space-based optical (SBO) telescope in 2005. The goal of the SBO study was to analyse how the existing knowledge gap in the space debris population in the millimetre and centimetre regime may be closed by means of a passive optical instrument. The SBO instrument was requested to provide statistical information on the space debris population in terms of number of objects and size distribution. The SBO instrument was considered to be a cost-efficient with 20 cm aperture and 6° field-of-view and having flexible integration requirements. It should be possible to integrate the SBO instrument easily as a secondary payload on satellites launched into low-Earth orbits (LEO), or into geostationary orbit (GEO). Thus the selected mission concept only allowed for fix-mounted telescopes, and the pointing direction could be requested freely. Since 2007 ESA focuses space surveillance and tracking activities in the Space Situational Awareness (SSA) preparatory program. Ground-based radars and optical telescopes are studied for the build-up and maintenance of a catalogue of objects. In this paper we analyse how the proposed SBO architecture could contribute to the space surveillance tasks survey and tracking. We assume that the SBO instrumentation is placed into a circular sun-synchronous orbit at 800 km altitude. We discuss the observation conditions of objects at higher altitude, and select an orbit close to the terminator plane. A pointing of the sensor orthogonal to the orbital plane with optimal elevation slightly in positive direction (0° and +5°) is found optimal for accessing the entire GEO regime within one day, implying a very good coverage of controlled objects in GEO, too. Simulations using ESA’s Program for Radar and Optical Observation Forecasting (PROOF) in the version 2005 and a GEO reference population extracted from DISCOS revealed that the proposed pointing scenario provides low phase angles together with low angular velocities of the objects crossing the field-of-view. Radiometric simulations show that the optimal exposure time is 1–2 s, and that spherical objects in GEO with a diameter of below 1 m can be detected. The GEO population can be covered under proper illumination nearly completely, but seasonal drops of the coverage are possible. Subsequent observations of objects are on average at least every 1.5 days, not exceeding 3 days at maximum. A single observation arc spans 3° to 5° on average. Using a simulation environment that connects PROOF to AIUB’s program system CelMech we verify the consistency of the initial orbit determination for five selected test objects on subsequent days as a function of realistic astrometric noise levels. The initial orbit determination is possible. We define requirements for a correlator process essential for catalogue build-up and maintenance. Each single observation should provide an astrometric accuracy of at least 1”–1.5” so that the initially determined orbits are consistent within a few hundred kilometres for the semi-major axis, 0.01 for the eccentricity, and 0.1° for the inclination.     

Response of the Earth’s lower ionosphere to the Ground Level Enhancement event of December 13, 2006

In this study we analyze the Ground Level Enhancement Event No 70 observed on December 13, 2006, by correlating the observations from two research topics: Cosmic rays and Very Low Frequency (VLF < 30 kHz) wave propagation, as two ground based techniques for the detection of solar proton events, and their impact on the lower ionosphere. The observations have been endorsed from recordings of worldwide network ground based Neutron Monitors as well as by satellite data from the satellites GOES 12 (www.swpc.noaa.gov) and Pamela (www.pamela.roma2infn.it).