agronomy Assessment of the Short-Term Impact of Anaerobic Digestate on Soil C Stock and CO 2 Emissions in Shallow Water Table Conditions

: Anaerobic digestion has been recently proposed as a more sustainable energy supply chain able to strengthen the existing security of supply provisions. Anaerobic digestate (AD) is the by-product of the anaerobic digestion process and presents organic fertilizer characteristics, but its agronomic usage has been hindered by both agronomic and environmental concerns. In this study, the impact of agricultural liquid (LD) and solid (SD) AD and beef manure (Man) was evaluated on some agro-environmental characteristics. First, the AD fractions functional groups were characterized by FT-IR, and then LD and SD performances were evaluated on soil organic carbon (SOC) stock and carbon dioxide (CO 2 ) emissions and compared to beef manure (Man). The fertilizer impact was also considered with the presence or absence of a shallow water table. Results showed that SD increased of ca. 3 t ha − 1 the SOC stock in the 0–20 cm soil proﬁle with respect to LD, while it was comparable to Man despite the lower C input (1.7 vs. 8.7 t C ha − 1) , due to the presence of more stable compounds (e.g., lignin). The CO 2 ﬂuxes were affected by both fertilizer type and water table level. In the absence of a water table, the CO 2 emissions (5.5 g CO 2 m 2 d − 1 ) were driven by carbon content and quality in the fertilizer, while the presence of a shallower water table hindered mineralization of stable SOC and, in turn, reduced emissions (4.4 g CO 2 m 2 d − 1 ). AD can be considered a beneﬁcial solution to both maintain soil fertility and, at the same time, give new insight into a circular economy model, although further investigation on GHG emissions is required.


Introduction
Biomethane production has been recently proclaimed by the European Commission as a renewable and low-carbon gas able to strengthen the resilience of the gas system [1]. According to the 2018/844/EC [2], biogas and biomethane will account for 32% of the renewable energy share of EU energy consumption after 2020. Biomethane is produced in anaerobic digestion plants (i.e., biogas plants) where the anaerobic digestion process takes place, producing biogas and anaerobic digestate as a byproduct. It has been estimated that approximately 180 million tons of anaerobic digestate are produced annually in the EU [3]. For these reasons, the treatment of digestate is receiving increasing attention to cope Some studies have been conducted in the area to test the effects of manure fertilization on nitrogen leaching [24] and N 2 O emissions [25]. However, to the best of our knowledge, a gap exists in the evaluation of the short-term effects of anaerobic digestate fertilization on SOC stock and soil CO 2 emission in relation to groundwater depth. For these reasons, this study aimed to assess the impact of solid (SD) and liquid digestate (LD) fractions on SOC stock and CO 2 emissions as a function of their functional groups as derived by FT-IR spectroscopy. A lysimeter study was conducted involving different water table levels and beef manure (Man) as a standard organic fertilizer. Our starting hypotheses were that (i) SD would show a similar impact on SOC and CO 2 fluxes due to its similar functional group to Man, and (ii) LD would exhibit lower potentialities than SD on SOC sequestration but also have less impact on soil C emissions due to lower C input.

Experimental Design
The experiment was conducted at the Experimental Farm of the University of Padova in northeastern Italy (45 •  The local climate is sub-humid, with yearly average temperature and rainfall of 12 • C and 800-850 mm, respectively. Annual rainfall is mostly concentrated during the autumn and spring months. Reference evapotranspiration (ETo) is 800 mm with a peak in July (4.5 mm d −1 ). The study site was originally set up in 1984 and consists of twenty drainable lysimeters, each with dimensions of 1 m × 1 m × 1.5 m (length × width × depth). The bottom of the lysimeter creates a funnel that is connected via an underground drainpipe (slope of 1‰) to an external tube and a valve regulating both the water table level and the leaching discharge. Each lysimeter was filled in 1984 with soil excavated from the adjacent experimental farm, preserving the original soil horizons. To facilitate drainage of water and prevent soil washout, a 15 cm thick layer of gravel (30-50 mm diameter) covers the bottom of each lysimeter. The soil is a fluvi-calcaric cambisol (CMcf), representing ca. 50% of the low Venetian plain. Over 30 years of operation, the soil has been managed with ordinary fertilization techniques (i.e., farmyard manure or slurry integrated with mineral fertilizer) and cropped with maize (Zea mays L.), wheat (Triticum aestivum L.), and sunflower (Helianthus annuus L.). Treatments consisted of a combination of three groundwater conditions-a free drainage (FD) and two shallow water table levels set at −120 cm (WT120) and −60 cm (WT60)-combined with three organic fertilization types (beef manure, Man; liquid digestate, LD; and solid digestate, SD), for a total of nine treatments.
Water input was regulated throughout the experiment by a plastic roof that automatically closed to cover the lysimeter site during rain events. The cover not only prevented natural uncontrolled rainfall into the lysimeters, but also protected the crop from extreme weather, such as hailstorms. A total water input of 415 mm was provided to the system by a series of irrigation events during the spring-summer growing season.
Water table conditions were set according to the communicating vessels principle. Each lysimeter, at the bottom, was directly connected to an external column of the same height as the lysimeter (150 cm) and a diameter large enough to avoid the capillary effect. The column was equipped with a water level sensor, allowing the water level to deviate ±10 cm from the set reference level.
The experiment followed a randomized block design with two replicates. Therefore, 18 lysimeters out of the 20 available were used. CS635 TDR (Time Domain Reflectometry) probes (Campbell Scientific Inc., Lincoln, NE, USA) were installed in the lysimeters at depths of 15, 30, and 60 cm. The TDR probes, which were calibrated in our laboratory, had an accuracy of ±1.5%. The water-filled pore space (WFPS, %) was calculated as the ratio of water content to total porosity. A thermocouple system monitored soil temperature at depths of 0-15 cm and 15-30 cm. Soil volumetric water content and temperature were obtained at the same time as CO 2 measurements. Tillage, fertilization, and seeding operations were conducted on the same day. Tillage was manually performed on each lysimeter by means of spring spading at 20 cm deep (mid-April) simulating a harrowing and was followed by fertilization and maize (Zea mays L.) sowing at a density of 8 plants m −2 . The digestate used in this experiment was collected from a biogas plant fed with beef manure and silage maize. The total Kjeldahl nitrogen (TKN) content ranged between 2.1 and 9.6 g kg −1 (Table 1). All the tested organic fertilizers were applied according to the Nitrate Directive at a 300 kg N ha −1 rate and were subsequently followed by 20 cm soil tillage. The resulting C input from fertilizers was 8.7, 1.7, and 6.7 t C ha −1 for Man, LD, and SD, respectively.

Emission of CO 2 from Soil
To monitor the emission of carbon dioxide (CO 2 ) from the soil, an automated closedchamber dynamic system was installed [26], with one chamber per lysimeter (12 total). Each chamber was constructed with a steel cylindrical base (20 cm ID, 8 cm H), a motorized steel lid with a neoprene cover over its inner surface, and a rubber O-ring to guarantee a tight seal. A vent reduced wind-induced pressure differences between the chamber headspace air and the surrounding atmosphere during chamber closure [26][27][28][29]. A steel collar (20 cm ID, 10 cm H) was installed in the soil at a 7 cm depth in the center row of each lysimeter between two maize plants, carefully oriented to prevent shadowing that could affect soil temperature change. The chamber was inserted into the collar to prevent lateral gas diffusion. The chamber isolated the area of soil on which it is located, and gas exchange occurred vertically between the area of soil and the chamber interior. Chambers were only removed from the site during tillage operations. Analysis of CO 2 gas emissions was performed using an infrared gas analyzer (IRGA). Each chamber underwent six measurements per day, one measurement every four hours (i.e., 0:00, 4:00, 8:00, 12:00, 16:00, and 20:00), during the entire experimental period. Each measurement deployed 70-90 s after the closure of the lid and consisted of measuring CO 2 air concentration increase every 1 or 1.6 s. A CR-1000X Datalogger (Campbell Scientific, Inc., USA) controlled all monitoring operations (chamber closure, activation and heating of the IRGA analyzer, air circulation pumps, and aperture or closure of chamber valves). The CO 2 fluxes across the soil atmosphere were calculated from the linear or non-linear time-dependent increase in headspace CO 2 concentration [28][29][30][31]. Due to the limited number of closed chambers, emissions were monitored in two contrasting treatments, SD and Man, which were factorially combined with three WT levels and replicated twice (12 chambers in total).

Soil Samplings and Laboratory Analysis
Composite soil samples were collected from 0-5, 5-10, and 10-20 cm soil layers the day before (t 0 ) and 180 days after (t 1 ) fertilization. The soil was air-dried, 2 mm sieved, and analyzed for soil organic carbon (SOC) with a vario MACRO elemental analyzer (Elementar Analysensysteme GmbH, Langenselbold, Germany) after an acid pre-treatment for inorganic C removal. Soil bulk density was measured in the 0-5, 5-10, and 10-20 cm layers with the core method [32] and used to calculate the SOC stock.
To perform FT-IR characterization of SOC, bulk soil samples were taken from the soil profile at 0-20 cm. Replicates from each sampling were mixed to obtain a pooled sample. SOC extraction was performed on air-dried soil with 0.1 M NaOH, as previously described by Nardi et al. [33]. In this study, the SOC extracted from −60 cm WT lysimeters was not listed due to the similarity of the FT-IR spectra to those at −120 cm WT. SOC samples were Agronomy 2022, 12, 504 5 of 14 named as follows: SOC free drainage (SOC FD) and SOC −120 cm WT (SOC WT120). The compositions of the applied Man, LD, and SD were also analyzed after being freeze-dried before extraction with 0.1 M NaOH. The extracted organic C was processed with the same method as described for soil. All extracts were lyophilized before analysis.
The FT-IR spectra were measured using an ALPHA FTIR spectrometer (Bruker Optics, Ettlingen, Germany) equipped with an attenuated total reflectance (ATR) sampling device containing a diamond crystal. The spectra were recorded from 4000 cm −1 to 400 cm −1 at a spectral resolution of 4 cm −1 with 64 scans co-added and averaged. A background spectrum of air was carried out before each series of measurements. Spectra were processed with Grams/386 spectroscopy software (version 6.00, Galactic Industries Corporation, Salem, NH, USA). The overlapping bands from 1800 to 1200 cm −1 were resolved using a peak fitting analysis. The best fitting parameters were resolved with a Gaussian function determined by minimization of the reduced chi-square (χ 2 ). Good agreement was obtained between the experimental and calculated profiles, with the coefficients of determination, R-squared, ranging from 0.998 to 0.988 and the standard error from 0.002 to 0.005. All data are expressed as area percentages.

Statistical Analysis
The SOC data was analyzed with three-way ANOVA, considering water table level, fertilizer type, time, and their interactions. A two-way ANOVA, considering water table level and fertilizer type was then applied for CO 2 flux analysis. Post hoc pairwise comparison of the least-squares means was performed with the Tukey method to adjust for multiple comparisons at the significance level of p ≤ 0.05. The CO 2 fluxes did not follow a normal distribution and were log-transformed before the ANOVA. Statistical analysis was performed using SAS (SAS Institute Inc., Cary, NC, USA) version 5.1.
Soil temperatures in the first layer (0-15 cm) changed by month from 17.2 • C in April to 26.5 • C in July, on average ( Figure 1). Shallow water table conditions (WT60) induced an average temperature decrease (1 • C, corresponding to ca. 4%) with respect to FD. A similar pattern was also observed in the 15-30 cm soil layer, where the average temperature ranged between 16.8 • C and 26.0 • C, and WT60 reduced soil temperature by 0.8 • C on average (−4%).

Soil Organic Carbon (SOC) Stock
The SOC stock was 17.5-32.5 t C ha −1 in the 0-20 cm soil profile with an average stock of 6.2, 6.4, and 12.5 t C ha −1 in the 0-5, 5-10, and 10-20 cm layers, respectively. The SOC stock showed similar results for SOC content with any fertilizer effect at 0-5 cm (6.2 t ha −1 , Agronomy 2022, 12, 504 7 of 14 on average) and a higher value in Man compared to LD (6.8 vs. 5.9 t ha −1 ) at 5-10 cm, as well as in SD compared to LD (13.5 vs. 11.4 t ha −1 ) at 10-20 cm (Table A2, Figure 3a). Considering the entire 0-20 cm soil profile, the SD and Man SOC stock was ca. 3 t ha −1 greater than LD.
The SOC stock decreased in all three studied layers from t 0 to t 1 with an overall depletion of 4.0 t ha −1 among the 0-20 cm soil profiles, with a greater magnitude in the 10-20 layer (−18%) than in the 5-10 (−14%) or 0-5 cm (−8%) layers (Figure 3a). The water table effect was ranked as follows: WT60 > WT120 > FD, with a clearer effect at 0-5 cm (Figure 3b). The SOC stock was 17.5-32.5 t C ha −1 in the 0-20 cm soil profile with an average stock of 6.2, 6.4, and 12.5 t C ha −1 in the 0-5, 5-10, and 10-20 cm layers, respectively. The SOC stock showed similar results for SOC content with any fertilizer effect at 0-5 cm (6.2 t ha −1 , on average) and a higher value in Man compared to LD (6.8 vs. 5.9 t ha −1 ) at 5-10 cm, as well as in SD compared to LD (13.5 vs. 11.4 t ha −1 ) at 10-20 cm (Table A2, Figure  3a). Considering the entire 0-20 cm soil profile, the SD and Man SOC stock was ca. 3 t ha −1 greater than LD.

FT-IR Characterization of Soil Organic Carbon (SOC)
The most striking structural variations in the functional groups were observed in the region between 1800-1200 cm −1 (spectra not shown). Similar functional groups were found in all spectra and are described as follows: the band at 1722-1680 cm −1 was due to H-bonded acid/ketone carbonyl groups and carbonyl in aromatic carboxylic acids [34]; the band at 1660-1640 cm −1 was mainly attributed to amide I and C=O stretching in ketones; the band at 1616-1509 cm −1 was highly characteristic of the aromatic ring and carboxylate ions [34]; the band at 1550-1535 cm −1 was mainly assigned to amide II, the aromatic ring, and carboxylate ions in the lignin; at lower frequencies, the band at 1458-1430 cm −1 was due to C-H stretching, whereas the presence of the band at 1451-1341 cm −1 was assigned to symmetric CH3 stretching of branched chains [34].
By applying the Gaussian curve fitting function in the region between 1800-1200 cm −1 , a quantitative estimation was obtained (Figure 4). The applied LD and SD differed in terms of area percentage for the examined functional groups. More specifically, the band at ~1720 cm −1 (H-bonded acid/ketone carbonyl groups) contributed 12% in LD and 4% in raw SD (Figure 4). In applied LD, amide I (1659 cm −1 ) and amide II (1542 cm −1 ) Different letters indicate significant differences between treatments according to the Tukey test with p ≤ 0.05. Lower case letters refer to the specific soil layer (0-5, 5-10, and 10-20 cm) while uppercase letters refer to the entire 0-20 cm soil profile. LD: liquid digestate; SD: solid digestate; Man: beef manure; WT60: −60 cm water table; WT120: −120 cm water table; FD: free drainage.

FT-IR Characterization of Soil Organic Carbon (SOC)
The most striking structural variations in the functional groups were observed in the region between 1800-1200 cm −1 (spectra not shown). Similar functional groups were found in all spectra and are described as follows: the band at 1722-1680 cm −1 was due to H-bonded acid/ketone carbonyl groups and carbonyl in aromatic carboxylic acids [34]; the band at 1660-1640 cm −1 was mainly attributed to amide I and C=O stretching in ketones; the band at 1616-1509 cm −1 was highly characteristic of the aromatic ring and carboxylate ions [34]; the band at 1550-1535 cm −1 was mainly assigned to amide II, the aromatic ring, and carboxylate ions in the lignin; at lower frequencies, the band at 1458-1430 cm −1 was due to C-H stretching, whereas the presence of the band at 1451-1341 cm −1 was assigned to symmetric CH 3 stretching of branched chains [34].
The composition of the functional groups in SOC-FD fertilized with LD showed a similar pattern to that of applied LD, while the percentages of the corresponding areas varied. In particular, the band area at 1717 cm −1 decreased by 1%. In contrast, amide I and amide II contributed 29% and 15%, respectively, and the aromatic ring (1599 cm −1 ) contributed 28%. Such values were higher than LD. The aliphatic groups at 1457 cm −1 , 1401 cm −1 , and 1349 cm −1 accounted for 1.3%, 18%, and 7%, respectively. For SOC-WT120, the area percentage differed from applied fertilizer and SOC-FD in LD. The band at 1717 cm −1 disappeared, while amide I and amide II accounted for 45% and 7.5%, respectively. The aromatic ring at 1582 cm −1 was 17% and decreased with respect to SOC-FD (30%). The appearance of the band at 1513 cm −1 in the lignin accounted for 3%. The band at 1397 cm −1 exhibited a marked increase up to 22%.
In SOC-FD fertilized with SD, there was a slight decrease in the percentage of band area at 1680 cm −1 (3%), amide I (24%), and amide II (17%) compared to applied SD (Figure 4). The aromatic ring at 1589 cm −1 accounted for 22%, and a new band at 1509 cm −1 was 5%. The appearance of the band at 1407 cm −1 accounting for 12% was particularly significant. In SOC-WT120, the band at 1719 cm −1 accounted for 5%, while amide I (1647 cm −1 ) and amide II (1542 cm −1 ) were 42% and 8%, respectively. The other bands at 1589 cm −1 and 1512 cm −1 decreased compared to SOC-FD.

Soil Carbon Dioxide (CO2) Fluxes
The CO2 fluxes responded to fertilizer treatment and crop growth as shown by their evolution over time in Figure 5. In SD, it is seen as a double bell-shaped curve characterized by progressive increments in CO2 fluxes up to values > 6 µmol CO2 m −2 s −1 ca. 2 weeks after the fertilizer application (Figure 5a)  In the applied Man ( Figure 4) the most important bands were those of amide I (1650 cm −1 ) and amide II (1524 cm −1 ), accounting for 21% and 4%, respectively. The most prominent band at 1583 cm −1 was 39%. In addition, the aliphatic groups at 1398 cm −1 and 1345 cm −1 accounted for 14% and 20%, respectively. In SOC-FD fertilized with Man, amide I and amide II increased, accounting for 43% and 8%, respectively. The band at 1584 cm −1 was 19% and decreased with respect to applied Man. Other relevant bands at 1511 cm −1 and 1399 cm −1 accounted for 2.5% and 20%, respectively. In SOC-WT120, amide I and amide II were 28% and 4%, respectively. They markedly decreased with respect to SOC-FD. The largest contribution was made by the band at 1589 cm −1 , which accounted for 41%. The aromatic ring in lignin (1512 cm −1 ) was 2.6% and did not differ from that found in SOC-FD. A minor contribution was made by the aliphatic groups.

Soil Carbon Dioxide (CO 2 ) Fluxes
The CO 2 fluxes responded to fertilizer treatment and crop growth as shown by their evolution over time in Figure 5. In SD, it is seen as a double bell-shaped curve characterized by progressive increments in CO 2 fluxes up to values > 6 µmol CO 2 m −2 s −1 ca. 2 weeks after the fertilizer application ( Figure 5a). Afterward, values decreased to a minimum < 3 m −2 s −1 on 22 May 2015 and increased again to values > 5 µmol CO 2 m −2 s −1 following the crop development. A progressive CO 2 flux decline (<3 µmol CO 2 m −2 s −1 ) from the end of July to the postharvest period coincided with the maize physiological stage. The total CO2 fluxes were in the 0.2-60.9 g CO2 m 2 d −1 range. The SD application produced lower fluxes at WT60 (4.4 g CO2 m 2 d −1 ) than at WT120 and FD (average of 5.9 g CO2 m 2 d −1 ) ( Table 2). Conversely, Man application yielded higher results in WT conditions (average of 5.8 g CO2 m 2 d −1 ) than in FD (4.0 g CO2 m 2 d −1 ) ( Table 2). Analyzing the 2-week span that followed fertilizer application, the total CO2 fluxes were not In Man, CO 2 fluxes exhibited a sharp increase just one day after fertilizer distribution, reaching peaks of 13, 7, and 7 µmol CO 2 m −2 s −1 , respectively, in WT60, WT120, and FD. Afterward, CO 2 fluxes slowly decreased to a minimum in mid-May, then increased again, mirroring the evolution observed in SD during the crop development and postharvest stages (Figure 5b).
The total CO 2 fluxes were in the 0.2-60.9 g CO 2 m 2 d −1 range. The SD application produced lower fluxes at WT60 (4.4 g CO 2 m 2 d −1 ) than at WT120 and FD (average of 5.9 g CO 2 m 2 d −1 ) ( Table 2). Conversely, Man application yielded higher results in WT conditions (average of 5.8 g CO 2 m 2 d −1 ) than in FD (4.0 g CO 2 m 2 d −1 ) ( Table 2). Analyzing the 2-week span that followed fertilizer application, the total CO 2 fluxes were not significantly affected by WT levels in SD, despite a tendency of higher value in FD (4.5 g CO 2 m 2 d −1 ) than in WT conditions (3.5 g CO 2 m 2 d −1 , on average). On the contrary, Man WT60 exhibited a higher CO 2 flux (8.2 g CO 2 m 2 d −1 ) than WT and FD (4.3 g CO 2 m 2 d −1 ) ( Table 2). Topsoil (0-15 cm) and subsoil (15-30 cm) monthly temperatures were significantly correlated to CO 2 fluxes for SD (r = 0.70 and r = 0.66, respectively), while no dependency was observed for Man. A positive correlation was also found between CO 2 fluxes and 60-90 cm WFPS (r = 0.53) for Man.

Discussion
Anaerobic digestion is an interesting solution for agriculture as its residue, the digestate, can be used as an organic fertilizer due to its nutrient content [35]. Indeed, digestate is generally characterized by a high content of recalcitrant organic molecules (i.e., lignin is not degraded) and nutrients (N and P) that are crucial due to their availability to plants after field application [36]. Greater lignin content in anaerobic digestate is expected since lignin is not converted into biogas during the anaerobic digestion of silage. Therefore, fertilization with digestate is a way to increase SOC and ensure the input of organic substances and nutrients into the soil.
The presented outcomes highlighted an appreciable SOC increase in the short-term, resulting from treatment with SD and Man in the layers between 5-10 and 10-20 cm. In addition, the impact was visible in shallow groundwater conditions, with a higher carbon accumulation in WT60 than in WT120 and FD (8.9 vs. 8.2 and 7.8 t C ha −1 , respectively). Conversely, no significant difference was observed for LD. The amount and quality of C input in the fertilizers, as well as water table conditions, played key roles in SOC stock dynamics. Comparable results were observed in SD and Man treatments despite the lower C input of the former (6.7 vs. 8.7 t ha −1 ). The SD impact on SOC stocks was most likely amplified by the recalcitrance nature of the C input. Indeed, in SD, the contribution of the aromatic ring (1513-1509 cm −1 ) as a lignin derivate was present in SOC-FD (27%) and SOC-WT120 (19%). The aerobic condition observed in FD could have increased the proportion of lignin derivate because of the enhanced degradation of more C labile forms (e.g., amino acids). Moreover, oxidative enzymatic activities and lignin breakdown could have been inhibited by high concentrations of N in the soil [37]. In addition to these effects, reactions between polyphenols and amino acids may also produce relatively toxic compounds and recalcitrant aromatic groups that may inhibit lignin breakdown [38]. The presence of proteins or amino acids in anaerobic digestate has been well-recognized [35].
In our experiment, the contribution of amide I was particularly relevant, as it was due almost entirely to C=O stretching motions of peptide bonds (about 80%); indeed, this spectral region represents the most sensitive secondary structural component of proteins to environmental effects. The amide II band, on the contrary, is much less sensitive. In particular, amide I in SOC-FD was 29% and 24%, respectively, for LD and SD, but under higher WFPS (i.e., WT), values increased to 45-42%. The presence of poorer mixing of the soil with the solid digestate could have hindered the degradation. Indeed, the presence of organic fertilizer "clouds" was previously related to reduced access to solid manure for soil decomposers [39]. In this study, this effect might have been exacerbated by the higher moisture content in soils with a shallow (e.g., −60 cm) water table.
The potential of SD for increasing SOC stock was previously highlighted by Slepetiene et al. [36], who tested the solid and liquid digestate fractions in eroded soils from Lithuania and found that solid digestate positively affected soil fertility, quality, sustainability, and durability. The underperformance of LD on SOC stocks was also described by Slepetiene et al. [36]. In our experiments, the poor results obtained from LD on SOC stock were mainly related to the low C input (1. 7 t C ha −1 ). Tambone et al. [40] deduced that, because of the presence of biochemically protected SOC, the fertilizing capacity of LD could not be overlooked. Indeed, according to Greenberg et al. [41], LD improved occluded particulate organic matter by up to 11% compared to mineral fertilizers.
The effect of the water table on digestate turnover was also confirmed by the CO 2 fluxes measured during the two weeks after fertilizer distribution. The CO 2 fluxes generally varied according to the seasons and the type of soil disturbance, e.g., soil tillage and type of soil covering [42,43]. Nevertheless, the first peak emissions results were dominated by SD distribution, with higher values in FD with WFPS levels averaging at 55%. These results confirmed the outcomes obtained by Schaufler et al. [21], who observed higher CO 2 emissions in croplands with WFPS levels between 20% and 60%. The second peak occurred during the vegetative corn stage and mainly was related to the root development. Roots can contribute to CO 2 emissions directly through respiration and indirectly by supplying C to the soil through exudation [22]. According to Rochette et al. [44], the contribution of roots and associated microorganisms in a maize field reached up to 45% of the total respiration, mainly concentrated in mid-summer. The total CO 2 fluxes in the seeding-harvesting period were lower in SD WT60 than in the other treatments. It is reasonable to assume that the lower SD mineralization in the shallow water table may have reduced the N availability and root growth.
The Man fertilization in shallow water table conditions showed an opposite behavior in the 2 weeks after application. Indeed, CO 2 emissions were negatively affected by water table depth, as demonstrated by Crézé and Madramootoo [22]. The authors observed that, after fertilization under free drainage and subsurface irrigation with a target water table depth of 0.75 m, CO 2 fluxes from subsurface-irrigated plots were 21% greater than those of free drainage. However, the authors were not able to find any significant correlation to seasonal changes in soil WFPS. On the contrary, in our study, the CO 2 fluxes were dependent on the WFPS of deeper layers.
The rapid response of CO 2 fluxes to Man distribution was probably due to the greater content of labile C forms in manure, boosted by the anaerobic conditions ensured by the shallow water table. Indeed, Man fertilization led to a higher contribution of amide I (43%) in SOC-FD that decreased to 28% under shallow water conditions, suggesting enhanced degradation processes at high WFPS values. As previously observed in the same experiment by Cocco et al. [25] these conditions might have fostered denitrification processes and, in turn, CO 2 emissions. As a consequence of labile compound degradation, the fraction of aromatic rings increased as a result of their high recalcitrance.
In the same experimental area, Maucieri et al. [45] observed high CO 2 emissions during the first hours after application of LD, followed by a rapid flux reduction that reached values similar to unamended soil after 48 h. The CO 2 peaks after fertilizer application might be the result of the rapid release of CO 2 dissolved in the organic fertilizers and the microorganism degradation of easily degradable C substrates. This was particularly visible when manure was applied on lysimeters with a 60 cm water table, probably due to the interaction between the greater content of labile C forms in manure and the anaerobic conditions ensured by the shallow water table.

Conclusions
Agro-biomass recycling is a smart alternative to improve soil carbon sequestration, as it is a valuable source of nutrients for crops and has a beneficial effect on the physicochemical properties of the soil. Direct application of anaerobic digestate, including solid and liquid fractions, could be considered a beneficial solution both to maintain soil fertility and, at the same time, give new insight into a circular economy model. Starting from the raw product, by mechanical separation it is possible to obtain two different fractions that can be used for fulfilling contrasting agronomic requirements. Currently, the use of liquid digestate fraction is strongly encouraged by the high cost of synthetic fertilizers. On the contrary, solid fraction digestate can be beneficially exploited in the carbon credits market.
Nevertheless, knowledge of GHG emissions from the soil both during and after digestate application requires more extensive studies.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Comparison of significance levels among the linear mixed-effect models analysis of the soil organic carbon content for each studied layer (0-5, 5-10, and 10-20 cm). Significant effects (p ≤ 0.05) are reported in bold.