Impact of atmospheric parameters on the atmospheric Cherenkov technique

Atmospheric density profiles as well as several light absorption and scattering processes depend on geographic position and are generally time-variable. Their impact on the atmospheric Cherenkov technique in general (imaging or non-imaging) is investigated. Different density profiles lead to differences in Cherenkov light density of up to 60%. Seasonal variations at mid-latitude sites are of the order of 15–20%. The quest for improved energy calibration of Cherenkov experiments also shows the need for improved transmission calculations, taking all relevant processes into account and using realistic profiles of absorbers. Simulations including the scattering mechanisms also reveal the relevance of Rayleigh and Mie scattering for atmospheric Cherenkov experiments. Refraction and the differences between treating the atmosphere in plane-parallel or spherical geometry are also investigated.     

Impact of aerosols and adverse atmospheric conditions on the data quality for spectral analysis of the H.E.S.S. telescopes

The Earth’s atmosphere is an integral part of the detector in ground-based imaging atmospheric Cherenkov telescope (IACT) experiments and has to be taken into account in the calibration. Atmospheric and hardware-related deviations from simulated conditions can result in the mis-reconstruction of primary particle energies and therefore of source spectra. During the eight years of observations with the High Energy Stereoscopic System (H.E.S.S.) in Namibia, the overall yield in Cherenkov photons has varied strongly with time due to gradual hardware aging, together with adjustments of the hardware components, and natural, as well as anthropogenic, variations of the atmospheric transparency. Here we present robust data selection criteria that minimize these effects over the full data set of the H.E.S.S. experiment and introduce the Cherenkov transparency coefficient as a new atmospheric monitoring quantity. The influence of atmospheric transparency, as quantified by this coefficient, on energy reconstruction and spectral parameters is examined and its correlation with the aerosol optical depth (AOD) of independent MISR satellite measurements and local measurements of atmospheric clarity is investigated.     

Image and non-image parameters of atmospheric Cherenkov events: a comparative study of their γ-ray/hadron classification potential in ultrahigh energy regime

In this exploratory simulation study, we compare the event-progenitor classification potential of a variety of measurable parameters of atmospheric Cherenkov pulses which are produced by ultrahigh energy γ-ray and hadron progenitors and are likely to be recorded by the TACTIC (TeV atmospheric Cherenkov telescope with imaging camera) array of atmospheric Cherenkov telescopes. The parameters derived from Cherenkov images include Hillas, fractal and wavelet moments, while those obtained from non-image Cherenkov data consist of pulse profile rise time and base width and the relative ultraviolet to visible light content of the Cherenkov event. It is shown by a neural-net approach that these parameters, when used in suitable combinations, can bring about a proper segregation of the two event types, even with modest sized data samples of progenitor particles.     

Detection of atmospheric Cherenkov radiation using solar heliostat mirrors

There is considerable interest world-wide in developing large area atmospheric Cherenkov detectors for ground-based gamma-ray astronomy. This interest stems, in large part, from the fact that the gamma-ray energy region between 20 and 250 GeV is unexplored by any experiment. Atmospheric Cherenkov detectors offer a possible way to explore this region, but large photon collection areas are needed to achieve low energy thresholds. We are developing an experiment using the heliostat mirrors of a solar power plant as the primary collecting element. As part of this development, we built a detector using four heliostat mirrors, a secondary Fresnel lens, and a fast photon detection system. In November 1994, we used this detector to record atmospheric Cherenkov radiation produced by cosmic ray particles showering in the atmosphere. The detected rate of cosmic ray events was consistent with an energy threshold near 1 TeV. The data presented here represent the first detection of atmospheric Cherenkov radiation using solar heliostats viewed from a central tower.     

An Overview of the Disposition of Solar Radiation in the Lower Atmosphere: Connections to the SORCE Mission and Climate Change

Solar radiation is the primary energy source for many processes in Earth's environment and is responsible for driving the atmospheric and oceanic circulation. The integrated strength and spectral distribution of solar radiation is modified from the space-based {Solar {Radiation and {Climate (SORCE) measurements through scattering and absorption processes in the atmosphere and at the surface. Understanding how these processes perturb the distribution of radiative flux density is essential in determining the climate response to changes in concentration of various gases and aerosol particles from natural and anthropogenic sources, as is discerning their associated feedback mechanisms. The past decade has been witness to a tremendous effort to quantify the absorption of solar radiation by clouds and aerosol particles via airborne and space-based observations. Vastly improved measurement and modeling capabilities have enhanced our ability to quantify the radiative energy budget, yet gaps persist in our knowledge of some fundamental variables. This paper reviews some of the many advances in atmospheric solar radiative transfer as well as those areas where large uncertainties remain. The SORCE mission's primary contribution to the energy budget studies is the specification of the solar total and spectral irradiance at the top of the atmosphere.     

Characteristics of the radiation emerging from the top of a Rayleigh atmosphere—II Total upward flux and albedo

The solution of the equation of radiative transfer in a medium exhibiting Rayleigh scattering, as developed by S. Chandrasekhar, has been used for an extensive series of computations(³) of the characteristics of the scattered and diffusely reflected radiation emerging from the top of an atmospheric model which corresponds in many respects to the sunlit portion of the earth's atmosphere. The first part of this two-part discussion dealt with the intensity, degree of polarization, plane of polarization and the neutral points of the emergent light as functions of sun elevation, direction in the downward hemisphere, optical thickness of the model atmosphere and reflectivity of the underlying surface. This second part is concerned with the upward flux obtained by an integration of the intensity over the entire hemisphere, for the incident radiation (a) being independent of wavelength or (b) having the spectral distribution of the extra-terrestrial solar radiation. Integration with respect to wavelength in the latter case, together with an approximation for the sphericity of the atmosphere, yields a value of 7.6 per cent for the earth's planetary albedo due to scattering by the clear atmosphere. An approximation for ozone absorption decreases the computed albedo to 6.9 per cent.     

Energetic Neutral Atoms (ENA) at Mars: Properties of the hydrogen atoms produced upstream of the martian bow shock and implications for ENA sounding technique around non-magnetized planets

We have studied the interaction of fast solar wind hydrogen atoms with the martian atmosphere by a three-dimensional Monte Carlo simulation. These energetic neutral hydrogen atoms, H-ENAs, are formed upstream of the martian bow shock. Both H-ENAs scattered and non-scattered from the martian atmosphere/exosphere were studied. The colliding H-ENAs were found to scatter both to the dayside and nightside. On the dayside they contribute to the so-called H-ENA albedo. On the nightside the heated and scattered hydrogen atoms were found also in the martian wake. The density, the energy distribution function and the direction of the velocity of H-ENAs on the nightside are presented. The present study describes a novel “ENA sounding” technique in which energetic neutral atoms are used to derive information of the properties of planetary exosphere and atmosphere in a similar manner as the solar wind photons are used to derive atmospheric densities by measuring the scattered UV light. A detailed study of the direction and energy of the scattered and non-scattered H-ENAs suggest that the ENA sounding is a method to study the interaction between the planetary atmosphere and the solar wind and to monitor the density, and likely also the magnetization, of the planetary upper atmosphere. Already present-day ENA instrument should be capable to detect the analyzed particle fluxes.     

Status and Results from the HEGRA Air Shower Experiment

The HEGRA experiment investigates the relativistic (non thermal) universe by measuring air showers initiated in the atmosphere by photons and charged particles from the cosmos, using the imaging atmospheric Cherenkov light technique (> 500 GeV) and the showerfront sampling technique (> 20 TeV). New results concerning the search and study of gamma ray sources above 500 GeV, the measurement of the energy spectrum and chemical composition of cosmic rays above 1 PeV and the search of TeV counterparts of Gamma Ray Bursts are briefly described. This revised version was published online in July 2006 with corrections to the Cover Date.     

Examining the feasibility of air shower core-location from polarization properties of associated Cherenkov light

Based on this exploratory investigation involving CORSIKA simulation code generated Cherenkov photons and a linearly polarized, hypothetical photon beam, we make a case here for exploiting polarization properties of atmospheric Cherenkov events for providing an independent method for locating air-shower cores by a TACTIC-like array of atmospheric Cherenkov telescopes. Preliminary results based on simulations indicate that for a 3 TeV γ-ray having ∼30% degree of polarization for its associated Cherenkov light at a core distance of ∼100 m, core location can be found with an error of ∼27 m. Deceased This revised version was published online in July 2006 with corrections to the Cover Date.     

Small-scale trench in the north polar region of Mars: Evolution of surface frost and ground ice concentration

In this paper we attempt to answer the question, how formation of a small-scale trench in the martian regolith affects local distribution of the subsurface ice. We are especially interested in the consequences of digging a trench to search for buried ice, as has been done during the Phoenix Mars Lander mission. However, the results may be also applicable for natural troughs, or cracks. We present results of simulations of diurnal exchange of water between the regolith and the atmosphere. Our model includes the heat and vapor migration in the regolith surrounding the trench, as well as formation of diurnal frost. We take into account scattering of light in the atmosphere and on the trench facets, as well as changes of atmospheric humidity on diurnal and seasonal time scales. Our calculations show, that the measurements of ice content in a sample obtained within one, or two days from the beginning of digging should not be affected. However, on somewhat longer time scale at the south facing site of the trench the regolith can be significantly depleted from ice. This effect should be taken into account if the excavation and taking samples from different depths will be performed in stages separated in time by a month, or more.     

Characteristics of the radiation emerging from the top of a rayleigh atmosphere—I Intensity and polarization

S. Chandrasekhar's solution of the equation of radiative transfer m a medium scattering according to Rayleigh's law has been used for an extensive series of computations of the scattered and diffusely reflected radiation directed outward from the top of a plane stratified model of the atmosphere. The results are discussed in two parts, the first of which deals with the intensity, degree of polarization, plane of polarization and the neutral points of the scattered and diffusely reflected radiation. Part II (see following article) will be a discussion of the upward flux, integrated over the hemisphere, as a function of wavelength of the incident radiation, and of the albedo of the atmospheric model for an energy distribution of the incident light similar to that of the extra-terrestrial solar spectrum. The results are compared with similar quantities for the diffusely scattered radiation from the sunlit sky.     

On the optical studies of the atmospheric water vapour from the surface of Mars

Remote observations of the atmospheric water vapour from the Mars orbit were usually carried out to study its global distribution and variability. Measurements of the water vapour abundance onboard the landers have recently become an important complement to the orbital sounding. Narrow-band filter photometry and spectroscopy of the solar radiation from the surface of the planet proved to be a powerful tool in the study of atmospheric water. The Imager for Mars Pathfinder (IMP) was the first instrument to measure its amount from the surface. The Surface Stereo Imager (SSI) onboard the Mars Polar Lander (MPL) was to follow but the spacecraft was lost at landing. Nevertheless significant expertise in the optical measurements of atmospheric H₂O was gained during these missions. This paper summarizes this experience emphasizing the radiative transfer aspects of the problem. The results of this study could be of importance for future missions to Mars.     

The influence of forward-scattered light in transmission measurements of (exo)planetary atmospheres

The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the planetary atmosphere, the amount of transmitted flux not only depends on the total optical thickness of the slant light path that is probed, but also on the amount of forward-scattering by the scattering particles. Here, we present results of calculations with a three-dimensional Monte Carlo code that simulates the transmitted flux during occultations or transits. For isotropically scattering particles, like gas molecules, the transmitted flux appears to be well-described by the total atmospheric optical thickness. Strongly forward-scattering particles, however, such as commonly found in atmospheres of Solar System planets, can increase the transmitted flux significantly. For exoplanets, such added flux can decrease the apparent radius of the planet by several scale heights, which is comparable to predicted and measured features in exoplanet transit spectra. We performed detailed calculations for Titan’s atmosphere between 2.0 and 2.8 μm and show that haze and gas abundances will be underestimated by about 8% if forward-scattering is ignored in the retrievals. At shorter wavelengths, errors in the gas and haze abundances and in the spectral slope of the haze particles can be several tens of percent, also for other Solar System planetary atmospheres. We also find that the contribution of forward-scattering can be fairly well described by modelling the atmosphere as a plane-parallel slab. This potentially reduces the need for a full three-dimensional Monte Carlo code for calculating transmission spectra of atmospheres that contain forward-scattering particles.     

One-dimensional variational retrieval of refractivity profiles by GPS occultation

In recent years, the technique of Global Positioning System/Low Earth Orbit (GPS/LEO) radio occultation has opened up a new route for exploring the earth atmosphere. With this technique, the vertical profiles of pressure, temperature, water vapor of the earth atmosphere can be retrieved from GPS/LEO occultation data, thus making the latter a potentially valuable data resource for meteorological and atmospheric sciences. The technique may effectively improve the quality of the meteorological profiles, and so improve the current numerical weather forecasting model. In the one-dimensional variational (1DVAR) assimilation retrieval technique in current use, the profiles of atmospheric water vapor and temperature, as well as the sea level pressure are retrieved by assimilating the atmospheric refractivity or bending angle profile from GPS/LEO occultation data. In the system independently developed at Shanghai Observatory, the 1DVAR assimilation of CPS occultation data is accomplished by the Levenberg-Marquardt method with the analytical data from the European Center for Medium Range Weather Forecasting (ECMWF) as the background field, the refractivity profile obtained from CHAMP occultation data as the observational quantity. The result was checked with sounding balloon data near the occultation points.     

An high resolution FDIRC for the measurement of cosmic-ray isotopic abundances

Measurements of the relative abundance of cosmic isotopes and of the energy dependence of their fluxes may clarify our present understanding on the confinement time of charged cosmic rays in the Galaxy. Experimental studies of these propagation clocks have been carried out by balloon and space missions at energies of a few 100 MeV/amu by means of detection techniques based on multiple dE/dx sampling, coupled with a measurement of the energy released in a thick absorber. At larger energies, the isotopic separation of light nuclei (as, for instance,⁹Be/¹⁰Be) can be achieved by combining a precise measurement of the particle’s rigidity with an high resolution determination of its velocity, via the observation of the Cherenkov effect in a radiator.In this paper, we propose the introduction - for the first time in a space experiment - of the DIRC technique (Detection of Internal Reflected Cherenkov light) for the identification of cosmic-ray isotopes. This type of detector has been successfully used in electron-positron colliders for particle identification and in particular for π-K separation. While for particles with unit charge the light yield is a limiting factor, in the case of a nucleus of charge Z the larger photostatistics (due to the Z² dependence of Cherenkov light emission) is the key to reach an adequate angular resolution to provide a mass discrimination for isotopes of astrophysical interest. We report on the early development phase of a DIRC prototype with a focussing scheme (FDIRC) to collect the Cherenkov light onto a detector plane instrumented with a Silicon PhotoMultiplier (SiPM) array.     

Effects of the Absorption of Solar XUV Radiation by the Earth's Upper Atmosphere at Altitudes of 100–500 km in the X-ray Solar Images Obtained Onboard the Coronas-I (TEREK Telescope) and Coronas-F (SPIRIT X-ray Complex) Satellites

We consider the effects of the absorption of solar XUV radiation by the Earth's atmosphere that were observed in the solar images obtained with the TEREK-K telescope onboard the Coronas-I satellite in May–June 1994 at low solar activity and with the SPIRIT instrumentation onboard the Coronas-F satellite in October–November 2001 at maximum solar activity. The solar images were recorded during the satellite occultation: in the 175- and 304-A spectral ranges onboard Coronas-I with the TEREK-K telescope and in the 175-, 304-, and 8.42-A ranges onboard Coronas-F with the SPIRIT instrumentation. Based on the XUV solar images obtained during atmospheric sounding, apart from the total absorption, we can determine the direction of the atmospheric density gradient and study the local absorption variations with altitude on spatial scales of less than 1 km. The described method can significantly supplement the data obtained in studies of the upper atmosphere by the methods of mass spectrometry, incoherent radar scattering, and the drag of orbital spacecraft.     

Sources of error in the photoclinometric determination of planetary topography: A reappraisal

The technique of photoclinometry has frequently been used to determine planetary topography without proper consideration of possible sources of error. Previous studies of error sources have been limited in extent and have overlooked the importance of factors such as atmospheric scattering and the choice of a surface photometric function. This paper adopts a thorough and more direct approach to error analysis, whereby known topography is compared with photoclinometric profiles derived from synthetic quantised reflectance scans.Instrumental and geometric sources of error are found to exert a minimal influence on profiles in practice, provided that sufficient care is taken in the selection of images and the extraction of scans from those images. Environmental factors — relating to the scattering properties of the surface and, if present, atmosphere — are far more important. It is found that a simple Lommel-Seeliger law is unlikely to be appropriate to the majority of planetary terrains, given its inability to model the effects of multiple scattering or unresolved macroscopic roughness. It is further demonstrated that a Minnaert function or combination of Lommel-Seeliger and Lambert laws may empirically compensate for the first of these phenomena but not the second; in this respect, Hapke's equation is a far superior model of surface optical properties. In the case of an atmosphere, the need to correct for scattering by aerosols or suspended dust becomes more acute as atmospheric opacity increases and as particle scattering becomes more forward-biased. To perform this correction, a model for the combined reflectance of surface and atmosphere must be used when deriving profiles.Two case studies — of a small impact crater on Triton and a dust-mantled basaltic lava flow on Mars - are presented here. Regarding the latter, the implications that errors in photoclinometric flow thickness measurements have for inferred lava rheology are examined. Conservative estimates of errors in yield strength and apparent viscosity easily exceed 100% when one of the simplest photometric models possible — a Lommel-Seeliger law — is used to derive a profile.In the light of these findings, strategies are suggested for improving the results obtained from photoclinometry in the future.