Source localization and geoacoustic inversion performance depends on the location of receiving phones, their spacing, their number, and the array aperture. Performance suffers when practical design issues limit the capabilities of the array. In this work, sparse arrays are augmented leading to the generation of virtual arrays for better sampling of the ocean and, thus, improved estimation performance. To that effect, Gaussian Processes are employed, which are shown to create high-fidelity field “measurements” at virtual, densely spaced hydrophones. Kernel functions are key building blocks in the implementation of Gaussian Processes, as they quantify the field coherence at neighboring spatial points. Functions of interest are the squared exponential, Matern, and plane wave kernels. We validate our method with application to synthetic data as well as data collected during the Seabed Characterization Experiment conducted in 2022. [Work supported by ONR.]