Carbon Sequestration in Soil Aggregates under Different Cropping Patterns of Bangladesh ()
1. Introduction
Global warming is a major threat to the environment and soil management practices are now believed to contribute significantly to changing environmental conditions [1]. Attention has been increasingly paid to soil organic carbon (SOC) pool and its dynamics in land use changes concerning terrestrial ecosystem carbon sink and the uprising atmospheric carbon dioxide [2]. The global soil carbon (C) pool of 2500 gigatons (Gt) includes about 1550 Gt of SOC and 950 Gt of soil inorganic carbon (SIC). The soil C pool is 3.3 times the size of the atmospheric pool (760 Gt) and 4.5 times the size of the biotic pool (560 Gt). The SOC pool to 1-m depth ranges from 30 tons/ha in arid climates to 800 tons/ha in organic soils in cold regions, and a predominant range of 50 to 150 tons/ha. The SOC pool represents a dynamic equilibrium of gains and losses. Conversion of natural to agricultural ecosystems causes depletion of the SOC pool by as much as 60% in soils of temperate regions and 75% or more in cultivated soils of the tropics. Severe depletion of the SOC pool degrades soil quality, reduces biomass productivity, and adversely impacts water quality, and the depletion may be exacerbated by projected global warming [1]. Carbon stock occurring in the topsoil is the most important, because this is the portion more influenced by external environmental and human factors and therefore the mainly susceptible to mineralization or synthesis processes. On the contrary, in deep layers, soil organic carbon is more stable and less liable to transformation. Soil surface layers account for the main part of the soil carbon stocks: on average, the 47% of the carbon stored in soils is held in the first 30 cm; about 2/3 are stored within a depth of 50 cm and the 80% within one meter. Soil carbon storing capability depends on many factors, such as pedoclimatic conditions, crop management practices and the starting soil carbon levels and on the interactions among them.
Urgency of meeting increased demand for agricultural produce is rapidly degrading soil quality and exacerbating degradation. Most agricultural soils have low soil organic matter (SOM) reserves due to fertility-mining practices and widespread problem of soil degradation. This decline is also attributed to removal of crop residue, and changes in cropping systems etc. Crop yield and use efficiency of input are also adversely affected by low levels of SOC pool [3].
Land use practices that result in a net accumulation of SOM in the soil are said to be C sequestering because they result in a net removal of C from the atmosphere [4]. Several studies in temperate and tropical regions reported that no-tillage practices substantially increase SOM storage and improve soil aggregation [5]. Soil aggregation can increase SOC storage by reducing loss by erosion and from mineralization. Soil organic matter can be physically protected from microbial mineralization through sorption to clay minerals and enclosure within soil aggregates [6]. Improvement in soil aggregation is an important factor influencing SOC sequestration in soils [7]. Soil aggregation is an important process of C sequestration and hence a useful strategy to mitigate the increase in concentration of atmospheric CO2 [8].
Tillage causes loss of soil organic matter through destruction of macro-aggregates and microbial mineralization of the physically protected SOM pool. Micro-aggregates and their associated SOC are stable to tillage, and this “passive or chemically protected” SOC pool represents the minimum level of SOC. Rice-wheat surface soils (0 - 15 or 0 - 20 cm) from Bangladesh increase in SOC as their silt + clay fraction increases. This variability may reflect differences in manure use, years under cultivation, tillage, sampling depth or length of time that soils are flooded [9]. Significant amounts of OC will not accumulate in soils in these environments without the protective effects of interactions with mineral surfaces (as provided by silt/clay) and the formation of aggregates. Hopelass et al. [10] reported that there was a trend of increase in concentration of SOC with decreasing aggregate size, but significant differences in these parameters in different aggregate size fractions were found only in few cases. The SOC concentration was higher in >0.25 mm than in <0.25 mm aggregates. The SOC sequestration rate by judicious use of inorganic fertilizer was the greatest in the grain-meadow rotation, while that by application of FYM was the greatest in the all grain rotation. Shrestha et al. [8] studied water stability of soil aggregates (WSA) and SOC associated with aggregates and primary particles of cultivated and forest soils. They concluded that microaggregates (<0.5 mm) were abundant (56% - 63%) in cultivated soils but cultivated soils contained higher amounts of clay but less clay-associated SOC than forest soils.
Bangladesh is a tropical agrarian country with high population density and astounding food demand. Overexploitation of agricultural soils causes severe soil degradation in Bangladesh. Depletion of SOM is a widespread problem on croplands and grazing lands in the country. Most soils have low levels of SOC contents, ranging from 8 to 10 g/kg. Low external input of organic amendment causes depletion of SOC pool because nutrients harvested in agricultural products are not replaced, and are made available through mineralization of SOM. In some cases, soil is burnt to release nutrients contained in SOM. It is estimated that SOC loss due to agricultural activities in Bangladesh between 1967 and 1995 was 16.2 Mg C/ha, with a range of 3.8 to 30.5 Mg C/ha [3]. To meet the increased food demand cropping intensity must be increased and cropping pattern should be changed. Rabbi et al. [11] reported that in silt loam soil highest percent %SOM was estimated under Fallow-Fallow-T. Aman (2.90%) and lowest under Vegetable-Sugarcane/Jute pattern (0.86%). The soils under Vegetable, Vegetable-Vegetable, Fallow-Jute-T. Aman, Fallow-Aus-Vegetable and Rabi-Jute-T. Aman patterns had relatively lower %SOM than other patterns. The %SOM under Sesame-Fallow-T. Aman pattern varied in different agricultural fields and it was 1.22% to 2.11%. In loam soils highest %SOM was obtained under Rabi-Jute-Turmeric (1.54%) and lowest under Boro-Fallow-T. Aman pattern (0.55%). Soils under Fallow-Fallow-T. Aman, Boro-T. Aman/Fish and Fallow-Jute-T. Aman patterns also had higher %SOM whereas Soils under Boro-Shrimp had lower %SOM in loam soil. The %SOM in silty clay loam soil under Boro-Fallow-T. Aman varied between 0.96% and 1.54%. The highest %SOM (2.05%) was obtained under Wheat-Jute-T. Aman pattern in silty clay soil and lowest (0.76%) in clay loam soil under Vegetable-Aus/Jute/Vegetable pattern. Soils under Fallow-Fallow-T. Aman in silt loam soil had higher %SOM than in loam soil under same pattern. The higher %SOM was obtained under Boro-Fallow-T. Aman pattern in silty clay loam soil than loam soils. It was reported that micro-aggregation of soils was also influenced by cropping patterns of the studied agricultural fields. The aggregate stability of soils of Ganges Floodplain is closely related with soil organic matter [12]. Therefore, the objectives of the research project were to evaluate the state of SOC sequestration in the agricultural soils under different cropping patterns and estimate the carbon stock of soils.
2. Methods and Materials
The study was conducted on agricultural soils of Jessore Sadar and Bagherpara Upazila under Jessore District, Phultala Upazila under Khulna District and the Nachole Upazila under Chapainawabganj District (Figures 1-5). General information’s of sampling sites are given in Table 1. Thirty-six soil samples (Core and Bulk sample) were collected from 0 - 100 cm soil depth and all soil samples were kept in sealed plastic bags. Marking and labelling was performed with a detailed description of the selected sampling site on both the soil-plastic bags, and preserved in plastic bags until arrival at the laboratory for sample preparation. Then samples were air dried by spreading the soils on separate sheet of papers. After drying in air, the larger aggregates were broken through gentle crushing with a wooden hammer. A portion of the crushed soils was passed through by using different size sieves. The ranges were >2000 µm, 2000 - 250 µm, 250 - 53 µm and <53 µm aggregate size range. The sieved soils were then preserved in plastic bags and labeled properly. These samples were later used for various chemical analyses.
2.1. Soil Physical Properties
The particle size analyses of the soils were carried out by combination of sieving and hydrometer method as described by Gee and Bauder [13]. Textural classes were determined by using Marshall’s Triangular Coordinator system. Bulk density of soil was determined by core method as described by Blake and Hartge
Table 1. Sampling locations, cropping patterns with physiography.
Figure 1. Sampling sites at Jessore sadar.
Figure 4. Sampling sites at Batiaghata upazila.
Figure 5. Sampling sites at Nachol upazila.
[14]. Particle density was determined by pycnometer meter as described by Blake and Hartge [14] and total porosity was calculated from bulk and particle density.
Aggregate Stability (Normalized Stability Index) by Wet Sieving of the Aggregates
The stability of aggregates was determined by the method as described by Six et al. [15]. For the determination of aggregate stability soil samples were air dried and crushed by a wooden hammer. The crushed soils were then sieved through 8 mm sieve. The air-dried soils that were passed through 8 mm sieve but retained on 2 mm sieve divided into 8 - 2 mm, 2 - 0.25 mm and 0.25 - 0.05 mm size fractions by hand sieving. For wet sieving with slaking pretreatment 10 grams air dried samples from each aggregate size fraction were submerged for 5 minutes on the top of smaller sieve of each size range prior to sieving. Soils were separated manually by moving the sieve 3 cm up and down under water with 50 repetitions during a period of 2 minutes. This manual separation technique was repeated for each size fractions. For wet sieving with wetted pretreatment the air-dried samples were adjusted to field capacity by soaking with water for overnight before submerging in water. The soils were then sieved for 2 minutes by the method as stated before. The amount of aggregates retained after sieving was oven dried at 105˚C for 24 hours and then weighed. The number of primary particles retained on the sieves during wet sieving was determined by sieving after dispersing the soils with 5% sodium hexametaphosphate. The weight of primary particles was recorded after oven drying at 105˚C for 24 hours.
The normalized stability index (NSI) of aggregates was calculated by the following [5].
NSI = 1 − [DL/DL(max)]
The whole soil disruption level (DL) was calculated as:
where, n = number of aggregate size classes.
i = 1 for the smallest size class.
The disruption level of a size class upon slaking (DLSi) was calculated by the following formula:
where,
DLSi = disruption level for each size class i; Pio = proportion of total sample weight in size class i before disruption (i.e., rewetted); Pi = proportion of total sample weight in size class i after disruption (i.e., slaked); Sio = proportion of sand with size i in aggregates of size i (=aggregate-sized sand) before disruption; Si = proportion of sand with size i in aggregates of size i in aggregates after disruption.
The whole soil DL(max) was calculated by the following formula:
The maximum disruption [DLSi(max)] was calculated with the following formula:
Pp = primary sand particle content with the same size as the aggregates size class after complete disruption of the whole soil.
2.2. Soil Chemical Properties
The soil pH was measured with the help of a glass electrode pH meter using soil water suspension ratio of 1:2.5 as described by Jackson [16]. Electrical conductivity (EC) of soil was estimated by EC meter (Jackson EC meter). Maintain the ratio with of soil to water of 1:5 and the result was converted to the ratio of 1:1 (soil:water) as suggested by USDA [17]. The soil organic carbon was determined by wet oxidation method of Walkley and Black’s method as described by Jackson [16]. The soil organic carbon associated with >2000 μm, 2000 - 250 μm, 250 - 53 μm and <53 μm aggregates was determined. Stock of soil organic carbon was evaluated by the method of Ellert et al. [18]. Active soil organic carbon was determined by the method as described by Weil et al. [19].
2.3. Statistical Analyses
Statistical analysis was performed by Microsoft Excel and SPSS window version 17.0 (SPSS Inc., Chicago, USA) to test the difference between means and correlation (statistical significance) of individual treatments [20].
3. Results and Discussion
3.1. Particle Size Distribution
The term particle size distribution of a soil refers to the percentage distribution of various sized particles in a given volume of soils. Particle size distribution is one of the most stable soil characteristics, being little modified cultivation or other practices. Although the usefulness of particle size analysis in practical agriculture has sometimes questioned, its indirect benefits have been extensive.
Samples studied were different in both of texture and cropping pattern. The percentage of clay was the highest (48%) under Rice-Colocasia cropping pattern at sampling spot 11 and was the lowest (15%) under Vegetables/Mustard-Fallow-T. Aman cropping pattern. Silt was the highest (68%) under Vegetables/Mustard-Fallow-T. Aman and the lowest (40%) under Wheat-Dhaincha-T. Aman cropping pattern with a mean value of 55.14. The percentage of sand was the highest (23%) under Mahogany garden and the lowest (7%) under both Ridge Gourd/Ash gourd-T. Aman and Rabi vegetables-T. Aman cropping patterns with a mean value of 15.38 (Table 2).
Table 2. Physical properties of soil samples.
The texture of the samples varied within three soil texture groups, Silt Loam, Silty Clay Loam and Silty Clay (Table 2). This result resembles to several findings by many other researchers. SRDI staff [21] found that on the young Ganges meander floodplain, soils of ridges and inter-ridges depressions are silt loam to silty clay; and on the old Ganges meander floodplain, soils of ridges and depressions are loamy to clay in texture. They found that most ridge soils in old Brahmaputra floodplain are silt loam to silty clay loam and in inter-ridge depressions they are mostly silty clay loam to silty clay subsoil. Joshua and Rahman [22] found that the soils of the Tista floodplain, in general, contained a high percentage of silt and this property appears to be characteristic for the floodplain.
3.2. Bulk Density, Particle Density and Porosity
Bulk density varied in the range of 1.18 to 1.84 g/cc with mean value of 1.40 g/cc, particle density in the range of 2.44 to 2.59 gm/cc, porosity in the range of 25.88% - 52.98% in soil samples under different cropping patterns (Table 2). The variations in the particle density, bulk density and porosity values under different cropping patterns were significant (p < 0.01). To a remarkable degree, increased organic matter can counteract the ill effects of too much clay or too much sand. Increasing the SOM content usually increases total porosity and therefore decreases bulk density [23]. At very high levels of SOM, additional OM has little further effect on soil aggregation and influences bulk density mainly because of its low particle density [24].
3.3. Normalized Stability Index (NSI)
In this research we found that the NSI varied from 0.51 to 0.89 under different cropping patterns (Appendix TableA1). The highest NSI value was found under Wheat-Fallow-T. Aman cropping pattern in Silty clay soils (sample No 15) and lowest value was found Vegetables/Mustard-Fallow-T. Aman cropping pattern in Silt loam soils (sample No 17).
Normalized Stability Index (NSI) of Silt Loam soils under different cropping patterns are presented in Figure 6. In Silt Loam soils the highest NSI value (0.87) was found under Cauliflower/Pumkin/Spinach-T. Aman (Sample No 2) and the lowest value (0.51) was obtained under Vegetables/Mustard-Fallow-T. Aman cropping pattern (Sample No 17).
Normalized Stability Index of Silty Clay Loam soils under different cropping patterns are presented in Figure 7. In S Silty Clay Loam soils the highest NSI value (0.87) obtained under Chickpea/mustard-T. Aman (Sample No 31) and the lowest value (0.59) obtained under Vegetables/Mustard-Fallow-T. Aman cropping pattern (Sample No 24).
Figure 6. Normalized stability index of silt loam soils under different cropping patterns.
Figure 7. Normalized stability index of silty clay loam soils under different cropping patterns.
Normalized Stability Index of Silty Clay soils under different cropping patterns are presented in Figure 8. In Silty Clay soils the highest NSI value (0.89) obtained under Wheat-Fallow-T. Aman (Sample No 15) and the lowest value (0.61) obtained under Lentil-Fallow-T. Aman cropping pattern (Sample No 13). Rabbi et al. [12] reported that aggregate stability of soils of Ganges River and Tidal Floodplains of Khulna region increased with increasing clay percentage of soils.
The value of NSI can vary between 0 to 1 [15]. The lower NSI indicated that the aggregates were not water stable. Earlier investigation by Rabbi et al. [12] indicated that water stability of aggregates of silt loam texture was low. The mineralogy of soils may play important role in aggregate stability [15]. However, in the present study mineral identification was not done. In the present study the NSI increased with decreasing cropping intensity. Soils under Fallow-Fallow-T. Amon cropping pattern had higher NSI. The cultivation of rice in winter and monsoon season also decreased NSI of soils. Six et al. [25] reported that conventional tillage with high cropping intensity caused maximum destruction of soil. It has been reported that stable aggregate of tilled soils is lower than that of no-tilled soils [26] due to aggregate breakdown in tillage process.
The cultivation of rice twice a year accelerates the aggregate destruction rate. Adiku et al., [27] concluded that soils under rice base cultivation were more prone to degradation. The inclusion of shrimp in land use markedly decreased the NSI may be due to application of saline water in soil. Rabbi et al. [11] reported that shrimp culture can deteriorate the tilth of soil. Cropping patterns have pronounced effect on soil organic matter content and hence affect the tilth of soil. Soil structure plays a dominant role in controlling microbial access to organic substrates. The labile organic material may be physically protected from
Figure 8. Normalized stability index of silty clay soils under different cropping patterns.
decomposition by its incorporation into soil aggregates [28]. Soil organic matter can be: 1) physically stabilized, or protected from decomposition, through micro aggregation, or 2) intimate association with silt and clay particles, and 3) can be biochemically stabilized through the formation of recalcitrant SOM compounds [5]. Hence, Soil aggregation is an important process of C sequestration [8].
3.4. Soil Chemical Properties
The pH of the soils under different cropping patterns varied from 5.99 to 7.85 (Table 3). The highest pH value was found in Ganges Meander Floodplain soils of Bahadurpur, Jessore (Sample No 9) under Vegetables/Mustard-Fallow-T. Aman cropping pattern and the lowest was found in soils of High Barind Tract under Chickpea + Linseed-T. Aman cropping pattern (Sample No 32).
The EC value of studied soils varied from 0.29 to 2.46 dS∙m−1 (Table 3). These results indicate that all the soils of the study area were non saline to slightly saline. The highest value of EC was obtained in soils of Ganges Tidal Floodplain under Bean-Fallow-T. Aman (Sample No 24) and lowest value of EC was obtained in soils of High Barind Tract under Rabi vegetables-T. Aman (Sample No 12).
3.5. Analyses of Different Carbon Forms of Studied Soil
3.5.1. Active Carbon
Active carbon of soil samples under different cropping patterns varied from 619.23 to 735.20 mg∙kg−1 with a mean value of 712.16 mg∙kg−1 (Appendix TableA1). The highest value was observed under Chickpea/mustard-T. Aman (Sample No 31) and the lowest value was found in case of Wheat-Fallow-T. Aman cropping
Table 3. pH and EC of soil samples.
pattern (Sample No 30). Active carbon contents are presented in Figures 9-11 in Silt Loam, Silty Clay Loam and Silty Clay soils, respectively.
Figure 9. Active C (mg/kg) in Silt loam soils.
Figure 10. Active C (mg/kg) in silty clay loam soils.
Figure 11. Active C (mg/kg) in silty clay soils.
The active fraction of soil organic matter consists of materials with a relatively high average C/N ratio (about 15 to 30) and short half-lives (half of these materials can be metabolized in a matter of a few months to a few years). Components include the living biomass, some of the fine particles of detritus (referred to as particulate organic matter, or POM), most of the polysaccharides, and other non humic substances, as well as some of the more labile (easily decomposed) fulvic acids. This active fraction provides most of the readily accessible food for the soil organisms and most of the readily mineralizable nitrogen. It is responsible for most of the beneficial effects on structural stability that lead to enhanced infiltration of water, resistance to erosion, and ease of tillage. The active fraction can be readily increased by the addition of fresh plant and animal residues, but it is also very readily lost when such additions are reduced or tillage is intensified. This fraction rarely comprises more than 10% to 20% of the total soil organic matter [29].
The content of soil organic carbon (SOC) can be affected by many factors, such as forest types [30], soil moisture, soil type, temperature and precipitation [31]. In our result, highest active carbon was under Chickpea/mustard/tomato-T. Aman. It may be the cause of addition of fresh plant and animal residues, soil moisture, soil type, soil organic carbon stock and organic carbon %.
3.5.2. Aggregate Associated Soil Organic Carbon (SOC)
Soil organic carbon associated with different size aggregates were the highest (3.14%) under Mungbean/Ash gourd-T. Aman (Sample No 20) and were the lowest (0.36%) under Cauliflower/Pumkin/Spinach-T. Aman cropping pattern (Sample No 2). Aggregate associated SOC contents are presented in Figures 12-14 in Silt Loam, Silty Clay Loam and Silty Clay soils, respectively.
Organic carbon content in aggregate size ranges >2000 µm (SOC1), 2000 - 250 µm (SOC2), 250 - 53 µm (SOC3), and <53 µm (SOC4) varied from 0.36% - 1.90%,
Figure 12. Aggregate associated C in silt loam soils under different patterns.
Figure 13. Aggregate associated C in silty clay loam soils under different patterns.
Figure 14. Aggregate associated C in silty clay soils under different patterns.
0.52% - 2.10%, 0.50% - 2.60% and 0.50% - 1.62%, respectively (Appendix TableA2). The percentages of SOC associated with <53 µm aggregates were higher than those of >2000 µm, 2000 - 250 µm and 250 - 53 µm, aggregates.
In Silt Loam soils (Figure 12), the highest SOC1 (1.23%) was found in Vegetables/Mustard-Fallow-T. Aman cropping pattern (Sample No 17) and the lowest SOC1 (0.36%) was found in Cauliflower/Pumkin/Spinach-T. Aman cropping pattern (Sample No 2). The highest SOC2 (1.59%) was found in Boro-Aus-T. Aman cropping pattern (Sample No 34) and the lowest SOC2 (0.51%) was found in Mustard-Fallow-T. Aman cropping pattern (Sample No 7). The highest SOC3 (1.69%) was found in Mahogany garden and Boro-Aus-T. Aman cropping pattern (Sample No 6 and 34, respectively) and the lowest SOC3 (0.49%) was found in Mustard-Fallow-T. Aman cropping pattern (Sample No 7). The highest SOC4 (1.69%) was found in Mahogany garden and Boro-Aus-T. Aman cropping pattern (Sample No 6 and 34, respectively) and the lowest SOC4 (0.49%) was found in Mustard-Fallow-T. Aman cropping pattern (Sample No 7).
In Silty Clay Loam soils (Figure 13), the highest SOC1 (1.90%), SOC2 (2.10%), SOC3 (2.60%) and SOC4 (3.14%) was found in Mungbean/Ash gourd-T. Aman cropping pattern (Sample No 20) and the lowest SOC1 (0.46%), SOC2 (0.51%), SOC3 (0.49%) and SOC4 (0.49%) was found in Chickpea/mustard-T. Aman cropping pattern (Sample No 31).
In Silty Clay soils (Figure 14), the highest SOC1 (1.31%) and SOC2 (1.29%) was found in Lentil-Fallow-T. Aman cropping pattern (Sample No 13) and the lowest SOC1 (0.58%) was found in Vegetables/Mustard-Fallow-T. Aman cropping pattern (Sample No 9). The lowest SOC2 (0.75%) was found in Wheat-Dhaincha-T. Aman cropping pattern (Sample No 16). The highest SOC3 (1.59%) was found in Bean-Falow-T. Aman cropping pattern (Sample No 12) and the lowest SOC3 (0.73%) was found in both Khesari-Fallow-T. Aman and Wheat-Dhaincha-T. Aman cropping pattern (Sample No 14 and 16 respectively). The highest SOC4 (1.69%) was found in Vegetables/Mustard-Fallow-T. Aman cropping pattern (Sample No 9) and the lowest SOC4 (0.79%) was found in Khesari-Fallow-T. Aman cropping pattern (Sample No 14).
High cropping intensity under conventional cultivation commonly has decrease in organic carbon content, because the amount of organic material returned to the soils is considerably lower and tillage enhance the decomposition of native soil organic matter [32]. Paddy soils with periodic submergence were slightly higher in organic carbon content than upland soils.
The average SOC associated with aggregates was in the order of SOC4 > SOC3 > SOC2 > SOC1. The fresh organic carbon that was derived from the crops were first incorporated to larger aggregates and then shunted to micro aggregates after decomposition by soil microbes and this process was stimulated by disturbance by conventional tillage [33]. The percentage of SOC associated with 2000 - 250 µm, 250 - 53 µm and <53 µm aggregates increased with increasing clay percentages of soils. Wiseman and Puttmann [34] described the importance of specific surface of clays rather than percentage of clays in SOC sorption. Wattel-Koekkoek et al. [35] showed that smectites have large sorptive capacity for SOC. The existing reports on clay minerals of Ganges River Floodplain of Bangladesh concluded that illite is the dominant clay mineral of this floodplain [36]. The surface area of illite is about 70 - 120 m2∙g−1. So, the capacity of illite to sorb SOC at may play an important role and it requires further research to conclude.
3.5.3. Soil Organic Carbon Stock (SOC Stock)
In order to estimate the potential of carbon sequestration in soils the original C stocks in soils need to be determined. Soil organic carbon stock of Silt Loam, Silty Clay Loam and Silty Clay soils under different cropping pattern ranged from 0.35 to 1.07 kg C m−2, 0.42 to 1.62 kg C m−2, 0.48 to 1.07 kg C m−2, respectively (Figures 15-17).
The highest organic carbon stock was found in Silty Clay Loam soil under Mungbean/Ash gourd-T. Aman cropping pattern (Sample no 4). The lowest SOC stock was found in Silt Loam soil under Cauliflower/Pumkin/Spinach-T. Aman Cropping pattern (Sample No 2).
Batjes [37] discussed the total soil C stock distribution by major ecological zones. Such zones show large differences in organic carbon storage, mainly in relation to temperature and rainfall. Soil C stocks down to 1 m depth range from about 4 kg∙m−2 (in the arid zone) to 21 - 24 kg∙m2 (in polar or boreal regions); with intermediate values of 8 to 10 kg∙m−2 in tropical zones. But our value ranged from 1.07 ± 0.01 to 0.35 ± 0.04 kg C m−2 (top soil). This lower range may be due to the erosion of top soil and intense cultivation. The contribution of the tropical regions to the global pool of soil carbon is 384 - 403 Pg C to 1 m and 616 - 640 Pg C to 2 m depth [38], compared to about 1500 Pg C to 1 m for the world (2736 - 2456 Pg to 2 m depth). The arid zone, which covers 40 percent of the global land surface, stocks only 5 percent (100 Pg) of the total. These agro-ecological
Figure 15. Soil organic carbon stock in silt loam soil.
Figure 16. Soil organic carbon stock in silty clay loam soil.
Figure 17. Soil organic carbon stock in silty clay soil.
zones, developed by FAO, can constitute a reference framework to evaluate and monitor soil C storage in soils.
3.6. Relationship among Soil Properties
Significant positive correlations were found between SOC stock and SOC1, SOC stock and SOC2, SOC stock and SOC3, SOC stock and SOC4 (Appendix TableA3). The SOC content in all four size ranges under this study was positively correlated with each other.
4. Conclusion
The distribution of aggregate size classes was influenced by land cultivation, as tillage destroyed especially aggregates > 1000 µm. Average SOC stocks in the soils under investigation were in order of Silty Clay Loam > Silty Clay > Silt Loam. The average SOC associated with aggregates was in the order of SOC4 > SOC3 > SOC2 > SOC1. The aggregate associated SOC concentration was greater for microaggregates (<250 µm) than macroaggregates (>250 µm) in the soils.
Data Availability
The data used to support findings of this study are available from the corresponding author upon request.
Acknowledgements
This study was made possible with funding from University Grants Commission of Bangladesh for climate change issue. The support rendered by Khulna University towards implemented of the study is acknowledged. We would like to express our sincere gratitude to our supervisor Associate Professor Md. Sadiqul Amin for his continuous supervision, guidance, inspiration and thoughtful suggestion during completion of the work. Individual efforts alone can never contribute in totally to a successful completion of any venture. We would be failing in our duty if we did not state our gratitude and appreciation to the following individuals who have made the valuable contribution towards the review work.
Appendix
Table A1. Descriptive statistics: bulk density, particle density, porosity, sand, silt, clay, NSI, pH, EC, OC, SOC1, SOC2, SOC3, SOC4, SOC stock and active carbon.
Table A2. Normalized stability index and different carbon forms.
Table A3. Correlations among selected soil properties.
Cell Contents: Pearson correlation; P-Value.