The upper ocean primary production measurements from the Hawaii Ocean Time series (HOT) at Station ALOHA in the North Pacific Subtropical Gyre showed substantial variability over the last two decades. The annual average primary production varied within a limited range over 1991-1998, significantly increased in 1999-2000 and then gradually decreased afterwards. This variability was investigated using a one-dimensional ecosystem model. The long-term HOT observations were used to constrain the model by prescribing physical forcings and lower boundary conditions and optimizing the model parameters against data using data assimilation. The model reproduced the general interannual pattern in the observed primary production, and mesoscale variability in vertical velocity was identified as a major contributing factor to the interannual variability in the simulation. Several strong upwelling events occurred in 1999, which brought up nitrate at rates several times higher than other years and elevated the model primary production. Our model results suggested a hypothesis for the observed interannual variability pattern of primary production at Station ALOHA: Part of the upwelled nitrate input in 1999 was converted to and accumulated as semilabile dissolved organic nitrogen (DON), and subsequent recycling of this semilabile DON supported enhanced primary productivity for the next several years as the semilabile DON perturbation was gradually removed via export.