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.

From planetesimals to planets: volatile molecules

Context. Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. Aims. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. Methods. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. Results. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01 ± 0.33 and 0.8 ± 0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2–3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

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.

Gas composition of the main volatile elements in protoplanetary discs and its implication for planet formation

Context. Direct observations of gaseous exoplanets reveal that their gas envelope has a higher C/O ratio than that of the host star (e.g., Wasp 12-b). This has been explained by considering that the gas phase of the disc could be inhomogeneous, exceeding the stellar C/O ratio in regions where these planets formed; but few studies have considered the drift of the gas and planet migration. Aims. We aim to derive the gas composition in planets through planet formation to evaluate if the formation of giant planets with an enriched C/O ratio is possible. The study focusses on the effects of different processes on the C/O ratio, such as the disc evolution, the drift of gas, and planet migration. Methods. We used our previous models for computing the chemical composition, together with a planet formation model, to which we added the composition and drift of the gas phase of the disc, which is composed of the main volatile species H2O, CO, CO2, NH3, N2, CH3OH, CH4, and H2S, H2 and He. The study focusses on the region where ice lines are present and influence the C/O ratio of the planets. Results. Modelling shows that the condensation of volatile species as a function of radial distance allows for C/O enrichment in specific parts of the protoplanetary disc of up to four times the solar value. This leads to the formation of planets that can be enriched in C/O in their envelope up to three times the solar value. Planet migration, gas phase evolution and disc irradiation enables the evolution of the initial C/O ratio that decreases in the outer part of the disc and increases in the inner part of the disc. The total C/O ratio of the planets is governed by the contribution of ices accreted, suggesting that high C/O ratios measured in planetary atmospheres are indicative of a lack of exchange of material between the core of a planet and its envelope or an observational bias. It also suggests that the observed C/O ratio is not representative of the total C/O ratio of the planet.

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.

The influence of absorbed solar radiation by Saharan dust on hurricane genesis

To date, the radiative impact of dust and the Sahar an air layer (SAL) on North Atlantic hurricane activity is not yet known. According to previous studies, dust stabilizes the atmosphere due to absorption of solar radiation but thus shifts convection to regions more conducive for hurricane genesis. Here we analyze differences in hurricane genesis and frequency from ensemble sensitivity simulations with radiatively active and inactive dust in the aerosol-climate model ECHAM6-HAM. We investigate dust burden and other hurricane-related variables and determine their influence on disturbances which develop into hurricanes (developing disturbances, DDs) and those which do not (nondeveloping disturbances, NDDs). Dust and the SAL are found to potentially have both inhibiting and supporting influences on background conditions for hurricane genesis. A slight southward shift of DDs is determined when dust is active as well as a significant warming of the SAL, which leads to a strengthening of the vertical circulation associated with the SAL. The dust burden of DDs is smaller in active dust simulations compared to DDs in simulations with inactive dust, while NDDs contain more dust in active dust simulations. However, no significant influence of radiatively active dust on other variables in DDs and NDDs is found. Furthermore, no substantial change in the DD and NDD frequency due to the radiative effects of dust can be detected.

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.

The role of seasonality of mineral dust concentration and size on glacial/interglacial dust changes in the EPICA Dronning Maud Land ice core

We present a record of particulate dust concentration and size distribution in subannual resolution measured on the European Project for Ice Coring in Antarctica (EPICA) Dronning Maud Land (EDML) ice core drilled in the Atlantic sector of the East Antarctic plateau. The record reaches from present day back to the penultimate glacial until 145,000 years B.P. with subannual resolution from 60,000 years B.P. to the present. Mean dust concentrations are a factor of 46 higher during the glacial (~850–4600 ng/mL) compared to the Holocene (~16–112 ng/mL) with slightly smaller dust particles during the glacial comparedto the Holocene and with an absolute minimum in the dust size at 16,000 years B.P. The changes in dust concentration are mainly attributed to changes in source conditions in southern South America. An increase in the modal value of the dust size suggests that at 16,000 years B.P. a major change in atmospheric circulation apparently allowed more direct transport of dust particles to the EDML drill site. We find a clear in-phase relation of the seasonal variation in dust mass concentration and dust size during the glacial (r(conc,size) = 0.8) but no clear phase relationship during the Holocene (0

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

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.