A 10-year reanalysis of the PacIOOS Hawaiian Island Ocean Forecast System was produced using an incremental strong constraint 4D-Variational data assimilation with the Regional Ocean Modeling System (ROMS). Observations were assimilated from a range of sources: satellite-derived sea surface temperature (SST), salinity (SSS), and height anomalies (SSHA); depth profiles of temperature and salinity from Argo floats, autonomous SeaGliders, shipboard conductivity-temperature-depth (CTDs); and surface HFR velocity measurements from high frequency radar (HFR). The performance of the state-estimate is examined against a free-running forecast showing an improved representation of the observations, especially the realization of HFR surface currents. EOFs of the increments made during the assimilation to the initial conditions and atmospheric forcing components are computed, revealing the variables that are influential in producing the state-estimate solution and the spatial structure the increments form.