Aftershock sequences of the magnitude M-W =6.4 Bingol earthquake of 1 May, 2003 (Turkey) are studied to analyze the spatial and temporal variability of seismicity parameters of the b value of the frequency-magnitude distribution and the p value describing the temporal decay rate of aftershocks. The catalog taken from the KOERI contains 516 events and one month's time interval. The b value is found as 1.49 +/- 0.07 with Mc =3.2. Considering the error limits, b value is very close to the maximum b value stated in the literature. This larger value may be caused by the paucity of the larger aftershocks with magnitude M-D >= 5.0. Also, the aftershock area is divided into four parts in order to detect the differences in b value and the changes illustrate the heterogeneity of the aftershock region. The p value is calculated as 0.86 +/- 0.11, relatively small. This small p value may be a result of the slow decay rate of the aftershock activity and the small number of aftershocks. For the fitting of a suitable model and estimation of correct values of decay parameters, the sequence is also modeled as a background seismicty rate model. Constant background activity does not appear to be important during the first month of the Bingol aftershock sequences and this result is coherent with an average estimation of pre-existing seismicity. The results show that usage of simple modified Omori law is reasonable for the analysis. The spatial variability in b value is between 1.2 and 1.8 and p value varies from 0.6 to 1.2. Although the physical interpretation of the spatial variability of these seismicity parameters is not straightforward, the variation of b and p values can be related to the stress and slip distribution after the mainshock, respectively. The lower b values are observed in the high stress regions and to a certain extent, the largest b values are related to Holocene alluvium. The larger p values are found in some part of the aftershock area although no slip occurred after the main shock and it is interpreted that this situation may be caused by the alluvium structure of the region. These results indicate that the spatial distribution in b and p values are generally related to the rupture mechanism and material properties of an aftershock area.