The aim of this study is to develop robust guidance laws for the control motion of an underwater autonomous vehicle (UAV) in a three-dimensional (3D) space. The control design is based on the use of Averaged Sub-Gradient (ASG) version of a class of dynamic integral sliding mode (ISM) method being sequentially applied to the subsystems of the complete model realizing, the so-called backstepping (or cascade) approach. The mathematical form of the UAV model induces a backstepping formulation for solving the tracking trajectory problem sequentially for the position, translation velocity, angular velocity and actuators (thrusters) dynamics. The solution of the trajectory tracking problem at each stage implements the ASG-version of the ISM method. This problem is treated as the optimization of a suitable convex (not obligatory strongly convex) cost functional, depending on the tracking error and reaching its minimal value at the origin of the error tracking space. This study shows that the minimization of the proposed functional leads to the optimal tracking regime under the presence of uncertainties in the mathematical model description. A numerical example proves the effectiveness of the suggested robust dynamic controller. The comparison between the obtained trajectory tracking results and the outcomes produced by a set of standard proportional integral derivative (PID) controllers, is presented. The proposed controller exhibits a better tracking of the reference trajectory compared with the PID version, showing a smaller mean square estimation for the tracking error.