Diurnal atmospheric forcing and temporal variations of the nutation amplitudes

A 29-year time-series of four-times-daily atmospheric effective angular momentum (EAM) estimates is used to study the atmospheric influence on nutation. The most important atmospheric contributions are found for the prograde annual (77 μas), retrograde annual (53 as), prograde semiannual (45 as), and for the constant offset of the pole (δψsinɛ₀=−86 as, δɛ=77 as). Among them only the prograde semiannual component is driven mostly by the wind term of the EAM function, while in all other cases the pressure term is dominant. These are nonnegligible quantities which should be taken into account in the new theory of nutation. Comparison with the VLBI corrections to the IAU 1980 nutation model taking into account the ocean tide contribution yields good agreement for the prograde annual and semiannual nutations. We also investigated time variability of the atmospheric contribution to the nutation amplitudes by performing the sliding-window least-squares analysis of both the atmospheric excitation and VLBI nutation data. Almost all detected variations of atmospheric origin can be attributed to the pressure term, the biggest being the in-phase annual prograde component (about 30 as) and the retrograde one (as much as 100200 as). These variations, if physical, limit the precision of classical modeling of nutation to the level of 0.1 mas. Comparison with the VLBI data shows significant correlation for the retrograde annual nutation after 1989, while for the prograde annual term there is a high correlation in shape but the size of the atmospherically driven variations is about three times less than deduced from the VLBI data. This discrepancy in size can be attributed either to inaccuracy of the theoretical transfer function or the frequency-dependent ocean response to the pressure variations. Our comparison also yields a considerably better agreement with the VLBI nutation data when using the EAM function without the IB correction for ocean response, which indicates that this correction is not adequate for nearly diurnal variations. Received: 10 September 1997 / Accepted: 5 March 1998     

Airborne gravimetry used in precise geoid computations by ring integration

Two detailed geoids have been computed in the region of North Jutland. The first computation used marine data in the offshore areas. For the second computation the marine data set was replaced by the sparser airborne gravity data resulting from the AGMASCO campaign of September 1996. The results of comparisons of the geoid heights at on-shore geometric control showed that the geoid heights computed from the airborne gravity data matched in precision those computed using the marine data, supporting the view that airborne techniques have enormous potential for mapping those unsurveyed areas between the land-based data and the off-shore marine or altimetrically derived data. Received: 7 July 1997 / Accepted: 22 April 1998     

Fast integration of low orbiter's trajectory perturbed by the earth's non-sphericity

A fast algorithm is proposed to integrate the trajectory of a low obiter perturbed by the earth's non-sphericity. The algorithm uses a separation degree to define the low-degree and the high-degree acceleration components, the former computed rigorously, and the latter interpolated from gridded accelerations. An FFT method is used to grid the accelerations. An optimal grid type for the algorithm depends on the trajectory's permissible error, speed, and memory capacity. Using the non-spherical accelerations computed from EGM96 to harmonic degree 360, orbit integrations were performed for a low orbiter at an altitude of 170 km. For a separation degree of 50, the new algorithm, together with the predict-pseudo correct method, speeds up the integration by 145 times compared to the conventional algorithm while keeping the errors in position and velocity below 10⁻⁴ m and 10⁻⁷ m/s for a 3-day arc. Received: 28 July 1997 / Accepted: 1 April 1998     

Reference frame stability and nonlinear distortion in minimum-constrained network adjustment

The aim of this paper is to investigate the influence of the minimum constraints (MCs) on the reference frame parameters in a free-net solution. The non-estimable part of these parameters (which is not reduced by the available data) is analysed in terms of its stability under a numerical perturbation of the constrained datum functionals. In practice, such a perturbation can be ascribed either to hidden errors in the known coordinates/velocities that participate in the MCs or to a simple change of their a priori values due to a datum switch on a different fiducial dataset. In addition, a perturbation of this type may cause a nonlinear variation to the estimable part of the reference frame parameters, since it theoretically affects the adjusted observations that are implied by the network’s nonlinear observational model. The aforementioned effects have an impact on the quality of a terrestrial reference frame (TRF) that is established via a minimum-constrained adjustment, and our study shows that they are both controlled through a characteristic matrix which is inherently linked to the MC system (the so-called TRF stability matrix). The structure of this matrix depends on the network’s spatial configuration and the ‘geometry’ of the datum constraints, while its main role is the filtering of any MC-related errors into the least-squares adjustment results. A number of examples with different types of geodetic networks are also presented to demonstrate the theoretical findings of our study.     

Cascaded image analysis for dynamic crack detection in material testing

Concrete probes in civil engineering material testing often show fissures or hairline-cracks. These cracks develop dynamically. Starting at a width of a few microns, they usually cannot be detected visually or in an image of a camera imaging the whole probe. Conventional image analysis techniques will detect fissures only if they show a width in the order of one pixel. To be able to detect and measure fissures with a width of a fraction of a pixel at an early stage of their development, a cascaded image analysis approach has been developed, implemented and tested. The basic idea of the approach is to detect discontinuities in dense surface deformation vector fields. These deformation vector fields between consecutive stereo image pairs, which are generated by cross correlation or least squares matching, show a precision in the order of 1/50 pixel. Hairline-cracks can be detected and measured by applying edge detection techniques such as a Sobel operator to the results of the image matching process. Cracks will show up as linear discontinuities in the deformation vector field and can be vectorized by edge chaining. In practical tests of the method, cracks with a width of 1/20 pixel could be detected, and their width could be determined at a precision of 1/50 pixel.     

Phase center modeling for LEO GPS receiver antennas and its impact on precise orbit determination

Most satellites in a low-Earth orbit (LEO) with demanding requirements on precise orbit determination (POD) are equipped with on-board receivers to collect the observations from Global Navigation Satellite systems (GNSS), such as the Global Positioning System (GPS). Limiting factors for LEO POD are nowadays mainly encountered with the modeling of the carrier phase observations, where a precise knowledge of the phase center location of the GNSS antennas is a prerequisite for high-precision orbit analyses. Since 5 November 2006 (GPS week 1400), absolute instead of relative values for the phase center location of GNSS receiver and transmitter antennas are adopted in the processing standards of the International GNSS Service (IGS). The absolute phase center modeling is based on robot calibrations for a number of terrestrial receiver antennas, whereas compatible antenna models were subsequently derived for the remaining terrestrial receiver antennas by conversion (from relative corrections), and for the GNSS transmitter antennas by estimation. However, consistent receiver antenna models for space missions such as GRACE and TerraSAR-X, which are equipped with non-geodetic receiver antennas, are only available since a short time from robot calibrations. We use GPS data of the aforementioned LEOs of the year 2007 together with the absolute antenna modeling to assess the presently achieved accuracy from state-of-the-art reduced-dynamic LEO POD strategies for absolute and relative navigation. Near-field multipath and cross-talk with active GPS occultation antennas turn out to be important and significant sources for systematic carrier phase measurement errors that are encountered in the actual spacecraft environments. We assess different methodologies for the in-flight determination of empirical phase pattern corrections for LEO receiver antennas and discuss their impact on POD. By means of independent K-band measurements, we show that zero-difference GRACE orbits can be significantly improved from about 10 to 6 mm K-band standard deviation when taking empirical phase corrections into account, and assess the impact of the corrections on precise baseline estimates and further applications such as gravity field recovery from kinematic LEO positions.     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

Comparing spatially varying coefficient models: a case study examining violent crime rates and their relationships to alcohol outlets and illegal drug arrests

In this paper, we compare and contrast a Bayesian spatially varying coefficient process (SVCP) model with a geographically weighted regression (GWR) model for the estimation of the potentially spatially varying regression effects of alcohol outlets and illegal drug activity on violent crime in Houston, Texas. In addition, we focus on the inherent coefficient shrinkage properties of the Bayesian SVCP model as a way to address increased coefficient variance that follows from collinearity in GWR models. We outline the advantages of the Bayesian model in terms of reducing inflated coefficient variance, enhanced model flexibility, and more formal measuring of model uncertainty for prediction. We find spatially varying effects for alcohol outlets and drug violations, but the amount of variation depends on the type of model used. For the Bayesian model, this variation is controllable through the amount of prior influence placed on the variance of the coefficients. For example, the spatial pattern of coefficients is similar for the GWR and Bayesian models when a relatively large prior variance is used in the Bayesian model.      

Evaluation of pre-CHAMP gravity models `GRIM5-S1' and `GRIM5-C1' with satellite crossover altimetry

 The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients) in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified, i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic. Received: 20 February 2001 / Accepted: 24 October 2001