Platinum group element abundances in a peat layer associated with the Tunguska event, further evidence for a cosmic origin

We have measured excesses of Pd, Rh, Ru, REE, Co, Sr, and Y in a peat column from the Northern peat bog of the 1908 Tunguska explosion site. Earlier, in this peat column the presence of an Ir anomaly at the event layers (30- depth) has been found (Planet Space Sci. 48 (1998) 179). In these layers, Pd, Rh, Ru, Co, Sr, and Y show pronounced anomalies of a factor 4-7 higher than the background value. In the event layers there are also good correlations between the siderophile platinum group elements (Pd, Rh, Ru) and Co, indicators of cosmic material, which imply they might have the same source, i.e. the Tunguska explosive body. The patterns of CI-chondrite-normalized REE in the event layers are much flatter than those in normal peat layers and different from those in the nearby traps. Furthermore, in these layers the patterns of CI-chondrite-normalized PGEs and the element ratios (e.g. C/Pd, C/Rh, and between some siderophile elements) give evidence that the Tunguska explosive body was more likely a comet, although we cannot exclude the possibility that the impactor could be a carbonaceous asteroid. We have estimated the total mass of a solid component of the explosive body up to 10³-10⁶ tons.     

Properties of energetic water-group ions in the extended pick-up region surrounding comet Giacobini-Zinner

The properties of energetic (65–95 keV) cometary water-group ions in the extended solar wind pick-up region surrounding comet Giacobini-Zinner are examined using data from the EPAS instrument on the ICE spacecraft. In the outer part of this region, extending from cometocentric distances of several hundred thousand to a few million kilometres (the limit of pick-up ion detectability), it is found that large modulations of the ion flux occur (with J_MAX/J_MIN 10²-10³) which are related to the direction of the magnetic field. It is also found that the ions stream in a direction which is intermediate between the directions of the solar wind flow and the E × B drift, and that ions are present at energies somewhat above the local pick-up energy. These properties indicate that the waves which are excited by the unstable “ring-beam” pick-up ion velocity distributions do result in significant scattering of the ions in this region, both in pitch angle and in energy, but that they have insufficient amplitude to scatter the ions into near isotropy in the solar wind frame. Closer to the comet (but still upstream from the bow shock), the ion flux modulations are considerably reduced in amplitude and the ions respond less to the E × B drift, indicating that the ions are scattered nearer to isotropy in this region. Inbound, this transition takes place relatively abruptly at a distance of 4 × 10⁵ km in association with an increase in the solar wind speed, after which the ion flux increases, and ceases to be modulated by the field direction, while the streaming direction is continuously antisolar and unmodulated by the direction of the E × B drift. Outbound, weak vestiges of the ring-beam ion anisotropy are present in the region immediately upstream from the bow shock (at −1 × 10⁵ km), but these become more marked at distances in excess of t4 × 10⁵ km, increasing gradually with increasing distance from the comet. It is shown that the evolution of the ion properties is qualitatively consistent with expectations based on quasi-linear diffusion of the ions by the magnetosonic waves observed during the encounter.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

The cosmic ray primary composition in the “knee” region through the EAS electromagnetic and muon measurements at EAS-TOP

The evolution of the cosmic ray primary composition in the energy range 10⁶–10⁷ GeV (i.e. the “knee” region) is studied by means of the e.m. and muon data of the Extensive Air Shower EAS-TOP array (Campo Imperatore, National Gran Sasso Laboratories). The measurement is performed through: (a) the correlated muon number (N_μ) and shower size (N_e) spectra, and (b) the evolution of the average muon numbers and their distributions as a function of the shower size. From analysis (a) the dominance of helium primaries at the knee, and therefore the possibility that the knee itself is due to a break in their energy spectrum (at E_k^He=(3.5±0.3)×10⁶ GeV) are deduced. Concerning analysis (b), the measurement accuracies allow the classification in terms of three mass groups: light (p,He), intermediate (CNO), and heavy (Fe). At primary energies E₀≈10⁶ GeV the results are consistent with the extrapolations of the data from direct experiments. In the knee region the obtained evolution of the energy spectra leads to: (i) an average steep spectrum of the light mass group (γ_p,He>3.1), (ii) a spectrum of the intermediate mass group harder than the one of the light component (γ_CNO2.75, possibly bending at E_k^CNO≈(6–7)×10⁶ GeV), (iii) a constant slope for the spectrum of the heavy primaries (γ_Fe2.3–2.7) consistent with the direct measurements. In the investigated energy range, the average primary mass increases from lnA=1.6–1.9 at E₀1.5×10⁶ GeV to lnA=2.8–3.1 at E₀1.5×10⁷ GeV. The result supports the standard acceleration and propagation models of galactic cosmic rays that predict rigidity dependent cut-offs for the primary spectra of the different nuclei. The uncertainties connected to the hadronic interaction model (QGSJET in CORSIKA) used for the interpretation are discussed.     

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.     

On the joint robust estimation of parameters and variance

A method of joint estimation of the parameters and the variance σ² in the linear model is presented in this paper. The M-estimators with finite rejection points are adopted for estimating the parameters. In order to estimate the variance σ², the Grubbs statistics and the Kurtosis test statistics are adopted to test the residual sequence {;r_i};. And the sample variance of {;r_i}; after discarding the outliers is taken as the estimation of σ². This method of estimating σ² is less computational demanding and more accurate, compared with the well-known method which takes 1.483 med_i |r_i| as the estimation of σ. The breakdown point of the above estimation of σ² is more than 20%, while the brekdown point is less than 10%, if only the Grubbs statistics is used.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

Finding of probable Tunguska Cosmic Body material:: isotopic anomalies of carbon and hydrogen in peat

Method of a search for traces of Tunguska Cosmic Body (TCB) material using layer-by-layer analysis of the isotopic composition of light elements in peat has been offered. Four peat columns sampled at the explosion epicentre indicated significant carbon and hydrogen isotopic effects in its near catastrophic layers. The shifts, opposite in direction, for carbon (Δ¹³C reaches +4.3‰) and hydrogen (ΔD reaches −22‰) cannot be attributed to any known terrestrial reasons (fall-out of terrestrial dust and fire soot; emission from the Earth of oil–gas streams; climate changes, humification of peat, and so on). Moreover, the isotopic effects are clearly associated with the area and with the time of the 1908 event. They are absent in the uppermost and the lowest peat layers and also in the control peat columns sampled at the remote places. Since calculated δ¹³C value for an admixture of carbon (+51–64‰) is very high, these effects may not be explained by contamination of peat with material similar to ordinary chondrites or achondrites, too. Such heavy carbon occurs in the most primitive CI and CM types of carbonaceous chondrites. However, C/Ir ratio in a cosmic admixture is 10,000 times as many as in CI chondrites that points to cometary nature of the TCB. The isotopic effects are in agreement with the increase of the Ir content observed in peat, but, at the same time, small content of Ir points to the low content of dust in the Tunguska comet that sharply differs it from Halleys comet.     

Density profiles of dark matter halos with anisotropic velocity tensors

We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ³ ∝ r⁻, and (ii) the velocity anisotropy parameter is given by the relation β_a(r)=β₁+2β₂ (r/r_*)/[1+(r/r_*)²] where β₁, β₂ are parameters and r_* equals twice the virial radius, r_vir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10^−2.5r_vir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β₁ and for increasing concentration of the system. On the other hand, slopes at r=r_vir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10⁻³r_vir and r_vir and connects kinematic and structural characteristics of spherical systems is described.     

The probability of detecting first generation stars of the Galaxy

First generation stars are the oldest stars that were formed in post-big bang, primitive gas, containing no elements heavier than carbon, with ages greater than 14 Gyr and having undergone no evolution so far. Observations over a long time have confirmed that, up to now, no stars with zero metallicity ([Fe/H]) or with [Fe/H]≤ −6 have been found in the Galaxy. To explain this absence, we shall make a theoretical calculation of the probability of detecting first generation stars using Tsuiimoto et al.'s model of chemical evolution of the galactic halo and assuming an initial mass function of the Miller-Scalo form. We use all the observational data on the halo stars to constrain the parameters. Our result is that, if the mass of the cloud that formed the stars is 10⁶–10⁷ M, then the probability of detecting first generation stars is 6.14×10⁻⁴–6.14×10⁻⁵.     

Generalized limits to the number of light particle degrees of freedom from big bang nucleosynthesis

We compute the big bang nucleosynthesis limit on the number of light neutrino degrees of freedom in a model-independent likelihood analysis based on the abundances of ⁴He and ⁷Li. We use the two-dimensional likelihood functions to simultaneously constrain the baryon-to-photon ratio and the number of light neutrinos for a range of ⁴He abundances Y_p = 0.225–0.250, as well as a range in primordial ⁷Li abundances from (1.6 to 4.1) ×10⁻¹⁰. For (⁷Li/H)_p = 1.6 × 10⁻¹⁰, as can be inferred from the ⁷Li data from Population II halo stars, the upper limit to N_ν based on the current best estimate of the primordial ⁴He abundance of Y_p = 0.238 is N_ν < 4.3 and varies from N_ν < 3.3 (at 95% C.L.) when Y_p = 0.225 to N_ν < 5.3 when Y_p = 0.250. If ⁷Li is depleted in these stars the upper limit to N_ν is relaxed. Taking (⁷Li/H)_p = 4.1 × 10⁻¹⁰, the limit varies from N_ν < 3.9 when Y_p = 0.225 to N_ν 6 when Y_p = 0.250. We also consider the consequences on the upper limit to N_ν if recent observations of deuterium in high-redshift quasar absorption-line systems are confirmed.     

Observation of a type of sporadic solar mass ejection and its physical interpretation with numerical simulation

Seven mediated and small ejective events on the sun observed at Ganyu Observing Station of Purple Mountain Observatory in 2000 are investigated. It is found that they were not accompanied by brightening. Their lengths were in the range 1–2.5×10⁴ km, their widths, 3–5×10³ km, and their lifetimes, 3–7 minutes. They were produced at places of weak magnetic fields and far away from large sunspots. These ejections are interpreted by numerical simulation with 1-D hydrodynamic equations of flow along magnetic arcs. As demonstrated by the results, they are different from the spicules and surges simulated by Suematsu et al. and Shibata et al. They are not matter with photospheric or chromospheric densities pushed by shock waves or rebound shock waves toward the solar corona, rather, they are ejections formed by continuous matter flows after magnetic reconnection. After evolving for about 5 minutes, they can attain a stationary hydrodynamic state.     

Nightside ionosphere of Venus

Venera 9, 10 measurements of the nightside ionospheric profile and the night airglow were used for investigating ionosphere formation processes. The upper ionospheric layer may be formed by HeI 584 Å radiation; the lower layer by meteorite ionization. Upper limits on the electron energy flux, <4 × 10⁸eV cm⁻² s⁻¹, the helium ion flux <10⁷ cm⁻² s⁻¹, the nitric oxide mixing ratio, <1.5 × 10⁻⁴ and the atomic sulphur mixing ratio, <10⁻⁶, are deduced for ionospheric altitudes.     

Comparison of high-energy galactic and atmospheric tau neutrino flux

We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for 10³E/GeV10¹¹. The intrinsic and oscillated tau neutrino fluxes from both sources are calculated. The intrinsic galactic ν_τ flux (E10³ GeV) is calculated by considering the interactions of high-energy cosmic-rays with the matter present in our galaxy, whereas the oscillated galactic ν_τ flux is coming from the oscillation of the galactic ν_μ flux. For the intrinsic atmospheric ν_τ flux, we extend the validity of a previous calculation from E10⁶ GeV up to E10¹¹ GeV. The oscillated atmospheric ν_τ flux is, on the other hand, rather suppressed. We find that, for 10³E/GeV5×10⁷, the oscillated ν_τ flux along the galactic plane dominates over the maximal intrinsic atmospheric ν_τ flux, i.e., the flux along the horizontal direction. We also briefly mention the presently envisaged prospects for observing these high-energy tau neutrinos.     

A prediction for three neutrino masses and mixings and similarity between quark and lepton mixings

If neutrinos have mass, we give reasons for a possible pattern of three (squaed) mass eigenvalues: m₁² (2.8−5.8) (eV)², m₂² 0.01 (eV)², m₃² (1.5−1) × 10⁻⁴ (eV)². The flavor states ν_μ and ν_e are mixtures of the eigenstates with m₂ and m₃ with a significant mixing, corresponding to an effective mixing angle of about 0.45. The ν_τ is nearly the state with m₁; the other two effective mixing angles are about an order of magnitude smaller than 0.45. There is a marked similarity to mixing in the quark sector.     

Galactic cosmic-ray anisotropy and its heliospheric modulation, inferred from the sidereal semidiurnal variations observed in the rigidity range 300–600 GV with multidirectional muon telescope at Sakashita underground station

The existence of sidereal semidiurnal variation of cosmic-ray intensity in a rigidity region 10²-10³ GV has been reported by many researchers, but there is no consensus of opinion on its origin. In this paper, using the observed semidiurnal variations in a rigidity range (300–600 GV) with 10 directional muon telescopes at Sakashita underground station (geog. lat. = 36°, long. = 138°E, DEPTH = 80 m.w.e.), the authors determine the magnitudes (η₁, η₂) and directions (a₁, a₂) of the first- and second-order anisotropies in the following galactic cosmic-ray intensity distribution (j)

jdp = j₀{1 + η₁P₁(cos χ₁) + η₂P₂(cos χ₂)}dp

, where P_nis the nth order spherical function and χ_n is the pitch angle of cosmic rays with respect to a_n. For the determination, the influence of cosmic-ray's heliomagnetospheric modulation, geomagnetic deflection and nuclear interaction with the terrestrial material and also of the geometric configuration of the telescopes are taken into account. Usually, the semidiurnal variation is produced by the second-order anisotropy. The present observation, however, requires also the first-order anisotropy which usually produces only the diurnal variation, but can produce also the semidiurnal variation as a result of the heliospheric modulation. The first- and second-order anisotropies are characterized with η₁) > 0 and η₂ < 0 have almost the same direction (a₁ a₂) specified by the right ascension ( 0.75 h) and declination (δ 50°S) and, therefore, they can be expressed, as a whole, by an axis-symmetric anisotropy of loss-cone type (i.e. deficit intensities in a cone). It is noteworthy that this anisotropy approximately coincides with that inferred from the air shower observation at Mt Norikura in the rigidity region 10⁴ GV.     

The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Isotopic-geochemical study of nitrogen and carbon in peat from the Tunguska Cosmic Body explosion site

Isotopic-geochemical investigations were carried out on peat samples from the 1908 Tunguska Cosmic Body (TCB) explosion area. We analyzed two peat columns from the Northern peat bog, sampled in 1998, and from the Raketka peat bog, sampled during the 1999 Italian expedition, both located near the epicenter of the TCB explosion area. At the depth of the “catastrophic” layer, formed in 1908, and deeper, one can observe shifts in the isotopic composition of nitrogen (up to Δ¹⁵N = +7.2‰) and carbon (up to Δ¹³C = +2‰) and also an increase in the nitrogen concentration compared to those in the normal, upper layers, unaffected by the Tunguska event. One possible explanation for these effects could be the presence of nitrogen and carbon from TCB material and from acid rains, following the TCB explosion, in the “catastrophic” and “precatastrophic” layers of peat. We found that the highest quantity of isotopically heavy nitrogen fell near the explosion epicenter and along the TCB trajectory. It is calculated that 200,000 tons of nitrogen fell over the area of devastated forest, i.e., only about 30% of the value calculated by Rasmussen et al. (1984). This discrepancy is probably caused by part of the nitrogen having dispersed in the Earth’s atmosphere. The isotopic effects observed in the peat agree with the results of previous investigations [Kolesnikov et al 1998a], [Kolesnikov et al 1998b], [Kolesnikov et al 1999] and [Rasmussen et al 1999] and also with the increased content of iridium and other platinoids found in the corresponding peat layers of other columns [Hou et al 1998] and [Hou et al 2000]. These data favor the hypothesis of a cosmochemical origin of the isotopic effects.     

The density and magnetic field of the dense cores of the molecular cloud L 1630

The density and magnetic field strength of the dense cores in the Orion B molecular cloud are derived from the observed radius and FWHM line width based on the model of a uniformly magnetic sphere. We obtain the average magnetic field strength of 110μG and the average density of 8 × 10⁴/cm³ for the 39 cores, which agree closely with the observations. The method for deriving the density and magnetic field strength is applicable to the cores with R>0.2pc.     

Traces of cometary material in the area of the Tunguska impact (1908)

A short overview of the studies of the authors and their colleagues performed over many years, which resulted in the discovery of traces of cometary matter in the peat at the epicenter of the Tunguska catastrophe in 1908, is given here. In the epicenter of the Tunguska cosmic body (TCB) explosion, the shifts in the isotopic composition of hydrogen and carbon relative to their values for the upper and lower layers of the same column were found in the catastrophic layers of peat grown up in 1908. These shifts cannot be attributed to any known terrestrial processes: the conservation of mineral and organic dust in peat, peat humification, the emission of hydrocarbon gases from the Earth, climate changes, and other physical and chemical processes. In the catastrophic layers of the control peat columns, the isotopic shifts are absent. The isotopic data agree well with the increased concentration of iridium and other platinum-group elements in the same peat layers, which is a reliable indicator of the presence of cosmic material in terrestrial objects. The cosmogenic character of the isotopic effects is confirmed by the presence of “dead” carbon (not containing radioactive ¹⁴C) in the catastrophic layers. To provide the shifts observed in the isotopic composition of carbon, cosmic carbon preserved in peat should be isotopically superheavy—from +50‰ to +60‰ according to calculations. Such isotopically heavy carbon is absent both on the Earth and in ordinary meteorites. It occurs only in individual mineral phases of CI carbonaceous chondrites, close to cometary dust in chemical composition, ratios of the content of iridium and other platinoids and rear-earth elements also points to the cometary nature of the TCB. In the near-catastrophic peat layers, the anomalous increase of the concentration of many volatiles was detected, which also suggests that the TCB was a cometary core. The studies of the content and the isotopic composition of nitrogen in the peat revealed traces of heavy acid rains induced by the flyby and explosion of the TCB.