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     

Precise Real-Time Positioning in WADGPS Networks

In recent years the importance of real-time positioning and navigation with the Global Positioning System (GPS) has grown rapidly. Starting from the establishment of differential GPS (DGPS) reference stations for marine and land navigation, new users and applications have emerged that resulted in a high demand for the establishment of a high-density network of reference stations around the world. Many countries have established their own DGPS service, which is either governmentally or commercially owned. These services are referred to as Local Area DGPS Systems (LADGPS). However, the costs for the establishment and maintenance of a dense network of reference stations are very high. Therefore Wide Area DGPS Systems (WADGPS) are being developed to overcome the main drawbacks of LADGPS. In this case, only a few reference stations are used to cover a large area, such s a continent like Europe. To achieve high positioning accuracies, real-time modeling of the main error sources for long-range baselines is required as errors in the satellite orbit and ionospheric refraction do not cancel entirely in double differencing. In this article, a real-time correction model based on the Kalman filter for WADGPS and networked LADGPS services is discussed and results of field tests in a WADGPS network in Europe are presented. ? 2000 John Wiley & Sons, Inc.     

Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt

 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001     

Physico-chemical phenomena in comets—III: The continuum of comet Burnham (1960 II)

The new model of the cometary head proposed in papers I and II is developed and applied to comet Burnham. It takes into account the likely existence of a halo of large icy particles surrounding the nucleus. These particles are steadily stripped from the nucleus by evaporating gases. Their terminal velocity and their rate of evaporation set the size of the halo. The existence of the icy halo influences in two ways the photometric characteristics of the coma. This paper establishes the photometric shape of the continuum as reflected by the icy grains, and compares it to the observed continuum of comet Burnham. Paper IV will compare the predictions of the model with the photometric profile of the molecular emission bands of C₂, in the same comet.     

Past,present and future vegetation-cloud feedbacks in the Amazon Basin

To begin exploring the underlying mechanisms that couple vegetation to cloud formation processes, we derive the lifting condensation level (LCL) to estimate cumulus cloud base height. Using a fully coupled land–ocean–atmosphere general circulation model (HadCM3LC), we investigate Amazonian forest feedbacks on cloud formation over three geological periods; modern-day (a.d. 1970–1990), the last glacial maximum (LGM; 21 kya), and under a future climate scenario (IS92a; a.d. 2070–2090). Results indicate that for both past and future climate scenarios, LCL is higher relative to modern-day. Statistical analyses indicate that the 800 m increase in LCL during the LGM is related primarily to the drier atmosphere promoted by lower tropical sea surface temperatures. In contrast, the predicted 1,000 m increase in LCL in the future scenario is the result of a large increase in surface temperature and reduced vegetation cover.     

柴北缘鱼卡地区中侏罗统石门沟组含煤层段沉积有机相分析

柴达木盆地北缘鱼卡地区中侏罗统石门沟组含煤段中发育有油页岩、煤、碳质泥岩和泥岩等富有机质细粒沉积.为了研究其沉积有机相的类型及煤和油页岩形成的控制因素,本文通过岩心观察、工业分析结合有机地球化学测试分析等方法将含煤段细粒沉积物划分为了类型A(三角洲平原沉积环境)、B(三角洲前缘-浅湖沉积环境)、C1(湖沼非油页岩亚相)及C2(湖沼油页岩亚相)四种类型,其中类型C1沉积物中发育煤,类型C2中发育油页岩,且煤较油页岩具有更高的含油率、水分、挥发分以及发热量值.有机质类型特征方面,沉积有机相类型A、B及C1沉积物有机质类型均为Ⅱ2-Ⅲ型、有机质来源均以陆源和混合来源两种有机质来源为主,而类型C2中主要为Ⅱ2型,以混合有机质来源为主;类型A及C2其沉积物有机质均处于未成熟阶段,而类型B和C1中有机质均处于未成熟-低成熟阶段;有机质保存条件方面,类型A沉积物主要形成于缺氧的淡水环境,类型B主要形成于贫氧-缺氧的淡水-半咸水环境,类型C1形成于贫氧-缺氧的淡水-半咸水环境,而类型C2则主要形成于缺氧的淡水-半咸水环境.其中类型B较类型A,类型C2较类型C1,其沉积物均形成于更为还原且盐度更高的水体环境中.石门沟组含煤段是煤和油页岩的形成层位,稳定的沉积环境、丰富的湖泊有机质来源、良好的保存条件及较少的陆源碎屑的稀释共同促进了类型C2中油页岩的形成,而丰富的陆源植物供给及良好的保存条件则促进了类型C1中煤层的形成.     

Urban seismic stations: soil–structure interaction assessment by spectral ratio analyses

In this work we present and discuss the results of ambient seismic noise analyses computed at four sites where seismic stations, managed by the INGV (Italian Institute for Geophysics and Vulcanology) and the DPC (Italian Department of Civil Protection), are installed inside buildings. The experiments were performed considering different types of installation: sensor located at the bottom of a school, directly installed on rock (case 1); sensor located at the bottom of a medieval fortress, built on an isolate hill, directly installed on rock (case 2); sensor installed on the foundations of a medieval fortress, built on an isolate hill (case 3); sensor installed on the foundations of a school, built on alluvial deposits (case 4). Since recent works proposed the use of spectral ratio techniques to study the dynamic characterization of buildings, ambient seismic-noise measurements were performed for each site close to the stations (at the base of the structures), at the top of the structures and outside the buildings. In order to check the source of vibrations both horizontal to vertical spectral ratio (HVNR) and standard spectral ratio (SSR) techniques were applied. For all stations the results from ambient seismic noise were compared to those obtained from earthquakes (HVSR). In order to detect preferential directions of amplification, for each site average HVNRs and HVSRs were computed considering one azimuth for each set of 5°. We obtain different results for different types of installation: in cases 1 and 2, where the sensors are directly installed on rock, the vibrations of the structure do not affect the noise measures performed close to the stations, which show flat HVNR in the whole frequency range: in both cases the eigenfrequency of the building is given by the HVNR calculated from the measures performed at the top of the structure. In cases 3 and 4, where the sensors are installed on the foundations of the considered structures, both the amplification peaks between 5 and 9 Hz (case 3) and between 5.5 and 7 Hz (case 4) include the contribution of the free oscillations of the buildings. In particular, in case 4, HVNRs performed outside building highlight possible soil–structure resonance effects in case of an earthquake.     

Statistical biases and errors inherent in photoclinometric surface slope estimation with natural light

Photoclinometry is the most common method used to obtain high-resolution topographic maps of planetary terrain. We derive the likelihood function of photoclinometric surface slope from (1) the probability distribution of the measured photon count of natural sunlight through a Charge-Coupled Device (CCD) including uncertainty due to camera shot noise, camera read noise, small-scale albedo fluctuation and atmospheric haze, and (2) a photometric model relating photocount to surface orientation. We then use classical estimation theory to determine the theoretically exact biases and errors inherent in photoclinometric surface slope and show when they may be approximated by asymptotic expressions for sufficiently high sample size. We show how small-scale albedo variability often dominates biases and errors, which may become an order of magnitude larger than surface slopes when surface reflectance has a weak dependence on surface tilt. We provide bounds on the minimum possible error of any unbiased photoclinometric surface slope estimate, and compute the sample sizes necessary to constrain errors within desired design thresholds.     

Thermodynamic controls on element partitioning between titanomagnetite and andesitic–dacitic silicate melts

Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO₂) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg²⁺, Mn²⁺, Co²⁺, Ni²⁺, Zn²⁺) and Cu²⁺/Cu⁺ with increasing fO₂ between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO₂ implies a reciprocal relationship between a(Fe²⁺O) and a(Fe³⁺O_1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO₂ in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO₂ in silicic, but less so in mafic bulk systems.