The nonlinear interaction of an intense electromagnetic beam with relativistic collisionless unmagnetized plasma is investigated by invoking the reductive perturbation technique, resting on the model of three-dimensional nonlinear Schrödinger (NLS) equation with cubic nonlinearity which incorporates the effects of self-focusing, self-phase modulation, and diffraction on wave propagation. Relying on the derived NLS equation, the occurrence of spatiotemporal modulational instability is investigated in detail.