The uncertainty in the empirical ground motion prediction models (GMMs) for any region depends on several parameters. In the present work, we apply an artificial neural network (ANN) to design a GMM of peak ground acceleration (PGA) for Kachchh, Gujarat, India, utilizing independent input parameters viz., moment magnitudes, hypocentral distances, focal depths and site proxy (in terms of average seismic shear-wave velocity from the surface to a depth of 30 m (Vs30)). The study has been performed using a PGA dataset consisting of eight engineering seismoscope records of the 2001 Mw7.7 Bhuj earthquake and 237 strong-motion records of 32 significant Bhuj aftershocks of Mw3.3–5.6 (during 2002–2008) with epicentral distances ranging from 1.0 to 288 km. We apply a feed-forward back propagation ANN method with 8 hidden nodes, which is found to be optimal for the selected PGA database and input–output mapping. The standard deviation of the error has been utilized to examine the performance of our model. We also test the ground motion predictability of our ANN model using real recordings of the 2001 Bhuj mainshock, two Mw5.6 Kachchh aftershocks and the 1999 Mw6.4 Chamoli mainshock. The standard deviation of PGA prediction error estimates in log10 units is found to be ± 0.2554. Also, the model predictability of our ANN model suggests a good prediction of the PGA for earthquakes of Mw5.6–7.7, which are occurring in Kachchh, Gujarat, India.