The spectral boundary integral method is popular for simulating fault, fracture, and frictional processes at a planar interface. However, the method is less commonly used to simulate off-fault dynamic fields. Here we develop a spectral boundary integral method for poroelastodynamic solid. The method has two steps: first, a numerical approximation of a convolution kernel and second, an efficient temporal convolution of slip speed and the appropriate kernel. The first step is computationally expensive but easily parallelizable and scalable such that the computational time is mostly restricted by computational resources. The kernel is independent of the slip history such that the same kernel can be used to explore a wide range of slip scenarios. We apply the method by exploring the short-time dynamic and static responses: first, with a simple source at intermediate and far-field distances and second, with a complex near-field source. We check if similar results can be attained with dynamic elasticity and undrained pore-pressure response and conclude that such an approach works well in the near-field but not necessarily at an intermediate and far-field distance. We analyze the dynamic pore-pressure response and find that the P-wave arrival carries a significant pore pressure peak that may be observed in high sampling rate pore-pressure measurements. We conclude that a spectral boundary integral method may offer a viable alternative to other approaches where the bulk is discretized, providing a better understanding of the near-field dynamics of the bulk in response to finite fault ruptures.