Twentieth Century Delta δC Variability in Surface Water Dissolved Inorganic Carbon Recorded by Coralline Algae in the Northern North Pacific Ocean and the Bering Sea

The oxygen isotopic composition and Mg/Ca ratios in the skeletons of long-lived coralline algae record ambient seawater temperature over time. Similarly, the carbon isotopic composition in the skeletons record delta(13)C values of ambient seawater dissolved inorganic carbon. Here, we measured delta(13)C in the coralline alga Clathromorphum nereostratum to test the feasibility of reconstructing the intrusion of anthropogenic CO(2) into the northern North Pacific Ocean and Bering Sea. The delta(13)C was measured in the high Mgcalcite skeleton of three C. nereostratum specimens from two islands 500 km apart in the Aleutian archipelago. In the records spanning 1887 to 2003, the average decadal rate of decline in delta(13)C values increased from 0.03% yr(-1) in the 1960s to 0.095% yr(-1) in the 1990s, which was higher than expected due to solely the delta(13)C-Suess effect. Deeper water in this region exhibits higher concentrations of CO(2) and low delta(13)C values. Transport of deeper water into surface water (i.e., upwelling) increases when the Aleutian Low is intensified. We hypothesized that the acceleration of the delta(13)C decline may result from increased upwelling from the 1960s to 1990s, which in turn was driven by increased intensity of the Aleutian Low. Detrended delta(13)C records also varied on 4-7 year and bidecadal timescales supporting an atmospheric teleconnection of tropical climate patterns to the northern North Pacific Ocean and Bering Sea manifested as changes in upwelling.

North Pacific Carbon Cycle Response to Climate Variability on Seasonal to Decadal Timescales

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