CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.

Optimizing Models of the North Atlantic Spring Bloom Using Physical, Chemical and Bio-Optical Observations from a Lagrangian Float

The North Atlantic spring bloom is one of the main events that lead to carbon export to the deep ocean and drive oceanic uptake of CO(2) from the atmosphere. Here we use a suite of physical, bio-optical and chemical measurements made during the 2008 spring bloom to optimize and compare three different models of biological carbon export. The observations are from a Lagrangian float that operated south of Iceland from early April to late June, and were calibrated with ship-based measurements. The simplest model is representative of typical NPZD models used for the North Atlantic, while the most complex model explicitly includes diatoms and the formation of fast sinking diatom aggregates and cysts under silicate limitation. We carried out a variational optimization and error analysis for the biological parameters of all three models, and compared their ability to replicate the observations. The observations were sufficient to constrain most phytoplankton-related model parameters to accuracies of better than 15 %. However, the lack of zooplankton observations leads to large uncertainties in model parameters for grazing. The simulated vertical carbon flux at 100 m depth is similar between models and agrees well with available observations, but at 600 m the simulated flux is larger by a factor of 2.5 to 4.5 for the model with diatom aggregation. While none of the models can be formally rejected based on their misfit with the available observations, the model that includes export by diatom aggregation has a statistically significant better fit to the observations and more accurately represents the mechanisms and timing of carbon export based on observations not included in the optimization. Thus models that accurately simulate the upper 100 m do not necessarily accurately simulate export to deeper depths.

Size-Fractionated Nitrogen Uptake Measurements in the Equatorial Pacific and Confirmation of the Low Si-High-Nitrate Low-Chlorophyll Condition

The equatorial Pacific Ocean is the largest natural source of CO(2) to the atmosphere, and it significantly impacts the global carbon cycle. Much of the large flux of upwelled CO(2) to the atmosphere is due to incomplete use of the available nitrate (NO(3)) and low net productivity. This high-nutrient low-chlorophyll (HNLC) condition of the equatorial upwelling zone (EUZ) has been interpreted from modeling efforts to be due to low levels of silicate ( Si( OH) 4) that limit the new production of diatoms. These ideas were incorporated into an ecosystem model, CoSINE. This model predicted production by the larger phytoplankton and the picoplankton and effects on air-sea CO(2) fluxes in the Pacific Ocean. However, there were no size-fractionated rates available for verification. Here we report the first size-fractionated new and regenerated production rates (obtained with (15)N - NO(3) and (15)N - NH(4) incubations) for the EUZ with the objective of validating the conceptual basis and functioning of the CoSINE model. Specifically, the larger phytoplankton ( with cell diameters > 5 mu m) had greater rates of new production and higher f-ratios (i.e., the proportion of NO(3) to the sum of NO(3) and NH(4) uptake) than the picoplankton that had high rates of NH(4) uptake and low f-ratios. The way that the larger primary producers are regulated in the EUZ is discussed using a continuous chemostat approach. This combines control of Si(OH)(4) production by supply rate (bottom-up) and control of growth rate ( or dilution) by grazing ( top-down control).

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.