Compositional data analysis of Holocene sediments from the West Bengal Sundarbans, India: Geochemical proxies for grain-size variability in a delta environment

This paper is part of a special issue of Applied Geochemistry focusing on reliable applications of compositional multivariate statistical methods. This study outlines the application of compositional data analysis (CoDa) to calibration of geochemical data and multivariate statistical modelling of geochemistry and grain-size data from a set of Holocene sedimentary cores from the Ganges-Brahmaputra (G-B) delta. Over the last two decades, understanding near-continuous records of sedimentary sequences has required the use of core-scanning X-ray fluorescence (XRF) spectrometry, for both terrestrial and marine sedimentary sequences. Initial XRF data are generally unusable in ‘raw-format’, requiring data processing in order to remove instrument bias, as well as informed sequence interpretation. The applicability of these conventional calibration equations to core-scanning XRF data are further limited by the constraints posed by unknown measurement geometry and specimen homogeneity, as well as matrix effects. Log-ratio based calibration schemes have been developed and applied to clastic sedimentary sequences focusing mainly on energy dispersive-XRF (ED-XRF) core-scanning. This study has applied high resolution core-scanning XRF to Holocene sedimentary sequences from the tidal-dominated Indian Sundarbans, (Ganges-Brahmaputra delta plain). The Log-Ratio Calibration Equation (LRCE) was applied to a sub-set of core-scan and conventional ED-XRF data to quantify elemental composition. This provides a robust calibration scheme using reduced major axis regression of log-ratio transformed geochemical data. Through partial least squares (PLS) modelling of geochemical and grain-size data, it is possible to derive robust proxy information for the Sundarbans depositional environment. The application of these techniques to Holocene sedimentary data offers an improved methodological framework for unravelling Holocene sedimentation patterns.

Effective grain size distribution analysis for interpretation of tidal-deltaic facies: West Bengal Sundarbans

Research over the past two decades on the Holocene sediments from the tide dominated west side of the lower Ganges delta has focussed on constraining the sedimentary environment through grain size distributions (GSD). GSD has traditionally been assessed through the use of probability density function (PDF) methods (e.g. log-normal, log skew-Laplace functions), but these approaches do not acknowledge the compositional nature of the data, which may compromise outcomes in lithofacies interpretations. The use of PDF approaches in GSD analysis poses a series of challenges for the development of lithofacies models, such as equifinal distribution coefficients and obscuring the empirical data variability. In this study a methodological framework for characterising GSD is presented through compositional data analysis (CODA) plus a multivariate statistical framework. This provides a statistically robust analysis of the fine tidal estuary sediments from the West Bengal Sundarbans, relative to alternative PDF approaches.

Moleculare composition of sediments from ODP Hole 116-717C

We describe the molecular composition of a portion of the solvent-soluble organic material (lipid extract), from three organic rich muds (samples 116-717C-22X-1, 80-86 cm, 116-717C-34X-3, 130-135 cm, and 116-717C-55X-1, 65-70 cm). These samples were taken from Hole 717C, located on the Bengal Fan at a position of 0°55.8'S and 81°23.4'E. Both the palaeoenvironmental and digenetic significance of these lipid distributions have been assessed and found to be consistent with their suspected origins, i.e., turbidites from the upper slope of the western Bay of Bengal and the Ganges-Brahmaputra delta.

Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Chemical and isotopic compositions of pore fluids from ODP Sites 169-1034 and 169-1037

We have analysed the concentrations of Li, K, Rb, Cs, and B, and the isotopic ratios of Li and B of a suite of pore fluids recovered from ODP Sites 1037 (Leg 169; Escanaba Trough) and 1034 (Leg 169S; Saanich Inlet). In addition, we have analysed dissolved K, Rb, and Cs concentrations for estuarine mixing of the Ganges-Brahmaputra river system. Together, these data sets have been used to assess the role of sediments in the marine geochemical cycles of the alkali elements and boron. Uptake onto clay minerals during estuarine mixing removes 20-30% of the riverine input of dissolved Cs and Rb to the oceans. Prior to this study, the only other recognised sink of Rb and Cs was uptake during low-temperature alteration of the oceanic crust. Even with this additional sink there is an excess of inputs over outputs in their modern oceanic mass balance. Pore fluid data show that Li and Rb are transferred into marine sediments during early diagenesis. However, modeling of the Li isotope systematics of the pore fluids from Site 1037 shows that seawater Li taken up during marine sedimentation can be readily returned to solution in the presence of less hydrated cations, such as NH4+. This process also appears to result in high concentrations of pore fluid Cs (relative to local seawater) due to expulsion of adsorbed Cs from cation exchange sites. Flux calculations based on pore fluid data for a series of ODP sites indicate that early diagenesis of clay sediments removes around 8% of the modern riverine input of dissolved Li. Although NH4+-rich fluids do result in a flux of Cs to the oceans, on the global scale this input only augments the modern riverine Cs flux by ~3%. Nevertheless, this may have implications for the fate of radioactive Cs in the natural environment and waste repositories.

Bulk and amino acid parameters of core SO188-342KL from the Northern Bay of Bengal for the last 18 ka

The Northern Bay of Bengal (NBoB) is a globally important region for deep-sea organic matter (OM) deposition due to massive fluvial discharge from the Ganges-Brahmaputra-Meghna (G-B-M) rivers and moderate to high surface productivity. Previous studies have focused on carbon burial in turbiditic sediments of the Bengal Fan. However, little is known about the storage of carbon in pelagic and hemipelagic sediments of the Bay of Bengal over millennial time scales. This study presents a comprehensive history of OM origin and fate as well as a quantification of carbon sediment storage in the Eastern Bengal Slope (EBS) during the last 18 ka. Bulk organic proxies (TOC, TIC, TN, d13CTOC, d15NTN) and content and composition of total hydrolysable amino acids (THAA) in a sediment core (SO188-342KL) from the EBS were analyzed. Three periods of high OM accumulation were identified: the Late Glacial (LG), the Bölling/Alleröd (B/A), and the Early Holocene Climatic Optimum (EHCO). Lower eustatic sea level before 15 ka BP allowed a closer connection between the EBS and the fluvial debouch, favoring high terrestrial OM input to the core site. This connection was progressively lost between 15 and 7 ka BP as sea level rose to its present height and terrestrial OM input decreased considerably. Export and preservation of marine OM was stimulated during periods of summer monsoon intensification (B/A and EHCO) as a consequence of higher surface productivity enhanced by cyclonic-eddy nutrient pumping and fluvial nutrient delivery into the photic zone. Changes in the THAA composition indicate that the marine plankton community structure shifted from calcareous-dominated before 13 ka BP to siliceous-dominated afterwards. They also indicate that the relative proportion of marine versus terrestrial OM deposited at site 342KL was primarily driven by relative sea level and enlarged during the Holocene. The ballasting effect of lithogenic particles during periods of high coastal proximity and/or enhanced fluvial discharge promoted the export and preservation of OM. The high organic carbon accumulation rates in the EBS during the LG (18-17 ka BP) were 5-fold higher than at present and comparable to those of glacial upwelling areas. Despite the differences in sediment and OM transport and storage among the Western and Eastern sectors of the NBoB, this region remains important for global carbon sequestration during sea level low-stands. In addition, the summer monsoon was a key promotor of terrestrial and marine OM export to the deep-ocean, highlighting its relevance as regulator of the global carbon budget.

A modeling study of the processes of surface salinity seasonal cycle in the Bay of Bengal

In response to the Indian Monsoon freshwater forcing, the Bay of Bengal exhibits a very strong seasonal cycle in sea surface salinity (SSS), especially near the mouths of the Ganges-Brahmaputra and along the east coast of India. In this paper, we use an eddy-permitting (similar to 25 km resolution) regional ocean general circulation model simulation to quantify the processes responsible for this SSS seasonal cycle. Despite the absence of relaxation toward observations, the model reproduces the main features of the observed SSS seasonal cycle, with freshest water in the northeastern Bay, particularly during and after the monsoon. The model also displays an intense and shallow freshening signal in a narrow (similar to 100 km wide) strip that hugs the east coast of India, from September to January, in good agreement with high-resolution measurements along two ships of opportunity lines. The mixed layer salt budget confirms that the strong freshening in the northern Bay during the monsoon results from the Ganges-Brahmaputra river discharge and from precipitation over the ocean. From September onward, the East India Coastal Current transports this freshwater southward along the east coast of India, reaching the southern tip of India in November. The surface freshening results in an enhanced vertical salinity gradient that increases salinity of the surface layer by vertical processes. Our results reveal that the erosion of the freshwater tongue along the east coast of India is not driven by northward horizontal advection, but by vertical processes that eventually overcome the freshening by southward advection and restore SSS to its premonsoon values. The salinity-stratified barrier layer hence only acts as a ``barrier'' for vertical heat fluxes, but is associated with intense vertical salt fluxes in the Bay of Bengal.

Os and Re concentrations and Os isotope ratios for DSDP and ODP metalliferous sediments from the Pacific Ocean (Table 1)

We report new 187Os/186Os data and Re and Os concentrations in metalliferous sediments from the Pacific to construct a composite Os isotope seawater evolution curve over the past 80 m.y. Analyses of four samples of upper Cretaceous age yield 187Os/186Os values of between 3 and 6.5 and 187Re/186Os values below 55. Mass balance calculations indicate that the pronounced minimum of about 2 in the Os isotope ratio of seawater at the K-T boundary probably reflects the enormous input of cosmogenic material into the oceans by the K-T impactor(s). Following a rapid recovery to 187Os/186Os of 3.5 at 63 Ma, data for the early and middle part of the Cenozoic show an increase in 187Os/186Os to about 6 at 15 Ma. Variations in the isotopic composition of leachable Os from slowly accumulating metalliferous sediments show large fluctuations over short time spans. In contrast, analyses of rapidly accumulating metalliferous carbonates do not exhibit the large oscillations observed in the pelagic clay leach data. These results together with sediment leaching experiments indicate that dissolution of non-hydrogenous Os can occur during the hydrogen peroxide leach and demonstrate that Os data from pelagic clay leachates do not always reflect the Os isotopic composition of seawater. New data for the late Cenozoic further substantiate the rapid increase in the 187Os/186Os of seawater during the past 15 Ma. We interpret the correlation between the marine Sr and Os isotope records during this time period as evidence that weathering within the drainage basin of the Ganges-Brahmaputra river system is responsible for driving seawater Sr and Os toward more radiogenic isotopic compositions. The positive correlation between 87Sr/86Sr and U concentration, the covariation of U and Re concentrations, and the high dissolved Re, U and Sr concentrations found in the Ganges-Brahmaputra river waters supports this interpretation. Accelerating uplift of many orogens worldwide over the past 15 Ma, especially during the last 5 Ma, could have contributed to the rapid increase in 187Os/186Os from 6 to 8.5 over the past 15 Ma. Prior to 15 Ma the marine Sr and Os record are not tightly coupled. The heterogeneous distribution of different lithologies within eroding terrains may play an important role in decoupling the supplies of radiogenic Os and Sr to the oceans and account for the periods of decoupling of the marine Sr and Os isotope records.

Osmium isotope composition of leachates and bulk sediment from ODP Holes in the Bengal Fan

Os isotopic compositions and OS and Re concentrations were measured in H2O2-H2SO4 leachates and bulk sediment samples from Holes 717C and 718C of ODP Leg 116 in the Bengal Fan. Os isotopic results indicate that, at the sediment surface, the leachable Os fraction is derived from seawater. In contrast, leachable Os from Ganges River sediments has 187Os/188Os ratios (Pegram et al., 1994, doi:10.1016/0012-821X(94)90172-4) much higher than the marine value. This difference suggests that the leachable radiogenic Os carried by the river sediments is completely released to the oceans prior to sediment deposition in the Fan. A simple calculation, assuming these sediments to be typical of those delivered by the Ganges-Brahmaputra river system, suggests that this process can account for a substantial part of the rise in the seawater Os isotopic ratio observed over the past 16 m.y. Bengal Fan leachate 187Os/188Os ratios increase with increasing depositional age, in contrast to the seawater Os isotopic ratio, which decreases with increasing age. Several lines of evidence suggest that, at the time of sediment burial, the leachate Os compositions most likely reflected the seawater values. Thus, the current divergence is probably the result of post-depositional processes. One such process, in situ radiogenic ingrowth of 187Os, can be excluded because the measured Re concentrations of these sediments are too low. Similarly, since most of the bulk rock Os isotopic ratios were lower than those of the associated leachates, the high leachate 187Os/188Os values cannot be explained by in situ sediment alteration. Instead, it is proposed that the increase with age results from radiogenic OS brought in by thermoconvective circulation from further upslope in the Fan. The ultimate source of this 187Os would then be alteration of radiogenic sediments or post-depositional radioactive decay of Re in sediments rich in organic carbon. Finally, the divergence between the results obtained on Bengal Fan sediments and those obtained in the open ocean (Pegram et al., 1992, doi:10.1016/0012-821X(92)90132-F) by the same leaching technique suggest that Os sediment leachate data must be interpreted with caution.

(Table 1) Water storage changes of the 33 world's largest river basins between 2002-2009

Global change in land water storage and its effect on sea level is estimated over a 7-year time span (August 2002 to July 2009) using space gravimetry data from GRACE. The 33 World largest river basins are considered. We focus on the year-to-year variability and construct a total land water storage time series that we further express in equivalent sea level time series. The short-term trend in total water storage adjusted over this 7-year time span is positive and amounts to 80.6 ± 15.7 km**3/yr (net water storage excess). Most of the positive contribution arises from the Amazon and Siberian basins (Lena and Yenisei), followed by the Zambezi, Orinoco and Ob basins. The largest negative contributions (water deficit) come from the Mississippi, Ganges, Brahmaputra, Aral, Euphrates, Indus and Parana. Expressed in terms of equivalent sea level, total water volume change over 2002-2009 leads to a small negative contribution to sea level of -0.22 ± 0.05 mm/yr. The time series for each basin clearly show that year-to-year variability dominates so that the value estimated in this study cannot be considered as representative of a long-term trend. We also compare the interannual variability of total land water storage (removing the mean trend over the studied time span) with interannual variability in sea level (corrected for thermal expansion). A correlation of ~0.6 is found. Phasing, in particular, is correct. Thus, at least part of the interannual variability of the global mean sea level can be attributed to land water storage fluctuations.

Mineralogy of northern Indian Ocean surface sediments and of suspended sediments of Indian rivers (Table 2 and 5)

The solid phases from surface sediments, atmospheric dusts, and rivers of the Indian Ocean environment have been analyzed for their clay minerals and quartz. Such data have been used to delimit the transport paths and sources of the detrital minerals in the oceanic deposits. Diagnostic in distinguishing fluvial and eolian inputs to the northern Indian Ocean is a combination of the clay mineral assemblages and of their geographic distributions. River borne solids are the primary components of the Bay of Bengal deposits. The eastern part receives its continental input through the Ganges-Brahmaputra river system, while drainage of the Indian Peninsula by rivers introduces solids to the western part. The former materials are characterized by high illite and chlorite in the clay mineral assemblages; the latter by montmorillonite. The winds over the Bay bear distinctive dust burdens based upon their directions. However, their contributions to the sediments are insignificant. The eastern sector of the Arabian Sea receives major contributions of continental debris from the rivers and the high montmorillonite levels clearly indicate a source in the Indian Peninsula. The rest of the Sea appears to receive most of its land-derived materials from the north, perhaps the desert regions of northern India and West Pakistan, and they are wind-borne. These materials are also transported to the equatorial regions of the Indian Ocean. A gradient in attapulgite, just north of the equator, may indicate an eolian contribution to the Arabian Sea from the African continent. The halogenated hydrocarbon pesticides were assayed in the southwest monsoon winds and enter the Bay of Bengal at levels of a half ton per month, an amount comparable to those introduced by other wind and river systems to the marine environment.

Physical properties of 39 sediment cores from the Bengal Fan and northeast Indian Ocean

We obtained sediment physical properties and geochemical data from 47 piston and gravity cores located in the Bay of Bengal, to study the complex history of the Late Pleistocene run-off from the Ganges and Brahmaputra rivers and its imprint on the Bengal Fan. Grain-size parameters were predicted from core logs of density and velocity to infer sediment transport energy and to distinguish different environments along the 3000-km-long transport path from the delta platform to the lower fan. On the shelf, 27 cores indicate rapidly prograding delta foresets today that contain primarily mud, whereas outer shelf sediment has 25% higher silt contents, indicative of stronger and more stable transport regime, which prevent deposition and expose a Late Pleistocene relic surface. Deposition is currently directed towards the shelf canyon 'Swatch of No Ground', where turbidites are released to the only channel-levee system that is active on the fan during the Holocene. Active growth of the channel-levee system occurred throughout sea-level rise and highstand with a distinct growth phase at the end of the Younger Dryas. Coarse-grained material bypasses the upper fan and upper parts of the middle fan, where particle flow is enhanced as a result of flow-restriction in well-defined channels. Sandier material is deposited mainly as sheet-flow deposits on turbidite-dominated plains at the lower fan. The currently most active part of the fan with 10-40 cm thick turbidites is documented for the central channel including inner levees (e.g., site 40). Site 47 from the lower fan far to the east of the active channel-levee system indicates the end of turbidite sedimentation at 300 ka for that location. That time corresponds to the sea-level lowering during late isotopic stage 9 when sediment supply to the fan increased and led to channel avulsion farther upstream, probably indicating a close relation of climate variability and fan activity. Pelagic deep-sea sites 22 and 28 contain a 630-kyear record of climate response to orbital forcing with dominant 21- and 41-kyear cycles for carbonate and magnetic susceptibility, respectively, pointing to teleconnections of low-latitude monsoonal forcing on the precession band to high-latitude obliquity forcing. Upper slope sites 115, 124, and 126 contain a record of the response to high-frequency climate change in the Dansgaard-Oeschger bands during the last glacial cycle with shared frequencies between 0.75 and 2.5 kyear. Correlation of highs in Bengal Fan physical properties to lows in the d18O record of the GISP2 ice-core suggests that times of greater sediment transport energy in the Bay of Bengal are associated with cooler air temperatures over Greenland. Teleconnections were probably established through moisture and other greenhouse-gas forcing that could have been initiated by instabilities in the methane hydrate reservoir in the oceans.

Discarge compilation from The Global River Discharge (RivDIS) Project

The Global River Discharge (RivDIS) data set contains monthly discharge measurements for 1018 stations located throughout the world. The period of record varies widely from station to station, with a mean of 21.5 years. These data were digitized from published UNESCO archives by Charles Voromarty, Balaze Fekete, and B.A. Tucker of the Complex Systems Research Center (CSRC) at the University of New Hampshire. River discharge is typically measured through the use of a rating curve that relates local water level height to discharge. This rating curve is used to estimate discharge from the observed water level. The rating curves are periodically rechecked and recalibrated through on-site measurement of discharge and river stage.