Up to now, laser damage growth on the exit surface of fused silica optics has been mainly considered as exponential, the growth coefficient depending essentially on fluence. From experiments with large beams carried out at 351 nm under nanosecond pulses, a statistical analysis is conducted leading to a refined representation of the growth. The effect of several parameters has been taken into account to describe precisely the growth phenomenon. The two principal parameters proved to be the mean fluence and the size of the damage sites. Nevertheless, contributions of other parameters have been estimated too: the number of neighbors around the damage site, the shot number, etc. From experimental results, a model smoothed on a statistical approach is developed that permits the description of a complete sequence of growth. To evaluate the relevance of the modeling approach, the occluded area estimated from modeling is compared with the ones experimentally measured. For this purpose, numerical growth methods have been developed too. It is shown that the approach outlined is appropriate for a more precise description of the growth.