In this paper, two classes of high-order numerical schemes on the time discretization for the solutions of two-dimensional nonlinear Sobolev equations are analyzed. The twolevel Newton linearized difference scheme had been established in the previous literature, in which the error estimates in L-2-norm and H1-norm are obtained, but that in L infinity -norm remains unsolved. The second difference scheme is established based on a threelevel linearized difference technique in temporal direction. The numerical discretization in spatial dimension for both schemes utilizes the compact difference operator. We prove that the orders of convergence for both difference schemes in L infinity-norm are O(tau(2)+ h(1)(4)+ h(2)(4)) based on the energy argument, where tau denotes the temporal step size, h(1) and h(2) denote the spatial step sizes in x-direction and y-direction, respectively. The numerical tests further verify the theoretical results and demonstrate the efficiency of both difference schemes. (c) 2022 IMACS. Published by Elsevier B.V. All rights reserved.