Climate variability drives significant changes in the physical state of the North Pacific, and there may be important impacts of this variability on the upper ocean carbon balance across the basin. We address this issue by considering the response of seven biogeochemical ocean models to climate variability in the North Pacific. The models' upper ocean pCO2 and air-sea CO2 flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO2 and its temperature- and non-temperature-driven components at three contrasting oceanographic sites capture the basic features found in observations (Takahashi et al., 2002, 2006; Keeling et al., 2004; Brix et al., 2004). However, particularly in the Western Subarctic Gyre, the models have difficulty representing the temporal structure of the total pCO2 seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO2 flux interannual variability (1s) in the North Pacific is smaller (ranges across models from 0.03 to 0.11 PgC/yr) than in the Tropical Pacific (ranges across models from 0.08 to 0.19 PgC/yr), and the time series of the first or second EOF of the air-sea CO2 flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO2 flux anomalies are correlated with the PDO, their magnitudes are small (up to ±0.025 PgC/yr (1s)). Flux anomalies are damped because anomalies in the key drivers of pCO2 (temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO2 anomalies.