Abstract: In this study, by applying cross-gradient constraints, we joint early arrivals and surface waves to invert the 2D P-wave velocity, S-wave velocity, Q_P and Q_S simultaneously for the near-surface area. In order to improve the efficiency of computation, we propose a method that employs 1D elastic forward modeling and applies 2D Tikhonov regularization to invert the 2D velocity structures and attenuation model. Synthetic tests and real data application illustrate that the method can greatly improve the computational efficiency, and is able to invert reliable velocity and attenuation models simultaneously.