Regional relative sea level rise is exacerbating flooding hazards in the coastal zone. In addition to changes in the ocean, vertical land motion (VLM) is a driver of spatial variation in sea level change that can either diminish or enhance flood risk. Here, we apply state-of-the-art interferometric synthetic aperture radar and global navigation satellite system time series analysis to estimate velocities and corresponding uncertainties at 30-m resolution in the New York City metropolitan area, revealing VLM with unprecedented detail. We find broad subsidence of 1.6 mm/year, consistent with glacial isostatic adjustment to the melting of the former ice sheets, and previously undocumented hot spots of both subsidence and uplift that can be physically explained in some locations. Our results inform ongoing efforts to adapt to sea level rise and reveal points of VLM that motivate both future scientific investigations into surface geology and assessments of engineering projects.