We consider the propagation of high frequency electromagnetic pulses in complex materials with nonlinear polarization. The physical problem is modeled by Maxwell鈥檚 equations in variational form, and well-posedness results are established with respect to probability distributions on the polarization parameters (in a Prohorov metric sense).