A dated cave sediment record of Pleistocene transgressions on Berry Head,Southwest England

Berry Head, a limestone headland in Torbay, southwest England, exhibits a series of subaerial marine erosion platforms and raised beaches spanning an altitude range of 97 m. Solution caves on the headland show preferred horizontal development at elevations that are correlated with the marine erosion platforms, and developed in a marine/freshwater mixing zone whose position was controlled by high sea-level still-stands. Corbridge Cave in Berry Head Quarry lies below the raised beaches in Torbay, and contains evidence of three marine transgressions in the form of fine-grained marine ponding deposits with a marine microfauna. Uranium-series dating of intercalated speleothems indicates that a transgression during oxygen isotope stage 5e reached an elevation of 5.8 m OD, while an earlier transgression (probably during stage 7) reached at least 7.2 m OD. These findings are used to constrain possible interpretations of the aminostratigraphy of raised beaches in Southwest Britain, and a correlation of the Unnamed Stage of Bowen, Sykes, Reeves, Miller, Andrews, Brew and Hare with oxygen isotope stage 5e is proposed.     

Regolith Surface Reflectance: A New Attempt to Model

The reflectance coefficient of the regolith layer of celestial bodies has been studied in relation to the physical properties of regolith particles (size, refractive index, and packing density) on the basis of an accurate numerical radiative-transfer algorithm for a semi-infinite flat layer. Using the geometric-optics approximation, we have found that a shape mixture of randomly oriented spheroids can successfully model the single-scattering phase function of independent soil grains. In order to take into account the effect of packing density in a regolith layer, the concept of the so-called static structure factor was used. The main effect of increasing packing density is to suppress the forward-scattering peak of the phase function and to increase the albedo of the reflecting surface. We also investigated the influence of fine dust on the reflected light. An addition of small particles not only increases the surface albedo, but also changes the brightness profile and enhances the backscattering. Although the problem of unique solution, which is inherent in the retrieval of the properties of a medium from the measurements of the intensity of light scattered by this media, cannot be removed in the proposed model, the procedure used here, in contrast to widely used approximations, allows us to fit observational data with a set of real characteristics of the regolith. Semiempirical approaches are able to fit the measurements well with a small number of free parameters, but they do not explicitly contain crucial physical characteristics of the regolith such as grain sizes or the refractive index. We compared the numerical solution of the radiative-transfer equation with the Hapke approximation, which is most often used by investigators. The errors introduced by the Hapke model are small only for near-isotropic scattering by isolated particles. However, independent regolith grains are known to scatter light mainly in the forward direction.     

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.     

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.     

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.     

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.     

Predicted interplanetary distribution of Lyman-σ intensity and polarization

Our current understanding of the spatial and temporal distribution of interplan etary neutral hydrogen is currently limited to a comparison of Lyman-σ photometric data with predictions of the solar backscattered radiation using theoretical models. In this paper, how the uncertainties in current model calculations could be reduced through the future use of polarization measurements made from interplanetary spacecraft is investigated. In particular, inquiry into how a mapping of the degree of linear polarization made from a spacecraft at various locations in the Solar System can improve knowledge of the interstellar wind parameters, number density, temperature, and velocity, is made. A polarization measurement can, in principle, be made with very high precision. In this regard, being a relative quantity, a polarization measurement can be made independent of instrumental calibration and long-term sensitivity degredation. Furthermore, the sky distribution of both intensity and polarization has been calculated using a variety of models for the neutral hydrogen. It is found that the polarization distribution over the sky is quite different from that of the intensity distribution. It is also showed that the maximum degree of polarization of the Laymam-σ line increases with heliocentric distance of the spacecraft, varying from 0 up to ～ 18% at 20 AU.     

The determination and analysis of the orbit of NIMBUS 1 ROCKET, 1964-52B: Part 1: The orbits

Orbital parameters of the satellite NIMBUS 1 ROCKET (1964-52B) have been determined at 285 epochs using about 14,500 optical and radar observations during the last 2233 days of its life. 1964-52B was in an orbit inclined at about 98.7° to the Equator and it decayed on 13 August 1974. The orbital elements were found using the RAE orbit refinement program PROP 6 giving average orbital accuracies in perigee height and inclination at 110 m and 0.002° respectively.