The new CMOS Tracking Camera used at the Zimmerwald Observatory

During the last years the use of tracking cameras for SLR observations became less important due to the high accuracy of the predicted orbits. Upcoming new targets like satellites in eccentric orbits and space debris objects, however, require tracking cameras again. In 2013 the interline CCD camera was replaced at the Zimmerwald Observatory with a so called scientific CMOS camera. This technology promises a better performance for this application than all kinds of CCD cameras. After the comparison of the different technologies the focus will be on the integration in the Zimmerwald SLR system.

Contribution of Starlette, Stella, and AJISAI to the SLR-derived global reference frame

The contribution of Starlette, Stella, and AJI-SAI is currently neglected when defining the International Terrestrial Reference Frame, despite a long time series of precise SLR observations and a huge amount of available data. The inferior accuracy of the orbits of low orbiting geodetic satellites is the main reason for this neglect. The Analysis Centers of the International Laser Ranging Service (ILRS ACs) do, however, consider including low orbiting geodetic satellites for deriving the standard ILRS products based on LAGEOS and Etalon satellites, instead of the sparsely observed, and thus, virtually negligible Etalons. We process ten years of SLR observations to Starlette, Stella, AJISAI, and LAGEOS and we assess the impact of these Low Earth Orbiting (LEO) SLR satellites on the SLR-derived parameters. We study different orbit parameterizations, in particular different arc lengths and the impact of pseudo-stochastic pulses and dynamical orbit parameters on the quality of the solutions. We found that the repeatability of the East and North components of station coordinates, the quality of polar coordinates, and the scale estimates of the reference are improved when combining LAGEOS with low orbiting SLR satellites. In the multi-SLR solutions, the scale and the Z component of geocenter coordinates are less affected by deficiencies in solar radiation pressure modeling than in the LAGEOS-1/2 solutions, due to substantially reduced correlations between the Z geocenter coordinate and empirical orbit parameters. Eventually, we found that the standard values of Center-of-mass corrections (CoM) for geodetic LEO satellites are not valid for the currently operating SLR systems. The variations of station-dependent differential range biases reach 52 and 25 mm for AJISAI and Starlette/Stella, respectively, which is why estimating station dependent range biases or using station-dependent CoM, instead of one value for all SLR stations, is strongly recommended.This clearly indicates that the ILRS effort to produce CoM corrections for each satellite, which are site-specific and depend on the system characteristics at the time of tracking,is very important and needs to be implemented in the SLR data analysis.

Impact of Earth radiation pressure on LAGEOS orbits and on the global scale

The indirect solar radiation pressure caused by reflected or re-emitted radiation by the Earth’s surface is an important non-gravitational force perturbing the orbits of geodetic satellites (Rubincam and Weiss, 1986; Martin and Rubincam, 1996). In the case of LAGEOS this acceleration is of the order of 15% of the direct solar radiation pressure. Therefore, Earth radiation pressure has a non-negligible impact not only on LAGEOS orbits, but also on the SLR-derived terrestrial reference frame. We investigate the impact of the Earth radiation pressure on LAGEOS orbits and on the SLR-derived parameters. Earth radiation pressure has a remarkable impact on the semi-major axes of the LAGEOS satellites, causing a systematic reduction of 1.5 mm. The infrared Earth radiation causes a reduction of about 1.0 mm and the Earth’s reflectivity of 0.5 mm of the LAGEOS’ semi-major axes. The global scale defined by the SLR network is changed by 0.07 ppb, when applying Earth radiation pressure. The resulting station heights differ by 0.5-0.6 mm in the solution with and without Earth radiation pressure. However, when range biases are estimated, the height differences are absorbed by the range biases, and thus, the station heights are not shifted.

A system design for space-based space surveillance

This paper presents the capabilities of a Space-Based Space Surveillance (SBSS) demonstration mission for Space Surveillance and Tracking (SST) based on a micro- satellite platform. The results have been produced in the frame of ESA’s "As sessment Study for Space Based Space Surveillance Demonstration Mission (Phase A) " performed by the Airbus DS consortium. Space Surveillance and Tracking is part of Space Situational Awareness (SSA) and covers the detection, tracking and cataloguing of spa ce debris and satellites. Derived SST services comprise a catalogue of these man-made objects, collision warning, detection and characterisation of in-orbit fragmentations, sub-catalogue debris characterisation, etc. The assessment of SBSS in an SST system architecture has shown that both an operational SBSS and also already a well - designed space-based demonstrator can provide substantial performance in terms of surveillance and tracking of beyond - LEO objects. Especially the early deployment of a demonstrator, possible by using standard equipment, could boost initial operating capability and create a self-maintained object catalogue. Unlike classical technology demonstration missions, the primary goal is the demonstration and optimisation of the functional elements in a complex end-to-end chain (mission planning, observation strategies, data acquisition, processing and fusion, etc.) until the final products can be offered to the users. The presented SBSS system concept takes the ESA SST System Requirements (derived within the ESA SSA Preparatory Program) into account and aims at fulfilling some of the SST core requirements in a stand-alone manner. The evaluation of the concept has shown that an according solution can be implemented with low technological effort and risk. The paper presents details of the system concept, candidate micro - satellite platforms, the observation strategy and the results of performance simulations for GEO coverage and cataloguing accuracy

Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 % ; in the absence of the polar data gap, SLR residuals agree by ≈96 % . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.

A comment on the article “A collinearity diagnosis of the GNSS geocenter determination” by P. Rebischung, Z. Altamimi, and T. Springer

Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) showed that the geocenter z -component estimated from observations of global navigation satellite systems (GNSS) is strongly correlated to a particular parameter of the solar radiation pressure (SRP) model developed by Beutler et al. (Manuscr Geod 19:367–386, 1994). They analyzed the forces caused by SRP and the impact on the satellites’ orbits. The authors achieved their results using perturbation theory and celestial mechanics. Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026, 2013) also deal with the geocenter determination with GNSS. The authors carried out a collinearity diagnosis of the associated parameter estimation problem. They conclude “without much exaggerating that current GNSS are insensitive to any component of geocenter motion”. They explain this inability by the high degree of collinearity of the geocenter coordinates mainly with satellite clock corrections. Based on these results and additional experiments, they state that the conclusions drawn by Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) are questionable. We do not agree with these conclusions and present our arguments in this article. In the first part, we review and highlight the main characteristics of the studies performed by Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) to show that the experiments are quite different from those performed by Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026,2013) . In the second part, we show that normal equation (NEQ) systems are regular when estimating geocenter coordinates, implying that the covariance matrices associated with the NEQ systems may be used to assess the sensitivity to geocenter coordinates in a standard way. The sensitivity of GNSS to the components of the geocenter is discussed. Finally, we comment on the arguments raised by Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026, 2013) against the results of Meindl et al. (Adv Space Res 51(7):1047–1064, 2013).

Earth Rotation and Gravity Field Parameters from Satellite Laser Ranging

We present the results from a simultaneous estimation of the gravity field, Earth rotation parameters, and station coordinates from combined SLR solutions incorporating up to nine geodetic satellites: LAGEOS-1/2, Starlette, Stella, AJISAI, Beacon-C, Lares, Blits and LARES. These solutions cover all three pillars of satellite geodesy and ensure full consistency between the Earth rotation parameters, gravity field coefficients, and geometry-related parameters. We address benefits emerging from such an approach and discuss particular aspects and limitations of the gravity field recovery using SLR data. The current accuracy of SLR-derived polar motion, by the means of WRMS w.r.t. IERS-08-C04 series, is at a level of 118-149 μas, which corresponds to 4 to 5 mm on the Earth’s surface. The WRMS of SLR-derived Length-of-Day, when the gravity field parameters are simultaneously estimated, is 56 μs/day, corresponding to about 26 mm on the ground, and the mean bias of SLR-derived Length-of-Day is 6.3 μs/day, corresponding to 3 mm.

GOCE precise orbit determination for the entire Mission – challenges in the final mission phase

The Gravity field and steady-state Ocean Circulation Explorer (GOCE), ESA’s first Earth Explorer core mission, was launched on March 17, 2009 into a sunsynchronous dusk-dawn orbit and eventually re-entered into the Earth’s atmosphere on November 11, 2013. A precise science orbit (PSO) product was provided by the GOCE High-level Processing Facility (HPF) from the GPS high-low Satellite-to-Satellite Tracking (hl-SST) data from the beginning until the very last days of the mission. We recapitulate the PSO procedure and refer to the results achieved until the official end of the GOCE mission on October 21, 2013, where independent validations with Satellite Laser ranging (SLR) measurements confirmed a high quality of the PSO product of about 2 cm 1-D RMS. We then focus on the period after the official end of the mission, where orbits could still be determined thanks to the continuously running GPS receivers delivering high quality data until a few hours before the re-entry into the Earth’s atmosphere. We address the challenges encountered for orbit determination during these last days and report on adaptions in the PSO procedure to also obtain good orbit results at the unprecedented low orbital altitudes below 224 km.

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

The impact of common versus separate estimation of orbit parameters on GRACE gravity field solutions

Gravity field parameters are usually determined from observations of the GRACE satellite mission together with arc-specific parameters in a generalized orbit determination process. When separating the estimation of gravity field parameters from the determination of the satellites’ orbits, correlations between orbit parameters and gravity field coefficients are ignored and the latter parameters are biased towards the a priori force model. We are thus confronted with a kind of hidden regularization. To decipher the underlying mechanisms, the Celestial Mechanics Approach is complemented by tools to modify the impact of the pseudo-stochastic arc-specific parameters on the normal equations level and to efficiently generate ensembles of solutions. By introducing a time variable a priori model and solving for hourly pseudo-stochastic accelerations, a significant reduction of noisy striping in the monthly solutions can be achieved. Setting up more frequent pseudo-stochastic parameters results in a further reduction of the noise, but also in a notable damping of the observed geophysical signals. To quantify the effect of the a priori model on the monthly solutions, the process of fixing the orbit parameters is replaced by an equivalent introduction of special pseudo-observations, i.e., by explicit regularization. The contribution of the thereby introduced a priori information is determined by a contribution analysis. The presented mechanism is valid universally. It may be used to separate any subset of parameters by pseudo-observations of a special design and to quantify the damage imposed on the solution.

Satellite laser ranging to GPS and GLONASS

Satellite laser ranging (SLR) to the satellites of the global navigation satellite systems (GNSS) provides substantial and valuable information about the accuracy and quality of GNSS orbits and allows for the SLR-GNSS co-location in space. In the framework of the NAVSTAR-SLR experiment two GPS satellites of Block-IIA were equipped with laser retroreflector arrays (LRAs), whereas all satellites of the GLONASS system are equipped with LRAs in an operational mode. We summarize the outcome of the NAVSTAR-SLR experiment by processing 20 years of SLR observations to GPS and 12 years of SLR observations to GLONASS satellites using the reprocessed microwave orbits provided by the center for orbit determination in Europe (CODE). The dependency of the SLR residuals on the size, shape, and number of corner cubes in LRAs is studied. We show that the mean SLR residuals and the RMS of residuals depend on the coating of the LRAs and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including the single-photon and multi-photon detection modes. We also show that the SLR observations to GNSS satellites are important to validate GNSS orbits and to assess deficiencies in the solar radiation pressure models. We found that the satellite signature effect, which is defined as a spread of optical pulse signals due to reflection from multiple reflectors, causes the variations of mean SLR residuals of up to 15 mm between the observations at nadir angles of 0∘ and 14∘. in case of multi-photon SLR stations. For single-photon SLR stations this effect does not exceed 1 mm. When using the new empirical CODE orbit model (ECOM), the SLR mean residual falls into the range 0.1–1.8 mm for high-performing single-photon SLR stations observing GLONASS-M satellites with uncoated corner cubes. For best-performing multi-photon stations the mean SLR residuals are between −12.2 and −25.6 mm due to the satellite signature effect.

SMARTnet: A Sensor for Monitoring the Geostationary Ring

As the number of space debris is increasing in the geostationary ring, it becomes mandatory for any satellite operator to avoid any collisions. Space debris in geosynchronous orbits may be observed with optical telescopes. Other than radar, that requires very large dishes and transmission powers for sensing high-altitude objects, optical observations do not depend on active illumination from ground and may be performed with notably smaller apertures. The detection size of an object depends on the aperture of the telescope, sky background and exposure time. With a telescope of 50 cm aperture, objects down to approximately 50 cm may be observed. This size is regarded as a threshold for the identification of hazardous objects and the prevention of potentially catastrophic collisions in geostationary orbits. In collaboration with the Astronomical Institute of the University of Bern (AIUB), the German Space Operations Center (GSOC) is building a small aperture telescope to demonstrate the feasibility of optical surveillance of the geostationary ring. The telescope will be located in the southern hemisphere and complement an existing telescope in the northern hemisphere already operated by AIUB. These two telescopes provide an optimum coverage of European GEO satellites and enable a continuous monitoring independent of seasonal limitations. The telescope will be operated completely automatically. The automated operations should be demonstrated covering the full range of activities including scheduling of observations, telescope and camera control as well as data processing.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

SMARTnet: First Experience of Setting Up a Telescope System to Survey the Geostationary Ring

Space debris in geostationary orbits may be detected with optical telescopes when the objects are illuminated by the Sun. The advantage compared to Radar can be found in the illumination: radar illuminates the objects and thus the detection sensitivity depletest proportional to the fourth power of the d istance. The German Space Operation Center, GSOC, together with the Astronomical Institute of the University of Bern, AIUB, are setting up a telescope system called SMARTnet to demonstrate the capability of performing geostationary surveillance. Such a telescope system will consist of two telescopes on one mount: a smaller telescope with an aperture of 20cm will serve for fast survey while the larger one, a telescope with an aperture of 50cm, will be used for follow-up observations. The telescopes will be operated by GSOC from Oberpfaffenhofen by the internal monitoring and control system called SMARTnetMAC. The observation plan will be generated by MARTnetPlanning seven days in advance by applying an optimized planning scheduler, taking into account fault time like cloudy nights, priority of objects etc. From each picture taken, stars will be identified and everything not being a star is treated as a possible object. If the same object can be identified on multiple pictures within a short time span, the trace is called a tracklet. In the next step, several tracklets will be correlated to identify individual objects, ephemeris data for these objects are generated and catalogued . This will allow for services like collision avoidance to ensure safe operations for GSOC’s satellites. The complete data processing chain is handled by BACARDI, the backbone catalogue of relational debris information and is presented as a poster.