$\newcommand\th\theta\newcommand\R{\mathbb R}$Suppose that we have a parametric setting, that is, the unknown distribution $P$ belongs to a known parametric family $(P_\th)$ of distributions parameterized by a sufficiently low-dimensional parameter $\th$; this may be the case if $P$ is Gaussian. Then you can get (say) a maximum likelihood estimate (MLE) $\hat\th_n$ of $\th$ and estimate $R$ by $$\hat R_n:=\int_H p_{\hat\th_n}\,ds,$$ where $p_{\hat\th_n}$ is the density of $P_{\hat\th_n}$. The MLE is usually consistent, so that you will have $\hat\th_n\to\th$ in probability (as $n\to\infty$). If now $p_\th$ is continuous on $\th$ uniformly on compact subsets of $\R^m$, you will get $p_{\hat\th_n}\to p_\th$ uniformly on compact subsets of $\R^m$. If, moreover, you can control the tails of the densities $p_\th$ so as to have their local uniform integrability (as you would have in the Gaussian case), then you will end up with the desired conclusion $\hat R_n\to R$ in probability.