Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

On the evolution of the oceanic component of the IPSL climate models from CMIP3 to CMIP5: A mean state comparison

This study analyses the impact on the oceanic mean state of the evolution of the oceanic component (NEMO) of the climate model developed at Institut Pierre Simon Laplace (IPSL-CM), from the version IPSL-CM4, used for third phase of the Coupled Model Intercomparison Project (CMIP3), to IPSL-CM5A, used for CMIP5. Several modifications have been implemented between these two versions, in particular an interactive coupling with a biogeochemical module, a 3-band model for the penetration of the solar radiation, partial steps at the bottom of the ocean and a set of physical parameterisations to improve the representation of the impact of turbulent and tidal mixing. A set of forced and coupled experiments is used to single out the effect of each of these modifications and more generally the evolution of the oceanic component on the IPSL coupled models family. Major improvements are located in the Southern Ocean, where physical parameterisations such as partial steps and tidal mixing reinforce the barotropic transport of water mass, in particular in the Antarctic Circumpolar Current) and ensure a better representation of Antarctic bottom water masses. However, our analysis highlights that modifications, which substantially improve ocean dynamics in forced configuration, can yield or amplify biases in coupled configuration. In particular, the activation of radiative biophysical coupling between biogeochemical cycle and ocean dynamics results in a cooling of the ocean mean state. This illustrates the difficulty to improve and tune coupled climate models, given the large number of degrees of freedom and the potential compensating effects masking some biases.

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.

Relation between Excitation Power Density and Er3+ Doping Yielding the Highest Absolute Upconversion Quantum Yield

The upconversion quantum yield (UCQY) is one of the most significant parameters for upconverter materials. A high UCQY is essential for a succesful integration of upconversion in many applications, such as harvesting of the solar radiation. However, little is known about which doping level of the rare-earth ions yields the highest UCQY in the different host lattices and what are the underlying causes. Here, we investigate which Er3+ doping yields the highest UCQY in the host lattices β-NaYF4 and Gd2O2S under 4I15/2 → 4I13/2 excitation. We show for both host lattices that the optimum Er3+ doping is not fixed and it actually decreases as the irradiance of the excitation increases. To find the optimum Er3+ doping for a given irradiance, we determined the peak position of the internal UCQY as a function of the average Er−Er distance. For this purpose, we used a fit on experimental data, where the average Er−Er distance was calculated from the Er3+ doping of the upconverter samples and the lattice parameters of the host materials. We observe optimum average Er−Er distances for the host lattices β-NaYF4 and Gd2O2S with differences <14% at the same irradiance levels, whereas the optimum Er3+ doping are around 2× higher for β-NaYF4 than for Gd2O2S. Estimations by extrapolation to higher irradiances indicate that the optimum average Er−Er distance converges to values around 0.88 and 0.83 nm for β-NaYF4 and Gd2O2S, respectively. Our findings point to a fundamental relationship and focusing on the average distance between the active rare-earth ions might be a very efficient way to optimize the doping of rare-earth ions with regard to the highest achievable UCQY.

Vitamin D status among children and adolescents on anticonvulsant drugs in southern Switzerland.

INTRODUCTION It is recognised that vitamin D status is often inadequate (<50 nmol/l) in epileptic children, mainly because some anticonvulsant drugs induce the enzymes responsible for its metabolism. The purpose of the present study was to address vitamin D status among children and adolescents treated with anticonvulsant drugs and control subjects who reside in southern Switzerland, a high solar radiation region. METHODS Between January and May 2013, total serum 25-hydroxyvitamin D was assessed by liquid chromatography-tandem mass spectrometry in 58 children and adolescents with epilepsy and 29 controls residing in southern Switzerland. Dark-skinned individuals, females wearing dress styles covering practically the whole body and subjects with body mass index ≥85th percentile for age and sex were excluded. RESULTS Concentration of serum 25-hydroxyvitamin D was similar in epilepsy patients (48 [37-62] nmol/l; median and interquartile range) and controls (53 [47-64] nmol/l). An inadequate serum 25-hydroxyvitamin D concentration was common both among patients (55%) and control subjects (34%). Serum 25-hydroxyvitamin D was significantly lower among patients treated with anticonvulsant drugs that induce the metabolism of vitamin D (30 [21-51] nmol/l) than among the remaining patients (51 [40-65] nmol/l) and controls. CONCLUSIONS The present study indicates a relevant tendency towards inadequate vitamin D status among children with and without anticonvulsant drug management who reside in southern Switzerland. This tendency is more prominent in patients treated with anticonvulsant drugs that induce the metabolism of 25-hydroxyvitamin D.

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).

Galileo Orbit and Clock Quality of the IGS Multi-GNSS Experiment

The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) aims at the data collection and analysis of all available satellite navigation systems. In particular the new global and regional satellite navigation systems are of interest, i.e., the European Galileo, the Chinese BeiDou, the Japanese QZSS as well as satellite based augmentation systems. This article analyzes the orbit and clock quality of the Galileo products of four MGEX analysis centers for a common time period of 20 weeks. Orbit comparisons of the individual analysis centers have a consistency at the 5–30 cm level. Day boundary discontinuities range from 4 to 28 cm whereas 2-day orbit fit RMS values vary between 1 and 7 cm. The accuracy evaluated by satellite laser ranging residuals is on the one decimeter level with a systematic bias of about −5 cm for all analysis centers. In addition, systematic errors on the decimeter level related to solar radiation pressure mismodeling are present in all orbit products. Due to the correlation of radial orbit errors with the clock parameters, these errors are also visible as a bump in the Allan deviation of the Galileo satellite clocks at the orbital frequency.

Terrestrial Radar Interferometric Measurement of Hillslope Deformation and Atmospheric Disturbances in the Illgraben Debris-Flow Catchment, Switzerland

We describe a method for rapid identification and precise quantification of slope deformation using a portable radar interferometer. A rockslide with creep-like behavior was identified in the rugged and inaccessible headwaters of the Illgraben debris-flow catchment, located in the Central Swiss Alps. The estimated volume of the moving rock mass was approximately 0.5 x 10(6) m(3) with a maximum daily (3-D) displacement rate of 3 mm. Fast scene acquisition in the order of 6 s/scene led to uniquely precise mapping of spatial and temporal variability of atmospheric phase delay. Observations led to a simple qualitative model for prediction of atmospheric disturbances using a simple model for solar radiation, which can be used for advanced campaign planning for short observation periods (hours to days).

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.

Pasture degradation modifies the water and carbon cycles of the Tibetan highlands

The Tibetan Plateau has a significant role with regard to atmospheric circulation and the monsoon in particular. Changes between a closed plant cover and open bare soil are one of the striking effects of land use degradation observed with unsustainable range management or climate change, but experiments investigating changes of surface properties and processes together with atmospheric feedbacks are rare and have not been undertaken in the world's two largest alpine ecosystems, the alpine steppe and the Kobresia pygmaea pastures of the Tibetan Plateau. We connected measurements of micro-lysimeter, chamber, 13C labelling, and eddy covariance and combined the observations with land surface and atmospheric models, adapted to the highland conditions. This allowed us to analyse how three degradation stages affect the water and carbon cycle of pastures on the landscape scale within the core region of the Kobresia pygmaea ecosystem. The study revealed that increasing degradation of the Kobresia turf affects carbon allocation and strongly reduces the carbon uptake, compromising the function of Kobresia pastures as a carbon sink. Pasture degradation leads to a shift from transpiration to evaporation while a change in the sum of evapotranspiration over a longer period cannot be confirmed. The results show an earlier onset of convection and cloud generation, likely triggered by a shift in evapotranspiration timing when dominated by evaporation. Consequently, precipitation starts earlier and clouds decrease the incoming solar radiation. In summary, the changes in surface properties by pasture degradation found on the highland have a significant influence on larger scales.

Rapid responses of high-mountain vegetation to early Holocene environmental changes in the Swiss Alps

1 The Early Holocene sediment of a lake at tree line (Gouillé Rion, 2343 m a.s.l.) in the Swiss Central Alps was sampled for plant macrofossils. Thin (0.5 cm) slices, representing time intervals of c. 50 years each from 11 800 to 7800 cal. year bp, were analysed and the data compared with independent palaeoclimatic proxies to study vegetational responses to environmental change. 2 Alpine plant communities (e.g. with Salix herbacea) were established at 11 600–11 500 cal. year bp, when oxygen-isotope records showed that temperatures increased by c. 3–4 °C within decades. Larix decidua trees reached the site at c. 11 350 cal. year bp, probably in response to further warming by 1–2 °C. Forests dominated by L. decidua persisted until 9600 cal. year bp, when Pinus cembra became more important. 3 The dominance of Larix decidua for two millennia is explained by dry summer conditions, and possibly low winter temperatures, which favoured it over the late-successional Pinus cembra. Environmental conditions were a result of variations in the earth's orbit, leading to a maximum of summer and a minimum of winter solar radiation. Other heliophilous and drought-adapted species, such as Dryas octopetala and Juniperus nana, could persist in the open L. decidua forests, but were out-competed when the shade-tolerant P. cembra expanded. 4 The relative importance of Larix decidua decreased during periods of diminished solar radiation at 11 100, 10 100 and 9400 cal. year bp. Stable concentrations of L. decidua indicate that these percentage oscillations were caused by temporary increases of Pinus cembra, Dryas octopetala and Juniperus nana that can be explained by increases in moisture and/or decreases in summer temperature. 5 The final collapse of Larix decidua at 8400 cal. year bp was possibly related to abrupt climatic cooling as a consequence of a large meltwater input to the North Atlantic. Similarly, the temporary exclusion of Pinus cembra from tree line at 10 600–10 200 cal. year bp may be related to slowing down of thermohaline circulation at 10 700–10 300 cal. year bp. 6 Our results show that tree line vegetation was in dynamic equilibrium with climate, even during periods of extraordinarily rapid climatic change. They also imply that forecasted global warming may trigger rapid upslope movements of the tree line of up to 800 m within a few decades or centuries at most, probably inducing large-scale displacements of plant species as well as irrecoverable biodiversity losses.

Optical Light Curve Observations to Determine Attitude States of Space Debris

The currently proposed space debris remediation measures include the active removal of large objects and “just in time” collision avoidance by deviating the objects using, e.g., ground-based lasers. Both techniques require precise knowledge of the attitude state and state changes of the target objects. In the former case, to devise methods to grapple the target by a tug spacecraft, in the latter, to precisely propagate the orbits of potential collision partners as disturbing forces like air drag and solar radiation pressure depend on the attitude of the objects. Non-resolving optical observations of the magnitude variations, so-called light curves, are a promising technique to determine rotation or tumbling rates and the orientations of the actual rotation axis of objects, as well as their temporal changes. The 1-meter telescope ZIMLAT of the Astronomical Institute of the University of Bern has been used to collect light curves of MEO and GEO objects for a considerable period of time. Recently, light curves of Low Earth Orbit (LEO) targets were acquired as well. We present different observation methods, including active tracking using a CCD subframe readout technique, and the use of a high-speed scientific CMOS camera. Technical challenges when tracking objects with poor orbit redictions, as well as different data reduction methods are addressed. Results from a survey of abandoned rocket upper stages in LEO, examples of abandoned payloads and observations of high area-to-mass ratio debris will be resented. Eventually, first results of the analysis of these light curves are provided.

Characterising the three-dimensional ozone distribution of a tidally locked Earth-like planet

We simulate the 3D ozone distribution of a tidally locked Earth-like exoplanet using the high-resolution, 3D chemistry climate model CESM1(WACCM) and study how the ozone layer of a tidally locked Earth (TLE) (ΩTLE = 1/365 days) differs from that of our present-day Earth (PDE) (ΩPDE = 1/1 day). The middle atmosphere reaches a steady state a symptotically within the first 80 days of the simulation. An upwelling, centred on the subsolar point, is present on the day side while a downwelling, centred on the antisolar point, is present on the night side. In the mesosphere, we find similar global ozone distributions for the TLE and the PDE, with decreased ozone on the day side and enhanced ozone on the night side. In the lower mesosphere, a jet stream transitions into a large-scale vortex around a low-pressure system, located at low latitudes of the TLE night side. In the middle stratosphere, the concentration of odd oxygen is approximately equal to that of the ozone [(Ox) ≈ (O3)]. At these altitudes, the lifetime of odd oxygen is ~16 h and the transport processes significantly contribute to the global distribution of stratospheric ozone. Compared to the PDE, where the strong Coriolis force acts as a mixing barrier between low and high latitudes, the transport processes of the TLE are governed by jet streams variable in the zonal and meridional directions. In the middle stratosphere of the TLE, we find high ozone values on the day side, due to the increased production of atomic oxygen on the day side, where it immediately recombines with molecular oxygen to form ozone. In contrast, the ozone is depleted on the night side, due to changes in the solar radiation distribution and the presence of a downwelling. As a result of the reduced Coriolis force, the tropical and extratropical air masses are well mixed and the global temperature distribution of the TLE stratosphere has smaller horizontal gradients than the PDE. Compared to the PDE, the total ozone column global mean is reduced by ~19.3 %. The day side and the night side total ozone column means are reduced by 23.21 and 15.52 %, respectively. Finally, we present the total ozone column (TOC) maps as viewed by a remote observer for four phases of the TLE during its revolution around the star. The mean TOC values of the four phases of the TLE vary by up to 23 %.

A new, fast and flexible radiative transfer method for Venus general circulation models

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation (0.1-5.5 mu m) and thermal radiation (1.7-260 mu m). The solar radiation calculation is based on the delta-Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative-convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models. (C) 2014 Elsevier Ltd. All rights reserved.