We consider the Cauchy problem for the Gross-Pitaevskii (GP) equation. Using the DBAR generalization of the nonlinear steepest descent method of Deift and Zhou we derive the leading order approximation to the solution of the GP in the solitonic region of space time $|x| < 2t$ for large times and provide bounds for the error which decay as $t to infty$ for a general class of initial data whose difference from the non-vanishing background possesss a fixed number of finite moments and derivatives. Using properties of the scattering map for (GP) we derive as a corollary an asymptotic stability result for initial data which are sufficiently close to the N-dark soliton solutions of (GP).