Vertical Distribution of Phosphorous Fractions and Bioavailability of the Nutrient in the Southern Indian Ocean

The study evaluated concentration, vertical distribution, and bioavailability of phosphorus fractions in the Indian Ocean core sediments. Sediment samples were collected from two sites in the Indian Ocean during January–February 2010 using gravity corer. Phosphorus fractions in sediments were analyzed by a modified sequential extraction procedure (SEDEX). The vertical profile revealed silt and clay as dominant grain size fractions in both cores. The enriched levels of organic carbon were observed in samples due to adsorption on fine grained sediments. Loosely bound P was the most dominant fraction, and its higher concentration in sediments was due to increased productivity and subsequent inputs of biogenic components. Detritus from marine organisms and lower oxygen concentrations (depleted Eh values) enhanced the production of authigenic P. The lower iron bound P fraction at specific depths was attributed to the nature of Fe minerals present within sediments. Decreased level of organic phosphate fraction was due to mineralization. The presence of Detrital P, the diagenetically stable form of P, might be due to inputs from the Himalayan, peninsular Indian, and Sri Lankan regions. The biogenic form was the dominant phosphate fraction throughout the cores suggesting the dissolution of fish debris. Estimated bioavailable phosphorous (BAP) ranged from 82.2 to 98.7% in Core I and from 90.21 to 92.98% in Core II, suggesting the fact that a major portion of the phosphate was bioavailable and hence the alterations in environmental conditions can release the BAP to water column resulting in nutrient enrichment. Pioneer study related to nutrient dynamics in the study area. The contribution of phosphorus from the continental region was established. Enhanced bioavailable P can lead to a eutrophic state under changing environmental conditions. Pioneer study related to nutrient dynamics in the study area. The contribution of phosphorus from the continental region was established. Enhanced bioavailable P can lead to a eutrophic state under changing environmental conditions.


Introduction
Phosphorus (P) serves as a critical factor in regulating biocommunity structure and functions and significantly influences the biogeochemical cycles of other elements in the environment. It is an essential component of several biological molecules (such as ATP, DNA, RNA, protein, enzymes, etc.) and a building block of the cells, cellular elements, bones, teeth, etc. (Tyrrell 1999;Ruttenberg 2003;Reddy et al. 2005;Serna and Bergwitz 2020). Phosphorus occurs primarily as phosphate in deposits of apatite (Ca 5 F (PO 4 ) 3 ), a mineral found in igneous, sedimentary, and metamorphic rocks; however, approximately 300 additional minerals containing phosphate (PO 4 3− ) have been described by Jahnke (1992). Oceanic sediment acts as the sink of P (Leote et al. 2014;Abu Hmeidan et al. 2018;Jiang et al. 2018;Dan et al. 2020;Maslukah et al. 2021) as well as a source when dissolved P is low (Jiang et al. 2018;Gurung et al. 2020). Phosphorus is usually released into oceanic and marine waters owing to rock weathering and continental erosion. Phosphorus is discharged into aquatic environments in either dissolved or suspended form by the rivers. Other causes that contribute to increasing phosphorus loadings in the environment include the deposition of suspended particles in the atmosphere, volcanic ash, and inorganic dust (Delaney 1998;Benitez-Nelson 2000;Bagalwa 2021).
Phosphorus is present in the sediments either in bioavailable or non-bioavailable form. "Bioavailable phosphate refers to the sum of immediately available phosphate and phosphate, which can be transformed into a functional form by alterations in environmental conditions and biogeochemical processes (Anschutz et al. 2007; Andrieux-Loyer et al. 2008;Wang et al. 2009;Kerr et al. 2010;Cong et al. 2014). The bioavailability of P in aquatic environment influences the primary production rate, species distribution, and ecosystem structure. Monitoring of P content in water and sediment is essential to control and avoid eutrophication of the coastal aquatic environment. Most of the P in the current Page 3 of 15 77 natural environment is present in particulate form and is not biologically available, thereby limiting primary production (Yi et al. 2019).
Different forms of P in the sediments can provide valuable information on the origin of P, the degree of pollution from anthropogenic activities, the bioavailability, and also the burial and digenesis of P in sediments (Andrieux and Aminot 1997;Schenau et al. 2000). The sequential extraction method (SEDEX), described in Ruttenberg (1992), has been recognized as a suitable methodology for fractionating different forms of P from sediments. Furthermore, the SEDEX method has high analytical sensitivity and can measure concentrations as low as 0.0005 wt % P, hence suitable for studying marine sediments. This method is the best option for biogeochemical research as it can separate authigenic carbonate fluorapatite from fluorapatite, which is an important distinction as authigenic CFAP represents an oceanic sink for reactive P and detrital FAP does not (Wang et al. 2013). The method separates the sedimentary phosphorus into P adsorbed onto grain surfaces (Ads-P); P associated with easily reducible iron and manganese oxides and/or oxyhydroxides (Fe-P); authigenic P (this includes authigenic carbonate fluorapatite, biogenic apatite, and CaCO 3 -associated P); detrital P (Det-P); and organic P (OP). Different P phases have different geochemical behaviors, and only certain forms of P (Ads-P, Fe-P, and OP) are bioavailable and released into the water column under various physicochemical processes, such as dissolution, desorption, and reduction. This process depends on the predominant form of P deposited in the sediment layer and its influence on physicochemical characteristics (Ruttenberg 1992;Samadi-Maybodi et al. 2013). However, the method the extraction process is very long and there are chances for partial release of organic phosphorus by CDB reagent (Wang et al. 2013).
The southern part of the Indian Ocean exhibit higher levels of nutrient and other chemical species, primarily due to Indonesian through Flow (ITF) inputs. ITF links the Pacific and Indian Oceans and offers a setup for modifying the stratification within these oceans and sea-air fluxes that impact climatic phenomena such as ENSO and the Asian Monsoon (Potemr 2005). Moreover, the southern part of the Indian Ocean is least investigated regarding nutrient bioavailability/ dynamics. It is crucial to gather information on P fraction composition to understand the bioavailability and dynamics of phosphorus in the study area and its influence on marine productivity. This investigation is intended to evaluate the fraction composition, vertical distribution, and bioavailability of phosphorous in core sediments collected from two distinct locations in the Southern Indian Ocean.

Study Area
The Indian Ocean has been the least studied compared to the Atlantic and the Pacific Ocean in terms of deep-sea benthic geochemistry. The present investigation was carried out in the Central Indian Ocean Basin (CIOB) under the program "Equatorial Indian Ocean Process-Study Dynamics and Biogeochemistry (EIOPS). Two sediment core samples were collected from this ITF zone. The study area and sampling locations are depicted below (Fig. 1). The second sampling site was a rocky area with manganese nodules. The sampling was very tough at this site, and the length of obtained core was petite with few manganese nodules at the top.

Sampling and Analysis
Sampling was carried out on board from January to February 2010, under the flagship of ORV Sagar Kanya. Sediment samples were collected using a gravity corer (10 cm internal diameter and 10 m long) from the above locations. Total 57 samples for Core I and three samples for Core II were obtained. Samples were kept in a deep freezer at − 20 °C prior to analyses. For chemical fractionation of phosphorus, sediments were firstly sieved through a mesh-size screen of 63 μm, and then 0.5 g freeze dried sediment of each sample in 50 ml extractant was analyzed.
The wet sieving determined the percentage composition of sand, silt, and clay (Krumbein and Pettijohn 1938). Total organic carbon content in sediments was estimated using CHNS Analyzer (Vario EL III CHNS Analyzer) after removing inorganic carbon with 10% HCl (Bouillon et al. 2004). The redox potential of the sediment was measured by the potentiometric method (Radojevic et al.1999). The sedimentary phosphate fractions were separated by the five-step process of the SEDEX method initially devised by Ruttenberg (1992) and later modified by Anderson and Delaney (2000). The fractionation scheme for phosphorous was based on sequential extraction followed by spectrophotometry (Babu and Nath 2005). In Fig. 2, the sequential extraction procedure is displayed in detail, by which phosphorus is classified as different forms of loosely-P, ironbound P, authigenic-P, detrital-P, and organic-P. In addition, the sediment sample used in a stage was then applied for the next stage. It is worth noting that the extracted solution was filtered through a 0.45 μm GF/C filter membrane before phosphorus measurement, and the pH of the extract was then reached to neutral point (Zhang et al. 2010;Lin et al.2013). In each phase and after extraction, samples were centrifuged at 4000 rpm for 20 min, and phosphorus content was measured according to the method of molybdenum blue/ ascorbic acid (Murphy and Riley 1962) at 880 nm wavelength using a UV-Vis spectrophotometer (Genesys 10 UV Thermo Spectronic). Potassium dihydrogen orthophosphate was used as the calibration standard for phosphate estimation. All the containers involved in phosphate determination were cleaned with a nutrient free detergent, rinsed with ultrapure water, soaked in 10% HCl overnight, and rinsed again with ultrapure water (Worsfold et al. 2005). Containers used in the subsampling of sediments (after slicing the core) were twice with the water of interest prior to sample collection. Analytical blanks were employed carefully to nullify the contribution of phosphate from reagents and water involved in the estimation. Analysis was carried out in triplicates to ensure the reproducibility of the data. The vertical distribution graph was plotted using Origin Pro 22. Pearson correlation and Principal Component Analysis were carried out using Origin Pro 22 to assess the relationship between different forms of phosphorus with sediment grain size and total organic carbon (TOC).

Results and Discussion
Distribution of sand, silt, and clay in core sediments is shown in Tables 1 and 2. Texture analysis revealed the dominance of finer particles (silt + clay) throughout the Core. The contents of sand, silt, and clay in Core I and Core II in relation to depth are shown in Figs. 3 and 4, respectively. Sand content varied from 0 to 1.48%, with an average of 0.08%. The maximum content for sand was noticed at a depth range of 90-100 cm. Downcore variation of sand followed a zigzag pattern from the surface to the bottom till it reached  to 100-110 cm. Steady state levels of clay were observed at a depth range of 130-220 cm. The highest and lowest clay percentage was recorded in the depth range of 500-510 cm (77.63%) and 120-130 cm (5.94%), respectively. In the case of Core II, an increasing trend for clay level was observed, and it ranged from 5.84 to 27.01% (average: 19.80 ± 0.93%). From the surface layer till 420 cm of the Core I exhibited a higher silty horizon with a prevalence of fine sediment fraction with both silt and clay silt ranging from 35.91 to 93.96%, whereas the bottom horizon (430-565 cm) was dominated by clay. Core II exhibited dominance of silt. The top horizon (0-10 cm) was enriched with silt (more than 90%), and the subsurface (20-30 cm) was of clayey silt. The Core II exhibited a considerably higher level of sand content compared to Core I, with a variation of 1.89-7.15%. Silt fraction varied from 66.35 to 91.77% and clay 5.84-27.08%. Fine sediments (silt and clay) were the dominant fractions in both the cores in general, which indicated the influence of hydrodynamic conditions. The unique depositional environment helps to trap organic matter. The sediments act as an excellent nutrition reservoir because they are enriched with organic matter. The high nitrogen and phosphorus content associated with organic matter is made available to the overlying water column through various physico-chemical processes and is utilized by aquatic plants for their growth. Moreover, organic molecules with their ability to make chelate compounds play a significant role in the retention of nutrients in sediments. Many workers have observed that hypoxic and anoxic conditions facilitate the release of phosphorus from sediments to the water column (Mort et al. 2010). It has been reported that lowering redox potential (E h ) increases the solubility of phosphorus (Ann et al. 2000). Anoxic sediment with low E h facilitates the reduction of insoluble ferric oxy-hydroxide compounds to soluble ferrous oxy-hydroxide compounds, increasing the sorption sites and thereby resulting in adsorption of a large amount of phosphorus. However, due to low bonding energy, desorption potential is substantially high. In contrast to aerobic sediments, adsorption is relatively less, but retention of phosphorus is more due to high bonding energy (Barik et al.2019).
Information on (TOC) content in sediments is vital to assess the role played by the organic fraction of sediments in the transport, deposition, and retention of metals and nutrients. The distribution of organic carbon in both the core samples is presented above in Tables 1 and 2, and below its variations along the length is shown in Fig. 5 for each Core. In Core I, the maximum TOC content was 0.89 %, and that of the minimum was 0.12% at depth ranges 10-20 cm and 290-300 cm, respectively. In core II, the level of TOC ranged from 0.15 to 0.12% at depth ranges 10-20 cm and 20-30 cm, respectively. Generally, the TOC content was found in sediments with high clay and silt content and lower E h values. Indeed, organic carbon, buried and conserved in marine sediments, depends on several parameters such as the type of organic matter, terrestrial or marine, as terrestrial is more refractory. Higher TOC content in fine-textured sediments compared to coarse-textured can be ascribed to the fact that differences in C input, rather than long-term decomposition dynamics, since fine-textured sediments tend to be more fertile than coarse-textured due to likely differences in water storage capacity. In an aquatic environment, the content of TOC increases as the texture becomes finer. Clay particles are believed to protect some of the more easily decomposable organic compounds from rapid microbial breakdown through encrustation and entrapment (Rakesh et al. 2020).

Phosphorous Fractions
Sequential extraction provides five phosphorus fractions viz., (1) loosely bound phosphorus (Lo-P) (L1), iron-bound phosphorus (Fe-P) (L2), Authigenic P (P auth ) (L3), Detrital (Det-P) (L4) and organic phosphorus (OP) (L5). The values of each phase of P obtained at corresponding depths of each Core are shown in Tables 3 and 4. The fractions of P show a wide range of variations, which are controlled by a number of factors such as the texture of sediments, organic matter content, the intensity of mineralization of organic matter in the sediment, and redox conditions in the sediments.
Loosely bound P Changes in physical, chemical, and biological conditions can induce the release and transformation of phosphorus in the sediments. Variation in the concentration of phosphorus fraction depends on the sources and biogeochemical processes. It may also depend on the rate of sedimentation, post sedimentation changes, tectonics, and the physical and chemical conditions of the depositional basin (Nath et al. 2005). Loosely-P, exchangeable P, or biogenic P (L1) are necessary inorganic phosphorus forms that are adsorbed by sediments, readily released, then enter into the water column above the surface sediment (Yang et al. 2016). From the surface of the Core, there is a distinct change in the trend of Lo-P, which gradually increased to the bottom in the Core I. The bottom layers after 30 cm exhibited relatively higher values for Lo-P (Fig. 6) in Core I. Nevertheless, in Core II, an overall decrease from top to bottom was noticed. The higher concentration may be due to the enhanced input of biogenic components into the sediments after deposition. Schenau et al. (2000) recorded a more remarkable preservation of fish debris at shallower depths (diagenetic process). This might be the reason for higher Lo-P in the deeper layers and other biogenic proxies (sedimentary calcium carbonate). Higher Lo-P may also reflect increased productivity in the area. A reverse trend was noticed in Core II, with higher productivity in the upper layers. During the Holocene and late Pleistocene, which marks the termination of the last glacial age with warmer climates, are more suitable for increased productivity. Higher Lo-P at the top few layers and decreased values at the bottom reflect variations in productivity. An enhanced rate of productivity is usually associated with higher P concentrations. The hydroxyapatite (Ca 10 (PO 4 ) 6 (OH) 2 ) 8 crystals found in the hard portions of marine fishes are more soluble than fluorapatite. Since the seawater is under saturated with respect to biogenic apatite (Atlas and Pytkowicz 1977), dissolution of fish debris takes place in the upper layer of the sediment (aqueous phase). Exchangeable or loosely sorbed P, carbonate related P, and biogenic P are all included in Leach 1. Figure 6 depicts the variation of this leach with depth. Lo-P in Core I was found to range between 2319.09 and 7628.11 ppm (average: 5189.86 ± 1.67 ppm). A steep downward trend was recorded up to 50 cm, followed by a rise up to 100 cm, a zigzag pattern until 500 cm, and finally an increase at 550 cm. The percentage contribution of Lo-P to TP in Core I (Table 5) ranged from 73.4 to 97.4%, with an average of 88.0%. In the case of Core II, the concentration of L1 was noted to be decreasing. The contribution of Lo-P to TP was from 79.47 to 87.17% (average: 84.29%) in Core II (Table 6).
Easily reducible or reactive iron bound P Variation in concentration of L2 is depicted in Fig. 7. Iron bound P in Core II ranged from 9.32 (0.2% of TP) to 130.23 ppm (4.1% of TP) (average: 64.189 ± 2.32 ppm). However, in the case of Core II, concentration decreases from surface to bottom. The percentage contribution of iron bound P to total P in Core II ranged from 2.24 to 9.21%, with an average of 5.90%. In an alkaline environment, Fe-P can be exchanged for OH − and other inorganic P molecules that are soluble. Since this kind of phosphorus promotes phytoplankton  growth and development, it may be used to assess their growth (Sundby et al. 1992;Slomp et al. 1996). The environmental conditions influence the concentration of Fe-P. The release of Fe-P into the aqueous phase reduces the pH, and the nutrient fraction is used up by primary producers (Andrieux-Loyer and Aminot 2001;Gurung et al. 2020). P associated with iron oxides has shown a gradual decrease for Core II and an increasing trend in the case of Core I. Initially, the phosphorus is dissolved, and O or OH − can adsorb it to become a Fe hydroxide (P-O-Fe) or a Feoxyhydroxide (P-OH-Fe) form (Mortimer 1941). The low Fe-P at specific depths may therefore be due to the nature of Fe minerals present within sediments (Slomp et al. 1996;Tamburini et al. 2002). In addition to this, under oxidizing conditions, ferric and manganic oxides and hydroxides are important adsorption sites for P (Krom and Berner 1980;Slomp et al. 1996). Ferric and manganic phosphate minerals, such as strengite (FePO 4 .2H 2 O) and trivalent Mn phosphate  (MnPO 3 .1.5H 2 O), can form and survive oxidizing conditions. These minerals exist under reduced conditions, resulting in the increased dissolution and release of soluble P into the water (Patrick et al. 1973;Emerson 1976;Emerson and Widmer 1978;Krom and Berner 1980;Moore and Coale 2009). The Core II was taken from the manganese nodule area and also recorded a higher value for Fe-P and TP due to the above reasons. Adsorption of Fe also depends on DO content, and it is seen that there is reduced P sorption at low oxygen conditions. Desorption of P is also greater in acidic conditions. Phosphate buried in sediments will be remobilized later by bacterial degradation. Phosphate present in the pore waters can either be formed by the reductive dissolution of oxides of Fe or by the phosphate released by the degeneration of organic matter (Filipelli and Delaney 1996). However, oxygen depletion and redox potential decrease can stimulate P release from the fraction bound with iron (Kowalczewska-Madura et al. 2019). Authigenic P Figure 8 depicts the variation in concentration of P auth (L3) in cores I and II. P auth in core-II varied from 29.21 to 700.12 ppm (average: 272.01 ± 3.73 ppm). However, the concentration of L3 in Core II was found to diminish with depth. P auth (Fig. 8) occurred in greater concentrations at the surface of each Core. Percentage contribution of P auth ranged from 1.30 to 8.70% (average: 5.14%) and 1.1 to 17.2% (Average: 8.4%) in Core II (Table 6) and Core I (Table 5), respectively. At the same time, a decrease in both detrital and loosely adsorbed P was also noted, suggesting transformations of phosphate from one phase to another. Most of the remobilized phosphate in the sediment-water interface is removed as authigenic carbonate fluorapatite and does not return to the ocean reservoir (Ruttenberg and Berner 1993;Ingall and Jahnke 1994). Detritus from marine organisms is likely to be the source of P auth (Sekula-Wood et al. 2012). Likewise, conditions with higher and lower oxygen concentrations are favorable for the production of P auth . However, low levels of organic matter and medium oxygen levels limit the formation of this fraction of phosphorus (Kraal et al. 2012;Tsandev et al. 2012). As the Lo-P content is mainly related to productivity and consistent relation between the biogenic and P auth suggests a climatic connection.
Detrital P L4 exhibited a zigzag pattern in distribution. The concentration of Det-P ( Fig. 9) was lower in both cores. In Core I, it ranged from 13.97 to 120.90 ppm (average: 47.57 ± 3.08 ppm). L4 decreased at the top and increased towards the bottom in Core II. However in Core I, from the surface up to 200 cm, a gradual decrease followed by a Fig. 8 Variation of P auth with respect to depth in Core I and Core II sudden increase was observed (Fig. 10). Det-P contributed 0.2 to 2.9% of total P in Core I (Table 5), with an average of 0.9%. Whereas, in Core II (Table 6), it ranged from 1.08 to 5.7% (average: 3.29%). Det-P is usually derived from weathering of continental rocks. A distinct change at 200 cm suggested a change in provenance type. The difference in P concentration may offer variation in sources with time which can arise due to shifting of the fan limits. The cessation of turbiditic sedimentation due to the abandonment of the active channel may also be a reason for the change in P sources. Essentially two sources, viz., the Himalayan and peninsular Indian and Sri Lankan, can be responsible for a change in detrital P content. Det-P is diagenetically stable, mainly derived from marine sediments and metamorphic rocks (Ruttenberg 1992;Meng et al. 2014). Detrital P has been rarely influenced by biogenic particles due to its mineralogical properties (Ruttenberg 1992) and is therefore recognized as a permanent sink of sediment bound phosphate (Meng et al. 2014).
Organic P There was no particular trend in OP concentration with depth in Core I, even though the upper layers of the Core exhibited a gradual decrease. In contrast, a significant gradual decline was present in Core II (Fig. 10). The concentration of OP ranged from 15.15 to 295.53 ppm (average: 84.71 ± 1.22 ppm) in Core I and from 11.392 to 64.480 ppm (average, 40.31 ± 0.68 ppm) in Core II. OP (Fig. 10) exhibited a decreasing trend with increasing depth in the case of Core II, while no particular trend was observed in Core I. The percentage contribution OP to TP ranged from 0.4 to 4.7% (average, 1.5%) and 0.5 to 2.8% (average: 1.3%) in Core I (Table 5) and Core II (Table 6) respectively. Organic matter in the sediments was susceptible to alteration, and thus the P bound to organic fraction can be converted to other phases. Hence the decline in organic phosphate fraction was due to the biodegradation of OP due to mineralization. The relatively low OP arises due to the increased rate of organic matter decomposition, which releases orthophosphate into pore waters to be captured via Fe oxyhydroxides or authigenic carbonate fluorapatite (Wang et al. 2013). During early diagenesis, the P released from OP degradation is captured through authigenic carbonate fluorapatite formation of oxyhydroxides of Fe or Mn (Filippelli and Delaney 1996;Katsaounos et al. 2007;Hou et al. 2009).
Total Phosphorus (TP) and Reactive P Total concentration of phosphorous in Core I and II are depicted below (Fig. 11). In Core I, TP initially increases and follows a zigzag pattern and then exhibits a decreasing trend with depth. The concentration of TP ranged from 2976.81 to 8056.35 ppm (average: 5866.88 ± 87.32 ppm) in Core I. Concentration of TP in Core II ranged from 1599.68 to 8047.75 ppm (average: 3963.12 ± 56.90 ppm) and exhibited a decreasing trend (Fig. 11). Content of reactant P ranged from 2901.89 to 8020.09 ppm (Average, 5819.30 ± 67.18 ppm) and 1550.36 to 7960.07 ppm in Core I and Core II, respectively (Fig. 12).

Statistical Analysis
Principal component analysis was used to identify and interpret the relations among phosphorus types, grain size, and TOC. When we analyzed core I, total variance explained by PC1 and PC2 (Fig. 13) was 63.41%, with PC1 accounting for 41.41% and PC2 accounting for 22.37% of the total variance. Lo-P was positively correlated with sand (r = 0.26) and silt (r = 0.20) (Fig. 13). P auth was positively correlated with silt (r = 0.21), sand (r = 0.24), and TOC (r = 0.43). Positive relation between TOC with P auth implies the remarkable role of TOC in dispersion of this class of phosphorus  Stone and English (1993). A significant positive correlation exists between the Fe-P and TOC (r = 0.46 p = 0.046) (Fig. 13), implying that TOC is an important factor in the concentration and distribution of the P species in the sediments. Fe-P significantly correlated with the P auth (r = 0.47, p < 0.05) and OP (r = 0.79, p < 0.02) similar to other investigations (Mao et al. 2021;Souza et al. 2022). Det-P was mostly concentrated in clay (r = 0.64, p < 0.02), comparable with the observations of Yang et al. 2016, who suggested that Det-P is probably associated with the transport of terrigenous material.
In the case of Core II, the total variance explained by PC1 and PC2 (Fig. 14) was 100%, with PC1 accounting for 66.33% and PC2 accounting for 33.67% of the total variance. Lo-P positively correlated with silt (r = 0.99). It could be due to the larger surface area of the finer particles providing a higher adsorption area for phosphate (Wang et al. 2006;Barik et al. 2019). P auth positively correlated with silt (r = 0.96). Iron bound P also exhibited a similar correlation to Lo-P and P auth . Meanwhile OP correlated positively with the silt (r = 0.65), and TOC (r = 0.18).

Assessment of Bioavailability of Phosphorous
It is crucial to identify the potential bioavailable P in surface sediments because BAP can provide information on eutrophication or offset P limitation (Andrieux and Aminot 1997;Yang et al. 2016;Li et al. 2018). BAP refers to the sum of immediately available phosphorus and potential phosphorus that can be transformed into a functional form by biogeochemical process (Sonzogni et al. 1982;Wang et al. 2009). Bioavailable P includes Ex-P, Fe-P, and OP (Kang et al. 2017;Bastami et al. 2018;Loh et al. 2020;Gu et al. 2021). BAP is the upper limit of phosphorus in the sediments; it can be released into the overlying water column with changing environmental conditions (Meng et al. 2014;Kang et al. 2017;Wang et al. 2017). When sedimentary BAP is released to the overlying waters, the primary productivity in the water body will be enhanced (Coelho et al. 2004;Hou et al. 2009;Kang et al. 2017). Remineralization of OP enhances phosphorus bioavailability (Colman et al. 2005;Joshi et al. 2015). Loosely bound phosphate is the most labile P form because it is loosely sorbed, exchangeable, and water soluble. Fe-P is pH and redox-sensitive, and during anoxic conditions, it serves as a source of internal phosphate loading. During aerobic organic matter decomposition, OP is released as phosphate. The released phosphorus is available for algal growth and can further sustain the eutrophication processes (Coelho et al. 2004). Sedimentary OP can transform into Ex-P and/or Fe-P and Ca-P when environmental conditions change (Joshi et al. 2015;Kang et al. 2017). Organic matter containing phosphate decomposes and transform into inorganic form under situations like low sediment deposition and oxic condition (Colman et al. 2005;Paytan and McLaughlin 2007;Mort et al. 2010;Asmala et al. 2017). The contribution of BAP to total P content in Core I ranged from 82.2 to 98.7% (average: 90.6%) (Fig. 15). In the case of Core II, BAP accounted for 90.2 to 93.0% of total P (average: 91.6%). The high BAP content in the core samples implied the fact that the sediments in the study area act as a sink for phosphorous, and the changing environmental conditions can trigger the release of phosphate to the overlying water column leading to eutrophication.

Conclusion
Texture compositions in Core I and II were dominated by fine-grained sediments. Phosphorous fractionation revealed the fact that > 80% of sedimentary phosphorus in the study area was composed of Lo-P. P auth was perceived as an abundant phosphorus fraction in the study area. Statistical analysis revealed that phosphorous fractions correlated positively with the silt, clay, and TOC. Positive correlations of silt with all P fractions (except Det-P) suggested the role of adsorption on fine grained sediments facilitating the accumulation of phosphorous forms in both cores. The interrelationships among Fe-P, P auth, and OP suggested similar sources and patterns of behavior. The phosphate fraction composition in core I followed the trend: loosely bound > authigenic > organic phosphorous > iron bound > detrital. Whereas in core II loosely bound > authigenic > iron bound > detrital > organic phosphorous. Phosphorus is a limiting component for oceanic productivity, and its concentration was influenced not only by continental weathering but also by the physical, chemical, and biological processes in the Ocean at the time of sediment deposition, as well as post-depositional diagenetic changes. The enhanced bioavailable phosphate levels in the core sediments implied the role of sediments to act as a sink for phosphate, and the drastic alteration in environmental conditions can trigger the nutrient release to the water column and subsequent eutrophication.
Funding Open Access funding provided by the Qatar National Library.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.