A new algorithm for estimating sea ice age (SIA) distribution based on the Eulerian advection scheme is presented. The advection scheme accounts for the observed divergence/convergence and freezing/melting of sea ice and predicts consequent generation/loss of new ice. The algorithm uses daily gridded sea ice drift and sea ice concentration products from the Ocean and Sea Ice Satellite Application Facility and produces gridded fractions of sea ice of a specific age. Thus, each grid cell of the output product contains a frequency distribution of sea ice age allowing to apply mean, median, weighted average or other statistical measures. The produced SIA maps and time series are compared to the National Snow and Ice Data Center SIA product. Several improvements related to the usage of the new ice drift product, constraining the algorithm by the observed ice concentration and preventing undersampling by the Eulerian scheme are demonstrated. Muliyear ice (MYI) concentration is computed as a sum of all multi-year ice fractions and compared to the MYI products based on passive and active microwave and SAR products.