Shape optimization of a staggered herringbone groove micromixer using three-dimensional Navier-Stokes analysis has been carried using Kriging model. The analysis of the degree of mixing is performed by the calculation of spatial data statistics. The calculation of the variance of the mass fraction at various nodes on a plane in the channel is used to quantify mixing. A numerical optimization technique with Kriging model is applied to optimize the shape of the grooves on a single wall of the channel. Three design variables, namely, the ratio of groove width to groove pitch, the ratio of the groove depth to channel height ratio and the angle of the groove, are selected for optimization. A mixing index is used as the objective function. The results of the optimization show that the mixing is very sensitive to the shape of the groove which can be used in controlling mixing in microdevices.