18 surface sediment samples collected from a north-south transect along the Indian Ocean have been analyzed for planktonic Foraminifera content. Among the other planktonic foraminiferal faunas, Globigerina bulloides was present substantially in all samples. Census data of G. bulloides were measured for different parameters (average size, mean proloculus size, coiling direction, and number of chambers) and a Q-mode cluster analysis was applied on these data. Samples were segregated into two homogeneous clusters, each reflecting particular environmental conditions. Two clusters are as follows: (1) Cluster A, comprised of 6 samples and characterized by the highest range of foraminiferal and ecological parameters, except sea surface temperature and salinity which shows the lowest range, and (2) Cluster B, comprised of 12 samples and characterized by the lowest range of foraminiferal parameters and ecological parameters, except sea surface temperature and salinity which shows the highest range. The study suggests that the ecological parameters are the governing factors for the morphological characteristics of planktonic foraminiferal species G. bulloides.
1. Introduction
Morphological variations of planktic Foraminifera have been extensively used to decipher paleoclimatic, paleoenvironmental, and paleoecological reconstructions [1–3].
Globigerina bulloides d’Orbigny, a spinose planktic foraminifer, is substantially present in temperate to subpolar water masses and is also characteristic of upwelling areas in lower latitudes [4–11]. In these upwelling regions, G. bulloides contributes maximum foraminiferal flux to the ocean floor [12, 13] and therefore provides important geochemical information for paleoceanographic studies [14–17].
Although earlier workers [18–21] discussed the general distribution of G. bulloides in Indian Ocean surface waters and sediments along with its relation to the change of ecological parameters, no attempt was made to study the morphological variations of this planktic foraminiferal species along a north-south transect in the Indian Ocean region in connection with the ecological parameters.
In the present work, the results of a biometric study of the latitudinal variation in test size, proloculus size, number of chambers, and coiling direction of G. bulloides from the surface sediments of a north-south transect of the Indian Ocean were analyzed. The major objective of this study is to link the patterns of morphological variation with the changes in physicochemical properties of the surface water in order to comprehend the ecological control on morphological characteristics in G. bulloides in the modern marine environment.
2. Study Area
The study area falls within the southwestern Indian Ocean basin. Previously, Ichiye [22], Wyrtki [23], and Gordon [24] gave a detailed account of the physical oceanography of the southern Indian Ocean. Two critical water-mass boundaries are present in this area: the Antarctic Polar Front, which separates Antarctic water mass from Subantarctic water mass, and the Subtropical Convergence, which separates Subantarctic from subtropical water masses (Figure 1). The Antarctic Convergence exhibits the sharpest water-mass boundary in the southern Indian Ocean waters in terms of the change in temperature and salinity whereas the gradients at the Subtropical Convergence (38°S–40°S) are less steep.
Location of sampling stations along a north-south transect in the southwestern Indian Ocean (NEC: North Equatorial Current; EEC: Eastward Equatorial Current; SEC: South Equatorial Current; STF: Subtropical Front; SAF: Subantarctic Front; ACC: Antarctic Circumpolar Current; PF: Polar Front).
During the pilot expedition to the Southern Ocean (PESO) aboard the Oceanic Research Vessel Sagar Kanya (199C and 200th cruises), a total of 18 surface sediment samples (0-1 cm) (including core top and grab samples) were collected along a north-south transect between 9.50°N to 45°S latitude and 80°E to 40°E longitude, covering tropical, subtropical, and Subantarctic waters in the southwestern Indian Ocean to undertake planktonic foraminiferal analysis (Figure 1, Table 1).
Average values (mean as open values and range in parenthesis) of the morphological parameters of G. bulloides for each sampling location.
Sample number
Latitude
Test size (µm) [range]
MPS (µm) [range]
Number of chambers [range]
Dextral %
SK199C/3
9.5045
268.00[192–348]
8.56[6–12]
10.33[9–12]
73.33
SK199C/4
9.4051
336.82[260–442]
15.72[10.4–6]
10.64[9–12]
54.54
SK199C/5
8.9917
241.50[180–324]
9.15[6–18]
9.63[8–12]
87.50
SK199C/6
8.1333
347.90[260–429]
14.55[13–26]
10.48[8–13]
66.66
SK199C/7
5.5121
233.33[204–276]
7.07[6–12]
10.11[9–14]
88.88
SK199C/13
−7.3648
260.14[216–312]
8.99[6–18]
10.71[9–13]
57.14
SK199C/14
−9.1790
243.20[204–312]
8.64[6–18]
10.40[8–13]
60.00
SK199C/15
−11.4243
238.00[156–325]
7.47[6–13]
9.80[8–12]
80.00
SK199C/17
−15.2785
225.35[156–336]
7.96[6–13]
11.10[8–14]
50.00
SK199C/19
−16.2677
245.05[182–312]
8.13[6.5–13]
10.05[8–13]
55.00
SK200/5
−28.3215
238.52[156–312]
8.20[6.5–13]
10.61[9–12]
52.17
SK200/9
−30.9142
271.00[195–442]
9.40[6.5–26]
10.62[8–12]
53.84
SK200/14
−36.1217
383.19[257.1–485.7]
11.89[7.1–21.4]
10.12[8–12]
73.53
SK200/15
−37.0000
259.50[180–348]
9.15[6–18]
9.79[8–13]
75.00
SK200/17
−39.0285
248.25[192–336]
7.43[6–18]
9.38[8–10]
75.00
SK200/19
−40.9813
363.05[192–520]
16.55[6–26]
10.24[8–12]
76.19
SK200/21
−43.1500
352.73[273–416]
16.73[10.4–26]
10.33[8–13]
86.66
SK200/33
−55.0065
360.29[260–520]
15.69[20.4–26]
9.57[8–12]
64.28
3. Hydrologic Setting
The sampling stations are characterized by several distinct water regions spreading over latitudinal segments: tropical, subtropical, transitional, and Subantarctic. They are divided into these waterbodies depending upon several zoogeographic provinces which are mainly influenced by ecological and climatological parameters [4, 25]. Southern Ocean can be divided into three prominent zones based on the water dynamics: the Western Boundary Current (WBC) zone (between 35° and 45°S), the Antarctic Circumpolar Current (ACC) zone (between 45° and 60°S), and the Seasonal Sea Ice (SSI) zone (between 60° and 75°S). The WBC zone comprises of mainly three currents at western boundary; they are the Agulhas Current, the Brazil/Malvinas Current, and the East Australia Current. Hydrographic conditions in Southern Ocean (SO) are mainly controlled by an eastward flowing Antarctic Circumpolar Current (ACC) [26].
The western part of the southern Indian Ocean acquires heat from the warm western boundary current [27]. The available data on the hydrological fronts and freshwater input along 62°E and 30°E sections [28] highlighted that the areas west of the Crozet Plateau and east of the Kerguelen-Amsterdam passage are two vital regions where the fronts diverse and converse. Agulhas Return Front (ARF), Southern Subtropical Front (SSTF), and Northern Subantarctic Front (NSAF) were designated as the combined front between 40°15′S and 43°S suggesting that the combined isotherms exhibit temperature variation from 19°C to 10°C and the combined isohalines exhibit a drop in salinity from 35.54 to 34.11 psu across ~3° latitude [26].
The Southern Subantarctic Front is situated between 47° and 48°S (between 6 and 7°C isotherms). The Polar Front (PF1) is present between 49° and 50°S (isotherms varied from 5 to 4°C). Southern Polar Front (PF2) is present between 52° and 54°S (temperature range 3-2°C). Antarctica Intermediate Water (AAIW) flows at ~1150 and ~1200 m water depth and the characteristics of this water mass were recorded as temperature of ~4.4°C, minimum salinity of ~34.42 psu, and density of ~27.24 kg m−3 in the northern front of subtropical zone [29].
Circumpolar Deep Water is characterized by different features, such as temperature of ~2°C, salinity of ~34.77 psu, and density of ~27.8 kg m−3. It flows at the water depth between 2000 and 3800 m north of 45°S and rises sharply to shallower depths south of the frontal zone. North Atlantic Deep Water (NADW) with higher salinities (~34.8 psu) transported from the South Atlantic to the southwestern corner of the Indian Ocean and Madagascar ridge blocks NADW to reach east of ~45°E [30]. Below the CDW, temperature and salinity decrease due to the influence of Antarctica Bottom Water (AABW). AABW is recorded between 49° and 56°S at the depth of 4100 to 4700 m with a temperature of ~−0.165 to −0.62°C, salinity of ~34.67 to 34.65 psu, and density of ~27.85 to 27.86 kg m−3 [26]. The Antarctic Circumpolar Current (ACC) reaches the ocean floor to mix with the North Atlantic Deep Water (NADW) as well as deep waters from the Indian Ocean and the Pacific Ocean. Thus, the mixture of these deep waters along with the Circumpolar Deep Waters (CDW) spreads to other oceans basins.
AABW and CDW enter the Indian Ocean in the west around Madagascar and East Africa and in the east along the Ninety East Ridge [31]. A very strong and deep overturning cell about 1800 m below near 32°S carries layers of warm near-surface water and cold deep water in opposite directions [30]. The Southern Ocean has also a unique role in the global scale overturning circulation caused by the circumpolar connection in the Southern Ocean. Water found at intermediate and abyssal depths at low latitudes rises towards the surface in the Southern Ocean. Deep water that upwells closer to Antarctica loses its heat after coming in contact with the cold air blowing off the continent and its salinity is eventually increased by brine released during sea ice formation. The dense water formed in this process gradually sinks near the continental margin of Antarctica and takes a return path to the north in deep currents flowing along the sea floor.
4. Methods
All the sediment samples were processed as per standard procedures. An appropriate amount of sediment (~5 gm) from each sample was dried overnight at 45°C. Dried sediment samples were soaked in water and subsequently treated with sodium hexametaphosphate in order to dissociate clay lumps. The treated sediments were sieved over 63 μm sieve and dried and transferred to plastic vials. While processing the sediment samples, utmost care was taken to prevent any possible breaking of the foraminiferal test, by using extremely low water pressure for removing the finer (>63 μm) fraction. The >63 μm fraction was dry sieved and, from the >125 μm fraction, an aliquot was taken by quartering and coning, to pick a minimum of 40 specimens of the planktic foraminiferal species G. bulloides.
After the picking procedure, morphological features such as coiling direction, shell size, number of chambers, and proloculus size were observed and measured under a stereo zoom microscope using a scale with divisions of 14 μm. The test size was measured along the last chamber to its diagonally opposite chamber through the apex; coiling direction was observed from dorsal side of the test, while the number of chambers was counted starting from the proloculus chamber to the last chamber observed on dorsal side of the test. The average values of these four parameters of all 40 individual specimens from each location were presented as mean values for each location (Table 1).
The data obtained was subjected to Q-mode cluster analysis. The analysis was performed using the unweighted pair group averaging method (Figure 2). The morphological parameters were chosen as variables to increase the precision of the analysis. The results of cluster analysis are plotted in the form of a two-dimensional hierarchy dendrogram wherein locations are presented along the x-axis while similarity level is plotted on y-axis (Figure 2). This dataset was compared with the modern hydrological settings (sea surface temperature, sea surface salinity, nutrients, and dissolved oxygen; mean values of 0 m, 50 m, 100 m, 150 m, and 200 m, considered the average value of 0–200 m water column) along the north-south transect in the present study area. This hydrological data was retrieved from World Ocean Atlas [32–35]. Moreover, the clusterwise comparative chart of morphological parameters and ecological parameters is tabulated to show the correlation between these two (Table 2).
Comparative morphological and ecological data for each cluster/subcluster for G. bulloides.
Clusters
Subclusters
Morphological parameters
Ecological parameters
Avg. test size (µm)
Avg. MPS (µm)
Avg. number of chambers
Avg. dextrality (%)
Avg. SST (°C)
Avg. SSS (‰)
Avg. nitrate (psu)
Avg. phosphate (psu)
Avg. dissolved O2 (mg/L)
A
357.33
15.19
10.23
70.31
14.61
34.77
4.14
4.07
5.34
B
247.65
8.34
10.21
67.32
20.55
35.03
1.89
2.39
4.82
B1
238.03
8.23
10.54
54.29
21.27
35.08
0.48
1.68
4.71
B2
240.27
7.78
9.73
82.85
21.11
34.77
0.58
3.10
4.87
B3
264.66
9.02
10.36
64.83
19.27
35.23
4.29
2.56
4.90
Results of the cluster analysis of the four morphological parameters (coiling direction, mean proloculus size (MPS), average test size, and number of chambers) analysed from the planktonic foraminiferal species Globigerina bulloides.
5. Results
The Q-mode cluster analysis classified the samples into two homogeneous clusters (A and B) under the linkage distance 40 (Figure 2). Cluster B was in turn subdivided into subclusters B1, B2, and B3 under the linkage distance 20. Each cluster and subcluster is characterized by a particular association of ecological parameters. The following are the relation of clusters and subclusters with ecological parameters.
Cluster A comprises a total of 6 samples, of which two fall between the latitudes 9.4051°N and 8.1333°N (tropical zone) with the remaining samples within the latitudes 36.12°S and 45°S (towards the subpolar zone). This cluster is characterized by the highest range of all morphological characteristics: average test size (range 336.82–383.19 μm; average 357.33 μm), mean proloculus size (range 11.89–16.73 μm; average 15.19 μm), number of chambers (range 9.57–10.64; average 10.23), and dextrality (range 54.54–76.19%; average 70.31%). Sea surface temperature (SST) and sea surface salinity (SSS) show the lowest range of values between 1.14–24.23°C (average 14.61) and 34.06–35.50 psu (average 34.77 psu), respectively. On the other hand, the nitrate, phosphate, and dissolved oxygen contents of the water show the highest values. The nitrate content ranges from 1.165 μmol to 7.81 μmol with an average of 4.14 μmol, while the phosphate content shows a range of 1.774 μmol to 8.423 μmol with an average of 4.07 μmol, and total nutrients show a range of 4.202 μmol to 11.547 μmol with an average of 8.21 μmol. The dissolved oxygen ranges from 3.886 mg/L to 7.733 mg/L with an average 5.34 mg/L (Table 2).
Cluster B (12 samples) comprises of samples from a wider range of the study area (9.5040°N to 39.03°S latitudes), which mainly falls within the tropical to subtropical zones. This cluster is characterized by the lowest range of all morphological characteristics: average test size (range 225.35–271.00 μm; average 247.65 μm), mean proloculus size (range 7.07–9.40 μm; average 8.34 μm), number of chambers (range 9.38–11.10; average 10.21), and dextral coiling (range 50.00–88.88%; average 67.32%). Two ecological parameters, namely, SST (range 15.63–23.85°C; average 20.55°C) and SSS (range 33.53–35.50 psu; average 35.03 psu), show the highest ranges. The nutrients and dissolved oxygen showed the lowest ranges. The total nutrient values varied between 1.851 μmol and 11.578 μmol (average 4.28 μmol), while the dissolved oxygen ranged from 4.431 mg/L to 5.612 mg/L (average 4.82 mg/L). The nitrate values ranged from 0.3 μmol to 9.105 μmol (average 1.89 μmol) while phosphate values varied from 1.421 μmol to 5.95 μmol (average 2.39 μmol). Cluster B is further subdivided into three subclusters, namely, B1, B2, and B3 (Table 2).
Subcluster B1 consists of four samples, located within the latitudes of 9.179°S and 28.32°S. This subcluster is characterized by the lowest range of average test size (range 225.35–245.05 μm; average 238.03 μm) and dextrality (range 50.00–60.00%; average 54.29%). On the other hand, it shows the highest range of number of chambers (range 10.05–11.10; average 10.54) and moderate range of mean proloculus size (range 7.96–8.64 μm; average 8.23 μm). In the case of ecological parameters, this cluster shows a lower range of values with the exception of SST and SSS. SST shows the highest values ranging between 19.76 and 22.76°C with an average of 21.27°C. In contrast, SSS ranges moderately within the cluster (range 34.92–35.42 psu; average 35.08 psu). The values of nitrate and phosphate varied from 0.38 to 0.7 μmol (average 0.48 μmol) and from 1.421 to 2.245 μmol (average 1.68 μmol), respectively, which collectively ranged from 1.851 to 2.945 μmol (average 2.17 μmol). Dissolved oxygen ranged from 4.543 mg/L to 4.936 mg/L, with an average of 4.71 mg/L.
Subcluster B2, represented by four samples, is distributed between latitudes 5.5121°N and 39.03°S. This subcluster is characterized by the lowest range of mean proloculus size (range 7.07–9.15 μm; average 7.78 μm), number of chambers (range 9.38–10.11; average 9.73), and the highest range of dextrality (range 75.00–88.88%; average 82.85%). On the other hand, it shows a moderate range of average test size (range 58.06–95.65%; average 82.00%) within Cluster B. In the case of ecological parameters, the SSS shows the lowest range (33.53–35.44 psu) with an average of 34.77 psu. SST and dissolved oxygen show medium ranges (15.63–23.85°C, average 21.11°C and 4.431–5.612 mg/L, average 4.87 mg/L, resp.). In the case of nutrients, nitrate shows a medium value (range 0.3–1.07 μmol, average 0.58 μmol), while phosphate shows the highest value (range 1.568–5.95 μmol; average 3.10 μmol) within Cluster B; collectively, the nutrient values range from 1.938 μmol to 7.02 μmol with an average of 3.68 μmol, which shows a medium value within Cluster B.
Subcluster B3 consists of four samples, located within the latitudes of 9.5045°N and 37°S. This subcluster is characterized by the highest range of average test size (range 259.50–271.00 μm; average 264.6 μm) and mean proloculus size (range 8.56–9.40 μm; average 9.02 μm) within Cluster B. On the other hand, it shows moderate range of number of chambers (range 9.79–10.71; average 10.36) and dextrality (range 53.84–75.00%; average 64.83%) within Cluster B. In the case of ecological parameters, this cluster shows the lowest range of SST (16.81–22.08°C; average 19.27°C) and highest range of SSS (34.94–35.50 psu; average 35.23) within Cluster B. The value of nitrate varies from 0.3 to 9.105 μmol (average 4.29 μmol), which is highest within Cluster B. Though phosphate shows a moderate range (1.687–4.271 μmol; average 2.56 μmol), collective nutrients show the highest range of values (1.987–11.578 μmol; average 6.85 μmol) within Cluster B. Dissolved oxygen shows a range of 4.495 mg/L to 5.423 mg/L, with an average of 4.90 mg/L, which is highest within Cluster B.
The correlation coefficient between the morphological variations (average test size, proloculus size, number of chambers, and dextrality) and the ecological parameters, namely, SST, SSS, nitrate content, phosphate content, total nutrient (nitrate + phosphate), and dissolved oxygen, was performed (Tables 3(a) and 3(b)). Some correlation was found moderately correlatable, namely, test size-SST, test size-nitrate (hence total nutrient), test size-dissolved oxygen, number of chambers-dissolved oxygen, proloculus size-SST, proloculus size-nitrate (hence total nutrient), and proloculus size-dissolved oxygen.
(a) Correlation coefficient equation of all morphological characters against the various ecological parameters. (b) Comparative square values of sample correlation coefficient (R2 values) between all morphological and ecological parameters.
SST
SSS
Nitrate
Phosphate
Total nutrient
Dissolved oxygen
Test size
y=(-0.058)x+35.30
y=(-0.002)x+35.41
y=0.031x-5.988
y=0.009x-0.032
y=0.040x-5.944
y=0.009x+2.499
Number of chambers
y=4.773x-3019
y=0.364x+31.22
y=0.232x+0.744
y=(-1.562)x+18.68
y=(-1.317)x+19.32
y=(-0.977)x+15.18
Proloculus size
y=(-0.926)x+28.41
y=(-0.056)x+35.54
y=0.389x-1.096
y=0.106x+1.52
y=0.492x+0.493
y=0.147x+3.58
Dextrality
y=(-0.066)x+23.11
y=(-0.013)x+35.87
y=0.026x+1.331
y=0.043x-0.259
y=0.072x+0.969
y=0.018x+3.935
SST
SSS
Nitrate
Phosphate
Total nutrient
Dissolved oxygen
Test size
0.289
0.029
0.329
0.073
0.352
0.322
Number of chambers
0.13
0.097
0.001
0.129
0.023
0.226
Proloculus size
0.292
0.14
0.205
0.039
0.209
0.341
Dextrality
0.02
0.109
0.011
0.078
0.054
0.062
6. Discussion and Conclusion
The foraminiferal data have been subjected so far to statistical analysis, namely, cluster analysis in different geographical regions and for various purposes [36–46].
G. bulloides is abundant in high southern latitude water masses and is at its peak in high northern latitudes, low latitude upwelling regions, and nutrient-rich environments [47–50]. In this study, cluster analysis on the morphological characteristics of G. bulloides is used to differentiate different water masses depending more on ecological parameters than on latitudinal gradients. In general, two major clusters (Clusters A and B) show that morphological characteristics are directly correlated to nutrients and dissolved oxygen of the ambient water mass and inversely correlated to sea surface temperature and sea surface salinity. From the combined results of all subclusters, it was observed that average test size depends directly on nitrate values and dissolved oxygen content of the ambient environment and is inversely related to the temperature. The average size of G. bulloides in Cluster A (357.33 μm) is greater than that in Cluster B (247.65 μm). Comparing the nitrate value, dissolved oxygen content, and temperature, we found that the average nitrate value and dissolved oxygen value of surface water in Cluster A (4.14 psu and 5.34 mg/L, resp.) are much higher than those in Cluster B (1.89 psu and 4.82 mg/L, resp.). In case of the three subclusters of Cluster B, the same trend was also observed (Table 1). Bé et al. [51] and Hecht [52] pioneered the study of the ecological influence on the adult test size of planktic foraminifers. In the North Atlantic region, “environmental optima” were defined for temperature and salinities for many planktonic species including G. bulloides [52]. Largest sizes of the same species were reported in the Indian Ocean at temperatures around 6-7°C [53, 54]. A positive correlation between general size and the frequencies of G. bulloides suggests that optimum growth occurs in areas of optimum environmental conditions rather than in more marginal environments where delayed reproduction (and hence greater size) might have occurred [53]. The upwelling assemblage, predominated by the species G. bulloides, has an affinity to high nutrients, because fertility is the defining characteristic for these assemblages [55].
Mean proloculus size (MPS) increases with the increase in SSS. In all subclusters under Cluster B, the average MPS shows a direct relation with the average SSS. Moreover, in our study, it is recorded that an increase in phosphate leads to a rise in dextrality and number of chambers. In Cluster A, the average dextrality and number of chambers show higher values (70.21% and 10.23, resp.) than those in Cluster B (67.32% and 10.21, resp.), which shows a direct relationship with the higher average phosphate value in Cluster A (4.07 psu) and lower value in Cluster B (2.39 psu). The dextrality is also showing an inverse dependency on the average SST. Cluster A shows higher average dextrality (70.31%) with lower average SST (14.61°C), whereas Cluster B shows lower average dextrality (67.32%) with higher average SST (20.55°C). An association between surface-water temperatures and coiling direction in living G. bulloides was reported from the southwestern Atlantic Ocean [56]. Malmgren and Kennett [53] observed a distinct relationship between the average surface-water temperature and coiling direction of G. bulloides. However, in both Antarctic Ocean and Indian Ocean, the surface-water temperature and the percentage of sinistral specimens are significantly negatively correlated.
Morphologically defined species of marine plankton often harbor a considerable level of cryptic diversity [57]. The results of present study are purely based on the cluster analyses of the morphological features of planktic species Globigerina bulloides on the assumption that each morphospecies of planktonic Foraminifera represents a genetically continuous species with a unique habitat. However, the possibility of having cryptic species in the species level population in the study material may not be ruled out completely. No doubt hidden genetic diversity among modern planktonic Foraminifera has significant repercussions on paleoproxies derived from their fossil shells. Nevertheless, Kucera and Darling [58] have compiled the genetic diversity and found 33 cryptic genetic types in 9 out of the 22 sequenced morphospecies of modern planktonic Foraminifera, implying that the total number of cryptic genetic types per morphospecies is not large and that most genetic types show a nonrandom pattern of distribution in the oceans [58]. Furthermore, Morard et al. [59] also pointed out that the cryptic genetic species of planktonic Foraminifera often exhibit narrower biogeographic distributions and ecological preferences than the respective morphospecies. In theory, it should therefore be possible to improve the resolution of the paleoceanographic reconstructions based on sediment assemblages of these species. Since many morphospecies show cosmopolitan distribution, an understanding of biogeographic and evolutionary processes at the level of genetic diversity requires global sampling [57]. Such an approach is beyond the scope of the present study.
Though the inferences drawn in the present study based on Q-mode cluster analysis clearly establish a correspondence between the ecological parameters of the ambient water masses and the morphological variables of the planktic foraminiferal species G. bulloides, more transects covering a wide geographical region need to be covered for arriving at a clearer conclusion.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful to the Director, Birbal Sahni, Institute of Paleobotany, Lucknow, and Ministry of Earth Science for providing permission to publish this paper. The authors are thankful to Ms. Rosyta Afonso for processing the samples. The authors are also thankful to Ms. Neelam Pereira for linguistic upgrade of this paper.
LewisK. B.JenkinsC.Geographical variation of Nonionellina flemingi1969151112HechtA. D.Intraspecific variation in Recent populations of Globigerinoides ruber and Globigerinoides trilobus and their application to paleoenvironmental analysis19744812171234KennettJ. P.HedleyR. H.AdamsC. G.Phenotypic variation in some recent and late Cenozoic planktonic foraminifera1976New York, NY, USAAcademic Press111170BéA. W. H.TolderlundD. S.FunnellB. M.RiedelW. R.Distribution and ecology of living planktonic foraminifera in surface waters of the Atlantic and Indian oceans1971Cambridge, UKCambridge University Press105149BéA. W. H.HutsonW. H.Ecology of planktonic foraminifera and biogeographic patterns of life and fossil assemblages in the Indian Ocean19772336941410.2307/1485406CurryW. B.OstermannD. R.GuptaM. V. S.IttekkotV.SummerhaysC. P.PrellW. L.EmeisK. C.Foraminiferal production and monsoonal upwelling in the Arabian Sea: evidence from sediment traps1992London, UKThe Geological Society93106Geological Society Special Publication No. 64PrellW. L.MurrayD. W.ClemensS. C.AndersonD. M.Evolution and variability of the Indian Ocean summer monsoon: evidence from the Western Arabian sea drilling program1992Geophysical Monograph no. 7010.1029/GM070p0447AndersonD. M.PrellW. L.A 300 kyr record of upwelling off Oman during the late Quaternary: evidence of the Asian Southwest monsoon19938219320810.1029/93PA002562-s2.0-0027384421NaiduP. D.MalmgrenB. A.Monsoon upwelling effects on test size of some planktonic foraminiferal species from the Oman Margin, Arabian Sea199510111712210.1029/94PA026822-s2.0-0028794438NaiduP. D.MalmgrenB. A.Relationship between late quaternary upwelling history and coiling properties of Neogloboquadrina pachyderma and Globigerina bulloides in the Arabian Sea1996261647010.2113/gsjfr.26.1.642-s2.0-0029665207ConanS.-H.BrummerG. J. A.Fluxes of planktic foraminifera in response to monsoonal upwelling on the Somalia Basin margin2000479–112207222710.1016/S0967-0645(00)00022-92-s2.0-0033934689SautterL. R.ThunellR. C.Seasonal succession of planktonic foraminifera: results from a four- year time-series sediment trap experiment in the Northeast Pacific198919425326710.2113/gsjfr.19.4.2532-s2.0-0024919525SautterL. R.SancettaC.Seasonal associations of phytoplankton and planktic foraminifera in an upwelling region and their contribution to the seafloor199218426327810.1016/0377-8398(92)90043-J2-s2.0-0026459419BardE.ArnoldM.MauriceP.DupratJ.MoyesJ.DuplessyJ.-C.Retreat velocity of the North Atlantic polar front during the last deglaciation determined by 14C accelerator mass spectrometry198832861337917942-s2.0-0023523726KallelN.LabeyrieL. D.Juillet-LeclercA.DuplessyJ.-C.A deep hydrological front between intermediate and deep-wafpr , masses in the glacial Indian Ocean1988333617465165510.1038/333651a02-s2.0-0023723812SautterL. R.ThunellR. C.Planktonic foraminiferal response to upwelling and seasonal hydrographic conditions: sediment trap results from San Pedro Basin, Southern California Bight199121434736310.2113/gsjfr.21.4.3472-s2.0-0026391639ChiyonobuS.MoriY.OdaM.Reconstruction of paleoceanographic conditions in the northwestern Pacific Ocean over the last 500kyr based on calcareous nannofossil and planktic foraminiferal assemblages201296-97293710.1016/j.marmicro.2012.07.0022-s2.0-84865441439BelyaevaN. V.Distribution of planktonic Foraminifera in the water and on the ocean floor in the Indian Ocean1964681283BlairD. G.The distribution of planktonic Foraminifera in deep-sea cores from the Southern Ocean, Antarctica196510BoltovskoyE.Living planktonic foraminifera at the 90°E meridian from the equator Lo the Antarctic196915237255BandyO. L.Variations in Globigerina bulloides d'Orbigny as indices of water masses19727194195IchiyeT.FairbridgeR. W.Indian ocean: oceanography1962New York, NY, USAVan Nostrand-Reinhold383396WyrtkiK.1971Washington, DC, USANational Science FoundationGordonA. L.HayesD. E.Introduction: physical oceanography of the southeast Indian Ocean1972Washington, DC, USAAmerican Geophysical Union39Antarctic Research Series 19SchottG.1942Hamburg, GermanyBoysenAnilkumarN.LuisA. J.SomayajuluY. K.BabuV. R.DashM. K.PednekarS. M.BabuK. N.SudhakarM.PandeyP. C.Fronts, water masses and heat content variability in the Western Indian sector of the Southern Ocean during austral summer 20042006631-2203410.1016/j.jmarsys.2006.04.0092-s2.0-33750433165GordonA. L.Oceanography—the brawniest retroflection200342190490510.1038/421904aParkY.-H.CharriaudE.FieuxM.Thermohaline structure of the Antarctic surface water/Winter Water in the Indian sector of the Southern Ocean1998171–452310.1016/S0924-7963(98)00026-82-s2.0-0032421756NowlinW. D.Jr.KlinckJ. M.The physics of the antarctic circumpolar current19862446949110.1029/RG024i003p00469TooleJ. M.WarrenB. A.A hydrographic section across the subtropical South Indian Ocean199340101973201910.1016/0967-0637(93)90042-22-s2.0-0027836265WarrenB. A.Indian Ocean, Deep circulation1980New York, NY, USAMcGraw-Hill227229ConkrightM. E.LevitusS.BoyerT.1994Washington, DC, USAU.S. Government Printing OfficeLevitusS.BoyerT.1994Washington, DC, USAU.S. Government Printing OfficeNOAA Atlas NESDIS 2LevitusS.BoyerT.1994Washington, DC, USAPrinting OfficeNOAA Atlas NESDIS 4LevitusS.BurgettR.BoyerT.1994Washington, DC, USAU.S. Government Printing OfficeNOAA Atlas NESDIS 3KauslerR. L.Quantitative re-evaluation of ecology and distribution of Recent foraminifera and ostracoda of Todos Santos Bay, Baja California, Mexico1966The University of Kansas Paleontological InstituteMellowJ. F.BuzasM. A.An application of cluster analysis as a method to determining biofacies196842747758UjiieH.NagaseK.Cluster analysis of living plankyonic foraminifera from the southeastern Indian OceanProceedings of the 2nd International Planktonic Conference1971Rome, Italy12511258ValentineJ. W.PeddicordR. G.Evolution of fossil assemblages by cluster analysis196711198211DriniaH.AntonarakouA.TsaparasN.Diversity and abundance trends of benthic foraminifera from the southern part of the Iraklion Basin, Central Crete200436772781NigamR.SarupriaJ. S.Cluster analysis and ecology of living benthonic foraminifera from the inner shelf off Ratnagiri, West Coast, India198122175180ReddyA. N.ReddyK. R.Cluster analysis and distribution of recent foraminifera, Araniar River estuary, Tamil Nadu, India1989337681CoulbournW. T.Relationship between the distribution of foraminifera and geologic structures of the Arctic Bight, South America198054696718UjiiéY.UjiiéH.TairaA.NakamuraT.OguriK.Spatial and temporal variability of surface water in the Kuroshio source region, Pacific Ocean, over the past 21,000 years: evidence from planktonic foraminifera200349433536410.1016/S0377-8398(03)00062-82-s2.0-0345149765TaldenkovaE.BauchH. A.StepanovaA.Benthic and planktic community changes at the North Siberian margin in response to Atlantic water mass variability since last deglacial times201296-97132810.1016/j.marmicro.2012.06.0072-s2.0-84868153932RasmussenT. L.ThomsenE.Changes in planktic foraminiferal faunas, temperature and salinity in the Gulf Stream during the last 30,000 years: influence of meltwater via the Mississippi River201233425410.1016/j.quascirev.2011.11.0192-s2.0-84856220503DuplessyJ. C.BéA. W. H.BlancP. L.Oxygen and carbon isotopic composition and biogeographic distribution of planktonic foraminifera in the Indian Ocean19813394610.1016/0031-0182(81)90031-62-s2.0-0019367265HemlebenC.SpindlerM.AndersonO. R.1989New York, NY, USASpringerThiedeJ.SuessE.ThiedeJ.Skeletal plankton and nekton in upwelling water masses off northwestern south America and northwest Africa1983New York, NY, USAPlenum Press183207ZhangB.Living planktonic foraminifera from eastern Arabian Sea198512789798BéA. W. H.HarrisonS. M.LottL.Orbulina universa d'Orbigny in the Indian Ocean19731915019210.2307/1485162HechtA. D.An ecologic model for test size variation recent planktonic foraminifera: application to the fossil record19766295311MalmgrenB.KennettJ. P.Biometric analysis of phenotypic variation in Recent Globigerina bulloides d'Orbigny in the Southern Indian Ocean1976132510.1016/0377-8398(76)90003-72-s2.0-0002345168SchmidtD. N.RenaudS.BollmannJ.SchiebelR.ThiersteinH. R.Size distribution of Holocene planktic foraminifer assemblages: biogeography, ecology and adaptation2004503-431933810.1016/S0377-8398(03)00098-72-s2.0-1242277290SchiebelR.BijmaJ.HemlebenC.Population dynamics of the planktic foraminifer Globigerina hulloides from the Eastern North Atlantic1997449-101701171310.1016/S0967-0637(97)00036-82-s2.0-0030689797BoltovskoyE.Note on the determination of absolute surface water paleotemperature by means of the foraminifer Globigerina bulloides d'Orbigny1973473-41521552-s2.0-001146302910.1007/BF02985703WeinerA. K. M.WeinkaufM. F. G.KurasawaA.DarlingK. F.KuceraM.GrimmG. W.Phylogeography of the tropical planktonic foraminifera lineage Globigerinella reveals isolation inconsistent with passive dispersal by ocean currents201493e9214810.1371/journal.pone.00921482-s2.0-84899714401KuceraM.DarlingK. F.Cryptic species of planktonic foraminifera: their effect on palaeoceanographic reconstructions2002360179369571810.1098/rsta.2001.09622-s2.0-0037090773MorardR.QuillévéréF.EscarguelG.de Garidel-ThoronT.de VargasC.KuceraM.Ecological modeling of the temperature dependence of cryptic species of planktonic Foraminifera in the Southern Hemisphere2013391133310.1016/j.palaeo.2013.05.0112-s2.0-84888198278