Application of biogas slurries from energy crops to arable soils and their impact on carbon and nitrogen dynamics

The use of renewable primary products as co-substrate or single substrate for biogas production has increased consistently over the last few years. Maize silage is the preferential energy crop used for fermentation due to its high methane (CH4) yield per hectare. Equally, the by-product, namely biogas slurry (BS), is used with increasing frequency as organic fertilizer to return nutrients to the soil and to maintain or increase the organic matter stocks and soil fertility. Studies concerning the application of energy crop-derived BS on the carbon (C) and nitrogen (N) mineralization dynamics are scarce. Thus, this thesis focused on the following objectives: I) The determination of the effects caused by rainfall patterns on the C and N dynamics from two contrasting organic fertilizers, namely BS from maize silage and composted cattle manure (CM), by monitoring emissions of nitrous oxide (N2O), carbon dioxide (CO2) and CH4 as well as leaching losses of C and N. II) The investigation of the impact of differences in soil moisture content after the application of BS and temperature on gaseous emissions (CO2, N2O and CH4) and leaching of C and N compounds. III) A comparison of BS properties obtained from biogas plants with different substrate inputs and operating parameters and their effect on C and N dynamics after application to differently textured soils with varying application rates and water contents. For the objectives I) and II) two experiments (experiment I and II) using undisturbed soil cores of a Haplic Luvisol were carried out. Objective III) was studied on a third experiment (experiment III) with disturbed soil samples. During experiment I three rainfall patterns were implemented including constant irrigation, continuous irrigation with periodic heavy rainfall events, and partial drying with rewetting periods. Biogas slurry and CM were applied at a rate of 100 kg N ha-1. During experiment II constant irrigation and an irrigation pattern with partial drying with rewetting periods were carried out at 13.5°C and 23.5°C. The application of BS took place either directly before a rewetting period or one week after the rewetting period stopped. Experiment III included two soils of different texture which were mixed with ten BS’s originating from ten different biogas plants. Treatments included low, medium and high BS-N application rates and water contents ranging from 50% to 100% of water holding capacity (WHC). Experiment I and II showed that after the application of BS cumulative N2O emissions were 4 times (162 mg N2O-N m-2) higher compared to the application of CM caused by a higher content of mineral N (Nmin) in the form of ammonium (NH4+) in the BS. The cumulative emissions of CO2, however, were on the same level for both fertilizers indicating similar amounts of readily available C after composting and fermentation of organic material. Leaching losses occurred predominantly in the mineral form of nitrate (NO3-) and were higher in BS amended soils (9 mg NO3--N m-2) compared to CM amended soils (5 mg NO3--N m-2). The rainfall pattern in experiment I and II merely affected the temporal production of C and N emissions resulting in reduced CO2 and enhanced N2O emissions during stronger irrigation events, but showed no effect on the cumulative emissions. Overall, a significant increase of CH4 consumption under inconstant irrigation was found. The time of fertilization had no effect on the overall C and N dynamics. Increasing temperature from 13.5°C to 23.5°C enhanced the CO2 and N2O emissions by a factor of 1.7 and 3.7, respectively. Due to the increased microbial activity with increasing temperature soil respiration was enhanced. This led to decreasing oxygen (O2) contents which in turn promoted denitrification in soil due to the extension of anaerobic microsites. Leaching losses of NO3- were also significantly affected by increasing temperature whereas the consumption of CH4 was not affected. The third experiment showed that the input materials of biogas plants affected the properties of the resulting BS. In particular the contents of DM and NH4+ were determined by the amount of added plant biomass and excrement-based biomass, respectively. Correlations between BS properties and CO2 or N2O emissions were not detected. Solely the ammonia (NH3) emissions showed a positive correlation with NH4+ content in BS as well as a negative correlation with the total C (Ct) content. The BS-N application rates affected the relative CO2 emissions (% of C supplied with BS) when applied to silty soil as well as the relative N2O emissions (% of N supplied with BS) when applied to sandy soil. The impacts on the C and N dynamics induced by BS application were exceeded by the differences induced by soil texture. Presumably, due to the higher clay content in silty soils, organic matter was stabilized by organo-mineral interactions and NH4+ was adsorbed at the cation exchange sites. Different water contents induced highest CO2 emissions and therefore optimal conditions for microbial activity at 75% of WHC in both soils. Cumulative nitrification was also highest at 75% and 50% of WHC whereas the relative N2O emissions increased with water content and showed higher N2O losses in sandy soils. In summary it can be stated that the findings of the present thesis confirmed the high fertilizer value of BS’s, caused by high concentrations of NH4+ and labile organic compounds such as readily available carbon. These attributes of BS’s are to a great extent independent of the input materials of biogas plants. However, considerably gaseous and leaching losses of N may occur especially at high moisture contents. The emissions of N2O after field application corresponded with those of animal slurries.

Asymptotic model for shape resonance control of diatomics by intense non-resonant light

We derive a universal model for atom pairs interacting with non-resonant light via the polarizability anisotropy, based on the long range properties of the scattering. The corresponding dynamics can be obtained using a nodal line technique to solve the asymptotic Schrödinger equation. It consists of imposing physical boundary conditions at long range and vanishing the wavefunction at a position separating the inner zone and the asymptotic region. We show that nodal lines which depend on the intensity of the non-resonant light can satisfactorily account for the effect of the polarizability at short range. The approach allows to determine the resonance structure, energy, width, channel mixing and hybridization even for narrow resonances.

Asymptotic model for shape resonance control of diatomics by intense non-resonant light: universality in the single-channel approximation

Non-resonant light interacting with diatomics via the polarizability anisotropy couples different rotational states and may lead to strong hybridization of the motion. The modification of shape resonances and low-energy scattering states due to this interaction can be fully captured by an asymptotic model, based on the long-range properties of the scattering (Crubellier et al 2015 New J. Phys. 17 045020). Remarkably, the properties of the field-dressed shape resonances in this asymptotic multi-channel description are found to be approximately linear in the field intensity up to fairly large intensity. This suggests a perturbative single-channel approach to be sufficient to study the control of such resonances by the non-resonant field. The multi-channel results furthermore indicate the dependence on field intensity to present, at least approximately, universal characteristics. Here we combine the nodal line technique to solve the asymptotic Schrödinger equation with perturbation theory. Comparing our single channel results to those obtained with the full interaction potential, we find nodal lines depending only on the field-free scattering length of the diatom to yield an approximate but universal description of the field-dressed molecule, confirming universal behavior.

Energy saving innovations, non-exhaustible sources of energy and long-run: what would happen if we run out of oil?

We formulate and solve a model of factor saving technological improvement considering three factors of production: labor, capital and energy. The productive activities have three main characteristics: first, in order to use capital goods firms need energy; second, there are two sources of energy: non-exhaustible and exhaustible; third, capital goods can be of different qualities and the quality of these goods can be changed along two dimensions -reducing the need of energy or changing the source of energy used in the production process. The economy goes through three stages of development after industrialization. In the first, firms make use of exhaustible energy and the e¢ ciency in the use of energy is constant. In the second stage, as the price of energy grows the e¢ ciency in its use is increased. In the third stage, the price of exhaustible sources is so high that firms have incentives to use non-exhaustible sources of energy. During this stage the price of energy is constant. In this set up, the end of the oil age has level effects on consumption and output but it does not cause the collapse of the economic system.

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Referencing of street-level flows measured during the DAPPLE 2004 campaign

Street-level mean flow and turbulence govern the dispersion of gases away from their sources in urban areas. A suitable reference measurement in the driving flow above the urban canopy is needed to both understand and model complex street-level flow for pollutant dispersion or emergency response purposes. In vegetation canopies, a reference at mean canopy height is often used, but it is unclear whether this is suitable for urban canopies. This paper presents an evaluation of the quality of reference measurements at both roof-top (height = H) and at height z = 9H = 190 m, and their ability to explain mean and turbulent variations of street-level flow. Fast response wind data were measured at street canyon and reference sites during the six-week long DAPPLE project field campaign in spring 2004, in central London, UK, and an averaging time of 10 min was used to distinguish recirculation-type mean flow patterns from turbulence. Flow distortion at each reference site was assessed by considering turbulence intensity and streamline deflection. Then each reference was used as the dependent variable in the model of Dobre et al. (2005) which decomposes street-level flow into channelling and recirculating components. The high reference explained more of the variability of the mean flow. Coupling of turbulent kinetic energy was also stronger between street-level and the high reference flow rather than the roof-top. This coupling was weaker when overnight flow was stratified, and turbulence was suppressed at the high reference site. However, such events were rare (<1% of data) over the six-week long period. The potential usefulness of a centralised, high reference site in London was thus demonstrated with application to emergency response and air quality modelling.

An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance

This paper presents the model SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes), which is a vertical (1-D) integrated radiative transfer and energy balance model. The model links visible to thermal infrared radiance spectra (0.4 to 50 μm) as observed above the canopy to the fluxes of water, heat and carbon dioxide, as a function of vegetation structure, and the vertical profiles of temperature. Output of the model is the spectrum of outgoing radiation in the viewing direction and the turbulent heat fluxes, photosynthesis and chlorophyll fluorescence. A special routine is dedicated to the calculation of photosynthesis rate and chlorophyll fluorescence at the leaf level as a function of net radiation and leaf temperature. The fluorescence contributions from individual leaves are integrated over the canopy layer to calculate top-of-canopy fluorescence. The calculation of radiative transfer and the energy balance is fully integrated, allowing for feedback between leaf temperatures, leaf chlorophyll fluorescence and radiative fluxes. Leaf temperatures are calculated on the basis of energy balance closure. Model simulations were evaluated against observations reported in the literature and against data collected during field campaigns. These evaluations showed that SCOPE is able to reproduce realistic radiance spectra, directional radiance and energy balance fluxes. The model may be applied for the design of algorithms for the retrieval of evapotranspiration from optical and thermal earth observation data, for validation of existing methods to monitor vegetation functioning, to help interpret canopy fluorescence measurements, and to study the relationships between synoptic observations with diurnally integrated quantities. The model has been implemented in Matlab and has a modular design, thus allowing for great flexibility and scalability.

Vibrational energy levels of water

Several quartic force fields and a full sextic anharmonic force field for H,O have been determined from high-quality ab initio calculations, the highest at the aug-cc-pVQZ CCSD(T) level of theory. These force fields have been used to determine vibrational excited state band origins up to 15 000 cm - ’ above the zero-point level, using both a perturbation-resonancea pproach and a variational approach. An optimisedq uartic force field hasb eeno btained by least squares refinement of our best ab initio results to fit the observed overtone levels of 5 symmetrically substituted isotopomers of water (Hi60, Hi70, HisO, D,O, and T,O) with an rms error of less than 10 cm-‘, using the perturbation-resonancem odel for the vibrational calculation. Predicatel east squaresr efinement was usedt o provide a loose constraint of the refined force field to the ab initio results. The results obtained prove the viability of the perturbation-resonancem odel for usei n larger molecular systemsa nd also highlight someo f its weaknesse

A comparison of energy use in conventional and integrated arable farming systems in the UK

The LINK Integrated Farming Systems (LINK-IFS) Project (1992-1997) was setup to compare conventional and integrated arable farming systems (IAFS), concentrating on practical feasibility and economic viability, but also taking into account the level of inputs used and environmental impact. As part of this, an examination into energy use within the two systems was also undertaken. This paper presents the results from that analysis. The data used is from the six sites within the LINK-IFS Project, spread through the arable production areas of England and from the one site in Scotland, covering the 5 years of the project. The comparison of the energy used is based on the equipment and inputs used to produce I kg of each crop within the conventional and integrated rotations, and thereby the overall energy used for each system. The results suggest that, in terms of total energy used, the integrated system appears to be the most efficient. However, in terms of energy efficiency, energy use per kilogram of output, the results are less conclusive. (C) 2003 Elsevier Science B.V. All rights reserved.

Comparison of energy evaluation systems and a mechanistic model for milk production by dairy cattle offered fresh grass-based diets

Grass-based diets are of increasing social-economic importance in dairy cattle farming, but their low supply of glucogenic nutrients may limit the production of milk. Current evaluation systems that assess the energy supply and requirements are based on metabolisable energy (ME) or net energy (NE). These systems do not consider the characteristics of the energy delivering nutrients. In contrast, mechanistic models take into account the site of digestion, the type of nutrient absorbed and the type of nutrient required for production of milk constituents, and may therefore give a better prediction of supply and requirement of nutrients. The objective of the present study is to compare the ability of three energy evaluation systems, viz. the Dutch NE system, the agricultural and food research council (AFRC) ME system, and the feed into milk (FIM) ME system, and of a mechanistic model based on Dijkstra et al. [Simulation of digestion in cattle fed sugar cane: prediction of nutrient supply for milk production with locally available supplements. J. Agric. Sci., Cambridge 127, 247-60] and Mills et al. [A mechanistic model of whole-tract digestion and methanogenesis in the lactating dairy cow: model development, evaluation and application. J. Anim. Sci. 79, 1584-97] to predict the feed value of grass-based diets for milk production. The dataset for evaluation consists of 41 treatments of grass-based diets (at least 0.75 g ryegrass/g diet on DM basis). For each model, the predicted energy or nutrient supply, based on observed intake, was compared with predicted requirement based on observed performance. Assessment of the error of energy or nutrient supply relative to requirement is made by calculation of mean square prediction error (MSPE) and by concordance correlation coefficient (CCC). All energy evaluation systems predicted energy requirement to be lower (6-11%) than energy supply. The root MSPE (expressed as a proportion of the supply) was lowest for the mechanistic model (0.061), followed by the Dutch NE system (0.082), FIM ME system (0.097) and AFRCME system(0.118). For the energy evaluation systems, the error due to overall bias of prediction dominated the MSPE, whereas for the mechanistic model, proportionally 0.76 of MSPE was due to random variation. CCC analysis confirmed the higher accuracy and precision of the mechanistic model compared with energy evaluation systems. The error of prediction was positively related to grass protein content for the Dutch NE system, and was also positively related to grass DMI level for all models. In conclusion, current energy evaluation systems overestimate energy supply relative to energy requirement on grass-based diets for dairy cattle. The mechanistic model predicted glucogenic nutrients to limit performance of dairy cattle on grass-based diets, and proved to be more accurate and precise than the energy systems. The mechanistic model could be improved by allowing glucose maintenance and utilization requirements parameters to be variable. (C) 2007 Elsevier B.V. All rights reserved.

Effect of energy and protein supplementation on phosphorus utilization in lactating dairy cows

Two experiments were undertaken in which grass silage was used in conjunction with a series of different concentrate types designed to examine the effect of carbohydrate source, protein level and degradability on total dietary phosphorus (P) utilization with emphasis on P pollution. Twelve Holstein-Friesian dairy cows in early to mid-lactation were used in an incomplete changeover design with four periods consisting of 4 weeks each. Phosphorus intake ranged from 54 to 80 g/day and faecal P represented the principal route by which ingested P was disposed of by cows, with insignificant amounts being voided in urine. A positive linear relationship between faecal P and P intake was established. In Experiment 1, P utilization was affected by dietary carbohydrate type, with an associated output of 3.3 g faecal P/g milk P produced for all treatments except those utilizing low degradable starch and low protein supplements, where a mean value of 2.8 g faecal P/g milk P was observed. In Experiment 2, where two protein levels and three protein degradabilities were examined, the efficiency of P utilization for milk P production was not affected by either level or degradability of crude protein (CP) but a significant reduction in faecal P excretion due to lower protein and P intake was observed. In general, P utilization in Experiment 2 was substantially improved compared to the Experiment 1, with an associated output of 1.8 g faecal P/g milk P produced. The improved utilization of P in Experiment 2 could be due to lower P content of the diets offered and higher dry matter (DM) intake. For dairy cows weighing 600 kg, consuming 17-18 kg DM/day and producing about 25 kg milk, P excretion in faeces and hence P pollution to the environment might be minimized without compromising lactational performance by formulating diets to supply about 68 g P/day, which is close to recent published recommended requirements for P.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Vibrational energy levels for the electronic ground state of the diazocarbene (CNN) molecule

The vibrational energy levels of diazocarbene (diazomethylene) in its electronic ground state, (X) over tilde (3) Sigma(-) CNN, have been predicted using the variational method. The potential energy surfaces of (X) over tilde (3) A" CNN were determined by employing ab initio single reference coupled cluster with single and double excitations (CCSD), CCSD with perturbative triple excitations [CCSD(T)], multi-reference complete active space self-consistent-field (CASSCF), and internally contracted multi-reference configuration interaction (ICMRCI) methods. The correlation-consistent polarised valence quadruple zeta (cc-pVQZ) basis set was used. Four sets of vibrational energy levels determined from the four distinct analytical potential functions have been compared with the experimental values from the laser-induced fluorescence measurements of Wurfel et al. obtained in 1992. The CCSD, CCSD(T), and CASSCF potentials have not provided satisfactory agreement with the experimental observations. In this light, the importance of both non-dynamic (static) and dynamic correlation effects in describing the ground state of CNN is emphasised. Our best theoretical fundamental frequencies at the cc-pVQZ ICMRCI level of theory, v(1) = 1230, v(2) = 394, and v(3) = 1420 cm(-1) are in excellent agreement with the experimental values of v(1) = 1235, v(2) = 396, and v(3) = 1419cm(-1) and the mean absolute deviation between the 23 calculated and experimental vibrational energy levels is only 7.4 cm(-1). It is shown that the previously suggested observation of the v(3) frequency at about 2847cm(-1) was in fact the first overtone 2v(3).

A case of four-coordinate silver(I) adopting a high energy planar geometry. Experiment and theory

The ligands PhL and MeL are obtained by condensing 2-formylpyridine with benzil dihydrazone and diacetyl dihydrazone, respectively, in 2: 1 molar proportion. With silver( I), PhL yields a double-stranded dinuclear cationic helicate 1 in which the metal is tetrahedral but MeL gives a cationic one-dimensional polymeric complex 2 where silver( I) is distorted square planar and the ligand backbone is nearly planar. In both complexes, metal: ligand ratio is 1: 1. Ab initio calculations on the ligands at the HF/6-31+G* level reveal that while PhL strongly prefers a helical conformation, MeL has a natural inclination to remain in a planar conformation. Density functional theory calculations on model silver( I) complexes show that formation of the linear polymer in the case of MeL is also an important factor in imposing the planar geometry of Ag(I) in 2.

2,2-Dialkyl-2H-benzimidazoles, the high energy tautomers of the corresponding 1,2-dialkyl-1H-benzimidazoles. Syntheses and their complexes with Cu(I) and Ag(I)

Reaction of Cu(1,2-phenylenediamine)(2)(ClO4)(2) with neat RR'=O (R = methyl and/or ethyl) (lives Cu(2,2-dialkyl-2H-benzimidazole)ClO4. demetallation of which by the action of aqueous ammonia yields Pure 2,2-dialkyl-2H-benzimidazoles. These are characterised by NMR. hi the X-ray crystal Structure, Ag(2,2-methyl-2H-benzimi-dazolc)NO3 is Found to be a spiral 1D coordination polymer where the 2H-benzimidazole acts as an N,N bridge between two Ag(I) centus. Although 2H-benzimidazoles are very unstable in the free state, they are quite stable in their Cu(I)(1) and Ag(I) complexes. The 1,2-tautomerisation in imidazole and benzimidazole have been Studied by means of transition state calculations at B3LYP/6-3 11 +G(2d,p)* level.