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 pCO₂ and air-sea CO₂ flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO₂ 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 pCO₂ seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO₂ flux interannual variability (1σ) 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 CO₂ flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO₂ flux anomalies are correlated with the PDO, their magnitudes are small (up to ±0.025 PgC/yr (1σ)). Flux anomalies are damped because anomalies in the key drivers of pCO₂ (temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO₂ anomalies.

Using atmospheric forcing data generated from a general circulation climate model, sixteen land surface schemes participating in the Project for the Intercomparison of Land-surface Parametrization Schemes (PILPS) were run off-line to equilibrium using forcing data from a GCM representative of a tropical forest and a mid-latitude grassland grid point. The values for each land surface parameter (roughness length, minimum stomatal resistance, soil depth etc.) were provided. Results were quality controlled and analyzed, focusing on the scatter simulated amongst the models. There were large differences in how the models' partitioned available energy between sensible and latent heat. Annually averaged, simulations for the tropical forest ranged by 79 W m-2 for the sensible heat flux and 80 W m-2 for the latent heat flux. For the grassland, simulations ranged by 34 W m-2 for the sensible heat flux and 27 W m-2 for the latent heat flux. Similarly large differences were found for simulated runoff and soil moisture and at the monthly time scale. The models' simulation of annually averaged effective radiative temperature varied with a range, between all the models, of 1.4 K for tropical forest and 2.2 K for the grassland. The simulation of latent and sensible heat fluxes by a standard 'bucket' models was anomalous although this could be corrected by an additional resistance term. These results imply that the current land surface models do not agree on the land surface climate when the atmospheric forcing and surface parameters are prescribed. The nature of the experimental design, it being offline and with artificial forcing, generally precludes judgements concerning the relative quality of any specific model. Although these results were produced de-coupled from a host model, they do cast doubt on the reliability of land surface schemes. It is therefore a priority to resolve the disparity in the simulations, understand the reasons behind the scatter and to determine whether this lack of agreement in decoupled tests is reproduced in coupled experiments.