Heterogeneity of dynamics plays a vital role in membrane function, but the methods for quantifying this heterogeneity are still being developed. Here we examine membrane dynamical heterogeneity via molecular simulations of a single-component dipalmitoylphosphatidylcholine (DPPC) lipid bilayer using the MARTINI force field. We draw upon well established analysis methods developed in the study of glass-forming fluids and find significant changes in lipid dynamics between the fluid (L-alpha), and gel (L-beta) phases. In particular, we distinguish two mobility groups in the more ordered L-beta phase: (i) lipids that are transiently trapped by their neighbors and (ii) lipids with displacements on the scale of the intermolecular spacing. These distinct mobility groups spatially segregate, forming dynamic clusters that have characteristic time (1-2 mu s) and length (1-10 nm) scales comparable to those of proteins and other biomolecules. We suggest that these dynamic clusters could couple to biomolecules within the membrane and thus may play a role in many membrane functions. In the equilibrium membrane, lipid molecules dynamically exchange between the mobility groups, and the resulting clusters are not associated with a thermodynamic phase separation. Dynamical clusters having similar characteristics arise in many other condensed phase materials, placing membranes in a broad class of materials with strong intermolecular interactions.