There was an evaluation × treatment interaction on fluxes of N2O (P < 0.05), with average emissions ranging from 2 g N2O-N ha-1day-1 to all treatments in day 0 (after feces and dung beetle application) and significantly greater peak with 46 and 80 g N2O-N ha-1day-1 on day 6 to pot with just dung and pot with dung + dung beetle species, respectively (Fig. 1). The fluxes of N2O from pot with dung beetle species were the greatest and differed significantly to pot with dung and pot with just soil over time (P < 0.001), except in day 2, when N2O flux was greater to just dung than dung + beetles.
Fluxes of N2O were relatively high over time for pot with dung + beetles and were greater than the ones observed for soil and soil + dung (P < 0.001). However, O. taurus and D. gazella had the least N2O flux compared with other beetle treatments over time (Fig. 2). The N2O flux from the pot with dung (treatment 2) increased over time, however, since day 6 the fluxes decreased considerably until day 24, and the values ranged from 45 to 2.9 g N2O-N ha-1day-1, respectively. The treatment 3 (O. taurus) and treatment 4 (D. gazella) showed the most depleted N2O-N among the dung beetle treatments in days 0, 1, and 12, with averages of -3, 12.3, 25.8 g N2O-N ha-1day-1 and − 1, 17.5, 23.5 g N2O-N ha-1day-1, respectively. Treatment 7 (O. taurus + D. gazella + P. vindex) presented the greatest pick of N2O-N in day 6 with average of 145.7 g N2O-N ha-1day-1; however, in day 12 the N2O flux decreased considerable, not differing (P > 0.05) from T3, and T4 Fig. 2). Treatment 1 (control pot with just soil) had the least N2O flux and did not show significant variation over time. Treatment 5 showed a progressive increase over time, with the greatest pick of N2O-N in day 12 and day 24.
There was an evaluation × treatment interaction (P < 0.05) on ammonia volatilization, which varied from a maximum of 6431 g NH3-N ha− 1 for T6 in day 2, decreasing over time, showing the least value in day 24, with average of 241 g NH3-N ha− 1. The treatment T3 and T4 presented the least values with 1536 and 1575 g of NH3-N ha− 1, respectively, when compared to other beetle treatments and T2 (Fig. 4). The T4 presented the most depleted NH3-N emission from day 6, 12, and 24 with average of 1526, 1048, and 245 g of NH3-N ha− 1 when compared to T2 and to other beetle treatments. The T5 showed a peak on day 6, which was greater (P < 0.05) than T1, T2, T3, and T4. The T1 presented the least NH3-N emission, and it did not significantly vary over time (Fig. 3).
T1: just soil, T2: soil + dung, T3: soil + dung + O. taurus (OT), T4: soil + dung + D. gazella (DG), T5: soil + dung + P. vindex (PV), T6: soil + dung + OT + DG, T7: soil + dung + OT + DG + PV. (*): Indicates significant difference at the 0.05 probability level among treatments in the same month, according to orthogonal contrast test.
At the end of the experiment, the dung removal efficiency was not statistically analyzed. However, according to a visual observation, the pot with just dung did not show any changes over time, losing humidity and getting a superficial crust (Fig. 4). On the other hand, the dung beetle species on the pot surface buried almost all the dung into the soil as much as the single and assemblage species.
T2: soil + dung, T3: soil + dung + O. taurus (OT), T4: soil + dung + D. gazella (DG), T5: soil + dung + P. vindex (PV), T6: soil + dung + OT + DG, T7: soil + dung + OT + DG + PV.
There was a treatment effect (P < 0.05) on soil nitrogen content. The two control treatments (soil and soil + dung) did not differ (P > 0.05) among them. The soil from pot with dung + beetle application (treatment 3, 4, 5, 6 and 7) presented greater nitrogen concentration compared with treatment 1 and 2 (Fig. 5).
T1: just soil, T2: soil + dung, T3: soil + dung + O. taurus (OT), T4: soil + dung + D. gazella (DG), T5: soil + dung + P. vindex (PV), T6: soil + dung + OT + DG, T7: soil + dung + OT + DG + PV. Different letters in lowercase indicate statistically significant differences among treatments, according to student test.
All treatments with dung beetle species resulted in taller pear millet than just dung on the pot surface, with both groups differing significantly (P < 0.10). Millet plants under dung beetle effect had 41.8 cm, significantly greater than millet plants with just dung and no beetle, showing 39.9 cm (Fig. 6).
*Different letters in lowercase indicate statistically significant differences at the 0.05 probability level among treatments in the same month, according to single degree of freedom polynomial contrast test.
There was a harvest × treatment interaction (P < 0.05) on HA of pear millet. Dung application had a positive effect on the HA. The greatest HA was observed in all treatments with dung application in the first harvest with average of 8 g of DM pot− 1 (pot area of 0.32 m2), greater (P < 0.05) than T1 (control with just soil), which averaged 5 g of DM pot− 1. In the second harvest, treatments did not differ among them, presenting the least DM values (Fig. 7
There was a treatment effect (P < 0.05) on N yield (Fig. 8). Treatments T3, T6, and T7 resulted in greater N yield for the pearl millet than T1, T2, T4, and T5. The T3, T6, and T7 also had more soil N available, with average of 0.34, 0.30, and 0.31 g N pot− 1, respectively.
T1: just soil, T2: soil + dung, T3: soil + dung + O. taurus (OT), T4: soil + dung + D. gazella (DG), T5: soil + dung + P. vindex (PV), T6: soil + dung + OT + DG, T7: soil + dung + OT + DG + PV. Different letters in lowercase indicate statistically significant differences among treatments, according to orthogonal contrast test.
In this experiment, we hypothesized that dung beetle increased the N2O-N and NH3-N volatilization because of faster dung N mineralization in the soil as a result of dung incorporation and breakdown of dung pat crust, enhancing availability of oxygen for nitrification to occur (height and HA; Figs. 7 and 8). On the other hand, to simplify the response of the treatments of this experiment on gas emissions and nutrient cycling variables, a principal component analysis was used (Fig. 9).