We address an interesting question raised by Dos Santos Ferreira, Kenig and Salo [11] about regions Rg ⊂ C for which there can be uniform L 2n n+2 → L 2n n−2 resolvent estimates for ∆g + ζ, ζ ∈ Rg, where ∆g is the Laplace-Beltrami operator with metric g on a given compact boundaryless Riemannian manifold of dimension n ≥ 3. This is related to earlier work of Kenig, Ruiz and the third author [19] for the Euclidean Laplacian, in which case the region is the entire complex plane minus any disc centered at the origin. Presently, we show that for the round metric on the sphere, S n , the resolvent estimates in [11], involving a much smaller region, are essentially optimal. We do this by establishing sharp bounds based on the distance from ζ to the spectrum of ∆ S n . In the other direction, we also show that the bounds in [11] can be sharpened logarithmically for manifolds with nonpositive curvature, and by powers in the case of the torus, T n = R n /Z n , with the flat metric. The latter improves earlier bounds of Shen [22]. The work of [11] and [22] was based on Hadamard parametrices for (∆g + ζ) −1 . Ours is based on the related Hadamard parametrices for cos t −∆g, and it follows ideas in [26] of proving L p -multiplier estimates using small-time wave equation parametrices and the spectral projection estimates from [25]. This approach allows us to adapt arguments in Bérard [3] and Hlawka [15] to obtain the aforementioned improvements over [11] and [22]. Further improvements for the torus are obtained using recent techniques of the first author [5] and his work with Guth [7] based on the multilinear estimates of Bennett, Carbery and Tao [2]. Our approach also allows us to give a natural necessary condition for favorable resolvent estimates that is based on a measurement of the density of the spectrum of −∆g, and, moreover, a necessary and sufficient condition based on natural improved spectral projection estimates for shrinking intervals, as opposed to those in [25] for unit-length intervals. We show that the resolvent estimates are sensitive to clustering within the spectrum, which is not surprising given Sommerfeld's original conjecture [32] about these operators.