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.     

A 59-second period in the very high energy gamma-ray emission from the Geminga pulsar

An analysis of our observations of the Geminga object with the GT-48 ground-based gamma-ray telescope has shown that its very-high-energy gamma-ray flux is modulated with a 59-s period. The 59-s period and its time derivative previously inferred from satellite data have been confirmed. According to our data, the period was 61.94 s in 1997 at MSD=50573. The statistical significance of this result is (1?4.5)×10^?4.     

Pulsar reflections within the Crab Nebula

Between 1997 August and October, the radio pulses from the Crab pulsar were followed by discrete moving echoes, which appear to be reflections from part of an ionized shell in the outer part of the Crab Nebula, crossing the line of sight to pulsar. Similar events have now been recognized in recordings from the past 30 yr, and it seems that the Nebula must contain a large number of ionized shell-like surfaces on a much finer scale than recognized hitherto.     

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.     

On the mathematical description of fine structure under the conditions of weak variability

The idea of the weak variability of the fields is introduced, which means that when they are recorded by some measuring instrument (MI), the output signal obtains increments of the order of one/several quanta of sensitivity at the discreteness interval. At small scales (microstructure), this field is described in terms of random discrete values which reflect the properties of both the field measured and the MI. Theoretical relationships which define the probabilities of increments at the compound intervals in terms of the probabilities of increments at the sampling intervals, are derived for the case of independent random increments at the sampling intervals of discreteness. The validity of this model is illustrated by respective computations using data of microstructure measurements of the temperature profile in the main thermocline in the Sargasso Sea.Translated by Mikhail M. Trufanov.     

The interannual leading modes of the extratropical variability in the Southern Hemisphere simulated by the ECHAM-4 atmospheric model

 An ensemble of twenty-three 14-year experiments conducted with the ECHAM-4 GCM has been examined to test the model's capability to simulate the principal modes of interannual variability. The integrations were performed under specified monthly SST between 1979–1993. The analysis was focused on the Southern Hemisphere (SH) extratropics. Empirical orthogonal functions analysis (EOF) using seasonal anomaly fields has been performed to isolate the principal modes that dominate the southern extratropical variability at the interannual time scale. Leading patterns of 500 hPa geopotential height (z500) have been compared with those estimated from the ECMWF re-analysis dataset. The model is able to adequately reproduce the spatial pattern of the annular mode, but it represents the temporal variations of the oscillation less satisfactorily. The model simulation of the Pacific South American (PSA) pattern is better, both in the shape of the pattern and in the temporal evolution. To verify if the capability of the model to adequately simulate the temporal oscillation of the propagating patterns is related to the increased influence of the tropical external forcing, covarying SST-atmospheric modes have been identified by singular value decomposition (SVD). In winter (July-August-September, JAS) the tropical SST variability is highly correlated with the ENSO mode. In summer (January-February-March, JFM) the strength of the teleconnections is related to strong westerly anomalies, disrupted by a meridional out of phase relation near to South America. The large size of the ensemble was exploited by comparing the time-varying model spread and degrees of freedom of the simulated extratropical circulation. Results show that when the extratropical circulation has a few degrees of freedom, the reproducibility is relatively low and the ensemble is governed by a fairly robust zonally symmetric structure of dispersion. Received: 9 May 2000 / Accepted: 30 January 2001     

Polarimetric observations of pulsating variable stars

Results are presented for polarimetric observations of 17 red giants and supergiants, of which nine are long-period Mira variables, five are semiregular variables (SR), and three are slowly fluctuating variables (Lb). Light polarization is detected for eight stars, seven of them for the first time.Translated fromAstrofizika, Vol. 39, No. 3, pp. 385–391, July–September, 1996.     

Contribution of chromospheric features to UV irradiance variability from spatially-resolved Ca II K spectroheliograms

We have digitized the Ca ii K spectroheliograms, observed at the National Solar Observatory at Sacramento Peak, for the period 1980 (maximum of solar cycle 21), 1985 (minimum of solar cycle 21), 1987 (beginning of the ascending phase of solar cycle 22), 1988 and 1989 (ascending phase and maximum of solar cycle 22), and 1992 (declining phase of solar cycle 22). A new method for analyzing the K spectroheliograms has been developed and applied to the K images for the time interval of 1992. Using histograms of intensity, we have segregated and measured the cumulative intensity and area of various chromospheric features like the plages, magnetic network and intranetwork elements. Also, the full width at half maximum (FWHM) derived from the histograms has been introduced as a new index for describing the chromospheric activity in the K-line. The full-disk intensity (spatial K index) has been derived from spatially-resolved K images and compared to the spectral K index derived from the line profiles for the full disk. Both the spatial K index and FWHM have been compared to the UV irradiance measured in the Mg ii h and k lines by the NOAA9 satellite and found that they are highly correlated with the Mg ii h and k c/w ratio.NRC Resident Research Associate, on leave from Indian Institute of Astrophysics, Bangalore 560034, India.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.