Exploiting the complementary character of CryoSat-2 and Soil Moisture and Ocean Salinity satellite sea ice thickness products, daily Arctic sea ice thickness estimates from October 2010 to December 2016 are generated by an Arctic regional ice-ocean model with satellite thickness assimilated. The assimilation is performed by a Local Error Subspace Transform Kalman filter coded in the Parallel Data Assimilation Framework. The new estimates can be generally thought of as combined model and satellite thickness (CMST). It combines the skill of satellite thickness assimilation in the freezing season with the model skill in the melting season, when neither CryoSat-2 nor Soil Moisture and Ocean Salinity sea ice thickness is available. Comparisons with in situ observations from the Beaufort Gyre Exploration Project, Ice Mass Balance Buoys, and the NASA Operation IceBridge demonstrate that CMST reproduces most of the observed temporal and spatial variations. Results also show that CMST compares favorably to the Pan-Arctic Ice-Ocean Modeling and Assimilation System product and even appears to correct known thickness biases in the Pan-Arctic Ice-Ocean Modeling and Assimilation System. Due to imperfect parameterizations in the sea ice model and satellite thickness retrievals, CMST does not reproduce the heavily deformed and ridged sea ice along the northern coast of the Canadian Arctic Archipelago and Greenland. With the new Arctic sea ice thickness estimates sea ice volume changes in recent years can be further assessed.