This article is dedicated to the study of an SPDE of the form Lu(t, x) = σ(u(t, x))Ż(t, x) t > 0, x ∈ O with zero initial conditions and Dirichlet boundary conditions, where σ is a Lipschitz function, L is a second-order pseudo-differential operator, O is a bounded domain in R d , andŻ is an α-stable Lévy noise with α ∈ (0, 2), α = 1 and possibly non-symmetric tails. To give a meaning to the concept of solution, we develop a theory of stochastic integration with respect to Z, by generalizing the method of [11] to higher dimensions and non-symmetric tails. The idea is to first solve the equation with "truncated" noiseŻ K (obtained by removing from Z the jumps which exceed a fixed value K), yielding a solution u K , and then show that the solutions u L , L > K coincide on the event t ≤ τ K , for some stopping times τ K ↑ ∞ a.s. A similar idea was used in [22] in the setting of Hilbert-space valued processes. A major step is to show that the stochastic integral with respect to Z K satisfies a p-th moment inequality, for p ∈ (α, 1) if α < 1, and p ∈ (α, 2) if α > 1. This inequality plays the same role as the Burkholder-Davis-Gundy inequality in the theory of integration with respect to continuous martingales.MSC 2000 subject classification: Primary 60H15; secondary 60H05, 60G60