In this paper, we study existence and uniqueness of strong as well as weak solutions for general time fractional Poisson equations. We show that there is an integral representation of the solutions of time fractional Poisson equations with zero initial values in terms of semigroup for the infinitesimal spatial generator L and the corresponding subordinator associated with the time fractional derivative. This integral representation has an integral kernel q(t, x, y), which we call the fundamental solution for the time fractional Poisson equation, if the semigroup for L has an integral kernel. We further show that q(t, x, y) can be expressed as a time fractional derivative of the fundamental solution for the homogenous time fractional equation under the assumption that the associated subordinator admits a conjugate subordinator. Moreover, when the Laplace exponent of the associated subordinator satisfies the weak scaling property and its distribution is self-decomposable, we establish two-sided estimates for the fundamental solution q(t, x, y) through explicit estimates of transition density functions of subordinators. ∞ 0 min{1, x}(−dw(x)) < ∞, a generalized time fractional derivative with weight w is defined by∂ w t f (t) := d dt t 0 whenever the right hand side is well defined. See [8, 13, 24, 27]. In particular, when w(s) = 1 Γ(1−β) s −β for β ∈ (0, 1) (where Γ(t) = ∞ 0 s t−1 e −s ds is the Gamma function), ∂ w t f is just the Caputo derivative of order β, i.e., ∂ β t f (t) := 1 Γ(1 − β) d dt t 0 Zhen-Qing Chen