This study developed a method to represent 2D seismic velocity fields using 2D Haar wavelet series. The method simplifies previous approaches using polynomials or 1D wavelet series. Seismic velocity fields from 5 geological models were approximated using the 2D Haar wavelet algorithm. While the number of coefficients was initially high, generating coefficient averages over subintervals reduced the total coefficients by up to 90%, simplifying the representations. The findings support expanding this work to seismic reflection travel time modeling, ray tracing, tomography, and migration using the wavelet parameterization.