We extend weak Drazin inverses of square matrices to rectangular, and investigate some properties of these inverses. Applications to systems of matrix equations are also given.