Observation of a thin Es-layer by the eiscat radar

High resolution E-region measurements carried out on 16 November 1983 using the EISCAT incoherent scatter radar are presented. The experiment was monostatic with a vertical radar beam, and it was based on a Barker-coded four-pulse code on one frequency channel and Barker-coded single pulses on three channels. The basic integration time was 15 s and the spatial resolution 450 m. The results reveal a short-lived but intense thin sporadic E-layer at 18:00–18:06 U.T. at an altitude of about 106 km. Both before and during the event, downward ion velocities of the order of 100 m s⁻¹ are observed above this height. A convergent null in the vertical ion speed is occasionally seen at the layer altitude. The layer occurrence is associated with auroral arcs drifting across the radar beam. Simultaneous observations of the STARE radar show an ionospheric electric field of 25–30 mV m⁻¹. The field always has a westward component, which is in accordance with the observed downward plasma flow. Most of the time when the layer is intense, the field points into the NW-sector. Theoretically, this field direction should create convergent vertical plasma motion. Therefore it is suggested that the observed E_s-layer is created by the action of the auroral electric field rather than by the wind shear mechanism.     

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.     

Deactivation of N2A Σumolecules in the aurora

Recent rocket observations of the N₂ V-K (Vegard-Kaplan) system in the aurora have been reinterpreted using an atmospheric model based on mass spectrometer measurements in an aurora of similar intensity at the same time of year. In contrast to the original interpretation, we find that population by cascade from the C³Π_u and B³Π_g states in the A³Σ_u⁺v=0,1 levels, as calculated using recently measured electron excitation cross sections, accurately accounts for the observed relative emission rates (IV-K/12PG_0.0). In addition there is no need to change the production rate of A ³ Σ _u⁺ molecules relative to that of C³Π_uv=0 as a function of altitude in order to fit the profile of the deactivation probability to the atmospheric model. Quenching of A ³ Σ _u⁺ molecules at high altitudes is dominated by atomic oxygen. The rate constants for the v=0 and v=1 levels are 8 × 10⁻¹¹ cm³ sec⁻¹ and 1.7 × 10⁻¹⁰ cm³ sec⁻¹ respectively, as determined using the model atmosphere mentioned above. Recent observations with a helium cooled mass spectrometer suggest that conventional mass spectrometer measurements tend to underestimate the atomic oxygen relative concentration. The rate coefficients may therefore be too large by as much as a factor of 3. Below 130 Km we find that it is possible to account for the deactivation in bright auroras by invoking large nitric oxide concentrations, similar to those recently observed mass spectrometrically and using a rate constant of 8 × 10⁻¹¹ cm³ sec⁻¹ for both the v=1 levels. This rate constant is very nearly the same as that measured in the laboratory (7 × 10⁻¹¹ cm³ sec⁻¹). Molecular oxygen appears not to play a significant role in deactivating the lower A ³ Σ _u⁺ levels.     

Cratering of icy targets by different impactors: Laboratory experiments and implications for cratering in the Solar System

Studies of impacts (impactor velocity about 5 km s⁻¹) on icy targets were performed. The prime goal was to study the response of solid CO₂ targets to impacts and to find the differences between the results of impacts on CO₂ targets with those on H₂O ice targets. The crater dimensions in CO₂ ice were found to scale with impact energy, with little dependence on projectile density (which ranged from nylon to copper, i.e., 1150-8930 kg m⁻³). At equal temperatures, craters in CO₂ ice were the same diameter as those in water ice, but were shallower and smaller in volume. In addition, the shape of the radial profiles of the craters was found to depend strongly on the type of ice and to change with impact energy. The impact speed of the data is comparable to that for impacts on many types of icy bodies in the outer Solar System (e.g., the satellites of the giant planets, the cometary nuclei and the Kuiper Belt objects), but the size and thus energy of the impactors is lower. Scaling with impact energy is demonstrated for the impacts on CO₂ ice. The issue of impact disruption (rather than cratering) is discussed by analogy with that on water ice. Expressions for the critical energy density for the onset of disruption rather than cratering are established for water ice as a function of porosity and silicate content. Although the critical energy density for disruption of CO₂ ice is not established, it is argued that the critical energy to disrupt a CO₂ ice body will be greater than that for a (non-porous) water ice body of the similar mass.     

Impact experiments in low-temperature ice

New results of low-velocity impact experiments in cubic and cylindrical (20 cm) water-ice targets initially at 257 and 81 °K are reported. Impact velocities and impact energies vary between 0.1 and 0.64 km/sec and 10⁹ and 10¹⁰ ergs, respectively. Observed crater diameters range from 7 to 15 cm and are two to three times larger than values found for equal-energy impacts in basaltic targets. Crater dimensions in ice targets increase slightly with increasing target temperatures. Crater volumes of strength-controlled ice craters are about 10 to 100 times larger than those observed for craters in crystalline rocks. Based on similarity analysis, general scaling laws for strength-controlled crater formation are derived and are applied to crater formation on the icy Galilean and Saturnian satellites. This analysis indicates that surface ages, based on impact-crater statistics on an icy crust, will appear greater than those for a silicate crust which experienced the same impact history. The greater ejecta volume for cratering in ice versus cratering in silicate targets leads to accelerated regolith production on an icy planet.     

Impacts onto H2O ice: Scaling laws for melting, vaporization, excavation, and final crater size

Shock-induced melting and vaporization of H₂O ice during planetary impact events are widespread phenomena. Here, we investigate the mass of shock-produced liquid water remaining within impact craters for the wide range of impact conditions and target properties encountered in the Solar System. Using the CTH shock physics code and the new 5-phase model equation of state for H₂O, we calculate the shock pressure field generated by an impact and fit scaling laws for melting and vaporization as a function of projectile mass, impact velocity, impact angle, initial temperature, and porosity. Melt production nearly scales with impact energy, and natural variations in impact parameters result in only a factor of two change in the predicted mass of melt. A fit to the π-scaling law for the transient cavity and transient-to-final crater diameter scaling are determined from recent simulations of the entire cratering process in ice. Combining melt production with π-scaling and the modified Maxwell Z-model for excavation, less than half of the melt is ejected during formation of the transient crater. For impact energies less than about 2 × 10²⁰ J and impact velocities less than about 5 km s⁻¹, the remaining melt lines the final crater floor. However, for larger impact energies and higher impact velocities, the phenomenon of discontinuous excavation in H₂O ice concentrates the impact melt into a small plug in the center of the crater floor.     

Laboratory impacts into dry and wet sandstone with and without an overlying water layer: Implications for scaling laws and projectile survivability

Abstract— Scaling laws describing crater dimensions are defined in terms of projectile velocity and mass, densities of the materials involved, strength of the target, and the local gravity. Here, the additional importance of target porosity and saturation, and an overlying water layer, are considered through 15 laboratory impacts of 1 mm diameter stainless steel projectiles at 5 km s^?1 into a) an initially uncharacterized sandstone (porosity ?17%) and b) Coconino Sandstone (porosity ?23%). The higher‐porosity dry sandstone allows a crater to form with a larger diameter but smaller depth than in the lower‐porosity dry sandstone. Furthermore, for both porosities, a greater volume of material is excavated from a wet target than a dry target (by 27–30%). Comparison of our results with Pi‐scaling (dimensionless ratios of key parameters characterizing cratering data over a range of scales) suggests that porosity is important for scaling laws given that the new data lie significantly beneath the current fit for ice and rock targets on a π_v versus π₃ plot (π_v gives cratering efficiency and π₃ the influence of target strength). An overlying water layer results in a reduction of crater dimensions, with larger craters produced in the saturated targets compared to unsaturated targets. A water depth of approximately 12 times the projectile diameter is required before craters are no longer observed in the targets. Previous experimental studies have shown that this ratio varies between 10 and 20 (Gault and Sonett 1982). In our experiments ?25% of the original projectile mass survives the impact.     

Electron excitation and auroral emission parameters

Auroral luminosities of the main emission lines in the aurora have been calculated for excitation by an isotopic primary electron flux with spectra of the form J(E) = AE exp (−E/E₁) + B(E₂)E exp (−E/E₁). The variation of emissions from O and N₂⁺ with height are shown, as are the variations of column integrated intensities and pertinent intensity ratios with the characteristic energy E₂, this leading to a method of estimating the electron spectrum from ground observation.     

Cooling curves and initial models for low-mass white dwarfs (<0.25 M⊙) with helium cores

We present a detailed calculation of the evolution of low-mass (<0.25 M_⊙) helium white dwarfs. These white dwarfs (the optical companions to binary millisecond pulsars) are formed via long-term, low-mass binary evolution. After detachment from the Roche lobe, the hot helium cores have a rather thick hydrogen layer with mass between 0.01 and 0.06 M_⊙. As a result of mixing between the core and outer envelope, the surface hydrogen content ( X _surf) is 0.5–0.35 , depending on the initial value of the heavy element Z and the initial secondary mass. We found that the majority of our computed models experience one or two hydrogen shell flashes. We found that the mass of the helium dwarf in which the hydrogen shell flash occurs depends on the chemical composition. The minimum helium white dwarf mass in which a hydrogen flash takes place is 0.213 M_⊙ ( Z =0.003), 0.198 M_⊙ ( Z =0.01), 0.192 M_⊙ ( Z =0.02) or 0.183 M_⊙ ( Z =0.03). The duration of the flashes (independent of chemical composition) is between a few ×10⁶ and a few ×10⁷ yr. In several flashes the white dwarf radius will increase so much that it forces the model to fill its Roche lobe again. Our calculations show that the cooling history of the helium white dwarf depends dramatically on the thickness of the hydrogen layer. We show that the transition from a cooling white dwarf with a temporarily stable hydrogen-burning shell to a cooling white dwarf in which almost all residual hydrogen is lost in a few thermal flashes (via Roche lobe overflow) occurs between 0.183 and 0.213 M_⊙ (depending on the heavy element value).     

Effective collision strengths for forbidden transitions among the 3s23p3 fine-structure levels of Cl IIIIII

Effective collision strengths for the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 3s²3p³ configuration of Cl  iii are presented. The calculation employs the multichannel R-matrix method to compute the electron-impact excitation collision strengths in a close-coupling expansion, which incorporates the lowest 23 LS target eigenstates of Cl  iii . These states are formed from the 3s²3p³, 3s3p⁴, 3s²3p²3d and 3s²3p²4s configurations. The Maxwellian-averaged effective collision strengths are presented graphically for all 10 fine-structure transitions over a wide range of electron temperatures appropriate for astrophysical applications [log  T (K)=3.3−log  T (K)=5.9]. Comparisons are made with the earlier seven-state close-coupling calculation of Butler & Zeippen, and in general excellent agreement is found in the low-temperature region where a comparison is possible [log  T (K)=3.3−log  T (K)=4.7]. However, discrepancies of up to 30 per cent are found to occur for the forbidden transitions which involve the ⁴S^o ground state level, particularly for the lowest temperatures considered. At the higher temperatures, the present data are the only reliable results currently available.     

Auroral N2(c′4Σu → aΠg emission

Inspection of recent spectra presented by Sivjee (1983) show evidence of the 0–4 and 0–5 bands of the N₂(c′₄¹Σ_u⁺ → a¹Π_g) Gaydon-Herman system. In conjunction with earlier spectra, it is now possible that this band system is a significant auroral component, with an intensity approx. 7% that of the N2 2P system. The absence in aurorae of the potentially far stronger N₂(c′₄¹Σ_u⁺ → X¹Π_g) system is discussed. It is that the O₂(A³Σ_u⁺ → X³Σ_g⁻) band system is indiscernible in Sivjee's auroral spectra, under conditio the foreground nightglow is expected to be clearly visible. On the other hand, at least one relatively strong O₂(A′³Δ_u → a¹Δ_g) band appears to be present in these spectra.     

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of Ne iv

Effective collision strengths for electron-impact excitation of the N-like ion Ne  iv are calculated in the close-coupling approximation using the multichannel R-matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 2s²2p³ ground-state configuration. The expansion of the total wavefunction incorporates the lowest 11 LS eigenstates of Ne  iv , consisting of eight n  = 2 terms with configurations 2s²2p³, 2s2p⁴ and 2p⁵, together with three n  = 3 states of configuration 2s²2p²3s. We present in graphical form the effective collision strengths obtained by thermally averaging the collision strengths over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log T (K) = 3.6 to log T (K) = 6.1 (the range appropriate for astrophysical applications). Comparisons are made with the earlier, less sophisticated close-coupling calculation of Giles, and excellent agreement is found in the limited temperature region where a comparison is possible [log T (K) = 3.7 to log T (K) = 4.3]. At higher temperatures the present data are the only reliable results currently available.     

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of S x

Effective collision strengths for electron-impact excitation of the N-like ion S  x are calculated in the close-coupling approximation using the multichannel R -matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 2s²2p³ ground configuration. The total (e⁻+ion) wavefunction is expanded in terms of the 11 lowest LS eigenstates of S  x , and each eigenstate is represented by extensive configuration-interaction wavefunctions. The collision strengths obtained are thermally averaged over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log  T (K)=4.6–6.7 (the range appropriate for astrophysical applications). The present effective collision strengths are the only results currently available for these fine-structure transition rates.     

The absorption of radio waves in the ionosphere with respect to geomagnetic current control and other proposed models—a critical survey

It is shown that the arguments advanced by Shaw⁽¹⁾ to demonstrate that the absorption of radio waves in the ionosphere is controlled by the currents causing geomagnetic variations are unsound. Further the method used by Bandyopadhyay⁽²⁾ in deducing the nondeviative absorption leads to too high a proportion of this absorption in the total. The two D-regions model proposed by Rumi⁽³⁾ is also unsatisfactory in several respects. In all three papers, error arises because of the neglect of the deviative absorption in E-region. The reason for this neglect may be because of the resemblance between the frequency variation of E-region deviative absorption and that of the non-deviative absorption, except in the immediate vicinity of ƒ₀E.     

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.     

Gamma ray spectrum of the crab nebula in the multi TeV region

The THEMISTOCLE array of 18 Cherenkov detectors which has a 3 TeV gamma energy threshold, has detected a signal from the Crab nebula at a 5.8 standard deviation level. Information on the energy spectrum is obtained in the range 3–15 TeV. The integrated flux can be fitted with the form, Φ (> E) = (3.7 ± 0.5) × 10^-12 (E/5)-^{−1.5 ± 0.20} cm⁻² s⁻¹ (E in TeV) compatible with the extrapolation of results at lower energies. The Crab signal is used to measure the angular resolution of the multi-telescope technique. The value obtained is 2.3 mr (0.15°) in agreement with the results of simulations, and confirms the interest of this new method for multi-TeV gamma-ray detection.     

Ultra-high-resolution observations of circumstellar K i and C2 around the post-AGB star HD 56126

We have used the Ultra-High-Resolution Facility (UHRF) at the AAT, operating at a resolution of 0.35 km s⁻¹ (FWHM), to observe K  i and C₂ absorption lines arising in the circumstellar environment of the post-AGB star HD 56126. We find three narrow circumstellar absorption components in K  i , two of which are also present in C₂. We attribute this velocity structure to discrete shells resulting from multiple mass-loss events from the star. The very high spectral resolution has enabled us to resolve the intrinsic linewidths of these narrow lines for the first time, and we obtain velocity dispersions ( b -values) of 0.2–0.3 km s⁻¹ for the K  i components, and 0.54±0.03 km s⁻¹ for the strongest (and best defined) C₂ component. These correspond to rigorous kinetic temperature upper limits of 211 K for K  i and 420 K for C₂, although the b -value ratio implies that these two species do not co-exist spatially. The observed degree of rotational excitation of C₂ implies low kinetic temperatures ( T _k≈10 K) and high densities ( n ≈10⁶ to 10⁷ cm⁻³) within the shell responsible for the main C₂ component. Given this low temperature, the line profiles then imply either mildly supersonic turbulence or an unresolved velocity gradient through the shell.     

Flux limits on ultra high energy neutrinos with AMANDA-B10

Data taken during 1997 with the AMANDA-B10 detector are searched for a diffuse flux of neutrinos of all flavors with energies above 10¹⁶ eV. At these energies the Earth is opaque to neutrinos, and thus neutrino induced events are concentrated at the horizon. The background are large muon bundles from down-going atmospheric air shower events. No excess events above the background expectation are observed and a neutrino flux following E⁻², with an equal mix of all flavors, is limited to E²Φ(10¹⁵ eV < E < 3 × 10¹⁸ eV)  0.99 × 10⁻⁶ GeV cm⁻² s⁻¹ sr⁻¹ at 90% confidence level. This is the most restrictive experimental bound placed by any neutrino detector at these energies. Bounds to specific extraterrestrial neutrino flux predictions are also presented.     

N2 emission in the twilight

About a year's observations of the N₂⁺ band (3914 Å) at Kitt Peak (latitude 32°) are reported. Morning intensities are the same throughout the year, but there is a strong winter maximum in the evening. It is suggested that the additional ionization is produced by photoelectrons from the magnetic conjugate point. Heights are estimated by the zenith-horizon method, which gives 235 km for the constant component and 350 km during the evening enhancement. The intensity variation through twilight is therefore entirely due to changes of the N₂⁺ concentration; each ion scatters light at a constant rate. The rotational distribution resembles that for a temperature of 1600°K, much higher than the temperature of the atmosphere. It is suggested that part of the ions may be produced by charge transfer from metastable O⁺(²D). N₂⁺ concentrations resulting from photoionization are calculated; they give a fair account of the observed horizon intensities, but not the zenith. Non-local electrons from higher in the atmosphere are suggested as a possible extra source; alternatively, the zenith measurements may be perturbed by scattered horizon light. The band intensity in the nightglow cannot be measured; the upper limit is 1 R.     

The Morphology of Small Craters and Knobs on the Surface of Jupiter's Satellite Callisto

Using the images of Callisto's surface acquired at 15-km resolution by the Galileo spacecraft during its C21 orbit, we studied the morphology of craters with diameters of less than 1–2 km and knobs. By analogy with other regions of Callisto that have been studied, these craters and knobs are thought to be formed by the sublimation degradation of the rims of larger craters that are also present in the region under study. The small craters closely resemble similar-sized lunar craters and, by analogy with the latter, are also divided into morphological classes. The depths of 42 craters of different morphological classes are estimated using shadow lengths visible in the craters. The fractions of the craters of different classes in the subpopulation are determined as a function of the crater diameter. Evidence has been obtained that larger craters degrade at a slower rate than smaller ones. The mean thickness of the mantle of dark material (40 m) is estimated from the sizes of the craters ejecting the blocks of the basement ice material. The shape of the knob shadows shows that the knobs are heights of mostly conical form with slopes whose steepness is close to the angle of repose. Analysis has shown that the observed landforms and material units of the region under investigation have been formed during two successive stages of the geologic history of Callisto. Large craters, knobs, and the mantle of dark material were formed mostly at the end of the period of heavy meteorite bombardment. The leading processes of this period are impact cratering, the sublimation of Callisto's crustal ice with the accumulation of residual non-icy material, and downslope mass movement. The next stage, which continues until the present time, involved the formation of the subpopulation of small (<1–2 km) craters. This formation was accompanied by the impact reworking of the upper portion of the dark mantle. The key processes occurring at this stage are impact cratering and downslope mass movement. The mean intensity of resurfacing at this stage is much lower than at the preceding stage.