Existence of mild solutions for a system of partial differential equations with time-dependent generators

We give sufficient conditions for global existence of positive mild solutions for the weak coupled system: \begin{eqnarray*} \frac{\partial u_{1}}{\partial t} &=&\rho_{1}t^{\rho_{1}-1}\Delta_{\alpha_{1}}u_{1}+t^{\sigma_{1}}u_{2}^{\beta_{1}},\ \ u_{1}\left(0\right) =\varphi_{1}, \\ \frac{\partial u_{2}}{\partial t} &=&\rho_{2}t^{\rho_{2}-1}\Delta_{\alpha_{2}}u_{2}+t^{\sigma_{2}}u_{1}^{\beta_{2}},\ \ u_{2}\left(0\right) =\varphi_{2}, \end{eqnarray*} where $\Delta_{\alpha_{i}}$ is a fractional Laplacian, $0<\alpha_{i}\leq 2,\ \beta_{i}>1,\ \rho_{i}>0,\ \sigma_{i}>-1\$are constants and the initial data $\varphi_{i}$ are positive, bounded and integrable functions.