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.

Holocene vegetation dynamics in the northeastern Rub' al-Khali desert, Arabian Peninsula: a phytolith, pollen and carbon isotope study

The Holocene vegetation history of the Arabian Peninsula is poorly understood, with few palaeobotanical studies to date. At Awafi, Ras al-Khaimah, UAE, a 3.3 m lake sediment sequence records the vegetation development for the period 8500 cal. yr BP to similar to3000 cal. yr BP. delta(13)C isotope, pollen and phytolith analyses indicate that C3 Pooid grassland with a strong woody element existed during the early Holocene (between 8500 and 6000 cal. yr BP) and became replaced by mixed C3 and C4 grasses with a strong C4 Panicoid tall grass element between 5900 and 5400 cal. yr BP. An intense, arid event Occurred at 4100 cal. yr BP when the lake desiccated and was infilled by Aeolian sand. From 4100 cal. yr BP the vegetation was dominated by C4 Chloridoid types and Cyperaceae, suggesting an incomplete vegetation cover and Aeolian dune reactivation owing to increased regional aridity. These data outline the ecosystem dynamics and carbon cycling in response to palaeomon-soon and north-westerly variability during the Holocene. Copyright (C) 2004 John Wiley Sons, Ltd.

Dynamics of Dissolved Forms of Carbon and Inorganic Nitrogen in Small Watersheds of the Coastal Atlantic Forest in Southeast Brazil

Based on the fact that streamwater quality reflects landscape conditions, the objectives of this study were: to investigate nitrogen (N), carbon (C), and major ion concentrations in six streams crossing minimally disturbed Atlantic Forest areas, with similar geomorphological characteristics; to determine N and C fluxes in one of these pristine streams (Indaia); and assess the impact of human activity on the biogeochemistry of two other streams in the same region, crossing urbanized areas. The distribution pattern of carbon and inorganic nitrogen dissolved forms, as well as the major ion and biogenic gas concentrations in the streamwater, was similar in pristine streams, indicating that the C and N dynamics were determined by influence of some factors, such as climate, atmospheric deposition, geology, soil type, and land covering, which were analogous in the forested watersheds. The urban streams were significantly different from the pristine streams, showing low dissolved oxygen concentrations, high respiration rates, and high concentrations of carbon dioxide, dissolved inorganic nitrogen, dissolved inorganic carbon, and major ion. These differences were attributed to anthropogenic impact on water quality, especially domestic sewage discharge. Additionally, in the Indaia stream, it was possible to observe the importance of rainfall over temporal dynamics of dissolved carbon forms, and also, the obtained specific flux of dissolved inorganic nitrogen was relatively elevated (approximately 11 kg ha(-1) year(-1)). These results reveal the influence of human activity over the biogeochemistry of coastal streams and also indicate the importance N export of Atlantic Forest to the ocean.

On the unzipping mechanisms of carbon nanotubes: Insights from reactive molecular dynamics simulations

Unzipping carbon nanotubes (CNTs) is considered one of the most promising approaches for the controlled and large-scale production of graphene nanoribbons (GNR). These structures are considered of great importance for the development of nanoelectronics because of its dimensions and intrinsic nonzero band gap value. Despite many years of investigations some details on the dynamics of the CNT fracture/unzipping processes remain unclear. In this work we have investigated some of these process through molecular dynamics simulations using reactive force fields (ReaxFF), as implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. We considered multi-walled CNTs of different dimensions and chiralities and under induced mechanical stretching. Our preliminary results show that the unzipping mechanisms are highly dependent on CNT chirality. Well-defined and distinct fracture patterns were observed for the different chiralities. Armchair CNTs favor the creation of GNRs with well-defined armchair edges, while zigzag and chiral ones produce GNRs with less defined and defective edges. © 2012 Materials Research Society.

Controlled route to the fabrication of carbon and boron nitride nanoscrolls: A molecular dynamics investigation

Carbon nanoscrolls (graphene layers rolled up into papyrus-like tubular structures) are nanostructures with unique and interesting characteristics that could be exploited to build several new nanodevices. However, an efficient and controlled synthesis of these structures was not achieved yet, making its large scale production a challenge to materials scientists. Also, the formation process and detailed mechanisms that occur during its synthesis are not completely known. In this work, using fully atomistic molecular dynamics simulations, we discuss a possible route to nanoscrolls made from graphene layers deposited over silicon oxide substrates containing chambers/pits. The scrolling mechanism is triggered by carbon nanotubes deposited on the layers. The process is completely general and can be used to produce scrolls from other lamellar materials, like boron nitride, for instance. © 2013 American Institute of Physics.

MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Influence of nitrogen and soil physical characteristics on belowground carbon flux dynamics of woody plants

At ecosystem level soil respiration (Rs) represents the largest carbon (C) flux after gross primary productivity, being mainly generated by root respiration (autotrophic respiration, Ra) and soil microbial respiration (heterotrophic respiration, Rh). In the case of terrestrial ecosystems, soils contain the largest C-pool, storing twice the amount of C contained in plant biomass. Soil organic matter (SOM), representing the main C storage in soil, is decomposed by soil microbial community. This process produces CO2 which is mainly released as Rh. It is thus relevant to understand how microbial activity is influenced by environmental factors like soil temperature, soil moisture and nutrient availability, since part of the CO2 produced by Rh, directly increases atmospheric CO2 concentration and therefore affects the phenomenon of climate change. Among terrestrial ecosystems, agricultural fields have traditionally been considered as sources of atmospheric CO2. In agricultural ecosystems, in particular apple orchards, I identified the role of root density, soil temperature, soil moisture and nitrogen (N) availability on Rs and on its two components, Ra and Rh. To do so I applied different techniques to separate Rs in its two components, the ”regression technique” and the “trenching technique”. I also studied the response of Ra to different levels of N availability, distributed either in a uniform or localized way, in the case of Populus tremuloides trees. The results showed that Rs is mainly driven by soil temperature, to which it is positively correlated, that high levels of soil moisture have inhibiting effects, and that N has a negligible influence on total Rs, as well as on Ra. Further I found a negative response of Rh to high N availability, suggesting that microbial decomposition processes in the soil are inhibited by the presence of N. The contribution of Ra to Rs was of 37% on average.

Soil application of Meliaceae derivatives: effect on carbon and nitrogen dynamics in the soil-plant system

The effect of soil incorporation of 7 Meliaceae derivatives (6 commercial neem cakes and leaves of Melia azedarach L.) on C and N dynamics and on nutrient availability to micropropagated GF677 rootstock was investigated. In a first laboratory incubation experiment the derivatives showed different N mineralization dynamics, generally well predicted by their C:N ratio and only partly by their initial N concentration. All derivatives increased microbial biomass C, thus representing a source of C for the soil microbial population. Soil addition of all neem cakes (8 g kg-1) and melia leaves (16 g kg-1) had a positive effect on plant growth and increased root N uptake and leaf green colour of micropropagated plants of GF677. In addition, the neem cakes characterized by higher nutrient concentration increased P and K concentration in shoot and leaves 68 days after the amendment. In another experiment, soil incorporation of 15N labeled melia leaves (16 g kg-1) had no effect on the total amount of plant N, however the percentage of melia derived-N of treated plants ranged between 0.8% and 34% during the experiment. At the end of the growing season, about 7% of N added as melia leaves was recovered in plant, while 70% of it was still present in soil. Real C mineralization and the priming effect induced by the addition of the derivatives were quantified by a natural 13C abundance method. The real C mineralization of the derivatives ranged between 22% and 40% of added-C. All the derivatives studied induced a positive priming effect and, 144 days after the amendment, the amount of C primed corresponded to 26% of added-C, for all the derivatives. Despite this substantial priming effect, the C balance of the soil, 144 days after the amendment, always resulted positive.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

INTERACTIVE CONTROLS OF WATER TABLE POSITION AND PLANT FUNCTIONAL TYPES ON PEAT POREWATER CHARACTER IN NORTHERN BOG ECOSYSTEMS: IMPLICATIONS FOR CARBON CYCLING DYNAMICS

Northern wetlands, and particularly peatlands, have been shown to store around 30% of the world's soil carbon and thus play a significant role in the carbon cycle of our planet. Changes in climate are altering peatland hydrology and vegetation communities. These changes are possibly resulting in declines in the ability of peatlands to sequester carbon because losses through carbon oxidation and mineralization are likely to increase relative to C inputs from net primary production in a warmer, drier climate. However, the consequences of interactive effects of altered hydrology and vegetation on carbon storage are not well understood. This research evaluated the importance of plant species, water table, and their interactive effects on porewater quality in a northern peatland with an average pH of 4.54, ranging from 4.15 to 4.8. We assessed the effects of plant functional group (ericaceous shrubs, sedges, and bryophytes) and water table position on biogeochemical processes. Specifically, we measured dissolved organic carbon (DOC), total dissolved nitrogen (TDN), potential enzyme activity, organic acids, anions and cations, spectral indexes of aromaticity, and phenolic content. Our results indicate that acetate and propionate concentrations in the sedge-dominated communities declined with depth and water table drawdown, relative to the control and ericaceous treatments. DOC increased in the lowered water table treatments in all vegetation community types, and the peat porewater C:N ratio declined in the sedge-dominated treatments when the water table was lowered. The relationship between DOC and ferrous iron showed significant responses to vegetation type; the exclusion of Ericaceae resulted in less ferrous iron per unit DOC compared to mixed species treatments and Ericaceae alone. This observation was corroborated with higher mean oxidation redox potential profiles (integrating 20, 40, and 70 cm) measured in the sedge treatments, compared with the mixed and Ericaceae species treatments over a growing season. Enzymatic activities did not show as strong of a response to treatments as expected; the oxidative enzyme peroxidase and the hydrolytic enzyme phosphatase were the only enzymes to respond to water table, where the potential activity of both enzymes increased with water table drawdown. Overall, there were significant interactive effects between changes in vegetation and water table position on peat porewater composition. These data suggest that vegetation effects on oxidation reduction potentials and peat porewater character can be as important as water table position in northern bog ecosystems.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.