Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

MELTWATER FLUXES AT AN ARCTIC FOREST-TUNDRA SITE

Models of surface energy balance and snow metamorphism are utilized to predict the energy and meltwater fluxes at an Arctic site in the forest–tundra transition zone of north-western Canada. The surface energy balance during the melt period is modelled using an hourly bulk aerodynamic approach. Once a snowcover becomes patchy, advection from the bare patches to the snow-covered areas results in a large spatial variation in basin snowmelt. In order to illustrate the importance of small-scale, horizontal advection, a simple parameterization scheme using sensible heat fluxes from snow free areas was tested. This scheme estimates the maximum horizontal advection of sensible heat from the bare patches to the snow-covered areas. Calculated melt was routed through the measured snowcover in each landscape type using a variable flow path, meltwater percolation model. This allowed the determination of the spatial variability in the timing and magnitude of meltwater release for runoff. Model results indicate that the initial release of meltwater first occurred on the shallow upland tundra sites, but meltwater release did not occur until nearly two weeks later on the deep drift snowcovers. During these early periods of melt, not all meltwater is available for runoff. Instead, there is a period when some snowpacks are only partially contributing to runoff, and the spatial variation of runoff contribution corresponds to landscape type. Comparisons of melt with and without advection suggests that advection is an important process controlling the timing of basin snowmelt.     

VELOCITY DISPERSION OF SEISMIC WAVES*

In well velocity surveys made to calibrate Sonic (CV) Logs the calibration survey uses frequencies around 50 Hz whereas the Sonic Logging tool uses frequencies around 20 kHz. There thus exists the possibility of making a direct measure of velocity dispersion. In any one survey the disturbing factors, both instrumental and operational, will often mask any dispersive effect that might exist. Consequently this paper reports on a statistical analysis of the velocity differences resulting from calibration surveys and Sonic logs. Only Borehole Compensated Sonic Logs were used. Four areas were investigated: the North Sea, Abu Dhabi, Libya and Alaska. After rejecting logs and calibration records which were obviously in error there remained 424000 feet (about 130 km) of usable log distributed throughout 66 wells. The four areas were analysed separately and in no case was the estimated dispersion significantly different from zero. However, the mean values did correlate with lithology from (? 0.17 ± 0.18)% for the essentially carbonate section in Abu Dhabi to (+ 0.45 ± 0.25)% for the sand-shale section in Alaska, a positive sign meaning that the higher frequencies travelled faster. Except for Alaska the calibration surveys were made with a wall-clamp geophone, and for these areas amplitude measurements were made. After suitable corrections estimates of the absorption parameter Q were obtained. These varied from 20 to 200 with mean values of 63 for Libya, 70 for Abu Dhabi and 88 for the North Sea (excluding the Tertiary). If, as is usually assumed, the absorption mechanism is linear and is described by a Q which is independent of frequency, then these values would necessarily imply dispersion of several percent. As instanced above no such dispersion was observed. It is possible that the expected dispersion was compensated for by invasion of the mud filtrate into the borehole walls, but it is more likely that the absorption mechanism was substantially non-linear.     

The neotectonics and geodynamics of the Lower Oka Region (East European Craton)

The neotectonic structures of the Lower Oka (Nizhneokskii) Region formed under different geodynamic conditions. This is attested by the morphology, orientation, internal structure, and jointing of the structures. The Oka-Tsna arc formed under the effect of tension from an inner source on the one hand and stress from the Alpian belt on the other hand. The latitudinally-oriented structures of the northwestern slope of the Tokmovo arc emerged as a result of uplift and widening. Both types of structure are combined within the limits of the Oka-Murom trough, which is a geodynamically active zone.     

Observational study of the seasonality of the submonthly and intraseasonal signal over the tropics

Summary Atmospheric variability in outgoing long-wave radiation (OLR) and tropospheric relative vorticity (VOR) over the South American region was studied from 1979 to 1996 using the complex Morlet wavelet function. The analyses focus on spatial variation in intraseasonal and submonthly scales. Scalograms were used to measure submonthly intraseasonal oscillations in convection, which were found to be predominant in the tropical regions. However, 7-day and 15-day oscillations were observed at tropical and extratropical latitudes in spring and winter, indicating that transient disturbances play a more prominent role. Regarding VOR, tropical energy intensities were highest in the spring and summer, whereas subtropical and extratropical energy intensities were highest in the autumn and winter. The dynamics of the 25-day and 45-day VOR oscillations indicates a possible correlation with Rossby waves over the eastern tropical Pacific Ocean, mainly during the summer. During winter, the 7-day and 15-day VOR oscillations are more frequent at higher latitudes and are enhanced along storm tracks. It was also observed that convection amplitudes in the regions of maximum intensity change appreciably from year to year and from season to season, showing that the behavior of the submonthly and intraseasonal oscillations is nonperiodic and correlates strongly with El Ni?o/Southern Oscillation years. These results confirm the efficiency of wavelet analysis for time-scale studies of atmospheric variability.     

Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Lateral continuity of basement seismic reflections in 15 Ma ultrafast-spreading crust at ODP Site 1256

The Ocean Drilling Program (ODP) initiated drilling at Site 1256D in the Guatemala Basin, about 1,000 km off the East Pacific Rise to penetrate plutonic rocks, anticipated to be relatively shallow in this region, formed at an ultra-fast spreading rate. IODP Expedition E312 successfully drilled into gabbros at ~1,150 m in basement. Multi-channel seismic traces show weak laterally coherent sub-basement reflections at borehole depths. Synthetic reflectivity seismograms were computed using a Ricker wavelet and impedance profiles from borehole sonic logs. These seismograms show significant sub-basement amplitude peaks. A zero-offset vertical seismic profile, shot on E312, was processed to investigate the authenticity of these reflections and their relationship to borehole geology. A dual scheme of the median filtering and F–K dip filtering was used. Tests with synthetic seismograms indicate the approach is effective at reasonable SNR levels. Downgoing energy is clearly identified but negligible upgoing energy is visible over random noise. These results indicate that lava flows and igneous contacts in upper ocean crust have significant topography on lateral scales less than the Fresnel Zone (~300 m) due to igneous and tectonic processes.