Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

The Effect of Atmospheric Drag on Satellite Orbits During the Bastille Day Event

Geomagnetic storms driven by solar eruptions are known to have significant effects on the total density of the upper atmosphere in the altitude range 250–1000 km. This in turn causes a measurable effect on the orbits of resident space objects in this altitude range. We analyzed a sample of these orbits, both from sensor data and from orbital element sets, during the period surrounding the 14 July 2000 solar activity. We present information concerning the effects of this event on the orbits of resident space objects and how well accepted atmospheric models were able to represent it. As part of this analysis, we describe a technique for extracting atmospheric density information from orbital element sets. On daily time scales, the effect of geomagnetic activity appears to be more important than that of prompt radiation. However, the limitations in time and amplitude quantization of the accepted solar indices are evident. A limited comparison is also made with previous solar storm events.     

Long-term seismicity of Istanbul and of the Marmara Sea region

We examine the long-term seismicity of the Marmara Sea region in Turkey over a period of twenty centuries to attempt to answer two questions: (1) how well recorded in history is the earthquake activity of this region? and (2) does seismicity over this long period of time differ from that over the present century? We study this densely-populated and fast-developing region because it is one of the most tectonically active regions on the continents that appears today to be relatively quiescent.     

A viable mechanism for the production of energy in active galactic nuclei

It is suggested that a collapsing supermassive object, which acts as an ultra-high energy particle accelerator, is the precursor of an active galactic nucleus and that the gravitational energy released during the collapse of the object is locked in the quark-gluon plasma permeated by leptons into which the entire matter in the core of the object is converted as a result of the collapse. It is also pointed out that the collapse of the object to a space-time singularity is inhibited by Pauli's exclusion principle as well as by Heisenberg's uncertainty principle and that the object explodes, before it could collapse to a singularity, thereby releasing the enormous amount of energy locked in the quark-gluon plasma.     

Challenging a Paradigm: Do We Need Active and Inactive Areas to Account for Near-Nuclear Jet Activity?

We briefly describe an advanced 3D gas dynamical model developed for the simulation of theenvironment of active cometary nuclei. The model canhandle realistic nucleus shapes and alternative physical models for the gas and dust production mechanism.The inner gas coma structure is computed by solving self-consistently(a) near to the surface the Boltzman Equation(b) outside of it, Euler or Navier-Stokes equations.The dust distribution is computed from multifluid ``zero-temperature' Euler equations,extrapolated with the help of a Keplerian fountain model.The evolution of the coma during the nucleus orbital and spin motion,is computed as a succession of quasi-steady solutions. Earlier versions of the model using simple,``paedagogic' nuclei have demonstrated that the surface orographyand the surface inhomogeneity contribute similarly to structuring the near-nucleusgas and dust coma,casting a shadow on the automatic attribution of such structures to ``active areas'.The model was recently applied to comet P/Halley, for whichthe nucleus shape is available. In the companion paper of this volume,we show that most near-nucleus dust structuresobserved during the 1986 Halley flybys are reproduced, assuming that the nucleus is strictly homogeneous. Here, we investigate the effect of shape perturbations and homogeneityperturbations. We show that the near nucleus gas coma structure is robust vis-a-vissuch effects. In particular, a random distribution of active and inactive areaswould not affect considerably this structure, suggesting that such areas,even if present, could not be easily identified on images of the coma.     

Local potential analysis of MHD instability

In many astrophysical problems, the study of the stability of an atmosphere in the presence of a magnetic field is of importance. In most cases the MHD instabilities of atmospheres are studied by energy principle of Bernsteinet al. (1958). In this paper, a general method for studying the stability of a system subject to MHD equations of conditions has been proposed. This is based on the local potential concept put forward by Glansdorff and Prigogine (1964). The scheme for securing stability criteria has been demonstrated in two particular cases.     

Optimum merger time-scales and binary evolution

A study of galaxy mergers, on the basis of the collisional theory, using galaxy models without halos and considering the evolution of the proginator galaxies only from a time when the gravitational interaction between them is physically significant, indicates that most of the mergers are affected in 2 to 3 orbital periods for progenitors of comparable mass: shorter and longer time-scales being underabundant. These results have a bearing on the evolution of binary galaxies; indicating that once the relative orbit of a binary is circularized, the components will merge during the subsequent orbit or the next one (in a time-scale ~ 10⁸ years). These results are also indicative of the fact that binary evolution is very likely to cause a gradual evolution of the fundamental plane occupied by paired ellipticals from that of isolated ellipticals. After the merger, the remnant is very likely to define a fundamental plane with a slightly different slope.     

A Critical review of theoretical models of negatively polarized light scattered by atmosphereless solar system bodies

About a dozen physical mechanisms and models aspire to explain the negative polarization of light scattered by atmosphereless celestial bodies. This is too large a number for the reliable interpretation of observational data. Through a comparative analysis of the models, our main goal is to answer the question: Does any one model have an advantage over the others? Our analysis is based on new laboratory polarimetric and photometric data as well as on theoretical results. We show that the widely used models due to Hopfield and Wolff cannot realistically explain the phase-angle dependence of the degree of polarization observed at small phase angles. The so-called interference or coherent backscattering mechanism is the most promising model. Models based on that mechanism use well-defined physical parameters to explain both negative polarization and the opposition effect. They are supported by laboratory experiments, particularly those showing enhancement of negative polarization with decreasing particle size down to the wavelength of light. According to the interference mechanism, pronounced negative branches of polarization, like those of C-class asteroids, may indicate a high degree of optical inhomogeneity of light-scattering surfaces at small scales. The mechanism also seems appropriate for treating the negative polarization and opposition effects of cometary dust comae, planetary rings, and the zodiacal light.     

A laboratory test of the soil chemical submodels of two models of catchment acidification

Parameters for ion exchange selectivity and aluminium hydroxide dissolution in the soil chemical submodels used in applications of the Birkenes model and of MAGIC are compared and several discrepancies identified for organic soils. A laboratory column simulation of the soil chemical submodels is proposed and applied to soils from the Loch Dee area in Galloway. Experimental results were well predicted by a simplified version of MAGIC, with ion exchange selectivity parameters similar to those used in a previous simulation of one subcatchment of Loch Dee. The aluminium hydroxide dissolution parameter used previously was found to be too low for the organic soil materials, where a value of 106 predicted the experimental results more closely. The model developed also included a simple silicate weathering reaction to release base cations into the system. It is concluded that such simple laboratory simulations are useful for independent calibration of the soil chemical submodel of catchment models.