Not trivial to estimate this. We don't even have that one simple ground node as depicted in the equivalent circuit model.
Angelo Scuderi et al provide some useful equations for equivalent circuit modelling in "Integrated Inductors and Transformers" (CRC Press).
In this book, you will find : Rsub*Csub = epsr_sub * rho_sub (Pfost, Rein, Holztwarth 1996). However, when I tried to apply that for real world SiGe inductors with non-homogeneous substrate doping (EPI layer) that equation didn't work well.
My approach is to extract the equivalent circuit model parameters from EM simulated (or measured) inductor data.
https://muehlhaus.com/products/equivalent-circuit-model-fit-for-rfic-inductors