Multifrequency algorithms for precise point positioning: MAP3

We present the new MAP3 algorithms to perform static precise point positioning (PPP) from multifrequency and multisystem GNSS observations. MAP3 represents a two-step strategy in which the least squares theory is applied twice to estimate smoothed pseudo-distances, initial phase ambiguities, and slant ionospheric delay first, and the absolute receiver position and its clock offset in a second adjustment. Unlike the classic PPP technique, in our new approach, the ionospheric-free linear combination is not used. The combination of signals from different satellite systems is accomplished by taking into account the receiver inter-system bias. MAP3 has been implemented in MATLAB and integrated within a complete PPP software developed on site and named PCube. We test the MAP3 performance numerically and contrast it with other external PPP programs. In general, MAP3 positioning accuracy with low-noise GPS dual-frequency observations is about 2.5 cm in 2-h observation periods, 1 cm in 10 h, and 7 mm after 1 day. This means an improvement in the accuracy in short observation periods of at least 7 mm with respect to the other PPP programs. The MAP3 convergence time is also analyzed and some results obtained from real triple-frequency GPS and GIOVE observations are presented.     

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.     

From data to decision-making – experiences in Integrated Earth Observation (IEO)

Within the confines of Earth Observation Scientific Knowledge and Technology Transfer in Hungary (EKAT) titled ESA PECS project we gained forum for utilization of our earth observation (EO) knowledge and possibility for further development of experiences. The project aims the horizontal preparation of Hungary for ESA membership – for a dynamically developing sphere in the common part of science, technology and business. For the efficient transfer and dissemination, we preconceived target groups and principles to reach main take-up points and to open new forums of information exchange. Nowadays, user-friendly service, complex supply and real-validation process all have great importance. Joining units of suitable data choice, entire ingestion–processing chains and ergonomic interfaces can only offer a timely and professional solution for the challenges of industry and business sectors. It has to be clearly seen, that information and communication technologies (ICT) and information control becomes a more and more significant part of EO services.     

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

The ionosphere: effects,GPS modeling and the benefits for space geodetic techniques

The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.     

Progressive Transformation of a Flux Rope to an ICME

The solar wind conditions at one astronomical unit (AU) can be strongly disturbed by interplanetary coronal mass ejections (ICMEs). A subset, called magnetic clouds (MCs), is formed by twisted flux ropes that transport an important amount of magnetic flux and helicity, which is released in CMEs. At 1 AU from the Sun, the magnetic structure of MCs is generally modeled by neglecting their expansion during the spacecraft crossing. However, in some cases, MCs present a significant expansion. We present here an analysis of the huge and significantly expanding MC observed by the Wind spacecraft during 9 – 10 November 2004. This MC was embedded in an ICME. After determining an approximate orientation for the flux rope using the minimum variance method, we obtain a precise orientation of the cloud axis by relating its front and rear magnetic discontinuities using a direct method. This method takes into account the conservation of the azimuthal magnetic flux between the inbound and outbound branches and is valid for a finite impact parameter (i.e., not necessarily a small distance between the spacecraft trajectory and the cloud axis). The MC is also studied using dynamic models with isotropic expansion. We have found (6.2±1.5)×10²⁰ Mx for the axial flux and (78±18)×10²⁰ Mx for the azimuthal flux. Moreover, using the direct method, we find that the ICME is formed by a flux rope (MC) followed by an extended coherent magnetic region. These observations are interpreted by considering the existence of a previously larger flux rope, which partially reconnected with its environment in the front. We estimate that the reconnection process started close to the Sun. These findings imply that the ejected flux rope is progressively peeled by reconnection and transformed to the observed ICME (with a remnant flux rope in the front part).     

Spatial and temporal patterns of dry spell lengths in the Iberian Peninsula for the second half of the twentieth century

Summary Daily pluviometric records of 43 meteorological stations across the Iberian Peninsula have permitted a detailed analysis of dry spell patterns for the period 1951–2000 by distinguishing daily amount thresholds of 0.1, 1.0, 5.0 and 10.0 mm/day. The analyses are based on three annual series, namely the number of dry spells, N, the average dry spell length, L, and the extreme dry spell length, L _max. First, the statistical significance of local trends for the annual series of N, L and L _max has been investigated by means of the Mann-Kendall test and significant field trends have been established by means of Monte Carlo simulations. Clear signs of negative field trends are detected for N (1.0 and 10.0 mm/day) and L (0.1 mm/day). Second, the Weibull model fits well the empirical distributions of dry spell lengths for all the rain gauges, whatever the daily amount threshold, with a well ranged spatial distribution of their parameters u and k. On the basis of the Weibull distribution, return period maps for 2, 5, 10, 25 and 50 years have been obtained for dry spell lengths with respect to the four daily threshold levels. While for 0.1 and 1.0 mm/day the longest dry spells are expected at the south of the Iberian Peninsula, for 5.0 and 10.0 mm/day they are mostly detected at the southeast. Finally, the elapsed time between consecutive dry spells has been analysed by considering the same rain amount thresholds and different dry spell lengths at increasing intervals of 10 days. This analysis makes evident a significant negative field trend of the elapsed time between consecutive dry spells of lengths ranging from 10 to 20 days for daily amount thresholds of 1.0, 5.0 and 10.0 mm/day. Authors’ addresses: X. Lana, C. Serra, Departament de Física i Enginyeria Nuclear, ETSEIB, Universitat Politècnica de Catalunya, Av. Diagonal 647 planta 11, 08028 Barcelona, Spain; M. D. Marínez, Departament de Física Aplicada, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain; A. Burgue?o, Departament de Meteorologia i Astronomia, Universitat de Barcelona, 08028 Barcelona, Spain; J. Martín-Vide, L. Gómez, Grup de Climatologia, Universitat de Barcelona, 08028 Barcelona, Spain.     

Geometry and accuracy of reflecting points in bistatic satellite altimetry

The analysis of the time and space distribution of specular (reflecting) points in bistatic altimetry between GPS and CHAMP satellites or SAC-C (taken as examples) is extended from Wagner and Klokočník (2003 J. Geod 77: 128–138). We demonstrate a significantly higher number and density of reflecting points in bistatic altimetry in comparison with traditional monostatic altimetry. After an outline of our older accuracy assessment for the vertical position of the reflecting point, we add a new independent derivation and compare both approaches. We account for orbit errors of both the transmitters (GPS) and receiver (CHAMP) satellites, and the measurement (delay) error. We found that the accuracy of the vertical position of the reflecting point decreases only slowly with increasing off-nadir angle and that the orbit errors must be accounted for if decimeter and better accuracy is required. In this paper, we do not study errors such as state of the ocean, technical parameters of the receiving system, and atmospheric corrections.     

Diurnal Variation of Gravity Wave Activity at Midlatitudes in the Ionospheric F Region

We investigate the short-term fluctuations in the period range from 15 to 180 minutes in the electron density variations of the F region ionosphere. Electron density profiles obtained at the ionospheric stations of Pruhonice (49.9° N, 14.5° E) and Ebro (40.8° N, 0.5° E) at five minute time sampling have been used for this analysis. The diurnal changes of the activity of the acoustic gravity wave fluctuations (AGW) show a clear enhancement during and several hours after sunrise. The periods of such AGW's are about 60 to 75 minutes and these waves propagates vertically through the ionosphere from a source located at an altitude of 180-220 km. The most likely source for these events seems to be passage of the Solar terminator.     

Deep inside a neoproterozoic intra-oceanic arc: growth,differentiation and exhumation of the Amalaoulaou complex (Gourma,Mali)

We show here that the Amalaoulaou complex, in the Pan-African belt of West Africa (Gourma, Mali), corresponds to the lower and middle sections of a Neoproterozoic intra-oceanic arc. This complex records a 90–130-Ma-long evolution of magmatic inputs and differentiation above a subducting oceanic slab. Early c. 793 Ma-old metagabbros crystallised at lower crustal or uppermost mantle depths (25–30 km) and have geochemical characteristic of high-alumina basalts extracted from a depleted mantle source slightly enriched by slab-derived sedimentary components ((La/Sm)_N < 1; ε_Nd: +5.4–6.2; ⁸⁷Sr/⁸⁶Sr: 0.7027–0.7029). In response to crustal thickening, these mafic rocks were recrystallised into garnet-granulites (850–1,000°C; 10–12 kbar) and subject to local dehydration–melting reactions, forming trondhjemititic leucosomes with garnet–clinopyroxene–rutile residues. Slightly after the granulitic event, the arc root was subject to strong HT shearing during partial exhumation (detachment faults/rifting or thrusting), coeval with the emplacement of spinel- and garnet-pyroxenite dykes crystallised from a high-Mg andesitic parental magma. Quartz and hornblende-gabbros (700–660 Ma) with composition typical of hydrous volcanic rocks from mature arcs ((La/Sm)_N: 0.9–1.8; ε_Nd: +4.6 to +5.2; ⁸⁷Sr/⁸⁶Sr: 0.7028–0.7031) were subsequently emplaced at mid-arc crust levels (~15 km). Trace element and isotopic data indicate that magmas tapped a depleted mantle source significantly more enriched in oceanic sedimentary components (0.2%). Exhumation occurred either in two stages (700–660 and 623 Ma) or in one stage (623 Ma) with a final exhumation of the arc root along cold P-T path (550°C, 6–9 kbar; epidote–amphibolite and greenschist facies conditions) during the main Pan-African collision event (620–580 Ma). The composition of magmas forming the Cryogenian Amalaoulaou arc and the processes leading to intra-arc differentiation are strikingly comparable to those observed in the deep section of exposed Mezosoic oceanic arcs, namely the Kohistan and Talkeetna complex. This evolution of the Amalaoulaou oceanic arc and its accretion towards the West African craton belong to the life and closure of the Pharusian Ocean that eventually led to the formation of the Greater Gondwana supercontinent, a similar story having occurred on the other side of the Sahara with the Mozambique Ocean.