“…Spatiotemporal variation among values is incorporated through three-dimensional functions based on spline smoothing 22 – 27 , accounting for spatial and temporal autocorrelation and improving prediction and inference 28 , 29 . Similar structures have been widely used in the related spatially-varying coefficient models 30 , 31 and temporally-varying coefficient models 32 , including applications in forestry 33 , ecology 34 , and economics 35 , 36 . Spline-based methods, which are often used to estimate in SIVCMs 37 , 38 , are more robust for spatially correlated data and do not require the spatial variation to be specified by a functional form 39 , 40 .…”