Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

Methodology for the combination of sub-daily Earth rotation from GPS and VLBI observations

A combination procedure of Earth orientation parameters from Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI) observations was developed on the basis of homogeneous normal equation systems. The emphasis and purpose of the combination was the determination of sub-daily polar motion (PM) and universal time (UT1) for a long time-span of 13 years. Time series with an hourly resolution and a model for tidal variations of PM and UT1-TAI (dUT1) were estimated. In both cases, 14-day nutation corrections were estimated simultaneously with the ERPs. Due to the combination procedure, it was warranted that the strengths of both techniques were preserved. At the same time, only a minimum of de-correlating or stabilizing constraints were necessary. Hereby, a PM time series was determined, whose precision is mainly dominated by GPS observations. However, this setup benefits from the fact that VLBI delivered nutation and dUT1 estimates at the same time. An even bigger enhancement can be seen for the dUT1 estimation, where the high-frequency variations are provided by GPS, while the long term trend is defined by VLBI. The estimated combined tidal PM and dUT1 model was predominantly determined from the GPS observations. Overall, the combined tidal model for the first time completely comprises the geometrical benefits of VLBI and GPS observations. In terms of root mean squared (RMS) differences, the tidal amplitudes agree with other empirical single-technique tidal models below 4 μas in PM and 0.25 μs in dUT1. The noise floor of the tidal ERP model was investigated in three ways resulting in about 1 μas for diurnal PM and 0.07 μs for diurnal dUT1 while the semi-diurnal components have a slightly better accuracy.     

Preparing for climate change in Washington State

Climate change is expected to bring potentially significant changes to Washington State’s natural, institutional, cultural, and economic landscape. Addressing climate change impacts will require a sustained commitment to integrating climate information into the day-to-day governance and management of infrastructure, programs, and services that may be affected by climate change. This paper discusses fundamental concepts for planning for climate change and identifies options for adapting to the climate impacts evaluated in the Washington Climate Change Impacts Assessment. Additionally, the paper highlights potential avenues for increasing flexibility in the policies and regulations used to govern human and natural systems in Washington.     

Evaporation estimates from Nasser Lake, Egypt, based on three floating station data and Bowen ratio energy budget

Nasser Lake is located in a hyper-arid region in the south of Egypt. Evaporation is by far the most important factor in explaining the water losses from the lake. To obtain better management scenarios for Nasser Lake, an accurate estimation of the lake evaporation losses thus is essential. This paper presents an update of previous evaporation estimates, making use of local meteorological and hydrological data collected from instrumented platforms (floating weather stations) at three locations on the lake: at Raft, Allaqi, and Abusembel (respectively at 2, 75, and 280 km upstream of the Aswan High Dam). Results from six conventional evaporation quantification methods were compared with the values obtained by the Bowen ratio energy budget method (BREB). The results of the BREB method showed that there is no significant difference between the evaporation rates at Allaqi and Abusembel. At Raft, higher evaporation rates were obtained, which were assumed to be overestimated due to the high uncertainty of the Bowen ratio (BR) parameter. The average BR value at Allaqi and Abusembel was used to eliminate this overestimates evaporation. Variance-based sensitivity and uncertainty analyses on the BREB results were conducted based on quasi-Monte Carlo sequences (Latin Hypercube sampling). The standard deviation of the total uncertainty on the BREB evaporation rate was found to be 0.62 mm day^?1. The parameter controlling the change in stored energy, followed by the BR parameter, was found to be the most sensitive parameters. Several of the six conventional methods showed substantial bias when compared with the BREB method. These were modified to eliminate the bias. When compared to the BREB-based values, the Penman method showed most favorably for the daily time scale, while for the monthly scale, the Priestley–Taylor and the deBruin–Keijman methods showed best agreement. Differences in mean evaporation estimates of these methods (against the BREB method) were found to be in the range 0.14 and 0.36 mm day^?1. All estimates were based calculations at the daily time scale covering a 10-year period (1995–2004).     

Testing GISS-MM5 physics configurations for use in regional impacts studies

The Mesoscale Modeling System Version 5 (MM5) was one-way nested to the Goddard Institute for Space Studies global climate model (GISS GCM), which provided the boundary conditions for present (1990s) and future (IPCC SRES A2 scenario, 2050s) five-summer “time-slice” simulations over the continental and eastern United States. Five configurations for planetary boundary layer, cumulus parameterization, and radiation scheme were tested, and one set was selected for use in the New York City Climate and Health Project—a multi-disciplinary study investigating the effects of climate change and land-use change on human health in the New York metropolitan region. Although hourly and daily data were used in the health project, in this paper we focus on long-term current and projected mean climate change. The GISS-MM5 was very sensitive to the choice of cumulus parameterization and planetary boundary layer scheme, leading to significantly different temperature and precipitation outcomes for the 1990s. These differences can be linked to precipitation type (convective vs. non-convective), to their effect on solar radiation received at the ground, and ultimately to surface temperature. The projected changes in climate (2050s minus 1990s) were not as sensitive to choice of model physics combination. The range of the projected surface temperature changes at a given grid point among the model versions was much less than the mean change for all five model configurations, indicating relative consensus for simulating surface temperature changes among the different model projections. The MM5 versions, however, offer less consensus regarding 1990s to 2050s changes in precipitation amounts. All of the projected 2050s temperature changes were found to be significant at the 95th percent confidence interval, while the majority of the precipitation changes were not.     

Mapping the magnetosphere of PSR B1055−52

We present a geometric study of the radio and γ-ray pulsar B1055−52 based on recent observations at the Parkes radio telescope. We conclude that the pulsar's magnetic axis is inclined at an angle of 75° to its rotation axis and that both its radio main pulse and interpulse are emitted at the same height above their respective poles. This height is unlikely to be higher or much lower than 700 km, a typical value for radio pulsars.
It is argued that the radio interpulse arises from emission formed on open fieldlines close to the magnetic axis which do not pass through the magnetosphere's null (zero-charge) surface. However, the main pulse emission must originate from fieldlines lying well outside the polar cap boundary beyond the null surface, and farther away from the magnetic axis than those of the outer gap region where the single γ-ray peak is generated. This casts doubt on the common assumption that all pulsars have closed, quiescent, corotating regions stretching to the light cylinder.     

57-second oscillations in Nova Centauri 1986 (V842 Cen)

High-speed photometry in 2008 shows that the light curve of V842 Cen possesses a coherent modulation at 56.825 s, with sidebands at 56.598 and 57.054 s. These have appeared since this nova remnant was observed in 2000 and 2002. We deduce that the dominant signal is the rotation period of the white dwarf primary and the sidebands are caused by reprocessing from a surface moving with an orbital period of 3.94 h. Thus, V842 Cen is an intermediate polar (IP) of the DQ Herculis subclass, is the fastest rotating white dwarf among the IPs and is the third fastest known in a cataclysmic variable. As in other IPs, we see no dwarf nova oscillations, but there are often quasi-periodic oscillations in the range 350–1500 s. There is a strong brightness modulation with a period of 3.78 h, which we attribute to negative superhumps, and there is an even stronger signal at 2.886 h which is of unknown origin but is probably a further example of that seen in GW Lib and some other systems. We used the Swift satellite to observe V842 Cen in the ultraviolet and in X-rays, although no periodic modulation was detected in the short observations. The X-ray luminosity of this object appears to be much lower than that of other IPs in which the accretion region is directly visible.     

Melting and Sublimation of Planetary Ices Under Low Pressure Conditions: Laboratory Experiments with a Melting Probe Prototype

One possibility to explore the subsurface layers of icy bodies is to use a probe with a “hot tip", which is able to penetrate ice layers by melting. Such probes have been built and used in the past for the exploration of terrestrial polar ice sheets and may also become useful tools to explore other icy layers in the Solar System. Examples for such layers are the polar areas of Mars or the icy crust of Jupiter’s moon Europa. However, while on Earth a heated probe launched into an ice sheet always causes melting with subsequent refreezing, the behaviour of such a probe in a low pressure environment is quite different. We report on the results of some experiments with a simple “melting probe" prototype with two different kinds of hot tips in a vacuum environment. For one of the tips the probe moved into two types of ice samples: (i) compact water ice and (ii) porous water ice with a snow (firn) like texture. It was also found that the penetration behaviour was basically different for the two sample types even when the same kind of tip was used. While in the porous sample the ice was only subliming, the phase changes occurring during the interaction of the tip with the compact ice are much more complex. Here alternating phases of melting and sublimation occur. The absence of the liquid phase has severe consequences on the performance of a “melting probe" under vacuum conditions: In this environment we find a high thermal resistance between the probe surface and the underlying ice. Therefore, only a low percentage of the heat that is generated in the tip is used to melt or sublime the ice, the bulk of the power is transferred towards the rear end of the probe. This is particularly a problem in the initial phases of an ice penetration experiment, when the probe has not yet penetrated the ice over its whole length. In the compact ice sample, phases could be observed, where a high enough gas pressure had built up locally underneath the probe, so that melting becomes possible. Only during these melting periods the thermal contact between the probe and the ice is good and in consequence the melting probe works effectively.     

Comparison of Positron Emission Tomography and X-ray radiography for studies of physical processes in sandstone

The migration of moisture in cores of porous homogeneous sandstone of Ledian age (Belgian Eocene, Tertiary) is monitored with PET (Positron Emission Tomography) and micro-focus X-ray radiography. In the case of PET, a nuclear medical imaging technique, ⁵⁵Co-EDTA (Ethylenediamine Tetraacetic acid) and R-¹⁸F were used as water-soluble tracers. The X-ray projection method has evolved from the better-known medical technique and allows a fast and accurate determination of the two-dimensional transient moisture content profiles. Results indicate that both techniques can deliver important information concerning physical processes in situ.