Abstract. Let M • be a complete noncompact manifold and g an asymptotically conic manifold on M • , in the sense that M • compactifies to a manifold with boundary M in such a way that g becomes a scattering metric on M . A special case that we focus on is that of asymptotically Euclidean manifolds, where the induced metric at infinity is equal to the standard metric on S n−1 ; such manifolds have an end that can be identified with R n \ B(R, 0) in such a way that the metric is asymptotic in a precise sense to the flat Euclidean metric. We analyze the asymptotics of the resolvent kernel (P + k 2 ) −1 where P = ∆g + V is the sum of the positive Laplacian associated to g and a real potential function V that is smooth on M and vanishes to third order at the boundary (i.e. decays to third order at infinity on M • ). We show that on a blown up version of M 2 × [0, k 0 ] the resolvent kernel is conormal to the lifted diagonal and polyhomogeneous at the boundary, and we are able to identify explicitly the leading behaviour of the kernel at each boundary hypersurface. Using this we show that the Riesz transform of P is bounded on L p (M • ) for 1 < p < n, and that this range is optimal if V ≡ 0 or if M has more than one end. The result with V ≡ 0 is new even when M • = R n , g is the Euclidean metric and V is compactly supported. When V ≡ 0 with one end, the range of p becomes 1 < p < pmax where pmax > n depends explicitly on the first non-zero eigenvalue of the Laplacian on the boundary ∂M .Our results hold for all dimensions ≥ 3 under the assumption that P has neither zero modes nor a zero-resonance. In a follow-up paper we shall analyze the same situation in the presence of zero modes and zero-resonances.