Towards Efficient Modeling and Inference in Multi-Dimensional Gaussian Process State-Space Models