Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 1611–1626, doi:10.1175/JPO-D-12-0204.1. ; A new method is proposed for extrapolating subsurface velocity and density fields from sea surface density and sea surface height (SSH). In this, the surface density is linked to the subsurface fields via the surface quasigeostrophic (SQG) formalism, as proposed in several recent papers. The subsurface field is augmented by the addition of the barotropic and first baroclinic modes, whose amplitudes are determined by matching to the sea surface height (pressure), after subtracting the SQG contribution. An additional constraint is that the bottom pressure anomaly vanishes. The method is tested for three regions in the North Atlantic using data from a high-resolution numerical simulation. The decomposition yields strikingly realistic subsurface fields. It is particularly successful in energetic regions like the Gulf Stream extension and at high latitudes where the mixed layer is deep, but it also works in less energetic eastern subtropics. The demonstration highlights the possibility of reconstructing three-dimensional oceanic flows using a combination of satellite fields, for example, sea surface temperature (SST) and SSH, and sparse (or climatological) estimates of the regional depth-resolved density. The method could be further elaborated to integrate additional subsurface information, such as mooring measurements. ; JW and AM were supported by NASA (NNX12AD47G) and NSF (OCE 0928617). JLM was supported by the Office of Naval Research and the Office of Science (BER), U.S. Department of Energy under DE-GF0205ER64119. GRF is supported by OCE-0752346 and JHL by NORSEE (Nordic Seas Eddy Exchanges) funded by the Norwegian Research Council. ; 2014-02-01