Based on the concept of self-decomposability, we extend some recent multidimensional Lévy models built using multivariate subordination. Our aim is to construct multivariate Lévy processes that can model the propagation of the systematic risk in dependent markets with some stochastic delay instead of affecting all the markets at the same time. To this end, we extend some known approaches keeping their mathematical tractability, study the properties of the new processes, derive closed-form expressions for their characteristic functions and detail how Monte Carlo schemes can be implemented. We illustrate the applicability of our approach in the context of gas, power and emission markets focusing on the calibration and on the pricing of spread options written on different underlying commodities.