Sea Level Anomaly (SLA) and geostrophic currents in the subpolar Southern Ocean, including its sea-ice covered parts, from 2003 to 2021. The observational product relies on satellite altimetry, using specifically Envisat, Cryosat-2, Altika and Sentinel-3A missions. To retrieve physical signal in the ice-covered section of the Subpolar Southern Ocean, waveforms coming from leads are reprocessed and sea level measurements extracted from it. They are mapped along side sea level measurements from open ocean measurements to obtain a SLA field covering the whole Southern Ocean. SLA are referenced to the SLA field averaged over the whole period 2003-2021. A Mean Dynamic Topography from CNES/CLS is included to reconstruct the Absolute Dynamic Topography (ADT). ADT is the sum of SLA and MDT. Geostrophic currents derived from this ADT are included in the geostrophic_currents.nc file with u and v being respectively the eastward and northward velocities. The formal error of the mapping scheme is also included to give an estimate of the uncertainty. Dataset processing and description is available in Mosneron Dupin et al., Reconstructing Sea Level Anomalies in the open and ice-covered Southern Ocean from 2003 to 2021 [submitted].