Constrained optimization in seismic reflection tomography: a Gauss-Newton augmented Lagrangian approach