A kinetic model was developed, based on the power law of growth and Avrami's nuclei growth concepts, to describe the heterogeneous nature of pulping kinetics, taking into account the effects of chemical concentration and temperature. The general form of the kinetic equation is first order with a timedependent rate coefficient. The model was statistically tested using published data obtained from delignification of lignocellulosic materials (hemp woody core, giant reed, wheat straw, cottonwood, bamboo and flax fibres). The activation energies of the delignification range from 71 to 136 kJ/mol. The p-values obtained from the regression analysis are significantly small indicating that all the estimates of the model parameters were significant with very high levels of confidence. The correlation coefficients R 2 for these models range from 0.76 to 0.98.