The coastal ecosystem of the Gulf of Alaska (GOA) is especially vulnerable to the effects of ocean acidification and climate change that can only be understood within the context of the natural variability of physical and chemical conditions.Controlled by its complex bathymetry, iron enriched freshwater discharge, and wind and solar radiation, the GOA is a highly dynamic system that exhibits large inorganic carbon variability from subseasonal to interannual timescales. This variability is poorly understood due to the lack of observations in this expansive and remote region. To improve our conceptual understanding 5 of the system, we developed a new model set-up for the GOA that couples the three-dimensional Regional Oceanic Model System (ROMS), the Carbon, Ocean Biogeochemistry and Lower Trophic (COBALT) ecosystem model, and a high resolution terrestrial hydrological model. Here, we evaluate the model on seasonal to interannual timescales using the best available inorganic carbon observations. The model was particularly successful in reproducing observed aragonite oversaturation and undersaturation of near-bottom water in May and September, respectively. The largest deficiency of the model is perhaps its 10 inability to adequately simulate spring time surface inorganic carbon chemistry, as it overestimates surface dissolved inorganic carbon, which translates into an underestimation of the surface aragonite saturation state at this time. We also use the model to describe the seasonal cycle and drivers of inorganic carbon parameters along the Seward Line transect in under-sampled months. As such, model output suggests that a majority of the near-bottom water along the Seward Line is seasonally undersaturated with regard to aragonite between June and January, as a result of upwelling and remineralization. Such an extensive 15 period of reoccurring aragonite undersaturation may be harmful to CO 2 sensitive organisms. Furthermore, the influence of freshwater not only decreases aragonite saturation state in coastal surface waters in summer and fall, but simultaneously also decreases surface pCO 2 , thereby decoupling the aragonite saturation state from pCO 2 . The full seasonal cycle and geographic extent of the GOA region is undersampled, and our model results give new and important insights for months of the year and areas that lack in situ inorganic carbon observations. Seasonal and spatial variability in northern Gulf of Alaska surface-water iron concentrations driven by shelf sediment resuspension, glacial meltwater, a Yakutat eddy, and dust, https://doi.