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.

SOLAR RADIATION PRESSURE AND LOCAL INTERSTELLAR MEDIUM FLOW PARAMETERS FROM INTERSTELLAR BOUNDARY EXPLORER LOW ENERGY HYDROGEN MEASUREMENTS

Neutral hydrogen atoms that travel into the heliosphere from the local interstellar medium (LISM) experience strong effects due to charge exchange and radiation pressure from resonant absorption and re-emission of Lyα. The radiation pressure roughly compensates for the solar gravity. As a result, interstellar hydrogen atoms move along trajectories that are quite different than those of heavier interstellar species such as helium and oxygen, which experience relatively weak radiation pressure. Charge exchange leads to the loss of primary neutrals from the LISM and the addition of new secondary neutrals from the heliosheath. IBEX observations show clear effects of radiation pressure in a large longitudinal shift in the peak of interstellar hydrogen compared with that of interstellar helium. Here, we compare results from the Lee et al. interstellar neutral model with IBEX-Lo hydrogen observations to describe the distribution of hydrogen near 1 AU and provide new estimates of the solar radiation pressure. We find over the period analyzed from 2009 to 2011 that radiation pressure divided by the gravitational force (μ) has increased slightly from μ = 0.94 ± 0.04 in 2009 to μ = 1.01 ± 0.05 in 2011. We have also derived the speed, temperature, source longitude, and latitude of the neutral H atoms and find that these parameters are roughly consistent with those of interstellar He, particularly when considering the filtration effects that act on H in the outer heliosheath. Thus, our analysis shows that over the period from 2009 to 2011, we observe signatures of neutral H consistent with the primary distribution of atoms from the LISM and a radiation pressure that increases in the early rise of solar activity.

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.

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.

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.

Upconverter materials and upconversion solar-cell devices: simulation and characterization with broad solar spectrum illumination

Upconverter materials and upconverter solar devices were recently investigated with broad-band excitation revealing the great potential of upconversion to enhance the efficiency of solar cell at comparatively low solar concentration factors. In this work first attempts are made to simulate the behavior of the upconverter β-NaYF4 doped with Er3+ under broad-band excitation. An existing model was adapted to account for the lower absorption of broader excitation spectra. While the same trends as observed for the experiments were found in the simulation, the absolute values are fairly different. This makes an upconversion model that specifically considers the line shape function of the ground state absorption indispensable to achieve accurate simulations of upconverter materials and upconverter solar cell devices with broadband excitations, such as the solar radiation.

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

Infantile spasms: does season influence onset and long-term outcome?

To study whether onset of infantile spasms manifests seasonal variation, as previously reported, and whether any such seasonality is associated with treatment response and long-term outcome, data for 57 patients were retrospectively reviewed. The data were collected from hospital files and through a mail survey of children with infantile spasms born from 1980 to 2002 and monitored at the University Children's Hospital of Berne, Switzerland. The mean age at time of onset of infantile spasms was 7 months (range, 0.75-40), at diagnosis 8 months (range, 1-42) and at follow-up 11.3 years (range, 1-23 years). In 77% of participants, the etiology of infantile spasms was known (symptomatic); in the remaining 23% it was not known (nonsymptomatic). In contrast to previous findings, onset of infantile spasms was not associated with calendar month, photoperiod, or global solar radiation. Long-term prognosis was poor: 4 of the 57 (7%) children died; 49 (86%) had cognitive impairment and 40 (70%) had physical impairment; 31 (54%) had cerebral palsy, 37 had (65%) persistent seizures, and 9 (16%) had Lennox-Gastaut syndrome. Symptomatic infantile spasms were associated with worse cognitive outcome (P < 0.001), but treatment modality and overall duration of infantile spasms were not. There was no association of calendar month or photoperiod at onset with cognitive outcome or treatment response.

Climate and chemistry effects of a regional scale nuclear conflict

Previous studies have highlighted the severity of detrimental effects for life on earth after an assumed regionally limited nuclear war. These effects are caused by climatic, chemical and radiative changes persisting for up to one decade. However, so far only a very limited number of climate model simulations have been performed, giving rise to the question how realistic previous computations have been. This study uses the coupled chemistry climate model (CCM) SOCOL, which belongs to a different family of CCMs than previously used, to investigate the consequences of such a hypothetical nuclear conflict. In accordance with previous studies, the present work assumes a scenario of a nuclear conflict between India and Pakistan, each applying 50 warheads with an individual blasting power of 15 kt ("Hiroshima size") against the major population centers, resulting in the emission of tiny soot particles, which are generated in the firestorms expected in the aftermath of the detonations. Substantial uncertainties related to the calculation of likely soot emissions, particularly concerning assumptions of target fuel loading and targeting of weapons, have been addressed by simulating several scenarios, with soot emissions ranging from 1 to 12 Tg. Their high absorptivity with respect to solar radiation leads to a rapid self-lofting of the soot particles into the strato- and mesosphere within a few days after emission, where they remain for several years. Consequently, the model suggests earth's surface temperatures to drop by several degrees Celsius due to the shielding of solar irradiance by the soot, indicating a major global cooling. In addition, there is a substantial reduction of precipitation lasting 5 to 10 yr after the conflict, depending on the magnitude of the initial soot release. Extreme cold spells associated with an increase in sea ice formation are found during Northern Hemisphere winter, which expose the continental land masses of North America and Eurasia to a cooling of several degrees. In the stratosphere, the strong heating leads to an acceleration of catalytic ozone loss and, consequently, to enhancements of UV radiation at the ground. In contrast to surface temperature and precipitation changes, which show a linear dependence to the soot burden, there is a saturation effect with respect to stratospheric ozone chemistry. Soot emissions of 5 Tg lead to an ozone column reduction of almost 50% in northern high latitudes, while emitting 12 Tg only increases ozone loss by a further 10%. In summary, this study, though using a different chemistry climate model, corroborates the previous investigations with respect to the atmospheric impacts. In addition to these persistent effects, the present study draws attention to episodically cold phases, which would likely add to the severity of human harm worldwide. The best insurance against such a catastrophic development would be the delegitimization of nuclear weapons.

LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models

The African great lakes are of utmost importance for the local economy (fishing), as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E), East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP). Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model’s performance. Simulations are performed over the freshwater layer only (60 m) and over the average lake depth (240 m), since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

Assessment of detectability of neutral interstellar deuterium by IBEX observations

Context. The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. Aims: We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). Methods: Using both previous and the most recent determinations of the flow parameters of neutral gas in the local interstellar cloud (LIC) and an observation-based model of solar radiation pressure and ionization in the heliosphere, we simulated the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assessed the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulated the observations expected during an epoch of high solar activity. In addition, we calculated the expected counts of D atoms from the thin terrestrial water layer covering the IBEX-Lo conversion surface, sputtered by neutral interstellar He atoms. Results: Most D counts registered by IBEX-Lo are expected to come from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable to the expected number of sputtered D atoms registered during the same time intervals. Conclusions: The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, say, by making observations in the special deuterium mode of IBEX-Lo.