1 Introduction
The Cariaco Basin is a small tectonic depression located on the northern continental shelf of Venezuela. Its depth extends to 1400 m, but the basin is largely separated from the Caribbean Sea by the Tortuga Bank and is only connected to the open sea by two shallow, ca. 130 m deep, sills (Fig. 1). Due to seasonal upwelling, the modern Cariaco Basin has high primary productivity (520 gCorg.m−2.a−1 on average for 1996–2002; Müller-Karger et al., 2001). The high oxygen demand induced by organic matter (OM) decay, added to weak renewal of bottom waters, led to progressive oxygen depletion of the water column (Peterson et al., 1991). Currently, the water column is euxinic below 300 m, and the sediments are relatively organic-rich, with TOC contents generally ranging between 2 and 5% (Daumas et al., 1978; Lyons et al., 2003; Richards, 1975).
First studied in the 1950s, the Cariaco Basin has seen renewed interest during the last two decades because it offers one of the best low latitude paleoclimatic records of the Pleistocene-Holocene (Haug et al., 2001; Hughen et al., 1996a; Peterson et al., 1991; Peterson and Haug, 2006). Its annually laminated sediments have also yielded several studies focused on better calibration of radiocarbon dates (Hughen et al., 1998, 2004b, 2006). The Cariaco Basin is the object of a long term oceanographic study, as it represents an ideal site to better constrain the links between primary productivity, nutrient cycles and flux of OM to the sediment (Müller-Karger et al., 2001, 2010).
The sediments of the Cariaco Basin display a temporal color pattern (gray scale) that closely parallels temperature records from Greenland ice cores (Hughen et al., 1996b; Peterson et al., 2000). These color patterns are related to variation in the total organic carbon (TOC) content of the sediment, the organic-rich intervals being the darkest (Peterson et al., 2000). Compared to the Holocene, Preboreal (PB), and Bølling-Allerød (BA) intervals, lighter sediment color and lower TOC content characterize the Younger Dryas (YD) interval (Hughen et al., 1996b; Piper and Dean, 2002; Werne et al., 2000). The origin of these low TOC values during the YD remains inadequately known. Sediments of the YD are rich in biogenic opal and calcium carbonate (Piper and Dean, 2002), contributing to their light color. The accumulation rate for terrigenous debris was also very high during this interval (Peterson et al., 1995), leading several authors to invoke dilution by biogenic and/or detrital minerals as a dominant explanation for low TOC and authigenic pyrite concentrations during the YD (Dahl et al., 2004; Lyons et al., 2003; Piper and Dean, 2002; Werne et al., 2003).
In the present article, we demonstrate by routine calculations that the low TOC values of the YD are not primarily a product of dilution by biogenic materials. Consequently, other factors influencing OM accumulation in the sediments of the Cariaco Basin through time and potentially accounting for these low TOC values are discussed.
2 Material and methods
This study is based on sediment from core MD03-2625 of the IMAGES Program (10°40.65′ N, 64°58.24′ W, water depth 847 m) recovered during the PICASSO cruise of R/V Marion Dufresne in 2003 (Laj and Shipboard Scientific Party, 2004) and Site 1002 of the Ocean Drilling Program (ODP) (10°42.73′ N, 65°10.18′ W; water depth, 893 m; Shipboard Scientific Party, 1997) drilled in 1996. Together, the two cores captured excellent records of the Holocene and Preboreal periods and the Younger Dryas and Bølling-Allerød intervals, respectively. The sediment from both cores was sampled at an average interval of ∼10 cm (∼200 y). Carbonate content was determined using a Bernard Calcimeter (via acid digestion). The relative error is about 2%. Total organic carbon was determined by Rock-Eval pyrolysis using a Delsi Oil Show Analyzer (OSA) at the Université P. & M. Curie (Paris 6). Specifically, pyrolysis cycle 3 was used, which is dedicated to the analysis of recent OM (Espitalié et al., 1986). Relative error for the TOC determination is 5%.
Major elements in bulk sediment were analyzed at the Service d’analyse des roches et des minéraux of the Centre national de la recherche scientifique (Vandœuvre-lès-Nancy, France) using ICP-AES after fusion with LiBO2 and HNO3 dissolution. Precision and accuracy are both better than 1%. The biogenic silica content (opal) was calculated as follows: opal = SiO2–Al2O3 × (SiO2/Al2O3)detrital, where (SiO2/Al2O3)detrital is the ratio of the detrital fraction of the sediment, assumed to be constant over the study interval. This assumption is reasonable given the very strong positive correlation between Si and Al in samples of Holocene age (r = 0.992, n = 33, not shown) where no diatom silica was reported (Peterson et al., 1995; Piper and Dean, 2002). Consistent with a high abundance of clay minerals in the sediments of the Cariaco Basin (Clayton et al., 1999), most of the studied samples, including some containing diatoms, have a SiO2/Al2O3 ratio lower than the average shale value (3; Wedepohl, 1991). (SiO2/Al2O3)detrital (i.e., 2.59) was therefore taken as the minimum SiO2/Al2O3 value of the studied interval. As shown on Fig. 2, the opal contents thus obtained are in good agreement with opal contents measured by Piper and Dean (2002) on nearby core PL07-39PC. The terrigenous content (Terr) in weight percentage is calculated as follows: Terr = 100–%opal–%CaCO3–%TOC.
Due to the different sedimentation rates in the two cores, the results are discussed in terms of sediment age. The age model used for core ODP 1002 has been described previously (Haug et al., 2001; Hughen et al., 2004a). For core MD03-2625, the age model was determined by correlating the colorimetric parameter a*, describing sediment color in the range from red to green (Laj and Shipboard Scientific Party, 2004), with gray scale data for core PL07-56PC for which a revised age model has been generated (Hughen et al., 2000). Correlation of colorimetric data has been shown to be a particularly efficient method for sediments of the Cariaco Basin (Hughen et al., 1996b). For the correlation, core PL07-56PC was preferred over ODP Site 1002 because it is located closer to core MD03-2625 and is very well dated (Hughen et al., 2000). The software used for the correlation was AnalySeries (Paillard et al., 1996). The ages are presented as calibrated years before present (a cal. B.P.). The sedimentation rates obtained by this method range from 0.3 to 1.4 mm.a−1, which are within the typical range for the Cariaco Basin (Hughen et al., 1996a; Peterson et al., 1991).
3 Results
Temporal variations of TOC, CaCO3 and opal contents are presented in Fig. 2. Though less pronounced in our dataset than in other studies (Mertens et al., 2009a; Piper and Dean, 2002; Werne et al., 2000), previously noted temporal patterns are clearly recognizable. In particular, we note the relatively low TOC values during the YD (mean 2.4%), when compared to the BA and PB with mean values of 4.1% and 3.9%, respectively. Maximum TOC values are observed at the end of the PB-beginning of the Holocene. The CaCO3 content also shows a characteristic behavior during the YD (Fig. 2). While the CaCO3 content is approximately 20% during the BA and PB, it rises to 30% during the YD. Opal is present in the sediments of the BA to PB and disappears after the PB. Conversely to TOC and CaCO3, the opal content is not constant during the YD. Relatively low at the beginning of the YD (7%), the opal content increases during the course of the YD to reach maximum values (30%) at the beginning of the PB. The TOC is also presented on a carbonate-free basis (TOCCF; Fig. 2). Similarly to TOC values, TOCCF values are lower during the YD interval compared to the BA and PB (Fig. 2).
In the studied interval, hydrogen index (HI) values show moderate variations between 250 and 450 mg hydrocarbons (HC)/g TOC. Despite slight discrepancies observed between the two studied cores, the values generally agree with those previously obtained by Aycard (2004) on Hole ODP 1002 C (Fig. 2). The BA and Holocene intervals show similar and relatively constant values around 340 mg HC/g TOC. The YD begins with lower values that increase during the course of the YD, to reach a maximum at the YD-PB transition. Values then progressively decrease to Holocene values during the course of the PB.
Fig. 2 also shows the same dataset from core PL07-39PC (Piper and Dean, 2002; data courtesy Larry Peterson). These results are very similar to those obtained for core MD03-2625. Slight differences are seen between core PL07-39PC and ODP site 1002, the latter showing lower biogenic mineral content during the BA (Fig. 2). The dissimilarity probably reflects the different location of the sites relative to the upwelling center (Fig. 1). Recent modeling (Lane-Serff and Pearce, 2009) showed that due to lower sea level, the Cariaco upwelling cell was smaller and shifted to the east during the BA, so that ODP Site 1002 was relatively distant from the upwelling center compared to the sites of core PL07-39PC and MD03-2625 (Fig. 1). It is worth noticing the very similar TOC patterns for both cores despite the relatively different biogenic mineral contents.
4 Discussion
4.1 Role of dilution
The TOC content shows no correlation with opal content, while a slight positive tendency (r = 0.53, n = 35) and a negative correlation (r = −0.76, n = 35) are observed with the terrigenous and CaCO3 contents, respectively (Fig. 3). From these observations, the carbonate content, that varies antithetically to TOC content through the BA-YD-PB intervals (Fig. 2), appears as a diluting phase for the OM during the YD. The study of Mertens et al. (2009a) allowed identify the origin of the generally higher carbonate content of the sediment during the YD compared to PB and BA: from the lower relative abundance of pteropods in the sediment (“pteropod ratio”), it appears that the Cariaco water column was more aggressive towards carbonates during the BA and PB than during the YD. The higher CaCO3 content of the sediment during the YD therefore reflects a lower dissolution of carbonates and more alkaline conditions in the water column. TOC values calculated on a carbonate-free basis (TOCCF) allowing one to overcome dissolution problems, however remain low during the YD compared to the PB and BA intervals (Fig. 2). From this observation, it appears that if carbonate dilution contributes to the low TOC values observed during the YD, it is not the major cause.
4.2 A reduced organic flux during the Younger Dryas
Mass accumulation rates (MAR) could not be determined for core MD03-2625, but a sufficiently robust age model allowed MAR calculations for core PL07-39PC (Piper and Dean, 2002, as updated by L. Peterson, Pers. Comm.). These data show that, during the YD, the flux of all sediment fractions increased, including the flux of OM (Fig. 5; Table 1). This observation led several authors to conclude that the YD was the period of highest primary productivity in the basin over the last 14.7 ka (Peterson et al., 1995; Piper and Dean, 2002; Werne et al., 2000). Careful examination of MAR values, however, reveals that the increase in MAR for TOC during the YD was muted compared to the other sediment components (Fig. 5; Table 1).
Taux d’accumulation massique moyens (MAR, en g.cm−2.ka−1) des différentes fractions du sédiment de Cariaco dans la carotte PL07-39PC (L. Peterson, 2009 Comm. Pers.).
CaCO3 MAR | Opal MAR | Corg MAR | Terr MAR | |
PB | 6.27 | 4.49 | 0.93 | 13.58 |
YD | 19.93 | 10.70 | 1.34 | 24.45 |
BA | 7.83 | 6.28 | 1.11 | 12.16 |
MAR ratios | ||||
YD/BA | 2.54 | 1.70 | 1.21 | 2.01 |
YD/PB | 3.18 | 2.38 | 1.44 | 1.80 |
Modern trap studies show that mineral ballasting has a strong effect on the OM flux to the sediment in the Cariaco Basin, such that the TOC and total mineral fluxes in the water column are strongly correlated (Thunell et al., 2007; Fig. 4A). Consistent with sediment trap data, a positive correlation is observed between MAR values for TOC and total mineral inputs in core PL07-39 PC (Fig. 4A). The slope of this correlation is, however, lower than for present-day sediment traps (Fig. 4A). This lower slope might be related to the location of core PL07-39PC on the margins of the upwelling cell, while the sediment traps were located in the center of the cell (Fig. 1). However, this lower slope must also point to OM degradation within the sediment after deposition (Goñi et al., 2003). Sediments deposited during the YD yield a third correlation line, with a slope that is lower than for the rest of core PL07-39PC (Fig. 4A). This lower slope points to a decreased Corg/mineral flux ratio during the YD compared to the rest of the interglacial. This reduced Corg/mineral flux ratio partly results from the better preservation of carbonates during the YD, so the same data excluding the CaCO3 MAR are shown on Fig. 4B. Though there is more scatter, Fig. 4B also shows the decreased Corg/mineral flux ratio during the YD compared to the rest of the interglacial. From these data, it appears that during the YD the flux of OM to the sediment, though high, was lower than expected.
4.3 Organic matter degradation
A first explanation to the reduction of the flux of OM relatively to the mineral flux during the YD is that OM was more efficiently degraded, either in the water column or within the sediment. In the sediments of the Cariaco Basin where OM is mainly of autochthonous origin (Combaz and Pelet, 1978; Dahl et al., 2004; Wakeham, 1990; Werne et al., 2000), HI values will preferentially reflect the degree of degradation of the OM (Espitalié et al., 1986). HI values observed in the studied interval are in the range of moderately well-preserved marine OM (Type II), consistent with the fact that despite anoxic conditions, OM is degraded in the water column (Thunell et al., 2000). Increased OM degradation during the YD should be indicated by a marked decrease of HI values. Compared to the BA, no particular decrease of HI values characterizes the YD, indicating a similar degree of preservation/degradation than during the BA. The slight increase in HI values in the course of the YD and maximum HI values at the beginning of the PB, however, could indicate better preservation of OM during this interval. Recent palynological data (Mertens et al., 2009b) also support this interpretation: the high abundance of degradation-sensitive dinoflagellate cysts in the sediments deposited during the YD and the beginning of the PB indicates a low degree of OM degradation during these intervals. Increased degradation of dinocycts is observed during the BA, PB and Middle Holocene (Mertens et al., 2009b). Both palynological and geochemical data indicate that the low TOC values observed during the YD are not related to OM degradation.
4.4 Source organisms and Corg/mineral flux ratios
In a basin where the terrigenous flux (including its terrigenous OM content) is considered as constant, changes of the Corg/mineral flux ratio will directly result from changes in surface producers: organic-walled organisms (i.e. those that do not have mineralized tissues, mainly corresponding to dinoflagellates, green algae, cyanobacteria and bacteria) mostly export OM from surface waters and therefore have high Corg/mineral export ratios; conversely, carbonate plankton and diatoms export a mineral test in addition to OM, and therefore have lower Corg/mineral export ratios. All other things being equal, a change in surface producers from dinoflagellate-dominated to diatom-dominated should lead to a decreased OM content in the sediment, a mechanism that Tyson (1995) termed “planktonic autodilution”. This simple scheme, however, can be complicated by selective degradations or dissolutions occurring in the water column. In the Cariaco Basin, diatoms are dominant during periods of upwelling, while cyanobacteria and dinoflagellates dominate during non-upwelling seasons (Ferraz-Reyes, 1983). In the sediment traps data from the CARIACO time series (Thunell et al., 2007), the samples with the highest proportion of OM, and therefore the highest Corg/mineral ratio, rather show a low proportion of opal, suggesting that this OM mostly derives from cyanobacteria and dinoflagellates.
From these observations, the decreased Corg/mineral flux ratio observed during the YD interval can be explained by a relative decrease of the contribution from organic-walled organisms compared to diatoms. The recent palynological study performed by Mertens et al. (2009b) gives support to this hypothesis. Despite very good preservation, the concentration and MAR of dinoflagellate cysts markedly decrease during the YD (Mertens et al., 2009b). Conversely, the PB is marked by an increase in cyst concentration, despite a lower preservation of very sensitive cysts (Mertens et al., 2009b). Keeping in mind the fact that cysts in the sediment represent a biased picture of the dinoflagellate biomass that existed in the water column (Head, 1996; Zonneveld et al., 2007), the palynological data tend to indicate a decreased dinoflagellates productivity during the YD followed by an important increase during the PB (Mertens et al., 2009b; also see Zonneveld et al., 2007). From a high-resolution biomarker study, Dahl et al. (2004) also concluded to a marked decrease of dinoflagellates compared to diatoms during the YD, followed by a marked increase at the transition to the PB. Therefore, the marked decrease in TOC values during the YD is to be related to the important decrease in dinoflagellates proportion. In addition, it can be observed that compared to the BA, the PB shows a high increase in dinoflagellate cyst abundance (Mertens et al., 2009b) and in dinoflagellate biomarkers (Dahl et al., 2004), as well as by TOC values that also are generally higher than the BA. From these observations, organic-walled organisms and in particular dinoflagellates appear as an important factor controlling the TOC content of the sediments of the Cariaco Basin.
4.5 Implications
The contribution from organic-walled organisms to total OM can be simply modelled using MAR data from core PL07-39PC and estimating the Corg/mineral ratio of the different sediment fractions (Fig. 5; Supplementary material). The Corg/mineral ratios listed in Table 2 were adjusted to get MAR values compatible with the Corg MAR data calculated by L. Peterson (Pers. Comm.). Comparison of sediment trap for carbonate plankton (Klaas and Archer, 2002) and diatoms (Thunell et al., 2007), showed that the Corg/mineral ratio of diatoms is approximately six times that of carbonate plankton. Variations in the values (Table 2, Supplementary material) take account of the higher dissolution of CaCO3 during the BA, PB and Late Holocene (Mertens et al., 2009a), and of the contribution of diatoms during the Holocene, despite quasi absence of opal in the sediment during this interval (Fig. 2). Though this model is only indicative, the results are consistent with biomarker data of Dahl et al. (2004) showing a change in OM sources in the middle of the YD (Fig. 5). If organic-walled organisms contribute only ca 20% of the OM during the YD, they could account for approximately half the OM in sediments of the BA and 2/3 for the PB (Fig. 5). During the Holocene, organic-walled organisms could represent up to 80% of the OM (Fig. 5).
Rapport Corg/minéral des différentes sources organiques utilisées pour la modélisation (Détail dans le e).
Carbonate plankton | Diatoms | Terrestrial flux | |
Corg/mineral | 0.010–0.015 | 0.060–0.200 | 0.005 |
In high productivity marine settings such as in the Cariaco Basin, Santa Barbara Basin or Saanish Inlet, diatoms flourish in surface waters and the sediment is enriched in OM and in diatom frustules (Johnson and Grimm, 2001). Diatoms have been shown to play an important role in the quantitative and qualitative export of organic carbon to the sediment (Goutx et al., 2007; Ragueneau et al., 2000). It is therefore tempting to consider that the OM deposited in high productivity environments mainly derives from diatoms. The present study of the sediments of the Cariaco Basin however shows that the contribution of OM derived from organic-walled organisms is far from negligible and could represent more than 50% of the OM, despite abundant presence of diatoms in surface waters. Such conclusion is consistent with the observation of a decoupling between surface primary productivity and organic carbon flux in the Cariaco Basin (Thunell et al., 2007). This decoupling demonstrates that export of OM is not exclusively associated to the blooming of diatoms and implies that other organisms contribute to the flux of OM to the sediment. This conclusion is also supported by numerous studies that have shown that marine snow often contains organic particles of different origins associated to opal (Alldredge et al., 1998; Goutx et al., 2007; Kovač et al., 2005). The close correlation observed between Corg and opal settling flux in the Cariaco Basin (Thunnel et al., 2007) is not incompatible with a multiple origin of the OM and could simply indicate that OM of diverse origin agglomerated with diatoms prior to settling in the water column. In the case of the Cariaco Basin, the contribution from organic-walled OM finally appears as the controlling factor of the TOC content of the sediment. While diatoms are important to quantitatively export organic carbon and have an important role in the carbon cycle (Ragueneau et al., 2000), organic-walled organisms appear determinant in the formation of (very) organic-rich sediments and therefore in the formation of future petroleum source rocks. The formation of organic-rich sediment is now recognized as resulting from a subtle balance between primary productivity, OM preservation and dilution by mineral phases (Tyson, 2005). The sediments of the Cariaco Basin, an apparently simple system with high surface productivity and an anoxic water column, illustrate the subtlety of this balance, emphasizing the importance of the nature of surface producers. Organic-walled phytoplankton was prominent during the Proterozoic and most of the Phanerozoic; the development of phytoplankton with mineralized tissues is a “modern” feature, in particular since the major rise of diatoms during the Tertiary (Katz et al., 2004; Knoll et al., 2007). This also explains the frequent observation that recent organic-rich sediments from high productivity settings are not as OM-rich as ancient black-shales, pointing that modern conditions are not very good analogues for ancient OM accumulations (Tyson, 1995, 2005).
Acknowledgements
We wish to thank the ODP program for providing sediment samples from core 1002 as well as the IMAGES program and crew of the R/V Marion Dufresne. Financial support for this study was provided by the project “Enfouissement organique dans le Bassin de Cariaco : contrôles climatique et diagénétique” through the French national EVE-CYBER-LEFE program (INSU-CNRS). Larry Peterson is greatly acknowledged for providing MAR data from core PL07-39PC. J.R. Disnar, P. Martinez and A.-Y. Huc are acknowledged for their constructive comments on the manuscript.