*2.4. Statistical Analyses*

The analysis of bulk density data was performed using SAS 6.11 [31], a doubly repeated mixed model with time and depth as repeated factors [32]. This step was followed by analyses of effects sliced by the interacting factor and differences among least square means for significant effects. Analysis of effects was conducted starting from higher-order interactions according to Milliken and Johnson [33].
