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 located in Ooty, India records 4x10^9 muons daily. These muons are produced by interaction of primary cosmic rays (PCRs) in the atmosphere. The high statistics of muons enables GRAPES-3 to make precise measurement of various sun-induced phenomenon including coronal mass ejections (CME), Forbush decreases, geomagnetic storms (GMS) and atmosphere acceleration during the overhead passage of thunderclouds. However, the understanding and interpretation of observed data requires Monte Carlo (MC) simulation of PCRs and subsequent development of showers in the atmosphere. CORSIKA is a standard MC simulation code widely used for this purpose. However, these simulations are time consuming as large number of interactions and decays need to be taken into account at various stages of shower development from top of the atmosphere down to ground level. Therefore, computing resources become an important consideration particularly when billion of PCRs need to be simulated to match the high statistical accuracy of the data. During the GRAPES-3 simulations, it was observed that over 60% of simulated events don't really reach the Earth's atmosphere. The geomagnetic field (GMF) creates a threshold to PCRs called cutoff rigidity Rc, a direction dependent parameter below which PCRs can't reach the Earth's atmosphere. However, in CORSIKA there is no provision to set a direction dependent threshold. We have devised an efficient method that has taken into account of this Rc dependence. A reduction by a factor ~3 in simulation time and ~2 in output data size was achieved for GRAPES-3 simulations. This has been incorporated in CORSIKA version v75600 onwards. Detailed implementation of this along the potential benefits are discussed in this work.

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.

Aims. Forbush decrease (FD) is a transient decrease followed by a gradual recovery in the observed galactic cosmic ray intensity. We seek to understand the relationship between the FDs and near-Earth interplanetary magnetic field (IMF) enhancements associated with solar coronal mass ejections (CMEs). Methods. We use muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and magnitude > 0.25%. We use IMF data from ACE/WIND spacecrafts. We look for correlations between the FD profile and that of the one hour averaged IMF. We ask if the diffusion of high energy protons into the large scale magnetic field is the cause of the lag observed between the FD and the IMF. Results. The enhancement of the IMF associated with FDs occurs mainly in the shock-sheath region, and the turbulence level in the magnetic field is also enhanced in this region. The observed FD profiles look remarkably similar to the IMF enhancement profiles. The FDs typically lag the IMF enhancement by a few hours. The lag corresponds to the time taken by high energy protons to diffuse into the magnetic field enhancement via cross-field diffusion. Conclusions. Our findings show that high rigidity FDs associated with CMEs are caused primarily by the cumulative diffusion of protons across the magnetic field enhancement in the turbulent sheath region between the shock and the CME.

We seek to identify the primary agents causing Forbush decreases (FDs) observed at the Earth in high rigidity cosmic rays. In particular, we ask if such FDs are caused mainly by coronal mass ejections (CMEs) from the Sun that are directed towards the Earth, or by their associated shocks. We use the muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and can be reasonably well associated with an Earth-directed CME and its associated shock. We employ two models: one that considers the CME as the sole cause of the FD (the CME-only model) and one that considers the shock as the only agent causing the FD (the shock-only model). We use an extensive set of observationally determined parameters for both these models. The only free parameter in these models is the level of MHD turbulence in the sheath region, which mediates cosmic ray diffusion (into the CME, for the CME-only model and across the shock sheath, for the shock-only model). We find that good fits to the GRAPES-3 multi-rigidity data using the CME-only model require turbulence levels in the CME sheath region that are only slightly higher than those estimated for the quiet solar wind. On the other hand, reasonable model fits with the shock-only model require turbulence levels in the sheath region that are an order of magnitude higher than those in the quiet solar wind. This observation naturally leads to the conclusion that the Earth-directed CMEs are the primary contributors to FDs observed in high rigidity cosmic rays.

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.