Modelling the effects of nitrogen deposition and carbon dioxide enrichment on forest carbon balance

Atmospheric CO2 concentration ([CO2]) has increased over the last 250 years, mainly due to human activities. Of total anthropogenic emissions, almost 31% has been sequestered by the terrestrial biosphere. A considerable contribution to this sink comes from temperate and boreal forest ecosystems of the northern hemisphere, which contain a large amount of carbon (C) stored as biomass and soil organic matter. Several potential drivers for this forest C sequestration have been proposed, including increasing atmospheric [CO2], temperature, nitrogen (N) deposition and changes in management practices. However, it is not known which of these drivers are most important. The overall aim of this thesis project was to develop a simple ecosystem model which explicitly incorporates our best understanding of the mechanisms by which these drivers affect forest C storage, and to use this model to investigate the sensitivity of the forest ecosystem to these drivers. I firstly developed a version of the Generic Decomposition and Yield (G’DAY) model to explicitly investigate the mechanisms leading to forest C sequestration following N deposition. Specifically, I modified the G’DAY model to include advances in understanding of C allocation, canopy N uptake, and leaf trait relationships. I also incorporated a simple forest management practice subroutine. Secondly, I investigated the effect of CO2 fertilization on forest productivity with relation to the soil N availability feedback. I modified the model to allow it to simulate short-term responses of deciduous forests to environmental drivers, and applied it to data from a large-scale forest Free-Air CO2 Enrichment (FACE) experiment. Finally, I used the model to investigate the combined effects of recent observed changes in atmospheric [CO2], N deposition, and climate on a European forest stand. The model developed in my thesis project was an effective tool for analysis of effects of environmental drivers on forest ecosystem C storage. Key results from model simulations include: (i) N availability has a major role in forest ecosystem C sequestration; (ii) atmospheric N deposition is an important driver of N availability on short and long time-scales; (iii) rising temperature increases C storage by enhancing soil N availability and (iv) increasing [CO2] significantly affects forest growth and C storage only when N availability is not limiting.

Grazing incidence pumped Zr x-ray laser for spectroscopy on Li-like ions

X-ray laser fluorescence spectroscopy of the 2s-2p transition in Li-like ions is promising to become a widely applicable tool to provide information on the nuclear charge radii of stable and radioactive isotopes. For performing such experiments at the Experimental Storage Ring ESR, and the future NESR within the FAIR Project, a grazing incidence pumped (GRIP) x-ray laser (XRL) was set up at GSI Darmstadt using PHELIX (Petawatt High Energy Laser for heavy Ions eXperiments). The experiments demonstrated that lasing using the GRIP geometry could be achieved with relatively low pump energy, a prerequisite for higher repetition rate. In the first chapter the need of a plasma XRL is motivated and a short history of the plasma XRL is presented. The distinctive characteristic of the GRIP method is the controlled deposition of the pump laser energy into the desired plasma density region. While up to now the analysis performed were mostly concerned with the plasma density at the turning point of the main pump pulse, in this thesis it is demonstrated that also the energy deposition is significantly modified for the GRIP method, being sensitive in different ways to a large number of parameters. In the second chapter, the theoretical description of the plasma evolution, active medium and XRL emission properties are reviewed. In addition an innovative analysis of the laser absorption in plasma which includes an inverse Bremsstrahlung (IB) correction factor is presented. The third chapter gives an overview of the experimental set-up and diagnostics, providing an analytical formula for the average and instantaneous traveling wave speed generated with a tilted, on-axis spherical mirror, the only focusing system used up to now in GRIP XRL. The fourth chapter describes the experimental optimization and results. The emphasis is on the effect of the incidence angle of the main pump pulse on the absorption in plasma and on output and gain in different lasing lines. This is compared to the theoretical results for two different incidence angles. Significant corrections for the temperature evolution during the main pump pulse due to the incidence angle are demonstrated in comparison to a simple analytical model which does not take into account the pumping geometry. A much better agreement is reached by the model developed in this thesis. An interesting result is also the appearance of a central dip in the spatially resolved keV emission which was observed in the XRL experiments for the first time and correlates well with previous near field imaging and plasma density profile measurements. In the conclusion also an outlook to the generation of shorter wavelength XRL’s is given.

Noble-metal nanoparticles produced with colloidal lithography: fabrication, optical properties and applications

In this work, metal nanoparticles produced by nanosphere lithography were studied in terms of their optical properties (in connection to their plasmon resonances), their potential application in sensing platforms - for thin layer sensing and bio-recognition events -, and for a particular case (the nanocrescents), for enhanced spectroscopy studies. The general preparation procedures introduced early in 2005 by Shumaker-Parry et al. to produce metallic nanocrescents were extended to give rise to more complex (isolated) structures, and also, by combining colloidal monolayer fabrication and plasma etching techniques, to arrays of them. The fabrication methods presented in this work were extended not only to new shapes or arrangements of particles, but included also a targeted surface tailoring of the substrates and the structures, using different thiol and silane compounds as linkers for further attachment of, i.e. polyelectrolyte layers, which allow for a controlled tailoring of their nanoenvironment. The optical properties of the nanocrescents were studied with conventional transmission spectroscopy; a simple multipole model was adapted to explain their behaviour qualitatively. In terms of applications, the results on thin film sensing using these particles show that the crescents present an interesting mode-dependent sensitivity and spatial extension. Parallel to this, the penetrations depths were modeled with two simplified schemes, obtaining good agreement with theory. The multiple modes of the particles with their characteristic decay lengths and sensitivities represent a major improvement for particle-sensing platforms compared to previous single resonance systems. The nanocrescents were also used to alter the emission properties of fluorophores placed close to them. In this work, green emitting dyes were placed at controlled distances from the structures and excited using a pulsed laser emitting in the near infrared. The fluorescence signal obtained in this manner should be connected to a two-photon processes triggered by these structures; obtaining first insight into plasmon-mediated enhancement phenomena. An even simpler and faster approach to produce plasmonic structures than that for the crescents was tested. Metallic nanodiscs and nanoellipses were produced by means of nanosphere lithography, extending a procedure reported in the literature to new shapes and optical properties. The optical properties of these particles were characterized by extinction spectroscopy and compared to results from the literature. Their major advantage is that they present a polarization-dependent response, like the nanocrescents, but are much simpler to fabricate, and the resonances can be tailored in the visible with relative ease. The sensing capabilities of the metallic nanodiscs were explored in the same manner as for the nanocrescents, meaning their response to thin layers and to bio-recognition events on their surface. The sensitivity of these nanostructures to thin films proved to be lower than that of the crescents, though in the same order of magnitude. Experimental information about the near field extension for the Au nanodiscs of different sizes was also extracted from these measurements. Further resonance-tailoring approaches based on electrochemical deposition of metals on the nanodiscs were explored, as a means of modifying plasmon resonances by changing surface properties of the nanoparticles. First results on these experiments would indicate that the deposition of Ag on Au on a submonolayer coverage level can lead to important blue-shifts in the resonances, which would open a simple way to tailor resonances by changing material properties in a local manner.

Optical resonances of sphere-on-plane geometries

The optical resonances of metallic nanoparticles placed at nanometer distances from a metal plane were investigated. At certain wavelengths, these “sphere-on-plane” systems become resonant with the incident electromagnetic field and huge enhancements of the field are predicted localized in the small gaps created between the nanoparticle and the plane. An experimental architecture to fabricate sphere-on-plane systems was successfully achieved in which in addition to the commonly used alkanethiols, polyphenylene dendrimers were used as molecular spacers to separate the metallic nanoparticles from the metal planes. They allow for a defined nanoparticle-plane separation and some often are functionalized with a chromophore core which is therefore positioned exactly in the gap. The metal planes used in the system architecture consisted of evaporated thin films of gold or silver. Evaporated gold or silver films have a smooth interface with their substrate and a rougher top surface. To investigate the influence of surface roughness on the optical response of such a film, two gold films were prepared with a smooth and a rough side which were as similar as possible. Surface plasmons were excited in Kretschmann configuration both on the rough and on the smooth side. Their reflectivity could be well modeled by a single gold film for each individual measurement. The film has to be modeled as two layers with significantly different optical constants. The smooth side, although polycrystalline, had an optical response that was very similar to a monocrystalline surface while for the rough side the standard response of evaporated gold is retrieved. For investigations on thin non-absorbing dielectric films though, this heterogeneity introduces only a negligible error. To determine the resonant wavelength of the sphere-on-plane systems a strategy was developed which is based on multi-wavelength surface plasmon spectroscopy experiments in Kretschmann-configuration. The resonant behavior of the system lead to characteristic changes in the surface plasmon dispersion. A quantitative analysis was performed by calculating the polarisability per unit area /A treating the sphere-on-plane systems as an effective layer. This approach completely avoids the ambiguity in the determination of thickness and optical response of thin films in surface plasmon spectroscopy. Equal area densities of polarisable units yielded identical response irrespective of the thickness of the layer they are distributed in. The parameter range where the evaluation of surface plasmon data in terms of /A is applicable was determined for a typical experimental situation. It was shown that this analysis yields reasonable quantitative agreement with a simple theoretical model of the sphere-on-plane resonators and reproduces the results from standard extinction experiments having a higher information content and significantly increased signal-to-noise ratio. With the objective to acquire a better quantitative understanding of the dependence of the resonance wavelength on the geometry of the sphere-on-plane systems, different systems were fabricated in which the gold nanoparticle size, type of spacer and ambient medium were varied and the resonance wavelength of the system was determined. The gold nanoparticle radius was varied in the range from 10 nm to 80 nm. It could be shown that the polyphenylene dendrimers can be used as molecular spacers to fabricate systems which support gap resonances. The resonance wavelength of the systems could be tuned in the optical region between 550 nm and 800 nm. Based on a simple analytical model, a quantitative analysis was developed to relate the systems’ geometry with the resonant wavelength and surprisingly good agreement of this simple model with the experiment without any adjustable parameters was found. The key feature ascribed to sphere-on-plane systems is a very large electromagnetic field localized in volumes in the nanometer range. Experiments towards a quantitative understanding of the field enhancements taking place in the gap of the sphere-on-plane systems were done by monitoring the increase in fluorescence of a metal-supported monolayer of a dye-loaded dendrimer upon decoration of the surface with nanoparticles. The metal used (gold and silver), the colloid mean size and the surface roughness were varied. Large silver crystallites on evaporated silver surfaces lead to the most pronounced fluorescence enhancements in the order of 104. They constitute a very promising sample architecture for the study of field enhancements.

Charged colloidal suspensions in confined geometries

Im Zentrum dieser Arbeit steht das Verhalten von geladenen kolloidalen Suspensionen in eingeschränkten Geometrien. Es wurden verschiedene keilförmige Zellen verwendet, die eine kontinuierliche Variation der Abstände zwischen den Platten ermöglichen. In Zellen mit fluid geordneten Suspensionen bei niedrigen Salzkonzentrationen akkumulieren die kolloidalen Partikel in der Keilspitze und bilden kristallin geordnete Strukturen. Systematische Experimente zu diesem Akkumulationseffekt führen zu dem Schluss, dass es um eine elektrostatischer Fallensituation handeln muss, was durch ein einfaches theoretisches, von Löwen et al vorgeschlagenes Modell bestätigt wird. In Abhängig von der Zellhöhe lässt sich in den auftretenden kristallinen Strukturen eine charakteristische Abfolge erkennen. Diese Struktursequenz wurde schon zuvor in eingeschränkten Keilgeometrien beobachtet, jedoch ermöglichen die in unseren Experimenten realisierbaren kleinen Keilwinkel die Beobachtung neuer Strukturen. Einige dieser neuen Strukturen zeigen eine exotische Anordnung die keine atomare Entsprechung besitzen. Basierend auf experimentellen Beobachtungen schlagen wir Modelle für unterschiedliche Übergangsmechanismen zwischen den verschiedenen Strukturen vor, unter der physikalisch motivierten Vorraussetzung, dass sich die Partikel wie einem hohen Druck unterworfene harte Kugeln verhalten. Des Weiteren wurde eine Zelle mit variabler Höhe konstruiert, die zur Untersuchung des vollständigen Phasenverhaltens geladener, zwischen parallelen Platten eingeschlossener Kugeln dient. Die vorläufigen Ergebnisse werden mit theoretischen Prognosen verglichen.

Tunnelkontakte mit Heusler-Elektrode: Spinpolarisation von Co 2Cr 0,6Fe 0,4Al und Einfluss der Barriere

Ziel dieser Arbeit ist die Bestimmung der Spinpolarisation von der Heusler-Verbindung Co2Cr0,6Fe0,4Al. Dieses Ziel wurde durch die sorgfältige Präparation von Co2Cr0,6Fe0,4Al basierten Tunnelkontakten realisiert. Tunnelwiderstandsmessungen an Co2Cr0,6Fe0,4Al-basiertenrnTunnelkontakten ergaben einen Tunnelmagnetowiderstand von 101% bei 4 K. DieserrnTunnelmagnetowiderstand legt eine untere Grenze von 67% für die Spinpolarisation von Co2Cr0,6Fe0,4Al fest.rnrnCo2Cr0,6Fe0,4Al ist eine Heusler-Verbindung, der die Eigenschaften eines halbmetallischen Ferromagneten zugeschrieben werden. Ein halbmetallischer Ferromagnet hat an der Fermikante nur Elektronenspinzustände mit einer Polarisation. Als Folge davon können bei einem spinerhaltenden Tunnelprozess nur Elektronen einer Spinrichtung in den halbmetallischen Ferromagneten tunneln. Mit einem magnetischen Feld und einer durch einen Antiferromagneten fixierten Gegenelektrode, können an einem Tunnelkontakt mit einem spinpolarisierten Ferromagneten deshalb zwei Zustände, eine hohe und eine niedrige Tunnelleitfähigkeit, erzeugt werden. Daher finden spinpolarisierte Tunnelkontakte in Form von MRAM in der Datenspeicherung Verwendung. Bislang wurde jedoch keine Verbindung gefunden, der eine Spinpolarisation von 100% experimentell eindeutig nachgewiesen werden konnte. Für Co2Cr0,6Fe0,4Al lagen die höchsten gemessenen Spinpolarisationen um 50%.rnrnTunnelspektroskopie ist eine zuverlässige und anwendungsnahe Methode zur Untersuchung der Spinpolarisation. Inelastische Tunnelprozesse und eine reduzierte Ordnung an Grenzflächen bewirken einen reduzierten Tunnelmagnetowiderstand. Eine symmetriebrechende Barriere, wie amorphes AlOx, ist Voraussetzung für die Anwendung des Jullière-Modells zur Bestimmung der Spinpolarisation. Das Jullière-Modell verknüpft die Spin-aufgespaltenenrnZustandsdichten der Elektroden mit dem Tunnelmagnetowiderstand. Ohne einernsymmetriebrechende Barriere, zum Beispiel mit MgO als Isolatorschicht, können höhere Tunnelmagnetowiderstände erzwungen werden. Ein eindeutiger Rückschluss auf die Spinpolarisation ist dann jedoch nicht mehr möglich. Mit Aluminiumoxid-basierten Barrieren liefert die Anwendung des einfachen Jullière-Modells eine Untergrenze der Spinpolarisation.rnrnUm die Spinpolarisation von Co2Cr0,6Fe0,4Al durch Tunnelspektroskopie zu bestimmen, musste die Präparation der Tunnelkontakte verbessert werden. Dies wurde ermöglicht durch den Anbau einer neuen Sputterkammer mit besseren UHV-Bedingungen an ein bestehendes Präparationscluster. Co2Cr0,6Fe0,4Al wird mit Hilfe von Radiofrequenz-Kathodenzerstäuben deponiert. Die resultierenden Schichten verfügen nach ihrer Deposition über einen höheren Ordnungsgrad und über eine geordnete Oberfläche. Durch eine Magnesium-Pufferschicht war es möglich, auf diese Oberfläche eine homogene amorphe AlOx-Barriere zu deponieren. Als Gegenelektrode wurde CoFe als Ferromagnet mit MnFe als Antiferromagnet gewählt. Diese Gegenelektrode ermöglicht Tunnelmessungen bis hin zu Raumtemperatur.rnrnMit den in dieser Arbeit vorgestellten optimierten Analyse- und Präparationsmethoden ist es möglich, die Untergrenze der Spinpolarisation von Co2Cr0,6Fe0,4Al auf 67% anzuheben. Dies ist der bisher höchste veröffentlichte Wert der Spinpolarisation von Co2Cr0,6Fe0,4Al.rn

Measurement and modeling of cloud condensation nuclei in continental air

Atmospheric aerosol particles serving as cloud condensation nuclei (CCN) are key elements of the hydrological cycle and climate. Knowledge of the spatial and temporal distribution of CCN in the atmosphere is essential to understand and describe the effects of aerosols in meteorological models. In this study, CCN properties were measured in polluted and pristine air of different continental regions, and the results were parameterized for efficient prediction of CCN concentrations.The continuous-flow CCN counter used for size-resolved measurements of CCN efficiency spectra (activation curves) was calibrated with ammonium sulfate and sodium chloride aerosols for a wide range of water vapor supersaturations (S=0.068% to 1.27%). A comprehensive uncertainty analysis showed that the instrument calibration depends strongly on the applied particle generation techniques, Köhler model calculations, and water activity parameterizations (relative deviations in S up to 25%). Laboratory experiments and a comparison with other CCN instruments confirmed the high accuracy and precision of the calibration and measurement procedures developed and applied in this study.The mean CCN number concentrations (NCCN,S) observed in polluted mega-city air and biomass burning smoke (Beijing and Pearl River Delta, China) ranged from 1000 cm−3 at S=0.068% to 16 000 cm−3 at S=1.27%, which is about two orders of magnitude higher than in pristine air at remote continental sites (Swiss Alps, Amazonian rainforest). Effective average hygroscopicity parameters, κ, describing the influence of chemical composition on the CCN activity of aerosol particles were derived from the measurement data. They varied in the range of 0.3±0.2, were size-dependent, and could be parameterized as a function of organic and inorganic aerosol mass fraction. At low S (≤0.27%), substantial portions of externally mixed CCN-inactive particles with much lower hygroscopicity were observed in polluted air (fresh soot particles with κ≈0.01). Thus, the aerosol particle mixing state needs to be known for highly accurate predictions of NCCN,S. Nevertheless, the observed CCN number concentrations could be efficiently approximated using measured aerosol particle number size distributions and a simple κ-Köhler model with a single proxy for the effective average particle hygroscopicity. The relative deviations between observations and model predictions were on average less than 20% when a constant average value of κ=0.3 was used in conjunction with variable size distribution data. With a constant average size distribution, however, the deviations increased up to 100% and more. The measurement and model results demonstrate that the aerosol particle number and size are the major predictors for the variability of the CCN concentration in continental boundary layer air, followed by particle composition and hygroscopicity as relatively minor modulators. Depending on the required and applicable level of detail, the measurement results and parameterizations presented in this study can be directly implemented in detailed process models as well as in large-scale atmospheric and climate models for efficient description of the CCN activity of atmospheric aerosols.

Essays in Macroeconomics

The dissertation consists of four papers that aim at providing new contributions in the field of macroeconomics, monetary policy and financial stability. The first paper proposes a new Dynamic Stochastic General Equilibrium (DSGE) model with credit frictions and a banking sector to study the pro-cyclicality of credit and the role of different prudential regulatory frameworks in affecting business cycle fluctuations and in restoring macroeconomic and financial stability. The second paper develops a simple DSGE model capable of evaluating the effects of large purchases of treasuries by central banks. This theoretical framework is employed to evaluate the impact on yields and the macroeconomy of large purchases of medium- and long-term government bonds recently implemented in the US and UK. The third paper studies the effects of ECB communications about unconventional monetary policy operations on the perceived sovereign risk of Italy over the last five years. The empirical results are derived from both an event-study analysis and a GARCH model, which uses Italian long-term bond futures to disentangle expected from unexpected policy actions. The fourth paper proposes a DSGE model with an endogenous term structure of interest rates, which is able to replicate the stylized facts regarding the yield curve and the term premium in the US over the period 1987:3-2011:3, without compromising its ability to match macro dynamics.

A framework for cyber-physical system simulation

In these last years, systems engineering has became one of the major research domains. The complexity of systems has increased constantly and nowadays Cyber-Physical Systems (CPS) are a category of particular interest: these, are systems composed by a cyber part (computer-based algorithms) that monitor and control some physical processes. Their development and simulation are both complex due to the importance of the interaction between the cyber and the physical entities: there are a lot of models written in different languages that need to exchange information among each other. Normally people use an orchestrator that takes care of the simulation of the models and the exchange of informations. This orchestrator is developed manually and this is a tedious and long work. Our proposition is to achieve to generate the orchestrator automatically through the use of Co-Modeling, i.e. by modeling the coordination. Before achieving this ultimate goal, it is important to understand the mechanisms and de facto standards that could be used in a co-modeling framework. So, I studied the use of a technology employed for co-simulation in the industry: FMI. In order to better understand the FMI standard, I realized an automatic export, in the FMI format, of the models realized in an existing software for discrete modeling: TimeSquare. I also developed a simple physical model in the existing open source openmodelica tool. Later, I started to understand how works an orchestrator, developing a simple one: this will be useful in future to generate an orchestrator automatically.

Geometry-Aware Multiscale Image Registration via OBBTree-Based Polyaffine Log-Demons

Non-linear image registration is an important tool in many areas of image analysis. For instance, in morphometric studies of a population of brains, free-form deformations between images are analyzed to describe the structural anatomical variability. Such a simple deformation model is justified by the absence of an easy expressible prior about the shape changes. Applying the same algorithms used in brain imaging to orthopedic images might not be optimal due to the difference in the underlying prior on the inter-subject deformations. In particular, using an un-informed deformation prior often leads to local minima far from the expected solution. To improve robustness and promote anatomically meaningful deformations, we propose a locally affine and geometry-aware registration algorithm that automatically adapts to the data. We build upon the log-domain demons algorithm and introduce a new type of OBBTree-based regularization in the registration with a natural multiscale structure. The regularization model is composed of a hierarchy of locally affine transformations via their logarithms. Experiments on mandibles show improved accuracy and robustness when used to initialize the demons, and even similar performance by direct comparison to the demons, with a significantly lower degree of freedom. This closes the gap between polyaffine and non-rigid registration and opens new ways to statistically analyze the registration results.

THE RESONATE-AND-FIRE NEURON: TIME DEPENDENT AND FREQUENCY SELECTIVE NEURONS IN NEURAL NETWORKS

The means through which the nervous system perceives its environment is one of the most fascinating questions in contemporary science. Our endeavors to comprehend the principles of neural science provide an instance of how biological processes may inspire novel methods in mathematical modeling and engineering. The application ofmathematical models towards understanding neural signals and systems represents a vibrant field of research that has spanned over half a century. During this period, multiple approaches to neuronal modeling have been adopted, and each approach is adept at elucidating a specific aspect of nervous system function. Thus while bio-physical models have strived to comprehend the dynamics of actual physical processes occurring within a nerve cell, the phenomenological approach has conceived models that relate the ionic properties of nerve cells to transitions in neural activity. Further-more, the field of neural networks has endeavored to explore how distributed parallel processing systems may become capable of storing memory. Through this project, we strive to explore how some of the insights gained from biophysical neuronal modeling may be incorporated within the field of neural net-works. We specifically study the capabilities of a simple neural model, the Resonate-and-Fire (RAF) neuron, whose derivation is inspired by biophysical neural modeling. While reflecting further biological plausibility, the RAF neuron is also analytically tractable, and thus may be implemented within neural networks. In the following thesis, we provide a brief overview of the different approaches that have been adopted towards comprehending the properties of nerve cells, along with the framework under which our specific neuron model relates to the field of neuronal modeling. Subsequently, we explore some of the time-dependent neurocomputational capabilities of the RAF neuron, and we utilize the model to classify logic gates, and solve the classic XOR problem. Finally we explore how the resonate-and-fire neuron may be implemented within neural networks, and how such a network could be adapted through the temporal backpropagation algorithm.

Markers Models in Survival Analysis and Applications to Issues Associated with AIDS

Jewell and Kalbfleisch (1992) consider the use of marker processes for applications related to estimation of the survival distribution of time to failure. Marker processes were assumed to be stochastic processes that, at a given point in time, provide information about the current hazard and consequently on the remaining time to failure. Particular attention was paid to calculations based on a simple additive model for the relationship between the hazard function at time t and the history of the marker process up until time t. Specific applications to the analysis of AIDS data included the use of markers as surrogate responses for onset of AIDS with censored data and as predictors of the time elapsed since infection in prevalent individuals. Here we review recent work on the use of marker data to tackle these kinds of problems with AIDS data. The Poisson marker process with an additive model, introduced in Jewell and Kalbfleisch (1992) may be a useful "test" example for comparison of various procedures.

Dynamic modeling of active regeneration in catalyzed and non-catalyzed diesel particulate filters

Particulate matter (PM) emissions standards set by the US Environmental Protection Agency (EPA) have become increasingly stringent over the years. The EPA regulation for PM in heavy duty diesel engines has been reduced to 0.01 g/bhp-hr for the year 2010. Heavy duty diesel engines make use of an aftertreatment filtration device, the Diesel Particulate Filter (DPF). DPFs are highly efficient in filtering PM (known as soot) and are an integral part of 2010 heavy duty diesel aftertreatment system. PM is accumulated in the DPF as the exhaust gas flows through it. This PM needs to be removed by oxidation periodically for the efficient functioning of the filter. This oxidation process is also known as regeneration. There are 2 types of regeneration processes, namely active regeneration (oxidation of PM by external means) and passive oxidation (oxidation of PM by internal means). Active regeneration occurs typically in high temperature regions, about 500 - 600 °C, which is much higher than normal diesel exhaust temperatures. Thus, the exhaust temperature has to be raised with the help of external devices like a Diesel Oxidation Catalyst (DOC) or a fuel burner. The O2 oxidizes PM producing CO2 as oxidation product. In passive oxidation, one way of regeneration is by the use of NO2. NO2 oxidizes the PM producing NO and CO2 as oxidation products. The passive oxidation process occurs at lower temperatures (200 - 400 °C) in comparison to the active regeneration temperatures. Generally, DPF substrate walls are washcoated with catalyst material to speed up the rate of PM oxidation. The catalyst washcoat is observed to increase the rate of PM oxidation. The goal of this research is to develop a simple mathematical model to simulate the PM depletion during the active regeneration process in a DPF (catalyzed and non-catalyzed). A simple, zero-dimensional kinetic model was developed in MATLAB. Experimental data required for calibration was obtained by active regeneration experiments performed on PM loaded mini DPFs in an automated flow reactor. The DPFs were loaded with PM from the exhaust of a commercial heavy duty diesel engine. The model was calibrated to the data obtained from active regeneration experiments. Numerical gradient based optimization techniques were used to estimate the kinetic parameters of the model.

Study of particle movement in the nearshore region of Lake Superior with radioisotope tracers

Biogeochemical processes in the coastal region, including the coastal area of the Great Lakes, are of great importance due to the complex physical, chemical and biological characteristics that differ from those on either the adjoining land or open water systems. Particle-reactive radioisotopes, both naturally occurring (210Pb, 210Po and 7Be) and man-made (137Cs), have proven to be useful tracers for these processes in many systems. However, a systematic isotope study on the northwest coast of the Keweenaw Peninsula in Lake Superior has not yet been performed. In this dissertation research, field sampling, laboratory measurements and numerical modeling were conducted to understand the biogeochemistry of the radioisotope tracers and some particulate-related coastal processes. In the first part of the dissertation, radioisotope activities of 210Po and 210Pb in a variability of samples (dissolved, suspended particle, sediment trap materials, surficial sediment) were measured. A completed picture of the distribution and disequilibrium of this pair of isotopes was drawn. The application of a simple box model utilizing these field observations reveals short isotope residence times in the water column and a significant contribution of sediment resuspension (for both particles and isotopes). The results imply a highly dynamic coastal region. In the second part of this dissertation, this conclusion is examined further. Based on intensive sediment coring, the spatial distribution of isotope inventories (mainly 210Pb, 137Cs and 7Be) in the nearshore region was determined. Isotope-based focusing factors categorized most of the sampling sites as non- or temporary depositional zones. A twodimensional steady-state box-in-series model was developed and applied to individual transects with the 210Pb inventories as model input. The modeling framework included both water column and upper sediments down to the depth of unsupported 210Pb penetration. The model was used to predict isotope residence times and cross-margin fluxes of sediments and isotopes at different locations along each transect. The time scale for sediment focusing from the nearshore to offshore regions of the transect was on the order of 10 years. The possibility of sediment longshore movement was indicated by high inventory ratios of 137Cs: 210Pb. Local deposition of fine particles, including fresh organic carbon, may explain the observed distribution of benthic organisms such as Diporeia. In the last part of this dissertation, isotope tracers, 210Pb and 210Po, were coupled into a hydrodynamic model for Lake Superior. The model was modified from an existing 2-D finite difference physical-biological model which has previously been successfully applied on Lake Superior. Using the field results from part one of this dissertation as initial conditions, the model was used to predict the isotope distribution in the water column; reasonable results were achieved. The modeling experiments demonstrated the potential for using a hydrodynamic model to study radioisotope biogeochemistry in the lake, although further refinements are necessary.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.