We consider, for the first time, the effects of strong capacitive and inductive coupling between radio frequency superconducting quantum interference devices (rf SQUIDs) in an overlapping metamaterial geometry when driven by rf flux at and near their self-resonant frequencies. The equations of motion for the gauge-invariant phases on the Josephson junctions in each SQUID are set up and solved. Our model accounts for the high-frequency displacement currents through capacitive overlap between the wiring of SQUID loops. We begin by modeling two overlapping SQUIDs and studying the response in both the linear and nonlinear high-frequency driving limits. By exploring a sequence of more and more complicated arrays, the formalism is eventually extended to the N×N×2 overlapping metamaterial array, where we develop an understanding of the many ( 8N2−8N+3 ) resulting resonant modes in terms of three classes of resonances. The capacitive coupling gives rise to qualitatively new self-resonant responses of rf SQUID metamaterials, and is demonstrated through analytical theory, numerical modeling, and experiment in the 10–30 GHz range on capacitively and inductively coupled rf SQUID metamaterials.