Realisation of a celestial system based on VLBI geodynamics programs

The relative orientations of various VLBI celestial reference frames are evaluated on the basis of coordinate differences of common sources. It is shown that an accuracy better than 0.001 can be achieved. Possible regional deformations in the different catalogues are investigated; they are found to reach a few 0.001 in some restricted zones. The application of these studies to the realisation of a combined celestial reference frame consistent with the BIH Terrestrial System is outlined.     

Processing of high-frequency multibeam echo sounder data for seafloor characterization

Processing simultaneous bathymetry and backscatter data, multibeam echosounders (MBESs) show promising abilities for remote seafloor characterization. High-frequency MBESs provide a good horizontal resolution, making it possible to distinguish fine details at the water-seafloor interface. However, in order to accurately measure the seafloor influence on the backscattered energy, the recorded sonar data must first be processed and cleared of various artifacts generated by the sonar system itself. Such a preprocessing correction procedure along with the assessment of its validity limits is presented and applied to a 95-kHz MBES (Simrad EM 1000) data set. Beam pattern effects, uneven array sensitivities, and inaccurate normalization of the ensonified area are removed to make possible further quantitative analysis of the corrected backscatter images. Unlike low-frequency data where the average backscattered energy proves to be the only relevant feature for discriminating the nature of the seafloor, high-frequency MBES backscatter images exhibit visible texture patterns. This additional information involves different statistical distributions of the backscattered amplitudes obtained from various seafloor types. Non-Rayleigh statistics such as K-distributions are shown to fit correctly the skewed distributions of experimental high-frequency data. Apart from the effect of the seafloor micro-roughness, a statistical model makes clear a correlation between the amplitude statistical distributions and the signal incidence angle made available by MBES bathymetric abilities. Moreover, the model enhances the effect of the first derivative of the seafloor backscattering strength upon statistical distributions near the nadir and at high incidence angles. The whole correction and analysis process is finally applied to a Simrad EM 1000 data set.     

Applications géodésiques du système DORIS à l'Institut géographique national

IGN is in charge of the installation and maintenance of the DORIS orbit determination network. More recently, in collaboration with JPL, precise geodetic computations were performed. The goal of this paper is to recall the various historic contributions of IGN to the DORIS system in their international context and then to describe a new estimation technique developed for a multi-satellite mode, making full profit of a better modeling for satellites and ground clocks as well as tropospheric correction parameters. Derived geodetic results demonstrate a precision in the order of 1 cm for station positions. To cite this article: P. Willis et al., C. R. Geoscience 337 (2005).     

Review: Theory-guided machine learning applied to hydrogeology—state of the art,opportunities and future challenges

Hydrogeology Journal - Thanks to recent technological advances, hydrogeologists now have access to large amounts of data acquired in real time. Processing these data using traditional modelling...     

Sediment yield in South Africa — A preliminary geographical analysis

The study examines the spatial relationships between sediment yield and 15 independent environmental variables in 54 catchments in South Africa. Rooseboom's (1978) data on the sediment yield from the catchments were standardized for a single time period. Bivariate regression analyses reveal no simple relationships. Multivariate regression analyses conducted for the whole and various sub-areas of South Africa indicate that latitude and longitude are the primary variables affecting spatial variations in sediment yield. This may be as a result of latitude and longitude being surrogate variables reflecting variation in other environmental variables (e. g. geology, vegetation). Within the sub-areas, 43.4% to 97.8% of the variation in sediment yield is explained by the combined variation in a number of different environmental variables. This study highlights a need for the collection and analysis of more sediment yield data, which would allow the analyses to be refined, to predict sediment yields from ungauged catchments in South Africa.     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

A review of algebraic constraints in terrestrial reference frame datum definition

 Station coordinates are combined with velocities estimated by space geodesy techniques to produce the International Terrestrial Reference System. The input is sets of coordinates and velocities calculated by International Earth Rotation Service analysis centers using space geodesy techniques. The working reference system of individual analysis centers is generally conventionally defined. However, the implications of such processing can have an effect on the resulting combined set. The problem of datum definition as a function of coordinate combinations is reviewed. In particular, the problem of minimum constraints is clearly emphasized and the reference system effect is defined. The goal is to build a process that could be used generally to remove uncertainties in the underlying coordinate system without disturbing the underlying information with additional unnecessary information. Received: 25 January 1999 / Accepted: 20 September 2000     

Improvement of the transformation between ITRF and Doppler-Realized WGS84

The WGS84 (World Geodetic System 1984) reference system is, originally, mathematically defined from the NSWC-9Z2 (Naval Surface Weapons Center — 9Z2) reference system. The WGS84 associated realization, called in this paper WGS84-D, is a 1 meter consistency NNSS (US Navy Navigation Satellite System) Doppler realized reference frame. In contrast, the ITRF (IERS Terrestrial Reference Frame) is a 1 centimeter consistency reference frame realized through the most accurate techniques of Space Geodesy. This work intends to improve the transformation parameters between the WGS84-D and the ITRF through the use of both a NSWC-9Z2/Doppler realization and an extension of the ITRF network. A strong linear correlation was also modeled between the Doppler determined scale factor and the mean smoothed sunspot number, due to uncompensated ionospheric effects. This correction improved NSWC-9Z2 (i.e. WGS84) Doppler realization consistency. The uncertainty of adjusted transformation parameters between the ITRF and the WGS84-D is improved by a factor 2 over previous determinations.     

DORIS Contribution to ITRF2005

We examine the contribution of the Doppler Orbit determination and Radiopositioning Integrated by Satellite (DORIS) technique to the International Terrestrial Reference Frame (ITRF2005) by evaluating the quality of the submitted solutions as well as that of the frame parameters, especially the origin and the scale. Unlike the previous versions of the ITRF, ITRF2005 is constructed with input data in the form of time-series of station positions (weekly for satellite techniques and daily for VLBI) and daily Earth orientation parameters (EOPs), including full variance–covariance information. Analysis of the DORIS station positions’ time-series indicates an internal precision reaching 15 mm or better, at a weekly sampling. A cumulative solution using 12 years of weekly time-series was obtained and compared to a similar International GNSS Service (IGS) GPS solution (at 37 co-located sites) yielding a weighted root mean scatter (WRMS) of the order of 8 mm in position (at the epoch of minimum variance) and about 2.5 mm/year in velocity. The quality of this cumulative solution resulting from the combination of two individual DORIS solutions is better than any individual solution. A quality assessment of polar motion embedded in the contributed DORIS solutions is performed by comparison with the results of other space-geodetic techniques and in particular GPS. The inferred WRMS of polar motion varies significantly from one DORIS solution to another and is between 0.5 and 2 mas, depending on the strategy used and in particular estimating or not polar motion rate by the analysis centers. This particular aspect certainly needs more investigation by the DORIS Analysis Centers.     

Blind Submarine Seismic Deconvolution for Long Source Wavelets

In seismic deconvolution, blind approaches must be considered in situations where reflectivity sequence, source wavelet signal, and noise power level are unknown. In the presence of long source wavelets, strong interference among the reflectors contributions makes the wavelet estimation and deconvolution more complicated. In this paper, we solve this problem in a two-step approach. First, we estimate a moving average (MA) truncated version of the wavelet by means of a stochastic expectation-maximization (SEM) algorithm. Then, we use Prony's method to improve the wavelet estimation accuracy by fitting an autoregressive moving average (ARMA) model with the initial truncated wavelet. Moreover, a solution to the wavelet initialization problem in the SEM algorithm is also proposed. Simulation and real-data experiment results show the significant improvement brought by this approach.