Air showers recorded by the Haverah Park Array during the years 1974 - 1987 have been re-analysed. For the original estimate of the energy spectrum, a relationship between the ground parameter $\rho$(600) and the primary energy, as determined by Hillas in the 1970s, was used. Here we describe the energy spectrum obtained using the QGSJET98 interaction model in the CORSIKA Monte Carlo code, together with GEANT to simulate the detailed detector response to ground particles. A new energy spectrum in the range 3 $\times$ 10$^{17}$ eV to 4 $\times$ 10$^{18}$ eV is presented.

At the Haverah Park Array a number of air shower observables were measured that are relevant to the determination of the mass composition of cosmic rays. In this paper we discuss measurements of the risetime of signals in large area water-Cherenkov detectors and of the lateral distribution function of the water-Cherenkov signal. The former are used to demonstrate that the CORSIKA code, using the QGSJET98 model, gives an adequate description of the data with a low sensitivity, in this energy range, to assumptions about primary mass. By contrast the lateral distribution is sufficiently well measured that there is mass sensitivity. We argue that in the range 0.2-1.0 EeV the data are well represented with a bi-modal composition of 34+-2 % protons and the rest iron. We also discuss the systematic errors induced by the choice of hadronic model.

The arrival directions of extremely high energy cosmic rays (EHECR) above $4\times10^{19}$ eV, observed by four surface array experiments in the northern hemisphere,are examined for coincidences from similar directions in the sky. The total number of cosmic rays is 92.A significant number of double coincidences (doublet) and triple coincidences (triplet) are observed on the supergalactic plane within the experimental angular resolution. The chance probability of such multiplets from a uniform distribution is less than 1 % if we consider a restricted region within $\pm 10^{\circ}$ of the supergalactic plane. Though there is still a possibility of chance coincidence, the present results on small angle clustering along the supergalactic plane may be important in interpreting EHECR enigma. An independent set of data is required to check our claims.

We give an overview of the current status and scientific goals of VERITAS, a proposed hexagonal array of seven 10 m aperture imaging Cherenkov telescopes. The selected site is Montosa Canyon (1390 m a.s.l.) at the Whipple Observatory, Arizona. Each telescope, of 12 m focal length, will initially be equipped with a 499 element photomultiplier camera covering a 3.5 degree field of view. A central station will initiate the readout of 500 MHz FADCs upon receipt of multiple telescope triggers. The minimum detectable flux sensitivity will be 0.5% of the Crab Nebula flux at 200 GeV. Detailed simulations of the array's performance are presented elsewhere at this meeting. VERITAS will operate primarily as a gamma-ray observatory in the 50 GeV to 50 TeV range for the study of active galaxies, supernova remnants, pulsars and gamma-ray bursts.

In this Letter we examine the arrival directions of the most energetic cosmic rays (E > 2 * 10^19 eV) detected by several air shower experiments. We find that data taken by different air shower arrays show positive correlations, indicating a non--uniform arrival direction distribution. We also find that the events with energy $ > 4 * 10^19 eV exhibit a correlation with the general direction of the supergalactic plane, where a large number of potential sources is located. If confirmed by data from other experiments our results would support models for the extragalactic origin of the highest energy cosmic rays.