We perform a set of 36 nonprecessing black-hole binary simulations with spins either aligned or counteraligned with the orbital angular momentum in order to model the final mass, spin, and recoil of the merged black hole as a function of the individual black hole spin magnitudes and the mass ratio of the progenitors. We find that the maximum recoil for these configurations is Vmax = 526 ± 23 km s −1 , which occurs when the progenitor spins are maximal, the mass ratio is qmax = m1/m2 = 0.623 ± 0.038, the smaller black-hole spin is aligned with the orbital angular momentum, and the larger black-hole spin is counteraligned (α1 = −α2 = 1). This maximum recoil is about 80 km s −1 larger than previous estimates, but most importantly, because the maximum occurs for smaller mass ratios, the probability for a merging binary to recoil faster than 400 km s −1 can be as large as 17%, while the probability for recoils faster than 250 km s −1 can be as large as 45%. We provide explicit phenomenological formulas for the final mass, spin, and recoil as a function of the individual BH spins and the mass difference between the two black holes. Here we include terms up through fourth-order in the initial spins and mass difference, and find excellent agreement (within a few percent) with independent results available in the literature. The maximum radiated energy is E rad /m ≈ 11.3% and final spin α max rem ≈ 0.952 for equal mass, aligned maximally spinning binaries.