We have used the latter almost exclusively in our own work, in part because the latter approach
is embedded in a modeling environment (ANUDEM or the Topo to Raster tool in ArcGIS 10.x) that tackles the last three problems unwanted depressions, unresolved flow directions in flat terrain, and reconciliation of elevation and hydrography datasets) simultaneously.