Relating soil physical properties to other soil properties and crop yields

doi:10.21203/rs.3.rs-1627501/v1

Download PDF

Research Article

Relating soil physical properties to other soil properties and crop yields

https://doi.org/10.21203/rs.3.rs-1627501/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Purpose: Soil physical properties can be related to other soil properties and crop yields, but their evaluations as soil health indicators relating extensive soil properties and long-term crop yields need further exploration. We evaluated the long-term (14 and 36 yr) effect of cropping systems and N fertilizations on selected soil physical properties and related them to 66 soil physical, chemical, biological, and biochemical properties and mean crop yields across years at two dryland farming sites in the semiarid region of the northern Great Plains, USA.

Methods: Treatments were rotations of no-tillage and conventional tillage spring wheat (Triticum aestivum L.), barley (Hordeum vulgare L.), pea (Pisum sativum L.), and fallow with or without N fertilization. Soil samples collected in April 2019 were analyzed for physical, chemical, biological, and biochemical properties and mean crop yields determined.

Results: The average slake aggregate (ASA), wet aggregate stability index (WASI), and intact core available water holding capacity (IAWHC) were associated with most soil physical, chemical, biological, and biochemical properties and clay concentration was associated with nutrient concentrations. These parameters were also better related to mean crop yields across years than other soil physical properties.

Conclusions: Because of enhanced relationship with soil properties and crop yields and simple and inexpensive measurement, ASA can be used as a potential soil health indicator in dryland cropping systems in the semiarid regions.

Crop rotation

Crop production

Tillage

Nitrogen application

Soil characteristics

Soil quality.

Soil physical properties as important components of soil health influence water and nutrient movements, aeration, soil temperature, nutrient cycling, and root growth that affect crop yields and environmental quality (Blanco-Canqui and Benjamin 2013; Logsdon and Karlen 2004; Maddoni et al. 1999). For example, increased bulk density (BD) due to increased soil compaction results in decreased pore volume and water infiltration, increases aeration stress, lowers soil temperature and nutrient cycling, increases denitrification, losses mycorrhizal fungi, and reduces root growth (Logsdon and Karlen 2004). In contrast, increased soil aggregation enhances water and nutrient movements, reduces soil erosion, promotes C sequestration, favors microbial activity and abundance, and increases root growth and crop yields (Blanco-Canqui and Benjamin, 2013; Maddoni et al. 1999). Clay concentration is an important indicator of soil health that enhances the retention of soil water and nutrients (Norkaew et al. 2019). While increased soil water retention enhances crop yields (Logsdon and Karlen 2004; Stone and Schlegel 2010), reduced water infiltration capacity of the soil results in anaerobic condition that hampers nutrient cycling and root growth, thereby reducing crop production (Blanco-Canqui and Benjamin 2013; Stone and Schlegel 2010;).

Some soil physical properties are related to other soil properties. Numerous researchers (Blanco-Canqui and Banjamin 2013; Blanco-Canqui et al. 2011; McVay et al. 2006) reported that BD was negatively correlated, but soil aggregation, water retention, and water infiltration capacity were positively correlated to soil organic C (SOC). In contrast, Nemes et al. (2005) observed that soil hydraulic conductivity (SHC) was negatively correlated to SOC. Blanco-Canqui and Benjamin (2013) concluded that the relationship between soil physical properties and SOC are affected by soil and climatic conditions and management practices. Shaver et al. (2002) found that soil aggregation was related to clay concentration. Norkaew et al. (2019) showed that water-stable aggregation (WSA) was positively related to SOC, soil total N (STN), potential N mineralization (PNM), and microbially active C (MAC), but BD was negatively related to these parameters. Shestak and Busse (2005) found that BD was negatively correlated to macroporosity (MAP), but was not related to microbial properties. Stott et al. (2010) reported that soil aggregation and water infiltration capacity were positively correlated, but BD was negatively correlated to β-glucosidase (BG) activity. Similarly, several researchers (Dick et al. 1988; Martens et al. 1992) demonstrated that BD was positively correlated, but SHC was negatively correlated to activities of N-acetyl-β-glucosamidinase (NAG) and phosphomonoesterase (PME).

Soil physical properties also vary in their relationships to crop production. Sene et al. (1985) found that sand concentration and SHC were positively correlated, but soil aggregation was negatively correlated to crop yield. Shang et al. (2018) reported that soil water content and aggregate stability (AS) were positively related, but BD was negatively related to crop yield. Similarly, several researchers (Abu-Hamdeh 2008; Rubio et al. 2021) observed that BD was negatively related to crop yield. In contrast, Logsdon and Karlen (2004) showed that BD was not related to crop yield. Maddoni et al. (1999) showed that AS was an important indicator of soil health that related better to corn (Zea mays L.) yield than labile C and N fractions.

Although some soil physical properties are related to other soil properties and crop yields, an in-depth examination of their relationships with other soil properties and long-term crop yields are needed for their proper evaluation for potential soil health indicators. We chose six soil physical properties (clay concentration, BD, average slake aggregate [ASA], wet aggregate stability index [WASI], intact core available water capacity (IAWHC], and SHC) representing particle size, soil aggregation, compaction, water retention, and water infiltration capacity that have broad impact on soil health, such as water and nutrient movements, aeration, soil temperature, C sequestration, nutrient cycling, and microbial properties and related them to 66 other soil properties and long-term crop yields at two dryland farming sites in eastern Montana, USA. We hypothesized that these physical properties would relate to other soil properties and long-term crop yields and that some physical properties may be used as potential soil health indicators due to their close relationships with above parameters in dryland cropping system in the semiarid region. Our objectives were to: (1) determine if physical properties are related to other 66 other soil physical, chemical, biological, and biochemical properties, and (3) select which soil physical properties can be used as potential soil health indicators that are related to most soil properties and crop yields.

Experimental details

The long-term experiments occurred at two dryland farming sites (Froid and Sidney) in eastern Montana, USA. The 36-yr-old experiment at Froid (48^o 20’N, 104^o 29’W, elevation 617 m) had a Dooley sandy loam (fine loamy, mixed, frigid, Typic Argiboroll) soil with 645 g kg^− 1 sand, 185 g kg^− 1 silt, 170 g kg^− 1 clay, 14.9 g kg^− 1 SOC, and 6.2 pH at the 0–15 cm depth at the beginning of the experiment in 1983. The site had mean annual air temperature (30-yr average) of 8 ^oC and average annual precipitation of 357 mm. The 14-yr-old experiment at Sidney (48^o 33’N, 104^o 50’W, elevation 592 m) had a Williams loam (fine-loamy, mixed, superactive, frigid, Typic Argiustolls) soil with sand, silt, and clay concentrations of 350, 325, and 325 g kg^− 1 respectively, pH 7.2, and SOC 13.2 g kg^− 1 at the 0–20 cm depth at the initiation of the experiment in 2006. The site had mean annual air temperature of 8 ^oC and average annual precipitation of 340 mm. The distance between the two sites was 88 km.

Aase and Pikul (1995) described the treatments and crop management for the Froid site and Sainju and Alasinrin (2020) for the Sidney site. In short, treatments in Froid had four cropping systems which were fall and spring till continuous spring wheat (FSTCW), no-till continuous spring wheat (NTCWA), no-till spring wheat-barley (1984–1999) replaced by spring wheat-pea (2000–2019) (NTWPA), and spring till spring wheat-fallow (STWF), each replicated for four times. In each cropping system, each phase of the crop appeared in every year. Tillage in FSTCW occurred with a tandem disc in the fall and a field cultivator in the spring to a depth of 8 cm to prepare the seedbed. Similarly, tillage in STWF occurred with a field cultivator to a depth of 8 cm in the spring and during the fallow period as needed to control weeds. In other treatments, tillage was not applied. The control treatment was STWF which is a conventional cropping in the region. The size of the plot was 12 m × 30 m.

In Sidney, treatments were cropping system as the main plot and N fertilization as the split-plot factor arranged in a randomized block design with three replications. Cropping systems included conventional till barley/spring wheat-fallow (CTWF, traditional system), no-till barley/spring wheat-fallow (NTWF), no-till continuous barley/spring wheat (NTCWB), and no-till barley/spring wheat-pea (NTWPB). In each treatment, barley was grown from 2006 to 2011, which was replaced by spring wheat from 2012 to 2019. All phases of crops in the rotations occurred in each year. Tillage in CTWF was performed with a field cultivator to a depth of 8 cm to prepare seedbeds in the spring and during the fallow period to control weeds. Tillage was not applied to other treatments. Nitrogen fertilizer was broadcast to barley at 0 (N0) and 80 kg N ha^− 1 (N1) from 2006 to 2011 and to spring wheat at 0 (N0) and 100 kg N ha^− 1 (N1) from 2012 to 2019. As in Froid, N fertilization rates were adjusted to soil residual N to a depth of 60 cm determined in the autumn of the previous year. Therefore, N fertilization rates included both soil and fertilizer N. Pea did not receive N fertilizer. The size of the split plot was 12.0 × 6.0 m.

At both sites, spring wheat, barley, and pea were seeded at recommended seed rates with a no-till drill at a row spacing of 20 cm. A banded application of P and K fertilizers at 11 kg P ha^− 1 and 27 kg K ha^− 1, respectively, was done to each crop at seeding. Herbicides and pesticides were applied to each crop as needed. In mid-July to early August, crops were harvested with a combine from a swath of 11.0 × 1.5 m and grain yield was determined after oven drying subsamples at 65 ^oC for 7 d. After grain harvest, crop residues were returned to the soil.

Soil sampling and analysis of physical, chemical, biological, and biochemical properties were described in detail by Sainju et al. (2021). For determining BD, two core (3.5 cm inside diameter) samples to a depth of 7.5 cm were collected within a plot and oven dried at 110^o C for 24 hr. The BD for a treatment was calculated as the mean values of two cores by dividing the weight of oven-dried soil by the volume of the core. The SHC was determined with an infiltrometer using the two ponding head method after measuring the water flow rate to a depth of 10 cm (Reynolds and Elrick, 1990). Sand, silt, and clay concentrations were determined using the pipette method (Gee and Bauder 1986). The WASI was calculated as mean-weight diameter of aggregates that passed through 1.00, 0.25, and 0.053 m sieves after shaking in water for 10 min (Franzluebbers et al. 2000). For this, field-moist soil was oven-dried at 55 ^oC for 2–3 d and passed through a 4.75 mm sieve, out of which 19 g soil was sieved in a stack of sieves in water. The mean-weight diameter of aggregates was calculated by dividing the sum of the products of weight and diameter of each aggregate by the total weight of the aggregates. An image recognition software smartphone app “SLAKES” was used to determine ASA where non-sieved aggregates, about 4–10 mm in diameter, were dispersed in a petri dish and the picture of area covered by the aggregates was taken by a hanging smartphone above the petri dish (Fajardo et al. 1986). The areas covered by aggregates before and after immersing aggregates in water for 10 min were calculated by the software and then compared for measuring the slaking resistance of the aggregates. The IAWHC was determined as difference between volumetric water contents at 33 and 1500 kPa which were determined by using the pressure plate technique (Reynolds and Topp 2008).

Statistical analysis of data

The Pearson’s correlation analysis was used to correlate soil physical properties. To determine the associations among soil physical, chemical, biological, and biochemical properties, a multivariate regression analysis, e.g. principal component analysis (PCA) using the PROC FACTOR procedure of SAS (SAS Users’ Guide, 2019) was used to analyze the soil data. Data were loaded into two principal components (PC1 and PC2) that explained most variations among physical treatments due to treatments. Soil properties were considered positively associated to each other when they are at acute angles, negatively associated when at right angle, and not associated when at obtuse angle. Regression analysis was used to explore the relationships between soil physical properties and mean crop yields across years for individual and combined sites. For these, soil physical properties were used as independent variables and crop yields as dependent variables. Because data for crop yields were absent during the fallow periods in alternate years and also for the latest year and yields can be affected by long-term management, it was decided to use mean crop yield across years for each treatment rather than yield for the latest year to relate with physical properties. Mean crop yield was determined by dividing total yield by the number of years. Yield during the fallow period was considered zero for data analysis. Statistical significance was observed at P ≤ 0.05, unless mentioned otherwise.

Relationships among soil properties

Soil physical properties

At Froid, clay concentration was positively correlated to WASI (Table 1). The ASA was positively correlated to WASI and IAWHC, but negatively correlated to BD. The WASI was positively correlated to IAWHC and SHC. The BD was negatively correlated to IAWHC. At Sidney, WASI was positively correlated to ASA and IAWHC. Similarly, SHC was positively correlated to IAWHC.

Table 1

Correlation (r) among clay concentration (Clay), average slake aggregate (ASA), wet aggregate stability index (WASI), bulk density (BD), intact core available water holding capacity (IAWHC), and saturated hydraulic conductivity (SHC) at Froid (n = 16) and Sidney (n = 24), Montana.
Soil physical properties	Clay	ASA	WASI	BD	IAWHC	SHC
Froid
Clay (g kg^− 1)	---------	0.14	0.50*	-0.32	0.36	0.14
ASA	0.14	---------	0.66**	-0.54*	0.52*	0.28
WASI	0.50*	0.66**	---------	-0.38	0.51*	0.69**
BD (Mg m^− 3)	-0.32	-0.54*	-0.38	---------	-0.71**	-0.13
IAWHC (cm³ cm^− 3)	0.36	0.52*	0.51*	-0.71**	---------	0.14
SHC (cm hr^− 1)	0.14	0.28	0.69**	-0.13	0.14	---------
Sidney
Clay (g kg^− 1)	---------	-0.25	-0.24	0.05	0.34	0.17
ASA	-0.25	---------	0.67**	-0.22	0.38	0.11
WASI	-0.24	0.67**	---------	-0.36	0.46*	0.32
BD (Mg m^− 3)	0.05	-0.22	-0.36	---------	0.08	0.15
IAWHC (cm³ cm^− 3)	0.34	0.38	0.46*	0.08	---------	0.43*
SHC (cm hr^− 1)	0.17	0.11	0.32	0.15	0.43*	---------
* and ** Significant at P ≤ 0.05 and 0.01, respectively.

The PCA indicated that soil physical properties contributed 34% of the total variability in principal component 1 (PC1) and 20% variability in principal component 2 (PC2) at Froid (Fig. 1). Clay concentration was positively associated with volumetric water content at 0.3 and 33 kPa (VWC0.3 and VWC33, respectively), volumetric water content at water saturation (VWCFS), and MAP, but negatively associated with BD. The IAWHC was positively associated with VWC33 and WSA. The ASA and WASI were positively associated with WSA, AS, and volumetric water content at the field-moist soil (VWCFM). The SHC was positively associated with VWCFM and dry aggregate stability index (DASI), but negatively associated with stone content (SC). Soil physical properties that were not associated with clay concentration, ASA, WASI, IAWHC, SHC, and BD were sand concentration, repacked core available water holdinf capacity (RAWHC), volumetric water content at 10 kPa (VWC10), mesoporosity (MEP), volumetric water content at 1500 kPa (VWC1500), silt concentration, and total shrinkage (TS).

At Sidney, soil physical properties contributed 33% of the total variability in PC1 and 19% of the variability in PC2 (Fig. 1). The ASA was positively associated with AS and WSA and negatively associated with VWC33, VWC10, RAWHC, and VWC1500. The WASI was positively associated with MEP, VWC0.3, and AS, but negatively associated with VWC1500. The IAWHC was positively associated with DASI. The BD was positively associated with VWCFM, but negatively associated with VWC0.3. The clay concentration was positively associated with silt concentration, but negatively associated to sand concentration. Soil physical properties that were not associated with clay concentration, ASA, WASI, IAWHC, SHC, and BD were VWCFS, TS, MAP, RAWHC, and SC.

Soil chemical properties

In the PCA, soil chemical properties contributed 30% of the total variability in PC1 and 24% of the variability in PC2 at Froid (Fig. 2). The ASA and IAWHC were positively associated with Cd, inorganic P (IP), Mn, Pb, K, and Zn concentrations, but negatively associated with Na-absorption ratio (SAR). The WASI was positively associated with electrical conductivity (EC) and Zn, Mn, Pb, and K concentrations. The SHC was positively associated with Na concentration. Clay concentration was positively associated with EC and Cu, Co, B, As, Ba, and Ni concentrations, but negatively associated with organic P (OP) concentration. The BD was positively associated with SAR, but negatively with Cd concentration. Soil chemical properties that were not associated with ASA, BD, Clay, IAWHC, SHC, and WASI were pH, buffer pH (BpH), and Al, Fe, Ca, Mg, and S concentrations.

At Sidney, chemical properties contributed 41% of the total variability in PC1 and 15% of the variability in PC2 (Fig. 2). The ASA was positively associated with K, IP, and OP concentrations, but negatively associated with Ni and Mg concentrations. The WASI was positively associated with K concentration, but negatively associated with Ba concentration. The BD was positively associated with Ni and Mg concentrations. Clay concentration was positively associated with Ba and Cr concentrations. Soil chemical properties that were not associated with ASA, BD, Clay, IAWHC, SHC, and WASI were pH, BpH, EC, SAR, and Pb, As, Ca, Mg, Co, B, Cu, Na, Al, MN, S, Fe, Zn, Cd, OP, and IP concentrations.

Soil biological and biochemical properties

At Froid, the PCA indicated that soil biological and biochemical properties contributed 52% of the total variability in PC1 and 12% of the variability in PC2 (Fig. 3). The ASA and WASI were positively associated with phospholipid-derived fatty acid (PLFA), water-extractable C (WEC), CO₂ flush at 4-d incubation (CO₂-4d), KMnO₄-extractable C (POXC), soil organic matter (SOM), MAC, BG, and CO₂ flush at 1-d incubation (CO₂-1d). The IAWHC was positively associated with water-extractable total N (WETN), PNM, water-extractable organic N (WEON), CO₂-1d, and NO₃-N. Clay concentration was positively associated with PLFA. The SHC was positively associated with NO₃-N and WEON. In contrast, BD was negatively associated with PME, autoclaved citrate-extractable protein (ACEP), NAG, soil total C (STC), and STN. Soil biological and biochemical properties that were not associated with ASA, BD, Clay, IAWHC, SHC, and WASI were arylsulfatase (AST), NH₄-N, PME, NAG, ACEP, STC, and STN.

At Sidney, soil biological and biochemical properties contributed 34% of the total variability in PC1 and 16% of the variability in PC2 (Fig. 3). The ASA was positively associated with PNM, MAC, CO₂-1d, NAG, and SOM. The WASI was positively associated with SOM, NAG, and STC. The IAWHC was positively associated with PLFA. The SHC was positively associated with AST. In contrast, clay concentration was negatively associated with ACEP, WEON, and BG. The BD was negatively associated with NAG, SOM, and MAC. Soil biological and biochemical properties that were not associated with ASA, BD, Clay, IAWHC, SHC, and WASI were NO₃-N, WETN, NH₄-N, WEC, STN, CO₂-4d, and POXC.

Relationship between soil physical properties and crop yields.

Clay concentration was strongly (P ≤ 0.01) related to mean crop yield across years from the combined data from Froid and Sidney, but not for the individual sites (Fig. 4). Clay concentration accounted for 24% of the variability in crop yield. An increase in clay concentration by 1 g kg^− 1 increased crop yield by 0.03 Mg ha^− 1. The ASA was very strongly (P ≤ 0.001) related to mean crop yield at Froid and weakly (P ≤ 0.05) related at Sidney. The ASA accounted for 79% of variability in crop yield at Froid and 15% of the variability at Sidney. The WASI was strongly related to mean crop yield at Froid, weakly related at Sidney, and very strongly related to the combined data from Froid and Sidney. The WASI accounted for 50% of the variability in crop yield at Froid, 25% of the variability at Sidney, and 30% of the variability for the combined data from Froid and Sidney.

The BD was the only soil physical property that negatively related to mean crop yield (Fig. 5). The BD was strongly related to crop yield at Froid and very strongly related to the combined data from Froid and Sidney (Fig. 5). The BD accounted for 39% of the variability in crop yield at Froid and 24% variability for the combined data from Froid and Sidney. An increase in BD by 1 Mg m^− 3 decreased crop yield by 1.85 Mg ha^− 1 at Froid and 1.80 Mg ha^− 1 for the combined data from Froid and Sidney. The IAWHC was weakly related to mean crop yield at Froid and strongly related to the combined data from Froid and Sidney. The IAWHC accounted for 33% of the variability in crop yield at Froid and 22% variability for the combined data from Froid and Sidney. The SHC was not related to crop yields at both sites or for the combined data from Froid and Sidney.

Relationship among soil properties

The positive correlations among ASA, WASI, IAWHC, and SHC at Froid and Sidney suggest that well aggregated soils can also enhance soil water retention and infiltration capacities. As soil aggregation is related to SOC (Blanco-Canqui and Benjamin 2013), it is possible that increased SOC in enhanced aggregated soil helped to retain more water. Also increased soil aggregation may have increased macroporosity which increased SHC. In contrast, negative correlations among BD, ASA, and IAWHC at Froid suggest that increased soil compaction can reduce slaking resistance of soil aggregates and water holding capacity. Blanco-Canqui and Benjamin (2013) reported that well aggregated soils have lower BD, better consistency, and greater macropores and water retention capacity than poorly aggregated soils. Benjamin et al. (2007) found that soil aggregation was positively correlated to SHC, but negatively correlated to BD. Manrique and Jones (1991) reported that BD was negatively correlated to soil water content. Our result of positive correlation between clay concentration and WASI was similar to that observed by Shaver et al. (2002).

The positive associations among ASA, WASI, WSA, AS, VWC33, IAWHC, and VWCFM at Froid and Sidney suggest that increased soil aggregation enhanced aggregate stability, slaking resistance, and water holding capacity. Water-stable aggregates with enhance stability have higher slaking resistance and water holding capacity and are less erodible than aggregates with reduced stability (Blanco-Canqui and Benjamin 2013). The close association between IAWHC and VWC33 at Froid indicates that water holding capacity increased as soil water content increased. Similarly, the close association between SHC and VWCFM suggests that water infiltration rate increased as soil water content at saturation increased. The negative association between clay concentration and BD suggests that soils with high clay concentration are less dense than soils with less concentration. Similarly, the negative association between clay and sand concentrations at both sites suggests that clay concentration increased as sand concentration decreased.

The positive associations among ASA, WASI, and some nutrients at Froid and Sidney suggest that improved soil aggregation can favor retention and availability of nutrients probably by enhancing nutrient cycling. As soil aggregation is closely related to SOC (Blanco-Canqui and Benjamin 2013; Blanco-Canqui et al. 2011), it is likely that increased SOC that enhanced soil aggregation also favored nutrient cycling. Similarly, the positive associations among IAWHC and some nutrients indicate that increased water holding capacity of soil also may help to retain nutrients that may be leached out of soil profile. Increased clay concentration also may have favored retention of nutrients, as clay concentration was positively associated with nutrients at both sites. Negatively charged clay particles may have attracted cations, helping to retain positively charged nutrients. The BD and SHC had minor associations with soil chemical properties, as soil compaction and water infiltration capacity may be less important for affecting chemical properties and nutrient levels.

Greater number of soil chemical properties were associated with ASA, BD, Clay, IAWHC, SHC, and WASI at Froid than at Sidney (Fig. 2). This suggests that the longer duration of the experiment conducted in coarse-textured soil may have altered both soil physical and chemical properties that were more related to each other than the shorter duration of the experiment conducted in fine-textured soil. This indicates that long duration experiments may be more ideal to evaluate soil physical properties as potential soil health indicators relating to chemical properties than short duration experiments in dryland cropping systems.

The positive associations among ASA, WASI, and most biological and biochemical properties reveal that soil aggregation is related to SOM, microbial abundance and activity, nitrogen mineralization, and enzyme activities. The binding agents from soil organic matter and fungal growth enhance soil aggregation (Blanco-Canqui and Benjamin 2013), as fungi contribute more to soil aggregation than bacteria because of their greater abundance (Hattori and Hattori, 1976; Paul & van Veen, 1978). Some researchers (Blanco-Canqui and Benjamin 2013; Blanco-Canqui et al. 2011; McVay et al. 2006) reported that AS was positively correlated, but BD was negatively correlated to SOM. The increased SOM associated with improved soil aggregation enhances soil microbial biomass and activity and N mineralization (Blanco-Canqui and Benjamin 2013; Rivero et al. 2020). Soil aggregation has been shown to be positively correlated to C mineralization (Liebig et al. 1995; van Es and Karlen 2019), WEC (Veum et al. 2012), and BG activity (Stott et al. 2010). Norkaew et al. (2019) reported that WSA was positively related to SOM, STN, PNM, MAC, but BD was negatively related to these parameters. It is likely that soils get less dense, as SOM is increased, favoring soil biological and biochemical properties. The positive associations among IAWHC, WETN, PNM, WEON, CO₂-1d, NO₃-N, and PLFA suggest that increased soil water holding capacity can enhance N mineralization and availability by stimulating microbial biomass and activity. Rice and Havlin (1994) reported that soil water content was linearly related to PNM. Similarly, the positive associations among SHC, NO₃-N, and WEON indicate that increased water infiltration capacity may enhance N availability by increasing N mineralization through increased aeration. In contrast, the negative associations among clay concentration and ACEP, WEON, and BG indicates that increased clay concentration might inhibit N availability by reducing enzyme activity.

Relationship between soil physical properties and crop yields

The positive relationships among ASA, WASI, IAWHC, clay concentration, and mean crop yields across years suggest that improved soil aggregation, water holding capacity, and fine-textured soils increased crop yields probably by enhancing soil water retention and availability, aeration, nutrient cycling, and root growth, as soil aggregation promotes these variables (Blanco-Canqui and Benjamin 2013; Rivero et al. 2020). Several researchers (Maddoni et al. 1999; Shang et al. 2018) found that soil aggregation and stability as well as soil water content were positively related to crop yields. Our results are in contrast to those reported by Sene et al. (1985) who observed that soil aggregation was negatively related, but SHC was positively related to crop yield. In contrast, the negative relationship between BD and crop yield indicates that increased soil compaction reduced crop yield, likely due to reduced porosity, aeration, water and nutrient movements, and root growth. Many researchers (Abu-Hamdeh 2003; Rubio et al. 2021; Shang et al. 2018) reported that BD was negatively related to crop yield. However, Logsdon and Karlen (2004) demonstrated that BD had no effect on crop yield.

Implications of soil physical properties

The ASA, WASI, and IAWHC were associated with most soil physical, chemical, biological, and biological properties, followed by clay concentration (Figs. 1, 2, and 3). The BD and SHC were minimally associated with soil properties. The ASA, WASI, and IAWHC were also weakly to very strongly related to mean crop yields across years at Froid and Sidney (Figs. 4 and 5). Clay concentration, BD, and SHC were either not related or weakly to strongly related to crop yield.

These results suggest that ASA, WASI, and IAWHC can be regarded as potential soil health indicators that were associated with most soil physical, chemical, biological, and biochemical properties and related to mean crop yields across years. However, it takes long time to measure WASI and IAWHC. Soils need to be oven dried for 2–3 d and aggregates needs to be separated in a stack of sieves in water to measure WASI. Similarly, soil water contents need to be determined at various hydraulic pressures and soils oven dried for several days to determine IAWHC. In contrast, the determination of ASA is simple, rapid, and inexpensive where the slaking resistance of aggregates is measured by a smart phone loaded with SLAKES software that calculates areas covered by aggregates before after slaking in water for 10 min. Therefore, ASA can be considered as a potential soil health indicator under dryland cropping systems in arid and semiarid regions.

The ASA, WASI, IAWHC, BD, clay concentration, and SHC varied in their relationships to soil properties and mean crop yields across years under dryland cropping systems at Froid and Sidney, Montana, USA. The ASA, WASI, and IAWHC were associated with most soil physical, chemical, biological, and biochemical properties and also weakly to strongly related to mean crop yields across years at both sites. The BD, clay concentration, and SHC were less associated with soil properties and related to crop yields. As a result, ASA, WASI, and IAWHC can be used as potential soil health indicators. However, it takes long time to measure WASI and IAWHC. Because of its simplicity, rapidity, and inexpensiveness as well as strong relationships to other soil properties and crop yields, ASA may be considered as an attractive potential soil health indicator under dryland cropping systems in arid and semiarid regions.

AECP, autoclaved citrate-extractable protein; AS, aggregate stability; ASA, average slake aggregate; AST, arylsulfatase; BD, bulk density; BP, buffer pH; BG, β-glucosidase; CTWF, conventional till barley-fallow (2006-2011) replaced by conventional till spring wheat-fallow (2012-2019); DASI, dry aggregate stability index; EC, electrical conductivity; FSTCW, fall and spring till continuous spring wheat; IAWHC, intact core available water holding capacity; IP, inorganic P; MAC, microbially active C; MAP, macroporosity; MEP, mesoporosity; NAG, N-acetyl-β-glucosaminidase; OP, organic P; PCM, potential C mineralization; PLFA, phospholipid-derived fatty acid; PME, phosphomonoesterase; POXC, KMnO₄-extractable C; NTCW1, no-till continuous spring wheat; NTCW2, no-till continuous barley (2006-2011) replaced by no-till continuous spring wheat (2012-2019); NTWF, no-till barley-fallow (2006-2011) replaced by no-till spring wheat-fallow (2012-2019); NTWP1, no-till spring wheat-barley (1984-1999) replaced by no-till spring wheat-pea (2000-2019); NTWP2, no-till barley-pea (2006-2011) replaced by no-till spring wheat-pea (2012-2019); PNM, potential N mineralization; RAWHC, repacked core available water holding capacity; SAR, Na-absorption ratio; SC, stone content; SIC, soil inorganic C; SOC, soil organic C; STN, soil total N; STWF, spring till spring wheat-fallow; TS, total shrinkage; VWC, volumetric water content; VWCFM, volumetric water content in the field-moist soil; VWCFS, volumetric water content at water saturation; WASI, wet aggregate stability index; WEC, water-extractable C; WEN, water-extractable N; and WSA, water-stable aggregation.

Funding

This work was supported by the Foundation for Food and Agricultural Research (Grant no. 423926). Daniel Liptzin has received the funding.

Competing Interests

The authors have no relevant financial or non-financial interests to disclose.

Authors’ contribution

Conceptualization: D. Liptzin and U. Sainju; Methodology: D. Liptzin; Formal analysis and investigation: U. Sainju and S Dangi; Data acquisition: D. Liptzin, U. Sainju, and S. Dangi; Writing: U. Sainju; Writing - review and editing: U. Sainju, D. Liptzin, and S. Dangi; Funding acquisition: D. Liptzin; Resources: D. Liptzin; Supervision: D. Liptzin and U. Sainju.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request”.

Aase JK, Pikul JL Jr (1995) Crop and soil responses to long-term tillage practices in the northern Great Plains. Agron J 87:652–656
Abu-Hamdeh NH (2008) Compaction and subsoiling effects on corn growth and soil bulk density. Soil Sci Soc Am J 67:1213–1219
Benjamin JG, Mikha MM, Nielsen DC, Vigil MF, Calderon F, Henry WB (2007) Cropping intensity effects on physical properties of a no-till loam soil. Soil Sci Soc Am J 71:1160–1165
Blanco-Canqui H, Benjamin JG (2013) Impacts of soil organic carbon on soil physical behavior. In: Logsdon S, Berli M, Horn R (eds) Quantifying and modeling soil structure dynamics. Soil Science Society of America, Madison, WI, pp 11–40
Blanco-Canqui H, Mikha MM, Presley DR, Claasen MM (2011) Addition of cover crops enhances no-till potential for improving soil physical properties. Soil Sci Soc Am J 75:1471–1482
Dick RP, Myrold DD, Kerle EA (1988) Microbial biomass and soil enzyme activities in composted and rehabilitated skid trail soils. Soil Sci Soc Am J 52:513–516
Fajardo M, McBratney AB, Field DJ, Minasty B (2016) Soil slaking assessment using image recognition. Soil Tillage Res 163:119–123
Franzluebbers AJ, Wright SF, Studemann JA (2000) Soil aggregation and glomalin under pastures in the southern Piedmont, USA. Soil Sci Soc Am J 64:1018–1026
Gee GW, Bauder JW (1986) Particle-size analysis. In: Klute A (ed) Methods of soil analysis. Part. 1. Physical and mineralogical methods, 2nd edn. American Society of Agronomy and Soil Science Society of America, Madison WI, pp 383–411
Hattori T, Hattori Q (1976) The physical environment in soil microbiology: An attempt to extend the principals of microbiology to soil microorganisms. Crit Rev Microbio 4:423–461
Liebig MA, Jones AJ, Doran JW, Mielke LN (1995) Potential soil respiration and relationship to soil properties in ridge tillage. Soil Sci Soc Am J 59:1430–1435
Logsdon SD, Karlen DL (2004) Bulk density as a soil health indicator during conversion to no-tillage. Soil Tillage Res 78:143–149
Maddoni GA, Urricanet S, Gherse CH, Lavado RS (1999) Assessing soil quality in the Rolling Pampa using soil properties and maize characteristics. Agron J 91:280–287
Martens DA, Johnson JB, Frankenberger WT Jr (1992) Production and persistence of soil enzymes with repeated addition of organic residues. Soil Sci 153:53–61
McVay KA, Budde JA, Fabrizzi K, Mikha MM, Rice CW, Schlegel AJ, Peterson DE, Sweeney DW, Thompson C (2006) Management effects on soil physical properties in long-term tillage studies in Kansas. Soil Sci Soc Am J 70:434–438
Manrique LA, Jones CA (1991) Bulk density of soils in relation to soil physical and chemical properties. Soil Sci Soc Am J 55:476–481
Nemes A, Rawls WJ, Pachepsky YA (2005) Influence of organic matter on the estimation of saturated hydraulic conductivity. Soil Sci Soc Am J 69:1330–1337
Norkaew S, Miles RJ, Brandt DK, Anderson SP (2019) Effects of 130 years of selected cropping management systems on soil health properties for Sanborn field. Soil Sci Soc Am J 83:1479–1490
Paul EA, van Veen JA (1978) The use of tracers to determine the dynamic nature of organic matter. 11th Transaction International Congress of Soil Science, 19–27 June 1978. Edmonton, Canada, pp 61–102
Reynolds WD, Elrick DE (1990) Ponded infiltration from a single ring: I. Analysis of steady flow. Soil Sci Soc Am J 54:1233–1241
Reynolds WD, Topp GC (2008) Soil water desorption and imbibition: Tension and pressure techniques. In: Carter MR, Gregorich EG (eds) Soil sampling and methods of analysis, 2nd edn. Canadian Society of Soil Science. Taylor and Francis, Boca Raton, FL, pp 981–987
Rice CW, Havlin JL (1994) Integrating mineralizable nitrogen indices into fertilizer nitrogen recommendations. In: Havlin JL, Jacobsen JS (eds) Soil testing prospects for improving nutrient recommendation, vol 40. Soil Science Society of America, Madison, WI, pp 1–13
Rivero C, Tourn SN, Garcia GV, Videla CC, Dominguez GF, Studdert GA (2020) Nitrogen mineralization in anaerobiosis as indicator of soil aggregate stability. Agron J 112:592–607
Rubio V, Quincke A, Ernst O (2021) Deep tillage and nitrogen do not remediate cumulative soil deterioration effects of continuous cropping. Agron J 113:5584–5596
Sainju UM, Alasinrin SY (2020) Changes in soil chemical properties and crop yields with long-term cropping system and nitrogen fertilization. Agrosyst Geosci Environ 3:e20019. Doi: 10.1002/agg2.20019
Sainju UM, Liptzin D, Stevens WB (2021) How soil carbon fractions relate to soil properties and crop yields in dryland cropping systems? Soil Sci Soc Am J. DOI: 10.1002/saj2.20399
SAS User’ Guide (2019) SAS/Share 9.4: User’s guide. 2nd ed. SAS, Cary, NC. https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/shrref/titlepage.htm
Sene M, Vepraskas MJ, Naderman GC, Denton HP (1985) Relationship of soil texture and structure to corn yield response to subsoiling. Soil Sci Soc Am J 49:422–427
Shang J, Zhu Q, Zhang W (2018) Advancing soil physics for securing food, water, soil, and ecosystem services. Vadose Zone J 17:1–7
Shaver TM, Petersen GA, Ahuja AR, Westfall DC, Sherrod LA, Dunn G (2002) Surface soil physical properties after twelve years of dryland no-till management. Soil Sci Soc Am J 66:1296–1303
Shestak CJ, Busse MD (2005) Compaction alters physical, but not biological indices of soil health. Soil Sci Soc Am J 69:236–246
Stone LR, Schlegel (2010) Tillage and crop rotation phases effects on soil physical properties in the west-central Great Plains. Agron J 102:483–491
Stott DE, Andrews SS, Liebig MA, Wienhold BJ, Karlen DL (2010) Evaluation of β-glucosidase activity as a soil quality indicator for the Soil Management Assessment Framework. Soil Sci Soc Am J 74:107–119
van Es HM, Karlen DL (2019) Reanalysis validates soil health indicator sensitivity and correlation with long-term crop yields. Soil Sci Soc Am J 83:721–733
Veum KS, Goyle KW, Kremer R, Motavalli PP (2012) Relationships among water-stable aggregates and organic matter fractions under conservation management. Soil Sci Soc Am J 76:2143–2153

Datafile.xlsx

Download PDF

Version 1

posted

You are reading this latest preprint version

Relating soil physical properties to other soil properties and crop yields

Status:

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

Soil physical properties

Soil chemical properties

Soil biological and biochemical properties

Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Version 1