Provenance of Bengal Shelf Sediments: 2. Petrology and Geochemistry of Sand

The Bangladesh lowlands are traversed by the largest sediment flux on the planet. Detritus generated mostly in Himalayan highlands and conveyed through the Ganga–Brahmaputra rivers and Meghna estuary reaches the Bay of Bengal, where it forms a composite deltaic system. This study integrates the vast existing database on Ganga–Brahmaputra sediments of all grain sizes from clay to sand with new petrographic, mineralogical, and geochemical data on estuarine and shallow-marine sands. A large spectrum of compositional signatures was used to: (i) assess the relative supply of the Ganga and Brahmaputra rivers to estuarine and shelfal sediments; (ii) define the compositional variability of estuarine sediments and the impact exerted by hydraulic sorting and climate-related chemical weathering on provenance signals; (iii) define the compositional variability of shelf sediments and the potential hydrodynamic segregation of fast-settling heavy minerals in coastal environments and of slow-settling platy micas on low-energy outer-shelf floors; (iv) consider the potential additional mud supply from the western subaerial part of the delta formerly built by the Ganga River; and (v) draw a preliminary mineralogical comparison between fluvio-deltaic sediments and turbidites of the Bengal–Nicobar deep-sea fan, thus tracing sediment dispersal across the huge sedimentary system extending from Tibet to the equatorial Indian Ocean. All investigated mineralogical and geochemical parameters, as well as Sr and Nd isotope ratios and clay–mineral assemblages, showed a clear prevalence in sediment supply from the Brahmaputra (60–70%) over the Ganga (30–40%). Heavy-mineral suites and Sr and Nd isotope fingerprints of Bengal shelf sediments are nearly identical to those of the Brahmaputra River and Meghna estuary, also because the Brahmaputra carries almost twice as many Ca-plagioclase grains and heavy minerals including epidote than the Ganga, and these minerals control the large majority of the Sr and Nd budgets. The experience gained in modern settings can be directly extrapolated only to the recent past, because sediments older than the late Pleistocene and buried more than a few hundred meters begin to lose less durable ferromagnesian minerals by selective chemical dissolution, which makes quantitative estimates progressively less robust in more deeply buried older strata.


The Bengal Delta and the Shelf
The Ganga and Brahmaputra rivers join~200 km upstream of the coastline forming a~100 km-long river tract named Padma (the name of the Ganga in Bangladesh). After receiving the Meghna River, which drains the Shillong Plateau and the frontal part of the Indo-Burman Ranges, this huge river system debouches into the Bengal Sea in the~100 km long, tide-dominated Lower Meghna estuary with four channels named, from west to east, Tetulia, Shahbazpur, Hatia, and Sandwip ( Figure 1). To the west of the Meghna estuary lies the Sundarban mangrove forest, formed during progressive eastward migration of the Ganga during the late Holocene highstand ( [55]).
Until the late 18th century, as shown in the map by Rennell [56] (Figure 19.4 in [54]), the Ganga discharged through the Tetulia channel, whereas the Brahmaputra flowed to the east of the Madhupur terrace and across the Sylhet basin to reach the sea joint with the Meghna River ( Figure 1). This Old Brahmaputra (Burrampooter) course was progressively abandoned during catastrophic floods between the latest 18th century and the beginning of the 19th century, and the Ganga and Brahmaputra became eventually joined around 1830 ( [57,58]).
The Ganga-Brahmaputra delta lies on an active plate-boundary and sits on the continent-ocean transition at the northern edge of the Indian plate, which is subducting under the Himalaya in the north and under the Tripura (Indo-Burman) fold-belt in the east ( [5,59]). The delta has thus long been subjected to strong tectonic subsidence, which has determined an accumulation of~16 km of orogenic sediment since Paleocene collision between India and Asia ( [3,60,61]). Despite continuing subsidence ( [62]; Figure 9 in [63]), huge sediment supply was sufficient to halt transgression during rapid sea-level rise in the early Holocene, and to initiate deltaic progradation at rates notably higher than most other deltas worldwide since~11 ka, leading to seaward coastal advance by 100-300 km over a 250 km wide front after maximum sea level was reached at 8-7 ka ( [64]). In the last 6-5 ka, the lower delta plain west of the modern estuary was deposited as a Ganga-dominated delta, whereas Brahmaputra influence remained confined to the eastern area of the Meghna estuary and to the supratidal part of the delta (Figure 13 in [65]; Figure 9 in [55]). Active sedimentation has continued until today in the estuary area, where on average 7 km 2 of new land have been added each year since 1792 (4.4 km 2 /a since 1840) while the inactive western delta was undergoing tidal reworking and net erosion with as much as 3-4 km of shoreline retreat ( [66]; Figure 2b in [6]). The Ganga-Brahmaputra delta lies on an active plate-boundary and sits on the continent-ocean transition at the northern edge of the Indian plate, which is subducting under the Himalaya in the north and under the Tripura (Indo-Burman) fold-belt in the east ( [5,59]). The delta has thus long been subjected to strong tectonic subsidence, which has determined an accumulation of ~16 km of orogenic sediment since Paleocene collision between India and Asia ( [3,60,61]). Despite continuing subsidence ( [62]; Figure 9 in [63]), huge sediment supply was sufficient to halt transgression during rapid sealevel rise in the early Holocene, and to initiate deltaic progradation at rates notably higher than most other deltas worldwide since ~11 ka, leading to seaward coastal advance by 100-300 km over a 250 km wide front after maximum sea level was reached at 8-7 ka ( [64]). In the last 6-5 ka, the lower delta plain west of the modern estuary was deposited as a Ganga-dominated delta, whereas Brahmaputra influence remained confined to the eastern area of the Meghna estuary and to the supratidal part of the delta (Figure 13 in [65]; Figure 9 in [55]). Active sedimentation has continued until today in the estuary area, where on average 7 km 2 of new land have been added each year since 1792 (4.4 km 2 /a since 1840) while the inactive western delta was undergoing tidal reworking and net erosion with as much as 3-4 km of shoreline retreat ( [66]; Figure 2b in [6]). Offshore, a large subaqueous delta is prograding seaward between the Hatia trough in the east and the Swatch of No Ground in the west, the gigantic submarine canyon that connected the Ganga-Brahmaputra mouth with the Bengal Fan during lowstands ( Figure 1). The canyon deeply incised into the shelf, acts as a barrier for farther westward sediment transport and as a sink for approximately a third of fluvial sediment load ( [5]). Beyond the estuary mouth, sandy topsets are wide and extend to~100 km offshore between water depths of~5 m and 15-30 m b.s.l. The slope steepens in delta-front foresets between water depths of 20-30 and~80 m b.s.l., where grain size decreases from very fine sand to fine silt ( [67]). The outer shelf, extending between 80 and 150 m b.s.l. at the shelf edge, includes muddy prodelta bottomsets passing seaward to a discontinuous layer of palimpsest sand and mud containing biogenic debris, originally deposited in tidal-flat and delta-plain environments during the last glacial period and reworked during the early Holocene transgression ( [68,69]).
Tidal amplitude is~2 m in the Sundarbans to the west, where saline waters may penetrate inland by more than 100 km, and then increases eastward in the Tetulia and Shahbazpur channels to reach~5 m in the Hatia and Sandwip channels ( [65]). Tidal range decreases in the inner shelf from 0.8 m on topset beds to 0.5 m on foreset and bottomset beds and becomes negligible in the outer shelf ( [70]). Agents of sediment transport and resuspension in the inner shelf are not only tides but largely tropical cyclones, which generate high waves that sweep the shelf mobilizing large amounts of sediment redeposited as graded beds ( [68]). In the November 1970 typhoon, winds of 220 km/h and a storm surge of 6-10 m recorded as far as 200 km inland resulted in 500,000 casualties in Bangladesh ( [71,72]). Although the number and power of cyclones has seemingly decreased in the last decades, cyclones have continued to dominate sediment distribution across the shelf, with a main southwestward direction of transport toward the Swatch of No Ground ( [7]).
Sedimentation rates have been assessed at 2-4 cm/a on the delta plain in the Holocene ( [73]), increasing to a maximum of 5-10 cm/a on the delta front, and dropping to 0.3 cm/a on the prodelta farther offshore ( [68,74]). Over the last centuries, the decline in monsoon precipitation contributed to a progressive decrease in the percentage of fluvial sediment load stored in the eastern part of the subaqueous delta, from~22% in the 19th century, to~18% in the first half of the 19th century, and to 11-17% thereafter. The depocentre is apparently shifting toward the western part of the subaqueous delta and to the Swatch of No Ground canyon, thus increasing the export as turbidity currents to the deep-sea fan ( [67,75]).

Analytical Methods
To investigate provenance of Bengal shelf sand and quantify the compositional differences between fluvial and marine sediments we analyzed the petrographic, heavy-mineral, and geochemical compositions of samples SO188_307VC 0-20, SO188_316VC 0-100, and SO188_347VC 280-300 (named SO307, SO316, and SO347 for brevity throughout the article (Figure 1 These new results allowed comparison with mineralogical and geochemical data on Bengal shelf mud illustrated in the companion paper by [19], with the previously obtained extensive mineralogical dataset on fluvial bedload and suspended load of the Ganga-Brahmaputra system illustrated in [12,13], and with our own data on sediment samples from the Padma River (S3562, BR820, and vertical profile BR524-BR529), the Meghna River (BGP36 and BR826), the Meghna estuary (BR219 and vertical profiles BR439-BR448), the delta, and the shelf. Bengal delta sediments were cored on land~20 km NNE of Kolkata (S3623; luminesce age~7.764 years, Andy Carter written communication, 2006) and on Bhola Island between the Tetulia and Shahbazpur channels (SBZ#1), or collected from a tidal distributary channel (Pussur BR434, BR435) and from Katka beach in the Sundarbans (BR809). Samples SO93_63KL 3-7 and SO93_105KL 686-690 (named SO63 and SO105 for brevity throughout the article) were collected by piston corer on the Bengal shelf during Sonne cruise 93 in January 1994, at water depths of 66 and 80 m b.s.l. ( [77]). Previously obtained results, including multiple-window petrographic and heavy-mineral data obtained separately for all significant size classes of two selected sand samples from the Padma River and Meghna estuary (S3562 collected during the dry season and BR446 collected during the wet monsoon season), and multiple-window heavy-mineral data on Meghna estuary shallow suspended load (sample BR448) and Bengal shelf sandy silt (sample SO63) allowed us to evaluate the compositional variability related to markedly seasonal transport, as well as various sources of error during sampling, sample treatment, and analysis.
Statistical techniques used to illustrate our petrographic and heavy-mineral datasets include multidimensional scaling, which produces a map of points in which samples with similar mineralogical signature cluster closely together and dissimilar samples plot far apart ( [78,79]) and the compositional biplot ( [80]). The biplot is a very efficient statistical/graphical tool that allows discrimination among multivariate observations (points) while shedding light on the mutual relationships among variables (rays). The length of each ray is proportional to the variance of the corresponding variable in the dataset. If the angle between two rays is close to 0 • , 90 • or 180 • , then the corresponding variables are directly correlated, uncorrelated, or inversely correlated, respectively. Key compositional parameters for selected samples are summarized in Table 1. Full information on sampling locations (Table S1 and Google Earth TM file Bengalsand.kmz) and the complete petrographic (Table S2), heavy-mineral (Table S3), and geochemical datasets (Table S4) are provided in the Supplementary Materials.

Petrography
A quartered aliquot of each bulk sample was impregnated with araldite, cut into a standard thin section stained with alizarine red to distinguish dolomite and calcite, and analyzed by counting 400 or 450 points under the microscope (Gazzi-Dickinson method; [81]). Sand samples were classified according to the relative abundance of the three main groups of framework components (Q = quartz; F = feldspars; L = lithic fragments), considered where exceeding 10% QFL. According to standard use ( [82,83]), the less abundant component goes first, the more abundant last (e.g., in a litho-feldspatho-quartzose sand Q > F > L> 10% QFL; classification scheme after [84,85]). The average rank of metamorphic rock fragments was expressed by the metamorphic indices MI or MI*, ranging, respectively, from 0 (detritus from sedimentary and volcanic rocks) or 100 (detritus from very low-grade metamorphic rocks) to 500 (detritus from high-grade metamorphic rocks; [86]). The median grain size of sand samples was determined both by sieving and in thin section by ranking and visual comparison with standards of φ/4 classes prepared in our laboratory.

Geochemistry
Bulk-sediment chemical analyses were carried out at Service d'Analyse des Roches et des Minéraux (SARM-CRPG, Nancy, France). Sediment samples were first powdered in an agate mortar. Element concentrations were measured by inductively-coupled-plasma optical emission spectrometry (ICP-OES) and mass spectrometry (ICP-MS) on bulk aliquots of~100 mg of sediment after lithium metaborate fusion ( [93,94]). The relative uncertainty for major element concentration is <5%, except for Si, Al, and Fe for which it is~2% (Table S5). For full information on analytical procedures and geostandards used see [94]. Grain-size subclasses split by sieving at 0.50 φ were previously analyzed at ACME Laboratories Vancouver after lithium metaborate/tetraborate fusion and nitric acid digestion; major oxides and several minor elements were determined by ICP-ES, and trace elements by ICP-MS ( [12] [98,99]) in our samples and in the Upper Continental Crust standard (UCC; [102,103]). The comparison among these different chemical indices allowed us to assess how they differently respond to controls other than weathering, including source-rock mineralogy, hydraulic sorting, grain size, and quartz addition ( [104]).

Results
In this section, we illustrate the petrographic, heavy-mineral, and geochemical data newly obtained on Meghna estuary and Bengal shelf sands.

Petrography and Heavy Minerals in Bengal Shelf Sand
The siliciclastic fraction of the studied samples cored in the subaqueous delta decreases seaward in grain size from fine sand on topsets (SO307; Figure 2G) to micaceous silty sand (SO316) and very coarse silt in upper foresets (SO105; Figure 2H), whereas it may be associated with ooids and bioclasts up to 300 µm in diameter in the outer shelf (SO347; Figure 2I). Bioclasts, commonly occurring at the core of ooids, include mostly foraminifera (miliolids and rotaliids with less common textularids, bolivinids, and rare globigerinids) and fragments of gastropods or echinoids ( [19]). The chambers of foraminiferal tests are commonly filled with glaucony.

Petrography and Heavy Minerals in Bengal Shelf Sand
The siliciclastic fraction of the studied samples cored in the subaqueous delta decreases seaward in grain size from fine sand on topsets (SO307; Figure 2G) to micaceous silty sand (SO316) and very coarse silt in upper foresets (SO105; Figure 2H), whereas it may be associated with ooids and bioclasts up to 300 µm in diameter in the outer shelf (SO347; Figure 2I). Bioclasts, commonly occurring at the core of ooids, include mostly foraminifera (miliolids and rotaliids with less common textularids, bolivinids, and rare globigerinids) and fragments of gastropods or echinoids ( [19]). The chambers of foraminiferal tests are commonly filled with glaucony.

Figure 3.
Petrography and heavy minerals in the Bengal sediment system (data for Ganga-Brahmaputra sand and Bengal shelf mud after [12,19]; data for Bengal-Nicobar Fan turbidites after [20][21][22][23]). Although detrital modes are controlled by several factors including seasonal and annual variability, hydraulic sorting, and grain size, litho-feldspatho-quartzose Brahmaputra sand is discriminated from mostly feldspatho-litho-quartzose Ganga sand in QFL (A) and LmLvLs (B) triangular diagrams, as well as in biplots (C: petrography; D: heavy minerals; E: petrography and heavy minerals combined). Padma River and Meghna estuary sands range from almost pure Brahmaputra to mixed Ganga-Brahmaputra, whereas sediments of the subaerial delta are akin to the Ganga in the west and to the Brahmaputra in the east. Bengal shelf and fan sediments are either intermediate or close to the Brahmaputra. See Section 3 for an explanation of the biplot (drawn using CoDaPack software by [106]). Petrographic and mineralogical parameters are as in Table 1.

Geochemistry of Meghna Estuary and Bengal Shelf Sand
The estuary and shelf sands have quite similar major element concentrations (mostly SiO 2 Table S4 and [17]). Estuary sand tends to be richer in Ti, Y, REEs (rare earth elements), Th, Nb, Ta, Cr, Mn, Cd, Sn, and Bi, reflecting concentration of ultra-dense minerals including garnet, monazite, and Fe-Cr-Ti oxides. Shelf sand, instead, displays higher Cs, Cu, As, and LOI; K and Rb also tend to be slightly higher, reflecting greater abundance of mica.
Among estuary sands, the coarsest-grained sample BR1051 is slightly enriched in Si and depleted in other major (Al, Fe, Mg, Ca, Ti, P) and minor elements (YREEs, Nb, Cr, Co, Ni). Conversely, sample BR8230 is depleted in Si and enriched in Fe, Mg, Ca, Ti, P, Mn, Sr, Ba, YREE, Th, U, Zr, Hf, V, Nb, Ta, Cr, Co, Ni, Zn, Cd, Sb, and Bi, which is largely explained by the much higher concentration in heavy minerals and less quartz and feldspar (Table 1). Among shelf sands, sample SO347 is rich in ooids and bioclasts ( Figure 2I), which explains its anomalously high CaO and LOI (4.4% and 4.9%, respectively) and relatively high Sr. Instead, Al, Fe, Na, K, Rb, Ba, and Pb are low.

The Bengal Sediment System: Processes and Products
In this section, we integrate our new mineralogical and geochemical data on the Meghna estuary and Bengal shelf sand and mud ( [19]) with previously obtained data from the Ganga and Brahmaputra rivers and Bengal delta ( [8,9,[11][12][13][14][15]17,107,108]). We also illustrate the strong variability of sediment composition observed in various parts of the Bengal sediment system and discuss specifically the control exerted by grain size and diverse hydrodynamic processes. These include suspension sorting (i.e., the partitioning of detrital grains at different depths in the water column according to their settling velocity, which depends in turn on their size, density, and shape; [109,110]) and selective entrainment (i.e., the preferential removal of lower density and, therefore, coarser settling-equivalent grains from deposits undergoing strong localized erosion; [111,112]).
As shown by previous studies, geochemical parameters may not be particularly helpful to discriminate sediment sources in provenance analysis ( [113]), especially because of superposed grain-size and hydraulic-sorting control ( [84,114]). This is particularly true for the Ganga-Brahmaputra sediment system, where powerful tractive currents locally promote a full order-of-magnitude increase in the concentration of dense minerals relative to "neutral" composition (i.e., the composition that sediment would ideally have everywhere in the absence of hydrodynamic control).
If, on the one hand, bulk-sediment chemical composition offers little help for provenance discrimination, on the other hand it provides essential information to identify and quantify hydraulic-sorting effects by the precise assessment of elements such as REE, Th, Zr, or Hf. These elements do not enter the crystal lattice of common silicates and are thus progressively concentrated in residual melts until the latest stages of magmatic crystallization, when they form ultra-dense accessory phases such as monazite or zircon. Chemical data also allow us to highlight the differences in weathering conditions found in different parts of the Bengal sediment system.

Ganga and Brahmaputra Rivers
Ganga and Brahmaputra sands transported from the central and eastern Himalaya to the Bengal Sea are feldspatho-litho-quartzose to litho-feldspatho-quartzose metamorphiclastic, with subequal amounts of K-feldspar and plagioclase, low-to-high rank metamorphic lithic grains, biotite, muscovite, amphibole, epidote, and garnet ( Figure 3). This is the typical fingerprint of orogenic detritus produced by focused erosion of deeply exhumed neometamorphic axial belts at the core of high-relief collision orogens (Axial Belt Provenance; [115]).
Ganga sand is richer in quartz and contains limestone and dolostone grains from Himalayan sedimentary and metasedimentary rocks. Heavy minerals include common garnet, tourmaline, diopsidic clinopyroxene, kyanite, and sillimanite, documenting prominent contributions from amphibolite-facies metasedimentary rocks of the Greater Himalaya ( [10]).
Brahmaputra sand is richer in feldspars (particularly plagioclase), lacks carbonate grains, and includes minor mafic and ultramafic detritus from ophiolitic units exposed along the Yarlung-Tidding suture zone. Heavy minerals are dominated by blue-green hornblende and epidote derived from amphibolite-facies Himalayan units exposed around the Namche-Barwa syntaxis as well as from Trans-Himalayan batholiths exposed in the Gangdese Range and Mishmi Hills. Hypersthene and chloritoid derived, respectively, from high-grade and low-grade rocks, also occur ( [9]).
As it is the case for bedload sand, Ganga suspended load contains more quartz and significant calcite and dolomite, whereas Brahmaputra suspended load contains more Ca-plagioclase, very little dolomite, and no calcite. A moderately-rich epidote-amphibole-garnet tHM suite characterizes Ganga suspended load, whereas a rich amphibole-epidote suite characterizes Brahmaputra suspended load ( [13]).

Hydraulic Sorting
Because of selective-entrainment effects, heavy-mineral concentration varies by one order of magnitude and more, from as low as 1.6 to 14.7 for Ganga bedload and from 3.0 to 22.6 for Brahmaputra bedload. Such a strong compositional variability is highlighted by geochemical data, showing that the concentration of zirconium, chiefly hosted in ultra-dense zircon grains, ranges from 116 to~1470 ppm in Ganga sand and from 145 to~2750 ppm in Brahmaputra sand ( [12,107]). Even garnet placers may form in Brahmaputra bedload, reaching a heavy-mineral concentration of 73.5, a grain density of 3.65 g/cm 3 (i.e., much higher than eclogite), and Th and Zr concentrations of~350 ppm and~6500 ppm (i.e.,~33 times higher than the UCC standard; [12]).
Because of suspension sorting, quartz, feldspar, and heavy minerals steadily increase with depth in the river channel relative to slow-settling platy micas and clay. Among heavy minerals, fast-settling zircon and garnet tend to increase with channel depth relatively to less dense amphibole and epidote. Steady mineralogical trends are faithfully reflected by geochemical trends, with higher concentration in chemical elements preferentially hosted in micas or associated with clay and oxy-hydroxides (i.e., Al, K, Rb, Cs, Ba, and heavy metals) close to the surface, and progressive increase in elements preferentially hosted in denser minerals (e.g., Zr, Hf) and more equant tectosilicates (Si, Na, Ca) toward the bed ( Table 2).

Meghna Estuary and Subaerial Bengal Delta
Similar hydrodynamic processes characterize the Padma River and the Meghna estuary downstream, where fluvial and tidal currents locally cause strong erosion of river bars and redeposition in lower-energy estuarine tracts (Figure 6b in [6]). In litho-feldspatho-quartzose sand, tHM suites range from rich to extremely rich and contain abundant amphibole, epidote, and common garnet (12-22%; Table 1). Zr concentration varies from 218 to 2000 ppm, and also Y, REE, Th, Hf, Ta, Mn, and P-all preferentially hosted in dense and ultra-dense minerals (monazite, zircon, rutile, garnet, allanite, apatite)-may be enriched by factors between three and nine ( Figure 4A). In the shallow suspended load of the Meghna estuary, instead, the rich amphibole-epidote tHM suite contains < 3% garnet and only 163-169 Zr ppm.  Figure 2I), sample SO347 displays low concentration in heavy metals, relatively high Zr, and low Al, suggesting effective reworking and winnowing of palimpsest sand on the outer shelf.

Bengal Shelf
The same mineralogical trend observed from bedload to deep and shallow suspended load in fluvial and estuary sediments is reproduced from sandy topsets to silty foresets and bottomsets on the subaqueous delta, with progressive enrichment in platy phyllosilicates in distal settings at the expense of tectosilicates and heavy minerals ( Table 2). Among heavy minerals, slower-settling amphibole tends to increase in fine silt at the expense of fast-settling garnet, which is more common in sand. The notable enrichment in durable ZTR minerals, dearth of pyroxene, and the high epidote/amphibole ratio in sample SO347 hint at early diagenetic dissolution of palimpsest sediment originally deposited by the Ganga River.  Table S4). (A) Fluvial, estuary, and subaerial delta sediments compared to Bengal shelf sand. The Katka beach and some Padma and Meghna estuary samples are markedly enriched in elements preferentially hosted in heavy minerals (e.g., Y, REE, Th, Zr, Ta, Mn, and P), indicating that they are semi-placer lags generated by selective entrainment of lower density grains. Shelf sand has more Cs, W, Cu, and As than fluvial and estuary sands, elements all markedly enriched in tidal mud of the Pussur distributary. (B) Bengal shelf mud compared to Bengal shelf sand. Geochemical signatures are chiefly controlled by grain size and clay content. Besides abundant Ca associated with allochems ( Figure 2I), sample SO347 displays low concentration in heavy metals, relatively high Zr, and low Al, suggesting effective reworking and winnowing of palimpsest sand on the outer shelf.
Subaerial Bengal delta samples range from feldspatho-litho-quartzose to litho-feldspatho-quartzose and tend to be richer in quartz and much poorer in plagioclase than Meghna estuary sand. The tHM suites are only moderately rich, including 13-16% garnet, <3% kyanite, and ≤0.5% staurolite; Zr is <400 ppm. The Katka beach sand, instead, has a very rich tHM suite including 32% garnet, 4% kyanite, and 2.5% staurolite; Zr is 1344 ppm. This sample can thus be considered as a semi-placer lag formed by selective removal of low-density grains by storm waves during beach erosion.

Bengal Shelf
The same mineralogical trend observed from bedload to deep and shallow suspended load in fluvial and estuary sediments is reproduced from sandy topsets to silty foresets and bottomsets on the subaqueous delta, with progressive enrichment in platy phyllosilicates in distal settings at the expense of tectosilicates and heavy minerals ( Table 2). Among heavy minerals, slower-settling amphibole tends to increase in fine silt at the expense of fast-settling garnet, which is more common in sand. The notable enrichment in durable ZTR minerals, dearth of pyroxene, and the high epidote/amphibole ratio in sample SO347 hint at early diagenetic dissolution of palimpsest sediment originally deposited by the Ganga River.
The analyzed Bengal shelf sands display Zr concentrations between 361 and 581 ppm, higher than the UCC values by factors between approximately two and three but not as high as the Padma River, Meghna estuary, and Katka beach semi-placers (Zr 1344-2000; Table 1). This reflects prevailing depositional processes on the subaqueous delta and only mild concentration of fast-settling heavy minerals by winnowing of slow-settling grains such as phyllosilicates rather than by strong erosional processes.
The joint inspection of Figure 4A,B reveals another aspect of hydraulic control on the compositional variability of Bengal Shelf sediments associated with either selective winnowing or settling of clay. The analyzed shelf sands display higher concentrations than fluvial and estuarine bedload in heavy metals and elements such as Cs and As, largely adsorbed on clay or associated with oxy-hydroxides and organic matter (e.g., [116][117][118]). The concentration of these elements is even markedly higher in tidal muds of the Pussur distributary channel, and highest in fine silt of distal foresets and bottomsets of the subaqueous delta.
Besides abundant Ca associated with ooids and bioclasts ( Figure 2I), the outer shelf sample SO347 displays the lowest concentration in heavy metals, relatively high Zr, and low Al ( Figure 4B), suggesting effective winnowing of clay by storm-generated bottom currents. Palimpsest sand including ooids and bioclasts represents drowned beach barriers deposited during the last lowstand stage around 20 ka ( [119,120]), followed by reworking and sediment starvation during the Holocene transgression, as also indicated by growth of glaucony within foraminiferal chambers and lack of reported sites where ooids are forming at present in the Indian Ocean ( [121,122]). Protracted sediment starvation just in front of the largest deltaic sediment system on Earth indicates that throughout the Holocene sea-level rise sediment was systematically swept by tropical cyclones toward the Swatch of No Ground, and from there funneled toward the Bengal Fan ( [7]), thus systematically bypassing the outer shelf.

Chemical Weathering
Besides physical effects such as selective mechanical breakdown, which can be considered as negligible for the purpose of this study ( [123,124]), and hydraulic-sorting, agents capable of compositional modification include chemical weathering. Previous studies have highlighted the difference between the essentially physical Brahmaputra River system-characterized by a cold and arid climate in Tibet, ultra-rapid erosion across the eastern Himalayan syntaxis, and very high runoff in the narrow floodplain-and the Ganga river system, where storage time in the floodplain is longer and weathering notably more intense ( [15,125]). This is reflected by clay-mineral assemblages ( [126][127][128][129][130]), which are dominated by illite (70-80%) with subequal amounts of chlorite and kaolinite (~10% each) and negligible smectite (≤5%) in Brahmaputra suspended load, whereas smectite/montmorillonite is much more common (20-40%) in Ganga suspended load.
The analysis presented here adds to the thorough study by [17], who used mobile to immobile element ratios including K/Si and H 2 O + /Si as well as the abundance of detrital calcite to document an increase of chemical weathering in the Himalayan system since the Last Glacial Maximum. The use of diverse chemical indices allows us not only to clarify their different response to different controls (Table 1) but also to perceive the different qualitative character of weathering processes in the Brahmaputra and Ganga catchments ( Table 2).

The Meaning of Chemical Indices
The broadly homogeneous composition of Meghna estuary and Bengal shelf sediments reflects the continuous supply of enormous volumes of detritus dominantly derived from the Himalayan orogen via the Ganga and Brahmaputra rivers. Because, in this huge system, sediment provenance can be safely considered as substantially unvaried in the short term, modern sediment samples offer a way to test whether the so-called chemical indices of weathering are indeed mainly controlled by climate-related weathering, which at the very short time scale of a few years can be held as roughly invariant, or they are largely or even overwhelmingly controlled by diverse factors including source-rock lithologies, hydraulic sorting, quartz addition, and grain size. Such a critical assessment is needed to prevent that these indices are used indiscriminately as an act of faith.
Close inspection of Table 1 is instructive in this regard. Al/Si, a widely used proxy for grain size (e.g., [8]), correlates quite well with the WIP (r = 0.83; significance level = 0.1%), which confirms that the WIP is strongly dependent on grain size and quartz content and, therefore, can hardly be used as a measure of weathering in sediments ( [98,99,131]). The CIA also correlates with Al/Si (i.e., with grain size; r = 0.64) and best with α Al Ca (r = 0.89), a correlation that disappears if the CIX is calculated instead. The addition of carbonate grains, also formed as allochems within the sedimentary basin (e.g., Figure 2I), may lead to the underestimation of weathering, even markedly if not properly identified and corrected for.
An almost perfect correlation exists between α Sm Na and Zr (r = 0.98; significance level = 0.1%), which demonstrates the overwhelming effect of concomitant hydraulic concentration of ultra-dense minerals such as REE-rich monazite and zircon. Dominant hydraulic-sorting control on trace-element concentration confirms that α indices using as a reference non-mobile elements such as Ti, Sm, Nd or Th ( [101]), which are preferentially hosted in ultra-dense minerals (e.g., monazite, allanite, ilmenite, rutile), cannot be used to evaluate weathering in sediments deposited by tractive currents ( [98,99]). The hydraulic-sorting effect is efficiently sterilized by using, instead, the non-mobile element Al as a reference ( Table 1).
The study of African sediments has shown that α Al Na is the chemical index most sensitive to weathering, largely because of the high emigration rate of Na + from plagioclase under active leaching ( [132][133][134][135]). The use of α Al Na-and in general of α Al E indices instead of α Sm Na, α Ti Ca, α Th K, α Nd Sr, and α Th Ba indices-avoids undesired perturbations up to even more than one order of magnitude associated with the locally strong hydraulic concentration or depletion of heavy minerals preferentially hosting Ti or REE (Table 1), and is thus recommended.
No chemical index, however, can be claimed to be a truthful proxy for weathering if taken alone. The use of the CIX instead of the CIA circumvents the problems associated with correction for CaO not hosted in silicates, although the CIX is more affected than the CIA by selective-entrainment processes and reaches extreme values even in the absence of weathering in placer deposits strongly depleted in low-density alkali feldspars. Being strongly affected by quartz addition, the WIP is best used in conjunction with the CIA or the CIX to evaluate the extent of recycling rather than the intensity of weathering ( [98,131]).

Weathering Processes
The inability of the WIP to measure weathering is clearly exposed by inspection of Table 2, which shows that all indices increase steadily with decreasing grain size from bedload to deep and shallow suspended load, being minimum in sand and maximum in clay in all tracts of the Bengal sediment system ( [8,13,15,108]). Higher values of the WIP, however, indicate lower and not higher weathering intensity. The reason why the WIP is unduly lower for sand than for clay is that the WIP is affected principally by quartz content-which is markedly higher in bedload sand than in phyllosilicate-rich suspended load-rather than by chemical weathering.
All indices but α Al Ca concur to confirm more intense weathering in the Ganga catchment than in the Brahmaputra catchment. The main reasons why α Al Ca and α Al Na give opposite indications for the two rivers do not only reflect lithological differences in source terranes but also the distinct character of weathering in the two river systems. Carbonate rocks provide a much larger share of detritus to the Ganga River, largely via the Gandak tributary that drains Tethys Himalayan and Greater Himalayan carbonate and metacarbonate rocks in the upper course and some Lesser Himalayan metacarbonates in the middle course ( [10]). The Gandak River is estimated to contribute possibly even as much as half of total Ganga load ( [100]). Consequently, despite conspicuous chemical weathering in the floodplain, the Ca/Al ratio is higher in Ganga bedload and deep suspended load than in the average upper continental crust. In Brahmaputra sediments, instead, α Al Ca displays the same values and trend as α Al Na from bedload to deep and shallow suspended load (Table 2) because in this case both elements are chiefly controlled by plagioclase weathering. The Brahmaputra River carries a significant amount of carbonate grains as far as the Himalayan mountain front (~6%, 20% of which calcite and 80% dolomite), reduced to~1% (100% dolomite) within the first 100 km of transport across the foreland basin, and progressively reduced to zero farther downstream ( [9]). On the same samples studied for mineralogy, manometrical measurement of CO 2 after 3 h for calcite and 1 week for dolomite gave comparable results: 2.2% calcite and 2.2% dolomite at the mountain front, 0.2% calcite and 0.4% dolomite after 100 km of floodplain transport, and 0.003-0.2% calcite and 0.004-0.5% dolomite at the entry point in Bangladesh, where calcite is, however, still 0.15% and dolomite 1.2% in suspended load ([107). Near complete dissolution of carbonate grains in the Assam foreland-basin lowlands is a weathering effect ascribed to huge Brahmaputra discharge with high pCO 2 levels in river waters undersaturated with respect to carbonates ( [136,137]).
In the Meghna estuary, weathering indices are either similar to Brahmaputra sediments or intermediate between Brahmaputra and Ganga sediments. In the subaerial Bengal delta, tidal mud of the Pussur distributary results as weathered as Ganga-Brahmaputra suspended load. In the subaqueous delta, sediments are markedly depleted in Ca and Na as grain size decreases and clay-the product of chemical weathering-progressively increases from sandy topsets to muddy distal foresets and bottomsets.

Provenance
In this section, we draw inferences about provenance and sediment mixing along the routing system that includes the Meghna estuary and the Bengal shelf, and terminates in the Bengal-Nicobar Fan.

Meghna Estuary
The Ganga and Brahmaputra rivers drain different geological domains besides the Himalayan thrust belt and, thus, carry sediments with distinguishable mineralogical and isotopic fingerprints ( Figure 3). Ganga sediments, largely derived from Greater and Lesser Himalayan rocks, are, for instance, characterized by more negative ε Nd and higher 87 Sr/ 86 Sr than Brahmaputra sediments, which are partly derived from juvenile Trans-Himalayan magmatic rocks in south Tibet ( [100,107]). The relative supply from the two big rivers to the Padma River, the Meghna estuary, and the Bengal shelf is, however, difficult to calculate precisely because of the strong variability of sediment composition, largely caused by hydraulic-sorting effects.
In the Meghna estuary, detrital modes are closer to Brahmaputra sand than to Ganga sand (Figure 3), and tHM suites are virtually indistinguishable from Brahmaputra tHM suites ( Table 2). The close similarity of tHM suites in the Brahmaputra, Padma River, and Meghna estuary ( Figure 6) is, however, partly an effect of heavy-mineral concentration, which is 50% to 100% higher in Brahmaputra bedload and suspended load than in Ganga bedload and suspended load ( [12,13]). The same holds largely true also for Nd and Sr isotopic ratios, which in the Padma River, Meghna estuary, and Bengal shelf are virtually the same as in Brahmaputra sediments (Figures 7 and 8 in [17,55]). As highlighted in [138], most Nd and much Sr are hosted in heavy minerals, the abundance of which is markedly higher in Brahmaputra than in Ganga sediments, as discussed further in the following subsection.
The mineralogical and geochemical signatures of suspended load are intermediate between the Ganga and Brahmaputra rivers, but notably closer to the Brahmaputra (e.g., relatively low quartz/feldspar ratio, heavy-mineral assemblages dominated by hornblende and epidote with minor garnet, diopside, sillimanite, and kyanite). A few carbonate grains occur in Padma bedload, and suspended load has more CaO than in the Brahmaputra river at any water depth but more so in shallow suspended load, suggesting that contributions from the Ganga River, although invariably subordinate, is more significant for finer-grained fractions ( Figure 5). Greater Ganga contribution for cohesive mud fractions is also indicated by the anomalous increase in the Q/F ratio with decreasing grain size of Bengal shelf muds ( [19]).
Minerals 2019, 9, x FOR PEER REVIEW 20 of 29 Figure 6. Geochemical discrimination between Ganga and Brahmaputra provenance. Brahmaputra sediments contain more plagioclase and epidote and, therefore, more Na (A) and Sr (B). Subaerialdelta sediments are akin to the Ganga, whereas subaqueous-delta sediments are akin to the Meghna estuary and to the Brahmaputra for sand but closer to the Ganga for fine silt. Sample SO347 is as low in Na as Ganga sand, Al poor because of the winnowing of clay, and Sr rich because of common calcareous ooids and bioclasts. The size of symbols is roughly proportional to the grain size for subaerial and subaqueous delta samples.

From the Bengal Shelf to the Bengal-Nicobar Fan
Neogene Bengal Fan turbidites display litho-feldspatho-quartzose metamorphiclastic composition with hornblende-dominated heavy-mineral assemblages ( [20,21]), comparing closely with modern Brahmaputra and Meghna estuary sediments and, thus, indicating continuous supply from the Ganga-Brahmaputra sediment system and dominant Himalayan provenance since the Miocene ( [4,142,143]). The mineralogy of upper Quaternary silty turbidites cored at ODP Sites 717-719 ( [22]) is very close to that of Meghna estuary sediments (Figure 5), documenting a remarkable compositional homogeneity of the fluvio-deltaic to deep-sea-fan depositional system over ~3000 km, from Bangladesh to the equatorial Indian Ocean.
Additional detritus from the Indian subcontinent is, however, indicated in clay-rich intervals deposited during periods when turbidite channels shifted eastward ( [22]). This is confirmed by geochemical data on uppermost Quaternary sediments of the western Bay of Bengal ( [144]) and by the occurrence of olivine-which is rare in Meghna estuary sediments-in upper Pleistocene distalfan turbidites, suggesting supply from Deccan Trap basalts via the Godavari and Krishna rivers ( [145]).
Nicobar Fan turbidites also display feldspatho-quartzose to litho-feldspatho-quartzose metamorphiclastic detrital modes and mostly moderately rich amphibole-epidote-garnet tHM suites that compare closely with those of Meghna estuary and Bengal shelf sediments (Figure 3), confirming that they also belong to the colossal Bengal sediment system ( [23]). The abundance of feldspar and the very close similarity of amphibole-epidote tHM suites of Bengal-Nicobar Fan turbidites with Brahmaputra River and Meghna estuary bedload ( Figure 5) indicates a clear prevalence of Brahmaputra contribution lasting since the late Miocene and, thus, since the onset of rapid unroofing of the eastern Himalayan syntaxis ( [146,147]). Geochemical discrimination between Ganga and Brahmaputra provenance. Brahmaputra sediments contain more plagioclase and epidote and, therefore, more Na (A) and Sr (B). Subaerial-delta sediments are akin to the Ganga, whereas subaqueous-delta sediments are akin to the Meghna estuary and to the Brahmaputra for sand but closer to the Ganga for fine silt. Sample SO347 is as low in Na as Ganga sand, Al poor because of the winnowing of clay, and Sr rich because of common calcareous ooids and bioclasts. The size of symbols is roughly proportional to the grain size for subaerial and subaqueous delta samples.
Forward mixing calculations based on integrated petrographic and heavy-mineral data ( [139,140]) indicate that sand supply from the Brahmaputra river is at least twice that of the Ganga River, despite the drainage basin of the Ganga being notably greater than that of the Brahmaputra. Similar calculations based on clay-mineral assemblages (as reported in Table 1 in [130]) indicate that the Ganga River contributes ≤ 40% and the Brahmaputra River~60% of clay minerals fed into the Bay of Bengal ( [19]). Because the drainage area of the Brahmaputra is~60% that of the Ganga, such a broad agreement among estimates based on different compositional parameters and grain-size fractions would imply that average denudation rates are about three times higher in the Brahmaputra catchment than in the Ganga catchment. However, a larger part of the Ganga sediment flux is sequestered in the foreland basin and part of it is diverted southward before reaching the Brahmaputra confluence.
Tracing of Sr concentration throughout the Holocene succession of the Bengal delta indicates that the western part of the subaerial delta consists almost entirely of low-Sr Ganga sediments, whereas high-Sr Brahmaputra sediments have remained dominant in the eastern region including the modern estuary ( [55]). This is consistent with petrographic and heavy-mineral modes of sand cored in the western part of the subaerial Bengal delta ( Figure 2D) and collected at Katka beach ( Figure 2E), which are closer to Ganga sand than to Brahmaputra sand (Figure 3), and with the tHM suite of tidal mud in the Pussur distributary ( Figure 2F), which is indistinguishable from deep suspended load in the Ganga River ( Figure 6). that the western part of the subaerial delta consists almost entirely of low-Sr Ganga sediments, whereas high-Sr Brahmaputra sediments have remained dominant in the eastern region including the modern estuary ( [55]). This is consistent with petrographic and heavy-mineral modes of sand cored in the western part of the subaerial Bengal delta ( Figure 2D) and collected at Katka beach ( Figure 2E), which are closer to Ganga sand than to Brahmaputra sand (Figure 3), and with the tHM suite of tidal mud in the Pussur distributary ( Figure 2F), which is indistinguishable from deep suspended load in the Ganga River ( Figure 5).

Figure 5.
Multidimensional scaling map based on heavy-mineral data (plotted using the "provenance" package of [141]). The distance among samples is approximately proportional to the Kolmogorov-Smirnov dissimilarity of their tHM suites. Solid and dashed lines link closest and second-closest neighbors, respectively; the "stress" value of the configuration is 5.4, indicating a "good" fit ( [78]). Samples are separated chiefly by grain size along the y-dimension and by provenance along the x-dimension, which highlights the Ganga affinity of subaerial Bengal delta Figure 6. Multidimensional scaling map based on heavy-mineral data (plotted using the "provenance" package of [141]). The distance among samples is approximately proportional to the Kolmogorov-Smirnov dissimilarity of their tHM suites. Solid and dashed lines link closest and second-closest neighbors, respectively; the "stress" value of the configuration is 5.4, indicating a "good" fit ( [78]). Samples are separated chiefly by grain size along the y-dimension and by provenance along the x-dimension, which highlights the Ganga affinity of subaerial Bengal delta sediments, and the Brahmaputra affinity of Meghna estuary, Bengal shelf, and Bengal Fan sediments. As in Figure 3D, Brahmaputra affinity of estuary, shelf, and fan sediments partly reflects the higher heavy-mineral concentration in Brahmaputra than in Ganga sediments.

Minerals Controlling the Sr and Nd Budgets
In shallow to deep suspended load, heavy minerals are estimated to contribute 17% to 21% of total Sr (13% to 16% from epidote) for the Ganga, and 25% (~20% from epidote) for the Brahmaputra. Heavy minerals are estimated to contribute 53% to 76% of total Nd (43-46% to 60-65% from allanite + monazite + titanite) in both Ganga and Brahmaputra shallow to deep suspended load. In bedload, heavy minerals may contribute up to 28% of total Sr (20% from epidote) and up to 96% of total Nd (88% from allanite + monazite + titanite) for the Ganga, and even up to 84% of total Sr (66% from epidote) and almost all of Nd (93% from allanite + monazite + titanite) for the Brahmaputra ( [12,13]).
Most of the remaining Sr is hosted in Ca-plagioclase, which is at least twice as abundant in Brahmaputra than in Ganga bedload and suspended load, which explains why Brahmaputra sediments contain 50-60% more Sr than Ganga sediments in any size fraction ( Figure 5). In summary, the Sr and Nd budgets are chiefly dependent on Ca-plagioclase and heavy minerals (between 56% and 70% for the Ganga and between 70% and 96% for the Brahmaputra), which in any size fraction are more abundant in Brahmaputra than in Ganga sediments by factors between 1.5 and ≥2 (Table 2).

Bengal Shelf Sediments
The petrographic and chemical composition of shelf sand is broadly intermediate between Brahmaputra and Ganga bedload ( Figure 5), indicating that both rivers contribute large amounts of sediment. Partly because of the reasons discussed above, the tHM suite of shelf sediments ( Figure 6) and their Sr and Nd isotopic fingerprints (Figures 7 and 8 in [17]) are instead barely distinguishable from those of Brahmaputra and Meghna estuary sediments.
Silt deposited on the Bengal shelf has quartz/feldspar ratio and proportions among feldspar minerals very close to Brahmaputra silt, and a tHM suite virtually indistinguishable from suspended load in both the Brahmaputra River and Meghna estuary ( Figure 6). Fine cohesive mud deposited on the distal foresets and bottomsets of the subaqueous delta may be interpreted as suspended load largely delivered by the Meghna estuary mixed with clay largely supplied by tidal distributaries draining the Sundarban swamps ( [19]). Illite-dominated clay-mineral assemblages include subordinate chlorite, smectite, and kaolinite in very similar proportions as clay minerals in Meghna estuary mud ( [19,130]).
These estimates confirm that sediment of all size fractions deposited on the subaqueous Bengal delta is mostly derived from the Meghna estuary-with possibly significant contributions of mud delivered by tidal distributaries draining the western part of the inactive subaerial delta-and ultimately supplied between 30% and 40% by the Ganga River and between 60% and 70% by the Brahmaputra River. Sediment flux across the Bay of Bengal and subsequently funneled via the Swatch of No Ground to the Bengal Fan is thus presently Brahmaputra-dominated ( Figure 7).  . Sediment in transit across the Bengal system ( [5][6][7]). Note the predominance of Brahmaputra (magenta) over Ganga (cyan) supply to the Meghna estuary and subaqueous delta. The subaerial delta, built by the Ganga River and progressively abandoned during the late Holocene highstand is now undergoing tidal reworking and locally net erosion, thus contributing mostly mud originally deposited by the Ganga River to the western part of the subaqueous delta. The studied sand and sandy silt samples are color-coded according to grain size as in Figure 1.
A robust quantitative comparison among mineralogical assemblages, however, can be drawn only for sediments buried less than a few hundreds of meters, because the tHM suites of ancient strata are exposed to selective intrastratal dissolution of less durable minerals progressing with increasing age and burial depth ( [148][149][150]). Heavy-mineral studies of Bengal Fan turbidites have documented Figure 7. Sediment in transit across the Bengal system ( [5][6][7]). Note the predominance of Brahmaputra (magenta) over Ganga (cyan) supply to the Meghna estuary and subaqueous delta. The subaerial delta, built by the Ganga River and progressively abandoned during the late Holocene highstand is now undergoing tidal reworking and locally net erosion, thus contributing mostly mud originally deposited by the Ganga River to the western part of the subaqueous delta. The studied sand and sandy silt samples are color-coded according to grain size as in Figure 1.

From the Bengal Shelf to the Bengal-Nicobar Fan
Neogene Bengal Fan turbidites display litho-feldspatho-quartzose metamorphiclastic composition with hornblende-dominated heavy-mineral assemblages ( [20,21]), comparing closely with modern Brahmaputra and Meghna estuary sediments and, thus, indicating continuous supply from the Ganga-Brahmaputra sediment system and dominant Himalayan provenance since the Miocene ( [4,142,143]). The mineralogy of upper Quaternary silty turbidites cored at ODP Sites 717-719 ( [22]) is very close to that of Meghna estuary sediments (Figure 6), documenting a remarkable compositional homogeneity of the fluvio-deltaic to deep-sea-fan depositional system over~3000 km, from Bangladesh to the equatorial Indian Ocean.
Additional detritus from the Indian subcontinent is, however, indicated in clay-rich intervals deposited during periods when turbidite channels shifted eastward ( [22]). This is confirmed by geochemical data on uppermost Quaternary sediments of the western Bay of Bengal ( [144]) and by the occurrence of olivine-which is rare in Meghna estuary sediments-in upper Pleistocene distal-fan turbidites, suggesting supply from Deccan Trap basalts via the Godavari and Krishna rivers ( [145]).
Nicobar Fan turbidites also display feldspatho-quartzose to litho-feldspatho-quartzose metamorphiclastic detrital modes and mostly moderately rich amphibole-epidote-garnet tHM suites that compare closely with those of Meghna estuary and Bengal shelf sediments (Figure 3), confirming that they also belong to the colossal Bengal sediment system ( [23]). The abundance of feldspar and the very close similarity of amphibole-epidote tHM suites of Bengal-Nicobar Fan turbidites with Brahmaputra River and Meghna estuary bedload ( Figure 6) indicates a clear prevalence of Brahmaputra contribution lasting since the late Miocene and, thus, since the onset of rapid unroofing of the eastern Himalayan syntaxis ( [146,147]).
A robust quantitative comparison among mineralogical assemblages, however, can be drawn only for sediments buried less than a few hundreds of meters, because the tHM suites of ancient strata are exposed to selective intrastratal dissolution of less durable minerals progressing with increasing age and burial depth ( [148][149][150]). Heavy-mineral studies of Bengal Fan turbidites have documented the common occurrence of unstable pyroxene and olivine only in the upper part of the studied cores (i.e., Zone I of [22], dated as~0.5 Ma, burial depth ≤ 160 m b.s.f.). Ferromagnesian minerals including amphibole progressively decrease in the Pliocene and older sediments buried more than 250-300 m at DSDP Sites 211 and 218 ( [20,151]), and the relative abundance of more durable zircon, tourmaline, rutile, and garnet consequently increases (correlation coefficients with core depth 0.54 to 0.72, all significant at the 1% level). Amphibole prevails over epidote + garnet in Quaternary strata (Zones I and IIa, burial depth ≤ 200 m b.s.f.), is less than both epidote and garnet in middle Miocene strata (Zone V, burial depth 650-750 m b.s.f.), to eventually become negligible in lower Miocene strata (Zone VI, burial depth > 750 m b.s.f.; [22]).

Conclusions
This study adds new data and observations that, combined with extensive previous work, contributes to a better understanding and quantification of the enormous sediment flux that throughout the Neogene has transited across the Bengal basin and accumulated on Bengal Sea floors. Detritus mostly generated from rapid erosion of the active Himalayan orogen and entrained by the Ganga and Brahmaputra rivers, which joined two centuries ago to form the Meghna estuary, is partly stored in the subaqueous delta prograding onto the Bengal shelf, as in the subaerial delta formed by the Ganga River to the west and progressively abandoned during the late Holocene highstand.
From the Meghna estuary and the Bengal shelf all the way to the Bengal and Nicobar deep-sea fans, detritus is closer in composition to Brahmaputra sediments than to Ganga sediments, as testified by all mineralogical and geochemical parameters. Transparent heavy-mineral suites in sand and silt of the Meghna estuary, Bengal shelf, and Bengal-Nicobar Fans, as well as Sr and Nd isotope fingerprints ( [17]), are virtually indistinguishable from those of Brahmaputra bedload and suspended load. Such a predominance over the Ganga is overemphasized because the Brahmaputra carries between 50% and ≥100% more Ca-plagioclase, epidote and heavy minerals in general, which together control most of the Sr and Nd budget in the sediment ( [138]). A more balanced contribution from the two rivers is suggested by the most commonly cited estimates of their suspended loads as well as by bulk-sediment mineralogy and geochemistry, which converge to indicate that the Ganga may supply as much as 40% of the total sediment flux. We conclude that, although the Brahmaputra River supplies no less than two-thirds of the sand, Ganga contribution may be more significant for finer silt fractions and clay. This is because the Ganga River, having lower flow velocity and competence, entrains finer particles than the Brahmaputra at all channel depths ( [8,14]), and also because the subaerial delta actively formed by the Ganga until a few ka ago is presently undergoing tidal reworking and erosion ( [6,66]), thus contributing fine material to the western part of the active subaqueous delta ( [19]). Although denser minerals tend to be segregated in coastal and shallow-marine environments, and slow-settling platy micas and clay are effectively winnowed and accumulated offshore, the studied shelf sediments never display the strong hydraulic-sorting effects observed in fluvial, estuarine, and beach sediments, reflecting the prevalence of depositional processes in the prograding subaqueous delta.
Although the Ganga and Brahmaputra drain also different geological domains, the Indian shield and the Trans-Himalayan batholiths and suture zone, respectively, the two rivers receive a large amount of detritus from the same Himalayan rock units. As a consequence, their sediments have mineralogical and, therefore, geochemical signatures that, although distinct, are not different enough to allow very precise estimates of their relative detrital supply. The task is made particularly difficult by the superposed effects of physical and chemical processes, which can produce a marked distortion of provenance signals. Hydraulic-sorting bias tends to efface any potential for distinction based on the geochemistry of bedload. Together with clay-mineral assemblages, geochemical data instead provide precious information on weathering conditions, highlighting not only the different intensity but also the different processes of weathering in Ganga and Brahmaputra lowlands. Ganga sediments undergo stronger alteration during more prolonged storage in the wider foreland basin. The Brahmaputra, instead, is chiefly a physical system, and yet the huge runoff leads to virtually complete dissolution of carbonate grains in bedload, whereas very little calcite and some dolomite is preserved in suspended load. Chemical dissolution becomes a much harder obstacle in the quantitative provenance assessment of ancient sediments, because unstable ferromagnesian minerals are selectively dissolved during burial diagenesis and the original mineralogical and geochemical fingerprints are thus progressively blurred in older strata. Nonetheless, the experience obtained from modern settings remains a fundamental aid to investigate the relative role played by climate and rainfall versus tectonically sustained relief in controlling the intensity of erosive processes across an active orogenic belt such as the Himalaya.