Vening Meinesz inverse isostatic problem with local and global Bouguer anomalies

Inverse problems in isostasy will consist in making the isostatic anomalies to be zero under a certain isostatic hypothesis. In the case of the Vening Meinesz isostatic hypothesis, the density contrast is constant, while the Moho depth (depth of the Mohorovičić discontinuity) is variable. Hence, the Vening Meinesz inverse isostatic problem aims to determine a suitable variable Moho depth for a prescribed constant density contrast. The main idea is easy but the theoretical analysis is somewhat difficult. Moreover, the practical determination of the variable Moho depths based on the Vening Meinesz inverse problem is a laborious and time-consuming task. The formulas used for computing the inverse Vening Meinesz Moho depths are derived. The computational tricks essentially needed for computing the inverse Vening Meinesz Moho depths from a set of local and global Bouguer anomalies are described. The Moho depths for a test area are computed based on the inverse Vening Meinesz isostatic problem. These Moho depths fit the Moho depths derived from seismic observations with a good accuracy, in which the parameters used for the fitting agree well with those determined geophysically. Received: 4 February 1999 / Accepted: 4 October 1999     

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.     

Atmospheric, oceanic and hydrological contributions to seasonal variations in length of day

 The annual and semiannual residuals derived in the axial angular momentum budget of the solid Earth–atmosphere system reflect significant signals. They must be caused by further excitation sources. Since, in particular, the contribution for the wind term from the atmospheric layer between the 10 and 0.3 hPa levels to the seasonal variations in length of day (LOD) is still missing, it is necessary to extend the top level into the upper stratosphere up to 0.3 hPa. Under the conservation of the total angular momentum of the entire Earth, variations in the oceanic angular momentum (OAM) and the hydrological angular momentum (HAM) are further significant excitation sources at seasonal time scales. Focusing on other contributions to the Earth's axial angular momentum budget, the following data are used in this study: axial atmospheric angular momentum (AAM) data derived for the 10–0.3 hPa layer from 1991 to 1997 for computing the missing wind effects; axial OAM functions as generated by oceanic general circulation models (GCMs), namely for the ECHAM3 and the MICOM models, available from 1975 to 1994 and from 1992 to 1994, respectively, for computing the oceanic contributions to LOD changes, and, concerning the HAM variations, the seasonal estimates of the hydrological contribution as derived by Chao and O'Connor [(1988) Geophys J 94: 263–270]. Using vector representation, it is shown that the vectors achieve a close balance in the global axial angular momentum budget within the estimated uncertainties of the momentum quantities on seasonal time scales. Received: 6 April 2000 / Accepted: 13 December 2000     

The stochastic estimation of satellite clock correction information in WADGPS

1 IntroductionChnventional DGPS is limited by the range overwhich the differential corrections are valid due tothe rapid decorrelation of the error sources with in-creasing distance from the reference station to user.In wide area differential GPS (WWPS) errorsources in GPS measurements are medeled sepa-rately,on the basis of a limited number of referencestaitOns, to overcome this drawback. The main er-ror sources are regarded as broadcast ephemeris er-ror,atmospheric refraction and satel…     

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.     

Neutron stars in the multidimensional einstein theory of gravitation

The radius, mass, total number of baryons, and other parameters of static, spherically symmetric, superdense stars are calculated. A model with one Ricci-flat inner space of arbitrary dimensionality and the approximation p₁=?0.5ε + ap for additional components of the energy — momentum tensor are used (ε and ρ are the total energy density and the pressure of the stellar matter and a is a fitting parameter). In the case of white dwarfs, the results of the multidimensional theory do not depend on the dimensionality D of space-time for ?10 ? a ? 10 and coincide with the analogous data of the general theory of relativity (GTR). For neutron stars there is a dependence on D and a. For D>4, in particular, the greatest mass M_max of a neutron star as a function of a has a maximum at 3<a(D) ? 4, which exceeds the greatest mass M _max ⁰ =2.14 M_⊙ in the GTR. A comparison of theoretical results with observational data determines the allowable values of a. Data for PSR 1913 + 16 lead to 0.2 ≤ a ≤ 9.2 in the case of D=26, while the results of [P. C. Joss and S. A. Rappaport, Annu. Rev. Astron. Astrophys.,22, 537 (1984)] lead to the stricter limits 1 ≤ a ≤ 7.4.     

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.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

Acid mine drainage at Cerro Rico de Potosí I: unabated high-strength discharges reflect a five century legacy of mining

Intensive mining and processing of Ag, Sn, Pb and Zn ores have occurred in various locations within and around the city of Potosí, Bolivia since 1545. Surface and subsurface waters, stream sediments and soils are contaminated with various ecotoxic metals in the headwaters of the economically vital, yet highly impacted, upper Rio Pilcomayo watershed. Previous studies have documented downstream trace metal contamination, however, not addressed specific sources. The AMD discharges identified in this study help link downstream pollution to primary origins. The majority of AMD would be considered high-strength due to metal concentrations and acidity orders of magnitude greater than typical AMD. Discharges from both operating and abandoned portals as well as tailings-related deposits displayed a high degree of heterogeneity with total metal concentrations ranging from 0.11 to 7480, <0.022 to 889, <0.0006 to 65.3, <0.001 to 310, 0.12 to 72,100, 0.3 to 402, <0.012 to 34.8, and 0.24 to 19,600 mg/L of Al, As, Cd, Cu, Fe, Mn, Pb and Zn, respectively. Net acidity and pH ranged from −10 to 246,000 mg/L as CaCO₃ equivalent and 0.90–6.94 standard units, respectively. Data were gathered during two sampling events centered around the most extreme periods of the dry and wet seasons of one water-year. Loadings to local streams were marginally greater for most metals in the wet season. If observed loadings are historically representative, Cerro Rico AMD has contributed thousands of tonnes of ecotoxic metals to the upper Rio Pilcomayo over the last five centuries. Metal and hydrogen ion concentrations in the majority of AMD sampled were several orders of magnitude above discharge limits set by the Bolivian government, yet no action has historically or contemporarily been taken.