Operator-matched spatial regression

A standard spatial regression model $Y = X^\topβ+ U$, with covariates $X$ and a Gaussian random field $U$ with Matérn covariance, decouples covariate effects from spatial dependence. In several applications this is mechanistically implausible, as covariates likely drive the dynamics generating the residual covariance, and it is statistically problematic, since the generalised least-squares estimator of $β$ may be biased by spatial self-confounding. We propose operator-matched spatial regression, where covariates enter both pointwise and as deterministic forcings in a stochastic partial differential equation (SPDE) that generates the random effect, partitioning the apparent effect into local and operator-mediated components. The formulation arises as the steady state of a transport–relaxation equation forced by $X$, giving the parameters physical interpretations. We show that self-confounding is resolved by construction and that maximum likelihood inference for the regression coefficients, implementable through standard covariance-based methods, is calibrated when the operator parameters are known and, under expanding-domain asymptotics, when they are jointly estimated. A computationally efficient finite element implementation is introduced and justified via bounds on the discretisation error. Applications to temperature–altitude regression and air-quality monitoring show that the method matches or outperforms alternatives and produces physically meaningful results.