Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies.

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $\gamma$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $\gamma$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$\sigma$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

We seek to estimate the average level of MHD turbulence near coronal mass ejection (CME) fronts as they propagate from the Sun to the Earth. We examine the cosmic ray data from the GRAPES-3 tracking muon telescope at Ooty, together with the data from other sources for three well observed Forbush decrease events. Each of these events are associated with frontside halo Coronal Mass Ejections (CMEs) and near-Earth magnetic clouds. In each case, we estimate the magnitude of the Forbush decrease using a simple model for the diffusion of high energy protons through the largely closed field lines enclosing the CME as it expands and propagates from the Sun to the Earth. We use estimates of the cross-field diffusion coefficient $D_{\perp}$ derived from published results of extensive Monte Carlo simulations of cosmic rays propagating through turbulent magnetic fields. Our method helps constrain the ratio of energy density in the turbulent magnetic fields to that in the mean magnetic fields near the CME fronts. This ratio is found to be $\sim$ 2% for the 11 April 2001 Forbush decrease event, $\sim$ 6% for the 20 November 2003 Forbush decrease event and $\sim$ 249% for the much more energetic event of 29 October 2003.

Earth's bow shock is known to produce non-thermal electrons which are generally observed as a `spike' in their flux profile. Here, in this paper, we present an analysis of electron and whistler wave properties for a quasi-perpendicular shock crossing that is supercritical, but subcritical to the so-called whistler critical Mach number, M$^w_{\rm crit}$, above which whistler waves cannot propagate upstream. We have found that the amplitudes of whistler waves increased exponentially as a function of time prior to the shock encounter, while the suprathermal ($>$ 2 keV) electron flux similarly increased with time, although with differing $e$-folding time scales. Comparison of the electron energy spectrum measured within the ramp with predictions from diffusive shock acceleration theory was poor, but the variation of pitch angle distribution showed scattering of non-thermal electrons in the upstream region. While not finding a specific mechanism to account for the electron diffusion, we suggest that the whistlers seen probably account for the differences observed between this `gradual' event and the `spike' events seen at shocks with no upstream whistlers.