Modelling coupled physical-biogeochemical processes in ice-covered oceans

The last decades have seen a large effort of the scientific community to study and understand the physics of sea ice. We currently have a wide - even though still not exhaustive - knowledge of the sea ice dynamics and thermodynamics and of their temporal and spatial variability. Sea ice biogeochemistry is instead largely unknown. Sea ice algae production may account for up to 25% of overall primary production in ice-covered waters of the Southern Ocean. However, the influence of physical factors, such as the location of ice formation, the role of snow cover and light availability on sea ice primary production is poorly understood. There are only sparse localized observations and little knowledge of the functioning of sea ice biogeochemistry at larger scales. Modelling becomes then an auxiliary tool to help qualifying and quantifying the role of sea ice biogeochemistry in the ocean dynamics. In this thesis, a novel approach is used for the modelling and coupling of sea ice biogeochemistry - and in particular its primary production - to sea ice physics. Previous attempts were based on the coupling of rather complex sea ice physical models to empirical or relatively simple biological or biogeochemical models. The focus is moved here to a more biologically-oriented point of view. A simple, however comprehensive, physical model of the sea ice thermodynamics (ESIM) was developed and coupled to a novel sea ice implementation (BFM-SI) of the Biogeochemical Flux Model (BFM). The BFM is a comprehensive model, largely used and validated in the open ocean environment and in regional seas. The physical model has been developed having in mind the biogeochemical properties of sea ice and the physical inputs required to model sea ice biogeochemistry. The central concept of the coupling is the modelling of the Biologically-Active-Layer (BAL), which is the time-varying fraction of sea ice that is continuously connected to the ocean via brines pockets and channels and it acts as rich habitat for many microorganisms. The physical model provides the key physical properties of the BAL (e.g., brines volume, temperature and salinity), and the BFM-SI simulates the physiological and ecological response of the biological community to the physical enviroment. The new biogeochemical model is also coupled to the pelagic BFM through the exchange of organic and inorganic matter at the boundaries between the two systems . This is done by computing the entrapment of matter and gases when sea ice grows and release to the ocean when sea ice melts to ensure mass conservation. The model was tested in different ice-covered regions of the world ocean to test the generality of the parameterizations. The focus was particularly on the regions of landfast ice, where primary production is generally large. The implementation of the BFM in sea ice and the coupling structure in General Circulation Models will add a new component to the latters (and in general to Earth System Models), which will be able to provide adequate estimate of the role and importance of sea ice biogeochemistry in the global carbon cycle.

External drivers of biogeochemical cycles in a tropical montane forest in Ecuador

Successful conservation of tropical montane forest, one of the most threatened ecosystems on earth, requires detailed knowledge of its biogeochemistry. Of particular interest is the response of the biogeochemical element cycles to external influences such as element deposition or climate change. Therefore the overall objective of my study was to contribute to improved understanding of role and functioning of the Andean tropical montane forest. In detail, my objectives were to determine (1) the role of long-range transported aerosols and their transport mechanisms, and (2) the role of short-term extreme climatic events for the element budget of Andean tropical forest. In a whole-catchment approach including three 8-13 ha microcatchments under tropical montane forest on the east-exposed slope of the eastern cordillera in the south Ecuadorian Andes at 1850-2200 m above sea level I monitored at least in weekly resolution the concentrations and fluxes of Ca, Mg, Na, K, NO3-N, NH4-N, DON, P, S, TOC, Mn, and Al in bulk deposition, throughfall, litter leachate, soil solution at the 0.15 and 0.3 m depths, and runoff between May 1998 and April 2003. I also used meteorological data from my study area collected by cooperating researchers and the Brazilian meteorological service (INPE), as well as remote sensing products of the North American and European space agencies NASA and ESA. My results show that (1) there was a strong interannual variation in deposition of Ca [4.4-29 kg ha-1 a-1], Mg [1.6-12], and K [9.8-30]) between 1998 and 2003. High deposition changed the Ca and Mg budgets of the catchments from loss to retention, suggesting that the additionally available Ca and Mg was used by the ecosystem. Increased base metal deposition was related to dust outbursts of the Sahara and an Amazonian precipitation pattern with trans-regional dry spells allowing for dust transport to the Andes. The increased base metal deposition coincided with a strong La Niña event in 1999/2000. There were also significantly elevated H+, N, and Mn depositions during the annual biomass burning period in the Amazon basin. Elevated H+ deposition during the biomass burning period caused elevated base metal loss from the canopy and the organic horizon and deteriorated already low base metal supply of the vegetation. Nitrogen was only retained during biomass burning but not during non-fire conditions when deposition was much smaller. Therefore biomass burning-related aerosol emissions in Amazonia seem large enough to substantially increase element deposition at the western rim of Amazonia. Particularly the related increase of acid deposition impoverishes already base-metal scarce ecosystems. As biomass burning is most intense during El Niño situations, a shortened ENSO cycle because of global warming likely enhances the acid deposition at my study forest. (2) Storm events causing near-surface water flow through C- and nutrient-rich topsoil during rainstorms were the major export pathway for C, N, Al, and Mn (contributing >50% to the total export of these elements). Near-surface flow also accounted for one third of total base metal export. This demonstrates that storm-event related near-surface flow markedly affects the cycling of many nutrients in steep tropical montane forests. Changes in the rainfall regime possibly associated with global climate change will therefore also change element export from the study forest. Element budgets of Andean tropical montane rain forest proved to be markedly affected by long-range transport of Saharan dust, biomass burning-related aerosols, or strong rainfalls during storm events. Thus, increased acid and nutrient deposition and the global climate change probably drive the tropical montane forest to another state with unknown consequences for its functions and biological diversity.

Copper and zinc stable isotope ratios as tracers of biogeochemical processes, sources and transport of Cu and Zn in soils

Copper and Zn are essential micronutrients for plants, animals, and humans; however, they may also be pollutants if they occur at high concentrations in soil. Therefore, knowledge of Cu and Zn cycling in soils is required both for guaranteeing proper nutrition and to control possible risks arising from pollution.rnThe overall objective of my study was to test if Cu and Zn stable isotope ratios can be used to investigate into the biogeochemistry, source and transport of these metals in soils. The use of stable isotope ratios might be especially suitable to trace long-term processes occurring during soil genesis and transport of pollutants through the soil. In detail, I aimed to answer the questions, whether (1) Cu stable isotopes are fractionated during complexation with humic acid, (2) 65Cu values can be a tracer for soil genetic processes in redoximorphic soils (3) 65Cu values can help to understand soil genetic processes under oxic weathering conditions, and (4) 65Cu and 66Zn values can act as tracers of sources and transport of Cu and Zn in polluted soils.rnTo answer these questions, I ran adsorption experiments at different pH values in the laboratory and modelled Cu adsorption to humic acid. Furthermore, eight soils were sampled representing different redox and weathering regimes of which two were influenced by stagnic water, two by groundwater, two by oxic weathering (Cambisols), and two by podzolation. In all horizons of these soils, I determined selected basic soil properties, partitioned Cu into seven operationally defined fractions and determined Cu concentrations and Cu isotope ratios (65Cu values). Finally, three additional soils were sampled along a deposition gradient at different distances to a Cu smelter in Slovakia and analyzed together with bedrock and waste material from the smelter for selected basic soil properties, Cu and Zn concentrations and 65Cu and 66Zn values.rnMy results demonstrated that (1) Copper was fractionated during adsorption on humic acid resulting in an isotope fractionation between the immobilized humic acid and the solution (65CuIHA-solution) of 0.26 ± 0.11‰ (2SD) and that the extent of fractionation was independent of pH and involved functional groups of the humic acid. (2) Soil genesis and plant cycling causes measurable Cu isotope fractionation in hydromorphic soils. The results suggested that an increasing number of redox cycles depleted 63Cu with increasing depth resulting in heavier 65Cu values. (3) Organic horizons usually had isotopically lighter Cu than mineral soils presumably because of the preferred uptake and recycling of 63Cu by plants. (4) In a strongly developed Podzol, eluviation zones had lighter and illuviation zones heavier 65Cu values because of the higher stability of organo-65Cu complexes compared to organo-63Cu complexes. In the Cambisols and a little developed Podzol, oxic weathering caused increasingly lighter 65Cu values with increasing depth, resulting in the opposite depth trend as in redoximorphic soils, because of the preferential vertical transport of 63Cu. (5) The 66Zn values were fractionated during the smelting process and isotopically light Zn was emitted allowing source identification of Zn pollution while 65Cu values were unaffected by the smelting and Cu emissions isotopically indistinguishable from soil. The 65Cu values in polluted soils became lighter down to a depth of 0.4 m indicating isotope fractionation during transport and a transport depth of 0.4 m in 60 years. 66Zn values had an opposite depth trend becoming heavier with depth because of fractionation by plant cycling, speciation changes, and mixing of native and smelter-derived Zn. rnCopper showed measurable isotope fractionation of approximately 1‰ in unpolluted soils, allowing to draw conclusions on plant cycling, transport, and redox processes occurring during soil genesis and 65Cu and 66Zn values in contaminated soils allow for conclusions on sources (in my study only possible for Zn), biogeochemical behavior, and depth of dislocation of Cu and Zn pollution in soil. I conclude that stable Cu and Zn isotope ratios are a suitable novel tool to trace long-term processes in soils which are difficult to assess otherwise.rn

Modelling biogeochemical cycles in forest ecosystems: a Bayesian approach

Forest models are tools for explaining and predicting the dynamics of forest ecosystems. They simulate forest behavior by integrating information on the underlying processes in trees, soil and atmosphere. Bayesian calibration is the application of probability theory to parameter estimation. It is a method, applicable to all models, that quantifies output uncertainty and identifies key parameters and variables. This study aims at testing the Bayesian procedure for calibration to different types of forest models, to evaluate their performances and the uncertainties associated with them. In particular,we aimed at 1) applying a Bayesian framework to calibrate forest models and test their performances in different biomes and different environmental conditions, 2) identifying and solve structure-related issues in simple models, and 3) identifying the advantages of additional information made available when calibrating forest models with a Bayesian approach. We applied the Bayesian framework to calibrate the Prelued model on eight Italian eddy-covariance sites in Chapter 2. The ability of Prelued to reproduce the estimated Gross Primary Productivity (GPP) was tested over contrasting natural vegetation types that represented a wide range of climatic and environmental conditions. The issues related to Prelued's multiplicative structure were the main topic of Chapter 3: several different MCMC-based procedures were applied within a Bayesian framework to calibrate the model, and their performances were compared. A more complex model was applied in Chapter 4, focusing on the application of the physiology-based model HYDRALL to the forest ecosystem of Lavarone (IT) to evaluate the importance of additional information in the calibration procedure and their impact on model performances, model uncertainties, and parameter estimation. Overall, the Bayesian technique proved to be an excellent and versatile tool to successfully calibrate forest models of different structure and complexity, on different kind and number of variables and with a different number of parameters involved.

Process-based modelling of lichens and bryophytes and their role in global biogeochemical cycles

This thesis presents a process-based modelling approach to quantify carbon uptake by lichens and bryophytes at the global scale. Based on the modelled carbon uptake, potential global rates of nitrogen fixation, phosphorus uptake and chemical weathering by the organisms are estimated. In this way, the significance of lichens and bryophytes for global biogeochemical cycles can be assessed. The model uses gridded climate data and key properties of the habitat (e.g. disturbance intervals) to predict processes which control net carbon uptake, namely photosynthesis, respiration, water uptake and evaporation. It relies on equations used in many dynamical vegetation models, which are combined with concepts specific to lichens and bryophytes, such as poikilohydry or the effect of water content on CO2 diffusivity. To incorporate the great functional variation of lichens and bryophytes at the global scale, the model parameters are characterised by broad ranges of possible values instead of a single, globally uniform value. The predicted terrestrial net uptake of 0.34 to 3.3 Gt / yr of carbon and global patterns of productivity are in accordance with empirically-derived estimates. Based on the simulated estimates of net carbon uptake, further impacts of lichens and bryophytes on biogeochemical cycles are quantified at the global scale. Thereby the focus is on three processes, namely nitrogen fixation, phosphorus uptake and chemical weathering. The presented estimates have the form of potential rates, which means that the amount of nitrogen and phosphorus is quantified which is needed by the organisms to build up biomass, also accounting for resorption and leaching of nutrients. Subsequently, the potential phosphorus uptake on bare ground is used to estimate chemical weathering by the organisms, assuming that they release weathering agents to obtain phosphorus. The predicted requirement for nitrogen ranges from 3.5 to 34 Tg / yr and for phosphorus it ranges from 0.46 to 4.6 Tg / yr. Estimates of chemical weathering are between 0.058 and 1.1 km³ / yr of rock. These values seem to have a realistic order of magnitude and they support the notion that lichens and bryophytes have the potential to play an important role for global biogeochemical cycles.

Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces.

Endocrine changes and liver mRNA abundance of somatotropic axis and insulin system constituents during negative energy balance at different stages of lactation in dairy cows

The liver has an important role in metabolic regulation and control of the somatotropic axis to adapt successfully to physiological and environmental changes in dairy cows. The aim of this study was to investigate the adaptation to negative energy balance (NEB) at parturition and to a deliberately induced NEB by feed restriction at 100 days in milk. The hepatic gene expression and the endocrine system of the somatotropic axis and related parameters were compared between the early and late NEB period. Fifty multiparous cows were subjected to 3 periods (1=early lactation up to 12 wk postpartum, 2=feed restriction for 3 wk beginning at around 100 days in milk with a feed-restricted and a control group, and 3=subsequent realimentation period for the feed-restricted group for 8 wk). In period 1, plasma growth hormone reached a maximum in early lactation, whereas insulin-like growth factor-I (IGF-I), leptin, the thyroid hormones, insulin, and the revised quantitative insulin sensitivity check index increased gradually after a nadir in early lactation. Three days after parturition, hepatic mRNA abundance of growth hormone receptor 1A, IGF-I, IGF-I receptor and IGF-binding protein-3 (IGFBP-3) were decreased, whereas mRNA of IGFBP-1 and -2 and insulin receptor were upregulated as compared with wk 3 antepartum. During period 2, feed-restricted cows showed decreased plasma concentrations of IGF-I and leptin compared with those of control cows. The revised quantitative insulin sensitivity check index was lower for feed-restricted cows (period 2) than for control cows. Compared with the NEB in period 1, the changes due to the deliberately induced NEB (period 2) in hormones were less pronounced. At the end of the 3-wk feed restriction, the mRNA abundance of IGF-I, IGFBP-1, -2, -3, and insulin receptor was increased as compared with the control group. The different effects of energy deficiency at the 2 stages in lactation show that the endocrine regulation changes qualitatively and quantitatively during the course of lactation.

Chemistry and Analysis of Phytocannabinoids and Other Cannabis Constituents

The Cannabis plant and its products consist of an enormous variety of chemicals. Some of the 483 compounds identified are unique to Cannabis, for example, the more than 60 cannabinoids, whereas the terpenes, with about 140 members forming the most abundant class, are widespread in the plant kingdom. The term “cannabinoids” [note: “ ” represents a group of C21 terpenophenolic compounds found until now uniquely in Cannabis sativa L. (1). As a consequence of the development of synthetic cannabinoids (e.g., nabilone [2], HU-211 [dexanabinol; ref. (3), or ajulemic acid [CT-3; ref. 4]) and the discovery of the chemically different endogenous cannabinoid receptor ligands (“endocannabinoids,” e.g., anandamide, 2-arachidonoylglycerol) (5,6), the term ’“phytocannabinoids’” was proposed for these particular Cannabis constituents (7).