The usefulness of GRACE and GRACE Follow-On (GRACE/FO) data products for ocean studies remains limited in a few areas, where large earthquake-driven mass signals from the 2011 Japan earthquake and three major earthquakes in the Andaman-Sumatra area (in 2004, 2007, and 2012) obscure the vastly smaller signals of interest related to ocean dynamics. Furthermore, these signals can not be simply removed from the the GRACE/FO products via a model during post-processing, due to the large noise they introduce in the estimation of the GRACE/FO products. Here, we create a multi-earthquake model from GRACE/GRACE-Follow On data, estimating the mass changes as seen in the Center for Space Research (CSR) mascons due to the four major earthquakes mentioned. An iterative principal component analysis was used to compute piecewise linear approximations of the signals in these regions, which describe the co-seismic bias changes during the known month of each earthquake as well as post-seismic trends in local mass after/between quakes. This model is then applied as part of the estimation background field during a new iteration of CSR mascon processing designed for oceanic use, using tightened regularization to reduce noise in the near-earthquake regions. By removing the CSR-derived earthquake model and reducing month-to-month variability via the new regularization, this the resulting mascon series is more suitable for oceanic studies near Japan and Andaman Bay.