Structure in the EAS size spectrum: analysis of the CASA-MIA results

We present an analysis of the very recent data from the CASA-MIA extensive air shower array (M.A.K. Glasmacher, Ph.D. Thesis, University of Michigan (1998) in terms of a search for the structure in the shower size spectrum claimed by us from an analysis of “the world's data”. Our earlier claim is found to be supported to the extent that there is strong evidence for the existence of structure in the spectrum which cannot obviously be explained by the conventional Galactic Modulation Model. There is modest evidence for the structure being of “our” form and strong support for “our” mass composition when “corrected” to the interaction model advocated by us. None of the results are inconsistent with there having been a recent, nearby, single supernova.     

A “division” treatment of the sky background in CCD images

Considering the errors transmitted from all pre-treatments, I suggest a “division” method for removing the sky background including interference fringes in CCD images. Compared to the traditional method, which may be called a “subtraction” method, the “division” method is easier to carry out and gives more accurate brightness measurement especially for the fainter images.     

Rushing to equilibrium: a simple model for the collisional evolution of asteroids

A simple mathematical model for the evolution of a system of collisionally interacting bodies—such as the asteroid population—consists of two coupled, nonlinear, first-order differential equations for the abundances of “small” and “big” bodies. The model easily allows us to recover Dohnanyi's value for the exponent of the equilibrium mass distribution. Moreover, the model shows that any initial value for the ratio of “big” to “small” bodies rapidly relaxes to the equilibrium ratio, corresponding to the exponent, and that integrating the evolution equations backward in time—an attractive possibility to investigate the mass distribution of primordial planetesimals—leads to strong numerical instability.     

Response time of the high-latitude dayside ionosphere to sudden changes in the north-south component of the IMF

The time scale of the response of the high-latitude dayside ionospheric flow to changes in the North-South component of the interplanetary magnetic field (IMF) has been investigated by examining the time delays between corresponding sudden changes. Approximately 40 h of simultaneous IMF and ionospheric flow data have been examined, obtained by the AMPTE-UKS and -IRM spacecraft and the EISCAT “Polar” experiment, respectively, in which 20 corresponding sudden changes have been identified. Ten of these changes were associated with southward turnings of the IMF, and 10 with northward turnings. It has been found that the corresponding flow changes occurred simultaneously over the whole of the “Polar” field-of-view, extending more than 2° in invariant latitude, and that the ionospheric response delay following northward turnings is the same as that following southward turnings, though the form of the response is different in the two cases. The shortest response time, 5.5 ± 3.2 min, is found in the early- to mid-afternoon sector, increasing to 9.5 ± 3.0 min in the mid-morning sector, and to 9.5 ± 3.1 min near to dusk. These times represent the delays in the appearance of perturbed flows in the “Polar” field-of-view following the arrival of IMF changes at the subsolar magnetopause. Overall, the results agree very well with those derived by Etemadi et al. (1988, Planet. Space Sci. 36, 471) from a general cross-correlation analysis of the IMF B_z and “Polar” beam-swinging vector flow data.     

Simulation of precise orbit determination of lunar orbiters

Based on the ongoing Chinese lunar exploration mission, i.e. the “Chang'e 1” project, precise orbit determination of lunar orbiters is analyzed for the actual geographical distribution and observational accuracy of the Chinese united S-band (USB) observation and control network as well as the very long baseline interferometry (VLBI) tracking network. The observed data are first simulated, then solutions are found after including the effects of various error sources and finally compared. We use the space data analysis software package, GEODYN, developed at Goddard Space Flight Center, NASA, USA. The primary error source of the flight orbiting the moon is the lunar gravity field. Therefore, the (formal) error of JGL165P1, i.e. the model of the lunar gravity field with the highest accuracy at present, is first discussed. After simulating the data of ranging and velocity measurement as well as the VLBI data of the time delay and time delay rate, precise orbit determination is carried out when the error of the lunar gravity field is added in. When the orbit is determined, the method of reduced dynamics is adopted with the selection of appropriate empirical acceleration parameters to absorb the effect of errors in the lunar gravity field on the orbit determination. The results show that for lunar missions like the “Chang'e 1” project, that do not take the lunar gravity field as their main scientific objective, the method of reduced dynamics is a simple and effective means of improving the accuracy of the orbit determination of the lunar orbiters.     

Solving the G-dwarf problem with a three-component model of the chemical evolution of the Galaxy

A three-component chemical evolution model of the Galaxy is presented, which we believe will cast a new light on the G-dwarf problem. The model is based on a scenario of the Galaxy consisting of three major evolutionary phases: halo, thick disk and thin disk, separated by two short interludes of rapid collapse. The evolution of different stellar populations are treated separately, the combination of which yields an overall metallicity distribution function for the solar neighbourhood. We tested three different models using the same set of basic equations: the “prompt initial enrichment” (PIE) model, the “proportional yield” (PPY) model and the “collapse” (CLP) model. Best-fit parameters are derived. The results show that the different populations have remarkably different IMFs, while mass exchange has only minimally affected the chemical evolution in the solar vicinity, so that the solar vicinity can be regarded as a closed system, at least in the late stage of the Galactic evolution.     

An overview of sage I and II ozone measurements

The Stratospheric Aerosol and Gas Experiments (SAGE) I and II measure Mie, Rayleigh, and gaseous extinction profiles using the solar occultation technique. These global measurements yield ozone profiles with a vertical resolution of 1 km which have been routinely obtained for the periods from February 1979 to November 1981 (SAGE I) and October 1984 to the present (SAGE II). The long-term periodic behavior of the measured ozone is presented as well as case studies of the observed short-term spatial and temporal variability.

A linear regression shows annual, semi-annual, and quasi-biennial oscillation (QBO) features at various altitudes and latitudes which, in general, agree with past work. Also, ozone, aerosol, and water vapor data are described for the Antarctic springtime showing large variation relative to the vortex. Cross-sections in latitude and altitude and polar plots at various altitudes clearly delineate the ozone hole vertically and areally. Comparisons of vertical profiles are made from 1979 to 1988.

Although there is a three-year gap between the SAGE I and II measurements, the two data sets have been used to determine long-term changes in ozone. The intercomparison generally shows decreases in the upper stratosphere (25–50 km) of 4% or less from 1980 to 1986.     

Fluctuations in tidal (24-, 12-h) characteristics and oscillations (8-h-5-d) in the mesosphere and lower thermosphere (70–110 km): Saskatoon (52°N, 107°W), 1979–1981

An M.F. radar (2.2 MHz) operating at Saskatoon, Canada (52°N, 107°W) has been used to produce continuous wind data ( 80–110km) from September 1978–April 1981. The 24-, 12-h tidal oscillations reveal regular summer-winter transitions; in particular the semi-diurnal tide demonstrates strikingly regular and rapid equinoctial changes over the three years. The vernal and autumnal equinox changes are clearly different in morphology. Shorter term tidal fluctuations (2d τ 10d) are compared with mean winds and gravity wave amplitudes, as well as with satellite-derived stratospheric temperatures.

Spectral analysis of monthly data sets for 1980, from 90–105 km, reveal oscillations of the expected 8-, 12-, 24-h periods, but also of 10-, 16- and 2-, 5/6d. A modulation of the “2-d” wave by the 12-h wave is suggested as a possible cause of these surprisingly regular oscillations.     

The Cr and Sr sources in BOREXINO as tools for neutrino magnetic moment searches

We calculate the event rates induced by a ⁵¹Cr ν_e source and by a ⁹⁰Sr---⁹⁰Y source in BOREXINO through elastic scattering on electrons, assuming a nonzero neutrino magnetic moment μ_ν. We consider a source activity of about 2 MCi and estimate the solar ν (“source-off”) background for various oscillation scenarios. It is shown that values of μ_ν as low as 0.5 × 10⁻¹⁰μ_B ( 0.2 × 10⁻¹⁰μ_B) can be proved with the ⁵¹Cr source (⁹⁰Sr source) in about 100 days of data taking.     

The HNCO ring in the Sgr B2 region

The large-scale structure associated with the 2′N HNCO peak in Sgr B2 [Minh, Y.C., Haikala, L., Hjalmarson, Å., Irvine, W.M., 1998. ApJ 498, 261 (Paper I)] has been investigated. A ring-like morphology of the HNCO emission has been found; this structure may be colliding with the Principal Cloud of Sgr B2. This “HNCO Ring” appears to be centered at (l,b) = (0.7°,−0.07°), with a radius of 5 pc and a total mass of 1.0 × 10⁵ to 1.6 × 10⁶ M. The expansion velocity of the Ring is estimated to be 30–40 km s⁻¹, which gives an expansion time scale of 1.5 × 10⁵ year. The morphology suggests that collision between the Ring and the Principal Cloud may be triggering the massive star formation in the Sgr B2 cloud sequentially, with the latest star formation taking place at the 2′N position. The chemistry related to HNCO is not certain yet, but if it forms mainly via reaction with the evaporated OCN⁻ from icy grain mantles, the observed enhancement of the HNCO abundance can be understood as resulting from shocks associated with the collision between the Principal Cloud and the expanding HNCO Ring.     

Statistical analysis of correlation between the positions of quasars and the uppsala galaxies

The question of the association of quasars with galaxies is re-examined using 785 quasars at |b| 30, δ − 02°30′ in the new Hewitt and Burbidge Catalog (1980) and all galaxies in the “Uppsala General Catalogue of Galaxies” (Nilson, 1973). The results of the two-point cross-correlation function are presented and they show that there are on the average 0.153 ± 0.011 more galaxies within 10'.0 of a QSO than would be expected if the QSOs were distributed randomly. We find the marginal significance of the tendency for correlation between QSOs and galaxies to increase with increasing redshifts z or apparent magnitudes V. The nearest neighbor test has also been taken to analyse the same data as well as get further evidence for the apparent association between QSOs and galaxies.     

The dangerous border of the 5:2 mean motion resonance

After discovering that asteroid (2953) Vysheslavia, a Koronis family member, has a typical dynamical lifetime of the order of 10 Myr only against “falling” into the 5:2 Kirkwood gap and ending up in a hyperbolic orbit (Milani and Farinella, 1995. Icarus 114, 209–212), we decided to repeat the same study for two other asteroids, 1991 UA₂ and 1993 FR₅₈, lying also very close to the outer edge of the same resonance. The orbital elements of these bodies were not yet accurate enough and we appealed to observers, both professional and amateur, to obtain more astrometric data. Eight groups in three different countries carried out observations, and as a result the Minor Planet Center could derive updated, accurate orbits for both objects. Here we report on the preliminary results of our study. Whereas the long-term behavior of 1991 UA₂ looks fairly regular and stable, 1993 FR₅₈ exhibits a typical “stable chaos” behavior; nevertheless, it does not end up in the resonance within a time span of 50 Myr. On the other hand, its fictitious clones, placed just 3–4 × 10⁻³ AU closer to the resonance, undergo large irregular semimajor axis changes, fall into the resonance and escape on a timescale <1 Myr.     

High time resolution observations of solar radio emission at four frequencies: evidence for energy release in microflares

As evidence for energy release in microflares, high time resolution observations of solar radio emission obtained with our “synchronous observation system of solar radio radiation with high time resolution at four frequencies (1.42, 2.13, 2.84 and 4.26 GHz)” from December 1989 to April 1993 are presented in this paper. The observed events include weak ms spikes, “spike-likes”, fast pulsations as well as two kinds of newly discovered fast fine structures, i.e., microwave type III bursts and microwave patch-like structures. A statistical study of the duration of fast fine structures has been made and on its basis the various types of phenomena are illustrated with actual examples.     

Triple splitting with the F2-region of the ionosphere at high and mid-latitudes

The occurrence of the third (z-ray) component of the F2-trace on ionograms is investigated at high- and mid-latitudes. Diurnal variations show a systematic shift, with magnetic inclination, of the time of maximum occurrence. Seasonal variations show a winter maximum, and an inverse sunspot-cycle relationship exists. Maximum occurrence appears between a magnetic inclination of 70° and 80° with a fall-off either side.

Evidence is presented to suggest a z-ray association with “Spread-F” fronts, and a possible mechanism for the recording of the z-ray trace at the transmitter site is described. This involves longitudinal propagation of the o-mode at its normal reflection level, coupling at this point, and ultimate reflection for the z-ray mode as a result of sloping ionization contours belonging to “Spread-F” fronts extending in directions perpendicular to the magnetic meridian.

An association with V.L.F. emissions (“dawn-chorus”) is discussed.     

Catastrophic fragmentation as a stochastic process: sizes and shapes of fragments

It is rather difficult to understand theoretically and to analyse the experimental data concerning the mass and shape distributions of fragments created by catastrophic collisions. The fragmentation process is discussed as being a purely stochastical phenomenon; the size and shape distributions obtained in this way are compared with the results of laboratory experiments. The results are presented of some computer simulations of random volume fragmentation processes; they are a 3-D generalization of the numerical experiments described in Grady and Kipp (J. Appl. Phys. 58(3), 1210–1222, 1985). The features of the size distribution are discussed, comparing it with the expectations of the Mott-Linfoot and Grady-Kipp theories. In the literature the shape of fragments is defined in terms of the ratios B/A and C/A, where A, B, C are defined as the sizes of a fragment along three orthogonal axes. The definition of the shape of a fragment cannot be considered unique, since it is not obvious in which order to define the three axes when the fragments are not ellipsoidal. A few possible methods are introduced explicity, and the resulting differences are discussed. In this light, the shape results (the mean values and the distribution of the axial ratios) obtained in recent laboratory experiments are rediscussed and critically reviewed. For what concerns the stochastical modelling, the results of various simulations, corresponding to different assumptions regarding fragmentation properties are presented. It is shown that the main features of the shape distributions from laboratory experiments cannot be satisfactorily reproduced. Comparison of the results with the outcomes of the semiempirical fragmentation model by Paolicchi et al. (Icarus 121, 126–157, 1996), as well as with some results coming out from hydrodynamical simulations, shows how only a “global” and physical model, not a purely statistical one (neither global nor “local”), can afford to reproduce the observed data.     

Geminid Meteor shower activity 2003–2005 as observed by Gadanki radar

The distribution of meteor signals reflected from a backscatter radar is considered according to their duration. This duration time (T) is used to classify the meteor echoes and to calculate the mass index (S) of different meteoroids of shower plus sporadic background. Observational data on particle size distribution of the Geminid meteor shower are very scarce, particularly at low latitudes. In this paper the observational data from Gadanki radar (13.46°N, 79.18°E) have been used to determine the particle size distribution and the number density of meteoroids inside the stream of the Geminid meteor shower. The mean variation of meteor number density across the stream has been determined for three echo duration classes, T<0.4, T=0.4–1 and T>1 s. We are more interested in the appearance of echoes of various durations and therefore meteors of various masses in order to understand more on the filamentary structure of the stream. It is observed that the faint particle flux peaks earlier than the larger particles. We found a decreasing trend in the mass index values from the day of peak activity to the next observation days. The mass index profile was found to be U-shaped with a minimum value near the time of peak activity. The observed minimum s values are 1.64±0.05 and 1.65±0.04 in the years 2003 and 2005, respectively. The activity of the shower indicates the mass segregation of meteoroids inside the stream. Our results are best comparable with the “scissors” structure model of the meteoroid stream formation of Ryabova [2007. Mathematical modeling of the Geminid meteoroid stream. Mon. Not. R. Astron. Soc. 375, 1371–1380] by considering the asteroid 3200 Phaethon as an extinct comet.     

Observations of waves artificially stimulated by an electron beam inside a region with auroral precipitation

The POLAR 5 sounding rocket, launched from Andøya, Norway on 1 February, 1976 was of a “mother-daughter” configuration. An electron accelerator, mounted on the “daughter,” produced a pulsed electron beam with currents up to 130 mA and electron energies up to 10 keV. The waves, artificially stimulated by the injected electron beam, was studied using wave receivers, mounted on the “mother.” The receivers covered the frequency range from 0.1 kHz to 5 MHz.

In addition to the stimulated waves observed during beam injection, enhanced wave emissions were observed 10–20 ms after the end of beam injection. This emission seemed to be relatively independent of whether the electron beam is launched up or down along the geomagnetic field.

The high frequency emission observed after beam injection is found to be correlated with the passage through an auroral arc. In particular this emission is closely correlated with the flux of 4–5 keV auroral electrons.

The low frequency emissions observed after beam injection are concentrated in two bands below the lower hybrid frequency.

Different mechanisms for explaining the observed time delays between the beam injection and the observation of the emissions are discussed.     

Some aspects of Sporadic-E at mid-latitudes

Evidence is presented to suggest that the ƒ₀E_s value tends to be high at positions in the E_s-layer where trough slope-lines and crest slope-lines of F2-layer irregularities meet the E_s-layer. These slope-lines are drawn through the troughs and crests, respectively, of the characteristic kinks in the F2-layer ionization contours, which are associated with F2-layer irregularities.

A rotating-loop direction-finding system has allowed an estimation of the distribution of ionization, which gives rise to Sporadic-E echoes. Analysis of Sporadic-E occurrence, on two occasions, suggests that the reflecting surfaces are frontal in nature, the fronts having a separation from each other of some tens of kilometres, and probably existing as closed curves, with diameters of the order of several hundreds of kilometres. A possible association between these structures and the occurrence of the green line of the airglow, is discussed.

A distribution of ionization, which will give contours showing “clouds” of ionization at some frequencies and a ripple structure at other frequencies, is proposed, in an endeavour to explain the apparent dual nature of Sporadic-E occurrence.

The evidence seems to indicate that the mechanism operating at the E_s-layer level, producing the phenomenon of Sporadic-E, is the same as that which produces the F2-layer irregularities which are responsible for “Spread-F”.     

Space–time foam and cosmic-ray interactions

It has been proposed that propagation of cosmic-rays at extreme-energy may be sensitive to Lorentz-violating metric fluctuations (“foam”). We investigate the changes in interaction thresholds for cosmic-rays and gamma-rays interacting on the CMB and IR backgrounds, for a class of stochastic models of space–time foam. The strength of the foam is characterized by the factor (E/M_P)^a, where a is a phenomenological suppression parameter. We find that there exists a critical value of a (dependent on the particular reaction: a_crit3 for cosmic-rays, 1 for gamma-rays), below which the threshold energy can only be lowered, and above which the threshold energy may be raised, but at most by a factor of two. Thus, it does not appear possible in this class of models to extend cosmic-ray spectra significantly beyond their classical absorption energies. However, the lower thresholds resulting from foam may have signatures in the cosmic-ray spectrum. In the context of this foam model, we find that cosmic-ray energies cannot exceed the fundamental Planck scale, and so set a lower bound of 10⁸ TeV for the scale of gravity. We also find that suppression of p→pπ⁰ and γ→e⁻e⁺ “decays” favors values aa_crit. Finally, we comment on the apparent non-conservation of particle energy–momentum, and speculate on its re-emergence as dark energy in the foamy vacuum.