Upper limit on the cosmic-ray photon flux above 10 eV using the surface detector of the Pierre Auger Observatory

A method is developed to search for air showers initiated by photons using data recorded by the surface detector of the Auger Observatory. The approach is based on observables sensitive to the longitudinal shower development, the signal risetime and the curvature of the shower front. Applying this method to the data, upper limits on the flux of photons of 3.8×10^-3, above , are derived, with corresponding limits on the fraction of photons being 2.0%, 5.1%, and 31% (all limits at 95% c.l.). These photon limits disfavor certain exotic models of sources of cosmic rays. The results also show that the approach adopted by the Auger Observatory to calibrate the shower energy is not strongly biased by a contamination from photons.     

New estimates for Io eruption temperatures: Implications for the interior

The initial interpretation of Galileo data from Jupiter's moon, Io, suggested eruption temperatures . Tidal heating models have difficulties explaining Io's prodigious heat flow if the mantle is , although we suggest that temperatures up to may be possible. In general, Io eruption temperatures have been overestimated because the incorrect thermal model has been applied. Much of the thermal emission from high-temperature hot spots comes from lava fountains but lava flow models were utilized. We apply a new lava fountain model to the highest reported eruption temperature, the SSI observation of the 1997 eruption at Pillan. This resets the lower temperature limit for the eruption from 1600 to . Additionally, viscous heating of the magma may have increased eruption temperature by as a result of the strong compressive stresses in the ionian lithosphere. While further work is needed, it appears that the discrepancy between observations and interior models is largely resolved.     

Interpretation of the underground muon charge ratio

The MINOS experiment has observed a rise in the underground muon charge ratio r_μ=μ⁺/μ^-. This ratio can be related to the atmospheric production ratios of π⁺/π^- and K⁺/K^-. Our analysis indicates that the relevant variable for studying the charge ratio is , rather than . We compare a simple energy dependent parameterization of the rise in the charge ratio with more detailed previously published Monte Carlo simulations and an analytical calculation. We also estimate the size of two previously neglected effects in this context: the charge sign dependency of the dE/dx in rock, and the energy dependence of heavy primaries on the derived K⁺/K^- ratio.     

Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV , respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top–down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.     

When did cosmic acceleration start? How fast was the transition?

Cosmic acceleration is investigated through a kink-like expression for the deceleration parameter (q). The new parametrization depends on the initial (q_i) and final (q_f) values of q, on the redshift of the transition from deceleration to acceleration (z_t) and the width of such transition (τ). We show that although supernovae (SN) observations (Gold182 and SNLS data samples) indicate, at high confidence, that a transition occurred in the past (z_t > 0) they do not, by themselves, impose strong constraints on the maximum value of z_t. However, when we combine SN with the measurements of the ratio between the comoving distance to the last scattering surface and the SDSS + 2dfGRS BAO distance scale (S_k/D_v) we obtain, at 95.4% confidence level, and for (S_k/D_v+Gold182), and and for (S_k/D_v+SNLS), assuming q_i = 0.5 and q_f = −1. We also analyze the general case, q_f  (−∞, 0) finding the constraints that the combined tests (S_k/D_v+SNLS) impose on the present value of the deceleration parameter (q₀).     

Search for VHE gamma rays from SS433/W50 with the CANGAROO-II telescope

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate VHE gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect VHE gamma rays using the 10 m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 h for ON source, and 80.8 h for OFF source data. No significant excess of VHE gamma rays has been found at three regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are G for the center of the brightest X-ray emission region and G for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.     

Thinned simulations of extremely energetic showers in dense media for radio applications

This paper presents our results for full 3-D simulations of very-high to ultra-high energy electromagnetic cascades – and the associated coherent Cherenkov radiation – as might be produced by high-energy neutrino interactions in dense media. Using “thinning” techniques, we develop an algorithm based on the existing “ZHS” code, and demonstrate that the new “ZHS-thinned” code can produce fast and accurate results for showers up to . Using ZHS-thinned, we develop new parameterisations for the radiation from showers in ice, salt, and the lunar regolith, with a separate treatment of the megaregolith (deep regolith). Our parameterisations include for the first time a method to simulate fluctuations in shower length induced by the LPM effect. Our results, which avoid the pit-falls of scaling simulations from lower energies, allow improved calculations of the detection probability for experiments searching for high-energy neutrinos using the radio technique.     

Correlations between environment and other properties of galaxies from the Sloan Digital Sky Survey at fixed color

From the volume-limited Main galaxy sample of the Sloan Digital Sky Survey Data Release 6 (SDSS DR6), we construct three samples with g–r color bins , labeled S1–S3, to investigate how other properties of galaxies depend on environment at fixed color. For each sample, we measure the local three-dimensional galaxy density in a comoving sphere with radius equal to the distance to the 5th nearest galaxy for each galaxy, select about 5% galaxies and construct the two subsamples at both extremes of density. Our study suggests that the environmental dependence of luminosity is mainly due to the environmental dependence of galaxy color and the correlation between color and luminosity. In addition, we preferentially conclude that concentration index and morphologies are not strongly correlated with local density at fixed color, and that galaxy color is a galaxy property very predictive of the local environment. Because SDSS spectroscopy is incomplete for bright galaxies at very low redshifts, we also use a volume-limited Main galaxy sample with a lower redshift limit z = 0.05, which contains 94,954 galaxies at 0.05 < z < 0.089 with −22.40 < M_r < −20.16, and reach the same conclusions.Due to the bimodality of the u–r color distribution, we classify galaxies as ‘red’ and ‘blue’, respectively, and further subdivide the samples into star-forming galaxies and passive ones using Hα equivalent width, W₀(Hα). Results show that color and star formation activity of galaxies are galaxy properties very predictive of the local environment.     

Development of an -based spectral analysis system for the coded-aperture hard X-ray balloon payload EXITE2

We present the spectral analysis system for the second-generation energetic X-ray imaging telescope experiment (EXITE2) balloon payload. EXITE2 is an imaging hard X-ray telescope using a coded-aperture mask and a NaI/CsI phoswich detector operating in the energy range 20–600 keV. The instrument was flown on a high-altitude scientific balloon from Ft. Sumner, NM on 7–8 May, 1997. We describe the details of the EXITE2 spectral analysis system, with emphasis on those aspects peculiar to coded-aperture instruments. In particular, we have made our analysis compatible with the standard X-ray spectral fitting package by generating a response matrix in the appropriate format including all the effects of a coded-aperture system. The use of , which may be a first for coded-aperture data, permits great flexibility in the fitting of spectral models. The additional effects of our phoswich system, or any other detector-specific considerations, may be easily included as well. We test our spectral analysis using observations of the Crab Nebula, and find that the EXITE2 Crab spectrum is consistent with those recorded by previous instruments operating in this energy range.     

Climate change of the last millennium inferred from borehole temperatures: regional patterns of climatic changes in the Czech Republic — Part III

Accurate temperature–depth profiles may help to assess the temperature variations associated with the climate changes in the past. Ninety-eight ground surface temperature histories inverted from the temperature–depth borehole logs drilled on the territory of the Czech Republic [Bodri, L., ermák, V., 1995. Climate changes of the last millennium inferred from borehole temperatures: results from the Czech Republic — Part I. Global Planet. Change 11, pp. 111–125; Bodri, L., ermák, V., 1997. Climate changes of the last two millennia inferred from borehole temperatures: results from the Czech Republic — Part II. Global Planet. Change 14, pp. 163–173.] are used to reconstruct the regional patterns of the respective climate change. The climate was mapped for the following periods: 1100–1300 A.D. (Little Climatic Optimum), 1400–1500 A.D., 1600–1700 A.D. (main phase of the Little Ice Age), and for the most recent climate trend after year 1960. Comparison of the obtained maps with the meteorological observations and proxy climatic reconstructions confirmed good applicability of the “geothermal” paleoclimatic reconstructions for the regional studies.     

Cosmic ray composition estimation below the knee of the spectrum from the time structure of erenkov light in EAS

Monte Carlo simulations show that the pulse profile of erenkov photons measured near the core of an extensive air shower is sensitive to the secondary muon/electron ratio of the cascade. erenkov pulses can easily be measured with a single large area mirror viewed by a photomultiplier tube subtending a small field of view (1°). Even for such a simple experiment, exposed to EAS from a range of core locations and arrival directions, strong statistical differences are shown to exist between the pulse parameter distributions of primary protons and those of heavier primary particles. A range of primary energies can be investigated by varying the zenith angle of observations. In this paper, results from simulations of primaries in the energy range 20 TeV to 400 TeV are presented, although in principle the technique could be extended to include the knee of the spectrum. At the lower end of this energy range results can be compared to direct measurements of the composition, while measurements at the upper end can augment results from existing ground based experiments.     

Performance of the stereoscopic system of the HEGRA imaging air erenkov telescopes: Monte Carlo simulations and observations

Based on the Monte Carlo simulations we have studied the performance of the HEGRA system of imaging air erenkov telescopes (IACTs) in its present configuration of 4 IACTs as well as in its future final configuration of 5 IACTs. Here we present the results on the basic characteristics of the IACT system which are used in the standard data analysis procedure, i.e., the collection areas, the detection rates, the angular resolution, the energy resolution, and the γ/hadron-separation efficiency. By comparing several key Monte Carlo predictions with experimental results it is possible to check the accuracy of the simulations. The Monte Carlo results concerning hadron-nuclear showers are tested with the recorded cosmic ray events and the results concerning photon-induced showers are tested with a large data sample of γ-rays observed from BL Lac object Mkn 501 during its high flaring activity in 1997. Summarizing the simulations and current observations we give the basic recommendations of using the instrument and the major values of its sensitivity.     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

The Titan N/N and C/C isotopic ratios in HCN from Cassini/CIRS

We report the detection of H¹³CN and HC¹⁵N in mid-infrared spectra recorded by the Composite Infrared Spectrometer (CIRS) aboard Cassini, along with the determination of the ¹²C/¹³C and ¹⁴N/¹⁵N isotopic ratios. We analyzed two sets of limb spectra recorded near 13-15° S (Tb flyby) and 83° N (T4 flyby) at 0.5 cm⁻¹ resolution. The spectral range 1210-1310 cm⁻¹ was used to retrieve the temperature profile in the range 145-490 km at 13° S and 165-300 km at 83° N. These two temperature profiles were then incorporated in the atmospheric model to retrieve the abundance profile of H¹²C¹⁴N, H¹³CN and HC¹⁵N from their bands at 713, 706 and 711 cm⁻¹, respectively. The HCN abundance profile was retrieved in the range 90-460 km at 15° S and 165-305 km at 83° N. There is no evidence for vertical variations of the isotopic ratios. Constraining the isotopic abundance profiles to be proportional to the HCN one, we find at 15° S, and at 83° N, two values that are statistically consistent. A combination of these results yields a ¹²C/¹³C value equal to 75±12. This global result, as well as the 15° S one, envelop the value in Titan's methane (82.3±1) [Niemann, H.B., and 17 colleagues, 2005. Nature 438, 779-784] measured at 10° S and is slightly lower than the terrestrial inorganic standard value (89). The ¹⁴N/¹⁵N isotopic ratio is found equal to at 15° S and at 83° N. Combining the two values yields ¹⁴N/¹⁵N = 56 ± 8, which corresponds to an enrichment in ¹⁵N of about 4.9 compared with the terrestrial ratio. These results agree with the values obtained from previous ground-based millimeter observations [Hidayat, T., Marten, A., Bézard, B., Gautier, D., Owen, T., Matthews, H.E., Paubert, G., 1997. Icarus 126, 170-182; Marten, A., Hidayat, T., Biraud, Y., Moreno, R., 2002. Icarus 158, 532-544]. The ¹⁵N/¹⁴N ratio found in HCN is ∼3 times higher than in N₂ [Niemann, H.B., and 17 colleagues, 2005. Nature 438, 779-784], which implies a large fractionation process in the HCN photochemistry.     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

Parametric analysis of modeled ion escape from Mars

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors - that of the EUV flux R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , where ?(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O⁺, we obtain the values, 1.73 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ?(i) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.     

Observations of the inner coma of C/1995O1 (Hale-Bopp)—gas and dust production

Hale-Bopp (C/1995 O1) was the most productive recent comet observed in terms of gas and dust output. Since its discovery in 1995 at a distance of 7.14 AU from the Sun, the comet has been well observed, revealing the dynamics of a rare and large comet. Hale-Bopp showed strong emissions of the principle cometary gases CN, C₃, and C₂, as well as an abundance of dust. The production rates of these gases were found to be 1.45×10²⁸, 1.71×10²⁸, and , respectively, with dust production, in terms of Afρ, , as measured in the green continuum (5260 Å). The observations for this paper are presented in two groups spanning 10 days each, one group centered near 32 days prior to and the other 21 days after perihelion. The averages of dust and gas production rates show a slightly higher value for each prior to perihelion than after perihelion, consistent with a possible peak in production a few weeks prior to perihelion passage.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.