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.     

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     

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.     

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.