DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.     

The role of space related non-governmental organizations (NGOs) in capacity building

Non-governmental organizations (NGOs) play a unique role in international affairs, providing access to resources, expertise, and assistance to supplement State resources. Sometimes the diplomatic skills and unofficial access of NGOs to policymakers through Track Two diplomacy can move a previously stalled critical issue forward and assist policymakers from different countries to find common ground outside official channels. Because they work outside of official channels, they are not bound by State policy that may inhibit negotiations between States. Some also have a convening power that sometimes makes it possible for State representatives to meet discipline experts and each other for informal discussions on issues of mutual interest. Finally, NGOs can draw attention to issues that may be overlooked or avoided by State organizations.     

Computational fluid dynamics and material processing (buoyancy and thermocapillary convections in melts during directional crystallization)

Modelisation and solution of heat and mass transfer problems relevant for material processing are generally hard to be handled, as they often involve 3D unsteady flows, viscous mixtures, phase changes, moving liquid-solid fronts, deforming liquid-gas interfaces, etc.… For space applications, material processing benefits of reduced buoyancy convection but can be faced to a strongly increased complexity due to variable g, mainly in manned flight.

Computational techniques used to analyse fluid motions in material processing, accounting for free surface, crystallization front and bulk convection in melt, are reviewed with emphasis to directional crystallization. Hydrodynamics stability and bifurcation analysis are shown to be useful complementary tools for correlating data, and for a better understanding of the physical laws. This last point will be illustrated in the case of the onset of oscillations in metallic melts.     

Relative amplitude of the total electron content variations depending on geomagnetic activity

We developed a method of estimation of a relative amplitude dI/I of the total electron content (TEC) variations in the ionosphere as deduced from the data of the global GPS receivers network. To obtain statistically significant results we picked out three latitudinal belts provided in the Internet by the maximum number of GPS sites. They are high-latitudinal belt (50–80°N, 200–300°E; 59 sites), mid latitude belt (20–50°N, 200–300°E; 817 sites), and equatorial belt (±20°N, 0–360°E; 76 sites). The results of the analysis of the diurnal and latitudinal dependencies of dI/I and dI/I distribution probability for 52 days with different levels of geomagnetic activity are presented. It was found that on average the relative amplitude of the TEC variations varies within the range 0–10% proportionally to the value of the K_p geomagnetic index. In quiet conditions the relative amplitude dI/I of the TEC variations at night significantly exceeds the daytime relative amplitude. At high levels of magnetic field disturbances, the geomagnetic control of the amplitude of TEC variations at high and middle latitudes is much more significant than the regular diurnal variations. At the equatorial belt, on average, the amplitude of TEC variations in quiet and disturbed periods almost does not differ. The obtained results may be useful for development of the theory of ionospheric irregularities.     

Investigation into impact of tropical cyclones on the ionosphere using GPS sounding and NCEP/NCAR Reanalysis data

Ionospheric response to tropical cyclones (TCs) was estimated experimentally on the example of three powerful cyclones – KATRINA (23–31 August 2005), RITA (18–26 September 2005), and WILMA (15–25 October 2005). These TCs were active near the USA Atlantic coast. Investigation was based on Total Electron Content (TEC) data from the international network of two-frequency ground-based GPS receivers and the NCEP/NCAR Reanalysis data. We studied the spatial–temporal dynamics of wave TEC disturbances over two periods of ranges (02–20 min and 20–60 min). To select the ionospheric disturbances which were most likely to be associated with the cyclones, maps of TEC disturbances were compared with those of meteorological parameters.     

GTD Terrain Reflection Model Applied to ILS Glide Scope

The capability of calculating the reflection of electromagnetic signals from uneven terrain has many applications. One of these is the determination of instrument landing system (ILS) glide slope performance. For this application the wavelength is approximately 1 m, incidence angles are usually near grazing, and the fields are horizontally polarized, so that gross irregularities such as dropoffs and hills are more important than surface roughness. Past approaches used to calculate the ground reflections for this application have been three-dimensional physical optics models which were very cumbersome and time consuming and which neglected important diffraction and shadowing phenomenon; a two-dimensional physical optics model which was faster than the three-dimensional models but ignored many shadowing and transverse terrain variation effects; and a half-plane diffraction model which is applicable only to a specified type of terrain geometry. In this paper a terrain reflection model based on the geometrical theory of diffraction (GTD) is described which can accommodate any piecewise linear terrain profile, requires less computer time than the physical optics models, is capable of including transverse terrain effects, and determines the reflected fields with all important diffraction and blockage effects included.     

Type Pi 1–2 magnetic field pulsations

The relationships of type Pi (broadband) pulsations to various other substorm-related phenomena are reviewed. Several of the more popular mechanisms for the origin of Pi activity are discussed in the light of the observations. There is only one mechanism in sight that tentatively accounts for observed characteristics of Pi 1–2 activity at auroral oval and polar cap latitudes and that is the three-dimensional current loop mechanism. If two or more mechanisms are involved in the generation of Pi noise, then it is possible that the garden-hose overstability and/or a drift Alfvén wave mechanism operating in the plasma sheet contribute to the observed pulsations.The common feature of all Pi 1–2 events is not the presence of temporal precipitation pulsations but the presence of an E-region, suggesting that enhanced conductivity and E-region currents are required. Pi activity appears to be closely related to unsteady convection in progress. Pi data promise to provide useful information on convection and field-aligned and ionospheric currents.     

An optimum filter for direct A/D conversion

The author describes a design procedure to find the optimum finite-impulse-response (FIR) filter coefficients when the characteristics of the bandpass filter preceding the analog-to-digital (A/D) converter are known. A simulation of direct A/D conversion performance shows the improvement over the filter coefficients designed to minimize the maximum error over a specified band of frequencies. DC cancellation may be added to the process with very little change in the coefficients of the optimum FIR filter or its performance.<>     

The magnetotail and substorms

The tail plays a very active and important role in substorms. Magnetic flux eroded from the dayside magnetosphere is stored here. As more and more flux is transported to the magnetotail and stored, the boundary of the tail flares more, the field strength in the tail increases, and the currents strengthen and move closer to the Earth. Further, the plasma sheet thins and the magnetic flux crossing the neutral sheet lessens. At the onset of the expansion phase, the stored magnetic flux is returned from the tail and energy is deposited in the magnetosphere and ionosphere. During the expansion phase of isolated substorms, the flaring angle and the lobe field strength decrease, the plasma sheet thickens and more magnetic flux crosses the neutral sheet.In this review, we discuss the experimental evidence for these processes and present a phenomenological or qualitative model of the substorm sequence. In this model, the flux transport is driven by the merging of the magnetospheric and interplanetary magnetic fields. During the growth phase of substorms the merging rate on the dayside magnetosphere exceeds the reconnection rate in the neutral sheet. In order to remove the oversupply of magnetic flux in the tail, a neutral point forms in the near earth portion of the tail. If the new reconnection rate exceeds the dayside merging rate, then an isolated substorm results. However, a situation can occur in which dayside merging and tail reconnection are in equilibrium. The observed polar cap electric field and its correlation with the interplanetary magnetic field is found to be in accord with open magnetospheric models.