Evidence of stable foraminifera biomineralization during the last two climate cycles in the tropical Atlantic Ocean

: Planktonic foraminiferal biomineralization intensity, reflected by the weight of their shell calcite mass, affects global carbonate deposition and is known to follow climatic cycles by being increased during glacial stages and decreased during interglacial ones. Here we measure the dissolution state and the mass of the shells of the planktonic foraminifera species Globigerina bulloides from a Tropical Eastern North Atlantic site over the last two glacial-interglacial climatic transitions and we report no major changes in plankton calcite production with the atmospheric pCO2 variations. We attribute this consistency in foraminifera calcification to the climatic and hydrological stability of the tropical regions. However, we recorded increased shell masses midway through the penultimate deglaciation (Termination II). In order to elucidate the cause of the increased shell weights we performed δ18O, Mg/Ca and μCT measurements on the same shells from a number of samples surrounding this event. We find that shells of increased mass are internally contaminated by sediment infilling and that shell weights are respond to local hydrographic changes.


Introduction
Planktonic foraminifera are important marine calcifiers, and the ongoing change in the oceanic carbon system makes it essential to understand the influence of environmental factors on the biomineralization of their shells [1]. Although shell weight is a prominent, easily measured feature of foraminiferal tests that has direct implications for the carbon cycle and carbonate budgets, it is currently not being widely recorded or discussed because its exact environmental meaning is unclear. The amount of calcite deposited by planktonic foraminifera during calcification has been hypothesized to reflect a range of environmental factors. Through the years changes in planktonic foraminifera shell w eights have been linked to different biotic and abiotic parameters such as dissolution [2], carbonate ion concentration [3], optimum growth conditions [4], phosphate concentrations [5], temperature [6] or salinity [1] and thus the various studies have used foraminifera shell weights each time as a different proxy. Recent studies have shown that planktonic foraminifera can alter their shell mass according to ambient seawater density [7] and that the degree of this alteration in time is a function of latitude [8].
Planktonic organisms are able to biosynthesize out of equilibrium with their ambient environment by maintaining chemical gradients . However, as passive floaters they must always retain equilibrium with the seawater in order to remain afloat at certain depths. It can thus be argued that plankton physiology is more sensitive to the physical rather than the chemical characteristics of seawater. Since planktonic organisms lack active floatation devices their only inert way to counterbalance seawater buoyancy changes and remain at certain depths is to modify their shell mass. Shell biomineralization must thus be a function of both chemical and physical seawater properties. Although the effects of ocean chemistry on plankton have been extensively studied there is currently a lack of literature on the effects of physical oceanic properties such as buoyancy , density or pressure, which very likely affect foraminifera physiology and morphology [9]. Here we examine the shell mass of the planktonic foraminifera species Globigerina bulloides from a sediment core of the northeastern tropical Atlantic Ocean through the two most recent climatic cycles . After assessing the preservation of the foraminifera tests, we report consistent shell weights and thus steady foraminiferal calcification, independent of atmospheric pCO2, in agreement with a Pliocene Caribbean record [10]. We thus attribute this consistency to the stability of the hydrological conditions over time at the tropics, since the tropical environment is strongly associated with the notion of physical and chemical stability [11,12].
The stability in the foraminifera shell weights is briefly interrupted midway through Termination II, w here elevated weights are recorded. In order to understand this feature we used standard geochemical analyses (δ 18 O, Mg/Ca) for a set of samples that we combined with high resolution X-ray computed tomography (CT) to evaluate potential changes in the thickness of foraminifera shells [13] along with other biometric characteristics. The geochemical analyses confirmed the relationship between shell mass and water density, which is further supported by the μCT data that indicate clay contamination as the cause of the elevated shell weights. The μCT analysis allowed also the determination of cell volumes, volume normalized weights (i.e. shell density) or porosities and has proven a valuable tool in the study of foraminifera shells. The present record provides new evidence on the response of planktonic calcifiers to ocean acidification that will help to better constrain the role that shell mass variations ha ve on the sedimentary calcite budget and the carbon cycle.

Location and oceanographic setting
GeoB 8502 (19°13.27'N, 18°56.04'W) is a Tropical North Atlantic pelagic site at the lower reaches of the Cap Timiris Canyon approximately 250 km offshore the Mauritanian coast ( Figure 1). The canyon head abuts the Tamanrasset River System which, although not discha rging under presentday climate conditions, ranks among the largest river systems worldwide. The mouth of this potential river system is located off Cap Timiris but its flow pathways are at present covered by extensive Saharan sand dunes [14]. Presently, the Senegal River is the northernmost active drainage system of West Africa. Core GeoB 8502-2 was retrieved from 2956 m water depth on the lower Northwestern (NW) African continental rise and consists of levee sediments that are predominantly hemipelagic deposits.
Modern climate over the NW African margin is governed by the dynamics of the West African Monsoon, which is associated with the seasonal latitudinal shifts of the Intertropical Convergence Zone (ITCZ) [15]. In winter, the equatorward displacement of the ITCZ (5° N) causes a southward shift of dry subtropical air masses and is associated with the development of strong easterly Saharan Air Layer winds. The southward shift of the ITCZ and wind development causes dust transport from the Sahara [16]. The dust plume is generally located between 15 and 25° N along an E-W axis over the tropical Atlantic Ocean [17]. During boreal summer, dry subtropical air is shifted northward as the ITCZ is located around 20° N ( Figure 1). It represents the onset of the rainy season (summer monsoon) with heavy rainfall and changes in atmospheric circulation [18]. Trade winds from the southern hemisphere, loaded with water vapor, penetrate north to the West African continent. The moisture-laden air spreads over ocean and continent, permitting heavy rains over the area.
As part of the Eastern Boundary Current system, the Mauritanian upwelling region is one of the major upwelling areas in the Atlantic Ocean [19]. Along the NW African margin, the temporal dynamics of the coastal upwelling is driven basically by the intensity of the northeast trade-winds, itself dependent of the seasonal Intertropical Convergence Zone (ITCZ) migration [20,21] on a perennial basis, producing cold nutrient -rich surface waters with modern sea surface temperatures (SSTs) as low as 16C. The studied area is under the influence of the major return branch of the subtropical gyre, the Canary Current (CC), flowing southward along the north-west African coast, then becoming the North Equatorial Current (NEC) when turning southwestward and leaving the African continent. Further south the westward flowing North Equatorial Counter Current (NECC) is encountered, which transports low salinity water in the area and is known to show a strong seasonal cycle position with maximum velocities when the ITCZ is located at the northernmost position and weak velocities in northern spring [21].
The main water masses encountered in the upwelling region are the Tropical Surface Water (TSW), the North and South Atlantic Central Waters (NACW and SACW), and the Ant arctic Intermediate Water (AAIW) [21]. Due to the Cape Verde Frontal Zone (CVFZ), which is the transition boundary zone between NACW and SACW, the upwelling is fed by two different subsurface water masses depending on latitude [22]. The CVFZ is located at about 20N off Africa oriented southwestward to about 16N in the central tropical Atlantic. The front is associated with a convergence at the coast between the CC conveying NACW southward and a northward flow of SACW. With SACW and NACW occupying the same density range, the front is density-compensated and results in a multitude of intrusions, filaments and lenses [23]. The NACW is w armer and more saline compared to the SACW. Both central water masses appear in the permanent pycnocline between depths of 150 m and 600 m at temperatures greater than about 8 C below which lies the AAIW

Materials and methods
GeoB 8502-2 is 14.78m long with an average sedimentation rate of 6 cm/ky and extends back 200 ky to Marine Isotope Stage (MIS) 6. Samples were taken at a resolution of ~1700 years (10 cm sampling interval) by extracting a slice of material, 1cm in thickness, which corresponds approximately to an average of 170 years of sedimentation. All samples were freeze-dried, weighed, and washed through a 63-μm sieve. They were subsequently dry-sieved into several sieve fractions and non-fragmented G. bulloides shells from the 315-355 μm size fraction were picked for mass analysis. The very narrow size interval (40 μm) should be sufficient to overcome the greater proportion of natural size variability without further normalization [25]. Furthermore, this size fraction is widely used in paleoceanographic studies and can allow comparisons. However any record of shell weight is a composite signal of dissolution superimposed upon initial shell weight variability.
In order to assess whether the influence of dissolution is responsible for the observed shell weights the preservation state of the specimens underwent multiproxy assessment. Moreover to better understand the increase in G. bulloides shell mass during Termination II (T-II) selected specimens were analyzed both geochemically and by high resolution X-ray microcomputed tomography (XMCT). XMCT was used to inspect the interior and the internal structure of the foraminiferal tests, and apart from addressing test's integrity (manifestation of dissolution), we were also able to assess to what degree the increased recorded masses are the result of interference from shell inclusions or of changes in test thickness. Finally, the XMCT analysis lead to total shell volume estimates that allowed the calculation for the first time of volume norma lized shell weights or G. bulloides shell densities, presenting a more precise method of eliminating the contribution of shell size to shell weight.

Weight analysis
Where available ideally 50 (minimum 38) G. bulloides shells were weighed in a preweighed aluminum carrier using a Sartorius CP2P microbalance with a precision of ±1μg. Average shell masses were calculated by dividing the recorded mass by the total number of foraminifera weighed. Subsequently, for each sample each shell was weighed separately, in order to estimate standard shell mass deviations. As explained above, performing shell weight analyses on a narrow size fraction of foraminifera constrains the ontogenic stage of the specimens to a certain number of chambers and thus minimizes size related weight differences [e.g. 25,26,27]. The analytical error, estimated by triplicate measurements of 50 random specimens, ranged from 0.42 to 0.58 μg, which is in accordance with the error of the balance.

Dissolution assesment
Due to the immediate impact of dissolution on shell mass, a multiproxy approach was applied to determine the preservation state of the foraminiferal carbonate because indirect dissolution proxies may respond to influences other than dissolution alone. First, the conventional Fragmentation Index (F.I.) was counted in a sub-split of (~300) particles from the 315-355 μm sieve fraction. Subsequently, 30 of the weighed specimens were mounted for ultrastructural examination under the SEM using the G. Bulloides Dissolution Index (BDX′). As this species has a high dissolution susceptibility ranking [28], it is a reliable indicator of the foraminiferal carbonate preservation state. The BDX′ [29] is based on the direct assessment of corrosion of the specimen's surface and it consists of distinct dissolution stages that are determined by the decreasing preservation state of the four ultrastructural test features: pores, interpore space, spines, and ridges. As dissolution proceeds, pores are widened, the interpore areas are etched, and ridges and spines become denuded until the specimen is finally broken down. It thus consists of six stages worsening from 0 to 5 with severe dissolution starting above stage 3. The SEM analysis was performed with a Zeiss DSM 940A SEM in the Department of Geosciences at the University of Bremen.

X-ray microscopy (XRM)
For X-ray microscopic analysis, in total 27 specimens were scanned from 3 samples (9 specimens per sample) that correspond to the time interval of the increased mass event (132.2 ka), the penultimate glacial maximum (139.1 ka) prior to the event and a time interval of increased shell weights that followed (122.5 ka). Each batch of shells was poured into a quartz cylindrical carrier 1 mm in diameter. They were stabilized with diluted tragacanth glue and left to dry prior to scanning. The specimens were subsequently retrieved for the geochemical analyses. The micro-CT (μCT) scanning was carried out with a Zeiss Versa 520 at the X-ray Microscopy laboratory of London Natural History Museum. X-ray source and detector geometry were kept constant throughout the scans. Anode voltage was set at 100 kV, the X-ray tube current at 90 μΑ and the exposure time to 2 seconds at an optical magnification of 4X. By processing approximately 1024 images per sample a scan resolution voxel size of ∼1.2 µm 3 was typically achieved using this set up in order to maximize the number of specimens that could be analysed in a single scan. The images were combined to build a 3D rendering using Avizo 9.2 software which was also used for segmentation. The segmentation resulted in the separation of the tomographs into shell area, area occupied by clay infillings (dirt) and internal shell (protoplasm) voids ( Figure 2).

Figure 2:
Example of tomograph se gmentation using CT data visualization software , whe re (a) in yellow the foraminife ra te st is solely segmented, (b) in red the are a covered by dirt was also segmented, (c) dirt and inte rnal (chamber) voids we re considered together to calculate the area that was occupied by protoplasm (in blue) and (d) protoplasm and shell areas we re merged together (in cyan) to calculate total volume that a foraminife r occupies.
Subject to the degree of segmentation the X-ray microscopic analysis allows the determination and study of a number of biometric characteristics of the foraminifera shells, such as total shell volume ( Figure 2d) and thus shell density (volume normalized weight), calcite (test) volume and thus test density [30], calcite (test) surface (Figure 2a) etc. The ratio of calcite volume/calcite surface provides a linear unidimensional quantity in length units and can thus serve as a measure of average test thickness. In this study, in addition to shell density, i.e. the ratio of shell volume to shell mass, we use the "specific surface area", i.e. the ratio of test volume/test surface, as a measure of average test thickness [30] and the test density, i.e. the ratio of test volume to shell mass as indication of test porosity. Furthermore, by segmenting the area occupied by clay infillings we were able to calculate by volume percent the degree of contamination in weight measurements.

Geochemical analysis
In order to reconstruct the oceanographic regime of the time interval surrounding the increased shell mass event during T-II, 12 adjacent down core samples, covering the period of ~22,000 yrs between 119 and 140 ka, were geochemically analysed to determine shell Mg/Ca ratios for seawater temperature estimations and δ 18 Ο isotope (δ 18 Οshell). Mg-derived temperatures and foraminiferal δ 18 Οshell were used together to assess changes in seawater δ 18 O (δ 18 Osw) composition. Mg/Ca ratios were converted to temperatures using an equation for temperate ranges [31]. δ 18 Osw was calculated using a published paleotemperature equation [32], and a VPDB-to-SMOW δ 18 Ο conversion of 0.27 ‰ [33]. Salinity estimates were derived from δ 18 Osw using the Tropical Atlantic Salinity-δ 18 Osw relationship [34] and corrected for deglacial whole ocean salinity changes using a 120 m scaled sealevel curve [35] and an average ocean depth of 3800 m. Seawater densities were calculated from temperature and salinity estimates using the equation state of the water [36] at a (G. bulloides calcification) depth of 100 m. In order to test the correlation between shell mass and the different environmental and oceanographic data two-tailed regression analyses were performed using the Reduced Major Axis model, at n-1 degrees of freedom (n representing sample size).

Oxygen stable isotope determination
For oxygen isotopes determination, on average 5 G. bulloides tests from each sample were transferred in-to sample vials, crushed and then dried in an oven at 50 C. The vials were loaded into the carousel and analysed using a Thermo Kiel device attached to a Thermo MAT253 Mass Spectrometer in dual inlet mode at The Godwin Laboratory for Palaeoclimate Research, Department of Earth Sciences, University of Cambridge. The preparation system operates automatically analysing samples in sequence. 100% orthophosphoric acid is dropped onto the evacuated vial and reacts with the calcium carbonate sample. The evolved carbon dioxide is cryogenically dried and then admitted to the dual inlet mass spectrometer for isotopic analysis by comparison with a reference gas. Each run of 30 samples was accompanied by 10 reference carbonates and 2 control samples. The results are reported with reference to the international standard VPDB and the precision is b etter than +/-0.06 per mil for δ 13 C and +/-0.08 per mil for δ 18 Ο.

Mg/Ca determination
On average 15 G. bulloides tests from each sample were used for Mg/Ca determination. The specimens were cleaned using the standard protocol [37] omitting the reductive treatment and were analysed for Mg/Ca ratios by inductively coupled plasma optical emission spectrometry (ICP-OES) using an Agilent 5100 instrument at the Godwin Laboratory for Palaeoclimate Research, Department of Earth Sciences, University of Cambridge. Samples were dissolved in 0.1 M HNO3 and centrifuged to remove any undissolved material. After an initial run to determine Ca concentration, samples were diluted to constant [Ca] (100 ppm). Mg/Ca ratios were determined by the intensity ratio method of [38] using calibration standards prepared according to [39]. Cleaning efficiency and diagenetic effects were monitored by measuring Fe/Ca, Mn/Ca, Al/Ca, Si/Ca and Ba/Ca ratios.
Instrumental precision for the Mg/Ca ratios is ±0.51%, determined by replicate analyses of a standard solution containing Mg/Ca ratios of 1.3 mmol mol −1 and a Ca concentration of 100 ppm. Accuracy of Mg/Ca ratios has been established by interlaboratory calibration [40]. Sample heterogeneity for Mg/Ca ratios has been shown to be much greater than instrumental precision [37,41]. We estimate the reproducibility of planktonic foraminiferal Mg/Ca ratios is ∼8% from replicate analyses of G. bulloides picked from an Atlantic core-top sample.

Results
The record of G. bulloides shell mass attained from GeoB 8502-2 for the last 200 ky shows no distinct mode of variability (Figure 3a). In general, shell weights are stable, fluctuating only on a small scale (±1.1 μg, 1σ) around an average of 13.4 μg. They thus do not exhibit any glacial-interglacial cyclicity or follow the atmospheric pCO2 fluctuations (Figure 3d). Lower shell mass weights are found within MIS 6, while values almost increase progressively after MIS 5.2. Superimposed on this pattern is a broad maximum in shell weight centred on the MIS 5/6 boundary of T-II with a duration of approximately 2,300 yr, during which shell masses increase by 30% above average. The multiproxy carbonate preservation assessment of the studied core is shown in Figure 3b & 2c. The record of fragmentation strikingly parallels that of the atmospheric pCO2, indicating influence of carbonate chemistry changes on the study samples. Nonetheless, this influence does not prove sufficient to severely alter foraminifera shell mass. In general, foraminiferal carbonate is found to be well preserved (Figure 3b & 2c) since any severe dissolution takes place after BDX' values of 3 [29]. The BDX' and the F.I. assessments (Figure 3b and 2c respectively) are in agreement indicating that the preservation becomes better during the course of the last glacial period. However, regardless of the increase in foraminifera fragmentation after the Last Glacial Maximum, the G. bulloides specimens' ultrastructure was not found to be corroded. A similar disagreement between the dissolution proxies exists for MIS 6, where although the surface ultrastructure of G. bulloides specimens shows signs of corrosion, overall foraminiferal breakage is found to be limited. The deteriorated shell preservation state during MIS 6 as suggested by the BDX' might be adequate to explain the slightly decreased G. bulloides shell weight values of this period, but on the other hand the gradual decrease of the BDX' values towards the present may explain a gradual increase in mass. inte rglacial pe riods. Data are available in [43].
The results of the geochemical analyses of the samples surrounding T-II are summarized in Figure 4 together with a focused record of G. bulloides shell weights, atmospheric pCO2 and relative sea level. δ 18 O range from -0.09 to 0.71‰ and w ith the exception of the first sample in MIS 6, in general they are found to be more depleted within the glacial period. The most depleted value is found near MIS 6.2 (133.9 ka) after which valuesstart to increase and become heavier during the Eemian. Mg/Ca ratios range from 2.33 to 4.15 mmol/mol but are in general consistent around 2.8 mmol/mol during this period. However, they are slightly more depleted prior to the termination and exhibiting a spike with a ~46% increase in their values during the MIS 6.1 interstadial. This increase should indicate environmental signal as all the contamination indicators (Fe, Al etc) were l ow for this sample.
The results of the XRM analysis of the samples covering the 22 ky interval surrounding the 30% increase in shell weight event during T-II are summarized in Table 1. Shells are on average ~16% thicker during the time that where they are heavier (132.2 ka), with a mean thickness of 5.7 μm. They are also ~32% more voluminous/larger within the same sieve fraction during the same time but with the highest size variability. Apart from being larger and thicker they are also found to be more than three times contaminated by clay infillings (19% of their interior void) and thus the observed increased in weight is a result of both enhanced biomineralization and also contamination. Shell densities or volume normalized weights do not change significantly (up to 7% from the densest during the time of the highest masses to the lightest during t he interglacial) but values show high variabilitye. In terms of test density or porosity the heaviest and densest shells are only intermediate porous. Sample 865 at 139.1 ka should record the Penultimate Glacial Maximum (PGM) in the present archive, which in the northern hemisphere is placed at ~140 ka BP [44], since after that the sea level and atmospheric CO2 rise continuously (Figure 4). The low shell weights recorded during this period indicate reduced calcification that is further supported by the μCT scans that reveal small, thin and highly porous shells of low density. The scanned shells closest to the peak of the Eemian interglacial period were found to be of average mass, intermediate size, thickness and density but were also the least porous. The interior of the Eemian and PGM specimens were found to be filled with sediment only by 5% and 4% respectively w ith thus with minimum distortions to the weight measurements.  the me an shell mass of this pe riod and the shading the 1σ confide nce inte rval. Enlarged circles denote the sample s that have been X-ray microscopically analysed. The thickne ss of the ir e nvelope line is proportional to the ave rage shell wall thickne ss shown in Table 1; (b) δ 18 Oshell of the foraminiferal calcite ; (c) Mg/Ca ratios of the foraminife ral calcite; (d) re lative sea le vel re cord [45] is shown for comparison and (e ) the re cord of atmospheric CO2 from Vostok [42] is shown also for comparison.
The 1σ confide nce inte rval for the she ll mass plot is shown as transparent polygons. Numbe rs refer to Marine Isotopic Stages and substages. Geochemical and shell mass data are available in [43].
The geochemical analyses allow the reconstruction of the physical oceanographic parameters that characterize the period surrounding T-II at the study area and are shown in Figure 5a yrs between 139.1 and 135.7 ka and then a ~24% increase to ~38 psu after the PGM before decreasing back to 34 at the peak interglacial. The same applies with the reconstructed densities of the water column that reveal lighter waters during the glacial and denser waters after MIS 6.1. Following the previous study [7], shell weights were plotted against the reconstructed seawater density values and a significant (R 2 =0.43, p<0.05, n=12) relationship between these parameters was verified (Figure 5d), which becomes even more significant (R 2 =0.74, p<0.01, n=12) if the (upper right) 2 heaviest (most contaminated by clay infillings) samples are omitted. It is evident that these two samples of excessive shell weights deviate from the resulting regression line. If this regression line dictates shell mass according to shell's geochemical signal, then we can infer that these two samples are contaminated with 2-3 μg of sediment. The two uppe r right points are the increased shell mass sample s that de viate from the regression line by 2 and 3 μg re spectively. Enlarged circles de note the samples that have be en X-ray microscopically analysed. The thickne ss of the ir e nve lope line is proportional to the ave rage she ll wall thickness shown in Table 1. Data are available in [43] and XRF data from [46].

Discussion
The sieved based weight analysis of G. bulloides shells from the tropical Atlantic core GeoB8502-2 revealed a typical value of 13.4 μg and only small variability in the average mass throughout the last two climatic cycles, which does not follow the atmospheric pCO2 record. This stability in planktonic foraminifera shell masses despite changes in atmospheric pCO2 has also been reported in similar latitudes for the Pliocene [47] and may be the result of the environmental stability of the tropical regions. The tropics are generally an environment in which the physicochemical factors are not undergoing major changes through time [11] and this oceanographic stability can explain the observed stability in the foraminifera shell weight record that we report here. The incoming solar energy varies considerably from tropical to polar la titudes. At the equatorial region the average Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 10 September 2020 doi:10.20944/preprints202009.0222.v1 insolation is the highest with only minor peak energy changes, while at middle and high latitudes the regular periodicities in the Earth's orbit and tilt influence the amount and distribution of incoming energy and temporal changes are greater [48].. Test weights and range of test weight in North Atlantic Pliocene and Pleistocene G. bulloides are similar despite the higher pCO2 during the Pliocene [10]. However, the absolute G. bulloides shell mass values for specimens 300-350 μm in size reported here for a tropical Atlantic site are low compared to the North Atlantic records [47,49]. This is consistent with the influence of seawater densities on shell weights [7], as the lowest sea surface density values are found in the tropical waters, where salinity is lowest due to excess rainfall associated with the Intertropical Convergence Zone, and increases toward the poles [50]. At higher latitudes G. bulloides shell weight changes by around onethird over glacial-interglacial cycles [3,49,51] and similar variability in test weight is seen through the Pliocene glacial-interglacial cycles at a similar location [47]. At these higher latitudes that are more sensitive to (latitudinal) changes in solar radiation through time, foraminifera shell weights follow the pronounced hydrological changes between the climatic cycles, with the most sensitive latitudes for explaining the glacial/interglacial ice-dynamics around 65[52], leaving the tropical shell mass records almost invariable. The consistency in foraminifera shell mass may reflect hydrological stability in the study area and thus stability of the ITCZ locations during the last 200 ky.
Both the Fragmentation Index and Bulloides Dissolution Index preservation proxies agree that the biogenic carbonate preservation is in general good to very good. There is a clear match in the trends of both proxies for deterioration in the preservation of foraminifera specimens from MIS 2 to MIS 5.5. Increased specimen dissolution is also evident in the tomographs of the T-II sample [43]. This deteriorating trend however is not sufficient to cause obvious reductions in the weights of G. bulloides shells during this time interval. The two proxies disagree within MIS 6, where the degree of ultrastructural corrosion does not parallel the degree of specimens' fragmentation. This dis crepancy points either to the application limits of the proxies or is the result of different test architecture. The F.I. is an indirect measure of shell dissolution and may be the result of ecologic factors like shell initial thickness or different environmental ones (e.g. degree of bioturbation) but in this case it closely follows the atmospheric pCO2, implying influence of carbonate chemistry. On the other hand, the BDX' assesses specimen dissolution directly but it is more subjective than F.I. as a proxy. However, if specimens are smaller during the late MIS 6, as suggested by the CT analysis (Sample 865; Table 1), then they would have been less prone to fragmentation even at a higher degree of corrosion and this could explain the offset in the indices.
Carbonate dissolution in the equatorial Atlantic was traditionally believed to be greater during interglacial periods mainly due to lower carbonate contents in the glacial sequences that fluctuated in response to Quaternary climatic oscillations [53][54][55][56]. However, Broecker, et al. [57] showed, in a core from the north Atlantic, that net carbonate input actually increased in the last glacial but a proportionately higher detrital terrigenous influx diluted the carbonate in accordance with later studies [58][59][60] that showed that large glacial-interglacial differences in carbonate concentration are caused mainly by increased dilution by non-carbonate material during low sea levels. Based on the multiproxy dissolution assessment and the consistency of foraminifera shell weight, the present analysis shows that biogenic carbonate dissolution has overall been minimal in the eastern tropical north Atlantic, with a slight increase towards the peak interglacial. Slightly increased dissolution of fossil shells is to be expected during interglacials due to the bi omineralization of dissolution susceptible higher-Mg calcite. Together with findings from the western equatorial Atlantic [60], the present study indicates that the average carbonate productivity in the tropical Atlantic has been constant during the last two climatic cycles. Furthermore, the strong carbonate preservation in the study site suggests that the Cape Verde plateau has remained under the influence of the noncorrosive NADW during the past 200 ky, the carbonate chemistry of which closely followed the the atmospheric pCO2.
The homogeneity of the G. bulloides shell mass record is interrupted by a spike in the mass values during T-II, where weights increase by 30% above average. In order to explain this abnormality a series of geochemical and tomographic analyses were performed that point to sediment Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 10 September 2020 doi:10.20944/preprints202009.0222.v1 contamination as the cause of the increase of the measured weights. The first elevated masses are recorded at 132.2 ka during MIS 6.1, where high Mg/Ca temperatures are also recorded, and persisted for approx. 2,400 years. With the exception of MIS 6.1, when 20.3C was reconstructed, temperatures in the area are relatively stable around 16-17C and only slightly reduced during glacial times. For comparison, in the upwelling region of the area the full range of glacial to interglacial temperature changes amounted to more than 6C [61] and such variations in SST most probably result from changes in the intensity of coastal upwelling, and in turn the consequence of fluctuations in tradewind strength [62]. Peak warm stages, accordingly, correspond to a cessation of upwelling and trade winds. Indeed, temperatures of 25C are characteristic of the subtropical Atlantic outs ide the upwelling region [63]. T-II is defined by the fairly gradual δ 18 Οshell decrease between substages 6.2 and 5.5 as it has been manifested for the study area before [64]. Coupled δ 18 Οshell -Mg/Ca derived salinities and densities reveal that since temperatures are relatively stable during this time interval, it is salinity that dictates sea surface density (SSD) and that major hydrological changes in the area take place after MIS 6.1, when salinities return to normal marine conditions. During the late MIS 6 the geochemical data suggest salinities as low as 29 for ~3,500 years. The presence of a less-saline water surface layer will most often lead to a reduced production, as higher density differences will slow down the upwelling system. T-II contains two major meltwater pulses [65] (MWP-2A and MWP-2B) that on one hand can explain low salinities, which are centered on 139 ± 1 and 133 ± 1 ka [66] and coincided within uncertainties with two North Atlantic cooling episodes [67]. MWP-2A indicates an early phase of ice-sheet retreat. MWP-2B is more convincingly resolved and marks a steep ∼70 m sea-level rise (∼70% of the glacial-interglacial change) at rates of 28 ± 8 m ky −1 [65,66]. Although small scale perturbations in SSSs and SSDs such as during MIS 5.5 may have been explained by changes in the position of CVFZ and accordingly shifts in the influence of NACW and SACW, salinities in the order of 29 indicate influence of continental discharge, since similar values are found at sites closest to the river mouth in western Africa [68,69]. Fluvial delivery is also suggested by the increase in Ba 2+ core X-ray fluorescence (XRF) data (Figure 5b), which is used to infer freshwater inputs [69,70]. Seawater Ba concentrations at oceanic sites influenced by riverine runoff have a notably high inverse correlation to salinity , because dissolved Ba is high in riverine water [71]. Saharan aquifers were found recharged during glacials [72] and such humid periods trigger the reactivation of the Tamanrasset river system in Western Sahara [14] part of which is the Cap Timiris canyon.
Shell weights closely follow SSD changes (R 2 =0.43, p < 0.05, n=12) during the time interval analyzed (Figure 5d). The two samples that exhibit abnormally elevated shell masses within the record are also shown to be outliers in the shell mass -water density plot and according to the CT analysis performed in one them the specimens were 4 times more contaminated by clay impurities than the rest of the scanned samples. If these two contaminated samples are omitted the correlation between weighed shell mass and geochemically reconstructed seawater densities becomes more significant (R 2 =0.74, p < 0.01, n=10), while there is no significant correlation between shell weights and atmospheric pCO2 (R 2 =0.16 p < 0.10, n=12) during this period. By using the shell mass -density regression line to correct the elevated weights, a contamination of 2-3 μg is found. If these corrections are applied then the slope of the regression line (0.38) matches that of Zarkogiannis, et al. [7] (0.39) that appears to be characteristic for G. bulloides. Based on this relationship the low shell weights recorded prior to MIS 6.1 (Figure 4a) can be explained by the freshening of Mauritanian waters. Furthermore, during this time interval the XRM analysis showed that foraminifera shells have the lowest volume normalized weights and they are smaller, thinner and of highest porosity ( Table 1). The lowest porosities are found close to the peak interglacial, when all other shell cha racteristics are intermediate. The heaviest, thickest and densest shells are indeed recorded during the increased shell mass event (132.2 ka) when specimens are also largest but of intermediate porosity (test density). Increased shell wall thicknesses and shell and test densities are characteristics of increased biomineralization efforts and play a partial role in the observed increased shell masses.

Conclusions
Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 10 September 2020 doi:10.20944/preprints202009.0222.v1 Based on a record of planktonic foraminifera shell weights from a sediment core in the eastern tropical north Atlantic, we provide evidence of steady carbonate productivity in the equatorial Atlantic Ocean during the last two climatic cycles. A multiproxy dissolution assessmentfound foraminiferal carbonate to be well preserved, and the consistency in the measured shell masses is a reliable indication of stable foraminifera biomineralization strength independent of changes in atmospheric pCO2. We attribute this invariability in the mass of G. bulloides shells to the stability of the local hydrologic conditions, which also reveals a stability in the position of the ITCZ during the last two glacial-interglacial stages. Furthermore, the enhanced preservation of the foraminiferal record indicates that the Cape Verde plateau remained under the influence of the non-corrosive NADW during the study interval.
The geochemical analysis of the samples surrounding the termination revealed a freshening of the surface waters prior to PGM, which coincided with fresh water pulses or the reactivation of a sub-Saharan drainage system. Furthermore, the reported relationship between G. bulloides shell weight and ambient seawater densities was verified, reinforcing the use of foraminifera shell weight as a hydrological indicator. The X-ray tomographic analysis of selected samples allowed the determination of volume normalized weights and revealed that the spike observed in the G. bulloides shell weight record during T-II is an artifact of specimen contamination by clay infillings. Overall the tomographic analysis proved to be a very powerful micropaleontological tool for the determination of key shell characteristics (e.g. cell volume, shell thickness, porosity etc.), which are important for paleoclimatic studies.