Md. Noor E. Alam Siddique
Soil Resource Development Institute, Ministry of Agriculture, Rangpur 5402, Bangladesh
Lisa Lobry de Bruyn
School of Environmental and Rural Science, University of New England, Armidale, NSW 2350, Australia
Chris N. Guppy
School of Environmental and Rural Science, University of New England, Armidale, NSW 2350, Australia; cguppy@une.edu.au (C.N.G.); yosanai@une.edu.au (Y.O.)
Yui Osanai
School of Environmental and Rural Science, University of New England, Armidale, NSW 2350, Australia; cguppy@une.edu.au (C.N.G.); yosanai@une.edu.au (Y.O.)
Soil carbon; Soil pH; Paddy soils; Land inundation; Crop intensification; Fertilizers
Dinajpur district
Crop-Soil-Water Management
Carbon sequestration
2.1. Study Area Description Dinajpur district lies between 260040 north latitude and 890180 east longitude in the Northwestern region of Bangladesh. The total area of Dinajpur is 3438 km2. The cropped area was 2707 km2 in 2014–2015. Dinajpur is among one of the most intensively used agricultural areas in Bangladesh. The region has a humid, wet and hot subtropical climate with distinct summer, monsoon and winter seasons. The majority of soils are located in three physiographic types, i.e., Piedmont plain, Tista Floodplain and Barind Tract/Terrace. The soils of Piedmont plain and Tista Floodplain are noncalcareous grey soils (i.e., Gleysols) and Terrace is shallow grey soils (i.e., Planosols). The lands of Dinajpur possess three land types, based on flooding during the monsoon and/or flood season, which are HL (i.e., land above the normal flooding level) and MHL (i.e., land flooded up to 90 cm for at least two weeks) and the remainder is MLL (i.e., land flooded up to 90–180 cm for more than two weeks). The surface soil (i.e., 0–15 cm generally) texture is mainly loam and silt loam, but varies from silt loam to sandy loam and silty clay loam to clay loam [33]. Farmers throughout the year depend on monsoon rain and/or irrigation, and commonly practice double crops (i.e., Aman rice in monsoon and Boro rice in winter/dry season rice) or triple crops (i.e., Aman and Boro rice with rotation crop) in MHL and HL or single crop (i.e., Aman rice) in MLL areas. Land management is characterized by conventional farming, agrochemicals, irrigation in dry season and high yielding crop varieties. Puddling (i.e., wet-tillage) is a traditional soil management practice for land preparation before transplanting of seedlings in paddy fields, which involves plowing and harrowing of surface soil in water-saturated conditions. After paddy crop harvest, crop residues and stubbles are often removed for animal fodder and homestead purposes in Dinajpur. 2.2. Soil Legacy Data, Cropping Intensity and Fertilizer SOC and pH datasets from two separate periods (i.e., 1990s and 2010s) were extracted from legacy soil data that was collected for a national semi-detailed soil survey program, conducted by the Soil Resource Development Institute in Bangladesh. The first soil survey initiative was published in the Land and Soil Resource Utilization Guide (LSRUG) (recognized as Upazila Nirdeshika) for each upazila’s (sub-district) in Bangladesh. The second phase of the survey (i.e., follow-up soil survey) was published as reports with a soil polygon map (1:50,000). Dinajpur district consists of thirteen upazila, soil survey datasets were accessed from the primary survey (first phase) and follow-up survey (second phase). These soil datasets are referred to as “legacy data of 1990s and 2010s periods”. The potential contribution of increasing cropping intensity and fertilization on the trends of SOC and soil pH was explored from 1990 to 2019 in Dinajpur. The cropping intensity (%), which is [(gross cropped area (total cropping area sown once and/or more in a year)/Net-cropped area (total cropping area) × 100)]; and fertilizers (Nitrogen (N), Phosphorus (P), Potassium (K), Sulphur (S)) usage (kg/ha) data from 1990 to 2019 were collected from the District Agriculture Office and Bangladesh Bureau of Statistics (BBS). 2.3. Statistical Analysis For data analysis and visualization, the statistical software R for Windows version 1.2.5019 and Microsoft Excel were used. Statistical significance was taken at the p < 0.05 level. Data were tabulated and visualized for physiography, land types and the interactions between the two factors. The contents of SOC and soil pH in physiographic types, land types were visualized in line and bar graphs. To establish the relationships statistically among the variables, SOC and soil pH were considered as dependent variables, and soil sampling years (i.e., time periods 1990s and 2010s), physiography and land type were considered as independent variables. Three-way partial (unbalanced) factorial linear regression models were carried out for the analysis of variance. The partial factorial model was used because the land type category MLL is only present in the Floodplain and is absent in the Piedmont plain and Terrace physiography, so that the model was not constrained by the number of categories within variables. For comparing the means of SOC and soil pH between the 1990s and 2010s for each physiography across land types the Fisher’s Least Significance Difference (LSD) post hoc multiple comparison tests was applied. Prior to analysis, the box-cox transformation was applied to the dependent variables (i.e., SOC and soil pH), ANOVA assumptions were checked and the normality distributions were assessed through Q-Q plots.
Agronomy 2021, 11, 59. https://www.mdpi.com/journal/agronomy
Journal