The mathematical models and numerical simulations reported here are motivated by satellite observations of horizontal sea surface fluid motions that are closely coordinated with the vertical motion of waves or, after an approximation, an envelope of rapidly oscillating waves. This coordination of fluid movements with wave envelopes tends to occur when strong horizontal buoyancy gradients are present. The nonlinear models of this coordinated movement presented here may provide future opportunities for the optimal design of satellite imagery that could simultaneously capture the dynamics of both waves and currents directly.The models derived here appear first in their unapproximated form, then again with a slowly varying envelope (SVE) approximation using the WKB approach.The WKB wave-current-buoyancy interaction model derived here for a free surface with horizontal buoyancy gradients indicates that the mechanism for these correlations is the ponderomotive force of the slowly varying envelope of rapidly oscillating waves acting on the surface currents via the horizontal buoyancy gradient. In this model, the buoyancy gradient appears explicitly in the WKB wave momentum, which in turn generates density-weighted potential vorticity whenever the buoyancy gradient is not aligned with the wave-envelope gradient.