The temperature dependence, between 10 and 120 K, of electron spin-lattice relaxation at X-band was analyzed for a series of eight pyrrolate-imine complexes and for ten other copper(II) complexes with varying ligands and geometry including copper-containing prion octarepeat domain and S100 type proteins. The geometry of the CuN4 coordination sphere for pyrrolate-imine complexes with R=H, methyl, n-butyl, diphenylmethyl, benzyl, 2-adamantyl, 1-adamantyl, and tert-butyl has been shown to range from planar to pseudo-tetrahedral. The fit to the recovery curves was better for a distribution of values of T1 than for a single time constant. Distributions of relaxation times may be characteristic of Cu(II) in glassy solution. Long-pulse saturation recovery and inversion recovery measurements were performed. The temperature dependence of spin-lattice relaxation rates was analyzed in terms of contributions from the direct process, the Raman process, and local modes. It was necessary to include more than one process to fit the experimental data. There was a small contribution from the direct process at low temperature. The Raman process was the dominant contribution to relaxation between about 20 and 60 K. Debye temperatures were between 80 and 120 K. For samples with similar Debye temperatures the coefficient of the Raman process tended to increase as gz increased, as expected if modulation of spin-orbit coupling is a major factor in relaxation rates. Above about 60 K local modes with energies in the range of 260-360 K (180-250 cm-1) dominated the relaxation. For molecules with similar geometry, relaxation rates were faster for more flexible molecules than for more rigid ones. Relaxation rates for the copper protein samples were similar to rates for small molecules with comparable coordination spheres. At each temperature studied the range of relaxation rates was less than an order of magnitude. The spread was smaller between 20 and 60 K where the Raman process dominates, than at higher temperatures where local modes dominate the relaxation. Spin echo dephasing time constants, Tm, were calculated from two-pulse spin echo decays. Near 10 K Tm was dominated by proton spins in the surroundings. As temperature was increased motion and spin-lattice relaxation made increasing contributions to Tm. Near 100 K spin-lattice relaxation dominated Tm.